摘要

采用欧拉-拉格朗日法计算空气、水滴两相流动,根据质量守恒和能量守恒计算结冰量和溢流水量,溢流水的流动方向及不同方向的流量分配由当地空气的速度矢量、网格边界的法向量计算得到。采用非稳态方法、非结构网格计算了NACA0012平直翼、截面为GLC-305的后掠翼的结冰,并与文献中的实验结果进行了对比,结果显示最大结冰厚度与实验吻合较好,预测的冰形与实验冰形发展趋势一致。计算结果表明,提出的三维结冰计算方法是有效、可行的,为三维结冰计算提供了一条可选的途径。