Related papers: grim: A Flexible, Conservative Scheme for Relativi…
We present the first astrophysical application of the ECHO code in its recent version supplemented by a generalized Ohm law, namely a kinematic study of dynamo effects in thick accretion disks. High-order \emph{implicit-explicit}…
The role of radiation in general relativistic magnetohydrodynamic (GRMHD) accretion simulations is discussed through axisymmetric simulations of the evolution of an initial torus seeded with a weak magnetic field. The paper compares and…
Hyper-accreting discs occur in compact-object mergers and collapsars, and may power gamma-ray bursts (GRBs). We calculate the microscopic viscosity and resistivity of plasma in these discs, and discuss the implications for their global…
The correlated and coupled dynamics of accretion and outflow around black holes (BHs) are essentially governed by the fundamental laws of conservation as outflow extracts matter, momentum and energy from the accretion region. Here we…
We describe an implementation of a particle physics module available for the PLUTO code, appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a non-thermal component represented by relativistic charged…
Ideal systems like MHD and Euler flow may develop singularities in vorticity (w = curl v). Viscosity and resistivity provide dissipative regularizations of the singularities. In this paper we propose a minimal, local, conservative,…
We have developed a one-dimensional code to solve ultra-relativistic hydrodynamic problems, using the Glimm method for an accurate treatment of shocks and contact discontinuities. The implementation of the Glimm method is based on an exact…
Compressible magnetohydrodynamic (MHD) turbulence is ubiquitous in astrophysical phenomena ranging from the intergalactic to the stellar scales. In studying them, numerical simulations are nearly inescapable, due to the large degree of…
This paper presents an arbitrary h.o. accurate ADER DG method on space-time adaptive meshes (AMR) for the solution of two important families of non-linear time dependent PDE for compr. dissipative flows: the compr. Navier-Stokes equations…
A stable added-mass partitioned (AMP) algorithm is developed for fluid-structure interaction (FSI) problems involving viscous incompressible flow and compressible elastic solids. Deforming composite grids are used to effectively handle the…
For the first time ever, we derive equations governing the time-evolution of fully relativistic slim accretion disks in the Kerr metric, and numerically construct their detailed non-stationary models. We discuss applications of these…
An algorithm for simulating self-gravitating cosmological astrophysical fluids is presented. The advantages include a large dynamic range, parallelizability, high resolution per grid element and fast execution speed. The code is based on a…
Wormholes are a hypothetical object that connects disparate points in spacetime. It is a theoretically well-motivated black hole alternative and offers a potential observationally testable arena for probing strong-field gravity with…
Models of highly sub-Eddington accretion onto black holes commonly use a single fluid model for the collisionless, near-horizon plasma. These models must specify an equation of state. It is common to use an ideal gas with $p = (\gamma - 1)…
We have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite…
The constraint of incompressibility is often used to simplify the magnetohydrodynamic (MHD) description of linearized plasma dynamics because it does not affect the ideal MHD marginal stability point. In this paper two methods for…
General relativistic magnetohydrodynamic (GRMHD) simulations are providing influential models for black hole spin measurements, gamma ray bursts, and supermassive black hole feedback. Many of these simulations use the same initial…
(Abridged) We study stability of gas accretion in Active Galactic Nuclei. Our grid based simulations cover a radial range from 0.1 to 200 pc. Here, as in previous studies by our group, we include gas radiative cooling as well as heating by…
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…
The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms…