Related papers: Explicit-Implicit Scheme for Relativistic Radiatio…
Radiation-hydrodynamics (RHD) determines the bulk evolution and observable emission in a wide variety of high-energy astrophysical phenomena. Due to their complexity, RHD problems must usually be studied through numerical simulation. We…
In this paper we report on recent upgrades to our general relativistic radiation magnetohydrodynamics code, Cosmos++, including the development of a new primitive inversion scheme and a hybrid implicit-explicit solver with a more general…
Aims. The importance of radiation to the physical structure of protoplanetary disks cannot be understated. However, protoplanetary disks evolve with time, and so to understand disk evolution and by association, disk structure, one should…
We develop an asymptotic preserving scheme for the gray radiative transfer equation. Two asymptotic regimes are considered: one is a diffusive regime described by a nonlinear diffusion equation for the material temperature; the other is a…
We present a new multi-fluid, multi-temperature plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (non-linear, implicit) scheme for collisional low-temperature plasma. The particle transport is described using the…
We present a novel spectral solver for general relativistic magnetohydrodynamics on dynamical spacetimes. By combining a high order discontinuous spectral method on mapped Chebyshev Fourier grids, our scheme attains exponential convergence.…
In this paper we give an overview of Implicit-Explicit Runge-Kutta schemes applied to hyperbolic systems with stiff relaxation. In particular, we focus on some recent results on the uniform accuracy for hyperbolic systems with stiff…
We present a numerical solver for plasma dynamics simulations in Hall magnetohydrodynamic (HMHD) approximation in one, two and three dimensions. We consider both isotropic and anisotropic thermal pressure cases, where a general gyrotropic…
We present a hierarchical approach for enhancing the robustness of numerical solvers for modelling radiative MHD flows in multi-dimensions. This approach is based on clustering the entries of the global Jacobian in a hierarchical manner…
We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and…
In this paper we propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space…
The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate…
Modelling the interaction between ionizing photons emitted from massive stars and their environment is essential to further our understanding of galactic ecosystems. We present a hybrid Radiation-Hydrodynamics (RHD) scheme that couples an…
A module for the ZEUS-2D code is described which may be used to solve the equations of radiation hydrodynamics to order unity in v/c, in the flux-limited diffusion (FLD) approximation. In this approximation, the tensor Eddington factor f…
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…
This article studies a fully implicit finite volume scheme for transient nonlinear heat transport equations coupled by nonlocal interface conditions modeling diffuse-gray radiation between the surfaces of (both open and closed) cavities.…
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…
To make relevant predictions about observable emission, hydrodynamical simulation codes must employ schemes that account for radiative losses, but the large dimensionality of accurate radiative transfer schemes is often prohibitive.…
When modeling astrophysical fluid flows, it is often appropriate to discard the canonical magnetohydrodynamic approximation thereby freeing the magnetic field to diffuse with respect to the bulk velocity field. As a consequence, however,…
This is the first of two papers examining the critical collapse of spherically symmetric perfect fluids with the equation of state P = (Gamma -1)rho. Here we present the equations of motion and describe a computer code capable of simulating…