On the Performances of Standard and Kinetic Energy Preserving Time-Integration Methods for Incompressible-Flow Simulations

The effects of kinetic-energy preservation errors due to Runge–Kutta (RK) temporal integrators have been analyzed for the case of large-eddy simulations of incompressible turbulent channel flow. Simulations have been run using the open-source solver Xcompact3D with an implicit spectral vanishing viscosity model and a variety of temporal Runge–Kutta integrators. Explicit pseudo-symplectic schemes, with improved energy preservation properties, have been compared to standard RK methods. The results show a marked decrease in the temporal error for higher-order pseudo-symplectic methods; on the other hand, an analysis of the energy spectra indicates that the dissipation introduced by the commonly used three-stage RK scheme can lead to significant distortion of the energy distribution within the inertial range. A cost-vs-accuracy analysis suggests that pseudo-symplectic schemes could be used to attain results comparable to traditional methods at a reduced computational cost.