乙醇溶液液滴降壓蒸發過程傳熱傳質特性

王茉，劉璐，王鵬程，米夢龍，劉彥豐

（華北電力大學能源動力與機械工程學院，河北 保定 071003）

?

乙醇溶液液滴降壓蒸發過程傳熱傳質特性

王茉，劉璐，王鵬程，米夢龍，劉彥豐

（華北電力大學能源動力與機械工程學院，河北 保定 071003）

摘要：針對單個乙醇溶液液滴在降壓環境下蒸發的傳熱傳質過程建立了數學模型。模型基于液相的能量守恒和傳質擴散理論，利用經典拓展模型計算液滴的質量蒸發率，并引入活度系數考慮液滴表面的蒸氣分壓。采用液滴懸掛法進行實驗，分別記錄了乙醇溶液液滴和乙酸溶液液滴在降壓蒸發過程中的液滴內溫度變化。將實驗數據與計算結果對比，驗證了模型的有效性。通過模型計算獲得了液滴內部溫度分布以及濃度分布隨時間的變化。結果表明：快速降壓階段空氣流動較快，加之乙醇工質易揮發，液滴表面溫度下降迅速，液滴內部溫差和乙醇濃度梯度較大；壓力穩定后，空氣流速為零，液滴內部溫差和乙醇濃度梯度逐漸減小。由于液滴內部的熱擴散速率大于傳質擴散系數，內部溫度隨時間的變化比濃度隨時間的變化更快。

關鍵詞：液滴；蒸發；傳熱；傳質；數值模擬

降壓環境下液滴的蒸發過程，是一種伴隨相變的傳熱傳質過程，該過程能夠被廣泛應用于噴霧干燥、快速蒸餾和工業脫鹽等領域。實驗研究方面，OWEN和JALIL[1]采用液滴懸掛法記錄了閃蒸過程中的環境壓力和液滴的溫度變化。SATOH[2]、劉偉民[3]等實驗研究了純水液滴在閃蒸過程中的形態變化和溫度變化。理論研究方面，SHIN[4]和KIM[5]等建立了一種擴散控制模型對純水液滴快速蒸發時的傳熱傳質過程進行模擬。章學來等[6]通過建立集總熱力學模型，研究了液滴閃蒸和結冰過程中的溫度變化。

純水液滴降壓蒸發過程的實驗與理論研究已經得到了相當多的關注，但對于雙組分液滴蒸發的研究卻相對較少。駱騫[7]、高文忠[8]等實驗測得了鹽水（NaCl溶液）液滴在降壓環境下蒸發過程的溫度變化，并分析了其影響因素。劉偉民[9]測量了乙醇溶液液滴和煤油液滴降壓蒸發時的溫度變化。MUTHUNAYAGAM等[10]建立了蒸發擴散模型，研究降壓環境下鹽水液滴直徑和溫度的變化。在本文作者課題組之前的研究中，實驗研究了鹽水液滴的降壓蒸發過程，討論了液滴內部的傳質擴散和溫度梯度，理論分析了鹽水液滴溫度變化的影響因素[11-13]。

乙醇溶液液滴的蒸發過程比純水液滴和鹽水液滴更為復雜。因為在純水液滴和鹽水液滴中，都只有水分的蒸發，但是在乙醇溶液中，兩種組分的蒸發會同時進行，兩種組分的蒸發互相影響，增加了問題的復雜程度。一些學者[14-17]考慮了液滴內部溫度變化和液相傳質擴散以及活度系數的影響，建立了高溫環境下雙組分液滴加熱和蒸發的數學模型。然而，液滴內部的溫度與濃度分布以及液滴在降壓環境下的蒸發特性尚未提及。

本文針對兩種組分同時蒸發的乙醇溶液液滴在降壓環境下蒸發的傳熱傳質過程建立數學模型。模型考慮了液滴內部的傳熱傳質、液滴表面的蒸發和對流換熱，并與實驗數據對比驗證了模型的有效性。通過模型計算，獲得了液滴內部溫度分布以及濃度分布隨時間的變化。

1 乙醇溶液液滴降壓蒸發數學模型

以液相的能量守恒和傳質擴散為基礎，建立了乙醇溶液液滴降壓蒸發過程的數學模型。模型假設如下：①液滴始終為球形，僅考慮液滴徑向溫度和濃度的變化；②氣相的傳熱傳質為準穩態；③忽略液滴表面與環境的輻射換熱。

1.1 數學模型

在球坐標下，液相的能量守恒關系式如式(1)。

式中，χ的范圍在1～2.72之間，由式(3)估算[18]。

建立在式(2)和式(3)基礎上的有效熱導率模型考慮了內部環流對液滴熱傳導的加強作用，見式(4)～式(6)。

初始條件

兩個邊界條件

式中，qm和qh分別是由于液滴表面蒸發和對流引起的熱流密度，W/m2，其計算式將在下文中給出。

表示液相質量分數的關系式如式(7)。

式中，i=1，2 代表了液滴內的某種組分；Yli是液滴內部第i種組分的液相質量分數；Deff是液相的有效傳質擴散系數，m2/s，可由式(8)估算。

式中，Dl是液相傳質擴散系數，對于乙醇和水組成的雙組分液滴，Dl=1.2×10?9m2/s。與式(3)相似，系數χD考慮了液滴內部環流的強化傳質作用，其值可由式(9)確定[12]。初始條件

兩個邊界條件

式中，Yli s是液滴表面第i種組分的液相質量分數；iε是第i種組分在總蒸發率中所占分數；是質量蒸發率，kg/s，第i種組分的質量蒸發率為式(13)[12]。

式中，Hi是第i種組分的傳質系數；ρvis和ρvia分別為液滴表面和環境中第i種組分的蒸氣密度，kg/m3。

根據式(13)，由于蒸發造成的熱流密度qm可由式(14)、式(15)計算。

式中，Li是第i種組分的汽化潛熱，J/kg。

液滴表面由于對流引起的熱流密度qh可由式(16)計算。

式中，h是傳熱系數，W/(m2·K)，其值可由Nusselt數確定；Td,s和Ta分別為液滴表面溫度和環境溫度，K。

1.2 數值計算方法

由于蒸發過程中液滴尺寸逐漸減小，需要進行坐標變換以消除界面移動的影響。通過引入無量綱參數η對rs( t )進行坐標轉換，η的計算式如

式(17)。

則能量式(1)和傳質擴散方程式(7)有如式(18)、式(19)的形式。

本文中采用0.01ms的時間步長，液滴徑向的節點數為100。

2 實驗系統

圖1 實驗裝置系統圖1—真空泵；2—連通閥；3—真空罐；4—電磁閥；5—測試罐；6—加熱裝置；7—熱電偶；8—壓力傳感器；9—高速攝像機；10—NI數據采集系統；11—計算機

為驗證本文模型的準確性，采用液滴懸掛法進行相關實驗，將實驗數據與模擬結果進行對比。實驗系統如圖1所示，實驗過程如下：首先通過真空泵對真空罐進行抽氣至預定壓力，同時通過加熱裝置將測試罐內環境溫度加熱至預定工況，將液滴懸掛于熱電偶上并密封測試罐。開啟電磁閥（完全開啟時間0.5s），由于真空罐的容積是測試罐容積的800倍，測試罐內壓力將迅速下降，此時液滴發生閃蒸。同時，數據采集系統采集到的壓力和溫度信號、高速攝像機采集到的圖像數據傳輸至計算機。

實驗過程中，采用直徑為0.1mm的銅-康銅T形熱電偶測量液滴中心溫度，其絕對測量誤差為±0.2℃；采用德國Optronis 品牌高速攝像機CamRecord 450記錄液滴形態變化，其拍攝速度可達2000幀/s；Rosemount 3051CA絕對壓力變送器用于測量測試罐內環境壓力，其測量精度為0.075%；美國NI公司生產的數據采集系統用于采集和傳輸數據信號。

3 結果與分析

3.1 模型驗證

為了驗證上述模型的準確性，圖2和圖3分別比較了50%乙醇+50%水（質量分數，下同）的乙醇溶液液滴和75%乙酸+25%水的乙酸溶液液滴降壓蒸發過程中液滴內溫度變化的模擬值與實驗值。實驗數據由上述實驗系統得出，模擬計算中壓力下降引起的氣流速度ua按照文獻[11]中的方法進行計算。圖2中的初始工況和環境參數為：液滴的初始溫度為20℃，液滴初始直徑為1.8mm，環境溫度為34℃，環境壓力在前2.5s內從0.1MPa下降至2000Pa，2.5s后環境壓力保持在2000Pa。圖3的初始工況和環境參數為：液滴的初始溫度為26℃，液滴初始直徑為1.9mm，環境溫度為29℃，環境壓力在前2.5s內從0.1MPa下降至2000Pa，2.5s后環境壓力保持在2000Pa。

圖2 50%乙醇+50%水的乙醇溶液液滴溫度變化的模擬結果與實驗數據比較

圖3 75%乙酸+25%水的乙酸溶液液滴溫度變化的模擬結果與實驗數據比較

圖4 降壓過程中的環境壓力與空氣流速的變化

圖5 50%乙醇+50%水的乙醇溶液液滴的溫度分布隨時間的變化

在圖2中，由于實驗過程中熱電偶節點靠近液滴的中心，所以實驗測得的溫度值與模擬結果中的中心溫度更接近。中心溫度的模擬值與液滴溫度的實驗值相比較，其最大誤差僅為10.8%。從圖3可以看出，液滴溫度的模擬值與實驗測得的溫度值有相同的變化趨勢，且實驗值在中心溫度模擬值與表面溫度模擬值之間。但是，實驗過程中熱電偶節點更加靠近液滴的表面，所以表面溫度模擬值更加接近于液滴溫度的實驗值。由此可見，圖2和圖3中的模擬結果與實驗數據吻合較好，證明了模型的有效性。

3.2 降壓蒸發過程中液滴內部的溫度與濃度分布

圖4顯示了降壓過程中環境壓力與空氣流速隨時間的變化。從圖4中可以看出，環境壓力在2.5s內從0.1MPa下降至2000Pa，之后維持恒定。在t=0.83s時空氣流速ua達到最大值9.63m/s，之后空氣流速逐漸降低，直至t=4s時，由于環境壓力不再變化，空氣流速降為零。圖5為50%乙醇+50%水的乙醇溶液液滴溫度分布隨時間變化的示意圖，液滴初始溫度為18℃，初始直徑為1.9mm，環境溫度為32℃，環境壓力的變化與圖4所示一致。在最初的2s內，由于強烈的空氣流動，液滴蒸發迅速，使得液滴表面溫度迅速下降，造成了液滴內部相當大的溫差。空氣流速越快，液滴內部的溫度梯度越大。直至t=4s時，空氣流速降為零，液滴的蒸發速度減緩，液滴內部的溫差逐漸減小。值得一提的是，當t=10s時，液滴表面溫度升高。這主要是因為隨著蒸發的進行，環境空氣中的蒸氣質量分數增大，并且Spalding傳質系數減小，使得蒸發引起的熱流量減小。另一方面，隨著液滴溫度的降低，由于對流引起的熱流量將會增大。當由于對流引起的熱流量的增量大于因為蒸發引起的熱流量的減少量時，液滴表面溫度將會增加。

圖6顯示了乙醇溶液液滴中乙醇濃度分布隨時間的變化。在快速降壓階段，因為乙醇有較強的揮發性，液滴表面的乙醇濃度快速降低，而此時液滴內部的乙醇濃度未發生變化。當t=1s時，空氣流速達到最大值，液滴內部的濃度梯度也達到最大。之后隨著蒸發速率的降低，液滴內部的濃度梯度也開始下降，這意味著液滴內部的乙醇濃度逐漸減小，而表面的乙醇濃度逐漸增大。

圖6 50%乙醇+50%水的乙醇溶液液滴中乙醇濃度分布隨時間的變化

圖5和圖6所展示的降壓蒸發過程中液滴內部的溫度梯度和濃度梯度都非常有意義，從中可以看出，液滴內部溫度隨時間的變化比濃度隨時間的變化更快。這是由于液相的擴散系數D數量級大約是10?9m2/s，而溶液的熱擴散率α數量級大約是10?7m2/s。因此，液滴內部的傳質速率比熱擴散速率要慢。

4 結 論

本文針對兩種組分同時蒸發的乙醇溶液液滴在降壓環境下蒸發的傳熱傳質過程建立了數學模型。模型考慮了液滴內部傳熱傳質以及液滴表面的蒸發和對流換熱，并通過實驗數據對模型進行了驗證。通過模型計算，獲得了液滴內部溫度分布以及濃度分布隨時間的變化，獲得如下主要結論。

（1）快速降壓階段，空氣流動較快，液滴蒸發迅速，使得液滴表面溫度迅速下降，液滴內部溫差較大；壓力穩定后，空氣流速降為零，液滴的蒸發速度減緩，液滴內部的溫差逐漸減小。

（2）快速降壓階段，由于乙醇工質的易揮發性，液滴表面的乙醇濃度快速降低，而此時液滴內部的乙醇濃度未發生變化。之后隨著蒸發速率的降低，液滴內部的濃度梯度也開始下降。

（3）由于液滴內部的熱擴散速率大于傳質擴散系數，液滴內部溫度隨時間的變化比濃度隨時間的變化更快。

參 考 文 獻

[1] OWEN I，JALIL J M. Heterogeneous flashing in water drops[J]. International Journal of Multiphase Flow，1991，17（5）：653-660.

[2] SATOH I，FUSHINOBU K，HASHIMOTO Y. Freezing of a water droplet due to evaporation-heat transfer dominating the evaporationfreezing phenomena and the effect of boiling on freezing characteristics[J]. International Journal of Refrigeration，2002，25 （2）：226-234.

[3] 劉偉民，畢勤成，劉璐，等. 低壓閃蒸液滴形態和溫度變化的研究[J]. 工程熱物理學報，2007，28（6）：957-960.

