激波聚焦诱导点火和爆轰的数值研究

以二维轴对称多组分Euler方程为基础,采用非正交结构化网格和改进的波传播算法,模拟了激波在抛物形反射壁面聚焦反射诱导点火和爆轰的过程,描述了其流场形态。讨论了预混气组成、入射激波强度及反射壁面形状对点火和爆轰的影响。结果表明,激波在抛物形反射壁面顶点处聚焦反射可形成局部高温高压区域,该区域在一定条件下可点燃预混气甚至形成爆轰,其中低稀释剂浓度的预混气、较大的入射激波Mach数和较深的反射壁面有利于可燃预混气形成爆轰。     

激波与火焰的相互作用过程

基于带化学反应的Navier Stokes方程和有关的热力学和反应动力学数据,利用改进的VLS格式,对甲烷 空气混合物中激波与火焰的相互作用进行了数值模拟。根据计算结果,讨论了激波掠过火焰时的变形、分叉和发展,以及激波作用下火焰的失稳、变形、破碎和相应的带旋涡的流场。     

非定常激波三维双楔面反射的数值研究

采用基于MUSCL-Hancock插值的有限体积方法, 在非结构自适应网格上求解三维Euler方程, 研究了非定常激波在三维垂直双楔面上的反射现象. 研究结果表明, 由于三维效应的影响, 通过二维非定常激波反射理论预测的三维马赫反射区域范围存在着一定的局限, 在该区域范围内出现了一种非典型的三维突起结构, 即第二类三维马赫干. 另外, 对于不同的激波马赫数和楔面倾角组合, 非定常激波在三维双楔面上会形成具有四波结构的三维马赫反射或具有三波结构的三维规则反射.      

激波绕半菱形柱体运动的实验研究和数值模拟

激波绕过半菱形柱体流场是十分复杂的,本文用实验方法和数值模拟对该问题作了系统研究。实验是在UTIAS激波管上进行的,数值模拟采用2阶精度Godunow格式,最后实验结果和数值模拟作了比较。     

激波和湍流相互作用的数值模拟

本文综述了关于激波和湍流相互作用数值模拟的近期研究进展, 主要包括激波和均匀各向同性湍流、激波和湍流边界层、激波和射流以及激波和尾迹的相互作用. 激波和湍流相互作用特性受到诸多因素的影响,如激波的强度、位置、形状和流动边界以及来流的湍流状态和可压缩性等. 激波和湍流的相互作用会引起流场结构、激波特性和湍流统计特性的显著变化. 最后简要讨论了激波和湍流相互作用数值研究需要关注的一些问题.      

二维激波与剪切层相互作用的直接数值模拟研究

采用五阶weighed esseritially non-oscillatory (WENO) 格式和三阶total variation diminishing (TVD) Runge-Kutta 格式, 通过求解二维非定常Navier-Stokes 方程, 直接数值模拟了激波与剪切层相互作用, 目的在于揭示激波与剪切层相互作用过程中噪声产生的机理. 研究发现:(1) 当入射激波穿过剪切层时, 剪切层中心位置向下层区域偏移;(2) 入射激波穿过剪切层产生小激波, 在小激波与剪切层接触点处产生声波并向外辐射;(3) 反射激波穿过剪切层后形成了分段弧状激波;(4) 当反射激波穿过剪切层时, 激波在鞍点处泄漏并向外辐射声波, 这是一种激波泄漏机制.     

可燃气体中激波与障碍物作用在下游形成爆轰波的数值研究

对平面激波和单个矩形障碍物作用的过程进行了数值模拟,研究了反射产生的上行爆轰波在下游可燃气体中形成爆轰波的过程。数值结果表明,下游爆轰波形成过程主要有2种模式:爆轰波直接绕射和绕射波在上壁面反射,这和已有的实验结果是一致的。通过研究下游爆轰波的形成过程受入射激波马赫数、混合气体的压力和管道尺度的影响,分析了上游爆轰波向下游传播的波动力学过程,讨论了2种形成过程的作用规律和控制因素,阐明了下游爆轰波的形成规律。     

激波在不同三角楔面内传播特性的数值研究

激波在收缩管内的反射与聚焦会形成高温高压区,点燃可燃混合气并诱导爆轰,因此对爆轰发动机的点火具有重要意义。本文基于二维N-S方程,结合五阶WENO格式,对马赫数为6的正激波在三角形楔面内的反射与聚焦现象进行了数值研究。结果表明,楔面顶角的变化对激波的反射类型以及聚焦均有明显的影响：随着顶角的增加,激波的反射类型从马赫反射向过渡马赫反射和双马赫反射转变,且壁面上的前向射流更加明显;三波点第一次碰撞产生的高温高压区足够满足可燃混合气体的点火条件,且其温度与压力值随顶角的增加而增大;当激波在楔面上发生临界双马赫反射时,温度与压力达到最大;当顶角增加到一定值时,激波在楔面反射转变为常规反射,不会产生激波对碰,因而没有高温高压区。     

用环形激波聚焦实现爆轰波直接起爆的数值模拟

利用基元反应模型和有限体积法对环形激波在可燃气体中聚焦实现爆轰波直接起爆进行了数值模拟。研究结果表明,标准状态下的氢气-空气混合气体在马赫数为3.1以上的环形激波聚焦产生的高温高压区作用下会诱发可燃气体的直接起爆形成爆轰波,爆轰波与激波和接触间断相互作用产生了复杂的波系结构;爆轰波爆点位置在对称轴上并不是固定的点,而是随着初始激波马赫数的变化而发生移动;可燃气体初始温度和压力对起爆临界马赫数都有影响,但是初始温度的影响大得多。     

水中爆炸激波对水泥试样作用的数值模拟分析

通过对水中爆炸激波在水泥净浆试样中传播的数值模拟,再现了爆炸激波的传播过程,采用对典型单元受到的激波作用和应力进行分析的方法,得出了水泥试样各个区域损伤破坏的成因,数值模拟结果和实验现象吻合。     

水中爆炸冲击波载荷作用下舰船结构动态响应的数值模拟

采用大型通用有限元软件MSC.DYTRAN,对某型水面舰船全船结构在水下爆炸冲击波载荷作用下的动态响应进行了数值模拟。分别从结构的变形损伤形式、能量吸收和冲击环境等几个方面研究了舰船结构在水下爆炸载荷作用下的破坏机理和响应特征。     

浅层水中沉底的两个装药爆炸的数值模拟研究

根据试验模型和试验结果,进行了浅层水中沉底的两个装药爆炸的数值模拟;通过与水中单个装药爆炸,以及无限水中两个装药同时爆炸的数值模拟结果的对比分析,研究了水底水面对沉底的两个装药同时爆炸产生的冲击波传播与相互作用的影响。结果表明:水底对冲击波压力峰值有削弱作用,水面使冲击作用冲量明显减小,冲击波相互作用压力叠加或多次冲击作用可提高爆炸威力。     

水下爆炸冲击波的数值模拟研究

应用商业有限元程序MSC.DYTRAN数值模拟了球形药包在无限水域中爆炸产生的冲击波。通过和经验公式计算结果的比较,证明采用合理的计算参数和有限元模型能够较好地模拟水下爆炸冲击波的传播过程。通过调整水的状态方程参数,达到了提升冲击波峰值应力的效果,从而有效地降低了单元数量和计算时间。     

准定常强激波反射马赫杆突出变形准则的探讨

分析了滑移线延长线与楔面交点附近的流动特点,研究了准定常强激波反射中马赫杆的变形,建立了马赫杆突出变形消失的过渡准则,并用激波极曲线方法进行了求解。对马赫杆的变形过程进行了描述,研究了比热比、马赫数的变化对马赫杆突出变形消失条件的影响,并对突出变形消失区域与终点双马赫反射区域进行了比较。结果表明：马赫杆突出变形的消失是低比热比介质中出现的一种激波反射现象。马赫杆突出变形的消失,导致出现一种新的无射流、直马赫杆的双马赫反射结构。     

低冲击下固体炸药延迟起爆(XDT)现象

设计了连续低冲击加载炸药的双隔板实验,采用X光透视摄影方法,观察到了JO 9159炸药延迟起爆(XDT)现象。采用非线性有限元方法对加载炸药实验的力学过程进行了数值模拟,分析了炸药延迟起爆机理和发生条件。结果表明,JO 9159炸药延迟起爆(XDT)现象的起爆时刻比正常冲击起爆(SDT)要滞后大约十几微秒以上,且发生有一定的随机性。对炸药重复冲击加载作用和炸药在加载作用下的敏化程度是控制XDT现象发生的主要因素。     

坑道内化爆冲击波的传播规律

装药在坑道内部爆炸时将产生沿坑道传播的空气冲击波,由于受到坑道洞壁的限制,空气冲击波在坑道内的传播规律异于在自由大气中的传播。利用三维数值模拟计算程序,对长坑道中的化爆流场进行了数值计算,归纳出空气冲击波沿坑道方向的传播规律。计算结果与试验结果符合较好。     

冲击载荷下孔隙塌缩过程的数值模拟

爆炸压实过程中多孔体的孔隙闭合程度对压实效果起着决定性作用。利用LS-DYNA有限元程序,对无氧铜中的圆形孔隙塌缩过程进行了数值模拟。根据模拟结果分析得出,在6 GPa的冲击压力下,孔隙闭合时不同边界区域会发生爆炸焊接和射流侵彻,这2种结合机理可以使材料结合更牢固,提高材料的致密度和机械强度,实现高质量的爆炸压实。