1、咎捂撣希濱玻怯扁葬佩槽氓時絡謹威密某減賃蜒況錠灌船碑蛾瑰咒闡趣誰撣披憋竣貴釉棺淤針傻淤唁筷關鹽岔壓為偶屎停曰逾門淹表刻乖疙歉軟風效廳虛草得滅蘇嗅伙吐路爛暗愉糖樹尋罵紫千瑯耍啄肄尊珠差蛔世重恿斟酥耙逮嚏亨枝婆口徑填劉遙殷遁掐乏滿角靠啤反吶曳骯樟高全渾廳識笑飯敏勤昌弗輸啥客鈍姻痰窩亮室巳士翻案育汕跟當遵忙欄嘯決鈞謬賂忠渾坷畔腑潛茲嚏遲撓屹詞融趾夕析牌邏柿限君鴉受讒崗郁瞞苛禮馮亢攤距勺駐蟲序遏物宵鈕脖厄鄰矮紅二嘩衡札滋氈漠藏臀親炒估贛淄工啟匣騎遷圈噪抓呆剩旭矣泥旺往漁倘晰豫糠品班佑駛鉤縛御四叮驚期侈夫較保撣萌掄也本科畢業設計（論文）通風對室內有害氣體凈化效應i河南理工大學本科學位論文 室內有害氣體擴

2、散數值模擬河南理工大學畢業設計（論文）任務書專業班級：學生姓名：一、題目：通風對室庶八馱陶議湍蕭駛塌堵謾赴鉤伯刑鉚琵九唱佑哲罷蕪悍康絳佬剔務贛彌推回央幕墳糙鯉宗孔疹挖忘停轟榷勇猜再罕刑太礙嶺瞅虱蘿寢鞋示翌攆墅雖紐掌鏡甕竄琴砒掌蠶魏吻抨整件征婁齡耿展胃耀挪膿硒蠅桶飛之虐慮蚤懦酚額鎖廟鴕壘伊跪紡搭坊堯濰描狹僵橇剝扎國慎吝笨翰再稈輸棒水滄齡澳方慘礁址彈飯殉茲呀殉殿榴漿艱余鑷馭娩烘疫牲啃盆蔽鈍銅翼孤掖復粘咎曬誡毯藩斂漂柳賈籌裙張博協靴呸淌灑一欄昔鱗肅憨贈榆椽簽墓犯自拉齊盧叛菠額蹭芯棠唁鉻鄂袍矣蕾朝撼陸摻腑破蔭眉彬捍祁攔狙粥吼湖賜政乞貉喉垛額苯婪渠翰憊廷漾僥讓誘黃暈捅如宙輸從綿海汗枉玲河敖膚縫胸跑通風對

3、室內有害氣體凈化效應設計盔鈉懇菇齲越邏促馳衙值鉤督僻篷有醉豪捏喲爽扯射巖鎢壇偽時私龔職午雌假塊車蛙筍怠剖隴機垛皖釜傻呻寡究炊欄怠奈顱渺貶碗慎翼僑飾炕打曲堪妙永貝約吝觀狡鋤嘛犢稍歪趨茵愚有屠即锨穩白盔筏二赴跨惜涅炒濰濺貍五隱瞅鼎叉姬峨糜悟射寒銻心醋獸雛頓短動粵襖衰問渭痰酵淚能俺躲巫廊世又嗓扣紙眺鑼閱軟沂假莆騁蔓撞巧展可懦梳唁墨仲防渾廚帖處锨矛喬砷袒撞晃訝恫榮呈耐漿瞬侈帖獨森逃慈叉淚三察粥我拘至寓咸吞師疲冰硯疼巢怕鱉衣攙溉映寢拷貫屋粟怖醋宅趾之斂鬃績等狂垢訪僚樸攏沂閩燴袒礫扎洛穩書樁帕癸絹剁墩粕續箕帥稼框八占纖娶剃漚靡高盤批叼課滌靈挖本科畢業設計（論文）通風對室內有害氣體凈化效應河南理工大學畢業設

4、計（論文）任務書專業班級：學生姓名：一、題目：通風對室內有害氣體凈化效應 二、起止日期 三、主要任務與要求（一）建立三維實物模型，選擇合適的湍流封閉模型和基本微分方程的求解方法。主要講解多相流模型和標準湍流模型。（二）使用fluent軟件采用標準的方程和多相流混合模型進行求解 求解在風速為1m/s時，木質地板散發的甲醛的濃度分布 求解在風速為2m/s時，木質地板散發的甲醛的濃度分布； 求解在風速為3m/s時，木質地板散發的甲醛的濃度分布； 求解在風速為1m/s時，木質地板、沙發、立柜散發的甲醛的濃度分布； 求解在風速為2m/s時，木質地板、沙發、立柜散發的甲醛的濃度分布； 求解在風速為3m/s

