航空学报 > 2006, Vol. 27 Issue (2): 198-203

混合解析/数值方法及其在湍流数值模拟上的应用

杜涛1, 吴子牛2, 杨勇1   

  1. 1. 中国运载火箭技术研究院 研究发展中心, 北京 100076;2. 清华大学 工程力学系, 北京 100084
  • 收稿日期:2005-01-01 修回日期:2005-07-02 出版日期:2006-04-25 发布日期:2006-04-25

Mixed Analytical/Numerical Method and Its Application to Problems of Turbulent Flow Simulation

DU Tao1, WU Zi-niu2, YANG Yong1   

  1. 1. China Academy of Launch Vehicle Technology Research & Development Center, Beijing 100076, China;2. Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China
  • Received:2005-01-01 Revised:2005-07-02 Online:2006-04-25 Published:2006-04-25

摘要: 针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。

关键词: 计算流体力学, 混合解析/数值方法, k-ε湍流模型, 源项, 定常态解, 刚性

Abstract: The basic concept of mixed method is to split the original set of equations into a Part Differential Equation (PDE) for the convection-diffusion part, and an Ordinary Differential Equation (ODE) for the source term. The PDE is solved by a conventional numerical method, while the ODE is integrated analytically. It is shown that the mixed method improves the accuracy of source term numerical treatment, and is much more accurate than the traditional method. For the two equations turbulence model, the mixed method without time splitting is developed, which allows for convergence to steady state. The application of mixed method on turbulent flow simulation is studied. Numerical result shows that the mixed analytical/numerical method is much better that the conventional method. The mixed method improves the numerical stability, and accelerates the convergence.

Key words: CFD, mixed analytical/numerical method, k-ε turbulence model, source term, steady state, stiffness

中图分类号: