Related papers: cuHARM : a new GPU accelerated GR-MHD code and its…
The profound comprehension of the evolution and phenomenology of an Active Galactic Nucleus requires an accurate exploration of the dynamics of the magnetized gaseous disk surrounding the massive black hole in the centre. Many numerical…
The Event Horizon Telescope (EHT) has recently released high-resolution images of accretion flows onto two supermassive black holes. Our physical understanding of these images depends on accuracy and precision of numerical models of plasma…
Perturber objects interacting with supermassive black hole accretion disks are often invoked to explain observed quasi-periodic behavior in active galactic nuclei (AGN). We present global, 3D general relativistic magnetohydrodynamic (GRMHD)…
Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel…
We study the general-relativistic dynamics of matter being accreted onto and ejected by a magnetised and nonrotating neutron star. The dynamics is followed in the framework of fully general relativistic magnetohydrodynamics (GRMHD) within…
Gas falling into a black hole (BH) from large distances is unaware of BH spin direction, and misalignment between the accretion disc and BH spin is expected to be common. However, the physics of tilted discs (e.g., angular momentum…
It is widely accepted that quasars and other active galactic nuclei (AGN) are powered by accretion of matter onto a central supermassive black hole. While numerical simulations have demonstrated the importance of magnetic fields in…
In this paper, we develop a new parallel auxiliary grid algebraic multigrid (AMG) method to leverage the power of graphic processing units (GPUs). In the construction of the hierarchical coarse grid, we use a simple and fixed coarsening…
Many problems at the forefront of theoretical astrophysics require the treatment of magnetized fluids in dynamical, strongly curved spacetimes. Such problems include the origin of gamma-ray bursts, magnetic braking of differential rotation…
This work summarizes a program intended to unify three burgeoning branches of the high-energy astrophysics of relativistic jets: general relativistic magnetohydrodynamic (GRMHD) simulations of ever-increasing dynamical range, the…
There is an ongoing effort to develop tools that apply distributed computational resources to tackle large problems or reduce the time to solve them. In this context, the Alternating Direction Method of Multipliers (ADMM) arises as a method…
We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements…
We present numerical magnetohydrodynamic (MHD) simulations of a magnetized accretion disk launching trans-Alfvenic jets. These simulations, performed in a 2.5 dimensional time-dependent polytropic resistive MHD framework, model a resistive…
Recent numerical simulations of MHD jets from accretion disks are briefly reviewed with emphasis on the scaling law for jet speed and the role of magnetic reconnection in relation to time variability in accretion disks, jets, and flares. On…
Two-dimensional magnetohydrodynamic simulations are performed using the ZEUS-2D code to investigate the dynamics of a collapsar that generates a GRB jet, taking account of realistic equation of state, neutrino cooling and heating processes,…
We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an…
We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle…
Recent advances in general relativistic magnetohydrodynamic simulations have expanded and improved our understanding of the dynamics of black-hole accretion disks. However, current simulations do not capture the thermodynamics of electrons…
A finite-difference Micromagnetic simulation code written in MATLAB is presented with Graphics Processing Unit (GPU) acceleration. The high performance of Graphics Processing Unit (GPU) is demonstrated compared to a typical Central…
In this work we present a proof of concept of CUDA-capable, resistive, multi-fluid models of relativistic magnetohydrodynamics (RMHD). Resistive and multi-fluid codes for simulating models of RMHD suffer from stiff source terms, so it is…