CE/SE方法数值模拟炸药粉尘爆轰

采用CE/SE方法数值模拟悬浮在空气中的RDX炸药粉尘的两相爆轰过程.炸药颗粒在爆轰波阵面后的高温高速气流中加速并升温,释放能量支持爆轰波传播.数值模拟爆轰波管中的粉尘爆轰,得到爆轰波流场中的物理量分布,确定爆轰参数,数值结果与文献符合较好.数值模拟复杂通道中的炸药粉尘爆轰,预测了爆轰波的发展和传播过程以及爆轰波后的流场演化.数值结果表明CE/SE方法能成功模拟气体-固体两相爆轰,为粉尘爆轰的研究提供了新的数值预测手段.     

拐角空间中悬浮铝粉尘爆轰后效的三维数值模拟

建立三维的铝粉-空气两相爆轰计算模型,采用时-空守恒元解元（CE/SE）方法求解,并开发了悬浮铝粉尘爆轰的三维数值模拟程序.基于消息传递接口（MPI）技术实现了程序的并行化设计.通过对激波管问题以及爆轰管中铝粉-空气两相爆轰实验的模拟验证程序的可靠性.对拐角空间中左侧浓度为368 g·m^-3的铝粉-空气混合物两相爆轰及其在拐角空间右侧和下方空气域内形成的冲击波和温压效应开展数值模拟,获得复杂空间内爆轰波或冲击波的传播、反射以及绕射过程.结果表明：两相爆轰在离铝粉尘区域2 m远的空气域内产生的后效冲击波能达到2.66 MPa的固壁反射压力,火球燃烧范围会超出初始铝粉尘区域约0.8 m,并且造成初始铝粉尘区域附近1.5 m范围内空气的温度高达1 600 K.模拟程序可用于铝粉尘爆轰的后效研究,对工业安全及其防护具有指导意义.     

真实比热模型中铝粉尘两相爆轰波的数值研究

采用多流体模型对铝粉尘两相爆轰波进行数值模拟,研究颗粒能量计算方法对起爆和传播过程的影响.以前的固相颗粒能量的计算一般采用固定比热方法,本文采用随温度变化的真实比热.由于铝颗粒及其产物氧化铝的比热变化很大,模拟得到的爆轰波的速度、压力和波后参数变化和采用固定比热存在较大的差异.变比热计算得到的爆轰波压力、传播速度和实验结果更加接近,而固定比热的计算方法会对这些参数造成高估.对爆轰波的形成进行研究,发现起爆距离主要受起爆能量影响,但是相对于固定比热模型,采用变比热模型得到的起爆距离较短.     

RDX炸药爆轰产物物态方程

采用基于统计物理的爆轰产物物态方程模型,研究了不同初始密度下RDX炸药爆轰产物性质,发现炸药初始密度小于1.20g/cm3时,炸药爆轰后只有气相产物生成,随着炸药初始密度进一步增大,爆轰后有固相碳析出。炸药初始密度小于1.60g/cm3时,各气相产物按均匀混合处理,计算的爆速、爆压和实验值符合较好。随着炸药初始密度进一步增大,计算的爆速、爆压与实验值出现较大的偏差,经分析发现此处可能出现了超临界流体相分离现象,在物态方程模型中考虑了这种效应后,得到较好的结果。     

HMX炸药燃烧转爆轰数值模拟

以两相流模型为基础,气相产物状态方程采用基于统计物理的类CHEQ的计算结果,建立了HMX炸药的燃烧转爆轰数学模型。采用CE/SE方法模拟了颗粒度为125μm的HMX炸药的燃烧转爆轰过程,得到了爆轰参数及流场变化规律。模拟了装填密度对HMX炸药燃烧转爆轰的影响,并与实验进行了对比。数值模拟结果表明,在相同的点火条件下,爆轰成长距离在一定范围内随装填密度呈"U"形变化。     

气液爆轰波及冲击波传播的数值研究

采用高分辨率ＴＶＤ格式、ＭａｃＣｏｒｍａｃｋ格式和二维轴对称气液两相方程组数值模拟了气液爆轰的发展过程及冲击波误减规律,计算结果与国内外研究结果符合良好。     

弱点火条件下炸药燃烧转爆轰的数值模拟

考察颗粒炸药从传导燃烧到对流燃烧再到爆轰的过程.对装填密度为85%的HMX颗粒炸药的燃烧转爆轰过程进行数值模拟,分析传导燃烧、对流燃烧和爆轰的发展过程.点火早期燃烧速度很低,火焰面在8.16 ms之内只前进了不到0.2 mm;形成对流燃烧之后燃烧速度快速增加,只用了0.1 ms就形成了速度为8 165 m·s^-1的稳定爆轰.当炸药颗粒直径或点火压力减小时,形成稳定爆轰所需的时间增加.     

RDX激光点火的一维气相模型

为了研究含能材料RDX在激光作用下的点火特性,在Liau等人的工作基础上,考虑了物理参数随温度和组分浓度的变化,建立了RDX激光点火的气相模型.利用有限差分方法,计算了激光功率密度为400 W/cm2、环境压力为0.1 MPa时RDX的点火过程.结果表明:约90%的RDX以蒸发的形式进入气相,另外的10%在液相中发生分...     

一种以RDX为主炸药的低密度炸药

低密度炸药的特点就是密度低(一般都在0．08～0．8g／cm^3之间)，其爆速、爆压及威力均较低，有粉状、塑性、挠性和硬质等四种物理状态的制品。民用上广泛用于采矿及金属爆炸加工，军事上主要用于子母弹中心管抛撒药和其他抛撒器装药，也可用于装填普通水雷、浮雷和鱼雷等。研制这种炸药，一般可采取两种技术途径：(1)低密度泡沫炸药；(2)低密度塑料黏结炸药。低密度塑料黏结炸药具有成型工艺简单、操作安全，能在低比压下模压成型等优点，近几年已经做了大量的工作，已完成了以PETN为主炸药的JI-9016低密度炸药的研制，其密度为0．6g／cm^3，爆速为2390m/s，爆压为1．427MPa。     

