Abstract
High-order unstructured methods have become a popular choice for the simulation of complex unsteady flows. Flux reconstruction (FR) is a high-order spatial discretization method, which has been found to be particularly accurate for scale-resolving simulations of complex phenomena. In addition, it has been shown to provide sufficient dissipation for implicit large-eddy simulation (ILES). In conjunction with an FR discretization, an appropriate temporal scheme must be chosen. A common choice is explicit schemes due to their efficiency and ease of implementation. However, these methods usually require a small time-step size to remain stable. Recently, the development of optimal explicit Runge–Kutta (OERK) schemes has enabled stable simulations with larger time-step sizes. Hence, we analyze the fully-discrete properties of the FR method with OERK temporal schemes. We show results for first, second, third, fourth and eighth-order OERK schemes. We observe that OERK schemes modify the spectral behaviour of the semidiscretization. In particular, dissipation decreases in the region of high wavenumbers. We observe that higher-order OERK schemes require a smaller time step than the low-order schemes. However, they follow the dispersion relations of the FR scheme for a larger range of wavenumbers. We validate our analysis with simple advection test cases. It was observed that first and second-degree temporal schemes introduce a relatively large amount of error in the solutions. A one-dimensional ILES test case showed that, as long as the time-step size is not in the vicinity of the stability limit, results are generally similar to classical RK schemes.
Similar content being viewed by others
References
Huynh, H.T.: A flux reconstruction approach to high-order schemes including discontinuous Galerkin methods. In: 18th AIAA Computational Fluid Dynamics Conference, p. 4079 (2007)
Cockburn, B., Karniadakis, G.E., Shu, C.-W.: The development of discontinuous Galerkin methods. In: Cockburn, B., Karniadakis, G.E., Shu, C.-W. (eds.) Discontinuous Galerkin Methods. Lecture Notes in Computational Science and Engineering, pp. 3–50. Springer, Berlin (2000)
Cockburn, B., Lin, S.-Y., Shu, C.-W.: TBV Runge–Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems. J. Comput. Phys. 84(1), 90–113 (1989)
Wang, Z.J.: Spectral (finite) volume method for conservation laws on unstructured grids. Basic formulation: basic formulation. J. Comput. Phys. 178(1), 210–251 (2002)
Liu, Y., Vinokur, M., Wang, Z.J.: Discontinuous spectral difference method for conservation laws on unstructured grids. In: Computational Fluid Dynamics, pp. 449–454. Springer (2006)
Tsoutsanis, P., Antoniadis, A.F., Drikakis, D.: WENO schemes on arbitrary unstructured meshes for laminar, transitional and turbulent flows. J. Comput. Phys. 256, 254–276 (2014)
Ritos, K., Kokkinakis, I.W., Drikakis, D.: Physical insight into the accuracy of finely-resolved iLES in turbulent boundary layers. Comput. Fluids 169, 309–316 (2018)
Nogueira, X., Ramírez, L., Fernández-Fidalgo, J., Deligant, M., Khelladi, S., Chassaing, J.-C., Navarrina, F.: An a posteriori-implicit turbulent model with automatic dissipation adjustment for large Eddy simulation of compressible flows. Comput. Fluids 197, 104371 (2020)
Maulik, R., San, O., Behera, R.: An adaptive multilevel wavelet framework for scale-selective WENO reconstruction schemes. Int. J. Numer. Methods Fluids 87(5), 239–269 (2018)
Vermeire, B., Cagnone, J.-S., Nadarajah, S.: ILES using the correction procedure via reconstruction scheme. In: 51st AIAA aerospace sciences meeting including the new horizons forum and aerospace exposition, p. 1001 (2013)
Moura, R.C., Sherwin, S.J., Peiró, J.: Linear dispersion-diffusion analysis and its application to under-resolved turbulence simulations using discontinuous Galerkin spectral/HP methods. J. Comput. Phys. 298, 695–710 (2015)
Van den Abeele, K.: Development of high-order accurate schemes for unstructured grids. Ph.D. thesis in Vrije Universiteit Brussel (2009)
Vincent, P.E., Castonguay, P., Jameson, A.: Insights from von Neumann analysis of high-order flux reconstruction schemes. J. Comput. Phys. 230(22), 8134–8154 (2011)
Moura, R.C., Peiró, J., Sherwin, S.J.: Implicit LES approaches via discontinuous Galerkin methods at very large Reynolds. In: Direct and Large-Eddy Simulation XI, pp. 53–59. Springer (2019)
Stetter, H.J.: Analysis of Discretization Methods for Ordinary Differential Equations. Springer Tracts in Natural Philosophy. Springer, Berlin (1973)
Hairer, E., Nørsett, S.P., Wanner, G.: Solving Ordinary Differential Equations I: Nonstiff Problems. Springer, Berlin (1993)
Ketcheson, D.I., Ahmadia, A.J.: Optimal stability polynomials for numerical integration of initial value problems. Commun. Appl. Math. Comput. Sci. 7(2), 247–271 (2012)
Kubatko, E.J., Yeager, B.A., Ketcheson, D.I.: Optimal strong-stability-preserving Runge–Kutta time discretizations for discontinuous Galerkin methods. J. Sci. Comput. 60(2), 313–344 (2014)
Vermeire, B.C.: Paired explicit Runge–Kutta schemes for stiff systems of equations. J. Comput. Phys. 393, 465–483 (2019)
Yang, H., Li, F., Qiu, J.: Dispersion and dissipation errors of two fully discrete discontinuous Galerkin methods. J. Sci. Comput. 55, 552–572 (2013)
Vermeire, B.C., Vincent, P.E.: On the behaviour of fully-discrete flux reconstruction schemes. Comput. Methods Appl. Mech. Eng. 315, 1053–1079 (2017)
Castonguay, P., Vincent, P.E., Jameson, A.: A new class of high-order energy stable flux reconstruction schemes for triangular elements. J. Sci. Comput. 51(1), 224–256 (2012)
Diamond, S., Boyd, S.: CVXPY: a python-embedded modeling language for convex optimization. J. Mach. Learn. Res. 17(83), 1–5 (2016)
Agrawal, A., Verschueren, R., Diamond, S., Boyd, S.: A rewriting system for convex optimization problems. J. Control Decis. 5(1), 42–60 (2018)
Horn, R.A., Johnson, C.R.: Matrix Analysis. Cambridge University Press, Cambridge (2012)
Guo, W., Zhong, X., Qiu, J.-M.: Superconvergence of discontinuous Galerkin and local discontinuous Galerkin methods: eigen-structure analysis based on Fourier approach. J. Comput. Phys. 235, 458–485 (2013)
Cockburn, B., Shu, C.-W.: Runge–Kutta Discontinuous Galerkin methods for convection-dominated problems. J. Sci. Comput. 16(3), 173–261 (2001)
Hu, F.Q., Atkins, H.L.: Eigensolution analysis of the discontinuous Galerkin method with nonuniform grids: I. One space dimension. J. Comput. Phys. 182(2), 516–545 (2002)
Adams, N.A.: Implicit subgrid-scale modeling by adaptive local deconvolution. In: Friedrich, R., Geurts, B.J., Métais, O. (eds.) Direct and Large-Eddy Simulation V. ERCOFTAC Series, pp. 15–22. Springer, Dordrecht (2004)
Wang, Z.J., Li, Y.: A mathematical analysis of scale similarity. Commun. Comput. Phys. 21(1), 149–161 (2017)
Li, Y., Wang, Z.J.: A priori and a posteriori evaluations of subgrid stress models with the Burgers’ equation. In: 53rd AIAA Aerospace Sciences Meeting. American Institute of Aeronautics and Astronautics (2015)
Hesthaven, J.S., Warburton, T.: Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications. Springer, Berlin (2007)
Acknowledgements
We acknowledge the support of the Natural Sciences and Engineering Research Council of Canada (NSERC), [RGPAS-2017-507988, RGPIN-2017-06773]. This research was enabled in part by support provided by Calcul Quebec (www.calculquebec.ca), WestGrid (www.westgrid.ca), SciNet (www.scinethpc.ca), and Compute Canada (www.computecanada.ca) via a Resources for Research Groups allocation. We would also like to thank the anonymous reviewers whose comments improved the quality of this manuscript.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Data Availability Statement
Data relating to the results in this manuscript can be downloaded from the publication’s website under a CC BY-NC-ND 4.0 license.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Pereira, C.A., Vermeire, B.C. Fully-Discrete Analysis of High-Order Spatial Discretizations with Optimal Explicit Runge–Kutta Methods. J Sci Comput 83, 63 (2020). https://doi.org/10.1007/s10915-020-01243-8
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-020-01243-8