Doctor of Philosophy / We investigate the linear growth and vertical structure of the magnetorotational instability (MRI) in weakly ionised, stratified accretion discs. The magnetic field is initially vertical and perturbations have vertical wavevectors only. Solutions are obtained at representative radial locations from the central protostar for different choices of the initial magnetic field strength, sources of ionisation, disc structure and configuration of the conductivity tensor. The MRI is active over a wide range of magnetic field strengths and fluid conditions in low conductivity discs. For the minimum-mass solar nebula model, incorporating cosmic ray and x-ray ionisation and assuming that charges are carried by ions and electrons only, perturbations grow at 1 AU for B < 8G. For a significant subset of these strengths (200mG < B < 5 G), the growth rate is of order the ideal MHD rate (0.75 Omega). Hall conductivity modifies the structure and growth rate of global unstable modes at 1 AU for all magnetic field strengths that support MRI. As a result, at this radius, modes obtained with a full conductivity tensor grow faster and are active over a more extended cross-section of the disc, than perturbations in the ambipolar diffusion limit. For relatively strong fields (e.g. B > 200 mG), ambipolar diffusion alters the envelope shapes of the unstable modes, which peak at an intermediate height, instead of being mostly flat as modes in the Hall limit are in this region of parameter space. Similarly, when cosmic rays are assumed to be excluded from the disc by the winds emitted by the magnetically active protostar, unstable modes grow at this radius for B < 2 G. For strong fields, perturbations exhibit a kink at the height where x-ray ionisation becomes active. Finally, for R = 5 AU (10 AU), unstable modes exist for B < 800 mG (B < 250 mG) and the maximum growth rate is close to the ideal-MHD rate for 20 mG < B < 500 mG (2 mG < B < 50 mG). Similarly, perturbations incorporating Hall conductivity have a higher wavenumber and grow faster than solutions in the ambipolar diffusion limit for B < 100 mG (B < 10 mG). Unstable modes grow even at the midplane for B > 100 mG (B ~ 1 mG), but for weaker fields, a small dead region exists. When a population of 0.1 um grains is assumed to be present, perturbations grow at 10 AU for B < 10 mG. We estimate that the figure for R = 1 AU would be of order 400 mG. We conclude that, despite the low magnetic coupling, the magnetic field is dynamically important for a large range of fluid conditions and field strengths in protostellar discs. An example of such magnetic activity is the generation of MRI unstable modes, which are supported at 1 AU for field strengths up to a few gauss. Hall diffusion largely determines the structure and growth rate of these perturbations for all studied radii. At radii of order 1 AU, in particular, it is crucial to incorporate the full conductivity tensor in the analysis of this instability, and more generally, in studies of the dynamics of astrophysical discs.

Magnetohydrodynamic (MHD) flows in annular channels are of great current interest due to experimental search for the so-called magnetorotational instability (MRI) which is important for astrophysical applications (accretion disk physics, magnetic dynamo effect). <p>The main point of MRI experiments is to study the stability of liquid metal rotating in an external magnetic field. Two different types of fluid rotation are proposed: Taylor-Couette flow between rotating coaxial cylinders and electrically driven flow in transverse magnetic field. The implementation of MRI experiments and explanation of experimental results requires a theoretical study of the equilibrium and the stability of MHD flow in an annular channel. This is one of the main tasks of present thesis.<p>For study of equilibrium Taylor-Couette and electrically driven flows, a numerical code is developed which is based on the finite difference scheme with Jacobi iterations. The structure of flows is calculated for different parameters of the experiment. Effect of the inertia on the rotation profiles is investigated in detail. The approximate analytical expressions are obtained for radial profiles of rotation that can be used for optimization of the experimental device for MRI investigation. Equilibrium Taylor-Couette and electrically driven flows are compared from the perspective of experimental studies of MRI.<p>The spectral stability of electrically driven flow is studied by solving the eigen-value problem. This study is performed in the frames of both ideal and dissipative MHD models. It is shown that electrically driven flow can be destabilized through the mechanism of MRI if fluid velocity exceeds some instability threshold, which is determined by non-axisymmetric modes. The obtained results are compared with available experimental data.<p>A general variational method is developed for the stability study of MHD flows of ideal compressible fluids. It is shown that the linearized dynamics of such fluids has an infinite set of invariants. A necessary and sufficient stability criterion can be obtained after inclusion of one or several such invariants in analysis. An analytical example is presented to confirm the fruitfulness of the developed method.

