自适应多分辨率方法在反应多相流数值模拟中的应用

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

计算物理 ›› 2020, Vol. 37 ›› Issue (6): 639-652.DOI: 10.19596/j.cnki.1001-246x.8176

自适应多分辨率方法在反应多相流数值模拟中的应用

牛霄¹, 倪国喜¹, 马文铧²

1. 北京应用物理与计算数学研究所, 北京 100088;
2. 上海交通大学数学科学学院, 上海 200240

收稿日期:2019-11-22 修回日期:2020-03-11 出版日期:2020-11-25 发布日期:2020-11-25
作者简介:牛霄(1991-),男,博士研究生,主要从事多相反应流数值方法研究,E-mail:niuxiao17@gscaep.ac.cn
基金资助:
挑战计划（TZ2016002）资助项目

Application of Adaptive Multi-resolution Method in Numerical Simulation of Reactive Multiphase Flows

NIU Xiao¹, NI Guoxi¹, MA Wenhua²

1. Institute of Applied Physics and Computational Mathematics, Beijing 100088, China;
2. School of Mathematical Sciences, Shanghai Jiao Tong University, Shanghai 200240, China

Received:2019-11-22 Revised:2020-03-11 Online:2020-11-25 Published:2020-11-25

摘要/Abstract

摘要： 将自适应多分辨率方法应用到刚性气体状态方程CJ模型的数值模拟，采用多相流问题的守恒锐利界面格式，通过level-set方法和虚拟流体方法来追踪和处理界面，能够很好地处理时间尺度较大的界面交互问题.利用金字塔型数据结构和多分辨率自适应方法，提高算法的存储效率和计算效率.给出一维和二维的数值算例，证明该算法在反应多相流数值模拟中的稳定性和高效性.

关键词: 反应多相流, 自适应多分辨率方法, 刚性气体状态方程, CJ模型

Abstract: We applie adaptive multi-resolution method to numerical simulation of stiffened-gas equation of state based on CJ model. A conservation sharp interface method of multiphase flow, in which the interface problem is tracked and solved by level-set method and ghost fluid method, is adopted. This method handles long physical time scale interactions across interface well and reduces conservative associated errors. In addition, pyramid data structure and adaptive multi-resolution analysis are used to improve storage efficiency and calculation efficiency of the algorithm. Finally, numerical examples are given to demonstrate stability and effectiveness of the adaptive multi-resolution method in numerical simulation of reactive multiphase flows.

Key words: reactive multiphase flow, adaptive multi-resolution method, stiffened-gas equation of state, CJ model

中图分类号:

O381

牛霄, 倪国喜, 马文铧. 自适应多分辨率方法在反应多相流数值模拟中的应用[J]. 计算物理, 2020, 37(6): 639-652.

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.

参考文献

