潤濕特性對多孔介質(zhì)緩慢干燥過程影響的數(shù)值模擬

程 樹，吳 睿，崔國民，李延銘

（上海理工大學(xué) 能源與動力工程學(xué)院新能源所，上海 200093）

潤濕特性對多孔介質(zhì)緩慢干燥過程影響的數(shù)值模擬

程 樹，吳 睿，崔國民，李延銘

（上海理工大學(xué) 能源與動力工程學(xué)院新能源所，上海 200093）

采用二維孔隙網(wǎng)絡(luò)模型，對不同孔隙結(jié)構(gòu)下親水多孔介質(zhì)與憎水多孔介質(zhì)的緩慢干燥過程進(jìn)行模擬，考察潤濕特性的改變對二維多孔介質(zhì)干燥過程中液相分布、蒸發(fā)速率以及干燥總時間的影響。模擬結(jié)果表明，忽略重力與黏性力影響時，多孔介質(zhì)干燥過程只受毛細(xì)壓力控制，親水多孔介質(zhì)蒸發(fā)端為無規(guī)則形狀，而憎水多孔介質(zhì)蒸發(fā)端則根據(jù)孔隙尺寸分布與喉孔比的不同呈現(xiàn)不同的特性；干燥總時間與潤濕特性、蒸發(fā)速率、孔隙結(jié)構(gòu)相關(guān)；當(dāng)喉孔比相同時，緊湊多孔介質(zhì)的干燥總時間小于非緊湊多孔介質(zhì)；當(dāng)喉孔比較小時，親水多孔介質(zhì)的干燥總時間小于憎水多孔介質(zhì)；當(dāng)喉孔比較大時，憎水與親水多孔介質(zhì)的干燥總時間可以近乎相等。

多孔介質(zhì)；干燥；孔隙網(wǎng)絡(luò)模型；潤濕特性；喉孔比

多孔介質(zhì)干燥是工業(yè)生產(chǎn)中的一個重要環(huán)節(jié)，主要用于食品干燥、輕質(zhì)油回收、質(zhì)子交換膜燃料電池等方面［1-2］。在一些領(lǐng)域中，多孔介質(zhì)表現(xiàn)為親水特性，干燥過程類似于排水過程［3-4］。而在一些特殊領(lǐng)域如質(zhì)子交換膜燃料電池氣體擴(kuò)散層等，多孔介質(zhì)則表現(xiàn)為憎水特性，其干燥過程可類比為自吸過程。

孔隙網(wǎng)絡(luò)模型作為一種有效方法可應(yīng)用于研究多孔介質(zhì)的干燥過程。自Daian等［5］和Nowicki等［6］率先提出孔隙網(wǎng)絡(luò)模型以來，Prat等［7-8］結(jié)合入侵逾滲理論提出了多孔介質(zhì)干燥孔隙網(wǎng)絡(luò)模型。此后，重力［9-10］、黏性力［11-13］、熱效應(yīng)［14-16］、拓?fù)浣Y(jié)構(gòu)［17-21］等影響多孔介質(zhì)干燥的外界因素在建模過程中被考慮。另一方面，針對干燥過程中出現(xiàn)的液膜效應(yīng)［22-23］、收縮現(xiàn)象［24］進(jìn)行的研究進(jìn)一步完善了多孔介質(zhì)干燥的孔隙網(wǎng)絡(luò)模型。

近年來，Pillai等［19,25］對親水多孔介質(zhì)干燥過程進(jìn)行了研究。但對憎水多孔介質(zhì)的研究還較少。吳睿等［26］采用三維孔隙網(wǎng)絡(luò)模型對質(zhì)子交換膜燃料電池氣體擴(kuò)散層水分蒸發(fā)的干燥過程進(jìn)行模擬。近期，Wu等［27］采用三維孔隙網(wǎng)絡(luò)模型研究孔隙結(jié)構(gòu)對憎水多孔介質(zhì)蒸發(fā)過程的影響，指出孔隙結(jié)構(gòu)是影響憎水多孔介質(zhì)干燥特性的重要因素。Chapuis等［28］在研究相同孔隙結(jié)構(gòu)多孔介質(zhì)蒸發(fā)過程時發(fā)現(xiàn)，憎水多孔介質(zhì)蒸發(fā)端相對平整，而親水多孔介質(zhì)蒸發(fā)端為不規(guī)則狀態(tài)。Lenormand［29］指出在不同孔喉比的情況下，自吸過程的兩相流動特性截然不同。由此可見，潤濕性對多孔介質(zhì)的影響與其孔隙結(jié)構(gòu)密切相關(guān)。但潤濕性的變化對不同孔隙結(jié)構(gòu)多孔介質(zhì)干燥的影響尚未明確。

本工作對不同孔隙結(jié)構(gòu)多孔介質(zhì)潤濕特性的影響進(jìn)行研究，考察相位分布、蒸發(fā)速率、蒸發(fā)總時間對多孔介質(zhì)干燥特性的影響。

1 孔隙網(wǎng)絡(luò)模型

孔隙網(wǎng)絡(luò)模型是將多孔介質(zhì)抽象為理想的幾何形狀，復(fù)雜的孔隙空間由喉道及其相連的孔隙構(gòu)成，喉道代表較小的空間，孔隙代表喉道交接處較大的空間。采用二維方形孔隙網(wǎng)絡(luò)模型，其中，孔隙為正方形，連接兩個相鄰孔隙的為矩形截面的喉道，二維孔隙的網(wǎng)絡(luò)模型見圖1。

孔隙半徑為孔隙邊長的一半，取值區(qū)間為[rpmin，rpmax]：

式中，ψ為0～1之間的隨機(jī)數(shù)。

連接孔隙i和j的喉道半徑（rt）為喉道橫截面邊長的一半，表示為

式中，λ為喉孔比，是孔隙網(wǎng)絡(luò)模型的結(jié)構(gòu)參數(shù)，通常情況下λ≤1。相鄰孔隙之間的距離（l）為定值，設(shè)為25 μm。另外，孔隙網(wǎng)絡(luò)模型中孔隙網(wǎng)絡(luò)大小也是影響多孔介質(zhì)蒸發(fā)過程的一個重要因素，為了更為詳細(xì)地觀察蒸發(fā)過程中孔隙網(wǎng)絡(luò)中的兩相流動特性，采用較小的孔隙網(wǎng)絡(luò)（設(shè)定為20×20）。

圖1 二維孔隙的網(wǎng)絡(luò)模型Fig.1 2D pore-network models.

由于只考慮緩慢蒸發(fā)過程，因此兩相流動只受毛細(xì)力的控制，黏性力的影響則可忽略。在此情況下，多孔介質(zhì)蒸發(fā)過程的兩相流動可以視為準(zhǔn)靜態(tài)過程。蒸發(fā)過程中，每個孤立液體團(tuán)中門闕流動阻力最小的單元（孔隙或喉道）將被驅(qū)替。單元的門闕流動阻力與單元尺寸和潤濕性相關(guān)。以憎水多孔網(wǎng)絡(luò)為例，喉道門闕流動阻力（Pct）［29］可由式（3）計算。

式中，σ為表面張力，取值0.072 N/m；θ為接觸角，接觸角的大小與介質(zhì)的潤濕性有關(guān)，憎水多孔介質(zhì)θ取值60°，在潤濕性不變的情況下θ不變。

孔隙門闕流動阻力（Pcp）［27］可由式（4）計算。

式中，rpe為有效孔隙半徑；z為孔隙周圍被氣相占據(jù)的喉道數(shù)，z=1～3。親水Pcp的計算式與式（3）相同。

由于孔隙周圍被氣相占據(jù)的喉道數(shù)目不同，根據(jù)式（4）計算的Pcp有3種不同的結(jié)果，故喉道和與之相鄰孔隙之間的門闕流動阻力根據(jù)z的不同而改變。這與親水多孔網(wǎng)絡(luò)有著根本的不同。親水多孔網(wǎng)絡(luò)中，親水多孔介質(zhì)θ的取120°，計算孔隙阻力時均采用實際尺寸而非有效孔隙半徑［30］，由于喉道尺寸小于相鄰孔隙，故Pcp大于相鄰孔隙。一旦喉道被侵入則相鄰孔隙自動被侵入。

