Magnetorotational instability with smoothed particle hydrodynamics

Abstract

We present a thorough numerical study on the MRI using the smoothed particle magnetohydrodynamics method (SPMHD) with the geometric density average force expression (GDSPH). We perform shearing box simulations with different initial setups and a wide range of resolution and dissipation parameters. We show, for the first time, that MRI with sustained turbulence can be simulated successfully with SPH, with results consistent with prior work with grid-based codes. In particular, for the stratified boxes, our simulations reproduce the characteristic butterfly diagram of the MRI dynamo with saturated turbulence for at least 100 orbits. On the contrary, traditional SPH simulations suffer from runaway growth and develop unphysically large azimuthal fields, similar to the results from a recent study with mesh-less methods. We investigated the dependency of MRI turbulence on the numerical Prandtl number in SPH, focusing on the unstratified, zero net-flux case. We found that turbulence can only be sustained with a Prandtl number larger than $\sim$2.5, similar to the critical values of physical Prandtl number found in grid-code simulations. However, unlike grid-based codes, the numerical Prandtl number in SPH increases with resolution, and for a fixed Prandtl number, the resulting magnetic energy and stresses are independent of resolution. Mean-field analyses were performed on all simulations, and the resulting transport coefficients indicate no $\alpha$-effect in the unstratified cases, but an active $\alpha \Omega$ dynamo and a diamagnetic pumping effect in the stratified medium, which are generally in agreement with previous studies. There is no clear indication of a shear-current dynamo in our simulation, which is likely to be responsible for a weaker mean-field growth in the tall, unstratified, zero net-flux simulation.

References in corpus (21)

• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
• MHD simulations of the magnetorotational instability in a shearing box with zero net flux. II. The effect of transport coefficients
• Impact of dimensionless numbers on the efficiency of MRI-induced turbulent transport
• MHD simulations of the magnetorotational instability in a shearing box with zero net flux. I. The issue of convergence
• Magnetic field amplification in turbulent astrophysical plasmas
• Constrained Hyperbolic Divergence Cleaning for Smoothed Particle Magnetohydrodynamics
• The Turbulent Dynamo in Highly Compressible Supersonic Plasmas
• A self-sustaining nonlinear dynamo process in Keplerian shear flows
• On Self-Sustained Dynamo Cycles in Accretion Discs
• Nonlinear current helicity fluxes in turbulent dynamos and alpha quenching
• Aspect ratio dependence in magnetorotational instability shearing box simulations
• Simulations of Magnetorotational Turbulence with a Higher-Order Godunov Scheme
• Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds
• Magnetorotational instability driven dynamos at low magnetic Prandtl numbers
• Angular momentum transport and large eddy simulations in magnetorotational turbulence: the small Pm limit
• Statistical simulation of the magnetorotational dynamo
• The alpha effect with imposed and dynamo-generated magnetic fields
• Local simulations of MRI turbulence with meshless methods
• Exact shering box solutions of MHD flows with resistivity, viscosity and cooling
• The dual role of shear in large-scale dynamos

Cited by in corpus (9)

• Simulating the magnetorotational instability on a moving mesh with the shearing box approximation
• Simulating cold shear flows on a moving mesh
• On large-scale dynamos with stable stratification and the application to stellar radiative zones
• Turbulent Transport of Dust Particles in Protostellar Disks: The Effect of Upstream Diffusion
• Numerical dependencies of the galactic dynamo in isolated galaxies with SPH
• Three-dimensional simulations of the magnetorotational instability in eccentric disks
• Shedding light on the MRI driven dynamo in a stratified shearing box
• Unified treatment of mean-field dynamo and angular-momentum transport in magnetorotational instability-driven turbulence
• MHD in a cylindrical shearing box II: Intermittent Bursts and Substructures in MRI Turbulence