控制方程:
x和y的稳态动量,连续性,能量方程方程相对于 Pr 和 Ra 数是无量纲的。 有关等式的详细说明,请参阅“ Ouertatani等人,封闭环境中二维Rayleigh-Bénard对流的数值模拟。CRMecanique 336(2008)”。 最终结果也针对本文进行了验证。
边界条件(二维正方形):
u=v=0 在所有四个边无量纲温度 = 顶部 -0.5,底部 +0.5 域左侧和右侧的零温度通量
数值方法:
用于解决速度-压力耦合的SIMPLE算法。 请查看“Versteeg, Malalasekera:计算流体动力学介绍”教科书。 速度网格与压力网格错开。 温度与压力具有相同的网格。 在每次迭代中使用雅可比方法更新速度和温度,并使用五对角矩阵算法直接求解压力校正方程。 应选择适当的欠松弛因子进行收敛。 Jacobi 方法是效率最低的迭代方法,但它是可并行化的。 鼓励您
2021-12-22 17:45:28
68KB
matlab
1