基于Mie-Grüneisen状态方程的爆轰波运动数值模拟

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

计算物理 ›› 2021, Vol. 38 ›› Issue (1): 47-56.DOI: 10.19596/j.cnki.1001-246x.8190

基于Mie-Grüneisen状态方程的爆轰波运动数值模拟

吴宗铎¹, 严谨¹, 宗智², 庞建华³, 高云⁴

1. 广东海洋大学海洋工程学院, 广东湛江 524088;
2. 大连理工大学船舶工程学院, 辽宁大连 116024;
3. 广东海洋大学深圳研究院, 智能海洋工程中心, 广东深圳 518055;
4. 西南石油大学, 油气藏地质及开发工程国家重点实验室, 四川成都 610500

收稿日期:2020-01-02 修回日期:2020-07-07 出版日期:2021-01-25 发布日期:2021-01-25
作者简介:吴宗铎(1984-),男,博士,讲师,E-mail:wuzongduo0@aliyun.com
基金资助:
国家自然科学基金（11702066，5179030）资助项目

Numerical Simulation of Detonation Wave Motion Based on Mie-Grüneisen Equation of State

WU Zongduo¹, YAN Jin¹, ZONG Zhi², PANG Jianhua³, GAO Yun⁴

1. College of Ocean Engineering, Guangdong Ocean University, Zhanjiang, Guangdong 524088, China;
2. Department of Naval Architecture, Dalian University of Technology, Dalian, Liaoning 116024, China;
3. Ocean Intelligence Technology Center, Shenzhen Institute of Guangdong Ocean University, Shenzhen, Guangdong 518055, China;
4. State Key Laboratory of Oil and Gas Reservoir Gelogy and Exploration, Southwest Petroleum University, Chengdu, Sichuan 610500, China

Received:2020-01-02 Revised:2020-07-07 Online:2021-01-25 Published:2021-01-25

摘要/Abstract

摘要： 依据C-J（Chapman-Jouguet）理论，对爆轰问题中的气态爆轰产物和未反应炸药分别考虑不同的参考状态，并根据参考状态选用特定的Mie-Grüneisen状态方程。忽略化学反应过程，爆轰产物厚度为零的前导激波面以界面的形式存在。数值模拟中，爆轰波的演化分为波面传播以及与未反应介质相互作用两个部分。传播过程中，爆轰波的传播速度即恒定的爆速，爆轰产物在传播过程中瞬间形成，而相互作用过程则是通过Mie-Grüneisen多介质混合模型来计算爆轰波的持续冲击作用。借助于Mie-Grüneisen状态方程以及Mie-Grüneisen多介质混合模型，可以很好地模拟爆轰波的运动过程。对比理论参数及文献的计算结果发现，模拟结果具备较好的准确度。

关键词: 爆轰波, C-J理论, Mie-Grü, neisen方程, 爆速

Abstract: Based on C-J (Chapman-Jouguet) theory different reference states of gaseous detonation products and unreacted explosives in detonation problem are considered. According to these reference states, specific Mie-Grüneisen EOS (equation of states) is selected. As chemical reaction process is neglected, a zero-thickness section of guided shock wave exists as an interface in front of the detonation wave. In numerical simulation, evolution of detonation wave includes two parts:Propagation of wave section, as well as interaction with unreacted medium. In the propagation process, speed is defined as the constant detonation speed, and detonation products forms instantly. In the interaction process, Mie-Grüneisen mixture model is employed to simulate continuous impact of detonation wave. With Mie-Grüneisen EOS, as well as the Mie-Grüneisen mixture model, motion of detonation wave is simulated well. Comparing with related theoretical data and numerical results good performance was found.

Key words: detonation wave, C-J theory, Mie-Grüneisen equation, detonation speed

中图分类号:

O382

吴宗铎, 严谨, 宗智, 庞建华, 高云. 基于Mie-Grüneisen状态方程的爆轰波运动数值模拟[J]. 计算物理, 2021, 38(1): 47-56.

WU Zongduo, YAN Jin, ZONG Zhi, PANG Jianhua, GAO Yun. Numerical Simulation of Detonation Wave Motion Based on Mie-Grüneisen Equation of State[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2021, 38(1): 47-56.

参考文献

[1] CHAPMAN D L. On the rate of explosion in gases[J]. Philos Mag, 1899, 47:90-104.
[2] JOUGUET E J. On the propagation of chemical reactions in gases[J]. De Mathematiques Pures et Appliqquees, 1905, Appl 1:347-425.
[3] ZEL'DOVICH Y B. On the theory of the propagation of detonation in gaseous systems[J]. Journal of Experimental and Theoretical Physics, 1940, 10:543-568.
[4] Von NEUMANN J. Theory of detonation waves[M]//John von Neumann's collected works. New York:Macmilann, 1942:18-28.
[5] 张宝坪, 张庆民, 黄风雷. 爆轰物理学[M]. 第1版. 北京:兵器工业出版社, 2001.
[6] 王成, SHU C W. 爆炸力学高精度数值模拟研究进展[J]. 科学通报, 2015, 60(10):882-898.
[7] 姜宗林, 滕宏辉, 刘云峰. 气相爆轰物理的若干研究进展[J]. 力学进展,2012, 42(2):129-140.
[8] 滕宏辉, 刘云峰, 姜宗林. 一维脉冲爆轰波的能量释放规律研究[J]. 空气动力学学报, 2012, 30(6):731-736.
[9] TENG Honghui, JIANG Zonglin.Numerical investigation of one-dimensional overdriven detonation initiation[J]. Chinese Journal of Computational Physics, 2008, 25(1):58-64.
[10] QIN Q Y, ZHANG X B. Study on the transition patterns of the oblique detonation wave with varying temperature of the hydrogen-air mixture[J]. Fuel, 2020, 274(16):117827.
[11] NIU Xiao, NI Guoxi, MA Wenhua. Application of adaptive multi-resolution method in numerical simulation of reactive multiphase flows[J]. Chinese Journal of Computational Physics, 2020, 37(6):639-652.
[12] 林佳青, 赖耕, 盛万成. 压差方程Chapman-Jouguet燃烧模型爆轰波与激波的相互作用[J]. 应用数学与计算数学学报, 2017, 31(1):96-106.
[13] LIU Jingjing, ZENG Xianyang, NI Guoxi. Euler numerical methods for reactive flow with general equation of states in two dimensions[J]. Chinese Journal of Computational Physics, 2018, 35(1):29-38.
[14] 许爱国, 张广财, 应阳君. 燃烧系统的Boltzmann建模与模拟研究进展[J]. 物理学报, 2015, 64(18):184701:1-25.
[15] 邵宝力. 基于格子Boltzmann的多孔介质内多场耦合流动与传热模拟[D]. 大庆:东北石油大学, 2019:3-13.
[16] ZHANG Yudong, XU Aiguo, ZHANG Guangcai,et al. A highly-efficient discrete Boltzmann model for detonation[J]. Chinese Journal of Computational Physics, 2016, 33(5):505-515.
[17] LIN Chuandong, LUO K H. Mesoscopic simulation of nonequilibrium detonation with discrete Boltzmann method[J]. Combustion and Flame, 2018, 198:356-362.
[18] ZHANG Yudong, XU Aiguo, ZHANG Guangcai, et al. Kinetic modeling of detonation and effects of negative temperature coefficient[J]. Combustion and Flame, 2016, 173:483-492.
[19] 李诗尧, 于明. 固体炸药爆轰的一种考虑热学非平衡的反应流动模型[J]. 物理学报, 2018, 67(21):214704:1-12.
[20] 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.
[21] LEE E, FINGER M, COLLINS W. JWL equation of states coefficients for high explosives:UCID-16189[R]. USA:Lawrence Livermore Laboratory, University of California, 1973.
[22] 傅华, 李金河, 谭多望,等. 未反应炸药冲击绝热线的近似计算[J]. 高压物理学报,2008, 22(3):329-332.
[23] 王成, 徐文龙, 郭宇飞. 基于基因遗传算法和γ律状态方程的JWL状态方程参数计算[J]. 兵工学报, 2017, 38(suppl):167-173.
[24] SOUERS P C, HASELMAN L C. Detonation equation of state at LLNL:UCRL-ID-116113[R]. CA:Lawrence Livermore National Laboratory, 1994.
[25] ZHAO Yanhong, ZHANG Gongmu, ZHANG Qili, et al. Equation of state of detonation products for HMX explosive[J]. Chinese Journal of Computational Physics, 2017, 34(2):155-159.
[26] 李晓杰, 闫鸿浩, 王小红, 等. Mie-Grüneisen状态方程的物理力学释义[J]. 高压物理学报, 2014, 28(2):227-231.
[27] SHYUE Keh-Ming. A fluid-mixture type algorithm for compressible multicomponent flow with Mie-Grüneisen equation of states[J]. Journal of Computational Physics, 2001, 171(2):678-707.
[28] WU Zongduo, YAN Jin, ZONG Zhi, et al.A multi-component Mie-Grüneisen mixture model based on HLLC algorithm[J].Chinese Journal of Computational Physics, 2020, 37(1):55-62.
[29] WU Zongduo, ZONG Zhi. Numerical calculation of multi-component conservative Euler equations under Mie-Grüneisen equation of state[J]. Chinese Journal of Computational Physics, 2011, 28(6):803-809.
[30] MILLER G H, PUCKETT E G. A high-order Godunov method for multiple condensed phases[J]. Journal of Computational Physics, 1996, 128(1):134-164.
[31] LIU M B, LIU G R, ZONG Z, et al. Computer simulation of high explosive explosion using smoothed particle hydrodynamics methodology[J]. Computers and Fluids, 2003, 32(3):305-322.
[32] 武丹, 刘岩, 王健平. 圆柱形爆轰波的二维数值模拟[J]. 爆炸与冲击, 2015, 35(4):561-566.

基于Mie-Grüneisen状态方程的爆轰波运动数值模拟

Numerical Simulation of Detonation Wave Motion Based on Mie-Grüneisen Equation of State

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	吴宗铎, 严谨, 宗智, 庞建华. 扩散界面形式下一种节约时间步的质量分数混合模型[J]. 计算物理, 2022, 39(5): 510-520.
[2]	吴宗铎, 严谨, 宗智, 赵勇. 一种基于HLLC算法的Mie-Grüneisen多介质混合模型[J]. 计算物理, 2020, 37(1): 55-62.
[3]	刘晶晶, 曾现洋, 倪国喜. 二维多相反应流的欧拉数值方法[J]. 计算物理, 2018, 35(1): 29-38.
[4]	王建航, 陈方, 刘洪, 王兰, 郑忠华. 一种捕捉多向刚性爆轰波的激波定位随机投影方法[J]. 计算物理, 2014, 31(6): 648-658.
[5]	滕宏辉, 杨旸, 姜宗林. 真实比热模型中铝粉尘两相爆轰波的数值研究[J]. 计算物理, 2013, 30(1): 44-52.
[6]	吴宗铎, 宗智. Mie-Grüneisen状态方程下多介质守恒型欧拉方程组的数值计算[J]. 计算物理, 2011, 28(6): 803-809.
[7]	程俊霞. 高能炸药爆速曲率关系的-种计算方法[J]. 计算物理, 2011, 28(6): 817-824.
[8]	程俊霞. 曲率相关的爆轰波在非结构四边形网格上的数值方法[J]. 计算物理, 2011, 28(2): 199-206.
[9]	张磊, 袁礼. 龙格库塔间断有限元方法在计算爆轰问题中的应用[J]. 计算物理, 2010, 27(4): 509-517.
[10]	林文洲, 洪滔, 林忠, 王瑞利. 拉氏程序使用减敏模型模拟炸药绕爆问题[J]. 计算物理, 2010, 27(2): 181-189.
[11]	董刚, 范宝春. 气相爆轰模拟的化学反应加速算法[J]. 计算物理, 2009, 26(1): 27-34.
[12]	刘国昭, 张树道. 气相爆轰高阶中心差分-WENO组合格式自适应网格方法[J]. 计算物理, 2008, 25(4): 387-395.
[13]	滕宏辉, 姜宗林. 一维过驱动爆轰波形成的数值研究[J]. 计算物理, 2008, 25(1): 58-64.
[14]	常利娜, 姜宗林, 秦承森. 化学反应和球面聚心爆轰波的相互作用[J]. 计算物理, 2007, 24(3): 301-306.
[15]	常利娜, 张德良, 胡宗民, 姜宗林. 频散可控格式的一种推广形式及其在爆轰波马赫反射中应用[J]. 计算物理, 2005, 22(3): 189-196.