Latest ArticlesBy investigating the JWL equation of state of detonation products of condensed explosive, we present a method to determine JWL parameters by γ fitting. This approach does not require cylinder test and is more economical, secure, convenient and accurate than existing methods. Using this method, four kinds of common explosive, e.g., TNT, C-4, PETN and HMX have been studied. By comparing to the p-V curve of JWL equation of state given by cylinder test, we showed that the γ fitting has a high precision and meets the need of explosion mechanics application.
In order to investigate the configuration and mass distribution of collision bulging, dynamic behaviors of Pb plate driven by collision of head-on sliding detonation waves were characterized with high-speed frame photography and pulsed X-ray radiography. Experimental results showed that a bulging similar to a shaped charge jet were produced in the collision region, which appeared as a multi-zone structure across several layers. The tip of bulging with higher velocity had a characteristic of fragmentary and cavitation because of its speed gradient and the interaction with air, and its bulk density was much lower than the pre-shock initial density of Pb. Though the root of bulging had a lower velocity and a higher density relatively, it was also in non-dense state. It suggested that the collision of sliding detonation waves, strength of materials and shock melting might play a leading role in the bulging formation process.
The analytic expressions for equation of state (EOS) and thermodynamic properties have been derived for the generalized Morse (GM) potential fluids, by using the Ross variational perturbation theory and the analytic Percus-Yevick (PY) expression of radial distribution function of hard spheres. Extensive comparison of the numerical results with computer simulations shows that the precision of the analytic Ross theory is equivalent to or slightly better than the non-analytic modified Weeks-Chandler-Anderson (mWCA) theory, and is much better than the complicated optimized reference hyper-netted chain (RHNC) theory. This analytic expressions for equation of state have been applied to N2 fluid at ambient temperature and low pressure (below 1 GPa) with the parameters of the generalized Morse potential are obtained by fitting to experimental data, results of prediction of the pressure above 1 GPa validate the analytic EOS with good agreement within a wide range of pressure and temperature.
The nano ZnO and nano composite ZnO/SnO2 materials were prepared by sol-gel method, then the prepared specimens were treated by high pressure method of which the pressure is up to 6 GPa. X-ray diffraction (XRD) and transmission electron microscopy (TEM) and photoluminescence spectrometer were employed to study structures and luminescence properties of the specimens. The results indicate that the particle sizes of ZnO and ZnO/SnO2 decrease 5.9% and 26.3% respectively after high pressure treatment. And the photoluminescence spectra show that photoluminescence spectra intensity of nano ZnO and composite ZnO/SnO2 materials both reduce after high pressure treatment, but the reduced intensity of composite ZnO/SnO2 materials is less than the one of pure ZnO, and the possible reasons are discussed.
High pressure Raman spectra of n-octane were investigated by diamond anvil cell (DAC) at pressure up to 7 GPa and ambient temperature. The result shows that the liquid-solid phase transition of n-octane takes place around 0.8 GPa. Within the experimental pressure range, the Raman shift of n-octane moves to higher frequency, while the peak intensity decreases, and the width of peaks broaden with increasing pressure. A turning point is observed in the curve of pressure-dependent Raman shift when the pressure increases up to 6.8 GPa. This result indicates that a solid-solid phase transition for n-octane has been finished at this pressure. The pressure-induced phase transition pressure in alkanes decreases with the length of their chains increasing.
Thermodynamic properties and structure of water under superhigh pressure were investigated by Equilibrium Molecular Dynamics (EMD) simulations. The obtained heat of vaporization and self-diffusion coefficient are in good agreement with the experimental data. With the increase of pressure, the hydrogen bond strength between water molecules increases, while self-diffusion coefficient decreases. With the increase of temperature, the hydrogen bond strength between water molecules decreases, while self-diffusion coefficient increases. However in the supercooled water, self-diffusion coefficient increases with the increase of temperature.
To study the state parameters of methane-air mixture under shock compression, Aluminum target chamber filled with methane-air mixture were impacted by tungsten alloy projectiles. Projectiles were accelerated to about 5 km/s using a two-stage light gas gun. Methane-air mixture with different air ratioes of 0, 1%, 5%, and 10% were tested, respectively. The radiance signals from shocked methane-air gases were recorded by an instantaneous pyrometer system. The Hugoniot data of four mixed ratio gases were obtained. The experimental results indicated that the shock temperatures of mixtures increase with rising of the air mixed ratio. The mixtures emitted non-equilibrium thermal radiation in shock downstream. The ionization of methane-air gas with 10% air mixed ratio is calculated by Saha ionization function. Calculated results indicated that methane-air mixed gases with an air mixed ratio less than 10% doesn't ionize, and can protect pins effectively.
Detonation shock dynamics (DSD) is an approximation to the reactive Euler equations that allows numerically efficient tracking of curved detonation waves. The DSD parameters are the velocity-curvature relation and the boundary angle. A computer code was developed to facilitate the calibration of these parameters for JB-9014 insensitive high explosive using the generalized optics model of DSD. Calibration data were obtained from measurements of the detonation velocities and fronts in JB-9014 rate sticks at ambient temperature, with diameters of 10~30 mm. The steady-state detonation velocities and fronts predicted by these DSD parameters are in very good agreement with experiment.
A graded-density flier plate with gradual change impedance was designed to launch the secondary flier plate in hypervelocity on a light gas gun. The impact surface of graded-density flier plate was made of magnesium, the base surface was made of tungsten, and the middle part was made of the mixture of magnesium, copper and tungsten. A calculating model of hypervelocity launch was developed. The numerical simulations of hypervelocity launch with graded-density flier plate were conducted. The effects of change law of impedance, diameter of secondary flier plate and thickness of buffer on the secondary flier plate velocity were analyzed. The calculated results indicated that the maximum velocity of the secondary flier plate exceeded 16 km/s in the condition of optimizing the change law of impedance of graded-density flier plate, diameter of secondary flier plate, and thickness of buffer.
The method of 1D elastic-plastic hydrodynamic calculations, the mixing model of different material constituents, and the shock dynamics code of compressible fluids were verified through the overall comparisons of the calculated results with the experimental data from Lawrence Livermore National Laboratory (LLNL). A Mg-W system prepared for graded density impactor (GDI) for quasi-isentropic compression was proposed by means of numerical simulations. By considering the actual material fabrication, the technical restrictions of minimum thickness of discrete layer, and the limitation of maximum mixing rate of the two individual constituents, a model GDI of a novel Mg-W system was designed to actualize the required loading traces and anchor the structural characteristics of GDI. Numerical calculations show that the loading path and intensity are highly sensitive to the distribution of wave impedance of GDI. The experimental designs and predictions for different complicated loading/unloading processes can be realized through model calculations.