一、網格的劃分
網格生成對CFD計算至關重要,直接關系到CFD計算的成敗。常見的網格有結構化網格、非結構化網格與自適應笛卡爾網格。質量高的計算網格應滿足以下幾點:
(1)網格在求解域內的疏密變化應與被求變量的梯度變化相適應;
(2)整個求解域內的單元變化是光滑的進行;
(2)單元的扭曲度較小,單元的寬高比不大于5:1。
根據恒溫恒濕房間的幾何尺寸,采用非結構網格填充率不高,產生的網格數量要比結構網格的數量大得多,因而要求較大的計算機內存,計算時間長,所以本次模擬采用六面體結構網格。計算區域劃分步長不大于恒溫恒濕室實際尺寸的1/20。其他網格劃分規則統一規定為:流體窄縫(fluid gap)的最小網格數為4,固體邊的最小網格數為4,圓柱形與三角形的表面最小網格數均為4。所得網格數為253400個左右。
二、風口模型的簡化
空調通風房間一般都是通過送風口的送風射流來實現送風和室內空氣的混合,以達到空調和通風的目的。為了正確的預測室內的空氣流動,應將送風參數的詳細情況正確的描述成為入流邊界條件。而實際的風口幾何形狀十分復雜,種類也較多,如條縫風口,百葉風口、方型散流器等。目前主要采用間接描述的盒子類風口模型和直接描述類風口模型來簡化描述復雜的空調風口入流邊界條件。本文對于以上問題的處理辦法是將送風口的送風速度取為平均速度,以避免因考慮送風參數細節而導致網格過密的難題。因此,本論文的模擬計算采用傳統的CFD方法對風口進行描述,既將風口當作一個簡單的開口。
三、數值模擬方法
由于廠房和房內障礙物都簡化成六面體,所以模型建立在笛卡爾坐標系下,采用控制容積離散方法將微分方程離散成代數方程。將微分形式離散成代數形式的過程中采用的是上限差分格式,網格的劃分采用多重網格劃分法。
對整個控制方程的求解采用的是SIMPLE算法,在每一次迭代結束后進行下一輪迭代之前,所賦的新值采用欠松弛方法,以此加強非線性耦合方程迭代求解的穩定性。不同的變量所采用的松弛因子不同,對于速度變量松弛因子一般在0.3~0.4之間,本次模擬取0.4;壓力項的松弛因子一般0.6~0.7之間,本次模擬取0.7;紊流動能和紊流耗散能的松弛因子一般在0.25~0.35之間,本次模擬取0.35;其它變量的松弛因子取1.0。在求解的過程中,判斷各個變量收斂的標準也不相同,對于速度變量采用的是10-3,質量守恒項采用的是5×10-3,紊流動能和紊流耗散能采用的是5×10-3,濃度方程變量采用的是10-4,能量方程變量采用的是10-6,其殘差曲線如圖3-3所示。
在前述條件已定的情況下,決定變量是否收斂或者收斂快慢的兩個主要因素是網格的劃分和初值的給定。在最初的計算時,賦初值為零,網格為粗網格,將計算出來的結果作為下一次計算的初值。然后再細化網格,進行新一輪的計算。數值模擬后的殘差收斂曲線如圖3-3。
(2)它的外部大空間的溫度為23℃,因此邊界條件可以簡化處理,房間的地板、頂棚和四面墻都處理成絕熱面;
(3)廠房內空氣為輻射透明介質。
http://m.bag2do.cn