Photochemical initiation of gaseous detonations

Direct initiation of detonations in gaseous mixtures of C₂H₂-O₂, H₂-O₂ and H₂-Cl₂ in the pressure range of 10–150 torr using flash photolysis was studied. Similar to blast initiation using a concentrated powerful energy source, it was found that for photochemical initiation, there exists a certain threshold of flash intensity and energy for each mixture at any given initial pressure and composition below which a deflagration is formed. At the critical threshold, however, a fully developed detonation is rapidly formed in the immediate vicinity of the window of incident UV radiation. However, at super critical flash energies, the amplitude of the detonation formed decreases and combustion of the entire irradiated volume approaches a constant volume explosion. It was found that photo-chemical initiation requires both a certain minimum peak value of the free radical concentration generated by the photo-dissociation as well as an appropriate gradient of this free radical distribution. The minimum peak radical concentration permits rapid reaction rates for the generation of strong pressure waves, while the gradient is necessary for the amplification of the shock waves to a detonation. If the gradient is absent and the free radicals are uniformly distributed in the mixture, then the entire volume simply explodes as in a constant volume process. The present study reveals that the mechanism of photochemical initiation is one of proper temporal synchronization of the chemical energy release to the shock wave as it propagates through the mixture. In analogy to the LASER, the term SWACER is introduced to represent this mechanism of Shock Wave Amplication by Coherent Energy Release. There are strong indications that this SWACER mechanism is universal and plays the main role in the formation of detonations whenever a powerful concentrated external source is not used to generate a strong shock wave in the explosive.     

Hydrodynamic instabilities and transverse waves in propagation mechanism of gaseous detonations

The present study examines the role of transverse waves and hydrodynamic instabilities mainly, Richtmyer–Meshkov instability (RMI) and Kelvin–Helmholtz instability (KHI) in detonation structure using two-dimensional high-resolution numerical simulations of Euler equations. To compare the numerical results with those of experiments, Navier–Stokes simulations are also performed by utilizing the effect of diffusion in highly irregular detonations. Results for both moderate and low activation energy mixtures reveal that upon collision of two triple points a pair of forward and backward facing jets is formed. As the jets spread, they undergo Richtmyer–Meshkov instability. The drastic growth of the forward jet found to have profound role in re-acceleration of the detonation wave at the end of a detonation cell cycle. For irregular detonations, the transverse waves found to have substantial role in propagation mechanism of such detonations. In regular detonations, the lead shock ignites all the gases passing through it, hence, the transverse waves and hydrodynamic instabilities do not play crucial role in propagation mechanism of such regular detonations. In comparison with previous numerical simulations present simulation using single-step kinetics shows a distinct keystone-shaped region at the end of the detonation cell.     

A study of detonation in methane/air clouds

In the present paper, we study the problem of detonation in unconfined, gaseous mixtures of methane/oxygen/nitrogen. A numerical simulation approach is employed in which we use a one-dimensional (spherical symmetry), time-dependent computer model to simulate the coupled compressible fluid dynamics-chemical kinetics processes which occur upon direct initiation of detonation. We establish the magnitude of explosive yield of tetryl required to initiate detonation in mixtures of CH₄ + 2O₂ with varying degrees of nitrogen dilution, up to and including stoichiometric . The numerical simulations illustrate the features of direct initiation observed in many experimental investigations, e.g. shock-wave breakaway followed by detonation reestablishment via a quasi-steady, oscillatory flow regime which occurs before the establishment of a steadily propagating spherical detonation. Our results compare well with recent experimental data obtained by Bull et al. (1976) over the range of tetryl masses studied by them. We find that tetryl explosive masses in excess of 10⁷ grams would be required to initiate detonation in an unconfined, stoichiometric mixture.     

The failure of marginal detonations in expanding channels

