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

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

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

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

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

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

提高k-ωSST模型对翼型失速特性的模拟能力

采用CFD软件fluent,k-ω SST模型对S825翼型进行了二维数值模拟.针对k-ω SST模型对翼型分离失速特性模拟不准的情况,对分离区域内湍流的强非平衡输运特性进行了分析.研究表明:原始模型由于不能准确模拟分离区内湍流的强非平衡输运特性,导致其对翼型失速特性模拟失效;提出了通过修正模型系数a₁及β_∞^*,提高k-ω SST模型预测非平衡湍流输运特性的方法,从而提高对翼型失速特性的模拟精度.     

结冰模拟与AEDC方法有效性的验证

介绍了AEDC(Arnold Engineering Development Center)结冰缩比方法,并对其有效性问题,用数值仿真方法进行了验证.利用结冰仿真计算软件对NACA0012翼型及其1/2和1/3缩比模型、圆柱及其1/2缩比模型进行了结冰数值模拟,内容涉及流场、水滴撞击特性和冰生长计算等.缩比部件的结冰条件由AEDC方法确定.结果表明,翼型、圆柱与其缩比模型的结冰过程具有较好的相似性.说明AEDC方法科学合理,其确定的结冰缩比试验参数有效,同时也说明结冰模拟技术可用于结冰缩比相似性研究.     

改进SA模型对翼型分离流动的数值模拟

采用(CFD,Computational Fluid Dynamics)软件Fluent,选取(SA,Spalart-Allmaras)模型,对风力机翼型S825进行了二维数值模拟,并与实验数据进行比较.针对SA模型不能准确预测翼型尾部流动分离的情况,分析了分离区域内湍流的强非平衡输运特性.研究表明,SA模型中系数C_b1与湍流的非平衡输运特性密切相关,进而提出了修正模型系数C_b1改进SA模型的方法,结果表明该方法能够更准确地模拟风力机翼型分离流动.     

直升机前飞状态旋翼结冰风洞试验研究

为研究直升机旋翼在不同前飞状态的结冰情况,研制了直升机旋翼模型结冰试验台,在中国空气动力研究与发展中心结冰风洞4.8 m×3.2 m试验段中进行了旋翼结冰风洞试验。以某型直升机2 m直径旋翼缩比模型为试验对象,分析了旋翼转速、初始拉力系数对旋翼结冰的影响。结冰试验过程中保持旋翼操纵恒定,利用天平测量了旋翼拉力和扭矩性能的动态变化,并采用二维冰形切割及三维冰形扫描的方式分别测量了桨叶展向典型剖面的翼型及桨叶的整体结冰形态。试验得到了旋翼结冰关键数据,分析结果表明:旋翼桨叶结冰主要集中在桨叶前缘和下表面,结冰会在降低旋翼升力的同时增大旋翼扭矩和功率;小拉杆杆端轴承的积冰可能造成卡塞,变距拉杆上的积冰可能造成杆端轴承卡塞,从而使旋翼操纵失效。      

过冷大水滴动力学行为对翼型结冰的影响分析

过冷大水滴（SLD）是极端危险的飞行环境之一。因大粒径水滴独特的动力学行为变形破碎、飞溅反弹,传统结冰计算方法难以准确地反映SLD结冰情况。采用Navier-Stokes方法求解流场、Euler方法计算水滴撞击、Shallow Water模型模拟结冰,并与NASA实验结果进行了对比验证方法可信。结果表明：SLD动力学行为对结冰和冰形影响较大。其中,变形破碎改变了水滴运动轨迹和撞击范围,降低了水滴撞击极限,导致上下结冰极限减小2.83%、2.13%;飞溅降低了驻点附近水滴收集率,导致前缘积冰量减少8.09%;反弹显著降低了水滴撞击极限,导致上下结冰极限减小30.69%、20.01%;撞击后飞溅反弹二次水滴再入流场使得上下结冰极限增加6.14%、3.71%。同时,与干净翼型相比带冰翼型空气动力学性能严重退化,在相同迎角下,升力更小、阻力更大、气动效率更低。     

跨音速翼型的气动最优化设计方法

将计算流体动力学(CFD)与最优化技术相结合,通过数值求解欧拉方程,对翼型绕流流场作出了数值模拟.再利用几何、流动和最优化控制方程,反复迭代求得在一定约束条件下气动性能最优的翼型.本文以NACA0012为原始翼型,选取两种设计工况,都取得了满意的结果.     

多段翼积冰的数值模拟及防冰预测

为研究飞行器翼面积冰预测问题,通过对积冰形成机理的分析,提出了基于经典Messinger模型的积冰表面传质传热的改进模型;通过引入对多段翼空气流场的插值计算,建立了多段翼翼面积冰预测模型;考虑防冰系统向翼面控制体带入的等效热功率,提出了防冰条件下的积冰预测方法.通过对算例的分析对比,数值模拟计算结果与实验数据基本吻合,表明本文提出的结冰计算模型合理,能够用来进行积冰预测分析,为进一步研究翼面结冰后的飞行动力学问题及防除冰系统的设计奠定了基础.      

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

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

机翼防冰过程中冰脊问题的数值分析

对机翼防冰过程中冰脊的形成特点及冰脊对机翼气动特性的影响进行了计算分析.基于经典的Messinger结冰模型开发了多步结冰程序,对不同条件下过冷水滴的撞击结冰进行了热质耦合计算,计算结果和文献试验结果吻合较好,表明该热质耦合算法的正确性.在此基础上对不同环境温度、飞行速度和不同加热功率等条件下冰脊的生长特点和机翼的气动特性进行了计算分析.结果表明:在非霜冰条件下,冰脊主要在热防护极限外并紧挨着热防护极限的位置处形成和发展,而在霜冰条件下,冰脊主要在机翼下表面形成,但是在热防护区域内有显著的结冰出现,该条件下的结冰对机翼的气动特性具有较大的破坏性.     

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

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

旋翼结冰模型与纵列式直升机平衡特性分析

为了研究纵列式直升机结冰后的飞行动力学问题,基于CFD(ComputationalFluid Dynam ics)技术提出了旋翼结冰的理论模型;基于旋翼冰风洞试验数据,考虑桨叶表面附着冰脱落及桨叶表面局部温度对结冰的影响,提出了旋翼结冰的工程模型.建立了纵列式直升机结冰后的飞行动力学模型,研究了不同飞行速度条件下结冰对样机CH-47B平衡特性的影响,分析了环境温度、结冰时间、液态水含量以及平均水滴直径对旋翼扭矩的影响.计算结果的对比分析显示:两种结冰模型在计算精度上基本一致,表明结冰模型合理,能够用来进行结冰后的飞行动力学分析,为纵列式直升机结冰后的飞行动力学研究奠定了基础.     

一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点，考虑到算法的质量（收敛速度、稳定性、精度等），结合磁流体数值计算的特性，对三维球坐标磁流体动力学（MHD）方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacComack格式，发展了一类快捷的具有TVD特性的组合数值新方法，作为格式的检验，在一维情况下，将其与PPM格式进行了比较，对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果，然后将其诮到定态太阳风的数值模拟上，在不同等离子体β情形下，可得到理想的太阳风定态结构，为今后将此数值模式应用到具有复杂磁场位型或三维直实太阳风暴的数值模拟研究奠定了基础。