Magnetohydrodynamic (MHD) flows in annular channels are of great current interest due to experimental search for the so-called magnetorotational instability (MRI) which is important for astrophysical applications (accretion disk physics, magnetic dynamo effect). <p>The main point of MRI experiments is to study the stability of liquid metal rotating in an external magnetic field. Two different types of fluid rotation are proposed: Taylor-Couette flow between rotating coaxial cylinders and electrically driven flow in transverse magnetic field. The implementation of MRI experiments and explanation of experimental results requires a theoretical study of the equilibrium and the stability of MHD flow in an annular channel. This is one of the main tasks of present thesis.<p>For study of equilibrium Taylor-Couette and electrically driven flows, a numerical code is developed which is based on the finite difference scheme with Jacobi iterations. The structure of flows is calculated for different parameters of the experiment. Effect of the inertia on the rotation profiles is investigated in detail. The approximate analytical expressions are obtained for radial profiles of rotation that can be used for optimization of the experimental device for MRI investigation. Equilibrium Taylor-Couette and electrically driven flows are compared from the perspective of experimental studies of MRI.<p>The spectral stability of electrically driven flow is studied by solving the eigen-value problem. This study is performed in the frames of both ideal and dissipative MHD models. It is shown that electrically driven flow can be destabilized through the mechanism of MRI if fluid velocity exceeds some instability threshold, which is determined by non-axisymmetric modes. The obtained results are compared with available experimental data.<p>A general variational method is developed for the stability study of MHD flows of ideal compressible fluids. It is shown that the linearized dynamics of such fluids has an infinite set of invariants. A necessary and sufficient stability criterion can be obtained after inclusion of one or several such invariants in analysis. An analytical example is presented to confirm the fruitfulness of the developed method.

It is currently believed that angular momentum transport in accretion disks is mediated by magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability (MRI). More than 15 years after its discovery, an accretion disk model that incorporates the MRI as the mechanism driving the MHD turbulence is still lacking. This dissertation constitutes the first in a series of steps towards establishing the formalism and methodology needed to move beyond the standard accretion disk model and incorporating the MRI as the mechanism enabling the accretion process. I begin by presenting a local linear stability analysis of a compressible, differentially rotating flow and addressing the evolution of the MRI beyond the weak-field limit when magnetic tension forces due to strong toroidal fields are considered. Then, I derive the first formal analytical proof showing that, during the exponential growth of the instability, the mean total stress produced by correlated MHD fluctuations is positive and leads to a net outward flux of angular momentum. I also show that some characteristics of the MHD stresses that are determined during this initial phase are roughly preserved in the turbulent saturated state observed in local numerical simulations. Motivated by these results, I present the first mean-field MHD model for angular momentum transport driven by the MRI that is able to account for a number of correlations among stresses found in local numerical simulations. I point out the relevance of a new type of correlation that couples the dynamical evolution of the Reynolds and Maxwell stresses and plays a key role in developing and sustaining the MHD turbulence. Finally, I address how the turbulent transport of angular momentum depends on the magnitude of the local shear. I show that turbulent MHD stresses in accretion disks cannot be described in terms of shear-viscosity.

Doctor of Philosophy / We investigate the linear growth and vertical structure of the magnetorotational instability (MRI) in weakly ionised, stratified accretion discs. The magnetic field is initially vertical and perturbations have vertical wavevectors only. Solutions are obtained at representative radial locations from the central protostar for different choices of the initial magnetic field strength, sources of ionisation, disc structure and configuration of the conductivity tensor. The MRI is active over a wide range of magnetic field strengths and fluid conditions in low conductivity discs. For the minimum-mass solar nebula model, incorporating cosmic ray and x-ray ionisation and assuming that charges are carried by ions and electrons only, perturbations grow at 1 AU for B < 8G. For a significant subset of these strengths (200mG < B < 5 G), the growth rate is of order the ideal MHD rate (0.75 Omega). Hall conductivity modifies the structure and growth rate of global unstable modes at 1 AU for all magnetic field strengths that support MRI. As a result, at this radius, modes obtained with a full conductivity tensor grow faster and are active over a more extended cross-section of the disc, than perturbations in the ambipolar diffusion limit. For relatively strong fields (e.g. B > 200 mG), ambipolar diffusion alters the envelope shapes of the unstable modes, which peak at an intermediate height, instead of being mostly flat as modes in the Hall limit are in this region of parameter space. Similarly, when cosmic rays are assumed to be excluded from the disc by the winds emitted by the magnetically active protostar, unstable modes grow at this radius for B < 2 G. For strong fields, perturbations exhibit a kink at the height where x-ray ionisation becomes active. Finally, for R = 5 AU (10 AU), unstable modes exist for B < 800 mG (B < 250 mG) and the maximum growth rate is close to the ideal-MHD rate for 20 mG < B < 500 mG (2 mG < B < 50 mG). Similarly, perturbations incorporating Hall conductivity have a higher wavenumber and grow faster than solutions in the ambipolar diffusion limit for B < 100 mG (B < 10 mG). Unstable modes grow even at the midplane for B > 100 mG (B ~ 1 mG), but for weaker fields, a small dead region exists. When a population of 0.1 um grains is assumed to be present, perturbations grow at 10 AU for B < 10 mG. We estimate that the figure for R = 1 AU would be of order 400 mG. We conclude that, despite the low magnetic coupling, the magnetic field is dynamically important for a large range of fluid conditions and field strengths in protostellar discs. An example of such magnetic activity is the generation of MRI unstable modes, which are supported at 1 AU for field strengths up to a few gauss. Hall diffusion largely determines the structure and growth rate of these perturbations for all studied radii. At radii of order 1 AU, in particular, it is crucial to incorporate the full conductivity tensor in the analysis of this instability, and more generally, in studies of the dynamics of astrophysical discs.

