supersonic turbulent flow and inertial particle is carried out. For the compressible flow, an
eighth-order bandwidth optimization weighted essentially nonoscillatory scheme is used for
shock capturing, and the central finite difference scheme is used for the spatial discretization
of diffusion terms. The three-dimensional ghost zone immersed boundary method is adopted
for solid-fluid interface identification. These numerical schemes are integrated in a direct …