Numerical modelling of non-ideal detonation in ANFO explosives applying Wood-Kirkwood theory coupled with EXPLO5 thermochemical code

Ammonium nitrate and fuel oil(ANFO)based explosive is a classic example of non-ideal high explosives.Its detonation is characterized by a strong dependence of detonation parameters on explosive charge diameter,presence and characteristics of confinement,as well as incomplete consumption of explosive at the sonic point.In this work we propose a detonation model based on the Wood-Kirkwood(WK)theory coupled with the thermochemical code EXPLO5 and supplemented with reaction rate models.Our objective is to analyze the validity of the model for highly non-ideal ANFO explosives,with emphasis on effect of reaction rate models.It was found that both single-step and two-step pressure-based models can be calibrated to reproduce experimental detonation velocity-charge radius data of ANFO at radii significantly above the failure radius(i.e.for D/D_(id)>~0.6).Single-step pressure-based model,with the pressure exponent equal to 1.4,proved to be the most accurate,even in the vicinity of the failure radius.The impact of the rate models is most evident on temporal(and spatial)distribution of flow parameters in detonation driving zone,especially when it comes to the conversion and width of detonation driving zone.

Confinement Effect and Diameter Effect on the Non-ideal Detonation Parameters of Small Charge

An expanding model of the confinement of non-ideal detonation of small charge is established on the basis of the nozzle theory.Making use of the expanding model,the analytic relationship of small charge detonation velocity and the semi-empirical relationship of detonation pressure that both change with charge diameter and confinement condition are established.The detonation velocity and pressure of small charges are calculated and experimentally verified,and the detonation velocity deviation is less than 7% while the detonation pressure deviation is less than 9%.

Assessing the energy release characteristics during the middle detonation reaction stage of aluminized explosives

Afterburning behind the detonation front of an aluminized explosive releases energy on the millisecond timescale,which prolong the release of detonation energy and the energy release at different stages also shows significant differences.However,at present,there are few effective methods for evaluating the energy release characteristics of the middle reaction stage of such explosives,which can have a duration of tens to hundreds of microseconds.The present work demonstrates an approach to assessing the midstage of an aluminized explosive detonation based on a water push test employing a high degree of confinement.In this method,the explosive is contained in a steel cylinder having one end closed that is installed at the bottom of a transparent water tank.Upon detonation,the gaseous products expand in one direction while forcing water ahead of them.The resulting underwater shock wave and the interface between the gas phase products and the water are tracked using an ultra-high-speed framing and streak camera.The shock wave velocity in water and the expansion work performed by the gaseous detonation products were calculated to assess the energy release characteristics of aluminized explosives such as CL-20 and RDX in the middle stage of the detonation reaction.During the middle stage of the detonation process of these aluminized explosives,the aluminum reaction reduced the attenuation of shock waves and increased the work performed by gas phase products.A higher aluminum content increased the energy output while the presence of oxidants slowed the energy release rate.This work demonstrates an effective means of evaluating the performance of aluminized explosives.

Reheat effect on the improvement in efficiency of the turbine driven by pulse detonation

Due to the strong unsteadiness of pulse detonation,large flow losses are generated when the detonation wave interacts with the turbine blades,resulting in low turbine efficiency.Considering that the flow losses are dissipated into the gas as heat energy,some of them can be recycled during the expansion process in subsequent stages by the reheat effect,which should be helpful to improve the detonationdriven turbine efficiency.Taking this into account,this paper developed a numerical model of the detonation chamber coupled with a two-stage axial turbine,and a stoichiometric hydrogen-air mixture was used.The improvement in turbine efficiency attributable to the reheat effect was calculated by comparing the average efficiency of the stages with the efficiency of the two-stage turbine.The research indicated that the first stage was critical in suppressing the flow unsteadiness caused by pulse detonation,which stabilized the intake condition of the second stage and consequently allowed much of the flow losses from the first stage to be recycled,so that the efficiency of the two-stage turbine was improved.At a 95%confidence level,the efficiency improvement was stable at 4.5%—5.3%,demonstrating that the reheat effect is significant in improving the efficiency of the detonation-driven turbine.

Study on concentration distribution and detonation characteristics for non-axisymmetric fuel dispersal

The study of non-axisymmetric fuel dispersal and detonation can provide reference for the prevention of industrial cloud explosion accidents and the design of fuel air explosive(FAE).The concentration and detonation fields of 85 kg cylindrical and fan-shaped fuel are investigated by experiments and numerical simulations.A dynamic model of the whole process for fuel dispersal and detonation is built.The concentration distribution of fuel is used as the initial condition to calculate the detonation stage,thus solving the initial value problem of detonation field.The phase and component changes of fuel cloud at different locations are compared.The fuel cloud is divided into directions of 0°,90°,135°and 180°.The results show that the maximum cloud radius is 20.94 m in 135°and the minimum is 12.04 m in 0°.The diameter of the detonation fireball is 53.6 m,and the peak temperature is 3455 K.The highest peak overpressure is 3.44 MPa in 0°and the lowest is 2.97 MPa in 135°.The proportion of liquid phase in 0°is22.90%,and the fuel loss is 11.8% and 9% higher than that in 135°and cylindrical charge,respectively.The stable propagation distance of blast wave in 135°is 42.50% longer than 0°and 28.37% longer than cylindrical charge.

Atmospheric transmission algorithm for pulsed X-rays from high-altitude nuclear detonations based on scattering correction

In high-altitude nuclear detonations,the proportion of pulsed X-ray energy can exceed 70%,making it a specific monitoring signal for such events.These pulsed X-rays can be captured using a satellite-borne X-ray detector following atmospheric transmission.To quantitatively analyze the effects of different satellite detection altitudes,burst heights,and transmission angles on the physical processes of X-ray transport and energy fluence,we developed an atmospheric transmission algorithm for pulsed X-rays from high-altitude nuclear detonations based on scattering correction.The proposed method is an improvement over the traditional analytical method that only computes direct-transmission X-rays.The traditional analytical method exhibits a maximum relative error of 67.79% compared with the Monte Carlo method.Our improved method reduces this error to within 10% under the same conditions,even reaching 1% in certain scenarios.Moreover,its computation time is 48,000 times faster than that of the Monte Carlo method.These results have important theoretical significance and engineering application value for designing satellite-borne nuclear detonation pulsed X-ray detectors,inverting nuclear detonation source terms,and assessing ionospheric effects.

Research of detonation products of RDX/Al from the perspective of composition

Aluminized explosives exhibit excellent performance because the oxidation of aluminum(Al)powders enhances the pressure and temperature of detonation products.However,the equation of state(EOS)of detonation products has not been understood well.In the present study,we conducted long-time tests(approximately 1 ms)of a metal rod driven by detonation products of RDX,RDX/Li F,and RDX/Al.In addition,we used laser velocimetry(PDV)to measure the freesurface velocity of the rod.Thermochemical code DLCHEQ was successfully applied to the hydrodynamic program SSS to perform the roddriven test,and a novel method was established to study the EOS of detonation products from the perspective of composition.The reliability of DLCEHQ was validated by a small deviation(<10%)between the experimental rod free-surface velocity of RDX and the calculated results;the deviation was considerably less than that from the results obtained using the JWL EOS and ideal-gas EOS.The endothermic process and the reaction of Al powders(Al+H_(2)O+NO+CO_(2)→CO+H_(2)+N_(2)+Al_(2)O_(3))were analyzed by calculating the rod free-surface velocity of RDX/Li F and RDX/Al,respectively.The results of the present study demonstrated that the thermodynamic state of Al powders has notable influence on the EOS of aluminized detonation products,and the findings were compared with those of previous studies.First,the temperature equilibrium between Al powders and CHNO products is not always achieved,and the disequilibrium is more obvious when the reaction of Al powders is stronger.Second,the reaction rate of Al powders depends on pressure and Al content.Finally,the endothermic process of Al powders has a high contribution to the decrease in the work ability of RDX/Al instead of the gasconsumption mechanism of the Al reaction.More than half of the reaction heat of Al powders is used to heat itself,whereas the gas consumption during the reaction is negligible.

Deflagration to detonation transition in weakly confined conditions for a type of potentially novel green primary explosive:Al/Fe_(2)O_(3)/RDX hybrid nanocomposites

The properties of the combustion and deflagration to detonation transition(DDT)of Al/Fe_(2)O_(3)/RDX hybrid nanocomposites,a type of potentially novel lead-free primary explosives,were tested in weakly confined conditions,and the interaction of Al/Fe_(2)O_(3)nanothermite and RDX in the DDT process was studied in detail.Results show that the amount of the Al/Fe_(2)O_(3)nanothermite has a great effect on the DDT properties of Al/Fe_(2)O_(3)/RDX nanocomposites.The addition of Al/Fe_(2)O_(3)nanothermite to RDX apparently improves the firing properties of RDX.A small amount of Al/Fe_(2)O_(3)nanothermite greatly increases the initial combustion velocity of Al/Fe_(2)O_(3)/RDX nanocomposites,accelerating their DDT process.When the contents of Al/Fe_(2)O_(3)nanothermite are less than 20 wt%,the DDT mechanisms of Al/Fe_(2)O_(3)/RDX nanocomposites follow the distinct abrupt mode,and are consistent with that of RDX,though their DDT processes are different.The RDX added into the Al/Fe_(2)O_(3)nanothermite increases the latter's peak combustion velocity and makes it generate the DDT when the RDX content is at least 10 wt%.RDX plays a key role in the shock compressive combustion,the formation and the properties of the DDT in the flame propagation of nanocomposites.Compared with RDX,the fast DDT of Al/Fe_(2)O_(3)/RDX nanocomposites could be obtained by adjusting the chemical constituents of nanocomposites.

High-Order Bound-Preserving Finite Difference Methods for Multispecies and Multireaction Detonations

In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical reaction,say zi,is between 0 and 1,withΣz_(i)=1.Due to the lack of maximum-principle,most of the previous bound-preserving technique cannot be applied directly.To preserve those bounds,we will use the positivity-preserving technique to all the zi'is and enforceΣz_(i)=1 by constructing conservative schemes,thanks to conservative time integrations and consistent numerical fluxes in the system.Moreover,detonation is an extreme singular mode of flame propagation in premixed gas,and the model contains a significant stiff source.It is well known that for hyperbolic equations with stiff source,the transition points in the numerical approximations near the shocks may trigger spurious shock speed,leading to wrong shock position.Intuitively,the high-order weighted essentially non-oscillatory(WENO)scheme,which can suppress oscillations near the discontinuities,would be a good choice for spatial discretization.However,with the nonlinear weights,the numerical fluxes are no longer“consistent”,leading to nonconservative numerical schemes and the bound-preserving technique does not work.Numerical experiments demonstrate that,without further numerical techniques such as subcell resolutions,the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.

Memristor’s characteristics: From non-ideal to ideal

Memristor has been widely studied in the field of neuromorphic computing and is considered to be a strong candidate to break the von Neumann bottleneck. However, the non-ideal characteristics of memristor seriously limit its practical application. There are two sides to everything, and memristors are no exception. The non-ideal characteristics of memristors may become ideal in some applications. Genetic algorithm(GA) is a method to search for the optimal solution by simulating the process of biological evolution. It is widely used in the fields of machine learning, combinatorial optimization,and signal processing. In this paper, we simulate the biological evolutionary behavior in GA by using the non-ideal characteristics of memristors, based on which we design peripheral circuits and path planning algorithms based on memristor networks. The experimental results show that the non-ideal characteristics of memristor can well simulate the biological evolution behavior in GA.

Dynamic analysis of buried pipeline with and without barrier system subjected to underground detonation

Failure of pipe networks due to blast loads resulting from terrorist attacks or construction facilities, may cause economic loss, environmental pollution, source of firing or even it may lead to a disaster. The present work develops a closed-form solution of buried pipe with barrier system subjected to subsurface detonation. The solution is derived based on the concept of double-beam system. Euler Bernoulli's beams are used to simulate the buried pipe and the barrier system. Soil is idealized as viscoelastic foundation along with shear interaction between discrete Winkler springs(advanced soil model). The finite SineFourier transform is employed to solve the coupled partial differential equations. The solution is validated with past studies. A parametric study is conducted to investigate the influence of TNT charge weight, pipe material, damping ratio and TNT offset on the response of buried pipe with and without barrier system. Further a statistical analysis is carried out to get the significant soil and pipe input parameters. It is perceived that peak pipe displacements for both the cases(with and without barrier) are increases with increasing the weight of TNT charge and decreases with increasing the damping ratio and TNT offset. The deformation of pipe also varies with pipe material. Pipe safety against blast loads can be ensured by providing suitable barrier layer. The present study can be utilized in preliminary design stage as an alternative to expensive numerical analysis or field study.

Numerical simulation of flow field characteristics and the improvement of pressure oscillation of rotating detonation engine

Due to the inherent working mode of rotating detonation engine(RDE),the detonation flow field has the characteristics of pressure oscillation and axial kinetic energy loss,which makes it difficult to design nozzle and improve propulsion performance.Therefore,in order to improve the characteristics of detonation flow field,the three-dimensional numerical simulation of annular chamber and hollow chamber is carried out with premixed hydrogen/air as fuel in this paper,and then tries to combine the two chambers to weaken the oscillation characteristics of detonation flow field through the interaction of detonation flow field,which is a new method to regulate the detonation flow field.The results show that there are four states of velocity vectors at the outlet of annular chamber and hollow chamber,which makes RDE be affected by rolling moment and results in the loss of axial kinetic energy.In the external flow field of combined chamber,the phenomenon of cyclic reflection of expansion wave and compression wave on the free boundary is observed,which results in Mach disk structure.Moreover,the pressure monitoring points are set at the external flow field.The pressure signal shows that the high-frequency pressure oscillation at the external flow field of the combined chamber has been greatly weakened.Compared to the annular chamber,the relative standard deviation(RSD) has been reduced from 14.6% to5.6%.The results thus demonstrate that this method is feasible to adjust the pressure oscillation characteristics of the detonation flow field,and is of great significance to promote the potential of RDE and nozzle design.

Effect of spin on the instability of THz plasma waves in field-effect transistors under non-ideal boundary conditions

Terahertz(THz) radiation can be generated due to the instability of THz plasma waves in field-effect transistors(FETs). In this work, we discuss the instability of THz plasma waves in the channel of FETs with spin and quantum effects under non-ideal boundary conditions. We obtain a linear dispersion relation by using the hydrodynamic equation, Maxwell equation and spin equation. The influence of source capacitance, drain capacitance, spin effects, quantum effects and channel width on the instability of THz plasma waves under the non-ideal boundary conditions is investigated in great detail. The results of numerical simulation show that the THz plasma wave is unstable when the drain capacitance is smaller than the source capacitance;the oscillation frequency with asymmetric boundary conditions is smaller than that under non-ideal boundary conditions;the instability gain of THz plasma waves becomes lower under non-ideal boundary conditions. This finding provides a new idea for finding efficient THz radiation sources and opens up a new mechanism for the development of THz technology.

Study on the prediction and inverse prediction of detonation properties based on deep learning

The accurate and efficient prediction of explosive detonation properties has important engineering significance for weapon design.Traditional methods for predicting detonation performance include empirical formulas,equations of state,and quantum chemical calculation methods.In recent years,with the development of computer performance and deep learning methods,researchers have begun to apply deep learning methods to the prediction of explosive detonation performance.The deep learning method has the advantage of simple and rapid prediction of explosive detonation properties.However,some problems remain in the study of detonation properties based on deep learning.For example,there are few studies on the prediction of mixed explosives,on the prediction of the parameters of the equation of state of explosives,and on the application of explosive properties to predict the formulation of explosives.Based on an artificial neural network model and a one-dimensional convolutional neural network model,three improved deep learning models were established in this work with the aim of solving these problems.The training data for these models,called the detonation parameters prediction model,JWL equation of state(EOS)prediction model,and inverse prediction model,was obtained through the KHT thermochemical code.After training,the model was tested for overfitting using the validation-set test.Through the model-accuracy test,the prediction accuracy of the model for real explosive formulations was tested by comparing the predicted value with the reference value.The results show that the model errors were within 10%and 3%for the prediction of detonation pressure and detonation velocity,respectively.The accuracy refers to the prediction of tested explosive formulations which consist of TNT,RDX and HMX.For the prediction of the equation of state for explosives,the correlation coefficient between the prediction and the reference curves was above 0.99.For the prediction of the inverse prediction model,the prediction error of the explosive equation was within 9%.This indicates that the models have utility in engineering.

A numerical study of two-and three-dimensional detonation dynamics of pulse detonation engine by the CE/SE method

In this paper, the CE/SE method is developed to simulate the two- and three-dimensional flow-field of Pulse Detonation Engine (PDE). The conservation equations with stiff source terms for chemical reaction are solved in two steps. The detailed analysis of computational results of a PDE with a single detonation tube and a PDE with five detonation tubes are given in this paper. Complex wave systems are observed inside and outside a PDE. For a PDE with 5 detonation tubes, there is a big bow shock produced from a number of little shocks near the open ends of tubes. A lot of vortexes interact with shocks and a large expansion wave propagates forward and backward with respect to the PDE in a semi-oval shape.

Verification and validation of detonation modeling

The mathematical model used to describe the detonation multi-physics phenomenon is usually given by highly coupled nonlinear partial differential equations. Numerical simulation and the computer aided engineering (CAE) technique has become the third pillar of detonation research, along with theory and experiment, due to the detonation phenomenon is difficult to explain by the theoretical analysis, and the cost required to accredit the reliability of detonation products is very high, even some physical experiments of detonation are impossible. The numerical simulation technique can solve these complex problems in the real situation repeatedly and reduce the design cost and time stunningly. But the reliability of numerical simulation software and the serviceability of the computational result seriously hinders the extension, application and the self-restoration of the simulation software, restricts its independently innovational ability. This article deals with the physical modeling, numerical simulation, and software development of detonation in a unified way. Verification and validation and uncertainty quantification (V&V&UQ) is an important approach in ensuring the credibility of the modeling and simulation of detonation. V&V of detonation is based on our independently developed detonation multiphysics software-LAD2D. We propose the verification method based on mathematical theory and program function as well as availability of its program execution. Validation is executed by comparing with the experiment data. At last, we propose the future prospect of numerical simulation software and the CAE technique, and we also pay attention to the research direction of V&V&UQ.

Wave dynamic processes in cellular detonation reflection from wedges

When the cell width of the incident detonation wave （IDW） is comparable to or larger than the Mach stem height, self-similarity will fail during IDW reflection from a wedge surface. In this paper, the detonation reflection from wedges is investigated for the wave dynamic processes occurring in the wave front, including transverse shock motion and detonation cell variations behind the Mach stem. A detailed reaction model is implemented to simulate two-dimensional cellular detonations in stoichiometric mixtures of H2/O2 diluted by Argon. The numerical results show that the transverse waves, which cross the triple point trajectory of Mach reflection, travel along the Mach stem and reflect back from the wedge surface, control the size of the cells in the region swept by the Mach stem. It is the energy carried by these transverse waves that sustains the triple-wave-collision with a higher frequency within the over-driven Mach stem. In some cases, local wave dynamic processes and wave structures play a dominant role in determining the pattern of cellular record, leading to the fact that the cellular patterns after the Mach stem exhibit some peculiar modes.

Application of Two-dimensional Viscous CE/SE Method in Calculation of Two-phase Detonation

The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.

Formation and corrosion behavior of Fe-based amorphous metallic coatings prepared by detonation gun spraying

Amorphous metallic coatings with a composition of Fe48Cr15Mo14C15B6Y2 were prepared by detonation gun spraying process. Microstructural studies show that the coatings present a densely layered structure typical of thermally sprayed deposits with the porosity below 2%. Both crystallization and oxidation occurred obviously during spraying process, so that the amorphous fraction of the coatings decreased to 54% compared with fully amorphous alloy ribbons of the same component. Corrosion behavior of the amorphous coatings was investigated by electrochemical measurement. The results show that the coatings exhibit extremely wide passive region and low passive current density in 3.5% NaCl (mass fraction) and 1 mol/L HCl solutions, which illustrates excellent ability to resist localized corrosion.

Periodic oscillation and fine structure of wedge-induced oblique detonation waves

An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to capture the shock wave.The numerical results show that with the compression of the wedge wall the detonation wave front structure is divided into three sections：the ZND model-like strcuture,single-sided triple point structure and dual-headed triple point strucuture.The first structure is the smooth straight,and the second has the characteristic of the triple points propagating dowanstream only with the same velocity,while the dual-headed triple point structure is very complicated.The detonation waves facing upstream and downstream propagate with different velocities,in which the periodic collisions of the triple points cause the oscillation of the detonation wave front.This oscillation process has temporal and spatial periodicity.In addition,the triple point trace are recorded to obtain different cell structures in three sections.