宋粉紅 姜瑜輝 王志遠(yuǎn) 王忠旭 范晶

收稿日期：2023-06-08

基金項(xiàng)目：吉林省自然科學(xué)基金項(xiàng)目（20210101088JC）

第一作者：宋粉紅（1983-），女，教授，博士，研究方向?yàn)樾履茉撮_發(fā)與存儲(chǔ)、二氧化碳捕集、潤濕與電潤濕、微尺度傳熱傳質(zhì)等。E-mail： fenhongsong@neepu.edu.cn。

通信作者：范晶（1984-），女，教授，博士，研究方向?yàn)樘娲剂蠠嵛镄浴⒌推肺荒芾?、二氧化碳捕集等。E-mail： crystalfan@neepu.edu.cn。

文章編號(hào)：1673-5005（2024）02-0161-09??? doi：10.3969/j.issn.1673-5005.2024.02.018

摘要： CO2在粗糙固體壁面上的潤濕行為對CO2冷凝捕集技術(shù)具有重要意義。采用分子動(dòng)力學(xué)方法探究微結(jié)構(gòu)、勢能及溫度對CO2液滴潤濕的影響。結(jié)果表明：相較于平板壁面，微結(jié)構(gòu)壁面提高了CO2液滴的接觸角；在微結(jié)構(gòu)壁面上勢能參數(shù)影響CO2液滴的潤濕狀態(tài)，在欄柵結(jié)構(gòu)壁面上CO2液滴從Cassie-Baxter狀態(tài)轉(zhuǎn)向Wenzel狀態(tài)，而在三角和斜坡結(jié)構(gòu)壁面上CO2液滴均保持Wenzel狀態(tài)，3種壁面上CO2液滴的接觸角均隨勢能參數(shù)的增大而減小；此外隨著體系溫度升高，CO2分子克服界面能壘由液相轉(zhuǎn)向氣相，使得液滴尺寸變小，液滴的接觸角稍有減小。

關(guān)鍵詞：粗糙壁面； 潤濕性； CO2液滴； 接觸角； 分子動(dòng)力學(xué)

中圖分類號(hào)：TK 121??? 文獻(xiàn)標(biāo)志碼：A

引用格式：宋粉紅，姜瑜輝，王志遠(yuǎn)，等.CO2液滴在粗糙固體壁面上的潤濕特性［J］. 中國石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2024，48（2）：161-169.

SONG Fenhong， JIANG Yuhui， WANG Zhiyuan， et al. Wetting behaviors of CO2 droplets on? rough solid substrate［J］. Journal of China University of Petroleum （Edition of Natural Science），2024，48（2）：161-169.

Wetting behaviors of CO2 droplets on? rough solid substrate

SONG Fenhong， JIANG Yuhui， WANG Zhiyuan， WANG Zhongxu， FAN Jing

（School of Energy and Power Engineering， Northeast Electric Power University， Jilin 132012， China）

Abstract： Understanding the role of rough solid substrates on the CO2 droplet's wetting characteristics is of importance for CO2 capture by dropwise condensation. The molecular dynamic method was adopted to investigate the effect of microstructures， potential energy and temperature on the wetting behaviors of CO2 droplets. The results show that the contact angle of CO2 droplets is enhanced on the microstructured substrate in comparison with the smooth substrate. The wetting state of the CO2 droplets is affected by the potential energy parameters on the micro-structured substrate. On the solid surface with fence structures， the CO2 droplets change from the Cassie-Baxter state to the Wenzel state， while on the surface with the triangular and ramp structures CO2 droplets are always at the Wenzel state. With the increase of the potential energy parameters， all the corresponding contact angles of CO2 droplets on the three solid surfaces decrease. Moreover， as the temperature increases， the CO2 molecules overcome the interface energy barrier from the liquid phase to the gas phase， resulting in the decrease of the droplet size. The contact angle of droplets decreases very slightly with the increase of temperature.

Keywords： roughness wall; wettability; CO2 droplet; contact angle; molecular dynamics

CO2是導(dǎo)致全球氣候變暖的主要溫室氣體之一，影響人類生產(chǎn)生活，面對這些問題各國高度重視并采取舉措［1-2］。除了減少化石能源消耗進(jìn)行低碳能源轉(zhuǎn)型外［3-5］，碳捕集與封存技術(shù)是也是減緩溫室氣體排放的有效途徑［6-7］。捕集后的CO2可二次利用，例如用于制作環(huán)保、節(jié)能材料［8］，也可以用做制冷系統(tǒng)中安全的制冷劑［9］，還可以用于提高石油、天然氣采收率等［10-13］，給各國帶來經(jīng)濟(jì)、環(huán)境雙重效益。碳捕集技術(shù)是碳捕集與封存技術(shù)過程中關(guān)鍵的步驟，在過去幾十年里已經(jīng)提出了多種CO2捕集技術(shù)，包括物理/化學(xué)吸收法、吸附法、膜分離法、低溫蒸餾法等［14-18］，這些碳捕集技術(shù)已有實(shí)際應(yīng)用，但每項(xiàng)技術(shù)仍存在不足。例如溶劑吸收法中溶劑對設(shè)備有強(qiáng)腐蝕性，溶劑的排放污染環(huán)境［19-20］；吸附法面臨吸附劑發(fā)生降解、吸附率低、對CO2的選擇性差等問題［21］；膜分離法存在膜的選擇性低、滲透性差、污染環(huán)境等問題［22-23］。所以開發(fā)高效率、低成本的碳捕集技術(shù)成為研究重點(diǎn)。在潮濕空氣或霧中通過疏水或納米結(jié)構(gòu)表面［24-27］捕集水液滴的技術(shù)已得到應(yīng)用，學(xué)者從中受到啟發(fā)并提出CO2冷凝液滴捕集技術(shù)。在純CO2蒸氣環(huán)境中，通過降低固體壁面溫度使得CO2液滴停留在其壁面上，通過施加外力或傾斜固體基底等條件，液滴流動(dòng)到收集系統(tǒng)，捕集的CO2以液態(tài)形式存在，可直接運(yùn)輸。該捕集技術(shù)的特點(diǎn)是CO2直接以液態(tài)形式出現(xiàn)并且易流動(dòng)捕集，這對液滴的潤濕狀態(tài)要求較高，而通過調(diào)控壁面潤濕特性可以控制液滴的潤濕狀態(tài)。目前，關(guān)于CO2液滴在固體壁面上的潤濕特性研究有限。Snustad等［28］綜述了壁面疏液性和滴狀冷凝的進(jìn)展，對于CO2等低表面張力流體在固體壁面上的潤濕狀態(tài)可通過壁面結(jié)構(gòu)形貌、尺寸參數(shù)進(jìn)行調(diào)節(jié)。Wu等［29-30］研究了CO2液滴在光滑和柱狀結(jié)構(gòu)壁面上潤濕現(xiàn)象，勢能參數(shù)與壁面潤濕性息息相關(guān)。在光滑表面上，發(fā)現(xiàn)CO2液滴接觸角隨著勢能參數(shù)的增大而減小。在柱狀表面上CO2液滴的接觸角顯著增加，當(dāng)增加柱間距時(shí)，發(fā)現(xiàn)CO2液滴從Cassie-Baxter狀態(tài)過渡到Wenzel狀態(tài)，由此可知粗糙結(jié)構(gòu)可以提高接觸角，改變潤濕狀態(tài)，但是研究多集中于CO2液滴在單一粗糙結(jié)構(gòu)壁面，缺乏液滴在不同粗糙結(jié)構(gòu)上潤濕狀態(tài)的對比性。筆者建立CO2液滴在3種微結(jié)構(gòu)粗糙壁面的潤濕模型，討論勢能參數(shù)、壁面粗糙形貌及溫度對CO2液滴潤濕特性的影響。

1? 分子模型及細(xì)節(jié)

采用MD方法從分子角度分析探討CO2液滴在微結(jié)構(gòu)粗糙壁面的潤濕行為。圖1為模擬系統(tǒng)固體壁面微結(jié)構(gòu)及CO2液滴初始結(jié)構(gòu)示意圖，固體基底由銅原子按面心立方晶格（FCC）結(jié)構(gòu)排列組成，晶胞尺寸為0.361 nm，基底大小為20 nm×20 nm×1.5 nm，保證液滴可以完全在壁面上鋪展。為討論壁面微結(jié)構(gòu)對CO2液滴潤濕性的影響，建立結(jié)構(gòu)形狀為欄柵形、斜坡形和三角形3種微結(jié)構(gòu)粗糙壁面，各微結(jié)構(gòu)壁面正視圖如圖1所示。3種類型微結(jié)構(gòu)壁面包含的銅原子數(shù)分別為47980、46953、53280，粗糙度因子r均約為2.4，其中粗糙度因子r定義為固體基底真實(shí)表面積與表觀表面積之比。7187個(gè)CO2分子隨機(jī)分布組成立方體液滴置于基底中心位置，立方體尺寸為10 nm×10 nm×10 nm，換算直徑約為12 nm，液滴直徑大于10 nm可減小尺寸效應(yīng)的影響［31］。立方體CO2液滴下邊緣與固體基底銅表面間的間隙為1 nm。固體基底壁面銅原子間相互作用采用經(jīng)典的EAM勢函數(shù)［32］計(jì)算。

CO2分子采用單點(diǎn)粗?；Ｐ?，該模型已經(jīng)被證實(shí)能準(zhǔn)確預(yù)測CO2分子的氣液界面行為［33-34］，優(yōu)勢在于保證計(jì)算精度的情況下加快計(jì)算速度。CO2分子之間的相互作用采用CG SAFT-γ Mie勢函數(shù)［35］計(jì)算，表達(dá)式為

φMie（r）=εC（λa，λr）σrλr-σrλa.? （1）

其中

C（λa，λr）=λrλr-λaλrλaβ，

β=λaλr-λa.

式中，ε為勢能參數(shù)；σ為勢阱參數(shù);λa為分子的吸引指數(shù)；λr為分子的排斥指數(shù)。

由于CO2粗?；蟮闹睆脚cCu晶格常數(shù)相近，所以CO2分子與銅固體基底的相互作用也采用CG SAFT-γ Mie勢函數(shù)計(jì)算。表1為 CO2粗顆粒的勢函數(shù)參數(shù)及其與Cu原子的勢函數(shù)參數(shù)。

所有MD模擬均使用LAMMPS［36］軟件計(jì)算。選取NVT系綜，采用Nose-Hoover［37］恒溫?zé)嵩》刂企w系溫度，X、Y方向采用周期性邊界，Z方向選用非周期性邊界，上邊界采用鏡像邊界以防止氣相分子跑出模擬區(qū)域。采用Velocity-Verlet算法對牛頓運(yùn)動(dòng)方程求解，時(shí)間步長為10 fs，截?cái)喟霃綖?.5 nm。首先模擬體系內(nèi)CO2液滴在溫度T為223 K下弛豫5 ns達(dá)到氣液平衡狀態(tài)，然后與壁面作用并控制體系溫度為223 K，繼續(xù)模擬10 ns使得CO2液滴在各微結(jié)構(gòu)粗糙壁面上潤濕達(dá)到新的平衡狀態(tài)，每1000個(gè)時(shí)間步長采集軌跡樣本，統(tǒng)計(jì)計(jì)算CO2液滴在各微結(jié)構(gòu)粗糙壁面上的平衡接觸角、凹槽內(nèi)分子數(shù)、氣相分子數(shù)等參數(shù)分析討論CO2液滴在粗糙表面上的潤濕特性。最后通過改變勢能參數(shù)或體系溫度討論其對CO2液滴潤濕性的影響。

在系統(tǒng)平衡后，采集1 ns周期內(nèi)樣本統(tǒng)計(jì)密度云圖。首先對統(tǒng)計(jì)樣本中液滴的質(zhì)心位置進(jìn)行統(tǒng)一，以質(zhì)心為中心，在Y方向上截取厚度2 nm作為統(tǒng)計(jì)區(qū)域。沿X和Z方向?qū)⒔y(tǒng)計(jì)區(qū)域劃分為100 ×1×80個(gè)小長方體，每個(gè)長方體的體積為，統(tǒng)計(jì)小長方體內(nèi)的粒子總數(shù)。依據(jù)公式ρ=NM/（VNA）計(jì)算出每個(gè)小立方體中的密度，N為粒子總數(shù)，M為粒子物質(zhì)的量質(zhì)量，V為小立方體體積，NA= 6.02×1023。氣-液界面密度定義為0.25 g/cm3，刪除密度小于0.25 g/cm3的區(qū)域進(jìn)而得到含液滴輪廓的密度云圖。圖2給出了液滴潤濕狀態(tài)下的密度云圖及接觸角測量方法。

根據(jù)密度云圖統(tǒng)計(jì)液滴的接觸角時(shí)，在氣液界面上畫出三相接觸線（圖中為紅色的線），每個(gè)液滴的接觸線夾角測量五次后取平均值，得到最終接觸角。

2? 結(jié)果與分析

2.1? CO2液滴在各微結(jié)構(gòu)壁面的潤濕性

液滴在微結(jié)構(gòu)壁面的潤濕模型主要分為Cassie-Baxter狀態(tài)和Wenzel狀態(tài)。Cassie-Baxter模型認(rèn)為液滴不能完全潤濕凹槽，液滴與壁面之間存在空氣，固液實(shí)際接觸面積小于表觀面積。Wenzel模型認(rèn)為液滴完全充滿凹槽并潤濕壁面，固液實(shí)際接觸面積大于表觀面積。而液態(tài)CO2為低表面張力流體，溫度為-52.2 ℃時(shí)對應(yīng)16.54 mN/m，其黏性力較強(qiáng)，不利于液滴在微結(jié)構(gòu)壁面上保持Cassie-Baxter潤濕狀態(tài)，影響液滴的捕集。壁面結(jié)構(gòu)形貌是影響壁面潤濕性的重要因素，直接影響液滴在壁面上的潤濕狀態(tài)。以CO2與壁面銅原子勢能參數(shù)εCO2-Cu =0.01 eV為例，分析CO2液滴在3種微結(jié)構(gòu)粗糙壁面的潤濕性。

圖3為CO2液滴在各微結(jié)構(gòu)粗糙壁面上潤濕達(dá)到平衡狀態(tài)時(shí)某時(shí)刻快照。同時(shí)CO2液滴在各壁面潤濕的平衡接觸角分布如圖4所示。由圖3可以看出，CO2液滴在欄柵結(jié)構(gòu)壁面上處于Cassie-Baxter狀態(tài)，凹槽底部存在小間隙，液滴聚積在凹槽頂部，接觸角稍大；而在斜坡和三角結(jié)構(gòu)壁面上潤濕時(shí)CO2液滴處于Wenzel狀態(tài)，CO2分子鋪滿凹槽底部，此時(shí)粗糙壁面近似平板壁面，壁面微結(jié)構(gòu)對液滴潤濕的平衡狀態(tài)影響減弱，接觸角稍小。從圖4也可以看出，相比于光滑平板壁面［18］，CO2液滴在各微結(jié)構(gòu)壁面潤濕的接觸角均增加，且變化規(guī)律為θ欄柵>θ斜坡>θ三角>θ平板，即CO2液滴在欄柵結(jié)構(gòu)壁面的接觸角最大，其次為斜坡和三角結(jié)構(gòu)表面，均大于平板壁面上的接觸角。為確保得到的接觸角具有可靠性，分別測量5次后在圖4中給出誤差范圍，其中最大標(biāo)準(zhǔn)差為0.83，說明測量誤差范圍很小，數(shù)據(jù)可靠。

為進(jìn)一步分析CO2液滴在各微結(jié)構(gòu)粗糙壁面呈現(xiàn)不同潤濕狀態(tài)的原因，圖5為各微結(jié)構(gòu)粗糙壁面凹槽內(nèi)CO2分子密度隨時(shí)間變化，凹槽內(nèi)分子密度曲線一定程度上可以反映凹槽結(jié)構(gòu)對潤濕影響的差異性。從圖5中可以看出，當(dāng)液滴潤濕達(dá)到平穩(wěn)狀態(tài)時(shí)，各微結(jié)構(gòu)凹槽內(nèi)分子密度數(shù)值變化。欄柵結(jié)構(gòu)凹槽內(nèi)分子密度最小，三角結(jié)構(gòu)的最大，斜坡結(jié)構(gòu)凹槽內(nèi)分子密度介于二者之間。斜坡和三角微結(jié)構(gòu)凹槽的開口由大漸小，溝壑淺，CO2分子容易進(jìn)入凹槽，已進(jìn)入的分子受壁面強(qiáng)吸附作用黏附在凹槽內(nèi)部，當(dāng)微結(jié)構(gòu)凹槽被大量分子覆蓋后，對應(yīng)的粗糙壁面逐步趨近于平板壁面，導(dǎo)致CO2液滴的平衡接觸角稍小。與前兩類結(jié)構(gòu)壁面相反，欄柵微結(jié)構(gòu)凹槽間距密集，溝壑深，在一定程度上阻礙CO2分子滲入凹槽區(qū)域，部分CO2分子進(jìn)入凹槽對粗糙壁面的潤濕性影響小，CO2液滴接觸角較大。綜上所述，凹槽開口大，導(dǎo)致凹槽內(nèi)的分子密度大，此類凹槽結(jié)構(gòu)不易阻礙分子的滲入，液滴的接觸角稍??；而凹槽的溝壑深，間隙密集，導(dǎo)致凹槽內(nèi)的分子密度越少，此類凹槽結(jié)構(gòu)越能阻礙分子的滲入，相應(yīng)液滴接觸角稍大，此類型的微結(jié)構(gòu)壁面在CO2冷凝捕集技術(shù)上占有較大的優(yōu)勢。

2.2? 勢能參數(shù)對CO2液滴潤濕影響

CO2液滴在不同類型潤濕壁面上呈現(xiàn)出各異的潤濕狀態(tài)。若逐一建立不同類型壁面進(jìn)行模擬計(jì)算會(huì)消耗大量資源，通過調(diào)整勢能參數(shù)可以近似代替不同類型的潤濕壁面，進(jìn)而研究CO2液滴在各微結(jié)構(gòu)壁面的潤濕性。討論不同CO2與壁面間的勢能參數(shù)εCO2-Cu下CO2液滴在各微結(jié)構(gòu)壁面的潤濕規(guī)律。CO2液滴在3種微結(jié)構(gòu)壁面上潤濕達(dá)到平衡時(shí)對應(yīng)的接觸角隨勢能參數(shù)變化如圖6所示。很明顯，接觸角隨著勢能參數(shù)的增大而減小。當(dāng)勢能參數(shù)從0.007 eV增加到0.012 eV時(shí)，CO2液滴在欄柵形、三角形及斜坡形結(jié)構(gòu)壁面上的潤濕接觸角分別減小37.5°、61°和55°。

圖7為當(dāng)勢能參數(shù)變化時(shí)CO2液滴在各微結(jié)構(gòu)壁面上達(dá)到平衡狀態(tài)時(shí)的平衡狀態(tài)圖。隨著勢能參數(shù)增大，液滴在欄柵壁面上的狀態(tài)發(fā)生改變，從Cassie-Baxter狀態(tài)轉(zhuǎn)向Wenzel狀態(tài)（圖7（a））；而在三角和斜坡壁面上CO2液滴的狀態(tài)均處于Wenzel狀態(tài)（圖7（b）、（c））。當(dāng)CO2液滴處于Cassie-Baxter狀態(tài)時(shí)，分子在滲入凹槽之前聚集在凹槽附近區(qū)域上，同時(shí)在凹槽底部出現(xiàn)小間隙。此時(shí)液滴與表面的實(shí)際接觸面積小于真實(shí)表面積，壁面與液滴間相互作用力減弱，壁面的潤濕能力降低，接觸角較大。而當(dāng)CO2液滴處于Wenzel狀態(tài)時(shí)，液滴與粗糙壁面的實(shí)際接觸面積增加，粗糙壁面對液滴的作用力增強(qiáng)，使得接觸角減小。

為了進(jìn)一步揭示勢能參數(shù)對CO2液滴在微結(jié)構(gòu)粗糙壁面上潤濕的影響，圖8給出了在不同勢能參數(shù)ε下各微結(jié)構(gòu)壁面凹槽內(nèi)CO2分子數(shù)隨時(shí)間變化。由圖8可以看出，各微結(jié)構(gòu)壁面凹槽內(nèi)CO2分子數(shù)均隨著勢能參數(shù)增大而增加。這是因?yàn)殡S著勢能參數(shù)增加，即壁面與CO2分子間相互作用力增大，近壁面處的氣相分子先吸附在粗糙壁面，同時(shí)氣相與液相分子存在相互作用，驅(qū)使CO2液滴向表面移動(dòng)，當(dāng)液滴與壁面間距離小于截?cái)喟霃胶螅合喾肿优c壁面發(fā)生相互作用，最后液滴潤濕壁面。在此過程中大量CO2分子聚集壁面及凹槽內(nèi)，凹槽內(nèi)分子數(shù)不斷增加，當(dāng)CO2液滴處于平衡狀態(tài)時(shí)凹槽內(nèi)分子數(shù)保持相對穩(wěn)定數(shù)量。

當(dāng)勢能參數(shù)εCO2-Cu ≥0.01 eV時(shí)，模擬系統(tǒng)出現(xiàn)某瞬間凹槽內(nèi)CO2分子數(shù)驟增現(xiàn)象。這一現(xiàn)象可通過CO2液滴的動(dòng)態(tài)潤濕過程解釋。圖9給出了在εCO2-Cu =0.01 eV情況下，液滴在三角形微結(jié)構(gòu)壁面上的動(dòng)態(tài)潤濕過程。由圖9可知：當(dāng)t=1 ns時(shí)，CO2液滴懸浮在壁面上方，隨著模擬進(jìn)行，分子間的相互作用使得CO2液滴緩慢向壁面運(yùn)動(dòng)，質(zhì)心降低；當(dāng)t=3.5 ns時(shí)，CO2液滴底部到粗糙壁面距離小于截?cái)喟霃綍r(shí)，壁面對液滴的作用力增強(qiáng)，液滴下移加快；當(dāng)t=4.3 ns時(shí)，CO2液滴底端與粗糙壁面接觸，液滴質(zhì)心位置繼續(xù)下降，液滴在三角形微結(jié)構(gòu)壁面的潤濕面積增加；當(dāng)t=4.6 ns時(shí)，CO2液滴在粗糙壁面的潤濕面積約達(dá)38 nm2，凹槽內(nèi)CO2分子也隨之增加，因此在4.3～4.6 ns時(shí)間間隔內(nèi)凹槽內(nèi)CO2分子數(shù)呈現(xiàn)突然驟增的現(xiàn)象。當(dāng)CO2液滴在粗糙壁面潤濕達(dá)到動(dòng)態(tài)平衡時(shí)，粗糙壁面凹槽內(nèi)CO2分子數(shù)在某一數(shù)值附近上下波動(dòng)。

2.3? 溫度對液滴潤濕影響

體系溫度決定分子的微觀熱運(yùn)動(dòng)狀態(tài)，表征分子熱運(yùn)動(dòng)的劇烈程度，進(jìn)而影響液滴的形貌及潤濕狀態(tài)。以欄柵形微結(jié)構(gòu)壁面與CO2的勢能參數(shù)為εCO2-Cu =0.008和0.012 eV為例，即疏液壁面和親液壁面。選定T為223、233、243 和253 K四個(gè)體系溫度，探究溫度對CO2液滴在微結(jié)構(gòu)壁面潤濕行為的微觀機(jī)制。

在CO2液滴潤濕動(dòng)態(tài)過程中，體系內(nèi)CO2分子處于氣液共存狀態(tài)，溫度是影響氣液分子數(shù)變化的重要因素。統(tǒng)計(jì)氣液分子數(shù)可以揭示系統(tǒng)內(nèi)氣液動(dòng)態(tài)變化過程。根據(jù)模擬條件下液態(tài)CO2密度（ρ=1054 kg/m3），計(jì)算得出CO2分子的平均間距約為0.41 nm。在分子模擬中，當(dāng)CO2分子間的距離小于0.41 nm時(shí)，則認(rèn)定兩個(gè)CO2分子聚集一起形成團(tuán)簇。在采集樣本中統(tǒng)計(jì)壁面凹槽頂部以上空間內(nèi)形成的團(tuán)簇個(gè)數(shù)及團(tuán)簇內(nèi)包含的分子數(shù)，當(dāng)團(tuán)簇內(nèi)包含的分子數(shù)小于等于3 時(shí)，認(rèn)為是氣相分子并進(jìn)行累積得最終氣相分子數(shù)，反之累積為液相分子數(shù)。統(tǒng)計(jì)模擬中10 ns系統(tǒng)內(nèi)氣相和液相分子數(shù)，將每隔0.5 ns對應(yīng)的采集樣本計(jì)算出氣相分子和液相分子數(shù)并繪制隨時(shí)間變化曲線，見圖10。

由圖10可知，整體趨勢走向一致，隨著溫度升高，氣相分子數(shù)逐漸增加，液相分子數(shù)逐漸減少，但是趨勢變化幅度有差別。氣相分子數(shù)的起點(diǎn)不是零點(diǎn)，原因是模擬初始階段在溫度T為223 K下進(jìn)行弛豫，達(dá)到該溫度下的氣液平衡狀態(tài)，在此基礎(chǔ)上改變體系溫度，進(jìn)行其他溫度工況的模擬，因此圖中T為223 K的氣相和液相分子數(shù)變化曲線近似為直線。

氣相分子數(shù)和液相分子數(shù)為對立關(guān)系，以氣相分子數(shù)為例分析，對應(yīng)圖10（a）、（c）。對于疏液壁面，溫度階梯變化，氣相分子數(shù)對應(yīng)變化范圍大；對于親液壁面，氣相分子數(shù)對應(yīng)變化范圍小。溫度越高時(shí)，兩種類型壁面上的氣相分子數(shù)終值差距越明顯。

由于體系溫度升高，CO2分子熱運(yùn)動(dòng)劇烈，動(dòng)能增加，液相分子間的束縛力減弱，CO2液滴外層的CO2分子開始向氣相空間運(yùn)動(dòng)，使得氣相分子數(shù)不斷增加。在疏液壁面上，CO2液滴與壁面的作用力較小，當(dāng)受到溫度擾動(dòng)后，固液間作用力不足以將液相分子吸附在壁面上，大量液相分子轉(zhuǎn)為氣相分子。而液滴與親液壁面間固液作用力大，可有效約束CO2液滴外層的分子。

圖11為在不同溫度下CO2液滴在各類型壁面上潤濕達(dá)到平衡狀態(tài)時(shí)的平衡狀態(tài)圖，其對應(yīng)的平衡接觸角隨溫度的變化如圖12所示。可以看出，不同類型壁面上CO2液滴的平衡接觸角隨溫度升高而減小，但不論CO2液滴是在親液壁面還是疏液壁面上潤濕的接觸角降幅均較小。這是因?yàn)?，隨著溫度升高，分子熱運(yùn)動(dòng)加劇，而液態(tài)CO2分子間作用力較小，CO2分子極易克服界面能壘從液相變成氣相，僅少量液相分子在壁面吸引作用下在壁面鋪展，因而隨溫度升高CO2液滴潤濕形態(tài)變化不明顯。當(dāng)溫度從223 K升高至253 K時(shí)，疏液壁面上CO2液滴潤濕接觸角減小約5°，親液壁面上CO2液滴潤濕接觸角減小約7°。同時(shí)隨著溫度升高二氧化碳分子從液相轉(zhuǎn)向氣相使得CO2液滴尺寸明顯變小。

溫度對親疏液壁面上的二氧化碳?xì)庖悍肿訑?shù)變化規(guī)律影響明顯，隨著溫度升高，氣相分子數(shù)增加，液相分子數(shù)減小。此外在親疏液壁面上CO2液滴的接觸角隨溫度升高呈小幅度減小。

3? 結(jié)? 論

（1）微結(jié)構(gòu)壁面提高了CO2液滴的平衡接觸角，當(dāng)液滴處于Cassie-Baxter狀態(tài)，液滴與壁面的實(shí)際接觸面積較小，接觸角稍大；而液滴處于Wenzel狀態(tài)，液滴與壁面的實(shí)際接觸面積較大，接觸角稍小，通過對壁面微結(jié)構(gòu)的設(shè)計(jì)可以促進(jìn)壁面趨向疏液性。

（2）勢能參數(shù)可以控制CO2液滴在微結(jié)構(gòu)壁面上的潤濕狀態(tài)，隨著勢能參數(shù)增大，在欄柵結(jié)構(gòu)壁面上CO2液滴從Cassie-Baxter狀態(tài)轉(zhuǎn)向Wenzel狀態(tài)，而在三角和斜坡結(jié)構(gòu)壁面上CO2液滴均保持Wenzel狀態(tài)；欄柵壁面上CO2液滴處于Cassie-Baxter狀態(tài)，液滴與壁面的接觸面積小，接觸角大，有利于液滴的捕集。

（3）溫度對CO2液滴潤濕狀態(tài)影響較小，但明顯影響系統(tǒng)內(nèi)的氣液分子數(shù)變化；在親疏液兩種類型壁面上液滴的接觸角隨溫度升高稍有減小，CO2分子克服界面能壘從液相變成氣相，導(dǎo)致氣液分子數(shù)變化明顯，而僅少量液相分子在吸引力作用下在壁面上鋪展，液滴的接觸角變化幅度較小。

參考文獻(xiàn)：

［1］? ARESTA M， DIBENEDETTO A， ANGELINI A. Catalysis for the valorization of exhaust carbon： from CO2 to chemicals， materials， and fuels： technological use of CO2［J］. Chemical Reviews， 2014，114：1709-1742.

［2］? 呂建中.保障能源安全需做好傳統(tǒng)能源與新能源戰(zhàn)略協(xié)同［J］.世界石油工業(yè)，2022，29（2）：1-6.

L Jianzhong. Necessary to coordinate traditional and new energy strategies to ensure energy security［J］. World Petroleum Industry， 2022，29（2）：1-6.

［3］? 崔楊，安寧，付小標(biāo)，等.面向高比例新能源電力系統(tǒng)調(diào)峰需求的儲(chǔ)能容量配置方法綜述［J］.東北電力大學(xué)學(xué)報(bào)，2023，43（1）：1-8.

CUI Yang， AN Ning， FU Xiaobiao， et al. Overview of energy storage capacity allocation methods for high-proportion new energy power system peak shaving demand［J］. Journal of Northeast Electric Power University， 2023，43（1）：1-8.

［4］? 蘇娟，李沂，盧穎，等.大變局下的能源安全與低碳轉(zhuǎn)型：第五屆中國石油國際合作論壇綜述［J］.世界石油工業(yè)，2022，29（6）：13-18.

SU Juan， LI Yi， LU Ying， et al. Energy security and low-carbon transition in the changing world： overview of the 5th CNPC international cooperation forum［J］. World Petroleum Industry， 2022，29（6）：13-18.

［5］? 王宜林.踔厲邁向能源安全與能源轉(zhuǎn)型并舉的新征程［J］.世界石油工業(yè)，2023，30（2）：1-4.

WANG Yilin. Strive for a new journey towards energy security and energy transition［J］. World Petroleum Industry， 2023，30（2）：1-4.

［6］? 丁根.國際油服公司碳捕集、利用與封存業(yè)務(wù)的現(xiàn)狀與展望［J］.世界石油工業(yè)，2022，29（4）：28-33.

DING Gen. Current situation and prospects of carbon capture，utilization and sequestration business of international oilfield services companies［J］. World Petroleum Industry， 2022，29（4）：28-33.

［7］? BUDINIS S， KREVOR S， DOWELL N M， et al. An assessment of CCS costs， barriers and potential［J］. Energy Strategy Reviews， 2018，22：61-81．

［8］? SONG C.CO2 conversion and utilization： an overview［J］. ACS Symposium Series， 2002，809：1-30.

［9］? HUNT A J， SIN E H K， MARRIOTT R， et al. Generation， capture， and utilization of industrial carbon dioxide［J］. Chemistry Sustainability Energy Materials， 2010，3（3）：306-322.

［10］? 劉樹陽，孫寶江，宋永臣，等.CO2提高天然氣采收率的重力效應(yīng)與氣藏壓力影響模擬［J］.中國石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2020，44（3）：81-89.

LIU Shuyang， SUN Baojiang， SONG Yongchen， et al. Simulation on gravity effect and reservoir pressure influence analysis in CO2 enhanced gas recovery［J］．Journal of China University of Petroleum（Edition of Natural Science）， 2020，44（3）：81-89.

［11］? LI X， ELSWORTH D. Geomechanics of CO2 enhanced shale gas recovery ［J］. Journal of Natural Gas Science and Engineering， 2015，26：1607-1619.

［12］? 任韶然，杜翔睿，孫志雄，等.油田驅(qū)油用CO2泡沫劑體系研究進(jìn)展［J］.中國石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，46（4）：102-108.

REN Shaoran， DU Xiangrui， SUN Zhixiong， et al. State of the art in research on CO2 foam agents for oilfield applications ［J］.Journal of China University of Petroleum （Edition of Natural Science）， 2022，46（4）：102-108.

［13］? NICHOLAS A A， NAKLES D V， GORECKI C D， et al. CO2 storage associated with CO2 enhanced oil recovery： a statistical analysis of historical operations ［J］. International Journal of Greenhouse Gas Control， 2015（37）：384-397.

［14］? DERKS P W J， KLEINGELD T， AKEN C V. Kinetics of absorption of carbon dioxide in aqueous piperazine solutions［J］. Chemical Engineering Science， 2006，61：6837-6854.

［15］? BEN-MANSOUR R， HABIB M A， BAMIDELE O E， et al. Carbon capture by physical adsorption： materials， experimental investigations and numerical modeling and simulations-a review［J］. Applied Energy， 2016，161：225-255.

［16］? SREEDHAR I， VAIDHISWARAN R， KAMANI B M， et al. Process and engineering trends in membrane based carbon capture［J］. Renewable and Sustainable Energy Reviews， 2017，68：659-684.

［17］? GOTTLICHER G， PRUSCHEK R. Comparison of CO2 removal systems for fossil fueled power plants［J］. Energy Convers Manag， 1997，38：173-178.

［18］? YAVE W， CAR A， FUNARI S S， et al. CO2-philic polymer membrane with extremely high separation performance ［J］. Macromolecules， 2010，43（1）：326-333.

［19］? ROCHELLE G T. Thermal degradation of amines for CO2 capture ［J］. Current Opinion Chemical Engineering， 2012，1（2）：183-190.

［20］? ZOANNOU K S， SAPSFORD D J，GRIFFITHS A J. Thermal degradation of monoethanolamine and its effect on CO2 capture capacity ［J］. International Journal of Greenhouse Gas Control， 2013，17：423-430.

［21］? FAROOQ M， SAEED M A， IMRAN M， et al.CO2 capture through electro-conductive adsorbent using physical adsorption system for sustainable development［J］. Environment Geochemistry and Health， 2020，42：1507-1515.

［22］? BRUNETTI A， SCURA F， BARBIERI G， et al. Membrane technologies for CO2 separation ［J］. Journal Membrane Science， 2010，359（1）：15-25.

［23］? BERNARDO P， DRIOLI E， GOLEMME G. Membrane gas separation： a review/state of the art ［J］. Industrial & Engineering Chemistry Research， 2009，48：4638-4663.

［24］? LEI C， WANG S Y， XIANG X， et al. Mechanism of surface nanostructure changing wettability： a molecular dynamics simulation ［J］. Computational Materials Science， 2020，171：109223.

［25］? HAHMARDI A， TAMMISOLA O， CHINAPPI M， et al. Effects of surface nanostructure and wettability on pool boiling： a molecular dynamics study ［J］. International Journal of Thermal Sciences， 2021，167：106980.

［26］? 全學(xué)波，董佳奇，周健.疏水表面拓?fù)浣Y(jié)構(gòu)對其潤濕狀態(tài)影響的粗?；M［J］.化學(xué)學(xué)報(bào)，2014，72（10）：1075-1078.

QUAN Xuebo， DONG Jiaqi， ZHOU Jian. Effect of topology of hydrophobic surfaces on their wetting states by coarse-grained simulations［J］. Journal of Chemical，2014，72（10）：1075-1078.

［27］? 朱毅，吳曉敏，褚福強(qiáng).微觀結(jié)構(gòu)尺度對液滴潤濕狀態(tài)的影響［J］.工程熱物理學(xué)報(bào)，2019，40（11）：2632-2636.

ZHU Yi， WU Xiaomin， ZHU Fuqiang. Influence on droplet wetting state by microstructural scale［J］. Journal of Engineering Thermophysics， 2019，40（11）：2632-2636.

［28］? SNUSTAD I， ROE I T， BRUNSVOLD A， et al. A review on wetting and water condensation-perspectives for CO2 condensation［J］. Advance in Colloid Interface Science， 2018，256：291-304.

［29］? WU J， ERVIK A， SNUSTAD I， et al. Contact angle and condensation of a CO2 droplet on a solid surface［J］. The Journal of Physical Chemistry C， 2018，123（1）：443-451.

［30］? WU J， SNUSTAD I， ERVIK A， et al. CO2 wetting on pillar-nanostructured substrates［J］. Nanotechnology， 2020，31（24）：245403.

［31］? 曹炳陽，陳民，過增元.納米結(jié)構(gòu)表面浸潤性質(zhì)的分子動(dòng)力學(xué)研究［J］.高等學(xué)校化學(xué)學(xué)報(bào)，2005（2）：277-280.

CAO Bingyang， CHEN Min， GUO Zengyuan. Wettability of surface with nano-structures studied by using molecular dynamics simulation［J］. Chemical Journal of Chinese Universities. 2005（2）：277-280.

［32］? MISHIN Y， MEHL M J， PAPACONSTANTOPOULOS D A， et al. Structural stability and lattice defects in copper： Ab initio， tight-binding， and embedded-atom calculations［J］. Physical Review B， 2001，63：224106.

［33］? MLLER E A， JACKSON G. Force-field parameters from the SAFT-γ equation of state for use in coarse-grained molecular simulations［J］. Annual Review of Chemical and Biomolecular Engineering， 2014，5：405-427.

［34］? HOANG H， DELAGE-SANTACREU S， GALLIERO G. Simultaneous description of equilibrium， interfacial， and transport properties of fluids using a Mie chain coarse-grained force field［J］. Industrial & Engineering Chemistry Research， 2017，56：9213-9226.

［35］? AVENDANO C， LAFITTE T， GALINDO A， et al. SAFT-γ force field for the simulation of molecular fluids 1： a single-site coarse grained model of carbon dioxide［J］. The Journal of Physical Chemistry B， 2011，115（38）：11154-11169.

［36］? PLIMPTON S. Fast parallel algorithms for short-range moleculajr dynamics［J］. Journal of Computational Physics，1995，117（1）：1-19.

［37］? HNENBERGER P H. Thermostat algorithms for molecular mynamics simulations［M］.Advanced Computer Simulation， 2005：105-149.

（編輯? 沈玉英）