TY - JOUR
T1 - Assessment of localized artificial diffusivity scheme for large-eddy simulation of compressible turbulent flows
AU - Kawai, Soshi
AU - Shankar, Santhosh K.
AU - Lele, Sanjiva K.
N1 - Funding Information:
This work is supported by AFOSR-MURI (Grant FA9550-04-1-0387 ) and DOE-SciDAC (Grant DE-FC02-06-ER25787 ) programs. The first author also acknowledges financial support from the Basic Science Research Program of the Sumitomo Foundation. The present code is based on the extension to the code FDL3DI provided by Dr. M.R. Visbal, whom the authors thank for this. We gratefully acknowledge Dr. J. Larsson for providing the DNS and WENO data discussed in Section 3.2 and Prof. S. Pirozzoli for providing their DNS data in Section 3.3 . We are grateful to Dr. A. Mani for comments which led to improvements in the grid-dependent scaling of the artificial diffusivity. We are also grateful to our SciDAC project collaborators, and in particular Dr. J. Larsson, and Dr. A.W. Cook for various discussions on numerical issues for shock–turbulence interaction problems. We wish to thank the referees for their comments which helped to improve the paper.
PY - 2010/3/1
Y1 - 2010/3/1
N2 - The localized artificial diffusivity method is investigated in the context of large-eddy simulation of compressible turbulent flows. The performance of different artificial bulk viscosity models are evaluated through detailed results from the evolution of decaying compressible isotropic turbulence with eddy shocklets and supersonic turbulent boundary layer. Effects of subgrid-scale (SGS) models and implicit time-integration scheme/time-step size are also investigated within the framework of the numerical scheme used. The use of a shock sensor along with artificial bulk viscosity significantly improves the scheme for simulating turbulent flows involving shocks while retaining the shock-capturing capability. The proposed combination of Ducros-type sensor with a negative dilatation sensor removes unnecessary bulk viscosity within expansion and weakly compressible turbulence regions without shocks and allows it to localize near the shocks. It also eliminates the need for a wall-damping function for the bulk viscosity while simulating wall-bounded turbulent flows. For the numerical schemes used, better results are obtained without adding an explicit SGS model than with SGS model at moderate Reynolds number. Inclusion of a SGS model in addition to the low-pass filtering and artificial bulk viscosity results in additional damping of the resolved turbulence. However, investigations at higher Reynolds numbers suggest the need for an explicit SGS model. The flow statistics obtained using the second-order implicit time-integration scheme with three sub-iterations closely agrees with the explicit scheme if the maximum Courant-Friedrichs-Lewy is kept near unity.
AB - The localized artificial diffusivity method is investigated in the context of large-eddy simulation of compressible turbulent flows. The performance of different artificial bulk viscosity models are evaluated through detailed results from the evolution of decaying compressible isotropic turbulence with eddy shocklets and supersonic turbulent boundary layer. Effects of subgrid-scale (SGS) models and implicit time-integration scheme/time-step size are also investigated within the framework of the numerical scheme used. The use of a shock sensor along with artificial bulk viscosity significantly improves the scheme for simulating turbulent flows involving shocks while retaining the shock-capturing capability. The proposed combination of Ducros-type sensor with a negative dilatation sensor removes unnecessary bulk viscosity within expansion and weakly compressible turbulence regions without shocks and allows it to localize near the shocks. It also eliminates the need for a wall-damping function for the bulk viscosity while simulating wall-bounded turbulent flows. For the numerical schemes used, better results are obtained without adding an explicit SGS model than with SGS model at moderate Reynolds number. Inclusion of a SGS model in addition to the low-pass filtering and artificial bulk viscosity results in additional damping of the resolved turbulence. However, investigations at higher Reynolds numbers suggest the need for an explicit SGS model. The flow statistics obtained using the second-order implicit time-integration scheme with three sub-iterations closely agrees with the explicit scheme if the maximum Courant-Friedrichs-Lewy is kept near unity.
KW - Artificial diffusivity method
KW - Compact differences
KW - Compressible turbulence
KW - High-order methods
KW - Large-eddy simulation
KW - Shock capturing
UR - http://www.scopus.com/inward/record.url?scp=72449132905&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=72449132905&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2009.11.005
DO - 10.1016/j.jcp.2009.11.005
M3 - Article
AN - SCOPUS:72449132905
VL - 229
SP - 1739
EP - 1762
JO - Journal of Computational Physics
JF - Journal of Computational Physics
SN - 0021-9991
IS - 5
ER -