An experiment on deformation of flame under the effect of focusing shock wave reflection is performed with the help of multiple-spark camera to understand the flame instability of the deformation process. Methane and ...An experiment on deformation of flame under the effect of focusing shock wave reflection is performed with the help of multiple-spark camera to understand the flame instability of the deformation process. Methane and oxygen are mixed stoichiometrically to be used in the experiment. Based on Navier-Strokes equations,two-dimensional axisymmetric elementary reactions are numerically simulated. And the simulation results are solved by optical calculation. Shaded pictures by simulation fit well with experimental photos. Focusing reflecton shock waves can affect the flame,which accelerates the deformation of flame and renders violent burning in high-energy flammable gases behind waves. Therefore anticlockwise whirlpool appears. It clusters around the external surface of flame and has a tendency to develop toward the right. Finally,the whirlpool focuses on the right side of the flame,which involves the fresh unfired gases into the whirlpool circle,and consequently the head of mushroom cloud is formed. Meanwhile,when shock wave passes through the flame,the intensity of the shock waves on the axis is strengthened.展开更多
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 an...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.展开更多
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 bounda...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.展开更多
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 obstru...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.展开更多
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 comp...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.展开更多
An experimental study and a numerical simulation were conducted to investigate the mechanical and thermodynamic processes involved in the interaction between shock waves and low density foam. The experiment was done i...An experimental study and a numerical simulation were conducted to investigate the mechanical and thermodynamic processes involved in the interaction between shock waves and low density foam. The experiment was done in a stainless shock tube (80 mm in inner diameter, 10 mm in wall thickness and 5 360 mm in length). The velocities of the incident and reflected compression waves in the foam were measured by using piezo-ceramic pressure sensors. The end-wall peak pressure behind the reflected wave in the foam was measured by using a crystal piezoelectric sensor. It is suggested that the high end-wall pressure may be caused by a rapid contact between the foam and the end-wall surface. Both open-cell and closed-cell foams with different length and density were tested. Through comparing the numerical and experimental end-wall pressure, the permeability coefficients α and β are quantitatively determined.展开更多
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 t...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.展开更多
In this paper, by taking into account the thickness of the incident shock as well as the influence of the boundary layer, we point out that even in a regular reflection there should be present a contact discontinuity....In this paper, by taking into account the thickness of the incident shock as well as the influence of the boundary layer, we point out that even in a regular reflection there should be present a contact discontinuity. By using the smallest energy criterion, the inclined angle of this contact discontinuity can be determined for differen incident angle. Then, with this inclined contact discontinuity, together with the law of conservation of mass, the mechanism for the transition from a regular reflection to a Mach reflection or a von Neumann reflection becomes clear. The important roles played by the leftest point in the reflected shock polar are identified.展开更多
When the pressure ratio increases from the perfectly expanded condition to the third limited condition in which a normal shock is located on the exit plane, shock wave configurations outside the nozzle can be further ...When the pressure ratio increases from the perfectly expanded condition to the third limited condition in which a normal shock is located on the exit plane, shock wave configurations outside the nozzle can be further assorted as no shock wave on the perfectly expanded condition, weak oblique shock reflection in the regular reflection （RR） pressure ratio condition, shock reflection hysteresis in the dual-solution domain of pressure ratio condition, Mach disk configurations in the Mach reflection （MR） pressure ratio condition, the strong oblique shock wave configurations in the corresponding condition, and a normal shock forms on the exit plane in the third limited con- dition. Every critical pressure ratio, especially under regular reflection and Mach reflection pressure ratio conditions, is deduced in the paper according to shock wave reflection theory. A hysteresis phenomenon is also theoretically possible in the dual-solution domain. For a planar Laval nozzle with the cross-section area ratio being 5, different critical pressure ratios are counted in these con- ditions, and numerical simulations are made to demonstrate these various shock wave configurations outside the nozzle. Theoretical analysis and numerical simulations are made to get a more detailed understanding about the shock wave structures outside a Laval nozzle and the RR←→MR transition in the dual-solution domain.展开更多
文摘An experiment on deformation of flame under the effect of focusing shock wave reflection is performed with the help of multiple-spark camera to understand the flame instability of the deformation process. Methane and oxygen are mixed stoichiometrically to be used in the experiment. Based on Navier-Strokes equations,two-dimensional axisymmetric elementary reactions are numerically simulated. And the simulation results are solved by optical calculation. Shaded pictures by simulation fit well with experimental photos. Focusing reflecton shock waves can affect the flame,which accelerates the deformation of flame and renders violent burning in high-energy flammable gases behind waves. Therefore anticlockwise whirlpool appears. It clusters around the external surface of flame and has a tendency to develop toward the right. Finally,the whirlpool focuses on the right side of the flame,which involves the fresh unfired gases into the whirlpool circle,and consequently the head of mushroom cloud is formed. Meanwhile,when shock wave passes through the flame,the intensity of the shock waves on the axis is strengthened.
文摘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.
文摘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.
文摘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.
文摘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.
文摘An experimental study and a numerical simulation were conducted to investigate the mechanical and thermodynamic processes involved in the interaction between shock waves and low density foam. The experiment was done in a stainless shock tube (80 mm in inner diameter, 10 mm in wall thickness and 5 360 mm in length). The velocities of the incident and reflected compression waves in the foam were measured by using piezo-ceramic pressure sensors. The end-wall peak pressure behind the reflected wave in the foam was measured by using a crystal piezoelectric sensor. It is suggested that the high end-wall pressure may be caused by a rapid contact between the foam and the end-wall surface. Both open-cell and closed-cell foams with different length and density were tested. Through comparing the numerical and experimental end-wall pressure, the permeability coefficients α and β are quantitatively determined.
基金supported by the National Research Foundation(NRF)of Korea grant funded by the Korea government(MSIP)(NRF2016R1A2B3016436)
文摘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.
基金supported by the NNSF of China(11271323,91330105)the Zhejiang Provincial Natural Science Foundation of China(LZ13A010002)supported by a GRF grant(City U 11303015)from the Research Grants Council of Hong Kong SAR,China
文摘In this paper, by taking into account the thickness of the incident shock as well as the influence of the boundary layer, we point out that even in a regular reflection there should be present a contact discontinuity. By using the smallest energy criterion, the inclined angle of this contact discontinuity can be determined for differen incident angle. Then, with this inclined contact discontinuity, together with the law of conservation of mass, the mechanism for the transition from a regular reflection to a Mach reflection or a von Neumann reflection becomes clear. The important roles played by the leftest point in the reflected shock polar are identified.
基金supported by the National Natural Science Foundation of China (No. 10702009)
文摘When the pressure ratio increases from the perfectly expanded condition to the third limited condition in which a normal shock is located on the exit plane, shock wave configurations outside the nozzle can be further assorted as no shock wave on the perfectly expanded condition, weak oblique shock reflection in the regular reflection （RR） pressure ratio condition, shock reflection hysteresis in the dual-solution domain of pressure ratio condition, Mach disk configurations in the Mach reflection （MR） pressure ratio condition, the strong oblique shock wave configurations in the corresponding condition, and a normal shock forms on the exit plane in the third limited con- dition. Every critical pressure ratio, especially under regular reflection and Mach reflection pressure ratio conditions, is deduced in the paper according to shock wave reflection theory. A hysteresis phenomenon is also theoretically possible in the dual-solution domain. For a planar Laval nozzle with the cross-section area ratio being 5, different critical pressure ratios are counted in these con- ditions, and numerical simulations are made to demonstrate these various shock wave configurations outside the nozzle. Theoretical analysis and numerical simulations are made to get a more detailed understanding about the shock wave structures outside a Laval nozzle and the RR←→MR transition in the dual-solution domain.
