气泡成核过程的格子Boltzmann方法模拟米.PDFVIP

物理 学报 ActaPhys．Sin．Vo1．62，No．17(2013)176401 气泡成核过程的格子Boltzmann方法模拟米 曾建邦 ) 李隆键2)十 蒋方明)丰 1)(中国科学院广州能源研究所先进能源系统实验室，中国科学院可再生能源重点实验室，广州 510640) 2)(重庆大学动力工程学院，低品位能源利用技术及系统教育部重点实验室，重庆 400030) (2013年3月 19曰收到；2013年5月 l4曰收到修改稿) 利用精确差分格子Boltzmann模型探讨水在特定温度下的亚稳态及不稳定平衡态，获得等温相变过程中形成气 泡和液滴的条件，模型预测结果与理论解符合 良好．在该等温模型的基础上耦合能量方程，通过调节流体一壁面相互 作用力获得不同的气泡与固壁间接触角，从而建立了一种新的描述气液相变的格子 Boltzmann理论模型．利用该新 模型模拟不同流体 一壁面相互作用力下凹坑气泡成核过程，再现 了气泡成核过程中的三阶段特性；探讨了接触角、 曲率半径及气泡体积随气泡成核过程的变化关系，获得了与文献结果定性符合的曲率一气泡体积关系曲线． 关键词：格子Boltzmann方法，气泡成核过程，气液相变，接触角 PACS：64．60．qe，63．70．+h，68．35．Rh DoI：10．7498／aps．62．176401 变化)．Bank0 5)首次建立气泡成核理论，并提出 1引 言 沸腾表面凹坑成核条件；Griffith和Wallis0[J测出了 工程表面上水的静态接触角对气泡成核过程 的影 亚稳态两相流动过程研究，对相变成核机理 响，并定性地给出了气泡体积随曲率的变化关系； 和反应堆事故暂态过程研究十分有用，常用于解 Sato和 Matsumura[7】针对 Bankoff5[]成核理论进行 释反应堆失水事故喷放过程 中遇到的一些现象 1J． 改进，提出了新的预测沸腾壁面过热度 的公式；随 而相变成核又是沸腾传热学研究的重要 内容，更 后 Davis和 Anderson8[]将 Sato和 Matsumura[7]的 是沸腾传热学最基础的问题．因此，研究气泡成核 预测公式推广到非圆形气泡成核条件；Lorenz等 9【J 过程具有重要理论及工程意义．稍显遗憾 的是人 在Bankof~5]理论基础上引入新的凹坑半径条件； 们至今对气泡成核过程所基于的微尺度传输机理 Wang等 [1o]通过最小化单个沸腾系统的赫尔姆兹 还未完全弄清 此外，在实际工程 中，沸腾表面 自由能，改进了接触角条件；Mikic和Rohsenow[IlJ， 大多数包含 自然的或机械加工形成的凹坑、刻痕 Judd和 Hwang{12]以及 Dhir[131先后讨论了相邻活 等，呈半圆槽和沟槽形状，大小从 肉眼可见到微观 化核心问存在相互作用；Kenning和 Yan[14】在池沸 尺寸 I3J．Clark等 4【J用显微镜观察到活化凹坑尺 腾实验中发现气泡生长形成的过冷区域大小等于 寸大约 10lxm到 100btm，要描述如此小的范围内 该点上气泡脱离时的投影范围．Zhang和 Shoji[15J 的气泡成核过程，宏观数值求解方法是很难胜任 对人工 凹坑 实验 的活化核心 间相互作用进行 的．因此，学者们另辟蹊径，以求在理论上有所突 分类． 破，但至今无论是从理论还是实验上，只能探寻气 要数值模拟气泡成核过程，首先必须寻求适 泡成核条件，以及活化核心的相互作用，而无法 了 合