鄭福明，王暢

1 海軍裝備部，北京100841

2 中國艦船研究設(shè)計中心，湖北武漢430064

0 引 言

過冷沸騰傳熱與單相傳熱存在差異的主要原因體現(xiàn)在2 個方面：一是汽泡動力學效應(yīng)；二是由汽泡引起的湍流擾動。兩者的聯(lián)合作用導(dǎo)致過冷沸騰傳熱強度遠高于單相傳熱。近年來，研究人員又發(fā)現(xiàn)由窄矩形通道組成的換熱元件具有結(jié)構(gòu)緊湊、換熱性能優(yōu)越的特點，因此，其通過結(jié)合過冷沸騰與窄矩形通道強化換熱的優(yōu)點，設(shè)計出了多種適用于船舶領(lǐng)域的換熱設(shè)備。

已有研究表明，船用換熱設(shè)備與陸用換熱器在運行環(huán)境方面存在的差異主要在于船舶受搖擺運動的影響。因此，通過掌握船舶搖擺運動狀態(tài)對換熱設(shè)備熱工水力特性的影響規(guī)律及影響程度，便可以達到提高設(shè)備運行安全與可靠性的目的?？梢?，搖擺條件下不同類型換熱器內(nèi)的流動傳熱特性理論與試驗研究是傳熱傳質(zhì)領(lǐng)域的研究熱點之一。

秦勝杰和魏敬華等［1-2］分別通過理論分析及數(shù)值計算，研究了搖擺運動對過冷沸騰汽泡脫離點的影響，認為由搖擺運動引起的流量波動是導(dǎo)致?lián)Q熱特性發(fā)生變化的主要原因。幸奠川和鄢炳火等［3-4］分別針對搖擺狀態(tài)下單相強制循環(huán)流動特性進行了試驗及數(shù)值研究，同樣發(fā)現(xiàn)當搖擺運動引起流量發(fā)生波動時，將導(dǎo)致單相流動傳熱與定常流存在差異。王暢等［5-7］針對搖擺狀態(tài)下的流量變化特性進行了理論及試驗研究，發(fā)現(xiàn)對于強迫循環(huán)而言，搖擺運動并不足以導(dǎo)致單相流動發(fā)生波動。

顯然，現(xiàn)有搖擺狀態(tài)下的熱工水力特性研究多集中于單相流動傳熱特性及汽泡脫離行為，尚未發(fā)現(xiàn)有針對搖擺條件下過冷沸騰流動傳熱特性變化規(guī)律的相關(guān)研究報道。

本文將基于兩相分相流模型，通過建立靜止及搖擺條件下窄矩形通道內(nèi)的強迫循環(huán)過冷沸騰數(shù)學模型，探討搖擺運動對過冷沸騰傳熱特性的影響機理，并對理論分析結(jié)論進行試驗驗證。

1 靜止狀態(tài)受力分析

Yeoh 等［8］認為，豎直靜止流道內(nèi)的汽泡在脫離核化點前的受力如圖1 所示，汽泡同時受表面張力Fs、不穩(wěn)定曳力FDu、剪切升力FSL、液體壓力Fh、接觸壓力Fcp、流動方向穩(wěn)定曳力Fqs及浮力Fb的作用。

圖1 汽泡脫離受力分析Fig.1 Force analysis of bubble detaching

因此，汽泡在流道的動量守恒可由如下方程組成。

垂直于壁面方向：

平行于壁面方向：

式中：ρg為汽相密度；Vb為汽泡體積；ub為汽泡運動速度。

由文獻［9］可知，窄矩形通道內(nèi)的過冷沸騰汽泡行為與圓管內(nèi)的差異極大。如圖2 所示，脫離核化點的汽泡將在浮力Fb、曳力FD、由邊界層蒸發(fā)產(chǎn)生的外推力FI及由速度梯度引起的力FR的作用下滑移。圖中，Tsat為飽和溫度。

圖2 汽泡滑移示意圖Fig.2 Schematic of bubble sliding

Lee 和Mudawwar［10］以及Lin［9］等認為，在平行于流動方向，于豎直靜止狀態(tài)下聚合形成的大汽泡在沿加熱壁面滑移的過程中受到的浮力Fb和曳力FD為：

式中：ur為汽泡與其中心處流體的相對速度；ρl為液相密度；Cd為汽泡曳力系數(shù)；Db為汽泡脫離直徑；Lb為滑移汽泡的長度。

Weisman 和Pei［11］認為汽泡在滑移聚合而形成橢球形汽彈的過程中，由于受鄰近的汽彈影響，汽彈在沿加熱壁面滑移時于徑向的增長程度并不明顯，因此可以近似地認為，滑移汽彈的當量直徑等于汽泡的脫離半徑2rd，橢球形汽彈長邊與短邊長度之比約為3∶1。由此，有

在垂直于壁面方向，由汽泡底部的過熱微液層蒸發(fā)而產(chǎn)生的驅(qū)使汽泡遠離壁面的推力FI為

其中，

式中：q 為熱流密度；Vbs為蒸發(fā)速度，ilg為汽化潛熱。

此外，汽泡在近壁面流體的速度梯度的作用下會出現(xiàn)旋轉(zhuǎn)和翻滾，從而產(chǎn)生朝向壁面的力FR而限制汽泡的橫向運動。

式中：ul為液相流速；y 為流體距壁面的距離。

2 搖擺狀態(tài)受力分析

由文獻［6，12］可知，搖擺運動對流動傳熱特性的影響主要體現(xiàn)在2 個方面：一是搖擺運動將導(dǎo)致系統(tǒng)設(shè)備之間的空間相對位置發(fā)生變化，進而引起自然循環(huán)驅(qū)動壓頭發(fā)生改變；二是搖擺運動將引入附加慣性力場。

為便于分析搖擺運動的影響程度，本文采用分相流模型，分別針對搖擺狀態(tài)下過冷流動沸騰的液相及汽相變化規(guī)律進行研究。

2.1 搖擺運動對液相傳熱的影響

由于液相流動與傳熱特性相互耦合，當流動特性不發(fā)生改變時，傳熱特性也不發(fā)生變化。由文獻［6］可知，對于強迫循環(huán)流動而言，由空間位置改變引起的自然循環(huán)驅(qū)動力變化相對于循環(huán)泵的驅(qū)動壓頭而言極其微弱，因此，只需考慮由搖擺運動引起的附加慣性力對液相流動傳熱特性的影響。

根據(jù)力的作用方向，由搖擺運動引起的附加慣性力可以分解為平行于流動方向的附加慣性力和垂直于流動方向的附加慣性力。平行于流動方向的附加慣性力將直接影響流體的速度，而垂直于流動方向的附加慣性力則主要對微觀流場產(chǎn)生影響。

文獻［6］中的研究表明，當強迫循環(huán)驅(qū)動力大于10 倍的平行于流動方向的附加慣性力時，搖擺運動不足以引起流量出現(xiàn)波動，而在工程設(shè)計過程中，為了保證系統(tǒng)可靠性，循環(huán)泵的驅(qū)動壓頭通常都留有足夠的裕量，因此泵的驅(qū)動力遠大于平行于流動方向的附加慣性力。此外，文獻［13］采用CFD 方法針對垂直于流動方向的附加慣性力作用下的微觀流動傳熱特性進行了研究，發(fā)現(xiàn)盡管附加慣性力會誘導(dǎo)流道內(nèi)的二次流隨搖擺運動發(fā)生顯著變化，但相對于主流的湍流強度而言，其引起的流動傳熱特性變化極其微弱，不足以導(dǎo)致宏觀流動傳熱特性發(fā)生顯著變化。

因此，在強迫循環(huán)工況下，搖擺運動對過冷沸騰中液相傳熱特性的影響完全可以忽略。

2.2 搖擺運動對汽相傳熱的影響

當搖擺運動引起的附加慣性力作用于汽泡后，如圖3 所示，窄通道內(nèi)處于滑移狀態(tài)的汽泡除了受浮力Fb、曳力FD、由邊界層蒸發(fā)產(chǎn)生的外推力FI及由速度梯度引起的力FR的影響外，還會受到由搖擺運動引起的離心力Fce及切向力Fta的影響。

式中：ω 為搖擺角速度；β為搖擺角加速度；r 為汽泡距搖擺軸心的距離。

圖3 搖擺狀態(tài)下滑移汽泡受力分析Fig.3 Force analysis of sliding bubble in periodic force field in rolling motion

此外，由于流道的方向隨搖擺運動發(fā)生了變化，因此汽泡受到的浮力作用方向也隨之發(fā)生改變。浮力在流動方向及垂直于壁面方向上的分力分別為

因此，搖擺運動狀態(tài)下，汽泡在沿流動方向及垂直于流動方向上所受的合力可表示為：

通過對汽泡在不同方向上的受力進行分析（圖4，圖中tRo為搖擺周期），發(fā)現(xiàn)在不同搖擺角度和搖擺周期狀態(tài)下，汽泡在垂直于壁面方向上所受的合力Fy始終小于0，這表明在附加慣性力作用下，流道內(nèi)的汽泡仍處于貼附在壁面液膜上的滑移狀態(tài)，因此不會出現(xiàn)由于汽泡脫離壁面、進入過冷主流而引發(fā)的湮滅擾動現(xiàn)象。由文獻［9］可知，汽泡底部液膜的厚度δ 只與汽泡脫離直徑、熱流密度、質(zhì)量流速及系統(tǒng)壓力相關(guān)，而當搖擺運動未引起系統(tǒng)壓力及質(zhì)量流速出現(xiàn)顯著變化時，汽泡底部的液膜厚度δ 也不會發(fā)生明顯改變。因此，在垂直于壁面方向上，由搖擺運動引起的附加慣性力場對汽泡的影響可以忽略。

圖4 滑移汽泡在垂直于壁面方向的受力分析Fig.4 Force analysis of sliding bubble in direction of perpendicular to the wall

搖擺狀態(tài)下汽泡沿流動方向受到的附加慣性力、曳力及浮力如圖5 所示。顯然，即使在當前最劇烈的搖擺工況下，浮力Fb，z及曳力FD隨搖擺運動的變化幅度仍然極其微小，相對而言完全可以忽略。此外，盡管汽泡所受的離心力Fce發(fā)生了明顯的周期性變化，但由于離心力的數(shù)量級僅為浮力及曳力的百萬分之一，因而其影響相對于浮力及曳力而言完全可以忽略。

因此，當搖擺運動沒有引起流量出現(xiàn)周期性的變化時，其引入的附加慣性力場對過冷沸騰流動中汽泡行為的影響也可以完全忽略。

3 搖擺狀態(tài)下過冷沸騰傳熱

理論分析表明，當搖擺狀態(tài)下的質(zhì)量流量沒有出現(xiàn)明顯的波動時，其對過冷沸騰液相及汽相流動傳熱特性的影響均可以忽略。為了驗證理論分析的準確性，本文針對以上分析結(jié)論進行了試驗驗證，相關(guān)試驗系統(tǒng)的資料詳見文獻［7，14］。

如圖6 及圖7 所示，與理論分析的結(jié)論一致，當搖擺運動未引起窄通道內(nèi)的質(zhì)量流量G 出現(xiàn)明顯變化時，搖擺狀態(tài)下窄通道內(nèi)的過冷沸騰傳熱系數(shù)h 與豎直靜止狀態(tài)基本一致。

圖6 搖擺狀態(tài)下過冷沸騰質(zhì)量流量變化Fig.6 Variation of mass flux of sub-cooled boiling under rolling motion

圖7 搖擺狀態(tài)下過冷沸騰傳熱特性變化Fig.7 Variation of heat transfer characteristic of sub-cooled boiling under rolling motion

4 結(jié) 論

本文針對搖擺狀態(tài)下窄通道內(nèi)過冷沸騰傳熱特性進行了理論分析及試驗驗證，得到以下結(jié)論：

1）當搖擺運動未引起過冷沸騰質(zhì)量流量發(fā)生明顯波動時，只需考慮搖擺狀態(tài)下附加慣性力及浮力變化對汽泡行為的影響。

2）搖擺狀態(tài)下垂直于壁面方向的附加慣性力無法驅(qū)動汽泡脫離壁面進入過冷主流，汽泡在搖擺運動過程中始終處于貼附于加熱壁面滑移的狀態(tài)，因此不會出現(xiàn)由于汽泡進入過冷主流“冷凝—湮滅”引起的擾動現(xiàn)象。

3）搖擺狀態(tài)下平行于流動方向的浮力及曳力隨搖擺運動的變化幅度極其微弱，不足以引起系統(tǒng)內(nèi)質(zhì)量流量發(fā)生改變。盡管離心力隨搖擺運動發(fā)生了顯著的波動，但其量級相對于浮力及曳力完全可以忽略。

