Related papers: Explicit-Implicit Scheme for Relativistic Radiatio…
A new method for solving relativistic ideal hydrodynamics in (1+3)D is developed. Longitudinal and transverse radial flows are explicitly embedded and the hydrodynamic equations are reduced to a single equation for the transverse velocity…
In this thesis different numerical methods, as well as applications of the methods to a number of current problems in relativistic astrophysics, are presented. In the first part the theoretical foundation and numerical implementation of a…
We describe a numerical scheme for studying time-dependent, multifluid, magnetohydrodynamic shock waves in weakly ionized interstellar clouds and cores. Shocks are modeled as propagating perpendicular to the magnetic field and consist of a…
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…
In this paper, we consider numerical approximations for the viscous Cahn-Hilliard equation with hyperbolic relaxation. This type of equations processes energy-dissipative structure. The main challenge in solving such a diffusive system…
Upcoming large-scale-structure surveys can shed new light on the properties of dark energy. In particular, if dark energy is a dynamical component, it must have spatial perturbations. Their behaviour is regulated by the speed of sound…
Time-implicit schemes are attractive since they allow numerical time steps that are much larger than those permitted by the Courant-Friedrich-Lewy criterion characterizing time-explicit methods. This advantage comes, however, with a cost:…
Active galactic nuclei, x-ray binaries, pulsars, and gamma-ray bursts are all believed to be powered by compact objects surrounded by relativistic plasma flows driving phenomena such as accretion, winds, and jets. These flows are often…
We describe a numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference scheme on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a…
This paper develops the genuinely multidimensional HLL Riemann solver and finite volume scheme for the two-dimensional special relativistic hydrodynamic equations on Cartesian meshes and studies its physical-constraint-preserving (PCP)…
In this paper a new semi-implicit relaxation scheme for the simulation of multi-scale hyperbolic conservation laws based on a Jin-Xin relaxation approach is presented. It is based on the splitting of the flux function into two or more…
We report on recent upgrades to our general relativistic radiation-magnetohydrodynamics code, Cosmos++, which expands the two-moment, $\bf{M}_1$, radiation treatment from grey to multi-frequency transport, including Doppler and…
We present algorithms to solve relativistic hydrodynamics in 3+1--dimensional situations without apparent symmetry to simplify the solution. In simulations of heavy--ion collisions, these numerical schemes have to deal with the physical…
The description of general relativistic radiation hydrodynamics in spherical symmetry is presented in natural coordinate choices. For hydrodynamics, comoving coordinates are chosen, and the momentum phase space for the radiation particles…
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with…
For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we…
We describe compressible two-phase flows by a single-velocity six-equation flow model, which is composed of the phasic mass and total energy equations, one volume fraction equation, and the mixture momentum equation. The model contains…
We present a new explicit and stable numerical algorithm to solve the homogeneous heat equation. We illustrate the performance of the new method in the cases of two 2D systems with highly inhomogeneous random parameters. Spatial…
This paper presents approximation methods for time-dependent thermal radiative transfer problems in high energy density physics. It is based on the multilevel quasidiffusion method defined by the high-order radiative transfer equation (RTE)…
We present an analytic solution of the Riemann problem for the equations of relativistic hydrodynamics with the ultra-relativistic equation of state and non-zero tangential velocities. A 3-dimensional numerical code solving such equations…