摘要

针对爆炸冲击波数值模拟中经常要面临的强间断奇异性问题,在一维或多维空间中引入弧长参数,通过添加一个约束方程,致使原网格自适应移动变形而能够在强间断区域分布较多的网格点,以此提高大梯度物理量捕捉的分辨率。针对非均匀变形物理空间中格式构造困难的问题,提出了坐标变换,将物理空间中的控制方程映射至均匀正交的弧长计算空间。基于MUSCL格式编写了伪弧长算法相应的一维、二维程序,并对其开展人为解方法验证。将伪弧长算法(PALM)的数值结果同有限体积法进行比较,经过误差分析,验证了伪弧长算法在间断附近具有较高的分辨率。