Three-dimensional ghost-fluid large-scale numerical investigation on air explosion

被引:21
作者
Wang, Cheng [1 ]
Ding, JianXu [1 ]
Shu, Chi-Wang [2 ]
Li, Tao [1 ]
机构
[1] Beijing Inst Technol, State Key Lab Explos Sci & Technol, Beijing 100081, Peoples R China
[2] Brown Univ, Div Appl Math, Providence, RI 02912 USA
基金
美国国家科学基金会; 中国国家自然科学基金;
关键词
Air explosion; WENO scheme; RGFM; Local Riemann solver; Parallel computation; EFFICIENT IMPLEMENTATION; LEVEL; RECONSTRUCTION; INTERFACES; SIMULATION; FLOWS;
D O I
10.1016/j.compfluid.2016.07.015
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Based on the double shockwave approximation procedure, a local Riemann solver for strongly nonlinear equations of state (EOS) such as the Jones-Wilkins-Lee (JWL) EOS is presented, which has suppressed successfully numerical oscillation caused by high-density ratio and high-pressure ratio across the interface between explosion products and air. The real ghost fluid method (RGFM) and the level set method have been used for converting multi-medium flows into pure flows and for implicitly tracking the interface, respectively. A fifth order finite difference weighted essentially non-oscillatory (WENO) scheme and a third order TVD Runge-Kutta method are utilized for the spatial discretization and the time advance, respectively. An enclosed-type MPI-based parallel methodology for the RGFM procedure on a uniform structured mesh is presented to realize the parallelization of three-dimensional (3D) air explosion. The overall process of 3D air explosion of both TNT and aluminized explosives has been successfully simulated. The overpressures at different locations of 3D air explosion for both explosives mentioned above are monitored and analyzed for revealing the influence of aluminum powder combustion on the over-pressure of the explosion wave. Numerical results indicate that, due to aluminum powder afterburning, the attenuation of the explosion wave formed by aluminized explosives is slower than that caused by TNT. (C) 2016 Elsevier Ltd. All rights reserved.
引用
收藏
页码:70 / 79
页数:10
相关论文
共 28 条
[1]   The fast construction of extension velocities in level set methods [J].
Adalsteinsson, D ;
Sethian, JA .
JOURNAL OF COMPUTATIONAL PHYSICS, 1999, 148 (01) :2-22
[2]   Multi-material interface reconstruction on generalized polyhedral meshes [J].
Ahn, Hyung Taek ;
Shashkov, Mikhail .
JOURNAL OF COMPUTATIONAL PHYSICS, 2007, 226 (02) :2096-2132
[3]  
AMSDEN AA, 1966, LA3466 LOS AL SCI LA
[4]  
[Anonymous], 2002, Level Set Methods and Dynamic Implicit Surfaces
[5]   A partial differential equation approach to multidimensional extrapolation [J].
Aslam, TD .
JOURNAL OF COMPUTATIONAL PHYSICS, 2004, 193 (01) :349-355
[6]   Reconstruction of multi-material interfaces from moment data [J].
Dyadechko, Vadim ;
Shashkov, Mikhail .
JOURNAL OF COMPUTATIONAL PHYSICS, 2008, 227 (11) :5361-5384
[7]  
Fedkiw RP, 1999, J COMPUT PHYS, V152, P457, DOI 10.1006/jcph.1999.6136
[8]   A study of numerical methods for the level set approach [J].
Gremaud, Pierre A. ;
Kuster, Christopher M. ;
Li, Zhilin .
APPLIED NUMERICAL MATHEMATICS, 2007, 57 (5-7) :837-846
[9]   NUMERICAL CALCULATION OF TIME-DEPENDENT VISCOUS INCOMPRESSIBLE FLOW OF FLUID WITH FREE SURFACE [J].
HARLOW, FH ;
WELCH, JE .
PHYSICS OF FLUIDS, 1965, 8 (12) :2182-&
[10]   Differential equation based constrained reinitialization for level set methods [J].
Hartmann, Daniel ;
Meinke, Matthias ;
Schroeder, Wolfgang .
JOURNAL OF COMPUTATIONAL PHYSICS, 2008, 227 (14) :6821-6845