Recently micro shock tubes have been widely used in many engineering and industrial fields, but the characteristics of unsteady flow are not well known to date in micro shock tubes. Compared to conventional shock tube...Recently micro shock tubes have been widely used in many engineering and industrial fields, but the characteristics of unsteady flow are not well known to date in micro shock tubes. Compared to conventional shock tubes with macro scales, flows related to shock waves in micro shock tubes are highly complicated. Stronger viscous and dissipative interactions make shock wave dynamic behaviors significantly different from theoretical predictions. In the present study, a CFD work was applied to the unsteady compressible Navier-Stokes equations which were solved using a fully implicit finite volume scheme. The diaphragm pressure ratio and shock tube diameter were varied to investigate their effects on micro shock tube flows. Different wall boundary conditions were also performed to observe shock wave and contact surface propagation with no slip and slip walls. Detailed flow characteristics at the foot of shock wave and contact surface propagation were known from the present numerical simulations.展开更多
Explosive blast injury has become the most prevalent injury in recent military conflicts and terrorist attacks. The magnitude of this kind of polytrauma is complex due to the basic physics of blast and the surrounding...Explosive blast injury has become the most prevalent injury in recent military conflicts and terrorist attacks. The magnitude of this kind of polytrauma is complex due to the basic physics of blast and the surrounding environments. Therefore, development of stable, reproducible and controllable animal model using an ideal blast simulation device is the key of blast injury research. The present review addresses the modeling of blast injury and applications of shock tubes.展开更多
The thermochemical non-equilibrium phenomena encountered by hypersonic vehicles present significant challenges in their design.To investigate the thermochemical reaction flow behind shock waves,the non-equilibrium rad...The thermochemical non-equilibrium phenomena encountered by hypersonic vehicles present significant challenges in their design.To investigate the thermochemical reaction flow behind shock waves,the non-equilibrium radiation in the visible range using a shock tube was studied.Experiments were conducted with a shock velocity of 4.7 km/s,using nitrogen at a pressure of 20 Pa.To address measurement difficulties associated with weak radiation,a special square section shock tube with a side length of 380 mm was utilized.A high-speed camera characterized the shock wave’s morphology,and a spectrograph and a monochromator captured the radiation.The spectra were analyzed,and the numerical spectra were compared with experimental results,showing a close match.Temperature changes behind the shock wave were obtained and compared with numerical predictions.The findings indicate that the vibrational temperatures are overestimated,while the vibrational relaxation time is likely underestimated,due to the oversimplified portrayals of the non-equilibrium relaxation process in the models.Additionally,both experimental and simulated time-resolved profiles of radiation intensity at specific wavelengths were analyzed.The gathered data aims to enhance computational fluid dynamics codes and radiation models,improving their predictive accuracy.展开更多
This study presents a fragility curve to assess explosively induced damage to military vehicle tires based on shock tube experiments.To replicate lateral damage scenarios that may occur in real battlefield environment...This study presents a fragility curve to assess explosively induced damage to military vehicle tires based on shock tube experiments.To replicate lateral damage scenarios that may occur in real battlefield environments involving missile or bomb detonations,extreme overpressure conditions were generated using a shock tube.The influence of explosive charge mass on tire damage was quantitatively evaluated.Experimental results identified two critical failure thresholds:for loss of pressure,the threshold was 354 kPa peak overpressure and 3052 kPa·ms impulse;for rupture,the values were 485 kPa and 4237 kPa-ms,respectively.The same damage profile was reproduced through finite element analysis(FEA),verifying the reliability of the simulation.A Single Degree of Freedom(SDOF)model and Kingery-Bulmash(K-B)chart were employed to generate pressure-impulse data as a function of standoff distance.These data were applied to a finite element tire model using the BLAST ENHANCED keyword in LS-DYNA.The applied peak overpressures were identical to the experimental values with a 24%-27%difference in impulse.The simulation also captured recurring bead rim separation phenomenon,leading to internal pressure loss consistent with high-speed camera observations from the experiments.The resulting fragility curve clearly defines the threshold conditions for tire damage and provides a standardized damage assessment model applicable to various explosive charge masses and stand-off distances.The proposed model offers a quantitative basis for evaluating tire vulnerability,providing foundational reference data for defense applications.Specifically,the findings are expected to serve as a reliable source for weapon effects analysis and target vulnerability assessments involving wheeled military vehicles.展开更多
In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dim...In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dimensional unsteady two-phase flows. The two-phase shock wave and the flow field behind it in a dusty gas shock tube are calculated and the time-dependent change of the flow parameters for the gas and particle phase are obtained. The numerical results indicate that both the two methods can give the relaxation structure of the two-phase shocks with a sharp discontinuous front and that the GRP method has the advantages of less time-consuming and higher accuracy over the RCM method.展开更多
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.展开更多
The performance of combustion driver ignited by multi-spark plugs distributed along axial direction has been analysed and tested. An improved ignition method with three circumferential equidistributed ignitors at main...The performance of combustion driver ignited by multi-spark plugs distributed along axial direction has been analysed and tested. An improved ignition method with three circumferential equidistributed ignitors at main diaphragm has been presented, by which the produced incident shock waves have higher repeatability, and better steadiness in the pressure, temperature and velocity fields of flow behind the incidence shock, and thus meets the requirements of aerodynamic experiment. The attachment of a damping section at the end of the driver can eliminate the high reflection pressure produced by detonation wave, and the backward detonation driver can be employed to generate high enthalpy and high density test flow. The incident shock wave produced by this method is well repeated and with weak attenuation. The reflection wave caused by the contracted section at the main diaphragm will weaken the unfavorable effect of rarefaction wave behind the detonation wave, which indicates that the forward detonation driver can be applied in the practice. For incident shock wave of identical strength, the initial pressure of the forward detonation driver is about 1 order of magnitude lower than that of backward detonation.展开更多
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.展开更多
This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters n_e,v,ω,L,ω_b.The φ800 mm high temperature shock tube ...This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters n_e,v,ω,L,ω_b.The φ800 mm high temperature shock tube has been used to produce a uniform plasma.In order to get the attenuation of the electromagnetic wave through the plasma behind a shock wave,the microwave transmission has been used to measure the relative change of the wave power.The working frequency is f=(2～35)GHz(ω=2πf,wave length λ=15 cm～8 mm).The electron density in the plasma is n_e=(3×10^(10)～1×10^(14))cm^(-3).The collision frequency v=(1×10~8～6×10^(10))Hz.The thickness of the plasma layer L=(2～80)cm.The electron circular frequency ω_b=eBo/m_e,magnetic flux density B_0=(0～0.84)T.The experimental results show that when the plasma layer is thick (such as L/λ(?)10),the correlation between the attenuation coefficients of the electromagnetic waves and the parameters n_e,v,ω,L determined from the measurements are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma.When the plasma layer is thin(such as when both L and λ are of the same order),the theoretical results are only in a qualitative agreement with the experimental observations in the present parameter range, but the formula of the electromagnetic wave propagation theory in an uniform infinite plasma can not be used for quantitative computations of the correlation between the attenuation coefficients and the parameters n_e,v,ω,L.In fact,if ω<ω_p,v^2(?)ω~2,the power attenuations K of the electromagnetic waves obtained from the measurements in the thin-layer plasma are much smaller than those of the theoretical predictions.On the other hand,if ω>ω_p,v^2(?)ω~2(just v≈f),the measurements are much larger than the theoretical results.Also,we have measured the electromagnetic wave power attenuation value under the magnetic field and without a magnetic field.The result indicates that the value measured under the magnetic field shows a distinct improvement.展开更多
We propose a method of applying a static magnetic field to reduce the attenuation of the magnetic field component(SH) of low-frequency electromagnetic(LF EM) waves in dense plasma. The principle of this method is to a...We propose a method of applying a static magnetic field to reduce the attenuation of the magnetic field component(SH) of low-frequency electromagnetic(LF EM) waves in dense plasma. The principle of this method is to apply a static magnetic field to limit electron movement, thereby increasing the equivalent resistance and thus reducing the induced current and SH. We consider the static magnetic field acting on the plasma of the entire induced current loop rather than on the local plasma, where the induced current is excited by the magnetic field component of LF EM waves. Analytical expressions of SH suitable for magnetized cylindrical enveloping plasma are derived by adopting an equivalent circuit approach, by which SHis calculated with respect to various plasma parameter settings. The results show that SH can be reduced under a static magnetic field and the maximum magnetic field strength that mitigates blackout is less than 0.1 T. Experiments in which LF EM waves propagate in a shock-tubegenerated magnetized cylindrical enveloping plasma are also conducted. SH measured under the magnetic field(the magnetic field strength B0 acting on the magnetic field probe was about0.06 T) reduces at f=10 MHz and f=30 MHz when ne≈1.9×1013 cm-3, which is consistent with theoretical results. The verification of the theory thus suggests that applying a static magnetic field with a weak magnetic field has the potential to improve the transmission capacity of LF EM waves in dense plasma.展开更多
In shock tube experiments,the interaction between the reflected shock and boundary layer can induce shock bifurcation and weak ignition.The weak ignition can greatly affect the ignition delay time measurement in a sho...In shock tube experiments,the interaction between the reflected shock and boundary layer can induce shock bifurcation and weak ignition.The weak ignition can greatly affect the ignition delay time measurement in a shock tube experiment.In this work,two-dimensional simulations considering detailed chemistry and transport are conducted to investigate the shock bifurcation and non-uniform ignition behind a retlected shock.The objectives are to interpret the formation of shock bifurcation induced by the reflected shock and boundary layer interaction and to investigate the weak ignition and its transition to strong ignition for both hydrogen and dimethyl ether.It is found that the non-uniform reflection of the incident shock at the end wall produces a wedge-shaped oblique shock foot at the wall.The wedge-shaped structure results in strong interactions between reflected shock and boundary layer,which induces the shock bifurcation.It is demonstrated that the local high-temperature spots at the foot of the bifurcated shock is caused by viscous dissipation and pressure work.As the post-reflected shock temperature increases,the transition from weak ignition to strong ignition in a stoichiometric hydrogen/oxygen mixture is observed.The relative sensitivity of ignition delay time to the post-rellected shock temperature is introduced to characterize the appearance of weak ignition behind the reflected shock.Unlike in the hydrogen/oxygen mixture,weak ignition is not observed in the stoichiometric dimethyl-ether/oxygen mixture since it has a relatively longer ignition delay time and smaller relative sensitivity.展开更多
The shock tubes with area change are used in the free piston shock tunnels,owing to its higher driver effect.For optimized operation of this kind of shock tube,a computer program for fast simulation of transient hyper...The shock tubes with area change are used in the free piston shock tunnels,owing to its higher driver effect.For optimized operation of this kind of shock tube,a computer program for fast simulation of transient hypersonic flow is presented.The numerical modeling embodied within this code is based on a quasi-one-dimensional Lagrangian description of the gas dynamics.In this code,a mass-loss model is also applied by using Mirels′theory of shock attenuation.The simulation of particular condition for T4 free piston shock tunnel is conducted and compared with experimental measurements and numerical simulation.The results provide good estimate for shock speed and pressure obtained after shock reflection.展开更多
Ignition delay times for n-decane/O2/Ar mixtures were measured behind reflected shock waves using endwall pressure and CH＊ emission measurements in a heated shock tube. The initial postshock conditions cover pressure...Ignition delay times for n-decane/O2/Ar mixtures were measured behind reflected shock waves using endwall pressure and CH＊ emission measurements in a heated shock tube. The initial postshock conditions cover pressures of 0.09-0.26 MPa, temperatures of 1 227-1 536 K, and oxygen mole fractions of 3.9%-20.7% with an equivalence ratio of 1.0. The correlation formula of ignition delay dependence on pressure, temperature, and oxygen mole fraction was obtained. The current data are in good agreement with available low-pressure experimental data, and they are then compared with the prediction of a kinetic mechanism. The current measurements extend the kinetic modeling targets for the n-decane combustion at low pressures.展开更多
This report describes a new method for measuring the temperature of the gas behind the reflected shock wave in shock tube, corresponding to the reservoir temperature of a shock tunnel, based on the chemical reaction o...This report describes a new method for measuring the temperature of the gas behind the reflected shock wave in shock tube, corresponding to the reservoir temperature of a shock tunnel, based on the chemical reaction of small amount of CF4 premixed in the test gas. The final product C2F4 is used as the temperature indicator, which is sampled and detected by a gas chromatography in the experiment. The detected concentration of C2F4 is correlated to the temperature of the reflected shock wave with the initial pressure P-1 and test time tau as parameters in the temperature range 3 300 K < T < 5 600 K, pressure range 5 kPa < P1 <12 kPa and tau similar or equal to 0.4 ms.展开更多
The pyrolysis of cyclohexane,methylcyclohexane,and ethylcyclohexane have been studied behind reflected shock waves at pressures of 5 and10 bar and at temperatures of 930-1550 K for 0.05%fuel diluted by Argon.A single-...The pyrolysis of cyclohexane,methylcyclohexane,and ethylcyclohexane have been studied behind reflected shock waves at pressures of 5 and10 bar and at temperatures of 930-1550 K for 0.05%fuel diluted by Argon.A single-pulse shock tube(SPST)is used to perform the pyrolysis experiments at reaction times varying from 1.65 to 1.74 ms.Major products are obtained and quantified using gas chromatography analysis.A flame ionization detector and a thermal conductivity detector are used for species identification and quantification.Kinetic modeling has been performed using several detailed and lumped chemical kinetic mechanisms.Differences in modeling results among the kinetic models are described.Reaction path analysis and sensitivity analysis are performed to determine the important reactions controlling fuel pyrolysis and their influence on the predicted concentrations of reactant and product species profiles.The present work provides new fundamental knowledge in understating pyrolysis characteristics of cyclohexane compounds and additional data set for detailed kinetic mechanism development.展开更多
Shock tube experiments were carried out to investigate dynamic behavior of Ultra-high hardness(UHH)steel and Aluminium(Al) sheets of 0.8 mm thickness at 0.55, 0.9 and 1.18 MPa peak-over pressure.Experimental results s...Shock tube experiments were carried out to investigate dynamic behavior of Ultra-high hardness(UHH)steel and Aluminium(Al) sheets of 0.8 mm thickness at 0.55, 0.9 and 1.18 MPa peak-over pressure.Experimental results showed that center point deflection increases with an increase in peak-over pressure for Al sheets. However, UHH steel sheets showed negligible deformation when loaded at low peak-over pressures and showed sudden brittle failure at high peak-over pressures. Similar results were obtained by quasi-static testing, UHH steel failed abruptly while Al showed ductile behavior. Results from literature indicate that to protect structures against shock loading it is necessary that they dissipate energy via plastic deformation. The Al sheets were shown to deform plastically both in quasi-static and shock loading. Thus, hardness along with ductility is required to dissipate supersonic shock waves.展开更多
In this study, 75% and 96% argon diluent conditions were selected to determine the ig- nition delay time of stoichiometric mixture of C2Ha/O2/Ar within a range of pressures （1.3-：3.0 arm） and temperatures （1092-17...In this study, 75% and 96% argon diluent conditions were selected to determine the ig- nition delay time of stoichiometric mixture of C2Ha/O2/Ar within a range of pressures （1.3-：3.0 arm） and temperatures （1092-1743 K）. Results showed a logarithmic linear rela- tionship of the ignition delay time with the reciprocal of temperatures. Under both two diluent conditions, ignition delay time decreased with increased temperature. By multiple linear regression analysis, the ignition delay correlation was deduced. According to this correlation, the calculated ignition delay time in 96% diluent was found to be nearly five times that in 75% diluent. To explain this discrepancy, the hard-sphere collision theory was adopted, and the collision numbers of ethylene to oxygen were calculated. The total collision numbers of ethylene to oxygen were 5.99×10^30 s^-1cm^-3 in 75% diluent and 1.53×10^29 s^-1cm^-3 in 96% diluent （about 40 times that in 75% diluent）. According to the discrepancy between ignition delay time and collision numbers, viz. 5 times corresponds to 40 times, the steric factor can展开更多
The localized differential quadrature (LDQ) method is a numerical technique with high accuracy for solving most kinds of nonlinear problems in engineering and can overcome the difficulties of other methods (such as di...The localized differential quadrature (LDQ) method is a numerical technique with high accuracy for solving most kinds of nonlinear problems in engineering and can overcome the difficulties of other methods (such as difference method) to numerically evaluate the derivatives of the functions.Its high efficiency and accuracy attract many engineers to apply the method to solve most of the numerical problems in engineering.However,difficulties can still be found in some particular problems.In the following study,the LDQ was applied to solve the Sod shock tube problem.This problem is a very particular kind of problem,which challenges many common numerical methods.Three different examples were given for testing the robustness and accuracy of the LDQ.In the first example,in which common initial conditions and solving methods were given,the numerical oscillations could be found dramatically;in the second example,the initial conditions were adjusted appropriately and the numerical oscillations were less dramatic than that in the first example;in the third example,the momentum equation of the Sod shock tube problem was corrected by adding artificial viscosity,causing the numerical oscillations to nearly disappear in the process of calculation.The numerical results presented demonstrate the detailed difficulties encountered in the calculations,which need to be improved in future work.However,in summary,the localized differential quadrature is shown to be a trustworthy method for solving most of the nonlinear problems in engineering.展开更多
In this work,the theoretical analysis and experiment results investigating the influence of plasma-induced reflected wave variations on microwave transmission characterization are presented.Firstly,an analytical trans...In this work,the theoretical analysis and experiment results investigating the influence of plasma-induced reflected wave variations on microwave transmission characterization are presented.Firstly,an analytical transmission line model for transmission characterization of plasma in shock tube is derived and validated against full-wave simulation.Then,the theoretical analysis of transmission characterization based on a time-dependent reconstruction algorithm that takes into account the variations of reflected wave is presented and the influence of reflection variations under various states of plasma is also investigated.The unusual increase in the amplitude of transmitted wave is theoretically predicted and experimentally demonstrated as well.Finally,the experiment results are also presented to illustrate the effects of reflected wave variations in practical microwave transmission characterization of supersonic plasma excited in shock tube.展开更多
The purpose of this paper is to consider 1D Riemann shock tube to investigate the formation and propagation of compression waves leading to formation, propagation and reflection of 1D normal shocks using simplified ma...The purpose of this paper is to consider 1D Riemann shock tube to investigate the formation and propagation of compression waves leading to formation, propagation and reflection of 1D normal shocks using simplified mathematical models commonly used in the published work as well as using complete mathematical models based on Conservation and Balance Laws (CBL) of classical continuum mechanics and constitutive theories for compressible viscous medium derived using entropy inequality and representation theorem. This work is aimed at resolving compression waves, the shock structure, shock formation, propagation and reflection of fully formed shocks. Evolutions obtained from the mathematical models always satisfy differentiability requirements in space and time dictated by the highest order of the derivatives of the dependent variables in the mathematical models investigated. All solutions reported in this paper including boundary conditions and initial conditions are always analytic. Solutions of the mathematical models are obtained using the space-time finite element method in which the space-time integral forms are space-time variationally consistent ensuring unconditionally stable computations during the entire evolution. Solution for a space-time strip or slab is calculated and is time marched upon convergence to obtain complete evolution for the desired space-time domain, thus ensuring time accurate evolutions. The space-time local approximation over a space-time element of a space-time strip or slab is p-version hierarchical with higher-order global differentiability in space and time, i.e., we consider scalar product approximation spaces in which k = (k<sub>1</sub>, k<sub>2</sub>) are the order of the space in space and time and p = (p<sub>1</sub>, p<sub>2</sub>) are p-levels of the approximations in space and time. Model problem studies are presented for different mathematical models and are compared with solutions obtained from the complete mathematical model based on CBL and constitutive theories for viscous compressible medium to illustrate the deficiencies and shortcomings of the simplified and approximate models in simulating correct physics of normal shocks.展开更多
基金supported by the National Research Foundation of Korea(NRF)grant funded by the Korea government(MEST)(2011-0017506)
文摘Recently micro shock tubes have been widely used in many engineering and industrial fields, but the characteristics of unsteady flow are not well known to date in micro shock tubes. Compared to conventional shock tubes with macro scales, flows related to shock waves in micro shock tubes are highly complicated. Stronger viscous and dissipative interactions make shock wave dynamic behaviors significantly different from theoretical predictions. In the present study, a CFD work was applied to the unsteady compressible Navier-Stokes equations which were solved using a fully implicit finite volume scheme. The diaphragm pressure ratio and shock tube diameter were varied to investigate their effects on micro shock tube flows. Different wall boundary conditions were also performed to observe shock wave and contact surface propagation with no slip and slip walls. Detailed flow characteristics at the foot of shock wave and contact surface propagation were known from the present numerical simulations.
基金We gratefully appreciate Academician Wang ZG for critically reading and editing the manuscript. This work was supported by the grants from National Natural Science Foundation of China （No. 81201461 ）, Key Project of Medicine and Health of PLA （No. 08G098） and the Natural Science Foundation of Chongqing, China （No. CSTC2012jjA10107）.
文摘Explosive blast injury has become the most prevalent injury in recent military conflicts and terrorist attacks. The magnitude of this kind of polytrauma is complex due to the basic physics of blast and the surrounding environments. Therefore, development of stable, reproducible and controllable animal model using an ideal blast simulation device is the key of blast injury research. The present review addresses the modeling of blast injury and applications of shock tubes.
基金supported by the Key-Area Research and Development Program of Guangdong Province(Grant No.2021B0909060004)the National Natural Science Foundation of China(Grant Nos.12072355 and 92271117)the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDB0620202).
文摘The thermochemical non-equilibrium phenomena encountered by hypersonic vehicles present significant challenges in their design.To investigate the thermochemical reaction flow behind shock waves,the non-equilibrium radiation in the visible range using a shock tube was studied.Experiments were conducted with a shock velocity of 4.7 km/s,using nitrogen at a pressure of 20 Pa.To address measurement difficulties associated with weak radiation,a special square section shock tube with a side length of 380 mm was utilized.A high-speed camera characterized the shock wave’s morphology,and a spectrograph and a monochromator captured the radiation.The spectra were analyzed,and the numerical spectra were compared with experimental results,showing a close match.Temperature changes behind the shock wave were obtained and compared with numerical predictions.The findings indicate that the vibrational temperatures are overestimated,while the vibrational relaxation time is likely underestimated,due to the oversimplified portrayals of the non-equilibrium relaxation process in the models.Additionally,both experimental and simulated time-resolved profiles of radiation intensity at specific wavelengths were analyzed.The gathered data aims to enhance computational fluid dynamics codes and radiation models,improving their predictive accuracy.
基金part of the Agency for Defense Development(ADD)research project on Weapon lethality/effectiveness analysis technology for material targets and grant funded by the korean goverment(511225-912A03301)。
文摘This study presents a fragility curve to assess explosively induced damage to military vehicle tires based on shock tube experiments.To replicate lateral damage scenarios that may occur in real battlefield environments involving missile or bomb detonations,extreme overpressure conditions were generated using a shock tube.The influence of explosive charge mass on tire damage was quantitatively evaluated.Experimental results identified two critical failure thresholds:for loss of pressure,the threshold was 354 kPa peak overpressure and 3052 kPa·ms impulse;for rupture,the values were 485 kPa and 4237 kPa-ms,respectively.The same damage profile was reproduced through finite element analysis(FEA),verifying the reliability of the simulation.A Single Degree of Freedom(SDOF)model and Kingery-Bulmash(K-B)chart were employed to generate pressure-impulse data as a function of standoff distance.These data were applied to a finite element tire model using the BLAST ENHANCED keyword in LS-DYNA.The applied peak overpressures were identical to the experimental values with a 24%-27%difference in impulse.The simulation also captured recurring bead rim separation phenomenon,leading to internal pressure loss consistent with high-speed camera observations from the experiments.The resulting fragility curve clearly defines the threshold conditions for tire damage and provides a standardized damage assessment model applicable to various explosive charge masses and stand-off distances.The proposed model offers a quantitative basis for evaluating tire vulnerability,providing foundational reference data for defense applications.Specifically,the findings are expected to serve as a reliable source for weapon effects analysis and target vulnerability assessments involving wheeled military vehicles.
文摘In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dimensional unsteady two-phase flows. The two-phase shock wave and the flow field behind it in a dusty gas shock tube are calculated and the time-dependent change of the flow parameters for the gas and particle phase are obtained. The numerical results indicate that both the two methods can give the relaxation structure of the two-phase shocks with a sharp discontinuous front and that the GRP method has the advantages of less time-consuming and higher accuracy over the RCM method.
基金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.
基金State Science and Technology CommitteeNational Natural Foundation of Science of China (19082012)+1 种基金 Chinese Academy of SciencesProject of National High Technology of China.
文摘The performance of combustion driver ignited by multi-spark plugs distributed along axial direction has been analysed and tested. An improved ignition method with three circumferential equidistributed ignitors at main diaphragm has been presented, by which the produced incident shock waves have higher repeatability, and better steadiness in the pressure, temperature and velocity fields of flow behind the incidence shock, and thus meets the requirements of aerodynamic experiment. The attachment of a damping section at the end of the driver can eliminate the high reflection pressure produced by detonation wave, and the backward detonation driver can be employed to generate high enthalpy and high density test flow. The incident shock wave produced by this method is well repeated and with weak attenuation. The reflection wave caused by the contracted section at the main diaphragm will weaken the unfavorable effect of rarefaction wave behind the detonation wave, which indicates that the forward detonation driver can be applied in the practice. For incident shock wave of identical strength, the initial pressure of the forward detonation driver is about 1 order of magnitude lower than that of backward detonation.
文摘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.
基金The project supported by the National Natural Science Foundation of China (19982005 and 10032050) and the National Defense Science Foundation of China
文摘This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters n_e,v,ω,L,ω_b.The φ800 mm high temperature shock tube has been used to produce a uniform plasma.In order to get the attenuation of the electromagnetic wave through the plasma behind a shock wave,the microwave transmission has been used to measure the relative change of the wave power.The working frequency is f=(2～35)GHz(ω=2πf,wave length λ=15 cm～8 mm).The electron density in the plasma is n_e=(3×10^(10)～1×10^(14))cm^(-3).The collision frequency v=(1×10~8～6×10^(10))Hz.The thickness of the plasma layer L=(2～80)cm.The electron circular frequency ω_b=eBo/m_e,magnetic flux density B_0=(0～0.84)T.The experimental results show that when the plasma layer is thick (such as L/λ(?)10),the correlation between the attenuation coefficients of the electromagnetic waves and the parameters n_e,v,ω,L determined from the measurements are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma.When the plasma layer is thin(such as when both L and λ are of the same order),the theoretical results are only in a qualitative agreement with the experimental observations in the present parameter range, but the formula of the electromagnetic wave propagation theory in an uniform infinite plasma can not be used for quantitative computations of the correlation between the attenuation coefficients and the parameters n_e,v,ω,L.In fact,if ω<ω_p,v^2(?)ω~2,the power attenuations K of the electromagnetic waves obtained from the measurements in the thin-layer plasma are much smaller than those of the theoretical predictions.On the other hand,if ω>ω_p,v^2(?)ω~2(just v≈f),the measurements are much larger than the theoretical results.Also,we have measured the electromagnetic wave power attenuation value under the magnetic field and without a magnetic field.The result indicates that the value measured under the magnetic field shows a distinct improvement.
基金supported by National Natural Science Foundation of China (Nos. 61771370, and 11704296)
文摘We propose a method of applying a static magnetic field to reduce the attenuation of the magnetic field component(SH) of low-frequency electromagnetic(LF EM) waves in dense plasma. The principle of this method is to apply a static magnetic field to limit electron movement, thereby increasing the equivalent resistance and thus reducing the induced current and SH. We consider the static magnetic field acting on the plasma of the entire induced current loop rather than on the local plasma, where the induced current is excited by the magnetic field component of LF EM waves. Analytical expressions of SH suitable for magnetized cylindrical enveloping plasma are derived by adopting an equivalent circuit approach, by which SHis calculated with respect to various plasma parameter settings. The results show that SH can be reduced under a static magnetic field and the maximum magnetic field strength that mitigates blackout is less than 0.1 T. Experiments in which LF EM waves propagate in a shock-tubegenerated magnetized cylindrical enveloping plasma are also conducted. SH measured under the magnetic field(the magnetic field strength B0 acting on the magnetic field probe was about0.06 T) reduces at f=10 MHz and f=30 MHz when ne≈1.9×1013 cm-3, which is consistent with theoretical results. The verification of the theory thus suggests that applying a static magnetic field with a weak magnetic field has the potential to improve the transmission capacity of LF EM waves in dense plasma.
基金supported by the National Natural Science Foundation of China(Grant Nos.52006001,and 52176096).
文摘In shock tube experiments,the interaction between the reflected shock and boundary layer can induce shock bifurcation and weak ignition.The weak ignition can greatly affect the ignition delay time measurement in a shock tube experiment.In this work,two-dimensional simulations considering detailed chemistry and transport are conducted to investigate the shock bifurcation and non-uniform ignition behind a retlected shock.The objectives are to interpret the formation of shock bifurcation induced by the reflected shock and boundary layer interaction and to investigate the weak ignition and its transition to strong ignition for both hydrogen and dimethyl ether.It is found that the non-uniform reflection of the incident shock at the end wall produces a wedge-shaped oblique shock foot at the wall.The wedge-shaped structure results in strong interactions between reflected shock and boundary layer,which induces the shock bifurcation.It is demonstrated that the local high-temperature spots at the foot of the bifurcated shock is caused by viscous dissipation and pressure work.As the post-reflected shock temperature increases,the transition from weak ignition to strong ignition in a stoichiometric hydrogen/oxygen mixture is observed.The relative sensitivity of ignition delay time to the post-rellected shock temperature is introduced to characterize the appearance of weak ignition behind the reflected shock.Unlike in the hydrogen/oxygen mixture,weak ignition is not observed in the stoichiometric dimethyl-ether/oxygen mixture since it has a relatively longer ignition delay time and smaller relative sensitivity.
文摘The shock tubes with area change are used in the free piston shock tunnels,owing to its higher driver effect.For optimized operation of this kind of shock tube,a computer program for fast simulation of transient hypersonic flow is presented.The numerical modeling embodied within this code is based on a quasi-one-dimensional Lagrangian description of the gas dynamics.In this code,a mass-loss model is also applied by using Mirels′theory of shock attenuation.The simulation of particular condition for T4 free piston shock tunnel is conducted and compared with experimental measurements and numerical simulation.The results provide good estimate for shock speed and pressure obtained after shock reflection.
基金supported by the National Natural Science Foundation of China (91016002)
文摘Ignition delay times for n-decane/O2/Ar mixtures were measured behind reflected shock waves using endwall pressure and CH＊ emission measurements in a heated shock tube. The initial postshock conditions cover pressures of 0.09-0.26 MPa, temperatures of 1 227-1 536 K, and oxygen mole fractions of 3.9%-20.7% with an equivalence ratio of 1.0. The correlation formula of ignition delay dependence on pressure, temperature, and oxygen mole fraction was obtained. The current data are in good agreement with available low-pressure experimental data, and they are then compared with the prediction of a kinetic mechanism. The current measurements extend the kinetic modeling targets for the n-decane combustion at low pressures.
基金The project supported by the China Aerodynamics Project for Basic Researches(J13.5.2 ZK04)
文摘This report describes a new method for measuring the temperature of the gas behind the reflected shock wave in shock tube, corresponding to the reservoir temperature of a shock tunnel, based on the chemical reaction of small amount of CF4 premixed in the test gas. The final product C2F4 is used as the temperature indicator, which is sampled and detected by a gas chromatography in the experiment. The detected concentration of C2F4 is correlated to the temperature of the reflected shock wave with the initial pressure P-1 and test time tau as parameters in the temperature range 3 300 K < T < 5 600 K, pressure range 5 kPa < P1 <12 kPa and tau similar or equal to 0.4 ms.
文摘The pyrolysis of cyclohexane,methylcyclohexane,and ethylcyclohexane have been studied behind reflected shock waves at pressures of 5 and10 bar and at temperatures of 930-1550 K for 0.05%fuel diluted by Argon.A single-pulse shock tube(SPST)is used to perform the pyrolysis experiments at reaction times varying from 1.65 to 1.74 ms.Major products are obtained and quantified using gas chromatography analysis.A flame ionization detector and a thermal conductivity detector are used for species identification and quantification.Kinetic modeling has been performed using several detailed and lumped chemical kinetic mechanisms.Differences in modeling results among the kinetic models are described.Reaction path analysis and sensitivity analysis are performed to determine the important reactions controlling fuel pyrolysis and their influence on the predicted concentrations of reactant and product species profiles.The present work provides new fundamental knowledge in understating pyrolysis characteristics of cyclohexane compounds and additional data set for detailed kinetic mechanism development.
文摘Shock tube experiments were carried out to investigate dynamic behavior of Ultra-high hardness(UHH)steel and Aluminium(Al) sheets of 0.8 mm thickness at 0.55, 0.9 and 1.18 MPa peak-over pressure.Experimental results showed that center point deflection increases with an increase in peak-over pressure for Al sheets. However, UHH steel sheets showed negligible deformation when loaded at low peak-over pressures and showed sudden brittle failure at high peak-over pressures. Similar results were obtained by quasi-static testing, UHH steel failed abruptly while Al showed ductile behavior. Results from literature indicate that to protect structures against shock loading it is necessary that they dissipate energy via plastic deformation. The Al sheets were shown to deform plastically both in quasi-static and shock loading. Thus, hardness along with ductility is required to dissipate supersonic shock waves.
文摘In this study, 75% and 96% argon diluent conditions were selected to determine the ig- nition delay time of stoichiometric mixture of C2Ha/O2/Ar within a range of pressures （1.3-：3.0 arm） and temperatures （1092-1743 K）. Results showed a logarithmic linear rela- tionship of the ignition delay time with the reciprocal of temperatures. Under both two diluent conditions, ignition delay time decreased with increased temperature. By multiple linear regression analysis, the ignition delay correlation was deduced. According to this correlation, the calculated ignition delay time in 96% diluent was found to be nearly five times that in 75% diluent. To explain this discrepancy, the hard-sphere collision theory was adopted, and the collision numbers of ethylene to oxygen were calculated. The total collision numbers of ethylene to oxygen were 5.99×10^30 s^-1cm^-3 in 75% diluent and 1.53×10^29 s^-1cm^-3 in 96% diluent （about 40 times that in 75% diluent）. According to the discrepancy between ignition delay time and collision numbers, viz. 5 times corresponds to 40 times, the steric factor can
文摘The localized differential quadrature (LDQ) method is a numerical technique with high accuracy for solving most kinds of nonlinear problems in engineering and can overcome the difficulties of other methods (such as difference method) to numerically evaluate the derivatives of the functions.Its high efficiency and accuracy attract many engineers to apply the method to solve most of the numerical problems in engineering.However,difficulties can still be found in some particular problems.In the following study,the LDQ was applied to solve the Sod shock tube problem.This problem is a very particular kind of problem,which challenges many common numerical methods.Three different examples were given for testing the robustness and accuracy of the LDQ.In the first example,in which common initial conditions and solving methods were given,the numerical oscillations could be found dramatically;in the second example,the initial conditions were adjusted appropriately and the numerical oscillations were less dramatic than that in the first example;in the third example,the momentum equation of the Sod shock tube problem was corrected by adding artificial viscosity,causing the numerical oscillations to nearly disappear in the process of calculation.The numerical results presented demonstrate the detailed difficulties encountered in the calculations,which need to be improved in future work.However,in summary,the localized differential quadrature is shown to be a trustworthy method for solving most of the nonlinear problems in engineering.
基金supported by National Natural Science Foundation of China(Nos.62001095 and U20B2043)the Fundamental Research Funds for the Central Universities of China(No.ZYGX2018KYQD200)+1 种基金the National Science and Technology Major Project of the Ministry of Science and Technology of China(No.SQ2019YFA040012)the National Defense Basic Scientific Research Program of China(No.2020-JCJQ-ZD-072)。
文摘In this work,the theoretical analysis and experiment results investigating the influence of plasma-induced reflected wave variations on microwave transmission characterization are presented.Firstly,an analytical transmission line model for transmission characterization of plasma in shock tube is derived and validated against full-wave simulation.Then,the theoretical analysis of transmission characterization based on a time-dependent reconstruction algorithm that takes into account the variations of reflected wave is presented and the influence of reflection variations under various states of plasma is also investigated.The unusual increase in the amplitude of transmitted wave is theoretically predicted and experimentally demonstrated as well.Finally,the experiment results are also presented to illustrate the effects of reflected wave variations in practical microwave transmission characterization of supersonic plasma excited in shock tube.
文摘The purpose of this paper is to consider 1D Riemann shock tube to investigate the formation and propagation of compression waves leading to formation, propagation and reflection of 1D normal shocks using simplified mathematical models commonly used in the published work as well as using complete mathematical models based on Conservation and Balance Laws (CBL) of classical continuum mechanics and constitutive theories for compressible viscous medium derived using entropy inequality and representation theorem. This work is aimed at resolving compression waves, the shock structure, shock formation, propagation and reflection of fully formed shocks. Evolutions obtained from the mathematical models always satisfy differentiability requirements in space and time dictated by the highest order of the derivatives of the dependent variables in the mathematical models investigated. All solutions reported in this paper including boundary conditions and initial conditions are always analytic. Solutions of the mathematical models are obtained using the space-time finite element method in which the space-time integral forms are space-time variationally consistent ensuring unconditionally stable computations during the entire evolution. Solution for a space-time strip or slab is calculated and is time marched upon convergence to obtain complete evolution for the desired space-time domain, thus ensuring time accurate evolutions. The space-time local approximation over a space-time element of a space-time strip or slab is p-version hierarchical with higher-order global differentiability in space and time, i.e., we consider scalar product approximation spaces in which k = (k<sub>1</sub>, k<sub>2</sub>) are the order of the space in space and time and p = (p<sub>1</sub>, p<sub>2</sub>) are p-levels of the approximations in space and time. Model problem studies are presented for different mathematical models and are compared with solutions obtained from the complete mathematical model based on CBL and constitutive theories for viscous compressible medium to illustrate the deficiencies and shortcomings of the simplified and approximate models in simulating correct physics of normal shocks.
