固体炸药爆轰的一种考虑热学非平衡的反应流动模型

基于固体炸药爆轰过程中化学反应混合区内的固相反应物与气相生成物处于力学平衡状态及热学非平衡状态的事实,提出一种考虑热学非平衡效应的反应流动模型来描述固体炸药的爆轰流动现象.该爆轰流动模型的主要特点是,在反应混合物Euler方程和固相反应物质量守恒方程的基础上,通过附加一套关于固相反应物的组分物理量的流动控制方程来表达固相反应物与气相生成物之间的热学非平衡效应.根据反应混合区内固相反应物与气相生成物这两种化学组分保持各自内能守恒的混合规则,并借助它们具有压力相等的性质以及满足体积分数总和为1的条件,推导获得的附加方程有：固相反应物的内能演化方程、体积分数演化方程及反应混合物的压力演化方程.这样,建立的爆轰模型包括：反应混合物的质量守恒方程、动量守恒方程、总能量守恒方程、压力演化方程,以及固相反应物的质量守恒方程、内能演化方程、体积分数演化方程.对所获得的爆轰模型方程组采用一个时空二阶精度的有限体积法进行数值求解,典型爆轰问题算例结果表明本文提出的固体炸药爆轰模型是合理的.     

Studies on aluminum powder combustion in detonation environment

The combustion mechanism of aluminum particles in a detonation environment characterized by high temperature (in unit 10³ K), high pressure (in unit GPa), and high-speed motion (in units km/s) was studied, and a combustion model of the aluminum particles in detonation environment was established. Based on this model, a combustion control equation for aluminum particles in detonation environment was obtained. It can be seen from the control equation that the burning time of aluminum particle is mainly affected by the particle size, system temperature, and diffusion coefficient. The calculation result shows that a higher system temperature, larger diffusion coefficient, and smaller particle size lead to a faster burn rate and shorter burning time for aluminum particles. After considering the particle size distribution characteristics of aluminum powder, the application of the combustion control equation was extended from single aluminum particles to nonuniform aluminum powder, and the calculated time corresponding to the peak burn rate of aluminum powder was in good agreement with the experimental electrical conductivity results. This equation can quantitatively describe the combustion behavior of aluminum powder in different detonation environments and provides technical means for quantitative calculation of the aluminum powder combustion process in detonation environment.     

H32.685N16.8430O3.0和H8.154N4.5770O3.0炸药爆轰参数的数值模拟

用Gibbs自由能最小原理,通过解化学平衡方程组,求解HNO型炸药H32·685N16·8430O3·0(70/30Hydrazine/HydrazineNitrate)和H8·154N4·5770O3·0(21/79Hydrazine/HydrazineNitrate)爆轰产物系统的平衡组分,计算结果与用BKW和LJD方法计算的结果相近。以Ree修正的WCA状态方程并考虑Ross软球修正的硬球微扰理论作为爆轰气相产物的物态方程,用自编的程序对H32·685N16·8430O3·0和H8·154N4·5770O3·0炸药爆轰参数作了预言,爆轰CJ点的爆速、爆压和爆温的计算结果与实验值吻合得很好。     

炸药爆轰的连续介质本构模型和数值计算方法

假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性. 关键词： 炸药爆轰 本构模型 化学反应率方程 数值模拟     

Numerical studies on autoignition and detonation development from a hot spot in hydrogen/air mixtures

Detonation development inside spark ignition engines can result in the so called super-knock with extremely high pressure oscillation above 200?atm. In this study, numerical simulations of autoignitive reaction front propagation in hydrogen/air mixtures are conducted and the detonation development regime is investigated. A hot spot with linear temperature distribution is used to induce autoignitive reaction front propagation. With the change of temperature gradient or hot spot size, three typical autoignition reaction front modes are identified: supersonic reaction front; detonation development and subsonic reaction front. The effects of initial pressure, initial temperature, fuel type and equivalence ratio on detonation development regime are examined. It is found that the detonation development regime strongly depends on mixture composition (fuel and equivalence ratio) and thermal conditions (initial pressure and temperature). Therefore, to achieve the quantitative prediction of super-knock in engines, we need use the detonation development regime for specific fuel at specific initial temperature, initial pressure, and equivalence ratio.     

高能炸药爆速曲率关系的一种计算方法

研究平面爆轰波爆轰化学能和状态方程之间的关系,得到平面爆轰波关系式和CJ条件.基于柱坐标系下的一维流体力学方程组,采用坐标变换方法得到含曲率的爆轰波关系式及广义声速条件.在反应区内采用压力平衡和温度平衡的混合规则,在给定的复杂状态方程和复杂反应率下求解,得到PBX9404和PBX9502炸药的爆速和曲率的关系.已知炸药...     

Elastic properties of layered crystals

The characteristic features of the elastic properties of layered crystals and their dependence on temperature and pressure are analyzed. The relations between the elastic constants of hexagonal layered crystals are given. It is shown that the anomalous behavior of the elastic constants in the temperature region of a phase transition affects both the magnitude and sign of the thermal expansion coefficients of layered crystals. From analyzing the pressure and temperature dependences of the elastic constants, it is found that the anharmonicity of the bonding forces between the layers is much greater than the anharmonicity of the intralayer forces. The contribution from thermal expansion to the variations of the elastic constants with temperature is estimated.     

炸药爆轰产物Jones-Wilkins-Lee状态方程不确定参数

Jones-Wilkins-Lee (JWL)状态方程是一种不显含化学反应、由实验方法确定参数的半经验状态方程, 能比较精确地描述爆轰产物的膨胀驱动做功过程. 在JWL状态方程中有多个未知(不确定)参数需要确定. 传统的确定JWL状态方程参数的方法是“调参数”, 人为因素影响较大, 无法给出参数的不确定性信息. 本文利用贝叶斯分析方法研究了炸药的不确定参数, 该方法能够基于以往的认识、实验和模拟数据标定(calibration)不确定参数. 在本文结果中, 不确定参数的后验分布均值与文献结果相符合, 基于参数标定结果的数值模拟90%置信区间完全包含实验数据. 数值标定结果说明贝叶斯参数标定适用于确定样品炸药的JWL状态方程参数. 特别是, 在本文JWL状态方程参数标定过程中极大地减少了人为因素的影响.     

爆轰产物的冷比内能与冷压的函数表达式

 猛炸药爆轰产物的状态可以用两相的强排斥-平动物态方程（简称为两相的排平物态方程）很好地描述。以爆轰产物分两段的等熵曲线为参考曲线的两相的排平（k, γ）物态方程,已经用于爆轰参数和强爆轰参数的理论估算,所得理论值与实验值符合得很好。为了更方便地估算爆温,有必要给出描述分子间相互作用的比内能项与压力项（分别简称为冷比内能与冷压）。参照描述分子间相互作用的Morse势和Mie势的排斥项,给出了带待定参数A、m、n和l的冷比能项和冷压项,这样的物态方程被称为两相的排平（A, m, n, l）物态方程。用TNT的{D, ρ₀}实验数据组,确定了两相的排平（A, m, n, l）物态方程的参数n=1和l=1/3,因此,可将其简称为两相的排平（A, m）物态方程。它适用于所有的猛炸药的爆轰产物。用硝基甲烷的强爆轰参数{p, D, T}实验数据组对其所做的检验表明,两相的排平（A, m）物态方程是恰当的爆轰产物物态方程。     

A first-principles study on the structural and elastic properties and the equation of state of wurtzite-type cadmium selenide under pressure

A first-principles investigation on the crystal structural and elastic properties and the equation of state of wurtzite-type cadmium selenide (w-CdSe) has been conducted using the plane-wave pseudo-potential density functional theory and the quasi-harmonic Debye model. The elastic constants, the aggregate elastic moduli, the elastic anisotropy, and Poisson's ratio under pressure have been investigated. Our calculated equilibrium lattice constants, the elastic constants, and the aggregate elastic moduli at zero pressure are in good agreement with the experimental data and other theoretical results. The variations in the compressional and shear elastic wave velocities with pressure at zero temperature up to pressure 2.7 GPa have been studied; the computed Debye temperature at zero pressure and zero temperature is in reasonable agreement with the result of Bonello et al., In addition, the equation of state of w-CdSe in the pressure range of 0–2.7 GPa and up to a temperature of 900 K has also been obtained.     

Equation of State of Protoneutron Star Matter

We have calculated the equation of state of protoneutron star matter by using the lowest order constrained variational method. In our calculations, the modern Argonne potential (AV₁₈) together with its older model potential (AV₁₄) are used. It is found that the equation of state for high lepton fraction is stiffer than for low lepton fraction. It is seen that the increasing effect of pressure due to high lepton fraction and due to entropy are comparable. It is shown that the temperature and adiabatic index depend on the values of both entropy and lepton fraction.     

First-principles calculations of elastic,phonon and thermodynamic properties of W

We investigate the elastic properties, lattice dynamical, thermal equation of state and thermodynamic properties of bcc phase W under high pressure using density functional theory. The calculated high-pressure elastic constants of bcc phase W agree well with experimental and theoretical data. Under compression, the phonon dispersion curves of bcc phase W do not show any anomaly or instability. Our calculated zero-pressure phonon dispersion curves are in excellent agreement with experiments. Within the quasiharmonic approximation, we predict the thermal equation of state and other properties including the thermal expansion coefficient, adiabatic bulk modulus, specific heat at constant volume and entropy.     

硅酸镁钙钛矿弹性及热力学特性的第一性原理计算

基于密度泛函理论的第一性原理计算,结合准谐德拜模型研究了高压下硅酸镁钙钛矿的弹性及热力学特性. 计算得到的物态方程数据、热容、热膨胀系数等在宽广的温度和压力范围与实验结果及其他理论计算结果吻合. 根据有限应变理论计算了硅酸镁钙钛矿的弹性常数,并讨论了杨氏模量、泊松比、德拜温度、晶体各向异性随压力的变化.     

Thermal equation of state, and melting and thermoelastic properties of bcc tantalum from molecular dynamics

We performed molecular dynamics simulations with the extended Finnis-Sinclair (EFS) potential to investigate thermal equation of state (EOS), and melting and thermoelastic properties of tantalum. The agreement of the obtained thermal EOS with experiments at ambient conditions is reasonably good. The EFS potential with the two-phase method also reproduced very satisfyingly the high-pressure melting curve, excellently consistent with both the experiments of melting temperature at ambient pressure and shock melting at high pressure. From molecular dynamics simulations, we also obtained the thermoelastic properties of Ta for temperatures up to 3000 K at ambient pressure. Fully including anharmonic effects in molecular dynamics, our calculated elastic constants are in excellent agreement with experimental data. Shear modulus G decreases quickly with increasing temperature.     

Heat conduction in a one-dimensional gas of elastically colliding particles of unequal masses

We study the nonequilibrium state of heat conduction in a one-dimensional system of hard point particles of unequal masses interacting through elastic collisions. A BBGKY-type formulation is presented and some exact results are obtained from it. Extensive numerical simulations for the two-mass problem indicate that, even for arbitrarily small mass differences, a nontrivial steady state is obtained. This state exhibits local thermal equilibrium and has a temperature profile as predicted by kinetic theory. The temperature jumps typically seen in such studies are shown to be finite-size effects. The thermal conductivity appears to have a very slow divergence with system size, different from that seen in most other systems.     

一种简化维里型状态方程预测高温甲烷PVT关系

为满足描述爆轰环境下高温气体高温、中高压状态的需求,本文提出了一种基于Lennard-Jones(LJ)势能函数的简化维里型状态方程Han-Long(HL).应用HL状态方程计算了甲烷1000 K以上112组理论和实验数据,计算所得体积平均绝对偏差约为1%,最大误差为3.28%,远低于DMW状态方程和BS状态方程的计算偏差.采用HL状态方程计算了甲烷冲击试验的热力学数据,计算所得体积偏差均小于3%.结果表明,HL状态方程能够很好的描述高温甲烷的热力学状态.     

A particle level-set based Eulerian method for multi-material detonation simulation of high explosive and metal confinements

The multi-material numerical simulation for energetic system that consists of a high explosive charge and an inert confinement is carried out with an accurate and state-of-the-art Eulerian method. An explosively driven copper tube results in a state of extreme temperature and pressure, coupled to a high speed structural response of metal due to a detonating high explosive (HE). We use the experimentally tuned Ignition and Growth (or JWL++) rate equation for the HE while the elasto-plastic response of inert is modeled by the Mie–Gruneisen equation of state (EOS) and the Johnson–Cook strength model. A new particle level-set based reactive Ghost Fluid Method (GFM) that imposes exact boundary conditions at the material’s interface according to physical restraints is developed to simulate the multi-material detonation problem. Our calculations reproduce the experimental data of both unconfined and confined rate stick problems, suggesting that the method is suitable for detonation simulation of energetic systems.     

从头算方法研究面心立方铝在高温高压下的热力学状态方程

用从头算方法优化计算了面心立方铝的电子结构和总能,得到了它在零温下的状态方程和弹性性质.将得到的总能与晶格体积拟合到Debye模型,获得了非平衡态下的Gibbs自由能与温度、压力之间的关系,在此基础上计算了相应的热状态方程,利用Burakovsky-Preston-Silbar （BPS） 熔化模型计算了铝的熔化曲线.所有的电子结构和总能计算都是基于局域密度近似（LDA）和广义梯度近似（GGA）的平均得到的.计算得到的铝在高温、高压下的状态方程与一些热力学性质和熔化曲线同冲击波和静高压实验数据在225 G 关键词： 铝 热力学状态方程 从头算 熔化曲线     

Ab initio study of lattice dynamics and thermal equation of state of Ni

Lattice dynamics and thermal equation of state of fcc nickel have been studied in the framework of density functional perturbation theory. The influence of the GGA+U on the structure is considered. The calculated phonon dispersion curve accords excellently with the experimental data. Within the quasi-harmonic approximation, the thermal equation of state, thermal expansion coefficient, thermal pressure, bulk moduli and Debye temperature are well reproduced. The thermal properties confirm the available experimental data and are extended to a wider pressure and temperature range.     

均匀炸药冲击起爆和起爆后的行为

 介绍并分析了Campbell等人研究均匀炸药冲击起爆和起爆后行为所获得的实验结果,但不涉及其冲击起爆条件。Campbell等人的实验表明,足够强的冲击波进入硝基甲烷后,经过若干微秒的感应时间,爆轰发生在隔板与炸药间的界面处。这就是说,在均匀炸药中,足够强的冲击虽非瞬时但直接（指不经过其它过程,如爆燃）引发了爆轰。重新处理后的实验数据表明：硝基甲烷起爆后,爆轰波的净爆速小于正常爆速;当进入硝基甲烷的初始冲击波的有效压力p_eff由8.82 GPa升至12.14 GPa时,感应时间t_ind的实验值由3.06 μs降至0.705 μs。以两相的排平（A,m）物态方程描述爆轰产物,较为严格地重新推导了基于热起爆理论的估算感应时间t_ind的公式。在上述p_eff的变化范围内,t_ind的理论值则由248 μs降至0.99 μs,明显地高于实验值。这表明,热起爆理论不适于描述硝基甲烷的冲击起爆行为。从本质上讲,热起爆理论对均匀炸药的冲击起爆行为的描述,不符合物质运动的微观图像,因此,它不适于描述均匀炸药的上述行为。