壁面溫度對(duì)室內(nèi)空氣流動(dòng)過(guò)程的影響摘要：為了分析在進(jìn)風(fēng)位置、進(jìn)風(fēng)溫度以及射流方向不同的條件下，房間壁溫對(duì)室內(nèi)空氣流動(dòng)的影響，利用CFD計(jì)算、流動(dòng)可視化和定量測(cè)量相結(jié)合的辦法研究了不同工況下室內(nèi)空氣的流動(dòng)形態(tài)、速度以及模型整體換熱情況。結(jié)果表明：進(jìn)風(fēng)溫度變化時(shí)，壁面溫度對(duì)水平送風(fēng)方式下的空氣流動(dòng)形態(tài)影響較大，對(duì)垂直送風(fēng)方式影響較小；進(jìn)風(fēng)溫度的降低引起水平進(jìn)風(fēng)方式的自由射流區(qū)縮短，右側(cè)垂直進(jìn)風(fēng)方式下射流向右側(cè)壁面偏轉(zhuǎn)，同時(shí)浮升力引起的自然對(duì)流作用增強(qiáng)，換熱溫差加大，增加了圍護(hù)結(jié)構(gòu)漏熱量。關(guān)鍵詞：空調(diào)通風(fēng)；流動(dòng)形態(tài)；CFD模擬；可視化實(shí)驗(yàn)通風(fēng)空調(diào)房間內(nèi)的空氣溫度、流速以及污染物濃度等對(duì)建筑物能耗、室內(nèi)空氣品質(zhì)和人體健康至關(guān)重要。自20世紀(jì)40年代以來(lái)，研究人員運(yùn)用不同方法來(lái)預(yù)測(cè)機(jī)械通風(fēng)房間內(nèi)的氣流組織結(jié)構(gòu)[1-6]。ChenQinyan[3,4]利用CFD模擬研究了不同的通風(fēng)方式、風(fēng)口結(jié)構(gòu)對(duì)房間內(nèi)空氣流動(dòng)的影響；Fitzner[5]通過(guò)原型實(shí)驗(yàn)得出室內(nèi)氣流主要受風(fēng)口結(jié)構(gòu)、位置，空氣流速以及房間熱負(fù)荷的影響，而排風(fēng)口位置對(duì)其影響較??；Zhang[6]實(shí)驗(yàn)研究了水平送風(fēng)條件下，不同進(jìn)口風(fēng)速、溫度以及房間熱負(fù)荷對(duì)室內(nèi)空氣速度和湍流特性的影響。然而，在進(jìn)風(fēng)口位置、進(jìn)風(fēng)射流方向以及進(jìn)風(fēng)溫度不同的條件下，房間壁面溫度對(duì)室內(nèi)空氣流動(dòng)形態(tài)的影響規(guī)律未作具體研究。本文通過(guò)CFD模擬和可視化模型實(shí)驗(yàn)相結(jié)合的辦法研究了不同的風(fēng)口位置、射流方向以及進(jìn)風(fēng)溫度下，房間壁面溫度對(duì)室內(nèi)空氣流動(dòng)過(guò)程的影響。1實(shí)驗(yàn)?zāi)Ｐ图皵?shù)值模擬方法圖1實(shí)驗(yàn)?zāi)Ｐ徒Y(jié)構(gòu)示意意圖設(shè)定如圖1所示的坐標(biāo)系統(tǒng)，模型的在x、y、z方向的長(zhǎng)度分別為1.2m、0.6m和0.8m，為了便于進(jìn)行流場(chǎng)可視化研究，模型主體采用有機(jī)玻璃加工?？紤]到大部分房間只有一面墻為外墻，其余墻面都為內(nèi)墻，且內(nèi)墻可近似認(rèn)為是絕熱邊界條件，所以在本實(shí)驗(yàn)中只將左壁面（x=0m）設(shè)為加熱壁（模擬外墻），為了維持恒壁溫條件，提高溫度均勻性，設(shè)計(jì)了由薄膜加熱器和厚4mm導(dǎo)熱鋁板組成的復(fù)合加熱板，通過(guò)調(diào)節(jié)輸入電壓可以控制實(shí)驗(yàn)所需的壁溫，而模型的其他外表面通過(guò)粘貼厚度為100mm的保溫材料來(lái)模擬絕熱邊界條件。本文討論A、B兩種不同位置的進(jìn)風(fēng)方式，風(fēng)口尺寸都為240mm×24mm，且以z=0.4m的平面前后對(duì)稱。其中進(jìn)風(fēng)口A設(shè)在x=1.2m的平面，水平中心線距模型頂面0.1m；進(jìn)風(fēng)口B設(shè)置在y=0.6m的頂面，右沿距模型x=1.2m的平面0.1m；排風(fēng)口設(shè)在y=0.6m的頂面，左沿距x=0m的平面0.1m。通常情況下房間內(nèi)的平均流速小于0.15m/s，只有射流區(qū)域的流速達(dá)到相對(duì)容易進(jìn)行測(cè)量分辨的0.4m/s以上。為此，房間內(nèi)部的風(fēng)速、溫度測(cè)量點(diǎn)主要布置在射流方向上。對(duì)于A型進(jìn)風(fēng)方式，7個(gè)測(cè)量點(diǎn)布置在y=0.3m、z=0.4m的直線上，其中1點(diǎn)、7點(diǎn)x坐標(biāo)分別為1.088m、0.112m，其余5點(diǎn)x坐標(biāo)在1點(diǎn)和7點(diǎn)之間均勻布置；對(duì)于B型進(jìn)風(fēng)方式，5個(gè)測(cè)量點(diǎn)均勻布置在x=1.088m、z=0.4m的直線上。表1給出本實(shí)驗(yàn)的研究工況，加熱表面的溫度為40℃、風(fēng)量為0.0078m3/s。表1實(shí)驗(yàn)研究工況工況序號(hào)ⅠⅡⅢⅣ風(fēng)口位置AABB進(jìn)風(fēng)溫度/℃2214.52214.5本文采用FLUENT6.0商業(yè)軟件進(jìn)行計(jì)算，入流湍動(dòng)能kin為來(lái)流平均動(dòng)能的1％，湍流耗散εin可按照入口湍流雷諾數(shù)Ret以及湍流特征長(zhǎng)度之間的關(guān)系來(lái)選取[7]。利用SIMPLEC算法進(jìn)行壓力——速度解耦，PRESTO！格式處理壓力項(xiàng)，QUICK離散格式處理控制方程；采用RNGk-ε湍流模型，用壁面函數(shù)法考慮邊界對(duì)流動(dòng)的影響。在劃分網(wǎng)格時(shí)，盡可能采用正六面體網(wǎng)格單元，固體壁面相鄰的第一個(gè)網(wǎng)格寬度限定為為1mm，最寬一級(jí)的網(wǎng)格寬度≤8cm，在進(jìn)風(fēng)口及壁面附近區(qū)域加密網(wǎng)格，將相鄰兩網(wǎng)格間距的尺度比嚴(yán)格控制地控制在0.8～1.25。2、結(jié)果及分析2．1CFD計(jì)算及可視化實(shí)驗(yàn)結(jié)果工況Ⅰ、工況Ⅱ下的z方向中心截面（z=0.4m，下同）的CFD模擬結(jié)果、流動(dòng)可視化照片如圖2所示?？梢钥闯觯诿鏈囟葘?duì)模型內(nèi)空氣的流動(dòng)形態(tài)有明顯的影響。進(jìn)風(fēng)溫度較高時(shí)，自由射流區(qū)較長(zhǎng)，射流進(jìn)入到1/2的寬度（x=0.6m的平面）后，開(kāi)始向下彎曲，射流能夠沖刷到左側(cè)壁面，并且在0.3m<x<0.6m、0m<y<0.3m范圍內(nèi)形成一個(gè)逆時(shí)針渦；而進(jìn)風(fēng)溫度較低時(shí)，射流進(jìn)入1/4的寬度后，便開(kāi)始向下彎曲，射流沖刷到底面，不能形成明顯的渦。產(chǎn)生這種差別的原因在于當(dāng)進(jìn)風(fēng)溫度較高時(shí)，空氣受到方向向下的重力作用較?。环粗M(jìn)風(fēng)溫度較低時(shí)，空氣受到的重力作用較大。圖2-a工況ⅠⅠ下的流動(dòng)形形態(tài)的計(jì)算結(jié)結(jié)果圖2-b工況Ⅰ下下的流動(dòng)形態(tài)態(tài)的可視化結(jié)結(jié)果圖2-b工況Ⅱ下下流動(dòng)形態(tài)的的計(jì)算結(jié)果圖2-b工況Ⅱ下下的流動(dòng)形態(tài)態(tài)的可視化結(jié)結(jié)果圖3顯示了工況Ⅲ、工況Ⅳ下的CFD模擬結(jié)果以及流動(dòng)可視化照片。對(duì)于垂直送風(fēng)方式，壁面溫度對(duì)模型內(nèi)空氣的流動(dòng)形態(tài)影響較小。只是當(dāng)進(jìn)風(fēng)溫度的較低時(shí)，射流末端向右壁面（x=1.2m）偏轉(zhuǎn)量增大。由于射流左側(cè)的空氣卷吸量大于右側(cè)，左側(cè)空氣速度小于右側(cè)，左側(cè)壓力大于右側(cè)，使得射流向右側(cè)偏轉(zhuǎn)，當(dāng)溫度較低時(shí)，射流所受的熱驅(qū)動(dòng)力和重力較大，射流速度較高，左、右側(cè)壓力差較大，射流偏轉(zhuǎn)量增大。圖3-a工況Ⅲ下下流動(dòng)形態(tài)的的計(jì)算結(jié)果圖3-b工況Ⅲ下流動(dòng)形態(tài)的可視視化結(jié)果圖3-c工況Ⅳ下流動(dòng)形態(tài)的計(jì)算算結(jié)果圖3-b工況Ⅳ下流動(dòng)形態(tài)的可視視化結(jié)果2．2溫度、速度結(jié)果圖4給出了A型進(jìn)風(fēng)方式下各測(cè)點(diǎn)的溫度、速度的計(jì)算值和測(cè)量值。其中速度的測(cè)量值略高于計(jì)算值，由于數(shù)值計(jì)算時(shí)進(jìn)風(fēng)速度在進(jìn)風(fēng)口截面上設(shè)定為處處相同，等于射流平均速度；而在一般情況下，實(shí)驗(yàn)時(shí)的進(jìn)口風(fēng)速在進(jìn)風(fēng)口截面上呈現(xiàn)中心高，四周低的現(xiàn)象，造成射流中心速度高于平均速度。從圖中還可以看出，進(jìn)風(fēng)溫度較高時(shí)，沿測(cè)量線溫度的變化比較平緩，溫度分別升高3.2℃；溫度較低時(shí)，則變化比較劇烈，溫度升高7.9℃,同一進(jìn)風(fēng)速度下，進(jìn)口溫度相差7.5℃時(shí)，出口溫度相差1.8℃，說(shuō)明進(jìn)口溫度較低時(shí)，自然對(duì)流對(duì)換熱起了很重要的作用。在射流的前半?yún)^(qū)（x>600mm），溫度對(duì)各點(diǎn)速度的影響不大，對(duì)于射流的后半?yún)^(qū)（x<600mm）,進(jìn)風(fēng)溫度對(duì)各點(diǎn)速度有一定的影響。進(jìn)風(fēng)溫度較高時(shí)，速度較高，進(jìn)風(fēng)溫度降低后，速度有所下降，從圖2的可視化圖譜可以看出，進(jìn)口溫度降低后，5、6、7三個(gè)測(cè)點(diǎn)已經(jīng)偏離了主射流區(qū)，使得所測(cè)速度變小。圖4-a工況Ⅰ下下各測(cè)點(diǎn)溫度度、速度分布布圖4-b工況Ⅱ下各各測(cè)點(diǎn)溫度、速速度分布圖5給出了B型進(jìn)風(fēng)方式下各測(cè)點(diǎn)的溫度、速度的計(jì)算值和測(cè)量值。在進(jìn)口處速度的測(cè)量值高于計(jì)算值，由于在垂直射流的條件下，進(jìn)口風(fēng)速在進(jìn)風(fēng)口截面上不均勻性增加，射流中心速度遠(yuǎn)高于平均速度，從而與數(shù)值計(jì)算時(shí)設(shè)定的均勻速度截面相矛盾。數(shù)值計(jì)算時(shí)，為了消除風(fēng)口速度不均勻性帶來(lái)的計(jì)算誤差，需要將一定長(zhǎng)度的風(fēng)管納入計(jì)算區(qū)域。從實(shí)驗(yàn)結(jié)果可以看出，進(jìn)風(fēng)溫度較高時(shí)，沿測(cè)量線溫度的變化比較平緩，溫度升高0.5℃；溫度較低時(shí)，則變化比較劇烈，溫度升高4.2℃,同一進(jìn)風(fēng)速度下，進(jìn)口溫度相差7.5℃時(shí)，出口溫度相差1.7℃；B型送風(fēng)方式的射流速度大于A型送風(fēng)方式的射流速度，并且當(dāng)y<300mm時(shí)，速度發(fā)生突降。以上現(xiàn)象說(shuō)明進(jìn)口溫度較低時(shí)，自然對(duì)流對(duì)換熱起了很重要的作用。由于射流的慣性力與受到的熱驅(qū)動(dòng)力、重力方向一致，空氣從進(jìn)風(fēng)口出來(lái)后，會(huì)在熱驅(qū)動(dòng)力、重力的作用下沿射流方向產(chǎn)生加速；同時(shí)，由于壁面加熱產(chǎn)生的浮升力作用，使得射流左側(cè)的順時(shí)針渦速度提高，從而使得射流速度較高，同時(shí)射流在y<300mm的區(qū)域內(nèi)已經(jīng)發(fā)生偏轉(zhuǎn)，測(cè)點(diǎn)位置已經(jīng)偏離了主射流區(qū)，使得測(cè)點(diǎn)速度明顯降低。圖4-a工況Ⅲ下下各測(cè)點(diǎn)溫度度、速度分布布圖4-b工況Ⅳ下各各測(cè)點(diǎn)溫度、速速度分布3、結(jié)論當(dāng)進(jìn)風(fēng)溫度變化時(shí)，壁面溫度對(duì)水平送風(fēng)方式下的空氣流動(dòng)形態(tài)影響較大，對(duì)垂直送風(fēng)方式影響較小；數(shù)值模擬時(shí)，為了減小實(shí)際進(jìn)風(fēng)口處的進(jìn)風(fēng)速度的不均勻性帶來(lái)的誤差，需要將風(fēng)管納入計(jì)算區(qū)域；對(duì)于水平進(jìn)風(fēng)方式，進(jìn)風(fēng)溫度較高引起自由射流區(qū)較長(zhǎng)，并在0.3m<x<0.6m、0m<y<0.3m范圍內(nèi)形成一個(gè)逆時(shí)針渦；進(jìn)風(fēng)溫度較低引起自由射流區(qū)較短，射流沖刷到底面，不能形成明顯的渦，對(duì)于垂直進(jìn)風(fēng)方式，射流末端向右側(cè)彎曲；進(jìn)風(fēng)溫度的降低引起浮升力作用下的自然對(duì)流作用增強(qiáng)，模型整體換熱量提高。參考文獻(xiàn)(References）[1] TuveGL.Airvelocitiesinventilatingjets[J].Heating,PipingandAirConditioning，1953，25(1)：181-191[2]HaghighatF,LiY,MegriAC.Developmentandvalidationofazonalmodel–POMA[J].Buildingandenvironment，2001，36(9)：1039-1047[3] ChenQ.,vanderkoviJ.,MeyersA.Measurementsandcomputationsofventilationefficiencyandtemperatureefficiencyinaventilationroom.EnergyandBuildings,1988,12(2):85-99[4] ChenQ.ANumericalStudyofIndoorAirQualityandThermalComfortunderSixKindsofAirDiffusion[J].ASHRAETransactions，1992，98(1)：203-217[5]FitznerKF.Airflowexperimentsinfull-scaletestrooms[J].ASHRAETransactions，1981，87(2)：1143-1153[6]ZhangJS