焊接工藝多源散發(fā)的高溫顆粒污染物擴(kuò)散特性數(shù)值模擬

韓 坤,莊加瑋,2,刁永發(fā)*,任 美,張儷安

焊接工藝多源散發(fā)的高溫顆粒污染物擴(kuò)散特性數(shù)值模擬

韓 坤1,莊加瑋1,2,刁永發(fā)1*,任 美1,張儷安1

(1.東華大學(xué)環(huán)境科學(xué)與工程學(xué)院,上海 201620；2.常州大學(xué)石油工程學(xué)院,江蘇 常州 213164)

為探究熱工藝過程產(chǎn)生的熱氣流作用下顆粒的遷移規(guī)律,基于氣固兩相離散粒子模型(DPM)對(duì)多源浮射流伴生的高溫顆粒的擴(kuò)散特性進(jìn)行了數(shù)值研究,討論了兩相流運(yùn)動(dòng)過程中熱氣流與顆粒群的溫度和速度的瞬時(shí)變化情況.結(jié)果表明,對(duì)于473K￡0￡673K的高溫顆粒群,溫度衰減趨勢類似,顆粒群溫度分布呈中心對(duì)稱;當(dāng)5μm￡p￡20μm時(shí),顆粒與氣流之間的跟隨性隨著粒徑的增大而降低;熱羽流隨時(shí)間經(jīng)歷了獨(dú)立發(fā)展和相互合并的過程,羽流之間渦旋結(jié)構(gòu)的消失導(dǎo)致中心位污染源散發(fā)的p=10μm的顆粒更容易發(fā)生沉降.

高溫顆粒；空氣質(zhì)量；熱工藝；擴(kuò)散特性；熱羽流

鋼鐵、鑄造等行業(yè)存在大量的高溫?zé)峁に囘^程,例如澆注、補(bǔ)焊、切割、高溫物料傾倒等工藝過程會(huì)產(chǎn)生初始溫度高于附近環(huán)境溫度的高溫浮力射流[1-2].這些高溫浮射流通常攜帶較多的顆粒污染物,顆粒以一定的初速度散發(fā)到周圍環(huán)境中,隨著浮射流的運(yùn)動(dòng)與發(fā)展,顆粒不斷與周圍空氣摻混擴(kuò)散,這些高溫顆粒一旦被人體吸入,可能會(huì)對(duì)人體產(chǎn)生較為嚴(yán)重的傷害[3-5].

目前高溫顆粒在室內(nèi)流場中的擴(kuò)散規(guī)律尚未明確,只有基于具體的工藝,探討其擴(kuò)散過程與周圍環(huán)境參數(shù)的相互關(guān)系,才能為制定相應(yīng)的控制方法提供科學(xué)有效的依據(jù)[6-8].國內(nèi)外學(xué)者近年來針對(duì)工業(yè)熱工藝過程產(chǎn)生的高溫顆粒的流動(dòng)特性進(jìn)行了相關(guān)研究,包括熱源附近顆粒群的微觀形成機(jī)理,以及其所形成的煙羽流在空間內(nèi)的擴(kuò)散運(yùn)動(dòng)[9-11],在揭示高溫顆粒運(yùn)動(dòng)過程的動(dòng)力學(xué)特性等方面取得了一定的研究進(jìn)展.Li等[12]通過數(shù)值研究發(fā)現(xiàn)單個(gè)污染源散發(fā)的高溫亞微米氣溶膠粒子運(yùn)動(dòng)過程主要受空氣阻力、布朗力、薩夫曼升力和熱電泳力的影響.Zhuang等[13]的結(jié)果表明瞬時(shí)工藝局部污染源散發(fā)的熱量愈多,顆粒遷移過程獲得的最大平均速度越高.

現(xiàn)有工業(yè)建筑中關(guān)于熱工藝產(chǎn)生的熱氣流作用下顆粒群遷移規(guī)律的研究,大多局限于單一局部污染源,而同一工藝過程往往存在多個(gè)獨(dú)立的污染源,已有部分學(xué)者的研究表明多污染源對(duì)應(yīng)單源時(shí)的流場特征存在顯著差異[14-16]. Lai等[17]發(fā)現(xiàn)浮射流的合并會(huì)導(dǎo)致與單個(gè)浮射流相關(guān)的被動(dòng)標(biāo)量場重疊,導(dǎo)致混合后的特性與獨(dú)立自由浮射流截然不同.為此他們提出了一種多股浮射流的動(dòng)態(tài)相互作用的一般半解析模型,預(yù)測了理想環(huán)境中多個(gè)浮射流的動(dòng)態(tài)相互作用.多個(gè)熱源產(chǎn)生的羽流之間相互作用會(huì)減少單股羽流的輸運(yùn)空間,Yang等[18]為此研究了不同間距的兩個(gè)相等熱污染源產(chǎn)生的熱羽流的流動(dòng)特性,獲得了不同間距比的空氣速度和溫度分布. Huang等[19]考慮了浮射流間距、初始速度和溫度比等因素,給出了控制兩個(gè)污染源的有效通風(fēng)設(shè)計(jì)指導(dǎo).

雖然前期學(xué)者們對(duì)多個(gè)污染源產(chǎn)生的高溫浮射流的流動(dòng)特性開展了一定的工作,但是尚缺乏針對(duì)多源浮射流中伴生的顆粒污染物的擴(kuò)散分布研究,且現(xiàn)有研究往往忽略了熱氣流與顆粒的兩相耦合作用. 本文以焊接工藝這一典型熱工藝為例,基于氣固兩相離散粒子模型(DPM)對(duì)多源瞬時(shí)散發(fā)的高溫顆粒的遷移規(guī)律進(jìn)行研究,討論兩相流運(yùn)動(dòng)過程中熱氣流溫度和速度隨時(shí)間的變化,對(duì)初始溫度0、顆粒粒徑p等不同影響因素作用下高溫顆粒群隨時(shí)間的動(dòng)態(tài)分布進(jìn)行分析和評(píng)價(jià),為多工位作業(yè)環(huán)境中顆粒污染物的控制提供一定的理論依據(jù).

1 模擬方法

1.1 物理模型

圖1 物理模型

研究建立的物理模型如圖1所示,該模型為實(shí)際工業(yè)建筑的一部分.計(jì)算域尺寸長()×寬()×高()=12×10×8m3,三個(gè)污染源位于廠房底部,考慮作業(yè)符合實(shí)際安全規(guī)范的前提下,每兩個(gè)污染源之間設(shè)置間隔為1m,污染源距地面高度0.8m,顆粒群與氣流的混合物由直徑為0.4m的圓形出口釋放. 其中顆粒為車間常見的鐵基顆粒,顆粒密度ρ= 7800kg/m3,比熱容C*=460J/(kg·K),導(dǎo)熱系數(shù)λ=80W/(m·K).

研究過程涉及到的工況如表1所示:

表1 工況設(shè)置

1.2 氣固耦合兩相流輸運(yùn)模型

1.2.1 氣流流動(dòng)模型 在解決房間內(nèi)部流體流動(dòng)的問題方面,目前常用的數(shù)學(xué)模型包括標(biāo)準(zhǔn)模型和Realizable～湍流模型等[20-21].由于Realizable～湍流模型對(duì)標(biāo)準(zhǔn)模型中的湍流粘性和湍流耗散率進(jìn)行了修正,對(duì)于模擬射流和混合流的自由流動(dòng)具有更高的計(jì)算精度和更好的穩(wěn)定性,因此研究過程中選取了Realizable～湍流模型.其中對(duì)于非穩(wěn)態(tài)流動(dòng)問題,其輸運(yùn)方程表達(dá)式可以表示為以下形式:

湍流動(dòng)能:

式中:為時(shí)間;為流體密度;σ是湍流動(dòng)能對(duì)應(yīng)的Prandtl數(shù);為湍流耗散率;G和G分別是因平均速度梯度與熱浮力而產(chǎn)生的湍流動(dòng)能;M表示不可壓縮湍流中因脈動(dòng)膨脹而產(chǎn)生的全部耗散率;為湍流粘性系數(shù);S為用戶定義源項(xiàng).

湍流耗散率:

式中:1ε、2ε和3ε是經(jīng)驗(yàn)常數(shù);為表面張力系數(shù);σ為耗散率對(duì)應(yīng)的Prandtl數(shù);是湍流時(shí)間尺度與平均速度的比值.

1.2.2 氣溶膠傳輸模型 依據(jù)牛頓第二運(yùn)動(dòng)定律,單個(gè)顆粒運(yùn)動(dòng)過程的受力平衡方程為:

其中曳力系數(shù)D與顆粒的雷諾數(shù)p有關(guān),其大小可通過式(5)得出[22-23]:

采用相間耦合的隨機(jī)游走模型(Discrete Radom Walk,DRW)預(yù)測因?yàn)榱黧w湍流引起的顆粒擴(kuò)散行為.該模型基于隨機(jī)格式通過連續(xù)相的瞬時(shí)速度(u)來預(yù)測顆粒的流動(dòng)軌跡:

式中:ζ是服從正態(tài)分布的隨機(jī)數(shù).

顆粒群的流動(dòng)相態(tài)采用顆粒物的體積分?jǐn)?shù)進(jìn)行表征:

式中:N為顆粒數(shù)量;V為單顆粒體積;為氣固兩相混合物體積.研究設(shè)置的顆粒質(zhì)量流量為0.00005kg/ s,經(jīng)計(jì)算,任意工況下顆粒體積百分?jǐn)?shù)均小于10-6.根據(jù)文獻(xiàn)[24]提到的兩相流劃分形式,視作稀疏兩相流,忽略顆粒間的相互作用.

考慮到顆粒與氣流之間存在動(dòng)量和熱量的傳遞,因此選擇雙向耦合方法(Two-Way Coupling Method)求解各個(gè)時(shí)刻離散相和連續(xù)相之間的動(dòng)量和熱量交換[25-26].兩相流間的動(dòng)量交換量為:

式中:p為顆粒的質(zhì)量流量,kg/s;?為時(shí)間步長,s.

熱量交換可以表示為:

式中:pin和pout分別代表顆粒發(fā)生熱交換前后的溫度,K;r為參考溫度,K.

2 模擬結(jié)果與分析

2.1 模型驗(yàn)證

2.1.1 熱空氣流動(dòng)模型可靠性驗(yàn)證 Huang等[27]的煙霧可視化實(shí)驗(yàn)結(jié)果以及Zhuang等[13]的速度流線輪廓模擬結(jié)果被用于驗(yàn)證本文中熱空氣流動(dòng)數(shù)值方法的合理性. Huang等[27]設(shè)置實(shí)驗(yàn)所在的房間尺寸為10.30m×6.30m×3.50m,直徑為0.38m的浮射流噴口布置在距地面0.80m位置處,測試時(shí),環(huán)境壓力為101325Pa,浮射流初始溫度為303K(工況一);Zhuang等[13]建立了直徑=0.4m的面散發(fā)源數(shù)值模型,面污染源沿豎直向上方向朝周圍環(huán)境直接釋放熱空氣和固態(tài)金屬顆粒.

在0=303K的工況一條件下,對(duì)單源散發(fā)污染物過程中熱氣流的流動(dòng)進(jìn)行了數(shù)值模擬,文獻(xiàn)數(shù)據(jù)與本文得到的模擬結(jié)果對(duì)比如圖2所示.

可以看出,氣流運(yùn)動(dòng)過程中在其邊緣存在明顯的渦旋結(jié)構(gòu).對(duì)比了實(shí)驗(yàn)及模擬的氣流演化進(jìn)程,本文數(shù)值模擬得到的渦旋結(jié)構(gòu)的演化過程與文獻(xiàn)[13,27]的結(jié)果吻合度較高. 這說明,本研究中使用的Realizable-湍流模型可以準(zhǔn)確地預(yù)測熱氣流的流動(dòng).

2.1.2 高溫顆粒輸運(yùn)模型的驗(yàn)證 為了進(jìn)一步驗(yàn)證本研究中采用的顆粒輸運(yùn)模型的合理性,對(duì)初始溫度373K工況二下工業(yè)廠房內(nèi)單污染源散發(fā)的粒徑為10μm的顆粒污染物的運(yùn)動(dòng)過程進(jìn)行了數(shù)值模擬,并將數(shù)值計(jì)算結(jié)果與Duan等[28]的數(shù)據(jù)進(jìn)行了比較.

圖2 速度流線軌跡隨時(shí)間的變化

圖3 顆粒群平均溫度隨時(shí)間的變化

圖4 顆粒群平均速度隨時(shí)間的變化

由圖3可見,對(duì)于初始溫度為373K的顆粒群,顆粒群在較短時(shí)間內(nèi)與周圍環(huán)境產(chǎn)生了明顯的熱量和動(dòng)量交換,=3.0s時(shí)溫度已接近常溫,隨著熱氣流速度的變化,顆粒群運(yùn)動(dòng)過程受重力和流體曳力的共同影響,速度呈現(xiàn)先增大而后減小的趨勢.這些結(jié)果與Duan等[28]的數(shù)據(jù)基本吻合,說明所采用的數(shù)值模擬方法可以準(zhǔn)確預(yù)測瞬時(shí)熱工藝過程中的氣固兩相流運(yùn)動(dòng)過程,其可用于進(jìn)一步研究.

2.2 多污染源氣流場流動(dòng)特性

在工況三0=473K,0=0.1s,0=0.3m/s且p= 10μm條件下,分析了局部熱氣流運(yùn)動(dòng)過程中速度場與溫度場的動(dòng)態(tài)變化情況.

多源時(shí)氣流速度場與單源時(shí)最本質(zhì)的區(qū)別就是各污染源產(chǎn)生的羽流會(huì)隨時(shí)間的發(fā)展逐漸靠近.如圖5所示,多股羽流在上升過程中首先經(jīng)歷獨(dú)立的向上部空間發(fā)展,熱羽流在浮升力的作用下向上運(yùn)動(dòng)不斷卷吸周圍空氣,各個(gè)羽流運(yùn)動(dòng)過程中在其邊緣存在明顯的渦旋結(jié)構(gòu).當(dāng)>20s后,觀察到羽流與羽流之間的渦旋結(jié)構(gòu)消失,速度場融合,羽流在發(fā)展過程中相互合并,這是由于羽流之間動(dòng)量通量的不平衡,附近的羽流會(huì)對(duì)這種卷吸起一定程度上的抑制作用,即產(chǎn)生了康達(dá)效應(yīng)[29],使得相鄰的羽流逐漸靠近.

圖6反映了工況三條件下氣流溫度場隨時(shí)間的變化情況,熱氣流沿程不斷卷吸周圍空氣,使羽流質(zhì)量增加,而傳遞的熱量守恒,熱羽流不斷與卷吸而來的常溫氣流進(jìn)行能量交換,熱氣流在短時(shí)間內(nèi)與周圍環(huán)境交換了大量的熱量,整體溫度急劇下降,=5s時(shí)氣流核心溫度已下降135K左右.伴隨著羽流之間的合并和卷吸作用的減弱,氣流場溫度衰減速度變緩.

圖5 氣流速度場隨時(shí)間的變化

圖6 氣流溫度場隨時(shí)間的變化

2.3 多源散發(fā)高溫顆粒溫度演變特征

工況三條件下熱氣流運(yùn)動(dòng)過程中高溫顆粒的溫度變化如圖7所示.初始階段,高溫顆粒與周圍環(huán)境存在較大溫差,受熱交換的作用,高溫顆粒的溫度隨著時(shí)間的推移迅速減小,=10s時(shí)顆粒群已接近環(huán)境溫度.顆粒群外圍的顆粒溫度衰減較快,而中心區(qū)域的顆?？梢韵鄬?duì)保持更高的溫度,一方面是因?yàn)橥鈱宇w粒與周圍空氣的接觸面積相對(duì)較大,更利于熱量傳遞;另一方面是因?yàn)橥鈬w粒更容易與沿渦旋結(jié)構(gòu)的回流區(qū)卷吸而來的常溫氣流換熱.結(jié)合圖5分析,受空氣卷吸作用的影響,隨著熱氣流的運(yùn)動(dòng),兩側(cè)的顆粒群逐漸往中間區(qū)域逼近.由于羽流在發(fā)展過程中發(fā)生了相互合并,外圍兩側(cè)的顆粒群沿渦旋結(jié)構(gòu)的回流區(qū)域繼續(xù)向上往中間區(qū)域偏轉(zhuǎn)運(yùn)動(dòng),=150s時(shí),原本處在中間區(qū)域的顆粒群因?yàn)閮蓚?cè)渦旋結(jié)構(gòu)的消失,在重力的作用下產(chǎn)生了明顯的沉降.

圖7 顆粒群的溫度隨時(shí)間的變化

以上結(jié)果表明,顆粒的溫度變化與顆粒運(yùn)動(dòng)過程中同熱氣流的跟隨性緊密相聯(lián).熱工藝過程前期,各污染源散發(fā)的顆粒表現(xiàn)出的溫度相對(duì)較高,并且越接近熱氣流與顆粒群混合物核心區(qū)域兩相流流速及溫度越高,人體暴露風(fēng)險(xiǎn)也越大,需要預(yù)防其對(duì)作業(yè)人員的燒傷.對(duì)于多個(gè)工位共同作業(yè)的工藝過程后期,處在中間工位的作業(yè)人員更容易受到顆粒污染物的影響,需要采取措施加強(qiáng)這一區(qū)域的安全衛(wèi)生防護(hù).

以工況三(0=473K)、工況四(0=573K)、工況五(0=673K)為例,探討了初始溫度這一因素對(duì)于高溫顆粒進(jìn)入環(huán)境后不同區(qū)域顆粒的溫度分布情況,顆粒群不同區(qū)域在=3s時(shí)沿X方向的溫度分布如圖8所示.

從圖中可以看出,顆粒群初始溫度與外界環(huán)境溫差越大,顆粒群溫度衰減的幅度越高.不同初始溫度條件下的高溫顆粒的溫度衰減的變化趨勢類似,溫度分布呈中心對(duì)稱的趨勢,比較了當(dāng)=3s時(shí)顆粒群位于=1.5m和=2m區(qū)域的溫度分布,發(fā)現(xiàn)每個(gè)污染源對(duì)應(yīng)的顆粒群核心區(qū)域上方的溫度相對(duì)下方較高,這是因?yàn)槠鸪躅w粒上升過程中主要受周圍氣流的曳力,顆粒群下方的顆粒在擴(kuò)散過程中更容易與熱氣流沿程不斷卷吸來的常溫氣流換熱.

2.4 高溫顆粒擴(kuò)散過程的動(dòng)力學(xué)分析

為了明確污染源的位置分布對(duì)于顆粒群在空間內(nèi)遷移過程的影響,在初始條件0=473K,0=0.1s,0=0.3m/s,p=10μm的工況下,追蹤了不同位置污染源散發(fā)出的典型顆粒,分析了顆粒運(yùn)動(dòng)過程中瞬時(shí)速度及擴(kuò)散距離隨時(shí)間的變化.

如圖9所示,<2s時(shí),顆粒的速度逐漸增大,這一時(shí)段顆粒上升的動(dòng)力主要來源于熱氣流向上運(yùn)動(dòng)附加給顆粒的曳力,而后隨著周圍空氣的卷吸作用逐漸減弱,熱氣流對(duì)顆粒的曳力逐漸降低,5s<￡40s時(shí),顆粒群在慣性力的作用下繼續(xù)做減速上升運(yùn)動(dòng).>40s后,中間污染源散發(fā)的顆粒運(yùn)動(dòng)過程中主要受重力影響,兩側(cè)污染源散發(fā)的顆粒受外圍空氣卷吸作用的影響,在曳力和慣性力的作用下以較低的速度繼續(xù)向上運(yùn)動(dòng).

圖9 顆粒遷移過程速度隨時(shí)間的變化

圖10 顆粒垂直方向擴(kuò)散距離隨時(shí)間的變化

由圖10可見,在高溫顆粒擴(kuò)散初始階段,對(duì)于任一污染源散發(fā)的顆粒,顆粒在垂直方向的擴(kuò)散距離都會(huì)隨時(shí)間的推移逐漸增大;在20s<￡40s時(shí),周圍空氣的卷吸作用逐漸減弱,熱氣流對(duì)顆粒的曳力逐漸降低,兩側(cè)污染源散發(fā)的顆粒在慣性力的作用下繼續(xù)向上做減速運(yùn)動(dòng);結(jié)合圖9可以發(fā)現(xiàn),羽流與羽流在發(fā)展過程中相互合并,當(dāng)>40s后,渦旋結(jié)構(gòu)的消失使得曳力和慣性力難以將中間污染源散發(fā)的典型顆粒運(yùn)輸至Y=4m以上的位置,顆粒在重力作用下逐漸沉降.

顆粒徑向擴(kuò)散距離隨時(shí)間的變化如圖11所示,由于兩側(cè)污染源散發(fā)的顆粒更容易受周圍空氣卷吸作用的影響,顆粒隨著渦旋結(jié)構(gòu)的擴(kuò)散區(qū)沿水平方向運(yùn)輸,后因顆粒在慣性力和重力作用下在垂直方向做減速運(yùn)動(dòng)進(jìn)入渦旋結(jié)構(gòu)回流區(qū),兩側(cè)顆粒逐漸向中心軸區(qū)域靠攏.

圖11 顆粒徑向擴(kuò)散距離隨時(shí)間的變化

進(jìn)一步考察了粒徑這一因素對(duì)顆粒群遷移過程的影響,保持初始條件0= 473K,0= 0.1s,0= 0.3m/s不變,工況六(p=5μm)和工況七(p=20μm)條件下不同粒徑顆粒群運(yùn)動(dòng)過程中的速度分布情況如圖12所示.

結(jié)合圖5氣流速度場的變化可以發(fā)現(xiàn),顆粒和氣流之間的跟隨性會(huì)隨著粒徑的增大而逐漸變差,當(dāng)=20s時(shí),粒徑為20μm的顆粒群已有較多顆粒發(fā)生沉降.除此之外,受熱空氣速度自核心區(qū)域逐漸向外圍衰減的影響,位于顆粒群邊緣處的顆粒速度相對(duì)較低,脫離氣流運(yùn)動(dòng)概率較大.

圖12 不同粒徑顆粒群在t=20s時(shí)的速度分布

3 結(jié)論

3.1 多股熱羽流隨時(shí)間經(jīng)歷了獨(dú)立發(fā)展和相互合并的過程,羽流與羽流之間動(dòng)量通量的不平衡導(dǎo)致了相鄰羽流間渦旋結(jié)構(gòu)的消失,使得顆粒遷移過程受到的氣流曳力和慣性力難以將中間污染源散發(fā)的粒徑p=10μm的顆粒運(yùn)輸至高度4m以上的位置,這部分顆粒在重力的作用下更容易發(fā)生沉降.對(duì)于多個(gè)工位共同作業(yè)的熱工藝過程后期,處在中間工位的作業(yè)人員更容易受到顆粒污染物的影響,需要采取措施加強(qiáng)這一區(qū)域的安全衛(wèi)生防護(hù).

3.2 各個(gè)污染源散發(fā)的高溫顆粒遷移過程與顆粒同熱氣流之間的跟隨性緊密相聯(lián),高溫顆粒群在熱氣流的作用下會(huì)同時(shí)沿著垂直方向和徑向擴(kuò)散,5μm

3.3 兩側(cè)污染源散發(fā)的顆粒更容易受周圍空氣卷吸作用的影響,顆粒會(huì)隨著渦旋結(jié)構(gòu)的擴(kuò)散區(qū)沿水平方向運(yùn)輸,后在慣性力和重力作用下在垂直方向做減速運(yùn)動(dòng)進(jìn)入渦旋回流區(qū),兩側(cè)的顆粒逐漸向中心軸區(qū)域靠攏.

[1] Huang Y, Wang Y, Ren X, et al. Ventilation guidelines for controlling smoke, dust, droplets and waste heat: Four representative case studies in Chinese industrial buildings [J]. Energy and Buildings, 2016,128: 834-844.

[2] 魏夢舉,崔海航,張鴻雁,等.鐵水傾倒過程的熱羽流流場特性研究 [J]. 熱加工工藝, 2017,(17):114-118.

Wei M J, Cui H H, Zhang H Y, et al. Study on flow field of thermal plume during pouring process of molten iron [J]. Hot Working Technology, 2017,(17):114-118.

[3] 莊加瑋,刁永發(fā),楚明浩,等.可吸入顆粒物在不同阻塞性呼吸道內(nèi)運(yùn)動(dòng)與沉積特性 [J]. 中國環(huán)境科學(xué), 2021,41(7):3349-3359.

Zhuang J W, Diao Y F, Chu M H, et al. Transport and deposition characteristics of inhalable particulate matters in different obstructive airways [J]. China Environmental Science, 2021,41(7):3349-3359.

[4] Sharma M, Kataria K K, Suri N M, et al. Monitoring respirable dust exposure in fettling work environment of a foundry: A proposed design intervention [J]. Int. J. Saf. Secur. Eng, 2020,10:759-767.

[5] Rahul M, Sivapirakasam S P, Vishnu B R, et al. Health issue owing to exposure with welding fumes and their control strategies at the source-A review [J]. Materials Today: Proceedings, 2021,46:9239- 9245.

[6] Szczurek A, Maciejewska M, Wy?omańska A, et al. Normal and anomalous diffusion in fluctuations of dust concentration nearby emission source [J]. Physica A: Statistical Mechanics and its Applications, 2018,491:619-631.

[7] Wang H, Wang T, Liu L, et al. Numerical evaluation of the performances of the ventilation system in a blast furnace casthouse [J]. Environmental Science and Pollution Research, 2021:1-15.

[8] Bémer D, Régnier R, Subra I, et al. Ultrafine particles emitted by flame and electric arc guns for thermal spraying of metals [J]. Annals of occupational hygiene, 2010,54(6):607-614.

[9] Wang Y, Zhao T, Cao Z, et al. The influence of indoor thermal conditions on ventilation flow and pollutant dispersion in downstream industrial workshop [J]. Building and Environment, 2021,187:107400.

[10] Zhiyuan W, Bo L, Fali H, et al. Numerical simulation of smoke particle motion [C]//2017 International Conference on Advanced Mechatronic Systems (ICAMechS). IEEE, 2017:131-134.

[11] Wang Y, Duan M, Ren X, et al. Experimental study of dust emission: Comparison between high-temperature and ambient-temperature materials [J]. Powder Technology, 2016,301:1321-1329.

[12] Li C, Wang H, Yu C W, et al. Diffusion characteristics of the industrial submicron particle under Brownian motion and turbulent diffusion [J]. Indoor and Built Environment, 2021:1420326X21991055.

[13] Zhuang J, Diao Y, Shen H. Numerical investigation on transport characteristics of high-temperature fine particles generated in a transiently welding process [J]. International Journal of Heat and Mass Transfer, 2021,176:121471.

[14] Ichimiya K, Saiki H. Behavior of thermal plumes from two-heat sources in an enclosure [J]. International Journal of Heat and Mass Transfer, 2005,48(16):3461-3468.

[15] Tokuhiro A, Kimura N. An experimental investigation on thermal striping: Mixing phenomena of a vertical non-buoyant jet with two adjacent buoyant jets as measured by ultrasound Doppler velocimetry [J]. Nuclear Engineering and Design, 1999,188(1):49-73.

[16] Yannopoulos P C, Noutsopoulos G C. Interaction of vertical round turbulent buoyant jets—Part II: Superposition method [J]. Journal of Hydraulic Research, 2006,44(2):233-248.

[17] Lai A C H, Lee J H W. Dynamic interaction of multiple buoyant jets [J]. Journal of fluid mechanics, 2012,708:539-575.

[18] Yang C, Li A, Gao X, et al. Interaction of the thermal plumes generated from two heat sources of equal strength in a naturally ventilated space [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2020,198:104085.

[19] Huang Y, Lu K, Guo J, et al. Study on ventilation performance of lateral exhaust hood under the influence of two high-temperature buoyant jets [J]. Building and Environment, 2020,177:106849.

[20] Shaheed R, Mohammadian A, Gildeh H K. A comparison of standard–and realizable–turbulence models in curved and confluent channels [J]. Environmental Fluid Mechanics, 2019,19(2):543-568.

[21] Lateb M, Masson C, Stathopoulos T, et al. Comparison of various types of–models for pollutant emissions around a two-building configuration [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2013,115:9-21.

[22] Zhong K, Yang X, Kang Y. Effects of ventilation strategies and source locations on indoor particle deposition [J]. Building and environment, 2010,45(3):655-662.

[23] Limam K, El Hamdani S, Abadie M, et al. Numerical study of transport and particle deposition inside buildings [J]. Energy Procedia, 2017,139:430-436.

[24] Karadimou D P, Markatos N C. Modelling of two-phase, transient airflow and particles distribution in the indoor environment by Large Eddy Simulation [J]. Journal of Turbulence, 2016,17(2):216-236.

[25] Ireland P J, Desjardins O. Improving particle drag predictions in Euler–Lagrange simulations with two-way coupling [J]. Journal of Computational Physics, 2017,338:405-430.

[26] Yang Y, Wang Y, Song B, et al. Transport and control of droplets: A comparison between two types of local ventilation airflows [J]. Powder Technology, 2019,345:247-259.

[27] Huang Y, Wang W, Lu K, et al. Flow-field evolution and vortex structure characteristics of a high-temperature buoyant jet [J]. Building and Environment, 2021,187:107407.

[28] Duan M, Wang Y, Gao D, et al. Modeling dispersion mode of high- temperature particles transiently produced from industrial processes [J]. Building and Environment, 2017,126:457-470.

[29] Freire A P S, Miranda D D E, Luz L M S, et al. Bubble plumes and the Coanda effect [J]. International journal of multiphase flow, 2002, 28(8):1293-1310.

Numerical simulation of diffusion characteristics of high-temperature particles generated by multiple pollution sources in welding process.

HAN Kun1, ZHUANG Jia-wei1,2, DIAO Yong-fa1*, REN Mei1, ZHANG Li-an1

(1.College of Environmental Science and Engineering, Donghua University, Shanghai 201620, China；2.College of Petroleum Engineering, Changzhou University, Changzhou 213164, China)., 2022,42(6)：2552～2560

In order to explore the migration law of particles under the action of hot airflow in the process of hot working, the diffusion characteristics of high-temperature particles associated with multiple buoyant jets are numerically studied based on the gas-solid two-phase discrete particle model (DPM). The instantaneous changes of the temperature and velocity of the hot airflow and the particle group during the two-phase flow are discussed. The results show that for the high-temperature particle group with 473K￡0￡673K, the temperature attenuation trend is similar, and the temperature distribution of the particle group is centrosymmetric. When the particle size is between 5μm and 20μm, the followability between particles and airflow decreases as the particle size increases.Thermal plumes have undergone independent development and merging processes over time.The disappearance of the vortex structure between the plumes causes the particles with a particle size of 10mm emitted by the central pollution source to settle moreeasily.

high-temperature particles；air quality；heat process；diffusion characteristics；thermal plume

X513

A

1000-6923(2022)06-2552-09

韓 坤(1995-),男,江蘇泰州人,東華大學(xué)博士研究生,主要從事工業(yè)建筑室內(nèi)環(huán)境污染物控制研究.發(fā)表論文3篇.

2021-11-10

國家重點(diǎn)研發(fā)計(jì)劃(2018YFC0705300);中央高?；究蒲袠I(yè)務(wù)費(fèi)重點(diǎn)項(xiàng)目(2232017A-09)

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