热气防冰参数对机翼表面溢流结冰的影响研究

机翼热气防冰数值模拟中,根据N-S方程求解流场,用Euler法获得水滴撞击特性,通过蒙皮导热将求解得到的内外流场进行耦合传热并达到稳定后,开始模拟结冰,来进行机翼热气防冰及形成溢流结冰的数值计算。计算结果表明,热空气防冰数值模拟是可行和合理的。采用数值模拟方法对机翼热空气防冰过程进行了模拟,得到了由于引气温度不足或机翼热空气保护面积不同而导致的不同溢冰高度和位置。分析了供气温度、防冰区域对翼面溢流冰形成的影响。结果表明：供气温度直接正向影响热交换后蒙皮表面温度。供气温度越低,溢流冰形高度越高,对气动特性影响也越大;热防护区域范围对溢流结冰结果也会产生影响,热防护区域越大,冻结位置距离驻点越远,而且冰形高度越低,对气动特性影响也越弱。     

基于内外传热耦合的热气防冰系统仿真计算

基于内外传热耦合原理,建立了热气防冰系统的AMESim仿真模型,研究热气防冰系统在不同引气状态下,管路流量、压力及蒙皮温度的变化规律.采用所建立的仿真方法,计算飞机机翼热气防冰系统的内部流动特性和内外耦合传热特性,将计算结果与热气防冰系统流量分配试验结果进行对比.结果表明:管路流量、压力和蒙皮温度的仿真计算结果与试验结果最大误差为9%,验证了该仿真模型的正确性.在此基础上,对飞机短舱热气防冰系统进行了仿真,分析了飞行包线下系统内外的换热、温度变化、加热效率等关键参数的瞬态特性.仿真结果为热气防冰系统的设计、分析与优化提供依据.      

基于流-固耦合传热的热气防冰系统干空气飞行蒙皮温度场计算研究

以机翼热气防冰系统为研究对象,建立了包含热气防冰系统防冰腔内外流场对流换热和固体结构导热的三维稳态流-固耦合传热物理模型,对整个计算区域生成混合网格,边界条件为第三类边界条件,采用计算流体力学方法以 FLUENT 软件为工具,对干空气飞行状态下流-固耦合传热模型进行了求解,获得防冰腔蒙皮内外表面对流换热系数分布和温度场结果,并对计算结果进行了分析。结果表明：防冰腔铝合金蒙皮沿展向和厚度方向导热显著,温度分布较均匀,防冰引气温度为200℃时,防冰腔蒙皮内外表面上最高温度为101℃,最低温度为21℃,3 mm厚的蒙皮同一点处内外表面最大温差仅为4℃,防冰腔排气口处气体的平均温度为63℃。热气防冰系统蒙皮温度场计算方法和计算结果,能够为热气防冰系统干空气飞行试验设计和测试中温度传感器的选型与布置提供依据。     

基于水膜流动与耦合传热的翼型防/除冰计算

针对飞行器防/除冰过程中翼面上空气-水膜-冰层-机翼之间的耦合传质传热现象,建立了一种基于水膜流动与耦合传热模型的翼型防/除冰数值模拟方法。基于Myers水膜流动模型建立了防/除冰热载荷作用下翼面溢流水流动、积冰及内部温度分布的数值计算理论。对于翼型及冰层内的传热现象,利用焓理论及有限体积法建立了复杂多层结构传热的数值模拟方法,对于冰层相变过程,提出了一种基于焓理论的相变修正方法以考虑相变潜热对温度变化的影响。最终实现了翼型防/除冰过程的耦合计算,结果表明:通过结合不同界面处的传热边界条件和考虑了相变潜热效应的焓理论对水膜流动与翼型/冰层传热模型进行耦合求解,能够对翼型/冰层内温度分布进行准确计算,可实现对翼型防/除冰过程中溢流水流动及积冰特性的有效预测与分析。      

飞机热气防冰系统与冰脊预测的数值模拟

基于流固耦合传热的思想建立了一套飞机热气防冰系统的的数值模拟方法,并将其与积冰热力学模型结合起来,实现了热气防冰系统开启时的机翼积冰预测.采用格心格式有限体积法求解N-S方程获得防冰腔与外流场;通过欧拉法在外流场的基础上获得过冷水滴撞击特性;求解三维热传导偏微分方程获得蒙皮的传热特性;采用交接面插值的方法实现防冰腔到外流场的热量传递;建立了考虑三维溢流效应的积冰热力学模型并在此基础上开展了机翼冰脊的数值预测.数值模拟结果表明:热气防冰系统开启时加热机翼表面温度最高可达308K,加热区后的上下机翼表面均有冰脊形成,通过对结果的分析表明该方法是合理可行的.      

