物质点法在爆轰波碰撞炸药猛度试验中应用

使用常规有限元法对炸药爆轰过程进行数值计算时,网格大尺度畸变,往往导致求解精度低,甚至求解失败的情况.采用同时兼顾拉格朗日方法和欧拉方法优点的物质点法(MPM)对爆轰波碰撞炸药猛度试验进行数值计算,得出爆轰波碰撞过程及其在铅柱中传播的应力、应变分布.与试验结果进行对比得知,铅柱压缩误差仅为0.17%.由此验证了物质点法在炸药猛度试验应用中的准确性,为爆轰波传播、碰撞过程分析及炸药猛度测试提供了一种新方法.     

压装TNT二维稳态爆轰反应区的流场分析

用实验测得两种直径的压装TNT药柱二维定常爆轰波的爆速和前导冲击波形状,利用二维C-J条件和守恒关系式求得其压力和反应度沿流线的变化曲线,从而研究压装TNT二维定常爆轰波反应区的流场性态,这对于爆轰理论发展和精密装药设计都是很重要的.     

基于不同状态方程的聚能装药爆轰波形对比及验证

为了对典型JO-9159聚能装药的爆轰波形进行分析,分别通过JWL状态方程、JWL与γ律联合方程表征JO-9159聚能装药的爆轰产物状态,计算得到了基于不同状态方程的爆轰波传播轨迹。采用高速扫描相机获得了该装药在特定位置处的爆轰波形随时间的变化规律,通过仿真与试验结果对比表明:通过JWL与γ律联合方程计算得到的爆轰波传播轨迹与试验结果更为接近,说明该方程可以更准确地描述JO-9159聚能装药的爆轰产物状态。研究结果可为聚能装药爆轰波形分析及破甲战斗部设计提供参考。     

凝聚相炸药爆轰波冲击空穴塌陷过程的高精度数值模拟

利用高精度数值格式对凝聚相炸药爆轰波冲击空穴塌陷过程进行了研究. 采用RGFM和Level-set相结合的物质界面处理方法,克服了爆炸流场中高密度比、高压强比的物质界面容易引起数值振荡这一问题. 采用5阶WENO高精度有限差分格式进行空间离散,采用3阶TVD Runger-Kutta法进行时间离散,开发了高精度动态并行程序,并运用该程序对凝聚相炸药爆轰波冲击空穴后的塌陷过程进行了数值模拟,给出了PBX-9404炸药在不同入射压力下的起爆距离及冲击起爆过程中的压力历史曲线图,并分别与实验结果对比后进行了验证. 对凝聚相炸药中的空穴塌陷过程进行了高精度数值计算,给出了不同时刻的密度、压力和粒子速度云图,研究结果表明由于空穴的存在,爆轰波的传播过程中会在空穴处形成射流现象以及密度和能量的汇聚,进一步产生高温高压的热点.      

炸药殉爆实验和数值模拟

为解决炸药殉爆实验可以给出炸药殉爆条件,但不能得到炸药爆炸过程细节的问题,进行了固黑铝(GHL)炸药殉爆实验,通过观测残留炸药和见证板变形,判断被发炸药爆炸情况. 并采用非线性有限元计算方法对炸药殉爆实验进行了数值模拟计算. 计算模型中主要考虑了主发炸药爆炸冲击波在空气中的传播及其对被发炸药的冲击起爆. 用欧拉法描述主发炸药及周围空气介质,用拉格朗日法描述被发炸药和见证板. 通过数值模拟计算,分析了炸药殉爆过程中,被发炸药爆轰波的成长历程. 结果表明:被发炸药起爆点位于药柱下端,爆轰波先向下传播,使底部炸药先爆炸,然后转为向上传播起爆整个炸药柱;炸药底端压力不高,远低于炸药C-J爆压,对见证板的破坏作用较小.     

FAE街道爆炸的三维数值模拟与面向对象编程

采用气液两相方程,对燃料空气炸药（FAE）气液两相爆轰的爆轰波发展过程进行数值模拟,以此结果作为计算的初值,用三维有限体积TVD方法和面向对象编程技术,模拟FAE爆炸产生的冲击波在街道中的传播规律以及对建筑物的毁伤效应,计算结果表明,三维有限体积TVD方法描述了FAE爆炸冲击波的反射,绕射,传播等物理现象,在计算机模型的构造和管理中,C＋＋和面向对象编程技术体现出极高的优越性,对连续边界的处理,采用了垂直扰动法,得到了清晰的物理图像。     