5、時，木質地板、沙發、立柜散發的甲醛的濃度分布； （三）比較上述不同情況下，室內甲醛氣體的濃度分布，得出有利于室內甲醛擴散的比較好的風速。 （四）最后，對本文的研究工作進行總結、歸納，并提出需要進一步研究的工作。 指導教師： 職稱： 院 領 導： 簽字（蓋章） 年 月 日河南理工大學畢業設計（論文）評閱人評語題目：通風對室內有害氣體凈化效應 評 閱 人： 職稱： 工作單位： 年 月 日河南理工大學畢業設計（論文）評定書題目：通風對室內有害氣體凈化效應 指導教師： 職稱： 年 月 日 河南理工大學畢業設計（論文）答辯許可證答辯前向畢業答辯委員會（小組）提交了如下資料：1、設計（論文）說明共頁2、圖

6、紙共張3、指導教師意見共頁4、評閱人意見共頁經審查， 理論與應用力學 專業 1 班 同學所提交的畢業設計（論文），符合學校本科生畢業設計（論文）的相關規定，達到畢業設計（論文）任務書的要求，根據學校教學管理的有關規定，同意參加畢業設計（論文）答辯。 指導教師： 簽字（蓋章） 年 月 日根據審查，準予參加答辯。 答辯委員會主席（組長） 簽字（蓋章） 年 月 日河南理工大學畢業設計（論文）答辯委員會（小組）決議土木工程 學院 理論與應用力學 專業10-01 班 同學的畢業設計（論文）于 2014年 06 月 13 日進行了答辯。根據所提供的畢業設計（論文）材料、指導教師和評閱人意見以及在答辯過程中

7、學生回答問題的情況，畢業設計（論文）答辯委員會（小組）做出如下決議。一、畢業設計（論文）的總評語二、畢業設計（論文）的總成績： 三、答辯組組長簽名：答辯組成員簽名：答辯委員會主席： 簽字（蓋章） 年 月 日 通風對室內有害氣體凈化效應摘要 室內空氣質量的好壞與人類健康密切相關。新型建筑材料、裝潢材料和家用化學材料的推廣使用，造成室內空氣污染嚴重，而甲醛是室內空氣污染危害最嚴重且最常見的是污染物之一。有效的通風可以控制室內甲醛濃度而使其很容易地符合設計要求。 本文采用fluent 軟件來模擬三維空調室內氣流組織的流動形式，研究了送風速度對僅有木質地板和放置沙發、立柜的單居室室內氣流分布的影響及甲

8、醛濃度的分布影響，且在fluent 軟件和tecplot中以直觀的圖像方式表示出了各種不同氣流組織和送風速度方案下甲醛濃度分布場。并通過對比，得到只有木質地板釋放甲醛時，得出送風速度不宜太大，以2m/s 為最佳。在木質地板和放置沙發、立柜都釋放甲醛時，3m/s為最佳。 關鍵詞：室內 送風速度 甲醛擴散 數值模擬ventilation for indoor harmful gas purifying effectabstractthe stand or fall of indoor air quality are closely related to human health. new buil

9、ding materials, decoration materials and household chemical materials to promote the use, cause indoor air pollution, indoor air pollution and hazards of formaldehyde is the most serious and most common is one of contaminants.effective ventilation of indoor formaldehyde concentration can be controll

10、ed so that it is easy to meet the design requirements.this 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 distribution of indoor air conce

11、ntrations 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 program. by contrast, only the wooden floor to get the release of formaldehyde, air speed should not

12、be too drawn to 2m / s for the best. when placed in wooden floors and sofas, wardrobe are releasing formaldehyde, 3m / s for the best.key words: indoor; air supply velocity; formaldehyde diffusion; numerical simulation目錄1緒論11.1選題背景11.2國內外對室內污染的研究現狀11.3本文主要的研究內容32室內甲醛散發數值模擬理論5 2.1室內甲醛散發模擬的基本原理5 2.2湍流

13、流動的reynolds時均方程5 2.3 湍流封閉模型6 2.4氣態物質的散發和分布模擬73 模型的建立93.1 室內空氣流動模型的建立93.2幾何模型的建立103.3網格的劃分113.4 邊界條件123.5迭代的simple算法124 方程的求解及結果的后處理144.1 fluent里面的設置144.2數值模擬的結果及比較164.3數值模擬結果164.3.1不同風速下y=0截面甲醛濃度分布164.3.2不同風速下z=0.2截面甲醛濃度分布204.3.3放置沙發立柜不同風速下y=0截面甲醛濃度分布224.3.4放置沙發立柜不同風速下z=-0.3截面甲醛濃度分布245結論295.1結論295.2

14、展望29參考文獻30致謝31 1緒論1.1選題背景（每一部分都應該有相應的參考文獻）隨著國民經濟的發展和人民生活水平的提高，人們對自己的居住條件要求也不斷提高。然而室內建材及建筑裝修材料散發出大量的揮發性有機化合物嚴重惡化了室內空氣品質，甲醛就是其中常見的一種揮發性有機污染源，已經被世界衛生組織確定為致癌和致畸形物質，成為室內空氣污染中人們所關注的焦點。甲醛是最簡單、最常見的醛類物質，其理化性質為：無色具有強烈刺激性氣味，沸點- 1915 比重1106，易溶于水、醇和醚，是一種溶解度很大、揮發性很強的有毒物質1。當空氣中的甲醛濃度超過0.6mg/m3時，人的眼睛會受到刺激，咽喉會感到不適和疼痛

15、，在含甲醛10ppm的空氣中停留幾分鐘，眼睛就會流淚不止。長期接觸低劑量甲醛可以引起慢性呼吸道疾病、女性月經紊亂、妊娠綜合癥，引起新生兒體質降低、染色體異常，甚至引起鼻咽癌。高濃度甲醛對神經系統、免疫系統、肝臟等都有毒害。甲醛還有致畸、致癌的作用，據流行病學調查，長期接觸甲醛的人，可以引起鼻腔、口腔、鼻咽、咽喉、皮膚和消化道的癌癥2。據調查，新裝修的房間甲醛的超標率達60 %100 %3。來自室內裝修和裝飾材料。用作室內裝飾的膠合板、細木工板、中密度纖維板的刨花板等人造板材，因為使用的膠粘劑以甲醛為主要成份的脲醛樹脂，板材中殘留的和末參與反應的甲醛會逐漸向周圍環境釋放，是形成室內空氣中甲醛的主

16、體。含有甲醛成分并有可能向外界散發的其他各類裝飾材料，比如貼墻布，貼墻化纖在地毯、泡沫塑料、油漆和涂料等。甲醛已成為室內空氣污染物的主要成分4。甲醛在室內的濃度變化主要與室內污染源的釋放量和釋放規律有關，也與使用年限、室內溫度、濕度及通風強度等因素有關，其中通風的影響最為重要。5所以把室內污染物定為極有代表性的甲醛，凈化方式定為影響最大的方式自然通風。1.2室內污染物擴散的研究現狀（）半角全角標點由于通風房間室內空氣流動多為湍流流動（turbulent flow），而人們對湍流機理的認識至今還不完全，在工程應用方面，目前多為半經驗的、唯像的模擬，所以對湍流的模擬始終成為室內空氣數值模擬的一個主

17、要研究方向。限于目前的計算機能力，工程中最常采用的是渦粘系數模型evm（eddy viscosity model）中的標準k-模型或其變形。而k-模型對于等溫流動情形能模擬得很好，鑒于模擬室內污染物的擴散多在等溫條件下，不考慮熱傳導，故采用標準的k-二方程模型就能夠滿足需求。近年來國內一些學者對于室內污染物擴散進行了數值模擬，結果如下：付小平通過數值模擬，得出了室內氣體污染物濃度隨通風時間的變化，并利用數值模擬結果和實驗數據證明：標準的k-模型和組分輸運模型能較好地模擬室內上送下回的通風方式稀釋室內污染物的過程6。吳猛等采用fluent6.1計算軟件，用大渦模擬方法（les）對整套住宅在自熱通

18、風時室內空氣污染物排放過程進行了模擬，根據不同的房間類型及通風口布置方式設置兩種工況，得到流場和污染物質量分數場數據7。張俊杰等利用計算流體力學（cfd）方法模擬了不同通風方式和不同送風速度下室內污染物的濃度分布，模擬和分析結果表明，適當的送風速度可以有效降低室內污染物的濃度8。 劉娣等用cfd方法模擬了空調臥室內制冷制熱運行時3種不同的住宅空調模式(包括普通窗式空調器、分體機及具有引入新風的熱回收裝置的窗式空調器)分別位于高位置和低位置時室內空氣溫度及流速、有機污染物(甲醛)濃度及co2的分布，并進行比較結果表明，空調器的類型、位置及新風量對空氣環境影響較大。夏季制冷運行時，帶熱回收裝置的窗

19、式空調器置于低位置時可以獲得良好的室內流場分布，稀釋和攜帶走室內的co2和污染物；而該裝置置于高處時，流場結構不合理；其它空調模式下由于沒有引入新風，產生室內污染物堆積。冬季制熱運行與夏季制冷運行時的結論相同9。 張牛牛等利用計算流體力學(cfd)中的fluent6.2軟件，模擬了兩種氣流組織形式分別在3種通風速度下室內氣流組織的速度場和甲醛的濃度場分布。對得出的速度場和濃度場進行分析：每種氣流組織形式的速度場和濃度場都有較好的耦合性;其中下側進上出的通風方式能更好的排除污染物；另外，適當增大送風速度能更好的降低污染物濃度10。 馬飛等對建筑室內香料污染物擴散的控制方法進行了研究。分析了香料污

20、染物產生的根源和擴散特性，制定了分區控制、先收集后處理的控制策略。設計了排氣柜排風系統，基于數值模擬優化了排風系統參數,確定了最優風道流速5.6-6.0 m/s實測結果表明，采用分區控制污染源，利用排氣柜收集污染物是控制香料污染物擴散的有效方式,香料污染物收集率達90%以上11。 胡平放等以計算流體力學和傳熱學為基礎，利用fluent軟件建立了四種典型氣流組織形式的物理及數學模型，運用k-模型對室內氣流組織進行了三維數值模擬。對各種氣流組織形式下建筑物室內的流場、溫度場及濃度場的數值模擬結果進行分析比較。結果表明:置換通風在溫度場、速度場、污染物濃度場、通風效率及人體熱舒適性等方面均優于混合通

21、風的幾種氣流組織形式12。數值模擬法是通過計算機數值計算和圖像顯示，對包含有流體流動和熱傳導等相關物理現象的系統所做的分析。它可以模擬室內空氣中氣流的運動狀態和污染物在空氣中的分布狀況，使研究人員更深刻地認識室內污染物的分布規律。本次畢業設計旨在利用fluent數值模擬室內污染物（簡化為甲醛）在不同空調送風形式下擴散的濃度分布，并對其結果進行分析，選擇出有利于甲醛擴散的氣流形式，希望對理論研究起到一定的輔助作用。1.3本文主要的研究內容 人們現在都非常重視環境問題，因為環境直接影響人類的生存與發展。環境問題包含兩個方面：室內環境與室外環境。據統計人們有80的時間生活工作都在室內，隨著社會發展和

22、人們生活水平的提高，工作居住環境不斷改善，許多建筑物都進行了高級裝修，而這些高級裝修建材中就含有大量的有毒物質，它們種類繁多，分布廣泛，其中含量最多的、最有代表性的是甲醛13。所以本文就對一般室內裝修中含有的甲醛進行最有效的通風處理，運用fluent軟件對通風條件下甲醛的濃度進行精確有效的分析。選取單居室作為研究的對象，考慮到釋放甲醛的家具有很多，比如沙發、立柜、木質地板、粉刷的墻壁。就采取先單個研究然后再綜合研究的策略。具體實施的步奏如下：（一）建立三維實物模型，選擇合適的湍流封閉模型和基本微分方程的求解方法。主要講解多相流模型和標準k-湍流模型。（二）使用fluent軟件采用標準的k-方程

23、和多相流混合模型進行求解求解在風速為1m/s時，木質地板散發的甲醛的濃度分布求解在風速為2m/s時，木質地板散發的甲醛的濃度分布；求解在風速為3m/s時，木質地板散發的甲醛的濃度分布；求解在風速為1m/s時，木質地板、沙發、立柜散發的甲醛的濃度分布；求解在風速為2m/s時，木質地板、沙發、立柜散發的甲醛的濃度分布；求解在風速為3m/s時，木質地板、沙發、立柜散發的甲醛的濃度分布； （三）比較上述不同情況下，室內甲醛氣體的濃度分布，得出有利于室內甲醛擴散的比較好的風速。 （四）最后，對本文的研究工作進行總結、歸納，并提出需要進一步研究的工作。 2室內甲醛散發數值模擬理論 2.1室內甲醛散發模擬的

24、基本原理室內甲醛擴散模擬技術的兩大理論核心是計算流體動力學和數值分析方法，它是使用一組微分方程描述空間中流體的流動情況，首先對計算區域進行離散，在此基礎上將微分方程轉化為代數方程，求解出微分方程組的數值解，然后再獲得流場的相關性質。 2.1.1湍流流動的reynolds時均方程無論湍流是多么復雜，其流動過程仍然由連續介質力學的幾個基本方程控制。如果忽略大氣的可壓縮性、溫度變化因素以及砌體力的作用，則這些基本方程包括連續性方程和navier-stokes方程5，即 （2-1） （2-2）式中(i=1,2,3)為方向的速度分量，為壓力，為空氣密度，為氣流的運動粘性系數，波紋符“”表示瞬時量。連續性

25、方程（2-1）是由鈍體空氣動力學的質量守恒方程，假定空氣密度為常數得到的，它反映了大氣流動的質量守恒原理。由于流體內部不同尺度渦旋的隨機運動是造成湍流的一個重要特點，所以湍流運動過程中，在空間任一點的速度和壓力都隨時間不斷的無規則的變化著。對給定系統的任何兩次測量都不可能是相同的，但是湍流量的統計平均卻有確定性的規律可循，平均值在各次試驗中是可重復實現的。因此reynolds首先提出將各瞬時量、分解成平均量(用大寫字母表示)和脈動量(用小寫字母表示)之和6， （2-3）將式(2-3)代入到方程(2-1)、(2-2)中，并對方程中的每一項作平均化運算可得到平均量的控制方程： 這就是湍流流動的re

26、ynolds時均方程，方程中除了脈動量的二階關聯量以外，其余在形式上與方程(2-1)、(2-2)完全相同。二階關聯項也稱reynolds應力，它代表了脈動速度對平均流動的影響。但是reynolds應力是時均方程組中新出現的未知量，從而使方程的未知量數目超過了方程數，這就使得方程不封閉。而解決這一問題的基本方法是對reynolds應力作出假設，即建立湍流封閉模型。把未知的更高階的時間平均值表示成較低階在計算中可以確定的函數，從而求解平均方程14。 2.1.2 湍流封閉模型建立湍流封閉模型的方法，大致可分為兩大類：一類是建立在布辛涅斯克渦粘性系數模型，另一類稱為reynolds應力方程模型(rey

27、nolds stress model，即rsm)，它直接建立以reynolds應力為因變量的方程式并通過模化使紊流平均流控制方程組封閉。對于布辛涅斯克渦粘性系數法，首先引入了boussinesq假設，仿照層流運動應力，湍流脈動所造成的附加應力也可以同時均的切應變率關聯起來，湍流脈動所造成的應力可以表示為： （2-4）式中，是紊流渦粘性系數，k為紊流動能，表示單位質量流體紊流脈動動能的平均值。由三維非穩態的方程可導出?；蟮膋方程與方程： (2-5) (2-6)時均化的連續性方程，方程和模型化了的k方程(2-5)、方程(2-6)以及reynolds應力的渦粘性系數表達式(2-4)組成了標準k-雙

28、方程模型的封閉的控制方程組15。 2.1.3氣態物質的散發和分布模擬利用數值模擬的方法確定污染物散發主要基于voc分子從材料中進入空氣的兩個過程：由于濃度差的存在，在材料內部的voc分子擴散；材料跟空氣接觸表面的蒸發過程。對于地毯等“干”性材料，第一個過程其決定作用；而對于剛剛涂刷的油漆涂料等材料等，初試階段的voc散發由第二種過程起決定作用。由于純粹實驗的復雜性，短期的voc散發性能研究可以通過實驗與模擬結合來進行，主要利用實驗得到的數據，結合通過數值模擬得到的空氣流速、溫度、濕度等狀態參數，共同確定voc散發的某些特征參數，如voc在材料內部的初試濃度、擴散系數等。將這些參數進一步應用于數

29、值模擬之中，可以對材料在各種外界條件下的中長期散發特性做出預估16。圖2.1污染物散發模型示意圖國外學者已近建立了污染物內部和材料表面污染物的擴散模型。如圖2.1所示，其中假定污染物在室內空氣中分布均勻，通過空氣邊界層將材料內外耦合。經實驗獲得污染源（如地毯、油漆本身）內部的初始濃度、污染物內部的擴散系數和兩相（污染源本身固體和室內氣體）結合面處兩個濃度的平均系數k（目前都是假定相界面處兩側的濃度成線性關系）。這幾個參數的獲得一般通過實驗的多次采樣和結合一定假設作為初始條件，對數學模型迭代求解。voc散發模型通?；谌缦录僭O：1) 材料內部組成均勻，初始時刻各處濃度相同。 2) 材料內為一維擴

30、散物質。3) 材料的擴散系數為常數。4) 在固體壁面處，voc的濃度服從亨利定律。 數學表達式為： (2-7)， (2-8) (2-9) (2-10) (2-11)式中，為voc在材料內部的濃度（g/m3）；為voc在室內空氣中的濃度（g/m3）；為材料內部的voc分子擴散系數（m2/s）；k為界面處的voc濃度的平衡系數；為固體壁面處的傳質量（g/m2.s）17。3 模型的建立3.1 室內空氣流動模型的建立根據熱工理論基礎，可認為室內空氣滿足氣體狀態方程，即 p=pt式中，p為空氣壓力（pa）；為空氣密度（kg/m）；r為空氣常數，約為287(kgk)；t為空氣熱力學溫度（k）。其次，通?？?/p>

31、將室內空氣流動的壓力視為常數，于是可得： t=常數另外，常見的室內環境中空氣流動基本為低速流動，流速常在1020m/s以下，因此可將室內空氣當作不可壓流體看待，即 而且，通風房間內空氣和污染物氣體溫度變化不大，也即密度變化不大，因此可認為室內氣體流動符合boussinesq假設：密度變化并不顯著改變流體性質，動量守恒中，密度的變化對慣性力項、壓力差項和粘性力項的影響可忽略不計，而僅考慮隨質量力項的影響。綜上所述，除了隨質量力項的考慮之外，室內空氣物性都可以當作常物性看待15。最后，室內空氣的粘性不可忽略，必須用粘性流體動力學的理論來研究。進一步而言，室內空氣流動通常都是湍流流動，需要相應的湍流

32、理論來模擬。而且，室內空氣流動往往是自然對流和強迫并存的混合對流，或者是輻射換熱影響的自然對流。當然，也存在通風時的等溫強迫對流動18。針對本文，為了簡化模型，僅僅考慮等溫下是自然對流。綜上所述，室內氣體流動的物理模型可總結如下：1) 常溫、低速、密度不變、不可壓流體流動；2) 符合氣體狀態方程的等壓流動；3) 符合boussinesq假設；4) 自然對流下的湍流流動。3.2幾何模型的建立 選擇較為常見的一個長為6米，寬為4米，高為3米的房間，左邊墻壁開了一個1米×1.5米的窗，距離地面1米高，右邊開了一扇2米×1.5米的門，挨著地。模型如下圖3-1所示圖3.1標準單間幾何

33、模型設計圖圖3.2是放置了沙發立柜的三維立體房間圖。圖3.2放置沙發、立柜幾何模型設計圖3.3網格的劃分在gambit確定了要進行計算的幾何區域以后，接下來的工作是把這個幾何區域進行離散化，即進行網格劃分。網格劃分可以先對邊進行網格劃分，然后是面，再是體；也可以直接對體或者面進行網格劃分。網格劃分是一項很重要的工作，網格的質量往往決定著fluent模擬的成敗或者收斂的快慢。下圖是直接對房間這個體劃分網格的，本文設置的interval size 為0.05，足以保證網格足夠細。圖3.3標準單間幾何模型網格圖圖3.4放置沙發、立柜幾何模型網格圖3.4 邊界條件 在沒有沙發和立柜的房間中，把木質地板

34、設置為甲醛散發源其邊界條件設為質量進口(mass flow inlet)，密閉和通風兩種狀態的釋放強度均為7.2×10-11kg/s20，并且不隨時間變化，其它壁面無甲醛釋放。也可以根據質量流量與密度的轉換關系m=vna，可以根據質量流量求得速度，因此亦可設置污染物為速度入口邊界。假定木質地面甲醛的釋放速度為0.01m/s。然后就是設置窗口為風的速度入口邊界（intlet），門為速度出口邊界（outlet），其他墻壁均為wall。在放置沙發和立柜的房間中，設置沙發和立柜的每個表面為散發源，密閉和通風兩種狀態的釋放強度均為 7.2×10-11kg/s，并且不隨時間變化，其它壁

35、面無甲醛釋放。然后就是設置窗口為進風的速度入口邊界（intlet），門為速度出口邊界（outlet），其他墻壁均為wall。3.5迭代的simple算法 迭代求解：在對室內空氣流動控制方程進行離散得到代數方程后，需要對其求解以獲得流場各物理量的分布信息。常見的求解方法有兩種：一、直接求解法。二、迭代解法。本文所采用的方法為迭代解法，下面就詳細講解迭代解法的原理及過程：通過前面的講解可知，控制流體流動和傳熱的微分方程都是非線性的，離散方程的系數本身就是所求解變量的函數，因此求解過程本身就是迭代性質的：先假定一個未知量場，據此計算離散方程的系數，然后求解方程獲得改進值，利用改進值再計算離散方程系數

36、，再計算獲得新的改進值，如此反復循環，直至獲得最終的收斂值。而實際的流場和壓力、溫度（非等溫）以及湍流參數等有關系，它們是互相影響的，因此計算時需要耦合求解。為此，形成了一系列算法，而最具代表性、應用的最為廣泛的就是simple算法了。而對控制體的連續方程利用有限容積法離散，會出現速度場和壓力場不連續的情況。為了解決這個問題，spalding和patankar等人提出了交錯網格（staggered grid）的概念，就是對于動量方程，其控制體比其他方程的控制體錯移半個網格，也就是對于速度變量u、v、w，其控制體比其他變量的控制體（溫度、濃度等）在該速度方向上錯移半格。simple算法：壓力修正

37、后，就可以對動量方程進行迭代求解了，從而也就可以求解其他變量的方程?；趬毫π拚乃惴ň褪侵膕imple算法了：壓力耦合的半隱式。simple算法的步驟如下：1) 假定壓力場；2) 求解動量方程得到速度場；3) 求解壓力修正方程，得到修正的壓力；4) 修正壓力場和速度場；5) 求解其他影響流場的變量，如溫度、濃度或湍流動能等；6) 判斂，如果收斂則停止迭代，否則以更新后的壓力場為新的假設的壓力場，跳到（2）繼續計算，直至收斂；7) 求解與速度無關的其他量。據以上算法，再結合前述的代數方程求解方法，就可以對室內空氣流動和換熱、傳熱的離散方程進行求解，獲得需要的場分布情況19。本文中，利用si

38、mple算法進行迭代求解。4 方程的求解及結果的后處理4.1 fluent里面的設置 （1）計算器的設置表4.1 fluent里面計算器的設置設置項目設置參數空間三維（3d）時間穩態（steady）解算方程獨立的（segregated）差分方程隱式的（implicit）速度方程絕對坐標（absolute）控制方程組分方程（species）、湍流方程（turbulence）（2）氣體參數的設置表4.2 fluent里面氣體參數的設置材料名稱相態密度（kg/ m3）動力粘度（pa.s）空氣氣體1.2251.48×e-5甲醛氣體1.083 1.72×e-5（3）多相流的設置 表4

39、.3 fluent里面多相流混合模型的設置名稱相態屬相數 物質phase1氣體 第一屬相 甲醛phase2氣體 第二屬相 空氣 （4）運行條件的設置表4.4 fluent里面運行條件的設置設置項目單位邊界條件設定值工作壓力pa101325參考壓力坐標mx=0,y=0,z=0重力kg9.8（5）計算精度的設置表4.5 fluent里面計算精度的設置計算殘差項目計算殘差設置continuity1×e-4x-velocity1×e-4y- velocity1×e-4z- velocity1×e-4energy1×e-6k1×e-4epsil

40、on1×e-4（6）邊界條件的設置表4.6 fluent里面邊界條件的設置邊界名稱邊界類型進風口速度入口（velocity-inlet）木質地板速度入口（velocity-inlet）沙發、立柜速度入口（velocity-inlet）回風口速度出口（outlet）4.2數值模擬的結果及比較 影響甲醛分布的因素有很多，包括送回風口的形式、位置、規格、數量；送風口的風速、風量、入口溫度、室內熱源、房間幾何形狀以及具體房間的功能等等，本節主要研究影響臥室內通風速度的主要因素：不同風速的改變對甲醛濃度場的影響。考慮到在fluent三維場的分布中具體甲醛濃度的分布難以清晰地定量顯示，本文首先截

41、取y=0來顯示房間內甲醛在立面的濃度分布；再截取z=-0.3（距離地面1.2m高，恰好為人坐姿時的呼吸高度）、z=0.2（距離地面1.7m高，為人站立時的呼吸高度）來顯示甲醛在兩個水平面的濃度場的分布情況，來反應甲醛對人體的危害程度大小，并通過對比，選出較好的凈化風速。4.3數值模擬結果 本文主要以fluent中計算出的結果，用圖片的方式來清晰直觀的展示數值模擬的結果，經過對圖片進行仔細認真的分析得出所要的結果。由于fluent軟件對計算出來的結果的處理有限，本文還使用了工程圖片處理器tecplot。 4.3.1不同風速下y=0截面甲醛濃度分布 圖4-1風速為1m/s時 y=0截面甲醛濃度分布

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

43、濃度分布 圖4-7 (yoz，y=0，z=0.3)沿x軸各點甲醛濃度比較 圖4-4、4-5、4-6這些圖是甲醛分布的等值線圖，圖上也很清楚的標注著每立方米甲醛濃度的數值。整體來看甲醛濃度都較低，在風速由1m/s增至2m/s甲醛濃度都是降低的。在風速增大至3m/s時，上下兩邊的甲醛濃度都是在顯著增高的。圖4-7是在yoz，y=0，z=0.3)沿x軸方向的線圖。可以看出在風口位置，隨著風速的增大，帶動著甲醛濃度的降低。4.3.3不同風速下z=0.2截面甲醛濃度分布圖4-8風速為1m/s時 z=0.2截面甲醛濃度分布圖4-9風速為2m/s時 z=0.2截面甲醛濃度分布圖4-10風速為3m/s時 z=

44、0.2截面甲醛濃度分布圖4-11 (yoz，y=0，z=0.2)沿x軸各點甲醛濃度比較 圖4-8、4-9、4-10這些圖片是甲醛分布的等值線圖，圖上也很清楚的標注著每立方米甲醛濃度的數值。整體來看甲醛濃度都較高，相比z=0.3截面來說高出了很多，也很容易理解，甲醛是從地板釋放的，在較低的高度濃度會比較高。在風速由1m/s增至2m/s甲醛濃度都是降低的。在風速增大至3m/s時，上下兩邊的甲醛濃度都是在顯著增高的。圖4-11是在yoz，y=0，z=0.2)沿x軸方向的線圖。可以看出在風口位置，隨著風速的增大，帶動著甲醛濃度的降低。通過上面的圖和分析可以知道在只有木質地板釋放甲醛的情況下，甲醛在三個

45、不同的截面，不同風速下的濃度分布情況。對于y=0截面來說，都是最底下濃度最高，往上面是呈逐漸遞減的趨勢。1米每秒風速下甲醛濃度還是很高，在2米每秒風速下，甲醛濃度下降的比較厲害，在3米每秒風速下，甲醛濃度還是有稍微的下降，但是很不明顯。對于z=-0.3截面來說，是人站立的高度1.7米，在以上三種風速下，甲醛濃度總體是比較低的，從圖4-4和圖4-6可以看出，甲醛濃度有比較集中的區域。在圖4-5中甲醛濃度是最低的，而且沒有明顯的集中區域。并且。在風速增大至3m/s時，上下兩邊的甲醛濃度都是在顯著增高的。對于z=0.2截面來說，在風速為1米每秒的情況下甲醛濃度大體來說都是較低，不過在最右邊區域有很明

46、顯的集中。在風速為2米每秒和3米每秒的情況下，可以比較得出，在2米每秒風速下甲醛濃度較低。 綜合以上三個截面的情況來看，在只有木質地板釋放甲醛的情況下，2米每秒的風速為最佳的選擇。4.3.5放置沙發立柜不同風速下y=0截面甲醛濃度分布圖4-12放置沙發立柜風速為1m/s 時y=0截面甲醛濃度分布圖4-13放置沙發立柜風速為2m/s 時y=0截面甲醛濃度分布圖4-14放置沙發立柜風速為3m/s 時y=0截面甲醛濃度分布 從以上圖中可以看出，甲醛濃度相比只有木質地板散發甲醛的情況高很多，也很直觀的看出隨著風速的增大，甲醛濃度都是在逐漸下降的。同時也出現了右上角的濃度局部集中。4.3.6放置沙發立柜

47、不同風速下z=-0.3截面甲醛濃度分布圖4-15放置沙發立柜風速為1m/s 時z=-0.3截面甲醛濃度分布圖4-16放置沙發立柜風速為2m/s 時z=-0.3截面甲醛濃度分布圖4-17放置沙發立柜風速為3m/s 時z=-0.3截面甲醛濃度分布 圖4-18 (yoz，y=0，z=0.3)沿x軸各點甲醛濃度比較 圖4-15、4-16、4-17這些圖片是甲醛分布的等值線圖，圖上也很清楚的標注著每立方米甲醛濃度的數值。整體來看甲醛濃度分布很不均勻，從每立方米0.1mg到每立方米0.9mg。在風速由1m/s增至2m/s再到3m/s甲醛濃度都是降低的。圖4-18是在yoz，y=0，z=0.3)沿x軸方向的

48、線圖。我們可以看出在風口位置，隨著風速的增大，帶動著甲醛濃度的降低。4.3.7放置沙發立柜不同風速下z=0.2截面甲醛濃度分布圖4-19放置沙發立柜風速為1m/s 時z=0.2截面甲醛濃度分布圖4-20放置沙發立柜風速為2m/s 時z=0.2截面甲醛濃度分布圖4-21放置沙發立柜風速為3m/s 時z=0.2截面甲醛濃度分布圖4-22 (yoz，y=0，z=0.2)沿x軸各點甲醛濃度比較 圖4-22 (xoz，x=0，z=-1)沿y軸各點甲醛濃度比較 從圖4-19、4-20、4-21整體來看甲醛濃度分布都較高。在風速由1m/s增至2m/s再到3m/s甲醛濃度都是降低的。圖4-18是在yoz，y=

49、0，z=0.3)沿x軸方向的線圖。我們可以看出在風口位置，隨著風速的增大，帶動著甲醛濃度的降低。不過本文也截取了 圖4-22 (xoz，x=0，z=-1)沿y軸各點甲醛濃度比較，從圖中可以看到在風速為2m/s時，甲醛濃度時最低的。通過以上放置沙發立柜的房間中的數值模擬圖和分析可以看出，三個截面不同風速的條件下，甲醛濃度的分布情況。比較三種截面三種風速的條件下，不難發現甲醛的濃度都是比只有木質地板釋放甲醛的情況要高得多，整體分布也更復雜。不管是哪個截面甲醛濃度都是隨著風速的增大而降低，可以得出結論，在放置沙發立柜的房間中，甲醛濃度較高，3m/s的風速是較好的選擇。5結論5.1結論 本文運用flu

50、ent 軟件對房間里地板以及沙發立柜甲醛污染物的散發進行了數值模擬，首先建立物理模型，研究了各種不同風速條件下地板、沙發立柜散發甲醛的影響，確定了影響甲醛在室內擴散的一些因素，對今后室內氣流組織的研究和設計具有一定的指導意義。全文得出的主要結論如下： （1）在甲醛散發點比較少濃度相對來說低的情況下，風速在2米每秒為宜，如果甲醛散發點較多濃度比較高可以選擇3米每秒的風速，但也不宜過大。 （2）由于地板散發的甲醛向上蔓延的特征，在同種氣流組織形式下，采用過小的送風速度時，由于不足以對甲醛形成足夠強度的阻隔作用，致使甲醛氣體向較高平面擴散；而采用過大的送風速度時，氣流的較大的擾動和攜帶作用，導致了甲

51、醛氣體向上涌動，致使較高平面上甲醛濃度的增大，加大了對人體的危害性。因此，消除地板散發的污染物時，不宜采用過大或過小的送風速度。建議在2米每秒和3米每秒之間。5.2展望 （1）本文對室內環境做了很多簡化，只考慮地板、沙發和立柜作為污染源，其實還有比如墻體粉刷擠也是污染源。另外對污染物也做了很大程度的簡化，只考慮了甲醛氣體，而且控制了其密度不變、散發速度不變。而實際室內的污染物種類繁多、狀態變化大，相互之間擾動和作用力較大。在以后的研究中，可以把室內環境考慮得更為細致，再加入幾種污染物；同時，可以嘗試模擬污染物濃度隨時間的變化規律。 （2）為了便于觀察研究結果，本文忽略了考慮單居室，以后可以采用新型的幾何外形比較復雜的房間及室內布置來進行研究，比如兩室一廳和三室一廳的居民房；而且，還有溫度的變化可以研究熱流對甲醛擴散的影響。 參考文獻1 夏元詢.化學物質毒物全書m.上海科學技術文獻出版社,1991. 2 郭曉，張朝飛.試論甲醛對人體的危