一个高效的离散Boltzmann爆轰模型

计算物理 ›› 2016, Vol. 33 ›› Issue (5): 505-515.

• 研究论文 • 下一篇

一个高效的离散Boltzmann爆轰模型

张玉东¹, 许爱国^2,3, 张广财², 祝成民¹

1. 北京航空航天大学宇航学院, 北京 100191;
2. 北京应用物理与计算数学研究所计算物理重点实验室, 北京 100088;
3. 北京大学应用物理与技术研究中心和高能量密度物理数值模拟教育部重点实验室, 北京 100871

收稿日期:2015-07-11 修回日期:2015-12-09 出版日期:2016-09-25 发布日期:2016-09-25
通讯作者: 许爱国(1970-),男,博士,教授,主要研究复杂物理系统的建模、模拟与分析,E-mail:Xu_Aiguo@iapcm.ac.cn;祝成民(1982-),男,博士,副教授,主要研究火箭发动机新技术、计算流体力学等,E-mail:cmzhu@buaa.edu.cn
作者简介:张玉东(1992-),男,硕士,主要研究燃烧问题的离散Boltzmann建模与模拟,E-mail:zhangyudongzzu@163.com
基金资助:
国家自然科学基金（11475028）和计算物理重点实验室基金资助项目

A Highly-efficient Discrete Boltzmann Model for Detonation

ZHANG Yudong¹, XU Aiguo^2,3, ZHANG Guangcai², ZHU Chengmin¹

1. School of Astronautics, Beijing University of Aeronautics and Astronautics, Beijing 100191, China;
2. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China;
3. Center for Applied Physics and Technology, MOE Key Center for High Energy Density Physics Simulations, College of Engineering, Peking University, Beijing 100871, China

Received:2015-07-11 Revised:2015-12-09 Online:2016-09-25 Published:2016-09-25

摘要/Abstract

摘要： 构建一个适用于爆轰过程模拟的离散Boltzmann模型.该模型由一个离散Boltzmann方程和一个唯象反应率方程构成；在物理建模上，它等效于一个传统Navier-Stokes模型外加一个关于热动非平衡行为的粗粒化模型.与传统流体模型相比，它能够提供更多的动力学和动理学信息.该模型采用16个离散速度，相比于使用33个离散速度的模型具有更高的运算效率，模型中引入了额外自由度，通过调节额外自由度的数目，可以模拟各种不同比热比的爆轰.采用爆轰问题中的一些经典算例对所建立的模型进行数值验证.结果表明：该模型不仅对传统流体模型所能模拟的爆轰问题有效，而且能够用于一些传统流体模型不能描述的非平衡过程，有利于对爆轰问题的深入研究.

关键词: 爆轰, 离散Boltzmann模型, 唯象反应率模型, 非平衡信息

Abstract: A discrete Boltzmann model is presented for detonation simulation. It consists of a discrete Boltzmann equation and a phenomenological reaction rate equation. In physical modeling, it is equivalent to a traditional Navier-Stokes model and a coarse grained model for thermodynamic non-equilibrium behavior. Compared with traditional fluid model, it provides more dynamic information and more kinetic information as well. The model adopts 16 discrete velocities, which has high computational efficiency. Additional degrees of freedom are introduced so that it can be used to simulate various detonations with different specific heat ratios. Several classic examples of detonation are simulated to validate the model. It shows that besides detonation dynamics where traditional hydrodynamic model works, it can also be used to study non-equilibrium effects where traditional fluid model does not work.

Key words: detonation, discrete Boltzmann model, phenomenological reaction rate model, non-equilibrium information

中图分类号:

O38
O354.7

张玉东, 许爱国, 张广财, 祝成民. 一个高效的离散Boltzmann爆轰模型[J]. 计算物理, 2016, 33(5): 505-515.

ZHANG Yudong, XU Aiguo, ZHANG Guangcai, ZHU Chengmin. A Highly-efficient Discrete Boltzmann Model for Detonation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(5): 505-515.

参考文献

[1] 姜宗林,滕宏辉,刘云峰.气相爆轰物理的若干研究进展[J].力学进展,2012,42(2):129-140.
[2] 孙锦山,朱建士.理论爆轰物理[M]. 北京:国防工业出版社,1995:88-172,386-395.
[3] JU Y G. Recent progress and challenges in fundamental combustion research[J].Advances in Mechanics,2014, 44(1):1-72.
[4] 吴望一.流体力学[M]. 北京:北京大学出版社, 2006:156-161.
[5] 许爱国, 张广财, 李英骏,等. 非平衡与多相复杂系统模拟研究——Lattice Boltzmann动理学理论与应用[J]. 物理学进展, 2014, 34(3):136-167.
[6] XU A G, ZHANG G C, GAN Y B, et al, Lattice Boltzmann modeling and simulation of compressible flows[J].Front Phys, 2012,7(5):582-600.
[7] SUCCI S. The lattice Boltzmann equation for fluid dynamics and beyond[M]. Oxford:Clarendon Press, 2001: 77-176.
[8] GAN Y B, XU A G, ZHANG G C,et al. Lattice BGK kinetic model for high-speed compressible flows: Hydrodynamic and nonequilibrium behaviors[J]. Europhysics Letter, 2013, 103 (24003): 1-6.
[9] LIN C D, XU A G, ZHANG G C, et al. Polar-coordinate lattice Boltzmann modeling of compressible flows[J]. Physical Review E: Statistical Nonlinear & Soft Matter Physics,2014, 89(1):390-400.
[10] CHEN F, XU A G, ZHANG G C, et al. Two-dimensional multiple-relaxation-time lattice Boltzmann model for compressible and incompressible flows[J]. Front Phys, 2014,9(2):246-254.
[11] YAN B, XU A G, ZHANG G C, et al. Lattice Boltzmann model for combustion and detonation[J].Front Phys, 2013,8(1):94-110.
[12] LIN C D, XU A G, ZHANG G C, et al. Polar coordinate lattice Boltzmann kinetic modeling of detonation phenomena[J]. Commun Theor Phys, 2014,62(5):737-748.
[13] XU A G, LIN C D, ZHANG G C, et al. Multiple-relaxation-time lattice Boltzmann kinetic model for combustion[J]. Physical Review E: Statistical Nonlinear & Soft Matter Physics, 2014, 91(4):003300.
[14] 许爱国,张广财,应阳君. 燃烧系统的离散Boltzmann建模与模拟研究进展[J]. 物理学报,2015, 64(18) :184701.
[15] XU A G, ZHANG G C, GAN Y B. Discrete Boltzmann modeling of liquid-vapor system[J]. Eprint Arxiv, 2014.
[16] GAN Y B, XU A G, ZHANG G C,et al. Discrete Boltzmann modeling of multiphase flows: hydrodynamic and thermodynamic non-equilibrium effects[J]. Soft Matter, 2015, 11.
[17] 何雅玲,王勇,李庆.格子Boltzmann方法的理论及应用[M]. 北京:科学出版社,2011:48-49.
[18] MINORU W, MICHISA T. Two-dimensional thermal model of the finite-difference lattice Boltzmann method with high spatial isotropy.[J]. Physical Review E: Statistical Nonlinear & Soft Matter Physics, 2003, 67(3Pt2):210-215.
[19] 水鸿寿.一维流体力学差分方法[M].北京:国防工业出版社,1989:500-509.
[20] 张宝平,张庆明,黄风雷.爆轰物理学[M].北京:兵器工业出版社,2001:64-79.

一个高效的离散Boltzmann爆轰模型

A Highly-efficient Discrete Boltzmann Model for Detonation

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	吴宗铎, 严谨, 宗智, 庞建华, 高云. 基于Mie-Grüneisen状态方程的爆轰波运动数值模拟[J]. 计算物理, 2021, 38(1): 47-56.
[2]	岳军政, 董贺飞, 洪滔. 拐角空间中悬浮铝粉尘爆轰后效的三维数值模拟[J]. 计算物理, 2020, 37(1): 37-45.
[3]	李诗尧, 于明. 一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法[J]. 计算物理, 2019, 36(5): 505-516.
[4]	刘晶晶, 曾现洋, 倪国喜. 二维多相反应流的欧拉数值方法[J]. 计算物理, 2018, 35(1): 29-38.
[5]	赵艳红, 张弓木, 张其黎, 李琼. HMX炸药爆轰产物物态方程研究[J]. 计算物理, 2017, 34(2): 155-159.
[6]	董贺飞, 洪滔, 张晓立. 弱点火条件下炸药燃烧转爆轰的数值模拟[J]. 计算物理, 2016, 33(1): 15-22.
[7]	王建航, 陈方, 刘洪, 王兰, 郑忠华. 一种捕捉多向刚性爆轰波的激波定位随机投影方法[J]. 计算物理, 2014, 31(6): 648-658.
[8]	滕宏辉, 杨旸, 姜宗林. 真实比热模型中铝粉尘两相爆轰波的数值研究[J]. 计算物理, 2013, 30(1): 44-52.
[9]	董贺飞, 洪滔, 张晓立. CE/SE方法数值模拟炸药粉尘爆轰[J]. 计算物理, 2012, 29(4): 495-502.
[10]	程俊霞. 曲率相关的爆轰波在非结构四边形网格上的数值方法[J]. 计算物理, 2011, 28(2): 199-206.
[11]	张磊, 袁礼. 龙格库塔间断有限元方法在计算爆轰问题中的应用[J]. 计算物理, 2010, 27(4): 509-517.
[12]	林文洲, 洪滔, 林忠, 王瑞利. 拉氏程序使用减敏模型模拟炸药绕爆问题[J]. 计算物理, 2010, 27(2): 181-189.
[13]	董刚, 范宝春. 气相爆轰模拟的化学反应加速算法[J]. 计算物理, 2009, 26(1): 27-34.
[14]	刘国昭, 张树道. 气相爆轰高阶中心差分-WENO组合格式自适应网格方法[J]. 计算物理, 2008, 25(4): 387-395.
[15]	滕宏辉, 姜宗林. 一维过驱动爆轰波形成的数值研究[J]. 计算物理, 2008, 25(1): 58-64.