Radiation Magnetohydrodynamics In Global Simulations Of Protoplanetary Disks

Abstract

Our aim is to study the thermal and dynamical evolution of protoplanetary disks in global simulations, including the physics of radiation transfer and magneto-hydrodynamic (MHD) turbulence caused by the magneto-rotational instability. We develop a radiative transfer method based on the flux-limited diffusion approximation that includes frequency dependent irradiation by the central star. This hybrid scheme is implemented in the PLUTO code. The focus of our implementation is on the performance of the radiative transfer method. Using an optimized Jacobi preconditioned BiCGSTAB solver, the radiative module is three times faster than the MHD step for the disk setup we consider. We obtain weak scaling efficiencies of 70% up to 1024 cores. We present the first global 3D radiation MHD simulations of a stratified protoplanetary disk. The disk model parameters are chosen to approximate those of the system AS 209 in the star-forming region Ophiuchus. Starting the simulation from a disk in radiative and hydrostatic equilibrium, the magnetorotational instability quickly causes MHD turbulence and heating in the disk. For the disk parameters we use, turbulent dissipation heats the disk midplane and raises the temperature by about 15% compared to passive disk models. A roughly flat vertical temperature profile establishes in the disk optically thick region close to the midplane. We reproduce the vertical temperature profile with a viscous disk models for which the stress tensor vertical profile is flat in the bulk of the disk and vanishes in the disk corona. The present paper demonstrates for the first time that global radiation MHD simulations of turbulent protoplanetary disks are feasible with current computational facilities. This opens up the windows to a wide range of studies of the dynamics of protoplanetary disks inner parts, for which there are significant observational constraints.

Accepted to A&A

References in corpus (8)

• PLUTO: a Numerical Code for Computational Astrophysics
• Protoplanetary Disk Structures in Ophiuchus
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• Global MHD simulations of stratified and turbulent protoplanetary discs. I. Model properties
• Radiation hydrodynamics integrated in the code PLUTO
• Thermodynamics of an Accretion Disk Annulus with Comparable Radiation and Gas Pressure
• Surface Structure in an Accretion Disk Annulus with Comparable Radiation and Gas Pressure
• Global simulations of magnetorotational turbulence I: convergence and the quasi-steady state

Cited by in corpus (60)

• The Disk Substructures at High Angular Resolution Project (DSHARP) VI: Dust trapping in thin-ringed protoplanetary disks
• Global simulations of protoplanetary disks with net magnetic flux: I. Non-ideal MHD case
• Bifurcation of planetary building blocks during Solar System formation
• On the water delivery to terrestrial embryos by ice pebble accretion
• Vertical shear instability in accretion disc models with radiation transport
• Global Models of Planet Formation and Evolution
• Radiation Hydrodynamical Turbulence In Protoplanetary Disks: Numerical Models and Observational Constraints
• 3D Radiation Non-ideal Magnetohydrodynamical Simulations Of The Inner Rim In Protoplanetary Disks
• The efficiency of dust trapping in ringed proto-planetary discs
• Gas and dust dynamics in starlight-heated protoplanetary disks
• Constraints from Dust Mass and Mass Accretion Rate Measurements on Angular Momentum Transport in Protoplanetary Disks
• Spiral arms in the proto-planetary disc HD100453 detected with ALMA: evidence for binary-disc interaction and a vertical temperature gradient
• Radiation hydrodynamical models of the inner rim in protoplanetary disks
• Dust settling and rings in the outer regions of protoplanetary discs subject to ambipolar diffusion
• Ring formation and dust dynamics in wind-driven protoplanetary discs: global simulations
• Burst occurrence in young massive stellar objects
• Grand challenges in protoplanetary disc modelling
• Collapse of turbulent massive cores with ambipolar diffusion and hybrid radiative transfer I. Accretion and multiplicity
• On the filtering and processing of dust by planetesimals 1. Derivation of collision probabilities for non-drifting planetesimals
• On shocks driven by high-mass planets in radiatively inefficient disks. II. Three-dimensional global disk simulations
• Eccentricity excitation and merging of planetary embryos heated by pebble accretion
• Magnetic Turbulence and Thermodynamics in the Inner Region of Protoplanetary Discs
• A general hybrid radiation transport scheme for star formation simulations on an adaptive grid
• Thermodynamics of the dead-zone inner edge in protoplanetary disks
• Collapse of turbulent massive cores with ambipolar diffusion and hybrid radiative transfer II. Outflows
• Impact of convection and resistivity on angular momentum transport in dwarf novae
• On the episodic excursions of massive protostars in the Hertzsprung-Russell diagram
• Global MHD Simulations of Accretion Disks in Cataclysmic Variables (CVs): II The Relative Importance of MRI and Spiral Shocks
• A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
• Protoplanetary dust porosity and FU Orionis Outbursts: Solving the mystery of Earth's missing volatiles
• Temperature Fluctuations driven by Magnetorotational Instability in Protoplanetary Disks
• Turbulent and wind-driven accretion in dwarf novae threaded by a large scale magnetic field
• Radiation Hydrodynamics using Characteristics on Adaptive Decomposed Domains for Massively Parallel Star Formation Simulations
• The origin of the eccentricity of the hot Jupiter in CI Tau
• VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method
• Three-temperature radiation hydrodynamics with PLUTO: Tests and applications to protoplanetary disks
• Oscillatory migration of accreting protoplanets driven by a 3D distortion of the gas flow
• No self-shadowing instability in 2D radiation-hydrodynamical models of irradiated protoplanetary disks
• Radiation hydrodynamics including irradiation and adaptive mesh refinement with AZEuS. I. Methods
• Electric heating and angular momentum transport in laminar models of protoplanetary disks
• Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks I. Angular momentum transport and turbulent heating
• exoALMA XV: Interpreting the height of CO emission layer
• Radiation transport methods in star formation simulations
• Variability of the inner dead zone edge in 2D radiation hydrodynamic simulations
• Magnetohydrodynamics of protoplanetary discs
• The weakly nonlinear magnetorotational instability in a global, cylindrical Taylor-Couette flow
• Inner dusty regions of protoplanetary discs -- II. Dust dynamics driven by radiation pressure and disc winds
• Global simulations of magnetorotational turbulence II: turbulent energetics
• The impact of accretion heating and thermal conduction on the dead zone of protoplanetary disks
• Migration of gap-opening planets in 3D stellar-irradiated accretion disks
• Non-LTE radiation hydrodynamics in PLUTO
• Meridional Circulation driven by Planetary Spiral Wakes in Radiative and Magnetized Protoplanetary Discs
• Multidimensional half-moment multigroup radiative transfer. Improving moment-based thermal models of circumstellar disks
• Multi dust species inner rim in magnetized protoplanetary disks
• MHD in a cylindrical shearing box II: Intermittent Bursts and Substructures in MRI Turbulence
• The role of detailed gas and dust opacities in shaping the evolution of the inner disc edge subject to episodic accretion
• Resurgence of CO in a warm bubble around accreting protoplanets and its observability
• Observable scattered light features from inclined and non-inclined planets embedded in protoplanetary discs
• Starlight-driven flared-staircase geometry in radiation hydrodynamic models of protoplanetary disks
• Radiation-magnetohydrodynamics with MPI-AMRVAC using flux-limited diffusion