不同磁场构型下Richtmyer-Meshkov不稳定性的数值研究及动态模态分解

基于磁流体动力学,本文通过数值模拟对不同磁场构型下轻质气柱界面Richtmyer-Meshkov不稳定性的演化过程进行了研究.结果显示:磁场对波系演化影响甚微,但能抑制界面不稳定性发展,且横向磁场抑制效果更好.无磁场时,界面形成涡串, SF_6射流穿过下游界面;有磁场时,界面光滑无涡串.其中,横向磁场下界面更光滑, SF_6射流不再穿过界面.此外,由于Richtmyer-Meshkov不稳定性的作用,磁力线在气柱界面发生扭曲,且上游界面处磁力线扭曲程度更大,产生强洛伦兹力,使涡量分层明显;下游界面处,纵向磁场产生的洛伦兹力较横向磁场更小,涡层之间相互干扰.最后,本文将动态模态分解用于界面不稳定性研究,发现:磁场作用下界面仍存在小涡,且纵向磁场下扰动更多.第一模态的稳定涡结构能反映主要流场信息,第二到第四模态下的小涡频率依次增加,且无磁场、纵向和横向磁场的同一模态下,小涡频率依次减小.因而磁场能抑制小涡频率,且横向磁场抑制效果更好.     

小口径膛口射流噪声的数值模拟

为了研究小口径武器的膛口气动噪声特性,采用CFD-CAA耦合算法对7.62 mm枪的射流噪声场进行了数值模拟。由于膛口流场结构复杂,在目前的计算发展水平下还不足以采用CAA直接法,因而本文中采用混合方法,即首先采用CFD方法计算7.62 mm枪的膛口流场,然后利用所得结果,采用声学方程计算射流噪声,具体为膛口近场采用LES进行计算,远场声场采用FW-H声拟法计算。通过对比验证实验,验证了该计算方法的可行性。然后,对7.62 mm枪射流噪声进行了数值模拟,分析了噪声指向性,绘制了声压级云图。研究表明:在本文的计算条件下,射流噪声强度主要集中在近膛口区域;且射流最大噪声主要分布在与轴线方向成30°~60°范围内。     

不同初始磁场对激波冲击R22重气柱过程影响的数值模拟

基于磁流体动力(magneto-hydrodynamic,MHD)方程,采用CTU+CT方法,对在不同初始磁场作用下的平面入射激波与磁化R22重质气柱作用过程进行了数值研究。数值结果清晰地描述了不同初始磁场条件时激波诱导R22气柱界面不稳定性的过程,揭示了磁场控制界面不稳定性的机理。另外,还分析了磁感应强度对界面不稳定性的影响,发现在磁场较小时,涡层附着于界面,但随着磁感应强度的增大,平均涡量随之增大,涡层与界面逐渐分离,最终更好地抑制了界面不稳定性。同时,还发现平均涡度拟能随着磁感应强度的增大而减小,而垂直磁场比平行磁场更能降低平均涡度拟能,因而平均涡度拟能可较好地反映磁场对不稳定性的影响效果。     

磁控条件下激波冲击三角形气柱过程的数值研究

本文基于磁流体动力学方程组,在保证磁场散度为零的条件下,采用CTU+CT（corner transport upwind+constrained transport）算法,对有无磁场控制下激波与重质或轻质三角形气柱相互作用过程进行数值研究.结果表明：无论有无磁场,两气柱在激波冲击下均具有完全不同的波系结构和射流现象.其中,入射激波与重气柱发生常规折射,形成介质射流,而与轻气柱作用则发生非常规折射,形成反相空气射流.无磁场时,气柱在激波冲击下,产生Richtmyer-Meshkov和Kelvin-Helmholtz不稳定性,界面出现次级涡序列,重气柱上下角卷起形成主涡对,轻气柱空气射流穿过下游界面后形成偶极子涡.施加横向磁场后,次级涡序列、主涡对以及偶极子涡均消失.进一步研究表明,在磁场作用下,洛伦兹力将不稳定性诱导产生的涡量向界面两侧的Alfvén波上输运,减少界面涡量沉积,抑制界面卷起失稳.最终,涡量沿界面两侧形成相互远离的涡层,界面不稳定性得到控制.此外,定量分析表明磁场能加快两气柱上游界面的运动,抑制下游界面的运动,且对轻气柱的控制效果更好.     

后效期火药燃气加速弹丸的数值研究

 弹丸出膛后受高压火药燃气作用将继续加速并达到最大速度,后效期内弹丸增速大小及增速规律对于武器系统的设计和射击精度等具有重要的影响。采用基于动网格的有限体积方法,对不同装药量、弹丸质量和口径等9种条件下的发射过程进行了数值模拟,讨论了弹丸出膛后的受力状况和运动过程,对最大速度出现的时间、位置和增速量及其变化特点、影响因素等进行了详细的分析。计算结果显示,后效期弹丸增速率最大可超过3%,且加速过程受弹丸质量和装药量的影响很小。     

截锥体头型弹丸低速斜入水实验研究

为研究截锥体头型弹丸在低速斜入水条件下，头部直径大小对入水空泡及弹道特性的影响，基于高速摄像方法，开展不同截锥体头型弹丸低速倾斜入水对比实验，得到了截锥体头弹丸头部直径大小对入水空泡、运动速度、俯仰角的影响规律。实验结果表明:截锥体头弹丸头部直径越大，尾部越早与空泡下壁面发生碰撞；头部直径大小对空泡深闭合时间几乎没有影响；弹丸空泡随头部直径增大而增大；头部直径过大或者过小均不利于入水稳定性；弹丸速度低于临界值时呈上升趋势，高于临界值时将呈现下降趋势。     

膛口反应流并行数值模拟

采用轴对称多组分N-S方程对含有高速运动弹丸的膛口反应流进行了数值模拟.控制方程采用时间分裂方法并在大型计算机上采用MPI方法进行多核并行求解,其中对流项采用二阶AUSM+格式和MUSCL插值方法进行处理,燃气采用氢气-空气混合气,反应机理为9组分19步基元反应.对于弹丸引起的网格运动,采用嵌套网格法处理.并行验证算例与串行计算结果一致,采用20个CPU计算时效率为64％.根据数值结果详细讨论了发射过程中的气体动力学和化学动力学过程,并且通过对两种条件下的计算结果比较分析了化学反应对膛口流场发展的影响.结果表明,上述算法能够较为正确地模拟弹丸和化学反应对膛口流场的影响,并大大提高了计算速度.