以DSD理论和LS方法为基础的程序燃烧法

给出了二维弯曲爆轰波后产物流场计算方法。爆轰波阵面传播规律满足Detonation Shock Dynamics (DSD)理论并用level set (LS)方法计算,波阵面传播规律与波后流场的耦合通过程序燃烧法实现,反应进程变量可作为LS函数的函数给出。爆轰波从刚性细管道向粗管道传播产生绕射的二维计算结果表明,化学反应速率不影响波后流场分布,只影响反应区结构。此方法可用于钝感炸药的驱动计算问题。     

凝聚炸药非理想爆轰的数值模拟

本文对凝聚炸药非理想爆轰的数值模拟进行了系统介绍。对敏感、钝感炸药的理想与非理想爆轰的数值模拟方法、宏观反应速率方法及近年发展起来的ＤＳＤ（ＤｅｔｏｎａｔｉｏｎＳｈｏｃｋＤｙｎｔｉｍｉｃｓ）方法和波阵面追踪法,以及对这些方法的适应能力、存在的主要问题进行了分析与讨论。      

驻定在高速弹丸上的斜爆轰波

给出了锥形弹丸以 16 0 0～ 2 10 0m/s的速度射入预混的氢 空气混合介质中形成驻定斜爆轰波的试验结果 ,并讨论锥顶角、飞行体速度以及可燃介质当量比等因素对形成驻定斜爆轰波的影响。     

多孔颗粒火药床燃烧转爆轰的敏感性研究

本文以不同粒度硝胺火药RDU-4,硝胺火药RPB,3/1-S为研究对象,探讨了多孔颗粒火药床燃烧转爆轰的敏感性。研究结果表明,火药颗粒尺寸越小,硝胺火药燃烧转爆轰越易实现:硝基胍组分具有抑制燃烧转冲击的作用,就同一装填条件下燃烧转冲击的难易而言,以RPB最易,3/1-S次之,RDU-4最难     

低温下小尺度钝感炸药非理想爆轰实验研究

钝感炸药爆轰波法向速度Dn与当地曲率密切相关,为研究这种非理想行为,测量了－30 ℃下直径为10、12.5、15、30 mm的钝感炸药的定态爆轰速度和爆轰波形,并与24 ℃下的实验结果进行了比较。研究表明,低温下爆速随药柱直径的增加而增加;直径相同时,温度越高,爆轰波形越平坦;温度相同时,直径越大,爆轰波形越平坦;直径较小时,低温下的爆速小于常温下的爆速,相应的法向爆速与当地曲率的关系曲线Dn()位于常温下曲线的下方;直径较大时,低温下的爆速大于常温下的爆速,相应的Dn()关系曲线位于常温下曲线的上方。     

金属约束下定常非理想爆轰的理论研究

为了对金属约束条件下的定常非理想爆轰进行理论研究,对未反应炸药和爆轰产物采用JWL形式状态方程,对金属采用p(,T)形式状态方程。采用过爆轰前导冲击波的流线偏转角在影响域内沿高度线性变化的假设,并且由未反应炸药和金属的冲击波极曲线的交点确定爆轰冲击边缘角,则可从未反应炸药斜冲击波极曲线关系式求出爆轰前导冲击波阵面的形状。采用爆轰前导冲击波阵面之后的流线是直线的假设,则爆轰流动控制方程由偏微分方程变为沿流线的常微分方程,沿着所有流线求解便给出爆轰化学反应区内声速线与化学反应结束线的位置。理论分析同时给出约束金属折射冲击波后面的流体的流动状态。理论结果给出的爆轰化学反应区结构特征和约束金属内的流动状态特征与高精度数值模拟的结果符合良好,说明本文中给出的理论方法具有良好的合理性和适用性。     

粒状火药床燃烧转爆轰的研究

以多－４５火药为研究对象，采用离子探针测速和应变片测压的技术，对强约束（钢管）条件下火药颗粒床燃烧转爆轰全过程中的火焰烽和压缩波进行了测量。提出了粒状火药燃烧转爆轰的机理。还研究了不同实验条件下火药燃烧转爆轰的敏感性。     

高密度B炸药的燃烧转爆轰实验研究

采用电探针及压力传感器测试技术对密度为1.597 g/cm3的固体B炸药（TNT/RDX=40/60）的燃烧转爆轰性能进行研究。实验结果表明,在较强的约束条件下（45号钢管,内径20 mm,外径64 mm,长500 mm）,B炸药形成了DDT现象,诱导爆轰距离为295~310 mm。     

国外炸药燃烧转爆轰的实验研究

报导了国外炸药燃烧转爆轰（ＤＤＴ）的实验研究情况，主要涉及国外ＤＤＴ研究中采用的实验装置测试方法、铸状炸药的ＤＤＴ机理，多孔装药的ＤＤＴ机理和各种因素对ＤＤＴ过程的影响。     

含能材料密实床燃烧转爆轰的数值模拟

本文建立了含能材料密实床燃烧转爆轰的全粘性欧拉二维非定常两相反应流数学模型。将SIMPLE型的数值计算方法引入燃烧转爆轰的二维数值计算。对跨音速流动的处理和可压缩流体压力校正方程的建立提出了改进方法,并以无起爆药雷管作为算例。结果表明,本方法较好克服了二维两相流数值解的振荡现象。     

流动系统中爆轰波传播特性的数值模拟

