A further study of numerical errors in large-eddy simulations
Numerical errors in large-eddy simulations (LES) arise from aliasing and discretization errors, and errors in the subfilter-scale (SFS) turbulence model. Using a direct numerical simulation (DNS) dataset of stably stratified shear flow to perform a priori tests, we compare the numerical error from several finite difference schemes to the magnitude of the SFS force. This is an extension of Ghosals analysis [J. Comput. Phys. 125 (1996) 187] to realistic flow fields. By evaluating different grid resolutions as well as different filtergrid ratios, we provide guidelines for LES: for a second-order finite difference scheme, a filtergrid ratio of at least four is desired; for a sixth-order Pad scheme, a filtergrid ratio of two is sufficient.
large-eddy simulation les turbulence numerical error analysis