Critical conditions for detonation failure due to tube expansion have been observed in marginal detonations propagating in a in. (6.35 × 76.2 mm) channel. In these experiments, a well established marginal detonation propagating in the narrow channel entered a test section in which one of the narrow walls was inclined to the central axis at positive angles which ranged from 10° to 45°. Experiments were performed at pressures ranging from 60 to 200 torr (8 to 26.7 kPa) in stoichiometric hydrogen-oxygen mixtures diluted with 20, 50 and 70% argon. Smoke track records obtained on the surface which is the major dimension of the tube, were used to determine failure, incipient failure or self-sustenance of the entering wave.Because of the narrow tube used in the studies the incident waves were marginal in that their velocity was below the expected CJ (Chapman-Jouguet) value, their transverse wave spacing was larger than one would see in a large tube, and the transverse waves were of greater strength than in an ordinary detonation. All of these indicators of marginal behavior became progressively more pronounced as the pressure dropped from 200 torr (26.7 kPa) to the limit pressure of approximately 58 torr (7.73 kPa).The most interesting result of this experimental investigation is that the theoretical analyses predicted that simple one-dimensional opening of the tube should not show a pressure dependence to failure, while the experiments showed a definite decrease in the opening angle required for failure as initial pressure decreased. This behavior is related to the marginality of the incident waves, which is observed to increase smoothly with decreased pressure. It is postulated that detonation failure in the hydrogen-oxygen system occurs when the shock velocity at the end of the cell drops to about 0.60 of the CJ value due either to marginal behavior or to an expansion of the cross section of the tube.     

Buoyancy effects on smoldering combustion

A theoretical and experimental study is carried out to determine the effect of buoyancy on the rate of spread of a cocurrent smolder reaction through a porous combustible material. Since buoyant forces are proportional to the product g(_gi − _g), they can be controlled experimentally by varying either the gravitational acceleration, g, or the density difference, _gi − _g. The latter approach was followed in the present work. Measurements are performed of the smolder spread rate through porous α-cellulose (0.83 void fraction) as a function of the ambient air pressure. The experiments are carried out in a pressure vessel for ambient pressures ranging from 0.5 to 1.2 atm. The rate of spread was obtained from the temperature histories of thermocouples placed at fixed intervals along the fuel centerline. The smolder velocity was found to increase as the ambient pressure was increased. Extinction was found to occur when the buoyancy forces could not overcome the drag forces, indicating that at least for the present experimental conditions transport by diffusion cannot, by itself, support the spread of a smolder reaction. This conclusion is particularly important for outer space conditions where gravity and consequently buoyancy could be negligible. In the analysis, which assumes one-dimensional processes, the transport equations are solved to give the smolder spread rate as a function of the inlet oxygen mass flux. This mass flux is then estimated by balancing buoyancy and drag forces. Assuming that the smolder chemical reaction is only weakly dependent on pressure, the analysis finally predicts a smolder velocity dependence of the form v Y_oig²_gi Pa², i.e. is proportional to the ambient pressure squared. Good qualitative agreement is found between the theoretical predictions and the experimental results.     

A numerical study of the charge diameter effect on unsteady detonation

A numerical study for the unsteady detonation of an unconfined tetryl charge of small diameter, which is assumed to be homogeneous, was performed by using the two-dimensional Lagrangian hydrodynamic computer code, 2 DL, with the first order Arrhenius equation of reaction rate. Becker-Kistiakowsky-Wilson (BKW) and Kihara-Hikita (KH) equations of state have been applied to the detonation products.In the case of BKW, it is shown that the rarefaction waves propagating inward from the lateral surface make the reaction rate slow and give a curvature to the front. Then after an induction time, a strong initiation occurs in the reaction zone near the lateral surface and higher pressure zone moves to the axis. This higher pressure accelerates the detonation propagation near the lateral surface and the curvature of detonation front is reduced. Then, the reaction at the lateral surface again begins to decay by the rarefaction waves. Such a sequence of process is repeated periodically.The possibility of the occurrence of the strong initiation depends on the pressure and temperature in the shocked zone near the surface. In a small diameter charge, the delayed explosion becomes weaker near the surface, while its frequency increases. No shock interaction occurs because the direction of the particle flow is always divergent.In the case of KH equation of state, the temperature of detonation is higher than that obtained by BKW and the behaviour of instability becomes rather different from the previous result, i.e. in the axis the pressure oscillates repeating the overdriven and underdriven detonation similar with the case of BKW.     

