论文部分内容阅读
标准的格子Boltzmann方法(LBM)在计算规则区域内的流体流动和传热问题时有一定的优势,然而实际问题中流场区域往往具有不规则的复杂边界,在这类问题中直角坐标下的LBM方法不仅无法准确描述计算域边界,而且计算效率也不高,因此将贴体坐标引入到标准的LBM热模型和多相流模型中具有重要的现实意义。考虑到通用的插值格子Boltzmann方法(GILBM)在不规则区域中对流体流动计算的优势,本文将标准的LBM热模型与之相结合,建立了贴体坐标下的LBM热模型。为了验证改进热模型的适用性,首先用其计算了圆柱绕流问题,计算所得的流场流线分布、圆柱表面压力系数及平均努赛尔数都与文献中的结果吻合良好。而后又进一步计算了水平圆环中的自然对流,当Ra=5×104时,计算所得环形空间中的等温线分布与Kuehn等人实验拍摄的结果高度吻合,同样环形空间内外壁面上的局部当量导热系数的数值计算结果与实验数据的对比显示：除了在环形空间底部的计算误差较大,其他各处均吻合良好,平均相对误差在3%以下。以上对比结果证明了贴体坐标下的LBM热模型在具有大曲率复杂边界流场中应用的可行性。鉴于贴体坐标下的LBM热模型在非等温单相流体系统中计算的成功,本文又尝试将GILBM方法与Dong等提出的直角坐标下的复合LBM相变模型相结合,建立了相应贴体坐标下的复合LBM相变模型。在验证模型的可行性方面,分别模拟了无重力条件下气泡在过热液体中的生长及水平过热壁面上的生长跃离,所得结果分别与Mikic的解析解以及Mukherjee的实验结果相吻合。然后进一步用改进的相变模型研究了不同Ja数下单气泡沿过热曲面的运动特性,发现气泡上升过程经历了沿壁面滑移和脱离壁面两个阶段,且Ja数越大气泡生长越快,上升速度变化也越快。曲面导致的气泡形变会使得气泡速度变化产生波动,波动频率和幅度与Ja数成正相关。气泡并非一直生长变大,其在脱离壁面后,由于受到周围过冷液体冷凝开始变小。气泡对温度场的扰动程度也与Ja数成正相关,同时温度场的剧烈变化也反过来影响气泡的生长。