Related papers: An Extension of the Athena++ Code Framework for Ra…
In paper I, we showed that time-dependent general relativistic magnetohydrodynamic (GRMHD) numerical models of accretion disks, although being highly turbulent, have surprisingly simple electromagnetic properties. In particular, the…
In the past few years, the Event Horizon Telescope (EHT) has provided the first-ever event horizon-scale images of the supermassive black holes (BHs) (M87*) and Sagittarius A$^*$ (Sgr A*). The next-generation EHT project is an extension of…
We present models of Galactic Center emission in the vicinity of Sagittarius A* that use parametrizations of the electron temperature or energy density. These models include those inspired by two-temperature general relativistic…
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 numerical study of relativistic magnetohydrodynamics (MHD) plays a crucial role in high-energy astrophysics, but unfortunately is computationally demanding, given the complex physics involved (high Lorentz factor flows, extreme…
Accurate modeling of active galactic nucleus (AGN) feedback, especially due to relativistic jets, is crucial for understanding the cool-core problem in galaxy clusters. We present a new subgrid method to model accretion onto and feedback…
This paper explores plasma wave modes using the extended magnetohydrodynamics (XMHD) model, incorporating Hall drift and electron inertia effects. We utilize the geometric optics ansatz to study perturbed quantities, with a focus on…
Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization…
The general-relativistic magnetohydrodynamical (GRMHD) formulation for black hole-powered jets naturally gives rise to a stagnation surface, wherefrom inflows and outflows along magnetic field lines that thread the black hole event horizon…
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…
Maintaining conservation laws in the fully discrete setting is critical for accurate long-time behavior of numerical simulations and requires accounting for discrete conservation properties in both space and time. This paper derives…
We construct finite element methods for the incompressible magnetohydrodynamics (MHD) system that precisely preserve magnetic and cross helicity, the energy law and the magnetic Gauss law at the discrete level. The variables are discretized…
We present a formulation and numerical algorithm to extend the scheme for grey radiation magneto-hydrodynamics (MHD) developed by Jiang (2021) to include the frequency dependence via the multi-group approach. The entire frequency space can…
We present the first dynamical model of plasma accretion onto traversable wormholes by performing General Relativistic magneto-hydrodynamical (GRMHD) simulations of the flow on both sides of the wormhole. We evolve the ideal MHD equations…
The classical, relativistic thin-disk theory of Novikov and Thorne (NT) predicts a maximum accretion efficiency of 40% for an optically thick, radiatively efficient accretion disk around a maximally spinning black hole (BH). However, when a…
We present the numerical relativity module within AthenaK, an open source performance-portable astrophysics code designed for exascale computing applications. This module employs the Z4c formulation to solve the Einstein equations. We…
We describe a novel splitting approach to numerical relativistic magnetohydrodynamics (RMHD) designed to expand its applicability to the domain of ultra-high magnetisation (high-$\sigma$). In this approach, the electromagnetic field is…
Cylindrical coordinates are often used in computational fluid dynamics, in particular, when one considers gas flow accreting onto a central object. Although the cylindrical coordinates have several advantages in describing rotation, they…
Hot, diffuse, relativistic plasmas such as sub-Eddington black hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects…
We provide a systematic description of the steps necessary -- and of the potential pitfalls to be encountered -- when implementing a two-moment scheme within an Implicit-Explicit (IMEX) scheme to include radiative-transfer contributions in…