热气防冰系统中短小销钉强化传热数值模拟

热气防冰是目前飞机常用的防冰技术之一,为了提高热气防冰系统的防冰效率,提出使用短小销钉绕流强化传热的方式.利用计算流体力学(CFD,Computational Fluid Dynamics)的方法研究了不同结构阵列销钉绕流对双蒙皮防冰通道的强化传热效果和流动阻力的影响,主要包括:销钉高径比、销钉阵列型式以及槽道高度.研究结果表明:加入销钉特别是叉排销钉能够强化传热,但也会引起流动阻力的增大;槽道高度小于1 mm时,传热效果的提高不明显,无需销钉来提高传热.研究结果为热气防冰系统双蒙皮通道的设计提供了一定的理论指导.     

直升机旋翼积冰的数值模拟

为研究直升机旋翼桨叶积冰问题,通过对积冰形成机理和直升机飞行原理的分析,建立了基于CFD(Computational Fluid Dynamics)技术的直升机旋翼桨叶结冰预测的理论模型;采用欧拉法求解水滴控制方程,在经典Messinger模型的基础上提出了积冰表面质能平衡的改进模型;以UH-60型直升机和CH-47型直升机为样机,具体分析了积冰条件下直升机桨叶冰形预测的方法,为进一步研究直升机结冰后的飞行动力学问题及防除冰系统的设计奠定了基础,对提高直升机飞行安全性具有积极意义.     

基于润滑理论的二维积冰数值模拟

将基于润滑理论的水膜-冰层模型推广到在任意二维截面上建立的正交曲线坐标系,得到描述二维型线表面水膜流动和积冰的偏微分控制方程.采用一种隐式-显式有限差分格式建立控制方程组的代数形式,给出求解非稳态控制方程组的数值方法.为验证模型和数值解法的有效性,在典型的航空积冰和传输线积冰环境下对翼型和传输线表面的积冰算例进行数值模拟.将水膜-冰层模型的冰形计算结果与传统Messinger模型的模拟结果以及冰风洞试验结果进行对比,低温明冰条件下采用当前方法计算的翼型表面冰形接近于Messinger模型的模拟冰形曲线;相对高温条件下的计算结果比传统Messinger模型更为精确.低速积冰环境下Messinger模型难以模拟的传输线表面积冰同样可以采用水膜-冰层模型进行有效 预测.     

机翼热气防冰系统设计

利用热流体系统仿真分析平台Flowmaster和计算流体力学软件Fluent进行联合仿真,对机翼热气防冰系统进行了设计和性能分析.设计参数包括防冰系统笛形管热气喷口孔径、孔数及孔间距等参数.基于Flowmaster对所设计的供气管路进行供气流量分配计算.基于Fluent模拟防冰腔内流动与换热,计算了防冰腔的热效率.根据防冰热载荷及热效率结果预测每段缝翼防冰所需热气流量,引入参数供气余度来分析防冰系统热性能.所提出的热气防冰系统设计方法可为热气防冰系统的工程设计与优化提供重要的帮助.     

基于弱耦合的翼伞气动变形数值模拟

对不带气室翼伞和带气室翼伞的气动特性和结构变形进行三维数值模拟,分析气室对翼伞气动特性和结构变形的影响.流场方面,采用有限体积法求解Navier-Stokes控制方程,选用剪切应力输运(SST,Shear-Stress Transport)k-ω湍流模型;结构方面,假定翼伞有初始形状,使用插值方法映射传递流固交界面的压力数据,通过ANSYS计算翼伞伞衣的气动变形.结果表明:气室对阻力影响不大,对升力的影响主要表现在大迎角情况下;翼伞柔性和尺寸大小对其有利迎角的范围影响不大;翼伞的气动变形和最大主应力主要集中在气室中前部分,带气室翼伞由于有肋片加强的缘故,伞衣气动变形较小,一定程度上保证了翼伞的气动特性.     

基于弱耦合的翼伞气动变形数值模拟

对不带气室翼伞和带气室翼伞的气动特性和结构变形进行三维数值模拟,分析气室对翼伞气动特性和结构变形的影响.流场方面,采用有限体积法求解Navier-Stokes控制方程,选用剪切应力输运(SST,Shear-Stress Transport)k-ω湍流模型;结构方面,假定翼伞有初始形状,使用插值方法映射传递流固交界面的压力数据,通过ANSYS计算翼伞伞衣的气动变形.结果表明:气室对阻力影响不大,对升力的影响主要表现在大迎角情况下;翼伞柔性和尺寸大小对其有利迎角的范围影响不大;翼伞的气动变形和最大主应力主要集中在气室中前部分,带气室翼伞由于有肋片加强的缘故,伞衣气动变形较小,一定程度上保证了翼伞的气动特性.     

基于CFD的水收集系数及防冰表面温度预测

利用计算流体力学(CFD,Computational Fluid Dynamics)软件FLUENT,对二维机翼在结冰气象条件下飞行时的水收集系数,以及防冰系统工作并达到稳定后的防冰表面温度进行了计算.水收集系数计算采用FLUENT中的欧拉两相流模型以及用户自定义函数.特定飞行条件及热气防冰系统结构下的防冰表面平衡温度分布预测采用防冰腔内外热流耦合的方法.所得结果和理论分析一致,表明了利用FLUENT可以有效的进行防冰研究,包括水收集系数及防冰表面温度的预测,为今后进一步研究设计热气或电热防冰系统打下基础.      

联翼布局俯仰力矩非线性变化特性的数值模拟

