基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

基于润滑理论的二维积冰数值模拟

将基于润滑理论的水膜-冰层模型推广到在任意二维截面上建立的正交曲线坐标系,得到描述二维型线表面水膜流动和积冰的偏微分控制方程.采用一种隐式-显式有限差分格式建立控制方程组的代数形式,给出求解非稳态控制方程组的数值方法.为验证模型和数值解法的有效性,在典型的航空积冰和传输线积冰环境下对翼型和传输线表面的积冰算例进行数值模拟.将水膜-冰层模型的冰形计算结果与传统Messinger模型的模拟结果以及冰风洞试验结果进行对比,低温明冰条件下采用当前方法计算的翼型表面冰形接近于Messinger模型的模拟冰形曲线;相对高温条件下的计算结果比传统Messinger模型更为精确.低速积冰环境下Messinger模型难以模拟的传输线表面积冰同样可以采用水膜-冰层模型进行有效 预测.