飞机旋转表面结冰数值模拟

在Shallow-Water结冰热力学模型的基础上, 搭建了一套适用于三维旋转表面的非稳态结冰模型, 采用Jacobi迭代法和Gauss-Seidel迭代法对表面结冰进行非稳态迭代求解。用所提模型对简化后的旋转桨叶模型进行计算, 并与FENSAP软件结果进行对比, 验证了所提模型的准确性, 分析了转速、水滴直径和液态水含量等因素对旋转表面结冰冰形和水膜流动的影响。结果表明:随着转速的增加, 结冰范围和水膜覆盖范围偏移愈加明显；结冰范围和水膜覆盖范围随水滴直径增长逐渐增加, 水膜厚度也逐渐增大；结冰厚度和水膜厚度随液态水含量增长而相应增加, 水膜覆盖范围也明显变大。 

Numerical simulation of icing on aircraft rotating surfaces

Based on Shallow-Water icing thermodynamic theory, an unsteady icing model is developed to simulate icing on three-dimensional rotating surfaces. Jacobi and Gauss-Seidel iterative algorithms are adopted to numerically solve the governing equations of unsteady processes of surface icing. This method is used to calculate the simplified rotating blade model, and the results are compared with those of FENSAP, verifying the accuracy of the model. Then the effect of the rotation speed, droplet diameter and liquid water content on the rotating surface water film flow and icing shape is investigated. Results show that as the rotation speed increases, the deviation of the icing range and the water film coverage becomes more obvious. With the increase of the droplet diameter, the icing range and the water film coverage and thickness gradually increase. The icing and water film thicknesses increase as a result of the increase of the liquid water content, and the water film coverage also increases significantly.