联翼布局飞机具有优良的升阻特性,是下一代亚声速飞机优先选择的气动布局型式之一,但在某些情况下其俯仰力矩随迎角的增大会表现出较强的非线性变化特点.针对该问题,在Ma=0.75条件下,采用数值模拟方法对某亚声速联翼布局气动性能及其绕流流场进行研究,通过对各部件气动特性分析,结合不同前翼绕流流动状态下前/后翼绕流场特点及截面气动力分布特点,揭示了前翼对后翼绕流流场干扰是引起其俯仰力矩非线性变化的主要原因.计算结果表明:在一定迎角下,该联翼布局飞机前翼绕流发生分离,从而影响后翼绕流流场,引起后翼气动效率下降,导致全机俯仰力矩随飞机迎角发生非线性上扬,对该机飞行性能的提高带来严重影响.     

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

积冰几何特性对翼型性能影响的神经网络预测

积冰几何形状对翼型气动系数的影响是复杂的.采用BP(Back Propagation)神经网络的LM(Levenberg-Marguardt)学习算法,建立明冰的典型几何特性(冰角前缘半径、冰角高度和冰角位置)对翼型气动系数影响的神经网络,得到该3种几何参数对气动系数影响的规律;建立了典型冰形参数对最大升力系数影响的神经网络,该网络能很好的预测冰形参数对应的最大升力系数值;此外,建立了冰型位置对舵面铰链力矩系数影响的神经网络.仿真结果表明,BP神经网络仿真结果与实验值具有高度一致性,并能预测非实验值条件下的气动系数;翼型表面积冰位置变化对气动系数影响最大;铰链力矩系数在失速迎角达到之前就发生突变,可以更安全地用来预测失速的发生.      

积冰对飞机飞行性能的影响

针对飞机飞行过程中出现的积冰所带来的不利影响,提出一种求解积冰条件下飞行性能的思路.基于二维有限基本解法对翼型前缘的积冰形状进行预测;在计算飞行性能时,利用数值模拟的方法求解RNG(Renormalization Group) k-ε湍流模型下的雷诺平均N-S(Navier-Stokes)方程,对飞机全机流场进行解算,从而得出积冰状态下飞机的升力、阻力特性.最终结果表明:积冰严重影响着飞机的飞行性能,这与由相关实验得出的结论一致.该方法可应用于工程领域,对积冰条件下飞机的飞行性能进行快捷、有效地求解,为研究飞机防冰技术、提高飞行性能提供参考.      

积冰对飞机纵向操稳特性的量化影响

提出了一种根据气象条件、飞机特征、飞行条件预测积冰对飞机飞行动力学特性影响的估算方法.采用热力学平衡分析方法完成了积冰增长过程的数学描述,利用逐步线性回归方法建立了结冰飞机的气动模型,计算得到了结冰飞机的气动参数,计算结果与试验结果吻合较好,表明这种方法具有一定的工程实用性.并在此基础上研究了不同部位积冰对飞机纵向操稳特性的影响,通过数值仿真计算得出,结冰飞机气动特性的急剧恶化会导致其纵向操稳特性显著降低,严重时会危及其飞行安全.      

基于CE/SE法的翼型流场激波捕捉计算

在二维四边形结构网格下,针对时空守恒元和解元(CE/SE, space-time Conservation Element and Solution Element)方法捕捉激波需要双倍空间网格点的问题,重新对守恒元和解元进行了定义,计算点同时包含单元中心和网格节点,并推导得到了新的CE/SE方法的计算公式.以此为基础,结合非敏感克朗数(CNIS, Courant Number Insensitive Scheme)计算格式和当地时间步长方法,对激波翼型流场进行数值模拟,并与原方法及AGARD报告进行对比.结果表明,对原CE/SE方法的改进是有效的,可明显提高激波分辨能力,并且激波前后无明显的数值震荡发生,新方法适合应用于翼型流场中捕捉激波.     

翼型波浪水面巡航地面效应数值模拟

数值求解非定常不可压缩流动的雷诺平均N-S(Navier-Stokes)方程和标准k-ε湍流模型,用VOR(Volume Of Fluid)方法模拟波浪水面,模拟了NACA2410翼型在波浪水面上方飞行的流场.研究了余弦波浪水面的生成方法,选取合适的计算网格和时间步长,避免了在生成规则余弦波浪水面时由于数值耗散使波面形状衰减.比较了固壁波浪与水面波浪计算结果的差异,发现固壁波浪的结果更接近余弦曲线分布.研究了波浪等级对翼型气动性能的影响,分析计算结果发现:在规则的余弦水面波浪上方飞行,翼型气动力呈现周期性,给出了一个周期内气动力的变化过程,以及波长和波幅对气动力平均值和波动幅度的影响规律.      

基于有限元仿真的航空整体结构件变形校正

针对航空整体结构件变形复杂、校正效果难以定量预测的工程问题,基于有限元仿真技术,综合研究加工变形—校正过程,对航空整体结构件的变形校正有限元建模、校正过程数值模拟、校正参数确定、校正载荷施加形式及校正残余应力等关键技术进行了研究,并运用上述技术对某航空铝合金材料整体结构件进行了变形校正有限元数值模拟.结果表明,综合研究加工变形—校正过程,提取变形有限元模型作为初始校正模型,进行有限元加载—卸载分析,可很好地模拟校正过程,并获得实际需要的校正参数.      

超声速场中的反向喷流数值模拟

采用了高分辨率的N-S(Navier-Stockes) 方程数值模拟方法,对钝头体头部反向喷流进行研究.计算并分析了头部反向喷流现象中喷流马赫数、来流马赫数、攻角等因素对流场细致结构的影响,并对反向喷流减少阻力和减少气动加热的原理进行了深入分析和探讨.研究分析表明这些因素通过决定弓形激波、马赫盘的强度和位置,来影响回流区的大小和位置,从而使有反向喷流时阻力系数所受到的影响比无喷流时更大.