空气自由场中强爆炸冲击波传播二维数值模拟

姚成宝; 李若; 田宙; 郭永辉

doi:10.11883/1001-1455(2015)04-0585-06

空气自由场中强爆炸冲击波传播二维数值模拟

doi: 10.11883/1001-1455(2015)04-0585-06

姚成宝^{1, 2,},
李若^{1, 3},
田宙¹,
郭永辉¹

1.
北京大学数学科学学院，北京 100871
2.
西北核技术研究所，陕西西安 710024
3.
北京大学应用物理与技术研究中心，北京 100871

详细信息

作者简介:
姚成宝(1984-), 男, 博士研究生, 助理研究员, yaocheng@pku.edu.cn

中图分类号: O381
计量
- 文章访问数: 4426
- HTML全文浏览量: 615
- PDF下载量: 760
- 被引次数: 0
出版历程
- 收稿日期: 2013-12-04
- 修回日期: 2014-07-09
- 刊出日期: 2015-07-25

Two dimensional simulation for shock wave produced by strong explosion in free air

Yao Cheng-bao^{1, 2
,},
Li Ruo^{1, 3},
Tian Zhou¹,
Guo Yong-hui¹

1.
School of Mathematical Sciences, Peking University, Beijing 100871, China
2.
Northwest Institute of Nuclear Technology, Xi'an 710024, Shaanxi, China
3.
Center of Applied Physics and Technology, Peking University, Beijing 100871, China

摘要

摘要: 编写了适用于模拟具有高密度比、高压力比的强激波问题的二维柱对称多介质流体计算程序。利用有限体积方法求解流体的Euler方程组，采用level set方法捕捉爆炸产物与空气的运动界面，并通过求解物质界面两侧Riemann问题的精确解来计算爆炸产物与空气之间的数值通量。研制了三角形网格自适应技术来实现网格的自动加密和粗化，在保证捕捉激波峰值的前提下有效地提高了计算效率。利用计算程序对1 kt TNT当量的空气自由场强爆炸问题进行数值模拟，计算得到的峰值超压、冲击波到达时间等物理参数与点爆炸理论结果基本一致。
- 爆炸力学 /
- 空中强爆炸 /
- level set /
- Riemann问题 /
- 网格自适应 /
- 自由场 /
- 冲击波
Abstract: Aimed to simulate the propagation of blast wave with high density ratio and high pressure ratio produced by strong explosion in the air, a two dimensional numerical program is written in which the problem is treated as a two-medium compressible flow with sharp material interface in Eulerian grids. In this method, the finite volume method is used to solve the Euler equations, level set method is used to capture the moving interface, and the numerical flux across the interface is calculated by exactly solving the Riemann problem. Mesh adaption technique in triangle meshes is adopted to refine or coarsen the meshes which can both capture the peak overpressure and improve the computational efficiency. One kiloton nuclear charge of strong explosion in free air is simulated. The shock wave parameters, including peak overpressure, shock arrival time and so on, are in consistence with the point explosion theory, which show the accuracy and efficiency of the numerical methods.
- mechanics of explosion /
- strong explosion in air /
- level set /
- Riemann problem /
- mesh adaption /
- free air /
- shockwave

HTML全文

图 1 Riemann问题解的类型

Figure 1. Solution type of the Riemann problem

下载: 全尺寸图片幻灯片

图 2 网格自适应示意图

Figure 2. Mesh adaption

下载: 全尺寸图片幻灯片

图 3 不同时刻的冲击波压力等值线图

Figure 3. Pressure contours at different time

下载: 全尺寸图片幻灯片

图 4 峰值超压

Figure 4. Peak overpressure

下载: 全尺寸图片幻灯片

图 5 正压冲量

Figure 5. Overpressure impulse

下载: 全尺寸图片幻灯片

图 6 冲击波到时

Figure 6. Arrival time

下载: 全尺寸图片幻灯片

图 7 正压作用时间

