计算物理 ›› 2017, Vol. 34 ›› Issue (6): 631-640.DOI: 10.19596/j.cnki.1001-246x.7568

• •    下一篇

DG方法求解可压缩气固两相流动

段茂昌1,2, 蔚喜军2, 陈大伟2, 黄朝宝1,2, 安娜3   

  1. 1. 中国工程物理研究院研究生院, 北京 100088;
    2. 北京应用物理与计算数学研究所, 计算物理实验室, 北京 100088;
    3. 山东师范大学数学与统计学院, 山东 济南 250014
  • 收稿日期:2016-10-28 修回日期:2017-03-08 出版日期:2017-11-25 发布日期:2017-11-25
  • 通讯作者: 蔚喜军,E-mail:yuxj@iapcm.ac.cn
  • 作者简介:段茂昌(1988-),男,博士研究生,主要从事计算流体力学和两相流研究,E-mail:duanmaochang@163.com
  • 基金资助:
    国家自然科学基金(11571002)、中国工程物理研究院科学基金(2013A0202011,2015B0101021)和国防基础科研计划(B1520133015)资助项目

DG Method for Compressible Gas-Solid Two-phase Flow

DUAN Maochang1,2, YU Xijun2, CHEN Dawei2, HUANG Chaobao1,2, AN Na3   

  1. 1. Graduate School of China Academy of Engineering Physics, Beijing 100088, China;
    2. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China;
    3. School of Mathematics and Statistics, Shangdong Normal University, Jinan 250014, China
  • Received:2016-10-28 Revised:2017-03-08 Online:2017-11-25 Published:2017-11-25

摘要: 给出求解双向耦合可压缩气固两相流的间断有限元方法,对所得到的气固两相流方程组不需要采用分裂的方法离散,对气相、颗粒相方程及其对流部分和源项可以统一处理,两相都采用基于近似Riemann解的数值通量.数值模拟低压含尘激波管内的两相非平衡流动,并与平衡流、冻结流的结果进行比较.分析颗粒相的存在对气体运动的影响,及激波后松弛区域内两相间相互作用规律.发现颗粒质量比决定两相平衡后的最终状态,而颗粒直径决定两相流从非平衡到平衡的过渡过程,即不同尺寸颗粒对应的驰豫时间、松弛距离不同.结果表明:本文提出的计算方法对求解可压缩气固两相流是可行的,为研究复杂的气固两相流动问题奠定了基础.

关键词: 间断有限元, 稀疏气粒两相流, 可压缩气固两相流, 含尘气体激波管, 松弛效应

Abstract: We present DG method for solving two-way coupling compressible gas-solid two-phase flow. Equations of both phases are discreted simultaneously, including convection term and source term. Splitting technique to discretize governing equations separately is avoided. Numerical flux of both phases is based on approximate Riemann solver. Dusty-gas shock tube problem with particles in low pressure section is simulated. Comparisons of equilibrium flow and frozen flow are made. Influence of particles in gas and interaction rules between two phases in relaxation zone behind shock are studied. It found that mass ratio of particles determines last equilibrium state and particles diameter determines transition process of two-phase flow from nonequilibrium to equilibrium flow. Namely, different diameter particles correspond to different relaxation time and distance. It shows that the numerical method proposed is reliable. It lays a foundation for more complicated gas-solid two-phase flow problems.

Key words: DG method, dilute gas particle flow, compressible gas solid two phase flow, dusty gas shock tube, relaxation effect

中图分类号: