A method for compressible multimaterial flows with condensed phase explosive detonation and airblast on unstructured grids
A method for compressible multimaterial flows with condensed phase explosive detonation and airblast on unstructured grids
复制标题
DOI：
10.1016/j.compfluid.2015.01.006
复制
发表时间：
2015
影响因子：
2.8
通讯作者：
M. Price;V. Nguyen;O. Hassan;K. Morgan
中科院分区：
文献类型：
--
作者：
M. Price;V. Nguyen;O. Hassan;K. Morgan
本文提出了一种用一般形式的状态方程模拟可压缩多物质流动的有效方法，并应用于炸药爆轰和空气爆炸问题。在三维非结构化网格上，采用体积分数型方法，对不混溶流体的多物质流动进行了数值模拟。五方程准保守系统在空间离散使用基于边缘的有限体积法与二阶精确HLLC近似黎曼求解器和时间离散与显式多级龙格库塔方法。计算模型足够强大，可以处理具有强冲击的流动，同时足够通用，可以模拟具有不同状态方程和物理状态的材料。数值试验表明，该方法对强激波和界面相互作用的准确性。采用程序燃烧法描述了凝聚相爆轰中固体未反应炸药向反应气体的转化过程。燃烧模型的准确性进行了验证，通过比较与出版的数值计算结果的流动轮廓在爆震和近场空气爆炸。通过对半球形和板形装药爆轰过程的数值模拟，研究了装药形状对空气冲击波的影响。模拟预测的压力和冲量与已发表的实验数据进行了比较。
An efficient method for the simulation of compressible multimaterial flows with a general form of equation of state is presented for explosive detonation and airblast applications. Multimaterial flows are modeled with a volume-fraction type approach for immiscible fluids governed by the compressible Euler equations on three-dimensional unstructured grids. The five-equation quasi-conservative system is discretized in space using an edge-based finite volume approach with a second-order accurate HLLC approximate Riemann solver and temporal discretization with an explicit multistage Runge–Kutta method. The computational model is robust enough to handle flows with strong shocks, while being general enough to model materials with different equations of state and physical states. Numerical tests demonstrate the accuracy of the method for strong shock and interface interactions. A program burn method is implemented to describe the conversion of solid unreacted explosive to reacted gases in condensed phase detonations. The accuracy of the burn model is validated by comparison with published numerical results of flow profiles during detonation and for near-field airblast. Numerical simulations of hemispherical and plate-shaped explosive charge detonations are performed to investigate the influence of charge shape on airblast. The predicted pressure and impulse from simulation compare well with published experimental data.