简介:采用算子分裂方法将辐射流体力学方程组分裂为对流项和刚性源项,并设计了一种高效求解刚性源项方程组的数值方法.数值实验表明:该方法对时间步长不敏感,计算精度能够满足工程计算要求.在不对方程组做任意近似的情况下,数值给出了较长时间内火球物理参量的空间分布.该结果与实际强爆炸中的一些经验公式吻合较好.
强爆炸火球数值模拟中的一种算子分裂方法