- Email: [email protected]
Hydrodynamic dispersion of reactive solute in a Hagen–Poiseuille flow of a layered liquid
S1004-9541(16)30609-7 doi:10.1016/j.cjche.2017.03.005 CJCHE 761
To appear in: Received date: Revised date: Accepted date:
23 June 2016 1 March 2017 2 March 2017
Please cite this article as: Sudip Debnath, Apu Kumar Saha, B.S. Mazumder, Ashis Kumar Roy, Hydrodynamic dispersion of reactive solute in a Hagen-Poiseuille flow of a layered liquid, (2017), doi:10.1016/j.cjche.2017.03.005
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
ACCEPTED MANUSCRIPT !"#$%#"&'()* #)+,-$+)%& %. $-'*/)0- +%12/- )& ' !'3-&45%)+-2)116%7 %. ' 1'"-$-# 1)82)#
"#$!%&'#(& )* +!&,#'!&-./0 1!&-)(!2 3(/&-&4&# )* 5#.,()2)670 86!%&!2!0 9::;<=0 3(>-! ? @24A-!2 +#.,!(-./ B!?)%!&)%70 C,7/-./ !(> 8$$2-#> +!&,#'!&-./ D(-&0 3(>-!( E&!&-/&-.!2 3(/&-&4� F)2G!&!0 9;;H;I0 3(>-!
SC
RI
!
PT
92#), :-;&'/<='>> ?,2 @2('$ 9'<'' > AB9B C'D2(#-$; > ?+<)+ @2('$ E%"'
!"#$%&'$
TE
D
MA
NU
?& '&'1"+)+ %. /<- +%12/- #)+,-$+)%& )& /<- 1)82)# 6%7)&3 /<$%23< ' ,),- ;" (-'&+ %. ?$)+4A'$/%&F+ G(-/<%# %. (%(-&/+F> 2&#-$ /<- H%)&/ -..-*/ %. +%(- I&)/- ")-1# +/$-++ '&# )$$-0-$+);1- ';+%$,/)%& )&/% /<- 7'11 )+ ,$-+-&/-# )& /<)+ ,',-$B J<- 1)82)# )+ *%&+)#-$-# '+ ' /<$--41'"-$ 1)82)# 7<-$/<- *-&/-$ $-3)%& )+ K'++%& 1)82)# +2$$%2&#-# ;" L-7/%&)'& 1)82)# 1'"-$B ? +)3&)I*'&/ *<'&3.$%( ,$-0)%2+ (%#-11)&3 -M-$*)+-+ )& /<- +/2#" %. <"#$%#"&'()* #)+,-$+)%&> #)..-$-&/ (%1-*21'$ #)..2+)0)/" <'+ ;--& *%&+)#-$-# .%$ /<- #)..-$-&/ $-3)%& "-/ /% ;- *%&+/'&/B N%$ '11 /)(- ,-$)%#> I&)/#)..-$-&*- )(,1)*)/ +*<-(- <'+ ;--& '#%,/-# /% +%10- /<- )&/-3$'1 (%(-&/ -82'/)%& '$)+)&3 .$%( /<- 2&+/-'#" *%&0-*/)0- #)..2+)%& -82'/)%&B J<- ,2$,%+- %. /<- +/2#" )+ /% I&# /<- #-,-&#-&*" %. +%12/- /$'&+,%$/ *%-.I*)-&/+ %& ';+%$,/)%& ,'$'(-/-$> ")-1# +/$-++> 0)+*%+)/" $'/)%> ,-$),<-$'1 1'"-$ 0'$)'/)%& '&# )& '##)/)%& 7)/< 0'$)%2+ #)..2+)0)/" *%-.I*)-&/ )& #)..-$-&/ 1)82)# 1'"-$+B J<)+ O)&# %. +/2#" ('" ;- 2+-.21 .%$ 2&#-$+/'&#)&3 /<- #)+,-$+)%& ,$%*-++ )& /<- ;1%%# 6%7 '&'1"+)+B
AC CE P
F#7J)%>/K K'++%& 1)82)#> P)-1# +/$-++> ?M)'14#)+,-$+)%& *%-.I*)-&/> Q$$-0-$+);1- $-'*/)%&> :)..2+)0)/"B () *+$%,-.'$/,+
:)+,-$+)%& )+ /<- ,$%*-++ 7<)*< )+ 2/)1)D-# '+ '& -.I*)-&/ (-'&+ /% '**%(,1)+< /<- #)12/)%& %$ ()M)&3B Q& $-*-&/ #'"+> /<- +/2#" %& /<- 1%&3)/2#)&'1 #)+,-$+)%& <'+ 3')&-# '//-&/)%& #2- /% )/+ ,$'*/)*'1 )(,%$/'&*- )& /<- I-1# %. *<-()*'1 -&3)&--$)&3> ,<"+)%1%3)*'1 62)# #"&'()*+> -&0)$%&4 (-&/'1 +*)-&*-+> ;)%(-#)*'1 -&3)&--$)&3 -/*B .%$ ;%/< *%$,%$-'1 '&# -M/$'4*%$,%$-'1B 9,-*)I*'11" &2/$)-&/+> (-/';%1)* )/-(+> #$23+> -/*B 3-/ /$'&+,%$/-# /<$%23< ;1%%# 6%7> '+ ' *%&+-82-&*- %. #)..2+)0- '&# *%&0-*/)0- (-*<'&)+(+B J<- ,)%&--$ $-+-'$*<-$ )& /<)+ I-1# 7'+ J'"1%$ RS> TU '&# ' 0)/'1 $-+21/)&3 *%(()/(-&/ 7'+ /')1%$-# ;" ?$)+ RVU> 7<% .$'(-# '&# ,'$/)'11" +%10-# -82'/)%&+ .%$ /<- )&/-3$'1 (%(-&/+ %. /<- *1%2# %. *%&/'()&'&/B 9'&O'$'+2;$'('&)'& '&# W)11 RXU #-0-1%,-# ' #-$)0'/)0- -M,'&+)%& (-/<%# 0)DB G3-&-$'1)D-# #)+,-$+)%& (%#-1F> 7<)*< <%1#+ /$2- .%$ 1'$3- '+ 7-11 '+ +('11 /)(-+B ? &2(-$)*'1 +/2#" 7'+ ,-$.%$(-# ;" ?&'&/<'O$)+<&'& -/ '1B RYU /% +%10- /<*%&0-*/)0- #)..2+)%& -82'/)%& ;" ' +/'&#'$# I&)/- #)..-$-&*- (-/<%# 7<)*< 7-&/ %& /% +<%7 /<'/ /<- J'"1%$4?$)+ #)+,-$+)%& /<-%$" 7'+ ,-$/)&-&/ %&1" .%$ '#-82'/-1" -M,'&+)0- -+/)('/)%& %. /<#)(-&+)%&1-++ /)(-> )B-B> t > 0.8B A'$/%& RZU $-I&-# /<- ?$)+F+ (%#-1 '&# -M/-&#-# /<- 1-3)/)('*" /% *%0-$ -0-& +('11 '&# (%#-$'/- /)(-+B L'!-2 !>>%#//K [email protected] [92#), :-;&'/<= \
C%#$%-(& /4?'-&&#> &) M,-(#/# N)4%(!2 )* M,#'-.!2 L(6-(##%-(6
+!%., OO0 O;H9
ACCEPTED MANUSCRIPT
AC CE P
TE
D
MA
NU
SC
RI
PT
N%11%7)&3 #)..-$-&/ /-*<&)82-+> ' 1'$3- &2(;-$ %. '$/)*1-+ [R]> ^> _> S`> SS> STU\ <'0- ;--& .%2&# )& 1)/-$'/2$- #-'1)&3 7)/< 1%&3)/2#)&'1 #)+,-$+)%& )& +/-'#" '&# 2&+/-'#" 6%7+ 7)/<)& /
T
ACCEPTED MANUSCRIPT
NU
SC
RI
PT
.%2&# /% ;- #-*$-'+)&3 7)/< /<- )&*$-'+- )& /<- -M,'&+)%& %. /<- ,-$),<-$'1 1'"-$B Q& ' /7%462)#> ;1%%# 6%7 '&'1"+)+> 9'&O'$ RV]> V^U %;+-$0-# /<'/ /<- 0-1%*)/" #)+/$);2/)%& '&# 6%7 $'/- #-*$-'+;2/ /<- 7'11 +<-'$ +/$-++> 7)#/< %. /<- ,123 6%7 $-3)%& '&# 1%&3)/2#)&'1 )(,-#'&*- )&*$-'+- 7<-& /<- ")-1# +/$-++ )&*$-'+-+B f- <'0- +--& 0-$" .-7 7%$O+ %. #)+,-$+)%& *%&+)#-$)&3 ,-$),<-$'1 1'"-$ -..-*/+> +2*< '+ 9<2O1' '&# 5'$)<'$ RV_U> '&# 9<2O1' '&# W2,/' RX`UB J% /<- ;-+/ %. %2$ O&%71-#3-> /<-$- )+ <'$#1" '&" 7%$O )& /<- 1)/-$'/2$- *%&+)#-$)&3 /<- -..-*/ %. /<- ,-$),<-$'1 1'"-$ %& #)+,-$+)%& ,$%*-++> 7<-& *%$- $-3)%& )+ K'++%& 1)82)# '&# ' ,-$),<-$'1 1'"-$ %. ,1'+(' )+ L-7/%&)'& 1)82)#B C%$-%0-$> /<- ,$-+-&/ +/2#" ('&'3-+ /<- )$$-0-$+);1- ';+%$,/)%& '/ /<- ;%2&#'$" %. /<- ,),-B :2/% /<- .'*/ %. ,-$),<-$'1 1'"-$ ,$-+-&*-> /<- 0)+*%+)/" $'/)% ',,-'$+ .%$ /<- I$+/ /)(- )& /<- +/2#" %. <"#$%#"&'()* #)+,-$+)%&B J<- ,$-+-&/ (%#-1 ('" ;- 2+-.21 .%$ /<- ,<"+)%1%3)*'1 ;1%%# 6%7 '&'1"+)+ '+ /<- &'/2$- %. /<- -..-*/+ *%&+)#-$-# )& /<- +/2#"> /<- 7%$O )+ (2*< *1%+-$ /% /<'/ )& ' ,<"+)%1%3)*'1 +)/2'/)%& *%(,'$-# /% +-0-$'1 %/<-$ 7%$O+ )& 1)/-$'/2$-B 0) 1&$234&$/'&5 6,%4.5&$/,+
AC CE P
TE
r
D
MA
? 2&)4#)$-*/)%&'1> +/-'#" 1'()&'$> 'M)4+"((-/$)*> .211"4#-0-1%,-# 6%7 %. '& )&*%(,$-++);1-> /<$--41'"-$ 1)82)# /<$%23< ' *)$*21'$ ,),- %. $'#)2+ R )+ *%&+)#-$-#B Q/ )+ '++2(-# /<'/ /<- *%$$-3)%& )+ K'++%& 1)82)# -&*1%+-# ;" ' /<)& 1'"-$ %. L-7/%&)'& 1)82)#B
− ∂p ∂z
*+,-().&)
R Ro
%&''() !"#$ !"#$
Rp
Z
%&''() *+,-().&)
N)32$- Sg 9*<-('/)* #)'3$'( %. /<- +-/2, 2&#-$ *%&+)#-$'/)%&B
N)3BS +<%7+ /<- 6%7 3-%(-/$" 7)/< *"1)&#$)*'1 *%%$#)&'/- +"+/-( 7<-$- /<- 'M)'1 '&# $'#)'1 *%%$#)&'/-+ '$- $-,$-+-&/-# ;" z '&# r $-+,-*/)0-1" [/<- ;'$ #-&%/-+ #)(-&+)%&'1 82'&/)/)-+\B h& '**%2&/ %. -&%23< 'M)'1 *%0-$)&3 %. /<- /2;- '&# .%$ /<- 1%7 E-"&%1#+ &2(;-$ 6%7> )/ *'& ;+<%7& /<'/ /<- $'#)'1 0-1%*)/" )+ &-31)3);1" +('11 '&# *'& ;- #)+$-3'$#-#> /<2+ /<- .211" #-0-1%,-# 6%7 )+ %&1" )& /<- 'M)'1 #)$-*/)%&B A-*'2+- %. /<- 'M)4+"((-/$" '11 82'&/)/)-+ '$- )&#-,-&#-&/ %. θB
V
ACCEPTED MANUSCRIPT J<- 3%0-$&)&3 -82'/)%&+ %. (%/)%& .%$ /<- 'M)'1 6%7 '$- 3)0-& ;" 0 ≤ r ≤ Ro ,
[S\
−
∂p 1 ∂(r τ n ) = ∂z r ∂r
Ro ≤ r ≤ R.
[T\
PT
∂p 1 ∂(r τ c ) = ∂z r ∂r
−
NU
SC
RI
J<- +<-'$ +/$-++ )& /<- K'++%& $-3)%& Rc [= [0, Ro ]\ '&# L-7/%&)'& $-3)%& Rn [= [Ro , R]\ '$3)0-& ;" /<- *%&+/)/2/)0- $-1'/)%&+g 1 1 1 ∂u 2 c τ c = τ y 2 + (−µc ∂r ) 2 , τ c ≥ τ y .%$ Rp ≤ r ≤ Ro , ∂uc [V\ = 0, τ c ≤ τ y .%$ 0 ≤ r ≤ Rp , . ∂r ∂un τ n = −µn ∂r , τ y = 0 .%$ Ro ≤ r ≤ R,
TE
D
MA
f<-$- τ y )+ /<- ")-1# +/$-++B Q& i8B[V\ 7<-&-0-$ τ c ≤ τ y > /<- 0-1%*)/" 3$'#)-&/ 7)11 ;- D-$% '&# '+ ' *%&+-82-&*-> ,123 6%7 $-3)%& Rp [= [0, Rp ]\ ',,-'$+ 7)/<)& Rc B J<- ;%2&#'$"> +"((-/$" '&# ('/*<)&3 *%&#)/)%&+ .%$ +%10)&3 i8+B[S\ 4 [V\ '$-g c τ c )+ I&)/- '&# ∂u = 0 '/ r = 0, ∂r τ c = τ n '&# uc = un '/ r = Ro , . [X\ un = 0 '/ r = R,
AC CE P
e-/ 2+ &%7 *%&+)#-$ /<- /$'&+,%$/ '&# +,$-'#)&3 %. *<-()*'1 +,-*)-+ 7)/< /<- +%10-&/ 6%7)&3 /<$%23< /<- /2;- 7<)*< )+ +2,,%+-# /% ;- *%(,1-/-1" ()+*);1- 7)/< /<- 1)82)#B J<- +,-*)-+ '$'11-3-# /% *%&/')& )$$-0-$+);1- ';+%$,/)%& )&/% /<- 7'11B Q& /<- ,$-+-&/ +*-&'$)%> /<- /$'&+,%$/ -82'/)%& /<'/ 3%0-$&+ C(t, r, z) )+ /<- 2&+/-'#" *%&0-*/)0#)..2+)%& -82'/)%& '&# )+ 3)0-& ;" % & ∂C ∂C ∂ 2C D ∂ ∂C r +D 2, [Y\ + u(r) = ∂z r ∂r ∂r ∂z ∂t 7<-$- D )+ /<- (%1-*21'$ #)..2+)0)/" '&# 7- *%&+)#-$ D /% ;- #)..-$-&/> "-/ *%&+/'&/ )& /<- /<$-$-3)%&+ Rp > Rc '&# Rn B 9% 7- '++2( Dp )& Rp D= . [Z\ D )& Rc c Dn )& Rn
J<- )&)/)'1 '&# ;%2&#'$" *%&#)/)%&+ .%$ /<- /$'&+,%$/ i8B[Y\ '$- *%&+)#-$-# '+ C(0, r, z) = C 0 B(r)ψ(z), ∂C =0 ∂r
0
'/ r = 0,
∂C +β C =0 ∂r X
'/ r = R,
[]\ [^\ [_\
ACCEPTED MANUSCRIPT
PT
7<-$- B(r) '&# ψ(z) '$- &--#-# /% ;- +,-*)I-# /% +%10- i8B[Y\> C 0 )+ /<- )&)/)'1 *%&*-&/$'/)%& %. /<- +123 )&,2/B J<- ';+%$;)&3 ;%2&#'$" *%&#)/)%& '/ /<- 7'11 %. /<- /2;- )+ $-,$-+-&/-# ;" i8B[_\B b+)&3 /<- .%11%7)&3 #)(-&+)%&1-++ 82'&/)/)-+g 2 R dp , uc = uu0c , un = uun0 , C = CC , t = Dn2 t , r = Rr , z = D2n z , u0 = − 4µ dz 0 n R R u0 Rp τy τc τn Ro [S`\ τc = µ ( u0 ) , τn = µ ( u0 ) , R0 = R , Rp = R , τy = µ ( u0 ) , β = β R, n R n R n R P e = uD0nR , m = µµc , D∗ = DDn ,
RI
n
% & ∂C D∗ ∂ 2 C r + , ∂r P e2 ∂z 2
[SS\
NU
∂C ∂C D∗ ∂ + u(r) = ∂t ∂z r ∂r
SC
/<- i8+B [Y 4 _\> ('" ;- 7$)//-& '+g
MA
Dp )& Rp Dc )& Rc , D∗ = 1 )& Rn
[ST\
'1%&3 7)/< )&)/)'1 '&# ;%2&#'$" *%&#)/)%&+ '+g
(0 < r < 1)
[SV\
D
C(0, r, z) = B(r)ψ(z),
TE
∂C =0 ∂r
'/ r = 0,
[SX\
AC CE P
∂C = −βC '/ r = 1. [SY\ ∂r !-$- u0 )+ /<- /)(-4'0-$'3-# 'M)'1 0-1%*)/"> β )+ /<- ';+%$,/)%& ,'$'(-/-$ %$ I$+/ %$#-$ $-'*/)%& $'/- $-,$-+-&/)&3 /<- $'/- %. 1%++ %& /<- 7'11B*P e+)+ /<- 5-*1-/ &2(;-$ 7<)*< * (-'+2$-+ /<- $-1'/)0+ R2 R *<'$'*/-$)+/)* /)(- %. /<- #)..2+)%& ,$%*-++ D /% /<- *%&0-*/)%& ,$%*-++ u0 B !-$- m )+ /<0)+*%+)/" $'/)% 7<)*< )&#)*'/-+ /<- )&62-&*- %. K'++%& 0)+*%+)/" [µc \ *%(,'$-# /% /<- L-7/%&)'& 0)+*%+)/" [µn \B J<- $'#)2+ %. /<- *"1)&#-$ *%&/')&)&3 +%12/- /% /<'/ %. /<- -&/)$- /2;- )+ '++2(-# '+ 0'12- %&-> +2;+-82-&/1"> /<- +123 )&)/)'11" %**2,)-+ /<- -&/)$- *$%++4+-*/)%& %. /<- /2;-B f)/< /<- '++)+/'&*- %. ';%0- ;%2&#'$" *%&#)/)%&+ [X\ )& #)(-&+)%&1-++ .%$(> /<- 0-1%*)/" #)+4 /$);2/)%&+ '$- '*82)$-# '+g 0 . / * 3 3+ , , 2τ 4 1 y 0 ≤ r ≤ Rp , Ro2 − Rp2 + τy (Ro − Rp ) up (r) = 1 − Ro2 + Ro2 − Rp2 − m 3 - 1 , uc (r) = 1 − Ro2 + m
.
- 4 , 2 Ro − r 2 −
/
+ 2τy * 23 3 2 Ro − r + τy (Ro − r) 3
, un (r) = 1 − r2
0
Rp ≤ r ≤ Ro , Ro ≤ r ≤ 1.
[SZ\
f<-$- uc '&# un '$- $-+,-*/)0-1" /<- 0-1%*)/)-+ .%$ +<-'$ 6%7 )& K'++%& '&# L-7/%&)'& 1)82)#B J<*%&+/'&/ 0-1%*)/" up ',,-'$+ %&1" )& /<- ,123 6%7 $-3)%& <'0)&3 $'#)2+ Rp = τ2y B ?1+% Ro (≤ 1) )+ Y
ACCEPTED MANUSCRIPT
PT
/<- $'/)% %. /<- *-&/$'1 *%$- $'#)2+ '&# &%$('1 ,),- $'#)2+B f- &%7 *%&+)#-$ /<- .%11%7)&3 /7% '114/)(- ',,$%'*<-+ )& +/2#")&3 #)+,-$+)%& %. /<- $-'*/)0+%12/-g [)\ ?$)+4A'$/%& ',,$%'*< '&# [))\ 9'&O'$'+2;$'('&)'&4W)11 ',,$%'*
RI
7) !%/#89&%$,+ !::%,&'2
NU
−∞
SC
J<- pth (%(-&/ %. /<- #)+/$);2/)%& %. /<- +%12/- )& /<- #)$-*/)%& %. 6%7 *'& ;- #-+*$);-# ;" ?$)+4A'$/%& RZU '+> 1 +∞ (p) z p C(t, r, z)dz. [S]\ C (t, r) =
MA
9% 2+)&3 i8+B[S]\> /<- #)..2+)%& i8B[SS\ +2;H-*/ /% /<- )&)/)'1 '&# ;%2&#'$" *%&#)/)%&+ [SV 4 SY\ *'& ;- 7$)//-& '+ % & D∗ ∂C (p) ∂C (p) D∗ ∂ r = p u(r) C (p−1) + p(p − 1)C (p−2) , [S^\ − ∂t r ∂r ∂r P e2 7)/<
(0, r) =
D
C
(p)
2
1 Pe
.%$ p = 0 .%$ p > 0
AC CE P
TE
0 ∂C (p) . = 0 '/ r = 0, ∂r(p) ∂C = −βC (p) '/ r = 1, ∂r J<- *$%++4+-*/)%&'1 '0-$'3- %. /<- pth (%(-&/ %. /<- #)+/$);2/)%& %. +%12/- )+ 3)0-& ;"g 1 1 3 (p) 4 C (t) = 2 rC (p) (t, r)dr,
[S_\
[T`\
0
7<-$- /<- '&31- ;$'*O-/ #-&%/-+ /<- *$%++4+-*/)%&'1 (-'&B f)/< /<)+ #-I&)/)%&> i8B[S^\ '&# /<- *%&#)/)%&+ [S_\ ;-*%(->
4 D∗ 3 3 (p−2) 4 # 3 (p) 4 , C + 2βD∗ C (p) (t, 1) = p u(r)C (p−1) + p(p − 1) C #t P e2
[TS\
7)/< /<- )&)/)'1 *%&#)/)%&
3
C (p) (0)
4
1 Pe = 0
=
.%$ p = 0 .%$ p > 0
[TT\
J<- pth %$#-$ *-&/$'1 (%(-&/ ';%2/ (-'& %. /<- *%&*-&/$'/)%& #)+/$);2/)%& *'& ;- #-I&-# '+ 5 1 5 2π 5 +∞ r(z − zg )p C#r#θ#z 0 0 −∞ µp (t) = , [TV\ 5 1 5 2π 5 +∞ rC#r#θ#z 0 0 −∞ Z
ACCEPTED MANUSCRIPT 7<-$-
3 (1) 4 555 C zC#v = , zg = 5 5 5 hC (0) i C#v
NU
SC
RI
PT
4 3 $-,$-+-&/+ /<- G*-&/$%)#F %$ GI$+/ (%(-&/F %. /<- #)+/$);2/)%& %. +%12/- '&# C (0) )+ /<- /%/'1 ('++ %. /<- *<-()*'1 +,-*)-+ )& /<- 6%7)&3 +/$-'(B J<- -M,$-++)%&+ .%$ *-&/$'1 (%(-&/+ *'& ;- %;/')&-# .$%( i8B[TV\ '+> hC (2) i 2 µ2 (t) = C (0) − zg , h i (3) hC i 3 µ3 (t) = C (0) − 3zg µ2 − zg , . [TX\ h i (4) hC i 2 4 µ4 (t) = C (0) − 4zg µ3 − 6zg µ2 − zg , h i
AC CE P
TE
D
MA
Q/ 7'+ ?$)+ RVU 7<% .%2&# /<- ,<"+)*'1 -++-&*- %. -'*< %. /<- )&/-3$'1 (%(-&/+ %. *%&*-&/$'/)%& #-I&-# ;" i8B[TX\ 7<)1- -M,1')&)&3 /<- #)+,-$+)%& ,<-&%(-&' %& 1)82)#B N%$ )&+/'&*-> /<- D-$%/< (%(-&/ 3)0-+ /<- '33$-3'/- $'&3- 2&#-$ /<- #)+/$);2/)%& *2$0- 7<)*< *%$$-+,%&#+ /% /<- /%/'1 ('++ %. /<- +%12/-> /<- I$+/ (%(-&/ (-'+2$-+ /<- 1%*'/)%& %. /<- *-&/-$ %. 3$'0)/" %. /<- +123 (%0)&3 7)/< /<- (-'& 0-1%*)/" %. /<- 1)82)# '&# +-*%&# *-&/$'1 (%(-&/ $-,$-+-&/+ /<- 0'$)'&*- %. /<#)+/$);2/)%& ';%2/ /<- (-'& ,%+)/)%& 7<%+- $'/- %. *<'&3- 3)0-+ /<- #)+,-$+)%& *%-.I*)-&/B J<- $-'*/)%& *%-.I*)-&/ %$ -M*<'&3- *%-.I*)-&/ M0 *'& ;- I32$-# .$%( /<- 7'" /<'/ /<- /%/'1 ('++ #-*'"+ -M,%&-&/)'11" 7)/< /)(- '**%$#)&3 /%g %1 t & 3 0 4 M0 (t) = C (0) Exp M0 dt . [TY\ 0
?3')& /<- *-&/-$ %. ('++ (%0-+ '/ ' +,--# -82)0'1-&/ /% /<- '#0-*/)%& +,--#> /<- '#0-*/)%& *%-.I4 *)-&/ M1 )+ 3)0-& ;"> dzg M1 = . [TZ\ dt N%11%7)&3 ?$)+ RVU /<- ',,'$-&/ #)+,-$+)%& *%-.I*)-&/> Da > *'& ;- -M,$-++-# '+ .%11%7+g Da (β, τy , Ro , m, P e, t) =
1 dµ2 . 2 dt
[T]\
3
J<- *%-.I*)-&/ %. +O-7&-++ ν2 (= µ3 /µ22 ) '&# O2$/%+)+ ν3 (= µ4 /µ22 − 3) '$- '1+% )(,-$'/)0- -14 -(-&/+ .%$ (-'+2$)&3 /<- #-3$-- %. +"((-/$" '&# ,-'O-#&-++ %. /<- *%&*-&/$'/)%& #)+/$);2/)%& $-+,-*/)0-1"B !'0)&3 #)+*2++-# )& #-/')1 /<- ?$)+4A'$/%& ',,$%'*< /% #)+,-$+)%&> 7- &%7 (%0- %& /% /<- #-/')1+ %. /<- 9'&O'$'+2;$'('&)'&4W)11 ',,$%'*/55 &::%,&'2 N%11%7)&3 /<- (%#-1 #-0-1%,-# ;" 9'&O'$'+2;$'('&)'& '&# W)11 RXU> /<- *%&*-&/$'/)%& -M,'&#+ '+
]
ACCEPTED MANUSCRIPT
C(t, r, z) =
∞ 6
fi (r, t)
i=0
fk′ s
∂ i Cm , ∂z i
[T^\
RI
PT
7<-$'$- /<- /)(-4#-,-&#-&/ *%-.I*)-&/+ /% ;- #-/-$()&-# '&# /<- #)(-&+)%&1-++ '$-'4'0-$'3*%&*-&/$'/)%& Cm )+ 3)0-& ;" 5 2π 5 1 1 1 dθ 0 rCdr 0 Cm = Cm (t, z) = 5 2π 5 1 Crdr, [T_\ =2 dθ 0 rdr 0 0
SC
2+)&3 i8B[T_\ )& i8B[SS\ '1%&3+)#- i8B[T^\ '&# '1+% /<- ;%2&#'$" *%&#)/)%&+ [SV 4 SY\> /<- .%11%74 )&3 #)+,-$+)%& (%#-1 .%$ Cm 3)0-+ ∞
NU
∂Cm 6 ∂ i Cm = Mi (t) . ∂t ∂z i i=0
[V`\
MA
?**%$#)&3 /% 9'&O'$'+2;$'('&)'& '&# W)11 RXU> )/ 7)11 ;- +2.I*)-&/ /% /$2&*'/- /<- )&I&)/- +2( )& i8B[V`\ './-$ )&0%10)&3 M2 (t)B f<-$->
D
M0 (t) = −2βD∗ f0 (1, t), ∗
TE
M1 (t) = −2βD f1 (1, t) − 2
1
1
u(r) rf0 (r, t)dr,
[VS;\
0
D∗ − 2βD∗ f2 (1, t) − 2 M2 (t) = 2 Pe
AC CE P
[VS'\
1
1
u(r) r f1 (r, t)dr.
[VS*\
0
Q& 1)3 M0 (t) '$)+-+> 7<)*< *%$$-4 +,%&#+ /% /<- ';+%$,/)%& ,'$'(-/-$B ?+ /<- +%12/- 2&#-$3%-+ 7)/< '& )$$-0-$+);1- $-'*/)%& '/ /<%2/-$ ;%2&#'$"> )/ 7)11 #-,1-/- 7)/< /<- /)(- /<'/ ('O-+ &-3'/)0- -M*<'&3- *%-.I*)-&/B Q& +2*< ' *'+- -M*<'&3- *%-.I*)-&/ 7%21# ;- ,%+)/)0- '+ β )& /<- i8B[SY\ 7%21# /'O- %& ' &-3'/)0- +)3&B M1 (t) '&# M2 (t) $-1'/- /% /<- *%&0-*/)0- '&# #)+,-$+)%& *%-.I*)-&/+ $-+,-*/)0-1"B Q/ )+ &-*-++'$" /% (-&/)%& /<'/ M2 (t) )+ '&'1%3%2+ /% /<- #)+,-$+)%& *%-.I*)-&/ [Da \ %;/')&-# .$%( ?$)+4A'$/%& (-/<%# %. (%(-&/+B N%$ +%10)&3 /<- ';%0- *%-.I*)-&/+ %&- &--#+ /% I&# /<- *%$$-+,%&#)&3 .2&*/)%&+ fk F+B 9%> 7- 2/)1)Dk i8B[T]\ )& i8B[SS\ '1%&3 7)/< i8B[V`\ '&# -82'/)&3 /<- *%-.I*)-&/ %. ∂∂zCkm > /<- .%11%7)&3 ,'$/)'1 #)..-$-&/)'1 -82'/)%& %;/')&-# '+ k
D∗ ∂ ∂fk 6 Mi (t)fk−i (r, t) + u(r)fk−1 − + ∂t r ∂r i=0
% & ∂fk D∗ r − fk−2 = 0, ∂r P e2
7<-$- f−1 = f−2 = 0 '&# k = 0, 1, 2. N$%( i8+B[SV 4 SY\> /<- )&)/)'1 *%&#)/)%&+ .%$ Cm '&# fk '$1 1 B(r) rdr, Cm (z, 0) = 2ψ(z) 0
^
[VT\
[VV\
ACCEPTED MANUSCRIPT
fk (r, 0) =
2
1 .%$ k = 0 0 .%$ k = 1, 2
7
,
[VX\
'&# /<- ;%2&#'$" *%&#)/)%&+ .%$ /<- ,),- '$-
'1+% Cm (t, ∞) = '&#
1 0
.%$ k = 0, 1, 2
[VZ\
∂cm (t, ∞) = 0, ∂z
[V]\
.%$ k = 0, 1, 2.
[V^\
1 fk (t, r) rdr = δk0 2
MA
1
PT
at r = 1,
[VY\
SC
∂fk = −βfk (r, t) ∂r
.%$ k = 0, 1, 2
RI
at r = 0,
NU
∂fk =0 ∂r
?) @/+/$3 -/663%3+'3 ,6 $23 .+#$3&-A ',+B3'$/B3 -/66.#/,+ 3C.&$/,+ /+ $23 $2%3385&A3% %3D/,+ $, ,"$&/+ $%&+#:,%$ ',36E'/3+$#
AC CE P
TE
D
:2- /% /<- ,$-+-&*- %. *%(,1-M)/" )& /<- '&'1"/)*'1 +%12/)%&> )& /<)+ +-*/)%&> 7- <'0- &2(-$4 )*'11" *%(,2/-# /<- +%12/)%& %. /<- (%(-&/ i8B[S^\ 7)/< )&)/)'1 '&# ;%2&#'$" *%&#)/)%&+ [S_\B N%$ '11 /)(- ,-$)%#+> 7<-& /<- +%12/- )+ )&H-*/-# 2&).%$(1" %0-$ /<- *$%++4+-*/)%& %. /<- ,),-> 7#-0-1%,-# ' I&)/- #)..-$-&*- .%$(21' )& 0)-7 %. K$'&O4L)*<%1+%& )(,1)*)/ +*<-(-> 7<-$- /<- -&4 /)$- 7)#/< %. /<- ,),- <'+ ;--& #)0)#-# )&/% (M − 1) -82'1 ,'$/+ <'0)&3 1-&3/< ∆rB J<- 1-&3/<+ '&# /)(-+ '$- $-,$-+-&/-# ;" /<- 3$)# ,%)&/ j '&# i $-+,-*/)0-1"> +% /<'/ /<- 3-&-$'1 .%$(21'+ '$rj = (j − 1) × ∆r '&# ti = ∆t × (i − 1)B Q& /<)+ 7'"> .%$ j = 1, M *%$$-+,%&#+ /% /<- 'M)+ %. /<- ,),- [r = 0\ '&# /<- 7'11 [r = 1\ %. /<- ,),- $-+,-*/)0-1"> '1+% '+ i = 1 *%$$-+,%&#+ /% /<)&)/)'1 /)(- [t = 0\> 7<-$- ∆t '&# ∆r '$- /<- '23(-&/'/)%&+ %. t '&# r $-+,-*/)0-1"B C (p) (i, j) )&#)*'/-+ /<- 0'12- %. C (p) '/ /<- i4/< 3$)# ,%)&/ '1%&3 /<- t4'M)+ '&# j4/< 3$)# ,%)&/ '1%&3 r4'M)+B J<- $-+21/)&3 I&)/- #)..-$-&*- .%$(21' 1-'#+ /% ' +"+/-( %. '13-;$')* -82'/)%& 7)/< ' /$)4#)'3%&'1 *%-.I*)-&/ ('/$)M> Pj C (p) (i + 1, j + 1) + Qj C (p) (i + 1, j) + Rj C (p) (i + 1, j − 1) = Sj .
[V_\
f<-$- Pj > Qj > Rj '&# Sj '$- /<- ('/$)M -1-(-&/+B J<- I&)/- #)..-$-&*- .%$( %. /<- )&)/)'1 *%&#)/)%& )+> C (p) (1, j) = P1e , .%$ p = 0 = 0, .%$ p > 0 '&# /<'/ %. ;%2&#'$" *%&#)/)%&+ '$C (p) (i + 1, 0) = C (p) (i + 1, 2), '/ /<- 'M)+ %. ,),-> '&# C (p) (i + 1, M + 1) = C (p) (i + 1, M − 1) − 2∆rβC (p) (i + 1, M ), '/ /<- +2$.'*- %. /<- ,),- 7'11B ? C?Je?A *%#- <'+ ;--& #-0-1%,-# /% +%10- /<- /$)4#)'3%&'1 *%-.I*)-&/ ('/$)M $-&#-$-# .$%( i8B[V_\ ;" /<- (-/<%# %. J<%('+ '13%$)/<( [?&#-$+%& #& !2P RXSU\ 7)/< /<- <-1, %. ,$-+*$);-# )&)/)'1 '&# ;%2&#'$" *%&#)/)%&+B J<- *%(,2/'/)%& +/-,+ '$-g _
ACCEPTED MANUSCRIPT
MA
NU
SC
RI
PT
<$3:8(F J<- /)(- )&#-,-&#-&/ 'M)'1 0-1%*)/" u )+ *%(,2/-# I$+/ .$%( i8B[SZ\B <$3:80F J<- *%&*-&/$'/)%& C (p) )+ /<-& *'1*21'/-# .$%( i8B[S^\ '+ /<- 0'12- %. u(r) '/ -'*< %. /<3$)# ,%)&/+ '$- O&%7& .$%( <$3:8(B 4 3 <$3:87F N)&'11" /<- 0'12- %. C (p) )+ *'1*21'/-# .$%( i8B[T`\ ;" ',,1")&3 9)(,+%&F+ %&-4/<)$# $21-> 2+)&3 /<- 0'12-+ .%2&# .$%( <$3:8( '&# <$3:80 '/ /<- *%$$-+,%&#)&3 3$)# ,%)&/+B L2(-$)*'1 *'1*21'/)%&+ <'0- ;--& ,-$.%$(-# .%$ 0'$)%2+ 0'12-+ .%$ ';+%$,/)%& ,'$'(-/-$> ")-1# +/$-++> ,-$),<-$'1 1'"-$ /<)*O&-++> 0)+*%+)/" $'/)%> '&# #)+/)&*/ (%1-*21'$ #)..2+)0)/)-+ '/ #)..-$-&/ 1'"-$+ /% #)..-$-&/)'/- /<- )&#)0)#2'1 #)+,-$+)%& ,$%*-++-+B ?1/<%23< /<- ,$-+-&/ +*<-(- )+ 1)&-'$1" +/';1- .%$ '&" I&)/- 0'12-+ %. ∆t/(∆r)2 > 7<-$- ' (-+< +)D- (∆t = 0.00001, ∆r = 1/(M − 1)) '&# M = 100> 3)0-+ +'/)+.'*/%$" $-+21/+ .%$ #)..-$-&/ 0'12-+ %. ,'$'(-/-$+B Q& '11 /<- *'+-+> 7∂p | = 2 '&# P e = 103 B N$%( /<- '.%$-(-&/)%&-# +,'/)'1 '&# /-(,%$'1 #)+*$-/)D'/)%& <'0- /'O-& | ∂z ,'$'(-/-$+> 7- <'0- 32'$'&/--# %2$+-10-+ %. ' #-*-&/ $-82-+/ %. ,$-*)+)%& %. /<- $-+21/+ 7<)10'$")&3 /<- /)(- +/-, '&# 3$)# +,'*)&3+B ?*/2'11"> /<- +('11-$ /)(- )&/-$0'1 7)/< +2.I*)-&/1" +('11 +,'/)'1 #)+*$-/)D'/)%&+ <'+ ;--& 2/)1)D-# )& %$#-$ /% O--, /<- '**2$'*" %. /<- $-+21/+B !'0)&3 &2(-$)*'11" %;/')&-# /<- /$'&+,%$/ *%-.I*)-&/+ .$%( /<- ';%0- /7% '114/)(- ',,$%'*<-+> 7- &%7 %2/1)&- /<- ,$%*-#2$- %. %;/')&)&3 /<- (-'& +%12/- *%&*-&/$'/)%& #)+/$);2/)%& )& /<- 1'"-$-# 1)82)# (-#)2(B G) H&5'.5&$/,+ ,6 $23 +,+8>&.##/&+ 43&+ #,5.$3 ',+'3+$%&$/,+ .#/+D I3%4/$3 :,5A+,4/&5#
AC CE P
TE
D
J<- ;-<'0)%2$ %. *%&*-&/$'/)%& #)+/$);2/)%& ()3 )/ )+ ,%++);1- /% ',4 ,$%M)('/- /<- (-'& 'M)'1 *%&*-&/$'/)%& #)+/$);2/)%& Cm (z, t) %. /$'*-$+ 7)/<)& /<- 6%7 $-3)%& 7)/< /<- '++)+/'&*- %. !-$()/- ,%1"&%()'1 $-,$-+-&/'/)%& .%$ &%&4W'2++)'& *2$0-+ [C- ∞ 3 (0) 4 −η2 6 Cm (z, t) = C (t) e an (t)Hn (η), [X`\ n=0
√ hC (1) i 7<-$- η = (z − zg )/ 2µ2 , zg = C (0) '&# Hi > /<- !-$()/- ,%1"&%()'1+> +'/)+." /<- $-*2$$-&*h i $-1'/)%& 7)/< H0 (η) = 1.0 '+ Hi+1 (η) = 2ηHi (η) − 2iHi−1 (η), i = 0, 1, 2, . . .
J<- *%-.I*)-&/+ ai F+ '$a0 = 1/ (2πµ2 )1/2 , a1 = a2 = 0, a3 = 21/2 a0 ν2 /24, a4 = a0 ν3 /96 J<-$-.%$-> 3)0-& /<- +/'/)+/)*'1 ,'$'(-/-$+ )& i8B[TX\> /<- *%&*-&/$'/)%& #)+/$);2/)%& *'& ;- -+/)4 ('/-# .$%( i8B[X`\ '/ '&" 3)0-& 1%*'/)%& )& /<- 'M)'1 #)$-*/)%& '&# /)(-B J) K3#.5$# &+- L/#'.##/,+ J<- ,$-+-&/ '&'1"+)+ -++-&/)'11" #-'1+ 7)/< /<- +/$-'(47)+- #)+,-$+)%& %. +%12/- )& +%10-&/ 6%74 )&3 /<$%23< ' &'$$%7 ,),- 2&#-$ /<- -..-*/ %. ,-$),<-$'1 1'"-$ /<)*O&-++ 7)/< )$$-0-$+);1- $-'*/)%& %**2$ '/ /<- ;%2&#'$"B Q&1)&- 7)/< /<- #)+*2++)%&+ %2/1')# )& 9-*/)%&BS> /<- )(,%$/'&*- %. /<- -..-*/ %. ,-$),<-$'1 1'"-$ )& ;1%%# 6%7> /<- 0'12-+ %. ")-1# +/$-++ '&# ,-$),<-$'1 1'"-$ /<)*O&-++ '++2(' ,)0%/'1 +)3&)I*'&*-B ?+ )&#)*'/-# ;" C*:%&'1# RXVU> ;1%%# ,%++-++-+ '& ")-1# +/$-++ %. `B`X S`
ACCEPTED MANUSCRIPT
1
/06
/06
/04
/04
τ;<=/0/ /0/2 /0/4 /0/6
/06 /04
:/02
/
:/02
:/04
:/06
:/06
:/06
:/08
:/08
:/08
/01
/02
/03
/04
/05
u
/06
/07
/08
/09
1
:1
/
/01
/02
/03
AC CE P
/
TE
:/04
:1
/04
u
A<10/ 101 102 103
/02
D
:/02
/
>=B=?
/08
MA
/
!"<=/09/ /094 /098 10//
/02 r
/02 r
/08
NU
/08
1
>=@=?
r
>=&=?
SC
1
RI
PT
dyn/cm2 '/ <-('/%*$)/ 0'12- %. X`jB Q& *%&/$'+/ /% ,2;1)+<-# +/2#)-+ )& /<)+ #)$-*/)%&> /<- ,$-+-&/ +/2#" <'+ )&0-+/)3'/-# /<- ,'$'(-/$)* -..-*/ %. (%1-*21'$ #)..2+)0)/" '/ #)..-$-&/ 1)82)# 1'"-$+ 7<)*< <'+ <-1,-# /% ,$%H-*/ /<- )(,%$/'&*- %. 0)+*%+)/" $'/)% )& /<- ,$-+-&/ ,'$'#)3( %. +/2#"B Q& %$#-$ /% -+/';1)+< /<- #-+)$-# ,'$'(-/$)* -..-*/ ;%/< 82'1)/'/)0-1" '&# 82'&/)/'/)0-1"> 0'$)';1-+ 7-$- )&4 0-+/)3'/-# )& /<- .%11%7)&3 $'&3-+ g [)\ )$$-0-$+);1- ';+%$,/)%& ,'$'(-/-$ [β\ .$%( ` /% S``> [))\ ")-1# +/$-++ [τy \ .$%( ` /% `BT> [)))\ 0)+*%+)/" $'/)% [m\ $'&3-+ .$%( S /% SBT '&# [)0\ ,-$),<-$'1 1'"-$ 0'$)'/)%& [Ro \ 0'$)-+ .$%( `B^Z /% SB N)3BT #-(%&+/$'/-+ /<- 0-1%*)/" ,$%I1-+ 7)/< $-+,-*/ /% /<,),- $'#)2+ '/ #)..-$-&/ ,'$'(-/$)* 0'12-+ %. τy > Ro '&# mB
:/04
/05
/06
/07
/08
:1
/
/01
/02
/03
/04
u
/05
/06
/07
/08
N)32$- Tg 9/-'#" 6%7 0-1%*)/" ,$%I1-+ .%$ #)..-$-&/ 0'12-+ ")-1# +/$-++ [τy \> ,-$),<-$'1 1'"-$ 0'$)'/)%& [Ro \ '&# 0)+*%+4 )/" $'/)% [m\B
J% 0-$)." /<- '**2$'*" %. /<- &2(-$)*'1 +*<-(-> /<- .%11%7)&3 *'+-+ <'0- ;--& *%&+)#-$-#g [)\ Q& %$#-$ /% *%(,'$- /<- $-+21/+ %. #)+,-$+)%& %;/')&-# )& /<- ,$-+-&/ +/2#" 7)/< /<'/ -+/';4 1)+<-# )& /<- +/2#" %. 9'&O'$'+2;$'('&)'& '&# W)11 RXU /<- ,'$'(-/-$+ %. )&/-$,<'+- ('++ /$'&+4 .-$ 2&#-$ *%&+/'&/ ,$-++2$- 3$'#)-&/ <'+ ;--& *<'$'*/-$)D-# 7)/< /<- +'(- 0'12-+ %. Ro = 1> Dp = Dc = 1> m = 1 '&# Rp = 0.0, [.%$ L-7/%&)'& 62)#\B J<- ,$-+-&/ (%#-1 <'+ ;--& .2$/<-$ )&0-+/)3'/-# 7)/< 0'$")&3 Rp = 0.02, 0.04, 0.06, 0.08, 0.1, 0.2 7<)*< *%$$-+,%&#+ /% /<- #)+4 ,-$+)%& +/2#)-+ %. L'3'$'&) -/ '1B RXXU .%$ *%(,'$)+%& %. /<- (%#-1 0)';)1)/" )& /<- *'+- %. K'++%& 62)#B N)3BV $-,$-+-&/+ /<- ,1%/ %. #)+,-$+)%& *%-.I*)-&/ [M2 \ 0+ ';+%$,/)%& [β\B ?+ )+ -0)#-&/ .$%( /<- $-+21/+ %;/')&-#> /<- (%#-1 /$-&#+ *%(,1-/-1" *%&*2$ 82'1)/'/)0-1" 7)/< /<'/ %;/')&-# )& /<+/2#)-+ %. 9'&O'$'+2;$'('&)'& '&# W)11 RXU '&# 7)/< L'3'$'&) -/ '1B RXXUB !%7-0-$ /<- 82'&/)/'4 /)0- #)..-$-&*-+ $-'1)D-# )& /<- ,$-+-&/ +/2#" )+ '//$);2/-# /% #)..-$)&3 &%$('1)D'/)%& +*'1-+ '#%,/-# )& /<- ,$-+-&/ +/2#"B [))\ Q& %$#-$ /% 0'1)#'/- /<- ,$-+-&/ (%#-1 $-+,%&+- 7)/< /<%+- -+/';1)+<-# )& /<- +/2#" %. :'+< -/ '1B RTZU> /<- +)(21'/)%& ,'$'(-/-$+ %. Ro , Dp , Dc '&# m ;-'$)&3 /<-)$ )&#)0)#2'1 +)3&)I*'&*- <'# ;--& O-,/ '/ 2&)/"B J% /<)+ -&#> /<- (%#-1 $-+,%&+-+ %2/1')# )& N)3BX #)+,1'"+ *%(,1-/- '3$--(-&/ 7)/< /<%+- %;/')&-# )& /<- +/2#" %. +%12/- #)+,-$+)%& %. :'+< -/ '1B RTZU 7)/<%2/ *'/'1"/)* 7'11 $-'*/)%&> )B-B β = 0 [+-- N)3BX /%/'11" *%)&*)#-+ 7)/< N)3BT['\ )& RTZU\B
SS
ACCEPTED MANUSCRIPT :3
6
C=1/
Rp = 0.0 0.02 0.04 0.06 0.08 0.10 0.20
5
PT
3
RI
M2
4
SC
2
/ :2 1/
NU
1
:1
/
1/
1
1/ β
2
1/
1/
MA
N)32$- Vg k'$)'/)%& %. #)+,-$+)%& *%-.I*)-&/ [M2 \ 7)/< ';+%$,/)%& ,'$'(-/-$ [β\ .%$ 0'$)%2+ 0'12-+ %. ,123 $'#)2+ [Rp \ .%$ /<- *'+- 7<-& Ro = 1, m = 1 '&# Dp = Dc = 1B
D
1
TE
Rp = 0.0 0.02 0.04 0.06 0.08 0.10 0.15 0.20
/08 /07
AC CE P
! " 192 × M2 − 1/P e2
/09
/06 /05 /04 /03 /02 /01 /
/
/0/5
/01
/015
/02
/025
t
/03
/035
/04
/045
/05
, N)32$- Xg k'$)'/)%& %. #)+,-$+)%& *%-.I*)-&/ R192 × M2 − 1/P e2 U 7)/< /)(- [t\ .%$ 0'$)%2+ 0'12-+ %. ,123 $'#)2+ [Rp \ .%$ /<- *'+- 7<-& Ro = 1, m = 1, Dp = Dc = 1 '&# β = 0.
B
Q& %$#-$ /% 0'1)#'/- /<- $%;2+/&-++ %. /<- ,$-+-&/ (%#-1+> /<- *%&0-*/)%& #)..2+)%& (-*<'&)+( ,'$'#)3( 7'+ '##$-++-#> ;%/< /<$%23< /<- (-/<%# %. (%(-&/ RZU '+ 7-11 '+ 7)/< /<'/ %. 3-&-$'14 )D-# #)+,-$+)%& (%#-1 RXUB 9)(21'/)%& $-+,%&+- /% /<)+ -..-*/ <'+ ;--& -&*',+21'/-# )& N)3BYB J<0'12- %. /<- #)+,-$+)%& *%-.I*)-&/ [Da \ %;/')&-# .$%( /<- G(-/<%# %. (%(-&/F '&# /<- #)+,-$+)%& *%-.I*)-&/ [M2 \ '*82)$-# .$%( G3-&-$'1)D-# #)+,-$+)%& (%#-1F .%$ ;%/< L-7/%&)'& '&# K'++%& 1)84 2)# *'+- +/2#)-+ 2&#-$ /<- +'(- ,'$'(-/$)* -&0)$%&(-&/> ")-1#+ +/$)O)&3 82'&/)/'/)0- '&# 82'1)/'/)0$-+-(;1'&*-+ )&#-,-&#-&/ %. /<- (%#-1+ +% *<%+-& /% /<)+ -&#B
ST
ACCEPTED MANUSCRIPT
D & E
5
2
4
PT
*+,-().&)=".F#.G
%&''()=".F#.G
RI
3
SC
2
1
/ :4 1/
NU
Dispersion Coefficient
6
:3
:2
1/
:1
1/
/
1/
t
1
1/
1/
MA
N)32$- Yg k'$)'/)%& %. #)+,-$+)%& *%-.I*)-&/ 7)/< /)(- %;/')&-# ;" (%(-&/ (-/<%# [+%1)# 1)&-\ '&# 3-&-$'1)D-# #)+4 ,-$+)%& (%#-1 [#'+<-# 1)&-\B J<- *<%+-& ,'$'(-/-$ 0'12-+ '$-g m = 1> Dp = Dc = 1B N%$ /<- *'+- %. L-7/%&)'& 1)82)#> τy = 0, Ro = 1 '&# .%$ K'++%& 1)82)#> τy = 0.04, Ro = 1.
TE
D
N%$ /<- $-+/ %. %2$ '&'1"+)+> -..-*/ %. #)..-$-&/ ,'$'(-/-$+ %& /$'&+,%$/ *%-.I*)-&/+ <'0- ;--& -M'()&-# ;" ?$)+F+ (-/<%# %. (%(-&/+B
τ
>=&=?
/
1/
:1
1/
1
1/
>=@=?
;
/
/
1
2
3
4
5
t
6
7
8
9
1/
>=B=?
A<10/ 101 102
/0/4 /0/6
M0
M0
β
1
1/
M0
1
1/
AC CE P
9PHP LQ.,!(6# .)#*R.-#(& J<- 0'$)'/)%& %. -M*<'&3- *%-.I*)-&/ [M0 (t)\ 7)/< /)(- <'+ ;--& )&/-$,$-/-# /<$%23< N)3BZ 2&#-$ /<- )(,'*/ %. $-'*/)%& ,'$'(-/-$+ [β\> ")-1# +/$-++ [τy \ '&# 0)+*%+)/" $'/)% [m\B
1/
/
/
1
2
3
4
5
t
6
7
8
9
1/
1/
/
1
2
3
4
5
t
6
7
8
9
1/
N)32$- Zg J)(- '++-++(-&/ %. -M*<'&3- *%-.I*)-&/+ [M0 \ )& ' ,2$-1" K'++%& 1)82)# 6%7 .%$ /<- IM-# Ro = Dp = Dc = 1> ['\ 0'$)'/)%& %. ';+%$,/)%& ,'$'(-/-$ [β\ 7<-& m = 1.1 '&# τy = 0.04l [;\ 0'$)'/)%& %. ")-1# +/$-++ [τy \ 7<-& m = 1.1 '&# β = 1l [*\ 0'$)'/)%& %. 0)+*%+)/" $'/)% [m\ 7<-& τy = 0.04 '&# β = 1B
N)3BZ['\ *1-'$1" +<%7+ /<'/ /<- )&*$-'+- )& β ('O-+ '& '23(-&/'/)%& %. M0 B J<- )&*$-'+- )& M0 7)/< β *'& ;- '//$);2/-# /% /<- -M,'&#)&3 1%++ %. +%12/- '/ /<- 7'11 '+ β )&*$-'+-+B Q/ )+ /% ;&%/-# <-$- /<'/ M0 -M)+/+ #2- /% /<- ,$-+-&*- %. *<-()*'1 $-'*/)%& 7<)*< )+ )&#-,-&#-&/ %. τy '&# m '+ +--& )& N)3+BZ[;> *\B i8B[VS'\ .2$/<-$ -12*)#'/-+ /<'/ /<- ,-$),<-$'1 1'"-$ -..-*/ %& M0 #% &%/ SV
ACCEPTED MANUSCRIPT 0'$" /<- %2/*%(-+ '*82)$-# .%$ ,2$-1" K'++%& 1)82)# 6%7B
/048
>=&=?
/044
/046 /044
/0/4 /0/6
/042
/042
/038
/038 /036
/036
/034
/042 /04
/036 /034 /032
/03
/03
:3
1/
:2
1/
:1
t
1/
/
1/
:3
TE 1
1/
:3
:2
1/
1/
:1
t
1/
/
1/
1
:4
1/
AC CE P >=@=?
β
2
Da
105
1
/05
:3
1/
:2
1/
:1
t
1/
/
1/
1
1/
1/
:1
t
1/
/
1
1/
1/
C=1/
>=G=?
τ;
205
2
2
105
105
1
1
/05
/05
:3
1/
:2
1/
:1
t
1/
/
1/
1
1/
/ :4 1/
>=H=?
A<10/ 101 102
205
/0/4 /0/6
/ :4 1/
:2
1/
:3
C=1/
205
:3
1/
:3
C=1/
/ :4 1/
:4
1/
Da
/03 :4 1/
/038
/032
/032
Da
/044
D
/034
>=+=?
A<10/ 101 102
/046
MA
M1
M1
/04 /04
/048
>=B=?
NU
β
/046
τ;
M1
/048
SC
RI
PT
9POP M)(A#.&-)( !(> >-/$#%/-)( .)#*R.-#(&/ N%$ ,2$-1" K'++%& '&# ,-$),<-$'1 1'"-$ 1)82)# 6%7+> 0'$)'/)%&+ %. /<- *%&0-*/)0- *%-.I*)-&/ [M1 (t)\ '&# ',,'$-&/ #)+,-$+)%& *%-.I*)-&/ [Da (t)\ 7)/< $-+,-*/ /% /)(- [t\ '$- +<%7& )& N)3B] .%$ 0'$)%2+ -+/)('/)%&+ %. /<- $-'*/)%& ,'$'(-/-$ [β\> ")-1# +/$-++ [τy \ '&# 0)+*%+)/" $'/)% [m\B J<1)&-+ 7)/< '+/-$)+O <'0- ;--& +<%7& )& I32$-+ #2- /% K'++%& 1)82)# 6%7> 7<-$-'+ $-+/+ '$- #2- /% ,-$),<-$'1 1'"-$ 1)82)#B ? 7)#- $'&3- %. /)(- )+ *%&+)#-$-# /% /$'*O /<- '+"(,/%/)* ',,$%'*< %. /<#)+,-$+)%& *%-.I*)-&/B
:3
1/
:2
1/
:1
t
1/
/
1/
N)32$- ]g J)(- '++-++(-&/ %. /$'&+,%$/ *%-.I*)-&/+ [M1 '&# Da \ )& ,2$-1" K'++%& 1)82)# 6%7 [ '+/-$)+O 1)&-+\ .%$ /
5$%I1-+ )& /<- N)3+B]['> ;\ $-,$-+-&/ /<- ;-<'0)%2$ %. M1 '&# Da .%$ +'(- 0'12-+ %. β 0'$)'/)%&B N$%( N)3+B]['> ;\> )/ )+ %;+-$0-# /<'/ M1 )&*$-'+-+ 7<)1- Da #-*$-'+-+ (%&%/%&)*'11" 7)/< β 2&/)1 )/ $-'*<-+ /<-)$ $-+,-*/)0- +/-'#" +/'/- 0'12-+B 9)&*- /<- 7'11 $-'*/)%& -&<'&*-+ /<- /$'&+,%$/ %. +%12/- /%7'$#+ /<- 7'11> +2;+-82-&/1" M1 )&*$-'+-+ 7)/< βB N$%( N)3B]['\> )/ *'& ;- +')# /<'/ M1 $-'*<-+ -'$1)-$ /% )/+ '+"(,/%/)* 0'12- 7)/< /<- #-*$-'+- )& β> +)&*- 7-'O ';+%$,/)%& ('O-+ '& SX
1
1/
ACCEPTED MANUSCRIPT
AC CE P
TE
D
MA
NU
SC
RI
PT
-82)1);$)2( ;-/7--& /<- #)..2+)%& %. +%12/- '&# /<- 7'11 $-'*/)%& %**2$+ 7)/<)& +('11 /)(- .$'(-B Q& N)3B][;\> /<- ,<"+)*'1 ,-$*-,/)%& ';%2/ /<- #-*$-(-&/ %. Da 7)/< )&*$-'+)&3 β 1-'#+ /% '& )&*$-'+)& /<- 82'&/)/" %. (%1-+ %. $-'*/)0- ('/-$)'1 2&#-$3%)&3 *<-()*'1 $-'*/)%& %$ ';+%$,/)%& $-+21/)&3 )& ' *<'&3- )& /<- *%&*-&/$'/)%& #)+/$);2/)%& '*$%++ /<- /2;- 7)/< /<- +)(21/'&-%2+ *%&+-82-&*%. ' #$%, )& #)+,-$+)%& *%-.I*)-&/B J<- ,$-+-&/ +/2#" ")-1#+ '& 2&)82- %;+-$0'/)%& '+ -0)#-&/ )& /<- N)3+B]['> ;\ )+ /<'/ /<- ,$-+-&*- %. ' ,-$),<-$'1 1'"-$ ('O-+ 1'$3-$ 0'12-+ )& M1 '&# Da '+ *%(,'$-# /% +)&31- K'++%& 1'"-$ 1)82)# 6%7B J<)+ /$-&# )+ '1+% -0)#-&/ .$%( +2;+-82-&/ N)3B _B J<- N)3+B][*> #\ <'0- ;--& #$'7& /% -(,<'+)D- /<- -..-*/ %. τy %& M1 '&# Da )& %$#-$ /% <)3<1)3 #\> )/ )+ *1-'$1" -0)#-&/ /<'/ M1 '&# Da #-*$-'+-+ 7)/< )&*$-'+)&3 τy B J<)+ )+ '//$);2/-# /% /<- .'*/ /<'/ /<- )&*$-'+- )& τy [%$ )& %/<-$ 7%$#+> '& )&*$-'+- )& /<- &%&4L-7/%&)'& ,'$'(-/-$\ #-*$-'+-+ /<- 0-1%*)/" %. /<- 62)#B J<- )(,%$/'&*- %. ,-$),<-$'1 1'"-$ *'& ;- ',,$-*)'/-# .$%( /<- %;+-$0'/)%& %. /<- )&*$-'+-# 0'12%. /<- *%-.I*)-&/+ '+ *%(,'$-# /% K'++%& 1)82)# '1%&-B J<- $-'+%& ;-<)&# /<-+- )&*$-(-&/+ '$#2- /% L-7/%&)'& 1'"-$ -..-*/> 7<)*< <'+ 1-++-$ 0)+*%+)/" /<'& &%&4L-7/%&)'& K'++%& 1)82)#B k)+*%+)/" $'/)% [m\ )+ #-I&-# '+ /<- $'/)% ;-/7--& K'++%& 0)+*%+)/" [µc \ '&# L-7/%&)'& 0)+*%+4 )/" [µn \B Q& /<- ,$-+-&/ +/2#"> /<- ',,-'$'&*- %. 0)+*%+)/" $'/)% [m\ )+ #2- /% /<- ,-$),<-$'1 1'"-$B J<- (')& '#0'&/'3- %. 2+)&3 m %& /<- +/2#" %. #)+,-$+)%& )+ /<'/ )/ #)$-*/1" #)+*1%+-+ '+ /% <%7 +/$%&3 /<- )&62-&*- %. µc )+ *%(,'$-# /% µn B J<- /)(- -0'12'/)%& %. *%&0-*/)%& *%-.I*)-&/ [M1 \ '&# #)+,-$+)%& *%-.I*)-&/ [Da \ 2&#-$ /<- -..-*/ %. m <'0- ;--& ,$%H-*/-# )& N)3+B][-> .\B ?+ *'& ;- )&.-$$-# .$%( N)3+B][-> .\ > 7- *'& %;+-$0-# /<'/ M1 '&# Da )&*$-'+- 7)/< /)(- [t\ .%$ ' +<%$/ ,-$)%# ;-"%&# 7<)*< /<-" '//')&+ ' +/-'#" +/'/- /$-&#B ?3')&> '+ m )&*$-'+-+ /<- ('3&)/2#-+ %. M1 '&# Da #-*$-'+-B J<)+ <',,-&+ )& 1)3 7<-& m = 1.1 '&# t = 0 /<- )&*$-(-&/ )& M1 '&# Da '$- SBSV '&# S_^B] /)(-+ /<- 0'12- *%$$-+,%&#)&3 /% t = 10> ;2/ .%$ ' /<$--41'"-$ 1)82)# /<)+ )&*$-(-&/ )+ SBSY '&# T`` /)(-+ $-+,-*/)0-1"B f<)1- *%(,'$)&3 /<- K'++%& 1)82)# '&# /<- *%$$-+,%&#)&3 /<$--41'"-$ 1)82)# +/$-'(> )/ )+ %;+-$0-# /<'/> ;%/< M1 '&# Da )&*$-'+- (%$- )& /<- *'+- %. /<$--41'"-$ 1)82)# +/$-'( '+ *%(,'$-# /% /<'/ /<- K'++%& 1)82)# -&0)$%&(-&/B N$%( /<- )&0-+/)3'/)%&> )& N)3+B][-> .\> )/ *'& ;- *'1*21'/-# /<'/ .%$ /<- /)(- )&/-$0'1 R0, 10U '&# 0)+*%+)/" $'/)% [m = 1.1\> /<- 0'12- %. M1 )&*$-'+-+ )& /<- $'&3- [0.0365 − 0.0374\ 7<-$-'+ /<- 0'12- %. Da )&*$-'+-+ )& /<- $'&3- [1.49 × e−6 − 4.5 × e−5 \B 9PSP +#!( .)(.#(&%!&-)( C%)R2# N)3+B^ ,$-+-&/+ /<- +/2#" %. (-'& *%&*-&/$'/)%& #)+/$);2/)%& '*$%++ 'M)'1 #)+/'&*- .%$ +/-'#" &'/2$- %. +/$-'( '/ ' IM-# /)(- t = 1.5B N)3B^['\ #-'1+ ';%2/ /<- -..-*/ %. ';+%$,/)%& ,'$'(-/-$ [β\ %& (-'& *%&*-&/$'/)%& #)+/$);2/)%& 7<-$-'+ N)3B ^[;\ '&# N)3B^[*\ (-'+2$- /<- )&62-&*- #2/% τy '&# m $-+,-*/)0-1"B :2- /% *%&+/'&/ #-,1-/)%& %**2$$)&3 '/ /<- /2;- 7'11 )& ;%/< /<- *'+-+> /<,-'O %. /<- (-'& *%&*-&/$'/)%& #-*$-'+-+ '+ β )&*$-'+-+B Q/ ('" ;-*%(- 6'/ .%$ 1'$3- ';+%$,/)%&> /<2+ /<- #)+,-$+)%& *%-.I*)-&/ ('" &%/ <'0- ,$'*/)*'1 )(,%$/'&*- .%$ 1'$3- βB f<)1- %;+-$0)&3 N)3+B^[;> *\> )/ )+ .%2&# /<'/ '+ τy '&# m )&*$-'+-> /<- ,-'O %. /<- (-'& *%&*-&/$'/)%& #)+/$);2/)%& '1+% )&*$-'+-+B N$%( N)3+B^['> ;> *\> )/ *'& ;- *1-'$1" %;+-$0-# /<'/ ,-$),<-$'1 1'"-$ ,$-+-&*- ('O-+ .'+/-$ #-*$-(-&/+ )& ,-'O 1-&3/< '+ *%(,'$-# /% ,2$-1" K'++%& 1)82)# 6%7B
SY
ACCEPTED MANUSCRIPT :3
:4
C=1/
7
β
5
205
4
2
/0/4 /0/6
PT
3
3
105 2
SC
1
1
/05
/
:/02
:/015
:/01
:/0/5
/
/0/5
/01
z − zg
/015
/02
:4
C=1/
/
z − zg
/0/5
/01
/015
/02
/025
>=B=?
D
TE 3
AC CE P
Cm
4
:/025 :/02 :/015 :/01 :/0/5
A<10/ 101 102
6
5
/025
MA
7
NU
/ :/025
>=@=?
τ;
6
Cm
Cm
305
>=&=?
RI
4
C=1/
2
1
/ :/03 :/025 :/02 :/015 :/01 :/0/5
/
/0/5
z − zg
/01
/015
/02
/025
/03
N)32$- ^g k'$)'/)%& %. (-'& *%&*-&/$'/)%& #)+/$);2/)%&> '/ ' 3)0-& )&+/'&/ t = 1.5 .%$ K'++%& 1)82)# ['+/-$)+O 1)&-+\> τy = 0.04, Ro = 1 '&# 1'"-$-# 1)82)#> τy = 0.04, Ro = 0.94B h/<-$ #-+*$),/)%& '$- '+ )& N)3B]B
9P
' #)+*-$&);1- )(,'*/ <'+ ;--& %;+-$0-# /<$%23< N)3B_ '&# N)3BS`B J<- -M*<'&3- *%-.I*)-&/ [M0 \ )+ )&#-,-&#-&/ %. Ro > 7<)*< <'+ +<%7& )& N)3B_['\ '&# )/ *%(,1-/-1" '3$--+ %& /<- ,<"+)*+ %. /<- ,$-+-&/ (%#-1B N)3+B_[;> *\ #-,)*/ /<'/ /<- *%&0-*/)%& *%-.I*)-&/ [M1 \ '&# #)+,-$+)%& *%-.I*)-&/ [Da \ ;%/< #-*$-'+- '+ 0'12- %. Ro )&*$-'+-+ [%$ ,-4 $),<-$'1 1'"-$ /<)*O&-++ #-*$-'+-+\B J<- ,<"+)*'1 -M,1'&'/)%& ';%2/ +2*< O)&# %. ;-<'0)%2$ )+ /<'/ /<- 1'$3-$ 0'12- %. Ro ('O-+ /<- +('11-$ $-3)%& .%$ L-7/%&)'& 1)82)#> <-&*- /<- 0-1%*)/" %. /<1)82)# #-*$-'+-+ $-+21/)&3 )& /<- +('11-$ *%&0-*/)%& '&# #)..2+)%& $'/-B N)3BS` *1-'$1" )&#)*'/-+ /<'/ ,-'O %. /<- (-'& *%&*-&/$'/)%& #)+/$);2/)%& )&*$-'+-+ '+ Ro )&*$-'+-+B J<)+ )+ #2- /% +1%7 /$'*-$ +,$-'#)&3 %0-$ /<- -&/)$- $-3)%&B SZ
ACCEPTED MANUSCRIPT 1
1/
/044
I
>=@=?
I
>=&=?
(
/094 /098
/094 /098
/042
PT
/038
RI
M1
M0
/04
SC
/036
/034
/
/
1
2
3
4
5
6
t
7
8
9
:3
1/
:3
:2
1/
1/
I
TE
/05
/ :4 1/
:3
1/
:2
1/
:1
t
1/
/
1/
1
1/
N)32$- _g J)(- '++-++(-&/ %. /$'&+,%$/ *%-.I*)-&/+ [M0 > M1 '&# Da \ .%$ 0'$)%2+ 0'12-+ %. Ro .%$ /<- *'+- 7<-& τy = 0.04, m = 1.1, β = 1, Dp = Dc = 1B
S]
1
1/
1
AC CE P
Da
105
(
t
/
1/
D
2
:1
1/
>=B=?
MA
C=1/
/032 :4 1/
NU
1/
ACCEPTED MANUSCRIPT :4
C=1/
I
5
/094 /098
PT
3
RI
Cm
4
SC
2
/ :/03 :/025 :/02 :/015 :/01 :/0/5
NU
1
/
/0/5
z − zg
/01
/015
/02
/025
/03
MA
N)32$- S`g k'$)'/)%& %. (-'& *%&*-&/$'/)%& #)+/$);2/)%&> '/ ' 3)0-& )&+/'&/ t = 1.5 .%$ /<- *'+- 7<-& τy = 0.04, m = 1.1, β = 1, Dp = Dc = 1B
AC CE P
TE
D
9PTP L**#.& )* ')2#.42!% >-**4/-A-&7 )( >-/$#%/-)( .)#*R.-#(& J<- (%1-*21'$ #)..2+)0)/" %$ #)..2+)%& *%-.I*)-&/ [D\ )+ ' O)&-('/)* ,$%,-$/" %. /<- 62)# 7<)*< )(,1)-+ /<'/ /<- ('++ %. /<- +2;+/'&*- #)..2+-+ /<$%23< ' 2&)/ +2$.'*- )& ' 2&)/ /)(- '/ ' *%&*-&4 /$'/)%& 3$'#)-&/ %. 2&)/"B 9%> #)..2+)%& *%-.I*)-&/ [D\ $-,$-+-&/+ /<- +,$-'#)&3 ,$%*-++> /<- 1'$3-$ 0'12- %. D )&#)*'/-+ /<- $',)# +,$-'#)&3 %. /<- /$'*-$ *2$0- 7<-$-'+ +('11-$ D )(,1)-+ +1%7 +,$-'#4 )&3 %. +%12/- )& /<- 1)82)#B N)3BSS +<%7+ /<- -..-*/ %. D %& /<- ',,'$-&/ #)+,-$+)%& *%-.I*)-&/ [Da \B
S^
ACCEPTED MANUSCRIPT :3
:3
C=1/ 205
C=1/
>=&=?
D <1KD <1 J
B
J
/015 /02
DJ
105
1
RI
105
PT
2
Da
Da
2
SC
1
/05
/05
:3
1/
:2
1/
:1
t
/
1/
1
1/
1/
:3
C=1/
:3
:2
1/
1/
MA
3
/ :4 1/
NU
/ :4 1/
>=@=?
D
D
>=B=?
105
AC CE P
1
/05
/ :4 1/
:3
1/
:2
1/
:1
t
1/
/
1/
1
1/
N)32$- SSg k'$)'/)%& %. #)+,-$+)%& *%-.I*)-&/ [Da \ 7)/< /)(- .%$ 0'$)%2+ 0'12-+ %. #)..2+)0)/" .%$ /<- *'+- 7<-& τy = 0.04, Ro = 0.94, m = 1.1, β = 1. [;\ .%$ Dc = 0.7 [*\ .%$ Dp = 0.15B
S_
1
1/
TE
Da
2
/07/ /075
t
/
1/
D
205
:1
1/
ACCEPTED MANUSCRIPT :4
1/ 9 8
C=1/
D <1KD <1 J
B
DJ
PT
7
Cm
6 5
RI
4 3
SC
2
/ :/03 :/025 :/02 :/015 :/01 :/0/5
NU
1
/
/0/5
z − zg
/01
/015
/02
/025
/03
MA
N)32$- STg k'$)'/)%& %. (-'& *%&*-&/$'/)%& #)+/$);2/)%&> '/ ' 3)0-& )&+/'&/ t = 1.5 .%$ /<- *'+- 7<-& τy = 0.04, Ro = 0.94, m = 1.1, β = 1.
AC CE P
TE
D
N)3+BSS['4*\ *1-'$1" +<%7 /<'/ 1'$3-$ 0'12-+ )& #)..2+)0)/" ;-1%&3)&3 .$%( /<- $-3)%& Rc %$ ,123 $-3)%& Rp > #)+,-$+)%& *%-.I*)-&/ #-*$-'+-+B J<- $-'+%& .%$ /<)+ #-*$-(-&/ )+ #2- /% /<- .'+/ +,$-'#4 )&3 %. (%1'$ 62M 7<)*< *%(,1-/-1" '3$--+ %& /<- ,<"+)*+ %. /<- +"+/-(B N)3BST #-(%&+/$'/-# /<(-'& *%&*-&/$'/)%& #)+/$);2/)%& '3')&+/ /<- 'M)'1 #)+/'&*-B ?+ -0)#-&/ .$%( /<- N)3BST> 1'$3-$ #).4 .2+)0)/" ('O-+ *%&/)&2%2+ +,$-'#)&3 %. +%12/- )& /<- $'#)'1 #)$-*/)%&> '+ ' $-+21/> #-*$-'+- )& ,-'O <-)3 XZUB J<- )&0-+/)3'/)%& %. #)+,-$+)%& 2&#-$ /<- ,$%,%+-# (%#-1 *'& ;- #--(-# /% ,$%0)#- '& )&+)3 1%7 ,'$/)*1-+ '//$)/)%& '&# 1%7 ,),- -$%+)%& RX]UB A-*'2+- %. /<- ,$-+-&*- %. ,-$),<-$'1 1'"-$> #)..-$-&/ #)..2+)0)/" <'+ /% ;- *%&+)#-$-# '&# ' &-7 ,'$'(-/-$ 0)DB> 0)+*%+)/" $'/)% ',,-'$+ )& /<- +/2#" %. +%12/- #)+,-$+)%& 7<)*< )+ 0-$" )(,%$/'&/ .%$ /<- +/2#)-+ $-1'/-# /% (%$,<%1%3" RX^UB J<- $-+21/+ %. /<- ,$-+-&/ +/2#" *1-'$1" -+/';1)+<-+ )+ $)3
M) H,+'5.#/,+# N$%( /<- ,$-0)%2+1" +/'/-# #)+*%2$+-> )/ *'& ;- $-'#)1" *%&*12#-# /<'/ /<- 2/)1)/" %. /<- 2,4 #'/-# +*)-&/)I* (%#-1 #-+*$);-# )& /<)+ ,',-$ <'+ ' +)3&)I*'&*- )& /<- +/2#" %. +%12/- #)+,-$+)%&B N%$ '11 /)(- '&'1"+)+> /<- ?$)+4A'$/%&F+ G(-/<%# %. (%(-&/F <'+ ;--& ,-$.%$(-# %& /<- /)(- #-4 0-1%,-# /$'&+,%$/ *%-.I*)-&/+ ;" 82'1)/'/)0-1" '&# 82'&/)/'/)0-1"> 7<-$- ' /7%41)82)# (%#-1 <'+ T`
ACCEPTED MANUSCRIPT
RI
PT
;--& *%&+)#-$-# /<$%23< ' $)3)# ,),- 2&#-$ *'/'1"/)* 7'11 $-'*/)%&B ?3')&> /<- G3-&-$'1)D-# #)+,-$4 +)%& (%#-1F %. 9'&O'$'+2;$'('&)'& '&# W)11 <'0- ;--& -(,1%"-# '&# /<- $-+21/+ %;/')&-# .$%( ;%/< (-/<%#+ '$- )#-&/)*'1> 7<)*< -&+2$-+ /<- '2/<-&/)*)/" %. /<- $-+21/+ %;/')&-# )& /<)+ +/2#"B C'/<-('/)*'1 +/2#)-+ %& ;1%%# 6%7 7-$- (%+/1" +/2#)-# ;" K'++%& (%#-1> ;2/ /<- ,$-+-&*- %. ,-$),<-$'1 1'"-$ %& ;1%%# <'+ '1$-'#" ;--& '**%(,1)+<-# -M,-$)(-&/'11" ;" #)..-$-&/ '2/<%$+B J<)(,%$/'&*- %. ,-$),<-$'1 1'"-$ %& ;1%%# 6%7 <'+ ;--& (-&/)%&-# ;" .-7 '2/<%$+ 7<)1- +/2#")&3 $-+)+/'&*- /% 6%7> 7'11 +<-'$ +/$-++> ,123 *%$- $'#)2+ -/*B 9%> /<- +/2#" %. #)+,-$+)%& 2&#-$ +2*< O)&# %. +)/2'/)%&+ <'+ '& )(,%$/'&/ ,<"+)*'1 '**-,/'&*-B J<- .%11%7)&3 3-&-$'1 *%&*12+)%&+ *'& ;- ('#- .$%( /<- +/2#"g
NU
SC
[)\ N%$ '11 /)(- ,-$)%#> /<- ('3&)/2#- %. -M*<'&3- '&# *%&0-*/)%& *%-.I*)-&/ )&*$-'+- .%$ )&4 *$-'+)&3 ';+%$,/)%& $'/- '/ /<- ,),- ;%2&#'$" /<%23< 1'$3-$ ';+%$,/)%& ('O-+ +('11-$ #)+,-$4 +)%& *%-.I*)-&/ '&# $-'*<-+ /% )/+ +/-'#" 0'12- )& ' +<%$/ #2$'/)%& %. /)(-B [))\ J<- -M*<'&3- *%-.I*)-&/ )+ )&#-,-&#-&/ %. /<- ")-1# +/$-++> 0)+*%+)/" $'/)%> '&# ,-$),<-$'1 1'"-$ /<)*O&-++B
MA
[)))\ f<'/-0-$ ;- /<- +%10-&/ [,2$-1" K'++%& 1)82)# %$ /<$--41'"-$ 1)82)#\> /<- 0'12-+ %. *%&0-*/)%& '&# #)+,-$+)%& *%-.I*)-&/+ '17'"+ #-*$-'+- '+ ")-1# +/$-++ '&# 0)+*%+)/" $'/)% )&*$-'+-B
D
[)0\ J<- ,$-+-&*- %. ,-$),<-$'1 1'"-$ ('O-+ $)+- )& /<- 0'12-+ %. *%&0-*/)%& '&# #)+,-$+)%& *%-.4 I*)-&/+ '+ *%(,'$-# /% +)&31- 1'"-$ K'++%& 1)82)# 6%7B
AC CE P
TE
[0\ 9/$%&3-$ /<- ';+%$,/)%& <'0- (2*< 6'/ )& (-'& *%&*-&/$'/)%& #)+/$);2/)%& '&# '1+% /<- ,-'O %. /<- (-'& *%&*-&/$'/)%& #)+/$);2/)%& )&*$-'+-+ 7<-& ")-1# +/$-++ '&# 0)+*%+)/" $'/)% )&4 *$-'+-B [0)\ J<- *%&0-*/)%& '&# #)+,-$+)%& *%-.I*)-&/ ;%/< #-*$-'+- '+ /<- ,-$),<-$'1 1'"-$ /<)*O&-++ #-*$-'+-B [0))\ ?+ #)..2+)0)/" )&*$-'+-+> #)+,-$+)%& *%-.I*)-&/ #-*$-'+-+> *%&+-82-&/1" ,-'O %. /<- (-'& *%&*-&/$'/)%& #)+/$);2/)%& '1+% 3%-+ #%7&7'$#B [0)))\ J<- I&#)&3+ %. /<- +/2#" '$- '1+% )& '3$--(-&/ 7)/< /<%+- %. -M)+/)&3 7%$O+B
@.$.%3 #',:3 ,6 #$.-A J<- ,'$'(-/-$+ )&*12#-# )& /<- ,$-+-&/ +*%,- %. +/2#" <'0- ;--& +)&321'$1" (%/)0'/-# /% -(21'//<- '*/2'1 &'/2$- %. /<- ,<"+)%1%3)*'1 +)/2'/)%& '$)+)&3 )& ;1%%# 6%7 '&'1"+)+ '+ .'$ '+ ,%++);1- .$%( ' ('/<-('/)*'1 ,-$+,-*/)0-B !%7-0-$ '*/2'1 ,<"+)%1%3)*'1 *%&#)/)%&+ '//$);2/)&3 /% /<- *'+- +/2#" *'& ;- .2$/<-$ 0'$)-# /% -&<'&*- /<- I#-1)/" %. /<- -+/';1)+<-# (%#-1 )& .2$/<-$ $-+-'$*
!'=+,N53-D343+$ J<- '2/<%$+ 3$-'/1" '*O&%71-#3- /<- 0'12';1- +233-+/)%&+ '&# *%((-&/+ ;" 5$%.B 5$'#--, WB 9)##<-+<7'$> A'&3'1%$- b&)0-$+)/"> A'&3'1%$-> Q&#)'B h&- %. 2+ [92#), :-;&'/<\ )+ /<'&O.21 /% LQJ ?3'$/'1' 7<-$-'+ [?+<)+ O2('$ E%"\ )+ /<'&O.21 /% C!E:> Q&#)' .%$ I&'&*)'1 +2,,%$/ .%$ ,2$+2)&3 /<)+ 7%$OB f- /<'&O /<- $-.-$--+ .%$ 2+-.21 *%((-&/+ /<'/ $-I&-# /<- ,',-$ /% /<- ,$-+-&/ .%$(B
TS
ACCEPTED MANUSCRIPT
0'$)';1- #-I&-# ;" i8B[SV\
Hi
!-$()/- ,%1"&%()'1+
zg
*-&/$)%#
C (p)
pth (%(-&/ %. /<- #)+/$);2/)%& %. +%12/-
fi
ith #)+,-$+)%& .2&*/)%& [i = 0, 1, 2, 3, ...\
C
*%&*-&/$'/)%& %. /<- +%12/-
C0
)&)/)'1 *%&*-&/$'/)%& %. /<- +%12/-
Cm
(-'& *%&*-&/$'/)%& %. /<- +%12/-
D∗
(%1-*21'$ #)..2+)%& *%-.I*)-&/
Da
',,'$-&/ #)+,-$+)%& *%-.I*)-&/ .%$ (%(-&/+ (-/<%#
Dc
(%1-*21'$ #)..2+)0)/" .%$ K'++%& 6%7 $-3)%&
Dn
m M0 M1 M2
RI
SC
NU
MA
D
(%1-*21'$ #)..2+)0)/" .%$ L-7/%&)'& 6%7 $-3)%& (%1-*21'$ #)..2+)0)/" .%$ ,123 6%7 $-3)%&
AC CE P
Dp
PT
B
TE
O,43+'5&$.%3
0)+*%+)/" $'/)% -M*<'&3- *%-.I*)-&/ *%&0-*/)%& *%-.I*)-&/ ',,'$-&/ #)+,-$+)%& *%-.I*)-&/ .%$ 3-&-$'1)D-# #)+,-$+)%& (%#-1
Pe
5-*1-/ &2(;-$
R
$'#)2+ %. ,),-
r
$'#)'1 *%%$#)&'/-
Ro
$'/)% %. /<- *-&/$'1 *%$- $'#)2+ /% /<- &%$('1 ,),- $'#)2+
Rp
,123 *%$- $'#)2+
t
/)(-
u
'M)'1 0-1%*)/" '1%&3 /<- 6%7
u0
*-&/-$ 1)&- 0-1%*)/"
uc
'M)'1 0-1%*)/" .%$ K'++%& 6%7 $-3)%& TT
ACCEPTED MANUSCRIPT un
'M)'1 0-1%*)/" .%$ L-7/%&)'& 6%7 $-3)%&
up
'M)'1 0-1%*)/" .%$ ,123 6%7 $-3)%&
z
'M)'1 *%%$#)&'/-
0'$)';1- #-I&-# ;" i8B[SV\
θ
'&321'$ *%%$#)&'/-
δij
@$%&-*O-$ #-1/'
Rc
[0 ≤ r ≤ Ro \ K'++%& 1)82)# $-3)%&
Rn
[Ro ≤ r ≤ 1\ L-7/%&)'& 1)82)# $-3)%&
Rp
[0 ≤ r ≤ Rp \ ,123 6%7 $-3)%& 7)/<)& K'++%& 1)82)# $-3)%&
β
7'11 ';+%$,/)%& ,'$'(-/-$
µc
0)+*%+)/" %. /<- K'++%& 1)82)#
µn
0)+*%+)/" %. /<- L-7/%&)'&
µp
ν3 τc τn τy
∂p ∂z
SC
NU
MA
D
pth %$#-$ *-&/$'1 (%(-&/ +O-7&-++
AC CE P
ν2
RI
ψ
TE
PT
>%33=
O2$/%+)+ +<-'$ +/$-++ %. /<- K'++%& 1)82)# +<-'$ +/$-++ %. /<- L-7/%&)'& 1)82)# ")-1# +/$-++
,$-++2$- 3$'#)-&/
P) K363%3+'3# RSU WB J'"1%$> :)+,-$+)%& %. +%12;1- ('//-$ )& +%10-&/ 6%7)&3 +1%71" /<$%23< ' /2;-> 5$%*B EB 9%*B e%&#%&> 9-$B ? TS_ [S_YV\ S^ZmT`VB RTU WB J'"1%$> K%&#)/)%&+ 2&#-$ 7<)*< #)+,-$+)%& %. ' +%12/- )& ' +/$-'( %. +%10-&/ *'& ;- 2+-# /% (-'+2$- (%1-*21'$ #)..2+)%&> 5$%*B EB 9%*B e%&#%&> 9-$B ? TTY [S_YX\ X]VmX]]B RVU EB ?$)+> h& /<- #)+,-$+)%& %. ' +%12/- )& ' 62)# 6%7)&3 /<$%23< ' /2;-> 5$%*B EB 9%*B e%&#%&> 9-$B ? TVY [S_YZ\ Z]m]]B
TV
ACCEPTED MANUSCRIPT RXU EB 9'&O'$'+2;$'('&)'&> fB W)11> b&+/-'#" *%&0-*/)0- #)..2+)%& 7)/< )&/-$,<'+- ('++ /$'&+4 .-$> 5$%*B EB 9%*B e%&#%&> 9-$B ? VVV [S_]V\ SSYmSVTB
PT
RYU kB ?&'&/<'O$)+<&'&> fB W)11> ?B dB A'$#2<&> e'()&'$ #)+,-$+)%& )& *',)11'$)-+g 5'$/ QB ('/<-4 ('/)*'1 '&'1"+)+> ?1K h& /<- (-/<%# %. (%(-&/+ .%$ +%12/- #)+,-$+)%&> dB N12)# C-*
SC
RI
R]U dB E'&'> 5B C2$/<"> 9%12/- #)+,-$+)%& )& ,21+'/)1- *'++%& 62)# 6%7 )& ' /2;- 7)/< 7'11 ';+%$,4 /)%&> dB N12)# C-*
NU
R^U dB E'&'> 5B C2$/<"> b&+/-'#" +%12/- #)+,-$+)%& )& &%&4&-7/%&)'& 62)# 6%7 )& ' /2;- 7)/< 7'11 ';+%$,/)%&> 5$%*B EB 9%*B e%&#%&> 9-$B ? X]T [T`SZ\ T`SZ`T_XB R_U AB 9B C'D2(#-$> 9B 5'21> :)+,-$+)%& %. $-'*/)0- +,-*)-+ 7)/< $-0-$+);1- '&# )$$-0-$+);1- 7'11 $-'*/)%&+> !-'/ C'++ J$'&+.-$ X^ [T`ST\ _VVm_XXB
MA
RS`U AB 9B C'D2(#-$> 9B A'&#"%,'#<"'"> h& +%12/- #)+,-$+)%& .$%( '& -1-0'/-# 1)&- +%2$*- )& '& %,-&4*<'&&-1 6%7> dB i&3B C'/
D
RSSU 5B L'3'$'&)> AB 9-;'+/)'&> :)+,-$+)%& %. ' +%12/- )& ,21+'/)1- &%&4&-7/%&)'& 62)# 6%7 /<$%23< ' /2;-> ?*/' C-*
TE
RSTU AB 9B C'D2(#-$> 9B @B :'+> i..-*/ %. ;%2&#'$" $-'*/)%& %& +%12/- #)+,-$+)%& )& ,21+'/)1- 6%7 /<$%23< ' /2;-> dB N12)# C-*
AC CE P
RSVU EB ?$)+> h& /<- #)+,-$+)%& %. ' +%12/- )& ,21+'/)&3 6%7 /<$%23< ' /2;-> 5$%*B EB 9%*B e%&#%&> 9-$B ? TY_ [S_Z`\ V]`mV]ZB RSXU !B JB K211)&'&> h& /<- +%12/)%& %. (21/)*%(,%&-&/ #)+,-$+)%& ,$%;1-(+> K<-(B i&3B dB Y [S_]V\ T]YmT]_B RSYU 5B WB 9)##<-+<7'$> 9B C'&H2&'/<> b&+/-'#" *%&0-*/)0- #)..2+)%& 7)/< <-/-$%3-&-%2+ *<-()4 *'1 $-'*/)%& )& ' ,1'&-4,%)+-2)11- 6%7 %. ' ()*$%,%1'$ 62)#> Q&/B dB i&3B 9*)B V^ [T```\ ]ZYm ]^VB RSZU iB C'$2c+)n*45'1%O'> QB 5'cD'&)&> h& $-'*/)0- +%12/- /$'&+,%$/ /<$%23< ' *2$0-# ,),-> ?,,1B C'/ C%#-1)&3 %. +%12/- #)+,-$+)%& )& ' *)$*21'$ ,),- I11-# 7)/< ()*$%,%1'$ 62)#> C'/ kB :-0)3&-> KB k'& :2)H&> E)3%$%2+ 2,+*'1)&3 %. /<- $-'*/)0- 6%7 /<$%23< ' ,%$-> 2&#-$ #%()&'&/ ,-*1-/ '&# #'(O%<1-$ &2(;-$+> 9Q?C dB C'/ i..-*/ %. 7'11 %+*)11'/)%& %& #)+,-$+)%& )& 'M)4+"((-/$)* 6%7+ ;-/7--& /7% *%'M)'1 *"1)&#-$+> o?CC4dB ?,,1B C'/ AB 9B C'D2(#-$> h& *%&/'()&'&/ #)+,-$+)%& )& 2&+/-'#" 3-&-$'1)+-# *%2-//- 6%7> Q&/B dB i&3B 9*)B V] [S___\ SX`]mSXTVB TX
ACCEPTED MANUSCRIPT RTSU 9B A'&#"%,'#<"'"> AB 9B C'D2(#-$> b&+/-'#" *%&0-*/)0- #)..2+)%& )& ' ,21+'/)1- 6%7 /<$%23< ' *<'&&-1> ?*/' C-*
PT
RTTU ?B 9'$O'$> WB d'"'$'('&> J<- -..-*/ %. 7'11 ';+%$,/)%& %& #)+,-$+)%& )& %+*)11'/%$" 6%7 )& '& '&&212+g ',,1)*'/)%& /% ' *'/<-/-$)D-# '$/-$"> ?*/' C-*
RI
RTVU AB 9B C'D2(#-$> @B @B C%&#'1> h& +%12/- /$'&+,%$/ )& %+*)11'/%$" 6%7 /<$%23< '& '&&21'$ ,),- 7)/< ' $-'*/)0- 7'11 '&# )/+ ',,1)*'/)%& /% ' *'/<-/-$)D-# '$/-$"> pB dB C-*
SC
RTXU CB @B 9<'$,> 9<-'$4'23(-&/-# #)+,-$+)%& )& &%&4&-7/%&)'& 62)#+> ?&&B A)%(-#B i&3B TS [S__V\ X`]mXSYB
NU
RTYU 9B ?3$'7'1> WB d'"'$'('&> L2(-$)*'1 +)(21'/)%& %. #)+,-$+)%& )& /<- 6%7 %. ,%7-$ 1'7 62)#+ )& *2$0-# /2;-+> ?,,1B C'/
MA
RTZU EB :'+<> WB d'"'$'('&> @B C- 9<-'$ '23(-&/-# #)+,-$+)%& %. ' +%12/- )& ' *'++%& 62)# 6%7)&3 )& ' *%)/> ?&&B A)%(-#B i&3B T^ [T```\ V]VmV^YB
D
RT]U 5B WB 9)##<-+<7'$> 9B C'$O'&#-> 9B C'&H2&'/<> i..-*/ %. )&/-$,<'+- ('++ /$'&+.-$ %& 2&4 +/-'#" *%&0-*/)0- #)..2+)%&g 5'$/ Q> ,1'&-4,%)+-2)11- 6%7 %. ' ,%7-$41'7 62)# )& ' *<'&&-1> K<-(B i&3B K%((2&B S^` [T```\ S^]mT`]B
AC CE P
TE
RT^U 5B WB 9)##<-+<7'$> 9B C'$O'&#-> 9B C'&H2&'/<> i..-*/ %. )&/-$,<'+- ('++ /$'&+.-$ %& 2&4 +/-'#" *%&0-*//0- #)..2+)%&g 5'$/ QQ <'3-& ,%)+-2)11- 6%7 %. ' ,%7-$ 1'7 62)# )& ' /2;-> K<-(B i&3B K%((2&B S^` [T```\ T`_mTT_B RT_U 5B L'3'$'&)> WB 9'$%H'(('> WB d'"'$'('&> iM'*/ '&'1"+)+ %. 2&+/-'#" *%&0-*/)0- #)..2+)%& )& *'++%& 62)# 6%7 )& '& '&&212+m',,1)*'/)%& /% *'/<-/-$)D-# '$/-$"> ?*/' C-* CB :-0)11-> 521+'/)1- 6%7 %. &%&4&-7/%&)'& 62)#+ /<$%23< '$/-$)'1 +/-&%+-+> dB A)%(-* qB q2> @B E'(&'$)&-> 5B !%+O)&+> L2(-$)*'1 )&0-+/)3'/)%& %. ,<"+)%1%3)*'11" $-'14 )+/)* ,21+'/)1- 6%7 /<$%23< '$/-$)'1 +/-&%+)+> dB A)%(-* ?& -82'/)%& .%$ /<- 6%7 %. ;1%%#> ,1'+(' '&# +-$2( /<$%23< 31'++ *',)11'$)-+> L'/2$- S^V [S_Y_\ ZSVmZSXB RVVU 9B K<'$(> WB @2$1'&#> k)+*%(-/$" %. <2('& ;1%%# .%$ +<-'$ $'/-+ %. `4S``>``` +-*4 S> L'/2$T`Z [S_ZY\ ZS]mZS^B RVXU WB A231)'$-11%> dB 9-0)11'> k-1%*)/" #)+/$);2/)%& '&# %/<-$ *<'$'*/-$)+/)*+ %. +/-'#" '&# ,214 +'/)1- ;1%%# 6%7 )& I&- 31'++ /2;-+B> A)%$<-%1B ] [S_]`\ ^YmS`]B RVYU WB EB K%O-1-/> J<- $<-%1%3" %. <2('& ;1%%#> A)%(-* CB 9'M-&'> J7%41'"-$-# (%#-1 %. *'++%& 62)# 6%7 /<$%23< +/-&%/)* ;1%%# 0-++-1+g ',,1)*'/)%&+ /% /<- *'$#)%0'+*21'$ +"+/-(> dB A)%(-*
TY
ACCEPTED MANUSCRIPT RV]U :B 9'&O'$> ? /7%462)# (%#-1 .%$ ,21+'/)1- 6%7 )& *'/<-/-$)D-# ;1%%# 0-++-1+> Q&/B dB L%& e)&-'$ C-*
PT
RV^U :B 9'&O'$> bB e--> J7%462)# *'++%& (%#-1 .%$ ,21+'/)1- ;1%%# 6%7 /<$%23< +/-&%+-# '$/-$)-+g ? /<-%$-/)*'1 (%#-1> K%((2&B L%&1)&-'$ 9*)B L2(-$B 9)(21B SY [T`S`\ T`^ZmT`_]B RV_U dB 9<2O1'> EB 5'$)<'$> :)+,-$+)%& /<$%23< ' 6%7)&3 62)# 7)/< 0'$)';1- 0)+*%+)/" '&# #)..2+)%& *%-.I*)-&/> 5<"+B N12)#+ TV [S_^`\ ST]SmST]VB
SC
RI
RX`U dB 9<2O1'> 9B W2,/'> :)+,-$+)%& )& ' ;)&3<'( ,1'+/)* 62)#g i..-*/+ %. ,-$),<-$'1 1'"-$> E<-%1B ?*/' TS [S_^T\ ]SVm]S_B
NU
RXSU :B ?B ?&#-$+%&> dB KB J'&&-<)11> EB !B 51-/*<-$> K%(,2/'/)%&'1 62)# (-*<'&)*+ '&# <-'/ /$'&+.-$ <-()+,<-$-> L-7 P%$O[S_^X\B RXTU EB C- EB C-$+%&> AB C*K%"> !-$()/- ,%1"&%()'1 $-,$-+-&/'/)%& %. *<$%('/%3$',<" -12/)%& *2$0-+> dB K<$%('/%3$B ? ^^ [S\ [S_]X\ SmZB
MA
RXVU :B ?B C*:%&'1#> A1%%# 6%7 )& '$/-$)-+> i#7'$# ?$&%1#g e%&#%& [S_]X\B
D
RXXU 5B L'3'$'&)> WB 9'$%H'(('> WB d'"'$'('&> i..-*/ %. ;%2&#'$" ';+%$,/)%& )& #)+,-$+)%& )& *'++%& 62)# 6%7 )& ' /2;-> ?&&B A)%(-#B i&3B VT [T``X\ ]`Zm]S_B
TE
RXYU WB d'"'$'('&> ?B e'2/)-$> AB4CB !2&3> WB d'$$"> :B e'2$-&/> L2(-$)*'1 +*<-(- .%$ (%#-11)&3 %M"3-& /$'&+.-$ )& /2;21'$ %M"3-&'/%$+> C-#B A)%1B i&3B K%(,2/B S_ [S_^S\ YTXmYVXB
AC CE P
RXZU 9B C)##1-('&> J$'&+,%$/ ,<-&%(-&' )& /<- *'$#)%0'+*21'$ +"+/-(> d%<& f)1-" r 9%&+ [S_]T\B RX]U 9B 9<'21> !B @'1('&> J<$-- ,123+ (%#-1> 5%7#-$ J-*<&%1B T^V [T`SY\ Y]_mY_TB RX^U AB :B N'0)+> dB45B K<'1).%2M> J<- -..-*/ %. 0)+*%+)/" $'/)% %& /<- (%$,<%1%3" %. ,%1",$%,"4 1-&-a,%1"*'$;%&'/- ;1-&#+ #2$)&3 ,$%*-++)&3> 5%1"(B i&3B 9*)B T] [S_^]\ SY_SmSZ``B
(Q) !RRSOL*T L3%/B&$/,+ ,6 B35,'/$A :%,E53 /+ $23 $2%3385&A3% %3D/,+ J<- #)(-&+)%&1-++ (%(-&/2( .%$( %. i8+B [S\ 4 [X\ '$- g
1 2
1 2
∂(r τc ) = 4r ∂r
0 ≤ r ≤ Ro ,
[?S\
∂(r τn ) = 4r ∂r
Ro ≤ r ≤ 1,
[?T\
τc = τy + (−m ∂uc = 0, ∂r n , τn = − ∂u ∂r
∂uc 12 ) , ∂r
τc ≥ τy .%$ Rp ≤ r ≤ Ro , τc ≤ τy .%$ 0 ≤ r ≤ Rp , . τy = 0 .%$ Ro ≤ r ≤ 1,
TZ
[?V\
ACCEPTED MANUSCRIPT c τc )+ I&)/- '&# ∂u = 0 '/ r = 0, ∂r '/ r = Ro , . τc = τn '&# uc = un '/ r = 1, un = 0
[?X\
PT
9%10)&3 /<- ;%2&#'$" 0'12- ,$%;1-( %. i8B[?S\ 4 i8B[?X\
un (r) = (1 − r2 )
Ro ≤ r ≤ 1,
- 4 , 2 Ro − r 2 −
/
[?Z\
SC
+ 2τy * 23 3 Ro − r 2 + τy (Ro − r) 3
MA
.
NU
τc = 2r
- 1 , uc (r) = 1 − Ro2 + m
[?Y\
RI
τn = 2r
[?]\
0
Rp ≤ r ≤ Ro , [?^\
2
.
,
2
Ro2 − Rp −
4
TE
1 up (r) = 1 − Ro + m ,
D
J<- ,123 6%7 0-1%*)/" *%2&/-$,'$/ up (r) *'& ;- %;/')&-# ;" /'O)&3 r = Rp )& i8B[?^\ '+g /
3+ 2τy * 32 Ro − Rp2 + τy (Ro − Rp ) 3
0
0 ≤ r ≤ Rp ,
AC CE P
[?_\ i8+B[?_> ?^> ?Z\ '$- /<- 0-1%*)/" #)+/$);2/)%&+ )& /<- /<$--41'"-$+ #%*2(-&/-# )& i8B[SZ\B
T]