封闭腔内自然对流数值方法研究

摘要

在文献 [1]的基础上 ,将非定常流函数涡量方程的数值求解方法推广至非等距网格剖分 ,其中流函数一阶导数即速度项采用二阶精度公式 ,包含温度在内的离散方程组采用ADI迭代方法求得定常解 ,以封闭腔内自然对流为例 ,进行了不同瑞利数 (Ra)条件下数值试验 ,对Ra =10 6的计算进行了必要的处理 .计算结果表明 ,该数值方法推导简单 ,计算稳定 ,为采用K -ε模式计算封闭腔内层流到湍流的转捩打下基础 .