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

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

Calculation of airfoil anti-icing/deicing characteristics based on water film flow and coupled heat transfer

Numerical simulation methods for airfoil anti-icing/deicing system based on water film flow and coupled heat transfer model were established to solve the coupled mass and heat transfer phenomena among air/water film/ice layer/airfoil on the airfoil surface in the aircraft anti-icing/deicing process. The numerical computation methods for runback flow water, ice accretion on the airfoil surface as well as inner temperature distribution under the effect of the thermal load of anti-icing/deicing system were established based on the Myers water film flow model. For the heat transfer process in the airfoil and ice layer, the enthalpy theory and the finite volume method were used to establish the numerical simulation method of heat transfer in the complex multilayer structure. For the phase change process in the ice layer, the phase change correction method based on enthalpy theory was proposed to consider the effect of phase change latent heat on temperature change. Finally, coupled calculation of the airfoil anti-icing/deicing process was realized. The results showed that through the coupled calculation of the water film flow and airfoil/ice layer heat transfer model with combination of the heat transfer boundary conditions at different interfaces and the enthalpy theory considering the latent heat effect of phase change, the temperature distribution can be accurately computed in the airfoil/ice layer, meanwhile the effective prediction and analysis of the flow and ice accretion characteristics of the runback water during the airfoil anti-icing/deicing process can be realized.