· a c# r ˇp ˚ ˇ ˚ o˚ ˇ> ... e˘u˚4 ˘0 7 k e˘u˚4 ˘0 ˙b# r. ˇ ˆ ˘ ˇ ˘ cfd 0...

!

CFD

! "#$%

&'()+

* * "+,-# ../#+,-01$&#

2 34$54 %6+ % ! "#$%

789:4' +;$*< ' =/' /# ' <. )$3> %? @ +

0 $+ ;$*< A /$' 4B ' 8 ! C D?+,-E )<A$ $+ @=#F. ''/ - G 4 @;$*< *< 'E $:-

HCFD/$I%% C @7 '! ';$*< >?>J/'/# '( 3J >J K(MRF)' '/ E $@ =+RNG k-ε RSM % *< # I E C

/&$<8 /' +E $ =4 "G ! . K# L % # I E @;$*< % ' C /5 )/&$<8/' +/$I% CFDE *< 'H <M

RNG k-εRSM E $! % .!D$'1$J 9N<>/=+/# 'D' ' E 7 ' # I /1+ /'+RSM E +RNG k-ε@, ' ' + + O @;$*< % ' C/&$<8/' + # I E E $?>J/=+

D$P+Q> ;$*< D$ -7 ' R> +& <> 3$ST##0%/' @;$*< E $=4+C

/$I% @;$*< CFD >JK( 3 @A$ $+C

KEY WORDS: Gas centrifuge, CFD simulation, Hydrodynamic, Multiple Rotating Reference

Frame (MRF).

*0 I# !+E-mail: [email protected]

"#$%&'('

)

!"# $ %& '&()*+ ) , !"# $ -&

%& ./0123& 4 2 5$26& 12 7&

/* 2 ) , #'/&8+!9: %& /'H2S2CO2

;& <=>*+)?@ ,A B/+'/& C!9: CO22H2S <=>

D 2 /& C /+ '! E & % /) , C 2 .F' /+ '/&

BC)$(05 /&6'! G H/&;+; 4 I&6$ '2)]J@[,

$.F)'/+'/& C & /&8+!9: %& /& C'

5K),%& ;2F L&$#/)22('

6MH$7& N(FONP$ 2: > &K2 Q&8F /N3&;& : &# R

- 23/&2'S5 & 2 2+, %& 0 ' 2 3/&2'

BTUBV6 2ID>/''%& , 7& (F > NF 2 ># W7& 2 :

F<) 0 2 < KG> -(BR X$/,><22<3&ONP$

( : $ =< O: .& Y$3& ;/F&', 1"$+2;&O<D2F L&$ H$

;=/- &0)P%& L&$B# R;& 2/,

%& 3& KD7(U$ &H+ ;& - Z">8 ; -B'2[5$<\

;& 3/&2 H /# R M+G'2[53(# R:6L

12'C B5F,O:)9< NNP$''4:N0 &

M/$ =& Y$ KG + $ -(BR 7& ]$ :7& 7'^@_)H?` ,A7& -B'2[5

;& O > a C# R P /<&b ! G F<) F c ,7/S('

d<=& Y$O: >#H'3&L&$3&M/$< cN* & ,7($N 2 / ^J_

@e`"]f[0(';& 2POQ$12=)%& 3&P;Y$ ,M+

<&'KEU\# *g"KG O: P ;Y$ a C , @ee"

h5 2;'^`_]i[<\ ';2 ;& 2 #cP 3& 8 F !\< # RO: 2 < F ' ' &2 jk

4 $: H0 &95 , M+H'FF #5 F<) 2P'/'

c 0&Y ) H , lm$ 7/S(''F 2 %& Z$ P ' ';

),'^f_@e"]n[ lbF3&O: &2' k ,@ei"

6&-+2NV ^i_]o[O:<F'' &2k lpB &5:'\# *& ,

22-'^n_J"][ 2 &NC6-' %& ,<M+'

4 $ 3& ;& 4: &952 N O: $q, H';& r\#+

&) ;& r\# PG 2 >#OC 4';& BU2)! >s2 # H2 C

;& PG> ;& (DU 2 &OC;& 2! P/7&5! >/N&N)

2H N2s2 # Hs U(05),H't'2% ;B M+;& S5 - 3&

H&= &2;& >U= Y\/),

UVW Cohen UXW Harada

UYW Cloutman and Gentry UZW Ungurish and Israeli

U[W Park and Hyun U\W Hasegawa

&#*+$,-./0'1#,-/CFD

u; 2u)^@_Jn"]e[=)2%& 3&<;& ;'2[5M+

<O: 2P4H:12' <, ;& ; H'=) 7& -

B5s2 # H# :;& (05 X$t&2'> :\#/

5 7& & 0(' ;& H& OF ) ,7/S(' M+& M 7 s2 # H& 2;# R ^J_

%& <$H&;$ 2h) lm$ )';# R&:;+BC&6CvF))

OF7& 24'2 2 BC;& ;$ 2<$;#%& ,

=) CFD/5 & 6&5 2 3/&' C Y&/ S5 ')

[&2&'N2;& C/52/C')/&6' 52OB'-B&+,12=M2F'

C B52/CH ; 6C U CFDC&',3&7& /H/=FCFD)<$ $ONP$ >#0:

!\<)' 5 X$ &95, -&'2[5' #%& ;& C&/) 'w#

OC l c (' t 2 ';B '',7&7& /'2[5ID2 B ;& C $

3/&2' 2 /' 2 $N(F ]& ) lm$ 2%& CFD D ,

! "#

$%CFD%&'()

< (0D 7& ' : 2 -5 7L xD^`_&2)2 '

y& B$ 2 2 h P : ) 2 )<'xD 4:;# RL ^f_ 21B*

D ,

%&*(+,-.

xD ;& % 2 / -5 <7L &KG/),

^@_( ). v 0t

∂ρ+ ∇ ρ =

∂

^J_( )

( ) ( )v

. vv P . gt

∂ ρ+ ∇ ρ = −∇ + ∇ τ + ρ

∂

^`_C T v. T .qt

∂ ρ + ∇ = −∇

∂

P3$BCτ2/$0/$gρ

B 2 ),

^f_( )T 2v v .vI

3

τ = µ ∇ + ∇ − ∇

%&/,-. 012345

;& 07)$//M4M'';6U& B# R ;& %& 2 r"# '

22 24&) ,< '> a C 5 7 KEU ^0& & _ O: 2 )z# ;'

<')O: ::KEU, B'(0D7KEUD2 ::'

-7&4/:N0''&5&/',2 (0D! >h P L xD3& xU: '

:/;$N0'7&&5 # OV0, N0' 4))/:|' ;$

/ ;# R xD^n_(SRF) ./' 'S5 (: 7& $SRF0 7( ) 7&'40Y$$=P:7&/RN02)

&122:7&C ),

•;# RL xD/R(MRF)^o_

•PG'>"#(MPM)^_

PG&'>"#(MPM) 2O; $/(' 7C c 27':'2lh P

;$=P&8=N>,

UVW Glombok and Bill UZW Swirling and rotating flows

UYW Rotor U\W Single rotating reference frame

U[W Inertial reference frame U]W Multiple rotating reference frame

UXW Rrotating reference frame U^W Mixing plane model

"#$%&'('

2

7&&MRF) & 2 OD D 7&5 0 $=P ;=D

=)%& ,&3&6= : : ;# Rh P ) a C2 &:;+7

/0' , & 7& /MRF< '4: ;& :2%& 2:'7

: 7 l 2 ) O:' 2 7h P&),V6 (k5&7& $)z < 2'h P: 4:;) ,3&

&2 F KEUω

L xD3& =07>5O)# R@);B,

# R P2 N :2 ^& _ap& <$)W

î_ru r= ω×

r rv v u= −

aω = ω

< 'F &L xD '

O&=$ ;# R L xD 7) , rv

F

=0^F &;# RxD ) _ 2v

IN\F

^ F &7 xD ) _ 2ru

,# R F ^@_^h PxD.=F _/),

;# RL xD3& : 2 Y<p& <$ &KE),

^n_r. v 0t

∂ρ+ ∇ ρ =

∂

ô_( )

( )r

r

v. v v

t

∂ ρ+ ∇ ρ +

∂

( )v P . Fρ ω× = −∇ + ∇ τ +

^_( ) ( ) ( )r r

E. v H u . k T

t

∂ ρ+ ∇ ρ + ρ = ∇ ∇

∂

6708,/6

1B*; l y& B$ ./'cPF ' X2'<Y'Y),

$=H&0(' ''%& 1B*2<\7&k-RNG2

RSM l7C c ';2;& 1B*)%& ,

0"34&56789%2:;

=P ; ED /($D $ '7) 2 =D ODDk- ;& <2'4'

&/) 7&[&2';& 4'%& ) ,CQ 27& <$ 7/S(';

N( < CQ' RNG k-2 MD ;/(>OD;& $ 2B# R'F C&

z &+ , ' 51B*I\ 3/&2' k-<'=P &)W

ê_( ) ( )i

i

k ku

t x

∂ ρ ∂ ρ+ =

∂ ∂

k eff k b Mj j

kG G Y

x x

∂ ∂α µ + + − ρε − ∂ ∂

^@_( ) ( )i

i

u

t x

∂ ρε ∂ ρε+ =

∂ ∂

( )2

eff 1 k 2 b 2j j

C G C G C Rx x k k

ε ε ε ε ε

∂ ∂ε ε εα µ + + − − ∂ ∂

< '8Gk 3/ % $ 1B* \2;& 2]F 'Gb3/% $1B*\22/)G#2YM-)+r0=4M

4 $1B* N!"$ ~ &95,kα 2εα < lbO 5F.$ $k2ε,

/&/$(RSM)/$7<$ /&'< O:' !"$ ~ =P 2 Y()3&

-&CQ</,/<&;& <2'i ;& 2<'O:$KG<'/,

UVW Swirl Velocity

z$=P2P

# RP

KEU40

7

KEU40

B# R

&#*+$,-./0'1#,-/CFD

;RRSM1 ;& ]# /P l^@_# R ^J_2 X$' M ~ L& < 2'2

<3&x0: H 2+;& $'/ O: M KG S5 ,<' Y

; &KG/&/$=P)W

^@@_( ) ( )i j k i jk

u u u u ut x

∂ ∂′ ′ ′ ′ρ + ρ =

∂ ∂

( )i j k kj i ik jk

u u u p u ux

∂ ′ ′ ′ ′ ′− ρ + δ + δ ∂

( ) j ii j i j j k

k k k k

u uu u u u u u

x x x x

∂ ∂∂ ∂′ ′ ′ ′ ′ ′+ µ − ρ + −

∂ ∂ ∂ ∂

( ) jii j j j

j i

uug u g u p

x x

′∂′∂ ′ ′ρβ θ + θ + + − ∂ ∂

( )jik j m ikm i m jkm

k i

uu2 . 2 u u u u

x x

′∂′∂′ ′ ′ ′µ − ρΩ ε + ε

∂ ∂

N( ;B<R(2NPIB/'N( ! <R(2 &Cij&(/'N( 6( #+$2/';B.$ $

d(' $ /$ $ d 2/)G# $ 2 !"$2 # R N .$ $

'(DT,ijDL,ijPijGijφijεij2Fij&(/',

90)7%:;<:

7&=) <\ %& 3/&2' )H;Ch O0 ,&()

U ('-O);+KEJ2 2@;B),;('(0Dh#;& TUB

6 N#^Z$m Jf _ - 2), 2 H; s# 2 C8.$ $(052

\D s 7&5 2nZ$ 2 'fne`J2@eJ ) ,OC 7/S('2 8 &'

Z$%& 7&5'fi@2J@`C D , () I=>&~ )7&5 OC 2 &'

=<$ ;6&6Cc/),C ' X2]& )F=)2 BC O) )TUB$=P / 0&

7/S('2F Y'(r '('),

<%*&=+$,-&#;

(Y

%& Z$f

%& \Dn

N# \Dn

OCUQn

~ \D7&5OC'J

22 \D;& f

0"2=+$,-#&#;

& -('$/)ID&C:2]& )&+,

*+4c/$22;& BC<%& BC) 24c/$ 2 BC Y 6F I& >7&) , D 22 BC r ) 2 )

' 5(' 2 22 < 6),Y 22 BCPaf@ 4c/$) ,' 7/S('

2 #BC r )7&528 2 #)YPaf@i;& 2 C2 2 #Paf@

4c/$ 2 # (05 ) , r ) C =16X&'OC2N#p2pY7/S('2'

) C c %& , 5 UQ cm@ # R : ;/F 2 &

F %& )m/s`@),UVW Swirl UYW Rotation

"#$%&'('

>

<')O:2P4H:124:, >'upwind 2 2 2 =$ 00

<' ) ,-&4SIMPLEO: c/ $.F jPE$ 2 BCy) ,

$=P< 2 N@Ji@o) 00 N,S5 $=)H& ;2 C/''

)40Y$(0D7&/R%& ./,t'&8:27&528 2 #7&5OC7& &28OC

5 2=) > ) /, $)<& -('F 2D&6C B

) C& # $=P =) : ,5=)]$ %& 6 :O 2;# R'

)/ )) ,($=P =) O))`);B,

$=P ; : c/ ; 4' C 2H' =0DI0:<\ =)N<$xU'

M&)HH'FNY0=)2+ 2 , X$ 7C c ')&6 &

&'=)'& $6 &7:'), ;+ -('&@l :J $ '), -('&O:/=D ('' 0F (H:<$2& -

' 5 ) ,$=P ; $ 8 =)'CFD=) '

1 5^JfF 0'GHzfJ_)H,

= >

6? @A< A

=) <$ X$O:/ ) $=P'H'=)lm$P$ D',7&7& /

C D& $=P =) lU, >&+&6C)6&CD&=)

;+<$ BH'// X$,"Y H'+F=)=)'=)'

/= ]& )C &95HZ/R ,~FY ';O) $=P =)f )4, Z<)0'&62F ;7&5O4

0"$?1='"=#&;

# R : # & )'Y~ F ;0& : 2 7& /)2=) 7 K2$'

N#:&) ,&6C $ iF6' X$$=P=)<$'6R'Y~F

+ 2 ,7& /;$ C H7/R=) <$ Y $=P '@ifJ`i )H'

=)'#=)OY0,

+ /BC0D

ONP$ 2 &/)--R%& ;2 : # &8 (' ,(0D 7& 7& /

OC2 5';& F 2 BC ; D ,=) c/7& /= ]& )&'

/&#*+$,-.0'1#,-/CFD

@

<%2<$A1"=&$*$*&;

=)<$()

@Ji@oMesh1

@ifJìMesh2

JèoiJMesh3

0">"=B CD$*$*&;

^&# RF m/s`@7&528 2 #'BC

.$ $ f@i2f@ 522 BC 2 paf@_ )H,

H'=)7&KGOC2 5&'2F >L\YBC%& O)i;B),

;(' 5 $)# R lm$ &6C & 2 'YF :'# RN#:&&6 &&# RF #

),:7&.$ $7&# RF <m/s`@)7& B:7&F )2

Z<)' F & F<)KG #& $;# R'7N#& ,7&5 28

lm$NF%& 'OC';& PF ':& # R'YC&&6CN#F

7&F 2: 8,:'7&52&8%& Ds2 #' H vF# R

F &)' >#7('2):($>s2 # HOC2 56RYF ;B/',

OC2 5'O) BC'i 2 pnp ;BF O&C2 5I\'&BC;# R'

& $F<)KG 2 Y7& BN#'' I\ &H' )16'2[5'

7B5)?@eA , 'BID $'YBC ! >s2 #' HOC2 5BC'27&:'8

);BGE#> %& Ô)np_, :'! > H's2 #7&528s2 # H_Y\/'

BC4&ZQ 7&NF ; ;B )OC2 5'F ' 2BCÔ)njx_2s2 # H

)48OC2 5]2Q FL\YBCs2 # H 9O) np;B';&

B5BC4s2 # H+ 2,' F (MCB# R :);;

s2 # H2 :;& 2)! P/ Hs2 #2C&H =YFy\;& s2 # H

-# Y\/& 2 )^@_B5 s2 # H

27C NGC &+s2 # HF ' M C O)# R,7& /NFH&:'

! > BC4 H's2 ## R: # 28F B5 H&s2 # H)

'%& Nv:=&--RP : ) ,OC2 5 ZQ 7& F '

P1 'Q F;Z$2;' O) o ) 4 , ID 7C& :&'

OC2 5/ &=PF 2/F 'G= K4),;(''

;B/'--Rh#22\YP : X$ , /< 7&1 'h# /&5

&3&6;& :%& 7&5(') :8 2 6 2/ , 1 'OP&8

&/;& h#I&6$8(2'2)M/$:'/PF 6 &),

P'(S5C X$7&NF Z<) BC '; O) ) 4 $, $ H+ '$=P

h#22\YTUBôX_CBC UVW Wake

"#$%&'('

:

0"@EFGF$H0GIPaJ30G&#*K F$HIm/sJ;

0":EFG0G IPaJ3K /Im/sJL/( M/N$";

&#*+$,-./0'1#,-/CFD

O

Y=

Y=

Y=

^h#22_

JY=

@Y=

0"O&$'K &$-6Im/sJL/P,1( &$*$*;

"#$%&'('

!

p(oX x(J X

s(` X (f X

0"!G0&PAP,1<$H.1-G$*$*&;

\Y;/F/() '^O)^p,__7&/<

N# & &6 %& 7&5 2 8 : 2 #'BC lm$P$^.$ $ f@i 2f@5_

lm$P$22 3&6:=0 BBC 22 BCf@) /'# /) ,;& .$ $ 7&

7&5 2h#22:(8 &&6;& 6 (%& &,Z<)'f2`

27&O >;22 ) ,BCKG7& :' 7&5 2 8 =06 Y' B

7&528(6 ;& :vF2)'#)'#^O)2s,_

6$ 708

xU1B* ./ c 7C l'1B* ;& ; [&2;& B# R '

# &8(' ,1B*($ : OV0 2 ';2 2 m 2

xU1B*;6 ;& 3&6C /$c:"02$=PL/KD2D), =

3&xU1B*&$&./&2P2'' ($: 2 ') ' , ' l 1B*%& /&2'; &)

]& ) I=> N0 xU =) ),

&#*+$,-./0'1#,-/CFD

pF F<)OC2 5 xBCF<)OC2 5

0"PG A0GCDL-Q<R$1SRSMRNG k-;

2%& ;28# RF $'

k-RNG2RSM7&;& =)']& )&')H/= ,H'=)F

Z$ %& ;2 F<) BC 2 # Rf O)%& pe);B,;('

65U&/)BC2F '/B5>2)OC2 527&N#&&6$

'2,;('O)e);B 2 '(' F 2 BC L&$ 0& 2 H;+ 7 6R !"# 2 '&),

RNGk-7)\2N( '<O:CQε;& '/$DC& 'U,&

/&NF RSM/$7<$ O: /&'<'!"$~ =P 2Y()<3&

-&CQ=P;/'D;+2H'C'#8 ,7 6R!"# $ 7& /H' 2 1B* $=P RNGk-=/$=P/&6'' >#),

DE

$&I=> 20l #&/2\#<3&& X$/>

.$ $%& 222ppY````i`2`f`7&N

F r ) ;/ 4c/$) ,H'C2 5O' 0 2 + P $=>L\Y %& O)@ ) ;B ,;('

O)@&&6);B&6C;+8 $ %& p \#KG #

5/,8:7/S('&8OC &2OCP$/&5 lm$3&6& %& p 2pY

' ;+ ,C2 5O'kF C 0C2 5;B&O/',\YKG7&/B

0Y&/(kF 2)^;DI=> & ,_;(' U '7B5BC ) vP

%& #&&6/B;)3&62)& ( N# ' ' BC & ,7& /

&&67& B;# R)4'#2BC BC 6 =$ /& $ 0&BC =0

$N(F 7&5 0 X$ )' BC F<)0& X$ Y 2)9 'U 02 /R lm$

NF 7& 7& / '# 0 // 7<$ C2 5O/0';Y -&('=002',

=

L&/G%& ('O;c/6 -7&3/&2'&5:'

"#$%&'('

))

0")EFJ3+$,-<$H0GJ-C&T1L/+$,-F$H;

2/CCFD=)< KG)

,<'< O) ;& 4:' Y 2 / -5 %& 2 :

&MRF2:A' >727:7 l2)O:&h P),

Y&0H+''1B*RNG k-2RSM N( ' Z<) M BC 2 F &%

/R7H)RNG k-./$=PD 7&5 $=P /&6' (' ) ,H'

+C2 5KG$=PO2F &';B>L\YBC&'BC #'

F 2/BBC 2 F Z<) ' 2 ) & F<)KG& $;# R #N#

' 7& ,OC lm$ 7/S('2 BC H& ' F /B& # ;# R &

C D ,H'=)BC4 : H& tB5;B'/',F P'

$=POG:&)--R : X$h#22\YP,/<7&&3&6 ;& h# /&5 L>Y ( '

:%& 7&5826 2) :C7&kF%& 82/&),

F BGD

EU `kg/mρ

F m/sν

/$0/$Paτx)B

m/sg

BCPaP

3/&&0&2kg/m.sµ

=0F m/srν

bi

m/sgi

lb&0&2kg/m.sµeff

3/% 1B*J/s

Jmk

&#*+$,-./0'1#,-/CFD

)

/% !"$~ 1B*`/s

Jmε

< lbO 5F2kk , εα α

3/% $1B*`kg/m.sGk

/% $1B*\2Gb 2/)2 `

kg/m.s

-)+r0=1B*4 $ &95`kg/m.sYM

d1B*`kg/m.sDt,ij

d`kg/m.sDL,ij

/$$/$$`kg/m.sPij

$2/)2 `kg/m.sGij

!"$2N# R`kg/m.sεij

0'+BCB $`kg/m.sφij,1

L& BCB $`kg/m.sφij,2

&< $`kg/m.sφij,ω

;st

&2F $21/sω,Ω

&0/$I

& a

&$C&WJ`@@@`e@&$1 &95WJ@@@`eJ

+, [1] Golombok M., Bill K., Removal of CO2 from a Gas Stream Using An Experimental Centrifugr,

Ind. Eng. chem, 44, p. 4721 (2005).

[2] Van Wissen R., Golombok M., Brouwers J.J.H., Separation of Carbon Dioxide and Methane in

Continuous Countercurrent Gas Centrifuges, Chemical Engineering Science, 60, p. 4397

(2005).

[3] Cohen K., "The Theory of Isotope Separation as Applied to the LargeScale Production of

U 235", McGraw-Hill, New York, (1951).

[4] Cloutman, L.D.,Gentry R.A., "Numerical Simulation of the Countercurrent Flow in a Gas

Centrifuge", Technical Report, Los Alamos National Laboratory, (1983).

[5] Park J.K., Hyun J.M., Numerical Solutions for Thermally Driven Compressible Flows in a

Rapidly Rotating Cylinder, Fluid Dynamics Research, 6, p. 139 (1990).

[6] Harada I., A Numerical Study of Weakly Compressible Rotating Flows in a Gas Centrifuge,

Nucl. Sci. Eng., 73, p. 225 (1980).

[7] Ungarish M., Israeli M., Axisymmetric Compressible Now in a Rotating Cylinder with Axial

Convection, J. Fluid Mech., 154, p. 121 (1985).

[8] Kai T.,HasegawaK., Numerical Calculation of Flow and Isotop Separation for SF6 Gas

Centrifuge, Journal of Nuclear Science and Technology, 37, p. 153 (2000).

[9] Dongjun J.,Shi Z., CFD Simulation of 3D Flow Field in a Gas Centrifuge, 14th International

Conference on Nuclear Engineering, Florida, USA, (2006).

· a c# r ˇp ˚ ˇ ˚ o˚ ˇ> ... e˘u˚4 ˘0 7 k e˘u˚4 ˘0 ˙b# r. ˇ ˆ ˘ ˇ ˘ cfd 0...

Documents