1. Angular momentum transport and large eddy simulations in magnetorotational turbulence: the small Pm limit
- Author
-
Pierre-Yves Longaretti, S. Fromang, Geoffroy Lesur, Heloise Meheut, M. Joos, Astrophysique Interprétation Modélisation (AIM (UMR7158 / UMR_E_9005 / UM_112)), Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris Diderot - Paris 7 (UPD7)-Centre National de la Recherche Scientifique (CNRS), Institut de Planétologie et d'Astrophysique de Grenoble (IPAG ), Observatoire des Sciences de l'Univers de Grenoble (OSUG), Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Centre National de la Recherche Scientifique (CNRS)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Centre National de la Recherche Scientifique (CNRS)-Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS), The simulations presented in this paper were granted access to the HPC resources of CCRT under the allocation x2012042231 made by GENCI (Grand Équipement National de Calcul Intensif), European Project: 258729,EC:FP7:ERC,ERC-2010-StG_20091028,PETADISK(2011), European Project: 294110,EC:FP7:PEOPLE,FP7-PEOPLE-2011-CIG,HALLDISCS(2012), Astrophysique Interprétation Modélisation (AIM (UMR_7158 / UMR_E_9005 / UM_112)), Institut de Planétologie et d'Astrophysique de Grenoble (IPAG), Centre National d'Études Spatiales [Toulouse] (CNES)-Université Grenoble Alpes (UGA)-Centre National de la Recherche Scientifique (CNRS)-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire des Sciences de l'Univers de Grenoble (OSUG), Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National de la Recherche Scientifique (CNRS)-Université Grenoble Alpes (UGA)-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Joseph Fourier - Grenoble 1 (UJF)-Centre National de la Recherche Scientifique (CNRS), Centre National de la Recherche Scientifique (CNRS)-Institut national des sciences de l'Univers (INSU - CNRS)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Université Paris Diderot - Paris 7 (UPD7), Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS)-Observatoire des Sciences de l'Univers de Grenoble (OSUG), and Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Centre National de la Recherche Scientifique (CNRS)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Centre National de la Recherche Scientifique (CNRS)
- Subjects
Angular momentum ,FOS: Physical sciences ,Astrophysics ,Magnetohydrodynamic turbulence ,01 natural sciences ,magnetohydrodynamics (MHD) ,methods: numerical ,accretion ,Magnetorotational instability ,0103 physical sciences ,Magnetic pressure ,Magnetic Prandtl number ,010306 general physics ,010303 astronomy & astrophysics ,Solar and Stellar Astrophysics (astro-ph.SR) ,Earth and Planetary Astrophysics (astro-ph.EP) ,Physics ,Turbulence ,accretion disks ,turbulence ,protoplanetary disks ,Fluid Dynamics (physics.flu-dyn) ,Astronomy and Astrophysics ,Physics - Fluid Dynamics ,Computational physics ,Magnetic field ,Amplitude ,Astrophysics - Solar and Stellar Astrophysics ,Space and Planetary Science ,[PHYS.ASTR]Physics [physics]/Astrophysics [astro-ph] ,Astrophysics - Earth and Planetary Astrophysics - Abstract
Angular momentum transport in accretion discs is often believed to be due to magnetohydrodynamic turbulence mediated by the magnetorotational instability. Despite an abundant literature on the MRI, the parameters governing the saturation amplitude of the turbulence are poorly understood and the existence of an asymptotic behavior in the Ohmic diffusion regime is not clearly established. We investigate the properties of the turbulent state in the small magnetic Prandtl number limit. Since this is extremely computationally expensive, we also study the relevance and range of applicability of the most common subgrid scale models for this problem. Unstratified shearing boxes simulations are performed both in the compressible and incompressible limits, with a resolution up to 800 cells per disc scale height. The latter constitutes the largest resolution ever attained for a simulation of MRI turbulence. In the presence of a mean magnetic field threading the domain, angular momentum transport converges to a finite value in the small Pm limit. When the mean vertical field amplitude is such that {\beta}, the ratio between the thermal and magnetic pressure, equals 1000, we find {\alpha}~0.032 when Pm approaches zero. In the case of a mean toroidal field for which {\beta}=100, we find {\alpha}~0.018 in the same limit. Both implicit LES and Chollet-Lesieur closure model reproduces these results for the {\alpha} parameter and the power spectra. A reduction in computational cost of a factor at least 16 (and up to 256) is achieved when using such methods. MRI turbulence operates efficiently in the small Pm limit provided there is a mean magnetic field. Implicit LES offers a practical and efficient mean of investigation of this regime but should be used with care, particularly in the case of a vertical field. Chollet-Lesieur closure model is perfectly suited for simulations done with a spectral code., Comment: Accepted for publication in A&A
- Published
- 2015