abstract: Understanding the temperature structure of protoplanetary disks (PPDs) is paramount to modeling disk evolution and future planet formation. PPDs around T Tauri stars have two primary heating sources, protostellar irradiation, which depends on the flaring of the disk, and accretional heating as viscous coupling between annuli dissipate energy. I have written a "1.5-D" radiative transfer code to calculate disk temperatures assuming hydrostatic and radiative equilibrium. The model solves for the temperature at all locations simultaneously using Rybicki's method, converges rapidly at high optical depth, and retains full frequency dependence. The likely cause of accretional heating in PPDs is the magnetorotational instability (MRI), which acts where gas ionization is sufficiently high for gas to couple to the magnetic field. This will occur in surface layers of the disk, leaving the interior portions of the disk inactive ("dead zone"). I calculate temperatures in PPDs undergoing such "layered accretion." Since the accretional heating is concentrated far from the midplane, temperatures in the disk's interior are lower than in PPDs modeled with vertically uniform accretion. The method is used to study for the first time disks evolving via the magnetorotational instability, which operates primarily in surface layers. I find that temperatures in layered accretion disks do not significantly differ from those of "passive disks," where no accretional heating exists. Emergent spectra are insensitive to active layer thickness, making it difficult to observationally identify disks undergoing layered vs. uniform accretion. I also calculate the ionization chemistry in PPDs, using an ionization network including multiple charge states of dust grains. Combined with a criterion for the onset of the MRI, I calculate where the MRI can be initiated and the extent of dead zones in PPDs. After accounting for feedback between temperature and active layer thickness, I find the surface density of the actively accreting layers falls rapidly with distance from the protostar, leading to a net outward flow of mass from ~0.1 to 3 AU. The clearing out of the innermost zones is possibly consistent with the observed behavior of recently discovered "transition disks." / Dissertation/Thesis / Ph.D. Physics 2012

Cette thèse est consacrée à la vérification de l'analogie entre l'instabilité viscoélastique (VEI) et l'instabilité magnéto-rotationnel (MRI) dans un écoulement képlérien, afin de mieux comprendre le transport du moment dans les disques d'accrétion. Le discriminant de Rayleigh élasto-rotationnel est établi pour clarifier le rôle de l'élasticité dans le VEI. L'analyse de stabilité linéaire (LSA) avec le modèle d’Oldroyd-B est effectuée pour prédire les paramètres critiques des modes viscoélastiques. Il fait apparaître également l'influence de l'élasticité, la viscosité polymérique et d'autres paramètres de contrôle pour le VEI. Des expériences bien contrôlées avec des solutions aqueuses de polyoxyéthylène (POE) et de polyéthylène glycol (PEG) sont effectuées. Nous avons observé le mode stationnaire axisymétrique supercritique avec des solutions de faible élasticité et modes désordonnés sous-critiques avec des solutions de grande élasticité. Les formes et les valeurs critiques de ces modes sont en bon accord avec les prédictions théoriques de LSA. Selon l'analogie, le mode axisymétrique stationnaire est probablement l'analogue de MRI standard, tandis que le mode désordonné est probable que l'analogue de MRI hélicoïdale. La thèse contient aussi des résultats théoriques expérimentaux sur quatre autres régimes de rotation et un cas de limite d'élasticité infinie. / This thesis is devoted to the verification of the analogy between the viscoelastic instability (VEI) and the magnetorotational instability (MRI) in a Keplerian flow, in order to get better understanding of the momentum transportation in accretion disks.The elasto-rotational Rayleigh discriminant is deduced to clarify the role of the elasticity in the VEI. The linear stability analysis (LSA) with Oldroyd-B model is performed to predict critical parameters of viscoelastic modes, and it reveals the influence of the elasticity, polymer viscosity on the VEI. Experiments with well controlled aqueous solutions of polyoxyethylene (POE) and polyethylene glycol (PEG) are conducted. We have observed supercritical stationary axisymmetric mode with solutions of small elasticity and subcritical disordered modes with solutions of large elasticity. Both the flow patterns and the critical values of these modes are in good agreement with the LSA predictions. According to the analogy, the stationary axisymmetric mode is likely the analog of the standard MRI while the disordered mode is likely the analog of the helical MRI. The thesis contains also theoretical and experimental results with four other rotation regimes and the limit case of infinite elasticity.