![Page 1: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/1.jpg)
BGUWISAP - 20111
Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear
Theory
Edward Liverts
Michael Mond
Yuri Shtemler
![Page 2: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/2.jpg)
BGUWISAP - 20112
Accretion disks – A classical astrophysical problem
Original theory is for infinite axially uniform cylinders. Recent studies indicate stabilizing effects of finite small thickness. Coppi and Keyes, 2003 ApJ., Liverts and Mond, 2008 MNRAS have considered modes contained within the thin disk.
Weakly nonlinear studies indicate strong saturation of the instability. Umurhan et al., 2007, PRL.
Criticism of the shearing box model for the numerical simulations and its adequacy to full disk dynamics description. King et al., 2007 MNRAS, Regev & Umurhan, 2008 A&A.
Magneto rotational instability (MRI) Velikhov, 1959, Chandrasekhar, 1960 – a major player in the transition to turbulence scenario, Balbus and Hawley, 1992.
![Page 3: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/3.jpg)
BGUWISAP - 20113
Algebraic Vs. Spectral Growthspectral analysis
ni tntf t f e ( , ) ( )
r r
Spectral stability if nIm 0
- Eigenvalues of the linear operatorn
for all s n
Short time behavior
1000eR
f
f tG t Max
f(0)
( )( )
(0)
Transient growth for Poiseuille flow.
Schmid and Henningson, “Stability and Transition in Shear Flows”
![Page 4: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/4.jpg)
BGUWISAP - 20114
To investigate the linear as well as nonlinear development of small perturbations in realistic thin disk geometry.
Goals of the current research
Derive a set of reduced model nonlinear equations
Follow the nonlinear evolution.
Have pure hydrodynamic activities been thrown out too hastily ?
Structure of the spectrum in stratified (finite) disk. Generalization to include short time dynamics. Serve as building blocks to nonlinear analysis.
![Page 5: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/5.jpg)
BGUWISAP - 20115
Magneto Rotational Instability (MRI)
Governing equations
.00,v
,4
,v-v
,1v 22
Bt
BcjBBdtBd
Bjc
rcdtd s
tzbbzB ,,v,v,
![Page 6: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/6.jpg)
BGUWISAP - 20116
z
R0
H0
Thin Disk GeometryVertical stretching
z
0
0
HR
1z
Supersonic rotation:
1sM
r
![Page 7: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/7.jpg)
BGUWISAP - 20117
2 2
2
d Pdt n
V j B
Vertically Isothermal Disks
( , ) ( , ) ( ) P r n r T r
*2*
4 PB
Vertical force balance:2 2/2 ( )( , ) ( ) H rn r N r e
Radial force balance: 3/2( ) ( ) , (r)=V r r r r
Free functions: The thickness is determined by .
( )T r( )H r ( ) ( ) / ( )H r T r r
( )N r
![Page 8: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/8.jpg)
BGUWISAP - 2011
• I. Dominant axial component • II. Dominant toroidal component
Two magnetic configurations
0 ( , )
( , )z zB B r
B B r
0 0
( , )
( , )
z zB B r
B B r
( )B r
Magnetic field does not influence lowest order equilibrium but makes a significant difference in perturbations’ dynamics.
Free function: ( )zB r Free function: ( )B r
8
![Page 9: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/9.jpg)
BGUWISAP - 20119
Linear Equations: Case Iz
zBzB2 2/2 ( )( , ) ( ) H rn r N r e
Self-similar variable
( )
H r
0B Consider
2 2
2( ) ( ) ( )( )
( )
z
N r H r rrB r
2 2 2( ) ( ) ( )sH r r C r
![Page 10: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/10.jpg)
BGUWISAP - 20111010
2
( ) 12 ( ) 0( ) ( )
1 ( ) 1 02 ( ) ( )
( ) 0
ln( ) 0ln
r r
r
r r
r
V Brr Vt r n
V BrVt r n
B Vrt
B V dr Bt d r
The In-Plane Coriolis-Alfvén System
2 /2( ) n e
2 2
2( ) ( ) ( )( )
( )
z
N r H r rrB r
No radial derivatives Only spectral instability possible
Boundary conditions
, 0 rB
![Page 11: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/11.jpg)
BGUWISAP - 201111
The Coriolis-Alfvén System
n2 /2( ) n e
2( ) sec ( ) n h b
22 /2sec ( ) h b d e d
2
b
Approximate profile
Change of independent variable
tanh( ) b
1
2
Spitzer, (1942)
![Page 12: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/12.jpg)
BGUWISAP - 201112
The Coriolis-Alfvén Spectrum
,ˆ( )( ) 0 rK K V L L
2[(1 ) ]
d dd d
L Legendre operator
2 22
( ) 3 2 9 162 rK
b
ˆ( , , ) ( , ) , ( ) ( )f r t f r e r r i t
,ˆ ( ) r kV P ( 1) , 1,2,K k k k
3
![Page 13: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/13.jpg)
BGUWISAP - 2011
The Coriolis-Alfvén Spectrum
13
![Page 14: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/14.jpg)
BGUWISAP - 201114
The Coriolis-Alfvén Dispersion Relation
2 22
( ) 3 2 9 162 rK
b ( 1)k k
2 4 2 ( )ˆ ˆ ˆ ˆ ˆ( 6) 9(1 ) 0 , cr
r
2( 1)
3cr
b k k
Universal condition for instability:
ˆ 1 cr
Number of unstable modes increases with .
![Page 15: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/15.jpg)
BGUWISAP - 201115
The Coriolis-Alfvén Dispersion Relation-The MRI
Im , Re Im
( 1) 0.42cr k max ˆIm for 3
Im
![Page 16: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/16.jpg)
BGUWISAP - 2011
( ) 0( )
( ) ( ) 0
z
z
V nrt n
n r n Vt
ˆ ˆ( ) ( 1) 0n n
The Vertical Magnetosonic System
2 22
2 2
ˆ ˆ(1 ) 2 01
d n nd
ˆ( , , ) ( , )( ) ( )
f r t f r e
r r
i t
2( ) sec ( ) n h b
tanh( ) b
16
![Page 17: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/17.jpg)
BGUWISAP - 201117
Continuous acoustic spectrum
21 1 2 2( ) 1 ( ) ( )n c f c f
21
1( ) ( )1
f
22
1( ) ( )1
f
2 21
Define: 2( ) 1 ( )n g
2 2 22
2 2 21(1 ) 2 2 01
d g d g gd d
Associated Legendre
equation
![Page 18: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/18.jpg)
BGUWISAP - 201118
Continuous eigenfunctions
1( )f
2 ( )f
![Page 19: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/19.jpg)
BGUWISAP - 201119
Linear Equations: Case II z
B
( , )
( , )z zB B r
B B r
2 2/2 ( )( , ) ( ) H rn r N r e
Axisymmetric perturbations are spectrally stable.Examine non-axisymmetric perturbations.
( ) , - , , ( ) ( )
r drr tH r H r
![Page 20: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/20.jpg)
BGUWISAP - 201120
Two main regimes of perturbations dynamics
I. Decoupled inertia-coriolisand MS modes
II. Coupled inertia-coriolisand MS modes
A. Inertia-coriolis waves driven algebraically by acoustic modes
B. Acoustic waves driven algebraically by inertia-coriolis modes
![Page 21: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/21.jpg)
BGUWISAP - 201121
I.A – Inertial-coriolis waves driven by MS modes
1 1 0( ) ( )
[ ( ) ]1 0( )
z
z
z
bV nn
n Vnn
b V
2
2( ) ( )( )
( )sN c
B
The magnetosonic eigenmodes
ˆ( , , , ) ( , ) i ikf f e
22
2
ˆ ˆ( ) 1 ( ) ˆ( 2) 0( ) 1 ( ) 1
d b n db n bn d nd
![Page 22: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/22.jpg)
BGUWISAP - 201122
WKB solution for the eigenfunctions
Turning points at *
The Bohr-Sommerfeld condition
*2 2
*
1( ) ( 1) , ( ) ( * )2 4
d m
1MS m ˆ( , , , ) ( , ) i ikf f e
![Page 23: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/23.jpg)
BGUWISAP - 201123
The driven inertia-coriolis (IC) waves
1 0
0
2 { , , } { , , }
1 { , , }2
rz z
r z
V V L n V b L n V b
V V K n V b
IC operator Driving MS modes
2
ˆ[ ( ) ]1 1 1ˆˆ( ) ( ) exp[ ]( ) ( )1
MS zr MS
MS
i d vv D ik b i ikd
1( )
rr
dvbdS
ln( ) dDd
![Page 24: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/24.jpg)
BGUWISAP - 201124
II – Coupled regime
Need to go to high axial and azimuthal wave numbers1, zk k
22
212 2 ( ) ( ) {[ ( ) 1] ( ) }r z r
r r r r r zd b db dbk k k k b k k b
d dd
22
21 [(1 ) ( ) ]z
z z r rd b k k b k k bd
3( , ) ( )2r r rk k D k
![Page 25: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/25.jpg)
BGUWISAP - 201125
Instability
![Page 26: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/26.jpg)
BGUWISAP - 201126
Hydrodynamic limit
Both regimes have clear hydrodynamic limits. Similar algebraic growth for decoupled regime.Much reduced growth is obtained for the coupled regime.
Consistent with the following hydrodynamic calculations:Umurhan et al. 2008Rebusco et al. 2010Shtemler et al. 2011
![Page 27: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/27.jpg)
BGUWISAP - 201127
II – Nonlinear theory
),()1((4
1)1( 22
1
Ovv
vbb
bvv
Lvv
tr
zr
zr
Ar
),()1( 22
Obb
vvv
bbb
Lbb
tr
zr
zr
Ar
),(1
)1(81
122
22
O
vbbv
Lv
t z
rzs
z
![Page 28: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/28.jpg)
BGUWISAP - 201128
Weakly nonlinear asymptotic model
.)(,,32
1427
cc i
~A
Ginzburg-Landau equation3AA
tA
Landau, (1948)
Phase plane
![Page 29: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/29.jpg)
BGUWISAP - 201129
Numerical Solution
...)()()()()()(),( 3)3(
,2)2(
,1)1(,, yPtvyPtvyPtvytv zzzz
![Page 30: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/30.jpg)
BGUWISAP - 201130
Weakly nonlinear theory, range I
.)(,,32,04])3(2][2[ 14
27222cc i
)].()([)()(
)(),(1
)()( 2021874
32
176
yPyPtAbb
yPtAvv
c
crcr
)561)(()1()(
42
2
tCtBvz
).(10
),(94
),(358
3
32
4
AOBdtdC
AOACdtdB
AOABAdtdA
nmmn ndPP
1
1 122)()(
![Page 31: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/31.jpg)
BGUWISAP - 201131
Nonlinear theory, range II
)()()()(),()()()(
4422
3311
PtCPtCPtBPtBvz
).(720
),(56
51
),(1335
52245
104315
),(1749
52105
104385
),()175
8758(
33
4
331
2
3242
3
3242
1
321
AOBdtdC
AOBBdtdC
AOACCdtdB
AOACCdtdB
AOBBAAdtdA
![Page 32: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/32.jpg)
BGUWISAP - 201132
Generation of toroidal magnetic field
)]()()[( 20 PPtAb
...)()()()()()(),( 3)3(
2)2(
1)1( yPtbyPtbyPtbytb
rtdrdbb r
![Page 33: Spectral and Algebraic Instabilities in Thin Keplerian Disks: I – Linear Theory](https://reader036.vdocument.in/reader036/viewer/2022070421/56816354550346895dd3f95b/html5/thumbnails/33.jpg)
BGUWISAP - 201133
Conclusions Generation of toroidal magnetic field