Related papers: Fast Higher-Order Interpolation and Restriction in…
On the basis of a recently proposed strategy of finite element integration in time domain for partial differential equations with a singular source term, we present a fourth order algorithm for non-rotating black hole perturbations in the…
Motivated by the substantial instability of the fundamental and high-overtone quasinormal modes, recent developments regarding the notion of black hole pseudospectrum call for numerical results with unprecedented precision. This work…
Electromagnetic waves interacting with three--dimensional periodic structures occur in many applications of great scientific and engineering interest. These three dimensional interactions are extremely complicated and subtle, so it is…
General-relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a graphics processing unit (GPU) accelerated GRMHD code \hammer{} with multi-faceted optimizations…
With the forthcoming VLBI images of Sgr A* and M87, simulations of accretion flows onto black holes acquire a special importance to aid with the interpretation of the observations and to test the predictions of different accretion…
Mesh generation is essential for accurate and efficient computational fluid dynamics simulations. To resolve critical features in the flow, adaptive mesh refinement (AMR) is routinely employed in certain regions of the computational domain,…
A single-step high-order implicit time integration scheme with controllable numerical dissipation at high frequencies is presented for the transient analysis of structural dynamic problems. The amount of numerical dissipation is controlled…
The application of suitable numerical boundary conditions for hyperbolic conservation laws on domains with complex geometry has become a problem with certain difficulty that has been tackled in different ways according to the nature of the…
As an entry for the 2001 Gordon Bell Award in the "special" category, we describe our 3-d, hybrid, adaptive mesh refinement (AMR) code, Enzo, designed for high-resolution, multiphysics, cosmological structure formation simulations. Our…
Identifying weak gravitational wave signals in noise and estimating the source properties require high-precision waveform templates. Numerical relativity (NR) simulations can provide the most accurate waveforms. However, it is challenging…
Gravitational instabilities naturally give rise to multi-scale structure, which is difficult for traditional Eulerian hydrodynamic methods to accurately evolve. This can be circumvented by adaptively adding resolution (in the form of…
Tropical cyclones (TCs) are powerful, natural phenomena that can severely impact populations and infrastructure. Enhancing our understanding of the mechanisms driving their intensification is crucial for mitigating these impacts. To this…
(abridged) In the first two papers of this publication series, we present a comprehensive library of synthetic EHT observations and used this library to train and validate Bayesian neural networks for the parameter inference of accreting…
New implicit and implicit-explicit time-stepping methods for the wave equation in second-order form are described with application to two and three-dimensional problems discretized on overset grids. The implicit schemes are single step,…
According to the theory of general relativity, the relative acceleration of masses generates gravitational radiation. Although gravitational radiation has not yet been detected, it is believed that extremely violent cosmic events, such as…
Observations of supermassive black holes by the Event Horizon Telescope reveal significant inhomogeneities, most likely related to density and magnetic field perturbations. To model these features, we conduct high-resolution 2D…
High-precision regression of physical parameters from black hole images generated by General Relativistic Ray Tracing (GRRT) is essential for investigating spacetime curvature and advancing black hole astrophysics. However, due to…
The efficient condition assessment of engineered systems requires the coupling of high fidelity models with data extracted from the state of the system `as-is'. In enabling this task, this paper implements a parametric Model Order Reduction…
Binary black hole (BBH) mergers detected via gravitational waves are addressing key open questions in astrophysics, cosmology, and fundamental physics. Our scientific conclusions rely on extracting accurate source parameters, for which we…
We propose the use of automatic differentiation through the programming framework jax for accelerating a variety of analysis tasks throughout gravitational wave (GW) science. Firstly, we demonstrate that complete waveforms which cover the…