[4] SHIN H T，LEE Y P，JURNG J. Spherical-shaped ice particle production by spraying water in vacuum chamber[J]. Applied Thermal Engineering，2000，20（5）：439-454.

[5] KIM B S，SHIN H T，LEE Y P，et al. Study on ice slurry production by water spray[J]. International Journal of Refrigeration，2001，24（2）：176-184.

[6] 章學來，王為，李志偉，等. 靜止水滴真空閃蒸模型及實驗研究[J]. 工程熱物理學報，2012，33（8）：1419-1422.

[7] 駱騫，畢勤成，韓彥寧，等. 鹽水液滴降壓環境下蒸發過程[J]. 化工學報，2013，64（6）：2001-2006.

[8] 高文忠，時亞茹，韓笑生，等. 混合除濕鹽溶液液滴閃蒸機理[J]. 化工學報，2012，63（11）：3453-3459.

[9] 劉偉民. 真空閃蒸過程中換熱與流動機理的研究[D]. 西安：西安交通大學，2006.

[10] MUTHUNAYAGAM A E，RAMAMURTHI K R，PADEN J R. Modelling and experiments on vaporization of saline water at low temperatures and reduced pressures[J]. Applied Thermal Engineering，2005，25（5/6）：941-952.

[11] LIU L，MI M L. Theoretical investigation on rapid evaporation of a saline droplet during depressurization[J]. Microgravity Science and Technology，2014，25：295-302.

[12] 劉璐，王茉，劉琰，等. 鹽水液滴降壓蒸發析鹽過程傳熱傳質特性研究[J]. 化工學報，2015，66（7）：2426-2432.

[13] 劉璐，舒盼盼，王茉，等. 鹽水液滴降壓蒸發析鹽過程數值模擬研究[J]. 化工進展，2015，34（1）：49-53.

[14] SAZHIN S S，ELWARDANY A，KRUTITSKII P A，et al. A simplified model for bi-component droplet heating and evaporation[J]. International Journal of Heat and Mass Transfer，2010，53（21/22）：4495-4505.

[15] MAQUA C，CASTANET G，LEMOINE F. Bicomponent droplets evaporation：temperature measurements and modelling[J]. Fuel，2008，87：2932-2942.

[16] YARIN A L，BRENN G，RENSINK D. Evaporation of acoustically levitated droplets of binary liquid mixtures[J]. International Journal of Heat and Fluid Flow，2002，23：471-486.

[17] 田章福，陶玉靜，蘇凌宇，等. 多組分液滴蒸發過程的理論模型及試驗研究[J]. 推進技術，2006，27（6）：568-571.

[18] ABRAMZON B，SIRIGNANO W A. Droplet vaporization model for spray combustion calculations[J]. International Journal of Heat and Mass Transfer，1989，32：1605-1618.

研究開發

Heat and mass transfer characteristics of an ethanol solution droplet during depressurized evaporation process

WANG Mo，LIU Lu，WANG Pengcheng，MI Menglong，LIU Yanfeng
（School of Energy Power and Mechanical Engineering，North China Electric Power University，Baoding 071003，Hebei，China）

Abstract：A mathematical model was constructed to simulate the temperature change and the concentration variation of a single ethanol solution droplet during depressurized evaporation. The present model was based on energy conservation in liquid phase，as well as the species diffusion in liquid phase. Film model was used to calculate the droplet evaporation rate. The activity coefficient was introduced to calculate the partial vapour pressure at the droplet surface. Using the droplet suspension method，the droplet temperatures during depressurized evaporation of ethanol solution droplets and acetic acid solution droplets were recorded，respectively. The present model was verified by comparing the numerical results with the experimental data. By numerical calculations，the variations of temperature and concentration profile within a droplet were discussed. Results showed that during the fast pressure dropping stage，the droplet evaporation rate was rapid due to the strong air flow，leading to a fairly large temperature difference inside the droplet. The concentration gradient inside the droplet was also very great. When the ambient pressure maintained constant，the temperature difference inside the droplet diminished gradually，and the concentration gradient within the droplet was also decreased. Because the thermal diffusivity is greater than the mass transfer coefficient，thebook=718,ebook=71droplet internal temperature change with time is faster than the concentration change versus time.

Key words：droplet；evaporation；heat transfer；mass transfer；numerical simulation

基金項目：國家自然科學基金（51406055）、河北省自然科學基金（E2012502069）及中央高校基本科研業務費專項資金項目（2014MS106）。

收稿日期：2015-09-14；修改稿日期：2015-10-10。

DOI：10.16085/j.issn.1000-6613.2016.03.010

中圖分類號：TK 121

文獻標志碼：A

文章編號：1000–6613（2016）03–0717–05

第一作者：王茉（1990—），女，碩士研究生。聯系人：劉璐，講師，研究方向為相變傳熱傳質。E-mail 122446254@qq.com。