Accretion Disk Matt Coleman Institute for Advanced Study Boundary - - PowerPoint PPT Presentation

accretion disk
SMART_READER_LITE
LIVE PREVIEW

Accretion Disk Matt Coleman Institute for Advanced Study Boundary - - PowerPoint PPT Presentation

Accretion Disk Matt Coleman Institute for Advanced Study Boundary Layers mcoleman@ias.edu UNLV - Athena++ User Meeting - 3/20/19 The Boundary Layer P disk, ram > P , mag Where the disk meets the star MRI doesnt work How


slide-1
SLIDE 1

UNLV - Athena++ User Meeting - 3/20/19

Accretion Disk Boundary Layers

Matt Coleman Institute for Advanced Study mcoleman@ias.edu

slide-2
SLIDE 2

The Boundary Layer

  • Where the disk meets the star
  • MRI doesn’t work
  • How does angular momentum transport here work?

ℳ = vkep (R⋆) Cs ≫ 1

Mach number

Pdisk, ram > P⋆, mag

slide-3
SLIDE 3

The Boundary Layer

ℳ = ( 8π 3 G3σm4

u

k4

b

)

1/8

M3/2μ1/2 · M−1/8R−1/8τ−1/8 (CV hot) = 32 ( M 0.6M⊙ )

3/8

( μ 0.6)

1/2

( · M−9 R 9 Mm τ 104 )

−1/8

(CV cold) = 270 ( M 0.6M⊙ )

3/8

( μ 2 )

1/2

( · M−11 R 9 Mm τ 5 )

−1/8

(AM CVn hot) = 55 ( M 1.1M⊙ )

3/8

( μ 1.4)

1/2

( · M−9 R 4.7 Mm τ 5 × 104 )

−1/8

(AM CVn cold) = 330 ( M 1.1M⊙ )

3/8

( μ 4 )

1/2

( · M−12 R 4.7 Mm τ 2000 )

−1/8

(PPD) = 53 ( M M⊙ )

3/8

( μ 2 )

1/2

( · M−8 R 2R⊙ τ 75 )

−1/8

(FU Ori) = 5.3 ( M M⊙ )

3/8

( μ 0.6)

1/2

( · M−5 R 2R⊙ τ 6 × 105)

−1/8

Numeric considerations

h⋆ ∼ ℳ−2, dt ∼ h⋆ ℳ ∼ ℳ−3

slide-4
SLIDE 4

Supersonic Shear

Hydrodynamically Unstable

slide-5
SLIDE 5

Possible DNO

slide-6
SLIDE 6

One-armed Spiral

ϕ = ϕ0 − sign (r − rp) 2 R⋆ r + r R⋆ ( R⋆ rp )

3/2

− 3 R⋆ rp

slide-7
SLIDE 7

Azimuthal Wave Modes

slide-8
SLIDE 8
slide-9
SLIDE 9

Angular Momentum Transport

1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 −2 −1 CS ×10−5

CS 10 5 11 16 4 sum

1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 −0.5 0.0 ×10−4

CL CA CS

m9-0.90 t/2π = 299.9 − 399.9

CS = Stress Transport CA = Advective Transport CL = CS + CA = Total Transport

slide-10
SLIDE 10

2D Runs

Sim Class M Runs done Nr Nphi 2DM5 5 2 1024 1024 2DM6 6 4 1024 1024 2DM6HR 6 6 2048 2048 2DM7 7 1 2048 2048 2DM8 8 1 2048 2048 2DM9LR 9 7 2048 2048 2DM9 9 7 4096 4096 2DM9HR 9 1 8192 8192 2DM10 10 1 4096 4096 2DM11 11 1 4096 4096 2DM12 12 5 4096 4096 2DM13 13 1 8192 8192 2DM14 14 1 8192 8192 2DM15 15 1 8192 8192

slide-11
SLIDE 11

Dispersion Relation

Curves are single parameter fit.

5 10 15 20 25 30 m 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 Ωp

M=5, r0 =0.82 M=6, r0 =0.85 M=7, r0 =0.86 M=8, r0 =0.88 M=9, r0 =0.89

Ωp Ω (R⋆) = ℳ−2 + ( ℳR⋆ 2mr0 )

2

slide-12
SLIDE 12

3D Preliminary Results

slide-13
SLIDE 13

3D Preliminary Results

slide-14
SLIDE 14

Accretion Disk Boundary Layers

Summary

  • Stable to MRI
  • Supersonic shear AM Transport
  • Possible explanation of DNO
  • Fourer decomposition
  • Working on 3D

Matt Coleman Institute for Advanced Study mcoleman@ias.edu