氢氧爆轰波在激波管中的成长机制研究

 针对气相爆轰波成长机制研究,采用压力传感器和高速摄影技术,测试了氢氧混合气体在点火后的火焰波、前驱冲击波以及爆轰波的成长变化过程,计算了冲击波过程参数和气体状态参数,分析了火焰加速机制。实验结果表明,APX-RS型高速摄影系统可用于拍摄气相爆轰波的成长历程;氢氧爆轰波的产生是由于湍流火焰和冲击波的相互正反馈作用,导致反应区内多处发生局部爆炸,爆炸波与冲击波相互耦合,最终成长为定常爆轰波。     

矩形管内临界爆轰动力学数值分析

 对矩形管内临界爆轰动力学特征进行了数值分析。采用基元反应描述爆轰化学反应过程,采用二阶附加半隐的龙格-库塔法和5阶WENO格式求解二维反应欧拉方程。对于25%氩稀释化学计量比的氢氧预混气体,当管道宽度为30 mm、初温为300 K时,产生临界爆轰的预混气体初压为3.5 kPa。在此临界条件下,获得了临界爆轰胞格结构、沿壁面的速度和峰值压力曲线及流场波系演变特征。着重对比分析了矩形管内临界爆轰与普通爆轰在爆轰波速度、平均速度、胞格宽长比、横波结构、未反应气囊及旋涡结构之间的差异,深入认识了临界爆轰的不稳定性和化学反应动力学特征。     

气相爆轰波在分叉管中传播现象的数值研究

数值研究气相爆轰波在分叉管中的传播现象.用二阶附加半隐龙格-库塔法和5阶WENO格式求解二维欧拉方程,用基元反应描述爆轰化学反应过程,得到了密度、压力、温度、典型组元质量分数场及数值胞格结构和爆轰波平均速度.结果表明:气相爆轰波在分叉管中传播,分叉口左尖点的稀疏波导致诱导激波后压力、温度急剧下降,诱导激波和化学反应区分离,爆轰波衰减为爆燃波(即爆轰熄灭).分离后的诱导激波在垂直支管右壁面反射,并导致二次起爆.畸变的诱导激波在水平和垂直支管中均发生马赫反射.分叉口上游均匀胞格区和分叉口附近大胞格区的边界不是直线,其起点通常位于分叉口左尖点上游或恰在左尖点.水平支管中马赫反射三波点迹线始于右尖点下游.分叉口左尖点附近的流场中出现了复杂的旋涡结构、未反应区及激波与旋涡作用.旋涡加速了未反应区的化学反应速率.反射激波与旋涡作用并使旋涡破碎.反射激波与未反应区作用,加速其反应消耗,并形成一个内嵌的射流.数值计算得到的波系演变和胞格结构与实验定性一致.     

一维爆轰波不稳定性的数值模拟

 对当量比为1的乙烷与空气的混合气体的一维爆轰不稳定问题进行数值模拟，得到了不同大小的网格对爆轰不稳定问题数值模拟结果的影响。网格大小Δx由从ZND模型分析得到的导引长度L_in确定，网格大小从0.2L_in变到0.002 5L_in。随着网格变细，没有得到振幅趋于一致的解，每一种网格尺寸得到的解的振幅都互不相同。当网格大小为Δx=0.01L_in、0.005L_in时，得到有规则的爆轰激波阵面压力的振荡，振荡的模式是峰值一大一小的振荡。网格更细时，爆轰波的振荡在计算范围内由一些有规则的振荡和一些较不规则的振荡组成。但爆轰激波阵面压力振荡的波长最后趋于一致，为91~93 mm，与实验得到的胞格长度88 mm很接近。     

大型水平管道中环氧丙烷-铝粉-空气混和物爆燃转爆轰的实验研究

 在长为32.4 m、内径为0.199 m的大型长直水平管道中,对环氧丙烷-空气两相流云雾及环氧丙烷-铝粉-空气三相流云雾的爆燃转爆轰（DDT）过程进行了实验研究。对弱点火条件下多相混和物DDT过程的不同阶段特征进行了分析,对比研究了不同浓度时混和物的燃爆情况。结果表明：浓度为513 g/m³的环氧丙烷-空气混和物及浓度为237和643 g/m³的环氧丙烷-铝粉-空气混和物均能在管道中完成爆燃向爆轰的转变,进入自持爆轰阶段,其胞格尺寸分别为0.28和0.50 m。     

隔板深度对激波聚焦起爆影响的数值模拟研究

考虑几何结构参数对激波聚焦触发爆轰波的复杂影响,对H2/Air预混气的环形射流激波聚焦起爆现象开展了数值模拟研究,详细分析了不同隔板深度条件下的激波聚焦过程、流场演化特征以及爆轰波参数变化规律.研究结果表明,凹腔内激波聚焦诱导的局部爆炸以及隔板前缘处射流形成"卷吸涡"是引起爆轰波触发的两个重要机制,而隔板深度是影响环形...     

地下爆炸空腔压力和温度历程数值模拟

将气体的渗透、气体与空腔壁的辐射换热、气体与周围介质的热传导三种因素相结合,建立地下爆炸空腔稳定后气体温度和压力变化的数学模型,编制一维球对称有限体积数值模拟程序.用该程序对美国内华达试验场冲积土介质中一次地下核爆炸空腔温度和压力历程进行数值模拟,并对压力监测结果进行分析.