Related papers: A finite volume method for two-moment cosmic-ray h…
Rotationally supported, cold, gaseous disks are ubiquitous in astrophysics and appear in a diverse set of systems, such as protoplanetary disks, accretion disks around black holes, or large spiral galaxies. Capturing the gas dynamics…
We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum,…
In this paper, a second-order accurate method was developed for calculating fluid flows in complex geometries. This method uses cut-Cartesian cell mesh in finite volume framework. Calculus is employed to relate fluxes and gradients along…
In numerical magnetohydrodynamics (MHD), a major challenge is maintaining zero magnetic field-divergence (div-B). Constrained transport (CT) schemes can achieve this at high accuracy, but have generally been restricted to very specific…
We investigate the transport of spectrally resolved cosmic ray (CR) protons with kinetic energies between $1-100$ GeV within the dynamic, multiphase interstellar medium (ISM), using a two-moment CR fluid solver applied to a TIGRESS MHD…
We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++: an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of…
We study the propagation of cosmic rays (CRs) through a simulation of magnetohydrodynamic (MHD) turbulence at unprecedented resolution of $10{,}240^3$. We drive turbulence that is subsonic and super-Alfv\'enic, characterized by $\delta…
We present a numerical implementation of radiative transfer based on an explicitly photon-conserving advection scheme, where radiative fluxes over the cell interfaces of a structured or unstructured mesh are calculated with a second-order…
In this work we develop a class of high-order finite difference weighted essentially non-oscillatory (FD-WENO) schemes for solving the ideal magnetohydrodynamic (MHD) equations in 2D and 3D. The philosophy of this work is to use efficient…
We have developed a new, very efficient numerical scheme to solve the CR diffusion convection equation that can be applied to the study of the nonlinear time evolution of CR modified shocks for arbitrary spatial diffusion properties. The…
We describe a module for the Athena code that solves the gray equations of radiation hydrodynamics (RHD), based on the first two moments of the radiative transfer equation. We use a combination of explicit Godunov methods to advance the gas…
We consider the depth-integrated non-hydrostatic system derived by Yamazaki et al. An efficient formally second-order well-balanced hybrid finite volume finite difference numerical scheme is proposed. The scheme consists of a two-step…
The method of moments is widely used for the reduction of kinetic equations into fluid models. It consists in extracting the moments of the kinetic equation with respect to a velocity variable, but the resulting system is a priori…
Cosmic rays (CRs) have critical impacts in the multiphase interstellar medium (ISM), driving dynamical motions in low-density plasma and modifying the ionization state, temperature, and chemical composition of higher-density atomic and…
This paper develops the genuinely multidimensional HLL Riemann solver for the two-dimensional special relativistic hydrodynamic equations on Cartesian meshes and studies its physical-constraint preserving (PCP) property. Based on the…
Magnetohydrodynamics is an important tool to study the dynamics of plasma Space Physics. In this context, we introduce a three-dimensional magnetohydrodynamic solver with divergence-cleaning in the adaptive multiresolution CARMEN code. The…
Many problems at the forefront of theoretical astrophysics require the treatment of magnetized fluids in dynamical, strongly curved spacetimes. Such problems include the origin of gamma-ray bursts, magnetic braking of differential rotation…
We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2 and 3 spatial dimensions and different systems of coordinates. The code provides a multi-physics, multi-algorithm modular environment particularly oriented…
We consider the three-dimensional magnetohydrodynamics (MHD) equations in the presence of a spatially degenerate stochastic forcing as a model for magnetostrophic turbulence in the Earth's fluid core. We examine the multi-parameter singular…
We describe how analytic solutions for linear hydromagnetic waves can be used for testing cosmological magnetohydrodynamic (MHD) codes. We start from the comoving MHD equations and derive analytic solutions for the amplitude evolution of…