Boundary-layer transition prediction using a simplified correlation-based model

This paper describes a simplified transition model based on the recently developed correlation-based c ? Reht transition model. The transport equation of transition momentum thick-ness Reynolds number is eliminated for simplicity, and new transition length function and critical Reynolds number correlation are proposed. The new model is implemented into an in-house com-putational fluid dynamics (CFD) code and validated for low and high-speed flow cases, including the zero pressure flat plate, airfoils, hypersonic flat plate and double wedge. Comparisons between the simulation results and experimental data show that the boundary-layer transition phenomena can be reasonably illustrated by the new model, which gives rise to significant improvements over the fully laminar and fully turbulent results. Moreover, the new model has comparable features of accuracy and applicability when compared with the original c ? Reht model. In the meantime, the newly proposed model takes only one transport equation of intermittency factor and requires fewer correlations, which simplifies the original model greatly. Further studies, especially on separation-induced transition flows, are required for the improvement of the new model.     

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.     

Calibration of a γ-Re_θ transition model and its validation in low-speed flows with high-order numerical method

Based on the Reynolds-averaged Navier–Stokes(RANS)equations and structured grid technology,the calibration and validation of γ-Re_θ transition model is preformed with fifth-order weighted compact nonlinear scheme(WCNS),and the purpose of the present work is to improve the numerical accuracy for aerodynamic characteristics simulation of low-speed flow with transition model on the basis of high-order numerical method study.Firstly,the empirical correlation functions involved in the γ-Re_θ transition model are modified and calibrated with experimental data of turbulent flat plates.Then,the grid convergence is studied on NLR-7301 two-element airfoil with the modified empirical correlation.At last,the modified empirical correlation is validated with NLR-7301 two-element airfoil and high-lift trapezoidal wing from transition location,velocity profile in boundary layer,surface pressure coefficient and aerodynamic characteristics.The numerical results illustrate that the numerical accuracy of transition length and skin friction behind transition location are improved with modified empirical correlation function,and obviously increases the numerical accuracy of aerodynamic characteristics prediction for typical transport configurations in low-speed range.     

A novel hybrid method for aerodynamic noise prediction of high-lift devices

Aerodynamic noise of High-Lift Devices(HLDs) is one of the main sources of airframe noise, and has immediate impacts on the airworthiness certification, environmental protection and security of commercial aircraft. In this study, a novel hybrid method is proposed for the aerodynamic noise prediction of HLD. A negative Spalart-Allmaras(S-A) turbulence model based Improved Delayed Detached Eddy Simulation(IDDES) method coupling with AFT-2017b transition model is developed, in order to elaborately ...     

Gain self-scheduled H_∞ control for morphing aircraft in the wing transition process based on an LPV model

This article investigates gain self-scheduled H 1 robust control system design for a tailless fold- ing-wing morphing aircraft in the wing shape varying process. During the wing morphing phase, the aircraft’s dynamic response will be governed by time-varying aerodynamic forces and moments. Nonlinear dynamic equations of the morphing aircraft are linearized by using Jacobian linearization approach, and a linear parameter varying (LPV) model of the morphing aircraft in wing folding is obtained. A multi-loop controller for the morphing aircraft is formulated to guarantee stability for the wing shape transition process. The proposed controller uses a set of inner-loop gains to provide stability using classical techniques, whereas a gain self-scheduled H 1 outer-loop controller is devised to guarantee a specific level of robust stability and performance for the time-varying dynamics. The closed-loop simulations show that speed and altitude vary slightly during the whole wing folding process, and they converge rapidly after the process ends. This proves that the gain self-scheduled H 1 robust controller can guarantee a satisfactory dynamic performance for the morphing aircraft during the whole wing shape transition process. Finally, the flight control system’s robustness for the wing folding process is verified according to uncertainties of the aerodynamic parameters in the nonlinear model.     

Instability analysis and drag coefficient prediction on a swept RAE2822 wing with constant lift coefficient

