Numerical calculation of early fireball radiation spectrum in strong explosion
-
摘要: 基于强爆炸火球光辐射的多群辐射流体力学方法, 采用算子分裂方法将方程组分裂为对流项和刚性源项, 其中源项部分根据方程形式, 进一步分裂为各群内的单独求解。数值计算表明:该方法克服了直接求解过程中辐射与流体耦合所带来的强不稳定性, 时间步长大幅提高, 给出的火球光辐射能谱特征与已有规律一致。可为定量分析光辐射能谱特征提供有效手段。Abstract: On the basis of multi-group radiation hydrodynamics method of fireball radiation in strong explosion, operator splitting method is used to split the equations into convection items and source items, which are split into radiation groups due to the equation formation and solved individually.Numerical calculations show that the method used here overcomes strong instability when solving the equations directly because of the coupling items between radiation and fluid.In the meantime, the time step in the calculation is increased obviously.Fireball radiation spectrum is obtained in fine accordance with the result in the literature.
-
表 1 光子分群能量
Table 1. Photon energy of each group
g pe/eV 1 0.01~0.5 2 0.5~1.0 3 1.0~1.8 4 1.8~2.1 5 2.1~2.5 6 2.5~3.1 7 3.1~4.0 8 4.0~7.0 9 7.0~10.0 10 10.0~20.0 11 20.0~40.0 12 40.0~70.0 13 40~70 14 100~200 15 200~400 16 400~1 000 17 1 000~2 000 18 2 000~5 000 19 5 000~10 000 20 10 000~20 000 21 20 000~80 000 - - - - - - 表 2 不同高度空气初始状态参数
Table 2. Air state parameters at different altitude
h/km p/kPa ρ/(kg·m-3) T/k cs/(m·s-1) 0 103.3 1.225 297.25 343.6 10 26.5 0.414 225.90 299.5 20 5.5 0.089 219.20 295.1 -
[1] 乔登江.核爆炸物理概论[M].北京: 国防工业出版社, 2003: 225-262. [2] 屠琴芬.核爆炸火球的辐射流体力学计算中的几个问题[R].西安: 西北核技术研究所, 1986. [3] Pomraning G C. The equations of radiation hydrodynamics[J]. Astrophysical Radiation Hydrodynamics, 1986(188): 45-69. [4] Brode H L. Fireball phenomenology[R]. AD0612197, 1964. [5] Crowley B K, Glenn H D, Marks R E. An analysis of marvel: A nuclear shock-tube experiment[J]. Journal of Geophysical Research, 1971, 76(14): 3356-3374. doi: 10.1029/JB076i014p03356 [6] Marrs R E, Moss W C, Whitlock B. Thermal radiation from nuclear detonations in urban environments[R]. UCRL-TR-231593, 2007. [7] Lowrie R B, Edwards J D. Radiative shock solutions with grey nonequilibrium diffusion[J]. Shock Waves, 2008, 18(2): 129-143. [8] 陈健华, 王心正, 谢龙生, 等.均匀大气中的强爆炸一维辐射流体力学数值解[J].爆炸与冲击, 1981, 1(2): 37-49.Chen Jian-hua, Wang Xin-zheng, Xie Long-sheng, et al. A one-dimensional radiation hydrodynamic numerical solution for a strong explosion in uniform atmosphere[J]. Explosion and Shock Waves, 1981, 1(2): 37-49. [9] 田宙, 乔登江, 郭永辉.不同当量强爆炸早期火球现象的数值模拟[J].爆炸与冲击, 2009, 29(4): 408-412.Tian Zhou, Qiao Deng-jiang, Guo Yong-hui. Numerical simulation on early fireball phenomenology of strong explosions for different yields[J]. Explosion and Shock Waves, 2009, 29(4): 408-412. [10] 田宙, 乔登江, 郭永辉.不同高度强爆炸早期火球数值研究[J].兵工学报, 2009, 30(8): 1078-1083.Tian Zhou, Qiao Deng-jiang, Guo Yong-hui. Numerical investigation of early fireball of strong explosion for different altitudes[J]. Acta Armamentarii, 2009, 30(8): 1078-1083. [11] 田宙, 乔登江, 郭永辉.强爆炸早期火球现象的一维数值研究[J].计算物理, 2010, 27(1): 9-14.Tian Zhou, Qiao Deng-jiang, Guo Yong-hui. A one dimensional numerical study on early fireball in strong explosion[J]. Chinese Journal of Computational Physics, 2010, 27(1): 9-14. [12] Symbalisty E M D, Zinn J, Whitaker R W. Radflo physics and algorithms[R]. LA-12988-MS, 1995. [13] 高银军, 田宙, 刘峰, 等.强爆炸早期火球光辐射能谱的分群计算[J].四川兵工学报, 2011, 32(3): 21-24.Gao Yin-jun, Tian Zhou, Liu Feng, et al. Calculation of energy spectrum of early fireball radiation in strong explosion with multi-group method[J]. Journal of Sichuan Ordnance, 2011, 32(3): 21-24. [14] 闫凯.二维辐射流体动力学方程组的数值求解[D].西安: 西北核技术研究所, 2011. [15] 闫凯.二维辐射流体动力学方程组的数值求解[C]//第六届全国青年计算物理学术会议.太原, 2011. [16] Minerbo G N. Maximum entropy eddington factors[J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 1977, 20(6): 541-545. [17] 王文高, 张建泉.物质辐射不透明性Ⅰ[R].西安: 西北核技术研究所, 1978. [18] 王文高, 张建泉.物质辐射不透明性Ⅱ[R].西安: 西北核技术研究所, 1978.