欧拉坐标系下具有锐利相界面的可压缩多介质流动数值方法研究

摘要

可压缩多介质流动问题的数值模拟在国防和工业领域内均具有重要的研究价值,诸如武器设计、爆炸安全防护等,通常具有大变形、高度非线性等特点,是一项极具挑战性的研究课题.本文提出了一种基于Euler坐标系的非结构网格、具有锐利相界面的二维和三维守恒型多介质流动数值方法,可用于模拟可压缩流体和弹塑性固体在极端物理条件下的大变形动力学行为.利用分片线性的水平集函数重构出单纯形网格内分段线性的相界面,并在混合网格内构建出具有多种介质的相界面几何结构,理论上可以处理全局任意种介质、局部3种介质的多介质流动问题.利用传统的有限体积格式来计算单元边界上同种介质间的数值通量,并通过在相界面法向上求解局部一维多介质Riemann问题的精确解来计算不同介质间的数值通量,保证了相界面上的通量守恒.提出了一种非结构网格上的单元聚合算法,消除了由于网格被相界面分割成较小碎片、违反CFL条件,进而可能带来数值不稳定的问题.针对一维多介质Riemann问题、激波与气泡相互作用问题、浅埋爆炸问题、空中强爆炸冲击波和典型坑道内冲击波传播问题开展了数值模拟研究,将计算结果与相关的理论、实验结果进行比对,验证了数值方法的正确性和可靠性.