基于曲面重构的飞机三维结冰数值模拟计算

在飞机结冰数值模拟中,采用多时间步长法可以准确模拟实际结冰情况。首先给出了三维结冰数值模拟方法,通过求解N-S方程计算空气流场、采用欧拉法计算水滴轨迹、求解质量和能量守恒方程计算结冰增长。对结冰后的结果数据进行重组,建立新的拓扑关系,采用非均匀双三次B样条曲面重构算法重构冰形曲面,得到几何模型,实现飞机三维结冰的多时间步长数值模拟。比较了单时间步长和多时间步长的计算结果;并与冰风洞试验数据及LEWICE预测数据进行对比,分析了不同时间步长对结冰外形的影响,得到了表面水收集系数随着时间的变化规律。结果表明,随着时间步长的不断缩小,水收集系数在机翼上表面不断增加,影响冰形的增长。时间步长越小,计算结果越接近试验数据,验证了曲面重构算法的正确性,说明多时间步长在结冰数值模拟方法中的优越性。     

多时间步长结冰数值模拟方法研究

时间步长是影响飞机结冰数值模批精度的重要参数之一。本文通过求解雷诺时均N-S方程,湍流模型采用k-ε两方程模型,获得空气流场,求解水滴运动轨迹方程获得水滴撞击特性,基于Messinger热力学模型求解能量和质量守恒方程计算冰形。并采用扇形分区法,更新机翼前缘结冰区网格,保持网格拓扑结构不变,实现多时间步长结冰数值计算。比较了采用单时间步长法争多时间步长法对翼型表面结冰增长数值模拟的计算结果,并与冰风洞试验数据及LEWICE预测数据进行对比。在此基础上,计算分析了不同时间步长对部件表面的结冰冰形的影响。结果表明,只有采用多时间步长法进行飞机结冰数值模拟方是有效的,并且存在一个合适的时间步长,既满足计算精度要求,又能提高计算效率。