基于Myers模型的三维结冰数值仿真

三维结冰表面上的水膜流动和结冰增长是结冰计算模型应考虑的核心内容，其中广泛应用的是Myers模型。Myers模型考虑了空气剪切力和空气压力对结冰表面水膜流动的影响，以及冰层、水膜和空气之间的导热与对流传热对结冰速率的影响。本文在使用Myers模型进行结冰预测时，发现Myers模型对霜冰转化为明冰的判断标准存在缺陷，会在结冰极限处产生不合理的冰角。因此对Myers模型的结冰类型判断标准进行了修改，对机翼表面的结冰过程进行了更加准确的模拟，并应用了有效的离散算法计算水膜流动和结冰增长过程。对比了二维NACA0012翼型的单步法、多步法计算结果和实验结果。明冰结冰温度较低时，本文计算结果与实验结果吻合很好，明冰结冰温度较高时，本文对上冰角的计算与实验结果有一定差距。本文提供了三维GLC-305后掠翼的结冰计算结果和实验结果的对比，冰角厚度的计算结果略小于实验结果，但整体趋势一致。

Three-dimensional numerical simulation of icing using Myers model

Calculating the flow of thin water film and the accretion of ice layer is important for predicting of aircraft icing. The Myers model takes into consideration the water flow driven by air shear and pressure, and the ice accretion rate influenced by heat conduction and convection among ice layer, water film and ambient air. When the Myers model is used to predict aircraft icing, undesired ice horn will occur as the model's criterion for determining the ice type is simple. To solve this problem, the criterion for ice type determination is revised, and an effective numerical method is presented to solve the equation for water flow and ice accretion. The single-step and multi-step computational results of the two-dimensional NACA0012 airfoil are presented and compared with the experimental results. For glaze ice accretion at lower temperature, the computational results agree well with the experimental results. For glaze ice accretion at higher temperature, the computational results of the upper ice horn are slightly different from the experimental results. The computational results of ice accretion on the three-dimensional GLC-305 airfoil are also presented and compared with the experimental results. The computational results of the ice horn thickness is slightly less than the experimental results, but the computational results of the overall ice accretion tendency agree well the experimental results.