初始狀態(tài)時孔隙網(wǎng)絡(luò)完全充滿液態(tài)水，水蒸氣擴(kuò)散只通過不封閉的上界面進(jìn)行。在計算中做如下簡化假設(shè)：1）忽略黏性力和重力的作用；2）水蒸氣擴(kuò)散過程為準(zhǔn)靜態(tài)過程；3）整個干燥過程在恒溫條件下進(jìn)行；4）水蒸氣視為理想氣體。

由以上假設(shè)條件，位于模型開口界面處的喉道與外界環(huán)境相接觸，由菲克定律可得：

孔隙網(wǎng)絡(luò)模型內(nèi)部被氣相占據(jù)的孔隙j的水蒸氣蒸發(fā)量可由質(zhì)量守恒原理得到：

式中，qpitij為孔隙i到喉道ij（連接孔隙i與孔隙j的喉道）的水蒸氣蒸發(fā)量，可由菲克定律得到：

若孔隙或喉道兩者中有一項為液相，那么該項的水蒸氣濃度等于飽和水蒸氣濃度，即cpi=1或ctij=1。

通過以上分析，多孔介質(zhì)干燥過程的孔隙網(wǎng)絡(luò)模型模擬算法可歸納為：1）確定網(wǎng)絡(luò)模型中的液體團(tuán)；2）根據(jù)質(zhì)量守恒定理迭代計算孔隙網(wǎng)絡(luò)中每個被氣相占據(jù)的孔隙中水蒸氣的擴(kuò)散量；3）確定每個液體團(tuán)中具有最小門闕流動阻力的單元；4）通過計算液體團(tuán)中每個含有液相的單元和與其比鄰的氣相單元之間的擴(kuò)散量來確定每個液體團(tuán)的總蒸發(fā)量；5）確定每個液體團(tuán)中具有最小門闕流動阻力單元中液體被完全蒸發(fā)所需時間，并確定最小時間；6）根據(jù)第5）步中確定的最小時間所在的單元更新計算網(wǎng)絡(luò)模型的相分布；7）重復(fù)1）～6）步。

2 結(jié)果與討論

主要分析潤濕性對不同結(jié)構(gòu)多孔介質(zhì)干燥過程的影響。多孔介質(zhì)的潤濕性可由液相接觸角表征。多孔介質(zhì)結(jié)構(gòu)用孔隙分布和喉孔比來表征。構(gòu)建兩種多孔網(wǎng)絡(luò)：一為緊湊多孔網(wǎng)絡(luò)，其孔隙半徑分布區(qū)間較小，為5.5～6.5 μm；二為非緊湊多孔網(wǎng)絡(luò)，其孔隙半徑分布區(qū)間較大，為0～12 μm。此外，通過改變多孔網(wǎng)絡(luò)的喉孔比（0.1，0.8，1.0），以獲得潤濕性對不同孔隙結(jié)構(gòu)多孔介質(zhì)干燥過程的影響。

2.1 液相分布

多孔介質(zhì)干燥特性與其內(nèi)部液相分布相關(guān)。為了揭示潤濕性對多孔介質(zhì)干燥過程的影響，對比不同液相飽和度情況下多孔介質(zhì)干燥過程的液相分布。液相飽和度為含有液體的單元體積與孔隙網(wǎng)絡(luò)中所有單元體積之比。

緊湊多孔介質(zhì)的干燥過程見圖2。由圖2可知，喉孔比較小的憎水多孔介質(zhì)，蒸發(fā)過程呈現(xiàn)幾乎平整蒸發(fā)端，即蒸發(fā)過程逐行進(jìn)行。多孔介質(zhì)侵入過程只由毛細(xì)壓力控制，因此比較Pct和Pcp能進(jìn)一步解釋蒸發(fā)侵入過程。憎水多孔介質(zhì)侵入過程中，喉道和周圍被氣相占據(jù)喉道數(shù)目不同的孔隙的門闕流動阻力大小排序為：

式中，Pcp（1），Pcp（2），Pcp（3）為孔隙周圍被氣相占據(jù)喉道數(shù)目分別為1，2，3的Pcp。因此，對于喉孔比較小的憎水緊湊多孔介質(zhì)，門闕流動阻力小的表面喉道優(yōu)先被侵入，第一層喉道完全侵入后，位于第二層的孔隙才開始蒸發(fā)，而直到第二層孔隙蒸發(fā)完全后蒸發(fā)過程才在第三層喉道中進(jìn)行，故蒸發(fā)端呈平整下降。而當(dāng)喉孔比為1時，孔隙與喉道尺寸差距小且連接相鄰兩孔隙的喉道尺寸以較小喉道尺寸為基準(zhǔn)，故喉道與孔隙之間的門闕流動阻力大小順序為：Pct≥Pcp。此時，侵入競爭發(fā)生在整個多孔介質(zhì)空間不同尺寸的喉道與孔隙之間，其干燥相圖表現(xiàn)為毛細(xì)指徑現(xiàn)象。而對于親水多孔介質(zhì)，蒸發(fā)相圖均呈現(xiàn)為毛細(xì)指徑現(xiàn)象，且喉孔比的變化對蒸發(fā)過程的影響不明顯。

非緊湊多孔介質(zhì)的干燥過程見圖3。在非緊湊多孔介質(zhì)空間內(nèi)，尺寸相差較大的孔隙和喉道呈隨機(jī)分布。根據(jù)入侵逾滲理論，憎水多孔介質(zhì)蒸發(fā)過程選擇侵入門闕流動阻力最小的孔隙或喉道，由于這個孔隙和喉道可能出現(xiàn)在孔隙網(wǎng)絡(luò)模型中的任意一處，所以無論多孔介質(zhì)的潤濕性如何，干燥相圖都呈現(xiàn)出毛細(xì)指徑現(xiàn)象，而且相同飽和度下喉孔比的變化對蒸發(fā)過程的影響不大。

對比圖2和圖3中憎水多孔介質(zhì)的蒸發(fā)相圖可以發(fā)現(xiàn)：喉孔比較小的憎水緊湊多孔介質(zhì)蒸發(fā)相圖中只在最后一行邊界處存在孤立的孔隙；喉孔比較大的緊湊多孔介質(zhì)相圖中既存在孤立喉道又存在孤立的孔隙。而非緊湊多孔介質(zhì)中，喉孔比較小的多孔介質(zhì)干燥過程中只出現(xiàn)孤立的孔隙，此孤立孔隙可以出現(xiàn)在孔隙網(wǎng)絡(luò)中的任意位置；喉孔比較大的多孔介質(zhì)干燥過程中既有孤立的孔隙又存在孤立的喉道。由此可見，孔隙結(jié)構(gòu)不同的多孔介質(zhì)干燥過程中相位分布存在著差異，且喉孔比對憎水多孔介質(zhì)相位分布的影響更大。

圖2 緊湊多孔介質(zhì)的干燥過程Fig.2 Drying process of porous media with narrow pore size distribution(PSD). Points represented the pores dominated by liquid， and lines denoted throats， the same below.

圖3 非緊湊多孔介質(zhì)的干燥過程Fig.3 Drying process of porous media with wide PSD.

為解釋這一現(xiàn)象，取憎水孔隙網(wǎng)絡(luò)的某孤立小液團(tuán)部分進(jìn)行具體分析。無孤立喉道氣相侵入的原理見圖4。對于喉孔比小于1的多孔介質(zhì)，無論孔隙尺寸如何分布，與孔隙相連的喉道尺寸一定小于孔隙尺寸（見圖4a）。由于憎水多孔介質(zhì)的侵入過程發(fā)生在最小門闕流動阻力處，所以蒸發(fā)首先在喉道t1中進(jìn)行，此后按門闕流動阻力等級依次在t2，t3，t4中進(jìn)行。隨著t4的蒸發(fā)完全，孔隙網(wǎng)絡(luò)中出現(xiàn)了孤立的孔隙p1（p1為微元中最大的孔隙，見圖4b）。憎水緊湊多孔介質(zhì)當(dāng)蒸發(fā)進(jìn)行到最后一行時，位于邊界處的Pcp（1）大于其他位置的Pcp（2），故在蒸發(fā)相圖中為孤立孔隙，且只存在孤立孔隙。對于喉孔比等于1的多孔介質(zhì)，孔隙尺寸分布不同的多孔介質(zhì)存在著不同的現(xiàn)象。

圖4 無孤立喉道氣相侵入的原理Fig.4 Schematic diagram of gas introduction into un-isolated throats.

無孤立孔隙氣相侵入的原理見圖5。由于緊湊多孔介質(zhì)的相鄰孔隙的尺寸差距較小，假設(shè)p2和p3為微元中的最大孔隙，且rp2=rp3（見圖5a）則rp2=rp3=rt23，兩孔隙周圍與空氣接觸的喉道t1～t6先于孔隙p2，p3，t23被侵入（見圖5b）。之后，蒸發(fā)在孔隙p2和p3間隨機(jī)選擇，若p2被侵入則孔隙網(wǎng)絡(luò)中留有孔隙p3和喉道t23（見圖5c）。此時，孔隙p3周圍被氣相占據(jù)的喉道數(shù)為1，侵入過程隨之隨機(jī)發(fā)生在孔隙p3和喉道t23之間，若p3被侵入，則孔隙網(wǎng)絡(luò)中出現(xiàn)孤立的喉道（見圖5d），若t23被侵入，則出現(xiàn)孤立孔隙。另外，非緊湊多孔介質(zhì)，由于其孔隙尺寸取值范圍較大，上述兩種情況皆可能出現(xiàn)。

當(dāng)小液團(tuán)中只有一個最大尺寸孔隙時，蒸發(fā)過程如圖4所示，而小液團(tuán)中具有兩個相同尺寸的孔隙時，蒸發(fā)過程則如圖5所示。

圖5 無孤立孔隙氣相侵入的原理Fig.5 Schematic diagram of gas introduction into un-isolated pores.

2.2 蒸發(fā)速率

潤濕性與孔隙結(jié)構(gòu)的變化使得多孔介質(zhì)中液相分布趨勢隨之改變，而這一改變直接影響多孔介質(zhì)液相飽和度的變化趨勢。給出不同孔隙結(jié)構(gòu)多孔介質(zhì)的液體飽和度-時間曲線以揭示液體飽和度隨蒸發(fā)時間的變化趨勢。多孔介質(zhì)干燥過程中液體飽和度隨蒸發(fā)時間的變化關(guān)系見圖6。

從圖6可看出，無論多孔介質(zhì)的孔隙尺寸分布是否緊湊，小喉孔比多孔介質(zhì)的液相飽和度緩慢下降，而大喉孔比多孔介質(zhì)的液相飽和度則極速下降。

圖6 多孔介質(zhì)干燥過程中液相飽和度隨蒸發(fā)時間的變化趨勢Fig.6 Variation of S with evaporation time.

由液相飽和度與時間的關(guān)系可以得到蒸發(fā)速率的變化趨勢。在不同喉孔比條件下不同多孔介質(zhì)中的無量綱蒸發(fā)速率的結(jié)果見表1。從表1可看出，無量綱蒸發(fā)速率基本上隨喉孔比的增加而減小。通過計算一個微元網(wǎng)絡(luò)中瞬時蒸發(fā)速率來解釋多孔介質(zhì)蒸發(fā)速率的變化規(guī)律。微元網(wǎng)絡(luò)的蒸發(fā)原理見圖7。

表1 在不同喉孔比條件下緊湊和非緊湊多孔介質(zhì)中的無量綱蒸發(fā)速率Table 1 Non-dimensional evaporation rate(Mv) in the narrow and wide PSD networks at different λ

圖7是一個1×1的二維孔隙網(wǎng)絡(luò)，除上界面外其他3個界面皆封閉。假設(shè)初始狀態(tài)孔隙與喉道完全充滿液體，表面喉道中的液體與外界相連。隨著蒸發(fā)的進(jìn)行，蒸發(fā)端下降至孔隙中點位置。由菲克定律可知，初始狀態(tài)微元網(wǎng)絡(luò)的蒸發(fā)量（q）為：

隨著蒸發(fā)的進(jìn)行，當(dāng)孔隙中一半的液體已被蒸發(fā)，此時的孔隙中的蒸發(fā)量qpt1為：

喉道與外界的蒸發(fā)量qt13為：

由式（14）可知，J隨喉孔比的增大而減小。

圖7 微元網(wǎng)絡(luò)的蒸發(fā)原理Fig.7 Schematic diagram of the evaporation for microelement network.

2.3 干燥總時間

喉孔比對干燥總時間的影響見圖8。從圖8可知，在相同的喉孔比條件下，緊湊多孔介質(zhì)的干燥總時間小于非緊湊多孔介質(zhì)。多孔介質(zhì)干燥總時間本質(zhì)上由多孔介質(zhì)空隙體積與蒸發(fā)速率二者聯(lián)合競爭決定。

多孔介質(zhì)空隙體積為多孔介質(zhì)內(nèi)部孔隙與喉道體積之和。非緊湊多孔介質(zhì)的空隙體積大于緊湊多孔介質(zhì)，則在相同的喉孔比條件下，緊湊多孔介質(zhì)的所需干燥總時間比非緊湊多孔介質(zhì)更短。

由圖8可見，隨喉孔比的增大，緊湊多孔介質(zhì)和非緊湊多孔介質(zhì)的干燥總時間均呈現(xiàn)出先陡降再緩降最后趨于平衡的趨勢；隨喉孔比的增大，親水與憎水多孔介質(zhì)的干燥時間之間的差距逐漸縮小。其原因在于喉孔比較小時，親水多孔介質(zhì)在干燥初期開口上界面還留有部分與外界相連的喉道，這些表面喉道為干燥過程提供了較高的蒸發(fā)速率，從而縮短了蒸發(fā)時間，而憎水多孔介質(zhì)在干燥初期與外界相連的表面喉道已被蒸發(fā)，所以其蒸發(fā)速率較小，從而干燥完全所需時間較長。當(dāng)喉孔比較大時，親水與憎水多孔介質(zhì)干燥初期都存在部分表面喉道，整個干燥過程中的平均干燥速率相似，所以它們的干燥總時間接近。

圖8 喉孔比對干燥總時間的影響Fig.8 Effects of different λ on the total drying time.

由圖8還可見，當(dāng)喉孔比較小時，親水多孔介質(zhì)的干燥過程比憎水多孔介質(zhì)快，這一現(xiàn)象與文獻(xiàn)［29］給出的結(jié)論相同；當(dāng)喉孔比較大時，親水與憎水緊湊多孔介質(zhì)干燥總時間幾乎相等，其原因在于喉孔比為0.1～0.8時，憎水緊湊多孔介質(zhì)蒸發(fā)逐行進(jìn)行，表面喉道快速蒸發(fā)階段持續(xù)時間短暫，此后蒸發(fā)速率降低，因而所需干燥時間較長，而親水多孔介質(zhì)在干燥初期由于表面喉道蒸發(fā)階段持續(xù)時間較長，所以蒸發(fā)速率快。而在干燥后期，多孔介質(zhì)被突破后，主液體團(tuán)被撕裂成許多小液團(tuán)，從而引起了較高的蒸發(fā)速率。因此，親水多孔介質(zhì)的蒸發(fā)速率總大于憎水多孔介質(zhì)。當(dāng)喉孔比為0.9～1.0時，憎水多孔介質(zhì)的蒸發(fā)不再逐行進(jìn)行，而是呈現(xiàn)與親水多孔介質(zhì)類似的毛細(xì)指徑現(xiàn)象，且在達(dá)到突破前已經(jīng)被撕裂成許多小液團(tuán)。小液團(tuán)蒸發(fā)速率較快，造成憎水多孔介質(zhì)的蒸發(fā)速率非常接近于親水多孔介質(zhì)，從而使得兩種不同潤濕性多孔介質(zhì)的干燥時間相近。

3 結(jié)論

1） 憎水緊湊多孔介質(zhì)，喉孔比較小時蒸發(fā)端平整，喉孔比較大時蒸發(fā)端無規(guī)則形狀；不同喉孔比的親水緊湊多孔介質(zhì)蒸發(fā)端形狀不規(guī)則。非緊湊多孔介質(zhì)，潤濕性的改變對蒸發(fā)端的影響較小。

2）喉孔比較大的多孔介質(zhì)干燥過程中既存在孤立的孔隙又存在孤立的喉道。當(dāng)喉孔比較小時，憎水緊湊多孔介質(zhì)蒸發(fā)相圖中只在最后一行邊界處存在孤立的孔隙，非緊湊多孔介質(zhì)中孤立孔隙可出現(xiàn)在孔隙網(wǎng)絡(luò)中的任意位置。

3） 隨喉孔比的增大，親水和憎水的多孔介質(zhì)蒸發(fā)速率均減小。

4） 當(dāng)喉孔比相同時，緊湊多孔介質(zhì)的干燥總時間小于非緊湊多孔介質(zhì)；當(dāng)喉孔比較小時，親水多孔介質(zhì)的干燥總時間小于憎水多孔介質(zhì)；當(dāng)喉孔比較大時，憎水與親水多孔介質(zhì)的干燥總時間可以近乎相等。

［1］ 陳永利，曹立勇，何威，等. 分形多孔介質(zhì)的孔隙特性對氣體擴(kuò)散的影響［J］. 化工學(xué)報，2011， 62（11）：3024 -3029.

［2］ 馬亮，何榕. 分形多子L介質(zhì)中氣體的非穩(wěn)態(tài)擴(kuò)散［J］. 化工學(xué)報，2013， 64（9）：3139 - 3144.

［3］ Prat M. Recent Advances in Pore-Scale Models for Drying of Porous Media［J］.Chem Eng J，2002，86（1/2）：153 - 164.

［4］ Metzger T，Tsotsas E. Network Models for Capillary Porous Media［J］.Chem Ing Tech，2010，82（6）：869 - 879.

［5］ Daian J F，Saliba J. Determining a Representative Pore-Network for Moisture Sorption and Migration in Cement Mortar［J］.Int J Heat Mass Transfer，1991，34（8）：2081 - 2096.

［6］ Nowicki S C，Davis H T，Scriven L E. Microscopic Determination of Transport Parameters in Drying Porous Media［J］.Dry Technol，1992，10（4）：926 - 946.

［7］ Prat M. Percolation Model of Drying Under Isothermal Conditions Inporous Media［J］.Int J Multiphase Flow，1993，19（4）：691 - 704.

［8］ Tsimpanogiannis I N， Yortsos Y C， Poulou S， et al. Scaling Theory of Drying in Porous Media［J］.Phys Reν E：Stat，Nonlinear，Soft，Matter Phys，1999，59 （4）：4353 - 4365.

［9］ Laurindo J B，Prat M. Numerical and Experimental Network Study of Evaporation in Capillary Porous Media：Phase Distribution［J］.Chem Eng Sci，1996，51（23）：5171 - 5185.

［10］ Yiotis A G，Tsimpanogiannis I N，Stubos A K，et al. Pore Network Study of the Characteristic Speriods in the Drying of Porous Materials［J］.J Colloid Interface Sci，2006，297（2）：738 - 748.

［11］ Yiotis A G，Stubos A K，Boudouvis A G， et al. A 2-D Pore-Network Model of the Drying of Single Component Liquids in Porous Media［J］.Adν Water Res，2001，24（3/4）：439 -460.

［12］ Metzger T，Irawan A，Tsotsas E. Isothermal Drying of Pore Networks：Influence of Friction for Differentpore Structures［J］.Dry Technol，2007，25（1/3）：49 - 57.

［13］ Metzger T，Tsotsas E. Viscous Stabilization of Drying Front：Three-Dimensional Pore Network Simulations［J］.Chem EngRes Des，2008，86（7A）：739 - 744.

［14］ Plourde F，Prat M. Pore Network Simulations of Drying of Capillary Porous Media. Influence of Thermal Gradients［J］.Int J Heat Mass Transfer，2003，46（7）：1293 - 1307.

［15］ Surasani V K，Metzger T，Tsotsas E. Consideration of Heat Transfer in Pore Network Modeling of Convective Drying［J］.Int J Heat Mass Transfer，2008，51（9/10）：2506 - 2518.

［16］ Surasani V K，Metzger T，Tsotsas E. Influence of Heating Mode on Drying Behavior of Capillary Porous Media：Pore Scale Modeling［J］.Chem Eng Sci，2008，63（21）：5218 -5228.

［17］ Metzger T，Irawan A，Tsotsas E. Inf uence of Pore Structure on Drying Kinetics：A Pore Network Study［J］.AIChE J，2007，53（12）：3029 - 3041.

［18］ Surasani V K，Metzger T，Tsotsas E. Drying Simulations of Various 3D Pore Structure by a Nonisothermal Pore Network Model［J］.Dry Technol，2010，28（5）：615 - 623.

［19］ Pillai K M，Prat M，Marcoux M. A Study on Slow Evaporation of Liquids in a Dual-Porosity Porous Medium Using Square Network Model［J］.Int J Heat Mass Transfer，2009，52 （7/8）：1643 - 1656.

［20］ Segura L A，Toledo P G. Pore-Level Modeling of Isothermal Drying of Pore Networks：Effects of Gravity and Pore Shape and Size Distribution on Saturation and Transport Parameters［J］.Chem Eng J，2005，111（2/3）：237 - 252.

［21］ Yiotis A G，Tsimpanogiannis I N，Stubos A K. Three-Dimensional Modeling of the Evaporation of Volatile Hydrocarbons from Anisotropic Porous Media［J］.J Porous Media，2010， 13（3）：209 - 219.

［22］ Yiotis A G，Boudouvis A G，Stubos A K，et al. Effect of Liquid Films on the Drying of Porous Media［J］.AIChE J，2004，50（11）：2721 - 2737.

［23］ Yiotis A G，Stubos A K，Boudouvis A G，et al. Pore Network Modeling of Isothermal Drying in Porous Media［J］.Trans Porous Media，2005，58 （1/2）：63 - 86.

［24］ Thomas Metzger，Evangelos Tsotsas. Network Models for Capillary Porous Media：Application to Drying Technology［J］.Chem Ing Tech，2010，82 （6）：869 - 879.

［25］ Sghaier N，Prat M，Nasrallah S B. On the Infuence of Sodium Chloride Concentration on Equilibrium Contact Angle［J］.Chem Eng J，2006，122（1/2）：47 - 53.

［26］ 吳睿，朱恂，廖強(qiáng)，等. PEMFC氣體擴(kuò)散層干燥過程孔隙網(wǎng)絡(luò)模擬［J］. 工程熱物理學(xué)報，2011，32（1）：97 - 100.

［27］ Wu Rui，Cui Guomin，Chen Rong. Pore Network Study of Slow Evaporation in Hydrophobic Porous Media［J］.Int J Heat Mass Transfer，2014，68（1）：310 - 323.

［28］ Chapuis O，Prat M. Influence of Wettability Conditions on Slow Evaporation in Two-Dimensional Porousmedia［J］.Phys Reν E：Stat，Nonlinear，Soft，Matter Phys，2007，75（4）：046311.

［29］ Lenormand R. Liquids in Porous Media［J］.J Phys：Condens Matter，1990，2（1）：79 - 88.

［30］ Prat M. Discrete Models of Liquid-Vapour Phase Change Phenomena in Porous Media［J］.Reνue Générale Thermique，1998，37（11）：954 - 961.

（編輯 李治泉）

Numerical Simulation of Influence of Wetting Property on Slow Drying of Porous Media

Cheng Shu, Wu Rui, Cui Guomin, Li Yanming
（Institute of New Energy Science and Engineering，School Energy and Power Engineering，University of Shanghai for Science and Technology，Shanghai 200093，China）

The slow drying processes of hydrophilic and hydrophobic porous media with different pore structure were simulated by using a pore network model and the inf uences of wetting property on evaporation behavior were investigated. The simulation results indicated that when adhesion and gravity were ignored，capillary pressure dominated the drying processes; the drying front ends in the hydrophilic porous media kept irregular shape，but in the hydrophobic porous media the drying front ends showed different features according to the pore size distribution(PSD) and throat-pore ratio; the total drying time was affected by the evaporation rate，pore structure and wetting property. For the wide PSD network，the total drying time of the hydrophilic porous media is less than that of the hydrophobic porous media when the throat-pore ratio is low. For the narrow PSD，the total drying time of the hydrophobic porous media is similar to that of the hydrophilic porous media when the throat-pore ratio is high.

porous media；evaporation；pore network；wetting property；throat-pore ratio

1000 - 8144（2014）09 - 1035 - 09

TQ 021.4

A

2014 - 03 - 27；［修改稿日期］ 2014 - 06 - 06。

程樹（1988—），女，浙江省安吉市人，碩士生，電郵 shucheng1022@gmail.com。聯(lián)系人：吳睿，電話 021 -55270155，電郵 ruiwu1986@gmail.com。

國家自然科學(xué)基金青年自然基金項目（51306124）；上海市自然科學(xué)基金青年項目（13ZR1458300）；重慶大學(xué)教育部重點實驗室開放基金項目（LLEUTS-201305）；滬江基金研究基地專項（D14001）。