趙彥琳 金陶勝

摘要 [目的]研究甲苯在氣液界面表面微層中的遷移變化規(guī)律。[方法]使用計(jì)算流體力學(xué)（CFD）軟件對(duì)表面微層中甲苯氣體的遷移變化進(jìn)行數(shù)值模擬，通過(guò)對(duì)液相流速、氣相流速、氣相濃度的模擬比較，討論不同條件下的甲苯氣體遷移變化規(guī)律。[結(jié)果]液相流速的增大，影響表面微層的形成，從而減小了氣體的遷移通量；氣相流速的增大，提高了氣體的遷移通量，同時(shí)縮短了達(dá)到傳質(zhì)平衡的時(shí)間；在吸收開(kāi)始的最初階段，氣相濃度越高，表面微層中氣體的吸收也越快，隨著時(shí)間的變化，對(duì)氣體的吸收效率逐漸下降。[結(jié)論]研究所獲得的成果將有助于認(rèn)識(shí)有機(jī)污染物在有機(jī)氣-液界面的遷移機(jī)理。

關(guān)鍵詞 甲苯；氣液界面；表面微層；遷移變化；數(shù)值模擬

中圖分類號(hào) S181.3 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2014）20-06646-04

Study on the Numerical Simulation of the Migration Variation of Toluene in the Surface Microlayer of the Gasliquid Interface

ZHAO Yanlin et al

（Yangling Vocational & Technical College，Yangling， Shaanxi 712100）

Abstract [Objective] The research aimed to study the migration variation of toluene in the surface microlayer of the gasliquid interface. [Method] A kind of CFD （Computational Fluid Dynamics） software was used to simulate the migration variation of toluene in the surface microlayer. Trough the simulation and comparison of the liquid flow rate， gas flow rate and concentration of gas phase，the migration variation of toluene under different conditions were discussed.[Result] Higher flow rate of liquid phase would slow down the formation of the surface microlayer and the transfer flux would decrease obviously. Higher flow rate of gas phase could shorten the time of reaching transfer balance， at the same time， the transfer flux would increase. At the beginning of the absorption， higher concentration of gas phase can accelerate the absorption of toluene in the surface microlayer. With the time changing， the concentration of toluene reached steady state. [Conclusion] The study is helpful to understand the migration mechanism of the organic pollutant in the surface microlayer of the gasliquid interface.

Key words Toluene；Gasliquid interface；Surface microlayer；Transfer variation；Numerical simulation

大部分的有機(jī)廢氣毒性大，可經(jīng)人體呼吸道、皮膚和消化道吸收，長(zhǎng)期暴露在有機(jī)空氣中會(huì)給人體造成巨大的傷害[1]。因此，對(duì)有機(jī)廢氣的控制處理研究有至關(guān)重要的意義。目前國(guó)內(nèi)常見(jiàn)的處理方法為吸收法[2]，即利用液體吸收劑對(duì)有機(jī)廢氣的選擇性吸收或反應(yīng)，進(jìn)而達(dá)到去除廢氣的最終目的。經(jīng)研究發(fā)現(xiàn)這種吸收作用主要發(fā)生在氣液界面處的表面微層中，而長(zhǎng)久以來(lái)，國(guó)內(nèi)外學(xué)者針對(duì)氣液界面表面微層的研究多為水體界面[3-5]，對(duì)于有機(jī)液面的研究鮮有報(bào)道。另外，由于表面微層的厚度通常只有幾十微米到幾百微米之間[6]，常規(guī)的實(shí)驗(yàn)室很難達(dá)到相應(yīng)的實(shí)驗(yàn)操作條件，因此，筆者擬針對(duì)典型的甲苯吸收液進(jìn)行表面微層的性質(zhì)研究，用數(shù)值模擬的方法分析在該表面微層內(nèi)甲苯的濃度變化以及影響該濃度變化的因素，以便有助于認(rèn)識(shí)有機(jī)污染物在有機(jī)氣-液界面的遷移機(jī)理。

1 模型的建立

1.1 模型基本假設(shè)

基于吸收液對(duì)吸收質(zhì)的物理吸收過(guò)程，對(duì)模型作如下假設(shè)：①氣體吸收操作過(guò)程為穩(wěn)態(tài)，且液體不可壓縮，即液體密度不隨溫度和濃度的變化而變化；②僅僅氣相中的吸收質(zhì)組分被吸收液吸收，吸收液中的成分不向氣相中傳遞；③忽略氣相與液相間熱量傳遞；④氣體吸收過(guò)程是絕熱的，即吸收容器與環(huán)境之間沒(méi)有熱量交換。

1.2 吸收傳質(zhì)模型的選擇

為了更貼近現(xiàn)實(shí)操作，該研究利用多相流模型和物質(zhì)傳輸模型以及軟件提供的用戶自定義接口UDF，對(duì)甲苯氣體在表面微層中的遷移變化進(jìn)行模擬，以達(dá)到試驗(yàn)研究的目的。

1.3 幾何模型的構(gòu)建

為了通過(guò)數(shù)值模擬考察氣液兩相吸收的傳質(zhì)過(guò)程，將吸收容器簡(jiǎn)化為如圖1所示的二維結(jié)構(gòu)，模型尺寸為長(zhǎng)20 cm、寬10 cm。其中上層為氣相主體，下層為液相主體，每相寬分別為5 cm，兩相的交界處即為氣液交界面，其中左側(cè)分別為氣體、液體速度入口，右側(cè)分別為氣體、液體壓力出口。構(gòu)建好模型后需要對(duì)模型進(jìn)行網(wǎng)格劃分，網(wǎng)格長(zhǎng)度為0.2 cm，共劃分5 000個(gè)網(wǎng)格（圖2）。同時(shí)，由于要對(duì)表面微層作著重研究，因此在氣液交界面處增加邊界層網(wǎng)格，起到網(wǎng)格加密的作用，使氣液界面處的模擬研究更加準(zhǔn)確精密。其中第一層邊界層距氣液交界面距離為10 μm，增長(zhǎng)系數(shù)Growth factor設(shè)為1.1，共畫(huà)出25層邊界層，整個(gè)邊界層厚度約為1 mm（圖3）。

1.4 邊界條件及其他參數(shù)設(shè)定

氣體和液體入口均設(shè)為

速度入口（velocityinlet）；氣體和液體出口均設(shè)為壓力出口

（pressureoutlet）；氣液交界面屬性設(shè)為wall；湍流模型選用k-ε模型，非穩(wěn)態(tài)unsteady過(guò)程；設(shè)置錄像功能，記錄不同時(shí)間內(nèi)甲苯的傳質(zhì)變化；氣相主體材料為甲苯和空氣的混合物，其中甲苯的質(zhì)量分?jǐn)?shù)為10%；液相主體材料為甲苯和吸收溶劑的混合物，假設(shè)初始狀態(tài)時(shí)，甲苯的質(zhì)量分?jǐn)?shù)為0。

2 結(jié)果與分析

設(shè)置氣相進(jìn)口速度為0.1 m/s，液相進(jìn)口速度為0.01 m/s，氣液分配系數(shù)為2，在模型模擬過(guò)程中設(shè)定觀看錄像功能Animate。模擬結(jié)果發(fā)現(xiàn)，到12 s時(shí)傳質(zhì)過(guò)程達(dá)到平衡狀態(tài)，且在傳質(zhì)過(guò)程中，液相中（深度為0～5 cm）甲苯的濃度呈分層分布（圖4）。其中，達(dá)到甲苯傳質(zhì)平衡時(shí)（12 s時(shí)），距離模型左側(cè)進(jìn)口約15～18 cm處表面微層中（小于500 μm）

的濃度分布變化如圖5所示。從圖5可以看出，表面微層中也存在明顯的分層現(xiàn)象，且越接近界面處，甲苯質(zhì)量分?jǐn)?shù)濃度越大，即在表面微層中，甲苯具有明顯的富集現(xiàn)象，因此也可以將甲苯質(zhì)量分?jǐn)?shù)濃度差看成其積累遷移的主要?jiǎng)恿?。影響表面微層中甲苯遷移的因素可能包括氣相、液相近界面處的湍動(dòng)程度以及氣液質(zhì)量分配系數(shù)等，對(duì)于湍動(dòng)程度的模擬可以通過(guò)改變液相流速和氣相流速實(shí)現(xiàn)。取表面微層中一固定縱剖面為研究對(duì)象，假設(shè)以模型左下角點(diǎn)為原點(diǎn)（0，0），則?。?.18 m，h）縱剖面研究對(duì)象，其中深度h分別取值10、50、100、150、200、250、300、350、400、450、500 μm。

圖5 表面微層中甲苯的質(zhì)量分?jǐn)?shù)等值線

2.1 液相流速的影響

分別設(shè)定液相流速為0.005、0.010、0.020、0.030 m/s（氣相流速不變，仍為0.1 m/s，氣液分配系數(shù)為2），得到不同深度甲苯的質(zhì)量分?jǐn)?shù)變化（圖6）。從圖6可以看出，隨著流速的增大，甲苯的質(zhì)量分?jǐn)?shù)逐漸減小，這是因?yàn)榱魉佥^大時(shí)，液相近表面處湍動(dòng)程度較大，影響了液膜即表面微層的形成，從而導(dǎo)致遷移通量的減??；相反，液相流速較低時(shí)，甲苯在表面微層中的遷移通量較大，可以理解為甲苯此時(shí)在表面微層中得到了充分的吸收，有利于甲苯的去除。另外，甲苯質(zhì)量分?jǐn)?shù)在表面微層中隨深度呈明顯的線性關(guān)系，且擬合結(jié)果非常好，對(duì)于液相速度為0.005、0.010、0.020、0.030 m/s時(shí)甲苯隨深度的變化率分別為-9.239 46×10-10、-9.665 52×10-10、-1.241 27×10-9、-1.426 13×10-9，可見(jiàn)流速為0.030 m/s時(shí)甲苯質(zhì)量分?jǐn)?shù)隨深度的變化較大，即可認(rèn)為其在液相中的擴(kuò)散速率較大，同理，流速為0.020、0.010 m/s時(shí)，擴(kuò)散速率逐漸減小，流速為0.005 m/s時(shí)擴(kuò)散速率最小。這是因?yàn)榱魉佥^大時(shí)，液相的湍動(dòng)程度較大，進(jìn)而可以理解為物質(zhì)擴(kuò)散的推動(dòng)力較大，從而加快了甲苯向下的的擴(kuò)散速率。此外，液相流速的增加可以看成是液相吸收液流量的增大，即單位體積的液相吸收劑單位時(shí)間內(nèi)吸收的甲苯的量減小了，增大了液相中的濃度梯度，使傳質(zhì)過(guò)程向著更有利于擴(kuò)散的方向發(fā)展。

圖6 不同深度甲苯的質(zhì)量分?jǐn)?shù)變化

綜上考慮，液相流速增大時(shí)，液相中甲苯的擴(kuò)散速率較大，但同時(shí)影響了液膜即表面微層的形成，從而減小了甲苯的遷移通量，因此，要綜合考慮這兩方面的影響因素，選擇較為合適的液相流速，進(jìn)而達(dá)到最好的吸收傳質(zhì)效果。

2.2 氣相流速的影響

分別設(shè)置氣相流速為0.05、0.10、0.20、0.30 m/s，設(shè)定液相進(jìn)口流速為0.01 m/s，氣液分配系數(shù)為2，進(jìn)行模擬運(yùn)算，結(jié)果發(fā)現(xiàn)，當(dāng)流速為0.05 m/s時(shí)，達(dá)到迭代收斂需要21 s；流速為0.10 m/s時(shí)，達(dá)到迭代收斂需要12 s；流速為0.20 m/s時(shí)，需要9 s；流速為0.30 m/s時(shí)，僅需要6 s即可達(dá)到迭代收斂?？梢?jiàn)，隨著氣相流速的增大，迭代收斂加快，即從模擬運(yùn)算角度來(lái)看，氣相流速越大，達(dá)到氣液平衡越快，即氣相流速的增大，縮短了甲苯分子傳遞至界面的時(shí)間，進(jìn)而加快了整個(gè)吸收過(guò)程達(dá)到傳質(zhì)平衡。這是因?yàn)榱鲃?dòng)的氣相通過(guò)對(duì)氣液界面直接的剪切作用以及因此產(chǎn)生的表面波的破碎等，造成較明顯的湍流現(xiàn)象，進(jìn)而增強(qiáng)了氣相傳輸速度。從雙膜理論方面考慮，氣相流速的增大可以使氣相薄膜厚度變薄，使其氣相阻力減小，因而甲苯分子可以較快的傳遞至氣液交界面。從不同流速下的甲苯質(zhì)量分?jǐn)?shù)值（圖7）可以看出，流速較大時(shí)，表面微層中甲苯的含量值較高，隨著氣相流速的減小，甲苯的質(zhì)量分?jǐn)?shù)也明顯降低。

圖7 不同氣相流速甲苯的質(zhì)量分?jǐn)?shù)變化趨勢(shì)

綜上考慮，氣相流速的增大增加了甲苯的遷移通量，同時(shí)縮短了達(dá)到氣液平衡的時(shí)間，因此，在條件允許的情況下，可以適當(dāng)加大氣相流速，來(lái)提高氣液吸收效果。

2.3 氣相濃度對(duì)甲苯遷移的影響

關(guān)于進(jìn)口氣相濃度對(duì)甲苯在氣液界面遷移的影響，筆者考慮甲苯進(jìn)口質(zhì)量濃度分別為1 000、4 000、8 000 mg/m3，同時(shí)假設(shè)氣相進(jìn)口流速為0.10 m/s，液相進(jìn)口流速為0.01 m/s，氣液分配系數(shù)為2保持不變，取液相深度為100 μm處的界面為研究對(duì)象，模擬結(jié)果顯示，在初始階段，不同進(jìn)口氣相濃度下，甲苯的吸收效率增長(zhǎng)均很快，其中8 000 mg/m3的情況最為明顯，即氣相濃度越高，甲苯的吸收效率越高（圖8）。

隨著時(shí)間的推移，甲苯的吸收效率逐漸下降，這與液相吸收劑對(duì)甲苯的吸收容量有關(guān)，另外，氣相濃度的增加同時(shí)又會(huì)造成氣相阻力層變厚，傳質(zhì)阻力增大，因此傳質(zhì)會(huì)向著不利于吸收的方向發(fā)展。

喬運(yùn)娟也對(duì)甲苯進(jìn)氣濃度對(duì)吸收效果的影響進(jìn)行了研究[7]，

圖8 1～4 s時(shí)不同氣相濃度下甲苯質(zhì)量分?jǐn)?shù)濃度

圖9 不同進(jìn)氣濃度對(duì)吸收效率的影響

從圖9可以看出，吸收開(kāi)始30 min之前，進(jìn)氣濃度越高，吸收效率越高，這與該研究所模擬的結(jié)果較為相似；而30 min之后，甲苯的吸收效率出現(xiàn)明顯降低，但濃度為1 000 mg/m3的情況下，吸收效率的改變比較平緩，吸收效果比較穩(wěn)定。但在該研究模擬結(jié)果中，隨著時(shí)間的變化，不同濃度值下吸收效率的變化趨勢(shì)是相似的。出現(xiàn)這種差異的原因主要是喬運(yùn)娟在試驗(yàn)過(guò)程中，添加了乳化劑，達(dá)到一定吸收程度后，氣相濃度越大，對(duì)甲苯的吸收抑制作用也就越大，因此會(huì)出現(xiàn)吸收效率下降情況的不同，而該研究的模擬僅考慮了較為理想的物理吸收過(guò)程，并沒(méi)有考慮其他影響因素，如化學(xué)反應(yīng)、周?chē)h(huán)境條件等，這為之后的進(jìn)一步研究和工作提供了思路。

3 結(jié)論

通過(guò)流體力學(xué)計(jì)算軟件，構(gòu)建適合該試驗(yàn)的物理模型，研究了表面微層中甲苯的遷移規(guī)律，分析了液相流速、氣相流速對(duì)其的影響，主要研究結(jié)果如下。

（1）液相流速對(duì)表面微層中甲苯遷移變化的影響主要表現(xiàn)為當(dāng)流速增大時(shí)，液相中甲苯的擴(kuò)散速率隨之增大，但同時(shí)影響了液膜即表面微層的形成，從而減小了甲苯的遷移通量，因此，要綜合考慮這兩方面的影響因素，選擇較為合適的液相流速，進(jìn)而達(dá)到最好的吸收傳質(zhì)效果。

（2）氣相流速的增大加快了甲苯分子到達(dá)氣液界面，縮短了達(dá)到氣液平衡的時(shí)間，同時(shí)增加了甲苯的遷移通量，因此，在條件允許的情況下，可以適當(dāng)加大氣相流速，來(lái)提高氣液吸收效果。

（3）在吸收開(kāi)始的最初階段，氣相濃度越高，表面微層中甲苯的吸收效率增長(zhǎng)速度也越快，隨著時(shí)間的變化，甲苯的吸收效率出現(xiàn)下降。

參考文獻(xiàn)

[1]

吳慶輝.溶劑吸收法在鞋業(yè)“三苯”廢氣治理中的應(yīng)用[J].福建環(huán)境，1999，16（4）：13-15.

[2] 劉英.苯系物的新型吸收劑制備及性能研究[D].青島：山東科技大學(xué)，2010.

[3] 戴樹(shù)桂，黃國(guó)蘭，雷紅霞.水體表面微層的環(huán)境化學(xué)研究[J].環(huán)境化學(xué)，1994，13（4），287-295.

[4] HARDY J T，APTS C W，CRECELIUS EA et al.Sea-surface microlayer metals enrichments in an urban and rural bay[J].Estuariune，Coastal and Shef Science，1985，20：299-312.

[5] HARDY J T，APTS C W，CRECELIUS EA et al.The sea-surface microlayer： fate and residence times of atmospheric metals[J].Limnd Oceanogy，1985，30（1）：93-101.

[6] LISS PS，DUCE RA.The Sea Surface and Global Change[M].Cambridge：Cambridge University Press，1997.

[7] 喬運(yùn)娟.液膜法吸收有機(jī)廢氣試驗(yàn)研究[D].杭州： 浙江大學(xué)，2008.