叢曉春，王 成，曲京華

(山東科技大學 山東省土木工程防災減災重點實驗室，山東 青島 266590)

?

水平采樣直管內顆粒傳輸特性的研究

叢曉春，王 成，曲京華

(山東科技大學 山東省土木工程防災減災重點實驗室，山東 青島 266590)

以0.1～20 μm的單分散顆粒為研究對象，利用數值模擬的方法，對在層流狀態(tài)下，顆粒在內徑為4、7和10 mm，長度為0.2、0.5、1.0和1.5 m的水平采樣直管內的運動特性進行模擬，得到水平采樣直管內的流場及顆粒的濃度分布，通過對模擬結果的整理分析得到顆粒在水平采樣直管內的傳輸效率。數值模擬發(fā)現(xiàn)，粒子的損失主要來源于重力和擴散作用；顆粒的傳輸效率要受到粒徑、管長、管內徑的影響；在層流狀態(tài)下，0.1～1.0 μm顆粒的傳輸效率幾乎達到100%，不受管內徑、管長和重力、布朗擴散力的影響；2.5～20 μm顆粒的傳輸效率隨粒徑、管內徑、管長的增大而減小。分析可得，通過減少采樣管的輸送長度或者在流量一定的基礎上減小管內徑均可提高顆粒在水平采樣直管內的傳輸效率。通過對比分析模擬結果與理論公式，可以近似認為顆粒在水平采樣直管內的傳輸效率等于在重力和擴散力單獨作用下顆粒傳輸效率的乘積。

顆粒物；采樣管；數值模擬；傳輸效率

隨著經濟的高速發(fā)展，城市化、工業(yè)化規(guī)模的擴大，空氣中各種環(huán)境污染物的排放量迅速增加，以顆粒物尤其是細顆粒物污染為主[1]。據現(xiàn)有文獻可知，對于顆粒物運動特性的研究主要集中在工業(yè)管道和空調系統(tǒng)的通風管道[2-4]，有關顆粒物在采樣儀器輸送管路中的傳輸及沉降特性的研究很少。本文以層流狀態(tài)下的采樣輸送管路中的水平直管段為研究對象，探討顆粒物在采樣管內的運動及沉降特性，分析管內徑、管長、粒徑等基本參數對傳輸效率的影響，為正確評價實際的采樣結果和大氣顆粒物的監(jiān)測提供技術支持。

1 數值模擬參數

本文通過數值模擬0.1、0.5、1.0、2.5、5.0、10、20 μm的單分散顆粒在層流(2.83 L/min)狀態(tài)下，在外徑為6.35、9.53、12.70 mm，壁厚1～2 mm(即內徑分別為4、7、10 mm)的水平直管內的運動，來了解0.1～20 μm的顆粒物在采樣管中的沉降及傳輸特性。數學模型及數值模擬研究所需詳細參數如表1所示。

對于顆粒相，在本模擬中采用可溯源至美國標準技術研究院(NIST)的單分散球形標準粒子(PSL小球)，其密度設定為美國DUCK標準小球的密度，即1.05 g/cm3，濃度采用2 000 個/cm3，這樣則可以忽略顆粒在輸送過程中顆粒的凝聚作用，其相關參數見表2。

2 模擬結果與分析

2.1 速度場分布

以內徑為7 mm的水平采樣管的模擬結果為例進行管內速度場的分析。在數值模擬過程中，設置入口流速是均勻分布的，從圖1中可以看出，在采樣管入口截面，氣流均勻分布，只是管道外邊壁處的流速幾乎為零；在出口截面上，流速成明顯的同心圓分布，且趨于穩(wěn)定狀態(tài)，管道中心的流速最大，靠近管道外壁面的流速幾乎為零，此時的邊界層比較薄，管內流體分層嚴重。這種速度分布恰好與圓管內層流的流動規(guī)律及邊界層理論相符。

表1 數學模型的主要參數

表2 顆粒相的主要參數

2.2 顆粒濃度場分布

以內徑為7 mm的水平采樣管為例，對顆粒的濃度場進行分析，層流狀態(tài)下管內顆粒濃度(kg/m3)分布如圖2所示，而其他工況則以顆粒在管內的傳輸效率的形式進行描述，如圖3所示。

在采樣管的入口處，我們設定顆粒是均勻分布的。從圖2可以看出，在顆粒進入采樣管之后，顆粒的運動受到管內氣流的影響，粒徑為0.1～1.0 μm顆粒受管內氣流的影響較大，其隨動性較強，受重力的影響不大；粒徑為2.5～ 20 μm的顆粒受氣流的影響較小，相對來說，其隨動性較差，在跟隨氣流的運動過程中，顆粒受到自身的重力影響較大。

從圖3可知，在層流狀態(tài)下，隨著顆粒粒徑的增大，其在管內的傳輸效率逐漸減小，其機理是重力作用。對于0.1～1 μm的粒子來說，其傳輸效率幾乎達到100%，管內徑、管長和重力、布朗擴散力對其產生的影響較小，這與文獻[5]中的結論是相同的；對于2.5～20 μm的粒子而言，在同一管長的采樣管內，隨著管內徑的增大，其傳輸效率逐漸減小，這是由于管內徑較大，粒子更可能碰撞到管道壁面而被捕集；在同一內徑的采樣管內，管長越大，其傳輸效率越小，這種影響隨粒徑的增大而越明顯。例如，內徑為4 mm的采樣管內，當管長為0.2 m時，5 μm粒子的傳輸效率為98.7%，20 μm粒子的傳輸效率為80.9%；當管長為1.0 m時，5 μm粒子的傳輸效率降低為93.6%，下降了5.2%，而20 μm粒子的傳輸效率則降低為23.2%，下降了71.3%。Liu[6]也得到了相同的變化趨勢。

2.3 理論公式對比分析

在圓形水平管道的層流中，F(xiàn)uchs和Thomas[7]假設管內氣流分布呈拋物線形，這樣就獨立解決了有關重力沉降所帶來的問題，最先得出了層流狀態(tài)下粒子在圓形水平管道中重力沉降的傳輸效率，隨后Heyder和Gebhart[8]進一步的研究，得到了顆粒重力沉降的傳輸效率的一般關系式：

(1)

(2)

式中：Z—重力沉降參數；L—管道長度，m；d—管道內徑，m。

在布朗擴散的作用下，粒子會在管壁上發(fā)生沉積，顆粒擴散損失的傳輸效率可以定義如下：

(3)

(4)

(5)

式中：Vdiff-顆粒擴散沉積速度，也稱為質量傳遞系數；Sh-舍伍德數，無量綱的質量轉移系數；D-顆粒的擴散系數。

利用上述表達式計算顆粒擴散傳輸效率的理論結果與Gormley和Kennedy的層流管道顆粒擴散沉積實驗的結果差別甚微，能夠很好的吻合[9]。

本文分別在只考慮重力作用、只考慮布朗擴散作用以及同時考慮兩者共同作用的條件下，對不同粒徑的顆粒在長度為0.5 m，內徑為4 mm水平采樣直管內的運動及傳輸效率進行數值模擬，將得到的模擬結果與上述理論公式的計算結果進行如下比較，如表3所示。

從表3中看出，當單獨考慮重力、單獨考慮布朗擴散力的影響時，通過數值模擬得到的結果與理論計算結果比較相近；當同時考慮重力和布朗擴散力的共同影響時，通過數值模擬得到的結果與重力、布朗力單獨作用時理論傳輸效率乘積之間的偏差都在10%以內。因此，在層流狀態(tài)下，可以近似認為顆粒在水平直管中的傳輸效率是其在重力、布朗擴散力單獨作用下的傳輸效率之積。

表3 顆粒傳輸效率的模擬結果與理論公式計算結果比較(層流)

3 結語

1)在層流狀態(tài)下，水平圓形采樣管中流體的速度沿流動方向逐漸增大，靠近管道壁面的速度較小，管道中心處的速度最大。

2)在水平采樣管中，層流狀態(tài)下，0.1～1.0 μm顆粒的傳輸效率幾乎達到100%，不受管內徑、管長和重力、布朗擴散力的影響；2.5～20 μm顆粒的傳輸效率隨粒徑、管內徑、管長的增大而減小。通過減少采樣管的輸送長度或者在流量一定的基礎上減小管內徑，均可提高顆粒在水平采樣管內的傳輸效率。

3)在層流狀態(tài)下，顆粒在水平采樣管內總的傳輸效率近似于重力和擴散力單獨作用下顆粒傳輸效率的乘積。

[1]張 霞，孟琛琛，王麗濤，等.邯鄲市大氣污染特征及變化趨勢研究[J].河北工程大學學報：自然科學版，2015，32(4)：69-74.

[2]JIANG HAI，LU LIN，SUN KE.Experimental study and numerical investigation of particle penetration and deposition in 90°bent ventilation ducts[J].Building and Environment,2011(46)：2195-2202.

[3]張 錦，劉澤華.通風管道中氣溶膠顆粒沉降的探討[J].節(jié)能與環(huán)保，2009(6)：37-39.

[4]TANG YINGJIE，GUO BING，ANDREW R MCFARLAND.A Computational Fluid Dynamics Study of Particle Penetration through an Omni-Directional Aerosol Inlet[J].Aerosol Science and Technology，2010(44)：1049-1057.

[5]巴 倫，維勒克.氣溶膠測量原理、技術及應用[M].白志鵬，張燦等譯.北京：化學工業(yè)出版社，2007.

[6]LIU DELING.Evaluating Aerosol Transmission Efficiency in a Sampling System[J].Journal of the IEST，2013，56(2)：29-42.

[7]Fuchs N.A.Sampling of Aerosols[J].Atmos Environ，1975(9)：697-707.

[8]HEYDER J, GEBHART J.Gravitational deposition of particles from laminar aerosol flow through inclined circular tubes[J].Journal of Aerosol Science，1977(8): 289-295.

[9]GORMLEY P G，KENNEDY M.Diffusion from a stream flowing through a cylindrical tube[J].Royal Irish Academy，1949(52)：163-169.

(責任編輯王利君)

The study of transmission characteristics of particles in Straight horizontal sampling tube

CONG Xiaochun, WANG Cheng, QU Jinghua

(Shandong Provincial Key Laboratory of Civil Engineering Disaster Prevention and Mitigation, Shandong University of Science and Technology, Shandong Qingdao 266590, China)

This paper chooses 0.1 μm～20 μm monodisperse particles as the research object. We adopted numerical simulation methods to study the transmission characteristics of particles in the straight horizontal sampling tube under laminar flow. Diameters of tube are 4 mm, 7 mm and 10 mm. The lengths of tube are 0.2 m, 0.5 m, 1.0 m and 1.5 m. We got the flow field of the tube and the particle concentration distribution. Through the analysis of simulation results, we obtained the transport efficiency of particles in the straight horizontal sampling tube. Based on the numerical simulation of particle transmission characteristics within the straight horizontal sampling tube, it is found that the particle loss is mainly from gravity and diffusion. Simultaneously, the particle transport efficiency is influenced by particle size, tube length and inner diameter. Under laminar flow, the transport efficiency of 0.1 μm～1.0 μm particles is almost 100%, without the influence of inner diameter, tube length , gravity and Brownian diffusion force. The transport efficiency of 2.5 μm～20 μm particles decreases with the increasing of particle size, tube length and inner diameter. Through the analysis, it will improve the transport efficiency of particles by reducing the length of the sampling tube or the inner diameter on the basis of a flow.The simulation results compared with the theoretical formula, it can be approximated that the transport efficiency of particles is equal to the product of particle transport efficicency respectively under the action of gravity and diffusion.

particulate matter; sampling tube; numerical simulation; transport efficicency

2015-12-15

國家科技支撐計劃項目(2013BAK12B08)

叢曉春(1974-)，女，山東威海人，博士，副教授，從事顆粒物污染防治方面的研究。

1673-9469(2016)02-0078-04

10.3969/j.issn.1673-9469.2016.02.017

X851

A