气泡临界半径的解释.docxVIP

中国工程热物理学会传热传质学第十届年会编号：013140气泡临界半径的解释陆建峰彭晓峰清华大学热能工程系北京】00084电话：010mail：pxf-dte唾jtsinghua．edu．cn摘要从动力学的角度对胚泡临界半径的研究已经相当的深入，但从多方面加以全面 分析的却是较少。本文主要运用热流法和广义力法得出胚泡的临界半径，实际上也得出了对胚泡演化过程的正确数理描述。关键词杨氏方程，热流，广义力1．概述气泡临界半径是沸腾中的重要参量，对此用经典的杨氏方程和动力学方程已有较 满意的答案，但是动力法得出临界半径不是唯一的方法。柴立和等【1]力图通过演化方程的分析来得出临界半径的值，取得有意义的结果。本文在这些基础上，运用热流法t广义力法和能量法分析得出临界半径的值，经过对比得出同样的结果，可见对于临界半径和演化过程的分析有一定的有效’肚。’2．热力学方程法最初获得沸腾临界半径是通过热力学方程的方法，即由Young—Laplace方程，，r只=R=只+=(1)‘和Clausiuz．Clapeyron方程，—皇fL≥』卫(2)瓦(v”一V’)△zk求解得到r+％2箍(3)厂v艘一一卵这就是气泡能液体中存在的临界半径。 在用非线性分析法来分析气泡演化中，柴等[1】从热力学理论出发初步研究了液体重沸腾核化的非线性特征，导出液相中形成一个半径为r的球形气泡核心所引起的统 自由焓的变化为：△G：一兰里蛀些r3+4聊z(4)3r,其中右边第一项为体积焓增量，第二项为表面张力增量，看成梯度系统后导出非平衡 态的演化方程为生：一堂：塑生竖：一沏．r=ar2-br々一=—‘一a舾r(5)L3J西静I其中。：!堡k些，6：8yro(6)，；式中第一项为非线形项，第二项为线形项。从热力学稳定性导出方程的临界点为，·：o，r．诅：粤，=u，r=c。=——■—?=一(7)(7)P，n詹△Yk这种按梯度系统导出演化方程的不足在于物理意义不明确，虽然对临界半径的结 果没有影响，但直观上看方程左右两侧的量纲是不相同的，而且，柴等[1】没有给出式(5)成立的数学物理证明。实际上，从物理意义上进行分析和改进演化方程是一个重要的出发点。3．热流法热流法主要是通过考虑热流和胚泡成长的能量平衡建立演化方程而取得临界半 径。为简化起见，把胚泡假定为准静态的演化的球形系统．胚泡边界看成壁面，而与 主流区交换热，建立温度场方程。同时，分析中忽略气泡问的非线形作用和重力的影 响。在通常情况下，主流区的温度Z可以认为变化不太，而气泡壁面(相边界)温度 瓦可近似看成是相变温度，因此主流区与气泡壁的总温差可表示为AT=正一正(8)进一步考虑总温差在气泡壁的温度梯度下所需的传热特征长度．由于温度场的对 称性，特征长度应与气泡大小应成正比，设为Z，，如图l所示·图1分析模型 在上述的简化情况下，以热流为考虑对象，则J=2VT=怂r／(∥)=k／r(9)其中五为综合热传导系数(对流影响可以折算)k=怂r／Z(10)由能量守恒定理可得，自由能的增量应该等于热流量．即4nr2J：—aA-G(11)日由于AG主要与，相关联，其它因素与r相关性很小，同时AG不显含自变量r，则～aAG：y堂．纽+一aAG：堂．-加2：；．一aAG(12)at々aqt融氆ar磺静其中g，为包括r在内各种广义自变量。综合式(5)、(9)-(11)，得(一胛2+br)r=4∥2·纠，=4nIc-，=c-，(13)其中c=4nlc，口，b由式(6)给出。 进一步简化其结果则得(kI一屯叶，=1(14)其中k．=b／c．k2=a／c。引入t3=√b，式(13)演化为348r=l／(‘一t：r)=(i／t。)／(1_tr)=(11tI)【l+k3，+霹r2+·-·】(15)=I／七I+．吒／女I·r+酵／七l，+···由上面的展开(应该在无量纲后分解，但以后会看到t，r本身为无量纲数，因此 不影响物理本质)可知，必须‘·，1，否则方程将发散，这也就是热流条件下的半径临界值，一1鱼：!生(16)t3口nh丘△t。也即是胚泡半径的临界值，与熟力学方程法的结果相一致，但这是严格推导而得到的结果。4．广义力法前面用熟流法得§U,-F确的临界胚泡值，而且可以从物理上说明胚泡的演化，至少 在热流为控制因素时是较为成功的方法。在胚泡初始阶段惯性因素居于主导地位，这时热流法可能是欠妥的，下面我们可以换一个角度来考虑这个问题，即从广义力(动力系统)，等效力(等效势)的角度来导出初始态的演化方程。i．广义力法由Hamiton原理可得g．-等，毒一罢‘17)其中H为Hamiton作用量，此处可取为H=AG，利用第二式，并令(18)P=m’·一m。为有效质量。令m‘=p+·，3，其中p’为有效密度，此处已包含相变焙变等众多 物理效应