Abstract
Time integration of Fourier pseudo-spectral DNS is usually performed using the classical fourth-order accurate Runge–Kutta method, or other methods of second or third order, with a fixed step size. We investigate the use of higher-order Runge–Kutta pairs and automatic step size control based on local error estimation. We find that the fifth-order accurate Runge–Kutta pair of Bogacki & Shampine gives much greater accuracy at a significantly reduced computational cost. Specifically, we demonstrate speedups of 2x-10x for the same accuracy. Numerical tests (including the Taylor–Green vortex, Rayleigh–Taylor instability, and homogeneous isotropic turbulence) confirm the reliability and efficiency of the method. We also show that adaptive time stepping provides a significant computational advantage for some problems (like the development of a Rayleigh–Taylor instability) without compromising accuracy.
| Original language | English (US) |
|---|---|
| Pages (from-to) | 79-93 |
| Number of pages | 15 |
| Journal | International Journal for Numerical Methods in Fluids |
| Volume | 92 |
| Issue number | 2 |
| DOIs | |
| State | Published - Oct 15 2019 |
Bibliographical note
KAUST Repository Item: Exported on 2020-10-01Acknowledgements: This research used the resources of the Supercomputing Laboratory and Extreme Computing Research Center at the King Abdullah University of Science & Technology (KAUST) in Thuwal, Saudi Arabia. N.S. was supported by the KAUST Visiting Student Research Program. N.S. alsoacknowledges support from the Priority Programme SPP1881 Turbulent Superstructures of theDeutsche Forschungsgemeinschaft. M. M. acknowledges support from the 4DSpace StrategicResearch Initiative at the University of Oslo.