Optimized Runge-Kutta Methods with Automatic Step Size Control for Compressible Computational Fluid Dynamics
arXiv:2104.06836 · doi:10.1007/s42967-021-00159-w
Abstract
We develop error-control based time integration algorithms for compressible fluid dynamics (CFD) applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime. Focusing on discontinuous spectral element semidiscretizations, we design new controllers for existing methods and for some new embedded Runge-Kutta pairs. We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice. We compare a wide range of error-control-based methods, along with the common approach in which step size control is based on the Courant-Friedrichs-Lewy (CFL) number. The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances, while additionally providing control of the temporal error at tighter tolerances. The numerical examples include challenging industrial CFD applications.
References in corpus (8)
- Split Form Nodal Discontinuous Galerkin Schemes with Summation-By-Parts Property for the Compressible Euler Equations
- Relaxation Runge-Kutta Methods: Fully-Discrete Explicit Entropy-Stable Schemes for the Compressible Euler and Navier-Stokes Equations
- Entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations
- General Relaxation Methods for Initial-Value Problems with Application to Multistep Schemes
- Optimal stability polynomials for numerical integration of initial value problems
- Fully-Discrete Explicit Locally Entropy-Stable Schemes for the Compressible Euler and Navier-Stokes Equations
- On the robustness and performance of entropy stable discontinuous collocation methods for the compressible Navier-Stokes equations
- Elastic wave propagation in anisotropic solids using energy-stable finite differences with weakly enforced boundary and interface conditions