【摘 要】
：
DSD (Detonation Shock Dynamics)算法是钝感炸药爆轰过程数值模拟的重要手段,与基于唯象化学反应率的直接数值模拟方法相比,DSD算法具有计算量小的优势。本文介绍了在二维、三维DSD算法研究与程序研制方面所开展的工作与进展情况,包括DSD理论与数值方法的简要概述,以及通过若干典型算例对算法与程序的考核与验证。
【机 构】
：
北京应用物理与计算数学研究所,北京100094 北京应用物理与计算数学研究所,北京100094;北
【出 处】
：
2018第十二届全国爆炸力学学术会议
论文部分内容阅读
DSD (Detonation Shock Dynamics)算法是钝感炸药爆轰过程数值模拟的重要手段,与基于唯象化学反应率的直接数值模拟方法相比,DSD算法具有计算量小的优势。本文介绍了在二维、三维DSD算法研究与程序研制方面所开展的工作与进展情况,包括DSD理论与数值方法的简要概述,以及通过若干典型算例对算法与程序的考核与验证。
其他文献
温压炸药爆炸后，产物由高爆炸药组分分解产生的气体和未反应的固体颗粒组成。本文从气固两相流的角度，采用耦合的有限差分-物质点法，对温压炸药的爆炸问题进行数值模拟。其中，有限差分法计算气相运动，物质点法描述颗粒相的运动和反应过程，气相和颗粒相通过耦合算法实现相互作用。编制了一维计算程序，计算了含铝炸药的空爆问题，给出了冲击波的传播过程以及铝颗粒在流场中的运动反应过程。计算结果表明，该算法能够正确地描述
有化学反应的气相爆轰是一个非常复杂的过程。由于其高度的不可控性,造成经济或人员的损失,所以对爆轰的风险评估很有必要,而对于不同方式起爆的爆轰机理进行研究便成为了重要领域之一。通常来说,爆轰起爆可以分为两大类：非直接起爆和直接起爆。非直接起爆对应于反应混合物被一个弱能量源点燃,紧随其后的是一个爆燃转爆轰的过程(DDT)。对于直接起爆,其定义为通过一个小体积无约束的可燃混合物急速释放大量的能量瞬间形成
为了获得爆炸波入射至CL-20/TNT共晶含能材料后的初始冲击引爆细节,利用刚体活塞冲击压缩共晶含能材料l ps,随后撤销活塞以激发脉冲波的产生和传播.活塞的初始冲击速度up分别为4.0 km/s和6.0 krn/s,用以比较不同强度冲击波引发共晶含能材料的初始化学反应特征.在冲击波传播特征方面,压缩波入射至含能材料后,受压区材料发生化学反应并向左膨胀.up=4.0 km/s和up=6.0 km/
接触-碰撞问题广泛存在于实际工程,如冲压成形、汽车碰撞、弹体穿甲/侵彻等。大量的数值模拟计算表明,接触界面模拟是影响数值计算效率和精度的重要因素。因此,发展高效、精确、健壮的接触-碰撞算法具有重要工程意义。有限元方法中,接触-碰撞的计算可以分解为三个步骤：全局搜索、局部搜索和接触约束施加。通常,局部搜索和约束施加算法决定了接触计算的精度,而全局搜索决定了接触计算效率。为提高接触搜索效率,本文以接触
爆轰加载下膨胀的HR2柱壳中的破坏是一个复杂的过程。实验结果显示膨胀HR2柱壳中破坏包括多种损伤机制,如径向层裂、环向拉伸断裂和剪切带破坏等。本文中,我们采用了CHAP3D程序,对爆轰驱动膨胀的HR2厚壁圆筒中的多阶段破坏行为进行了数值模拟研究。作为一个材料不稳定过程,HR2柱壳内破坏的形成与材料的不均匀性息息相关。我们采用了一个满足Gauss分布的初始屈服强度的扰动来描述材料的初始不均匀性,同时
不同密度流体的扰动界面在冲击波或者加速度作用下会发生不稳定性增长，这种现象称为Richtmyer-Meshkov不稳定性和Rayleigh-Taylor不稳定性。如果不稳定性问题发生在金属材料界面上，相应地称为金属材料界面不稳定性。类似于流体界面不稳定性，金属材料界面不稳定性也会形成尖钉和气泡结构；但是和流体界面不稳定性的最大不同之处在于金属材料具有强度效应，而已有的研究表明金属材料强度能够抑制扰
为了更清晰地理解楔面诱导产生的斜爆轰波起爆区的波系结构特征，采用两步诱导-放热反应模型对其进行数值研究。重点关注来流马赫数、楔面角度、化学放热反应速率的影响，结果表明：高来流马赫数和大楔面角度时，起爆区主要是由斜激波、燃烧波系构成；低来流马赫数和小楔面角度时，起爆区主要是由斜激波、燃烧波系、一系列压缩波以及斜激波和斜爆轰波诱导出斜激波构成；更低来流马赫数时，起爆区结构发生突变，其构成主要是斜激波、
采用level-set方法计算爆轰波的传播时边界条件的处理非常重要。本文采用虚拟网格实现边界条件，采用显示方法更新边界节点的level-Set函数值。通过为角度边界节点选择已知程度更高的模板，并将所有边界点进行拓扑排序和强联通分量分解，使程序能计算形状较复杂的炸药。最后通过与精确解对比，定量验证本文程序(DSDLS)的精度。