趙英博 章文杰 貢桐丹 周建偉

（南京理工大學(xué)能源與動力工程學(xué)院 南京 210094）

0 引言

隨著我國人民生活水平的不斷提高，人們對室內(nèi)環(huán)境的要求也越來越高，各式各樣與室內(nèi)環(huán)境控制有關(guān)的空調(diào)系統(tǒng)應(yīng)運(yùn)而生，它們以新穎的理念和先進(jìn)的技術(shù)水平活躍在各種場合的建筑應(yīng)用中，輻射空調(diào)在其中具有舉足輕重的地位。

與傳統(tǒng)空調(diào)相比，輻射供冷空調(diào)因其獨(dú)特的供冷方式，以熱舒適性好、換熱效率高、運(yùn)行安靜和節(jié)能性高等諸多優(yōu)點(diǎn)[1]受到了人們的青睞，然而在系統(tǒng)運(yùn)行中也存在著冷表面易結(jié)露和單位面積換熱量不足等問題[2,3]。針對輻射空調(diào)冷表面易結(jié)露的問題，國內(nèi)外學(xué)者從供水控制、獨(dú)立送風(fēng)和輻射板優(yōu)化等方面做了大量的研究：袁玉潔[4]提出了輻射供冷加新風(fēng)復(fù)合系統(tǒng)中提高供水溫度、匹配新風(fēng)除濕系統(tǒng)等基礎(chǔ)設(shè)計(jì)措施是防止輻射供冷空調(diào)系統(tǒng)結(jié)露的有效舉措；Kang[2]對比并分析了兩種復(fù)合輻射冷卻空調(diào)系統(tǒng)對室內(nèi)熱舒適性的影響，并列出了置換通風(fēng)通過承擔(dān)室內(nèi)部分負(fù)荷以提高輻射板水溫等多項(xiàng)防結(jié)露的優(yōu)勢，Yin[5]等人設(shè)計(jì)了三塊0.16m 的輻射冷板，通過對比實(shí)驗(yàn)研究發(fā)現(xiàn)與金屬輻射冷板和純管板相比，石膏輻射冷板的傳熱性能最好，并且相同條件下其防結(jié)露性能最佳。

由此可見，國內(nèi)外對輻射冷板的防結(jié)露問題進(jìn)行了大量的研究，在增加能耗或改變能源輸入形式的情況下，很多方法的確能起到很好的防結(jié)露效果，本文所考慮的是在現(xiàn)有新風(fēng)承擔(dān)室內(nèi)全部濕負(fù)荷的情況下，不通過高能耗諸如置換通風(fēng)的手段，僅通過改變新風(fēng)的送風(fēng)方式來達(dá)到防結(jié)露目的。首先探究南京市夏季新風(fēng)貼附射流與輻射供冷復(fù)合空調(diào)系統(tǒng)理想的防結(jié)露工況點(diǎn)，隨后建立貼附射流送風(fēng)和普通散流器送風(fēng)的建筑模型，并對兩種送風(fēng)方式復(fù)合的輻射冷板空調(diào)系統(tǒng)防結(jié)露效果進(jìn)行模擬驗(yàn)證和對比分析。

1 模型及計(jì)算條件

1.1 物理模型

本文以南京市某辦公樓的標(biāo)準(zhǔn)辦公室為建模參考依據(jù)進(jìn)行模擬分析。辦公室的房間規(guī)格（長×寬×高）為5m×4m×2.5m。辦公室僅南墻為外墻，且南墻中間裝有規(guī)格（長×高）為2m×1.5m 的單層玻璃鋁合金窗，該窗戶設(shè)有密織布深色內(nèi)遮陽。辦公室東、西、北面均為內(nèi)隔墻，假設(shè)相鄰房間的室內(nèi)情況與該辦公室相同。

