王芳 王鵬浩 徐洪祥

摘?要：針對冬季東北地區(qū)房間通風性差及常規(guī)地板鋪設方式帶來的甲醛污染問題，使用C語言自定義編程，確定密閉條件下室內(nèi)地板散發(fā)甲醛的數(shù)學模型，利用Fluent求解及分析此條件下室內(nèi)甲醛濃度隨時間的變化規(guī)律，探討密閉條件下室內(nèi)甲醛濃度隨位置變化的擴散機理;研究自然通風時，室內(nèi)揮發(fā)性氣體的速度分布和甲醛的濃度分布。研究結(jié)果表明，采用自然通風方式可有效控制室內(nèi)甲醛污染，對改善室內(nèi)空氣品質(zhì)具有指導意義。

關(guān)鍵詞：自然通風;甲醛;數(shù)值模擬;濃度場

DOI：10.15938/j.jhust.2018.06.006

中圖分類號： TU834.1

文獻標志碼： A

文章編號： 1007-2683（2018）06-0029-06

Abstract：In order to solve the problem of formaldehyde pollution in northeastern winter，which caused by the poor ventilation and the conventional floor laying method?the mathematical model of indoor formaldehyde distribution was established by C language programming under enclosed condition.?By using Fluent under this condition?the changing rule and diffusion mechanism of indoor formaldehyde concentration was determined by analyzing its changes over time and distribution respectively.?The velocity distribution of volatile gas and the concentration distribution of formaldehyde were also studied under the natural ventilation mode.?The results showed that the concentration of formaldehyde could be kept at relatively low levels under the natural ventilation mode?and it is instructive to improve the indoor air quality.

Keywords：natural draft;formaldehyde;numerical simulation;concentration field

0?引?言

隨著人們生活水平和審美的提高，越來越多的人們會在室內(nèi)鋪設各種各樣的地板。地板具有節(jié)能、質(zhì)輕美觀和可調(diào)節(jié)溫濕度等特點[1]。但鋪設地板特別是價格便宜的中低檔木板會降低室內(nèi)空氣質(zhì)量，引起室內(nèi)的空氣污染。其中最具有代表性的空氣污染成分就是甲醛。學者李艷麗指出甲醛濃度在每立方米空氣中達到0.06～0.07mg/m3時，兒童就

會發(fā)生輕微氣喘。當室內(nèi)空氣中甲醛含量為0.1mg/m3時，就有異味和不適感;達到0.5mg/m3時，可刺激眼睛，引起流淚;達到0.6mg/m3，可引起咽喉不適或疼痛。濃度更高時，可引起惡心嘔吐，咳嗽胸悶，氣喘甚至肺水腫;達到30mg/m3時，會立即致人死亡[2]。學者李庭進行了民用住宅空氣污染擴散數(shù)值模擬的研究，其所研究的住宅全部鋪設地板[3]。大多學者對地板散發(fā)甲醛的研究存在一個慣性即認為室內(nèi)全部鋪設地板[4-6]。目前家庭裝修中，客廳、餐廳走動多，若用地板，耐磨性差易磨損，因此，客廳餐廳用地磚，臥室、書房用地板 [7]。前人已經(jīng)在利用自然通風控制室內(nèi)污染方面做了許多的研究[8-10]，本文所研究的內(nèi)容區(qū)別以往研究的地方是甲醛的釋放源為主臥和次臥的地板，房間其他地方不存在甲醛的釋放。

本文研究對象為哈爾濱的某一住宅，其具有冬季東北地區(qū)房間的共性即密封性好，房間漏風很小，可以忽略不計[11-12]。趙毅仁和孫世鈞學者提出了哈爾濱節(jié)能舉措，例如推廣節(jié)能窗等，這些舉措都在一定程度上增加了房間的密封性[13]。由于房間具有良好的保溫能力，本文忽略了室內(nèi)外的熱量傳遞。本文以甲醛作為代表性氣體污染物，采用數(shù)值模擬的方法分析房間密閉條件下和自然通風條件下甲醛的擴散規(guī)律，總結(jié)室內(nèi)污染物的控制方法，對冬季東北地區(qū)室內(nèi)空氣污染的控制具有指導意義[14-15]。

1?物理模型

模擬對象為哈爾濱某一住宅，住宅具體平面結(jié)果如圖1所示。整個戶型南北（北向為Y軸正方向）朝向，西部由南向北依次為主臥、衛(wèi)生間和次臥，東部由南向北依次為客廳、餐廳和廚房。為了盡可能的物理模型貼合實際，提高數(shù)值模擬的準確性，主臥和次臥設置有雙人床，客廳設置有茶幾，餐廳設置有餐桌。房間的幾何模型如圖2所示。

2?數(shù)學模型

室內(nèi)的空氣視為不可壓縮理想流體，利用理想氣體狀態(tài)方程求解室內(nèi)空氣密度，但室內(nèi)空氣的流動仍以不可壓縮流體的方法求解[16]。室內(nèi)空氣的流動要受到質(zhì)量守恒定律、動量守恒定律和質(zhì)量擴散定律的制約[17-19]。

流態(tài)的判別：密閉條件時室內(nèi)空氣的流態(tài)為層流。

通風條件時Re=4600，由于室內(nèi)結(jié)構(gòu)復雜，空氣流動過程中速度變化較大，所以室內(nèi)空氣的流態(tài)為湍流。

3?基本假設和求解條件

3.1?基本假設

1）室內(nèi)密閉，無室外空氣通過門縫或窗縫等滲入室內(nèi)。

2）忽略室內(nèi)內(nèi)熱源（照明、烹飪、人的散熱等）和室外通過房間的圍護結(jié)構(gòu)向室內(nèi)的傳熱對室內(nèi)甲醛擴散的影響。

3）忽略室內(nèi)固體（墻壁、屋頂、家具等）壁面的吸附作用。

4）甲醛的擴散為純物理過程，沒有化學發(fā)應，污染物不產(chǎn)生消耗。

5）甲醛在氣固交界面上的對流擴散只沿地板的外法線方向，不沿其他方向擴散，為一維質(zhì)擴散，遵循Fick定律。

6）甲醛與室內(nèi)空氣彼此之間也沒有化學反應，互不影響。甲醛為強還原劑，在微量堿性時還原性更強。在空氣中能緩慢氧化成甲酸。

7）氣固交界面上甲醛的傳質(zhì)是均勻的及臥室內(nèi)地板散發(fā)甲醛的單位面積的質(zhì)量流量相同。

3.2?求解條件

邊界條件：

1）通風條件時窗戶一設置為速度進口，速度的大小取冬季哈爾濱室外的平均風速V=2m/s。

2）通風條件時新風溫度取冬季哈爾濱室外的平均溫度T=253.15K。

3）通風條件時窗戶二和窗戶三為壓力出口，壓力等于室外大氣壓力P=101.325kg。

4）通風條件時室內(nèi)地板供暖溫度為T=293.15K。

5）密閉條件時臥室一與臥室二地板甲醛的散發(fā)率R采用二階動態(tài)擴散模型[20]R=aeb+ced。

當t<5h時，甲醛釋放速率衰減的較快a=18660，b=-1.16t，c=0;

當5h≤t≤10h時，甲醛釋放速率衰減變慢，a=16950，b=-1.157t，c=102.6，d=-0.003552t;

當t>5h時，甲醛的釋放速率趨于穩(wěn)定，a=33.63，b=-0.007705t，c=68.48，d=-0.0009882t（單位μg/m2·h）。

本文利用c語言將甲醛的散發(fā)率變?yōu)樵错棓U散函數(shù)，作為邊界條件嵌入地板條件中。通風條件時臥室一和臥室二的地板以恒定速率釋放甲醛R=1.83×10-9μg/m2·h。

初始條件：

1）室內(nèi)壓強為一個標準大氣壓。

2）室內(nèi)甲醛濃度為零。

3）室內(nèi)溫度T=293.15K。

4?密閉條件下甲醛濃度計算結(jié)果及分析

哈爾濱地區(qū)冬季室外嚴寒其冬季采暖室外計算溫度為-26℃。人們?yōu)榱吮Ｅ?，防止室外冷空氣通過窗縫或門縫等滲入室內(nèi)，一般都會采取各種各樣的保溫措施，例如給門窗加裝和更換密封條、玻璃貼保溫膜和涂刷保溫涂料等使冬季東北地區(qū)住宅具有密封性好特點，在進行密閉條件下甲醛擴散的數(shù)值模擬時忽略窗戶的漏風。東北地區(qū)在冬季，室內(nèi)外溫度差較大，但由于住宅墻壁和窗戶的導熱熱阻大，所以在短時間內(nèi)忽略室外溫度變化和室內(nèi)微弱的自然對流對室內(nèi)溫度的影響即認為室內(nèi)恒溫，室內(nèi)外無熱量的傳遞，室內(nèi)的傳質(zhì)方式為質(zhì)擴散。我國《民用建筑工程室內(nèi)環(huán)境污染控制規(guī)范》中，對室內(nèi)空氣采樣前要求的密閉時間是1h。本文研究室內(nèi)密閉1h時，甲醛濃度在水平方向和高度方向上的分布，更好的掌握密閉條件下臥室地板釋放的甲醛在室內(nèi)的三維擴散過程。監(jiān)控面分別為Z=10mm，X=4800mm。直線line-1的兩端點為（2000，4000，0）、（2000，4000，2700）（單位mm）處于臥室一中，直線line-2的兩端點為（4000，4000，0）、（4000，4000，2700）（單位為mm）處于客廳中。下圖中，顏色由藍至紅代表污染物濃度由低到高。

圖3顯示的是水平面Z=10mm，的甲醛濃度場，主要研究甲醛在房間水平方向上的擴散。如圖3所示，當臥室一和臥室二的門打開時地板散發(fā)的甲醛會從臥室一和臥室二不斷向其他房間擴散，隨著距臥室一門與臥室二門的距離越遠，甲醛的濃度越小，濃度梯度也越小。根據(jù)質(zhì)擴散菲克定律，當房間密閉時，甲醛擴散的決定性因素為濃度梯度，主要方式為分子運動，當濃度梯度為零時，甲醛在房間的質(zhì)量通量密度也為零。由于臥室一門和臥室二門對向，臥室一與臥室二的向外擴散的甲醛會在兩門之間發(fā)生碰撞，在碰撞交界處甲醛的濃度梯度為零而在其他地方仍然存在較大的濃度梯度，進而改變了其擴散方向，最終導致臥室一地板散發(fā)的甲醛主要向客廳和餐廳擴散，臥室二地板散發(fā)的甲醛主要向衛(wèi)生間、餐廳和廚房擴散。擴散現(xiàn)象是氣體分子的內(nèi)遷移現(xiàn)象。從微觀上分析是大量氣體分子做無規(guī)則熱運動時，分子之間發(fā)生相互碰撞的結(jié)果。由于不同空間區(qū)域的分子密度分布不均勻，分子發(fā)生碰撞的情況也不同。這種碰撞迫使密度大的區(qū)域的分子向密度小的區(qū)域轉(zhuǎn)移，最后達到均勻的密度分布 。甲醛擴散過程中會受到室內(nèi)餐桌、茶幾和墻壁等物體的阻擋，在餐桌、茶幾和墻壁的周圍甲醛的濃度相對較低，主要是因為甲醛分子在運動過程中會撞擊墻壁和家具，家具和墻壁會給甲醛分子施加反作用力，其改變了甲醛分子的運動方向進而改變了甲醛的擴散方向最終導致在家具和墻壁周圍的甲醛的濃度相對較低。

圖4顯示的是截面x=4800mm的甲醛濃度場，主要研究甲醛在房間高度方向上的擴散，由于擴散時間較短，房間上部大部分甲醛濃度為零。對比圖3和圖4明顯看出水平方向上甲醛的濃度梯度大于高度方向上的甲醛濃度梯度，導致甲醛在室內(nèi)水平方向上的擴散速率大于在高度方向上的擴散速率，造成這種現(xiàn)象的原因為甲醛在高度方向的擴散受到了重力的抑制。如圖4所示甲醛從地板逐漸向上擴散，出現(xiàn)了明顯的分層現(xiàn)象，甲醛的濃度隨著高度的增加，在同一水平面上趨于一致。圖5很好的驗證了這一現(xiàn)象，大約在0～150mm的范圍內(nèi)兩曲線的間隔逐漸縮小，當大于150mm時兩直線基本重合。主要原因為甲醛在高度方向的擴散速率較慢，隨著高度的增加房間甲醛的濃度衰減的較快，加之甲醛在房間的各個高度都存在擴散速率相對較快的水平方向上的擴散，所以在房間內(nèi)隨著高度的增加水平面的甲醛濃度差別逐漸減小。

5?通風條件下甲醛濃度計算結(jié)果及分析

東北地區(qū)冬季室內(nèi)外溫差較大，進行自然通風時通風方式的選擇要考慮人的熱舒適性。由于主臥和次臥裝修所用的地板為甲醛的釋放源，東北冬季的主導風向為由南向北，選用窗戶一為房間的進風口，窗戶二和窗戶三為房間的出風口，當房間進行自然通風時，既可以降低房間內(nèi)甲醛的濃度又兼顧了人的熱舒適性。監(jiān)控對象為Z=1350mm。

為了研究通風條件下甲醛濃度場，首先研究室內(nèi)氣體的速度分布。如圖6和圖7所示室外空氣通過窗戶一進入主臥后，射流斷面的速度場從射流中心開始逐漸向邊界衰減，并沿射程不斷變化，流量沿程增加，射流直徑加大，但在各斷面上的總動量保持不變?？梢詫⑸淞鞣譃檩S心速度不變的起始段和軸心速度變小的主體段。其形成原因是由于射流邊界與周圍介質(zhì)間的紊流動量交換，周圍空氣不斷卷入，射流不斷擴大，形成了非等溫受限射流。對比圖6，圖8和圖9可以看出溫度場，濃度場與速度場存在相似性，只是射流邊界比速度分布的邊界有所擴大，其主要原因就是射流與周圍空氣的摻混的結(jié)果。由于低壓和墻壁約束的共同作用，主臥內(nèi)在射流的東西兩側(cè)形成了漩渦，由于射流的出口及主臥室門在臥室的東側(cè)所以射流的左側(cè)的漩渦遠遠大于右側(cè)的漩渦。室外新鮮空氣進入主臥后會攜帶臥室一內(nèi)的含甲醛氣體通過主臥的門進入室內(nèi)其他房間，由于衛(wèi)生間的壓強較小，室內(nèi)主流會向衛(wèi)生間分流，分流在慣性的作用下貼附衛(wèi)生間的墻壁流動導致混合空氣在衛(wèi)生間內(nèi)形成渦流。室內(nèi)主流由于主臥門的作用導致過流斷面變小，流速增大，卷吸作用增強。室內(nèi)主流在剛出主臥門的時候會卷吸客廳和餐廳的部分空氣，隨后室內(nèi)主流流速減小并被房間的墻壁一份為二，一部分流向次臥另一部分部分流向廚房。流向廚房的分流，當?shù)竭_廚房時流速減小，在低壓強的作用下會向餐廳分流，再加上墻壁的約束，使餐廳與客廳的交接處形成了整個房間內(nèi)最大的漩渦，客廳其他部分的混合氣體的流動主要是在此漩渦的帶動作用下形成。次臥的分流主要是沿著房間的分隔壁流向窗戶二，次臥的其他位置在低壓和墻壁約束的作用會生成漩渦。廚房的主流主要是沿著東墻流向窗戶三，廚房的其他位置在低壓和墻壁約束的作用下也會生成漩渦。

如圖9所示各房間內(nèi)甲醛濃度差異較大。雖然客廳和衛(wèi)生間的氣體流速較小，但客廳和衛(wèi)生間的甲醛濃度相對于主臥和次臥分布均勻，不存在濃度過高的死角。其主要為室內(nèi)質(zhì)量的傳遞方式從質(zhì)擴散變?yōu)榱速|(zhì)量傳遞速率更快的對流擴散，主臥和次臥地板生成的甲醛主要是跟隨室內(nèi)主流流動，所以主臥和次臥向其他房間傳遞轉(zhuǎn)移的甲醛量較小，主臥和次臥內(nèi)局部存在漩渦，其氣體的較難得到更新加之處于甲醛的污染源內(nèi)進而導致了臥室一和臥室二存在污染物濃度過高的死角。臥室二的漩渦占據(jù)了臥室二的較大的一部分的空間，自然通風結(jié)束后，隨著甲醛的擴散方式有對流傳質(zhì)變?yōu)橘|(zhì)擴散，臥室二的甲醛平均濃度將會達到最高，所以剛?cè)胱》块g時，盡量住在臥室一而不是臥室二，可以采用其他方式輔助降低臥室二的甲醛濃度如擺放綠色植物等。

PMV指標代表了對同一環(huán)境絕大多數(shù)人的冷熱感覺，因此可用PMV指標預測熱環(huán)境下人體的熱反應。由于人與人之間生理的差別，故用預期不滿意百分率（PPD）指標來表示對熱環(huán)境的不滿意百分數(shù)。通過表1與圖8可以看出，進行自然通風時客廳的熱舒適性較好，甲醛濃度遠小于國家標準值，所以窗戶一進風，窗戶二與窗戶三出風，陽臺關(guān)閉的通風方式比較合理。

6?結(jié)?論

1）在重力的作用下，甲醛在房間水平方向上的擴散速率大于房間高度上的擴散速率，進而導致了在高度方向上，甲醛濃度分層明顯并且隨著高度的增加房間在水平面上甲醛的濃度差別逐漸變小。

2）臥室一和臥室二存在甲醛濃度過高的死角，特別是臥室二，渦旋占據(jù)了房間較大一部分區(qū)域，自然通風除甲醛的效果不好，采用其他方法除甲醛如在室內(nèi)擺放綠色植物、活性炭等。

3）當采取窗戶一進風，窗戶二和窗戶三排風的自然通風方式時可以較為有效降低室內(nèi)各個房間的甲醛濃度，并且客廳的熱舒適相對較好但主臥與次臥存在甲醛濃度較高的死角，仍然需要采用其他方式去除甲醛。

參 考 文 獻：

[1]?靳俊杰，唐中華，高理福，等.熱舒適性情況下的低溫地板輻射供暖節(jié)能研究[J].四川建筑科學研究，2013，39（4）： 370-377.

[2]?李艷麗，尹詩，黃寶妍.室內(nèi)甲醛污染來源及其對人體危害[J].佛山科學技術(shù)學院學報，2003，21（1）： 49-52.

[3]?李庭.民用單元住宅空氣污染擴散數(shù)值模擬研究[D].哈爾濱：哈爾濱工業(yè)大學，2015：20-54.

[4]?駱偉，劉榮華，邱汗峰.自然通風作用下室內(nèi)甲醛濃度分布數(shù)值模擬分析[J].制冷與空調(diào)，2013，27（5）：504-508.

[5]?張牛牛，薛敏.室內(nèi)污染物排放的數(shù)值模擬分析[J].江蘇工業(yè)學院學報，2009，21（2）：27-30.

[6]?馬飛，羅海亮，謝慧，等.建筑室內(nèi)污染物擴散模擬與控制[J].北京科技大學學報，2007，29（2）：16-18.

[7]?彭沖.實木拼接地板制造工藝及性能研究[D].廣州：華南農(nóng)業(yè)大學，2016：1-20.

[8]?曹國慶，張益昭.?通風與空氣過濾對控制室內(nèi)生物污染的影響研究[J]. 土木建筑與環(huán)境工程，2009，31（1）：130-135.

[9]?謝偉樊，越勝，黃亦平，等.自然通風室內(nèi)外PM10、PM2.5、PM1.0污染特征[J].潔凈與空調(diào)技術(shù)，2015，（3）：58-61.

[10]肖楚璠.自然通風條件下室內(nèi)甲醛擴散模擬研究[J].潔凈與空調(diào)技術(shù).2013（1）：15-18.

[11]張姝，張希浩.嚴寒地區(qū)既有居住建筑節(jié)能改造策略研究[J].應用科技，2017，44（4）：82-86.190-191.

[12]陸書齋，呂平.哈爾濱市住宅建筑外墻保溫技術(shù)的研究[J].山西建筑，2015，41（35）：190-191.

[13]趙毅仁，孫世鈞.哈爾濱建筑節(jié)能的探討[J].黑龍江科技信息，2003（11）：198.

[14]MOHAMED M F，KING S，BEHNIA Metal.?The Effects of ?Balconies on the Natural Ventilation Performance of Cross-ventilated High-rise Buildings[J]. journal of green ?building，2014，9（2）：145-160.

[15]王韌剛，吳志湘，岳斌佑.?利用CFD技術(shù)預測工業(yè)廠房通風降溫效果的必要性及注意事項[J]. 潔凈與空調(diào)技術(shù)，2013，（4）：22-24.

[16]郭雅迪，張人偉，劉曰帥，等.?基于FLUENT模擬的選煤廠煤倉瓦斯超限治理研究[J]. 安全與環(huán)境工程，2014，（3）：148-153.

[17]龔波.教學樓風環(huán)境和自然通風教室數(shù)值模擬研究[D].成都：西南交通大學，2005：47-49.

[18]王樂.辦公建筑自然通風研究的CFD分區(qū)分析方法探索[D].天津：天津大學，2005：8-16.

[19]陶文銓.數(shù)值傳熱學[M]. 西安：西安交通出版社，2001，16-17.

[20]嚴彥，王光學，楊旭.?木質(zhì)人造板材甲醛釋放規(guī)律的研究[J].環(huán)境科學學報，2003，23（1）：134-137.

（編輯：王?萍）