计算物理 ›› 2010, Vol. 27 ›› Issue (4): 509-517.

• 研究论文 • 上一篇    下一篇

龙格库塔间断有限元方法在计算爆轰问题中的应用

张磊1,2, 袁礼1   

  1. 1. 中国科学院数学与系统科学研究院计算数学所, 科学与工程计算国家重点实验室, 北京 100190;
    2. 中国矿业大学(北京)理学院, 北京 100083
  • 收稿日期:2009-04-21 修回日期:2009-08-02 出版日期:2010-07-25 发布日期:2010-07-25
  • 作者简介:张磊(1977-),男,北京,讲师,博士生,从事问断有限元方法及在流动计算中的应用研究.
  • 基金资助:
    国家973计划(2005CB321703);国家自然科学基金(10531080,10729101)资助项目

Runge-Kutta Discontinuous Galerkin Method for Detonation Waves

ZHANG Lei1,2, YUAN Li1   

  1. 1. LSEC and Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China;
    2. School of Science, China University of Mining and Technology(Beijing), Beijing 100083, China
  • Received:2009-04-21 Revised:2009-08-02 Online:2010-07-25 Published:2010-07-25

摘要: 构造求解带源项守恒律方程组的龙格库塔间断有限元(RKDG)方法,并分别结合源项的Strang分裂法和无分裂法数值求解模型守恒律方程和反应欧拉方程.为了和有限体积型WENO方法进行比较,设计计算源项的WENO重构格式.对一维带源项守恒律的计算表明,对于非刚性问题,RKDG方法比有限体积型WENO方法的误差更小;对于刚性问题,RKDG方法对于间断面位置的捕捉更为精确.对于一二维爆轰波问题的计算结果表明,RKDG方法对爆轰波结构的分辨和爆轰波位置的捕捉能力更强.

关键词: 龙格库塔间断有限元方法, 爆轰波, 反应Euler方程, 刚性源项

Abstract: A Runge-Kutta discontinuous Galerkin(RKDG) method for conservation law with source term is shown.The method is implemented with Strang split or unsplit methods,and is applied to solve one-dimensional conservation law with source term as well as one and two-dimensional detonation wave problems.In order to compare with the fifth-order finite volume WENO method,a special reconstruction method is proposd to calculate integration of the source term with high-order spatial accuracy.Numerical tests in one dimension show that the RKDG method has smaller errors than WENO method for nonstiff problems and is more accurate in capturing position of discontinuity in stiff problems.Numerical simulations of detonation waves demonstrate that the RKDG method is more effcient in resolving detailed structure of detonation waves and location of detonation front.

Key words: Runge-Kutta discontinuous Galerkin method, detonation wave, reactive Euler equations, stiff source term

中图分类号: