Related papers: H-AMR: A New GPU-accelerated GRMHD Code for Exasca…
Adaptive Mesh Refinement (AMR) with subcycling in time enables different grid levels to advance using their own time steps, ensuring finer grids employ smaller steps for accuracy while coarser grids take larger steps to improve…
(Abridged) We present one of the first physically-motivated two-dimensional general relativistic magnetohydrodynamic (GRMHD) numerical simulations of a radiatively-cooled black-hole accretion disk. The fiducial simulation combines a…
We present axisymmetric numerical simulations of radiatively inefficient accretion flows onto black holes combining general relativity, magnetohydrodynamics, self-consistent electron thermodynamics, and frequency-dependent radiation…
This work provides a proof of concept for the computation of pure gluonic amplitudes in quantum chromodynamics (QCD) on graphics processing units (GPUs). The implementation relies on the Berends-Giele recursion algorithm and, for the first…
Structured Adaptive Mesh Refinement (SAMR) is a popular numerical technique to study processes with high spatial and temporal dynamic range. It reduces computational requirements by adapting the lattice on which the underlying differential…
Large language models (LLMs) have been widely applied but face challenges in efficient inference. While quantization methods reduce computational demands, ultra-low bit quantization with arbitrary precision is hindered by limited GPU Tensor…
GPU accelerators have had a notable impact on high-performance computing across many disciplines. They provide high performance with low cost/power, and therefore have become a primary compute resource on many of the largest supercomputers.…
We present a set of eleven two-temperature, radiative, general relativistic magnetohydrodynamic (2TGRRMHD) simulations of the black hole M87* in the magnetically arrested (MAD) state, surveying different values of the black hole spin $a_*$.…
Accretion onto supermassive black holes (BHs) can launch relativistic jets that inject energy and momentum into their surroundings. Understanding how such feedback shapes large-scale accretion is key to bridging observations from galactic…
Rapid growth in scientific data and a widening gap between computational speed and I/O bandwidth make it increasingly infeasible to store and share all data produced by scientific simulations. Instead, we need methods for reducing data…
We study the general-relativistic dynamics of matter being accreted onto and ejected by a magnetised and nonrotating neutron star. The dynamics is followed in the framework of fully general relativistic magnetohydrodynamics (GRMHD) within…
General-relativistic radiative transfer (GRRT) calculations coupled with the calculation of geodesics in the Kerr spacetime are an essential tool for determining the images, spectra and light curves from matter in the vicinity of black…
We present a comprehensive analysis of super-Eddington black hole accretion simulations that solve the GRMHD equations coupled with angle-discretized radiation transport. The simulations span a range of accretion rates, two black hole…
Smoothed Particle Hydrodynamics (SPH) is essential for modeling complex large-deformation problems across various applications, requiring significant computational power. A major portion of SPH computation time is dedicated to the Nearest…
The remarkable results by the Event Horizon Telescope collaboration concerning the emission from M87* and, more recently, its polarization properties, require an increasingly accurate modeling of the plasma flows around the accreting black…
In this work, we introduce GRChombo: a new numerical relativity code which incorporates full adaptive mesh refinement (AMR) using block structured Berger-Rigoutsos grid generation. The code supports non-trivial "many-boxes-in-many-boxes"…
Clawpack is a library for solving nonlinear hyperbolic partial differential equations using high-resolution finite volume methods based on Riemann solvers and limiters. It supports Adaptive Mesh Refinement (AMR), which is essential in…
Understanding the dynamics of low angular momentum accretion flow around black holes (BHs) is essential for probing extreme plasma behavior in strong gravity, where shock formation can naturally produce variability signatures. In this…
This work presents a GPU thread mapping approach that allows doing fast parallel stencil-like computations on discrete fractals using their compact representation. The intuition behind is to employ two GPU tensor-core accelerated thread…
Efficient simulation of plasmas in various contexts often involves the use of meshes that conform to the intrinsic geometry of the system under consideration. We present here a description of a new magnetohydrodynamic code, Gamera (Grid…