Grid-based magnetohydrodynamic (MHD) simulations have proven invaluable for the study of astrophysical accretion disks. However, the fact that angular momentum transport in disks is mediated by MHD turbulence (with structure down to very small scales) raises the concern that the properties of the modeled accretion disks are affected by the finite numerical resolution of the simulation. By implementing an orbital advection algorithm into the Athena code in cylindrical geometry, we have performed a set of global (but unstratified) Newtonian disk simulations extending up to resolutions previously unattained. We study the convergence of these models as a function of spatial resolution and initial magnetic field geometry. The usual viscosity parameter (α) or the ratio of thermal-tomagnetic pressure (β) are found to be poor diagnostics of convergence, whereas the average tilt angle of the magnetic field in the (r, φ)-plane is a very good diagnostic of convergence. We suggest that this is related to the saturation of the MHD turbulence via parasitic modes of the magnetorotational instability. Even in the case of zero-net magnetic flux, we conclude that our highest resolution simulations (with 32-zones and 64-zones per vertical scale height) have achieved convergence.Our global simulations reach resolutions comparable to those used in local, shearing box models of MHD disk turbulence. We find that the saturation predictors derived from local simulations correspond well to the instantaneous correlations between local flux and stress found in our global simulations. However, the conservation of magnetic flux implicit in local models is not realized in our

Orbiting disks may exhibit bends due to a misalignment between the angular momentum of the inner and outer regions of the disk. We begin a systematic simulational inquiry into the physics of warped disks with the simplest case: the relaxation of an unforced warp under pure fluid dynamics, i.e. with no internal stresses other than Reynolds stress. We focus on the nonlinear regime in which the bend rate is large compared to the disk aspect ratio. When warps are nonlinear, strong radial pressure gradients drive transonic radial motions along the disk's top and bottom surfaces that efficiently mix angular momentum. The resulting nonlinear decay rate of the warp increases with the warp rate and the warp width, but, at least in the parameter regime studied here, is independent of the sound speed. The characteristic magnitude of the associated angular momentum fluxes likewise increases with both the local warp rate and the radial range over which the warp extends; it also increases with increasing sound speed, but more slowly than linearly. The angular momentum fluxes respond to the warp rate after a delay that scales with the square-root of the time for sound waves to cross the radial extent of the warp. These behaviors are at variance with a number of the assumptions commonly used in analytic models to describe linear warp dynamics.

During geomagnetic storms the intensities of the outer radiation belt electron population can exhibit dramatic variability. Deep depletions in intensity during the main phase are followed by increases during the recovery phase, often to levels that significantly exceed their prestorm values. To study these processes, we simulate the evolution of the outer radiation belt during the 17 March 2013 geomagnetic storm using our newly developed radiation belt model (Conservative Hamiltonian Integrator for Magnetospheric Particles) based on test particle and coupled 3‐D ring current and global magnetohydrodynamic (MHD) simulations, and driven solely with solar wind and F10.7 flux data. Our approach differs from previous work in that we use MHD information to identify regions of strong, bursty, and azimuthally localized earthward convection in the magnetotail where test particles are then seeded. We validate our model using in situ Van Allen Probe electron intensities over a multiday period and show that our model is able to reproduce meaningful qualitative and quantitative agreement. Analysis of our model enables us to study the processes that govern the transition from the prestorm to poststorm outer belt. Our analysis demonstrates that during the early main phase of the storm the preexisting outer belt is largely wiped out via magnetopause losses, and subsequently, a new outer belt is created during a handful of discrete, mesoscale injections. Finally, we demonstrate the potential importance of magnetic gradient trapping in the transport and energization of outer belt electrons using a controlled numerical experiment.

Efficient simulation of plasmas in various contexts often involves the use of meshes that conform to the intrinsic geometry of the system under consideration. We present here a description of a new magnetohydrodynamic code, Gamera (Grid Agnostic MHD for Extended Research Applications), designed to combine geometric flexibility with high-order spatial reconstruction and constrained transport to maintain the divergence-free magnetic field. Gamera carries on the legacy of its predecessor, the LFM (Lyon-Fedder-Mobarry), a research code whose use in space physics has spanned three decades. At the time of its initial development the LFM code had a number of novel features: eighth-order centered spatial differencing, the Partial Donor Cell Method limiter for shock capturing, a non-orthogonal staggered mesh with constrained transport, and conservative averaging-reconstruction for axis singularities. A capability to handle multiple ion species was also added later. Gamera preserves the core numerical philosophy of LFM while also incorporating numerous algorithmic and computational improvements. The upgrades in the numerical schemes include accurate grid metric calculations using high-order Gaussian quadrature techniques, high-order upwind reconstruction, non-clipping options for interface values, and improved treatment of axis singularities. The improvements in the code implementation include the use of data structures and memory access patterns conducive to aligned vector operations and the implementation of hybrid parallelism, using MPI and OMP. Gamera is designed to be a portable and easy-to-use code that implements multi-dimensional MHD simulations in arbitrary non-orthogonal curvilinear geometries on modern supercomputer architectures.