Related papers: GRMHD in axisymmetric dynamical spacetimes: the X-…
Astrophysical plasmas in relativistic spacetimes, such as black hole accretion flows, are often weakly collisional and require kinetic modeling to capture non-local transport and particle acceleration. However, the extreme scale separation…
Various radio observations have showed that the hot atmospheres of galaxy clusters are magnetized. However, our understanding of the origin of these magnetic fields, their implications on structure formation and their interplay with the…
We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and…
We analyze the cosmic-ray magnetohydrodynamic (CR MHD) equations to improve the numerical simulations. We propose to solve them in the fully conservation form, which is equivalent to the conventional CR MHD equations. In the fully…
We present two general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of magnetically arrested disks (MADs) around non-spinning ($a_*=0$) and spinning ($a_*=0.9$) supermassive black holes (BHs). In each simulation, the…
The detection of gravitational waves has opened a new era for astronomy, allowing for the combined use of gravitational wave and electromagnetic emissions to directly probe the physics of compact objects, still poorly understood. So far,…
We present a new two-dimensional numerical code called Nada designed to solve the full Einstein equations coupled to the general relativistic hydrodynamics equations. The code is mainly intended for studies of self-gravitating accretion…
A new formulation of time-dependent Relaxed Magnetohydrodynamics (RxMHD) is derived variationally from Hamilton's Action Principle using microscopic conservation of mass, and macroscopic conservation of total magnetic helicity, cross…
Magnetic field strengths inferred for relativistic outflows including gamma-ray bursts (GRB) and active galactic nuclei (AGN) are larger than naively expected by orders of magnitude. We present three-dimensional relativistic…
This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not admit a…
We describe a method of implementing the axisymmetric evolution of general-relativistic hydrodynamics and magnetohydrodynamics through modification of a multipatch grid scheme. In order to ease the computational requirements required to…
We present results of models of the physical space and parameters of the accretion disk of Sagittarius A*, as well as simulations of its emergent spectrum. This begins with HARM, a 2D general relativistic magneto-hydrodynamic (GRMHD) model,…
This paper resolves the global regularity problem for the three-dimensional compressible magnetohydrodynamics (MHD) equations in the three-dimensional whole space, in the presence of a background magnetic field. Motivated by geophysical…
Recent observations by the Event Horizon Telescope (EHT) of supermassive black holes M87* and Sgr A* offer valuable insights into their spacetime properties and astrophysical conditions. Utilizing a library of model images (~2 million for…
The study of the electrodynamics of static, axisymmetric and force-free Kerr magnetospheres relies vastly on solutions of the so called relativistic Grad-Shafranov equation (GSE). Different numerical approaches to the solution of the GSE…
This paper presents both a numerical method for general relativity and an application of that method. The method involves the use of harmonic coordinates in a 3+1 code to evolve the Einstein equations with scalar field matter. In such…
We model the structure and evolution of black hole accretion disks, and their neighboring regions, using numerical simulations. The numerics is governed by the equations of general relativistic magneto-hydrodynamics (GRMHD). In particular,…
We describe a conservative, shock-capturing scheme for evolving the equations of general relativistic magnetohydrodynamics. The fluxes are calculated using the Harten, Lax, and van Leer scheme. A variant of constrained transport, proposed…
We describe an implicit general relativistic hydrodynamics code. The evolution equations are formulated in comoving coordinates. A conservative finite differencing of the Einstein equations is outlined, and artificial viscosity and…
We find exact and explicit solutions of the axisymmetric MHD equations of a self-gravitating polytropic gas. These solutions are able to describe a flat (uniform density) subsonic internal core, contracting homologously, of a collapsing…