Related papers: GRAthena++: puncture evolutions on vertex-centered…
One of the most promising sources of gravitational radiation is coalescence of binary neutron stars or black holes. In order to study gravitational radiation at the merging phase of coalescing binary neutron stars which is called the last…
We present an update of the General-relativistic multigrid numerical (Gmunu) code, a parallelized, multi-dimensional curvilinear, general relativistic magnetohydrodynamics code with an efficient non-linear cell-centred multigrid (CCMG)…
Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in…
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…
With the advent of GRAVITY+, the upgrade to the beam combiner GRAVITY at the Very Large Telescope Interferometer (VLTI), fainter and higher redshift active galactic nuclei (AGNs) are becoming observable, opening an unprecedented opportunity…
We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and…
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…
Feedback from supernovae (SNe) is an essential mechanism that self-regulates the growth of galaxies, and a better model of SN feedback is still needed in galaxy formation simulations. In the first part of this paper, using an Eulerian…
Numerical methods for radiative transfer play a key role in modern-day astrophysics and cosmology, including study of the inhomogeneous reionization process. In this context, ray tracing methods are well-regarded for accuracy but notorious…
We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code which solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including…
In this thesis, we use numerical relativity to investigate gravitational waves from binary black holes in extensions of GR. We first study spherically symmetric gravitational collapse in cubic Horndeski theories of gravity. By varying the…
Active galactic nuclei (AGN) are accreting supermassive black holes co-evolving with their host galaxies through a complex interplay of feeding and feedback. In this chapter, we first discuss AGN fuelling in galaxies, both in interacting…
The new era of multimessenger astrophysics requires the capability of studying different aspects of the evolution of compact objects. In particular, the merger of neutron star binaries is a strong source of gravitational waves and…
Many astrophysical systems of interest to numerical relativity-such as rapidly rotating stars, black hole accretion disks, and core-collapse supernovae-exhibit near-symmetries. These systems generally consist of a strongly gravitating…
We have developed a new computer code, RELDAFNA, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite…
To overcome the communication bottlenecks observed in state-of-the-art parallel time-dependent flow solvers at extreme scales, an asynchronous computing approach that relaxes communication and synchronization at a mathematical level was…
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)…
Modern gravitational-wave science demands increasingly accurate and computationally intensive numerical relativity (NR) simulations. The Python-based, open-source NRPy framework generates optimized C/C++ code for NR, including the complete…
Emission and absorption line features are important diagnostics for the physics underlying extragalactic astronomy. The interpretation of observed signatures involves comparing against forward modeled spectra from galaxy formation…
A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the…