Fluid mechanics of pulse detonation thrusters

The advantages of constant volume combustion cycle as compared to constant pressure combustion in terms of thermodynamic efficiency have focused the researches of advanced propulsion on detonation engines. The paper gives coverage of efforts undertaken during past decades in adjusting detonations for propulsion applications, and highlights new challenges in studying fluid flow dynamics relevant to onset of detonation.     

煤油温度对于爆震波形成影响的实验研究

在内径为30mm的脉冲爆震发动机模型上，以煤油为燃料，以空气为氧化剂，成功地进行了两相爆震实验，获得了充分发展的脉冲爆震波。测试了在化学恰当比，不同爆震频率及燃油温度下的爆震波压力，并对其变化进行了分析。通过分析实验结果发现，在化学恰当比下，爆震频率不变时，煤油温度的升高明显促进了爆震的形成，在内径小于混合物胞格尺寸的爆震管内，可以形成充分发展的两相脉冲爆震波。     

Detonation in two-layer systems

Detonation in heterogeneous systems involving parallel layers of two different substances was investigated. It was assumed that the denser medium can undergo fast chemical reaction without mixing with the other medium, and the detonation propagates along the two layers. The process of initiation is associated with a shock wave advancing through the medium of lesser density. An idealized gasdynamic model of the phenomenon is proposed in order to evaluate detonation parameters of the two-layer system on the basis of the characteristic features of its constituents. As calculation shows, in the absence of mixing between the two layers, the detonation is capable of propagating at a higher velocity than in the case when the constituents forming the layers are mixed.     

Analysis of Low-Frequency Microaccelerations onboard the Foton-11Satellite

The method and the results of investigating the low-frequency component of microaccelerations onboard the Foton-11satellite are presented. The investigation was based on the processing of data of the angular velocity measurements made by the German system QSAM, as well as the data of measurements of microaccelerations performed by the QSAM system and by the French accelerometer BETA. The processing was carried out in the following manner. A low-frequency (frequencies less than 0.01 Hz) component was selected from the data of measurements of each component of the angular velocity vector or of the microacceleration, and an approximation was constructed of the obtained vector function by a similar function that was calculated along the solutions to the differential equations of motion of the satellite with respect to its center of mass. The construction was carried out by the least squares method. The initial conditions of the satellite motion, its aerodynamic parameters, and constant biases in the measurement data were used as fitting parameters. The time intervals on which the approximation was constructed were from one to five hours long. The processing of the measurements performed with three different instruments produced sufficiently close results. It turned out to be that the rotational motion of the satellite during nearly the entire flight was close to the regular Eulerian precession of the axially symmetric rigid body. The angular velocity of the satellite with respect to its longitudinal axis was about 1 deg/s, while the projection of the angular velocity onto the plane perpendicular to this axis had an absolute value of about 0.2 deg/s. The magnitude of the quasistatic component of microaccelerations in the locations of the accelerometers QSAM and BETA did not exceed 5 × 10^–5–10^–4m/s²for the considered motion of the satellite.     

High speed gas flows in explosions of cavitated explosives

A possibility of attaining steady flow of detonation products with specific energy much larger than the specific chemical energy of explosive is demonstrated in the case when a cylindrical charge of explosive is fitted with an evacuated cavity. Simple estimates and results of numerical analysis of the process are presented. Steady process may be considered to occur under the following assumptions: (1) effects arising due to jet interaction with cavity walls are negligible; (2) the detonation process is steady. In the case of limited explosive lengths these assumptions have been shown to be correct.When the cavity is filled with gas or liquid, a variety of steady and non-steady flow regimes is possible, depending on the density of the filling medium. One well-known case is that of flow with irregular reflection of shock waves at the cavity axis accompanied with the formation of Mach intersections. Another interesting flow regime is observed to occur in the case of low density filling medium (liquid hydrogen, for example). In this case the filling medium is driven by a “detonation piston” at constant velocity, equal to the velocity of detonation, forming a uniform growing column of hot shock-compressed matter, specific energy of which exceeds by one order of magnitude the specific energy of the explosive. Obviously, the walls of the vessel containing hydrogen must be able to withstand radial loads for a sufficiently long time (20 μ sec).The relative merits of these methods in comparison to others in high speed gas-dynamics is discussed.     

