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

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

γ-Re_(θt)-CF转捩模型在Spalart-Allmaras湍流模型中的推广及验证

横流(CF)不稳定性是三维流动中诱发转捩的一项非常重要的因素,考虑到γ-Reθt-CF转捩模型对流向Tollmien-Schlichting波和横流波不稳定性引起转捩的判定均是完全基于当地变量,且Spalart-Allmaras(SA)湍流模型计算效率高,因而将γ-Reθt-CF转捩模型与SA湍流模型相结合,并将其引入开源Standford University Unstructured(SU2)计算流体力学分析平台。为了考察和验证模型的预测精度,分别使用原始γ-Reθt模型和γ-Reθt-CF-SA模型,对NLF(2)-0415后掠翼型和标准6∶1椭球模型进行了转捩预测数值模拟。算例结果表明,γ-Reθt-CF-SA模型的计算结果与试验数据吻合程度远远优于原始γ-Reθt模型,γ-Reθt-CF-SA模型能正确地预测出三维流动中的横流不稳定性引起转捩的现象。     

γ-Re_θ模式应用于高速边界层转捩的研究

应用γ-Reθ转捩模式对超声速、高超声速边界层转捩进行数值模拟。γ-Reθ模式通过求解关于当地雷诺数和间歇因子两个输运方程给出转捩起始位置和转捩区长度等信息。本文对马赫数3.5至7范围内的四种算例进行计算,研究了来流雷诺数、攻角变化和头部钝化半径等关键参数的变化对γ-Reθ模式预测流动转捩性能的影响。给出了壁面摩擦阻力系数、热流值与实验值的对比以及壁面附近间歇因子等值线分布等计算结果。γ-Reθ模式能够正确预测攻角和钝化半径变化时转捩位置和流动参数的变化趋势;在较大雷诺数时计算结果与实验值吻合很好。γ-Reθ模式对于超声速或高超声速边界层转捩的模拟仍需修正和改进。     

转捩模型在翼型气动性能预测 中的应用研究

边界层转捩位置的准确预测对于提高飞行器气动性能的预测精度具有重要意义.选取与k-ω SST湍流模型相耦合的γ-Reθt模型,以零压力梯度平板为研究对象,通过求解基于有限体积法的雷诺平均N-S方程验证该模型自动捕捉流动转捩的准确性;将该模型应用于传统有压力梯度的NACA 0012翼型的流场特性和气动性能的研究中,并与原始k-ω SST模型的计算结果及全湍流试验数据进行比较.结果表明:远场边界距离对翼型阻力系数有较大的影响;与无转捩模型相比,γ-Reθt转捩模型对翼型阻力系数的预测精度有一定程度的提高;对于二维模型,γ-Reθt转捩模型难以捕捉翼型表面的三维效应和非定常分离特性.     

Gamma-Theta经验转捩模型在DES中的应用

与RANS湍流模型相比,DES湍流模型的边界层外自由流部分存在LES区,该区域的涡粘性会有所降低。为了研究了这一特征对γ-Reθt转捩模型的影响,以平板实验为算例,选取基于k-w SST湍流模型的DES模型和DDES模型,检验γ-Reθt转捩模型在DES模型和DDES模型中的表现。结果表明：DES模型和DDES模型使转捩位置推迟,SST—DES模型由于对网格过度敏感而没有触发转捩,SST—DDES模型使转捩位置后移较少;DDES模型比DES模型更适合γ-Reθt转捩模型;转捩经验关系需要进行修正后才能将γ-Reθt应用于DDES模型。     

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

γ-Re_θ转捩模型的标定与应用

简单介绍了γ-Reθ转捩模型,并在mbns3d内部软件上实现。应用零压力梯度平板试验数据对关联函数进行了标定,得到了一组关联函数。对A-aerofoil迎角13.1°及13.3°进行了验证计算,结果表明该模型可以比较准确模拟转捩效应,提高了升力、阻力计算精度。随后研究了NLR7301两段翼型,通过应用环量修正及转捩模型提高了阻力计算精度,阻力最大误差控制到4%以内。     

基于新的经验关联公式的γ-Re_θ转捩模型在高超声速流动中的应用

