摘要
采用算子分裂方法将辐射流体力学方程组分裂为对流项和刚性源项,设计一种高效求解刚性源项方程组的数值方法.数值实验表明:该方法对时间步长不敏感,计算精度能够满足工程计算要求.在不对方程组做任何近似的情况下,数值给出了较长时间内火球冲击波阵面压力峰值及阵面位置的变化,结果与实际强爆炸中的经验公式吻合较好.
An operator splitting method is developed for numerical simulation of strong explosion fireball.In the approach the system is split into two parts.One part couples radiation and fluid in a hyperbolic subsystem.Another parabolic part evolves radiation diffusion and source-sink terms.We propose a numerical method for stiff source term.It shows that the method maintains efficiency and satisfies the needs of engineering.At last,spatial distributions of parameters in nuclear explosion reflecting fireball evolution are given.The calculation agrees well with empirical formulas.
出处
《计算物理》
CSCD
北大核心
2013年第3期379-386,共8页
Chinese Journal of Computational Physics
关键词
火球
辐射流体力学
刚性源项
算子分裂方法
fireball
radiation hydrodynamics
stiff source term
operator splitting method
