Related papers: Implementation of the Shearing Box Approximation i…
We present a new magnetohydrodynamic-particle-in-cell (MHD-PIC) code integrated into the Athena++ framework. It treats energetic particles as in conventional PIC codes while the rest of thermal plasmas are treated as background fluid…
We derive a relativistic magnetohydrodynamics (RMHD) theory for a dilute electron-ion gas governed by the Boltzmann-Vlasov equation, using the method of moments. This yields an extended MHD framework beyond standard astrophysical…
One of the unitary forms of the quantum mechanical time evolution operator is given by Cayley's approximation. A numerical implementation of the same involves the replacement of second derivatives in Hamiltonian with the three-point…
We present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit…
We test a new "hybrid" scheme for simulating dynamical fluid flows in which cylindrical components of the momentum are advected across a rotating Cartesian coordinate mesh. This hybrid scheme allows us to conserve angular momentum to…
We consider an incompressible magnetohydrodynamics (MHD) model in which the classical first-order time derivatives in the momentum and magnetic induction equations are replaced by variable-order Caputo time-fractional derivatives. This…
ADER (Arbitrary DERivative in space and time) methods for the time-evolution of hyperbolic conservation laws have recently generated a fair bit of interest. The ADER time update can be carried out in a single step, which is desirable in…
We perform a suite of two- and three-dimensional magnetohydrodynamic (MHD) simulations with the Athena code of the non-driven Kelvin-Helmholtz instability in the subsonic, weak magnetic field limit. Focusing the analysis on the non-linear…
Aims. We present a new radiation hydrodynamics code, called "ARK-RT" which uses a two-moment model with the M1 closure relation for radiative transfer. This code aims at being ready for high-performance computing, on exascale architectures.…
We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type…
We present a new approach to Eulerian computational fluid dynamics that is designed to work at high Mach numbers encountered in astrophysical hydrodynamic simulations. The Eulerian fluid conservation equations are solved in an adaptive…
We propose a new class of finite element approximations to ideal compressible magnetohydrodynamic equations in smooth regime. Following variational approximations developed for fluid models in the last decade, our discretizations are built…
We develop a Schwinger--Keldysh effective theory for quantum-interference corrections in a two-dimensional electron system in the hydrodynamic regime. Starting from the clean hydrodynamic fixed point, we introduce a minimal random-friction…
Digital quantum simulation relies on Trotterization to discretize time evolution into elementary quantum gates. On current quantum processors with notable gate imperfections, there is a critical tradeoff between improved accuracy for finer…
We study a system of Maxwell's equations that describes the time evolution of electromagnetic fields with an additional electric scalar variable to make the system amenable to a mixed finite element spatial discretization. We demonstrate…
Magnetorotational turbulence provides a viable mechanism for angular momentum transport in accretion disks. We present global, three dimensional (3D), MHD accretion disk simulations that investigate the dependence of the turbulent stresses…
We derive expressions for the local ideal magnetohydrodynamic (MHD) equations for a warped astrophysical disc using a warped shearing box formalism. A perturbation expansion of these equations to first order in the warping amplitude leads…
Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows which cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus…
In the framework of chiral kinetic theory (CKT), we consider a system of right- and left-handed Weyl fermions out of thermal equilibrium in a homogeneous weak magnetic field. We show that the Lorentz invariance implies a modification in the…
We develop a numerical hydrodynamics code using a pseudo-Newtonian formulation that uses the weak field approximation for the geometry, and a generalized source term for the Poisson equation that takes into account relativistic effects. The…