Swept wing is widely used in civil aircraft,whose airfoil is chosen,designed and optimized to increase the cruise speed and decrease the drag coefficient.The parameters of swept wing,such as sweep angle and angle of attack,are determined according to the cruise lift coefficient requirement,and the drag coefficient is expected to be predicted accurately,which involves the instability characteristics and transition position of the flow.The pressure coefficient of the RAE2822 wing with given constant lift coefficient is obtained by solving the three-dimensional Navier-Stokes equation numerically,and then the mean flow is calculated by solving the boundary layer(BL) equation with spectral method.The cross-flow instability characteristic of boundary layer of swept wing in the windward and leeward is analyzed by linear stability theory(LST),and the transition position is predicted by eNmethod.The drag coefficient is numerically predicted by introducing a laminar/turbulent indicator.A simple approach to calculate the lift coefficient of swept wing is proposed.It is found that there is a quantitative relationship between the angle of attack and sweep angle when the lift coefficient keeps constant;when the angle of attack is small,the flow on the leeward of the wing is stable.when the angle of attack is larger than 3°,the flow becomes unstable quickly;with the increase of sweep angle or angle of attack the disturbance on the windward becomes more unstable,leading to the moving forward of the transition position to the leading edge of the wing;the drag coefficient has two significant jumping growth due to the successive occurrence of transition in the windward and the leeward;the optimal range of sweep angle for civil aircraft is suggested.     

A DISSIPATION-DISPERSION TURBULENCE THEORY AND ITS MODIFIED k-ε MODEL

Based on the perturbation equation of the Navier-Stokes equation and the instability analysis of strongly nonlinear flow. a set of new control equations is presented featuring the interaction between dissipation and dispersion of turbulence. Three extra terms repr(?)nting respectiveiy the effects of extra strains, rotation and dispersion on the turbulent stresses introduced into the momentum equation. The defects of the original k-(?) model are removed mainly by introducing new (?) terms. This modified k-(?) model has advantages of simplicity, self-consistency and clarity of physics It is abie to handle complex flows and provides new perspectives for looking at a diversity of turbulence problems, such as inter(?)ttency, coberent structure,anisotropy, turbulent energy inversion and the cut-off of turbulence spectra.     

Design and validation of a variable camber wing structure

Variable camber wing technology is one of the important development trends of green aviation at present.Through smooth,seamless,continuous and adaptive change of wing camber,the aerodynamic performance is improved in achieving increase in lift and reduction in resistance and noise.Based on the aerodynamic validation model CAE-AVM,Chinese Aeronautical Establishment(CAE) has carried out the design and validation of a variable camber wing,proposed an aerodynamic deformation matrix for the leading a...     

利用气动弹性的柔性振翼模型建立及分析(英文)

Making use of modal characteristics of the natural vibration of flexible structure to design the oscillating wing aircraft is proposed. A series of equations concerning the oscillating wing of flexible structures are derived. The kinetic equation for aerodynamic force coupled with elastic movement is set up, and relevant formulae are derived. The unsteady aerodynamic one in that formulae is revised. The design principle, design process and range of application of such oscillating wing analytical method are elaborated. A flexible structural oscillating wing model is set up, and relevant time response analysis and frequency response analysis are conducted. The analytical results indicate that adopting the new-type driving way for the oscillating wing will not have flutter problems and will be able to produce propulsive force. Furthermore, it will consume much less power than the fixed wing for generating the same lift.     

Transonic wing stall of a blended flying wing common research model based on DDES method

Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation(DDES) and unsteady Reynolds-averaged Navier-Stokes(URANS) equations methods based on the shear stress transport(SST) turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10~6. A joint time step/grid density study is performed based on power spectrum density(PSD) analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x_(cr)= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.     

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

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

梯形翼风洞试验模型数值模拟技术

基于雷诺平均Navier-Stokes (RANS)方程和结构网格技术,采用二阶空间离散精度的MUSCL格式,并结合k-ω剪切应力输运(SST)两方程湍流模型和γ-Re_θ转捩模型,研究了梯形翼风洞试验模型中前缘缝翼、后缘襟翼连接装置对气动特性的影响。简要介绍了本文采用的计算方法;介绍了梯形翼的风洞试验模型及风洞试验结果;在网格收敛性研究的基础上,采用"全湍流"方式和转捩模型研究了梯形翼试验模型连接装置对气动特性的影响。通过与不带连接装置的计算结果的对比,采用"全湍流"模拟方式,计算模型中考虑试验模型的连接装置引起升力系数下降、阻力系数下降、低头力矩减小以及失速迎角提前;通过与试验数据的对比,进一步考虑转捩影响可以提高梯形翼风洞试验模型气动特性的计算结果与试验结果的吻合程度,梯形翼风洞试验模型失速迎角附近的气动特性数值模拟技术还需要进一步的研究。     

