Computational study of shock‐wave interaction with solid obstacles using immersed boundary methods

Abstract

SUMMARYIn this study, an immersed boundary (IB) method based on a direct forcing is coupled with a high-order weighted-essentially non-oscillatory (WENO) scheme to simulate fluid-solid interaction (FSI) problems with complex geometries. The IB is a general simulation method for FSI, whereas the WENO is an efficient scheme for fluid flow simulations and shock waves, and both of them work on regular cartesian grids. The effectiveness and the accuracy of the coupled scheme are first analyzed on well-documented supersonic test problems for a wide range of Mach numbers. The results are in good agreement with both analytical and experimental data. A comprehensive analysis of the interaction of the moving shock through an array of cylinder matrix is then conducted by varying the number of cylinders in the matrix block while keeping the same opening passage. The relaxation length between two adjacent columns of cylinders is kept identical to study uniquely the effect of surface-to-volume ratio of the obstacle matrix. It is shown that the configuration with higher surface-to-volume ratio produces more post-shock flow instabilities downstream of the matrix block. The complex shock/shock and shock/vortex interactions are well resolved by the present computation. It is being observed that after the passage of the shock through the cylinder matrix, eddies of different length scales are generated, but the later stage of shock/vortex and shocklet/vortexlet interactions are different for the two cases. The analysis of the PSD of the total kinetic energy globally conforms to Richardson's inviscid cascade. An intermittent peaked PDF of downstream instantaneous vorticity field is obtained in the limit of Re ! 1. The baroclinic production of vorticity is found to be feeble as previously founded by Sun and Takayama (J. Fluid Mech. 2003; 478:237-256).

Numerical Study of Shock Wave Attenuation in Two-Dimensional Ducts Using Solid Obstacles: How to Utilize Shock Focusing Techniques to Attenuate Shock Waves

“…Chaudhuri et al [15,16] provided a detailed numerical study of shock wave propagation through different arrays of solid obstacles and its degree of attenuation. Obstacles of cylindrical, square and triangular shape were placed inside a shock tube using array-matrix arrangements in both non-staggered and staggered columns.…”

Analysis of shock-wave diffraction over double cylindrical wedges. Part II: Vorticity generation

“…Chaudhuri et al [14] used an immersed boundary (IB) method to study the interaction of the moving shock through an array of cylinder matrix. Their analysis confirmed earlier findings of Sun and Takayama [5], where the baroclinic production of the vorticity is found to be feeble.…”

A new formulation of a spray dispersion model for particle/droplet-laden flows subjected to shock waves

“…In the past, there have been many theoretical and experimental studies investigating the physics of the interaction of droplets or solid particles/obstacles with shock waves (Carrier 1958;Rudinger 1964;Olim et al 1990;Geng et al 1994;Chaudhuri et al 2012Chaudhuri et al , 2013Balakrishnan & Bellan 2017;Mouronval et al 2019;Gai et al 2020). Commonly, the particles are assumed to be at rest before they meet the shock at a given velocity.…”