脉冲爆轰推进装置中波过程的实验研究与分析

研究了脉冲爆轰推进装置(PDE)中各种波的过程与相互作用,包括爆轰波的建立与传播、爆轰波后的Tay-lor稀疏波以及因开口处气流膨胀产生的反射稀疏波等.在实验中引入了预燃腔从而实现了氢气/空气混合气体的快速起爆.针对实验结果,采用CJ理论以及不带化学反应的多组分Euler方程对整个PDE装置的波运动过程进行了分析,其结果显示出口端反射稀疏波足以排出爆轰管内的爆轰产物,使爆轰管内的压力降低到较低的水平.最后给出了一定燃料/氧化剂的混合气体和初始压力下估算该PDE的推力性能以及最大可能的工作频率的方法.     

环形激波和爆轰波会聚过程的气体动力学特性

首先应用修正的激波CCW关系式分析了激波和爆轰波会聚过程及其化学反应的影响, 得到了初始Mach数相同条件下激波聚焦效果比爆轰波强的定性结论. 然后求解Euler方程和基元化学反应模型对环形激波和爆轰波的会聚过程进行了数值模拟, 分析了会聚过程的气体动力学特性、波系结构的差异及其产生原因. 通过对波绕射、会聚和聚焦点附近状态参数变化的对比分析, 揭示了爆轰波化学反应在绕射和会聚不同阶段所起的作用, 验证了理论分析结果.     

方爆轰管内压力波和火焰加速过程的试验

为探究缓燃向爆轰转捩过程中压力波与火焰加速的过程,在60 mm×60 mm×2 000 mm方爆轰管内,用乙炔(C2H2)和空气混合物进行了单次爆轰性能研究.试验在光滑爆轰管和加扰流器的爆轰管内分别进行,在爆轰管上按200 mm的间距同截面布置压力传感器与离子探针,用来监控压力波和火焰各自的信号.根据压力波和火焰位置、速度随时间的变化规律来分析压力波和火焰加速过程.结果表明:在光滑和加入扰流器的爆轰管内火焰的加速传播起主导作用;光滑管内峰值压力波速度大于初始压力波和火焰速度,而加入扰流片后,过爆前初始压力波、峰值压力波和火焰的速度基本上同步发展,过爆后初始压力波速度短暂跃升后最终三者耦合在一起形成爆轰波.     

用修正的KHT状态方程预后炸药爆轰性能

预报炸药的爆轰性能。方法修正部分气态爆轰产物的ＫＨＴ状态方程参数及游离Ｃ的Ｃｏｗａｎ状态方程参数数,计算ＣＨＮＯ炸药的爆轰参数及圆筒试验过程。结果计算结果与实验结果相吻合。     

气相爆轰波贯穿胞格变化研究

根据气相爆轰波胞格结构的规则特则特性建立一种以子单元分解胞格结构的新方法,根据该子单元的性质和斜冲击波关系,求解胞格结构中的三波点对撞问题,并推导贯穿胞格时,爆轰速度的相对波动幅度仅依赖于胞格几何形状的关系,然后,采用爆炸波衰减模型,计算气相爆轰波贯穿胞格时的衰减过程,计算结果与实验符合得较好。     

气相爆轰波与空气冲击波在斜劈上Mach反射的特征参数比较

以传统的烟迹技术记录胞格结构,证实了气相爆轰波 在管道中与斜劈要互作用发生Mach反射,得到了系统的实验结果,在此基础上详细地分析了爆轰波反射后流场的特征参数（三波咪这线与楔面的夹角X与斜劈楔角X与斜劈楔角及气体初压P0的关系等）,并与空气冲击波在斜劈上Mach反射后流场的特征参数进行了比较,初步探讨了参数偏差的原因,这一研究有助于比较深入,全面地了解爆轰波的前导冲击波和横波在Mach反射中所起的作用。     

气相爆轰波在衬多孔钢板管道中的衰减和恢复

对气相爆轰波在衬有多孔钢板管道中的传播进行了实验研究. 结果表明, 多孔钢板对气相爆轰波的传播有衰减作用;同时它所产生的紊流作用又能使已衰减的爆轰波在一定程度上得到恢复.     

