Related papers: Implementation of the Shearing Box Approximation i…
The effect of large magnetic Prandtl number $\text{Pm}$ (the ratio of viscosity to resistivity) on the turbulent transport and energetics of the magnetorotational instability (MRI) is poorly understood, despite the realization of this…
We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and Xin.…
We study a truncated accretion disk using a well-resolved, semi-global magnetohydrodynamic simulation that is evolved for many dynamical times (6096 inner disk orbits). The spectral properties of hard state black hole binary systems and…
We investigate the evolution of black holes on orbits with small inclinations ($i < 2^\circ$) to the gaseous discs of active galactic nuclei. We perform 3D adiabatic hydrodynamic simulations within a shearing frame, studying the damping of…
We present global magnetohydrodynamic (MHD) simulations of accretion disks with a strong toroidal magnetic field using an equation of state that fixes the gas thermal scale height. The disk forms from the inflow of a rotating magnetized gas…
We analyze the long-term evolution of hierarchical triple systems in Newtonian gravity to second order in the quadrupolar perturbation parameter, and to sixth order in $\epsilon = a/A$, the ratio of the semimajor axes of the inner and outer…
The hydrodynamic response of the inviscid small shearing box model of a midplane section of a rotationally supported astrophysical disk is examined. An energy functional ${\cal E}$ is formulated for the general nonlinear problem. It is…
Previous work has shown that interactions between a central binary system and a circumbinary disc (CBD) can lead to the binary orbit either shrinking or expanding, depending on the properties of the disc. In this work, we perform…
We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in…
The traveltime of compressional (P) and shear (S) waves have proven essential in many applications of earthquake and exploration seismology. An accurate and efficient traveltime computation for P and S waves is crucial for the success of…
We present a constrained transport (CT) algorithm for solving the 3D ideal magnetohydrodynamic (MHD) equations on a moving mesh, which maintains the divergence-free condition on the magnetic field to machine-precision. Our CT scheme uses an…
We examine the effect of accuracy of high-order spectral element methods, with or without adaptive mesh refinement (AMR), in the context of a classical configuration of magnetic reconnection in two space dimensions, the so-called…
We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented…
We study super-Eddington accretion flows onto black holes using a global three dimensional radiation magneto-hydrodynamical simulation. We solve the time dependent radiative transfer equation for the specific intensities to accurately…
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve…
We construct a structure-preserving finite element method and time-stepping scheme for compressible barotropic magnetohydrodynamics (MHD) both in the ideal and resistive cases, and in the presence of viscosity. The method is deduced from…
We investigate variational principles for the approximation of quantum dynamics that apply for approximation manifolds that do not have complex linear tangent spaces. The first one, dating back to McLachlan (1964) minimizes the residuum of…
A method based on separated integration to estimate anharmonic corrections to energy and vibration of molecules in a second-order diagrammatic vibrational many-body Green's function formalism has already been presented. A severe bottleneck…
In this paper we describe a numerical method designed for modelling different kinds of astrophysical flows in three dimensions. Our method is a standard explicit finite difference method employing the local shearing-box technique. To model…
We describe local shearing box simulations of turbulence driven by the magnetorotational instability (MRI) in a collisionless plasma. Collisionless effects may be important in radiatively inefficient accretion flows, such as near the black…