短舱防冰系统三维内外流耦合计算方法

为提高防冰性能评估效率，更好地支撑发动机短舱防冰系统设计工作，基于松散耦合方法开展了适用于短舱进气道笛形管防冰系统的三维内外流耦合仿真计算方法研究。将三维全尺寸短舱进气道计算模型分割为内、外流计算域，并以防冰腔外表面温度和换热系数作为计算域之间的交互数据，以实现内外流耦合迭代。根据干、湿空气条件下防冰表面流动传热特性，考虑湿空气条件下由过冷水滴撞击产生的质量与能量源项，同时引入短舱进气道周向非对称特性导致的三维溢流效应影响，总结了干、湿空气条件下短舱进气道防冰系统内外流松散耦合迭代策略。计算结果表明：在干空气条件下防冰内外流松散耦合迭代仅需3个轮次即可收敛；在湿空气条件下，受防冰表面溢流影响，内外流耦合迭代4个轮次后趋于收敛。     

液态水含量对防冰表面水膜流动换热的影响

为了研究液态水含量对防冰表面水膜流动换热的影响,基于机翼防冰表面水膜及空气相互作用机理,并考虑水膜表面传热传质过程,建立了水膜与空气的流动换热模型,得到溢流水膜及空气边界层流动换热的积分控制方程,通过对比文献试验结果验证了模型的准确性.在此基础上,比较了不同液态水含量条件下防冰表面水膜厚度及主要热流量的分布情况.结果表明:液态水含量对水膜沿表面厚度分布有明显影响,而对换热过程中各项热流的影响主要集中在水滴撞击区域,加热热流与散热热流随液态水含量的增加呈现相反的变化趋势.     

某型飞机发动机短舱热气防冰系统性能数值模拟

使用三维内外强固传热耦合方法计算校核发动机短舱热气防冰系统的性能,并分析发动机进气流量对蒙皮表面温度的影响.内、外部表面传热系数计算均采用纯三维的CFD方法,在内、外部网格数据交互时使用了距离加权反比插值法.通过计算获得发动机短舱的局部水收集系数、蒙皮表面温度的分布情况、各处溢流水量,并由此判定此防冰系统性能是否达到要求.分析表明此发动机短舱热气防冰系统符合防冰性能要求;当发动机进气流量增大时,蒙皮表面温度下降,且溢流水量增加.     

二维机翼防冰表面溢流水膜破裂数值模拟

为了准确地预测机翼防冰系统热载荷,对二维机翼防冰表面溢流水流动状态进行了仿真计算.依据剪切力驱动模型,分别推导了水膜破裂前后的速度分布;引入能量最低原理,推导了溪状流模型,运用数值方法对水膜破裂的临界厚度进行了求解;使用基于Fluent软件的用户自定义函数方法编程并加载计算模块,分别以水膜模型和溪状流模型计算得到了溢流水的流动特性.结果表明:在一定结冰条件下,引入溪状流模型后,溢流水的溢流范围在两个不同算例中分别增大了33%和17%,这会很大程度影响防冰系统热载荷的分布.      

三维内外热耦合计算热气防冰系统表面温度

提出了一种三维内外强固传热耦合计算热气防冰系统表面温度的方法并进行算例计算.外部表面传热系数计算采用附面层积分方法,内部热气流动与换热采用CFD(计算流体动力学)方法,利用FLU-ENT的用户自定义函数UDFs(user-defined functions)实现蒙皮外表面热载荷计算及热流边界条件的自动加载,使迭代自动推进.对表面温度及外部热载荷进行亚松弛来稳定迭代计算,用强固传热耦合迭代来加快计算速度.所研究的三维热气防冰系统表面温度计算新方法可应用于防冰系统的热性能验证及系统优化.      

溢流条件下飞机结冰过程的传热特性研究

基于液/固相变和液膜流动的基础理论,对飞机结冰过程的传热特性进行了研究,建立了溢流与液/固相变耦合的结冰传热模型,并采用该模型对来流各参数对冰层生长特性的影响进行了分析.研究结果表明,结冰表面液态水的溢流将促进液/固相变过程,并进而提高冰层生长速率.气动剪切力是影响溢流效应的主要因素,来流速度越高,气动剪切力越大,溢流效应也越明显;反之,则溢流效应越微弱.在溢流条件下,来流参数中来流温度和速度是影响冰层生长速率的主要原因.来流速度越高或温度越低,则冰层生长速率越大,反之则生长速率越小.比较而言,液态水含量和水收集系数的变化对冰层生长速率的影响相对较小.      

来流速度对防冰表面溢流水流动换热的影响

为研究来流速度对防冰表面溢流水流动形态及换热的影响,基于空气-水两层相互作用的质量、动量和能量守恒,建立防冰表面溢流水水膜流动换热及破裂的数学模型,分析了防冰表面溢流水在不同来流条件下的流动形态和表面换热情况.计算分析表明:来流速度增加时,防冰表面相同位置处的连续水膜厚度减小,水膜破裂位置随之延后;较高来流速度条件下,破裂处水膜厚度稍有增加,使得破裂后形成的溪流厚度和宽度增大;作为主要的表面散热项,连续水膜表面蒸发及对流换热热流均随来流速度的增加而增大.此外,由水膜破裂引起的表面溢流水流态变化对防冰表面蒸发热流有一定影响.     

Effects of flow parameters on thermal performance of an inner-liner anti-icing system with jets impingement heat transfer

The multiple jets impingement heat transfer is widely applied in the wing anti-icing system. It is challenging to apply the similarity criterion to carry out the anti-icing experiments due to the complex flow and heat transfer behavior. In the present study, the full-scale slat model is used to carry out anti-icing experimental researches in a 2 m × 3 m icing wind tunnel of China Aerodynamics Research and Development Center. The effects of icing parameters Liquid Water Content (LWC) and Median Volume Diameter (MVD) and hot air parameters (mass flow rate and temperature) on the thermal performance of an inner-liner anti-icing system with jets impingement heat transfer are studied. The effects of the experimental parameters are analyzed in detail by combining impingement and evaporation heat transfer mechanisms. The impingement hot air mass flow rate dramatically affects the heat transfer performance of the impingement stagnation region within the range of the experimental parameters. The temperature of impingement hot air and that of wing skin are approximately linear correlated. The experimental results show the effects of LWC and MVD on water film formation and runback ice accretion. The formation of water film is analyzed by an analytical method based on the wing skin temperature difference of dry and wet air conditions.     

高温下加热/冷却管内参与性介质的辐射与对流耦合换热研究

本文通过数值模拟方法研究了高温下加热/冷却管内层流充分发展流动的参与性介质热起始段的辐射与对流耦合换热。采用离散坐标法求解圆柱坐标系下的非线性积分微分辐射传递方程,采用有限差分法求解辐射与对流的耦合换热能量方程。通过计算,分析了不同壁温及入口温度下,介质导热系数、壁面发射率、介质散射反照率及散射相函数等因素对局部努谢尔数、壁面热流密度及有效加热/冷却长度的影响,并同时采用努谢尔数和壁面热流密度评价了耦合换热中辐射对对流换热的影响。      

冰脱落对电热除冰传热特性的影响研究

冰脱落对电热除冰传热特性的影响是电热除冰过程中的重要特征.采用有限体积法和混合网格技术,基于热焓方法将传统的相变传热模型应用到电热除冰计算中,研究了冰脱落对电热除冰传热特性的影响,主要研究了加热模式和加热单元间隔对除冰结构传热特性的影响.研究发现:无间隔除冰单元是有间隔除冰单元的特殊情况,其蒙皮表面温度历程在本质上是一致的.冰脱落后,无间隔除冰单元蒙皮表面温度开始剧烈上升然后逐渐缓慢上升.而对于有间隔除冰单元,随着间隔上加热功率的降低,温升过程逐渐缓和.单元间隔的加热功率降低到一定程度时,表面温度出现温降过程.当传递到蒙皮表面的热量和对流换热带走的热量差逐渐趋于稳定时,表面温度又开始逐渐稳定的上升.     

多影响因素作用下碳氢燃料跨临界过程换热恶化的数值研究

为深入理解多影响因素作用下碳氢燃料跨临界过程换热恶化的特性,基于开源计算软件OpenFOAM对超临界RP-3的流动换热过程进行数值模拟。采用广义对应状态法则对碳氢燃料替代模型的物性进行计算,湍流模型选用SST(shear stress transport)k-ω湍流模型。与实验数据比较,热流密度为300~400kW/m2内的计算壁温平均误差小于3%。研究分析了换热恶化机理,讨论压力、进口温度、热流密度与质量流量之比对RP-3换热特性的影响。结果表明:拟临界温度附近RP-3热物性的剧烈变化是强制对流下发生换热恶化的主要原因;提高压力、降低热流密度与质量流量之比或减小进口温度是避免流体在拟临界温度附近发生换热恶化的有效措施;提出了换热恶化预测关联式,为主动再生冷却技术提供参考。