王曉春 徐萌 蘇釗 丁爽 龔慧峰 吳濤

（1.青島海爾智能技術(shù)研發(fā)有限公司，山東青島 266101；2.數(shù)字化家電國家重點實驗室，山東青島 266101；3.海爾智家股份有限公司，山東青島 266101；4.青島海爾空調(diào)器有限總公司，山東青島 266101）

掛機空調(diào)室內(nèi)機送風葉輪為貫流風機，其特點為沿著軸線方向送風相對均勻、噪聲低、結(jié)構(gòu)緊湊。貫流風機理論并不復(fù)雜，但由于其結(jié)構(gòu)的獨特性及高轉(zhuǎn)速特性，導(dǎo)致使用CFD 研究三維計算規(guī)模較大[1]。

對于制冷系統(tǒng)仿真而言，有效地獲取截面的風量分配是獲取換熱量、制冷劑壓降計算的關(guān)鍵因素之一，如何快速有效地獲取表面風量分布決定著制冷系統(tǒng)仿真應(yīng)用的效率。已有的研究中大多探索了三維貫流風機或者二維貫流風機流動特性[2,3]，通過仿真結(jié)合試驗的方法[4]研究風機葉片優(yōu)化對風量、噪聲的影響因素[5,6]，并沒有闡述三維模型及二維模型計算的差異性。

本文在前人仿真試驗對標模型的基礎(chǔ)上，利用焓差室風量測試數(shù)據(jù)做依據(jù)，分析了不同轉(zhuǎn)速等條件下，總風量仿真的可靠性及三維和簡化二維模型總風量、風量分布情況的可行性，為以二維風量分配計算在制冷系統(tǒng)仿真業(yè)務(wù)開展提供了依據(jù)。

1 仿真模型建立及可靠性驗證

1.1 仿真模型

為了確保仿真與試驗的一致性，在三維模型中，保留了完整的葉輪結(jié)構(gòu)、進風格柵、出風導(dǎo)葉片等。翅片采用多孔介質(zhì)模型替代，其慣性阻力系數(shù)和粘性阻力系數(shù)由風洞實驗室測試數(shù)據(jù)（如表1 所示）處理得到慣性阻力系數(shù)為116 kg/m4，粘性阻力系數(shù)314 kg/m3·s。

表1 風洞P-Q 測試數(shù)據(jù)

湍流模型采用k-e，標準壁面函數(shù)，非定常仿真動靜交界面采用Motion 網(wǎng)格技術(shù)處理，時間步長設(shè)定為1e-05 秒。為了捕捉到光滑的葉輪表面信息，研究了表面網(wǎng)格控制尺度在0.15 mm～0.5 mm，各尺度下葉輪表面網(wǎng)格形態(tài)如圖1，模型選擇0.15 mm 作為最小表面網(wǎng)格尺寸。

圖1 不同尺度表面網(wǎng)格對葉片表面光滑度影響

1.2 網(wǎng)格無關(guān)性驗證及精度確保

本文以高頻轉(zhuǎn)速1 300 rpm、非定常三維仿真開展網(wǎng)格無關(guān)性研究，受限于計算資源及效率，探索了5 種網(wǎng)格數(shù)量的流量計算值差異，圖2 為不同網(wǎng)格尺度下的計算結(jié)果，采用6 130 萬網(wǎng)格與試驗值834 m3/h 相比，誤差為2.73%，模型處于較高精度，驗證了計算結(jié)果的可靠性。

圖2 網(wǎng)格無關(guān)性研究

后續(xù)以此網(wǎng)格配置為基礎(chǔ)開展相關(guān)研究。圖3 為該網(wǎng)格尺度下葉輪的Y+分布情況。

圖3 Y+分布

2 結(jié)果對比分析及應(yīng)用

2.1 三維仿真不同轉(zhuǎn)速可靠性分析

基于上述三維非定常模型及多孔介質(zhì)模型，計算了900 rpm、1 300 rpm、1 606 rpm 三種轉(zhuǎn)速下的總風量、蒸發(fā)器表面風量分配情況，表2 為相應(yīng)的仿真與試驗統(tǒng)計信息，整體偏差在8%以內(nèi)，判定仿真計算有效。

表2 不同轉(zhuǎn)速三維瞬態(tài)仿真試驗對比

本文研究的目的是論證服務(wù)于制冷系統(tǒng)仿真的換熱器表面分配用二維替代三維進行CFD計算的可靠性，為達成這一目標，將換熱器表面按照圖4 虛切線的形式將換熱器進風側(cè)拆分為17段，葉輪展向作為一整段取數(shù)統(tǒng)計，換熱器進風表面風速分布如圖5 所示，葉輪展向分布基本一致，論證了后續(xù)二維界面計算的可行性。

圖4 風量分段統(tǒng)計示意圖

圖5 橫向進風風速分布云圖

統(tǒng)計不同轉(zhuǎn)速下流經(jīng)各段的風量占比，統(tǒng)計結(jié)果如圖6 所示。計算結(jié)果表明，不同轉(zhuǎn)速下分流分配基本一致，在做整機系統(tǒng)仿真時可以只做額定轉(zhuǎn)速的風量分配，其他轉(zhuǎn)速參照額定轉(zhuǎn)速分配，進而提高仿真計算效率，統(tǒng)計還發(fā)現(xiàn)，掛機三段式蒸發(fā)器，前中后的總風量占比分別約為25%、45%、30%。

圖6 不同轉(zhuǎn)速下三維仿真計算風量分配分布圖

2.2 二維替代三維可行性分析

在2.1 節(jié)的基礎(chǔ)上，補充計算了二維仿真計算，并將結(jié)果分別與試驗值、三維結(jié)果做對比，表3 為相關(guān)統(tǒng)計數(shù)據(jù)。

表3 不同轉(zhuǎn)速二維與三維分析數(shù)據(jù)

結(jié)果表明二維仿真相較實測，不同轉(zhuǎn)速下（900 rpm、1 300 rpm、1 606 rpm）總風量差異-7.58%、1.04%、4.53%；二維相較三維，不同轉(zhuǎn)速下（900 rpm、1 300 rpm、1 606 rpm）風量差異0.37%、-1.74%、-2.23%，二維仿真模型計算結(jié)果與三維相比差異不大，可以很好地預(yù)測總風量。從速度場截面看（如圖7 所示），風速分布位置、偏心渦位置也沒有差異。

圖7 不同仿真模式下的截面速度云圖

圖8 展示了不同轉(zhuǎn)速下二維與三維關(guān)于17 段風量分配的仿真計算結(jié)果，相較三維，二維計算的各段風量占比偏差統(tǒng)計值在±0.7%以內(nèi)，可用二維替代三維做風量分配計算。

圖8 不同轉(zhuǎn)速二維與三維仿真計算的風量分配

綜上所述，使用二維仿真替代三維進行總風量及風量分配計算具備可行性。

2.3 結(jié)果應(yīng)用

本研究的初衷是為空調(diào)一維制冷系統(tǒng)仿真在掛機應(yīng)用中提供一種簡便可靠的風量及分布方法，關(guān)于制冷系統(tǒng)仿真模型搭建本文不做具體介紹，只將研究結(jié)果進行應(yīng)用效果展示。

利用上述研究方法，在1.5 P φ7 蒸發(fā)器上進行分流設(shè)計，通過制冷劑環(huán)路的溫度分布進行應(yīng)用效果驗證。換熱器分流設(shè)計的理論設(shè)計圖與管路焊接實物圖如圖9 所示，制冷劑上下兩路從進口到出口各布置5 個溫度點。判斷分流設(shè)計優(yōu)劣的一個重要指標是換熱器不同流路出口溫度分布的一致性，制冷劑環(huán)路不同位置的溫度分布對比如表4 所示。仿真數(shù)據(jù)顯示，在該工況下，只有蒸發(fā)器出口的每一路最后一只管出現(xiàn)過熱，且上下兩路出口的溫度分別為14.00℃和14.20℃，處于較好的分流設(shè)計狀態(tài)。將該分流設(shè)計制作樣件測試，試驗數(shù)據(jù)也是每一路最后一只管處于過熱狀態(tài)，出口溫度高度一致，差值僅為0.15℃，很好地驗證了應(yīng)用二維仿真替代三維進行風量及分布計算的可行性。通過優(yōu)化分流，實現(xiàn)了換熱器換熱量從3 450 W 提升到3 563 W，換熱量提升3.3%。

圖9 某型號換熱器分流設(shè)計圖

表4 制冷劑環(huán)路不同位置試驗與仿真溫度對比（單位：℃）

3 結(jié)束語

本文通過不同轉(zhuǎn)速下的仿真與試驗對比，從總風量、分布進行仿真試驗對比，論證了二維替代三維計算風量及進風表面分布的可靠性，得出以下結(jié)論：

（1）調(diào)試的模型及網(wǎng)格策略在900 rpm、1 300 rpm、1 606 rpm 三種轉(zhuǎn)速下的總風量及分配的仿真試驗對比偏差在8%以內(nèi)，仿真計算精度可靠。

（2）二維仿真模型計算結(jié)果與三維相比差異不大，偏差在3%以內(nèi)，可以很好地預(yù)測總風量，17 段風量分配的仿真計算結(jié)果表明，相較三維，二維計算的各段風量占比偏差統(tǒng)計值在±0.7%以內(nèi)，可用二維替代三維做風量分配計算。

應(yīng)用二維仿真替代三維進行風量及分布計算的方法，通過制冷系統(tǒng)各路分流的表面溫度分布情況驗證了方法的可靠性。通過優(yōu)化分流，實現(xiàn)了換熱器換熱量從3 450 W 提升到3 563 W，換熱量提升3.3%。

本文僅對進風總風量進行了試驗對比研究，且從實際產(chǎn)品應(yīng)用的角度判定了二維CFD 計算風量分配在制冷系統(tǒng)仿真中的可靠性。受限于篇幅，出風口二維及三維的仿真對比分析未在本文展開討論。后續(xù)可以通過PIV 等測試方法，完善進出風口仿真與試驗的可靠性研究。