机翼前缘结冰对大飞机纵向模态特性的影响

结冰导致飞机气动特性恶化，进而影响飞行品质。针对机翼前缘结冰条件下飞机全包线范围飞行品质评估问题，提出了基于自适应拟配初值的等效系统拟配方法，构建了机翼前缘结冰构型，并通过数值模拟得到背景飞机机翼前缘结冰气动数据。建立了飞机系统模型，进行升降舵倍脉冲操纵，进而得到飞机的响应数据。分析数据特征计算了自适应拟配初值，在飞机全包线范围内干净构型及不同结冰严重程度条件下进行了等效系统拟配，获得了纵向短周期飞行模态特性参数，进而得出相应的飞行品质等级。仿真结果表明：结冰会对飞机模态特性造成影响，使得飞行品质降低，严重时可能导致飞行品质发生降级。     

机翼展向不同部位结冰对飞机气动力特性影响研究

机翼展向不同部位结冰对飞机气动力特性的影响规律是机翼防除冰系统设计需要考虑的重要因素之一。通过风洞试验方法,将机翼不同部位的模拟冰型加装在飞机模型上进行常规测力试验,研究巡航构型、着陆构型下的机翼展向不同部位结冰后的升力特性、阻力特性、俯仰力矩特性的变化规律。同时通过数值计算的手段,分析机翼不同部位结冰的流场细节特征。研究结果表明,机翼中段结冰对飞机气动力特性影响最为严重,翼根和翼尖结冰影响较小,研究结果为制定高效合理的防除冰系统提供技术依据。     

盒式机翼布局气动特性研究

采用数值模拟的方法研究了盒式机翼布局的气动特性。在数值模拟过程中,控制方程为雷诺平均N-S方程,采用基于非结构网格的有限体积方法进行求解。空间离散采用了Osher逆风格式,时间离散则采用五步Runge-Kutta格式。紊流模型采用了S-A方程模型。对盒式机翼布局外形在不同状态下的流场进行了数值模拟,给出了每个状态下的气动参数,并进行了分析,结果显示盒式机翼布局具有较好的气动特性。     

基于欧拉方法的旋翼翼型结冰数值模拟及参数影响研究

从研究二维翼型结冰问题入手,基于水滴流场的概念,采用基于气液两相流的欧拉方法求解水滴流场,得到水滴收集特性,再结合热力学模型数值模拟旋翼翼型结冰过程.通过计算结果和冰风洞试验数据的比较,验证了该方法正确有效.同时研究了气象参数对结冰的影响.基于本文建立的方法,计算分析了飞行速度、平均水滴直径、水滴浓度、温度等不同参数对水滴收集效率、冰型、结冰量、结冰范围的影响,深入地探讨了翼型结冰机理和结冰形成条件,得出气象参数影响特征.     

半模机翼振动对气动性能影响的风洞试验研究

为了验证翼型振动测压试验结果和有关振动数值模拟文献的计算结果,设计完成了对称翼型NA-CA0012和层流翼型NACA64-210两种半模机翼的低速风洞振动试验.通过专门研制的半模振动模型,选取5种激振方式,用直接测力法得到了模型在静态和不同激振方式下升力特性变化曲线,分别研究了振动频率、雷诺数、自然转捩和固定转捩、数据采集方式等参数对机翼气动特性的影响.试验结果表明:半模机翼振动的气动特性与二元翼型和数值计算情况有所不同,其影响因不同翼型构型、不同采集方式和翼面不同流动模式等会产生不同效果.并对产生原因进行初步探讨.     

水滴多尺寸分布对水滴撞击特性和结冰增长的影响

水滴直径是影响水滴撞击特性的重要参数,而撞击特性在很大程度上影响结冰冰形,因此,在结冰数值模拟中,如何合理选取水滴分布模型,对于保证计算精度有重要意义.本文基于求解雷诺时均N-S方程获得的流场结果,采用拉格朗日法计算水滴撞击特性,并采用Messinger模型进行结冰增长模拟,研究了水滴直径、水滴多尺寸分布对水滴撞击特性和结冰冰形的影响.结果表明,水滴的多尺寸分布对于结冰数值模拟结果具有重要影响,采用Langmuir D分布能够满足数值模拟精度要求.     

民用飞机适航用临界冰形的确定及验证