Figure 7. Positive time duration

下载: 全尺寸图片幻灯片

图 8 不同爆心距的冲击波超压历程曲线

Figure 8. Overpressure history at different radius

下载: 全尺寸图片幻灯片

参考文献(15)

[1]	奥尔连科Л П.爆炸物理学[M].孙承纬, 译.北京: 科学出版社, 2013, 456-463.
[2]	张洪武, 何扬, 张昌权.空中爆炸冲击波地面载荷的数值模拟[J].爆炸与冲击, 1992, 12(2): 188-198. http://www.bzycj.cn/article/id/10725 Zhang Hong-wu, He Yang, Zhang Chang-quan. Numerical simulation on ground surface loading of shock wave from air explosion[J]. Explosion and Shock Waves, 1992, 12(2): 188-198. http://www.bzycj.cn/article/id/10725
[3]	岳鹏涛, 徐胜利, 彭金华. FAE爆炸波对地面目标作用的三维数值模拟[J].爆炸与冲击, 2000, 20(2): 195-201. http://www.bzycj.cn/article/id/10262 Yue Peng-tao, Xu Sheng-li, Peng Jin-hua. 3D numerical simulations on the interaction between FAE blast waves and ground targets[J]. Explosion and Shock Waves, 2000, 20(2): 195-201. http://www.bzycj.cn/article/id/10262
[4]	北京工业学院.爆炸及其作用[M].北京: 国防工业出版社, 1979.
[5]	Liu T G, Khoo B C, Yeo K S. Ghost fluid method for strong shock impactingon material interface[J]. Journal of Computational Physics, 2003, 190(2): 651-681. https://www.sciencedirect.com/science/article/pii/S0021999103003012
[6]	Fedkiw R P, Aslam T, Merriman B, et al. A non-oscillatory eulerian approach to interfaces in multimaterial flows: The ghost fluid method[J]. Journal of Computational Physics, 1999, 152(2): 457-492. https://www.sciencedirect.com/science/article/pii/S0021999199962368
[7]	Colella P, Glaz H M. Efficient solution algorithms for the Riemann problem for real gases[J]. Journal of Computational Physics, 1985, 59(2): 264-289. https://www.sciencedirect.com/science/article/pii/0021999185901469
[8]	乔登江.强爆炸物理概论(上册)[M].北京: 国防工业出版社, 2003: 51-55.
[9]	杨秀敏.爆炸冲击现象数值模拟[M].合肥: 中国科学技术大学出版社, 2010: 122-125.
[10]	Osher S, Sethian J A. Fronts propagating with curvature dependent speed: Algorithms based on Hamilton-jaccobi formulations[J]. Journal of Computational Physics, 1988, 79(1): 12-49. https://www.sciencedirect.com/science/article/pii/0021999188900022
[11]	Sussman M, Semerka P, Osher S. A level set approach for computing solutions to incompressible two-phase flow[J]. Journal of Computational Physics, 1994, 114(4): 146-159. https://www.sciencedirect.com/science/article/pii/S0021999184711557
[12]	张军, 赵宁, 任登凤, 等. level set方法在自适应Catersian网格上的应用[J].爆炸与冲击, 2008, 28(5): 156-161. Zhang Jun, Zhao Ning, Ren Deng-feng, et al. Application of the level set method on adaptive cartesian grids[J]. Explosion and Shock Waves, 2008, 28(5): 156-161.
[13]	Eleuteuo F T. Riemann solvers and numerical methods for fluid dynamics[M]. First Edition. Springer, 2009: 16-123.
[14]	刘儒勋, 舒其望.计算流体的若干新方法[M].北京: 科学出版社, 2003: 121-130.
[15]	Di Y N, Li R, Tang T, et al. level set calculations for incompressible two-phase flows on a dynamically adaptive grid[J]. Journal of Scientific Computing, 2007, 31(1/2): 75-98. doi: 10.1007/s10915-006-9119-3