Related papers: The Lagrangian Numerical Relativity code SPHINCS_B…
We discuss a new algorithm to generate multi-scale initial conditions with multiple levels of refinements for cosmological "zoom-in" simulations. The method uses an adaptive convolution of Gaussian white noise with a real space transfer…
The precise measurement of neutron star (NS) spins can provide important insight into the formation and evolution of compact binaries containing NS. While traditional methods of NS spin measurement rely on pulsar observations, gravitational…
A new lattice based scheme for numerical relativity will be presented. The scheme uses the same data as would be used in the Regge calculus (eg. a set of leg lengths on a simplicial lattice) but it differs significantly in the way that the…
In this work numerical methods for solving Einstein's equations are developed and applied to the study of inhomogeneous cosmological models. A two-dimensional computer code is described which implements two advanced numerical methods:…
We study the behavior of Lagrangian perturbative solutions. For a spherical void model, the higher order the Lagrangian perturbation we consider, the worse the approximation becomes in late-time evolution. In particular, if we stop to…
This is the first in a series of papers on the construction and validation of a three-dimensional code for general relativistic hydrodynamics, and its application to general relativistic astrophysics. This paper studies the consistency and…
We present the first quantitative comparison of two independent general-relativistic hydrodynamics codes, the Whisky code and the SACRA code. We compare the output of simulations starting from the same initial data and carried out with the…
In this thesis different numerical methods, as well as applications of the methods to a number of current problems in relativistic astrophysics, are presented. In the first part the theoretical foundation and numerical implementation of a…
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…
We present a numerical analysis of the cosmological evolution of scalar field dark matter (SFDM) in the Boltzmann code $\texttt{CLASS}$, based on a dynamical system analysis of previous works. We show a detailed study of the evolution of…
Building efficient, accurate and generalizable reduced order models of developed turbulence remains a major challenge. This manuscript approaches this problem by developing a hierarchy of parameterized reduced Lagrangian models for…
We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989)…
Direct numerical simulation of subsonic turbulence with smoothed particle hydrodynamics (SPH) has traditionally been hampered by zeroth-order (E0) errors, inaccurate gradient evaluations, and excessive numerical dissipation. We demonstrate…
Black hole-neutron star (BHNS) binary mergers are candidate engines for generating both short-hard gamma-ray bursts (SGRBs) and detectable gravitational waves. Using our most recent conformal thin-sandwich BHNS initial data and our fully…
Physics-Informed Neural Networks (PINNs) have gained widespread popularity for solving inverse and forward problems across a range of scientific and engineering domains. However, most existing PINN frameworks are limited to the Eulerian…
While the use of numerical general relativity for modeling astrophysical phenomena and compact objects is commonplace, the application to cosmological scenarios is only just beginning. Here, we examine the expansion of a spacetime using the…
We derive the equation of matter density perturbations on sub-horizon scales for a general Lagrangian density f(R, phi, X) that is a function of a Ricci scalar R, a scalar field phi and a kinetic term X=-(nabla phi)^2/2. This is useful to…
We experiment with several new modifications for the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of the Einstein field equations and demonstrate how these modifications affect the stability of numerical black hole evolution…
There has been interest in recent years to assess the ability of astrophysical hydrodynamics codes to correctly model the Kelvin-Helmholtz instability. Smoothed particle hydrodynamics (SPH), in particular, has received significant…
We describe the implementation and testing of a smoothed particle hydrodynamics (SPH) code that solves the equations of radiation hydrodynamics in the flux-limited diffusion (FLD) approximation. The SPH equations of radiation hydrodynamics…