航空学报 > 1998, Vol. 19 Issue (6): 67-70

有限元法求解QPNS方程模拟高速旋成体气动力和热

张国富, 花锡青   

  1. 南京航空航天大学601 教研室, 南京, 210016
  • 收稿日期:1996-12-17 修回日期:1998-04-06 出版日期:1998-12-25 发布日期:1998-12-25

SOLUTION OF QPNS EQUATIONS TO SIMULATE AERODYNAMIC FORCES AND HEATING FOR ROTATION BODIES AT HIGH SPEED USING FINITE ELEMENT METHOD

Zhang Guofu, Hua Xiqing   

  1. Faculty 601, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016
  • Received:1996-12-17 Revised:1998-04-06 Online:1998-12-25 Published:1998-12-25

摘要: 从FNS方程中略去流向粘性导数项导出熵变量形式QPNS方程(拟抛物化N-S方程)其离散解能自动满足熵增不等式。利用GLS(伽辽金/最小二乘法)有限元法构造出QPNS方程的弱解式,所形成的整体方程组用GMRES(广义最小余数法)求解。采用空间推进求解QPNS方程,能明显降低计算机内存和运算时间。算例给出典型旋成体的物面压强、剪应力、热流量分布以及流场中密度、温度分布。

关键词: GLS有限元法, QPNA方程, 旋成体

Abstract: The QPNS equations in the entropy variable are derived from the FNS equations by neglecting streamwise viscous derivatives. The discrete solution always satisfies the entropy production inequality. The QPNS equations are solved with a spacemarching method. The step back strategy is used to provide the data of the initial plane. The Galerkin/least squares finite element method is employed to formulate the weak solution of QPNS equations. The nonsymmetric linear equation systems are solved by GMRES algorithm. The nodal block diagonal preconditioning technique is used to accelerate the convergence rate efficiently. The present method markedly reduces the amount of computer storage and run time. It shows a good stability as well. Numerical examples are the supersonic flows past the parabolic body and tangent ogive body. The results include the pressure, viscous stress, heat flux on the body surface, and the density, temperature in the flow field of windward and leeward.

Key words: GLS finite element method, QPNS equations, rotation body

中图分类号: