Related papers: Performance-Portable Binary Neutron Star Mergers w…
We present the extension of the differentiable hydrodynamics code, diffhydro, enabling scalable PDE-constrained inference and integrated hybrid physics-ML models for a wide range of astrophysical applications. New physics additions include…
We present \texttt{SACRA-2D}, a new MPI and OpenMP parallelized, fully relativistic hydrodynamics (GRHD) code in dynamical spacetime under axial symmetry with the cartoon method using the finite-volume shock-capturing schemes for…
Efficient hybrid DFT simulations of solid state materials would be extremely beneficial for computational chemistry and materials science, but is presently bottlenecked by difficulties in computing Hartree-Fock (HF) exchange with plane wave…
We perform high-resolution magnetohydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese supercomputer K. The neutron stars and merger remnants are covered by a grid spacing of 70\,m, which yields…
We investigate close encounters by stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGN) as a potential formation channel of binary black holes (BBHs). We perform a series of 2D isothermal viscous…
The early post-merger phase of a binary neutron-star coalescence is shaped by characteristic rotational velocities as well as violent density oscillations and offers the possibility to constrain the properties of neutron star matter by…
We implement general relativistic hydrodynamics in the moving-mesh code AREPO. We also couple a solver for the Einstein field equations employing the conformal flatness approximation. The implementation is validated by evolving isolated…
We present the numerical methods and GPU-accelerated implementation underlying a Total Lagrangian finite element framework for finite-deformation flexible multibody dynamics, introduced in the companion paper [1]. The framework supports…
We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a…
The powerful high-energy phenomena typically encountered in astrophysics invariably involve physical engines, like neutron stars and black hole accretion disks, characterized by a combination of highly magnetized plasmas, strong…
The second paper of this series presents two robust entropy stable shock-capturing methods for discontinuous Galerkin spectral element (DGSEM) discretizations of the compressible magneto-hydrodynamics (MHD) equations. Specifically, we use…
We investigate the influence of magnetic fields on the evolution of binary neutron-star (BNS) merger remnants via three-dimensional (3D) dynamical-spacetime general-relativistic (GR) magnetohydrodynamic (MHD) simulations. We evolve a…
Axisymmetric magnetohydrodynamics (MHD) can be invoked for describing astrophysical magnetized flows and formulated to model stellar magnetospheres including main sequence stars (e.g. the Sun), compact stellar objects [e.g. magnetic white…
We investigate the thermodynamics of close encounters between stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGN), during which binary black holes (BBHs) may form. We consider a suite of 2D viscous…
We present a co-scaling grid formalism and its implementation in the magnetohydrodynamics code Athena++. The formalism relies on flow symmetries in astrophysical problems involving expansion, contraction, and center-of-mass motion. The grid…
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 a GPU-accelerated version of the real-space SPARC electronic structure code for performing hybrid functional calculations in generalized Kohn-Sham density functional theory. In particular, we develop a batch variant of the…
Simulations of relativistic hydrodynamics often need both high accuracy and robust shock-handling properties. The discontinuous Galerkin method combines these features --- a high order of convergence in regions where the solution is smooth…
We present new initial data codes for constructing stationary, axisymmetric equilibrium models of differentially rotating neutron stars in full general relativity within the Frankfurt University/KADATH (FUKA) suite of initial data codes.…
We present the first-ever moving-mesh general relativistic hydrodynamics solver for static spacetimes as implemented in the code, MANGA. Our implementation builds on the architectures of MANGA and the numerical relativity Python package…