气泡成核过程的格子boltzmann方法模拟-物理学报.pdfVIP

物理学报 Acta Phys. Sin. Vol. 62, No. 17 (2013) 176401 气泡成核过程的格子Boltzmann 方法模拟* 1 2† 1‡ 曾建邦 李隆键 蒋方明 1) ( 中国科学院广州能源研究所先进能源系统实验室, 中国科学院可再生能源重点实验室, 广州 510640 ) 2) ( 重庆大学动力工程学院, 低品位能源利用技术及系统教育部重点实验室, 重庆 400030 ) ( 2013 年3 月19 日收到; 2013 年5 月14 日收到修改稿) 利用精确差分格子Boltzmann 模型探讨水在特定温度下的亚稳态及不稳定平衡态, 获得等温相变过程中形成气 泡和液滴的条件, 模型预测结果与理论解符合良好. 在该等温模型的基础上耦合能量方程, 通过调节流体-壁面相互 作用力获得不同的气泡与固壁间接触角, 从而建立了一种新的描述气液相变的格子Boltzmann 理论模型. 利用该新 模型模拟不同流体- 壁面相互作用力下凹坑气泡成核过程, 再现了气泡成核过程中的三阶段特性; 探讨了接触角、 曲率半径及气泡体积随气泡成核过程的变化关系, 获得了与文献结果定性符合的曲率-气泡体积关系曲线. 关键词: 格子Boltzmann 方法, 气泡成核过程, 气液相变, 接触角 PACS: 64.60.qe, 63.70.+h, 68.35.Rh DOI: 10.7498/aps.62.176401 变化). Bankoff5 首次建立气泡成核理论, 并提出 1 引言 沸腾表面凹坑成核条件; Griffith 和Wallis6 测出了 工程表面上水的静态接触角对气泡成核过程的影 亚稳态两相流动过程研究, 对相变成核机理 响, 并定性地给出了气泡体积随曲率的变化关系; 和反应堆事故暂态过程研究十分有用, 常用于解 Sato 和Matsumura7 针对Bankoff5 成核理论进行 1 释反应堆失水事故喷放过程中遇到的一些现象 . 改进, 提出了新的预测沸腾壁面过热度的公式; 随 而相变成核又是沸腾传热学研究的重要内容, 更 后Davis 和Anderson8 将Sato 和Matsumura7 的 是沸腾传热学最基础的问题. 因此, 研究气泡成核 预测公式推广到非圆形气泡成核条件; Lorenz 等9 过程具有重要理论及工程意义. 稍显遗憾的是人 在Bankoff5 理论基础上引入新的凹坑半径条件; 们至今对气泡成核过程所基于的微尺度传输机理 10 Wang 等 通过最小化单个沸腾系统的赫尔姆兹 2 11 还未完全弄清 . 此外, 在实际工程中, 沸腾表面 自由能, 改进了接触角条件; Mikic 和Rohsenow , 大多数包含自然的或机械加工形成的凹坑、刻痕 Judd 和Hwang12 以及Dhir13 先后讨论了相邻活 等, 呈半圆槽和沟槽形状, 大小从肉眼可见到微观 化核心间存在相互作用; Kenning 和Yan14 在池沸 3 4 尺寸 . Clark 等