Turbulence spectra characteristics of high order schemes for direct and large eddy simulation

Abstract The comparative resolution of the high wavenumber portion of compressible turbulence energy spectrum by some high order numerical schemes is presented in this paper. Included in this are the essentially nonoscillatory (ENO) schemes, the weighted essentially nonoscillatory schemes (WENO), and the compact differencing schemes. The governing equations are the Navier–Stokes equations and the objective is to identify the numerical scheme that best represents the physics of compressible turbulence. Mach numbers M 1 values of 0.1, 0.5 and 0.7 are studied. The compact differencing schemes need filters for numerical stability. It is found in this work that a parameter in the filter scheme provides some flexibility for controlling the physical turbulence energy transfer rate at high wavenumbers, vis-a-vis the numerical dissipation at those scales. Although the ENO schemes do not require filters for numerical stability, the present study shows that the addition of filters improves the energy transfer process at high wavenumbers. Without filtering, with relatively coarse grids, numerical turbulence caused by stencil adaptation persists. This limits the useful wavenumber resolution range of the ENO schemes. The WENO schemes do not require the stabilizing filters but the results tend to be slightly more dissipative. Finally, at low Mach numbers, the current compact differencing and filter scheme formulation gives better results but as the Mach number increases the relative suitability of the ENO method increases.

[1]  Joe F. Thompson,et al.  Numerical grid generation , 1985 .

[2]  Miguel R. Visbal,et al.  High-Order Schemes for Navier-Stokes Equations: Algorithm and Implementation Into FDL3DI , 1998 .

[3]  Xiaodan Cai,et al.  Uniform mean scalar gradient in grid turbulence: Asymptotic probability distribution of a passive scalar , 1996 .

[4]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[5]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[6]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[7]  Xiaodan Cai,et al.  Thermodynamic behavior in decaying, compressible turbulence with initially dominant temperature fluctuations , 1996 .

[8]  D. Gaitonde,et al.  Practical aspects of high-order accurate finite-volume schemes for electromagnetics , 1997 .

[9]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[10]  D. Fyfe,et al.  Economical evaluation of Runge-Kutta formulae , 1966 .

[11]  F. Ladeinde,et al.  A parallelized eno procedure for direct numerical simulation of compressible turbulence , 1996 .

[12]  Eckart Meiburg,et al.  A numerical study of the convergence properties of ENO schemes , 1990 .

[13]  W. Matthaeus,et al.  Low Mach number two-dimensional hydrodynamic turbulence - Energy budgets and density fluctuations in a polytropic fluid , 1992 .

[14]  S. Osher,et al.  Some results on uniformly high-order accurate essentially nonoscillatory schemes , 1986 .

[15]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[16]  Datta Gaitonde,et al.  High-order accurate methods for unsteady vortical flows on curvilinear meshes , 1998 .

[17]  Chi-Wang Shu.  Numerical experiments on the accuracy of ENO and modified ENO schemes , 1990 .

[18]  Xiaodan Cai,et al.  Advection by polytropic compressible turbulence , 1995 .

[19]  Gordon Erlebacher,et al.  High-order ENO schemes applied to two- and three-dimensional compressible flow , 1992 .

[20]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[21]  T. Passot,et al.  Numerical simulation of compressible homogeneous flows in the turbulent regime , 1987, Journal of Fluid Mechanics.