Related papers: Implementation of the Shearing Box Approximation i…
This paper presents a new method to approximate the time-dependent convection-diffusion equations using conforming finite element methods, ensuring that the discrete solution respects the physical bounds imposed by the differential…
We present the new general-relativistic magnetohydrodynamics (GRMHD) capabilities of the Einstein Toolkit, an open-source community-driven numerical relativity and computational relativistic astrophysics code. The GRMHD extension of the…
We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume…
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…
Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfv\'{e}n waves, which interact and rapidly…
We present a systematic shearing-box investigation of MRI-driven turbulence in a weakly collisional plasma by including the effects of an anisotropic pressure stress, i.e. anisotropic (Braginskii) viscosity. We constrain the pressure…
Time-centered, hence second-order, methods for integrating the relativistic momentum of charged particles in an electromagnetic field are derived. A new method is found by averaging the momentum before use in the magnetic rotation term, and…
This paper investigates the various spherically symmetric wormhole solutions in the presence of tidal forces and applies numerous methods, such as test particle orbital dynamics, ray-tracing and microlensing. We make the theoretical…
We analyze the wave equation in mixed form, with periodic and/or Dirichlet homogeneous boundary conditions, and nonconstant coefficients that depend on the spatial variable. For the discretization, the weak form of the second equation is…
We obtain equations of motion for the boost-invariant expansion of a system of chiral particles. Our analysis is based on the Boltzmann equation for left- and right-handed massless particles in the relaxation time approximation. We assume…
We present high-resolution, three-dimensional hydrodynamic simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent medium on galactic scales, taking M87* as a typical case. The simulations use a new…
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…
Development of a new framework for derivation of order-by-order hydrodynamics from Boltzmann equation is necessary as the widely used Anderson-Witting formalism leads to violation of fundamental conservation laws when the relaxation-time…
We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a…
We report on the second phase of our study of slightly rotating accretion flows onto black holes. We consider magnetohydrodynamical (MHD) accretion flows with a spherically symmetric density distribution at the outer boundary, but with…
The Energy Conserving Descent (ECD) algorithm was recently proposed (De Luca & Silverstein, 2022) as a global non-convex optimization method. Unlike gradient descent, appropriately configured ECD dynamics escape strict local minima and…
We present a novel approach to solve the advection-diffusion equation under arbitrary transporting fields using a quantum-inspired 'Schrodingerisation' technique for Hamiltonian simulation. Although numerous methods exist for solving…
We propose several continuous data assimilation (downscaling) algorithms based on feedback control for the 2D magnetohydrodynamic (MHD) equations. We show that for sufficiently large choices of the control parameter and resolution and…
Second-order optimization methods are among the most widely used optimization approaches for convex optimization problems, and have recently been used to optimize non-convex optimization problems such as deep learning models. The widely…
The 2D time dependent solution of thin accretion disk in a close binary system have been presented on the equatorial plane around the Schwarzschild black hole. To do that, the special part of the General Relativistic Hydrodynamical(GRH)…