Application of approximate dispersion-diffusion analyses to under-resolved Burgers turbulence using high resolution WENO and UWC schemes

• Published in: Journal of computational physics 2021-06, Vol.435, p.110246, Article 110246
• Main Authors: Solán-Fustero, P., Navas-Montilla, A., Ferrer, E., Manzanero, J., García-Navarro, P.
• Format: Article
• Language: English
• Subjects: Burgers' turbulence; Computational fluid dynamics; Computational physics; Diffusion; Dispersion-diffusion analysis; High-order schemes; Implicit Large Eddy Simulation; Integrators; Numerical dissipation; Runge-Kutta method; Turbulence; Von Neumann; Weighted essentially non-oscillatory WENO
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2021.110246

This paper presents a space-time approximate diffusion-dispersion analysis of high-order, finite volume Upwind Central (UWC) and Weighted Essentially Non-Oscillatory (WENO) schemes. We perform a thorough study of the numerical errors to find a-priori guidelines for the computation of under-resolved turbulent flows. In particular, we study the 3-rd, 5-th and 7-th order UWC and WENO reconstructions in space, and 3-rd and 4-th order Runge-Kutta time integrators. To do so, we use the approximate von Neumann analysis for non-linear schemes introduced by Pirozzoli. Moreover, we apply the “1% rule” for the dispersion-diffusion curves proposed by Moura et al. [41] to determine the range of wavenumbers that are accurately resolved by each scheme. The dispersion-diffusion errors estimated from these analyses agree with the numerical results for the forced Burgers' turbulence problem, which we use as a benchmark. The cut-off wavenumbers defined by the “1% rule” are evidenced to serve as a good estimator of the beginning of the dissipation region of the energy cascade and they are shown to be associated to a similar level of dissipation, with independence of the scheme. Finally, we show that WENO schemes are more diffusive than UWC schemes, leading to stable simulations at the price of more dissipative results. It is concluded both UWC and WENO schemes may be suitable schemes for iLES turbulence modeling, given their numerical dissipation level acting at the appropriate wavenumbers. •A space-time dispersion-diffusion analysis of UWC-RK and WENO-RK is presented.•The “1% rule” determines the wavenumbers that are accurately resolved.•Cut-off wavenumbers are obtained for each scheme and different CFL numbers.•The results are supported by the forced Burgers turbulence problem.•UWC and WENO may be suitable for iLES, being WENO more dissipative than UWC.