Nonlinear MRI. Jeremy Goodman. CMPD/CMSO Winter School 7-12 January 2008, UCLA
|
|
- Octavia Brianne Green
- 5 years ago
- Views:
Transcription
1 Nonlinear MRI Jeremy Goodman CMPD/CMSO Winter School 7-12 January 2008, UCLA
2 Questions for the nonlinear regime How does the linear instability saturate? i.e., What nonlinear mode-mode interactions brake growth? What are the consequences for the mean flow? i.e., What does the final state look like?
3 The final state depends on what s driving the mean flow Taylor-Couette flow is contained by pressure & usually driven viscously by rotation of the boundaries. Viscous driving is weak at large Reynolds number Re=VL/ Mean flow can change drastically Disks are driven by gravity, which is stronger than magnetic & thermal energies Mean flow can t change much Dissipation is balanced instead by radial accretion
4 Taylor-Couette flow weight & pulley drive inner cylinder fluid-filled gap between cylinders motor drives outer cylinder G.I. Taylor s (1936) apparatus (B=0) Nominal parameters of the Princeton MRI experiment (B0)
5 Mechanical energy is minimized by uniform rotation Energy: Angular momentum: E = L = N p i U(r 1,...,r N ) i=1 N r i p i i=1 Minimize E at fixed L: ( E L)= 0 p i = r i p i A system of minimal energy rotates uniformly with angular velocity (=).
6 TC flow can rotate uniformly L is constant, neglecting viscous & magnetic interactions with the container and coils Centrifugal force is balanced by radial pressure gradient: 2 r = P 1 ( r c J B), if any r...in cylindrical coordinates z,r, with = e z P(r max ) P(r min ) 2 ( 2 ) 2 r max 2 r min This pressure drop is supported by the container. Absent external driving (torques), TC flow will relax to uniform rotation
7 Disks cannot rotate uniformly Boundaries are free (P0 at boundaries) except inner edge, sometimes (e.g. stellar surface) Orbital energy dominates: (r) 2 >>P/ c s 2 This means disks are thin: z<< r Centrifugal force is balanced mainly by gravity: 2 r GM * r 2 (at z << r) r 3/2 ("keplerian") Dissipation causes little change in, but a radial drift toward the star (accretion), i.e. toward lower gravitational potential
8 MRI saturation in disks Since the available free energy is so large, saturation occurs not by reduction in free energy but by secondary instabilities that feed on the main MRI modes. dynamo? MRI linear instability primary MRI modes 2 ndary instabilities These 2 ndary may include Kelvin-Helmholtz dissipation tearing modes magnetic Rayleigh-Taylor & Parker instability (if sat ~1) This has been studied mostly via analysis and simulations of a local model for a small part of a disk called the shearing box. Main goal is to determine the rate of momentum transport & dissipation at saturation ( parameter)
9 The Shearing Box: A local, corotating Cartesian approximation x r r 0 y r [ 0 0 (t)]; 0 (t) (r 0 ) z z Dimensions of the box are typically comparable to vertical scale height h c s r 0, c s sound speed
10 Shearing-box equations of motion dv dt = v t + v v = 2ẑ v 4Ax ˆx 2 zẑ 1 P + 1 c J B + 2 v d dt = v B t = ( v B )+ 2 B 2A r d dr r=r0 : shear rate 4Ax ˆx 2 zẑ : tidal field : kinematic viscosity dimensions : L 2 T 1 [ ] : magnetic diffusivity ditto
11 Standard equilibrium state of the shearing box Gas or plasma pressure dominates P >> B 2 /8, i.e. >>1 Radial gradients (of pressure, etc.) are negligible J 0 With these assumptions, 2ẑ v 4Ax ˆx v 0 = 2Axŷ 1 0 P 0 z = 2 z P 2 z+ 1 z ẑ 0 (r 0 ) + d ( ) dr r0 r r 0
12 Recap of linear theory Consider the (most important) special case, B 0 = B 0 ẑ = constant, V A B 0 4, (v,,p, B) 1 exp( t + ikz) : ( +k 2 )( +k 2 )+ kv A ( ) ( + A) +k 2 ( ) 2 + 4A( kv A ) 2 = 0 In ideal MHD (=0=), weak field is always unstable: ~ A at k ~/V A presuming d 2 /dr = 4A < 0, as in disks. Thus growth occurs on small scales (<< h) for weak fields ( >>1) For 0, modes at k > ( A /) 1/2 are stabilized. For =0 but 0, unstable range of k is unaffected but is reduced.
13 Incipient turbulence: Parasitic instabilities
14 Exact Fourier modes On scales l << c s /, perturbations are incompressible: v = v 0 + v 1, v 1 0 = 0 + 1, 1 0, etc. For a single mode ( v,p,b) 1 exp(ik r) + c.c., nonlinear terms vanish (next slide). Single Fourier modes can grow to arbitrary amplitude Assuming that boundary conditions allow single modes; The shearing box does, with shearing-periodic boundaries; TC flow does not. Though exact, sufficiently large-amplitude single modes are not stable.
15 Proof of exactness v 1 = 0 = B 1 k v 1 = 0 = k B 1 ; ( ) v 1 + c.c. = 0. v 1 v 1 v 1 ikv 1 + c.c. = ikv 1 Similarly, v 1 B 1 = 0 = B 1 v 1 = B 1 B 1. Also, 1 P 1 is linear because = 0. More generally, v 1 (r,t) = a n (t)exp(ik n r) n so v 1 v 1 ik m a n + c.c. ( ) a m 0 in general unless k m k n = 0. Modes at parallel wavenumbers can be linearly superposed provided >>1, i.e. in the incompressible limit Also true in non-ideal MHD, since 2 v 1 & 2 v 1 are linear
16 2D (x,z) simulations in ideal MHD Hawley & Balbus 1992, ApJ 400, 595 z Above: growth of magnetic energy in several runs. Below: power versus (k x,k z ) x or r v y, 1 at t=2.4, 3.1, 3.6, & 4.1 orbits. =4000, box size=0.25c s / Cartoon of evolution of the magnetic field.
17 Parasitic modes Goodman & Xu 1994, ApJ 432, 213 Consider stability of a single large-amplitude mode v 1 (r,t) = V h exp(st)sin(kz) v 0 = 2Ax ˆx, B 1 (r,t) = B h exp(st)cos(kz) B 0 = B 0 ẑ. V h, B h are constant, horizontal (xy) vectors; V h B h = 0 in ideal MHD. (1) Linearize about (1) as basic state, neglecting (v 0,B 0 ) Also neglect v because we seek 2 nd 1, B 1 ary (= parasitic ) growth rates ~ Kv 1 >> S Describe 2 nd ary mode by its lagrangian displacement : v = v 0 + v 1 + v /2, B = B 0 + B 1 + B 2 v 2 = d 2 dt t + v 1 2, and in ideal MHD, B 2 = 2 B 1 ( ) Since basic state is independent of x, y, and (by assumption) t, and periodic in z, seek a mode of Bloch form (like an electron wavefunction in a crystal): 2 X(z)exp(it + ik x x + ik y y), where X is quasiperiodic: X(z + ) = X(z)exp(ik z ) with 2 and K k K z 2, K 2
18 Parasitic modes (continued) It can be shown that in ideal MHD, the vertical displacement satisfies 1 z z z k 2 h z = 0, where k h (k x,k y ), the horizontal part of k, and (2) ( k h v 1 ) 2 (k h B 1 )2 4 = (z) Growth rates Im()=O(b), where b=b 1 /B 0 Fastest growth in 3D (k y 0) with k h B 1 =0 These are Kelvin-Helmholtz modes No perturbed magnetic tension since k h B 1 =0 slower modes with k h B 1 0 that are almost discontinuous at z where B 1 (z)0 Proto-tearing modes(?) But =0 in this analysis Non-ideal (0, 0) is feasible but has not been done! Im()/b vs. k z /K for k h B 1 =0 z vs. zk for a k h B 1 0 mode
19 Fully developed MRI turbulence Distinguish net-flux (i.e., <B z > 0) cases shearing-box boundary conditions preserve <B z > persistent instability & turbulence are guaranteed but at what level, and how does this depend on <B z >,,,...?...from zero-net-flux (<B z >=0) cases turbulence must act as a magnetic (fluctuation?) dynamo to sustain the field against reconnection & Ohmic dissipation impossible in 2D: Cowling anti-dynamo theorem even in 3D ideal MHD, persistent turbulence isn t guaranteed
20 Turbulent momentum transport: The parameter Let ij = ji turbulent stress tensor: flux of i th component of momentum in the j th direction (discounting transport by the mean flow). Dimensions: ij = [ force area]= [ energy volume]. Associated dissipation rate = ij i v j energy/volume/time [ ] In disks, this is dominated by r r r 2A xy in shearing box. xy = v x v y + B xb y 4 where v i v i v i. Reynolds stress Maxwell stress Shakura & Sunyaev (1973, before MRI was understood) postulated r ( xy )= P dimensionless viscosity parameter
21 Measured values of Numerical simulations of MRI varies with large-scale field, dissipation terms Protostellar disks based on disk masses, temperatures, accretion rates, and lifetimes Cataclysmic variables based on models of dwarf nova outbursts AGN ( QSOs, quasars, Seyferts...) direct observational constraints are few to none?
22 What we used to think we knew about MRI turbulent transport If <B z >0 ( nonzero net flux ), then 0.1 iff Re m >1, and Lu >1 *, and V A2 < c 2 s (i.e., plasma >1) Necessary only that <B z >0 within the computational domain, usually having size ~ h << r log <B 2 > time in orbits If <B z > =0 ( zero net flux ), then 0.1 iff Re m >10 4 : magnetic [energy] dynamo no external support required * Re m =VL/; Lu = V A L/ Fleming, Stone, & Hawley 2000
23 10-1 Zero net flux (closeup view) P mag /P ~ 10-4 orbits Fleming, Stone, & Hawley 2000
24 <B z > in net-flux cases scaled scaled <B z > Figure above is taken from Pessah, Chan, & Psaltis (2007), but the scaling with <B z > was noted already by Hawley, Gammie, & Balbus (1995).
25 Saturation depends on grid scale when <B z > = 0 N x =64 N x =128 N x =256 log P B t [2/] log k Level & dominant lengthscale of the turbulence decrease grid scale, not domain size. - Not what one expects from an inertial cascade. Calculations above (Fromang & Papaloizou 2007) were carried out without explicit or, hence dissipation is purely numerical.
26 increases with magnetic Prandtl number Pm = / With sufficiently large explicit viscosity () and diffusivity (), saturation should be independent of grid resolution. increases with Pm in simulations with <B z > 0 (Lesur & Longaretti 2007) Pm < in protostellar disks =0 when <B> = 0 unless Pm > 1 (Fromang, Papaloizou, Lesur, & Heinemann 2007)
27 Pm << 1 in most astrophysical & laboratory systems Protostellar disks: Pm ~ is very uncertain due to nonthermal ionization Cataclysmic-variable (white-dwarf) disks: Pm ~ 1 Quasar disks: Pm ~ [broad-line region] Liquid metals (e.g. Na, Ga): Pm ~
School and Conference on Analytical and Computational Astrophysics November, Angular momentum transport in accretion disks
2292-13 School and Conference on Analytical and Computational Astrophysics 14-25 November, 2011 Angular momentum transport in accretion disks Gianluigi Bodo Osservatorio Astronomico, Torino Italy Angular
More informationThe Magnetorotational Instability
The Magnetorotational Instability Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics March 10, 2014 These slides are based off of Balbus & Hawley (1991), Hawley
More informationNonmodal Growth and the Unstratified MRI Dynamo
MPPC general meeting, Berlin, June 24 Nonmodal Growth and the Unstratified MRI Dynamo Jonathan Squire!! and!! Amitava Bhattacharjee MRI turbulence Observed accretion rate in disks in astrophysical disks
More informationGlobal magnetorotational instability with inflow The non-linear regime
Global magnetorotational instability with inflow The non-linear regime Evy Kersalé PPARC Postdoctoral Research Associate Dept. of Appl. Math. University of Leeds Collaboration: D. Hughes & S. Tobias (Dept.
More informationAmplification of magnetic fields in core collapse
Amplification of magnetic fields in core collapse Miguel Àngel Aloy Torás, Pablo Cerdá-Durán, Thomas Janka, Ewald Müller, Martin Obergaulinger, Tomasz Rembiasz Universitat de València; Max-Planck-Institut
More informationHydrodynamic Turbulence in Accretion Disks
Hydrodynamic Turbulence in Accretion Disks Banibrata Mukhopadhyay, Niayesh Afshordi, Ramesh Narayan Harvard-Smithsonian Center for Astrophysics, 6 Garden Street, MA 38, USA Turbulent viscosity in cold
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Rough breakdown of MHD shocks Jump conditions: flux in = flux out mass flux: ρv n magnetic flux: B n Normal momentum flux: ρv n
More informationGlobal Simulations of Black Hole Accretion. John F. Hawley Department of Astronomy, University of Virginia
Global Simulations of Black Hole Accretion John F. Hawley Department of Astronomy, University of Virginia Collaborators and Acknowledgements Julian Krolik, Johns Hopkins University Scott Noble, JHU Jeremy
More informationGlobal Magnetorotational Instability with Inflow
Global Magnetorotational Instability with Inflow Evy Kersalé PPARC Postdoctoral Research Associate Dept. of Appl. Maths University of Leeds Collaboration: D. Hughes & S. Tobias (Appl. Maths, Leeds) N.
More informationviscous stress tensor: T r!
Accretion in a rotating system is only possible if matter looses its angular momentum! viscous stress tensor: T r! mass ang. mom. Viscosity could do this, but molecular viscosity is far too low. Thus one
More informationMomentum transport from magnetic reconnection in laboratory an. plasmas. Fatima Ebrahimi
Momentum transport from magnetic reconnection in laboratory and astrophysical plasmas Space Science Center - University of New Hampshire collaborators : V. Mirnov, S. Prager, D. Schnack, C. Sovinec Center
More informationThe Physics of Collisionless Accretion Flows. Eliot Quataert (UC Berkeley)
The Physics of Collisionless Accretion Flows Eliot Quataert (UC Berkeley) Accretion Disks: Physical Picture Simple Consequences of Mass, Momentum, & Energy Conservation Matter Inspirals on Approximately
More informationAccretion Disks I. High Energy Astrophysics: Accretion Disks I 1/60
Accretion Disks I References: Accretion Power in Astrophysics, J. Frank, A. King and D. Raine. High Energy Astrophysics, Vol. 2, M.S. Longair, Cambridge University Press Active Galactic Nuclei, J.H. Krolik,
More informationSmall-Scale Dynamo and the Magnetic Prandtl Number
MRI Turbulence Workshop, IAS, Princeton, 17.06.08 Small-Scale Dynamo and the Magnetic Prandtl Number Alexander Schekochihin (Imperial College) with Steve Cowley (Culham & Imperial) Greg Hammett (Princeton)
More informationStellar Magnetohydrodynamics
Stellar Magnetohydrodynamics 12 December 2006 by JJG Magnetohydrodynamics (MHD) is the generalization of hydrodynamics to electrically conducting fluids. Major astrophysical applications include stellar
More informationMagnetized Accretion Discs and the Dwarf Nova Phenomenon
and the Dwarf Nova Phenomenon 1, E. Breedt 1,2 and M. Erwee 1 1 Department of Physics, University of the Free State, PO Box 339, Bloemfontein, 9300, South Africa 2 Department of Physics and Astronomy,
More informationNIMROD simulations of dynamo experiments in cylindrical and spherical geometries. Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest,
NIMROD simulations of dynamo experiments in cylindrical and spherical geometries Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest, 1 Introduction Two experiments, Madison Plasma Couette Experiment
More informationProtoplanetary Disks: Gas Dynamics
Spring School on Planet Formation, IASTU, Beijing, 05/2014 Protoplanetary Disks: Gas Dynamics Xuening Bai Institute for Theory and Computation, Harvard-Smithsonian Center for Astrophysics Topics to be
More informationThe Physics of Fluids and Plasmas
The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the
More informationSimulations of magnetic fields in core collapse on small and large scales
Simulations of magnetic fields in core collapse on small and large scales Miguel Ángel Aloy Torás, Pablo Cerdá-Durán, Thomas Janka, Ewald Müller, Martin Obergaulinger, Tomasz Rembiasz CAMAP, Departament
More informationLarge scale magnetic fields and Dynamo theory. Roman Shcherbakov, Turbulence Discussion Group 14 Apr 2008
Large scale magnetic fields and Dynamo theory Roman Shcherbakov, Turbulence Discussion Group 14 Apr 2008 The Earth Mainly dipolar magnetic field Would decay in 20kyr if not regenerated Declination of the
More informationAttilio Ferrari. CIFS, Università di Torino. 12th Agile Workshop, May 8, 2014
Attilio Ferrari CIFS, Università di Torino 12th Agile Workshop, May 8, 2014 Plasma processes of astrophysical relevance Highly nonlinear (relativistic) physics Huge extension of physical parameters Scalability?
More informationTurbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields
Turbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields Jungyeon Cho Dmitry Shapovalov MWMF Madison, Wisconsin April 2012 The Large Scale Dynamo The accumulation of magnetic
More informationThe role of magnetic fields for planetary formation
fields for esimal Leiden Observatory, Leiden University IAU Symposium 259: Cosmic Magnetic Fields: from s, to Stars and Galaxies Puerto Santiago, Tenerife, November 2008 Star and The four stages of low
More informationTurbulent transport and its effect on the dead zone in protoplanetary discs
Astronomy & Astrophysics manuscript no. ilgner nelson4-final c ESO 2008 March 3, 2008 Turbulent transport and its effect on the dead zone in protoplanetary discs Martin Ilgner & Richard P. Nelson Astronomy
More informationSYMMETRIES, SCALING LAWS, AND CONVERGENCE IN SHEARING-BOX SIMULATIONS OF MAGNETO-ROTATIONAL INSTABILITY DRIVEN TURBULENCE
C 011. The American Astronomical Society. All rights reserved. Printed in the U.S.A. doi:10.1088/0004-637x/739//8 SYMMETRIES, SCALING LAWS, AND CONVERGENCE IN SHEARING-BOX SIMULATIONS OF MAGNETO-ROTATIONAL
More informationThe Turbulent Magnetic Prandtl Number of MHD Turbulence in Disks
The Turbulent Magnetic Prandtl Number of MHD Turbulence in Disks Xiaoyue Guan and Charles F. Gammie 1 Astronomy Department, University of Illinois, 1002 West Green St., Urbana, IL 61801, USA ABSTRACT The
More informationThe MRI in a Collisionless Plasma
The MRI in a Collisionless Plasma Eliot Quataert (UC Berkeley) Collaborators: Prateek Sharma, Greg Hammett, Jim Stone Modes of Accretion thin disk: energy radiated away (relevant to star & planet formation,
More informationTwo Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch
1 Two Fluid Dynamo and Edge-Resonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of Wisconsin-Madison, Madison,
More informationAccretion Disks Angular momentum Now let s go back to black holes. Black holes being what they are, something that falls into one disappears without
Accretion Disks Angular momentum Now let s go back to black holes. Black holes being what they are, something that falls into one disappears without a peep. It might therefore seem that accretion onto
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationVortices in accretion discs: formation process and dynamical evolution
Vortices in accretion discs: formation process and dynamical evolution Geoffroy Lesur DAMTP (Cambridge UK) LAOG (Grenoble) John Papaloizou Sijme-Jan Paardekooper Giant vortex in Naruto straight (Japan)
More informationFluctuation dynamo amplified by intermittent shear bursts
by intermittent Thanks to my collaborators: A. Busse (U. Glasgow), W.-C. Müller (TU Berlin) Dynamics Days Europe 8-12 September 2014 Mini-symposium on Nonlinear Problems in Plasma Astrophysics Introduction
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationAccretion disks. AGN-7:HR-2007 p. 1. AGN-7:HR-2007 p. 2
Accretion disks AGN-7:HR-2007 p. 1 AGN-7:HR-2007 p. 2 1 Quantitative overview Gas orbits in nearly circular fashion Each gas element has a small inward motion due to viscous torques, resulting in an outward
More informationAxisymmetric Shearing Box Models of Magnetized Disks
Axisymmetric Shearing Box Models of Magnetized Disks Xiaoyue Guan and Charles F. Gammie Astronomy Department, University of Illinois, Urbana, IL 61801 ABSTRACT The local model, or shearing box, has proven
More information36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson
36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through
More informationPAPER 57 DYNAMICS OF ASTROPHYSICAL DISCS
MATHEMATICAL TRIPOS Part III Monday, 9 June, 2014 1:30 pm to 3:30 pm PAPER 57 DYNAMICS OF ASTROPHYSICAL DISCS Attempt no more than TWO questions. There are THREE questions in total. The questions carry
More informationAXISYMMETRIC SHEARING BOX MODELS OF MAGNETIZED DISKS
The Astrophysical Journal Supplement Series, 174:145Y157, 2008 January # 2008. The American Astronomical Society. All rights reserved. Printed in U.S.A. AXISYMMETRIC SHEARING BOX MODELS OF MAGNETIZED DISKS
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationChallenges of Predictive 3D Astrophysical Simulations of Accreting Systems. John F. Hawley Department of Astronomy, University of Virginia
Challenges of Predictive 3D Astrophysical Simulations of Accreting Systems John F. Hawley Department of Astronomy, University of Virginia Stellar Evolution: Astro Computing s Early Triumph Observed Properties
More informationarxiv: v1 [astro-ph.sr] 24 Oct 2011
Astrophysics and Space Science DOI 10.1007/s - - - The role of Hall diffusion in the magnetically threaded thin accretion discs Mohsen Shadmehri 1 arxiv:1110.5142v1 [astro-ph.sr] 24 Oct 2011 c Springer-Verlag
More informationGLOBAL MAGNETOROTATIONAL INSTABILITY WITH INFLOW. II. THE NONLINEAR DEVELOPMENT OF AXISYMMETRIC WALL MODES
The Astrophysical Journal, 638:382 390, 2006 February 10 # 2006. The American Astronomical Society. All rights reserved. Printed in U.S.A. GLOBAL MAGNETOROTATIONAL INSTABILITY WITH INFLOW. II. THE NONLINEAR
More informationAccretion Disks: angular momentum Now let s go back to black holes. Black holes being what they are, something that falls into one disappears without
Accretion Disks: angular momentum Now let s go back to black holes. Black holes being what they are, something that falls into one disappears without a peep. It might therefore seem that accretion onto
More informationInterfacial Dynamics in Protoplanetary Accretion Disks BERN 2017
DEPARTMENT OF EARTH SCIENCES, RAYMOND AND BEVERLY SACKLER FACULTY OF EXACT SCIENCES Interfacial Dynamics in Protoplanetary Accretion Disks BERN 2017 Author: Ron Yellin-Bergovoy Supervisors: Prof. Eyal
More information[2 marks] Show that derivative of the angular velocity. What is the specific angular momentum j as a function of M and R in this Keplerian case?
k!! Queen Mary University of London M. Sc. EXAM I N AT1 0 N ASTMOOS Angular Momentum and Accretion in Astrophysics Fkiday, 26th May, 2006 18:15-19:45 Time Allowed: lh 30m This paper has two Sections and
More informationTheory of star formation
Theory of star formation Monday 8th 17.15 18.00 Molecular clouds and star formation: Introduction Tuesday 9th 13.15 14.00 Molecular clouds: structure, physics, and chemistry 16.00 16.45 Cloud cores: statistics
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationSimulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions
Chapter 1 Earth Science Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Project Representative Yozo Hamano Authors Ataru Sakuraba Yusuke Oishi
More informationACTIVE GALACTIC NUCLEI: FROM THE CENTRAL BLACK HOLE TO THE GALACTIC ENVIRONMENT
Julian H. Krolik ACTIVE GALACTIC NUCLEI: FROM THE CENTRAL BLACK HOLE TO THE GALACTIC ENVIRONMENT PRINCETON UNIVERSITY PRESS Princeton, New Jersey Preface Guide for Readers xv xix 1. What Are Active Galactic
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationDynamics of Astrophysical Discs
Dynamics of Astrophysical Discs 16 lectures, 3 example classes http://www.damtp.cam.ac.uk/user/hl278/dad.html Henrik Latter e-mail: hl278@cam.ac.uk 16 lectures Tu. Th. 10 Course Outline Office: F1.19 hl278@cam.
More informationThe total luminosity of a disk with the viscous dissipation rate D(R) is
Chapter 10 Advanced Accretion Disks The total luminosity of a disk with the viscous dissipation rate D(R) is L disk = 2π D(R)RdR = 1 R 2 GM Ṁ. (10.1) R The disk luminosity is half of the total accretion
More informationThe large-scale magnetic field in protoplanetary disks
The large-scale magnetic field in protoplanetary disks Jérôme Guilet MPA, Garching Max-Planck-Princeton center for plasma physics In collaboration with Gordon Ogilvie (Cambridge) 1/24 Talk outline 1) Impacts
More informationSW103: Lecture 2. Magnetohydrodynamics and MHD models
SW103: Lecture 2 Magnetohydrodynamics and MHD models Scale sizes in the Solar Terrestrial System: or why we use MagnetoHydroDynamics Sun-Earth distance = 1 Astronomical Unit (AU) 200 R Sun 20,000 R E 1
More information28. THE GRAVITATIONAL INTERCHANGE MODE, OR g-mode. We now consider the case where the magneto-fluid system is subject to a gravitational force F g
8. THE GRAVITATIONAL INTERCHANGE MODE, OR g-mode We now consider the case where the magneto-fluid system is subject to a gravitational force F g =!g, where g is a constant gravitational acceleration vector.
More informationProtoplanetary Disks. Ernst Dorfi WS 2012
Protoplanetary Disks Ernst Dorfi WS 2012 Disks and star formation Disks of gas and dust are observed around young stars for about 107 years Disks form because of angular momentum (conservation), collapse
More informationCataclysmic variables
Cataclysmic variables Sander Bus Kapteyn Astronomical Institute Groningen October 6, 2011 Overview Types of cataclysmic stars How to form a cataclysmic variable X-ray production Variation in outburst lightcurve,
More informationMagnetorotational Instability: Plans at MPA, IPP and TU Berlin
Magnetorotational Instability: Plans at MPA, IPP and TU Berlin Wolf-Christian Müller 1, Ewald Müller 2, Thomas Janka 2, Frank Jenko 3 1 Max-Planck-Institut für Plasmaphysik, Garching 2 Max-Planck-Institut
More informationSelected Topics in Plasma Astrophysics
Selected Topics in Plasma Astrophysics Range of Astrophysical Plasmas and Relevant Techniques Stellar Winds (Lecture I) Thermal, Radiation, and Magneto-Rotational Driven Winds Connections to Other Areas
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationDedalus: An Open-Source Spectral Magnetohydrodynamics Code
Dedalus: An Open-Source Spectral Magnetohydrodynamics Code Keaton Burns University of California Berkeley Jeff Oishi SLAC National Accelerator Laboratory Received ; accepted ABSTRACT We developed the magnetohydrodynamic
More informationStabilization of sawteeth in tokamaks with toroidal flows
PHYSICS OF PLASMAS VOLUME 9, NUMBER 7 JULY 2002 Stabilization of sawteeth in tokamaks with toroidal flows Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College
More informationarxiv:astro-ph/ v2 19 Nov 2006
APJ IN PRESS Preprint typeset using LATEX style emulateapj v. 05/04/06 THREE DIMENSIONAL COMPRESSIBLE HYDRODYNAMIC SIMULATIONS OF VORTICES IN DISKS YUE SHEN, JAMES M. STONE AND THOMAS A. GARDINER Department
More informationRADIATION FROM ACCRETION ONTO BLACK HOLES
RADIATION FROM ACCRETION ONTO BLACK HOLES Accretion via MHD Turbulence: Themes Replacing dimensional analysis with physics MRI stirs turbulence; correlated by orbital shear; dissipation heats gas; gas
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationEquilibrium Structure of Radiation-dominated Disk Segments
Equilibrium Structure of Radiation-dominated Disk Segments Shigenobu Hirose The Earth Simulator Center, JAMSTEC, Japan collaborators Julian Krolik (JHU) Omer Blaes (UCSB) Reconstruction of standard disk
More informationChapter 7. Basic Turbulence
Chapter 7 Basic Turbulence The universe is a highly turbulent place, and we must understand turbulence if we want to understand a lot of what s going on. Interstellar turbulence causes the twinkling of
More informationAccretion Disks Black holes being what they are, something that falls into one disappears without a peep. It might therefore seem that accretion onto
Accretion Disks Black holes being what they are, something that falls into one disappears without a peep. It might therefore seem that accretion onto a black hole would release no energy. It isn t the
More informationJet Stability: A computational survey
Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,
More informationarxiv: v2 [astro-ph.sr] 20 Nov 2018
Draft version November 21, 2018 Typeset using LATX twocolumn style in AASTeX61 ACTIV MODS AND DYNAMICAL BALANCS IN MRI-TURBULNC OF KPLRIAN DISKS WITH A NT VRTICAL MAGNTIC FILD D. Gogichaishvili, 1 G. Mamatsashvili,
More informationOverview spherical accretion
Spherical accretion - AGN generates energy by accretion, i.e., capture of ambient matter in gravitational potential of black hole -Potential energy can be released as radiation, and (some of) this can
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More informationDisc-Planet Interactions during Planet Formation
Disc-Planet Interactions during Planet Formation Richard Nelson Queen Mary, University of London Collaborators: Paul Cresswell (QMUL), Martin Ilgner (QMUL), Sebastien Fromang (DAMTP), John Papaloizou (DAMTP),
More informationarxiv: v1 [astro-ph.he] 10 Jan 2012
Accretion Disc Theory since Shakura and Sunyaev Andrew King 1 arxiv:1201.2060v1 [astro-ph.he] 10 Jan 2012 ABSTRACT I briefly review the progress of accretion disc theory since the seminal paper of Shakura
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationTurbulence and Transport The Secrets of Magnetic Confinement
Turbulence and Transport The Secrets of Magnetic Confinement Presented by Martin Greenwald MIT Plasma Science & Fusion Center IAP January 2005 FUSION REACTIONS POWER THE STARS AND PRODUCE THE ELEMENTS
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationarxiv:astro-ph/ v2 5 Aug 1997
Dissipation of a tide in a differentially rotating star Suzanne Talon Observatoire de Paris, Section de Meudon, 92195 Meudon, France and arxiv:astro-ph/9707309v2 5 Aug 1997 Pawan Kumar Institute for Advanced
More informationThe Effects of Anisotropic Transport on Dilute Astrophysical Plasmas Eliot Quataert (UC Berkeley)
The Effects of Anisotropic Transport on Dilute Astrophysical Plasmas Eliot Quataert (UC Berkeley) in collaboration with Ian Parrish, Prateek Sharma, Jim Stone, Greg Hammett Hydra A w/ Chandra Galactic
More informationFundamentals of Turbulence
Fundamentals of Turbulence Stanislav Boldyrev (University of Wisconsin - Madison) Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas What is turbulence? No exact definition.
More informationFrancesco Califano. Physics Department, University of Pisa. The role of the magnetic field in the interaction of the solar wind with a magnetosphere
Francesco Califano Physics Department, University of Pisa The role of the magnetic field in the interaction of the solar wind with a magnetosphere Collaboration with M. Faganello & F. Pegoraro Vien na,
More information3D Reconnection of Weakly Stochastic Magnetic Field and its Implications
3D Reconnection of Weakly Stochastic Magnetic Field and its Implications Alex Lazarian Astronomy Department and Center for Magnetic Self- Organization in Astrophysical and Laboratory Plasmas Collaboration:
More informationHigh-Energy Astrophysics
Oxford Physics: Part C Major Option Astrophysics High-Energy Astrophysics Garret Cotter garret@astro.ox.ac.uk Office 756 DWB Michaelmas 2011 Lecture 7 Today s lecture: Accretion Discs Part I The Eddington
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More informationHartmann Flow in a Rotating System in the Presence of Inclined Magnetic Field with Hall Effects
Tamkang Journal of Science and Engineering, Vol. 13, No. 3, pp. 243 252 (2010) 243 Hartmann Flow in a Rotating System in the Presence of Inclined Magnetic Field with Hall Effects G. S. Seth, Raj Nandkeolyar*
More informationResistive MHD, reconnection and resistive tearing modes
DRAFT 1 Resistive MHD, reconnection and resistive tearing modes Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 6 May 18 1. Introduction
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationLecture 14: Viscous Accretion Disks. HST: μm image of the A0e star AB Aurigae
Lecture 14: Viscous Accretion Disks HST: 0.2-1 μm image of the A0e star AB Aurigae OB1a (10 Myr) Subgroups in the Orion OB1 Associations OB1b (5 Myr) OB1c (2.5 Myr) OB1d (1 Myr) The Orion OB1 association
More informationIdeal Magnetohydrodynamics (MHD)
Ideal Magnetohydrodynamics (MHD) Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 1, 2016 These lecture notes are largely based on Lectures in Magnetohydrodynamics
More informationarxiv: v1 [astro-ph.sr] 5 May 2016
Under consideration for publication in J. Fluid Mech. 1 arxiv:1605.01462v1 [astro-ph.sr] 5 May 2016 Linear and nonlinear responses to harmonic force in rotating flow Xing Wei 1,2 1 Institute of Natural
More informationPlasma Physics for Astrophysics
- ' ' * ' Plasma Physics for Astrophysics RUSSELL M. KULSRUD PRINCETON UNIVERSITY E;RESS '. ' PRINCETON AND OXFORD,, ', V. List of Figures Foreword by John N. Bahcall Preface Chapter 1. Introduction 1
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationThe Turbulent Magnetic Prandtl Number of MHD Turbulence in Disks
The Turbulent Magnetic Prandtl Number of MHD Turbulence in Disks Xiaoyue Guan and Charles F. Gammie 1 arxiv:0903.3757v1 [astro-ph.he] 22 Mar 2009 Astronomy Department, University of Illinois, 1002 West
More informationChapter 5 Thin Accretion Disks
Chapter 5 Thin Accretion Disks 1. Disk formation with oche lobe overflow In a non-rotating frame, the components of stream velocity at the L 1 point along and perpendicular to the instantaneous line of
More informationarxiv:astro-ph/ v1 31 Aug 1998
A&A manuscript no. (will be inserted by hand later) Your thesaurus codes are: (3.3.4;..; 7..; 7.9.; 8..3) ASTRONOMY AND ASTROPHYSICS On the treatment of the Coriolis force in computational astrophysics
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationEART164: PLANETARY ATMOSPHERES
EART164: PLANETARY ATMOSPHERES Francis Nimmo Last Week - Dynamics Reynolds number, turbulent vs. laminar flow Velocity fluctuations, Kolmogorov cascade Brunt-Vaisala frequency, gravity waves Rossby waves,
More information