Related papers: A Tetrad-First Approach to Robust Numerical Algori…
We analyze inexact Riemannian gradient descent (RGD) where Riemannian gradients and retractions are inexactly (and cheaply) computed. Our focus is on understanding when inexact RGD converges and what is the complexity in the general…
We compare a particular selection of approximate solutions of the Riemann problem in the context of ideal relativistic magnetohydrodynamics. In particular, we focus on Riemann solvers not requiring a full eigenvector structure. Such solvers…
A main disadvantage of many high-order methods for hyperbolic conservation laws lies in the famous Gibbs-Wilbraham phenomenon, once discontinuities appear in the solution. Due to the Gibbs-Wilbraham phenomenon, the numerical approximation…
Accretion of magnetized gas on compact astrophysical objects such as black holes has been successfully modeled using general relativistic magnetohydrodynamic (GRMHD) simulations. These simulations have largely been performed in the Kerr…
The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms…
In this paper, we are concerned with the recovery of the geometric shapes of inhomogeneous inclusions from the associated far field data in electrostatics and acoustic scattering. We present a local resolution analysis and show that the…
In this paper, we propose a novel approach to hyperspectral image super-resolution by modeling the global spatial-and-spectral correlation and local smoothness properties over hyperspectral images. Specifically, we utilize the tensor…
The relativistic magnetohydrodynamics (RMHD) set of equations has recently seen increased use in astrophysical computations. Even so, RMHD codes remain fragile. The reconstruction can sometimes yield superluminal velocities in certain parts…
We describe the implementation of sophisticated numerical techniques for general-relativistic magnetohydrodynamics simulations in the Athena++ code framework. Improvements over many existing codes include the use of advanced Riemann solvers…
Computing geodesics for Riemannian manifolds is a difficult task that often relies on numerical approximations. However, these approximations tend to be either numerically unstable, have slow convergence, or scale poorly with manifold…
We present a new numerical method of special relativistic resistive magnetohydrodynamics with scalar resistivity that can treat a range of phenomena, from nonrelativistic to relativistic (shock, contact discontinuity, and Alfv\'en wave).…
The Riemann problem, and the associated generalized Riemann problem, are increasingly seen as the important building blocks for modern higher order Godunov-type schemes. In the past, building a generalized Riemann problem solver was seen as…
We investigate the revised Deser-Woodard model of nonlocal gravity involving four auxiliary scalar fields, introduced to explain the standard cosmological background expansion history without fine-tuning issues. In particular, we simplify…
Many astrophysical systems of interest to numerical relativity-such as rapidly rotating stars, black hole accretion disks, and core-collapse supernovae-exhibit near-symmetries. These systems generally consist of a strongly gravitating…
We construct a new Godunov type relativistic hydrodynamics code in Milne coordinates, using a Riemann solver based on the two-shock approximation which is stable under the existence of large shock waves. We check the correctness of the…
A Continuous Galerkin method-based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial…
General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole shadows from the…
We present two-dimensional general relativistic radiative magnetohydrodynamical simulations of accretion disks around non-rotating stellar-mass black hole. We study the evolution of an equilibrium accreting torus in different grid…
We present a simple method for applying excision boundary conditions for the relativistic Euler equations. This method depends on the use of Reconstruction-Evolution methods, a standard class of HRSC methods. We test three different…
We revisit the definition of transverse frames and tetrad choices with regards to its application to numerically generated spacetimes, in particular those from the merger of binary black holes. We introduce the concept of local and…