沈　希，張　舒，徐　鳴，黃躍進，顧江萍

(浙江工業(yè)大學(xué) 機械工程學(xué)院，浙江 杭州 310014)

保護器測試系統(tǒng)恒溫箱的CFD研究與實驗驗證

沈希，張舒，徐鳴，黃躍進，顧江萍

(浙江工業(yè)大學(xué) 機械工程學(xué)院，浙江 杭州 310014)

摘要：為了測試壓縮機保護器在特定溫度的運行情況，需要對溫度進行必要的控制.針對保護器溫度特性檢測系統(tǒng)中的恒溫箱均勻度差的問題，建立了恒溫箱流體域CFD數(shù)值模型，得到了恒溫箱的溫度場、速度場模型.恒溫箱內(nèi)5個測溫點溫度數(shù)據(jù)與模擬數(shù)據(jù)基本吻合，驗證了該數(shù)值模型的正確性.分析結(jié)果表明：送風(fēng)溫差適當(dāng)?shù)慕档陀欣跍囟染鶆蛐缘母纳?，同時送風(fēng)速度及恒溫箱結(jié)構(gòu)也對恒溫箱均勻度有很大的影響.

關(guān)鍵詞：恒溫箱；溫度場；均勻度；數(shù)值模擬

CFD modelling and experimental verification on thermostat

of protector testing system

SHEN Xi, ZHANG Shu, XU Ming, HUANG Yuejin, GU Jiangping

(College of Mechanical Engineering, Zhejiang University of Technology, Hangzhou 310014, China)

Abstract:In order to test the operating conditions of the compressor protector under specific temperature protection， necessary control of the temperature is needed. Aiming at problem of average degree of the temperature characteristic test system of protector，the CFD numerical model of thermostat is established to get temperature field and velocity field .Using Delphi for upper machine, data collected from five points on the test surface matched the simulated data. It verifies the correctness of the numerical model. Analysis results show that appropriate reducing of supply air temperature is conducive to the improvement of average degree and the supply air’s velocity and the structure of thermostat also have great influence on incubator evenness.

Keywords:thermostat; temperature field; average degree; numerically simulation

目前，大部分壓縮機選配的熱保護器都是圓盤型過熱過流保護裝置.此結(jié)構(gòu)中的電觸點，利用碟片的動作控制電路通道.保護器串聯(lián)到一個單相電動機公共端子，電流流過保護器的發(fā)熱絲和雙金屬片產(chǎn)生了熱量.雙金屬片的溫度既受電機繞組電流的影響也受壓縮機殼體溫度的影響，保護器兼有過熱和過流的保護功能.因為保護器的溫度特性關(guān)乎制冷設(shè)備的安全性，所以動作溫度和恢復(fù)溫度便成了極其重要的參數(shù)[1].在恒溫箱的溫度測試系統(tǒng)設(shè)計上，已有很多成果.高峰以提高溫度控制的精度，減小溫度控制的過渡過程時間，降低測控系統(tǒng)的硬件成本為目標(biāo)，進行了優(yōu)化設(shè)計[2]，張雷杰等通過CFD軟件對恒溫箱內(nèi)的兩種氣流組織進行了模擬和實驗驗證[3]，周曉設(shè)計出了無線傳感網(wǎng)絡(luò)用于采集溫度并上傳[4].在數(shù)值模擬方面，張建勝采用了帶旋流修正的RealizableK—ε湍流模型對風(fēng)壓進行了分析[5]，童躍平等對離心泵的內(nèi)部流場流動情況進行了兩種模型的對比分析[6].

在以往的溫度測試中，我們主要關(guān)注的是測試過程和測試結(jié)果是否滿足規(guī)范要求，對于影響測試結(jié)果的溫度均勻度等因素往往有些忽視[7].因此，筆者以一臺保護器測試裝置的恒溫箱作為研究對象，采用fluent6.3軟件，建立恒溫箱流場三維參數(shù)化實體模型，對恒溫箱的溫度場與速度場進行了分析，對數(shù)值模擬的結(jié)果與實驗進行驗證.