为提高γ-Reθ转捩模型在高超声速流动下的转捩预测精度,参考部分高超声速风洞的实验数据对该模型进行改进,在原始经验关联公式的基础上提出了转捩动量厚度雷诺数基于马赫数、湍流度、压力梯度因子的经验关联公式;选取高超声速平板绕流、高超声速圆锥绕流等算例对提出的经验关联公式进行了验证,并将原始经验关联公式和提出的经验关联公式计算得到的重要参数进行对比,分析了提出的公式的作用机理。结果表明:基于低速不可压流动经验关联公式的γ-Reθ转捩模型在高超声速流动中计算得到的流场中Reθt分布偏低(一般只有真实值的40%~60%),进而导致边界层内Reθt偏低而过早达到转捩触发准则,与实验值相比所预测的位置可能提前10%(以参考长度为基准)甚至更多。新的经验关联公式有效地改善了这一问题,计算结果与实验值符合得更好。     

不同襟翼偏角梯形翼构型气动特性数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用二阶离散精度的单调迎风格式(MUSCL),结合剪切应力输运(SST)两方程湍流模型和γ-Reθ转捩模型,研究了梯形翼高升力构型襟翼偏角变化对气动特性的影响。主要目的是进一步确认Trisonic Platform(TRIP)软件模拟高升力梯形翼不同襟翼偏角引起微小气动特性变化的能力。首先,简要介绍了采用的计算方法;其次,介绍了两种襟翼偏角的梯形翼模型及风洞试验;最后,在网格收敛性研究的基础上,采用全湍流和转捩两种方式模拟了梯形翼构型不同襟翼偏角对气动特性的影响。与试验数据的对比结果表明,采用全湍流和转捩两种方式均可以较好地模拟不同襟翼偏角对气动特性的影响量,采用γ-Reθ转捩模拟方式可以提高梯形翼构型气动特性的模拟精度,失速迎角附近的气动特性模拟需要进一步研究。     

[γ-Reθt]转捩模型在跨声速涡轮叶栅中的应用

为了解涡轮叶栅在跨声速条件下的流动特性和准确预测涡轮叶栅外换热情况,对γ-Reθt转捩模型,依据零压力梯度时自由流湍流度衰减实验结果,给出了一种直接估算来流粘性比的方法,以保证叶片前缘附近具有正确的自由流湍流度分布,提高换热预测准确度;同时减少试算次数。对MARK II与VKI两种叶栅通道跨声速工况下的流动换热情况使用CFX软件,选取层流模型、SST k-ω模型以及缺省粘性比和设定合理粘性比的γ-Reθt转捩模型进行了数值模拟验证,计算结果与实验数据的对比表明:转捩模型优于其他模型;而采用本文方法给定进口粘性比,能准确预测转捩位置,同时显著改善γ-Reθt转捩模型对不同来流湍流度下涡轮叶栅表面换热的预测精度;当入口湍流度较高,相比采用缺省粘性比情况,压力面上换热系数的相对误差降低30%以上,控制在7%左右。     

模型的转捩流动计算分析简

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

Conjugate heat transfer investigations of turbine vane based on transition models

The accurate simulation of boundary layer transition process plays a very important role in the prediction of turbine blade temperature field. Based on the Abu-Ghannam and Shaw (AGS) and c-Re h transition models, a 3D conjugate heat transfer solver is developed, where the fluid domain is discretized by multi-block structured grids, and the solid domain is discretized by unstructured grids. At the unmatched fluid/solid interface, the shape function interpolation method is adopted to ensure the conservation of the interfacial heat flux. Then the shear stress transport (SST) model, SST & AGS model and SST & c-Re h model are used to investigate the flow and heat transfer characteristics of Mark II turbine vane. The results indicate that compared with the full turbulence model (SST model), the transition models could improve the prediction accuracy of temperature and heat transfer coefficient at the laminar zone near the blade leading edge. Compared with the AGS transition model, the c-Re h model could predict the transition onset location induced by shock/boundary layer interaction more accurately, and the prediction accuracy of temperature field could be greatly improved.     

