Massively parallel simulation of shock/boundary-layer interactions

摘要

Massively parallel direct numerical simulation (DNS) of turbulent flows, including shock/boundary-layerinteractions (SBLI), has been carried out in this study. The code uses high order central finite differences for evaluating spatial derivatives in the 3D compressible Navier-Stokes (NS) equations and explicit multi-stage Runge-Kutta algorithm for the time advancement. To detect discontinuities a shock-capturing method of the total variation diminishing (TVD) family, combined with the artificial compression method (ACM), is employed. Some recently developed techniques like entropy splitting of Euler terms and stable boundary treatment are also implemented. The code has been parallelized using the message passing interface (MPI) library and demonstrated its capability to achieve an efficient simulation. Validations have been carried out on a wide range of flow problems, including some well-documented channel and boundary-layer flows. A simulation has been carried out for a viscous flow over a curved bump, a configuration with simple geometry but exhibiting complicated flow physics, in which a shock wave appears under certain flow conditions and the flow separates due to strong adverse pressure gradients and shock/boundary-layer interactions. Simulation results are given and comparisons with theory and existing DNS data are promising.