Abstract:An improved ghost-cell immersed boundary method coupled with a high-order finite difference scheme is currently proposed to numerically simulate the reflected shock/boundary layer interaction in a wavy-wall tube. The main improvement of current algorithm involves the particular treatment of the solid boundary reconstruction procedure by using both ghost points in solid domain and forcing points in fluid domain together. Owing to this kind of improvement operation, the amount of points used to recover the solid boundary increases greatly and will inevitably promote the accuracy of solid boundary reconstruction, which will play a significant role in the current high-order accurate numerical simulation. A typical validation case of 2-dimensional shock/circular cylinder interaction is conducted and demonstrate that the current improved numerical method can be used to capture the flow details such as shock bifurcation and small vortex very well and improve the numerical precision greatly. The improved method is then applied to numerically investigate the effects of boundary layer on the propagation of the shock in a wavy-wall tube. It is observed that, at low Mach number, the existence of boundary layer has little effect on the shock propagation. As the Mach number increases to 1.9, some certain complicated flow features will occur due to the enhanced boundary layer effects. In the early stage, the flow structures are nearly the same as those without the boundary layers. In the transition stage, a high-speed shear band behind the leading shock is observed and finally evolves into a worm-like structure, whose instability is illustrated by qualitative assessment on its flow features. In the developing stage, a remarkable shock bifurcation phenomenon is present and the worm-like structures will start shedding regularly, which will induce a series of scattering-like waves and asymmetrical flow fields. These characteristics demonstrate that the bifurcated shock waves located in upper and lower walls of the wavy-wall shock tube act as a primary trigger to gradually destabilize the shocked flow.

The Numerical Simulation of Radiative Shocks. I. The Elimination of Numerical Shock Instabilities Using a Local Oscillation Filter

The numerical simulation of two-dimensional fluid flow with strong shocks

Numerical Simulations of Shock/obstacle Interactions Using an Improved Ghost-Cell Immersed Boundary Method

Radiative cooling instability in 1D colliding flows

Numerical Simulation of an Idealised Richtmyer-Meshkov Instability Shock Tube Experiment

A shock-stable numerical scheme accurate for contact discontinuities: Applications to 3D compressible flows

A discussion on numerical shock stability of unstructured finite volume method: Riemann solvers and limiters

Numerical stability analysis of shock-capturing methods for strong shocks II: high-order finite-volume schemes

Continuum Formulation For Non-Equilibrium Shock Structure Calculation

Valid Physical Processes from Numerical Discontinuities in Computational Fluid Dynamics

Temperature-reducing shocks in optically-thin radiative MHD -- analytical and numerical results

Grid-converged Solution and Analysis of the Unsteady Viscous Flow in a Two-dimensional Shock Tube

The Importance of Heat Flux in Quasi-Parallel Collisionless Shocks

Numerical Study of the Slip Line Instabilities in Shock-Wavy Wall Reflection

Assessment of a high-order shock-capturing central-difference scheme for hypersonic turbulent flow simulations

An exact Riemann solver based solution for regular shock refraction

Numerical Simulation of Planar Moving Shock Interacting with a Single Row of Water Columns and Multi-material Interfaces

Shock Wave Formation in Radiative Plasmas

Simulating the long-term evolution of radiative shocks in shock tubes

Radiation and Heat Transport in Divergent Shock-Bubble Interactions

Vorticity Suppression by Particle Lag Effects in Shock-Driven Multiphase Instability