舒朝暉，張 強，段亞雄

貫流風機葉輪參數(shù)的模擬研究

舒朝暉，張 強，段亞雄

（華中科技大學，湖北武漢 430074）

貫流風機葉輪參數(shù)對風機性能有著決定性的影響，基于某企業(yè)提供的原型葉輪，通過二維穩(wěn)態(tài)模擬探討了葉輪前后緣半徑Ri與Ro、葉片傾角a以及葉片彎度角γ 等3個重要葉片參數(shù)對空調用貫流風機流場和性能的影響。結果表明，葉片前后緣半徑比對葉片壁面邊界層分離存在較大影響，當Ri/Ro＞1時，通常能夠較好地抑制邊界層分離，得到較大的出口流量；對于原型葉片而言，其傾角a的最佳取值范圍介于25°到30°之間，且隨a增大，葉輪入口流場波動減?。粡澏冉铅迷酱?，葉輪對氣流做功越多，γ 為90°時，取得最大流量值716.47m3/h，相比原葉輪，流量提升4.04%。

貫流風機；葉輪參數(shù)；流場特性；風機性能

1 前言

貫流風機（crossflowfan）作為空氣輸送裝置被廣泛應用于空調領域，尤其是應用于分體壁掛式空調中，其送風性能和噪聲特性備受關注。由于其流場極其復雜，目前仍然沒有通用的設計理論。

葉輪是貫流風機最關鍵的送風部件，其結構很大程度上決定了貫流風機的性能。前人對貫流風機葉輪結構做了一系列研究［1～6］。2003年，ALazzaretto總結前人的研究成果，提出并定義了葉輪和外殼的詳細形狀和位置參數(shù)，確定了能夠獲得較好性能和效率的各參數(shù)的取值范圍［6］。ALazzaretto提出的葉輪特征參數(shù)包括葉輪軸向長度與外徑之比L/D2、葉輪內外徑之比D1/D2、葉片數(shù)Z、葉片入口角b1和出口角b2，主要特征參數(shù)的取值范圍為 0.7 ≤ D1/D2≤ 0.85；20°≤β2≤45°；70°≤β1≤130°。其后，張師帥等以葉片出入口角、葉片中部厚度、葉片端部厚度、葉片斜扭角度以及葉片圓周分布角等為主要設計參數(shù)，實現(xiàn)了貫流風機葉輪幾何建模的參數(shù)化［7］。另外，很多學者研究了葉片扭曲角、葉片斜度、葉片分布以及葉片邊緣形狀（如鋸齒邊）等葉輪參數(shù)對貫流風機性能的影響［8～17］。

前人對葉輪橫截面參數(shù)的研究多集中在葉輪入口角β1和出口角β2上，鮮有學者研究葉片兩端厚度對性能的影響。另外，由于氣流兩次貫穿葉片，氣流相對葉片并沒有固定的入口角和出口角，所以，從葉片幾何設計的角度出發(fā)，本文參考Mcnally繪制雙圓弧葉片串列葉柵的形位參數(shù)，設定了一組葉片特征參數(shù)：外徑D2、弦長C、前后緣半徑Ri與Ro、最大厚度TM、彎度角γ和傾角a。結合某風機廠家提供的空調用貫流風機原型葉輪，本文維持參數(shù)D2和C不變，通過二維數(shù)值模擬的方法，分別探討葉片前后緣半徑Ri與Ro、彎度角γ和傾角a的改變對貫流風機性能的影響，為生產實踐提供參考［18］。

2 貫流風機特征及其物理模型

本文研究的貫流風機位于壁掛式空調室內機的下游，主要由蝸殼、蝸舌以及前彎葉片葉輪構成（如圖1）。與軸流風機和離心風機不同的是，貫流風機氣流由徑向進入，兩次貫穿葉片后從出口流出，其流場最典型的特點是，葉輪尾緣會產生大量的脫落渦，這些渦不斷破碎、擴散和組合，最終在蝸舌附近（圖1中虛線方框）形成一個較穩(wěn)定的偏心渦，其位置、大小和強度決定著貫流風機的流場特性。

圖1 空調室內機結構示意

葉片橫截面型線和各參數(shù)如圖2所示，圖（a）為由葉片角β1、β2和內外徑定義的等厚葉片（T為厚度），其中角β1與β2的頂點分別為葉片中弧線與葉輪內外圓的交點。圖（b）為由本文設定的特征參數(shù)確定的葉片。

圖2 葉片橫截面示意

本文的物理模型由企業(yè)提供的原型機簡化而來，如圖1所示，葉片型線與圖2（b）相同。蝸舌間隙4.6mm，蝸殼間隙3.7mm，葉片不等距分布，葉輪主要特征參數(shù)列于表1中。

表1 企業(yè)提供的貫流風機葉輪主要特征參數(shù)

葉輪轉速為1205r/m in，葉輪外緣速度為5.94m/s，馬赫數(shù)Ma＜0.02，以葉片弦長為特征長度的雷諾數(shù)Re=4838，可視為不可壓湍流模型。且葉片都為直葉片，貫流風機流場具有良好的軸向一致性。綜上，實際模型可簡化為二維不可壓湍流模型。

3 數(shù)值模型

3.1 網格劃分

計算網格如圖3所示，計算域為半徑為10D2的半圓區(qū)域，外圍圓環(huán)區(qū)域（圖中標記A，B）采用結構網格，其它區(qū)域采用非結構網格（葉輪區(qū)域D為三角形網格），并對葉輪、蝸舌和蝸殼附近進行加密處理。在進行網格無關性驗證后，最終選取的全局網格數(shù)約為19.5萬（改變葉片參數(shù)后葉輪區(qū)域網格會微小變化），其中葉輪區(qū)域D約10萬，外圍結構網格數(shù)為15520，換熱器區(qū)域C為7920，葉輪壁面邊界層為15層。轉子圓環(huán)寬度d以及圓環(huán)內外圓與鄰近葉尖的間歇對流場模擬結果影響較大，未見對其大小設置作嚴格規(guī)定的相關理論。本文中，如無特殊說明，則區(qū)域D的內外徑保持不變，分別為71.582mm和96.152mm，即分別與企業(yè)提供的原型葉輪內外葉尖間隔1mm。

圖3 網格示意

3.2 模擬條件設置

邊界條件設置：室內機外圍進出口（圖3（a））分別設置為壓力進出口，表壓均為0Pa；圖3中交界邊1，2，5和6設置為Interface邊界；交界邊3，4，7，8以及蝸舌加密區(qū)的兩端邊線設置為Interior邊界；其它邊界均設置為壁面Wall；換熱器區(qū)域C設為多孔介質區(qū)域；葉輪區(qū)域D采用多坐標參考系MRF，即設置為Fram eMotion；葉片壁面設置為旋轉壁面。

其它條件設置：空氣密度根據(jù)通風機標準進氣狀態(tài)設定為1.2kg/m3，選取RNGk-e湍流模型，壁面采用Scalab leWall模型，壓力速度耦合采用SIMPLEC格式。壓力采用標準差分格式，其它空間離散設置為二階迎風格式。

4 模擬結果及分析

4.1 前、后緣半徑的影響

原葉輪前、后緣半徑分別為Ri=0.218mm、Ro=0.268mm，為了探究前、后緣半徑的影響，保持其它葉片橫截面特征參數(shù)不變，分別改變Ri，Ro（D1會隨其微小變化，具體值將列出）進行二維穩(wěn)態(tài)計算，不同條件的性能參數(shù)變化情況見圖4。

圖4 不同前、后緣半徑情況的性能參數(shù)變化

葉輪功率：

式中P——葉輪功率，W

T——單位長度的葉輪轉矩，N·m

ω——葉輪旋轉角速度，rad/s

全壓效率：

式中Pe——風機有效功率，Pe=ΔPeQ

ΔPe——風機全壓，Pa

Q——出口風量（注意軸向長度），m3/h

各序號對應的前后緣半徑值如表2所示（為方便表述，令常數(shù)a=0.218mm，a～1.23a表示0.218～0.268mm），其中序號4對應原葉片，同時 附上對應的葉片角β1，β2以供參照。

表2 不同前后緣半徑的流量變化和對應的葉片角

由圖4和表2可得，固定Ri=a時，當Ro從2.1a（此時后緣直徑已十分接近于葉片最大厚度TM）逐漸減小到a時（序號1到5），風機出口流量、葉輪功率和全壓效率均小幅增大，說明這一過程中葉輪對氣流做功增多且流動損失相對減??；固定Ro=a時，當Ri從a逐漸增大到2.1a時（序號5到9），風機流量逐漸減小，但變化程度十分微小，葉輪功率基本不變，全壓效率呈微小的減小趨勢，說明流量的微小變化主要由流動損失的增大所致。

總體而言，隨葉片前后緣半徑變化，3個性能參數(shù)均先微小增大（序號1到5）后基本不變（序號5到9），都在序號5處取得最大值。流量最大值相對于流量極低值（序號1）增大2.42%，相對于原風機（序號4）增大0.58%。

圖5，6依次給出了a～1.8a，a～a和2.1a～a 3種葉片入口和出口局部渦量云圖。

圖5 葉片入口局部渦量云圖

圖6 葉片出口局部渦量云圖

由圖可知，隨葉片前緣半徑增大（Ri/Ro增大），入口葉片吸力面上的渦量更加貼近壁面，吸力面附近尾渦區(qū)稍有減?。▍⒁妶D中橢圓標記），這是因為，根據(jù)翼型理論，葉片入口端大出口端?。▓A頭尖尾）能夠有效減小壁面逆壓梯度，從而抑制邊界層的分離或使分離點后移，使得尾渦減小，流動阻力與損失減??；而出口葉片上的渦量分布則剛好與入口情況相反，只是沒有入口明顯，仔細觀察出口葉片尾緣（圖中圓形標記），可以發(fā)現(xiàn)，隨前緣半徑增大，葉片尾緣附近的渦量范圍略有增大且尾緣附近渦量強度很大。即隨葉片前緣半徑增大（或后緣半徑減?。?，貫流風機葉輪進口流道的流阻減小而出口流道的流阻增大。因此，理論上存在一組最佳的前后緣半徑值，使得整個葉輪流道的流阻最小，從而使得風機流量最大。由以上9組模擬結果可知，當Ri≥Ro時，一般都能獲得較大的出口風量，且此時全壓效率和葉輪功率基本不變。

綜上可知，葉片前、后緣半徑Ri，Ro對貫流風機性能的影響較小。一般而言，可使用Ri略小于Ro的葉片來適當降低葉輪功率，而為獲得較大的流量值和全壓效率最好使得Ri≥Ro。

4.2 葉片傾角α的影響

由于葉片傾角改變時，葉輪內徑發(fā)生較大變化，為保證轉子區(qū)域D的網格質量，本節(jié)把區(qū)域D（圖3）的內徑由71.582mm改為70mm。不同葉片傾角時貫流風機性能參數(shù)變化情況如圖7所示，虛線對應原葉片。不同情況的具體流量值和葉片角示于表3中。

圖7 不同傾角α時的性能曲線

表3 不同葉片傾角的流量變化和對應的葉片角

需要特別說明的是，當傾角為15°、20°和23°時，風機出口流量模擬值存在比較明顯的波動，流量值無法收斂，推測這3種情況下風機內流場十分紊亂。為便于對照，本文選取從8000步到10000步的模擬結果的流量平均值作為這3種情況對應的參考流量。下文將對傾角為20°的流場予以分析。另外，對于葉片彎度角為45°和60°的情況，本文也采取了相同的處理方式，后文將不再說明。

由圖7和表3可知，傾角α由25°增大到30°的過程中，風機各性能參數(shù)的變化都十分微?。◣缀醪蛔儯．敠劣?0°變?yōu)?5°時，各性能參數(shù)均明顯減小。

圖8給出了不同傾角時貫流風機內部的靜壓云圖和流線分布。

圖8 不同傾角靜壓云圖和流線分布

結合圖7和圖8可知，隨傾角α增大，偏心渦位置越靠近蝸舌，這是因為傾角越大，渦核下游附近氣流以更大的負攻角通過葉片流道，下游葉片壓力面產生的脫落渦更多，對偏心渦的誘導作用更大。而一般渦心越靠近蝸舌，出口回流區(qū)域越小，風機流量趨于增大；偏心渦位置向蝸舌移動使得入口氣流方向發(fā)生明顯改變，氣流在入口葉片吸力面產生的脫落渦減少，因而入口流場紊亂程度明顯減小，流動損失減少，風量和效率趨于增大。另外，渦核壓力隨傾角增大而變大，這導致渦核對流場的誘導作用減小，從而使風機整體流量趨于減小?？偟膩碚f，可以推測，傾角為15°、20°和23°時，葉輪入口流場的紊亂性對整個風機流場產生了強烈擾動，使得出口流量無法收斂；傾角由25°增大到35°的過程，一方面由于偏心渦靠近蝸舌使得風機流量趨于增大，另一方面由于渦心壓力變大導致流量趨于減小，當前者趨勢大于后者時，風機流量增大，反之就減小。同時，由于傾角越大，風機葉輪入口流場波動越小，所以在流量變化不大時，可以選擇較大的葉片傾角。

4.3 葉片彎度角γ

圖9給出了不同葉片彎度角時貫流風機各性能參數(shù)的變化情況，虛線對應原葉片，不同情況的具體流量值和葉片角示于表4中。

圖9 不同彎度角γ時的性能曲線

由圖9和表4知，出口流量隨彎度角γ 增大先增大后減小，當γ＞75.49°后出口流量增大趨勢明顯減小，當γ=90°時取得流量最大值716.47m3/h，相對原葉輪（γ=75.49°），風機流量提升4.04%；葉輪功率隨γ 增大近于線性增大，當γ＞90°后增大幅度減小，γ=90°時的葉輪功率相比原葉輪功率增大14.87%；全壓效率隨γ 增大先逐漸增大，當γ＞75.49°后基本不變，當γ＞90°后，開始減??；γ=90°時的全壓效率與原葉輪基本相同。這說明，隨彎度角增大流量的增大主要得益于葉輪對氣流做功增多；當γ＞75.49°后流量的增量遠小于葉輪功率的增量且全壓效率基本不變，說明這一過程流動損失并沒有明顯變化，葉輪對氣流做功主要使得出口氣流的全壓大幅提升而對氣體輸送量的提升并不十分明顯。

圖10不同彎度角速度云圖和流線分布

圖10 給出了彎度角為45°、90°、100°3種情況的流線圖和速度云圖。由圖可知，葉輪內的低速區(qū)主要出現(xiàn)在偏心渦區(qū)域和葉輪入口部分湍動嚴重的葉片尖端。隨γ增大，葉輪入口低速區(qū)域和旋渦區(qū)域明顯減小，結合流線圖可以看出這主要得益于葉輪入口尾渦減?。徊⑶胰~輪出口高速區(qū)域明顯增大，說明彎度角的增大提升了葉輪對氣流的做功能力，這與圖9中性能曲線的變化情況相符合。另外，結合流線圖和速度云圖可以明顯看出，彎度角的改變使偏心渦渦心位置發(fā)生了變化，一方面，隨彎度角增大，渦心位置逐漸向右移動，越來越靠近蝸舌，出口回流區(qū)域明顯減??；另一方面，彎度角越大，渦心位置越靠近葉輪軸心，偏心渦上游循環(huán)區(qū)域明顯增大。因此，存在一個最佳的葉片彎度角γ，使循環(huán)區(qū)域盡可能小，從而使得風機出口流量盡可能大。由模擬結果可知，在不考慮風機效率的情況下，這一最佳彎度角應在90°附近取得。渦心位置的改變以及入口尾渦減小的原因與“4.2葉片傾角α的影響”中的分析相同。

5 結論

（1）葉片前后緣半徑Ri與Ro、彎度角γ和傾角α等葉片參數(shù)對貫流風機流場有著重要影響，它們與葉片表面邊界層的分離息息相關，一定程度上決定了偏心渦的位置、大小和強度。

（2）葉片前后緣半徑比Ri/Ro主要影響著葉片壁面的邊界層分離以及葉輪尾渦區(qū)域的大小，一般而言，Ri/Ro值越大，葉輪總體尾渦區(qū)域越小，流動損失越小。

（3）葉片傾角α 越大，偏心渦越靠近蝸舌，出口回流區(qū)域越小，但渦心壓力越來越大，不利于渦核對氣流的吸引誘導。對于原型葉輪，當α≤23°時，葉輪入口渦脫落摻混十分嚴重，貫流風機流場無法穩(wěn)定，出口流量不能收斂；當α＞30°時，風機流量和全壓效率均明顯下降；對企業(yè)提供的貫流風機而言，最佳葉片傾角范圍為25°≤α≤30°。

（4）葉片彎度角γ越大，葉輪對氣流做功越多。當γ＞75.49°后，葉輪對氣流做功主要使得出口氣流的全壓大幅提升而對氣體輸送量的提升并不明顯。就本文而言，若考慮風機的綜合性能，最佳彎度角應在75.49°附近取得；若只是為了獲得較大的出口流量，則彎度角γ可取為90°，此時取得最大流量值 716.47m3/h。

［1］埃克B，著.沈陽鼓風機研究所，譯.通風機（第一版）［M］.北京：機械工業(yè)出版社，1983.

［2］Yam afu jiK.Studiesontheflowo fcross-flowim pe-llers（1streport，analyticalstudy）［J］.Bu lletino ftheJSME，1975，18（123）：1018-1025.

［3］Yam afu jiK.Studiesontheflowo fcross-flowim pe-llers（2ndreport，analyticalstudy）［J］.Bu lletino ftheJSME，1975，18（126）：1425-1431.

［4］ Lazzaretto A，Toffolo A，Martegan，A D. A systematic experimental approach to cross-flow fan design［J］. Journal of Fluids Engineering，2003，125：684-693.

［5］ Lazzarotto L，Lazzaretto A，Martegani A D，et al. On cross-flow fan similarity：effects of casing shape ［J］. Journal of Fluids Engineering，2001，123：523-531.

［6］ Lazzaretto A. A criterion to define cross-flow fan design parameters［J］. Journal of Fluids Engineering，2003，125（4）：680-683.

［7］張師帥，羅亮.空調用貫流風機葉輪幾何建模的參數(shù)化［J］.風機技術，2006（5）：14-16.

［8］劉敏，王嘉冰，吳克啟.數(shù)值模擬不等距葉片對貫流風機的影響［J］.工程熱物理學報，2007，28（2）：211-214.

［9］田杰，歐陽華，李游,等.基于交錯葉輪技術的橫流風機氣動聲學特性研究［J］.機械工程學報，2010，46（3）：97-102.

［10］朱宗銘，鄔烈鋒，譚華兵，等.貫流風機性能參數(shù)實驗與扭曲角優(yōu)化［J］.機械設計與制造，2012（1）：203-205.

［11］王銀姣，盧劍偉，江斌，等.利用CFD技術研究葉片斜度對貫流風機性能的影響［J］.合肥工業(yè)大學學報：自然科學版，2012，35（7）：882-887.

［12］ Yamakawa H. Development of high efficiency and low noise cross-flow fans for room air conditioner indoor unit［C］.ASME 2014 4th Joint US-European Fluids Engineering Division Summer Meeting collocated with the ASME 2014 12th International Conference on Nanochannels，Microchannels，and Minichannels. American Society of Mechanical Engineers，2014：V01BT10A024.

［13］陳陽，黃友根，王軍，等.傾斜蝸舌對前彎離心風機蝸舌區(qū)域流動影響的數(shù)值研究［J］.流體機械，2015，43（9）：26-30.

［14］丁學亮，葉學民，李春曦.軸流風機葉片展向結構變化對性能影響的數(shù)值分析［J］.流體機械，2015,43（9）：43-49.

［15］賈志彬，陳強，汪軍，等.微穿孔消聲器在小型高速離心風機中的應用研究［J］.流體機械，2015,43（8）：6-10.

［16］池優(yōu)陽，阮競蘭.糧油機械產品的低碳設計制造［J］.包裝與食品機械，2015,33（5）：55-58.

［17］鞏桂芬，陳寧.螺旋橫向運動式去皮裹包果蔬產品的機械及力學研究［J］.包裝與食品機械，2016,34（5）：33-37.

［18］ Crouse J E，Nally W D M. Fortran program for computing coordinates of circular arc single and tan -dem turbomachinery blade sections on a plane［Z］. 1970.

Numerical Investigations on the Impeller Parameters of the Cross Flow Fan

SHU Zhao-hui，ZHANG Qiang，DUAN Ya-xiong
（Huazhong University of Science and Technology，Wuhan 430074，China）

Impeller parameters have great effect on the performance of the cross flow fan. Based on the blade provided by a company，the influences of several important blade parameters such as leading radiu（Ri），rear radiu（Ro），blade tilt angle（a） and camber angle（γ） on the flow field and performance of the cross flow fan，has been studied by two-dimensional steady simulation.The results show that，the ratio of Riand Ro has great effect on boundary-layer separation of the blade surface，the boundary layer separation is inhibited and higher outlet flow rate is obtained when Ri/ Ro＞1. Blade tilt angle a of the prototype blade can be chosen in the ranges 25°≤a≤30°，and in which range，the fluctuation of the impeller inlet flow decreases with a increases. The greater the camber angle γ，the more the impeller works on the airflow. Moreover，when γ is equal to 90°，the flow rate achieves the maximum of 716.47 m3/h，enhanced 4.04% compared to the original impeller.

cross flow fan；impeller parameters；flow field characteristics；fan performance

TH43

A

10.3969/j.issn.1005-0329.2017.11.005

1005-0329（2017）11-0022-07

2016-11-24

2017-04-19

廣東省科技計劃項目（2013B090600061）

舒朝暉（1970-），男，博士，副教授，主要從事制冷空調及其節(jié)能技術、高效熱交換設備與技術等方面的研究,通訊地址：430074湖北武漢市華中科技大學能源與動力工程學院化工裝備與技術研究所，E-m ail：shuzhaohui@hust.edu.cn。