華鳳皎,亢燕銘,鐘 珂 (東華大學(xué)環(huán)境科學(xué)與工程學(xué)院,上海 201620)

湍流效應(yīng)對(duì)雨滴表面粒子捕集過程的影響

華鳳皎,亢燕銘*,鐘 珂 (東華大學(xué)環(huán)境科學(xué)與工程學(xué)院,上海 201620)

為考察湍流運(yùn)動(dòng)對(duì)捕集效率的影響,本文在不同湍流強(qiáng)度下分別模擬了 4種直徑雨滴對(duì)不同粒徑粒子的捕集過程.根據(jù)湍流場(chǎng)中粒子捕集過程的特征,通過引入碰撞概率的概念,重新定義了適用于湍流狀態(tài)下捕集效率的計(jì)算公式,并在此基礎(chǔ)上分析研究了湍流捕集效應(yīng)對(duì)捕集效率的影響.結(jié)果表明,所有粒徑粒子的碰撞概率分布特征都會(huì)受到背景湍流強(qiáng)度的影響,尤其對(duì) 5μm及以下粒子影響明顯.湍流捕集效應(yīng)對(duì)捕集效率的影響與粒子密度和粒徑直接相關(guān),與雨滴尺寸無(wú)關(guān).

氣溶膠粒子；雨滴；湍流強(qiáng)度；捕集效率

近幾年來(lái),霾天氣在我國(guó)多數(shù)大城市頻繁出現(xiàn),已成為熱點(diǎn)環(huán)境問題.顯然,大氣的自然凈化過程是霾的重要清除機(jī)制之一.在自然凈化過程中,氣溶膠粒子主要通過干沉積和濕清除從大氣中被移除.盡管干沉積的貢獻(xiàn)率更高、更有效,但濕清除在降雨過程中單位時(shí)間內(nèi)的清除效率更高.在多雨地區(qū),濕清除對(duì)大氣環(huán)境的凈化作用不可忽略.目前在全球和區(qū)域大氣氣溶膠輸送模式及空氣質(zhì)量評(píng)估模式等的研究中[1-4],尤其關(guān)注氣溶膠濕清除率,并認(rèn)為其計(jì)算的準(zhǔn)確性將影響到這些模式對(duì)氣溶膠輸送過程估計(jì)的精度.

Slinn公式是目前常用的雨滴對(duì)粒子捕集效率的半經(jīng)驗(yàn)公式[5],但它是在以雨滴表面繞流為層流的基礎(chǔ)上得到的,而大氣運(yùn)動(dòng)具有明顯的湍流特征,大氣湍流強(qiáng)度甚至可能達(dá)到 50%[6,7].背景大氣流場(chǎng)的湍動(dòng)效應(yīng)通過影響粒子軌跡來(lái)改變雨滴對(duì)粒子的捕集效果.盡管最近幾年來(lái),已有不少研究者指出[8-11]應(yīng)該考慮這種湍流擾動(dòng)對(duì)雨滴捕集氣溶膠粒子的影響,但相應(yīng)的研究結(jié)果目前還很少見到.通過分析背景湍流捕集效應(yīng)對(duì)雨洗過程中粒子軌跡的干擾和畸變,得出慣性粒子在湍流狀態(tài)下被雨滴捕獲的效率將大于層流時(shí)的捕集效率,但沒有對(duì)湍流效應(yīng)產(chǎn)生的具體影響做進(jìn)一步研究[12].

計(jì)算流體力學(xué)(CFD)方法是求解湍流運(yùn)動(dòng)的常用手段,通過耦合求解流動(dòng)方程和粒子運(yùn)動(dòng)方程,可得到湍流狀態(tài)下雨滴對(duì)粒子的捕集效率.與實(shí)測(cè)相比,CFD模擬可節(jié)省大量實(shí)驗(yàn)成本,并極大地縮短問題求解時(shí)間從而降低研究成本.因此,本文將采用經(jīng)過實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證后的CFD模型進(jìn)行上述問題的求解.由于不同湍流強(qiáng)度下,氣流速度的脈動(dòng)幅度不同,進(jìn)而導(dǎo)致粒子隨湍流脈動(dòng)隨機(jī)游走的過程中運(yùn)動(dòng)軌跡差別較大,最終造成湍流對(duì)粒子捕集效率的作用受湍流特性和粒子特性的影響很大,故本文將針對(duì)不同尺寸的雨滴和粒子來(lái)分析湍流強(qiáng)度對(duì)捕集過程和捕集效率的具體影響,以便為修正清除效率估計(jì)式(即Slinn公式)提供必要的科學(xué)依據(jù).

1 數(shù)值方法與驗(yàn)證

1.1 物理模型和計(jì)算工況

實(shí)際大氣中雨滴以終端沉降速度u0下落,為了簡(jiǎn)化問題,將下落雨滴視為靜止液體球,氣流以恒定速度 u0流經(jīng)雨滴.考慮雨滴球體的對(duì)稱性,取雨滴的1/4進(jìn)行建模,如圖1所示.計(jì)算域入口邊界與雨滴前向駐點(diǎn)間的距離取為雨滴直徑的5倍,計(jì)算域出口與雨滴后向駐點(diǎn)的距離取為 8倍直徑.為避免計(jì)算區(qū)域端部效應(yīng)的影響,雨滴表面與計(jì)算域其余邊界的距離均取5倍直徑,將緊鄰雨滴表面處的網(wǎng)格做加密處理,網(wǎng)格最小尺寸為1.0μm.計(jì)算區(qū)域與網(wǎng)格處理詳見文獻(xiàn)[12-13].

雨滴終端速度[1]u0(dw)如下:

式中:dw為雨滴直徑,cm.針對(duì)本文研究的雨滴直徑,表1給出了不同雨滴對(duì)應(yīng)的終端速度.

湍流強(qiáng)度是衡量氣流湍流脈動(dòng)大小的參數(shù),氣溶膠粒子隨氣流運(yùn)動(dòng)過程中受到大氣湍流脈動(dòng)和氣流速度的影響,從而增加粒子運(yùn)動(dòng)軌跡的隨機(jī)脈動(dòng)性.為研究湍流強(qiáng)度對(duì)雨滴-粒子碰撞概率和捕集效率的影響,模擬了 5種背景湍流強(qiáng)度I0情況下dw=800μm雨滴分別對(duì)5種粒徑粒子(dp=2.5、3.5、5、7和10μm)進(jìn)行捕集;在研究雨滴直徑對(duì)湍流捕集效應(yīng)的影響時(shí),選取相同的湍流強(qiáng)度I0=15%,模擬4種雨滴直徑(dw=100、200、 400和800μm)對(duì)5種粒徑粒子(dp=2.5、3.5、5、7和 10μm)的捕集;討論粒子密度對(duì)湍流捕集效應(yīng)的影響時(shí),模擬湍流強(qiáng)度為 I0=50%,雨滴直徑dw=800μm 時(shí),5種粒徑粒子(dp=2.5、3.5、5、7和10μm)密度分別為ρp=1500和2800kg/m3時(shí)被雨滴捕集的效率.

圖1 計(jì)算模型示意Fig.1 Schematic diagram of the computational model

表1 不同直徑雨滴對(duì)應(yīng)的終端沉降速度Table 1 Terminal velocities corresponding to different diameter raindrops

1.2 計(jì)算方法與邊界條件

本文采用標(biāo)準(zhǔn) k?ω模型[14-15]對(duì)雨滴附近湍流流場(chǎng)進(jìn)行模擬,其中壁面湍流特征用標(biāo)準(zhǔn)壁面函數(shù)描述,采用有限容積法離散控制方程,對(duì)離散方程的差分采用二階迎風(fēng)格式.假定計(jì)算域入口流速均勻,入口邊界設(shè)為 inlet,出口邊界設(shè)為outflow,即該流出區(qū)域的質(zhì)量流率和流入?yún)^(qū)域的質(zhì)量流率相等.將雨滴視作硬球,即表面無(wú)滑移,計(jì)算域的壁面設(shè)為有滑移壁面.粒子視作實(shí)心球,并假設(shè)粒子間無(wú)凝并行為.為了與經(jīng)典實(shí)驗(yàn)數(shù)據(jù)比較并驗(yàn)證數(shù)值模型以獲得簡(jiǎn)明結(jié)果,在數(shù)值程序中,僅考慮重力、布朗力和Saffman升力.控制方程如下:

連續(xù)性方程

動(dòng)量守恒方程

湍流動(dòng)能k和擴(kuò)散速率ω的輸運(yùn)方程

粒子運(yùn)動(dòng)的控制方程

式中:Fg是重力項(xiàng);FB為布朗力;FL為Saffman升力,表達(dá)式如下:

式中:ζ是均值為0、方差為1的高斯隨機(jī)函數(shù);Δt是顆粒的時(shí)間步長(zhǎng);S0是光譜強(qiáng)度函數(shù);K=2.594; dij為流體變形速率張量.

1.3 模型的合理性驗(yàn)證

已有不少研究者對(duì)于層流條件下球形捕集體上的粒子慣性捕集進(jìn)行了實(shí)驗(yàn)研究[17-21].為驗(yàn)證本文采用程序的可靠性,以下以雨滴直徑dw=800μm 為例,模擬了層流下慣性粒子在雨滴表面的捕集效率.圖2給出了模擬計(jì)算結(jié)果與已有文獻(xiàn)實(shí)測(cè)數(shù)據(jù)的比較(圖中E為捕集體對(duì)粒子的捕集效率,St為粒子的Stokes數(shù)[22]).

圖2 球形捕集體上粒子慣性捕集效率的比較Fig.2 Comparison of calculated efficiencies of aerosols on a sphere by impaction vs previous results

由圖2可見,在只考慮慣性碰撞機(jī)制時(shí),捕集效率的模擬結(jié)果與文獻(xiàn)中給出的一系列實(shí)測(cè)值基本一致.數(shù)值結(jié)果還表明,與 Slinn公式[5]求得的慣性效率相比,本文計(jì)算曲線更接近實(shí)測(cè)結(jié)果,而Langmuir和Blodgett的結(jié)果與所有實(shí)驗(yàn)值都偏離較遠(yuǎn).因此,本文采用的數(shù)值程序可以用于對(duì)粒子運(yùn)動(dòng)軌跡的計(jì)算和雨滴對(duì)粒子捕集效率的模擬和估計(jì).

2 結(jié)果與討論

2.1 背景湍流強(qiáng)度對(duì)雨滴附近湍流場(chǎng)的影響

圖3給出了背景湍流強(qiáng)度分別為I0=25%和50%時(shí),dw=800μm 雨滴附近流場(chǎng)的湍流強(qiáng)度及流線的分布特征.

從上圖中容易看出,雨滴上半游的氣流湍流強(qiáng)度自入口向遠(yuǎn)場(chǎng)擴(kuò)散,當(dāng)雨滴對(duì)氣流運(yùn)動(dòng)產(chǎn)生阻礙作用時(shí),氣流擾動(dòng)劇烈,湍流強(qiáng)度顯著增強(qiáng).此外,盡管湍流強(qiáng)度不同,但雨滴附近的湍流強(qiáng)度等值線輪廓和流線相似,表明平均流場(chǎng)特征不受背景湍流強(qiáng)度的影響.

圖3 雨滴附近流場(chǎng)湍流強(qiáng)度及流線分布Fig.3 Distributions of the turbulent intensity and streamlines near a raindrop

圖4 雨滴迎風(fēng)面附近湍流強(qiáng)度與背景流場(chǎng)湍流強(qiáng)度的相關(guān)性比較Fig.4 Comparison of correlations between turbulent intensity near windward side of a raindrop and background turbulent intensity

圖4給出了雨滴迎風(fēng)面附近(距離雨滴迎風(fēng)面中心0.01dw處)投影面積上的平均湍流強(qiáng)度Idw與背景湍流強(qiáng)度I0的相關(guān)性,可以看到,雨滴附近的湍流強(qiáng)度隨著背景湍流強(qiáng)度的增大而增加.當(dāng)背景湍流強(qiáng)度較小時(shí),雨滴附近湍流強(qiáng)度遠(yuǎn)大于背景湍流強(qiáng)度,且隨著背景湍流強(qiáng)度的增大,二者逐漸接近,直至背景湍流強(qiáng)度大于雨滴附近湍流強(qiáng)度.不同的背景湍流強(qiáng)度將造成不同的局部氣流擾動(dòng),從而直接導(dǎo)致雨滴對(duì)粒子的捕集效果發(fā)生變化.

2.2 雨滴捕集粒子的湍流捕集效應(yīng)

層流狀態(tài)下,粒子的運(yùn)動(dòng)軌跡與釋放位置一一對(duì)應(yīng),因此存在如圖1所示的臨界位置y0,即粒徑為dp的粒子極限軌跡對(duì)應(yīng)的釋放位置.所有在中心對(duì)稱軸線的距離大于 y0的位置處的釋放粒子,其運(yùn)動(dòng)軌跡將全部繞過雨滴逃逸,小于y0位置釋放的粒子軌跡全部碰撞到雨滴表面被捕集.為此,層流狀態(tài)時(shí)雨滴對(duì)粒子的捕集效率 EL(dw,dp)的定義為:

式中:dw為雨滴直徑,μm;dp為粒子直徑,μm;N為被捕集的粒子數(shù),個(gè);N和N0分別為被捕集的粒子數(shù)和投影面積上釋放的總粒子數(shù),個(gè).

但是在湍流流場(chǎng)中由于受到氣流脈動(dòng)的影響,粒子軌跡將不具有層流狀態(tài)時(shí)的唯一性,而具有一定的隨機(jī)性,即在層流狀態(tài)下無(wú)法被捕集的粒子,在湍流背景場(chǎng)中卻可能被雨滴捕獲,而在層流狀態(tài)下一定會(huì)被捕集的粒子在湍流背景下也有可能不會(huì)被捕集.為此,本文引入碰撞概率ε(dw,dp,Y)來(lái)分析湍流下捕集效率的微觀圖像,其中無(wú)量綱釋放距離 Y=y/y0.碰撞概率定義為從位置y處發(fā)出的粒徑為dp的粒子與直徑為dw雨滴發(fā)生碰撞而被捕集的可能性.

圖5給出了不同背景湍流強(qiáng)度下dw=800μm的雨滴對(duì)不同粒徑粒子的碰撞概率沿?zé)o量綱釋放位置Y(Y=y/y0)的分布圖.由圖可以看到,碰撞概率隨著粒子無(wú)量綱釋放位置Y減小而增大.這是因?yàn)榱Ｗ舆\(yùn)動(dòng)軌跡首先受到慣性力控制,湍流在此基礎(chǔ)上通過速度脈動(dòng)影響粒子運(yùn)動(dòng)軌跡,故靠近雨滴中心對(duì)稱軸線釋放的粒子的碰撞概率較大.對(duì)于大粒子,碰撞概率都在Y小到一定值時(shí)增加到 100%,而對(duì)于小粒子,慣性力作用較弱,湍流干擾對(duì)粒子運(yùn)動(dòng)軌跡影響非常大,使得即使在Y=0處釋放的粒子碰撞概率也無(wú)法達(dá)到100%.

由圖5還可以看到,隨著背景湍流強(qiáng)度的增大,湍流擾動(dòng)作用相對(duì)于慣性力的作用增強(qiáng),導(dǎo)致碰撞概率達(dá)到100%對(duì)應(yīng)的臨界位置Yc不斷減小,使碰撞概率隨粒子釋放位置 Y的變化曲線更加偏離層流狀態(tài)時(shí)的特征.

圖5 湍流強(qiáng)度對(duì)雨滴-粒子碰撞概率的影響Fig.5 Influence of turbulent intensities on the collision probabilities between a raindrop and aerosols

上述模擬計(jì)算結(jié)果表明,在湍流脈動(dòng)作用下,從雨滴上游任意位置y處發(fā)出的粒子,碰撞到雨滴表面的概率都介于 0～100%.這就意味著,考慮了湍流捕集效應(yīng)后,層流下的捕集效率預(yù)測(cè)公式(12)不再適用.需要根據(jù)湍流狀態(tài)下粒子與雨滴的碰撞概率特征重新定義適用于湍流場(chǎng)的捕集效率.

2.3 湍流脈動(dòng)對(duì)捕集效率的影響

由于不同的釋放位置 Y對(duì)應(yīng)的碰撞概率以及對(duì)應(yīng)的雨滴投影面積均不同,本文將雨滴的整個(gè)圓形投影區(qū)域分成若干個(gè)捕集微元環(huán)dS,如圖6所示.

圖6 湍流捕集效率計(jì)算示意Fig.6 Schematic diagram of calculating turbulent collection efficiency

釋放位置 y+dy對(duì)應(yīng)的投影圓環(huán)面積為2πydy,球形捕集體周圍流場(chǎng)關(guān)于軸對(duì)稱的特性使得Y值與碰撞概率ε(dw,dp,Y)為一一對(duì)應(yīng)關(guān)系.根據(jù)對(duì)稱性原理,每個(gè)捕集微元環(huán)對(duì)應(yīng)的碰撞概率可以認(rèn)為是一個(gè)定值.因此,湍流場(chǎng)中粒子的捕集效率可表示為:

通過數(shù)值模擬計(jì)算可以得到湍流場(chǎng)中在不同Y值釋放的粒子與雨滴的碰撞概率,為便于利用數(shù)值結(jié)果計(jì)算雨滴對(duì)粒子的捕集效率,將上式改寫為:式中:y1,y2,…… yn分別為第1,2,…n個(gè)捕集微元環(huán)外圓半徑,μm;ε1,ε2,……εn分別為 y1,y2,……yn位置處對(duì)應(yīng)的碰撞概率.

根據(jù)圖 5的數(shù)值模擬結(jié)果和式(14),可以計(jì)算得到不同湍流強(qiáng)度時(shí)粒子的捕集效率.圖7給出了5種湍流強(qiáng)度下dw=800μm雨滴對(duì)不同粒徑粒子的捕集效率.

由圖 7可知,較小粒子(dp=2.5～5μm)捕集效率隨著湍流強(qiáng)度的增大而升高,并且粒徑越小,湍流對(duì)捕集效率的增大作用越顯著.而較大粒子盡管碰撞概率曲線(圖5)受湍流的影響程度隨湍強(qiáng)的增大而增大,但捕集效率卻幾乎不受湍流強(qiáng)度影響.相同的湍流強(qiáng)度下,捕集效率均隨著粒徑的增大而增加.

圖7 不同湍流強(qiáng)度下粒子捕集效率的比較Fig.7 Comparison of collection efficiencies of aerosols with different turbulent intensities

2.4 雨滴尺寸對(duì)湍流捕集效應(yīng)的影響

由于5μm以上大粒子的捕集效率受湍流影響不明顯,因此,本文針對(duì)2.5和5μm粒子研究雨滴尺寸對(duì)上述湍流捕集效應(yīng)的響應(yīng).圖8分別給出了背景湍流強(qiáng)度I0=15%時(shí),dp=2.5和5μm粒子與不同直徑雨滴的碰撞概率曲線.

圖8的結(jié)果表明,dp=2.5μm的粒子在大雨滴上的碰撞概率略高于在小雨滴上的碰撞概率,這是由于dp=2.5μm粒子在靠近雨滴中軸線釋放時(shí),受布朗擴(kuò)散作用影響較大,所以,被不同直徑雨滴捕集的概率隨機(jī)性亦較大.但總體上看,雨滴尺寸對(duì)粒子的碰撞概率影響不大.原因是阻礙物周圍流場(chǎng)特征主要受其形狀影響,而本文研究的各種尺寸的雨滴均視為球形.因此,不同尺寸雨滴附近流場(chǎng)湍流特征也基本相同,故捕集效率對(duì)湍流效應(yīng)的響應(yīng)便近似相同.

圖8 雨滴直徑對(duì)湍流捕集效應(yīng)的影響Fig.8 Influence of raindrop diameters on turbulent collection effect

2.5 粒子特性對(duì)湍流捕集效應(yīng)的影響

由圖5可以看出,在湍流流場(chǎng)中粒徑dp越大,受到的慣性作用越顯著,與雨滴發(fā)生碰撞的概率ε (dw,dp,Y)越接近層流下兩點(diǎn)垂直式的碰撞概率分布,表明湍流捕集效應(yīng)隨著粒徑增大而減小.

圖9比較了ρp=1500和2800kg/m3粒子在背景湍流強(qiáng)度為I0=50%時(shí)與dw=800μm雨滴的碰撞概率.由圖可知,粒子密度越小,湍流捕集效應(yīng)越明顯,但不同粒徑粒子受密度的影響程度不同.圖9(a)表明,當(dāng)dp=2.5μm,ρp=1500kg/m3時(shí),雨滴與粒子的碰撞概率很低,并且與釋放位置Y相關(guān)性很小.而ρp=2800kg/m3時(shí),碰撞概率明顯增高,且與釋放位置相關(guān)性很強(qiáng).這是因?yàn)榱Ｗ邮艿降膽T性力隨密度增大而增大,密度較小的dp=2.5μm粒子運(yùn)動(dòng)軌跡幾乎完全受控于湍流脈動(dòng),慣性力相對(duì)于湍流作用力較弱.由圖9(b)可知,對(duì)于大粒子,雖然較小密度會(huì)導(dǎo)致慣性力減小、碰撞概率受到的湍流捕集效應(yīng)增大,但由于粒徑較大,保證了粒子受到的流場(chǎng)慣性力不會(huì)過小,因此,小密度和大密度粒子碰撞概率分布規(guī)律均與釋放位置 Y密切相關(guān).

圖9 粒子密度對(duì)湍流捕集效應(yīng)的影響Fig.9 Influence of particle densities on turbulent collection effect

圖10 層流和湍流狀態(tài)下捕集效率的比較Fig.10 Comparison of the collection efficiencies between laminar and turbulent flows

圖10給出了I0=50%和層流下dw=800μm雨滴對(duì)不同密度粒子的捕集效率.可以看到在湍流作用下,大密度粒子和小密度粒子的捕集效率都有所提高.這是因?yàn)橄鄬?duì)于層流,盡管湍流捕集效應(yīng)降低了Y較小位置釋放粒子的碰撞概率,但增大了Y較大位置釋放粒子的碰撞概率,而后者對(duì)應(yīng)的圓環(huán)面積[圖 6(b)]大于前者,使得雨滴捕集體有效捕集面積增大,最終提高了粒子的捕集效率.

為了明確湍流對(duì)不同密度和不同粒徑粒子捕集效率的影響,定義粒子捕集效率的湍流增加率λ為:

式中:EL,ET分別為該粒徑粒子對(duì)應(yīng)的層流和湍流捕集效率.

圖11分別給出了dw=800μm的雨滴對(duì)不同密度粒子捕集效率的湍流增加率.由圖11可以看到,粒子的密度和粒徑dp越小,湍流脈動(dòng)對(duì)捕集效率的增加率λ越大.湍流對(duì)dp=10μm的2種密度粒子捕集效率的提高幅度分別為 8.3%和 9.5%.但對(duì)于dp=2.5μm的2種密度粒子,提高幅度高達(dá)72%和120%.

圖11 不同粒子密度下捕集效率的湍流增加率Fig.11 Turbulent increase percentages of collection efficiency with different particle densities

3 結(jié)論

3.1 由于湍流脈動(dòng)對(duì)粒子運(yùn)動(dòng)軌跡的隨機(jī)干擾,雨滴對(duì)不同位置釋放的粒子的捕集效果與層流完全不同.為此,原有的層流狀態(tài)下粒子捕集效率的計(jì)算公式不再適用于湍流場(chǎng),本文通過定義碰撞概率的概念,重新定義了適用于湍流狀態(tài)下的捕集效率計(jì)算公式.

3.2 所有粒徑粒子的碰撞概率分布特征都會(huì)受到背景湍流強(qiáng)度的影響,但捕集效率受到湍流強(qiáng)度明顯影響的粒徑范圍為dp≤5μm.

3.3 由于不同大小的球形阻礙物對(duì)周圍流場(chǎng)的干擾作用相似,因此,雨滴尺寸對(duì)碰撞概率分布規(guī)律的影響可以忽略.

3.4 隨著粒子密度和粒子直徑的減小,粒子的碰撞概率分布將越分散,湍流脈動(dòng)對(duì)捕集效率的增大作用越顯著.

[1] Wang X, Zhang L, Morna M D. Uncertainty assessment of current size-resolved parameterizations for below - cloud particle scavenging by rain [J]. Atmospheric Chemistry & Physics, 2010, 10(12):5685-5705.

[2] Wang X, Zhang L, Morna M D. On the discrepancies between theoretical and measured below-cloud particle scavenging coefficients for rain – a numerical investigation using a detailed one-dimensional cloud microphysics model [J]. Atmospheric Chemistry & Physics, 2011,11(22):11859-11866.

[3] Henzing J S, Olivie D J L, van Velthoven P F J. A parameterization of size resolved below cloud scavenging of aerosols by rain [J]. Atmospheric Chemistry & Physics, 2006, 6(11):3363-3375.

[4] Sportisse B. A review of parameterizations for modeling dry deposition and scavenging of radionuclides [J]. Atmospheric Environment, 2007,41(13):2683-2698.

[5] Slinn W G N. Precipitation scavenging. In: Atmospheric Sciences and Power Production 1979. Washington, D.C.: Division of Biomedical Environmental Research, U.S. Department of Energy, 1983.

[6] 高會(huì)旺,顧 明,王仁磊,等.北黃海海域大氣湍流強(qiáng)度特征及風(fēng)速標(biāo)準(zhǔn)差相似性分析 [J]. 中國(guó)海洋大學(xué)學(xué)報(bào), 2009,39(4):563-568.

[7] 徐玉貌,周朝輔,李振華,等.廣州市近地層大氣的湍流微結(jié)構(gòu)和譜特征 [J]. 大氣科學(xué), 1993,17(3):338-348.

[8] Feng J. A 3-mode parameterization of below-cloud scavenging of aerosols for use in atmospheric dispersion models [J]. Atmospheric Environment, 2007,41(32):6808-6822.

[9] Pinsky M B, Khain A P. Turbulence effects on droplet growth and size distribution in clouds - A review [J]. Journal of Aerosol Science, 1997,28(7):1177-1214.

[10] Andronache C, Gronholm T, Laakso L, et al. Scavenging of ultrafine particles by rainfall at a boreal site: observations and model estimations [J]. Atmospheric Chemistry & Physics, 2006, 6(12):4739-4754.

[11] Khain A P, Pinsky M. B. Turbulence effects on the collision kernel, II: Increase of the swept volume of colliding drops [J]. Quarterly Journal of the Royal Meteorological Society, 1997, 123(542):1543-1560.

[12] 華鳳皎,劉曉燕,亢燕銘. 雨滴表面慣性粒子捕集的湍流捕集效應(yīng)數(shù)值分析 [J]. 中國(guó)環(huán)境科學(xué), 2013,33(9):1585-1590.

[13] Hosseini S A, Vahedi Tafreshi H. Modeling particle-loaded single fiber efficiency and fiber drag using ANSYS–Fluent CFD code [J]. Computers & Fluids, 2012,66:157-166.

[14] Longest P W, Vinchurkar S. Validating CFD predictions of respiratory aerosol deposition: Effects of upstream transition and turbulence [J]. Journal of Biomechanics, 2007,40:305-316.

[15] Geng T, Longest P W, et al. Development of a stochastic individual path (SIP) model for predicting the tracheobronchial deposition of pharmaceutical aerosols: Effects of transient inhalation and sampling the airways [J]. Journal of Aerosol Science, 2011,42:781-799.

[16] Wilcox D C. Reassessment of the scale-determining equation for advanced turbulence model [J]. AIAA Journal, 1988,26(11):1299-1310.

[17] Waldenmaier M. Measurements of inertial deposition of aerosol particles in regular arrays of spheres [J]. Journal of Aerosol Science, 1999,30(10):1281-1290.

[18] Ranz W E, Wong J B. Impaction of dust and smoke particles on surface and body collectors [J]. Industrial & Engineering Chemistry, 1952,44(6):1371-1380.

[19] Walton W, Woolcock A. The suppression of airborne dust by water spray [J]. International Journal of Air Pollution, 1960,3: 129-153.

[20] Schmidt M, Loffler F. Experimental investigations on two- phase flow past a sphere using digital Particle – image – velocimetry [J]. Experiments in Fluids, 1993,14(5):296-304.

[21] Langmuir I, Blodgett K B. Mathematical investigation of water droplet trajectories [J]. Atmospheric Phenomena, 1961,1062:335-347.

[22] Friedlander S K. Smoke, Dust, and Haze: Fundamentals of Aerosol Dynamics. Oxford University Press, USA, 2000, p71.

Influence of turbulent effect on the collection process of aerosols by raindrops.

HUA Feng-jiao, KANG Yan-ming*, ZHONG Ke
(School of Environmental Science and Engineering, Donghua University, Shanghai 201620, China). China Environmental Science, 2017,37(1)：13～20

To explore the influence of turbulent motion on collection efficiency, the collecting processes of different diameter aerosols captured by four diameter raindrops were estimated by employing numerical simulation with different turbulent intensities, respectively. Considering the characteristics of scavenging process in the turbulent flow field, the concept of collision probability was introduced to give the new definition of computational formula of collection efficiency in turbulent flow condition. On the basis of the concept, the influence of turbulent effect on the collection efficiency was discussed and analyzed. The numerical results show that the collision probability distributions for particle of different diameters are affected by background turbulent intensity. Furthermore, collection efficiency is strongly affected by turbulent intensity for dp≤ 5μm. The influence of turbulent effect on collection efficiency depends on the density and diameter of particle, which was independent of the raindrop diameter.

aerosol particle；raindrop；turbulent intensity；collection efficiency

X513,P426.6

A

1000-6923(2017)01-0013-08

華鳳皎(1989-),女,山東煙臺(tái)人,東華大學(xué)博士研究生,主要從事城市大氣環(huán)境與室內(nèi)空氣品質(zhì)研究.發(fā)表論文3篇.

2016-03-10

國(guó)家自然科學(xué)基金資助項(xiàng)目(41275157);上海市教委科研創(chuàng)新重點(diǎn)項(xiàng)目(14ZZ073)

* 責(zé)任作者, 教授, ymkang@dhu.edu.cn