微米粗糙结构表面上汽泡沸腾成核的热力学分析 董丽宁.doc

中国工程热物理学会 传热传质学 学术会议论文 编号：113428 董丽宁1，全晓军1，郑平1，王蓉2 上海交通大学机械与动力工程学院200240） （2. 东风汽车股份有限公司汽车分公司，湖北441004） （Tel.: 021 E-mail: pingcheng@sjtu.edu.cn） 摘要：基于Gibbs自由能和有效能的变化，对微米粗糙结构的加热面上进行热力学分析， 关键词：非均相成核微米粗糙结构亲疏水性 1. 引言 近年来，随着微加工工艺的蓬勃发展，各种微结构的加热表面广泛应用于反应器和换热设备中。其中，表面粗糙度和表面亲/疏水性对沸腾成核和换热的影响作用受到广泛的关注[1]。粗糙表面上的微结构有助于气体的滞留和增大核化点密度，从而降低初始过热度和核态沸腾过热度，进而提高沸腾换热系数[2,3]。 当液相温度场呈均匀分布时，光滑固-液界面上的非均相汽泡成核可以使用经典的均相成核理论[4,5]。然而，对于具有微结构的表面，由于汽泡形状发生变化，光滑表面的非均相成核理论不再适用。Cole[6]基于热力学理论研究了不同形状参数的凸面和凹面上的非均相成核半径及成核功，但是没有具体分析表面形状及特性对成核的影响。Fletcher[7]指出了凸面上非均相成核过程的尺寸效应，而没有研究凹槽的尺寸效应对成核的影响。Blander和Katz[8]对均相和分均相汽泡成核从理论和实验两方面进行了总结，提出接触角大于68°时，汽泡倾向于非均相成核。 本文从Gibbs自由能和有效能的变化出发，对加热表面上工凹槽或凸起的非均相沸腾成核进行了热力学分析，考虑了表面特性对成核的影响，推导出适用于微米粗糙度表面的非均相成核功公式，并分析了表面粗糙结构及亲/疏水性对成核的影响作用。 2. 非均相成核功 (1) 当成核开始后，系统的有效能由液相Ψl、汽相Ψv以及汽-液界面Ψi三部分组成，表示为 (2) 其中， (3) (4) (5) 成核前后有效能的变化即成核功， (6) 其中，mT为系统的总质量，mv为蒸汽质量；自由能定义为，之差即为自由能的变化。’s方程可知汽-液界面的表面张力有如下关系，。 对于具有微米粗糙度的过热表面，假设其表面形貌由微米尺度的凸起和凹穴构成，其曲率半径R的尺寸在微米范围之内。汽泡胚核在凸面和凹面上的沸腾成核模型分别如图1(a)和(b)所示。凸面和凹面上的沸腾汽泡形状不同，汽-液、汽-固界面表面积Alv、Asv以及汽泡体积V也有不同的表达方式。 (7a,b) (8a,b) (9a,b) 其中，θ为接触角，r为汽泡半径，R为表面粗糙结构的曲率半径，α、β、d分别为汽泡形状的几何参数（见图1）；下标+和-分别代表了凸面和凹面的情况。将上述公式(7)-(9)代入公式(6)，可以得到适用于不同结构表面的成核功公式。由于公式形式较为复杂，对其进行简化，得出通用公式形式如下， (10) 其中，，；和分别为与表面结构有关的无量纲参数，是表面曲率半径与汽泡半径之比x和接触角θ的函数，对于不同类型的表面，其表达方式有所区别： (i) 对均相成核（）， (11a) (ii) 对于平面（）， (11b) (iii) 对凸面和凹面， (11c) 3. 结果分析及讨论 由上述公式的推导过程可以看出，成核功与表面的结构形态有关，不同结构的表面上（平面、凸面或凹面），成核功的公式形式也有所不同。过热表面上的沸腾汽泡成核受到接触角θ的影响；此外，曲率半径R与汽泡半径r的比值x也会影响非均相成核功。标准大气压下（P=101325Pa）水的物性参数如下：饱和温度Tsat=373K，密度ρ=0.598kg/m3，汽化潜热hfg=2257000J/kg，表面张力σlv=0.059N/m，气体常数Rg=461J/kg?K。将上述参数以及不同的曲率半径R和接触角θ值，分别代入公式(10)-(11)，可以得到不同表面条件下的均相及非均