Robust design of NLF airfoils

 A robust optimization design approach of natural laminar airfoils is developed in this paper. First, the non-uniform rational B-splines (NURBS) free form deformation method based on NURBS basis function is introduced to the airfoil parameterization. Second, aerodynamic characteristics are evaluated by solving Navier-Stokes equations, and the γ-Re_θt transition model coupling with shear-stress transport (SST) turbulent model is introduced to simulate boundary layer transition. A numerical simulation of transition flow around NLF0416 airfoil is conducted to test the code. The comparison between numerical simulation results and wind tunnel test data approves the validity and applicability of the present transition model. Third, the optimization system is set up, which uses the separated particle swarm optimization (SPSO) as search algorithm and combines the Kriging models as surrogate model during optimization. The system is applied to carry out robust design about the uncertainty of lift coefficient and Mach number for NASA NLF-0115 airfoil. The data of optimized airfoil aerodynamic characteristics indicates that the optimized airfoil can maintain laminar flow stably in an uncertain range and has a wider range of low drag.     

γ-Reθ转捩模型在二维低速问题中的应用

介绍了基于当地关联的γ-Reθ转捩模型.通过Spalart建议的环境源项法对自由流湍流度衰减的有效控制,实现了该转捩模型的实用化.利用中国空气动力研究与发展中心的TRIP软件平台上实现的该转捩模型,对S809翼型、NLR7301两段翼型和30P30N三段翼型3个低速翼型绕流问题进行了计算,将不同迎角下的计算结果与试验结...     

基于改进CST参数化方法和转捩模型的翼型优化设计

为提高翼型优化设计效率,增大设计空间,采用B样条基函数替代传统的形状类别函数(CST)方法中的Bezier多项式,增强了对翼型参数化表达的局部控制能力并提高了翼型局部表达精度。为了确保翼型在优化设计过程中的几何光顺特性和代理模型的准确性,采用小波分解技术提出了多分辨率翼型的局部光顺方法。采用基于本征正交分解(POD)的流场数值代理模型,并结合γ-Reθt转捩模型实现了快速准确的气动力与流动转捩预测。采用小波技术光顺的CST翼型参数化建模、POD流场数值计算代理模型以及γ-Reθt转捩模型,结合遗传算法建立了完整的翼型气动优化设计系统。针对低速层流翼型与超临界翼型进行优化设计,优化设计后的翼型升阻比分别提高了47.42%和45.85%,且对改进前后参数化建模方法的优化性能进行了对比,结果表明本文构建的翼型气动优化设计系统具备很高的优化效率。     

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

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

Quantification of parametric uncertainty in γ-Reθ model for typical flat plate and airfoil transitional flows

Owing to the lack of physical knowledge of boundary layer transition, the γ-Re_θ transition model introduces closure parameters, which increase the uncertainty of transition prediction. The objective of this work is to quantify the uncertainties of closure parameters in the quantities of interests and identify the key parameters. The six closure parameters in the uncertainty intervals are used as input variables, and the uncertainties of the output results are propagated by a stochastic expansion based on the point-collocation nonintrusive polynomial chaos method. The relative contribution of each parameter to uncertainty is evaluated by the Sobol index. The computational cases include natural and bypass transitional flows on zero-pressure-gradient flat plates, and subsonic and transonic flows around airfoils. For most cases, c_e2, c_a2, and c_a1 dominate the uncertainty, and the influence of σ_θt is also significant when the history effects of flow are evident. The contribution of parameters in airfoils is more complex than that in flat plates. The transonic airfoil case shows that flow separation dramatically changes the distribution of Sobol indices, which poses a challenge to the accurate prediction of transition. Generally, c_e2 and c_a2 are the key parameters of the γ-Re_θ model.     

利用遗传算法和LMI设计固定结构H2/H∞飞行控制律

 在设计飞翼式无人机(UAV)的横航向飞行控制系统时,为了使无人机具有较好的动态特性和阵风抑制能力,同时又便于工程实现,提出了固定结构的H₂/H_∞控制律设计方法。对于由此遇到的双线性矩阵不等式(BMI)问题,先用线性矩阵不等式(LMI)方法得到控制律参数和H₂/H_∞性能指标的映射关系,再用此映射关系作为适应度函数,用改进的遗传算法求解使H₂/H_∞性能最优的控制律参数。仿真结果表明,使用固定结构的H₂/H_∞控制方法的无人机动态响应迅速平滑,在侧风干扰下的滚转角振荡幅值仅是原经典控制方法的一半。     

基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.