宇航学报 ›› 2016, Vol. 37 ›› Issue (9): 1135-1141.doi: 10.3873/j.issn.1000-1328.2016.09.014

• 环境试验与器件 • 上一篇    下一篇

增强高超声速化学反应流数值计算的稳定性方法

王羽,蔡晋生,屈崑,粟虹敏   

  1. 西北工业大学航空学院,西安 710072
  • 收稿日期:2016-03-18 修回日期:2016-06-12 出版日期:2016-09-15 发布日期:2016-09-25
  • 基金资助:

    国家自然科学基金(11472221)

A Method for Improving the Stability of Numerical Calculation in Hypersonic Chemically Reacting Flows

WANG Yu, CAI Jin sheng, QU Kun, SU Hong min   

  1. School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China
  • Received:2016-03-18 Revised:2016-06-12 Online:2016-09-15 Published:2016-09-25

摘要:

针对高超声速化学反应流数值模拟中包含源项对角化雅可比矩阵的数值格式稳定性很差的问题,在对采用上下对称的高斯-赛德尔(LU-SGS)格式的离散方程进行线性稳定性分析的基础上,提出一种依据流动状态对流场的块网格进行局部加密的方法,以增强流场数值计算稳定性。在高超声速流场中,对包含激波后亚声速高温区域的块网格给予加密,而其他区域中的块网格可以适当进行网格粗化或保持不变。算例结果表明,相比于初始网格,局部加密后的改进网格并不会产生过多的网格量,而且能够采用较大的柯朗-弗里德里奇-列维(CFL)数以达到数值求解过程中稳定计算和快速收敛的目的。

关键词: 高超声速, 对角化, 稳定性分析, 网格加密

Abstract:

In numerical simulations of hypersonic flows with chemical reactions, the stability of the numerical scheme including the diagonalized Jacobian of chemical source becomes terrible. In this paper, according to the linear stability analysis of the discrete equation with the lower-upper symmetric Gauss-Seidel (LU-SGS) scheme, a local block grid refinement method is suggested based on the flow conditions for improving the stability. In the hypersonic chemically reacting flows, the refinement should be carried out for the block grid in the areas of subsonic state and high temperature behind the shock wave, and the block grid in other areas may be coarsened or invariable. The numerical results demonstrate that the improved grid has a stable computation and rapid convergence with a larger Courant-Friedrichs-Lewy (CFL) number, but cannot produce overmuch grid cells compared with initial grid.

Key words: Hypersonic, Diagonalization, Stability analysis, Grid refinement

中图分类号: