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

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

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

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

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

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

考虑水膜蒸发的三维明冰积冰数值研究

明冰积冰过程中冰层表面会覆盖一层薄水膜,水膜在冰层表面流动和传热,水膜蒸发会引起质量和能量传递,对积冰相变产生影响。为了研究水膜蒸发对积冰的影响,在已有三维积冰模型的基础上,考虑了明冰积冰过程中水膜的蒸发,建立了考虑水膜蒸发的三维积冰数学模型,并开发了相应的积冰计算程序。选取NASA提供的典型算例,模拟了明冰积冰条件下NACA0012翼型的积冰过程,并将计算结果与NASA试验和数值结果进行对比。结果表明水膜的平均蒸发量为1.5g/(m~2·s),蒸发吸收的热量占总传热的35%。与不考虑水膜蒸发的计算结果相比,本文得到的积冰量和冰形都与试验结果具有更好的一致性,说明三维明冰计算中考虑水膜蒸发是有必要和合理的。     

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

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

机翼电加热防冰表面内外传热的耦合计算

分析了机翼防冰表面的传热传质现象,对电加热防冰表面内外传热耦合计算进行了研究.采用附面层积分法对机翼表面传热系数进行计算.根据传热传质比拟得到了防冰表面蒸发量.建立表面微元的能量与质量守恒差分方程组,利用高斯-赛德尔迭代法求解,将外部防冰载荷、蒙皮弦向导热和内部加热热流三者耦合,采用亚松弛系数稳定迭代过程,最后得到平衡表面温度、溢流水和溢流结冰的范围.计算结果和文献中试验以及其他数值仿真结果比较,吻合较好,表明了模型及计算方法的正确性.      

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

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

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

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

某型飞机发动机短舱防冰系统设计计算

介绍了某型飞机发动机短舱防冰系统的结构形式和防冰原理,按照可利用热量应大于所需热量的原则,通过计算完全蒸发和湿状态下的热载荷,确定了发动机短舱防冰严酷状态及所需热流量.基于机翼水收集系数和数据相似原则,利用欧拉法两相流计算方法,确定了短舱防冰系统水收集系数.按照在湿状态下允许有一定结冰量存在的设计原则,根据发动机的吞冰能力,利用VSAERO/ICE工具对发动机短舱防冰系统进行了计算,计算结果表明,在严酷状态下,防冰系统能够满足允许结冰厚度要求.      

迎风面三维积冰过程中水膜流动的计算方法

对影响迎风面三维积冰过程的冰层表面水膜流动进行分析,建立描述其流动的数学模型,引入在流动前锋下游预设薄水膜的方法处理其在表面干区的流动推进过程,提出超过临界厚度的水膜被吹除的方法处理角状冰下游局部水膜被吹离表面的现象.通过对翼型-平板结构中翼型上的冰形对比和水膜流动规律分析,验证了该水膜流动数学模型的合理性和计算方法的可行性.研究表明:与Messinger模型相比,积冰模型中引入水膜流动后,可以较好地模拟三维明冰积冰,同时明冰极限的模拟结果也更接近试验值;冰层表面的水膜厚度约为10-5m,流动速度约为10-2m/s.      

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.     

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

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

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

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

旋转帽罩电加热防冰瞬态过程研究

基于欧拉两相流理论,对旋转帽罩水滴撞击特性进行了数值模拟,提出了旋转部件表面对流换热计算方法,并基于改进Messinger结冰模型开发了旋转帽罩电热防冰计算程序,对旋转帽罩瞬态防冰过程进行了数值模拟研究.结果表明:旋转运动对旋转帽罩表面对流换热起到增强作用,且在供给相同加热热流密度时转速越大,防冰表面温度越低;防冰系统启动阶段,旋转帽罩表面会发生结冰和冰脱落现象,考虑结冰过程后系统的响应时间缩短;当电加热功率相同时,周期电加热防冰方式更为节能;当电加热能耗相同时,周期电加热方式系统响应更快.     

某型飞机机翼防冰系统计算分析

主要针对某型飞机机翼的热气防冰系统计算分析,得到水滴直径变化对撞击极限的影响,飞行马赫数变化对机翼表面换热系数的影响,分析了不同飞行高度湿表面和干表面的温度分布.结果表明水滴撞击区随着水滴直径增加而增大;机翼表面的换热系数随飞行马赫数的增加而增加;在相同计算条件下,干表面温度比湿表面温度要高.对多个典型截面以及其在不同飞行状态的计算结果表明,在给定的计算条件下,4km及7km时防冰系统工作都是有效的,7km时表面部分位置湿表面温度低于0℃.      

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

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

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

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