Related papers: GENESIS: A high-resolution code for 3D relativisti…
We construct a new relativistic viscous hydrodynamics code optimized in the Milne coordinates. We split the conservation equations into an ideal part and a viscous part, using the Strang spitting method. In the code a Riemann solver based…
We describe the CRASH (Center for Radiative Shock Hydrodynamics) code, a block adaptive mesh code for multi-material radiation hydrodynamics. The implementation solves the radiation diffusion model with the gray or multigroup method and…
Many astrophysical systems can only be accurately modelled when the behaviour of their baryonic gas components is well understood. The residual distribution (RD) family of partial differential equation (PDE) solvers produce approximate…
Current GPU-accelerated supercomputers promise to enable large-scale simulations of turbulent flows. Lattice Boltzmann Methods (LBM) are particularly well-suited to fulfilling this promise due to their intrinsic compatibility with highly…
With increasing engineering demands, there need high order accurate schemes embedded with precise physical information in order to capture delicate small scale structures and strong waves with correct "physics". There are two families of…
GenASiS (General Astrophysical Simulation System) is a new code being developed initially and primarily, though by no means exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This…
We have built a code to obtain the exact solutions of Riemann problems in ideal magnetohydrodynamics (MHD) for an arbitrary initial condition. The code can handle not only regular waves but also switch-on/off rarefactions and all types of…
The paper describes a new upwind conservative numerical scheme for special relativistic resistive magnetohydrodynamics with scalar resistivity. The magnetic field is kept approximately divergence free and the divergence of the electric…
We adopt a two-moment formalism, together with a reference-metric approach, to express the equations of relativistic radiation hydrodynamics in a form that is well-suited for numerical implementations in curvilinear coordinates. We…
We implement advanced Riemann solvers HLLC and HLLD \cite{Mignone:2005ft,MUB:2009} together with an advanced constrained transport scheme \cite{Gardiner:2007nc} in a numerical-relativity neutrino-radiation magnetohydrodynamics code. We…
A two dimensional hydrochemical hybrid code, KM2, is constructed to deal with astrophysical problems that would require coupled hydrodynamical and chemical evolution. The code assumes axisymmetry in cylindrical coordinate system, and…
Computation of gas dispersal in urban places or hilly grounds requires a large amount of computer time when addressed with conventional multidimensional models. Those are usually based on two-phase flow or Navier-Stokes equations. Different…
We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary…
We present a five-wave Riemann solver for the equations of ideal relativistic magnetohydrodynamics. Our solver can be regarded as a relativistic extension of the five-wave HLLD Riemann solver initially developed by Miyoshi and Kusano for…
Many systems of current interest in relativistic astrophysics require a knowledge of radiative transfer in a magnetized gas flowing in a strongly-curved, dynamical spacetime. Such systems include coalescing compact binaries containing…
The equations of relativistic hydrodynamics are transformed so that steps forward in time preserves local simultaneity. In these variables, the space-time coordinates of neighboring points on the mesh are simultaneous according to co-moving…
We show that in hydrodynamic simulations for relativistic heavy-ion collisions, strong resonance decay calculations can be performed with fewer species of particle resonances while preserving good accuracy in single particle spectra and…
The characteristic decomposition for GRMHD in the comoving frame of the fluid has been known for a long time. However, it has not been known in the coordinate frame of the simulation and in terms of the conserved variables evolved in…
Heavy-ion collisions at center-of-mass energies between 1 and 100 GeV/nucleon are essential to understand the phase diagram of QCD and search for its critical point. At these energies the net baryon density of the system can be high, and…
A cosmological hydrodynamic code is described, which includes a routine to compute cosmic ray acceleration and transport in a simplified way. The routine was designed to follow explicitly diffusive acceleration at shocks, and second-order…