Some aspects of blast from fuel-air explosives

Numerical solutions of the blast wave flow from a spherical explosive charge were obtained using the artificial viscosity technique as employed by Oppenheim. The flow is treated as adiabatic and inviscid and ideal equations of state are used for reactants, products and the surrounding air environment. Differences are noted in the peak pressure, static impulse and dynamic impulse resulting from three representative types of idealized initiation: (1) centrally initiated, self-similar Chapman-Jouguet detonation, (2) edge initiated spherical implosion and (3) constant volume energy release followed by sudden venting to the environment. These are compared to the ideal point blast with counterpressure of the same total energy release. In addition, numerical solutions are presented for centrally initiated, stoichiometric hydrogen-oxygen mixtures surrounded by air for detonation and for deflagration according to an empirically determined, non-steady flame velocity. The greater energy transfer to the environment in the case of detonation is demonstrated.     

Cylindrical heterogeneous detonation waves

Further experimental studies of blast wave initiated cylindrical heterogeneous (liquid fuel drops, gas oxidizer) detonation waves are described. A pie-shaped shock tube, used for these studies, was altered in certain ways so as to improve the modeling of cylindrical waves. These modifications, along with some operational aspects, are briefly discussed. The breech of the facility, where the blast wave is generated by an explosive, became distorted with usage. Results are presented which show that lower detonation velocities are realized with the damaged breech (other conditions being the same). A photographic and pressure switch wave time of arrival study was made to ascertain the wave shape. Photographs are shown which show that the waves, blast as well as detonation, are close to cylindrical. However, in some cases there is appreciable distortion of the wave front by debris ahead of the wave. Presumably this debris comes from the blasting cap used to ignite the condensed explosive. A series of experiments was conducted using kerosene drops of 388 μm diameter dispersed in air through use of a large number of hypodermic needles. Radial fuel void regions were established by cutting off the fuel flow to a number of needles. Preliminary results relating to the effect of the size of the cloud gap on detonation velocity, quenching, and the initiator energy levels required for detonation are discussed.     

Low-velocity detonations in cast and pressed high explosives

A model of low velocity detonations in charges of cast and pressed high explosives confined in metal tubes with thin walls is suggested. An analytical solution is obtained and velocities and pressures of stationary LVD waves are calculated as functions of the parameters of the wall confinement. The fraction of the total chemical energy transferred to maintain the shock wave and the average rate of reaction in the wave are estimated and the dynamics of transient predetonation processes are analysed. An explanation of the causes of the secondary shock wave formation in the transition from deflagration to normal detonations is suggested and the delay of this secondary wave is calculated and compared with experimental data.     

Effects of a hot jet on detonation initiation and propagation in supersonic combustible mixtures

The effects of a hot jet on detonation initiation and propagation in supersonic combustible mixtures has been studied with two-dimensional numerical simulations with the open-source program AMROC that uses a block-structured adaptive mesh refinement method. Results indicate that the hot jet could ignite the detonation effectively in supersonic combustible mixtures like a pneumatic ramp. After the realization of the detonation initiation, the hot jet can still play an important role on the detonation propagation during its continuous ejection. For a hot jet with certain diameter, it can result in an overdriven detonation with almost constant overdrive degree. After the shutdown of the hot jet, the stable CJ detonation combustion was realized finally in the supersonic combustible mixtures. With the re-ejection of the hot jet, the failed detonation could be reinitiated quickly. Through the control of the re-ejection of the hot jet, it plays a key role not only in the initiation process, but also in the subsequent continuous detonation combustion period.     

