NS方程中耦合转捩判断的风力机翼型气动力数值模拟

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合Michel经验判据和Chen-Thy-son转捩模型对风力机翼型S809进行计算。由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用Michel经验判据分析判断层流边界层的解得到转捩点的位置,使用Chen-Thyson转捩模型得到从层流到湍流的转捩过渡区,这样随着流场的迭代求解,求解器自动更新判断出转捩点位置。对S809翼型进行全湍流和耦合转捩判断的RANS方程计算,可以看到考虑转捩判断后得到的升阻力系数与实验值吻合较好,验证了方法的可行性。     

耦合转捩自动判断的机翼粘性绕流计算与优化设计

具有边界层转捩自动判断功能的三维雷诺平均Navier-Stokes方程(RANS)方程求解器是自然层流机翼设计的必要手段。将一种基于线性稳定性理论的NTS/NCF转捩判断方法与RANS求解器耦合,以提高求解器的计算精度。通过求解三维稳定性方程,然后使用包络线方法计算对应Tollmien-Schlichting不稳定性的NTS,使用固定β方法计算对应横流不稳定性的NCF,然后根据转捩阀值(NTS)tr和(NCF)tr判断转捩位置。在上述工作的基础上,开展了自然层流机翼平面形状优化的气动减阻设计研究。     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

基于双eN方法的翼身组合体流动转捩自动判断

发展翼身组合体复杂外形流动转捩自动判断方法,对高亚声速民机自然层流（NLF）机翼设计具有重要意义。使用多块结构化网格和三维雷诺平均Navier-Stokes（RANS）方程求解器,耦合边界层方程求解和基于线性稳定性理论（LST）的完全双e^N方法,发展了一套可同时计及Tollmien-Schlichting波和横流不稳定性扰动诱导转捩的翼身组合体流动转捩自动判断方法。对DLR-F4翼身组合体绕流进行了转捩自动判断,将得到的转捩位置与试验结果进行比较,验证了所发展方法的正确性。使用上述方法对配置自然层流机翼的中短程民机翼身组合体外形进行了数值模拟,并将结果与单独机翼的转捩位置进行了对比,结果表明机身三维位移效应增强了自然层流后掠机翼边界层的横流不稳定性强度,导致翼根转捩位置提前至前缘区。     

基于响应面法的低速翼型气动优化设计

响应面方法相较于其它直接优化方法有其高效、实用的优势,此前的研究更多地将响应面方法用于超音速和跨音速翼型的减阻优化设计中。本文将此方法应用于低速翼型优化设计中,进行了基于RANS(Reynolds-Aver-aged Navier-Stokes)方程和自由转捩预测耦合求解的低速翼型气动优化设计。通过计算附面层方程得到附面层参数并用en方法计算转捩位置,并考虑了T-S波和层流分离造成的转捩。RANS方程计算中,使用了转捩过渡区模型,以保证附面层外边界压力分布的精度。RANS方程和转捩预测迭代进行至转捩位置收敛。在响应面模型计算中,使用不含二阶交叉项的二阶多项式模型,减少了构造模型所需的计算量;合理的选择设计空间保证了构造的响应面模型具有较高的精度。使用三个设计点的多目标优化设计,保证了设计的合理性。通过对NACA64(1)-112翼型优化计算结果表明,本文的方法可以有效地进行低速翼型的气动优化,各设计点上转捩位置也得到了改善,有较好的工程实用前景。     

基于eN-数据库方法复杂构型飞机转捩预测

为探索边界层转捩对大型运输机在起降条件和有较大层流区的巡航条件下的气动力精确计算问题，通过在三维RANS求解器中引入eN-数据库方法来预测飞行器表面的转捩位置，并探索转捩对气动力的影响规律。方法与目前流行的基于间歇因子控制方程的转捩预测方法相比，具有计算效率高、易于工程应用、且考虑TS不稳定性转捩因素的特点。在此基础上，通过计算NASA梯形翼来分析起降构型条件下气动力受转捩影响的规律，并通过计算DLR－F6翼身组合体来探索三维构型在巡航条件下的气动力精度。使用eN-数据库转捩判断方法的计算结果与实验值吻合较好，验证了所构建的基于RANS求解器的eN-数据库转捩预测方法的有效性，并为大型运输机气动力精确计算提供了分析工具。     

RANS方程和附面层方程耦合求解转捩位置的方法

研究了基于线性稳定性分析的RANS（Reynolds-Averaged Navier-Stokes）方程和附面层方程耦合求解转捩位置的方法。RANS方程计算得附面层外边界的速度分布作为附面层方程计算的输入,求解附面层方程得到稳定性分析需要的附面层参数。由于附面层的解在分离点处的奇异性,用交互式附面层方程求解方法。用基于线性稳定性分析的e^N方法计算转捩位置,并考虑了Tollmien-Schliching波和层流附面层分离造成的转捩。迭代求解RANS方程和附面层方程,直到转捩位置收敛。RANS求解中使用转捩过渡区模型,避免了点转捩模型引入的数值扰动。通过对NA-CA0012翼型的计算并和实验结果及Xfoil计算的结果进行比较,吻合较好,有较好的工程实用价值。     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。     

采用γ-Re_(θt)模型的转捩流动计算分析

为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Reθt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序(HUNS3D)。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Reθt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明:γ-Reθt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

Numerical research on airfoil transition delay by alternative current dielectric barrier discharge actuation

A Dielectric Barrier Discharge (DBD) plasma actuator can create a body force which locally accelerates the base flow leading to an attenuation of broadband disturbance to delay the transition. In this study, numerical simulation on an NLF0416 airfoil is conducted to investigate transition delay and drag reduction by a DBD plasma actuator. To simulate plasma’s effect more accurately, boundary-layer data is acquired from Reynolds Averaged Navier Stocks (RANS) equations instead of laminar boundary layer equations, although RANS equations need a much finer boundary-layer grid, and the linear stability analysis method is used to analyze the boundary layer and get the transition point. In this study, the influences of different actuation intensities and positions are investigated, and results show that if the actuation intensity is stronger and the actuation position is closer to the base transition point, more drag reduction can be obtained. However, the efficiency of plasma transition delay is really low. For example, when the actuation voltage is 16 kV, the actuation frequency is 1 kHz, and the main Mach number is 0.1, the saved power due to drag reduction is about 5.09 W, but the power consumed is about 32.61 W, and the efficiency is just 15.6%.     

基于双eN方法的短舱层流转捩影响因素

发展自然层流短舱对提升现代民机的经济性和环保性具有重要意义，而对影响短舱层流转捩的因素进行研究有助于更好地开展短舱的层流设计。本文基于线性稳定性分析方法，将双e^N方法同RANS方程求解器耦合，建立了一套可同时计算T-S（Tollmien-Schlichting）波和横流（CF）驻波诱导转捩的流动转捩预测方法，通过标准椭球算例验证了该方法的正确性，进而研究了来流马赫数、雷诺数、湍流度以及迎角对短舱转捩的影响。结果表明：马赫数和迎角会带来压力梯度的明显改变从而引起转捩位置发生变化；而在此构型的高雷诺数工况下，雷诺数和湍流度对转捩位置影响相对较小，转捩位置移动的区域不超过短舱长度的5%。因此在设计阶段，在高雷诺数条件下保持层流设计要尽量避免较大的逆压梯度，保持顺压梯度。     

Validation of the RANS-simulation of laminar separation bubbles on airfoils

This paper presents RANS simulations of the low-Reynolds-number flow past an SD7003 airfoil at $\mathit{Re}=6\times {10}^4$, where transition takes place across a laminar separation bubble. The transition prediction procedure using an approximate envelope method as well as a linear stability solver is discussed. The numerical results are validated against PIV- and force measurements obtained in several wind- and watertunnels and are also compared to XFOIL results. Good agreement is found within the operational range of the airfoil.     

Wing aeroelasticity analysis based on an integral boundary-layer method coupled with Euler solver

An interactive boundary-layer method, which solves the unsteady flow, is developed for aeroelastic computation in the time domain. The coupled method combines the Euler solver with the integral boundary-layer solver(Euler/BL) in a ‘‘semi-inverse" manner to compute flows with the inviscid and viscous interaction. Unsteady boundary conditions on moving surfaces are taken into account by utilizing the approximate small-perturbation method without moving the computational grids. The steady and unsteady flow calculations for the LANN wing are presented. The wing tip displacement of high Reynolds number aero-structural dynamics(HIRENASD) Project is simulated under different angles of attack. The flutter-boundary predictions for the AGARD445.6 wing are provided. The results of the interactive boundary-layer method are compared with those of the Euler method and experimental data. The study shows that viscous effects are significant for these cases and the further data analysis confirms the validity and practicability of the coupled method.     

基于γ-Reθt转捩模型的翼型气动特性计算精度改进研究

将结合SSTk-ω湍流模型的γ-Reθt转捩模型与二维雷诺平均NS（RANS）方程数值求解相耦合,发展了考虑转捩影响的翼型气动特性计算方法。得到的计算结果与没有考虑转捩影响的计算结果的比较表明,与实验更加吻合,计算精度得到了明显提高。