翟玉玲,鐘曉暉,勾昱君,曹樹星,周樹光

(河北聯(lián)合大學(xué)冶金與能源學(xué)院,河北 唐山063009)

0 引 言

冷庫作為現(xiàn)代物流的基礎(chǔ)設(shè)施之一,呈快速發(fā)展的趨勢。發(fā)達國家如日本、美國、芬蘭、加拿大等國的制冷行業(yè)尤為完善。特別是作為世界上最大速凍食品生產(chǎn)國的日本,其低溫庫(-20℃)的數(shù)量占冷庫總數(shù)的80%以上[1]。自從20世紀(jì)60年代中期我國開始建造第一座冷庫起,至今冷庫網(wǎng)絡(luò)遍布全國各地,主要分布在以肉類、水果、蔬菜為主產(chǎn)區(qū)的城市。

隨著生活水平的提高,人們對冷凍食品的需求量和質(zhì)量不斷提高。但是冷庫的冷凍冷藏能力越大,對能源的消耗也越大。在如今以倡導(dǎo)環(huán)保節(jié)能型為口號的社會中,冷庫的節(jié)能問題越來越受到關(guān)注。但是冷庫的設(shè)計也存在很多不合理的地方,不合理的設(shè)計導(dǎo)致庫內(nèi)氣流組織的不平衡,造成能源的浪費[2]～[8]。比如冷庫門長時間的開啟,增加冷庫的能耗。一般具有冰凍保鮮能力的冷庫溫度都在0℃以下,與室外環(huán)境溫度差別較大。在開關(guān)門過程中,由于室內(nèi)外冷熱空氣的熱濕、熱質(zhì)交換,將導(dǎo)致冷庫大門處易于結(jié)霜,難以開啟。數(shù)值模擬因其具有理論性和實踐性的雙重特點,為現(xiàn)代科學(xué)中許多復(fù)雜流動和傳熱問題提高了一種行之有效的方法,被廣泛應(yīng)用于實踐工程中。本論文應(yīng)用CFD模擬技術(shù)分別模擬冷庫在開關(guān)門各10s、20s、30s三種情況下庫門附近的速度場和溫度場。

1 物理模型及數(shù)學(xué)模型的建立

1.1 物理模型

選取唐山市某小型冷庫為研究對象,該庫的內(nèi)部尺寸為長×寬×高=4.5m×3.3m×3.0m。采用單機吊頂式翅片管式蒸發(fā)器制冷,由于蒸發(fā)器的形狀對冷庫內(nèi)的流場和溫度場影響不大,因而將蒸發(fā)器模型簡化為一長方體,該蒸發(fā)器的尺寸為長×寬×高=0.7m×0.3m× 0.5m,置于冷庫天花板側(cè)面;回風(fēng)口位于冷庫下方,尺寸為長×寬×高=0.5m×0.3m×0.5m,送風(fēng)口和回風(fēng)口位置見圖1;冷庫外墻采用保溫隔熱層;冷庫門位于中間,尺寸為長×高=1.2m×2.0m。將冷庫視為三維模型,該模型立體結(jié)構(gòu)圖如圖1所示。

圖1 冷庫的立體結(jié)構(gòu)示意圖

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

建立數(shù)學(xué)模型的目的在于建立反應(yīng)求解問題的基本控制方程和相應(yīng)的定解條件。具體地說就是建立反應(yīng)問題的各個量之間的微分方程,數(shù)值模擬的基礎(chǔ)就是建立正確合理的數(shù)學(xué)模型,脫離了這個基本點,數(shù)值模擬無從而談。本論文模擬的開關(guān)門工況有外界動力強制庫內(nèi)氣體流動,空氣流動呈混亂無序狀態(tài),流場的雷諾數(shù)數(shù)量級大,可達到106,屬于大空間強制對流換熱問題?？刹捎霉こ躺蠌V泛使用的標(biāo)準(zhǔn)k-ε模型求解流動及換熱問題,此模型共包含的5個控制方程,分別是連續(xù)性方程、動量守恒方程、能量守恒方程、湍動能k方程及湍動耗散率ε方程。對于近壁處由于空氣流動緩慢,呈貼壁狀態(tài),其流動有可能出現(xiàn)層流現(xiàn)象,所以對這一區(qū)域的氣流采用壁面函數(shù)法處理。上述5個控制方程表述如下[9]:

連續(xù)性方程:

能量守恒方程:

動量守恒方程:

標(biāo)準(zhǔn)k-ε模型:

其中,

模型各常數(shù)值:C1ε=1.44,C2ε=1.92,Cμ=0.09,σk=1.0,σε=1.3

2 邊界條件的設(shè)定

一個完整數(shù)學(xué)描述應(yīng)該包括控制方程組及相應(yīng)的定解條件,只有合理的邊界條件才有可能計算出流場的解。對于非穩(wěn)態(tài)問題,必須給定初始條件,它是所研究對象開始時刻各物理量在空間上的分布情況,而對于穩(wěn)態(tài)問題,則無須給定。邊界條件是在求解區(qū)域上對所求變量或某一導(dǎo)數(shù)隨空間和時間的變化規(guī)律。

冷庫門附近的流體由于風(fēng)壓和氣壓的共同作用,將產(chǎn)生強烈的熱質(zhì)交換,特別是當(dāng)有貨物進出冷庫時,冷庫門需要長時間的開啟,將破壞冷庫內(nèi)部的氣體流場,嚴(yán)重影響冷藏品的質(zhì)量。本文分別模擬開關(guān)門10s、20s、30s三種情況下庫門附近的速度場和溫度場。為簡化模型,將冷庫視為空庫處理。

該工況下設(shè)置:

(1)入口邊界:蒸發(fā)器入口采用速度邊界,V=2.5m/s,t=-15℃;開門時,大門處的流體處于自然對流狀態(tài),可采用速度入口,V=0.2m/s,t=30℃;關(guān)門時,庫門設(shè)為壁面條件。

(2)蒸發(fā)器出口邊界:采用outflow邊界,出口壓力未知,風(fēng)速由軟件計算得出。

(3)壁面:采用第一類邊界條件。壁面溫度取庫內(nèi)初始溫度,其值為30℃。粘性流動中,壁面可默認(rèn)為無滑移邊界條件。外壁與周圍環(huán)境無熱量交換,視為絕熱。

3 網(wǎng)格劃分

網(wǎng)格是劃分幾何模型的基礎(chǔ),其質(zhì)量直接影響到計算效率和結(jié)果的精確度。對于復(fù)雜的幾何模型,網(wǎng)格的生成既耗時又易出錯。據(jù)不完全統(tǒng)計,對于大型且復(fù)雜的幾何模型,網(wǎng)格的選擇和生成時間占實際CFD計算時間的一半以上。因此,網(wǎng)格劃分是CFD模擬的關(guān)鍵部分,必須給予足夠的關(guān)注。網(wǎng)格又分為結(jié)構(gòu)網(wǎng)格和非結(jié)構(gòu)網(wǎng)格,在GAMBIT軟件里建立模型,無論是結(jié)構(gòu)化網(wǎng)格還是非結(jié)構(gòu)化網(wǎng)格,其生成過程都需遵循圖2所示的步驟。本文的冷庫模型相對簡單,所以采用結(jié)構(gòu)網(wǎng)格劃分,網(wǎng)格單元選取3D的六面體網(wǎng)格。對局部區(qū)域如進風(fēng)口、出風(fēng)口及大門處加密。冷庫的體網(wǎng)格圖如圖3所示。

圖2 GAMBIT建模步驟

圖3 冷庫體網(wǎng)格圖

4 求解過程及結(jié)果分析

4.1 求解過程

無論是氣體流動問題、傳熱問題,還是穩(wěn)態(tài)問題、非穩(wěn)態(tài)問題,CFD求解過程都可按照圖4所示的總體計算流程計算。對所選的控制方程組按照控制體積法(Control Volume Method,簡稱CVM)進行離散,流場的求解方法可選用基于壓力(Pressure-Based Method)的解法,開關(guān)門工況屬于非穩(wěn)態(tài)穩(wěn)問題,可采用壓力的隱式算子分割算法(Pressure Implicitwith Sp litting of Operators,簡稱PISO),這是一種基于無迭代的瞬態(tài)計算程序,采用預(yù)測-修正-再修正三步,它的精度依賴于所選的時間步長,所以在預(yù)測修正過程中,使用的時間步長越小,可取得的計算精度越高,但所用的時間也越長。考慮到氣體流動傳熱,啟動能量方程,激活傳熱機制。對于本文設(shè)定工況下的簡單氣體流動,在開始計算時松弛系數(shù)可選較大值,這可大大改善求解器的性能,加快收斂速度。

4.2 結(jié)果分析

由于三維計算結(jié)果在FLUENT里無法清楚顯示冷庫內(nèi)的速度場及溫度場,因此本文選取幾個具有代表性的切面,以便分析結(jié)果。選取代表性切面x=2.25 m(包含冷庫門的界面)。圖5～7分別是開門10 s、20 s、30 s的速度場,圖中的淺色區(qū)域速度為負(fù)值,說明室外空氣隨著開門時間的推進,慢慢滲透進冷庫內(nèi)部,最后致使庫內(nèi)的冷熱氣體嚴(yán)重混亂。圖8～10分別是開門10 s、20 s、30 s的溫度場,對應(yīng)的溫度場變化很大,說明在x方向上氣體進行了強烈的熱質(zhì)交換。當(dāng)室外風(fēng)進入庫內(nèi)時,由于冷空氣的密度比熱空氣大,室外風(fēng)上揚,首先影響冷庫上部的流場,隨著時間的推進的逐步影響整個流體空間,但是由于蒸發(fā)器的作用,因此室內(nèi)深部的溫度場變化不大。

圖4 CFD工作流程圖

圖11～13是關(guān)門10 s、20 s、30 s的速度場,圖內(nèi)所示速度均為正值,說明庫內(nèi)關(guān)閉后已無室外空氣進入,蒸發(fā)器入口開始對熱空氣進行降溫,此階段庫內(nèi)的氣流仍混亂,說明恢復(fù)到穩(wěn)態(tài)狀態(tài)還需一段時間。

圖14～16為關(guān)門10 s、20 s、30 s的溫度場,從圖中可以看到,從關(guān)門瞬間起庫內(nèi)的溫度下降速率很快,這是因為蒸發(fā)器風(fēng)速是室外風(fēng)速的十幾倍。由于蒸發(fā)器位于冷庫上部,關(guān)門后冷庫上部的溫度場恢復(fù)較快,隨著時間的推進冷庫內(nèi)部溫度場逐步恢復(fù)到穩(wěn)態(tài)狀態(tài)。

5 小 結(jié)

本文利用FLUENT模擬了非穩(wěn)態(tài)狀態(tài)下開啟冷庫門的工況,分析模擬結(jié)果發(fā)現(xiàn)冷庫門開啟時對庫內(nèi)的速度場和溫度場擾動劇烈,尤其是對庫內(nèi)附近的影響更為嚴(yán)重。因為存在嚴(yán)重的熱濕交換,為了保證庫門處的溫度場的穩(wěn)定及避免大門處結(jié)霜、結(jié)冰現(xiàn)象,可在門處加上空氣幕,部分阻隔室外風(fēng)的進入。

[1] 中華工程資格考試網(wǎng)[EB/OL].http://www.100gczg.com,2010-1-8.

[2] 謝晶,吳天.小型冷庫開關(guān)門過程溫度場的數(shù)值模擬[J].上海水產(chǎn)大學(xué)學(xué)報,2006,15(3):333～339.

[3] 楊磊,汪小旵.冷藏庫預(yù)冷降溫過程中溫度場的數(shù)值模擬與試驗研究[J].西北農(nóng)林科技大學(xué)學(xué)報:自然科學(xué)版,2008,36(9):219～223.

[4] 劉妍玲,張巖,王世清等.果蔬擺放形式對冷庫內(nèi)氣流場分布影響的研究[J].青島農(nóng)業(yè)大學(xué)學(xué)報:自然科學(xué)版,2008,25(1),24～27.

[5] 謝晶,瞿曉華,徐世瓊.冷藏庫內(nèi)氣體流場數(shù)值模擬及驗證[J].農(nóng)業(yè)工程學(xué)報,2005,21(2):11～16.

[6] H u Z H,Sun DW.CFD simu lation of heat andm oisture transfer for predicting cooling rate and w eight lossof cooked ham during air -b last chilling process[J].Journal of Food Engineering,2000,46(3):189～197.

[7] Cotella G,,M anzan M,Com ini G.CFD sim ulation of refrigerated disp lay cabinets[J].Internat ional Journal of Refrigerat ion, 2001,24(3):250～260.

[8] Guyot S,MarneyN,Sanoner P,et al.Direct thiolysis on crude applematerials for high-performance liquid chromatog raphy characterization and quantification of polyphenols in cider app le tissues and juices[J].M ethods Enzym ol、2001(335):57～70.

[9] 王福軍.計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004:122～123,212.