1恒溫箱流場數(shù)值模型研究

1.1基本控制方程

假設(shè)流體為不可壓縮的牛頓型流體且考慮常物性和宏觀熱能守恒的，恒溫箱內(nèi)流體的流動和傳熱必須滿足下面3個方程和Boussinesq模型[8-9].

質(zhì)量守恒方程為

(1)

動量守恒方程(運動方程)為

(2)

(3)

(4)

能量守恒方程為

(5)

Boussinesq模型為

(ρ-ρ0)g≈-ρ0β(T-T0)g

(6)

式中：u，v，w分別為流體在x，y，z三個方向上的速度分量；ρ為流體密度；μ為動力粘度；Su，Sv，Sw分別為動量守恒方程的廣義源項；T為流體溫度；λ為流體的導(dǎo)熱系數(shù)；CP為比熱容；ST為流體的內(nèi)熱源及粘性作用流體機械能轉(zhuǎn)化為熱能的部分，簡稱粘性耗散項；ρ0為流體的常密度；T0為操作溫度；β為熱膨脹系數(shù).

1.2計算模型

對于模型的建立先參考實際的恒溫箱畫出模型如圖1所示，再導(dǎo)入到前處理軟件gambit中，采用非結(jié)構(gòu)化的三角形網(wǎng)格,如圖2所示.建模所采用的恒溫箱的具體尺寸由于篇幅太長不再贅述，計算工況如表1所示.

圖1　恒溫箱模型圖Fig.1　Modle of thermostat

圖2　恒溫箱網(wǎng)格圖Fig.2　Grid of thermostat

名稱參數(shù)名稱參數(shù)流體空氣熱膨脹系數(shù)0.0187墻壁鋼鐵重力加速度/(m·s-2)9.81初始溫度/℃45比熱容/(J·kg-1·℃-1)1012葉輪轉(zhuǎn)速/(rad·s-1)293加熱管功率/(kW·m-3)769

1.3模型計算與邊界條件的簡化

由于恒溫箱尺寸跨度比較大，離心葉輪的結(jié)構(gòu)相對復(fù)雜，所以在物理量基本準(zhǔn)確的前提下，模型的建立可進行如下簡化假設(shè)：1) 忽略上下孔板結(jié)構(gòu)；2) 忽略離心葉輪的扇葉的具體尺寸；3) 假設(shè)流體流動和傳熱過程是均一、穩(wěn)定的；4) 流體為牛頓流體、不可壓縮、各向同性且連續(xù)[10]；5) 不考慮箱體壁厚；6) 只考慮加熱情況，不考慮制冷情況.

計算邊界條件的簡化：1) 由于模型中的流到屬于靜網(wǎng)格，扇葉屬于動網(wǎng)格，所以在考慮采用MRF多重參考系統(tǒng)進行耦合(Multi-reference.Frame)；2) 空氣出口采用出流(outflow)；3) 加熱管加熱采用共軛傳熱，加熱管區(qū)域設(shè)置為熱源.

本模型屬于K—ε湍流，計算中選用模型進行湍流模擬.求解的條件用非耦合求解、隱式算法、三維空間、非定常流動與絕對速度；控制方程的離散采用一階迎風(fēng)式，速度與壓力的耦合采用SIMPLE算法，其他保持默認設(shè)置[11-12].

2計算結(jié)果與分析

圖3是計算260 s后恒溫箱Y=0 m中截面等速度分布圖，從圖3中可以看出：在恒溫腔體的左側(cè)存在著渦流區(qū)域，這是由于腔體左側(cè)聚集了能量較低的流體，而腔體流道中存在逆壓梯度，這樣在左側(cè)就產(chǎn)生了氣體的回流.就因為這樣渦流區(qū)域的存在，影響了氣流向左側(cè)的移動，從而對恒溫箱內(nèi)氣流溫度與速度分布的均勻性產(chǎn)生影響.

圖3　Y=0 m等速度分布圖Fig.3　The Velocity profile of Y=0 m

圖4是Z=0 m的中截面速度矢量圖，由圖4可以看出：葉輪流場的矢量圖明顯為非對稱，盡管恒溫箱總體在結(jié)構(gòu)上有對稱性，但在圖4下方正方形缺口處存在一個需要給觸點接線的金屬板結(jié)構(gòu)，該結(jié)構(gòu)破壞了恒溫箱整體的對此性.在靠近金屬板的一側(cè)產(chǎn)生了速度較大的射流，這使得氣流在循環(huán)過程剛開始出現(xiàn)了不對稱.從上面的分析可以預(yù)見，要想優(yōu)化恒溫箱溫度均勻性，可以將觸點板移至外部，從而改變氣流組織條件.

圖5是測試平面的溫度圖，由仿真結(jié)果可以看出：溫度右邊的要明顯大于左邊，且根據(jù)多組數(shù)據(jù)觀測，右下角與中間區(qū)域位置的溫度最高.產(chǎn)生這樣的原因就是因為氣流從頂部回到葉輪進氣口過程中，產(chǎn)生了比較大的渦流，使得測試板左邊空氣不能正常的進行換熱，即產(chǎn)生了溫度死區(qū).

圖4　Z=0 m速度矢量圖Fig.4　The velocity vector diagram of Z=0 m

圖5　測試平面溫度圖Fig.5　The image of test surface temperature

3實驗驗證與分析

實驗過程中采用研華工控機為控制單元，PCL-730卡為I/O卡，PCL-813卡為A/D采樣卡，溫度傳感器采用PT100，上位機軟件由Delphi編寫，采用access 2007數(shù)據(jù)庫對需要采集的數(shù)據(jù)進行自動存儲，可以避免人工的查詢與記錄.

測試平面溫度傳感器布線方式如圖6所示，在模型中進行多次仿真然后與結(jié)果相對比，圖7為工況一的實驗與模擬對比圖，經(jīng)過多次實驗對比發(fā)現(xiàn)，測試點4處的模擬值與實際值偏差較大，測試平面的總體溫度趨勢與實際的溫度趨勢吻合的較好，故該模型可用于實際恒溫箱內(nèi)部溫度場的預(yù)測.為了研究葉輪的轉(zhuǎn)速，加熱管功率等在恒溫箱內(nèi)的影響，在表1仿真基礎(chǔ)上設(shè)置幾組邊界條件進行模擬如表2所示.由工況一、二可以看到，當(dāng)增加葉輪轉(zhuǎn)速，即增大送風(fēng)量可明顯改善溫度均勻性；由工況一、三可以看出，當(dāng)減小加熱管功率，即減小送風(fēng)溫差也能明顯改善溫度場的均勻性.

在實驗中，送風(fēng)溫差取決于加熱管與制冷機組熱平衡后的熱負荷，送風(fēng)量取決于風(fēng)機的轉(zhuǎn)速.采用熱平衡方式，即制冷機組常開，用加熱管平衡多余的制冷量可以減小送風(fēng)溫差；還可以對加熱管開關(guān)量進行更加精細的PID連續(xù)調(diào)節(jié)，這樣可以減小溫度波動從而增加溫度均勻性；對于增加送風(fēng)量，在試驗標(biāo)準(zhǔn)允許條件下提高風(fēng)速，可以加強空氣在箱內(nèi)流動，消除死區(qū)，從而提高溫度均勻度.

圖6　溫度探頭布置圖Fig.6　The layout of temperature sensor

圖7　工況一的實驗與模擬對比圖Fig.7　The experimental contrast figure of condition 1

邊界條件工況一工況二工況三葉輪轉(zhuǎn)速/(r·s-1)300100300加熱管功率/(kW·m-3)800800400觀測面最大溫差/℃2.44.62.0

4結(jié)論

根據(jù)實際恒溫箱的具體尺寸、工況建立了模型，對其進行了數(shù)值模擬和實驗驗證，利用fluent軟件對恒溫箱的整體進行了數(shù)值計算，揭示了箱體內(nèi)速

度場和溫度場的分布特點.因為網(wǎng)格生成的難度，所以模型中簡化了上下孔板的結(jié)構(gòu)，過程中也沒有考慮扇葉的具體尺寸.溫度實驗值與模擬值最大偏差在5%內(nèi)[13]，總體的模擬結(jié)果與實驗測得吻合良好，溫度分布符合實際情況.由仿真結(jié)果可以看出：減小送風(fēng)溫差和提高風(fēng)循環(huán)速度都可以提高恒溫箱內(nèi)的溫度均勻性，但從實際角度考慮，減小送風(fēng)溫差的方法更加經(jīng)濟可行，這就需要在溫度控制上從原來的PID控制轉(zhuǎn)為更加精細的控制方法.由于恒溫箱無法做到結(jié)構(gòu)上完全對稱，在靠近離心葉輪的一方產(chǎn)生了速度較大的射流，氣流進入測試環(huán)境后，左邊氣流速度大于右邊，出現(xiàn)溫度不均勻，為了彌補結(jié)構(gòu)上的缺陷，可在恒溫箱上方添加擋風(fēng)板和交叉條，從而增加結(jié)構(gòu)的對稱性.

參考文獻：

[1]中國國家標(biāo)準(zhǔn)化管理委員會.GB 4706.17—2010家用和類似用途電器的安全 電動機-壓縮機的特殊要求[S].北京：中國標(biāo)準(zhǔn)出版社，2010.

[2]高峰.高精度恒溫箱測控系統(tǒng)的優(yōu)化設(shè)計[D].合肥：合肥工業(yè)大學(xué),2010.

[3]張雷杰,張敏,楊樂.兩款恒溫箱的氣流組織模擬與實驗驗證[J].安徽農(nóng)業(yè)科學(xué),2009,37(18):8752-8764.

[4]周曉,李杰,邊裕挺.基于無線傳感網(wǎng)絡(luò)的環(huán)境溫度監(jiān)測系統(tǒng)設(shè)計[J].浙江工業(yè)大學(xué)學(xué)報,2013，41(4)：440-443.

[5]張建勝,周鋒,盧成原.典型低矮雙坡屋蓋平均風(fēng)壓的數(shù)值模擬[J].浙江工業(yè)大學(xué)學(xué)報,2014，42(3)：302-306.

[6]童躍平,張淑佳,李賢華.標(biāo)準(zhǔn)K—ε模型與RSM模型在離心泵三維模擬中的比較[J].浙江工業(yè)大學(xué)學(xué)報,2008，36(6)：678-681.

[7]陳謀義.氣候環(huán)境試驗箱(室)溫度均勻性問題及其改進[J].環(huán)境技術(shù),2001(5):27-30.

[8]古新,董其伍,王珂.基于三維實體模型的管殼式換熱器殼程流場和溫度場數(shù)值研究[J].壓力容器,2008,25(5):1-5.

[9]陶文銓.數(shù)值傳熱學(xué)[M].西安：西安交通大學(xué)出版社,1988.

[10]付磊,付麗婭,唐克倫,等.基于FLUENT的管殼式換熱器殼程流場數(shù)值模擬研究[J].四川理工學(xué)院學(xué)報,2012(3):17-21.

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

[12]韓占忠,王敬,蘭小平.FLUENT流體工程仿真計算實例與應(yīng)用[M].北京：北京理工大學(xué)出版社,2004.

[13]胡宗,楊文杰,朱柏山.某會展中心報告廳空調(diào)氣流組織的CFD模擬分析[J].重慶建筑,2012(11):58-60.

(責(zé)任編輯：劉巖)

中圖分類號：TP23

文獻標(biāo)志碼：A

文章編號：1006-4303(2015)03-0237-04

作者簡介：沈希(1966—),男,浙江杭州人,教授,主要從事熱物理量的測量和控制研究,E-mail:sx@zjut.edu.cn.

基金項目：國家自然科學(xué)基金資助項目(51076143)

收稿日期：2014-12-19