Related papers: Implementation of the Shearing Box Approximation i…
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to…
We describe a new Godunov algorithm for relativistic magnetohydrodynamics (RMHD) that combines a simple, unsplit second order accurate integrator with the constrained transport (CT) method for enforcing the solenoidal constraint on the…
This study aims to develop second-order relativistic viscous magnetohydrodynamics (MHD) derived from kinetic theory within an extended relaxation time approximation (momentum/energy dependent) for the collision kernel. The investigation…
We present a structure-preserving discretization of the hybrid magnetohydrodynamics (MHD)-driftkinetic system for simulations of low-frequency wave-particle interactions. The model equations are derived from a variational principle,…
We develop a new relativistic radiation hydrodynamics code based on the Monte-Carlo algorithm. In this code, we implement a new scheme to achieve the second-order accuracy in time in the limit of a large packet number for solving the…
We present the implementation of a three-dimensional, second order accurate Godunov-type algorithm for magneto-hydrodynamic (MHD), in the adaptive-mesh-refinement (AMR) cosmological code {\tt CHARM}. The algorithm is based on the full…
We describe a new algorithm to solve the time dependent, frequency integrated radiation transport (RT) equation implicitly, which is coupled to an explicit solver for equations of magnetohydrodynamics (MHD) using {\sf Athena++}. The…
We show how the local approximation of astrophysical discs, which is the basis for the well known model of the shearing box, can be used to study many aspects of the dynamics of warped discs. In the local model, inclination of the orbit of…
We employ a variety of numerical simulations in the local shearing box system to investigate in greater depth the local hydrodynamic stability of Keplerian differential rotation. In particular we explore the relationship of Keplerian shear…
We present a fourth-order accurate finite volume method for the solution of ideal magnetohydrodynamics (MHD). The numerical method combines high-order quadrature rules in the solution of semi-discrete formulations of hyperbolic conservation…
In numerical models of thin astrophysical disks that use an Eulerian scheme, gas orbits supersonically through a fixed grid. As a result the time step is sharply limited by the Courant condition. Also, because the mean flow speed with…
We explore the thermodynamic and global structural properties of a local patch of an accretion disk whose parameters were chosen so that radiation pressure and gas pressure would be comparable in magnitude. Heating, radiative transport, and…
This paper presents a series of stratified shearing-box simulations of collisionless accretion disks in the recently developed framework of kinetic magnetohydrodynamics (MHD), which can handle finite non-gyrotropy of a pressure tensor.…
We present here both analytical and numerical results of hydrodynamic stability investigations of rotationally supported circumstellar flows using the shearing box formalism. Asymptotic scaling arguments justifying the shearing box…
We numerically evolve turbulence driven by the magnetorotational instability (MRI) in a 3D, unstratified shearing box and study its structure using two-point correlation functions. We confirm Fromang and Papaloizou's result that shearing…
We present a divergence-free semi-implicit finite volume scheme for the simulation of the ideal magnetohydrodynamics (MHD) equations which is stable for large time steps controlled by the local transport speed at all Mach and Alfv\'en…
Magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability can provide diffusive transport of angular momentum in astrophysical disks, and a widely studied computational model for this process is the ideal, stratified,…
Magnetically driven mass accretion in protoplanetary disks plays a crucial role in understanding disk evolution and planet formation. However, the $\alpha$ prescription lacks a direct connection to physical processes, and no systematic…
Due to the prevalence of magnetic fields in astrophysical environments, magnetohydrodynamic (MHD) simulation has become a basic tool for studying astrophysical fluid dynamics. To further advance the precision of MHD simulations, we have…
Axisymmetric incompressible modes of the magneto-rotational instability (MRI) with a vertical wavenumber are exact solutions of the non-linear local equations of motion for a disk (shearing box). They are referred to as "channel solutions".…