二维多相反应流的欧拉数值方法

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

计算物理 ›› 2018, Vol. 35 ›› Issue (1): 29-38.DOI: 10.19596/j.cnki.1001-246x.7594

二维多相反应流的欧拉数值方法

刘晶晶, 曾现洋, 倪国喜

北京应用物理与计算数学研究所, 北京 100088

收稿日期:2016-12-01 修回日期:2017-05-16 出版日期:2018-01-25 发布日期:2018-01-25
作者简介:刘晶晶(1990-),女,汉,硕士,主要从事多相反应流数值方法研究,E-mail:liujingjing0618@126.com
基金资助:
挑战计划（TZ2016002）及国家自然科学基金（11171154，11671050和11771050）资助项目

Euler Numerical Methods for Reactive Flow with General Equation of States in Two Dimensions

LIU Jingjing, ZENG Xianyang, NI Guoxi

Institute of Applied Physics and Computational Mathematics, Beijing 100088, China

Received:2016-12-01 Revised:2017-05-16 Online:2018-01-25 Published:2018-01-25

摘要/Abstract

摘要： 本文研究二维多相反应流的数值模拟方法.对于非完全流体，考虑两种状态方程：刚性气体状态方程和JWL（Jones-Wilkins-Lee）型的状态方程.针对这些非完全流体，研究一种混合单元物理量重构法，它由单元内的已知物理量出发，应用单元内不同相的物理量之间的关系：混合密度与各自密度之间的关系、及混合内能与各物质内能之间的关系，再由温度与压强平衡条件，建立关于单变量的非线性方程，对所得方程通过"移动跟踪法"求解，得到单元内的压力、温度及声速等.再利用HLLC解法器，求解各个单元边界数值通量，从而构造一种高效的数值方法.数值结果表明：该算法能清晰地捕捉爆轰波的结构特征，胞格排列有序，三波点的特征明显.

关键词: 多相反应流, 爆轰波, Godunov格式, HLLC解法器, 移动跟踪法

Abstract: We present an efficient method to simulate reactive flow for general equation of states in two dimensions. Two kinds of nonideal equation of states for compressible fluid model coupling with reactive rate equation are concerned. The important aspect is to deal with mixture of different phase in one cell, which inevitably happens in Euler method for reactive flows. Physical variables such as pressure,velocity and speed of sound in each cell are reconstructed which results in nonlinear algebra equations. They are used to obtain flux by HLLC Riemann solver. Numerical examples of detonation in two dimensions with different equation of states demonstrate accuracy and robustness of the method.

Key words: reactive flow, detonation wave, HLLC solver, Godunov scheme, general equation of states

中图分类号:

刘晶晶, 曾现洋, 倪国喜. 二维多相反应流的欧拉数值方法[J]. 计算物理, 2018, 35(1): 29-38.

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.

参考文献

[1] DAVIS S F. Simplified second order Godunov type methods[J]. SIAM Journal of Scientific and Statistical Computing, 1988, 9:445-473.
[2] DORING W. Uberder detonation vergang in gasen[J]. Ann Phys, 1943,43(5):421-436.
[3] MARTELOT S. NKONGA B, SAUREL R. Liquid and liquid-gas flows at all speeds[J].J Comput Phys, 2013,255:53-82.
[4] PAPALEXANDRIS M V, LEONARD A, DIMOTAKIS P E. Unsplit schemes for hyperbolic conservation laws with source terms in one space dimension[J]. J Comput Phys, 1997,134(1):31-61.
[5] WESCOTT B L, STEWART D S, DAVIS W C.Equation of state and reaction rate for condensed phase explosives[J]. Journal of Applied Physics, 2005,98:053514.
[6] ZELDOVICH Ya B. On the theory of the propagation of detonation in gaseous systems[J]. Sov Phys JETP, 1940,10(5):542-568.
[7] BOURLIOUX A, MAJDA A, VICTOR ROYTBURD. Theoretical and numerical structure for unstable one-dimensional detonations[J]. Siam J Appl Math, 1991,51(2):303-343.
[8] TRUESDELL C. Rational thermodynamics[M]. New York:McGraw-Hill, 1984.
[9] TORO E F. Riemann solvers and numerical methods for fluid dynamics:A practical introduction[M]. Spronger Berlin Heidelberg, 2009.
[10] TORO E F, CASTRO C, LEE B J. A novel numerical flux for the 3D Euler equations with general equation of state[J]. J Compute Phys, 2015, 303:80-94.
[11] 李俊红,张亮,俞继军,等.高超声速可压缩流中粗糙壁热流研究[J].计算物理,2017,34(2):165-174.
[12] 蒋华,董刚,陈霄.激波与火焰面相互作用数值模拟的GPU加速[J].计算物理,2016,33(1):23-29.
[13] AZARENOK B N, TANG T. Second-order Godunov-type scheme for reactive flow calculations on moving meshes[J], J Comput Phys, 2005,206:48-80.
[14] BEAR M R, NUNZITO J W. A two-phase mixture theory for the deflagration-to-detonation transition(DDT) in reactive granular materials[J]. Int J Multiphase Flow, 1986,12(6):861-889.
[15] WANG C, JIANG Z, GAO Y. Half-cell law of regular cellular detonation[J]. Chinese Physics Letters, 2008, 25(10):3704-3707.
[16] FICKETT W, DAVIS W C. Detonation:Theory and experiment[M]. Dover Publications Inc, Mineola, NY, 1979.
[17] FICKETT W, WOOD W C. Flow calculations for pulsating one-dimensional detonations[J]. Phys Fluid, 1966,9(5):903-916.
[18] SAUREL R, PETIPAS F, BERRY R. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures[J]. Journal of Computational Physics, 2009, 228:1678-1712.
[19] NI G X, JIANG S, WANG S H. A remapping-free, efficient Riemann-solvers based, ALE method for multi-material fluids with general EOS[J]. Computers and Fluids, 2013,71:19-27.
[20] TORO E F. Riemann solvers and numerical methods for fluid dynamics[M]. 3rd ed.Springer-Verlag Berlin Heidelberg, 2009.
[21] ZENG X Y, XIAO M, NI G X.An efficient numerical method for reactive flow with general equation of states (one dimensional)[J]. Int J Numer Meth Fluids, DOI:10.1002/fld.4233.
[22] 曾现洋.多相反应流的数值模拟方法[D].中国工程物理研究院学位论文,2016.
[23] 曾现洋,倪国喜.振动NACA0012翼型的移动网格数值模拟[J].计算物理,2016,33(3):266-272.
[24] VON NEUMANN J. Theory of detonation waves[R].Institute for advanced study Princeton N J, National Defense Research, 1942.

二维多相反应流的欧拉数值方法

Euler Numerical Methods for Reactive Flow with General Equation of States in Two Dimensions

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	吴宗铎, 严谨, 宗智, 庞建华, 高云. 基于Mie-Grüneisen状态方程的爆轰波运动数值模拟[J]. 计算物理, 2021, 38(1): 47-56.
[2]	戢美璇, 蔚喜军, 宋明阳. 一维含化学反应流体力学方程组的ALE间断有限元方法[J]. 计算物理, 2020, 37(1): 1-9.
[3]	王建航, 陈方, 刘洪, 王兰, 郑忠华. 一种捕捉多向刚性爆轰波的激波定位随机投影方法[J]. 计算物理, 2014, 31(6): 648-658.
[4]	滕宏辉, 杨旸, 姜宗林. 真实比热模型中铝粉尘两相爆轰波的数值研究[J]. 计算物理, 2013, 30(1): 44-52.
[5]	程俊霞. 曲率相关的爆轰波在非结构四边形网格上的数值方法[J]. 计算物理, 2011, 28(2): 199-206.
[6]	张磊, 袁礼. 龙格库塔间断有限元方法在计算爆轰问题中的应用[J]. 计算物理, 2010, 27(4): 509-517.
[7]	林文洲, 洪滔, 林忠, 王瑞利. 拉氏程序使用减敏模型模拟炸药绕爆问题[J]. 计算物理, 2010, 27(2): 181-189.
[8]	董刚, 范宝春. 气相爆轰模拟的化学反应加速算法[J]. 计算物理, 2009, 26(1): 27-34.
[9]	刘国昭, 张树道. 气相爆轰高阶中心差分-WENO组合格式自适应网格方法[J]. 计算物理, 2008, 25(4): 387-395.
[10]	滕宏辉, 姜宗林. 一维过驱动爆轰波形成的数值研究[J]. 计算物理, 2008, 25(1): 58-64.
[11]	郑建国, 马东军, 孙德军, 尹协远. Van der Waals状态方程可压缩多流体流动的PPM目的[J]. 计算物理, 2007, 24(3): 287-294.
[12]	常利娜, 姜宗林, 秦承森. 化学反应和球面聚心爆轰波的相互作用[J]. 计算物理, 2007, 24(3): 301-306.
[13]	吴昊, 沈智军. Burgers方程的声速点故障问题[J]. 计算物理, 2006, 23(4): 396-404,440.
[14]	常利娜, 张德良, 胡宗民, 姜宗林. 频散可控格式的一种推广形式及其在爆轰波马赫反射中应用[J]. 计算物理, 2005, 22(3): 189-196.
[15]	胡宗民, 牟乾辉, 张德良, 姜宗林. 爆轰波在弯管内传播过程数值分析[J]. 计算物理, 2004, 21(5): 408-414.