王國強,邢超,趙玉軍
(陜西重型汽車有限公司,陜西 西安 710200)
在STAR-CCM+中風(fēng)扇的三種不同分析方法
王國強,邢超,趙玉軍
(陜西重型汽車有限公司,陜西 西安 710200)
文章通過三種不同的方法對風(fēng)扇進行了模擬,并用流量和壓降對分析結(jié)果進行評價。分析結(jié)果表明:三種方法均能較準(zhǔn)確的模擬風(fēng)扇,但在流場和壓力分布上有一定的差別。
風(fēng)扇;流場;壓降
CLC NO.:U467.3Document Code:AArticle ID:1671-7988 (2017)06-81-03
風(fēng)扇模型在計算分析中的運用,可使分析結(jié)果更貼合實際,但代價就是模型的前處理、計算模型的復(fù)雜性增加、計算資源和時間的增加。fan interface的提出,可以避免實際風(fēng)扇模型給計算帶來的諸多不便,可以使計算的收斂更快,但fan interface也有其自身的不足。fan interface 是對fan momentum source的一種簡單替代,它通過計算流過軸流風(fēng)扇的流量或速度模擬風(fēng)扇的壓降,然而卻不能像 fan momentum source那樣增加旋轉(zhuǎn)。究竟這兩者在計算結(jié)果上有啥差異,本文將通過對比分析進行探討。
1.1 “fan momentum source”介紹
為了更好地模擬軸流風(fēng)扇的性能,在STAR-CCM+中,使用了源項“fan momentum source”,“fan momentum source”使用了“actuator disk”理論來理論近似風(fēng)扇的工作過程,在做計算分析的過程中,可以不需要風(fēng)扇的實際模型。其理論分析過程如下:
為了得到風(fēng)扇作用到流動介質(zhì)上的力,本文先從風(fēng)扇前后介質(zhì)流動情況的速度三角形來陳述,如下圖所示:
圖1 風(fēng)扇葉片進出口速度三角形
符號列表:
Vin:介質(zhì)進口速度
Vaxial:介質(zhì)軸向速度
ε:葉片安裝角
β:葉片弦線與軸向的夾角
ω:風(fēng)扇旋轉(zhuǎn)角速度
r:葉片半徑
Vout,absolute:介質(zhì)出口絕對速度
Vout,relativ:介質(zhì)出口相對速度
介質(zhì)以Vin進入風(fēng)扇,并假設(shè)Vin= Vaxial,風(fēng)扇在θ方向上以ω的角速度進行旋轉(zhuǎn),且風(fēng)扇葉片本身不扭轉(zhuǎn)。
由連續(xù)性方程可得Vaxial.in=Vaxial.out,由動量守恒方程可知:作用于物體上所有力的總和等于通過物體表面的凈動量。基于此,在θ方向上有:
這是葉片作用于介質(zhì)的切向力,而總的力的方向垂直于葉片表面,故葉片作用于介質(zhì)的軸向力為:
如果將風(fēng)扇置于x-y平面中,則Fz=Faxial。為了將切向力在直角坐標(biāo)系下進行表述,現(xiàn)假設(shè)有一點P在直角坐標(biāo)下下的坐標(biāo)為(xp,yp,zp),c點風(fēng)扇的中心點,其坐標(biāo)為(xc,yc,zc),如下圖所示:
圖2 直角坐標(biāo)系下風(fēng)扇位置
為風(fēng)扇的壓升,Ψ為體積流量,C1和 C2是由載荷曲線上兩個工作點確定的常數(shù)。
1.2 風(fēng)扇定律
2.1 模型通用部分介紹
模型全長12米,進出口面積相等,直徑均為1.3米;進出口段長度相等,且直管段部分的長度為5.1米;風(fēng)扇段長度為0.125米。
圖3 風(fēng)扇的進氣段與出氣段
2.2 風(fēng)扇段不同建模方法
2.2.1 具有實際風(fēng)扇模型
風(fēng)扇順時針旋轉(zhuǎn),葉片數(shù)為9片,直徑為704.012mm,覆環(huán)直徑635mm。見圖4。
圖5 風(fēng)扇的計算域
圖6 Fan Interface面
圖4 實際風(fēng)扇模型
2.2.2 無實際風(fēng)扇模型,并啟用Moment Source Option項。
在此方法中,并無風(fēng)扇實際模型,僅保留風(fēng)扇的計算域,風(fēng)扇的特性由特性曲線給出。通過Tables和Reference Frames將風(fēng)扇特性加入到計算域中,并將風(fēng)扇域中 Momentum Source Option設(shè)置為Fan類型。見圖5。
2.2.3 無實際風(fēng)扇模型,無風(fēng)扇計算域,不啟動Moment Source Option項,僅用interface模擬風(fēng)扇性能。
在此方法中,僅用interface替代風(fēng)扇。將interface的類型設(shè)置為Fan,并通過Fan Curve將風(fēng)扇特性應(yīng)用于interface中。見圖6。
2.3 風(fēng)扇特性
當(dāng)風(fēng)扇轉(zhuǎn)速為3000rpm時,風(fēng)扇特性如下表所示:
表1 風(fēng)扇流量與壓升對應(yīng)關(guān)系
計算分析是在風(fēng)扇轉(zhuǎn)速一定(3000rpm)、進口流量一定的情況下進行分析計算的。并用風(fēng)扇壓升和風(fēng)扇后流線對不同方法進行評價。
進口邊界類型:Mass Flow Inle 進口流量:8.4kg/s 即6.99m3/s
出口邊界類型:Pressure Outlet 出口壓力:0Pa (參考壓力:101325Pa)
3.1 壓力流線圖
圖7 實際風(fēng)扇模型壓力流線圖
圖8 風(fēng)扇域模型壓力流線圖
圖9 Fan interface壓力流線圖
3.2 三種方式下的壓升對比
表2 三種方式下風(fēng)扇的壓升對比
從流線圖和表2上看:
(1)由于實際風(fēng)扇模型有9片葉片,在每個葉片的的后部都有明顯的尾流;
(2)相對于實際風(fēng)扇模型,風(fēng)扇域模型由于沒有風(fēng)扇實體模型,所以后部流線的分布和旋轉(zhuǎn)都比較均勻;
(3)在實際的分析過程中,F(xiàn)an interface模型模型簡單,計算效率高,但從流線來看,該模型并未考慮氣流的旋轉(zhuǎn),僅僅模擬了通過風(fēng)扇的壓升狀況;
(4)從表2來看,三種方法對壓升的分析均比較準(zhǔn)確,實際風(fēng)扇模型最好,風(fēng)扇域模型次之,F(xiàn)an interface模型的最大偏差也只有-6.7%,也基本能滿足分析要求。
風(fēng)扇模型在汽車熱管理分析或其它涉及到散熱、冷卻、通風(fēng)等方面都有廣泛的應(yīng)用,對不同的分析,考慮的著眼點不同,導(dǎo)致對分析結(jié)論的側(cè)重點不同,其中也涉及到模型的復(fù)雜程度和計算資源的配置情況,本文對風(fēng)扇模型的三種處理方法,希望能對計算結(jié)果的評價及計算分析過程中模型的選擇提供幫助。當(dāng)然,最好的模型還是實際風(fēng)扇模型。
[1] 李明,李明高.STAR-CCM+與流場計算[M].北京:機械工業(yè)出版社,2011.
[2] 林建忠,流體力學(xué)M].北京:清華大學(xué)出版社,2005.
Three Different Methods of Fan Simulation in STAR-CCM+
Wang Guoqiang, Xingchao, Zhao Yujun
( Shaanxi Heavy Duty Automobile Co., Ltd, Shaanxi Xi'an 710200 )
By means of three different methods simulation of fan is performed in STAR-CCM+. Its results were compared in according with flow rate and pressure drop. It could be found that the three methods can simulate the fan accurately, however,some differences in the distribution of flow field and pressure drop can be seen.
fan; flow field; pressure drop
U467.3
A
1671-7988 (2017)06-81-03
王國強,工程師,就職于陜西重型汽車有限公司汽車工程研究院,研究方向為汽車氣動與熱管理分析。
10.16638/j.cnki.1671-7988.2017.06.027