Related papers: The Lagrangian Numerical Relativity code SPHINCS_B…
We present a code for solving the coupled Einstein-hydrodynamics equations to evolve relativistic, self-gravitating fluids. The Einstein field equations are solved in generalized harmonic coordinates on one grid using pseudospectral…
Relativistic numerical cosmology is most often based either on the exact solutions of the Einstein equations, or perturbation theory, or weak-field limit, or the BSSN formalism. The Silent Universe provides an alternative approach to…
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…
We present a new three-dimensional general relativistic hydrodynamics code, the Whisky code. This code incorporates the expertise developed over the past years in the numerical solution of Einstein equations and of the hydrodynamics…
We examine the Lagrangian divergence of the displacement field, arguably a more natural object than the density in a Lagrangian description of cosmological large-scale structure. This quantity, which we denote \psi, quantifies the…
We have written and tested a new general relativistic magnetohydrodynamics (GRMHD) code, capable of evolving MHD fluids in dynamical spacetimes with adaptive-mesh refinement (AMR). Our code solves the Einstein-Maxwell-MHD system of coupled…
We extend a new finite element code, Einstein PHG (iPHG), to solve the evolution part of Einstein equations in first-order GH formalism. This paper is the third one of a systematic investigation of applying adaptive finite element method to…
We provide a description of the SNSPH code--a parallel 3-dimensional radiation hydrodynamics code implementing treecode gravity, smooth particle hydrodynamics, and flux-limited diffusion transport schemes. We provide descriptions of the…
Recent observations of compact binary systems have provided evidence for black holes with masses below standard expectations. When paired with neutron stars (NSs), such low-mass black holes (BHs) have the potential to be detected by future…
In this article, we propose a novel Stabilized Physics Informed Neural Networks method (SPINNs) for solving wave equations. In general, this method not only demonstrates theoretical convergence but also exhibits higher efficiency compared…
We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement…
A new numerical scheme to solve the Einstein field equations based upon the generalized harmonic decomposition of the Ricci tensor is introduced. The source functions driving the wave equations that define generalized harmonic coordinates…
Smoothed Particle Hydrodynamics is a multidimensional Lagrangian method of numerical hydrodynamics that has been used to tackle a wide variety of problems in astrophysics. Here we develop the basic equations of the SPH scheme, and we…
This work extends our previous study of isospin symmetric quarkyonic matter to quarkyonic neutron matter which might be relevant to the inner cores of neutron stars. The vector-isovector $\rho$ mesons are introduced to the model mainly to…
Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations -- the impact of massive neutrinos should therefore be included in these…
In the extreme violence of merger and mass accretion, compact objects like black holes and neutron stars are thought to launch some of the most luminous outbursts of electromagnetic and gravitational wave energy in the Universe. Modeling…
We present an open-source update to the spherically-symmetric, general-relativistic hydrodynamics, core-collapse supernova (CCSN) code GR1D (O'Connor & Ott, 2010, CQG, 27, 114103). The source code is available at http://www.GR1Dcode.org. We…
Luminous Red Novae (LRNe) are enigmatic transient events distinguished by a rapid rise in luminosity, a plateau in luminosity, and spectra which become redder with time. The best-observed system before, during, and after the outburst is…
One of the computational challenges of cluster formation simulations is resolving individual stars and simulating massive clusters with masses of more than $10^4 M_{\odot}$ without gravitational softening. Combining direct $N$-body code…
Within a cosmological context, we study the behaviour of collisionless particles in the weak field approximation to General Relativity, allowing for large gradients of the fields and relativistic velocities for the particles. We consider a…