可膨胀过热水系统多气核演化过程与临界过热度分析3.PDF

第 58 卷 第 11 期 2009 年 11 月 物 　理 　学 　报 Vol . 58 ,No . 11 ,November ,2009 ( ) 10003290200958 11 759607 ACTA PHYSICA SINICA 2009 Chin . Phys. Soc . 可膨胀过热水系统多气核演化过程 与临界过热度分析 厉思杰 　白博峰 ( 西安交通大学动力工程多相流国家重点实验室 ,西安 　710049) (2009 年 2 月 7 日收到 ;2009 年 3 月 10 日收到修改稿) 　　使用分子动力学方法 ,研究过热水系统均质沸腾核化过程. 采用 Langevin 动力学方法控制体积可变系统的温 度与压力 ,更好地模拟了沸腾实际物理过程. 得到了液相系统体积连续膨胀 、分子间距逐渐增大 ,最终稳定在汽相 的现象学规律. 当过热温度较高时 ,亚稳态液相系统可能在局部形成不同大小的近球形区域 :气核 ,这些气核是不 稳定的 ,处于不断演化之中. 通过分析分子所受引力与斥力的共同作用 ,得到了气核形成与消亡以及多个气核融合 的机理. 比较了模拟结果与经典沸腾理论的差异 ,提出了气核生长是比气泡生长更为微观过程的认识. 通过研究不 同过热度系统 ,得到了模拟条件下水的均质沸腾临界过热温度为 535 K 左右 ,在已有的实验结果范围内. 关键词 : 过热水系统 , 分子动力学 , 气核 , 临界过热度 PACC : 0570F , 6400 , 4460 　　　 1 引 言 根据经典沸腾理论 ,系统达到饱和状态后可能 并不立即沸腾 ,可以保持在液相且处于亚稳定的过 热状态 ;当系统温度低于临界过热温度时 ,虽然系统 处于热力学亚稳定状态 ,但各个状态量仍然处于不 断的涨落中 ,系统内部的局部瞬时分子数密度及热 力学能可能会偏离系统平均值 ,使局部分子密度下 降而首先形成微小气泡[ 1] ; 当系统温度达到或超过 临界过热温度时 ,液相系统是不稳定的 ,经核化沸腾 图 1 　纯物质 PV 图中过热液体 、过饱和蒸汽区以及范德瓦尔斯 达到稳定的汽相状态. 纯物质沸腾过程各区域如 PV 等温线 图 1 所示. ( - 9 ) 分子动力学验证两个方面研究了快速加热条件下均 在气泡形成初期 ,其直径在 nm 10 m 量级 ,核 化时间在 p s ( 10 - 12 s) 量级 , 因此需要从微观角度对 质核化规律及固体壁面对沸腾过程的影响[ 10 ] . Dirk 沸腾气泡形成进行研究. 分子动力学从分子之间相 Zahn 使用 256 个水分子组成的系统研究沸腾原理 , 互作用的微观角度进行研究 , 已被广泛用于热力学 得出了液体沸腾起始于 ( 10 - 10 m) 量级的大量不稳 的诸多领域[2 —7 ] . Maruyama 等研究了5488 个氩分子 定空核的结论[ 11] . Hisashi Okumura 等采用恒定压力 在固体壁面上的气泡核化过程 ,得到了贴壁气泡的 系统 ,并局部加热系统内的分子 ,观察到了被加热分 分子形貌[8 ] . Nagayamad 等用分子动力学方法得到了 子驱逐周围分子使分子间距扩大形成气核随后又迅 纳米通道内气泡的形成过程 ,并与经典气泡动力学 速被周围分子压缩消失的现象[ 12 ] . 李　等通过模拟 理论进行了比较[9 ] . Carey 等从热力学理论分析与 氩分子系统 ,将核化过程分为四部分 : 团聚段 、成核