航空学报 > 2024, Vol. 45 Issue (S1): 730598-730598   doi: 10.7527/S1000-6893.2024.30598

光滑型TENO非线性加权六阶低耗散WCNS-CU6-ST格式

吴文昌1, 韩省思2(), 闵耀兵1, 燕振国1, 马燕凯1   

  1. 1.空天飞行空气动力科学与技术全国重点实验室,绵阳 621000
    2.南京航空航天大学 能源与动力学院,南京 210016
  • 收稿日期:2024-04-24 修回日期:2024-05-13 接受日期:2024-07-09 出版日期:2024-08-07 发布日期:2024-08-05
  • 通讯作者: 韩省思 E-mail:xshan@nuaa.edu.cn
  • 基金资助:
    国家自然科学基金(52376114);国家科技重大专项(2017-III-0005-0029)

Sixth-order low-dissipation WCNS-CU6-ST scheme based on smooth TENO nonlinear weight

Wenchang WU1, Xingsi HAN2(), Yaobing MIN1, Zhenguo YAN1, Yankai MA1   

  1. 1.State Key Laboratory of Aerodynamics,Mianyang 621000,China
    2.College of Energy and Power Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China
  • Received:2024-04-24 Revised:2024-05-13 Accepted:2024-07-09 Online:2024-08-07 Published:2024-08-05
  • Contact: Xingsi HAN E-mail:xshan@nuaa.edu.cn
  • Supported by:
    National Natural Science Foundation of China(52376114);National Science and Technology Major Project(2017-III-0005-0029)

摘要:

空天飞行器的流动环境涉及激波与湍流等复杂流动现象,流场中激波间断与多尺度湍流结构共存使得高阶精度数值模拟面临巨大挑战。为了提升对复杂湍流流动问题模拟的精度,在加权紧致非线性格式(WCNS)框架下,发展了六阶中心/五阶迎风线性混合插值方法;同时考虑到高分辨率的激波捕捉特性,TENO格式中的尺度分离策略被用于发展非线性权函数,针对数值格式的数值收敛性,进一步发展了光滑型S-TENO非线性权函数。通过近似色散关系(ADR)、一维流动问题、三维槽道湍流和跨声速翼型自适应湍流模拟的算例测试,发现线性模板采用六阶中心/五阶迎风的混合插值方法,能够显著降低格式的数值耗散,从而改善槽道湍流的预测精度;同时光滑型S-TENO非线性权函数不仅保持了原TENO格式非线性权函数良好的色散耗散特性,并在跨声速翼型流动的激波/湍流边界层干扰模拟算例中,提升了数值收敛性,数值算例结果表明本文提出的基于光滑型S-TENO权函数的六阶低耗散WCNS-CU6-ST格式能够同时具备数值耗散低与收敛良好的数值特性,有助于进一步提升对激波/湍流等复杂流动问题模拟的准确性。

关键词: 高阶精度格式, 有限差分方法, 非线性权函数, 激波捕捉, 湍流模拟, 数值耗散, 数值收敛性

Abstract:

The flow environment of spacecraft involves complex flow phenomena such as shock waves and turbulence. The coexistence of shock wave discontinuities and multiscale turbulent structures in the flow field poses significant challenges to high-order accurate numerical simulations. To effectively reduce numerical dissipation in turbulence simulations, a sixth-order central/fifth-order upwind linear hybrid interpolation was developed under the framework of Weighted Compact Nonlinear Scheme (WCNS). Considering the high-resolution shock capturing characteristics, the scale separation strategy in the TENO scheme is used to develop nonlinear weight functions. For the numerical convergence of the numerical scheme, a smooth S-TENO nonlinear weight function is further developed. Through Approximate Dispersion Relations (ADR), one-dimensional convection problems, and self-adaptive turbulence eddy simulation of three-dimensional channel, it was found that the linear stencil using a sixth-order central/fifth-order upwind linear hybrid interpolation stencil can significantly reduce the numerical dissipation of the scheme, thereby improving the prediction accuracy of channel flow and transonic flow past an airfoil. The proposed smooth S-TENO nonlinear weight not only maintains the good dispersion and dissipation properties of the original TENO nonlinear weight, but also significantly improves the numerical convergence in the cases of shock wave/turbulent boundary layer interaction when transonic flow past an airfoil, thus improving the accuracy of computational results. The numerical results indicate that the optimized sixth-order low-dissipation WCNS-CU6-ST scheme based on smooth S-TENO weight can achieve both low numerical dissipation and good convergence, which is helpful to improve the accuracy of simulations in complex shock/turbulence problems.

Key words: high-order accurate numerical scheme, finite differencing method, nonlinear weight function, shock capturing, turbulence simulation, numerical dissipation, numerical convergence

中图分类号: