甘陽(yáng)陽(yáng)，李志生

（廣東工業(yè)大學(xué) 土木與交通工程學(xué)院，廣東 廣州 510006）

近年來(lái)，室內(nèi)空氣品質(zhì)問題備受關(guān)注。PM2.5是室內(nèi)空氣的主要污染源，可進(jìn)入人體呼吸系統(tǒng)，對(duì)人體健康造成很大威脅。廚房烹飪是PM2.5的重要來(lái)源，高溫烹調(diào)油和食物能釋放各種有害的顆粒物和有機(jī)化合物，包括揮發(fā)性有機(jī)物、油脂、烴類、醛酮類和多環(huán)芳香烴[1-2]等。國(guó)內(nèi)外對(duì)于廚房顆粒物的研究主要集中在了排放特征[3-5]、油煙散發(fā)實(shí)測(cè)[6-7]、濃度影響因素分析[8-10]等方面，關(guān)于廚房?jī)?nèi)氣流組織、溫度和PM2.5濃度分布的數(shù)值模擬研究不足。深入研究廚房PM2.5的運(yùn)輸、分布特征及其控制策略，可以有效防治廚房PM2.5污染，減少人員暴露時(shí)間。

盡管現(xiàn)有的抽油煙機(jī)已經(jīng)可以控制大部分烹飪煙氣，但由于烹飪油煙散發(fā)快速且劇烈，廚房?jī)?nèi)依然存在空氣質(zhì)量不良和熱舒適性差的現(xiàn)象[11]。利用空氣幕射流[12]作為動(dòng)力進(jìn)行補(bǔ)風(fēng)，可將PM2.5輸送到抽油煙機(jī)排風(fēng)口，幫助其排除；同時(shí)，利用射流產(chǎn)生的風(fēng)幕阻隔、控制PM2.5的擴(kuò)散。Simone等[13]研究表明，預(yù)測(cè)平均投票數(shù)/不滿意者百分?jǐn)?shù)(predicted mean vote/predicted percentage of dissatisfied, PMV/PPD)指數(shù)并不適用于廚房熱環(huán)境。Zhang等[14]研究了空氣幕的風(fēng)量和噴射角對(duì)煙氣控制的影響。Zhou等[15]結(jié)合實(shí)驗(yàn)和數(shù)值模擬，研究了空氣幕對(duì)CO2的控制效果，結(jié)果表明：空氣幕射流速度在0～1 m/s，4個(gè)射流角都為90°時(shí)，控制效果最佳。

本文研究對(duì)象為典型民居廚房，面積適中，窗戶在烹飪區(qū)側(cè)后方，區(qū)別于前人研究較多的側(cè)面通風(fēng)?；谟?jì)算流體動(dòng)力學(xué)(Computational fluid dynamics，CFD)的數(shù)值模擬在預(yù)測(cè)空氣運(yùn)動(dòng)和顆粒物模擬方面較實(shí)驗(yàn)研究具有優(yōu)勢(shì)，并在流場(chǎng)分布及其優(yōu)化控制研究上得到廣泛應(yīng)用[16-17]。因此，為對(duì)廚房PM2.5污染物和余熱進(jìn)行有效控制，采用空氣幕送風(fēng)方式，使用Fluent 16.0軟件對(duì)其控制效果進(jìn)行模擬分析，以期改善廚房?jī)?nèi)的空氣質(zhì)量和熱舒適性。通過氣流組織分析揭示了空氣幕的作用原理與效果，量化分析了3種空氣幕射流速度的降溫幅度和PM2.5排除率。研究結(jié)論可對(duì)廚房PM2.5控制策略提供參考，為空氣幕送風(fēng)系統(tǒng)的研究提供模擬數(shù)據(jù)和理論依據(jù)。

1 研究方法

1.1 物理模型

民居廚房尺寸為2.0 m×2.2 m×2.4 m(長(zhǎng)×寬×高)，內(nèi)部布局為L(zhǎng)型，廚房門、窗尺寸分別為1.96 m×0.69 m×0.1 m和1.27 m×0.75 m×0.15 m。廚房物理模型在實(shí)際測(cè)量的基礎(chǔ)上參考相關(guān)文獻(xiàn)[15,18-19]進(jìn)行了適當(dāng)簡(jiǎn)化，結(jié)果如圖1所示。廚房?jī)?nèi)有大氣嵌入式雙眼灶，灶簡(jiǎn)化為兩個(gè)尺寸為0.24 m×0.22 m×0.07 m的污染源。抽油煙機(jī)為側(cè)吸式，雙風(fēng)機(jī)結(jié)構(gòu)，擁有兩檔標(biāo)稱排風(fēng)量15 m3/min和17 m3/min。排風(fēng)口直徑為0.22 m，側(cè)面布置；油煙分離板、油槽等均被簡(jiǎn)化。空氣幕放置在烹飪?cè)粗車?，位于抽油煙機(jī)正下方，如圖2所示?？諝饽挥?個(gè)長(zhǎng)槽和2個(gè)短槽組成，長(zhǎng)度分別為0.84 m和0.39 m，槽寬均為0.03 m。人體高度168 cm，高度劃分遵循以下比例：身體頭部部分18%，上半身40%，下半身42%；與操作臺(tái)外圍的水平距離為0.2 m。

1.2 網(wǎng)格劃分

圖 1 廚房灶臺(tái)實(shí)物及簡(jiǎn)化物理模型Fig.1 Physical diagram of kitchen’s hearth & simplified physical model of kitchen

圖 2 空氣幕安裝示意圖Fig.2 Air curtain installation diagram

網(wǎng)格劃分是數(shù)值模擬的前提和重要組成部分，劃分質(zhì)量將直接決定計(jì)算的精度、時(shí)間和收斂性。本文采用靈活性較高的非結(jié)構(gòu)化網(wǎng)格，利用ICEM CFD軟件進(jìn)行劃分。為確保模擬結(jié)果的合理和準(zhǔn)確，對(duì)網(wǎng)格質(zhì)量和數(shù)量進(jìn)行了管理和控制，具體如下：

(1) 模型表面和內(nèi)部空間分別使用三角形和四面體網(wǎng)格單元。(2) 對(duì)排風(fēng)口、空氣幕、門和窗附近速度梯度較大的區(qū)域進(jìn)行網(wǎng)格局部加密。(3) 對(duì)顆粒物散發(fā)源即爐灶及附近區(qū)域、溫度梯度較大區(qū)域和人體進(jìn)行網(wǎng)格局部加密。(4) 通過不斷調(diào)整設(shè)置因子和最大允許數(shù)值排除網(wǎng)格質(zhì)量不好的單元。(5) 控制網(wǎng)格參數(shù)縱橫比和扭曲率等處在良好或可接受范圍內(nèi)。

此外，為確保網(wǎng)格數(shù)量差異對(duì)計(jì)算結(jié)果沒有影響，對(duì)廚房模型進(jìn)行了網(wǎng)格獨(dú)立性檢驗(yàn)。從計(jì)算精度、時(shí)間和對(duì)計(jì)算器的要求3個(gè)方面綜合考慮，在數(shù)值模擬中選擇了427 896網(wǎng)格數(shù)量的方案。網(wǎng)格劃分情況如圖3所示。

圖 3 廚房網(wǎng)格劃分示意圖Fig.3 Kitchen interior meshing chart

1.3 數(shù)學(xué)模型

引入雷諾平均模擬法，將瞬時(shí)Navier-Stokes方程對(duì)時(shí)間進(jìn)行統(tǒng)計(jì)平均，連續(xù)相模型選用常見的RNG k-ε，組建包括連續(xù)性方程、能量方程、動(dòng)量方程、k方程和ε方程在內(nèi)的一整套封閉方程組。近壁面區(qū)的流動(dòng)計(jì)算采用壁面函數(shù)法。顆粒相模型選擇離散相模型(DPM)。DPM基于歐拉-拉格朗日方法，對(duì)大量顆粒軌道進(jìn)行統(tǒng)計(jì)分析，計(jì)算得到顆粒群的運(yùn)動(dòng)特征。但由于逐一計(jì)算單個(gè)顆粒的運(yùn)動(dòng)軌跡，計(jì)算量較大，要求顆粒相的體積分?jǐn)?shù)小于10%。顆粒相在X方向的運(yùn)動(dòng)控制方程如式(1)所示。

式中， up、u 分 別為顆粒相、空氣的速度，m/s；t為時(shí)間，s； FD為 斯托克斯阻力系數(shù)，即弛豫時(shí)間，s-1； ρp、ρ 分別為顆粒相、空氣的密度，kg/m3； gx為重力加速度，m/s2； Fx為顆粒相所受其他單位質(zhì)量力，m/s2。

1.4 工況及邊界條件

模擬工況共設(shè)置4組，工況1、2、3的空氣幕射流速度分別為0.3 m/s、0.6 m/s、0.9 m/s，工況4作為沒有使用空氣幕的對(duì)照組。由于單、雙灶的差異主要體現(xiàn)在PM2.5的散發(fā)數(shù)量級(jí)上，對(duì)氣流組織影響相對(duì)不大，故為減少變量，4種工況啟用的均是左邊的爐灶和抽油煙機(jī)。邊界條件設(shè)置參考了相關(guān)文獻(xiàn)[15,18-21]。廚房維護(hù)結(jié)構(gòu)為固體壁面邊界，關(guān)閉房門，傳熱系數(shù)均為1.5 w/(m2·k)。窗戶保持開啟，為壓力出口邊界。假設(shè)人體站立并做少量機(jī)械功，設(shè)為恒溫31 ℃的壁面邊界。操作臺(tái)定義為絕熱邊界。考慮油煙分離板等的影響，抽油煙機(jī)排風(fēng)口設(shè)為速度入口邊界，方向?yàn)樨?fù)值，風(fēng)速為2.5 m/s。將爐灶定義為污染源和熱源，設(shè)為速度入口邊界，速度為0.6 m/s，溫度為50 ℃?？諝饽煌瑯釉O(shè)為速度入口邊界，射流角度為90°，溫度為25 ℃。

實(shí)際灶眼以圓形居多，但由于烹飪期間食物的PM2.5散發(fā)量遠(yuǎn)大于灶眼，故可以忽略灶眼形狀的影響；同時(shí)為與現(xiàn)有文獻(xiàn)所使用的方形散發(fā)源作對(duì)照，并假設(shè)使用平底鍋，故將散發(fā)源簡(jiǎn)化為面源[18-19,21]。不考慮室外顆粒物通過門、窗或維護(hù)結(jié)構(gòu)的滲透效應(yīng)。顆粒物粒徑設(shè)為2.5 μm，材料為默認(rèn)無(wú)煙煤細(xì)顆粒，密度為950 kg/m3，散發(fā)量為3×10-8kg/s[21-22]，接近煎炒烹飪方式。顆粒相與連續(xù)相的耦合采用雙向耦合方式。連續(xù)相迭代10次后，對(duì)顆粒物進(jìn)行一次軌跡計(jì)算，采用穩(wěn)態(tài)追蹤方式。速度出口和壓力出口的DPM邊界條件設(shè)為逃逸，其余設(shè)置為反射。

2 結(jié)果與討論

2.1 對(duì)廚房?jī)?nèi)氣流組織的影響

