Analysis of shock wave reflection from fixed and moving boundaries using a stabilized particle method

In the present paper, the efficiency of an enhanced formulation of the stabilized corrective smoothed particle method （CSPM） for simulation of shock wave propagation and reflection from fixed and moving solid boundaries in compressible fluids is investigated. The Lagrangian nature and its accuracy for imposing the boundary conditions are the two main reasons for adoption of CSPM. The governing equations are further modified for imposition of moving solid boundary conditions. In addition to the traditional artificial viscosity, which can remove numerically induced abnormal jumps in the field values, a velocity field smoothing technique is introduced as an efficient method for stabilizing the solution. The method has been implemented for one- and two-dimensional shock wave propagation and reflection from fixed and moving boundaries and the results have been compared with other available solutions. The method has also been adopted for simulation of shock wave propagation and reflection from infinite and finite solid boundaries.

Shock wave mitigation using zig-zag structures and cylindrical obstructions

The present study focuses on the mitigation of shock wave using novel geometric passages in the flow field.The strategy is to produce multiple shock reflections and diffractions in the passage with minimum flow obstruction,which in turn is expected to reduce the shock wave strength at the target location.In the present study the interaction of a plane shock front(generated from a shock tube)with various geometric designs such as,1)zig-zag geometric passage,2)staggered cylindrical obstructions and 3)zigzag passage with cylindrical obstructions have been investigated using computational technique.It is seen from the numerical simulation that,among the various designs,the maximum shock attenuation is produced by the zig-zag passage with cylindrical obstructions which is then followed by zig-zag passage and staggered cylindrical obstructions.A comprehensive investigation on the shock wave reflection and diffraction phenomena happening in the proposed complex passages have also been carried out.In the new zig-zag design,the initial shock wave undergoes shock wave reflection and diffraction process which swaps alternatively as the shock front moves from one turn to the other turn.This cyclic shock reflection and diffraction process helps in diffusing the shock wave energy with practically no obstruction to the flow field.It is found that by combining the shock attenuation ability of zig-zag passage(using shock reflection and diffraction)with the shock attenuation ability of cylindrical blocks(by flow obstruction),a drastic attenuation in shock strength can be achieved with moderate level of flow blocking.

A COMPARISON OF FOUR RECENT NUMERICAL SCHEMES GIViNG HIGH RESOLUTION OF SHOCK WAVE AND CONCENTRATED VORTEX

In recent ten years high resolution difference schenies for the computation of thefull unsteady Eulerian system of equations for invisid compressible gas finds celebratedprogress. This paper tests furtherly, by a complex two-dimensional unsteady problem,four recent schemes. to them attentions are paid. The test problem is the initial stageof a two-dimensional diffraction and reflection of a plane shock wave, impinging on arectangular obstacle. At whose top side there are two sharp corners, near which flow.parameters finds severe variation. There is occurrence of expansion fan with a centerand also concentrated vortices. To simulate them well, the schemes should have goodadaptivity. The special shock Mach number M,=2.068 is so chosen, that at this M,the partical velocity behind impinging shock in fixed coordinate system is just equal tothe speed of sound there, this condition also occurs along a curve in the region ofexpansion fan with a center at the corner. This can clarify the computational featureof different schemes in case,when one of the eigenvalues is just zero. Zero eigenvaluemay spoil some schemes locally. Graphical visualization of the computational resultsmay, show features of the tested schemes about the shock wave resolution, schemeviscosity, expansion wave and the ability. to simulate the process of the generation ofunsteadv concentrated vortex.

THE INTERACTION OF A SHOCK WAVE WITH THE BOUNDARY LAYER IN A REFLECTED SHOCK TUNNEL

The influence of a nontotal reflection on the interaction of a reflected shock wave with the boundary layer in a reflected shock tunnel has been investigated. The calculating method of the velocity, the temperature and the Mach number profiles in the boundary layer in reflected shock fixed coordinates has been obtained. To account for equilibrium real gas effects of nitrogen, the numerical results show that the minimum Mach number in the boundary layer has been moved from the wall into the boundary layer with the increasing of the incident shock Mach number. The minimum Mach number, the shock angle in the bifurcated foot and the jet velocity along the wall to the end plate are reduced owing to the Increasing of the area of nozzle throat. The numerical results are in good agreement with measurements.

A study on characteristics of shock train inside a shock tube

Shock tubes are devices which are used in the investigation of high speed and high temperature flow of compressible gas. lnside a shock tube, the interaction between the reflected shock wave and boundary layer leads to a complex flow phenomenon. Initially a normal shock wave is formed in the shock tube which migrates toward the closed end of the tube and that in turn leads to the reflection of shock. Due to the boundary layer interaction with the reflected shock, the bifurcation of shock wave takes place. The bifurcated shock wave then approaches the contact surface and shock train is generated. Till date only a few studies have been conducted to investigate this shock train phenomenon inside the shock tube. For the present study a computational fluid dynamics （CFD） analysis has been performed on a two dimensional axi-symmetric model of a shock tube using unsteady, compressible Navier-Stokes equations. In order to investigate the detailed characteristics of shock train, parametric studies have been performed by varying different parameters such as the shock tube length, diameter, pressure ratio used inside the shock tube.

Numerical Study of Hydrogen Peroxide Thermal Decomposition in a Shock Tube

Hydrogen peroxide(H_2O_2) has its significance during the combustion of heavy hydrocarbons in the internal combustion(IC) engines. Owing to its importance the measurements of H_2O_2 dissociation rate have been reported mostly using the shock tube apparatus. These types of experimental measurements are although quite reliable but require high cost. On the other hand, numerical simulations provide low cost and reliable solutions especially using computation fluid dynamics(CFD) software. In the current study an experimental shock tube flow is modeled using open access platform OpenFOAM to investigate the thermal decomposition of H_2O_2. Using two different convective schemes, limited Linear and upwind, the propagation of shock wave and resultant dissociation reaction are simulated. The results of the simulations are compared with the experimental data. It is observed that the rate constant measured using the simulation data deviates from the experimental results in the low temperature range and approaches the experimental values as the temperature is raised.