不同介电常数工质的气泡受力分析.docxVIP

不同介电常数工质的气泡受力分析 1 强化传热机理 由于能有效地加强核态沸腾交换，党中央的作用引起了人们的注意。关于电场强化换热,Allen等综述了现有电场强化单相以及相变换热的大量的实验结果,并且归纳总结了雷诺数,努塞尔数和沸腾换热中的临界热流密度值等准则数的拟合计算式。Yagoobi等具体的分析了电场强化传热的机理,认为强化的主要原因是受电场力的影响。即对于单相传热,电场库仑力使流体粒子运动加剧,有利于强化传热;而对于相变换热(特别是核态沸腾换热)电场的介电电泳力引起核化气泡的变形和移动,从而搅动了热边界层,换热加强。Zaghdoudi等将气泡假定为进行蠕动流的小球,针对n-戊烷工质计算了气泡的受力关系。本文利用Zaghdoudi提出的理论模型对n-戊烷以及液体和气体相对介电常数较大的制冷工质R123和R134a进行了计算,定量分析了电场的法向力、切向力和合力对沸腾气泡变形的影响。 2 本构模型的静电学特性,主要描述利用guass 本文考虑的是在各向同性电介质中,加热面接地,电极板接正电压这样一个系统,如图1所示。并且假设沸腾气泡为球形,气泡表面在液体中受到电场力、表面张力、粘性力作用,而气泡整体受浮力、重力、介电电泳力作用。考虑到气泡内外工质的介电常数不同,从而在外电场作用下,气泡表面有极化电荷产生,其中由于气泡界面内外工质介电常数的变化趋势对于气泡的上下部正好相反,因此气泡的上部感应出正电荷,下部感应出负电荷,但带电总量为零。 电场作用下流体粒子所受的电场力可写为 其中,为所受的电场力,为电场强度,q为自由电荷密度,ε为工质的介电常数,ρ为极化电荷密度。 利用Guass定理将式(1)表达成Maxwell应力张量形式为: 其中下标j取l时表示为液相,取v时表示为气相。s代表气液界面,为各相的单位法向分量。从而作用于气液界面的合力可写为: 将气泡界面上的条件代入式(3)整理得:(其中分别为液、气两相的切向分量) 由静电学理论,两种不同介质界面的静电场边界条件为: 其中qes为界面的极化电荷量。 代入式(5)中的边界条件以及介质中的场强分布可得气液界面上法向电应力fsn和切向电应力fst分别为: 式(6)和式(7)中E0为电场强度,θ指气泡界面上的点所在半径同z轴的夹角,ε0,εrl,εrv分别为真空介电常数、液体相对介电常数和气体相对介电常数,X为汽液电导率之比;Y是液汽介电常数比值。即有:X=σev/σel,Y=εl/εv=εrl/εrv。将法向力和切向力分别沿水平方向和垂直方向分解为fsnx,fsnz,fstx,fstz,如图1所示。 计算n-戊烷气泡界面所受的电场合力如图2所示。计算中所取的各参数值分别为:E0=20kV/cm,ε0=8.85×10-12F/m,εrl=1.8(n-戊烷),X=0.25,Y=1.78。 3 麻黄内开发内资源量的影响 如图2(a)所示本文计算了从-180°～180°范围内气泡界面所受到的电应力作用,其中在-90°～90°范围内本文的计算结果同Zaghdoudi的相同,如图2(b)所示。另外从图中显示:在气泡顶部(即θ=0°)和底部(即θ=180°)时,气泡分别受到向上和向下的最大垂直力作用,从而使气泡产生纵向拉伸变形,同Herman的实验结果一致,即电场作用下沸腾气泡沿场强方向变长。另外沸腾气泡在-35°和-145°时受到向左的最大水平力作用,在35°和145°时受到向右的最大水平力作用。所以当气泡对称时,此两水平力平衡,不会对气泡产生整体运动。但是在此期间瞬时的平衡极易破坏,一旦破坏,则不平衡的水平力使气泡产生沿接触表面的横向移动。 4 最优条件下充放电场常数之比 另外本文还选取了液体和气体相对介电常数比值较大的R123和R134a工质进行计算。图3是工质R123的计算结果,所取参数为E0=20 kV/cm,ε0=8.85×10-12F/m,εrl=3.46(R123),X=0.25,Y=3.46。图4是制冷工质R134a的计算结果,所取参数为E0=20 kV/cm,ε0=8.85×10-12F/m,εrl=9.5(R134a),x=0.25,Y=9.5。 从图3来看,在-90°处R123气泡受到向右的最大水平力作用;在90°处R123气泡受到向左的最大水平力作用,并且最大水平分力的绝对值比戊烷的对应值大。而R123气泡受到的垂直力沿气泡表面呈周期性的变化。图4显示R134a仍在±90°处受到最大水平力作用,其方向与R123相同,但是最大水平分力的绝对值比R123更大。而在±60°处R134a受到向下的最大垂直分力作用,在±135°处R134a受到向上的最大垂直分力作用。由此可见对于不同工质,液体同气体的相对介电常数之比(Y=εrl/εrv)越大,则沸腾中形成的气泡受到电场力的水平分力