圖4是以空氣幕作為流線出發(fā)點(diǎn)的室內(nèi)流線云圖，可以看出：工況1、2、3情況下，空氣幕射流氣流對(duì)烹飪區(qū)域產(chǎn)生了很好的包裹效應(yīng)。氣流垂直均勻向上運(yùn)動(dòng)，一部分氣流被排風(fēng)口所吸入，另一部分則繼續(xù)向上運(yùn)動(dòng)到達(dá)天花，在廚房?jī)?nèi)產(chǎn)生繞流或經(jīng)窗戶流出。工況2、3的流線相較工況1更均勻，但有不同程度的繞流。因此可見，空氣幕提供了均勻、穩(wěn)定的流場(chǎng)；良好的氣流組織可以使PM2.5和熱流得到更好的控制。主要原因是：(1)空氣幕射流可以使室內(nèi)外壓差減少，降低窗戶進(jìn)氣量，使烹飪區(qū)內(nèi)部獲得穩(wěn)定低壓區(qū)域，有利于抽油煙機(jī)排出更多油煙。(2)由于PM2.5受空氣分布影響很大，所以空氣幕射流能夠有效阻隔PM2.5的擴(kuò)散。

2.2 對(duì)廚房?jī)?nèi)溫度分布的影響

工況1～4下廚房?jī)?nèi)的溫度分布見圖5，可以看出：與工況4相比，使用了空氣幕的工況1、2、3的溫度明顯更低。溫度分布以操作臺(tái)為界線大致分為兩部分。一部分是操作臺(tái)上部，煙氣主體在烹飪區(qū)域內(nèi)部能達(dá)到50 ℃，隨著空氣幕射流速度的增加，煙氣上部截面逐漸變窄；在烹飪區(qū)外圍，工況1、2、3的溫度約在28～35 ℃，而工況4在36 ℃左右。另一部分是操作臺(tái)外部，人體周圍的空氣溫度較低，約在25～29 ℃；呼吸區(qū)域的溫度約在29 ℃，小于工況4的31 ℃。此外，工況3相比工況2溫度更高，可能的原因是空氣幕射流速度偏高，排風(fēng)口難以將所有氣流排出，致使烹飪煙氣外泄，溫度升高。因此可見，空氣幕很好地控制了高溫氣體的蔓延，改善了人體呼吸區(qū)域的熱舒適性。

圖 4 室內(nèi)流線云圖Fig.4 Interior streamline chart

2.3 對(duì)廚房?jī)?nèi)PM2.5濃度分布的影響

工況1～4下廚房?jī)?nèi)PM2.5濃度分布見圖6?？梢钥闯觯弘m然抽油煙機(jī)和空氣幕有效控制了大部分PM2.5，但依然有部分PM2.5在散發(fā)過程中會(huì)向外逃逸，進(jìn)而向周圍擴(kuò)散、懸浮或沉降。工況1、4的PM2.5分別在天花、地板和氣流渦旋區(qū)有積聚。將工況1、2、3與工況4對(duì)比，可以看出空氣幕顯著降低了廚房?jī)?nèi)的PM2.5濃度，射流氣流也有效阻隔了PM2.5向外的擴(kuò)散。在烹飪區(qū)上部，工況4的PM2.5濃度可達(dá)到5×10-7kg/m3，工況1、2、3的PM2.5濃度約在1×10-7～2.5×10-7kg/m3；而在操作臺(tái)外部，工況4的PM2.5濃度約在1×10-7～4×10-7kg/m3，工況2、3的PM2.5濃度基本在1.5×10-7kg/m3以下，說(shuō)明使用空氣幕后，PM2.5濃度明顯降低。

圖 5 廚房?jī)?nèi)的溫度體繪制圖Fig.5 Temperature volume rendering in kitchen

圖 6 廚房?jī)?nèi)的PM2.5質(zhì)量濃度體繪制圖Fig.6 PM2.5 mass concentration volume rendering in kitchen

此外，工況2、3相比工況1對(duì)PM2.5的擴(kuò)散有更好的阻隔，特別是在抽油煙機(jī)上部和人體前方。雖然工況3在天花附近較工況2有更嚴(yán)重的積聚，但依然可以看出：更高的空氣幕射流速度導(dǎo)致了更低的PM2.5濃度；同時(shí)，有效降低了烹飪?nèi)藛T呼吸區(qū)的暴露量。

工況1～4下廚房?jī)?nèi)平均溫度、PM2.5平均質(zhì)量濃度及其控制效果見表1。從表1可知：在0.3 m/s、0.6 m/s和0.9 m/s 三種射流速度下，三者的降溫幅度分別為1.07、2.28和1.64 ℃；PM2.5的排除幅度分別為7.84×10-8kg/m3、1.33×10-7kg/m3和1.34×10-7kg/m3。與工況4相比，使用了空氣幕的工況1、2、3的廚房平均降溫幅度在1～2 ℃，PM2.5濃度排除率提高了44%～75%。以上結(jié)果說(shuō)明空氣幕的使用有效改善了廚房?jī)?nèi)的熱環(huán)境和空氣質(zhì)量。

表 1 不同工況下廚房?jī)?nèi)的平均溫度和PM2.5平均質(zhì)量濃度Table 1 Average temperature & PM2.5 concentration in the kitchen under different cases

從表1還可以得知，當(dāng)空氣幕射流速度為0.3 m/s時(shí)，PM2.5的排除率和降溫效果最差。但是，在0.6 m/s和0.9 m/s的情況下，PM2.5的排除率已相差無(wú)幾；而0.6 m/s情況下的降溫效果還要優(yōu)于0.9 m/s。這是因?yàn)檫^高的射流速度會(huì)對(duì)烹飪煙氣的運(yùn)動(dòng)產(chǎn)生擾動(dòng)，可能導(dǎo)致抽油煙機(jī)難以將所有氣流排出，從而導(dǎo)致PM2.5的排除率降低。顯然，空氣幕射流速度的選取并不是越大越好，存在控制效果最佳的射流速度。

此外，當(dāng)送風(fēng)速度大于0.3 m/s時(shí)還會(huì)對(duì)人造成吹風(fēng)感[23]。在0.3 m/s、0.6 m/s和0.9 m/s 這3種風(fēng)速條件下，風(fēng)速越小則吹風(fēng)感越低；而從能耗角度考慮，送風(fēng)速度越高，意味著空氣幕的能耗也越高。為降低廚房?jī)?nèi)污染物濃度，應(yīng)同時(shí)提高抽油煙機(jī)的排風(fēng)量，僅僅依靠提高空氣幕射流速度是不合理的。綜上所述，本文認(rèn)為空氣幕射流速度的最佳值為0.6 m/s是比較合理的。

3 結(jié)論

本文提出運(yùn)用空氣幕與抽油煙機(jī)結(jié)合的方式對(duì)廚房PM2.5進(jìn)行有效控制。模擬得出了廚房?jī)?nèi)的氣流組織、溫度場(chǎng)和PM2.5濃度場(chǎng)。具體結(jié)論如下：(1) 氣流組織對(duì)PM2.5分布影響很大，空氣幕射流氣流對(duì)烹飪區(qū)域產(chǎn)生了很好的包裹效應(yīng)，提供了更加均勻、穩(wěn)定的流場(chǎng)，可以阻隔PM2.5的擴(kuò)散；同時(shí)減少窗戶的進(jìn)氣量，有利于抽油煙機(jī)工作。(2) 使用了空氣幕的工況1、2、3與沒有使用空氣幕相比，無(wú)論是在溫度還是PM2.5的控制上都有顯著效果，可以改善廚房?jī)?nèi)的熱環(huán)境和空氣質(zhì)量。廚房平均降溫幅度分別為1.08 ℃、2.28 ℃、1.64 ℃，PM2.5排除率分別為44%、74%、75%。(3) 從降溫、控制效果、能耗和吹風(fēng)感綜合考慮，選擇0.6 m/s作為空氣幕的最佳射流速度。