造雪機風扇位置布局研究

白鵬程，金向陽，2，王莎莎，薛加磊，3，林晶，2，董子昂

（1.哈爾濱商業大學 輕工學院，哈爾濱 150028；2.虛擬制造技術福建省高校重點實驗室，福建 泉州 362000；3.江蘇大學 流體機械工程技術研究中心，江蘇 鎮江 212000）

0 引言

伴隨著科技的發展，人造雪已成為滑雪場、雪雕制作的首選，較自然降雪由于有產量大、不限場地降雪等優點，已被各大雪場廣泛采用。隨著對造雪機的要求逐漸增高，高效率造雪機的開發成為被研究的課題。影響造雪機降雪品質的因素主要有兩個方面：1）造雪機核子器，核子器主要影響成雪率，是造雪機的核心之一；2）造雪機風扇氣動系統，氣動系統主要影響成雪品質和機器的能耗，是造雪機另一核心。因此研究造雪機的氣動系統對造雪機的優化有指導意義。

目前國內外一些學者對于風扇氣動系統的研究已有了一定的進展和成果，吳曉林等[1]研究了發動機冷卻風扇的布局問題，以風量和噪聲程度為評價標準，獲得了最優的風扇布局；唐釗[2]在多因素對風扇氣動性能影響的研究中引入正交試驗的方法，將風扇葉型的參數作為試驗因子，仿真分析了不同參數的風扇性能，在提高風扇氣動性能的指標下，得到參數的最優組合；Kohri等[3]在探究不同的湍流模型及網格尺寸對風扇流場仿真的影響時，采用了計算流體力學（CFD）仿真方法，對比了仿真結果，最后得到優化的模型和網格精度；Peter Gullberg和Lennart Lodahl[4]為研究風扇在不同轉速下的氣動變化規律，引入旋轉參考系法（MRF）對風扇模型進行仿真分析，并總結了基于MRF模型的簡單“速度校正”，探究了這種研究方法在仿真計算中的局限性和適用性。

造雪機氣動系統中，除了風扇布局的問題，風扇的性能和結構參數、風腔的結構參數、電動機的性能參數等都對氣動系統有較大的影響，本文在以造雪機原參數結構的基礎上運用CFD計算方法，在仿真軟件中模擬風扇的流場特性和氣動性能，并將仿真結果與原參數模型對比，結果顯示，最優布局下風扇的氣動性能明顯提升。

1 CFD理論基礎

計算流體力學運用離散數值模擬方法，模擬流體的無黏繞流和黏性流動特性，可以對流體運動的微分方程數值求解，得到流體的流場分布。

1）質量守恒方程。單位時間內流體中微元質量的增加量，等于同單位時間間隔流入該微元體的質量。微分形式的不可壓流體質量守恒方程可表述為式中，u，v，w分別為（x，y，z）坐標方向上的速度矢量。

2）動量守恒方程。單位時間微元中流體的動量變化率等于該流體微元上受到外界作用力的和[5]。基于此定律可以得到（x，y，z）方向的動量守恒方程，它規定了氣流里壓強的變化、速度的變化及與徹體力之間的關系。此方程稱為納維-斯托斯克（Navier-Stokes，N-S）方程：

式中：T為熱力學溫度，K；cp為定壓比熱容；k為流體的導熱系數；ST為黏性耗散項，是流體黏性作用轉化為熱能的部分[6]，對于不可壓縮流體，若熱交換很小，則可不引入此方程[7]。

2 仿真建模與計算

在三維軟件中建立造雪機風扇的結構模型，并導入計算軟件中劃分網格，設置計算條件和仿真參數，得到仿真計算結果。通過結果分析造雪機風扇的流場特性。

2.1 仿真模型的建立

在三維軟件中建立造雪機風扇的參數模型，風扇的葉片數為13片，葉片安裝角度為29°，葉片寬度為80 mm，輪轂比為0.43，距離入風口的距離為100 mm，如圖1所示。

圖1 風扇與風腔模型

將建立好的模型導入仿真計算軟件中，為研究造雪機風腔內的流場特性，建立圓周直徑為800 mm的圓柱，代表風扇內部的流場，主要模擬葉片周圍的流場特性，稱為內流域，在風腔的基礎上建立大的外部流場，模擬風扇在風腔內的流場特性，稱為外流域，內外流域一起模擬可以更好地觀察風扇的氣動性能。建立的風扇流域如圖2所示。

圖2 風扇的內外流域

將內外流域進行網格劃分，因風腔的區域面積較大，風腔內的整體流域采用非結構化網格劃分，網格大小為20 mm，在此基礎上對風扇的內流域進行網格劃分，網格大小為10 mm，為使計算精度滿足要求，風扇的整體結構進行局部細化[8]，輪轂及扇葉的尺寸大小為5 mm，內部旋轉域與外流域的交界面需要保持均勻的過渡[9]，網格大小設置為15 mm，網格劃分效果如圖3所示。

圖3 網格劃分結果

2.2 邊界求解設置

對于不同特點的旋轉機械仿真方法，主要包括3種：1）旋轉參考系法（MRF）。這種方法適用于參考系旋轉的穩態計算，因為只涉及參考系旋轉，網格未動，所以特點是計算快，求解穩態計算問題時非常適用。2）滑移網格法（RBM）。這種方法適用于求解瞬態計算的精細過程，在計算過程中，網格旋轉，對細節問題可以加以捕捉。3）網格重疊法（OVERSET）。這種方法適用于對網格邊界精度要求高的情況，因為使用了內外兩層網格嵌套，所以計算精度更容易滿足要求。

為了更好地捕捉旋轉細節，仿真計算采用了滑移網格法。在Fluent中設置計算參數，內外流域相交的交界面建立3個interface對，風腔及固定部件設置為固壁條件[10]，湍流模型選用k-ε湍流模型，進口條件設為Pressureinlet，出口條件為Pressure-outlet，選用滑移網格運動，轉速設為2960 r/min，時間步長設為0.02 s，迭代次數為20，仿真計算后對結果進行分析。

3 結果分析與比較

3.1 結果分析

仿真結果完成后，在后處理軟件中處理得到風扇旋轉的速度云圖、速度流線圖和速度矢量圖。圖4可以直觀地表示風扇旋轉域速度的變化，葉尖處因距離旋轉中心距離最遠，因此線速度最大，越靠近旋轉中心，線速度逐漸減小，在輪轂處的速度最小，這也證明了仿真的結果符合實際。

圖4 內流域速度云圖

圖5直觀地反映了風扇的流場特性，在外流域部分，靠近風扇的出口處，有部分流線雜亂無章，呈現紊流狀態；在整體流域上，風扇旋轉產生的速度場大多為層流狀態；風腔入口處速度平穩，經過風扇的旋轉作用，風向呈現螺旋向前，在出口處有波狀的擺動現象，流速越快，擺動風振幅和頻率越大，為過渡流。圖6直觀地反映局部細節的流場特性，圖6（a）所示為在扇葉周圍的速度矢量，氣體流速在從輪轂到葉尖的方向逐漸增大，在葉尖和風扇的邊緣處達到極值，在風扇的輪轂和風腔的邊緣處，氣體流速較小，風扇產生的風圍繞扇葉邊緣處形成了一個“風圈”，在圖6（b）中顯示為流速最大的部分。

圖5 速度流線圖

圖6 速度矢量圖

3.2 結果比較

為研究風扇在風腔的布局最優問題，以原模型風扇的旋轉中心為原點，建立向風扇葉尖處的X軸與Y軸，以風扇入風口為負方向，出風口為正方向建立Z軸的坐標系，討論風扇位置對風腔出口的風速、總風量和靜壓效率的影響。

3.2.1 軸向位置的影響

風扇在風腔中的軸向距離，表示了風扇在風腔中位置，為了對比3種不同情況下風扇的流場特性，建立3種不同軸向位置的風扇三維模型，將模型進行仿真后，處理數據得到3種情況的速度比較，如圖7所示。

圖7 不同軸向位置速度對比

在X的正半軸，風扇的軸向坐標Z=-50 mm時，最大風速約為45 m/s，越靠近風扇的輪轂風速越低，當軸向坐標處于原點時，風扇的最大風速并未有明顯的變化，而當軸向位置移動到50 mm處時，風扇的最大風速略微提升，但整體提升不大，這說明這3種風扇的位置對于風扇風速的影響并不明顯，提取3個位置的風扇的風量和靜壓效率，如表1所示。

表1 不同軸向位置風扇風量與靜壓效率

將3種位置的風扇風量和靜壓效率作成折線圖，可以更直觀地對比出3種位置的差別，如圖8、圖9所示。

圖8清楚地反映了3種位置下風量的變化情況，軸向距離Z由-50 mm變為原點時，風量由52 719.86 m3/h降到了49 028.37 m3/h，風量下降了約7%，由原點位置變為50 mm時，風量繼續下降到47 650.19 m3/h，可以判斷出風扇軸向位置越靠近風腔的入口處，風量的提升越大。

圖8 不同軸向位置風量對比

靜壓效率是表征風扇性能好壞的重要參數，圖9中風扇的位置Z由-50 mm到原點再到50 mm 變化的過程中，風扇的靜壓效率也是逐漸在減小，風扇在風腔的軸向位置不僅影響風量，對風扇的靜壓效率也產生較大影響，越靠近風腔的入口處，風扇的靜壓效率越高。

圖9 不同軸向位置靜壓效率對比

3.2.2 徑向位置的影響

造雪機風扇在風腔中的徑向位置，表示風扇與風腔間隙的大小，取風扇與風腔的間隙為60 mm、100 mm和140 mm 建立造雪機氣動系統，將3種模型進行仿真，對比3種位置下風速的大小，如圖10所示。

圖10 不同徑向位置的風扇風速

間隙為60 mm時，風扇的最大風速明顯高于其余兩種間隙的風扇風速，最高風速約為46 m/s。間隙增大到100 mm和140 mm，最大風速并無差別，風扇與風腔的間隙越小，可以產生的風速越快。通過計算將3種位置的風扇風量和靜壓效率制成表，如表2所示。

表2 不同徑向位置風扇風量與靜壓效率

以同樣的方式建立不同徑向位置的風量和靜壓效率折線圖，進行對比分析，如圖11、圖12所示。

從圖11可以很直觀地看出間隙在由60 mm增大到100 mm時，風量約由48 843 m3/h增長到了49 028 m3/h，而間隙繼續增加到140 mm時，風量則下降到了48 708 m3/h，證明了風扇與風腔的間隙并不是越小越好，當間隙過小時，會減小風扇的風量，同時間隙過小也會給扇葉和風腔的制造精度帶來更高的要求，并增大風扇的噪聲[12]。

圖11 不同徑向位置風扇風量對比

從圖12可以觀察出徑向間隙的增大，帶來的是風扇靜壓效率的降低，風扇與風腔的間隙越大，風扇的靜壓效率越低，這是因為較大的間隙帶來的是更多的葉尖處的回流[13]，導致了風扇靜壓效率的降低，因此需要選擇合理的間隙距離。

圖12 不同徑向位置風扇靜壓效率對比

4 結語

為探究造雪機風扇在風腔中的位置對風扇氣動性能的影響，提升造雪機的工作效率，以造雪機氣動系統為研究對象，對造雪機中的風扇建立三維數模，運用計算流體力學理論原理，對不同位置的風扇模擬仿真，對比風速、風量和靜壓效率的影響。結論如下：1）運用CFD計算造雪機風扇在風腔中的工作，模擬風扇的流場特性，證明了數值仿真計算的可靠性；2）造雪機風扇在風腔中的軸向位置對風速影響不明顯，但越靠近風腔入口處，風量和靜壓效率越高；3）造雪機風扇在風腔中的徑向位置越小，風扇產生的風速越高，但風量也會略微降低，合理選擇風腔的間隙，能使風扇風量和靜壓效率達到最優值。