基于伪弧长移动网格算法的爆炸与冲击多介质问题数值模拟  

Numerical Simulation of Explosion and Shock Problems of Multi-material Flows by Pseudo-arc Length Moving Mesh Algorithm

作　　者：赵金庆 马天宝[1] ZHAO Jinqing;MA Tianbao(State Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology,Beijing 100081,China)

机构地区：[1]北京理工大学爆炸科学与技术国家重点实验室,北京100081

出　　处：《兵工学报》2020年第S02期200-210,共11页Acta Armamentarii

基　　金：国家自然科学基金项目(11822203、11532012);爆炸科学与技术国家重点实验室自主课题项目(YBKT18-01)。

摘　　要：在爆炸与冲击问题中,大量物理现象都是由多种介质相互作用产生的。基于伪弧长移动网格(PALM)算法,结合五方程模型对爆炸与冲击多介质问题进行数值模拟,提出多介质伪弧长移动网格(MM-PALM)算法。该算法包括PALM算法与人工界面压缩技术。爆炸与冲击多介质问题采用由Mie-Grüneisen一般状态方程表示的非守恒形式五方程模型描述;控制方程采用2阶MUSCL格式与3阶TVD Runge-Kutta方法进行数值离散。利用PALM算法网格自适应的优点,降低爆炸与冲击问题数值计算过程中强冲击间断处的奇异性,并实现对爆炸冲击波阵面的高效精确捕捉。在多介质爆炸与冲击问题计算过程中加入人工界面压缩技术,以保持多介质之间的精确界面描述。一维激波管问题和二维水下爆炸问题的研究结果表明,MM-PALM算法是有效的,能够精确捕捉冲击波阵面与多介质间的界面。A large number of physical phenomena in explosion and impact are generated by the interaction of various materials.A multi-material pseudo arc-length algorithm is proposed to numerically simulate the explosion and impact of multi-material flows,which is described by non-conserved five-equation model with the Mie-Grüneisen equation of state.Multi-material pseudo arc-length moving mesh(MM-PALM)algorithm consists of pseudo arc-length moving mesh scheme and interface-sharpening technique.The five-equation non-conserved system is discretized in space by the finite volume approach with a second-order accurate MUSCL scheme and the temporal discretization with a multistage TVD Runge-Kutta method.Pseudo arc-length moving mesh scheme is used to reduce the singularity of strong impact discontinuity and capture the wave-front of explosive shock wave accurately.The interface-sharpening technique is adopted to maintain a precise interface between multi-material flows.The research results of one-dimensional shock tube problem and two-dimensional underwater explosion problem show that the MM-PALM algorithm is feasible,and research can accurately capture the wave front of explosive shock wave and describe the interface between multi-material flows.

关 键 词：爆炸与冲击 伪弧长移动网格算法 多介质流体 人工界面压缩技术 数值模拟 

分 类 号：O389[理学—流体力学]