首页 >

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

杨小权1, 程苏堃1, 杨爱明1,2, 孙刚1   

  1. 1. 复旦大学 力学与工程科学系, 上海 200433;
    2. 上海大学 上海市应用数学和力学研究所, 上海 200072
  • 收稿日期:2012-04-06 修回日期:2012-05-27 出版日期:2013-04-25 发布日期:2013-04-23
  • 通讯作者: 杨爱明,Tel.: 021-55665534 E-mail: amyang@fudan.edu.cn E-mail:amyang@fudan.edu.cn
  • 作者简介:杨小权 男, 博士研究生。主要研究方向: 计算流体力学和设计空气动力学。 E-mail: 09110290008@fudan.edu.cn;杨爱明 男, 博士, 副教授, 硕士生导师。主要研究方向: 计算流体力学和设计空气动力学。 Tel: 021-55665534 E-mail: amyang@ fudan.edu.cn
  • 基金资助:

    国家自然科学基金(11172070);上海市自然科学基金(09ZR1411400);复旦大学研究生创新基金(EYH2126022)

Time Spectral Method for Numerical Simulation of Unsteady Viscous Flow over Oscillating Airfoil and Wing

YANG Xiaoquan1, CHENG Sukun1, YANG Aiming1,2, SUN Gang1   

  1. 1. Department of Mechanics and Mechanical Engineering Science, Fudan University, Shanghai 200433, China;
    2. Shanghai Institute of Applied Mathematics and Mechanics, Shanghai University, Shanghai 200072, China
  • Received:2012-04-06 Revised:2012-05-27 Online:2013-04-25 Published:2013-04-23
  • Supported by:

    National Natural Science Foundation of China (11172070); Shanghai Municipal Natural Science Foundation (09ZR1411400); Graduates' Innovation Foundation (EYH2126022) *Corresponding author. Tel.: 021-55665534 E-mail: amyang@fudan.edu.cn

摘要:

传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。

关键词: 周期性非定常流动, 时间谱方法, 双时间方法, 俯仰振荡翼型/机翼, 多重网格, 数值模拟

Abstract:

The traditional dual time stepping formula is comparatively inefficient in solving unsteady problems due to its repeated transitional iteration stepping. In view of the characteristics of periodical unsteady flow, a high efficiency time spectral method based on discrete fourier transformation is hereby developed for such type of flow over an oscillating airfoil and wing. Roe's flux difference splitter scheme is utilized to discretize the convective terms of the time-space coupled Reynolds-averaged Navier-Stokes (RANS) equations. Time spectral method is used to treat the physical time derivative terms while the implicit LU-SGS (Lower-Upper Symmetric Gauss-Seidel) scheme deals with pseudo time stepping. Meanwhile, the physi-cal time terms of Spalart-Allmaras's one equation turbulence model, which acts as our turbulence model, is also discretized by the time spectral method in consideration of the turbulence's time-space coupling effect. In order to further improve efficiency, methods like local time step and multigrid algorithm are used to accelerate the convergence of RANS equation's solution. In the numerical examples section, the periodical unsteady flow field over an oscillating NACA0012 airfoil and Lann wing is simulated numerically. Results show that the time spectral method possesses obvious advantages in reducing computational cost and improving calculation accuracy as compared with the conventional dual time stepping formula.

Key words: time-periodic unsteady flow, time spectral method, dual time stepping method, oscillating airfoil/wing, multigrid, numerical simulation

中图分类号: