畢業(yè)論文--通風對室內有害氣體凈化效應

1、<p>　　本科畢業(yè)設計（論文）</p><p>　　通風對室內有害氣體凈化效應</p><p><b>　　年 月 日</b></p><p>　專業(yè)名稱：理論與應用力學</p><p>　年級班級：力學10-01班</p><p>　學生姓名：</p><p&

2、gt;　指導教師：</p><p>　　畢業(yè)設計（論文）任務書</p><p>　　二、起止日期 2014 年 4 月 1日 至 2014 年 6 月 10日</p><p><b>　　三、主要任務與要求</b></p><p>　　（一）建立三維實物模型，選擇合適的湍流封閉模型和基本微分方程的求解方法。主要講解多相

3、流模型和標準湍流模型。</p><p>　?。ǘ┦褂胒luent軟件采用標準的方程和多相流混合模型進行求解</p><p>　　求解在風速為1m/s時，木質地板散發(fā)的甲醛的濃度分布</p><p>　　求解在風速為2m/s時，木質地板散發(fā)的甲醛的濃度分布；</p><p>　　求解在風速為3m/s時，木質地板散發(fā)的甲醛的濃度分布；</

4、p><p>　　求解在風速為1m/s時，木質地板、沙發(fā)、立柜散發(fā)的甲醛的濃度分布；</p><p>　　求解在風速為2m/s時，木質地板、沙發(fā)、立柜散發(fā)的甲醛的濃度分布；</p><p>　　求解在風速為3m/s時，木質地板、沙發(fā)、立柜散發(fā)的甲醛的濃度分布；</p><p>　?。ㄈ┍容^上述不同情況下，室內甲醛氣體的濃度分布，得出有利于室內甲醛

5、擴散的比較好的風速。</p><p>　　（四）最后，對本文的研究工作進行總結、歸納，并提出需要進一步研究的工作。</p><p>　　指導教師： 職稱： </p><p>　　院 領 導： 簽字（蓋章）</p><p><b>　　年 月 日</b>&

6、lt;/p><p>　　通風對室內有害氣體凈化效應</p><p><b>　　摘要</b></p><p>　　室內空氣質量的好壞與人類健康密切相關。新型建筑材料、裝潢材料和家用化學材料的推廣使用，造成室內空氣污染嚴重，而甲醛是室內空氣污染危害最嚴重且最常見的是污染物之一。有效的通風可以控制室內甲醛濃度而使其很容易地符合設計要求。</p&g

7、t;<p>　　本文采用Fluent 軟件來模擬三維空調室內氣流組織的流動形式，研究了送風速度對僅有木質地板和放置沙發(fā)、立柜的單居室室內氣流分布的影響及甲醛濃度的分布影響，且在Fluent 軟件和Tecplot中以直觀的圖像方式表示出了各種不同氣流組織和送風速度方案下甲醛濃度分布場。并通過對比，得到只有木質地板釋放甲醛時，得出送風速度不宜太大，以2m/s 為最佳。在木質地板和放置沙發(fā)、立柜都釋放甲醛時，3m/s為最佳。&l

8、t;/p><p>　　關鍵詞：室內 送風速度 甲醛擴散 數(shù)值模擬</p><p>　　Ventilation for indoor harmful gas purifying effect</p><p><b>　　Abstract</b></p><p>　　The stand or fall of indoor

9、 air quality are closely related to human health. New building materials, decoration materials and household chemical materials to promote the use, cause indoor air pollution, indoor air pollution and hazards of formalde

10、hyde is the most serious and most common is one of contaminants.Effective ventilation of indoor formaldehyde concentration can be controlled so that it is easy to meet the design requirements.</p><p>　　This

11、paper uses the Fluent software to simulate the flow in the form of three-dimensional air-conditioned indoor air, studied the wood floors and air velocity only place sofa, distribution on single bedroom closet distributio

12、n of indoor air concentrations of formaldehyde, and in the Fluent software and the Tecplot in an intuitive way of showing a picture of formaldehyde concentration distribution field under various airflow and air velocity

13、program. By contrast, only the wooden floor to get the </p><p>　　Key words: indoor; air supply velocity; formaldehyde diffusion; numerical simulation</p><p><b>　　目錄</b></p>&l

14、t;p><b>　　1緒論1</b></p><p><b>　　1.1選題背景1</b></p><p>　　1.2國內外對室內污染的研究現(xiàn)狀1</p><p>　　1.3本文主要的研究內容3</p><p>　　2室內甲醛散發(fā)數(shù)值模擬理論5</p><p>

15、;　　2.1室內甲醛散發(fā)模擬的基本原理5</p><p>　　2.2湍流流動的Reynolds時均方程5</p><p>　　2.3 湍流封閉模型6</p><p>　　2.4氣態(tài)物質的散發(fā)和分布模擬7</p><p><b>　　3 模型的建立9</b></p><p>　　3.1 室

16、內空氣流動模型的建立9</p><p>　　3.2幾何模型的建立10</p><p>　　3.3網(wǎng)格的劃分11</p><p>　　3.4 邊界條件12</p><p>　　3.5迭代的SIMPLE算法12</p><p>　　4 方程的求解及結果的后處理14</p><p>　　4

17、.1 Fluent里面的設置14</p><p>　　4.2數(shù)值模擬的結果及比較16</p><p>　　4.3數(shù)值模擬結果16</p><p>　　4.3.1不同風速下y=0截面甲醛濃度分布16</p><p>　　4.3.2不同風速下z=0.2截面甲醛濃度分布20</p><p>　　4.3.3放置沙發(fā)立

18、柜不同風速下y=0截面甲醛濃度分布22</p><p>　　4.3.4放置沙發(fā)立柜不同風速下z=-0.3截面甲醛濃度分布24</p><p><b>　　5結論29</b></p><p><b>　　5.1結論29</b></p><p><b>　　5.2展望29</b

19、></p><p><b>　　參考文獻30</b></p><p><b>　　致謝31</b></p><p><b>　　1緒論</b></p><p>　　1.1選題背景（每一部分都應該有相應的參考文獻）</p><p>　　隨著國民經(jīng)

20、濟的發(fā)展和人民生活水平的提高，人們對自己的居住條件要求也不斷提高。然而室內建材及建筑裝修材料散發(fā)出大量的揮發(fā)性有機化合物嚴重惡化了室內空氣品質，甲醛就是其中常見的一種揮發(fā)性有機污染源，已經(jīng)被世界衛(wèi)生組織確定為致癌和致畸形物質，成為室內空氣污染中人們所關注的焦點。</p><p>　　甲醛是最簡單、最常見的醛類物質，其理化性質為：無色具有強烈刺激性氣味，沸點- 1915℃ 比重1106，易溶于水、醇和醚，是一種溶解

21、度很大、揮發(fā)性很強的有毒物質[1]。當空氣中的甲醛濃度超過0.6mg/m3時，人的眼睛會受到刺激，咽喉會感到不適和疼痛，在含甲醛10ppm的空氣中停留幾分鐘，眼睛就會流淚不止。長期接觸低劑量甲醛可以引起慢性呼吸道疾病、女性月經(jīng)紊亂、妊娠綜合癥，引起新生兒體質降低、染色體異常，甚至引起鼻咽癌。高濃度甲醛對神經(jīng)系統(tǒng)、免疫系統(tǒng)、肝臟等都有毒害。甲醛還有致畸、致癌的作用，據(jù)流行病學調查，長期接觸甲醛的人，可以引起鼻腔、口腔、鼻咽、咽喉、皮膚和消

22、化道的癌癥[2]。</p><p>　　據(jù)調查，新裝修的房間甲醛的超標率達60 %～100 %[3]。來自室內裝修和裝飾材料。用作室內裝飾的膠合板、細木工板、中密度纖維板的刨花板等人造板材，因為使用的膠粘劑以甲醛為主要成份的脲醛樹脂，板材中殘留的和末參與反應的甲醛會逐漸向周圍環(huán)境釋放，是形成室內空氣中甲醛的主體。含有甲醛成分并有可能向外界散發(fā)的其他各類裝飾材料，比如貼墻布，貼墻化纖在地毯、泡沫塑料、油漆和涂料等。

23、甲醛已成為室內空氣污染物的主要成分[4]。</p><p>　　甲醛在室內的濃度變化主要與室內污染源的釋放量和釋放規(guī)律有關，也與使用年限、室內溫度、濕度及通風強度等因素有關，其中通風的影響最為重要。[5]</p><p>　　所以把室內污染物定為極有代表性的甲醛，凈化方式定為影響最大的方式——自然通風。</p><p>　　1.2室內污染物擴散的研究現(xiàn)狀（）半角全角

24、標點</p><p>　　由于通風房間室內空氣流動多為湍流流動（Turbulent Flow），而人們對湍流機理的認識至今還不完全，在工程應用方面，目前多為半經(jīng)驗的、唯像的模擬，所以對湍流的模擬始終成為室內空氣數(shù)值模擬的一個主要研究方向。限于目前的計算機能力，工程中最常采用的是渦粘系數(shù)模型EVM（eddy viscosity model）中的標準k-ε模型或其變形。而k-ε模型對于等溫流動情形能模擬得很好，鑒于模

25、擬室內污染物的擴散多在等溫條件下，不考慮熱傳導，故采用標準的k-ε二方程模型就能夠滿足需求。近年來國內一些學者對于室內污染物擴散進行了數(shù)值模擬，結果如下：</p><p>　　付小平通過數(shù)值模擬，得出了室內氣體污染物濃度隨通風時間的變化，并利用數(shù)值模擬結果和實驗數(shù)據(jù)證明：標準的k-ε模型和組分輸運模型能較好地模擬室內上送下回的通風方式稀釋室內污染物的過程[6]。</p><p>　　吳猛等

26、采用Fluent6.1計算軟件，用大渦模擬方法（LES）對整套住宅在自熱通風時室內空氣污染物排放過程進行了模擬，根據(jù)不同的房間類型及通風口布置方式設置兩種工況，得到流場和污染物質量分數(shù)場數(shù)據(jù)[7]。</p><p>　　張俊杰等利用計算流體力學（CFD）方法模擬了不同通風方式和不同送風速度下室內污染物的濃度分布，模擬和分析結果表明，適當?shù)乃惋L速度可以有效降低室內污染物的濃度[8]。 </p><

27、;p>　　劉娣等用CFD方法模擬了空調臥室內制冷制熱運行時3種不同的住宅空調模式(包括普通窗式空調器、分體機及具有引入新風的熱回收裝置的窗式空調器)分別位于高位置和低位置時室內空氣溫度及流速、有機污染物(甲醛)濃度及CO2的分布，并進行比較．結果表明，空調器的類型、位置及新風量對空氣環(huán)境影響較大。夏季制冷運行時，帶熱回收裝置的窗式空調器置于低位置時可以獲得良好的室內流場分布，稀釋和攜帶走室內的CO2和污染物；而該裝置置于高處時，流

28、場結構不合理；其它空調模式下由于沒有引入新風，產(chǎn)生室內污染物堆積。冬季制熱運行與夏季制冷運行時的結論相同[9]。</p><p>　　張牛牛等利用計算流體力學(CFD)中的Fluent6.2軟件，模擬了兩種氣流組織形式分別在3種通風速度下室內氣流組織的速度場和甲醛的濃度場分布。對得出的速度場和濃度場進行分析：每種氣流組織形式的速度場和濃度場都有較好的耦合性;其中下側進上出的通風方式能更好的排除污染物；另外，適當增

29、大送風速度能更好的降低污染物濃度[10]。</p><p>　　馬飛等對建筑室內香料污染物擴散的控制方法進行了研究。分析了香料污染物產(chǎn)生的根源和擴散特性，制定了分區(qū)控制、先收集后處理的控制策略。設計了排氣柜排風系統(tǒng)，基于數(shù)值模擬優(yōu)化了排風系統(tǒng)參數(shù),確定了最優(yōu)風道流速5.6-6.0 m/s實測結果表明，采用分區(qū)控制污染源，利用排氣柜收集污染物是控制香料污染物擴散的有效方式,香料污染物收集率達90%以上[11]。&l

30、t;/p><p>　　胡平放等以計算流體力學和傳熱學為基礎，利用FLUENT軟件建立了四種典型氣流組織形式的物理及數(shù)學模型，運用k-ε模型對室內氣流組織進行了三維數(shù)值模擬。對各種氣流組織形式下建筑物室內的流場、溫度場及濃度場的數(shù)值模擬結果進行分析比較。結果表明:置換通風在溫度場、速度場、污染物濃度場、通風效率及人體熱舒適性等方面均優(yōu)于混合通風的幾種氣流組織形式[12]。</p><p>　　數(shù)

31、值模擬法是通過計算機數(shù)值計算和圖像顯示，對包含有流體流動和熱傳導等相關物理現(xiàn)象的系統(tǒng)所做的分析。它可以模擬室內空氣中氣流的運動狀態(tài)和污染物在空氣中的分布狀況，使研究人員更深刻地認識室內污染物的分布規(guī)律。本次畢業(yè)設計旨在利用Fluent數(shù)值模擬室內污染物（簡化為甲醛）在不同空調送風形式下擴散的濃度分布，并對其結果進行分析，選擇出有利于甲醛擴散的氣流形式，希望對理論研究起到一定的輔助作用。</p><p>　　1.3

32、本文主要的研究內容</p><p>　　人們現(xiàn)在都非常重視環(huán)境問題，因為環(huán)境直接影響人類的生存與發(fā)展。環(huán)境問題包含兩個方面：室內環(huán)境與室外環(huán)境。據(jù)統(tǒng)計人們有80％的時間生活工作都在室內，隨著社會發(fā)展和人們生活水平的提高，工作居住環(huán)境不斷改善，許多建筑物都進行了高級裝修，而這些高級裝修建材中就含有大量的有毒物質，它們種類繁多，分布廣泛，其中含量最多的、最有代表性的是甲醛[13]。所以本文就對一般室內裝修中含有的甲醛

33、進行最有效的通風處理，運用Fluent軟件對通風條件下甲醛的濃度進行精確有效的分析。</p><p>　　選取單居室作為研究的對象，考慮到釋放甲醛的家具有很多，比如沙發(fā)、立柜、木質地板、粉刷的墻壁。就采取先單個研究然后再綜合研究的策略。具體實施的步奏如下：</p><p>　?。ㄒ唬┙⑷S實物模型，選擇合適的湍流封閉模型和基本微分方程的求解方法。主要講解多相流模型和標準k-ε湍流模型。&

34、lt;/p><p>　?。ǘ┦褂肍luent軟件采用標準的k-ε方程和多相流混合模型進行求解</p><p>　?、偾蠼庠陲L速為1m/s時，木質地板散發(fā)的甲醛的濃度分布</p><p>　?、谇蠼庠陲L速為2m/s時，木質地板散發(fā)的甲醛的濃度分布；</p><p>　?、矍蠼庠陲L速為3m/s時，木質地板散發(fā)的甲醛的濃度分布；</p>

35、<p>　?、芮蠼庠陲L速為1m/s時，木質地板、沙發(fā)、立柜散發(fā)的甲醛的濃度分布；</p><p>　?、萸蠼庠陲L速為2m/s時，木質地板、沙發(fā)、立柜散發(fā)的甲醛的濃度分布；</p><p>　?、耷蠼庠陲L速為3m/s時，木質地板、沙發(fā)、立柜散發(fā)的甲醛的濃度分布；</p><p>　?。ㄈ┍容^上述不同情況下，室內甲醛氣體的濃度分布，得出有利于室內甲醛擴散的

36、比較好的風速。</p><p>　?。ㄋ模┳詈?，對本文的研究工作進行總結、歸納，并提出需要進一步研究的工作。</p><p>　　2室內甲醛散發(fā)數(shù)值模擬理論</p><p>　　2.1室內甲醛散發(fā)模擬的基本原理</p><p>　　室內甲醛擴散模擬技術的兩大理論核心是計算流體動力學和數(shù)值分析方法，它是使用一組微分方程描述空間中流體的流動情況，

37、首先對計算區(qū)域進行離散，在此基礎上將微分方程轉化為代數(shù)方程，求解出微分方程組的數(shù)值解，然后再獲得流場的相關性質。</p><p>　　2.1.1湍流流動的Reynolds時均方程</p><p>　　無論湍流是多么復雜，其流動過程仍然由連續(xù)介質力學的幾個基本方程控制。如果忽略大氣的可壓縮性、溫度變化因素以及砌體力的作用，則這些基本方程包括連續(xù)性方程和Navier-Stokes方程[5]，即

38、</p><p><b>　?。?-1）</b></p><p><b>　?。?-2）</b></p><p>　　式中(i=1,2,3)為方向的速度分量，為壓力，為空氣密度，為氣流的運動粘性系數(shù)，波紋符“～”表示瞬時量。</p><p>　　連續(xù)性方程（2-1）是由鈍體空氣動力學的質量守恒方程，