4）對于窄通道內(nèi)的強迫循環(huán)流動而言，船舶搖擺運動不足以導(dǎo)致過冷沸騰傳熱特性發(fā)生顯著變化。

［1］秦勝杰，高璞珍.搖擺運動對過冷沸騰流體中汽泡受力的影響［J］.核動力工程，2008，29（2）：20-23.QIN Shengjie，GAO Puzhen. Effect of rolling motion on forces acting on bubbles in sub-cooled boiling flow［J］.Nuclear Power Engineering，2008，29（2）：20-23.

［2］魏敬華，潘良明，徐建軍，等.附加慣性力作用下豎直矩形流道內(nèi)過冷流動沸騰的數(shù)值模擬［J］. 化工學報，2011，62（5）：1239-1245.WEI Jinghua，PAN Liangming，XU Jianjun，et al. Nu?merical simulation of subcooled flow boiling in vertical rectangular channel under additional inertial forces［J］.CIESC Journal，2011，62（5）：1239-1245.

［3］幸奠川，閻昌琪，劉洋，等.搖擺條件下矩形通道內(nèi)單相水強制循環(huán)流動特性研究［J］.核動力工程，2011，32（1）：112-116.XING Dianchuan，YAN Changqi，LIU Yang，et al. Re?search on single-phase forced circulation in rectangular channel flow characteristics under rolling conditions［J］.Nuclear Power Engineering，2011，32（1）：112-116.

［4］鄢炳火，顧漢洋，楊燕華，等.搖擺條件下矩形管內(nèi)湍流流動特性的LES 研究［J］. 核動力工程，2011，32（1）：108-111，116.YAN Binghuo，GU Hanyang，YANG Yanhua，et al.LES research on flowing characteristics of turbulent flow in rectangular tubes in rolling motion［J］. Nuclear Power Engineering，2011，32（1）：108-111，116.

［5］王暢，高璞珍，譚思超，等.搖擺條件下強迫循環(huán)流量脈動特性分析［J］.核動力工程，2012，33（3）：56-60.WANG Chang，GAO Puzhen，TAN Sichao，et al.Anal?ysis on characteristic of forced circulation flow rate pul?sation under rolling condition［J］. Nuclear Power Engi?neering，2012，33（3）：56-60.

［6］WANG C，WANG S W，WANG H，et al. Investiga?tion of flow pulsation characteristic in single-phase forced circulation under rolling motion［J］. Annals of Nuclear Energy，2014，64：50-56.

［7］WANG C，LI X H，WANG H，et al. Experimental study of friction and heat transfer characteristics of pulsating flow in rectangular under rolling motion［J］. Progress in Nuclear Energy，2014，71：73-81.

［8］YEOH G H，CHEUNG S C P，TU J Y，et al.Fundamen?tal consideration of wall heat partition of vertical sub?cooled boiling flows［J］. International Journal of Heat and Mass Transfer，2008，51（15/16）：3840-3853.

［9］LIN K W，LEE C H，HOURNG L W，et al. A theoreti?cal and experimental study on subcooled flow boiling at high heat flux［J］. W?rme-Und Stoffübertragung，1993，29（5）：319-327.

［10］LEE C H，MUDAWWAR I. A mechanistic critical heat flux model for subcooled flow boiling based on lo?cal bulk flow conditions［J］. International Journal of Multiphase Flow，1988，14（6）：711-728.

［11］WEISMAN J，PEI B S. Prediction of critical heat flux in flow boiling at low qualities［J］. International Jour?nal of Heat and Mass Transfer，1983，26（10）：1463-1477.

［12］WANG C，LI X H，WANG H，et al. Experimental study on single-phase heat transfer of natural circula?tion in circular pipe under rolling motion condition［J］. Nuclear Engineering and Design，2014，273：497-504.

［13］王暢，高璞珍，譚思超.搖擺運動對單相流動傳熱特性的影響［C］//中國工程熱物理學會傳熱傳質(zhì)學學學術(shù)會議.東莞，2012.WANG Chang，GAO Puzhen，TAN Sichao. Effect of rolling motion on single-phase flow and heat transfer characteristics under rolling motion［C］//Heat and Mass Transfer Conference. Dongguan，China,2012.

［14］WANG C，WANG H，WANG S W，et al. Experimen?tal study of boiling incipience in vertical narrow rect?angular channel［J］. Annals of Nuclear Energy，2014，66：152-160.