該辦公室采用新風(fēng)貼附射流和輻射供冷復(fù)合的空調(diào)系統(tǒng)，考慮到輻射冷板的結(jié)露問題，濕負(fù)荷全部由新風(fēng)承擔(dān)，氣流組織方式為上送下回，條縫形送風(fēng)口的尺寸（長×寬）為2000mm×50mm，排風(fēng)口的尺寸（長×寬）為250mm×200mm，與之形成對照的普通形式的新風(fēng)散流器送風(fēng)與輻射供冷復(fù)合空調(diào)系統(tǒng)送風(fēng)口與排風(fēng)口尺寸（長×寬）均為250mm×200mm。

辦公室內(nèi)設(shè)備負(fù)荷情況如下：照明設(shè)備為兩盞熒光燈，有效發(fā)熱功率為36W/盞；雙人辦公室配備兩臺電腦，有效發(fā)熱功率為120W/臺。辦公室內(nèi)人員設(shè)備簡化為兩個尺寸（長×寬×高）為0.4m×0.2m×1.2m 的六面體，折合發(fā)熱功率為281W。辦公室的人員設(shè)備及送排風(fēng)口位置如圖1所示。

圖1 辦公室模型示意圖Fig.1 Schematic diagram of the office model

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

1.2.1 模型假設(shè)

為使計(jì)算在保證精確性的情況下避免模擬過程的復(fù)雜化，對模型進(jìn)行如下理想化假設(shè)：

（1）由于辦公室內(nèi)空氣處于常溫、常壓、低速的狀態(tài)，因此假設(shè)辦公室內(nèi)的空氣近似為理想氣體，不考慮空氣調(diào)節(jié)過程中細(xì)微溫度變化帶來空氣相關(guān)參數(shù)的變化，計(jì)算中相關(guān)參數(shù)設(shè)置為常數(shù)；

（2）假設(shè)輻射冷表面溫度均勻恒定，冷板的邊界條件設(shè)置為恒溫邊界；

（3）除南墻以外的其他內(nèi)墻，假設(shè)相鄰辦公室及走廊的室內(nèi)環(huán)境與研究房間完全一致，沒有來自相鄰房間的熱量傳遞，因此內(nèi)墻設(shè)為絕熱邊界；

（4）簡化人體設(shè)備模型為兩個尺寸（長×寬×高）為0.4m×0.2m×1.2m 的穩(wěn)定體熱源。

1.2.2 控制方程

求解過程中采用的控制方程包括動量方程、能量方程和連續(xù)性方程，動量守恒方程如式（1）-（3）所示：

其中，τ為時間，s；u、v、ω分別表示x、y、z軸方向的速度，m/s；μ為動力粘性系數(shù)，取1.81×10-5N·s/m ；g為重力加速度，取9.8m/s 。

能量守恒方程如式（4）所示：

其中，t為溫度，℃；c為比熱容，kJ/kg·℃；q為單位體積的熱源發(fā)熱功率，J/s·m3。

連續(xù)性方程如式（5）所示：

考慮內(nèi)熱源的影響，輻射模型采用可以計(jì)算散射介質(zhì)并能考慮邊界發(fā)射率影響的Do 模型。

1.3 計(jì)算條件

1.3.1 計(jì)算參數(shù)

南京市夏季室外和室內(nèi)設(shè)計(jì)參數(shù)如表1所示，其中，室內(nèi)相對濕度的選擇主要考慮輻射供冷空調(diào)的特點(diǎn)，室內(nèi)的相對濕度越小，則輻射空調(diào)的冷板結(jié)露可能性越小。

表1 相關(guān)設(shè)計(jì)及計(jì)算參數(shù)Table 1 Related design and calculation parameters

新風(fēng)量選取30m3/h·人，則雙人辦公室的新風(fēng)總量為60m3/h，人體的散熱和散濕量選取按照普通事務(wù)工作下成年男子的活動情況，即散熱量60W/人，散濕量110g/h·人。

外墻為240mm 厚的磚墻，內(nèi)外兩側(cè)有水泥砂漿抹面，其傳熱系數(shù)為1.94W/（m ·K）；窗體為單層玻璃鋁合金窗，傳熱系數(shù)為6.64W/(m ·K)，窗戶設(shè)有深黃密織布材料的內(nèi)遮陽。

1.3.2 邊界條件

首先根據(jù)以上計(jì)算參數(shù)，確定房間的負(fù)荷情況，采用冷負(fù)荷系數(shù)法，應(yīng)用HDY-SMAD 空調(diào)負(fù)荷計(jì)算及分析軟件V4.0 進(jìn)行計(jì)算得到結(jié)果如下：夏季最大總冷負(fù)荷發(fā)生時刻為14:00，總冷負(fù)荷Q=1483W，其中新風(fēng)冷負(fù)荷Qx=615W，輻射冷板承擔(dān)負(fù)荷Qf=868W，濕負(fù)荷W=210g/h。

輻射冷板的表面溫度可由文獻(xiàn)《綠色建筑視野下輻射冷量計(jì)算方法》[6]的ASHRAE handbook 理論計(jì)算公式獲得：

其中，QF為單位面積的輻射傳熱量，W/m ；TEP為輻射體表面的平均熱力學(xué)溫度，K；TPJ為室內(nèi)非加熱表面的面積加權(quán)平均熱力學(xué)溫度，K。

采用室內(nèi)空調(diào)設(shè)計(jì)溫度替代其他表面的加權(quán)平均溫度，根據(jù)上述公式可獲得輻射冷板表面溫度TEP=292.7K，即19.7℃。

在新風(fēng)貼附射流和輻射供冷復(fù)合的空調(diào)系統(tǒng)中，由于室內(nèi)全部濕負(fù)荷由新風(fēng)承擔(dān)，據(jù)此獲得獨(dú)立新風(fēng)的送風(fēng)狀態(tài)點(diǎn)。具體分析計(jì)算如下：

其中，dl為室內(nèi)新風(fēng)處理至機(jī)器露點(diǎn)的含濕量，g/kga；G為室內(nèi)新風(fēng)的質(zhì)量，kg。

根據(jù)獲得的新風(fēng)處理至機(jī)器露點(diǎn)的含濕量dl可以獲得該狀態(tài)點(diǎn)的溫度tl，取風(fēng)管設(shè)備溫升為1℃，則可確定新風(fēng)的送風(fēng)溫度ti。

求得含濕量dl=6.51g/kga，此時的狀態(tài)點(diǎn)溫度tl=9.23℃，考慮風(fēng)管設(shè)備溫升，取1℃，則新風(fēng)送風(fēng)溫度to=10.23℃。

根據(jù)送回風(fēng)風(fēng)量和送回風(fēng)口尺寸，可確定模擬的送回風(fēng)口空氣流速，將貼附射流送風(fēng)口設(shè)置為速度入口邊界條件（velocity inlet），回風(fēng)口為壓力出口邊界條件（pressure outlet），出入口邊界條件參數(shù)如下表所示。

表2 送回風(fēng)設(shè)計(jì)參數(shù)Table 2 Design parameters of supply air and return air

透過玻璃窗進(jìn)入的太陽輻射得熱形成的逐時冷負(fù)荷按照下列公式計(jì)算：

其中，CLc為透過玻璃窗進(jìn)入的太陽輻射得熱形成的逐時冷負(fù)荷，W；C clC為透過無遮陽標(biāo)準(zhǔn)玻璃太陽輻射冷負(fù)荷系數(shù)；Cz為外窗綜合遮擋系數(shù)；Cw為外遮陽修正系數(shù)；Cn為內(nèi)遮陽修正系數(shù)；C s為玻璃修正系數(shù)；Djmax為夏季日射得熱因數(shù)最大值；FC為窗玻璃凈面積，m 。

使用Ansys Fluent 軟件完成模擬內(nèi)容，計(jì)算過程中進(jìn)行的設(shè)置情況如下：求解器采用壓力基求解器，壓力速度耦合采用半隱式鏈接壓力方程的SIMPLE 算法，動量方程和能量方程采用計(jì)算穩(wěn)定性較好的一階迎風(fēng)格式，控制方程的離散格式采用標(biāo)準(zhǔn)離散模型，湍流模型根據(jù)氣流特征選擇k-epsilon 雙方程中的標(biāo)準(zhǔn)k-ε模型，輻射模型采用可以計(jì)算散射介質(zhì)并能考慮邊界發(fā)射率影響的的Do 模型。當(dāng)能量殘差值低于10-6，其余監(jiān)測變量的殘差值低于10-3時，認(rèn)為計(jì)算收斂。

2 模擬結(jié)果及分析

2.1 室內(nèi)熱舒適性分析

室內(nèi)空調(diào)設(shè)計(jì)的根本目的是保證室內(nèi)環(huán)境的舒適性，因此在討論輻射供冷空調(diào)的防結(jié)露效果和能耗狀況前應(yīng)先保證該環(huán)境的舒適程度符合標(biāo)準(zhǔn)規(guī)定，根據(jù)上述邊界條件，本研究對室內(nèi)0.1m（腳踝處）、1.1m（人處于坐姿狀態(tài)時溫度感知位置）和1.7m（人處于站姿狀態(tài)時溫度感知位置）的室內(nèi)溫度及風(fēng)速分布進(jìn)行分析，以確保所選擇計(jì)算邊界條件的合理性，同時根據(jù)輻射冷板的表面溫度分布情況對防結(jié)露效果進(jìn)行具體分析，室內(nèi)不同高度的空氣溫度及速度分布情況如圖2所示。

圖2 不同高度空氣溫度模擬結(jié)果Fig.2 Simulation results of air temperature at different aititudes

根據(jù)圖2（a）-（c）可以看出，高度為0.1m處（人體腳踝位置）的溫度分布受南墻的輻射和室內(nèi)人員設(shè)備熱源的影響較大，加上距離輻射冷板和貼附射流最遠(yuǎn)，受南墻和室內(nèi)熱源的影響較大，溫度分布并不均勻，但絕大部分區(qū)域的溫度區(qū)間在23℃-27℃之間，所處溫度范圍仍在可接受的舒適范圍內(nèi)；1.1m 處和1.7m 處的溫度分布較為均勻，且溫度大多位于23℃-25℃之間，所處溫度范圍在人體熱舒適范圍內(nèi)，此外1.1m 處和1.7m 處靠近南墻的位置出現(xiàn)了溫度陡降的區(qū)域，結(jié)合后面的室內(nèi)空氣流速分布圖推測，是新風(fēng)射流在受到熱源加熱的抬升效果后將自身的低溫空氣和沿程卷吸的輻射冷板降溫的空氣輸送到南墻位置，被南墻阻擋后跌落形成的局部冷域；此外人員設(shè)備模型周圍的空氣溫度較高，是因?yàn)閷⒃撃Ｐ驮O(shè)置為體熱源的緣故。

根據(jù)《民用建筑供暖通風(fēng)空氣調(diào)節(jié)設(shè)計(jì)規(guī)范》[7]的相關(guān)規(guī)定，室內(nèi)風(fēng)速應(yīng)不大于0.3m/s，由圖3（a）-（c）速度分布云圖可以看出，由于輻射供冷的均勻性特點(diǎn)和新風(fēng)射流的風(fēng)量小、出風(fēng)速度低的原因，室內(nèi)各高度的速度分布絕大部分在0.04m/s 以下；根據(jù)室內(nèi)環(huán)境國際標(biāo)準(zhǔn)ISO7730 的標(biāo)準(zhǔn)要求，工作區(qū)內(nèi)（即0.1m-1.1m 高度區(qū)間內(nèi)坐姿辦公的空間區(qū)域）的垂直溫度差不應(yīng)大于3℃，由圖4可以看出，除前文分析的冷域部分，模擬辦公環(huán)境的2m 以下的絕大部分區(qū)域溫度分布較為均勻，溫度在296K-298K 之間，除去所假設(shè)熱源的集中影響，認(rèn)為辦公環(huán)境的熱舒適性符合要求，因此所選擇的送回風(fēng)參數(shù)及輻射冷板的表面溫度等設(shè)計(jì)參數(shù)是合理的。

圖3 不同高度空氣流速模擬結(jié)果Fig.3 Simulation results of air velocity at different altitudes

圖4 z=1.1m 處溫度分布云圖Fig.4 Cloud map of temperature distribution at z=1.1m

2.2 輻射冷板結(jié)露情況對比分析

分析圖5（a）可知，輻射冷板表面假設(shè)為溫度均勻的恒溫邊界，且模擬結(jié)果在比19.7℃稍高的溫度下進(jìn)行輻射傳熱，該結(jié)果可以由貼附射流在冷板表面加強(qiáng)了對流換熱得到解釋，此時的輻射冷板下方空氣溫度分布如圖5（b）所示，可以看出該高度區(qū)間的溫度分布受貼附射流的影響較大，從風(fēng)口向外送風(fēng)的過程中，射流的溫度逐漸提高，并在接觸南墻后向四周墻壁擴(kuò)散，假設(shè)室內(nèi)空氣組分均勻，除輻射冷板與送風(fēng)口接觸位置的溫度區(qū)域受10.4℃的新風(fēng)影響外，其余絕大部分溫度區(qū)域在19℃以上，高于室內(nèi)狀態(tài)點(diǎn)露點(diǎn)溫度13.3℃，溫差在6℃以上，因此輻射冷板在設(shè)定工況下不會結(jié)露。且從溫度分布的占比情況來看，14℃-17℃以及17℃-19℃的溫度區(qū)域主要受低溫的新風(fēng)貼附射流影響，由于新風(fēng)的露點(diǎn)溫度更低，為7.8℃，且19℃-21℃的溫度區(qū)間占貼壁處總體的面積很小，絕大部分的溫度區(qū)域分布在21℃-25℃之間，可以認(rèn)為貼壁處空氣的結(jié)露可能性極低。

圖5 新風(fēng)貼附射流冷板表面及貼壁處空氣分布Fig.5 Air distribution on the surface of the cold plate attached to the jet and adhering to the wall

分析圖6（a）可知，普通送風(fēng)方式下的輻射冷板表面由于受新風(fēng)的氣流組織影響較小，因此冷板表面溫度穩(wěn)定在19.7℃的設(shè)定工況下運(yùn)行，且此時輻射冷板下方貼壁處的空氣溫度分布如圖6（b）所示，從圖中可以看到，新風(fēng)對貼壁處空氣溫度的影響在房間的后半部分開始出現(xiàn)，且沿新風(fēng)輸送的方向核心溫度呈下降趨勢，最低降到了17℃左右，該溫度同樣高于室內(nèi)狀態(tài)點(diǎn)的露點(diǎn)溫度13.3℃，在設(shè)定工況下也不會結(jié)露，然而，從溫度分布的占比情況來看，新風(fēng)影響的低溫核心區(qū)（17℃-19℃）面積相比較貼附射流房間的低溫區(qū)面積要大很多，且絕大部分的溫度分布于19℃-23℃，和貼附射流房間比較貼壁處的整體平均溫度下降2℃左右。因此，當(dāng)辦公室內(nèi)濕負(fù)荷陡增，例如室內(nèi)人員突然增加時，普通送風(fēng)方式的輻射空調(diào)房間結(jié)露可能性將大大提高。

圖6 普通送風(fēng)冷板表面及貼壁處空氣溫度分布Fig.6 Air temperature distribution on the surface of the common air supply cold plate and the wall

3 結(jié)論

總結(jié)以上對比分析，在南京市氣象條件下，本文模擬的辦公室在新風(fēng)承擔(dān)室內(nèi)全部濕負(fù)荷時，無論是貼輻射流送風(fēng)還是普通送風(fēng)方式對輻射冷板都有防結(jié)露的作用，但是新風(fēng)貼附射流送風(fēng)強(qiáng)化了輻射冷板表面的對流換熱，使貼壁處整體溫度比普通送風(fēng)方式的貼壁處溫度更高，更加遠(yuǎn)離室內(nèi)露點(diǎn)溫度，由此在室內(nèi)濕負(fù)荷突然增大的情況下結(jié)露風(fēng)險(xiǎn)更小。