Related papers: GRMHD in axisymmetric dynamical spacetimes: the X-…
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…
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to…
Motivated by the desire for highly accurate numerical computations of compact binary spacetimes in the era of gravitational wave astronomy, we reexamine hyperbolicity and well-posedness of the initial value problem for popular models of…
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…
In this work, we introduce GRChombo: a new numerical relativity code which incorporates full adaptive mesh refinement (AMR) using block structured Berger-Rigoutsos grid generation. The code supports non-trivial "many-boxes-in-many-boxes"…
We describe a magnetohydrodynamic (MHD) constrained energy functional for equilibrium calculations that combines the topological constraints of ideal MHD with elements of Taylor relaxation. Extremizing states allow for partially chaotic…
The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate…
This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes.…
We report a new implementation for axisymmetric simulation in full general relativity. In this implementation, the Einstein equations are solved using the Nakamura-Shibata formulation with the so-called cartoon method to impose an…
This paper develops novel and robust central discontinuous Galerkin (CDG) schemes of arbitrarily high-order accuracy for special relativistic magnetohydrodynamics (RMHD) with a general equation of state (EOS). These schemes are provably…
I have developed two numerical codes to investigate the dynamics of collapsars. One is two-dimensional MHD code that are performed using the Newtonian (ZEUS-2D) code where realistic equation of state, neutrino cooling and heating processes…
We present a series of fully three-dimensional, dynamical-spacetime general relativistic magnetohydrodynamics (GRMHD) simulations of core-collapse supernovae (CCSNe) for a progenitor of zero-age-main-sequence (ZAMS) mass $25\, M_\odot$. We…
We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M$_\odot$ star with spectral neutrino transport. Our focus is to…
Higher order finite volume schemes for magnetohydrodynamics (MHD) and relativistic magnetohydrodynamics (RMHD) are very valuable because they allow us to carry out astrophysical simulations with very high accuracy. However, astrophysical…
Axisymmetric numerical simulations of rotating stellar core collapse to a neutron star are performed in the framework of full general relativity. The so-called Cartoon method, in which the Einstein field equations are solved in the…
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.…
The lack of a positive-definite and conserved energy is a serious obstacle in the black hole stability problem. In this work, we will show that there exists a positive-definite and conserved Hamiltonian energy for axially symmetric linear…
The current status of numerical solutions for the equations of ideal general relativistic hydrodynamics is reviewed. With respect to an earlier version of the article the present update provides additional information on numerical schemes…
This paper resolves the global regularity problem for the three-dimensional incompressible magnetohydrodynamics (MHD) equations in the upper half-space with slip boundary conditions, in the presence of a background magnetic field. Motivated…
We here present a new version of the publicly available general relativistic magnetohydrodynamic (GRMHD) code $\texttt{Spritz}$, which now includes an approximate neutrino leakage scheme able to handle neutrino cooling and heating. The…