开裂式方向舵在变前掠翼布局中的操纵性能研究

 针对变前掠翼(VFSW)无尾布局的横航向操纵,设计了开裂式方向舵(SR)操纵面,采用Navier-Stokes控制方程的有限体积离散方法以及剪切应力输运(SST)湍流模型,计算了变前掠翼中平直翼和典型的前掠翼布局开裂式方向舵的操纵性能,并对其操纵效率进行了比较,分析了其流场形态。计算结果表明:右侧开裂式方向舵打开后,平直翼时迎角对偏航力矩的影响较小,而舵偏量对偏航影响显著,利于偏航;前掠翼的偏航力矩随迎角的增加有所波动,但在小迎角时较为稳定,偏航作用随舵偏量增加而增强;右侧滑对滚转作用强于左侧滑,偏航作用低于左侧滑。经比较,平直翼的偏航作用明显强于前掠翼,平直翼和前掠翼的滚转和偏航作用均具有耦合性,但平直翼的耦合效应弱于前掠翼。     

基于边界层相似性解的放大因子输运模型

为了使二维线性稳定性理论能够适应现代CFD求解技术,通过求解Falkner-Skan边界层相似性方程获得各个形状因子下的相似速度型。基于线性稳定性理论,对各个速度型进行稳定性分析从而获得对应不同速度型的扰动放大因子包络线。最后使用标量输运方程的形式实现包络近似方法中放大因子的当地化求解,并结合原始γ-Re_θt转捩模型中的间歇因子输运方程,实现了自然转捩和分离泡转捩的建模。使用该输运模型对S&K平板、S809翼型、NLR7301翼型和DLR-F5机翼进行转捩预测,结果均与试验结果吻合较好,验证了该模型构建的合理性和可行性。     

低压涡轮叶栅流动中转捩模型的校验及改进

为评估并提高现有转捩模型的预测精度,采用计算流体力学(CFD)软件FLUENT 12.1,选取层流模型、全湍流模型、剪切应力输运(SST)低雷诺数模型、k-kl-ω模型以及γ-Reθ模型对低压涡轮叶栅T106-EIZ进行数值模拟,通过与实验数据的对比校验了后3种模型对于转捩以及相关参数的模拟能力,并对结果以及模型的作用机理进行分析,校验结果表明所有模型都不能准确地预测分离流转捩以及尾迹诱导转捩.选取预测效果较好的γ-Reθ模型进行了修正,提出通过修改间歇因子输运方程中的参数Ca1和Ca2的方法来修正该模型,结果表明该方法可以提高模拟精度.     

翼伞平面形状对翼伞气动性能的影响

对带气室的展弦比为3的不同平面形状翼伞模型的流场进行了三维定常数值模拟,详细考察了平面形状对翼伞气动性能的影响.运用有限体积法对三维坐标系下不可压雷诺时均Navier-Stokes(RANS)方程进行了直接求解,采用剪切应力输运(SST)k-ω二方程湍流模型进行湍流模拟.数值模拟得出的原始翼伞的气动性能参数与试验数据在...     

横流不稳定性转捩预测模型

由于Langtry和Menter提出的γ-Re_θt边界层转捩模型只能预测流向的边界层转捩现象,因此继续改进该转捩模型使其具有横流不稳定性转捩的预测能力显得非常必要。通过对经典Falkner-Skan-Cooke (FSC)三维边界层相似解的理论分析和数值求解,结合Thwaites压力梯度因子与当地后掠角构建的函数关系来求解复杂构型的当地Hartree压力梯度因子β_H以及形状因子H₁₂,采用由试验数据标定的C1准则关系式获得横流转捩位移厚度雷诺数,从而建立能够对复杂构型进行横流不稳定性转捩预测的转捩判据。应用所建模型对30°前缘后掠角的ONERA-M6机翼和变前缘后掠角的DLR-F5机翼以及标准6:1椭球标模进行了横流不稳定转捩数值模拟,计算结果显示转捩位置均与试验数据吻合较好,证明了该模型的合理性和实用性。     

基于放大因子与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%.      

大展弦比机翼模型设计对翼型流场气动特性的影响

采用SST两方程湍流模型,通过求解非定常Navier-Stokes方程,模拟了大展弦比机翼风洞模型振动条件下的翼型流场,总结了翼型不同振动状况下的流场和气动力特点,分析了模型设计中的不同振动情况对风洞试验结果的影响。研究结果表明：在大展弦比机翼风洞模型的设计中,将翼型的重心设计在机翼的弹性轴之后,对风洞试验的精度较为有利。此结论对大展弦比机翼的风洞实验模型设计有指导意义。     

巡飞器无舵偏情况下的气动数值分析

采用SSTk-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程实现对巡飞器无舵偏情况下湍流流场的数值模拟。对巡飞器低雷诺数绕流流场进行了计算,分析了其在无舵偏情况下气动特性随迎角、侧滑角的变化情况。仿真结果表明,巡飞器采用充气式机翼后,具有很好的失速性能与较好的大迎角稳定性。