对民用飞机结冰适航取证用临界冰形确定及验证进行研究。在民用飞机结冰适航取证过程中,首先需要确定临界冰形,为结冰后的性能和操稳评估提供冰形输入。首先确定临界冰形的判断标准,制定临界冰形确定的工作思路;然后采用SADRICE结冰数值软件开展结冰冰形参数的敏感性分析,获取临界结冰条件;最后通过冰风洞试验对参数敏感性结论进行验证。结果表明数值模拟参数敏感性分析结论与冰风洞试验结果一致,数值模拟所得临界冰形上冰角与冰风洞试验一致,下冰角存在一定差异,但不影响冰形临界性判断结论。参数敏感性分析所得临界结冰条件以及冰风洞试验所得临界冰形正确,可作为结冰适航取证临界结冰条件及临界冰形。临界冰形确定思路、方法可为其他民用飞机结冰适航取证提供参考。     

内融冰盘管动态蓄冰特性数值模拟与实验研究

采用焓法建立了管外结冰过程的二维动态数学模型,数值模拟了冰层沿换热管径向及长度方向的动态生长规律,并在此基础上分析了载冷剂进口温度、流量变化和不同结冰时间等对冰层动态生长及蓄冷特性的影响,并进行了实验验证。实验结果表明:在结冰初期的2 h内冰层厚度及动态蓄冷量数值模拟结果与实验数据误差较大,而在蓄冰过程结束时二者吻合较好,误差在10%以内;沿载冷剂流动方向冰层厚度逐渐减小,数值模拟与实验结果基本吻合。数值模拟结果与实验数据基本相符,验证了建立的数学模型的合理性。     

基于大型结冰风洞的航空发动机结冰与防冰试验技术

航空发动机结冰和防冰过程复杂,数值计算无法对其进行准确模拟,因此试验研究是发动机研究过程中必不可少的手段。我国首座大型结冰风洞已具备开展飞机翼段结冰试验的能力,有必要进一步发展航空发动机结冰与防冰试验技术,以满足下一步我国航空发动机型号设计与适航取证的需求。依托3 m×2 m大型结冰风洞,发展了进气模拟技术和热气供气技术,提出了一套结冰风洞试验流程及方法,并针对某型航空发动机进气部件开展了结冰风洞验证试验。结果表明:试验可真实模拟发动机内外流耦合和压气机引气防冰状态,且提出的结冰风洞试验流程及方法合理可行,实现了试验动态过程监测及进气道内流场压力测量,为下一步我国航空发动机结冰防护系统设计与安全适航符合性验证提供了技术支撑。     

基于预处理的低速小展弦比机翼数值模拟

以雷诺平均的Navier-Stokes 方程为基础, 应用非结构混合网格技术,基于有限体积的时间推进预处理方法,对小展弦比机翼的低速流动进行了数值模拟.其中对流项的离散应用了基于Roe的Riemann近似解迎风格式,时间离散应用了基于LU-SGS的隐式格式.为描述湍流流动,采用了Spalart-Allmaras一方程湍流模型.为了验证本文方法对低速小展弦比机翼粘性流动的数值分析效果, 对几个不同展弦比的机翼进行了数值模拟.本文的数值结果与文献的实验数据基本一致, 表明本文方法能够准确模拟小展弦比机翼在低速流动时机翼的展弦比对其气动特性的影响.     

Numerical Study on Ridge Ice Accretion and Its Effect Under Thermal Ice Protection

Under the condition of thermal anti-icing,the liquid water on the leading edge of the airfoil that would flow to the downstream non-protective zone will produce ridge ice,thus endangering flight safety.Based on the existing three-dimensional(3 D)icing model which considers the water film flow on the ice layer,an icing model with thermal boundary condition is introduced.With the boundary conditions of none anti-icing and thermal anti-icing,glaze ice accretion and ridge ice accretion are simulated on a simplified airfoil of unmanned aerial vehicle(UAV),and then the lift coefficient and drag coefficient are calculated and compared with the smooth airfoil under the same conditions.The results show that the lift-drag ratio obviously decreases after glaze ice occurred on the leading edge under the condition of none anti-icing;and that after setting the condition of anti-icing heat flux in the impingement area,the glaze ice on the leading edge becomes thinner and the ridge ice occurs in the non-protective zone,so the airfoil with this icing characteristic gets a lower lift-drag ratio.     

随积冰历程的机翼蒙皮载荷实验研究

积冰改变了翼型的气动外形和绕流流场,使得机翼气动载荷分布产生动态变化。蒙皮作为气动载荷的承受及传递对象,会在气动载荷的动态作用下产生不同的振动响应。以某大弯度翼型为研究对象,提取了典型积冰增长过程中尾缘上下蒙皮振动特征,采用载荷谱方法研究积冰全历程的蒙皮振动及流场变化特性,并分析了不同材质的蒙皮在结冰不同阶段的响应及结构稳定性。结果表明:不同材质蒙皮对冰致脱体涡及后缘分离涡具有不同的载荷感知特性;刚性蒙皮载荷谱及其能量相对集中,柔性蒙皮载荷相对分散;随着积冰增多、冰角增长,脱体涡主频逐渐前移,冰角表面小冰枝引起宽幅高频振动;前缘脱体涡与尾迹掺混造成翼型后缘绕流载荷能量增加。     