39、假定空氣密度為常數(shù)得到的，它反映了大氣流動的質量守恒原理。</p><p>　　由于流體內部不同尺度渦旋的隨機運動是造成湍流的一個重要特點，所以湍流運動過程中，在空間任一點的速度和壓力都隨時間不斷的無規(guī)則的變化著。對給定系統(tǒng)的任何兩次測量都不可能是相同的，但是湍流量的統(tǒng)計平均卻有確定性的規(guī)律可循，平均值在各次試驗中是可重復實現(xiàn)的。因此Reynolds首先提出將各瞬時量、分解成平均量(用大寫字母表示)和脈動量(用小

40、寫字母表示)之和[6]，</p><p><b>　?。?-3）</b></p><p>　　將式(2-3)代入到方程(2-1)、(2-2)中，并對方程中的每一項作平均化運算可得到平均量的控制方程：</p><p>　　這就是湍流流動的Reynolds時均方程，方程中除了脈動量的二階關聯(lián)量以外，其余在形式上與方程(2-1)、(2-2)完全相同。

41、二階關聯(lián)項也稱Reynolds應力，它代表了脈動速度對平均流動的影響。但是Reynolds應力是時均方程組中新出現(xiàn)的未知量，從而使方程的未知量數(shù)目超過了方程數(shù)，這就使得方程不封閉。而解決這一問題的基本方法是對Reynolds應力作出假設，即建立湍流封閉模型。把未知的更高階的時間平均值表示成較低階在計算中可以確定的函數(shù)，從而求解平均方程[14]。</p><p>　　2.1.2 湍流封閉模型</p>

42、<p>　　建立湍流封閉模型的方法，大致可分為兩大類：一類是建立在布辛涅斯克渦粘性系數(shù)模型，另一類稱為Reynolds應力方程模型(Reynolds Stress Model，即RSM)，它直接建立以Reynolds應力為因變量的方程式并通過?；刮闪髌骄骺刂品匠探M封閉。</p><p>　　對于布辛涅斯克渦粘性系數(shù)法，首先引入了Boussinesq假設，仿照層流運動應力，湍流脈動所造成的附加應力也可

43、以同時均的切應變率關聯(lián)起來，湍流脈動所造成的應力可以表示為：</p><p><b>　?。?-4）</b></p><p>　　式中，是紊流渦粘性系數(shù)，k為紊流動能，表示單位質量流體紊流脈動動能的平均值。</p><p>　　由三維非穩(wěn)態(tài)的方程可導出?；蟮膋方程與方程：</p><p><b>　　(2-5

44、)</b></p><p><b>　　(2-6)</b></p><p>　　時均化的連續(xù)性方程，方程和模型化了的k方程(2-5)、方程(2-6)以及Reynolds應力的渦粘性系數(shù)表達式(2-4)組成了標準k-ε雙方程模型的封閉的控制方程組[15]。</p><p>　　2.1.3氣態(tài)物質的散發(fā)和分布模擬</p>

45、<p>　　利用數(shù)值模擬的方法確定污染物散發(fā)主要基于VOC分子從材料中進入空氣的兩個過程：①由于濃度差的存在，在材料內部的VOC分子擴散；②材料跟空氣接觸表面的蒸發(fā)過程。對于地毯等“干”性材料，第一個過程其決定作用；而對于剛剛涂刷的油漆涂料等材料等，初試階段的VOC散發(fā)由第二種過程起決定作用。由于純粹實驗的復雜性，短期的VOC散發(fā)性能研究可以通過實驗與模擬結合來進行，主要利用實驗得到的數(shù)據(jù)，結合通過數(shù)值模擬得到的空氣流速、溫度

46、、濕度等狀態(tài)參數(shù)，共同確定VOC散發(fā)的某些特征參數(shù)，如VOC在材料內部的初試濃度、擴散系數(shù)等。將這些參數(shù)進一步應用于數(shù)值模擬之中，可以對材料在各種外界條件下的中長期散發(fā)特性做出預估[16]。</p><p>　　圖2.1污染物散發(fā)模型示意圖</p><p>　　國外學者已近建立了污染物內部和材料表面污染物的擴散模型。如圖2.1所示，其中假定污染物在室內空氣中分布均勻，通過空氣邊界層將材料內

47、外耦合。經(jīng)實驗獲得污染源（如地毯、油漆本身）內部的初始濃度、污染物內部的擴散系數(shù)和兩相（污染源本身固體和室內氣體）結合面處兩個濃度的平均系數(shù)K（目前都是假定相界面處兩側的濃度成線性關系）。這幾個參數(shù)的獲得一般通過實驗的多次采樣和結合一定假設作為初始條件，對數(shù)學模型迭代求解。</p><p>　　VOC散發(fā)模型通?；谌缦录僭O：</p><p>　　材料內部組成均勻，初始時刻各處濃度相同。

48、 </p><p>　　材料內為一維擴散物質。</p><p>　　材料的擴散系數(shù)為常數(shù)。</p><p>　　在固體壁面處，VOC的濃度服從亨利定律。 </p><p><b>　　數(shù)學表達式為：</b></p><p><b>　　(2-7)</b></p&

49、gt;<p>　　， (2-8)</p><p><b>　　(2-9)</b></p><p><b>　　(2-10)</b></p><p><b>　　(2-11)</b></p><p

50、>　　式中，為VOC在材料內部的濃度（μg/m3）；為VOC在室內空氣中的濃度（μg/m3）；為材料內部的VOC分子擴散系數(shù)（m2/s）；K為界面處的VOC濃度的平衡系數(shù)；為固體壁面處的傳質量（μg/m2.s）[17]。</p><p><b>　　3 模型的建立</b></p><p>　　3.1 室內空氣流動模型的建立</p><p>

