基于转捩/尺度适应模型与FW-H声学方程的气动噪声数值模拟研究

采用计算流体力学(CFD)同"声类比"相结合的方法进行噪声模拟,利用CFD数值模拟 MD30P30N多段翼型失速攻角附近流场以及其气动特性来校核近场精度,进一步通过结合可穿透数据面的FW-H声学方程进行气动噪声分析。为了准确捕捉近场流场信息,为噪声预测提供可靠的声源精度,本文基于k-ω剪应力输运(SST)湍流模型,建立了尺度适应分离流(SAS)模型,并采用γ-Re_θt转捩模型耦合k-ω SST湍流模型建立了边界层转捩数值模拟技术;充分利用尺度适应模型在边界层表现为雷诺平均(RANS)方法这一特点,将γ-Re_θt转捩模型与SAS模型结合,建立针对包含转捩、分离现象复杂流场的数值模拟技术。文中以RANS为控制方程,分别采用全湍流k-ω SST模型、γ-Re_θtSAS转捩/分离流模型对多段翼进行近场数值模拟计算,结合可穿透数据面的FW-H声学方程进行气动噪声预测,在分析其对流场及气动噪声影响的基础上,得出了几点结论。     

绕振荡水翼流动及其转捩特性的数值计算研究

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明:与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.      

后掠机翼低速流动转捩位置的升华法测量

报导了在西北工业大学低湍流度风洞中采用升华法测量后掠翼上表面流动转捩位置及其随攻角变化而改变情况,并参照后掠翼压力分布的实验结果进行分析.本文涉及的升华法,是根据边界层的层流区和湍流区流态差别来判断边界层转捩的一种测量方法.升华法的测量结果表明:在小攻角范围内,转捩位置逐渐提前,当达到某一攻角时,转捩快速接近前缘.这一结果和压力分布的实验结果相吻合,说明升华法能够较准确地显示转捩线位置.由于影响后掠翼转捩位置的因素较多,诸如边界层内横流的不稳定性、壁面的干扰、后掠翼的三维流动效应等,致使出现了攻角为4°时转捩位置呈现出一条斜线的现象.     

k-ω-γ模式对转捩影响因素的预测性能研究

高超声速边界层转捩的准确预测对飞行器的防热、减阻至关重要,而影响高超声速边界层转捩的因素众多.从模式角度出发研究边界层转捩的影响因素,采用k-ω-γ 转捩模式对5°圆锥的边界层转捩进行了数值分析,计算了不同头部钝度、来流雷诺数和湍流度情况下的边界层转捩,并与实验结果进行了对比. 研究结果表明:k-ω-γ 转捩模式基本能够反映头部钝度、来流雷诺数、来流湍流度对高超声速圆锥边界层转捩的影响规律,但对转捩后的热流峰值预测不准;从模式构造角度分析发现,雷诺数越高或头部钝度越小,层流区边界层越薄,k-ω-γ 转捩模式中第一、第二模态时间尺度增大,因此转捩起始位置提前;来流湍流度越大,等效脉动动能初值越大,导致层流区发展过程中等效脉动动能越大,因此转捩易于发生.      

两种湍流模型在高速旋转翼身组合弹箭中的对比研究

弹箭设计、弹道计算和稳定性研究都需要准确预测旋转弹箭的马格努斯力和力矩,国内针对旋转弹箭气动特性的数值模拟工作集中在旋成体上,对带翼外形进行完全时间相关的非定常研究鲜有见到;国外虽然有对带翼外形开展研究,但以验证方法为主,对湍流模型在复杂外形弹箭旋转中的研究未曾见到.采用完全时间相关的非定常N-S方程,对带翼弹箭开展计算,对比了一方程SA(Spalart-Allmaras)湍流模型和两方程k-ωSST(shear-stress-transport)湍流模型对马格努斯效应产生的影响,并分析了旋转导致的边界层和涡非对称畸变,以及周向压力分布和剪切应力分布非对称畸变.结果表明:旋转引起的物面流场参数变化主要体现在弹体中后部,SA和SST湍流模型预测的全弹马格努斯特性与阿诺德工程发展中心(Arnold Engineering Development Center,AEDC)实验及陆军研究实验室(Army Research Laboratory,ARL)的计算结果一致性很好,对动导数而言两湍流模型计算精度相当.两湍流模型计算的弹体左侧流场参数差异比右侧大,分析认为正向旋转使左侧壁面速度方向与来流速度相反,相互阻碍使气流脉动效应更强.壁面附近湍流黏性系数SA结果大于SST结果,y=0截面物面压力SA结果小于SST结果、最大相差6%,摩阻系数SA结果大于SST结果、最大相差35%.SA对旋转产生的分离抑制作用强于SST.     

曲率对机翼边界层二次失稳影响

采用非线性抛物化扰动方程(NPSE)计算了层流机翼NLF0415(2)在几个工况下横流涡的非线性发展, 应用Floquet理论分析了横流涡的二次失稳. 较系统地分析了后掠机翼在多种参数下表面曲率对其流动稳定性的影响. Haynes证明曲率对于横流在线性稳定性计算(LST)和非线性抛物化扰动方程计算中(NPSE)都起着非常明显的稳定作用; 然而, 该文计算结果表明, 曲率在二次失稳计算中影响不大.      

两种湍流模型在高速旋转翼身组合弹箭中的对比研究1）

弹箭设计、弹道计算和稳定性研究都需要准确预测旋转弹箭的马格努斯力和力矩,国内针对旋转弹箭气动特性的数值模拟工作集中在旋成体上,对带翼外形进行完全时间相关的非定常研究鲜有见到;国外虽然有对带翼外形开展研究,但以验证方法为主,对湍流模型在复杂外形弹箭旋转中的研究未曾见到.采用完全时间相关的非定常N-S方程,对带翼弹箭开展计算,对比了一方程SA（Spalart-Allmaras）湍流模型和两方程k-!SST（shear-stress-transport）湍流模型对马格努斯效应产生的影响,并分析了旋转导致的边界层和涡非对称畸变,以及周向压力分布和剪切应力分布非对称畸变.结果表明：旋转引起的物面流场参数变化主要体现在弹体中后部,SA和SST湍流模型预测的全弹马格努斯特性与阿诺德工程发展中心（Arnold Engineering Development Center,AEDC）实验及陆军研究实验室（Army Research Laboratory,ARL）的计算结果一致性很好,对动导数而言两湍流模型计算精度相当.两湍流模型计算的弹体左侧流场参数差异比右侧大,分析认为正向旋转使左侧壁面速度方向与来流速度相反,相互阻碍使气流脉动效应更强.壁面附近湍流黏性系数SA结果大于SST结果,y=0截面物面压力SA结果小于SST结果、最大相差6%,摩阻系数SA结果大于SST结果、最大相差35%.SA对旋转产生的分离抑制作用强于SST.     

绕椭球的低速流动研究

理解和预测绕椭球的流动对指导飞行器和潜艇等交通工具的设计具有很强的工程意义. 近年来, 针对椭球绕流开展了大量的实验和数值模拟研究. 对有攻角下椭球绕流分离的定性描述和定量研究, 促进了对三维分离的辨识和拓扑研究. 文章对流场特性进行了分析, 介绍了分离对气动力、噪声、尾迹的影响, 以及实验条件对流动的影响. 上述定常流动与非定常机动过程之间存在明显差异, 非定常机动过程不能作为定常或准定常问题处理, 在机动过程中, 分离出现明显延迟, 气动力出现明显变化. 随后介绍了数值模拟在求解绕椭球流动中的进展, 当前求解雷诺平均的N-S方程湍流模式仍然是解决绕椭球大范围分离流动的主要工程方法, 大涡模拟和分离涡模拟等也逐渐得到了广泛应用. 受限于计算能力, 直接数据模拟只能用于较低雷诺数, 在高雷诺数流动中还不适用. 非定常机动过程的数值模拟较定常状态, 与实验结果的差距要大一些. 最后, 介绍了对椭球绕流场转捩的研究进展, 对T-S转捩与横流转捩的机理和辨识已经较为准确, 数值模拟结果与实验结果基本相符, 但对再附转捩的认识还不够清晰, 尤其是迎风面, 因此椭球绕流转捩的研究还需要依靠实验.      

绕振荡水翼流动及其转捩特性的数值计算研究简

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明：与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.     

大迎角横向流动分离控制器设计

大迎角下横向流动引起的外翼过早分离是后掠翼飞机常见的一大问题,对此提出在机翼下翼面安装分离流控制器来抑制流动分离的方案.二维计算中采用Gamma-Theta转捩模型,三维研究使用SST湍流模型,研究发现分离流控制器的设计使流过上翼面的高能气流通量增加,这些高能流体具有更强的抵抗流向逆压梯度的能力,同时阻碍了边界层向外流动,从而改善了外翼的失速性能.通过参数化设计所得的分离流控制器与原始构型相比失速迎角增加,最大升力系数显著提高,阻力系数显著减小,达到了设计目标.     

Aerodynamic analysis and validation by using three turbulence models for narrow trench configuration

Three turbulence models SST Gamma Theta, k-ω and k-ε which all found in ANSYS CFX were used. Velocity contours, pressure coefficient profiles and turbulence levels contours were discussed. Results indicate that three models calculated the cross flow boundary layer with different thicknesses. This leads to a difference in calculation each the momentum of the cross flow fluid closer to the jet exit and cooling performance. At M = 0.5 from range X/D = 5:40, the average centerline effectiveness values in case SST Gamma Theta model and k-ω model were 21 and 9.2 % larger than k-ε model, respectively. While at high blowing ratios 1 and 1.4, the effectiveness values in case SST Gamma Theta model and k-ω model were 16, 24.4 % and 32, 46 % less than k-ε model, respectively. k-ω model shows a larger backflow regions than others which a maximum negative value of u/u_∞ = ?0.198 is reached in zone 4.     

A numerical investigation of CO_2 dilution on the thermochemical characteristics of a swirl stabilized diffusion flame

The turbulent combustion flow modeling is performed to study the effects of CO₂ addition to the fuel and oxidizer streams on the thermochemical characteristics of a swirl stabilized diffusion flame. A flamelet approach along with three well-known turbulence models is utilized to model the turbulent combustion flow field. The k-ω shear stress transport(SST) model shows the best agreement with the experimental measurements compared with other models. Therefore, the k-ω SST model is used...     

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

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

Robust design of natural laminar flow supercritical airfoil by multi-objective evolution method

A transonic, high Reynolds number natural laminar flow airfoil is designed and studied. The γ-θ transition model is combined with the shear stress transport(SST)k-w turbulence model to predict the transition region for a laminar-turbulent boundary layer. The non-uniform free-form deformation(NFFD) method based on the non-uniform rational B-spline(NURBS) basis function is introduced to the airfoil parameterization.The non-dominated sorting genetic algorithm-II(NSGA-II) is used as the search algorithm, and the surrogate model based on the Kriging models is introduced to improve the efficiency of the optimization system. The optimization system is set up based on the above technologies, and the robust design about the uncertainty of the Mach number is carried out for NASA0412 airfoil. The optimized airfoil is analyzed and compared with the original airfoil. The results show that natural laminar flow can be achieved on a supercritical airfoil to improve the aerodynamic characteristic of airfoils.     

Multi-layer analytic solution for k-ω model equations via a symmetry approach

Despite being one of the oldest and most widely-used turbulence models in engineering computational fluid dynamics (CFD), the k-ω model has not been fully understood theoretically because of its high nonlinearity and complex model parameter setting. Here, a multi-layer analytic expression is postulated for two lengths (stress and kinetic energy lengths), yielding an analytic solution for the k-ω model equations in pipe flow. Approximate local balance equations are analyzed to determine the key parameters in the solution, which are shown to be rather close to the empirically-measured values from the numerical solution of the Wilcox k-ω model, and hence the analytic construction is fully validated. The results provide clear evidence that the k-ω model sets in it a multilayer structure, which is similar to but different, in some insignificant details, from the Navier-Stokes (N-S) turbulence. This finding explains why the k-ω model is so popular, especially in computing the near-wall flow. Finally, the analysis is extended to a newly-refined k-ω model called the structural ensemble dynamics (SED) k-ω model, showing that the SED k-ω model has improved the multi-layer structure in the outer flow but preserved the setting of the k-ω model in the inner region.     

弯度翼型大迎角时分离涡流场的多平衡态现象

基于kω的SST两方程湍流模型,在时间域求解雷诺平均Navier-Stokes方程,模拟弯度翼型大迎角时的分离流动。通过给翼型施加一定形式的扰动,重点关注了翼型弯度对大迎角分离涡流场平衡态转移的影响。研究结果表明:与相同厚度20%以上的对称翼型相比,2%弯度的翼型出现分离涡流场平衡态转移的起始迎角变小2°左右,迎角区间变宽约1°;在厚度相对较小的NACA2416翼型上也发现上述分离涡平衡态转移现象。由此说明翼型弯度在一定程度上促使了分离涡平衡态的转移。     

Benchmark PIV database for the validation of CFD simulations in a transitional cavity flow

We present an experimental benchmark database for the transitional cavity flow. The database is obtained by planar Particle Image Velocimetry measurements at the median plane of the cavity model, for Reynolds numbers between 6300 and 19,000 based on the cavity height. A detailed uncertainty analysis of the experimental results is performed via the correlation statistics method for PIV uncertainty quantification and linear error propagation.The experimental results are compared to two-dimensional Reynolds-Averaged Navier Stokes (RANS) numerical simulations with different turbulence models. It is shown that, when the standard k-ω turbulence model is employed, the discrepancy between numerical simulations and experimental results exceeds the uncertainty of the latter. Conversely, RANS simulations with the SST k-ω turbulence model agree well with the experimental data in terms of time-averaged flow properties; however, the turbulent kinetic energy results present significant discrepancies at all considered Reynolds numbers. The data presented in this paper is made available for open-access download via the 4TU.ResearchData repository with DOI: https://doi.org/10.4121/14061233.     

Applications of EPSE method for predicting crossflow instability in swept-wing boundary layers

The nth-order expansion of the parabolized stability equation(EPSEn) is obtained from the Taylor expansion of the linear parabolized stability equation(LPSE) in the streamwise direction. The EPSE together with the homogeneous boundary conditions forms a local eigenvalue problem, in which the streamwise variations of the mean flow and the disturbance shape function are considered. The first-order EPSE(EPSE1) and the second-order EPSE(EPSE2) are used to study the crossflow instability in the swept NLF(2)-0415 wing boundary layer. The non-parallelism degree of the boundary layer is strong. Compared with the growth rates predicted by the linear stability theory(LST),the results given by the EPSE1 and EPSE2 agree well with those given by the LPSE.In particular, the results given by the EPSE2 are almost the same as those given by the LPSE. The prediction of the EPSE1 is more accurate than the prediction of the LST, and is more efficient than the predictions of the EPSE2 and LPSE. Therefore, the EPSE1 is an efficient e~N prediction tool for the crossflow instability in swept-wing boundary-layer flows.