張旭 康釗 胡松濤 朱輝 陳紀(jì)敏 張鳳嬌

1.海信空調(diào)有限公司 山東青島 266100；2.青島理工大學(xué)環(huán)境與市政工程學(xué)院 山東青島 266520

0 引言

全球性溫室效應(yīng)的發(fā)展帶來了一系列氣候問題，尤其是近年來愈演愈烈的熱浪事件（heat waves），給人們的健康和生命安全帶來了嚴(yán)重威脅[1]。在全球變暖和極端氣候事件頻發(fā)的背景下，空調(diào)對(duì)于人們的日常生活和生產(chǎn)尤為重要，尤其是環(huán)境脆弱性人群，如老人、兒童、病人等。因此，空調(diào)系統(tǒng)的研究和優(yōu)化依然是領(lǐng)域內(nèi)研究的重點(diǎn)。

壁掛式空調(diào)作為常見的家用空調(diào)，被很多普通家庭所使用。在夏季和冬季空調(diào)的使用率一般會(huì)急劇上升，所以空調(diào)房間的氣流組織問題越來越引起人們的關(guān)注。壁掛式空調(diào)的氣流組織效果對(duì)空調(diào)性能（如熱舒適、節(jié)能性）起著決定性的作用，因此通過優(yōu)化空調(diào)氣流組織的仿真來提高空調(diào)性能具有重要意義。

空調(diào)器氣流組織模擬仿真是預(yù)測和優(yōu)化空調(diào)器熱舒適性和節(jié)能性的重要手段，在早期設(shè)計(jì)階段可以進(jìn)行仿真預(yù)測并優(yōu)化設(shè)計(jì)方案，提高研發(fā)效率，節(jié)省試驗(yàn)成本。如早在1998年Murakami[2]等學(xué)者便開始采用計(jì)算流體動(dòng)力學(xué)（CFD）技術(shù)進(jìn)行了空調(diào)熱負(fù)荷的計(jì)算，并在后續(xù)研究中對(duì)動(dòng)態(tài)熱負(fù)荷進(jìn)行了數(shù)值模擬，為后續(xù)CFD技術(shù)在空調(diào)行業(yè)的發(fā)展應(yīng)用提供了基礎(chǔ)。近年來，有學(xué)者對(duì)我國東部地區(qū)日益增多的盒式空調(diào)通風(fēng)方式進(jìn)行了調(diào)查，并通過CFD技術(shù)分析了其氣流組織形式同時(shí)簡化了模型[3-4]。在空調(diào)的熱舒適性方面，有學(xué)者研究了壁掛式空調(diào)制熱的舒適問題，通過改變空調(diào)的送風(fēng)溫度來改變熱舒適，并提出雖然較高的送風(fēng)溫度可以提高系統(tǒng)的制熱能力，但是由于高風(fēng)量引起的氣流組織失穩(wěn)，室內(nèi)熱舒適的改善效果有限[5-8]。

在送風(fēng)角度方面國內(nèi)外學(xué)者也開展了一系列研究工作，發(fā)現(xiàn)壁掛式空調(diào)器的送風(fēng)導(dǎo)流板的設(shè)計(jì)較大地影響著空調(diào)的熱舒適性，尤其是在制冷工況下。劉丙磊等對(duì)單層和雙層導(dǎo)流板進(jìn)行了對(duì)比，發(fā)現(xiàn)雙層導(dǎo)流板可以有效地降低送風(fēng)噪聲[9]。此外對(duì)于導(dǎo)流片結(jié)構(gòu)和葉輪轉(zhuǎn)速的數(shù)值模擬發(fā)現(xiàn)，導(dǎo)流板一方面可以改變主流方向和分散射流速度軸心區(qū)，使室內(nèi)流場和溫度場更加均勻；另一方面當(dāng)導(dǎo)流板角度不合適時(shí)可能阻礙送風(fēng)[10-11]。

通過總結(jié)現(xiàn)有文獻(xiàn)資料發(fā)現(xiàn)，針對(duì)貫流風(fēng)扇室內(nèi)掛機(jī)目前行業(yè)內(nèi)的熱舒適性仿真方法主要有三維整機(jī)仿真、二維風(fēng)道仿真、等效出風(fēng)仿真三種，但均存在室內(nèi)流場分布誤差較大的問題（統(tǒng)一采用垂直風(fēng)向，室內(nèi)風(fēng)速，溫度失真），使仿真數(shù)據(jù)缺乏參考性。目前在空調(diào)器熱舒適性仿真領(lǐng)域，如何在保留貫流風(fēng)機(jī)、完整風(fēng)道、實(shí)驗(yàn)室尺寸等比例建模的前提下高精度預(yù)測空調(diào)器的熱舒適性能參數(shù)（包括風(fēng)速、溫度、濕度）是行業(yè)難題。

綜上所述，目前針對(duì)貫流風(fēng)機(jī)空調(diào)器室內(nèi)氣流組織缺乏一套成熟的仿真方法。上述三種相關(guān)的仿真方法均未能計(jì)算出高精度的室內(nèi)熱舒適性流場分布情況，其主要原因在于計(jì)算域的邊界條件理論研究欠缺；網(wǎng)格劃分問題需系統(tǒng)優(yōu)化；算法適配問題亟待解決，例如：換熱器邊界溫度、風(fēng)速不均分布的附加，濕度的附加、耗散的影響等。因此，本研究從以上問題出發(fā)，采用計(jì)算流體力學(xué)技術(shù)對(duì)某型號(hào)壁掛式貫流風(fēng)機(jī)空調(diào)器的室內(nèi)氣流組織進(jìn)行全尺度的數(shù)值模擬，旨在解決傳統(tǒng)模擬方法的一系列局限和問題，為新產(chǎn)品的開發(fā)和優(yōu)化提供參考。

1 研究方法

模擬對(duì)象為某標(biāo)準(zhǔn)熱舒適性實(shí)驗(yàn)室，該實(shí)驗(yàn)室的結(jié)構(gòu)如圖1所示，其中門為木制門，窗為雙層玻璃窗，墻體為標(biāo)準(zhǔn)的混凝土磚墻。其中輻射板僅在該實(shí)驗(yàn)室進(jìn)行輻射供冷/供熱試驗(yàn)時(shí)使用。本研究對(duì)空調(diào)掛機(jī)與房間進(jìn)行全尺度等比例建模，并導(dǎo)入CFD軟件進(jìn)行仿真計(jì)算。

圖1 模擬對(duì)象

1.1 網(wǎng)格劃分方案

在網(wǎng)格劃分過程中，對(duì)于房間內(nèi)的流體域采用六面體網(wǎng)格進(jìn)行劃分，但是由于房間空調(diào)器遠(yuǎn)端不需要重點(diǎn)觀察，因此房間流體域的網(wǎng)格間距按照1.05倍的比例逐漸變大（如圖2所示）?？紤]到貫流風(fēng)機(jī)為整個(gè)數(shù)值模擬的動(dòng)力源，該部分的流場涉及到整個(gè)仿真的準(zhǔn)確性，需要重點(diǎn)關(guān)注。因此在網(wǎng)格劃分過程中，對(duì)于貫流風(fēng)機(jī)和換熱翅片進(jìn)行網(wǎng)格加密處理。由于貫流風(fēng)機(jī)的邊界極其復(fù)雜，因此采用四面體網(wǎng)格以減少網(wǎng)格劃分工作量，而換熱器部分采用六面體網(wǎng)格（如圖3所示）。對(duì)于空調(diào)器內(nèi)除了貫流風(fēng)機(jī)外的其他流體域，采用六面體網(wǎng)格進(jìn)行劃分。最后生成網(wǎng)格數(shù)量為690萬，網(wǎng)格質(zhì)量滿足CFD軟件的計(jì)算要求。

圖2 房間內(nèi)流體域的網(wǎng)格劃分方案

圖3 換熱器和貫流風(fēng)機(jī)的網(wǎng)格劃分方案

1.2 湍流模型與離散方法

本數(shù)值模擬主要針對(duì)受限空間內(nèi)部三維穩(wěn)態(tài)流動(dòng)問題?？紤]到計(jì)算效率問題排除RSM、LES模型，最終選擇SST k-Ω湍流模型。SST k-Ω模型同時(shí)適用于轉(zhuǎn)捩流和低雷諾數(shù)流動(dòng)，可以針對(duì)性地求解本問題中的貫流風(fēng)機(jī)旋流、室內(nèi)空間的渦旋以及室內(nèi)空間的部分區(qū)域（如空調(diào)遠(yuǎn)端、角落等）的低速流體運(yùn)動(dòng)，適用于本數(shù)值模擬的對(duì)象特征[12]。其次，對(duì)于本數(shù)值模擬針對(duì)的穩(wěn)態(tài)問題，模擬過程中選擇壓力速度半隱式耦合算法SIMPLEC來進(jìn)行求解[13]。SIMPLEC算法與經(jīng)典的SIMPLE算法的基本思路一致，僅在通量修正方法上有所改進(jìn)，因而加快了計(jì)算的收斂速度。

1.3 邊界條件

本模擬主要針對(duì)制冷工況，邊界條件設(shè)置如下：對(duì)于貫流風(fēng)機(jī)采用旋轉(zhuǎn)邊界，分別設(shè)置三個(gè)模擬工況下的轉(zhuǎn)速為600 r/min、1200 r/min和1500 r/min。

為了簡化網(wǎng)格劃分和計(jì)算，換熱器處采用多孔介質(zhì)模型來代替翅片管。多孔介質(zhì)設(shè)置模塊中采用慣性阻力系數(shù)和粘性阻力系數(shù)法，其中Y、Z方向慣性阻力系數(shù)和粘性系數(shù)為同一數(shù)量級(jí)。X方向上的慣性阻力系數(shù)和粘性阻力系數(shù)假設(shè)為無窮大，因?yàn)閾Q熱器翅片流道極大地限制了氣流在X方向上的擴(kuò)散。

將換熱器滲透率設(shè)置為1（全部通過），通過測試該制冷系統(tǒng)下的出風(fēng)口溫度指定通過換熱器的流體（空氣）溫度為17.3℃。設(shè)置壁面為常熱流密度（第二類熱邊界條件），熱流密度值源自于試驗(yàn)過程中的測量值：四周圍護(hù)壁面的熱流密度折算為30 W/m2，地面設(shè)置為絕熱壁面。此外，在參考條件中，對(duì)空氣密度、重力加速度進(jìn)行了相應(yīng)設(shè)置，并采用Boussinesq假設(shè)來精確模擬溫度驅(qū)動(dòng)的浮升力對(duì)流場和溫度場的作用[14]。

2 模擬結(jié)果

2.1 貫流風(fēng)機(jī)流場

圖4為不同導(dǎo)流葉片（送風(fēng)角度）情況下，空調(diào)器內(nèi)部的二維流場分布情況。由圖可知，在無導(dǎo)流片、3導(dǎo)流片和4導(dǎo)流片的情況下，空調(diào)器內(nèi)部流場分布相似，并且在貫流風(fēng)機(jī)的蝸舌處均出現(xiàn)了偏心渦。偏心渦是貫流風(fēng)機(jī)原理的重要特征，這說明空調(diào)內(nèi)部流場的數(shù)值模擬結(jié)果具有較好的精確度，可以保證動(dòng)力源處流場的精度，從而保證整個(gè)房間內(nèi)部氣流組織的準(zhǔn)確度。

圖4 貫流風(fēng)機(jī)流場

其次，從圖4可以看出，空調(diào)器出口處（送風(fēng)口）的流場因?qū)Я髌霈F(xiàn)一定的差異。對(duì)于無導(dǎo)流片的情況，空調(diào)器出口流體出現(xiàn)上揚(yáng)現(xiàn)象，易引起回風(fēng)短路，而對(duì)于3導(dǎo)流片和4導(dǎo)流片的情況，上揚(yáng)現(xiàn)象受到明顯抑制，并且導(dǎo)流片越多，抑制效果越顯著。這印證了空調(diào)器送風(fēng)口處設(shè)置導(dǎo)流片的必要性和重要性：出口導(dǎo)流葉片可以在一定程度上抑制短路的產(chǎn)生，可有效地將渦旋對(duì)稱線“下壓”，從而使得空調(diào)機(jī)正下方的渦旋負(fù)壓區(qū)較大，通過“拖曳”作用抑制回風(fēng)短路現(xiàn)象產(chǎn)生。

2.2 600 r/min的流場和溫度場

圖5為貫流風(fēng)機(jī)600 r/min時(shí)的三維數(shù)值模擬結(jié)果。由圖5 a) 可知，出風(fēng)口流體的水平射流距離達(dá)到了房間的50%左右的距離，可以保證處于房間中間位置的人員有一定的吹風(fēng)感以滿足熱舒適需求。但是流場的分層現(xiàn)象并不顯著，這可能導(dǎo)致整個(gè)室內(nèi)溫度場的集中程度較高，如圖5 b) 所示。室內(nèi)溫度出現(xiàn)了一定的分層現(xiàn)象，但是靠近壁面處具有較厚的熱邊界層。尤其在空調(diào)器遠(yuǎn)端的壁面處溫度分層不明顯，可能出現(xiàn)局部溫度偏高，影響熱舒適度。因此低風(fēng)量情況雖然可以抑制回風(fēng)短路，但是室內(nèi)溫度場的分布并不理想，適用于冷負(fù)荷較小的情況。最后從圖5 c) 可以直觀地發(fā)現(xiàn)，整個(gè)空調(diào)房間以出風(fēng)射流線為中心，形成了上下兩個(gè)接近堆疊的大渦旋，有效抑制了出風(fēng)口回風(fēng)短路的現(xiàn)象，但帶來的后果是溫度分層現(xiàn)象不理想，可能影響房間壁面附近的熱舒適性。

圖5 600 r/min時(shí)的模擬結(jié)果

2.3 1200 r/min的流場和溫度場

圖6 為貫流風(fēng)機(jī)轉(zhuǎn)速為1200 r/min時(shí)的氣流組織和溫度分布情況。由圖6 a) 可知，隨著貫流風(fēng)機(jī)轉(zhuǎn)速的增大，送風(fēng)的水平射流距離增大到房間遠(yuǎn)端處，并且出現(xiàn)了較明顯的速度分層現(xiàn)象。這說明隨著送風(fēng)速度的增加，室內(nèi)氣流的擾動(dòng)，即湍流度進(jìn)一步增大，利于更均勻的溫度場的發(fā)展。圖6 b) 印證了這一現(xiàn)象，由圖6 b) 可知，送風(fēng)量的增大可提高整體房間內(nèi)的溫度均勻度，同時(shí)也使得溫度分層現(xiàn)象更加明顯。即使在空調(diào)器的遠(yuǎn)端也有著較舒適的空氣溫度和送風(fēng)速度，有利于保證空調(diào)房間的整體熱舒適性。圖6 c) 為1200 r/min的轉(zhuǎn)速下的流線分布情況?？梢园l(fā)現(xiàn)送風(fēng)速度的增加在一定程度上破壞了低送風(fēng)速度下渦旋的對(duì)稱性，但是由于射流距離的增加，渦旋的尺度變得更大。

圖6 1200 r/min時(shí)的模擬結(jié)果

2.4 1500 r/min的流場和溫度場

圖7為貫流風(fēng)機(jī)轉(zhuǎn)速1500 r/min時(shí)的氣流組織和溫度分布情況。由圖7 a) 可知，送風(fēng)速度的增大進(jìn)一步引起了射流距離的增加，并且速度分層現(xiàn)象更加顯著。這表明送風(fēng)速度的增加引起了更大的湍動(dòng)能和耗散現(xiàn)象，有利于房間內(nèi)出現(xiàn)更加均勻的溫度。圖7 b) 為轉(zhuǎn)速1500 r/min送風(fēng)下的溫度分布情況，可以發(fā)現(xiàn)空調(diào)器遠(yuǎn)端壁面處出現(xiàn)了明顯的溫度分層現(xiàn)象，整個(gè)空調(diào)房間內(nèi)人員活動(dòng)區(qū)域的溫度場更加理想。因此，大送風(fēng)量更加適用于房間冷負(fù)荷較大的工況，如空調(diào)房間的壁面為外墻，具有較高的太陽輻射得熱等工況。圖7 c) 為轉(zhuǎn)速1500 r/min的送風(fēng)下室內(nèi)流體的流線圖。可以發(fā)現(xiàn)送風(fēng)速度的增加破壞了渦旋的對(duì)稱性，左下側(cè)大渦旋底部向房間遠(yuǎn)端擴(kuò)散，而右上側(cè)大渦的空間被壓縮。

圖7 1500 r/min時(shí)的模擬結(jié)果

3 模擬結(jié)果準(zhǔn)確度驗(yàn)證

為了驗(yàn)證模擬結(jié)果的準(zhǔn)確度，進(jìn)行了相同邊界條件下中風(fēng)速工況（轉(zhuǎn)速1200 r/min）下的對(duì)標(biāo)試驗(yàn)，并將該工況下的出口質(zhì)量流量、溫度場的試驗(yàn)數(shù)據(jù)和模擬數(shù)據(jù)進(jìn)行了對(duì)比。

3.1 質(zhì)量流量

本次仿真監(jiān)測了空調(diào)出風(fēng)口的流量，用以檢驗(yàn)仿真結(jié)果與實(shí)際情況的一致性。通過CFD軟件提取空調(diào)器出風(fēng)口的面積加權(quán)平均質(zhì)量流量為0.2057 kg/s。根據(jù)該模擬值計(jì)算得到穩(wěn)態(tài)工況下空調(diào)器的出口體積流量為V=604 m3/h，與風(fēng)量臺(tái)實(shí)測風(fēng)量578 m3/h的相對(duì)誤差約為4.5%。

3.2 溫度場

對(duì)標(biāo)試驗(yàn)過程中，在實(shí)驗(yàn)室長、寬、高三個(gè)方向上分別布置了7×6×5只熱電偶，共210只，實(shí)時(shí)連續(xù)獲取穩(wěn)態(tài)制冷工況下室內(nèi)的溫度分布（如圖8所示）。模擬過程中，在CFD軟件中同樣設(shè)置了對(duì)應(yīng)的210個(gè)點(diǎn)，用于獲取穩(wěn)態(tài)制冷時(shí)室內(nèi)各個(gè)位置的溫度分布情況，如圖9所示。

圖8 實(shí)驗(yàn)室測點(diǎn)分布圖

圖9 模擬過程中設(shè)置的溫度點(diǎn)

本研究中分別獲得210個(gè)溫度實(shí)測值與模擬值。兩者的誤差計(jì)算方法為：（模擬值－實(shí)測值）/實(shí)測值。210個(gè)測點(diǎn)的總體誤差分布如圖10所示。由圖10可知，模擬結(jié)果和試驗(yàn)實(shí)測結(jié)果偏差度基本分布在5%左右，但是也存在局部數(shù)據(jù)異常點(diǎn)，數(shù)量為4個(gè)，約占210個(gè)溫度點(diǎn)數(shù)量的2%，幾乎可以忽略不計(jì)，具體的結(jié)果對(duì)比如表1所示。

表1 1200 r/min工況下的熱電偶溫度測點(diǎn)與對(duì)應(yīng)位置模擬溫度點(diǎn)的相對(duì)誤差值（單位：%）

圖10 模擬值與實(shí)測值誤差分布

表1給出了所有210個(gè)溫度點(diǎn)的實(shí)測值與模擬值的誤差情況，由表可知，制冷工況、1200 r/min工況條件下，數(shù)值模擬和試驗(yàn)測試所獲得的210個(gè)溫度點(diǎn)中僅有4個(gè)溫度異常點(diǎn)（誤差超過10%的溫度點(diǎn)）出現(xiàn)。通過核查試驗(yàn)結(jié)果，其中3處異常點(diǎn)可能是熱電偶發(fā)生故障欠校準(zhǔn)導(dǎo)致的，具體表現(xiàn)為該處熱電偶實(shí)測值異常高。此外，46個(gè)溫度點(diǎn)的誤差處于5%～10%之間，其余160個(gè)測點(diǎn)的誤差均處于5%以內(nèi)。

4 結(jié)論

本研究在三種不同的送風(fēng)速度下（轉(zhuǎn)速），采用CFD方法研究了帶有貫流風(fēng)機(jī)空調(diào)器房間內(nèi)部的氣流組織和溫度分布情況。研究表明，采用CFD二維仿真技術(shù)計(jì)算空調(diào)內(nèi)部及空調(diào)房間內(nèi)全尺度流場時(shí)，觀察到了貫流風(fēng)機(jī)靠近蝸舌的葉片處出現(xiàn)明顯的偏心渦，并且經(jīng)后續(xù)三維仿真驗(yàn)證，風(fēng)量誤差控制在5%以內(nèi)，說明使用本數(shù)值模擬方法具有較好的仿真精度。

其次發(fā)現(xiàn)，無導(dǎo)流葉片的二維全尺度風(fēng)機(jī)自循環(huán)流場仿真極易發(fā)生回風(fēng)短路現(xiàn)象，而適當(dāng)增加導(dǎo)流葉片數(shù)量可有效地將渦旋對(duì)稱線“下壓”，從而使得空調(diào)機(jī)正下方的渦旋負(fù)壓區(qū)較大，通過“拖曳”作用抑制回風(fēng)短路現(xiàn)象產(chǎn)生，但最終不能避免此現(xiàn)象產(chǎn)生。這可能與二維風(fēng)機(jī)自循環(huán)仿真模式下氣流分布缺少一個(gè)維度的流動(dòng)補(bǔ)償、容易上下波動(dòng)、回吸現(xiàn)象被明顯放大有關(guān)。

此外，在三維仿真時(shí)，室內(nèi)流場出現(xiàn)兩個(gè)對(duì)稱的大型渦旋，且渦旋的對(duì)稱性受到送風(fēng)速度的直接影響，送風(fēng)速度的大小也影響到了室內(nèi)的溫度均勻性，在高冷負(fù)荷時(shí)需要大送風(fēng)速度來改變渦旋對(duì)稱性，使得人員活動(dòng)區(qū)域具有較均勻的溫度，從而保證熱舒適性。

最后，通過對(duì)比1200 r/min工況下試驗(yàn)測試數(shù)據(jù)和模擬結(jié)果，發(fā)現(xiàn)210個(gè)溫度點(diǎn)中有160個(gè)的誤差低于5%，誤差異常點(diǎn)的數(shù)量僅有4個(gè)，約占總數(shù)的2%，因此本文的模擬方法具有較高的精度，其結(jié)果可為空調(diào)產(chǎn)品的開發(fā)和優(yōu)化提供參考，具有一定的應(yīng)用意義，并且高精度的氣流組織模擬方法有利于空調(diào)熱舒適性機(jī)理研究工作的深入開展。