弹丸速度对超爆轰冲压加速器性能影响

为研究超爆轰模态冲压加速器的推进性能，采用混合的Roe/HLL（Harten, Lax, Van Leer）格式，结合自适应网格加密技术(AMR )与沉浸边界法(IBM )，数值模拟了弹丸速度高于预混可燃气体C-J爆速的冲压加速器流场，揭示了弹丸速度对流场结构与推力的影响。结果表明当弹丸速度在一定范围时，斜爆轰波可驻定在弹丸肩部或头部，在弹丸尾部形成高压区加速弹丸，并且，斜爆轰波驻定在弹丸头部推力更高，稳定工作的速度范围 更宽 。     

返回舱再入过程密封舱气体泄漏计算研究

为分析返回舱再入过程中密封舱漏孔内外压差,并对漏孔变流量充气过程进行研究,采用离散化分析方法将返回舱再入过程分成若干个阶段,针对容积为14 m~3的密封舱和面积为10 cm~2的漏孔,计算并获得了密封舱内外压差、漏孔质量流率、漏孔流速等参数在50~5 km范围内随高度下降的变化规律。结果表明:在高度5 km开伞时刻,漏孔质量流率达到最大值0.134 kg/s,舱内外压差趋近于最大值,约20 172 Pa;返回舱下降过程中漏孔流速在148.4~181.5 m/s之间,处于亚声速区;漏孔气体流速与漏孔面积大小无关,仅与漏孔内外压力及漏孔进口空气密度有关。以上研究结果可为密封舱结构强度设计、伞舱弹伞设计提供参考。     

Shock compression of liquid hydrogen in various experimental geometries

In connection with the use of cryogenic liquids in high-speed gas dynamics and high-pressure physics, shock-wave processes in liquid hydrogen were investigated under plane, cylindrical and hemispherical loading.The plane loading of liquid hydrogen consisted of a multicyclic, nearly isentropic compression. A transducer employing a contact electrical effect was used to record this multicyclic compression process between a rigid wall and a flyer, resulting in a sequence of shock steps of decreasing amplitude, whose integrated action is equivalent to the isentropic compression of liquid hydrogen up to 500 kbar.The cylindrical loading was generated by detonating a high-explosive charge enclosing a cylindrical cavity along its axis that was filled with liquid hydrogen. Under these conditions shock velocities up to 13.7 km/sec were recorded, and pressure in the shock-compressed hydrogen reached 90 kbar. The formation of a boundary layer and expansion of the cylindrical cavity limited further pressure increases in the column of compressed liquid and lead to a decrease in the flow velocity. The observed increase in detonation velocity is associated with the influence of the channel wave on the detonation regime in the neighboring explosive layers.Under hemispherical loading, an increase in the converging shock velocity from 6 to 20 km/sec was recorded. The final pressure reached 210 kbar, and the total specific energy exceeded 200 kJ/g. During the release of the shock-compressed hydrogen into air at 0.1 torr, shock waves with velocities exceeding 50 km/sec were obtained.     

固体火箭发动机撞击靶板安全性数值分析

为研究固体火箭发动机撞击安全性,建立了固体火箭发动机撞击靶板的计算模型,模型中发动机的推进剂装药采用点火增长反应速率方程.采用非线性有限元流体动力学方法,对发动机径向撞击靶板过程进行了数值模拟,分析了不同撞击速度下发动机中推进剂装药的反应情况.计算结果表明,发动机径向撞击靶板爆炸的临界速度范围为150～200 m/s;低强度多次撞击过程中推进剂会发生延迟爆轰情况.     

超声速斜爆震发动机起爆过程研究综述

对超声速斜爆震发动机的起爆方式进行了比较分析,对起爆发展和稳定特性的研究历程和发展现状进行了综述,对相关的研究方法和技术进行了概括,提出了利用先进光学测量技术,结合激光诱导荧光技术对超声速斜爆震发动机起爆过程进行实验研究的设想。