单气泡沿过热曲面运动的格子Boltzmann方法模拟

发布时间：2017-07-14 11:02

  本文关键词：单气泡沿过热曲面运动的格子Boltzmann方法模拟

  更多相关文章： 复杂流场 贴体坐标 GILBM LBM相变模型 气泡生长运动

【摘要】：标准的格子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数成正相关,同时温度场的剧烈变化也反过来影响气泡的生长。
【关键词】：复杂流场 贴体坐标 GILBM LBM相变模型 气泡生长运动
【学位授予单位】：大连理工大学
【学位级别】：硕士
【学位授予年份】：2015
【分类号】：TK124
【目录】：

• 摘要4-5
• Abstract5-9
• 1 绪论9-15
• 1.1 研究背景及意义9-10
• 1.2 非等温系统中气泡动力学特性的研究现状10-12
• 1.3 LBM在复杂流场中的应用现状12-13
• 1.4 本文的主要工作13-15
• 2 格子Boltzmann方法的理论和模型15-32
• 2.1 格子Boltzmann方法的基本原理15-17
• 2.2 计算边界的处理17-20
• 2.2.1 周期边界处理方法17-18
• 2.2.2 反弹边界处理方法18-19
• 2.2.3 充分发展边界处理方法19
• 2.2.4 非平衡态外推处理方法19-20
• 2.3 热格子Boltzmann模型20-25
• 2.3.1 热模型简介20-22
• 2.3.2 封闭方腔自然对流22-25
• 2.4 格子Boltzmann多组分多相流模型25-31
• 2.4.1 多组分多相流模型简介25-26
• 2.4.2 Zheng模型26-28
• 2.4.3 Laplace定律验证28-31
• 2.5 本章小结31-32
• 3 贴体坐标下的LBM理论32-43
• 3.1 贴体坐标下的网格生成技术32-35
• 3.1.1 贴体坐标的基本概念32-33
• 3.1.2 贴体网格的生成方法33-35
• 3.2 贴体坐标下的格子Boltzmann方法35-41
• 3.2.1 通用的插值格子Boltzmann方法(GILBM)35-39
• 3.2.2 GILBM实施步骤39
• 3.2.3 45°斜方腔的顶盖驱动流39-41
• 3.3 本章小结41-43
• 4 单相复杂流场传热问题的数值计算43-59
• 4.1 贴体坐标下的热格子Boltzmann模型43-44
• 4.2 复杂流场传热问题的计算与模型验证44-58
• 4.2.1 圆柱绕流44-50
• 4.2.2 水平同心环形空间自然对流50-55
• 4.2.3 水平偏心环形空间自然对流55-58
• 4.3 本章小结58-59
• 5 过热曲面上气泡运动特性的数值研究59-73
• 5.1 贴体坐标下的复合LBM相变模型59-63
• 5.1.1 直角坐标下的复合LBM相变模型59-61
• 5.1.2 贴体坐标下的复合LBM相变模型的建立61-63
• 5.2 贴体坐标下的复合LBM相变模型的验证63-66
• 5.2.1 无重力条件下过热液体中气泡的生长63-64
• 5.2.2 单气泡在水平壁面上的生长和脱离64-66
• 5.3 过热曲面上单气泡运动特性的模拟计算66-72
• 5.3.1 数值计算条件66-67
• 5.3.2 数值计算结果67-72
• 5.4 本章小结72-73
• 6 结论与展望73-75
• 6.1 结论与创新73-74
• 6.2 展望74-75
• 参考文献75-81
• 攻读硕士学位期间发表学术论文情况81-82
• 致谢82-83

【共引文献】

中国期刊全文数据库 前4条

1 赵西增;;A CIP-Based Numerical Simulation of Free Surface Flow Related to Freak Waves[J];China Ocean Engineering;2013年06期

2 田辉;李国君;;改进粒子水平集法及其应用[J];计算物理;2013年06期

3 李维仲;冯玉静;董波;张宁;;直角和贴体坐标系下两种格子Boltzmann方法的比较研究[J];计算力学学报;2014年03期

4 王德;周成;徐国进;;上方侧注式双辊铸轧连铸工艺的数值模拟[J];特种铸造及有色合金;2013年12期

中国重要会议论文全文数据库 前1条

1 冯玉静;董波;李维仲;;使用改进的插值格子Boltzmann方法(GILBM)模拟复杂流动[A];高等学校工程热物理第十九届全国学术会议论文集[C];2013年

中国博士学位论文全文数据库 前4条

1 李婷婷;弹性动力学变形的格子Boltzmann方法研究[D];吉林大学;2013年

2 娄钦;微通道内气液两相界面动力学的格子Boltzmann方法研究[D];华中科技大学;2013年

3 夏盛勇;三氧化二铝液滴碰撞机理及模型研究[D];西北工业大学;2015年

4 蒋成刚;仿生超疏水表面湿润性研究及应用[D];大连理工大学;2014年

中国硕士学位论文全文数据库 前8条

1 辜坚;一种自适应浸入边界法数值模型及其粘性绕流数值模拟研究[D];华南理工大学;2013年

2 韩光华;越浪式波能发电装置的能量转换系统设计研究[D];中国海洋大学;2013年

3 林南盛;海洋结构物运动与载荷的SI-IB-VOF理论研究[D];华南理工大学;2014年

4 卞恩杰;格子Boltzmann方法的改进及应用[D];上海理工大学;2013年

5 金效兴;低微重力场下池内沸腾气泡动力学及传热数值模拟[D];上海大学;2014年

6 常麒;Coflowing液—液两相形成机制以及过相界面传质计算模型的研究[D];天津大学;2013年

7 辛琦;熔融铝液遇水碎化分析及爆炸冲击波能量转化[D];北京理工大学;2015年

8 赖彬彬;IBM-VOF数值模拟方法在船舶水动力分析运用[D];华南理工大学;2015年

，

本文编号：540784