氣體離心機(jī)二維供料射流流場(chǎng)數(shù)值模擬

氣體離心機(jī)二維供料射流流場(chǎng)的數(shù)值模擬

姜東君 　王黎明 　曾實(shí)

(清華大學(xué)工程物理系  北京100084)

    摘要：為研究氣體離心機(jī)供料射流的流動(dòng)狀況,采用數(shù)值求解二維Navier-Stokes方程組的方法對(duì)其進(jìn)行了模擬,計(jì)算得到了不同供料角度下的流場(chǎng)分布。結(jié)果表明,供料射流對(duì)離心機(jī)流場(chǎng)在徑向和軸向上的分布都有影響;供料角度為0°時(shí),氣體在軸向上反復(fù)膨脹、壓縮,形成膨脹波和激波,對(duì)徑向的影響較小;供料角度為45°和90°時(shí),射流對(duì)軸向的影響減弱,對(duì)徑向流動(dòng)的影響顯著增強(qiáng)。由計(jì)算結(jié)果可知,離心機(jī)環(huán)流的計(jì)算中應(yīng)當(dāng)對(duì)供料條件進(jìn)行修正。

    關(guān)鍵詞：供料射流;離心機(jī);二維Navier-Stokes方程組;激波

    中圖分類號(hào):O35文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):1000-0054(2008)08-1287-03

    為了提高氣體離心機(jī)的分離性能,通常在其內(nèi)部引入二次環(huán)流[1]。供取料驅(qū)動(dòng)是較為常用的驅(qū)動(dòng)方式,由于取料口位于離心機(jī)側(cè)壁附近,因此關(guān)于取料對(duì)二次環(huán)流影響的研究比較深入,對(duì)取料器附近流場(chǎng)的研究也比較多[2-4],而對(duì)位于中心區(qū)域的供料射流的研究相對(duì)較少。由于供料射流影響離心機(jī)內(nèi)部流場(chǎng)的粘性邊界,因而對(duì)環(huán)流的影響也必須加以考慮。在以前的研究中,供料驅(qū)動(dòng)條件往往是根據(jù)供料量和一定的假設(shè)在粘性邊界上給出的,如采用均勻亞聲速流動(dòng)條件[5],或者采用滿足一定函數(shù)分布的源項(xiàng)[6]。以上各種假設(shè)不能夠反映真實(shí)的流動(dòng)狀況,因此,有必要對(duì)氣體離心機(jī)供料射流進(jìn)行研究。

    離心機(jī)內(nèi)的供料射流同工程中遇到的射流問(wèn)題有較大不同,工程中射流的周圍為均勻流場(chǎng),在沒(méi)有供料射流時(shí),離心機(jī)內(nèi)的流場(chǎng)滿足旋轉(zhuǎn)剛體假設(shè)[1],即壓強(qiáng)和密度在徑向上呈指數(shù)分布,速度為線性分布,當(dāng)工作氣體射入離心機(jī)后,將受到強(qiáng)旋流場(chǎng)的影響,一方面受徑向壓強(qiáng)梯度的影響,在徑向上的膨脹波被壓縮;另一方面受到旋轉(zhuǎn)的影響而具有角向動(dòng)量。

    本文通過(guò)數(shù)值求解二維軸對(duì)稱Navier-Stokes方程組,模擬離心機(jī)內(nèi)的二維軸對(duì)稱供料射流流場(chǎng),并對(duì)不同供料角度下的流場(chǎng)進(jìn)行了分析,為實(shí)際離心機(jī)供料射流流場(chǎng)的研究打下了基礎(chǔ)。

    1　計(jì)算模型

    在本文的計(jì)算中,離心機(jī)參數(shù)按照Iguacu模型[5]給定,其轉(zhuǎn)子半徑為6cm,長(zhǎng)徑比為8,側(cè)壁線速度為600m/s,即旋轉(zhuǎn)角速度為104rad/s,側(cè)壁壓強(qiáng)為13.3kPa。采用二維軸對(duì)稱供料模型,如圖1所示。本文計(jì)算了3種供料模型,其區(qū)別在于射流角度的不同,如圖1a、1b和1c所示,射流角度α分別為0°、45°和90°。

    2　計(jì)算方法

    2.1　控制方程控制方程為二維軸對(duì)稱Navier-Stokes方程組:

    由于在射流出口為超聲速,流動(dòng)為紊流狀態(tài),在計(jì)算中選用Realizablek-ε湍流模型。

    2.2　計(jì)算方法

    在計(jì)算中采用有限體積法對(duì)控制方程進(jìn)行離散,計(jì)算中使用二階迎風(fēng)格式,隱式迭代。

    3　邊界條件和初始條

    件離心機(jī)處于高速旋轉(zhuǎn)狀態(tài),因此在計(jì)算中采用旋轉(zhuǎn)坐標(biāo)系,旋轉(zhuǎn)角速度為1×104rad/s。

    計(jì)算中用到的邊界如圖1所示,具體設(shè)置如下。

    供料入口:給定來(lái)流的壓強(qiáng)、速度和溫度。

    轉(zhuǎn)子側(cè)壁:旋轉(zhuǎn)壁面,給定壁面的旋轉(zhuǎn)速度(600m/s)和溫度(300K)。

    壓力出口:按照剛體旋轉(zhuǎn)條件[1]設(shè)定背壓:

    底部端蓋:旋轉(zhuǎn)的壁面條件,給定壁面的旋轉(zhuǎn)速度和溫度。

    在計(jì)算中,以無(wú)供料時(shí)的離心機(jī)流場(chǎng)(即旋轉(zhuǎn)剛體假設(shè))作為初場(chǎng),壓強(qiáng)分布見式(1),密度為

    4　計(jì)算結(jié)果及分析

    計(jì)算中使用的氣體為UF6,其相對(duì)分子量為352。入口條件給定為500Pa,溫度為300K,速度為0。在此供料條件下,計(jì)算得到了射流角度分別為0°、45°和90°3種供料模型的流場(chǎng)。

    圖2給出了3種供料模型下的等壓線分布,相應(yīng)的等馬赫數(shù)線分布如圖3所示。在供料角度為0°情況下,當(dāng)氣體進(jìn)入離心機(jī)后,迅速膨脹,在徑向上受到離心機(jī)壓力梯度的影響,被迅速壓縮;在軸向上存在壓縮-膨脹-再壓縮-再膨脹的過(guò)程。當(dāng)供料角度為45°時(shí),氣體在徑向上的膨脹壓縮過(guò)程較為明顯,軸向上的分布已經(jīng)明顯減弱。在供料角度為90°的情況下,氣體在進(jìn)入離心機(jī)內(nèi)之前,在管口處形成了激波,經(jīng)過(guò)激波后,氣體的速度明顯降低,可以看到,形成的膨脹-壓縮波基本上是以出口為中心的圓弧,在軸向上呈對(duì)稱分布,與供料角度為45°的情況相比,對(duì)離心機(jī)側(cè)壁處的影響減弱。在圖3中可以觀察到同樣的流動(dòng)現(xiàn)象。

    對(duì)于不同的供料模型,圖4分別給出了在半高處(x=0.24m)徑向上的Ma分布。從圖中可以看到,在r=0.04～0.06m之間,在0°供料角條件下,Ma基本保持線性分布;供料角為45°和90°兩種情況下,Ma的分布基本相同,也呈線性分布,但其分布曲線位于供料角為0°的曲線下方。這表明當(dāng)供料角度增大后,徑向速度增大,從而對(duì)離心機(jī)流場(chǎng)施加了更大的影響。

    供料射流對(duì)離心機(jī)流場(chǎng)的最主要影響是對(duì)近側(cè)壁處環(huán)流的干擾,因此環(huán)流計(jì)算內(nèi)邊界處,軸向和徑向速度在軸向上的分布是本文所關(guān)心的。圖5和圖6給出了半徑r=0.048m處軸向上的徑向和軸向速度分布曲線,從中可以觀察到,對(duì)于文中計(jì)算的3種供料條件,徑向速度在軸向上的分布比較復(fù)雜,在3種供料角度下,都存在不均勻分布。但一般離心機(jī)環(huán)流計(jì)算經(jīng)常是將供料假設(shè)為一股均勻的亞聲速流動(dòng),這是應(yīng)當(dāng)修正的。在供料角為45°和90°時(shí),軸向速度則基本上以半高處為中心對(duì)稱分布。

    5　結(jié)論和討論

    本文通過(guò)求解二維軸對(duì)稱Navier-Stokes方程組模擬了二維離心機(jī)射流流場(chǎng),采用有限體積法離散,使用二階隱式迎風(fēng)格式迭代計(jì)算,得到了不同供料角度下的流場(chǎng)分布。對(duì)計(jì)算結(jié)果的分析表明:

    1)射流對(duì)離心機(jī)流場(chǎng)在徑向上和軸向上的分布都有影響;

    2)在供料角度為0°時(shí),氣體在軸向上反復(fù)膨脹、壓縮,對(duì)徑向上的流場(chǎng)影響較小;

    3)在供料角為45°和90°時(shí),射流對(duì)軸向的影響減弱,對(duì)徑向流動(dòng)的影響顯著增強(qiáng);

    4)無(wú)論在何種供料條件下,供料射流所造成的流動(dòng)參數(shù)(包括徑向速度和軸向速度)在軸向上的變化,在離心機(jī)環(huán)流計(jì)算中都應(yīng)當(dāng)被考慮。

    本文計(jì)算的是二維情況下供料射流對(duì)離心機(jī)流場(chǎng)的影響,實(shí)際的供料模型更為復(fù)雜,需要進(jìn)行三維數(shù)值模擬。