地面效应下的不同翼型亚声速气动特性分析

采用数值方法研究了亚声速地面效应条件下不同翼型的气动特性,进一步以Ma=0.5来流工况为例,研究了翼型参数和飞行高度对气动特性的影响。计算结果表明在Ma为0.5、迎角为6°的地效情况下,翼型弯度减小,更容易在翼型前缘产生激波阻力;翼型下翼面后缘弯度增大使得后缘压力更高,升力系数和低头力矩相应增大;随着飞行高度的减小,地效作用加强,翼型下翼面压力增大,下翼面的升力增量大于上翼面吸力损失,机翼升力系数和升阻比增加越来越显著。     

基于合成射流的旋翼翼型动态失速控制研究

针对直升机旋翼工作环境下来流速度和迎角(Angle of attack,AoA)耦合引起的动态失速问题,建立了基于合成射流的旋翼动态失速控制的数值模拟方法。采用运动嵌套网格方法,通过对翼型的平移和旋转实现变来流速度-变迎角的耦合。以积分形式的雷诺平均N-S方程为主控方程,空间离散使用Roe格式,时间离散为隐式LU-SGS方法,以OA209翼型为研究对象,在翼型上表面放置合成射流激振器,开展了射流位置、动量系数、无量纲频率以及偏角等参数对轻度失速、深度失速下翼型动态失速控制的研究。研究发现,轻度失速下,射流位置靠近气流分离点时(20%c附近,c为翼型弦长),对逆压梯度引起的轻度失速控制效果最佳。深度失速下气流分离点虽在5%c之前,但射流位于前缘分离泡后端(10%c附近)时控制效果较好。大迎角需要较大的动量系数才能有效控制。射流频率对涡结构的尺寸和数量会产生一定影响,能改变气动特性波动幅度。较小的射流偏角对轻度失速的控制更有效,而深度失速则需要较大的偏角。     

直升机新型旋翼翼型气动特性与布局分析

利用结构网格计算流体力学(Computational fluid dynamics ,CFD)的翼型气动特性分析方法开展旋翼翼型气动特性计算。通过RAE282，NACA0012，OA212，OA207等翼型压力分布、升力和阻力等特性计算结果与试验结果的对比分析，验证了计算方法的准确性，并进一步完成了HF系列旋翼翼型的气动特性计算。基于翼型的气动特性，采用时间步进自由尾迹的旋翼气动性能分析方法开展旋翼桨叶翼型的气动布局优化设计，对悬停和前飞条件下的旋翼开展计算分析，得到两种条件下的旋翼气动特性。而后通过本文建立的优化方法开展旋翼翼型布局优化设计。     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围.在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术.通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式.研究表明:在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值.     

Highly-Efficient Aerodynamic Optimal Design of Rotor Airfoil Using Viscous Adjoint Method

In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.     

以风轮气动性能为目标的风力机翼型优化设计

提出了一种优化风力机翼型气动特性的设计方法。在现有风力机翼型的基础上,通过改变其最大厚度位置至后缘部分的外形进行优化计算。采用N-S方程计算翼型的流场及气动特性,以提高风力机切向力同时减小其法向力为目的,用二维线性插值的方法进行多目标优化。以DU 93W 210翼型为例,针对风力机翼型气动性能本身的矛盾性,根据实际需求提出几种方案加以优化。结果表明,在不同的性能需求下,本文所提出的设计方法能取得较好的效果。     

电热除冰的热力耦合特性及其对冰层的影响研究

采用计算和实验相结合的方法,初步研究了电热除冰过程中的热力耦合特性及其对冰层的影响.在电加热条件下,耦合外部气动力载荷的作用,采用有限元方法计算了不同热流密度下表面冰层和蒙皮间界面法向和切向应力的分布,比较了加热/不加热条件下界面法向和切向应力分布的差别,研究了冰层最大主应力随热流密度的变化规律.研究发现,电加热条件下,在表面冰层融化前,热力耦合特性将造成冰层内部应力的显著增加,从而 造成冰层局部区域的破裂,加速冰层的破坏.同时,设计的原理性实验结果验证了热力耦合特性对冰层的破坏影响.研究结果对于电热除冰理论和除冰技术的发展有现实意义.     

前缘缝翼构型平尾直升机气动特性分析

使用数值模拟方法研究了前缘缝翼翼型气动特性,并将该前缘缝翼翼型应用于直升机平尾,分析了该直升机全机气动特性,最后通过风洞试验对数值模拟结果进行验证。结果表明,前缘缝翼构型平尾能有效改善直升机的纵向静稳定性。