Related papers: IllinoisGRMHD: An Open-Source, User-Friendly GRMHD…
We assemble the equations of general relativistic magnetohydrodynamics (MHD) in 3+1 form. These consist of the complete coupled set of Maxwell equations for the electromagnetic field, Einstein's equations for the gravitational field, and…
We report on a new 3D numerical code designed to solve the Einstein equations for general vacuum spacetimes. This code is based on the standard 3+1 approach using cartesian coordinates. We discuss the numerical techniques used in developing…
The nascent field of gravitational-wave astronomy offers many opportunities for effective and inspirational astronomy outreach. Gravitational waves, the "ripples in space-time" predicted by Einstein's theory of General Relativity, are…
Heliospheric plasmas require multi-scale and multi-physics considerations. On one hand, MHD codes are widely used for global simulations of the solar-terrestrial environments, but do not provide the most elaborate physical description of…
We propose an explicit-implicit scheme for numerically solving Special Relativistic Radiation Hydrodynamic (RRHD) equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, 0th and 1st moment…
I present the publicly available code GLADIS (GLobal Accretion Disk Instability Simulation) developed in my reserach group over the years 2002-2017. It can be freely downloaded and modified by the users via the link from the Astrophysics…
We present preliminary results in our long-term project of studying the evolution of matter in a dynamical spacetime. To achieve this, we have developed a new code to evolve axisymmetric initial data sets corresponding to a black hole…
We present MHDuet, an open source evolution code for general relativistic magnetohydrodynamics with neutrino transport. The code solves the full set of Einstein equations coupled to a relativistic, magnetized fluid with an M1 neutrino…
We introduce a new GPU-accelerated general-relativistic magneto-hydrodynamic (GR-MHD) code based on HARM which we call cuHARM. The code is written in CUDA-C and uses OpenMP to parallelize multi-GPU setups. Our code allows us to run high…
Machine learning techniques offer a precious tool box for use within astronomy to solve problems involving so-called big data. They provide a means to make accurate predictions about a particular system without prior knowledge of the…
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 flexible spectral code for the study of general relativistic magnetohydrodynamics is presented. Aiming at investigating the physics of slowly rotating magnetized compact stars, this new code makes use of various physically motivated…
A wide variety of astrophysical phenomena involve the flow of turbulent magnetized gas with relativistic velocity or energy density. Examples include gamma-ray bursts, active galactic nuclei, pulsars, magnetars, micro-quasars, merging…
We investigate the ideal and incompressible magnetohydrodynamic (MHD) equations in three space dimensions for the development of potentially singular structures. The methodology consists in implementing the four-fold symmetries of the…
We present a new computational method for smoothly matching general relativistic ideal magnetohydrodynamics (MHD) to its force-free limit. The method is based on a flux-conservative formalism for MHD and its force-free limit, and a vector…
This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference Convex ENO…
Many astrophysical simulations involve extreme dynamic range of timescales around 'special points' in the domain (e.g. black holes, stars, planets, disks, galaxies, shocks, mixing interfaces), where processes on small scales couple strongly…
X-ray reflection spectroscopy is currently one of the leading techniques for studying the inner part of accretion disks around black holes, measuring black hole spins, and even testing fundamental physics in strong gravitational fields.…
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 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…