[1] CHAPMAN D L. On the rate of explosion in gases[J]. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1899, 47(284):90-104.
[2] JOUGUET E. On the propagation of chemical reactions in gases[J]. Journal De Mathematiques Pures Et Appliquees, 1905, 1:347-425, Continued in 1906, 2:5-85.
[3] TRYGGVASON G, BUNNER B, ESMAEELI A. A front-tracking method for the computations of multiphase flow[J]. Journal of Computational Physics, 2001, 169(2):708-759.
[4] GLIMM J, LI X, XU Z, ZHAO N. Conservative front tracking with improved accuracy[J]. SIAM Journal on Numerical Analysis, 2003, 41(5):1926-1947.
[5] TAI Y C, NOELLE S, GRAY J M N T. Shock-capturing and front tracking methods for granular avalanches[J]. Journal of Computational Physics, 2002, 175(2):269-301.
[6] SHIRANI E, ASHGRIZ N, MOSTAGIMI J. Interface pressure calcuation based on conservarion of momentum for front capturing methods[J]. Journal of Computational Physics, 203(1), 2005:154-175.
[7] HIRT C W, NICHOLS B D. Volume of fluid method for the dynamics of free boundaries[J]. Journal of Computational Physics, 1981, 39(1):201-225.
[8] YABE T, XIAO F, UTSUMI T. The constrained interpolation profile method for multiphase analysis[J]. Journal of Computational Physics, 2001, 169(2):556-593.
[9] ADALSTEINSSON D, SETHIAN J A. The fast construction of extension velocites in level set methods[J]. Journal of Computational Physics, 2002, 175(2):269-301.
[10] OSHER S, FEDKIW R, PIECHOR K. Level set methods:An overview and some recent results[J]. Journal of Computational Physics, 2001, 169(2):463-502.
[11] YUAN B. Level set parallel highly accurate evolution based on GPU cluster[J]. Chinese Journal of Computational Physics, 2016, 33(3):253-265.
[12] DING S R, XU S L, LU J F, ZHANG M P. Moving shock interacting with cylindrical water columns and a nozzle flow field:rGFM and level-set method[J]. Chinese Journal of Computational Physics, 2019, 36(2):165-174.
[13] FEDKIW R, ASLAM T, MERRIMAN B, OSHER S. A non-oscillatory Eulerian approach to interfaces in multimaterial flows[J]. Journal of Computational Physics, 1999, 152(2):457-492.
[14] HU X Y, KHOO B C, ADAMS N A, HUANG F L. A conservative interface method for compressible flows[J]. Journal of Computational Physics, 2006, 219(2):553-578.
[15] BERGER M J, OLIGER J. Adaptive mesh refinement for hyperbolic partial differential equations[J]. Journal of Computational Physics, 1984, 53(3):484-512.
[16] MINIATI F, COLELLA P. Block structured adaptive mesh and time refinement for hybrid, hyperbolic+N-body systems[J]. Journal of Computational Physics, 2007, 227:400-430.
[17] MACNEICE P, OLSON K M, MOBARRY C. PARAMESH:A parallel adaptive mesh refinement community toolkit[J]. Journal of Computational Physics, 2000, 126:330-354.
[18] HARTEN A. Adaptive multiresolution schemes for shock computations[J]. Journal of Computational Physics, 1994, 115:319-338.
[19] ROUSSEL O, SCHNEIDER K, TSIGULIN A. A conservative fully adaptive multiresolution algorithm for parabolic PDEs[J]. Journal of Computational Physics, 2003, 188:493-523.
[20] DOMINGUES M O, GOMES S M, ROUSSEL O. Extended generalized Lagrangian multipliers for magnetohydrodynamics using adaptive multiresolution methods[J].ESAIM:Proceedings, 2013, 34:95-107.
[21] GONG X F, YANG J M, ZHANG S D. A parallel SPH method with background grid of adaptive mesh refinement[J]. Chinese Journal of Computational Physics, 2016, 33(2):183-189.
[22] GAO X, GROTH C P T. A parallel adaptive mesh refinement algorithm for predicting turbulent non-premixed combusting ows[J]. International Journal of Computational Fluid Dynamics, 2006, 20(5):349-357.
[23] BIHARI B L, HARTEN A. Multiresolution schemes for the numerical solution of 2-D conservation laws I[J]. Siam Journal on Scientific Computing, 2006, 18(2):315-354.
[24] BIHARI B L, SCHWENDEMAN D. Multiresolution schemes for the reactive Euler equations[J]. Journal of Computational Physics, 1999, 154(1):197-230.
[25] HEJAZIALHOSSEINI B, ROSSINELLI D, BERGDORF M. High order finite volume methods on wavelet-adapted grids with local time-stepping on multicore architectures for the simulation of shock-bubble interactions[J]. Journal of Computational Physics, 2010, 229(22):8364-8383.
[26] HAN L H, HU X Y, ADAMS A. Adaptive multi-resolution method for compressible multi-phase ows with sharp interface model and pyramid data structure[J]. Journal of Computational Physics, 2014, 262:131-152.
[27] PAN S C, HU X Y, ADAMS A. High-resolution method for evolving complex interface networks[J]. Computer Physics Communications. 2018, 225:10-27.
[28] HARLOW F H, AMSDEN A. Numerical calculation of almost incompressible flow[J]. Journal of Computational Physics, 1968, 3:80-93.
[29] KAPILA A K, SCHWENDEMAN D W, BDZIL J B. A study of detonation diffraction in the ignition-and-growth model[J]. Combustion Theory and Modelling, 2007, 11(5):781-822.
[30] TENG Z H, CHORIN A J, LIU T P. Riemann problems for reacting gas with applications to transition[J]. SIAM Journal on Numerical Analysis, 1982, 42(5):946-981.
[31] DOMINGUES M O, GOMES S M, ROUSSEL O, et al. An adaptive multiresolution scheme with local time stepping for evolutionary PDEs[J]. Journal of Computational Physics, 2008, 227(8):3758-3780.
[32] FEDKIW R P, ASLAM T, XU S J. The ghost fluid method for deflagration and detonation discontinuities[J]. Journal of Computational Physics, 1999, 154:393-427.

