To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitro...To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.展开更多
A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simu...A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simulation models for high-speed piston impacts on explosive supercells were established.Simulations were also performed to investigate shock-induced reactions of various high-energy explosives.The fraction of reacted explosive molecules in an initial supercell changed linearly with the propagation distance of the shock-wave front.The corresponding slope could be used as a reaction rate for a specific shock-loading velocity.Reaction rates that varied with the shock-loading pressure exhibited two-stage linearities with different slopes.The two inflection points corresponded to the initial and accelerated reactions,which respectively correlated to the thresholds of shock-induced ignition and detonation.Therefore,the ignition and detonation critical pressures could be determined.The sensitivity could then be a quantitative prediction of the critical pressure.The accuracies of the quantitative shock sensitivity predictions were verified by comparing the impact and shock sensitivities of common explosives and the characteristics of anisotropic shock-induced reactions.Molecular dynamics simulations quantitatively predict and rank shock sensitivities by using only crystal structures of the explosives.Overall,this method will enable the design and safe use of explosives.展开更多
Chemical inclusions significantly alter shock responses of crystalline explosives in macroscale gap experiments but their microscale dynamics origin remains unclear.Herein shock-induced energy localization,overall phy...Chemical inclusions significantly alter shock responses of crystalline explosives in macroscale gap experiments but their microscale dynamics origin remains unclear.Herein shock-induced energy localization,overall physical responses,and reactions in a-1,3,5-trinitro-1,3,5-triazinane(a-RDX)crystal entrained various chemical inclusions were investigated by the multi-scale shock technique implemented in the reactive molecular dynamics method.Results indicated that energy localization and shock reaction were affected by the intrinsic factors within chemical inclusions,i.e.,phase states,chemical compositions,and concentrations.The atomic origin of chemical-inclusions effects on energy localization is dependent on the dynamics mechanism of interfacial molecules with free space volume,which includes homogeneous intermolecular compression,interfacial impact and shear,and void collapse and jet.As introducing various chemical inclusions,the initiation of those dynamics mechanisms triggers diverse decay rates of bulk RDX molecules and hereby impacts on growth speeds of final reactions.Adding chemical inclusions can reduce the effectiveness of the void during the shock impacting.Under the shockwave velocity of 9 km/s,the parent RDX decay rate in RDX entrained amorphous carbon decreases the most and is about one fourth of that in RDX with a vacuum void,and solid HMX and TATB inclusions are more reactive than amorphous carbon but less reactive than dry air or acetone inclusions.The lessdense shocking system denotes the greater increases in local temperature and stress,the faster energy liberation,and the earlier final reaction into equilibrium,revealing more pronounced responses to the present intense shockwave.The quantitative models associated with the relative system density(RD_(sys))were proposed for indicating energy-localization mechanisms and evaluating initiation safety in the shocked crystalline explosive.RD_(sys)is defined by the density ratio of defective RDX to perfect crystal after dynamics relaxation and reveals the global density characteristic in shocked systems filled with chemical inclusions.When RD_(sys)is below 0.9,local hydrodynamic jet initiated by void collapse dominates upon energy localization instead of interfacial impact.This study sheds light on novel insights for understanding the shock chemistry and physical-based atomic origin in crystalline explosives considering chemical-inclusions effects.展开更多
基金supported by the National Natural Science Foundation of China (22275018)the Project of State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology)(Grant No.QNKT20-04)。
文摘To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.
基金supported by the National Natural Science Foundation of China(Grant No.11832006).
文摘A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simulation models for high-speed piston impacts on explosive supercells were established.Simulations were also performed to investigate shock-induced reactions of various high-energy explosives.The fraction of reacted explosive molecules in an initial supercell changed linearly with the propagation distance of the shock-wave front.The corresponding slope could be used as a reaction rate for a specific shock-loading velocity.Reaction rates that varied with the shock-loading pressure exhibited two-stage linearities with different slopes.The two inflection points corresponded to the initial and accelerated reactions,which respectively correlated to the thresholds of shock-induced ignition and detonation.Therefore,the ignition and detonation critical pressures could be determined.The sensitivity could then be a quantitative prediction of the critical pressure.The accuracies of the quantitative shock sensitivity predictions were verified by comparing the impact and shock sensitivities of common explosives and the characteristics of anisotropic shock-induced reactions.Molecular dynamics simulations quantitatively predict and rank shock sensitivities by using only crystal structures of the explosives.Overall,this method will enable the design and safe use of explosives.
基金the financial support from National Natural Science Foundation of China(Grant Nos.11872119,12172051,and 11972329)Natural Science Foundation of Hubei Province(Grant No.2021CFB120)。
文摘Chemical inclusions significantly alter shock responses of crystalline explosives in macroscale gap experiments but their microscale dynamics origin remains unclear.Herein shock-induced energy localization,overall physical responses,and reactions in a-1,3,5-trinitro-1,3,5-triazinane(a-RDX)crystal entrained various chemical inclusions were investigated by the multi-scale shock technique implemented in the reactive molecular dynamics method.Results indicated that energy localization and shock reaction were affected by the intrinsic factors within chemical inclusions,i.e.,phase states,chemical compositions,and concentrations.The atomic origin of chemical-inclusions effects on energy localization is dependent on the dynamics mechanism of interfacial molecules with free space volume,which includes homogeneous intermolecular compression,interfacial impact and shear,and void collapse and jet.As introducing various chemical inclusions,the initiation of those dynamics mechanisms triggers diverse decay rates of bulk RDX molecules and hereby impacts on growth speeds of final reactions.Adding chemical inclusions can reduce the effectiveness of the void during the shock impacting.Under the shockwave velocity of 9 km/s,the parent RDX decay rate in RDX entrained amorphous carbon decreases the most and is about one fourth of that in RDX with a vacuum void,and solid HMX and TATB inclusions are more reactive than amorphous carbon but less reactive than dry air or acetone inclusions.The lessdense shocking system denotes the greater increases in local temperature and stress,the faster energy liberation,and the earlier final reaction into equilibrium,revealing more pronounced responses to the present intense shockwave.The quantitative models associated with the relative system density(RD_(sys))were proposed for indicating energy-localization mechanisms and evaluating initiation safety in the shocked crystalline explosive.RD_(sys)is defined by the density ratio of defective RDX to perfect crystal after dynamics relaxation and reveals the global density characteristic in shocked systems filled with chemical inclusions.When RD_(sys)is below 0.9,local hydrodynamic jet initiated by void collapse dominates upon energy localization instead of interfacial impact.This study sheds light on novel insights for understanding the shock chemistry and physical-based atomic origin in crystalline explosives considering chemical-inclusions effects.
