一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法

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

计算物理 ›› 2019, Vol. 36 ›› Issue (5): 505-516.DOI: 10.19596/j.cnki.1001-246x.7913

• • 下一篇

一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法

李诗尧¹, 于明²

1. 中国工程物理研究院研究生院, 北京 100083;
2. 北京应用物理与计算数学研究所, 北京 100094

收稿日期:2018-06-25 修回日期:2018-09-28 出版日期:2019-09-25 发布日期:2019-09-25
通讯作者: 于明,E-mail:yu_ming@iapcm.ac.cn
作者简介:李诗尧(1994-),男,硕士研究生,主要研究计算流体力学和凝聚炸药爆轰的相关算法
基金资助:
国家自然科学基金（11772066）及核科学挑战专题（TZZT2017-A1-F）资助项目

A Cell-centered Lagrangian Scheme Based on Characteristics Theory for Condensed Explosive Detonation

LI Shiyao¹, YU Ming²

1. Graduate School of China Academy of Engineering Physics, Beijing 100083, China;
2. Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

Received:2018-06-25 Revised:2018-09-28 Online:2019-09-25 Published:2019-09-25

摘要/Abstract

摘要： 提出一种数值模拟凝聚炸药爆轰问题的单元中心型Lagrange方法.利用有限体积离散爆轰反应流动方程组，基于双曲型偏微分方程组的特征理论获得离散网格节点的速度与压力，获得的网格节点速度与压力用于更新网格节点位置以及计算网格单元边的数值通量.以这种方式获得的网格节点解是一种"真正多维"的理论解，是一维Godunov格式在二维Riemann问题的推广.有限体积离散得到的爆轰反应流动的半离散系统使用一种显-隐Runge-Kutta格式来离散求解：显式格式处理对流项，隐式格式处理化学反应刚性源项.算例表明，提出的单元中心型Lagrange方法能够较好地模拟凝聚炸药的爆轰反应流动.

关键词: 凝聚炸药爆轰, 单元中心型Lagrange方法, 特征理论, 化学反应率

Abstract: We present a cell-centered Lagrangian scheme for numerical simulation of condensed-explosive detonation. It utilizes finite volume method to discrete detonation equations. Velocity and pressure of grid nodes are obtained with characteristics theory of hyperbolic partial differential equations. They are used to update position of grid nodes and calculate numerical flux of grid cells. Solution of grid nodes obtained by characteristics theory is a "genuinely multi-dimensional" theoretical solution, which is a generalization of one-dimensional Godunov scheme in two-dimensional Riemann problem. Semi-discrete system of detonation equations obtained by finite volume scheme is solved with an implicit-explicit Runge-Kutta method. Convection terms are explicitly treated, and stiff source terms of chemical reactions are implicitly treated. The numerical schemes select ZND model for detonation, JWL equations of state for unreacted explosives and detonation products, and use Ignition-Growth model to simulate evolution process in reaction zone. Typical numerical results show that the cell-centered Lagrangian scheme simulates condensed-explosive detonation problems well.

Key words: condensed-explosive detonation, cell-centered Lagrangian scheme, characteristics theory, chemical reaction rate

中图分类号:

O381

李诗尧, 于明. 一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法[J]. 计算物理, 2019, 36(5): 505-516.

LI Shiyao, YU Ming. A Cell-centered Lagrangian Scheme Based on Characteristics Theory for Condensed Explosive Detonation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 36(5): 505-516.

参考文献

[1] MADER C L. Numerical modeling of explosives and propellants[M]. 2nd ed. New York:CRC Press, 1998.
[2] BENSON D J. Computational methods in Lagrangian and Eulerian hydrocodes[J]. Computational Methods of Applied Mechanics and Engineering, 1992, 99:235-394.
[3] VON NEUMANN J, RICHTMYER R D. A method for the numerical calculations of hydrodynamical shocks[J]. Journal of Applied Physics, 1950, 21:232-238.
[4] WILKINS M L. Calculation of elastic-plastic flow:Methods in computational physics[M]. New York:Academic Press, 1964.
[5] JONAS Z. Introduction to hydrocodes[M]. Amsterdam:Elsevier Press, 2004.
[6] CARAMANA E J, ROUSCULP C L, BURTON D E. A compatible, energy and symmetry preserving Lagrangian hydrodynamics algorithm in three-dimensional Cartesian geometry[J]. Journal of Computational Physics, 2000, 157:89-119.
[7] GODUNOV S K. Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics[J]. Matematicheskii Sbornik, 1959, 47(3):271-306.
[8] ADDESSIO F L, CARROLL D E, DUKOWICZ J K, et al. CAVEAT:A computer code for fluid dynamics problems with large distortion and internal slip[R]. Los Alamos National Laboratory, UC-32, 1988.
[9] CHENG J, SHU C W. A high order ENO conservative Lagrangian type scheme for the compressible Euler equations[J]. Journal of Computational Physics, 2007, 227:1567-1596.
[10] DESPRÉS B, MAZERAN C. Lagrangian gas dynamics in two dimensions and Lagrangian systems[J]. Archive for Rational Mechanics and Analysis, 2005, 178:327-372.
[11] MAIRE P H. A cell-centered Lagrangian scheme for two-dimensional compressible flow problems[J]. SIAM Journal on Scientific Computing, 2007, 29(4):1781-1824.
[12] MAIRE P H, BREIL J. A second-order cell-centered Lagrangian scheme for two-dimensional compressible flow problems[J]. International Journal for Numerical Methods in Fluids, 2007(56):1417-1423.
[13] MAIRE P H. A high-order cell-centered Lagrangian scheme for two-dimensional compressible fluid flows on unstructured meshes[J]. Journal of Computational Physics, 2009, 228:2391-2425.
[14] MAIRE P H, LOUBERE R. Staggered Lagrangian discretization based on cell-centered Riemann solver and associated hydrodynamics scheme[J]. Communications of Computational Physics, 2011, 10(4):940-978.
[15] ROE P L. Discrete models for the numerical analysis of time-dependent multidimensional gas dynamics[J]. Journal of Computational Physics, 1986,63:458-476.
[16] MORTON K W, ROE P L. Vorticity-preserving Lax-Wendroff-type schemes for the system wave equation[J]. SIAM Journal on Scientific Computing, 2001, 23(1):170-192.
[17] LUNG T B, ROE P L. Toward a reduction of mesh imprinting[J]. International Journal for Numerical Methods in Fluids, 2014, 76:450-470.
[18] PRASAD P. Nonlinear hyperbolic waves in multi-dimensions[M]. New York:Chapman & Hall/CRC, 2001.
[19] LUKÁČOVÁ-MEDVID'OVÁ M, SAIBERTOVÁ J, WARNECKE G. Finite volume evolution Galerkin methods for nonlinear hyperbolic systems[J]. Journal of Computational Physics, 2002, 183:533-562.
[20] SUN Y T, REN Y X. The finite volume local evolution Galerkin method for solving the hyperbolic conservation laws[J]. Journal of Computational Physics, 2009, 228:4945-4960.
[21] SUN Y T, REN Y X, YU M, ZHANG S D. A finite volume method for 2D inviscous Lagrangian hydrodynamics based on characteristics theory[J]. Chinese Journal of Computational Physics, 2011, 28(1):19-26.
[22] SUN Y T, JIA Z P, YU M, REN Y X. A second order Lagrangian scheme based on characteristics theory for two-dimensional compressible flows[J]. Chinese Journal of Computational Physics, 2012, 29(6):791-798.
[23] MENIKOFF, RALPH. JWL equation of state[R]. LA-UR-15-29536, Los Alamos, 2015.
[24] TARVER C M, MCGUIRE E M. Reactive flow modeling of the interaction of TATB detonation waves with inert materials[C]. The 12th Symposium on Detonation, San Diego, 2002.
[25] CHERTOCK A, CUI S, KURGANOV A, WU T. Steady state and sign preserving semi-implicit Runge-Kutta method for ODES with stiff damping term[J]. SIAM Journal on Numerical Analysis, 2015, 53(4):2008-2029.
[26] MADER C L, ZUMBRO J D, FERM E N. Proton radiographic and numerical modeling of colliding, diverging PBX-9502 detonations[C]. The 12th Symposium on Detonation, San Diego, 2002.

一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法

A Cell-centered Lagrangian Scheme Based on Characteristics Theory for Condensed Explosive Detonation

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 3

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	孙宇涛, 贾祖朋, 于明, 任玉新. 基于特征理论的二维可压缩流动的二阶拉氏算法[J]. 计算物理, 2012, 29(6): 791-798.
[2]	孙宇涛, 任玉新, 于明, 张树道. 基于特征理论的二维无粘Lagrange流体力学有限体积法[J]. 计算物理, 2011, 28(1): 19-26.
[3]	楼建锋, 于恒. CTVD格式数值计算非均质炸药爆轰问题[J]. 计算物理, 2005, 22(4): 358-364.