基于带化学反应的二维Euler方程,采用氢气/空气的9组分19步基元反应简化模型,对充有当量 比的氢气/空气预混气体的矩形爆轰流场中爆轰波的传播过程进行了数值模拟,讨论了均匀来流对爆轰波传 播的影响。数值结果表明,在均匀来流的影响下,上游方向上的燃烧强度大于C-J爆速,下游方向上的燃烧强 度小于C-J爆速;上游方向传播的爆轰波的阵面压力大于下游方向传播的爆轰波的阵面压力。所以,经典的 C-J爆轰理论并不适用于流动系统中爆轰波传播特性的研究。     

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.The project supported by the National Natural Science Foundation of China (59906005), the Teaching and Research Award Program for Outstanding Young Teachers in High Education Institutions of MOE, China     

混合炸药爆速预报的新方法

为了提高混合炸药爆速尤其是含铝混合炸药爆速的计算精度,基于体积加权法思路构造一种新的混合炸药爆速计算公式,对近50种常见的混合炸药配方进行了计算,并与实验值进行对比,结果很好吻合,最大误差低于3%,平均误差低于1%,与BKW和Urizar常用方法相比,计算精度有显著提高。

     

A double-front structure of detonation wave as the result of phase transitions

Using thermochemical code calculations, we show that the nanographite–nanodiamond phase transition, which may occur in the detonation products of a number of carbon containing explosives, can affect the detonation properties and can cause a specific detonation regime with some unusual peculiarities. Among them, we first note the failure of the Chapman–Jouguet condition and the presence of the sonic plane, where the Mach number is equal to unity, in a detonation product expansion wave at a lower pressure than that at the Chapman–Jouguet point. The peculiarities of this detonation regime are demonstrated by the example of TNT, HNS, and RDX. The computed detonation velocities are in excellent agreement with experiments over a wide range of initial charge densities for all of the investigated explosives. The results of this work allow one to explain, e.g., contradictory experimental data on the detonation pressure and on the length of the reaction zone for TNT. We believe that some other solid–solid, solid–liquid, and liquid–liquid phase transformations in the detonation products may also cause a detonation regime with the same features as shown here for the nanographite–nanodiamond transition. We suggest a computational study that should facilitate proposing detonation experiments strongly arguing in favor of the model presented. PACS 47.40.-x; 47.40.Rs; 64.70.-p; 64.70.Kb; 05.70.-a; 05.70-.CeThis paper was based on the work that was presented at the 19th International Colloquium on the Dynamics of Explosions and Reactive Systems, Hakone, Japan, July 27–August 1, 2003.     

三维气相爆轰动态并行计算程序设计与开发

在三维气相爆轰数值研究中,网格精度和计算域的规模导致网格数占有非常庞大的计算资源,进而给数值模拟带来了极大的挑战。本文针对这一难题,采用5阶WENO格式对带化学反应Euler方程组进行空间离散,基于MPI（MessagePassingInterface）并行模式开发了高精度动态并行代码,并对爆轰波在带有障碍物的三维方形管道中的传播过程进行计算。计算结果表明,高精度动态并行计算能够很好的模拟三维气相爆轰波在大尺寸管道中的传播,不仅提高了计算效率,而且提高了爆轰波阵面的分辨率。与高精度静态并行相比,高精度动态并行计算减少了界面数据通信时间,从而进一步提高了计算效率。因此,高精度动态并行程序为探究三维气相爆轰新的物理机制提供有效的手段。     

驻定斜爆轰波并行数值模拟

采用多组分化学反应Euler方程组对驻定在高速飞行弹丸上的斜爆轰波流场进行了数值模拟。计算中分别采用TVD格式和基元反应模型,并基于并行编程模型MPI（message passing interface）实现了非结构网格上的并行计算,对流项和化学反应项用时间分裂法进行处理。计算结果表明并行计算能有效地提高计算速度,扩展计算规模,为进一步研究超驱爆轰推进技术奠定基础。     

爆轰法合成碳包覆纳米铜颗粒

以硝酸铜和柠檬酸制成的干凝胶为主要反应物,加入油酸有机碳源和黑索今炸药,在氮气保护气氛下在爆炸容器中引爆,成功地合成了碳包覆纳米铜颗粒。分别采用X射线衍射、透射电镜对产物形貌特征进行表征。结果表明,在爆轰产物中富含碳包覆纳米铜颗粒,产物呈圆球体,具有完好的核壳结构形貌,颗粒粒径在10~40nm 之间,外层碳壳结构主要由无定型碳和石墨构成。并对爆轰法合成碳包覆纳米铜颗粒的形成机理进行了分析。     

Diffraction of a planar strong detonation wave by a moving three-dimensional thin body

An analytic solution is obtained for the diffraction of a planar strong detonation wave by a three-dimensional thin body moving in the opposite direction. The planform and the thickness distribution of the body can be arbitrary and the speed of the body can be either supersonic or subsonic relative to the undisturbed stream ahead of the wave or to that behind the wave. The solution is a generalization of the previous solution of Ting and Gunzburger for the shock diffraction.     

松弛方法在计算凝聚炸药爆轰问题中的应用

利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。     

An approximate analytical solution of the cylindrical detonation wave generated by the linear explosion

The solution of the cylindrical detonation wave generated by the linear explosion was obtained by numerical method in ref.[1].In this paper,when the ratio of specific heatγ》1 by using the enlargement coordinate method,the first-order analytical solutions are obtained.The perturbation parameter isε=1/γ~2 .The correction of these solutions is checked at the end of this paper.