董 凡,成永剛

(西安交通大學(xué)網(wǎng)絡(luò)信息中心,陜西 西安 710049)

1 引言

冷卻空調(diào)系統(tǒng)設(shè)計不合理,會導(dǎo)致制冷設(shè)備出現(xiàn)局部過熱的現(xiàn)象[1]。為了保證冷卻空調(diào)系統(tǒng)運行過程中的溫度,需要調(diào)低整個空調(diào)機房的溫度,導(dǎo)致大部分的空調(diào)系統(tǒng)在運行過程中并未達到滿負荷運行的狀態(tài),不僅會造成大量的能源浪費,而且無法徹底解決空調(diào)局部過熱的現(xiàn)象[2]。在設(shè)計冷卻空調(diào)系統(tǒng)的氣流組織時,如果缺乏合理性,就會升高冷卻空調(diào)系統(tǒng)的能耗,因此對冷卻空調(diào)系統(tǒng)氣流組織進行數(shù)值模擬研究,有效處理空調(diào)設(shè)備產(chǎn)生的熱量,從而提高冷卻空調(diào)系統(tǒng)在運行中的制冷效率,在減少冷卻空調(diào)能耗的前提下,設(shè)置冷卻空調(diào)系統(tǒng)數(shù)據(jù)中心,保證冷卻空調(diào)系統(tǒng)在運行中的安全性和穩(wěn)定性[3]。

周鑫濤等人[4]提出了一種基于流體力學(xué)的氣流組織模擬方法,在不同氣流組織形式下,采集了深海載人平臺氣流組織數(shù)據(jù),分析出空氣中溫度的分布規(guī)律,結(jié)果表明,在深海載人平臺密閉艙室的熱源氣流部分不均勻的情況下,外側(cè)的送風(fēng)氣流組織不能滿足密閉艙室的氣流需求,通過更改氣流組的輸送方式,滿足氣流組織要求。趙福云等人[5]為了分析工作臺送風(fēng)口對電子潔凈室內(nèi)空氣環(huán)境的影響,結(jié)合流體力學(xué)的計算方法,對iso6模式下的電子潔凈室空氣環(huán)境數(shù)值進行模擬,根據(jù)送風(fēng)口的個數(shù)以及所處位置的不同,分析工作臺送風(fēng)口對電子潔凈室內(nèi)空氣環(huán)境的影響。實驗結(jié)果表明,處于電子潔凈室正前方的工作臺送風(fēng)口可以有效抑制熱氣流的影響,工作臺送風(fēng)口可以根據(jù)電子潔凈室的位置作出相應(yīng)的位置調(diào)整,送風(fēng)口應(yīng)設(shè)置在電子潔凈室正交正對的位置,方可保證送風(fēng)口的風(fēng)速,有助于電子潔凈室的空氣環(huán)境的流動性,通過數(shù)值模擬的結(jié)果,可以為電子潔凈室的送風(fēng)口的位置設(shè)置提供合理化建議。

基于以上研究背景,本文利用CFD仿真軟件對冷卻空調(diào)系統(tǒng)的氣流組織進行數(shù)值模型研究,從而確定影響冷卻空調(diào)系統(tǒng)氣流組織的因素。

2 冷卻空調(diào)系統(tǒng)氣流組織數(shù)值模擬

2.1 控制方程

控制方程包括質(zhì)量守恒、動量守恒和能量守恒[6],在冷卻空調(diào)系統(tǒng)氣流組織的數(shù)值模擬中,可以解決氣流組織的流動和傳熱問題[7]。本文以數(shù)據(jù)中心冷卻空調(diào)系統(tǒng)為研究對象,假設(shè)數(shù)據(jù)中心的密封性良好[8],那么控制方程為:

質(zhì)量守恒

(1)

其中,vx、vy和vz分別表示速度矢量在x、y和z方向的分量。

動量守恒

(2)

其中,ρ表示空氣的密度,p表示空氣中的大氣壓強,?表示運動黏度。

能量守恒

(3)

其中,T表示熱力學(xué)溫度,λ表示導(dǎo)熱系數(shù),Sh表示內(nèi)熱源,Φ表示耗散函數(shù)。

根據(jù)質(zhì)量守恒方程、動量守恒方程和能量守恒方程的表達式,得到了冷卻空調(diào)系統(tǒng)氣流組織數(shù)值模擬的控制方程。

2.2 構(gòu)建冷卻空調(diào)系統(tǒng)模型

將CFD仿真軟件應(yīng)用到冷卻空調(diào)系統(tǒng)模型的構(gòu)建中時,分為多孔介質(zhì)模型和多孔階躍模型[9],相比于多孔介質(zhì)模型而言,多孔階躍模型的結(jié)構(gòu)更加簡單,在計算過程中也更加穩(wěn)定[10]。多孔地板通常比較薄,本文利用多孔階躍模型進行建模[11],那么地板的壓力損失為

Δp=K(0.5ρv2)

(4)

式中,K表示流動阻力系數(shù),表示v地板附近的速度。

根據(jù)式(4)的結(jié)果,計算地板的流動阻力,公式為

(5)

其中,F表示地板的孔隙率。

以上根據(jù)地板的壓力損失,構(gòu)建冷卻空調(diào)系統(tǒng)模型,接下來利用CFD仿真軟件對冷卻空調(diào)系統(tǒng)氣流組織進行數(shù)值模擬。

2.3 CFD仿真

利用CFD仿真軟件對冷卻空調(diào)系統(tǒng)氣流組織進行數(shù)值模擬時,將仿真過程劃分為三個階段:先對冷卻空調(diào)系統(tǒng)模型進行處理[12],轉(zhuǎn)換成CFD仿真軟件可以識別的數(shù)據(jù)模型,再利用求解器讀取前處理的氣流組織數(shù)據(jù)[13],經(jīng)過計算之后,將氣流組織數(shù)據(jù)輸出,結(jié)合后處理軟件處理求解器輸出的氣流組織數(shù)據(jù)[14],并通過繪制圖表展示可視化結(jié)果,本文利用CFD仿真軟件建立數(shù)據(jù)中心冷卻空調(diào)系統(tǒng)模型[15],并將其導(dǎo)入到前處理器中,根據(jù)氣流組織數(shù)據(jù)的分析,輸出結(jié)果。

3 結(jié)果分析

3.1 架空地板高度對冷卻空調(diào)系統(tǒng)氣流組織的影響

當(dāng)冷卻空調(diào)系統(tǒng)采用地板下送風(fēng)形式時,架空地板高度影響冷卻空調(diào)機房的層高,另一方面,架空地板高度影響冷卻空調(diào)系統(tǒng)的送風(fēng)速度以及送風(fēng)的均勻性,從而影響到冷卻空調(diào)的整體制冷效果。因此,本文對不同架空地板高度的氣流組織數(shù)值進行模擬,通過選取200mm、400mm、600mm、800mm、1000mm五種不同的架空地板高度,分析架空地板高度對冷卻空調(diào)氣流組織數(shù)值的影響。

3.1.1 風(fēng)量分析

風(fēng)量大小直接反映冷卻空調(diào)系統(tǒng)的氣流組織到達機房各個區(qū)域的位置。根據(jù)模擬結(jié)果,統(tǒng)計不同架空地板高度情況下的出風(fēng)量情況,表1為不同架空地板高度下風(fēng)量大小分布的區(qū)間。

表1 地板格柵出風(fēng)量的分布情況

根據(jù)表1的風(fēng)量分布區(qū)間,繪制出冷卻空調(diào)系統(tǒng)氣流組織區(qū)間圖,如圖1所示。

圖1 出風(fēng)量分布區(qū)間圖

從圖1中可以看出,當(dāng)架空地板高度在200mm到1000mm之間時,冷卻空調(diào)系統(tǒng)氣流組織的風(fēng)量大小區(qū)間呈縮短趨勢。

3.1.2 壓力場分析

架空地板高度下的壓力場會影響冷卻空調(diào)系統(tǒng)氣流組織分布的均勻性,進而影響到冷卻空調(diào)的制冷效果。綜上所述風(fēng)量對冷卻空調(diào)系統(tǒng)氣流組織數(shù)值的影響,指出地板間隙的風(fēng)量的大小對空調(diào)內(nèi)側(cè)區(qū)域的影響,之所以會出現(xiàn)這種現(xiàn)象因為架空地板高度下的壓力場分布不均,受架空地板高度的影響,其壓力場的分布規(guī)律也不一樣。

分別對地板間隙的靜壓力與架空地板高度下的壓力場進行分析,根據(jù)模擬結(jié)果,不同架空地板高度下的地板間隙靜壓力與壓力場分布區(qū)間的統(tǒng)計情況如表2所示。

表2 地板格柵與平面的靜壓分布情況

根據(jù)表2的風(fēng)量分布區(qū)間,繪制了冷卻空調(diào)系統(tǒng)氣流組織的壓力場分布區(qū)間圖,如圖2所示。

圖2 靜壓分布區(qū)間圖

從圖2可以看出,當(dāng)架空地板高度在200mm到1000mm時,靜壓大小呈縮短趨勢,說明隨著架空地板的高度升高,地板間隙的靜壓力與架空地板高度氣流組織的壓力場的差異越來越小。

3.1.3 溫度場分析

在完成架空地板高度的出風(fēng)量與壓力場的分析后,可以看出架空地板高度對冷卻空調(diào)機房系統(tǒng)的氣流組織具有一定的影響作用。然而,架空地板高度的出風(fēng)量和壓力場對冷卻空調(diào)系統(tǒng)氣流組織數(shù)值的影響都會在制冷溫度上反映出來,只有保證冷卻空調(diào)系統(tǒng)的機房溫度適宜,才能保證冷卻空調(diào)系統(tǒng)的正常運行,因此通過對溫度場分析,可以獲取如何改善冷卻空調(diào)系統(tǒng)氣流組織的方法。為了分析不同的架空地板高度下的溫度場情況,分別截取不同高度下的溫度分布,并將五個不同的架空地板高度下的最高溫度值和冷卻空調(diào)系統(tǒng)機房的最高溫度進行統(tǒng)計,如表3所示。

根據(jù)表3得出不同架空地板高度下,冷卻空調(diào)系統(tǒng)氣流組織的最高溫度變化規(guī)律,如圖3所示。

由表3可以看出,升高架空地板的高度,會降低數(shù)據(jù)中心的溫度。當(dāng)架空地板高度為200mm時,冷卻空調(diào)系統(tǒng)機房會出現(xiàn)氣流組織分布不均勻的情況,此時的最高溫度為44℃。當(dāng)架空地板的高度為600mm時,冷卻空調(diào)系統(tǒng)機房的溫度會有所下降,因此可以總結(jié)出,隨著地板架空高度的逐漸升高,冷卻空調(diào)機房的溫度會隨之下降,此時的冷卻空調(diào)系統(tǒng)的冷卻效果會有所提高。

3.2 地板送風(fēng)格柵開孔率對冷卻空調(diào)系統(tǒng)氣流組織的影響

冷卻空調(diào)系統(tǒng)的冷卻目標(biāo)是保證冷卻空調(diào)系統(tǒng)的入口溫度不超過最高入口的溫度。冷卻空調(diào)系統(tǒng)機房的溫度過高,是因為地板送風(fēng)格柵的冷氣不足造成的,最終導(dǎo)致熱空氣倒流,形成局部溫度升高的現(xiàn)象,最終對冷卻空調(diào)系統(tǒng)的制冷效果產(chǎn)生影響。因此,對地板送風(fēng)格柵的開孔率進行優(yōu)化,是保證冷卻空調(diào)系統(tǒng)冷卻效率的重要環(huán)節(jié)。本文分別將地板送風(fēng)格柵的開孔率設(shè)置為20%,35%,50%及組合布置形式,并分別進行4組試驗以驗證地板送風(fēng)格柵開孔率對冷卻空調(diào)系統(tǒng)氣流組織的影響。

3.2.1 風(fēng)量分析

風(fēng)量大小可以直接反映冷卻空調(diào)系統(tǒng)氣流組織數(shù)值的大小,根據(jù)模擬的結(jié)果來看,通過統(tǒng)計地板送風(fēng)格柵的出風(fēng)量,表4表示不同地板送風(fēng)格柵開孔率影響下的風(fēng)量大小分布情況。

表4 地板格柵出風(fēng)量的分布情況

根據(jù)表4的風(fēng)量分布區(qū)間,繪制了冷卻空調(diào)系統(tǒng)氣流組織的區(qū)間圖,如圖4所示。

圖4 出風(fēng)量分布區(qū)間圖

由圖4可以看出,地板送風(fēng)格柵的開孔率在20%到50%區(qū)間時,風(fēng)量呈逐漸增大的趨勢,說明隨著地板送風(fēng)格柵開孔率減少,冷卻空調(diào)系統(tǒng)氣流組織的出風(fēng)量逐漸減小。

3.2.2 壓力場分析

分析地板送風(fēng)格柵與冷卻空調(diào)系統(tǒng)氣流組織壓力場,根據(jù)模擬結(jié)果,得到不同地板送風(fēng)格柵開孔率與壓力場的關(guān)系,如表5表示。

表5 地板格柵與平面的靜壓分布情況

從表5可以看出壓力場的分布區(qū)間,并繪制區(qū)間圖,如圖5所示。

圖5 靜壓分布區(qū)間圖

由圖5可以看出,當(dāng)?shù)匕逅惋L(fēng)格柵開孔率為20%時,冷卻空調(diào)系統(tǒng)氣流組織的壓力場值較大,隨著地板送風(fēng)格柵開孔率的增大,壓力場的值隨之減小,在組合布置形式下的靜壓力所對應(yīng)的板格柵開孔率為20%-35%,滿足冷卻空調(diào)系統(tǒng)的送風(fēng)量要求。

3.2.3 溫度場分析

為了分析不同地板送風(fēng)格柵開孔率下,冷卻空調(diào)系統(tǒng)氣流組織的溫度場情況,將各個截面的最高溫度值與冷卻空調(diào)機房的最高溫度值進行統(tǒng)計,結(jié)果如表6所示。

表6 冷卻空調(diào)系統(tǒng)氣流組織在不同界面的最高溫度分布情況

根據(jù)表6的結(jié)果,繪制不同地板送風(fēng)格柵開孔率下溫度場的變化情況,如圖6所示。

圖6 溫度變化折線圖

從表6和圖6可以看出,當(dāng)?shù)匕逅惋L(fēng)格柵的開孔率控制在20%-35%范圍內(nèi)時,隨著地板送風(fēng)格柵的開孔率的增加,溫度場的溫度明顯下降,說明冷卻空調(diào)系統(tǒng)的制冷效率明顯提高,冷卻效果良好,在此區(qū)間范圍內(nèi)可以滿足模擬要求。

4 結(jié)束語

本文提出了基于CFD仿真的冷卻空調(diào)系統(tǒng)氣流組織數(shù)值模擬研究,發(fā)現(xiàn)架空地板高度和地板送風(fēng)格柵開孔率會影響冷卻空調(diào)系統(tǒng)氣流組織的冷卻效果。但是本文的研究還存在很多不足,在今后的研究中,希望可以考慮到封閉冷/熱通道對冷卻空調(diào)系統(tǒng)氣流組織的影響,進一步改善冷卻空調(diào)系統(tǒng)氣流組織的冷卻效果。