Bayesian uncertainty analysis of SA turbulence model for supersonic jet interaction simulations

The Reynolds Averaged Navier-Stokes(RANS) models are still the workhorse in current engineering applications due to its high efficiency and robustness. However, the closure coefficients of RANS turbulence models are determined by model builders according to some simple fundamental flows, and the suggested values may not be applicable to complex flows, especially supersonic jet interaction flow. In this work, the Bayesian method is employed to recalibrate the closure coefficients of Spalart-Allma...     

TRIP2．0软件湍流模块的开发与确认

基于TRIP2.0_SOLVER数值模拟软件,开发了湍流模型方法库,目前该CFD软件中集成的湍流模型主要包括B-L代数湍流模型、SA一方程模型和SST两方程模型.为了考核不同的湍流模型工程适用性,本文采用对接网格技术,通过有限体积法数值离散三维任意坐标系下的RANS方程组,应用LU分解、MUSCL差分格式和低雷诺数SA和SST两种湍流模型,数值模拟了二维NLR-7301两段翼型和三维DLR-F6翼身组合体的绕流流场,计算与试验比较的内容包括了表面的压力系数分布、典型站位的速度型和气动特性曲线等内容.通过计算与试验的比较,在本文的计算范围内,采用两种湍流模型均可以得到与试验结果相吻合的压力分布和升力曲线,但在边界层内的速度型、粘性阻力和力矩特性等方面,不同的湍流模型具有明显的差异.     

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

基于SA和SSTγ-■模型对边界层转捩预测的比较

利用作者共同研发的In-house代码TRANS3D平台,在SA和SST两种常用湍流模型框架下构建了γ-■转捩模型,并以二维低速S809层流翼型和三维小后掠跨音速F5层流机翼为对象,比较了两种不同湍流模型构架下γ-■转捩模型预测的气动力特性和流场分布。结果表明:两种转捩模型均能预测航空转捩计算中常见的分离流转捩和自然转捩类型,明显改善了中低雷诺数流场下的预测精度,但由于选取基准湍流模型的不同,基于SA和SST的γ-■转捩模型在流动细节上依然存在着一些差异。     

风力机标模非定常数值模拟中的影响因素研究

基于自主研发的“亚跨超 CFD 软件平台”(TRIP3.0)，采用刚性运动网格技术和动态拼接网格技术，开发了针对旋转类机械的非定常求解模块。本文开展了 NREL Phase VI 风力机标模非定常数值模拟中的影响因素研究，影响因素主要包括不同时间步长的影响、不同湍流模型的影响、刚性动网格技术和动态拼接网格技术的影响三个方面。本文数值模拟采用的控制方程为雷诺平均 N-S 方程，采用有限体积法离散控制方程，离散方程的时间方向采用“双时间步”方法求解，空间方向无粘项离散采用 MUSCL-Roe 格式，湍流模型采用 SA 和 SST 湍流模型，并引入多重网格和并行计算技术加速求解。数值模拟结果表明：时间步长、湍流模型、网格方案等因素主要影响风力机叶片吸力面的流动结构，进而影响吸力面的压力分布，而对压力面的流动结构和压力分布基本没有影响；采用刚性运动网格结合 SA 湍流模型所得到的压力分布更接近实验值。     

湍流模型系数不确定度对翼型绕流模拟的影响

使用非嵌入式多项式混沌方法研究了湍流模型系数的不确定度对RAE2822跨声速翼型绕流模拟的影响。计算中关注了数值模拟的积分量（升力系数、阻力系数）和局部量（壁面压力、摩擦系数和空间马赫数分布）的不确定度量化结果。首先，从单输入变量入手，研究卡门常数的不确定度对数值模拟的影响。然后，同时考虑Spalart-Allmaras模型中9个参数的不确定度带来的影响。通过多项式混沌展开，得到系统输出对不确定输入变量的响应，由此可以得到输出的统计特性，包括平均值、方差和极值等信息。最后，在多变量不确定度量化过程中，通过Sobol指标来量化每个输入变量的不确定度对输出不确定度的贡献程度。本文计算只考虑了RAE2822跨声速翼型模拟的单一计算状态，影响规律是否可以推及其他工况和算例需要进一步检验。     

A new simpler rotation/curvature correction method for Spalart-Allmaras turbulence model

 A new and much simpler rotation and curvature effects factor, which takes the form of Richardson number suggested by Hellsten originally for SST k-x model, is presented for Spalart and Shur's rotation and curvature correction in the context of Spalart-Allmaras (SA) turbulence model. The new factor excludes the Lagrangian derivative of the strain rate tensor that exists in the SARC model, resulting in a simple, efficient and easy-to-implement approach for SA turbulence model (denoted as SARCM) to account for the effects of system rotation and curvature, techniquely. And then the SARCM is tested through turbulent curved wall flows: one is the flow over a zeropressure- gradient curved wall and the other is the channel flow in a duct with a U-turn. Predictions of the SARCM model are compared with experimental data and with the results obtained using original SA and SARC models. The numerical results show that SARCM can predict the rotation-curvature effects as accurately as SARC, but considerably more efficiently. Additionally, the accuracy of SARCM might strongly depend on the rotation-curvature model constants. Suggesting values for those constants are given, after some trials and errors.     

SA和SST湍流模型对高超声速边界层强制转捩的适应性

凸起物是高超声速流动中常用的一种人工转捩装置。采用高阶精度算法模拟了高超声速进气道压缩面上的强制转捩流动,转捩装置为一排高度为1mm的钻石型凸起物和斜坡型凸起物。考察了Spalart-Allmaras(SA)模型和剪切应力输运(SST)湍流模型对该问题的适应性。在考察过程中通过丰富的算例分析了网格规模、可压缩修正和空间离散格式等对计算结果的影响。在层流区,计算能与试验取得非常一致的结果。但在湍流区,计算得到的热流通常高出试验数据。经分析发现其原因是强制转捩的湍流边界层与自然转捩的湍流边界层在涡结构上存在较大差别,使得湍流模型的效果较差。针对高超强制转捩湍流涡结构丰富的特点,对SST湍流模型进行了修改。计算结果表明,该修改方法对提高热流精度具有一定效果。     

30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。