供料位置對離心機(jī)供料射流流場影響的數(shù)值模擬

張雯成，姜東君，曾 實(shí)

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

供料位置對離心機(jī)供料射流流場影響的數(shù)值模擬

張雯成，姜東君，曾 實(shí)

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

為研究供料口位置對離心機(jī)供料射流流場的影響，采用直接模擬蒙特卡羅（DSMC）方法對離心機(jī)中心射流區(qū)域進(jìn)行模擬。選擇Iguassu離心機(jī)模型和變徑硬球碰撞模型，得到了不同的軸向和徑向供料位置下徑向供料模型的供料射流流場參數(shù)分布情況。通過對模擬結(jié)果進(jìn)行分析可得：在徑向供料條件下，供料口附近區(qū)域的流場受供料條件的影響顯著。由于供料射流的膨脹作用，導(dǎo)致了供料口附近低溫區(qū)的出現(xiàn)。在靠近出流邊界處，由于較高的角向速度，形成較大的壓力梯度，使氣體的徑向速度降低，導(dǎo)致溫度明顯升高。軸向位置的變化導(dǎo)致流體參數(shù)的軸向平移。而隨著供料口徑向位置的增大，供料口附近的膨脹范圍減小，低溫區(qū)減弱，即該區(qū)域的最低溫度升高。而近出流邊界處，隨供料口徑向位置的增大，溫度相應(yīng)降低，徑向速度同時降低。

氣體離心機(jī)；稀薄流區(qū)；數(shù)值模擬

氣體離心機(jī)是一中空的轉(zhuǎn)子，內(nèi)部充入UF6同位素混合物氣體。在轉(zhuǎn)子高速旋轉(zhuǎn)形成的強(qiáng)離心場作用下，充入的氣體產(chǎn)生分離效應(yīng)。在無內(nèi)部環(huán)流的情況下，其內(nèi)部氣體的壓強(qiáng)在徑向上呈指數(shù)分布［1］。按側(cè)壁壓強(qiáng)13.3kPa計算［2］，對于Iguassu模型離心機(jī)［3］，努森數(shù)Kn的范圍為4.33×10-4～4.67×107?？梢妼τ谠撃Ｐ碗x心機(jī)，流場在徑向上跨越了自由分子流區(qū)、過渡區(qū)、滑移流動區(qū)和連續(xù)介質(zhì)區(qū)等4個流區(qū)，其中自由分子流區(qū)、過渡區(qū)和滑移流動區(qū)統(tǒng)稱為稀薄流區(qū)。

對于離心機(jī)中心區(qū)域的稀薄流區(qū)，由于連續(xù)性假設(shè)不成立，Navier-Stokes方程不再適用。因此，在離心機(jī)的流場計算中，常忽略中心的稀薄區(qū)域，僅以滿足Navier-Stokes方程的連續(xù)流區(qū)為研究對象，并引入黏性邊界，人工給定供料條件。但這一方法并不符合實(shí)際的流動情況，所以需對離心機(jī)內(nèi)部的稀薄流區(qū)進(jìn)行模擬，來得到供料射流的具體結(jié)構(gòu)。

國內(nèi)外學(xué)者均對該領(lǐng)域進(jìn)行過研究。Roblin等［4］采用DSMC（Direct Simulation Monte-Carlo）方法，對徑向供料模型進(jìn)行了模擬，得出了流場中的速度分布情況。姜東君等［5］采用CFD方法，同樣模擬了徑向供料的Iguassu模型。目前，該方面的研究主要集中在徑向供料的流場分布和供料速度等參數(shù)對流場的影響，對供料位置等其他參數(shù)研究得較少，因此需進(jìn)行進(jìn)一步的研究。

本文采用DSMC方法，對高速旋轉(zhuǎn)離心機(jī)內(nèi)部的供料射流區(qū)域進(jìn)行模擬，針對徑向供料模型，研究其不同的軸向和徑向供料位置對供料射流流場的影響。

1 DSMC方法簡介

DSMC方法是目前模擬稀薄流區(qū)較為合適的方法。該方法直接對流動進(jìn)行模擬，由Bird［2］首先提出。其思想是將分子的運(yùn)動和碰撞解耦，并用模擬分子（宏分子）代替真實(shí)的氣體分子（每一個模擬分子代表大量的真實(shí)分子）。在每個時間間隔內(nèi)，氣體分子無碰撞地按照自身速度和當(dāng)前位置自由運(yùn)動。然后根據(jù)分子的位置，運(yùn)用概率的方法來判斷分子與分子之間或分子與壁面之間是否發(fā)生碰撞。根據(jù)分子間的碰撞模型和分子與壁面間的碰撞模型，可計算得到分子碰撞后的狀態(tài)。流場中宏觀量（如速度、密度、溫度等）可通過對網(wǎng)格中的模擬分子進(jìn)行統(tǒng)計平均得到。在DSMC方法中，分子間的碰撞模型是比較關(guān)鍵的問題，本文中所采用的是變徑硬球模型。

2 模型及條件

2.1 離心機(jī)模型

本文采用Iguassu模型離心機(jī)［3］，其主要參數(shù)列于表1。

表1 Iguassu模型離心機(jī)［3］主要參數(shù)Table 1 Parameters of Iguassu model centrifuge［3］

2.2 計算模型

計算模型如圖1所示，采用徑向供料模型，氣體由沿徑向的供料口進(jìn)入離心機(jī)內(nèi)部，在供料口處給定供料射流的徑向速度、溫度和密度，不考慮供料射流的軸向和角向速度。由于離心機(jī)徑向上的密度變化很大，因此越靠近側(cè)壁，所需的模擬分子數(shù)就越多。受計算機(jī)計算能力的限制，目前尚無法開展對連續(xù)流區(qū)的DSMC模擬。因此，計算域選擇為供料口附近區(qū)域，徑向范圍為0.005～0.05m，軸向范圍為-0.12～0.12m。

圖1 計算模型Fig.1 Calculation model

2.3 氣體參數(shù)

本文中所選用的氣體介質(zhì)為UF6，其相關(guān)參數(shù)列于表2。

表2 UF6氣體相關(guān)參數(shù)Table 2 Physical properties of UF6

2.4 邊界條件及初始條件

模擬中，無供料條件下的模擬區(qū)域初始狀態(tài)設(shè)為均勻流，而在徑向供料條件下則設(shè)為真空。

邊界條件的設(shè)定如圖2所示。在無供料條件下，模擬區(qū)域的下邊界為固體壁面，而在徑向供料條件下，模擬區(qū)域的下邊界為入流邊界，供料量為12mg／s，供料氣體的溫度為300K，馬赫數(shù)為1。供料管固定，溫度為300K。模擬區(qū)域的上邊界為出流邊界，氣體的分子數(shù)密度為1.37×1021m-3，溫度為300K，旋轉(zhuǎn)線速度為500m／s。模擬區(qū)域的左右邊界設(shè)為旋轉(zhuǎn)固體壁面，溫度為300K，旋轉(zhuǎn)角速度為1×104rad／s。

圖2 邊界條件Fig.2 Boundary condition

3 計算程序

本文所用程序是在Bird［2］提供的G2程序基礎(chǔ)上進(jìn)行的再度開發(fā)，使其能適應(yīng)離心機(jī)流場的特殊需求。對于軸對稱流場，均勻的網(wǎng)格劃分會導(dǎo)致靠近軸線的網(wǎng)格體積小于靠近出流邊界的網(wǎng)格體積，再加上離心機(jī)中心區(qū)域極為稀薄的特點(diǎn)，靠近軸線區(qū)域的網(wǎng)格內(nèi)分子數(shù)會遠(yuǎn)少于靠近出流邊界處的。如果每個模擬分子所代表的真實(shí)分子數(shù)N為固定值，就會導(dǎo)致中心區(qū)域內(nèi)的分子數(shù)過少，從而碰撞數(shù)過少，不利于計算的收斂和宏觀量的統(tǒng)計平均。

為了解決這一問題，計算中設(shè)定模擬分子所代表的真實(shí)分子數(shù)N與半徑r的關(guān)系為：N＝krn。其中，k為常數(shù)。通過實(shí)際模擬計算，選取n＝3。

4 模擬結(jié)果及分析

本文首先對無供料情況進(jìn)行模擬，在此情況下，流體狀態(tài)存在解析解，在半徑方向上氣體壓強(qiáng)呈指數(shù)分布，速度呈線性分布v＝Ωr。圖3示出了無供料條件下z＝0處沿徑向的流動參數(shù)分布情況。可看出，密度、角向速度分布與解析解符合得較好，從而驗(yàn)證了本文所編DSMC計算程序的正確性。

圖3 無供料情況下z＝0處的流體參數(shù)分布Fig.3 Distribution of gas properties on z＝0on condition of no feed

固定徑向供料位置r＝0.005m，計算不同的軸向供料位置z為0、-0.01和-0.02m3種情況下的流場分布情況。對于實(shí)際離心機(jī)，軸向供料位置不可能位于接近端蓋處。因?yàn)榻咏松w處，受擋板、支臂等影響較大，在該處供料，離心機(jī)的分離功率較低［6］，對實(shí)際工作無意義。因此，本文選取了較為靠近中心處的3個軸向供料位置。圖4為在供料口軸向位置z為0時，流場參數(shù)的分布情況?？煽闯觯┝峡诟浇鼌^(qū)域流場受徑向供料的影響極為顯著。由于供料射流的膨脹作用，導(dǎo)致了供料口附近低溫區(qū)的出現(xiàn)。在近出流邊界處，由于較高的角向速度，形成很大的壓力梯度，使氣體的徑向速度降低，導(dǎo)致溫度明顯升高。而圖5示出3種不同的軸向供料位置下，r＝0.03m處的流體參數(shù)分布情況?？煽闯?，在供料口附近的一定軸向范圍內(nèi)，軸向位置的變化對流場的影響不甚明顯，僅使得參數(shù)的值隨供料口軸向位置的平移而平移。

圖4 供料口軸向位置z＝0時的流場參數(shù)分布Fig.4 Distribution of gas properties on condition of feed on z＝0

圖5 r＝0.03m處的流體參數(shù)分布Fig.5 Distribution of gas properties on r＝0.03m

固定軸向供料位置z＝0，計算不同的徑向供料位置0.005、0.010和0.025m的3種情況下的流場分布情況。圖6為3種情況下r＝0.05m處的流體參數(shù)分布情況?？煽闯?，隨供料口徑向位置的增大，由于軸向可擴(kuò)散的范圍減小，供料口附近的膨脹范圍減小，低溫區(qū)減弱。而近出流邊界處，隨供料口徑向位置r的增大，溫度相應(yīng)降低，徑向速度同時降低。

圖6 r＝0.05m處的流體參數(shù)分布Fig.6 Distribution of gas properties on r＝0.05m

5 結(jié)論

本文采用DSMC方法，對氣體離心機(jī)供料射流區(qū)域進(jìn)行了模擬，給出了無供料和徑向供料兩種條件下該區(qū)域的流體參數(shù)分布情況。本文實(shí)現(xiàn)了DSMC方法對離心機(jī)供料射流區(qū)域的模擬，并通過對模擬結(jié)果進(jìn)行分析，得到了以下結(jié)論：

1）在無供料條件下，離心機(jī)中心區(qū)域流場基本符合等溫剛體模型。其中，氣體壓強(qiáng)沿徑向呈指數(shù)分布，角向速度呈線性。

2）在徑向供料條件下，供料口附近區(qū)域的流場受供料條件的影響顯著。由于供料射流的膨脹作用，導(dǎo)致了供料口附近低溫區(qū)的出現(xiàn)。在近出流邊界處，由于較高的角向速度，形成較大的壓力梯度，使氣體的徑向速度降低，導(dǎo)致溫度明顯升高。

3）供料口附近一定軸向范圍內(nèi)，軸向位置的變化對流場的影響不甚顯著，僅使得參數(shù)的值隨供料口軸向位置的平移而平移。

4）隨供料口徑向位置r的增大，供料口附近的膨脹范圍縮小，低溫區(qū)減弱，即該區(qū)域的最低溫度升高。而近出流邊界處，隨供料口徑向位置r的增大，溫度相應(yīng)降低，徑向速度同時降低。因此，供料口的徑向位置r對離心機(jī)的流場影響相對較為顯著。

通過以上分析可看出，供料口的位置變化會給黏性邊界處的參數(shù)分布帶來較大影響。該影響會進(jìn)一步導(dǎo)致離心機(jī)分離功率的變化。因此，進(jìn)一步的工作應(yīng)將射流區(qū)域的流場與連續(xù)流區(qū)的流場進(jìn)行耦合計算，其中，連續(xù)流區(qū)的流場采用CFD的方法進(jìn)行計算，結(jié)合外部連續(xù)流區(qū)的環(huán)流進(jìn)行綜合分析，得出供料射流的位置對離心機(jī)分離功率等指標(biāo)的影響。

［1］ 張存鎮(zhèn).離心分離理論［M］.北京：原子能出版社，1987.

［2］ BIRD G A.Molecular gas dynamics and the direct simulation of gas flows［M］.Oxford：Clarendon Press，1994.

［3］ DONEDDU F，ROBLIN P，WOOD H G.Optimization studies for gas centrifuges［J］.Separation Science and Technology，2000，35：1 207-1 221.

［4］ ROBLIN P，DONEDDU F.Direct Monte-Carlo simulation in a gas centrifuge［C］∥Rarefied Gas Dynamics：22ndInternational Symposium.Sydney，Australia：American Institute of Physics，2001：169-173.

［5］ 姜東君，曾實(shí).氣體離心機(jī)二維供料射流的數(shù)值研究［J］.原子能科學(xué)技術(shù)，2008，42（增刊）：41-44.

JIANG Dongjun，ZENG Shi.2Dnumerical study of feed-jet flow in gas centrifuge［J］.Atomic Energy Science and Technology，2008，42（Suppl.）：41-44（in Chinese）.

［6］ 劉冰，姜東君，曾實(shí).供料條件對氣體離心機(jī)分離性能影響的數(shù)值研究［J］.原子能科學(xué)技術(shù)，2012，46（10）：1 153-1 157.

LIU Bing，JIANG Dongjun，ZENG Shi.Numerical investigation on influence of feed condition on separation performance of gas centrifuge［J］.Atomic Energy Science and Technology，2012，46（10）：1 153-1 157（in Chinese）.

Numerical Simulation of Influence of Feed Position on Flow Field

ZHANG Wen-cheng，JIANG Dong-jun，ZENG Shi
（Department of Engineering Physics，Tsinghua University，Beijing100084，China）

In order to study the influence of feed position on the counter-current in gas centrifuge，direct simulation Monte-Carlo（DSMC）method was chosen to simulation the rarefied region of gas centrifuge.By setting Iguassu centrifuge model and variable hard sphere（VHS）model，the flow distributions of different feed positions were acquired.The analysis of the calculation results notes the following flow phenomena：In case of radial feed，the region near the feed outlet is largely influenced by the feed condition.The expansion of jet results in the hypothermic region.Near the boundary of effluent，due to the fast peripheral velocity，pressure gradient is large，which results in the slow radial velocity and the high temperature.The axial translation of the jet results in the axial translation of the flow distributions.With the increase of the radius of feed，the region of expansion near the feed reduces and the lowest temperature of the hypothermic region rises.Near the boundary of effluent，with the increase of the radius of feed，thetemperature and the radial velocity fall.

gas centrifuge；rarefied region；numerical simulation

TL25

A

1000-6931（2014）02-0331-05

10.7538／yzk.2014.48.02.0331

2012-11-11；

2013-01-05

國家自然科學(xué)基金資助項(xiàng)目（11075088）；中俄（NSFC-RFBR）合作項(xiàng)目資助（11311120046）

張雯成（1989—），女，上海人，碩士研究生，核燃料循環(huán)與材料專業(yè)