北京航空航天大学学报 ›› 2004, Vol. 30 ›› Issue (07): 679-681.

• 论文 • 上一篇    下一篇

高温化学非平衡计算中源项处理方法

马辉, 王发民   

  1. 中国科学院力学研究所, 北京 100080
  • 收稿日期:2003-02-28 出版日期:2004-07-31 发布日期:2010-09-21
  • 作者简介:马 辉(1973-),女,黑龙江齐齐哈尔人,博士生, ma hui@iapcm.ac.cn.

Calculation method of source terms in high temperature chemical nonequilibrium flows

Ma Hui, Wang Famin   

  1. Institute of Mechanics, Chinese Academy of Sciences, Beijing 100080, China
  • Received:2003-02-28 Online:2004-07-31 Published:2010-09-21

摘要: 高超声速三维化学非平衡绕流流场的数值模拟中,在计算初期容易发生组元密度出现负值的非物理现象,另外源项的刚性是影响计算稳定性和收敛速度的主要原因。根据源项线性化理论,采用2种源项线性化方法处理化学反应源项.一种为根据流场内化学反应物理规律构造的线性化方法,该方法能抑制在计算过程中组元密度出现负值现象,提高了计算稳定性,加速了收敛速度.同时采用了时间预处理矩阵的线性化方法,较好地解决了非平衡化学反应与流场耦合的刚性问题。证明了所构造源项线性化方法相容性.数值实验表明该方法有效地避免了源项计算中密度出负的问题,加快了计算的收敛速度,从而提高了计算效率.

Abstract: In the numerical simulation of hypersonic chemical nonequilibrium flows, species density may become negative at the beginning of calculation, chemical stiffness is the main reason that affects the calculation’s stability and convergence. Based on the basic theory of source term’s linearazition, two linearazition methods were chosen. One is based on the physical rules of chemical reaction flows. The method’s advantage is that can improve the calculation’s stability and can avoid the species’ density to become negative during the calculation. The other method is the time preconditioner method, which has the function of diminishing the stiffness of the source term. The method was proved to be well approximates the partial differential equations. The calculation result shows that this method can effectively avoid negative density’s appearance and accelerate the calculation’s convergence. So the method effectively improves the calculation’s efficiency.

中图分类号: 


版权所有 © 《北京航空航天大学学报》编辑部
通讯地址:北京市海淀区学院路37号 北京航空航天大学学报编辑部 邮编:100191 E-mail:jbuaa@buaa.edu.cn
本系统由北京玛格泰克科技发展有限公司设计开发