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.

Study of the Shock Wave–Turbulent Boundary Layer Interaction Using a 3D Von Kármán Length Scale

A Turbulence Characteristic Length Scale for Compressible Flows

Research of Compressible Turbulence Model in Shock Wave/Boundary-Layer Interaction Flow at a Compression Corner

One-Dimensional Non-Equilibrium Shock Wave Structure Analysis by Simplified Conventional Burnett Equations

A Scale-Adaptive Turbulence Model Based on the k-Equation and Recalibrated Reynolds Stress Constitutive Relation

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

Direct numerical simulation of a supersonic turbulent boundary layer over a compression–decompression corner

A New Compressibility Modification K -Ε Turbulence Model with Shock Unsteadiness Effect

Numerical investigations on interactions between 2D/3D conical shock wave and axisymmetric boundary layer at Ma=2.2

Amplification and structure of streamwise-velocity fluctuations in compression-corner shock-wave/turbulent boundary-layer interactions

A modified k-ε turbulence model for compressible flow application

Three-Dimensional Lattice Boltzmann Model for High-Speed Compressible Flows

Assessment of a Novel K–ω Turbulence Model for Transonic Centrifugal Impeller Simulations

Reynolds-Averaged Numerical Simulations of Conical Shock-Wave/Boundary-Layer Interactions

Numerical Studies of Scaling Exponents in Compressible Wall-Bounded Turbulence

A K- Ɛ Turbulence Model Considering Compressibility in Three-Dimensional Transonic Turbulent Flow Calculation

Applicability Investigation about Several Common Schemes and Turbulence Models in High Speed Flow

Modeling of Turbulent Kinetic Equation in Converging-Reflecting Spherical Shock-Turbulence Interaction

Hybrid LES/RANS Simulation of Shock/Turbulent Boundary-Layer Interactions

Development of K-R Turbulence Model for Wall -Bounded Flows

Simulation of three-dimensional compressible decaying isotropic turbulence using a redesigned discrete unified gas kinetic scheme