近临界马赫数条件下楔形诱导斜爆轰数值模拟

在爆轰推进技术的有关研究中,形成稳定的爆轰波对于爆轰发动机产生稳定推力至关重要.通过理论推导确定了斜爆轰角与来流马赫数间的曲线关系,并给出临界马赫数.进一步采用一步总包反应模型,开展近临界马赫数条件下有、无黏性扩散的斜爆轰大尺度数值模拟.针对临界马赫数近CJ特征的斜爆轰,有、无黏性两种情况下的模拟结果表明,近临界马赫数条件下斜爆轰存在CJ爆轰上行现象,爆轰驻定后,具有全局波动特征.      

基于流固耦合算法的多孔铜爆炸压实数值模拟

本文利用LS-DYNA有限元程序及流固耦合算法,对多孔铜的爆炸压实过程进行数值模拟研究。建立了多孔铜爆炸压实的三维四分之一模型,模型只由空气、炸药和多孔铜三部分组成。炸药、空气与多孔铜的相互作用通过流固耦合算法来实现。根据得到的压力云图,对球面爆轰波和滑移爆轰的传播,多孔铜侧壁压力变化等问题进行了分析。结果表明,爆炸压实过程中多孔铜侧壁压力随滑移爆轰的向下传播而逐渐增加,当夹角增加到90度后,角度及压力不再变化,形成稳定的滑移爆轰,多孔铜侧壁压力达到最大值,约为炸药爆压的85%。     

激光支持等离子体爆轰波流场研究

采用高速纹影摄影法记录激光支持等离子体爆轰波流场的演化,得到流场冲击波阵面形状和尺寸随时间变化规律。建立激光支持等离子体爆轰波流场演化的物理模型,采用计算流体力学方法模拟计算激光击穿空气后,激光作用时间内及激光作用结束后等离子体流场演化过程。结果表明,纹影摄影法能有效观测到激光支持等离子体爆轰波的传播情况。等离子体爆轰波在激光作用初期为椭球形,随着时间的增加逐渐转变为球形。计算的等离子体爆轰波流场密度分布图与实验结果基本吻合。     

塑料导爆管内腔气体介质对传爆可靠性的影响

利用变色塑料导爆管、氮气、氧气和真空系统对管腔内介质影响导爆管传爆可靠性做了研究。氧气有利于传爆，氮气不利于传爆。当导爆管内腔接近真空状态时（－0．0987MPa以下），导爆管仍可传爆，证明导爆管可以作为航天火工品使用。进一步证实导爆管传爆的必要条件是超前冲击波、管道效应以及支持爆轰波波阵面继续传播的足够能量，而导爆和内腔气体介质不是导爆管传爆的必要条件。     

火焰聚心和激波聚焦诱导爆震波

对几种不同结构形式的爆震发生器进行了数值模拟,来研究激波聚焦和火焰聚心现象、气动特性及其机理。数值计算采用多组分理想气体详细的化学反应机理、二维轴对称非定常流动Navier-Stokes方程来模拟流体动力学和化学动力学过程。数值计算发现用较低的点火能量对射流火焰燃烧器中可燃混合物点火,层流火焰在狭窄管壁作用下加速,射流火焰在轴线上汇聚过程,有利于激波的加强。强激波加速火焰,在多重激波与火焰反复作用下,激波和火焰面之间出现热点,热点迅速放大并形成压力很高的过驱爆震波,而后衰减为稳定的爆震波,不同的激波聚焦腔爆震波的形成过程不同。对数值计算的结果分析得出了起爆距离和稳定爆震距离,为进一步试验提供参考。     

初始条件影响气体非稳定爆轰波在弯管中传播特性的实验研究

对丙烷 氧气 空气的预混气体非稳定爆轰波通过圆形截面弯管时传播特性包括弯管角度以及预混气体的初始浓度和初始压力的影响因素进行了初步的实验研究 .结果表明 ,初始条件对可燃气体非稳定爆轰波经过弯管时爆轰参数有很大影响 .随着弯管角度的增加 ,弯管后非稳定爆轰波强度先增加后减小 ,弯管前非稳定爆轰波强度变化很小 ;随着预混气体的初始浓度和初始压力的增加 ,弯管前后非稳定爆轰波强度增加 .这一实验研究结果对安全使用管道阻火器具有重要的实际意义