EI、Scopus 收录
中文核心期刊
杜一鸣, 高正红, 舒博文, 邱福生, 宋辰星. k-ω SST湍流模式三维激波分离流动修正. 力学学报, 2022, 54(6): 1485-1501. DOI: 10.6052/0459-1879-22-065
引用本文: 杜一鸣, 高正红, 舒博文, 邱福生, 宋辰星. k-ω SST湍流模式三维激波分离流动修正. 力学学报, 2022, 54(6): 1485-1501. DOI: 10.6052/0459-1879-22-065
Du Yiming, Gao Zhenghong, Shu Bowen, Qiu Fusheng, Song Chenxing. Three-dimensional shock separated flow corrections. DOI: 10.6052/0459-1879-22-065
Citation: Du Yiming, Gao Zhenghong, Shu Bowen, Qiu Fusheng, Song Chenxing. Three-dimensional shock separated flow corrections. DOI: 10.6052/0459-1879-22-065

k-ω SST湍流模式三维激波分离流动修正

THREE-DIMENSIONAL SHOCK SEPARATED FLOW CORRECTIONS OF k-ω SST MODEL

  • 摘要: “激波−边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.

     

    Abstract: Shock/boundary-layer separation is a typical turbulence non-equilibrium flow in the field of aeronautical aerodynamics. Accurate simulation of shock separated flow is of great significance for aerodynamic performance evaluation and optimization of transonic aircraft. The definition of eddy viscosity coefficient in conventional eddy-viscosity turbulence models (EVM), however, is not suitable for non-equilibrium flow. The Bradshaw assumption introduced by k-ω SST turbulence model for this purpose instead restricts the generation of Reynolds stress when applied to three-dimensional flow with strong adverse pressure gradient and large separation, which results in the invalidity of k-ω SST model as well as other commonly used EVMs for this kind of flow. Moreover, the existing nonlinear constitutive relation of Reynolds stress cannot effectively improve the simulation accuracy. To this end, two shock separated flow correction methods respectively based on Bradshaw assumption and length scale are proposed for k-ω SST model. The former correction relaxes the limitation of Reynolds stress generation by increasing Bradshaw constant. While based on the concept of turbulence length scale, the latter correction constructs a modified function for the dissipation term of the ω equation by using the mixing length theory, the generation-dissipation ratio of turbulent kinetic energy and a newly defined ratio of length scale to improve the modeling length scale in three-dimensional shock separated flow. The two methods both get better simulation results for the transonic flow of ONERA M6 wing at high angle of attack than those of Reynolds stress model. Further Reynolds stress analysis reveals that the concept of "major Reynolds-stress component" in three-dimensional shock separated flow is no longer tenable since the magnitude of each Reynolds-stress component is close. The grid convergence analysis, and verifications on other angles of attack and the wall-function law of turbulent boundary layer on the flat plate further confirm the validity, applicability and universality of the proposed correction methods.

     

/

返回文章
返回