A three-dimensional two-phase detonation model for aluminum dusts/air mixtures is built. A 3D numerical simulation program for detonation of suspended Al dusts is developed with space-time conservation element and solution element (CE/SE) method. Moreover, program parallelization is realized based on message passing interface (MPI) technique. Reliability of the program is demonstrated with simulation of shock tube problem and two-phase detonation experiment of Al particles/air mixtures in a detonation tube. Two-phase detonation of 368 g·m^-3 Al particles/air mixtures in left space of a corner and its aftereffect on air in the right and lower space of the corner are investigated. Propagation, reflection and diffraction of detonation wave or shock wave in complex space are obtained. It shows that air shock wave generated by two-phase detonation can reach a reflection pressure of 2.66 MPa at the solid wall 2 m far away from the Al dusts region. Fireball range exceeds initial Al dusts region by about 0.8 m. Air temperature within a range of 1.5 m near the initial Al dusts region is more than 1 600 K. The program can be used for aftereffect study of Al dusts detonation. It provides guidance for industrial safety and protection.

By introducing conductive burning process into classical deflagration-to-detonation transition (DDT) model, transition process from low speed conductive burning to convective burning to detonation was proposed. Transition process in HMX granular bed with 85% loading density was simulated. Development of conductive burning, convective burning and detonation was analyzed. In early stage combustion propagation rate is very slow. It propagates no more than 0. 2mm within 8. 16ms. After onset of convective burning, it tooks 20 ms to form a steady detonation with a velocity of 8 165 m·s^-1. Time to form detonation increases with decrease of particle diameter and ignition pressure.

Two-phase detonation of RDX particles suspended in air was numerically studied with CE/SE method.Behind leading shock front of detonation,explosive particles are accelerated and heated by the gas flow.Energy is released to support propagation of detonation wave.Dust detonation in a shock tube was numerically simulated.Distribution of physical quantity behind leading shock front was calculated.Parameters of detonation were obtained and they agree well with those in a reference.Dust detonation in a complex channel was numerically simulated.It shows that CE/SE method simulates gas-solid two-phase detonation successfully.

Theoretical study and numerical simulations are made on size distribution characteristic of initial liquid drops during aerosol explosive dispersal process.With thermodynamic consistent concept,the maximization of entropy generation in the initial breakup process,constraint equations of the problem are derived,and probability distribution function of initial drop size is obtained.Numerical predictions are in agreement with air-blast annular nozzle experimental data and Samirant experimental data.

High-energy explosive detonation diffraction is simulated with a 2D Lagrangian unstructured mesh detonation fluid program. With EOS of JWL and trinomial reaction rate, numerical result by Lagrangian program agrees with experimental result qualitatively. Phenomena of whirlpool is well simulated. Desensitization is added to reaction rate in simulation of dead zone in insensitive high explosives detonation diffraction. The dead zone is captured. The result agrees with reported data qualitatively.

A parallel three-dimensional smoothed particle hydrodynamics cede CSPH3D based on MPI is introduced, including computational scheme, parallel method, data structure, flow chart, and skills in programming. A simulation of three-dimensional ejection and penetration model shows that the CSPH3D code preferably simulates these phenomenon and the parallel efficiency is high. The parallel efficiency of a 1527402 particles micro-ejection model and a 1454225 particles penetration model achieves 80% with 100 processors.