邱國棟，蔡偉華，吳志勇，姜益強，姚 楊（.哈爾濱工業(yè)大學(xué)市政環(huán)境工程學(xué)院，50090哈爾濱；.哈爾濱工業(yè)大學(xué)能源科學(xué)與工程學(xué)院，5000哈爾濱）

Lee相變傳質(zhì)方程中傳質(zhì)系數(shù)取值的分析

邱國棟1，蔡偉華2，吳志勇1，姜益強1，姚 楊1
（1.哈爾濱工業(yè)大學(xué)市政環(huán)境工程學(xué)院，150090哈爾濱；2.哈爾濱工業(yè)大學(xué)能源科學(xué)與工程學(xué)院，150001哈爾濱）

針對Lee相變傳質(zhì)方程中傳質(zhì)系數(shù)主要依靠經(jīng)驗取值的問題，本文嘗試提出一種確定傳質(zhì)系數(shù)的方法，給出了評價傳質(zhì)系數(shù)合理性的2個指標——潛熱份額（潛熱換熱量占總換熱量的份額）和飽和溫差（流體溫度與飽和溫度的差值），建立了穩(wěn)態(tài)分析模型，通過理論推導(dǎo)，得出了潛熱份額和飽和溫差的理論表達式，據(jù)此分析了傳質(zhì)系數(shù)對模擬結(jié)果的影響，已有特定條件下的理論解驗證了本文分析結(jié)果的正確性.分析結(jié)果表明傳質(zhì)系數(shù)越大，計算精度越高.基于分析結(jié)果解釋了不同文獻對該系數(shù)取值差異很大的原因，并給出了傳質(zhì)系數(shù)的通用取值方法，最后推薦了常見工況模擬中傳質(zhì)系數(shù)的合理取值范圍.

相變；傳質(zhì)；冷凝；沸騰；數(shù)值模擬

在采用計算流體動力學(xué)（CFD）方法模擬冷凝或沸騰過程時，相變模型的準確性對整個冷凝或沸騰的模擬結(jié)果至關(guān)重要，這關(guān)系到能否準確地模擬出溫度場、干度變化、氣液相分布及整個換熱過程的合理性.在現(xiàn)有相變模型中，Lee［1］所提出的相變傳質(zhì)方程是在Hertz Knudsen方程［2-3］的基礎(chǔ)上推導(dǎo)而來，具有形式簡單、易于計算、可靠性高等優(yōu)點，因此得到了廣泛應(yīng)用［4-10］.但是Lee相變傳質(zhì)方程中存在一個難以確定的系數(shù)，該系數(shù)通常是根據(jù)研究者的經(jīng)驗來取值，缺乏理論依據(jù)，最終導(dǎo)致該系數(shù)在不同文獻中的取值差異非常大［5，6，9-12］.由于傳質(zhì)系數(shù)對模擬結(jié)果的準確性具有重要影響，因此本文針對傳質(zhì)系數(shù)的取值方法進行分析.

1 傳質(zhì)系數(shù)問題

關(guān)于單位體積冷凝或沸騰傳質(zhì)速率的計算，Lee［1］給出了結(jié)構(gòu)簡單、應(yīng)用廣泛的相變傳質(zhì)方程，即Lee方程，其形式為

式中：S為液相連續(xù)性方程中的質(zhì)量源項，kg／（m3·s）；r為相變傳質(zhì)系數(shù)，1／s；αl、αg分別是液相和氣相的體積分數(shù)；ρl、ρg分別是液相和氣相的密度，kg／m3；T、Ts分別是流體溫度和飽和溫度，K.

Lee方程涉及界面濃度（即單位體積內(nèi)的相界面面積）的計算，由于氣液兩相流存在復(fù)雜的流型，要計算出所有流型下的界面濃度十分困難，目前只能計算出某些特定流型下的界面濃度. ANSYS14.0 FLUENT理論手冊［13］給出了細泡狀流下的界面濃度計算公式其過程如下：

式中，Ai為界面面積，m2；Vcell為單元格體積，m3；d為氣泡（或液滴）直徑，m.

由此得到泡狀流沸騰過程r的表達式為

式中：β為調(diào)節(jié)系數(shù)；M為摩爾質(zhì)量，kg／mol；R為通用氣體常數(shù)，8 314 J／（mol·K）；hfg為氣化潛熱，J／kg；Ts為流體飽和溫度，K.

盡管式（3）、（4）是在泡狀流沸騰和霧狀流冷凝的情況下推導(dǎo)出的相變傳質(zhì)系數(shù)表達式，但表達式中的粒徑d和調(diào)節(jié)系數(shù)β仍難以確定．而在其他流型下因界面濃度難以計算，使得其傳質(zhì)系數(shù)r更加難以用理論分析來確定．因此Lee相變傳質(zhì)方程在使用時，相變傳質(zhì)系數(shù)一般取為經(jīng)驗常數(shù).目前不同文獻對該相變傳質(zhì)系數(shù)的取值差異非常大［5，6，9-12］，為0.1～107s-1.且這些文獻只給出了經(jīng)驗值，并未對其合理性進行深入分析.

2 傳質(zhì)系數(shù)對結(jié)果的影響

從理論上導(dǎo)出傳質(zhì)系數(shù)較為困難，本文不考慮哪些因素會影響傳質(zhì)系數(shù)，而是分析傳質(zhì)系數(shù)的不同取值會對相變換熱的模擬結(jié)果產(chǎn)生什么影響，最終通過分析其影響規(guī)律來反推傳質(zhì)系數(shù)的合理取值.

該方法需要一個能夠評價傳質(zhì)系數(shù)合理性的指標，從能量方程（5）出發(fā)來分析.

式中：ρ為流體的密度，kg／m3；h為流體比焓，J／kg；u為流體速度，m／s；λeff為考慮湍流的流體有效導(dǎo)熱系數(shù)，W／（m2·K）.等號右邊第一項為顯熱換熱量，第二相為潛熱換熱量，由于S與r成正比，所以r的大小決定了潛熱換熱量占總換熱量的比重.定義潛熱換熱量占總換熱量的比重為潛熱份額RL．常見的單質(zhì)冷凝（或沸騰）換熱一般認為是定壓相變換熱，在兩相區(qū)流體的溫度變化很小，換熱以相變換熱為主，潛熱份額RL應(yīng)接近1．故潛熱份額RL可作為評價r取值是否合理的指標之一.

很顯然潛熱份額RL隨r的增加而增加，r越大，RL越接近1，換熱過程越接近純相變換熱．理論上當r→+∞時RL→1，即當r→+∞時換熱過程才能保證是純相變換熱．這表明實際冷凝（或沸騰）模擬中，r盡可能地取較大的值，從而可保證模擬得到的溫度場和干度變化接近理論值，此時模擬得到的相界面溫度等于飽和溫度（即界面邊界條件［4，6］）．如果r取值過小，RL接近0，將會導(dǎo)致?lián)Q熱過程接近單相換熱，對于換熱量較大的模擬工況，它使得冷凝過程的氣相溫度遠低于沸點溫度，或沸騰過程的液相溫度遠高于沸點溫度，這樣的過程不僅不合理，而且因相變量太小導(dǎo)致流體干度與實際差異較大，模擬得到的流型和換熱特性也與實際差異較大．這表明兩相態(tài)流體溫度與飽和溫度的差值ΔT（飽和溫差）也可作為評價r取值是否合理的指標之一．

上述分析表明，r取較大的值可以保證任何情況下的模擬結(jié)果與實際接近，具有通用性．但是，r太大會帶來能量方程計算易發(fā)散的問題，r取多大會帶來能量方程易發(fā)散的問題與模擬問題的復(fù)雜性、網(wǎng)格質(zhì)量和算法中松弛因子取值有關(guān)，例如文獻［6］模擬一根水平直管內(nèi)冷凝過程，結(jié)構(gòu)簡單，網(wǎng)格質(zhì)量好，不發(fā)散的r最大可取到750 000；文獻［12］模擬一根蛇形管內(nèi)沸騰過程，結(jié)構(gòu)復(fù)雜，網(wǎng)格質(zhì)量稍低，不發(fā)散的r僅為100，因此r并非越大越好．由于精度和收斂性的矛盾，要求r應(yīng)該在保證精度的前提下盡量取較小的值．為此需要定量分析潛熱份額RL與傳質(zhì)系數(shù)r的數(shù)學(xué)關(guān)系．為了便于理論分析，建立穩(wěn)態(tài)分析模型，取一個圓柱形微元體，如圖1所示，直徑為D，長度為L，進口溫度為Tin=Ts，干度為xin，出口溫度為Tout，干度為xout，質(zhì)量流率為G．

圖1 微元段幾何模型

由于本模型不關(guān)心每個面的換熱情況，只關(guān)心換熱總量，故不妨假設(shè)換熱只在圓周面上進行，面熱流密度為q，對于冷凝過程，顯熱換熱量為

式中cp，m為氣液平均比熱；L很小，微元體內(nèi)平均干度近似等于入口干度，因此氣、液相平均比熱可按下式計算：

潛熱換熱量為

總換熱量為

聯(lián)立式（9）和（10）可得潛熱份額RL為

聯(lián)立式（8）、（9）和（10）可得飽和溫差ΔT為

將式（12）代入式（11）可得

由式（13）可看出，RL隨r的增加而增加，當r大到一定值時，再增加r對RL的影響很?。墒剑?2）可看出，ΔT隨r的增加而減小，當r大到一定值時，再增加r對ΔT的影響很小．這說明當r滿足精度要求時，再增加r的意義不大．

3 傳質(zhì)系數(shù)取值方法

依據(jù)r取值合理性的兩個評價指標，可得合理的r應(yīng)同時滿足如下要求：

1）模擬的換熱過程以相變換熱為主，即RL值接近1；

2）兩相態(tài)的流體溫度與飽和溫度接近，即飽和溫差ΔT接近0．

這兩個指標具有相關(guān)性，當r足夠大時能夠同時滿足，但是在某些情況下并不會同時滿足，例如當長度L或熱流密度q較小時，由式（12）可以看出，即使r較小也可能使ΔT接近0．不同的研究任務(wù)對上述兩個指標的要求不同，應(yīng)根據(jù)研究目的需要，對這兩個指標設(shè)定合理的精度．從式（12）～（15）可以看出，影響RL和ΔT的參數(shù)有r、ρg、ρl、x、hfg、L、cp，m、G、q、D和Ts；對于某一個特定的工況，除了r外其他參數(shù)均為已知．按上述合理的r的取值要求，既保證RL值接近1（如0.95＜RL＜1），同時保證ΔT接近0（如ΔT＜0.1 K），依據(jù)式（12）～（15）可以得到一個合理的r值．

由上述分析可以看出，合理的r并非一個定值，而是一個范圍，這可能是不同文獻對r取值不同的原因．本文選取常見介質(zhì)：水，繪制出不同干度x，不同長度L，不同直徑D，不同熱流密度q，不同對比壓力Pr，不同流量G時潛熱份額RL和飽和溫差ΔT隨r變化的曲線，如圖2～圖7所示．圖中當某些參數(shù)變化時，其他未標出的參數(shù)均為默認值，默認值如下：干度x=0.5，長度L=1 m，直徑D=0.014m，熱流密度q=10 kW·m-2，對比壓力Pr=0.2，流量G=600 kg／（m2·s）.

圖2 不同流量G下RL和ΔT隨r的變化曲線

圖3 不同對比壓力Pr下RL和ΔT隨r的變化曲線

圖4 不同干度x下RL和ΔT隨r的變化曲線

圖5 不同長度L下RL和ΔT隨r的變化曲線

圖6 不同管徑D下RL和ΔT隨r的變化曲線

圖7 不同熱流密度q下RL和ΔT隨r的變化曲線

由圖2～7可看出，不同工況下水的潛熱份額RL和飽和溫差ΔT隨r的變化趨勢基本一致，均表現(xiàn)為RL隨r的增加而增加，ΔT隨r的增加而減?。啥葂，長度L，對比壓力Pr，流量G，直徑D，熱流密度q對RL影響較小或無影響．當r取值較小時，這些參數(shù)對ΔT的影響較大，在常見的取值范圍中，相對來說長度L對ΔT的影響更大一些；但是當r較大時，這些因素對ΔT的影響均變的較小．各因素敏感性分析表明當r大到一定值以后，潛熱份額RL和飽和溫差ΔT幾乎不受其他參數(shù)影響，基本為一個常數(shù)．這意味著對于常見工況，存在一個較通用的r取值范圍．由圖2～7可看出，從精度角度看，常見工況的r推薦值在104以上較為合適，該推薦的范圍之所以能通用主要是因為在該范圍內(nèi)，其他參數(shù)的影響均可忽略．由于計算收斂性受多種因素影響（如模型復(fù)雜性，網(wǎng)格質(zhì)量，計算松弛因子取值等），且部分因素與研究者主觀經(jīng)驗有關(guān)，故從精度角度得到的r取值其收斂性難以預(yù)判．實際應(yīng)用時，可以先從精度分析推薦的范圍中取一個r值試算，如果能量方程不收斂，則降低r取值，如果r取值明顯小于104，則需用r兩個評價指標來評價其合理性．

從泡狀流沸騰的傳質(zhì)系數(shù)理論式（3）可以估算出常見工況r的量級大于106（估算數(shù)據(jù)如下：氣泡（或液滴）直徑d取10-4m，β取1，摩爾質(zhì)量M取0.018 kg／mol，通用氣體常數(shù)R取8 314 J／（mol·K），潛熱hfg取5×105J／kg，Ts取300 K，得r=1.02×106），從本文的分析可以看出該量級使得潛熱份額RL十分接近1，且飽和溫差ΔT十分接近0，可認為是理想的等溫相變過程，與理論公式（3）的前提假設(shè)是一致的．這也驗證了本文分析的合理性.由此還可以看出，雖然理論分析得到的r使得相變過程與理想情況一致，但是由此得到的r并不一定便于在CFD中使用.這是因為由此得到的r太大，使得能量方程容易發(fā)散.這表明花太多精力去研究傳質(zhì)系數(shù)理論解意義不大.

4 結(jié) 論

針對在冷凝或沸騰過程模擬研究中廣泛使用的相變傳質(zhì)Lee方程，建立了穩(wěn)態(tài)分析模型，分析了該方程中的傳質(zhì)系數(shù)對模擬結(jié)果的影響，解決了現(xiàn)有文獻針對該傳質(zhì)系數(shù)的取值多依靠經(jīng)驗而缺乏理論依據(jù)的問題，獲得以下結(jié)論：

1）傳質(zhì)系數(shù)r會影響潛熱份額RL以及飽和溫差ΔT，進而影響溫度場、干度變化和氣液相分布等模擬結(jié)果．

2）傳質(zhì)系數(shù)r增加，潛熱份額RL隨之增加，飽和溫差ΔT隨之減小，當r大到一定值以后，再增加r對RL和ΔT的影響很?。?/p>

3）傳質(zhì)系數(shù)r取較大的值可以使相變換熱接近理論值，從精度角度看，常見工況的r的合理取值范圍推薦為104以上．

［1］LEE W H.A pressure iteration scheme for two-phase flow modeling［J］.Multiphase Transport Fundamentals. Reactor Safety，Applications，1980，1：407-431.

［2］DOBSON M K，CHATO JC.Condensation in smooth horizontal tubes［J］.Journal of Heat Transfer-Transactions of ASME，1998，120：193-213.

［3］KNUDSEN M.The kinetic theory of gases.Somemodern aspects［M］.London：Methuen and Company，1934.

［4］DA RIVA E，DEL COL D.Numerical simulation of condensation in aminichannel［C］／／Proceedings of the ASME Micro／Nanoscale Heat and Mass Transfer International Conference.Shanghai：American Society of Mechanical Engineers，2010：139-145.

［5］BORTOLIN S，DA RIVA E，DELCOLD.Condensation in a square minichannel：application of the vofmethod［J］.Heat Transfer Engineering，2014，35（2）：193-203.

［6］DA RIVA E，DEL COL D.Numerical simulation of laminar liquid film condensation in a horizontal circular minichannel［J］.Journal of Heat Transfer-Transactions of the ASME，2012，134（5）：1-8.

［7］DA RIVA E，DEL COL D.Effect of gravity during condensation of r134a in a circular minichannel［J］. Microgravity Science and Technology，2011，23（1）：87-97.

［8］DA RIVA E，DELCOLD，CAVALLINIA.Modelling of condensation in a circular minichannel by means of the vofmethod［C］／／Proceedings of the Asme International Heat Transfer Conference-2010.Washington，DC：American Society ofMechanical Engineers，2010：205-213.

［9］DA RIVA E，DEL COL D，GARIMELLA S V，et al. The importance of turbulence during condensation in a horizontal circular minichannel［J］.International Journal of Heat and Mass Transfer，2012，55（13／14）：3470-3481.

［10］DE SCHEPPER SC K，HEYNDERICKX G J，MARIN G B.Modeling the evaporation of a hydrocarbon feedstock in the convection section of a steam cracker［J］.Computers and Chemical Engineering，2009，33：122-132.

［11］WU H L，PENG X F，YE P，et al.Simulation of refrigerant flow boiling in serpentine tubes［J］. International Journal of Heat and Mass Transfer，2007，50：1186-1195.

［12］YANG Z，PENG X F，YE P.Numerical and experimental investigation of two phase flow during boiling in a coiled tube［J］.International Journal of Heat and Mass Transfer，2008，51：1003-1016.

［13］ANSYS Inc.ANSYS14.0 FLUENT theory guid［Z］. Canonsburg：ANSYS Inc.2011.

［14］CHISHOLM D.Research note：void fraction during twophase flow［J］.Journal Mechanical Engineering Science，1973，15（3）：1959-1982.

（編輯 趙麗瑩）

Analysis on the value of coefficient ofmass transfer w ith phase change in Lee′s equation

QIU Guodong1，CAIWeihua2，WU Zhiyong1，JIANG Yiqiang1，YAO Yang1
（1.School of Municipal＆Environmental Engineering，Harbin Institute of Technology，150090 Harbin，China；2.School of Energy Science and Engineering，Harbin Institute of Technology，150001 Harbin，China）

To solve the problem that the mass transfer coefficient in Lee＇s equation need to be obtained depending on the experience，amethod was proposed，and two indexes evaluating the rationality of the mass transfer coefficient-latentheat share（the share of latent heat flux in total heat flux）and saturation temperature difference（difference between fluid temperature and saturation temperature）were presented.A steady state model was established to analyze this problem and the expression of latent heat share and saturation temperature difference were obtained by the theoretical guidance.By the expression，the effect ofmass transfer coefficient on the simulation was analyzed.The analysis resultwas validated by the existing theoretical result in particular cases.The results show that the calculation accuracy increases with the increase of the coefficient. Based on these results，the reason for large difference of the coefficient value in different literatures was explained，themethod of how to get the coefficient value was presented，and the value range of easy-to-use was suggested.

phase change；mass transfer；condensation；boiling；numerical simulation

TB61+1

A

0367-6234（2014）12-0015-05

2014-04-08.

教育部博士點基金項目（博導(dǎo)類）（20102302110045）.

邱國棟（1985—），男，博士研究生；姜益強（1973—），男，教授，博士生導(dǎo)師；姚 楊（1963—），女，教授，博士生導(dǎo)師.

姜益強，jyq7245＠sina.com.