航空学报 > 2003, Vol. 24 Issue (2): 111-115

混合解析/数值方法在含高频振荡点源项的欧拉方程数值计算上的应用

杜涛   

  1. 清华大学工程力学系 北京100084
  • 收稿日期:2001-10-23 修回日期:2002-10-09 出版日期:2003-04-25 发布日期:2003-04-25

Mixed Analytical/Numerical Method Applied to Euler Equations with High Frequence Oscillating Source Terms

DU Tao   

  1. Tsinghua University; Beijing 100084; China
  • Received:2001-10-23 Revised:2002-10-09 Online:2003-04-25 Published:2003-04-25

摘要: 通过模型方程分析看到了一个重要的现象, 如果源项涉及的时间尺度远小于对流项的时间尺度, 那么基于对流时间尺度作为步长的传统数值方法, 即使源项相对于墩流项是个小量, 也会导致平均尺度上错误的结果。为了克服这种困难, 采用时间分裂方法, 把方程分裂成含对流项部分的偏微分方程( PDE) 和包含源项的常微分方程(ODE) 基础上, PDE 使用传统的数值方法, ODE 用解析的方法求解。该混合方法在数值格式时间步长小于平均流动时间尺度时, 得到正确解, 而与点源所隐含的时间尺度无关。把这个方法应用在含源流动的Euler 方程的计算中, 计算了翼型振荡问题, 取得了理想的结果。

关键词: 源项, 刚性问题, 混合解析/数值方法

Abstract: A model equation shows that the conventional method gives very wr ong results for the time step just basedon the advection part, if the time scale of the source term is much smaller than that of the advection par t. Even ift he sour ce term is a small quantity compared with the advection par t, the large scale structure is completely wrong.In or der to over come the numerical difficulty, consider the operator splittingmethod. After splitting, the mean flowpart ( advection part) will be treated using any suitable numerical method, while the sour ce step ( ODE) is solvedanalytically. With such a met hod, one can use a time step just based on the t ime scale of the advection par t to cap2ture the cor rect large scale, independently of how small the time scale of the sour ce term is. This method is appliedto t he Euler equations with an oscillation source term. The satisfying results are given.

Key words: source terms, stiffness problem, mixed analytical/n umerical scheme

中图分类号: