强爆炸火球辐射流体自适应网格高精度数值模拟

doi:10.19596/j.cnki.1001-246x.8230

计算物理 ›› 2021, Vol. 38 ›› Issue (2): 146-152.DOI: 10.19596/j.cnki.1001-246x.8230

强爆炸火球辐射流体自适应网格高精度数值模拟

李康^1,2, 李守先², 刘娜^1,2

1. 中物院高性能数值模拟软件中心, 北京 100088;
2. 北京应用物理与计算数学研究所, 北京 100088

收稿日期:2020-05-08 修回日期:2020-08-01 发布日期:2021-09-29
作者简介:李康(1988-),男,副研究员,博士,主要从事多介质辐射流体高精度数值模拟研究,E-mail:kanglimech@hotmail.com
基金资助:
国家自然科学基金（11971070，U1830139）资助项目

High-precision Numerical Simulation of Strong Explosion Fireball with Adaptive Mesh

LI Kang^1,2, LI Shouxian², LIU Na^1,2

1. Software Center for High Performance Numerical Simulation, China Academy of Engineering Physics, Beijing 100088, China;
2. Institute of Applied Physics and Computational Mathematics, Beijing 100088, China

Received:2020-05-08 Revised:2020-08-01 Published:2021-09-29

摘要/Abstract

摘要： 提出一种自适应网格方法，应用于基于Euler方法的强爆炸辐射流体高精度数值求解。通过与Zinn数值结果对比，验证该方法的正确性。研究自适应网格对冲击波和光辐射输出模拟精度的影响，对比不同网格尺度下的计算耗时。在相同的条件下，使用自适应网格与均匀网格加密3倍得到的冲击波超压分布、光辐射输出演化接近，计算效率提升约8.5倍。由此可知，提出的自适应网格方法可用于强爆炸问题的高精度数值模拟，显著地提升模拟精度和模拟效率。

关键词: 辐射流体动力学, 火球, Euler方法, 高精度模拟, 自适应网格

Abstract: A method of adaptive mesh is proposed and simulation of strong explosion radiation hydrodynamics based on Euler method is carried out. To demonstrate feasibility of the approach, comparisons with Zinn's numerical results are made. It shows that profiles of over-pressure and visible output power with adaptive mesh is close to those with uniform mesh with three times grid number. Meanwhile, the time consumed decrease to 1/8.5. It shows that the adaptive mesh method is feasible in high-precision numerical simulation and accuracy and efficiency are increased prominently.

Key words: radiation hydrodynamics, fireball, Euler method, high-precision numerical method, adaptive mesh

中图分类号:

李康, 李守先, 刘娜. 强爆炸火球辐射流体自适应网格高精度数值模拟[J]. 计算物理, 2021, 38(2): 146-152.

LI Kang, LI Shouxian, LIU Na. High-precision Numerical Simulation of Strong Explosion Fireball with Adaptive Mesh[J]. Chinese Journal of Computational Physics, 2021, 38(2): 146-152.

参考文献

[1] 乔登江. 核爆炸火球物理[J]. 物理学进展, 1982, 3(2):236-267.
[2] 乔登江. 核爆炸物理概论[M]. 北京:国防工业出版社, 2003.
[3] 孙景文. 高空核爆炸和美苏高空核试验的评述[R]. 中国核科技报告, 1999.
[4] 王成, 舒其望. 爆炸力学高精度数值模拟研究进展[J]. 科学通报, 2015, 60(10):882-898.
[5] JOHNSEN E, COLONIUS T. Implementation of WENO schemes in compressible multicomponent flow problems[J]. Journal of Computational Physics, 2006, 219:715-732.
[6] LI L, LIU X, LI X, et al. On statistical spectral property of nonlinear weights in high order schemes[J]. Chinese Journal of Computational Physics, 2019, 36(2):127-140.
[7] ZHAO F, PAN L, WANG S. Weighted essentially non-oscillatory schemes on unstructured quadrilateral meshes[J]. Chinese Journal of Computational Physics, 2018, 35(5):525-534.
[8] BRODE H L, ASANO A, PLEMMONS M, et al. A program for calculating radiation flow and hydrodynamic motion[R]. The RAND Corporation, RM-5187-PR, 1967.
[9] ZINN J. A finite difference scheme for time-dependent spherical radiation hydrodynamics problems[J]. Journal of Computational Physics, 1973, 13:569-590.
[10] ZINN J, SUTHERLAND C D. Special numerics for a nuclear fireball model[R]. Los Alamos National Laboratory, LA-9413-MS, UC-34, 1982.
[11] 毛孝勇, 王敏, 李茂生, 等. 空爆X光在大气中转换的数值模拟研究[R]. 第三卷,中国科学技术进展报告, 2013.
[12] 陈建华, 王心正, 谢龙生, 等. 均匀大气中的强爆炸一维辐射流体力学数值解[J]. 爆炸与冲击, 1981, 2:37-49.
[13] WANG X, SUI W. Two-dimension radiation hydrodynamics calculation of the high-altitude fireball[J]. Chinese Journal of Computational Physics, 1987, 4(2):159-168.
[14] 田宙, 乔登江, 郭永辉. 不同当量强爆炸早期火球现象的数值模拟[J]. 爆炸与冲击, 2009, 29(4):408-412.
[15] 闫凯, 田宙, 郭永辉, 等. M1近似应用于强爆炸火球辐射输运数值模拟[J]. 爆炸与冲击, 2014, 34(2):241-246.
[16] TIAN Z, QIAO D, GUO Y. A one-dimensional numerical study on early fireball in strong explosion[J]. Chinese Journal of Computational Physics, 2010, 27(1):8-14.
[17] YAN K, LI R, TIAN Z, et al. An operator splitting method for numerical simulation of strong explosion fireball[J]. Chinese Journal of Computational Physics, 2013, 30(3):379-386.
[18] 高银军, 闫凯, 田宙, 等. 强爆炸早期火球光辐射能谱的数值计算[J]. 爆炸与冲击, 2015, 35(3):289-295.
[19] 高银军, 田宙, 闫凯, 等. 强爆炸光辐射脉冲辐照特征与爆炸当量的相关性[J]. 爆炸与冲击, 2017, 37(3):549-553.
[20] LIU X, XU X, ZHANG A. A fast box set subtranction algorithm for parallel structured adaptive mesh refinement applications[J]. Chinese Journal of Computational Physics, 2017, 34(5):563-573.
[21] TANG H, TANG T. Adaptive mesh methods for one-and two-dimensional hyperbolic conservation laws[J]. SIAM J NUMER ANAL, 2003, 41(2):487-515.
[22] SYMBALISTY E M, ZINN J, et al. Physics and algorithms[R]. Los Alamos National Laboratory, LA-12988-MS, UC-700, 1995.

强爆炸火球辐射流体自适应网格高精度数值模拟

High-precision Numerical Simulation of Strong Explosion Fireball with Adaptive Mesh

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	徐建军, 史卫东, 李兴伟, 舒适. 块结构自适应网格上任意区域和任意界面的数值积分[J]. 计算物理, 2017, 34(1): 10-18.
[2]	王路, 徐江荣, 刘保银. 基于Lagrange方法封闭的两相湍流场方程模型[J]. 计算物理, 2016, 33(3): 305-310.
[3]	宫翔飞, 杨基明, 张树道. 使用AMR型背景网格的并行SPH方法[J]. 计算物理, 2016, 33(2): 183-189.
[4]	孙文俊, 范征峰. 辐射流体力学的分子动理学自适应加密方法(英文)[J]. 计算物理, 2015, 32(3): 277-292.
[5]	舒适, 岳孝强, 周志阳, 徐小文. SAMR网格上扩散方程有限体格式的逼近性与两层网格算法[J]. 计算物理, 2014, 31(4): 390-402.
[6]	闫凯, 李若, 田宙, 郭永辉, 曹渊. 强爆炸火球数值模拟中的算子分裂方法[J]. 计算物理, 2013, 30(3): 379-386.
[7]	田宙, 乔登江, 郭永辉. 强爆炸早期火球现象的一维数值研究[J]. 计算物理, 2010, 27(1): 8-14.
[8]	刘国昭, 张树道. 气相爆轰高阶中心差分-WENO组合格式自适应网格方法[J]. 计算物理, 2008, 25(4): 387-395.
[9]	马天宝, 郝莉, 宁建国. Euler多物质流体动力学数值方法中的界面处理算法[J]. 计算物理, 2008, 25(2): 133-138.
[10]	柏劲松, 陈森华, 李平. 多介质流界面高精度自适应欧拉算法[J]. 计算物理, 2003, 20(2): 95-101.
[11]	康红文, 柳崇健, 徐祥德. 一维无迭代自适应网格技术在Burgers方程中的应用[J]. 计算物理, 2002, 19(6): 553-556.
[12]	康红文, 王鹏云, 徐祥德. 自适应网格技术及其在二维动力锋生成过程模拟上的应用研究[J]. 计算物理, 2000, 17(S1): 137-148.
[13]	周春华. 非匹配网格上广义Stokes问题的区域分裂解法及事后误差估算[J]. 计算物理, 2000, 17(4): 372-380.
[14]	徐涛, 水鸿寿. 一种二维自适应网格构造方法及其实现[J]. 计算物理, 1999, 16(1): 66-76.
[15]	朱广圣, 吴锤结. 应用自适应网格方法研究着火和燃烧反应流动[J]. 计算物理, 1997, 14(6): 846-850.