51、;　　根據(jù)熱工理論基礎，可認為室內空氣滿足氣體狀態(tài)方程，即</p><p><b>　　P=ρPT</b></p><p>　　式中，P為空氣壓力（Pa）；ρ為空氣密度（kg/m）；R為空氣常數(shù)，約為287(kgK)；T為空氣熱力學溫度（K）。</p><p>　　其次，通?？蓪⑹覂瓤諝饬鲃拥膲毫σ暈槌?shù)，于是可得：</p>&

52、lt;p><b>　　ρT=常數(shù)</b></p><p>　　另外，常見的室內環(huán)境中空氣流動基本為低速流動，流速常在10~20m/s以下，因此可將室內空氣當作不可壓流體看待，即</p><p>　　而且，通風房間內空氣和污染物氣體溫度變化不大，也即密度變化不大，因此可認為室內氣體流動符合Boussinesq假設：密度變化并不顯著改變流體性質，動量守恒中，密度的變

53、化對慣性力項、壓力差項和粘性力項的影響可忽略不計，而僅考慮隨質量力項的影響。綜上所述，除了隨質量力項的考慮之外，室內空氣物性都可以當作常物性看待[15]。</p><p>　　最后，室內空氣的粘性不可忽略，必須用粘性流體動力學的理論來研究。進一步而言，室內空氣流動通常都是湍流流動，需要相應的湍流理論來模擬。而且，室內空氣流動往往是自然對流和強迫并存的混合對流，或者是輻射換熱影響的自然對流。當然，也存在通風時的等溫

54、強迫對流動[18]。針對本文，為了簡化模型，僅僅考慮等溫下是自然對流。</p><p>　　綜上所述，室內氣體流動的物理模型可總結如下：</p><p>　　常溫、低速、密度不變、不可壓流體流動；</p><p>　　符合氣體狀態(tài)方程的等壓流動；</p><p>　　符合Boussinesq假設；</p><p>　　

55、自然對流下的湍流流動。</p><p>　　3.2幾何模型的建立</p><p>　　選擇較為常見的一個長為6米，寬為4米，高為3米的房間，左邊墻壁開了一個1米×1.5米的窗，距離地面1米高，右邊開了一扇2米×1.5米的門，挨著地。模型如下圖3-1所示</p><p>　　圖3.1標準單間幾何模型設計圖</p><p>　

56、　圖3.2是放置了沙發(fā)立柜的三維立體房間圖。</p><p>　　圖3.2放置沙發(fā)、立柜幾何模型設計圖</p><p><b>　　3.3網(wǎng)格的劃分</b></p><p>　　在GAMBIT確定了要進行計算的幾何區(qū)域以后，接下來的工作是把這個幾何區(qū)域進行離散化，即進行網(wǎng)格劃分。網(wǎng)格劃分可以先對邊進行網(wǎng)格劃分，然后是面，再是體；也可以直接對體或

57、者面進行網(wǎng)格劃分。網(wǎng)格劃分是一項很重要的工作，網(wǎng)格的質量往往決定著Fluent模擬的成敗或者收斂的快慢。下圖是直接對房間這個體劃分網(wǎng)格的，本文設置的interval size 為0.05，足以保證網(wǎng)格足夠細。</p><p>　　圖3.3標準單間幾何模型網(wǎng)格圖</p><p>　　圖3.4放置沙發(fā)、立柜幾何模型網(wǎng)格圖</p><p><b>　　3.4 邊

58、界條件 </b></p><p>　　在沒有沙發(fā)和立柜的房間中，把木質地板設置為甲醛散發(fā)源其邊界條件設為質量進口(mass flow inlet)，密閉和通風兩種狀態(tài)的釋放強度均為7.2×10-11kg/s[20]，并且不隨時間變化，其它壁面無甲醛釋放。也可以根據(jù)質量流量與密度的轉換關系m=ρVnA，可以根據(jù)質量流量求得速度，因此亦可設置污染物為速度入口邊界。假定木質地面甲醛的釋放速度為0.

59、01m/s。然后就是設置窗口為風的速度入口邊界（intlet），門為速度出口邊界（outlet），其他墻壁均為wall。</p><p>　　在放置沙發(fā)和立柜的房間中，設置沙發(fā)和立柜的每個表面為散發(fā)源，密閉和通風兩種狀態(tài)的釋放強度均為 7.2×10-11kg/s，并且不隨時間變化，其它壁面無甲醛釋放。然后就是設置窗口為進風的速度入口邊界（intlet），門為速度出口邊界（outlet），其他墻壁均為wa

60、ll。</p><p>　　3.5迭代的SIMPLE算法 </p><p>　　迭代求解：在對室內空氣流動控制方程進行離散得到代數(shù)方程后，需要對其求解以獲得流場各物理量的分布信息。常見的求解方法有兩種：一、直接求解法。二、迭代解法。本文所采用的方法為迭代解法，下面就詳細講解迭代解法的原理及過程：</p><p>　　通過前面的講解可知，控制流體流動和傳熱的微分方程都

61、是非線性的，離散方程的系數(shù)本身就是所求解變量的函數(shù)，因此求解過程本身就是迭代性質的：先假定一個未知量場，據(jù)此計算離散方程的系數(shù)，然后求解方程獲得改進值，利用改進值再計算離散方程系數(shù)，再計算獲得新的改進值，如此反復循環(huán)，直至獲得最終的收斂值。</p><p>　　而實際的流場和壓力、溫度（非等溫）以及湍流參數(shù)等有關系，它們是互相影響的，因此計算時需要耦合求解。為此，形成了一系列算法，而最具代表性、應用的最為廣泛的就

62、是SIMPLE算法了。而對控制體的連續(xù)方程利用有限容積法離散，會出現(xiàn)速度場和壓力場不連續(xù)的情況。為了解決這個問題，Spalding和Patankar等人提出了交錯網(wǎng)格（Staggered Grid）的概念，就是對于動量方程，其控制體比其他方程的控制體錯移半個網(wǎng)格，也就是對于速度變量U、V、W，其控制體比其他變量的控制體（溫度、濃度等）在該速度方向上錯移半格。</p><p>　　SIMPLE算法：壓力修正后，就可

63、以對動量方程進行迭代求解了，從而也就可以求解其他變量的方程?；趬毫π拚乃惴ň褪侵腟IMPLE算法了：壓力耦合的半隱式。SIMPLE算法的步驟如下：</p><p><b>　　假定壓力場；</b></p><p>　　求解動量方程得到速度場；</p><p>　　求解壓力修正方程，得到修正的壓力；</p><p>

64、;　　修正壓力場和速度場；</p><p>　　求解其他影響流場的變量，如溫度、濃度或湍流動能等；</p><p>　　判斂，如果收斂則停止迭代，否則以更新后的壓力場為新的假設的壓力場，跳到（2）繼續(xù)計算，直至收斂；</p><p>　　求解與速度無關的其他量。</p><p>　　據(jù)以上算法，再結合前述的代數(shù)方程求解方法，就可以對室內空氣流

65、動和換熱、傳熱的離散方程進行求解，獲得需要的場分布情況[19]。</p><p>　　本文中，利用SIMPLE算法進行迭代求解。</p><p>　　4 方程的求解及結果的后處理</p><p>　　4.1 Fluent里面的設置 </p><p><b>　?。?）計算器的設置</b></p><p

66、>　　表4.1 Fluent里面計算器的設置</p><p>　　（2）氣體參數(shù)的設置</p><p>　　表4.2 Fluent里面氣體參數(shù)的設置</p><p><b>　?。?）多相流的設置</b></p><p>　　表4.3 Fluent里面多相流混合模型的設置</p><p>　

67、?。?）運行條件的設置</p><p>　　表4.4 Fluent里面運行條件的設置</p><p>　?。?）計算精度的設置</p><p>　　表4.5 Fluent里面計算精度的設置</p><p>　　（6）邊界條件的設置</p><p>　　表4.6 Fluent里面邊界條件的設置</p>&l

68、t;p>　　4.2數(shù)值模擬的結果及比較</p><p>　　影響甲醛分布的因素有很多，包括送回風口的形式、位置、規(guī)格、數(shù)量；送風口的風速、風量、入口溫度、室內熱源、房間幾何形狀以及具體房間的功能等等，本節(jié)主要研究影響臥室內通風速度的主要因素：不同風速的改變對甲醛濃度場的影響?？紤]到在Fluent三維場的分布中具體甲醛濃度的分布難以清晰地定量顯示，本文首先截取y=0來顯示房間內甲醛在立面的濃度分布；再截取z=

69、-0.3（距離地面1.2m高，恰好為人坐姿時的呼吸高度）、z=0.2（距離地面1.7m高，為人站立時的呼吸高度）來顯示甲醛在兩個水平面的濃度場的分布情況，來反應甲醛對人體的危害程度大小，并通過對比，選出較好的凈化風速。</p><p><b>　　4.3數(shù)值模擬結果</b></p><p>　　本文主要以Fluent中計算出的結果，用圖片的方式來清晰直觀的展示數(shù)值模擬

70、的結果，經(jīng)過對圖片進行仔細認真的分析得出所要的結果。由于Fluent軟件對計算出來的結果的處理有限，本文還使用了工程圖片處理器Tecplot。 </p><p>　　4.3.1不同風速下y=0截面甲醛濃度分布</p><p>　　圖4-1風速為1m/s時 y=0截面甲醛濃度分布</p><p>　　圖4-2 風速為2m/s時 y=0截面甲醛濃度分布</p

71、><p>　　圖4-3風速為3m/s時 y=0截面甲醛濃度分布</p><p>　　從以上圖中可以看出，甲醛濃度集中于較低的高度，也很直觀的看出隨著風速的增大，下部分甲醛整體濃度都是在下降的，其中當風速由1m/s增加至2m/s時，甲醛濃度下降的很明顯，當再增加到3m/s時，下降的程度是很小的。特別的對于風口以上區(qū)域，在風速3m/s時，甲醛整體濃度都高于風速為2m/s。</p>&

72、lt;p>　　4.3.2不同風速下z=-0.3截面甲醛濃度分布 </p><p>　　圖4-4風速為1m/s 時z=-0.3截面甲醛濃度分布</p><p>　　圖4-5風速為2m/s 時z=-0.3截面甲醛濃度分布</p><p>　　圖4-6風速為3m/s 時z=-0.3截面甲醛濃度分布&l

73、t;/p><p>　　圖4-7 (YOZ，y=0，z=–0.3)沿X軸各點甲醛濃度比較</p><p>　　圖4-4、4-5、4-6這些圖是甲醛分布的等值線圖，圖上也很清楚的標注著每立方米甲醛濃度的數(shù)值。整體來看甲醛濃度都較低，在風速由1m/s增至2m/s甲醛濃度都是降低的。在風速增大至3m/s時，上下兩邊的甲醛濃度都是在顯著增高的。圖4-7是在YOZ，y=0，z=–0.3)沿X軸方向的線圖。

74、可以看出在風口位置，隨著風速的增大，帶動著甲醛濃度的降低。</p><p>　　4.3.3不同風速下z=0.2截面甲醛濃度分布</p><p>　　圖4-8風速為1m/s時 z=0.2截面甲醛濃度分布</p><p>　　圖4-9風速為2m/s時 z=0.2截面甲醛濃度分布</p><p>　　圖4-10風速為3m/s時 z=0.2截面甲醛濃

75、度分布</p><p>　　圖4-11 (YOZ，y=0，z=0.2)沿X軸各點甲醛濃度比較</p><p>　　圖4-8、4-9、4-10這些圖片是甲醛分布的等值線圖，圖上也很清楚的標注著每立方米甲醛濃度的數(shù)值。整體來看甲醛濃度都較高，相比Z=0.3截面來說高出了很多，也很容易理解，甲醛是從地板釋放的，在較低的高度濃度會比較高。在風速由1m/s增至2m/s甲醛濃度都是降低的。在風速增大至

76、3m/s時，上下兩邊的甲醛濃度都是在顯著增高的。圖4-11是在YOZ，y=0，z=0.2)沿X軸方向的線圖?？梢钥闯鲈陲L口位置，隨著風速的增大，帶動著甲醛濃度的降低。</p><p>　　通過上面的圖和分析可以知道在只有木質地板釋放甲醛的情況下，甲醛在三個不同的截面，不同風速下的濃度分布情況。</p><p>　　對于y=0截面來說，都是最底下濃度最高，往上面是呈逐漸遞減的趨勢。1米每秒風

77、速下甲醛濃度還是很高，在2米每秒風速下，甲醛濃度下降的比較厲害，在3米每秒風速下，甲醛濃度還是有稍微的下降，但是很不明顯。</p><p>　　對于z=-0.3截面來說，是人站立的高度1.7米，在以上三種風速下，甲醛濃度總體是比較低的，從圖4-4和圖4-6可以看出，甲醛濃度有比較集中的區(qū)域。在圖4-5中甲醛濃度是最低的，而且沒有明顯的集中區(qū)域。并且。在風速增大至3m/s時，上下兩邊的甲醛濃度都是在顯著增高的。&l

78、t;/p><p>　　對于z=0.2截面來說，在風速為1米每秒的情況下甲醛濃度大體來說都是較低，不過在最右邊區(qū)域有很明顯的集中。在風速為2米每秒和3米每秒的情況下，可以比較得出，在2米每秒風速下甲醛濃度較低。</p><p>　　綜合以上三個截面的情況來看，在只有木質地板釋放甲醛的情況下，2米每秒的風速為最佳的選擇。</p><p>　　4.3.5放置沙發(fā)立柜不同風速下

79、y=0截面甲醛濃度分布</p><p>　　圖4-12放置沙發(fā)立柜風速為1m/s 時y=0截面甲醛濃度分布</p><p>　　圖4-13放置沙發(fā)立柜風速為2m/s 時y=0截面甲醛濃度分布</p><p>　　圖4-14放置沙發(fā)立柜風速為3m/s 時y=0截面甲醛濃度分布</p><p>　　從以上圖中可以看出，甲醛濃度相比只有木質地板散發(fā)

80、甲醛的情況高很多，也很直觀的看出隨著風速的增大，甲醛濃度都是在逐漸下降的。同時也出現(xiàn)了右上角的濃度局部集中。</p><p>　　4.3.6放置沙發(fā)立柜不同風速下z=-0.3截面甲醛濃度分布</p><p>　　圖4-15放置沙發(fā)立柜風速為1m/s 時z=-0.3截面甲醛濃度分布</p><p>　　圖4-16放置沙發(fā)立柜風速為2m/s 時z=-0.3截面甲醛濃度分

81、布</p><p>　　圖4-17放置沙發(fā)立柜風速為3m/s 時z=-0.3截面甲醛濃度分布</p><p>　　圖4-18 (YOZ，y=0，z=–0.3)沿X軸各點甲醛濃度比較</p><p>　　圖4-15、4-16、4-17這些圖片是甲醛分布的等值線圖，圖上也很清楚的標注著每立方米甲醛濃度的數(shù)值。整體來看甲醛濃度分布很不均勻，從每立方米0.1mg到每立方米0

82、.9mg。在風速由1m/s增至2m/s再到3m/s甲醛濃度都是降低的。圖4-18是在YOZ，y=0，z=–0.3)沿X軸方向的線圖。我們可以看出在風口位置，隨著風速的增大，帶動著甲醛濃度的降低。</p><p>　　4.3.7放置沙發(fā)立柜不同風速下z=0.2截面甲醛濃度分布</p><p>　　圖4-19放置沙發(fā)立柜風速為1m/s 時z=0.2截面甲醛濃度分布</p><

83、;p>　　圖4-20放置沙發(fā)立柜風速為2m/s 時z=0.2截面甲醛濃度分布</p><p>　　圖4-21放置沙發(fā)立柜風速為3m/s 時z=0.2截面甲醛濃度分布</p><p>　　圖4-22 (YOZ，y=0，z=0.2)沿X軸各點甲醛濃度比較</p><p>　　圖4-22 (XOZ，x=0，z=-1)沿Y軸各點甲醛濃度比較</p>&l

84、t;p>　　從圖4-19、4-20、4-21整體來看甲醛濃度分布都較高。在風速由1m/s增至2m/s再到3m/s甲醛濃度都是降低的。圖4-18是在YOZ，y=0，z=–0.3)沿X軸方向的線圖。我們可以看出在風口位置，隨著風速的增大，帶動著甲醛濃度的降低。不過本文也截取了 圖4-22 (XOZ，x=0，z=-1)沿Y軸各點甲醛濃度比較，從圖中可以看到在風速為2m/s時，甲醛濃度時最低的。</p><p>　

85、　通過以上放置沙發(fā)立柜的房間中的數(shù)值模擬圖和分析可以看出，三個截面不同風速的條件下，甲醛濃度的分布情況。</p><p>　　比較三種截面三種風速的條件下，不難發(fā)現(xiàn)甲醛的濃度都是比只有木質地板釋放甲醛的情況要高得多，整體分布也更復雜。不管是哪個截面甲醛濃度都是隨著風速的增大而降低，可以得出結論，在放置沙發(fā)立柜的房間中，甲醛濃度較高，3m/s的風速是較好的選擇。</p><p><b&

86、gt;　　5結論</b></p><p><b>　　5.1結論</b></p><p>　　本文運用Fluent 軟件對房間里地板以及沙發(fā)立柜甲醛污染物的散發(fā)進行了數(shù)值模擬，首先建立物理模型，研究了各種不同風速條件下地板、沙發(fā)立柜散發(fā)甲醛的影響，確定了影響甲醛在室內擴散的一些因素，對今后室內氣流組織的研究和設計具有一定的指導意義。</p>

87、<p>　　全文得出的主要結論如下：</p><p>　?。?）在甲醛散發(fā)點比較少濃度相對來說低的情況下，風速在2米每秒為宜，如果甲醛散發(fā)點較多濃度比較高可以選擇3米每秒的風速，但也不宜過大。</p><p>　?。?）由于地板散發(fā)的甲醛向上蔓延的特征，在同種氣流組織形式下，采用過小的送風速度時，由于不足以對甲醛形成足夠強度的阻隔作用，致使甲醛氣體向較高平面擴散；而采用過大的送風

88、速度時，氣流的較大的擾動和攜帶作用，導致了甲醛氣體向上涌動，致使較高平面上甲醛濃度的增大，加大了對人體的危害性。因此，消除地板散發(fā)的污染物時，不宜采用過大或過小的送風速度。建議在2米每秒和3米每秒之間。</p><p><b>　　5.2展望</b></p><p>　?。?）本文對室內環(huán)境做了很多簡化，只考慮地板、沙發(fā)和立柜作為污染源，其實還有比如墻體粉刷擠也是污染

89、源。另外對污染物也做了很大程度的簡化，只考慮了甲醛氣體，而且控制了其密度不變、散發(fā)速度不變。而實際室內的污染物種類繁多、狀態(tài)變化大，相互之間擾動和作用力較大。在以后的研究中，可以把室內環(huán)境考慮得更為細致，再加入幾種污染物；同時，可以嘗試模擬污染物濃度隨時間的變化規(guī)律。</p><p>　?。?）為了便于觀察研究結果，本文忽略了考慮單居室，以后可以采用新型的幾何外形比較復雜的房間及室內布置來進行研究，比如兩室一廳和

90、三室一廳的居民房；而且，還有溫度的變化可以研究熱流對甲醛擴散的影響。 </p><p><b>　　參考文獻</b></p><p>　　[1] 夏元詢.化學物質毒物全書[M].上?？茖W技術文獻出版社,1991. </p><p>　　[2] 郭曉，張朝飛.試論甲醛對人體的危害及預防措施[J].河南農(nóng)業(yè)出版社

91、 2012，6(2):35-36</p><p>　　[3] 李延紅，蘇瑾，楊濱等.裝飾材料中甲醛對居室空氣污染和健康危害研究[M]. 勞動衛(wèi)生，2001，18(1):25-27</p><p>　　[4] 張會群，涂明揚.室內空氣中甲醛的危害及防治[J]廣西輕工業(yè)·資源與環(huán) 境,2008，2(2):

92、68-69</p><p>　　[5] 胡麗慧,葉茂,山桂云.學生寢室甲醛污染研究[J].生態(tài)·環(huán)保,2006,1:39</p><p>　　[6] 付小平，王遠成，姜作福. 上送下回通風方式稀釋室內污染物的數(shù)值模擬研究節(jié)能[J]. 2008，8:21-23.</p><p>　　[7] 吳猛，胡桂林，樊建人. 住宅室內空氣污染物排放的三維大渦數(shù)值模擬[J

93、]. 浙江大學學報（工學版），2006，40(6):987-989.</p><p>　　[8] 張俊杰，任雁秋. 室內污染物換氣過程的數(shù)值模擬[J]. 包頭鋼鐵學院學報，2006，25(3):266-268.</p><p>　　[9] 劉娣，湯廣發(fā)，趙福云. 空調臥室空氣環(huán)境的數(shù)值模擬[J]. 湖南大學學報，2006，3(33):28-32.</p><p>　

94、　[10] 張牛牛，薛敏. 室內污染物排放的數(shù)值模擬分析[J]. 江蘇工業(yè)學院學報，2009，21(2): 27-30.</p><p>　　[11] 馬飛，羅海亮，謝慧，等. 建筑室內污染物擴散模擬與控制[J]. 北京科技大學學報，2007，29(2): 16-20.</p><p>　　[12] 胡平放，蔡芬. 氣流組織形式對室內空氣環(huán)境影響的數(shù)值模擬[J]. 華中科技大學學報(城市科

95、學版)，2006，23(2): 28.</p><p>　　[13] 嚴玉梅. 撲捉人類健康的殺手) 甲醛 [J]. 廣東建材, 2006 ，3( 7): 43</p><p>　　[14] 趙彬，李先庭，彥啟森. 暖通空調氣流組織數(shù)值模擬的特殊性[J]. 暖通空調, 2004，34(11): 122.</p><p>　　[15] 李先庭，趙彬. 室內空氣流動數(shù)值

96、模擬[M]. 機械工業(yè)出版社，2009.</p><p>　　[16] 朱紅鈞，林元華，謝龍漢. FLUENT12流體分析及工程仿真[M]. 清華大學出版社，2011.</p><p>　　[17] 于勇，姜連田. FLUENT入門與進階教程[M].北京理工大學出版社，2008.</p><p>　　[18] 李進良，李承曦，胡仁喜. FLUENT6.3流場分析[M

97、]. 化學工業(yè)出版社，2009.</p><p>　　[19] 溫正. FLUENT流體計算應用教程[M]. 清華大學出版社，2013.</p><p>　　[20] GB/T 18883-2002.《室內空氣質量標準》.</p><p><b>　　致謝</b></p><p>　　本次論文是在**老師的悉心教導下