自适应多分辨率方法在反应多相流数值模拟中的应用

Application of Adaptive Multi-resolution Method in Numerical Simulation of Reactive Multiphase Flows

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 14

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	肖敏, 徐喜华, 倪国喜. 基于任意拉格朗日—欧拉型移动网格的反应流广义黎曼问题方法[J]. 计算物理, 2020, 37(2): 127-139.
[2]	岳军政, 董贺飞, 洪滔. 拐角空间中悬浮铝粉尘爆轰后效的三维数值模拟[J]. 计算物理, 2020, 37(1): 37-45.
[3]	李诗尧, 于明. 一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法[J]. 计算物理, 2019, 36(5): 505-516.
[4]	赵艳红, 张弓木, 张其黎, 李琼. HMX炸药爆轰产物物态方程研究[J]. 计算物理, 2017, 34(2): 155-159.
[5]	滕宏辉, 杨旸, 姜宗林. 真实比热模型中铝粉尘两相爆轰波的数值研究[J]. 计算物理, 2013, 30(1): 44-52.
[6]	林文洲, 洪滔, 林忠, 王瑞利. 拉氏程序使用减敏模型模拟炸药绕爆问题[J]. 计算物理, 2010, 27(2): 181-189.
[7]	张德良, 王景焘, 王刚. 高阶精度CE/SE算法及其应用[J]. 计算物理, 2009, 26(2): 211-220.
[8]	刘国昭, 张树道. 气相爆轰高阶中心差分-WENO组合格式自适应网格方法[J]. 计算物理, 2008, 25(4): 387-395.
[9]	童慧峰, 唐志平, 张凌. 由强激光产生的烧蚀等离子体随时间演化的数值模拟[J]. 计算物理, 2007, 24(6): 667-672.
[10]	郑建国, 马东军, 孙德军, 尹协远. Van der Waals状态方程可压缩多流体流动的PPM目的[J]. 计算物理, 2007, 24(3): 287-294.
[11]	常利娜, 姜宗林, 秦承森. 化学反应和球面聚心爆轰波的相互作用[J]. 计算物理, 2007, 24(3): 301-306.
[12]	王昌建, 徐胜利, 朱建士. 气相爆轰波在分叉管中传播现象的数值研究[J]. 计算物理, 2006, 23(3): 317-324.
[13]	楼建锋, 于恒. CTVD格式数值计算非均质炸药爆轰问题[J]. 计算物理, 2005, 22(4): 358-364.
[14]	田宙, 郭永辉, 郝保田. 气液爆轰波及冲击波传播的数值研究[J]. 计算物理, 2000, 17(S1): 131-136.