01背景介紹
螺桿式壓縮機(jī)具有尺寸小、重量輕、易維護(hù)等特點(diǎn),是制冷壓縮機(jī)中發(fā)展較快的一種機(jī)型。隨著工作可靠性的不斷提高,螺桿式壓縮機(jī)在中等制冷量范圍內(nèi)已逐漸替代往復(fù)式壓縮機(jī),并占據(jù)了離心式壓縮機(jī)的部分市場(chǎng)。
壓縮機(jī)內(nèi)部的溫度變化對(duì)其壓縮性能和效率有著至關(guān)重要的影響。固體受熱膨脹會(huì)導(dǎo)致金屬部件發(fā)生過(guò)度磨損,從而造成泄露間隙的改變,進(jìn)而對(duì)產(chǎn)品的性能產(chǎn)生影響,為了準(zhǔn)確預(yù)測(cè)壓縮機(jī)的熱力學(xué)特性,兼顧產(chǎn)品研發(fā)周期和經(jīng)濟(jì)性,通常采用CFD技術(shù)對(duì)壓縮機(jī)的流固共軛傳熱(CHT)問(wèn)題進(jìn)行研究。但是,由于固體的傳熱速度要比壓縮氣體的慢得多,如果利用CFD技術(shù)直接對(duì)固體結(jié)構(gòu)和氣體介質(zhì)進(jìn)行耦合傳熱模擬,可能需要計(jì)算足夠多的壓縮機(jī)旋轉(zhuǎn)數(shù)后才能得到一個(gè)穩(wěn)定的CHT解,那么模擬的運(yùn)行時(shí)間可能就變得不切實(shí)際。因此,尋求一種先進(jìn)的方法解決流固共軛傳熱問(wèn)題迫在眉睫。
02難點(diǎn)分析
目前,利用CFD技術(shù)進(jìn)行螺桿壓縮機(jī)流固共軛傳熱分析存在以下問(wèn)題:
間隙設(shè)計(jì)是螺桿壓縮機(jī)中的重要問(wèn)題,間隙必須足夠小以提高容積效率,同時(shí)又必須有足夠的間隙來(lái)防止轉(zhuǎn)子干涉。對(duì)傳統(tǒng)CFD軟件而言,如果預(yù)留的最小間隙與實(shí)際情況一致,會(huì)導(dǎo)致動(dòng)網(wǎng)格生成失敗,或使網(wǎng)格質(zhì)量下降、網(wǎng)格總數(shù)急劇增加,因此對(duì)于螺桿壓縮機(jī)的CFD模擬幾乎成為不可能的任務(wù);
由于流體側(cè)和固體側(cè)傳熱相互影響,進(jìn)行溫度場(chǎng)計(jì)算時(shí),無(wú)法確定流固交界面的邊界條件,因此需要考慮流固共軛傳熱,將流固之間難以確定的邊界條件,轉(zhuǎn)換成耦合計(jì)算的內(nèi)部邊界,使計(jì)算更符合實(shí)際工況;
模型需考慮流固共軛傳熱,流體為理想氣體;固體包括機(jī)殼、陽(yáng)轉(zhuǎn)子和陰轉(zhuǎn)子等部件,需要仿真軟件具有高效的前處理能力,可以快速實(shí)現(xiàn)建模與計(jì)算;
由于螺桿壓縮機(jī)運(yùn)動(dòng)的復(fù)雜性,需要構(gòu)建高質(zhì)量的網(wǎng)格并設(shè)置動(dòng)網(wǎng)格,這對(duì)于傳統(tǒng)CFD軟件是一個(gè)較大的挑戰(zhàn);
需要同時(shí)考慮穩(wěn)態(tài)、瞬態(tài)以及傳熱等,對(duì)于CFD軟件的求解能力要求較高;
流體域計(jì)算與固體域計(jì)算時(shí)間相差較大,如何更好地實(shí)現(xiàn)流固共軛傳熱仿真,對(duì)于傳統(tǒng)CFD仿真難度較大。
03Simerics-MP+解決方案
鑒于上述螺桿壓縮機(jī)流固共軛傳熱分析幾何模型、物理現(xiàn)象和運(yùn)動(dòng)的復(fù)雜性,對(duì)CFD分析軟件的選取提出了較高的要求。Simerics-MP+作為專(zhuān)業(yè)的運(yùn)動(dòng)機(jī)械CFD模擬專(zhuān)家,在螺桿壓縮機(jī)流固共軛傳熱領(lǐng)域具備獨(dú)特的優(yōu)勢(shì):
Simerics-MP+具備螺桿壓縮機(jī)轉(zhuǎn)子網(wǎng)格導(dǎo)入接口,可直接將螺桿壓縮機(jī)轉(zhuǎn)子部件的網(wǎng)格快速導(dǎo)入并完成動(dòng)網(wǎng)格設(shè)置;
Simerics-MP+引進(jìn)一種先進(jìn)的共軛傳熱求解的新方法—混合時(shí)間尺度耦合法(Mixed Timescale Coupling method),這種方法可以解決熱量傳播過(guò)程中的時(shí)間尺度問(wèn)題,快速獲取可靠結(jié)果;
Simerics-MP+內(nèi)置的應(yīng)變應(yīng)力求解器可進(jìn)行單向流固耦合求解,將轉(zhuǎn)子壁面溫度反饋到固體結(jié)構(gòu)進(jìn)行熱應(yīng)力和熱膨脹變形的求解;
Simerics-MP+具有高效的求解功能,其求解器基于傳統(tǒng)的CFD求解器進(jìn)行了優(yōu)化,相較于傳統(tǒng)的CFD工具,其求解速度更快。結(jié)合強(qiáng)大的網(wǎng)格技術(shù),Simerics-MP+可以順利進(jìn)行系統(tǒng)級(jí)的CFD分析工作。目前Simerics-MP+在系統(tǒng)級(jí)分析領(lǐng)域已積累較多經(jīng)驗(yàn)。
04雙螺桿壓縮機(jī)案例
本文以N35無(wú)油雙螺桿壓縮機(jī)為研究對(duì)象,陽(yáng)轉(zhuǎn)子的運(yùn)行速度為6000rpm到14000rpm,陽(yáng)轉(zhuǎn)子直徑為127.45mm,陰轉(zhuǎn)子直徑為120.02mm,兩個(gè)轉(zhuǎn)子之間的中心距離為93.00mm。轉(zhuǎn)子的長(zhǎng)徑比為1.6,陽(yáng)轉(zhuǎn)子的包角為285.0度。同時(shí)考慮流體域和固體域,以實(shí)現(xiàn)兩者之間的共軛傳熱問(wèn)題。
4.1 流體模型
采用專(zhuān)業(yè)的網(wǎng)格生成軟件SCORG進(jìn)行雙螺桿轉(zhuǎn)子部分流體域的網(wǎng)格劃分。在SCORG中根據(jù)不同的旋轉(zhuǎn)角度創(chuàng)建了一系列的轉(zhuǎn)子網(wǎng)格文件。通過(guò)Simerics-SCORG網(wǎng)格接口將轉(zhuǎn)子網(wǎng)格文件讀入求解器,而流體域的吸入端和排出端則使用Simerics二叉樹(shù)非結(jié)構(gòu)化網(wǎng)格生成器來(lái)劃分網(wǎng)格。所有的流體域均通過(guò)交互面MGI進(jìn)行連接,流體域網(wǎng)格總數(shù)約為145萬(wàn),流體模型幾何及網(wǎng)格如圖1所示
4.2 固體模型
固體模型包括三個(gè)域:機(jī)殼、陽(yáng)轉(zhuǎn)子和陰轉(zhuǎn)子,網(wǎng)格為基于二叉樹(shù)算法的笛卡爾網(wǎng)格,網(wǎng)格數(shù)共約40萬(wàn),固體模型幾何和網(wǎng)格如圖2所示
4.3 邊界條件設(shè)置
對(duì)于流體域仿真而言,流固交界面設(shè)置為固定溫度邊界,溫度值從固體模型模擬結(jié)果映射而來(lái);對(duì)于固體域仿真而言,流固交界面設(shè)置為固定的熱通量邊界,熱通量值由流體模型的模擬結(jié)果映射得到,仿真分析原理如圖3所示。
就邊界條件而言,機(jī)殼外表面設(shè)置為熱對(duì)流邊界,機(jī)殼外表面在300K環(huán)境溫度下的對(duì)流換熱系數(shù)假設(shè)為10W/m2*K。流體域和固體域的初始溫度設(shè)置為300K,其余邊界條件如表1所示。
本案例所用硬件配置為2.20GHz E5-2630 v4 雙CPU的工作站,采用20核并行計(jì)算,流體模型中陽(yáng)轉(zhuǎn)子每旋轉(zhuǎn)一周用時(shí)1h,相較于流體模型,固體模型的模擬時(shí)間可忽略不計(jì)。
4.4 結(jié)果與討論
陽(yáng)轉(zhuǎn)子旋轉(zhuǎn)大約5周后,結(jié)果開(kāi)始呈現(xiàn)出周期性。圖4所示為流體和轉(zhuǎn)子之間的瞬時(shí)熱通量和平均熱通量。最終,轉(zhuǎn)子壁面最大瞬時(shí)熱通量約為400W,從機(jī)殼外表面帶走的平均熱通量約為100W。
陽(yáng)轉(zhuǎn)子、陰轉(zhuǎn)子和機(jī)殼的最終平均溫度分別為345.3K,349.0K和329.6K。圖5(a)所示為轉(zhuǎn)子橫截面溫度分布,圖5(b)所示為轉(zhuǎn)子表面溫度分布。圖例范圍為300K到400K,品紅色代表高溫,藍(lán)色代表低溫。固體內(nèi)部溫度從入口至出口呈分層分布,溫度由低到高。
為了說(shuō)明流體和固體之間共軛傳熱的效果,本文將流固交界面設(shè)為絕熱邊界,其余設(shè)置保持不變,以進(jìn)行比較。圖6為交界面為絕熱邊界狀態(tài)下轉(zhuǎn)子在5個(gè)不同的曲軸轉(zhuǎn)角下的溫度分布計(jì)算結(jié)果,由圖可見(jiàn)瞬時(shí)溫度不再呈現(xiàn)出分層分布的特點(diǎn),在每個(gè)轉(zhuǎn)角下具有相似的溫度分布。而且,溫度范圍明顯更高。這意味著由于金屬具有極大的熱慣性,轉(zhuǎn)子的表面溫度實(shí)際上更溫和,更均勻,并且呈現(xiàn)為分層分布,這與絕熱壁面假設(shè)模型有很大的區(qū)別。
圖7所示為共軛傳熱(CHT)條件下5個(gè)曲軸轉(zhuǎn)角下轉(zhuǎn)子的壓力云圖,圖例范圍為1bar到2.5bar,品紅色代表高壓,藍(lán)色代表低壓。與預(yù)想的一致,每個(gè)轉(zhuǎn)角下的壓力水平相似。當(dāng)腔內(nèi)流體從入口流至出口時(shí),由于體積減小而引起壓力升高。與溫度分布不同,對(duì)于考慮或不考慮共軛傳熱(CHT)的情況,轉(zhuǎn)子表面的壓力分布幾乎相同。這意味著共軛傳熱(CHT)對(duì)壓縮機(jī)性能的影響可能很小。
表2比較了考慮和不考慮共軛傳熱(CHT)情況下的氣體質(zhì)量流量和轉(zhuǎn)子功率。從表中可以看出,兩種情況下氣體質(zhì)量流量和轉(zhuǎn)子功率的差異小于1%,與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比發(fā)現(xiàn),兩種情況下的質(zhì)量流量仿真值均高于試驗(yàn)結(jié)果約4%至5%,這種差異可能是由于間隙尺寸的誤差導(dǎo)致的。功率的預(yù)測(cè)值與實(shí)驗(yàn)數(shù)據(jù)的差異約為1%。因此,共軛傳熱對(duì)本文所用螺桿壓縮機(jī)性能的影響較小。
Simerics-MP+具備單向流固耦合求解能力,根據(jù)固體溫度模擬結(jié)果,對(duì)螺桿壓縮機(jī)固體熱應(yīng)力和熱膨脹變形進(jìn)行求解。圖8所示為由固體熱膨脹引起的轉(zhuǎn)子徑向位移,圖例范圍為0到50微米,品紅色代表大位移,藍(lán)色代表小位移。由圖可見(jiàn),徑向最大位移約為50微米。但值得注意的是,當(dāng)前的熱膨脹變形是單向耦合模擬所得,熱膨脹結(jié)果并沒(méi)有反饋到流體模擬中。
05結(jié)論與展望
5.1 結(jié)論
本文采用混合時(shí)間尺度耦合法(Mixed Timescale Coupling method)進(jìn)行螺桿壓縮機(jī)的流固共軛傳熱CFD分析,該方法可以有效解決時(shí)間尺度差異的問(wèn)題,模擬運(yùn)行時(shí)間與不考慮CHT的模擬時(shí)間相當(dāng),計(jì)算結(jié)果快速精確;
為了驗(yàn)證CHT對(duì)于螺桿壓縮機(jī)性能的影響,本文分別進(jìn)行了考慮和不考慮CHT條件的兩種工況模擬,結(jié)果表明CHT對(duì)壓縮機(jī)性能的影響很小,因此不考慮CHT的模擬是可以接受的;
根據(jù)固體溫度分布結(jié)果,利用軟件內(nèi)置的應(yīng)變應(yīng)力求解器進(jìn)行了單向耦合的熱應(yīng)力和熱膨脹模擬,結(jié)果穩(wěn)定可靠。
5.2 展望
未來(lái)將完成CHT模擬的集成封裝,以進(jìn)一步減少用戶(hù)的設(shè)置工作;
在熱膨脹模擬中將進(jìn)一步與流體模擬所得的幾何改變量相耦合,自動(dòng)考慮運(yùn)行過(guò)程中由CHT引起的所有變化,以提高模擬精度。