sailfish.solvers.srhd_1d

One-dimensional relativistic hydro solver supporting homologous mesh motion.

The solver configuration is:

  • Planar cartesian, or spherical coordinate system in 1d

  • Possible homologous radial expansion in spherical coordinates mode

  • Four conserved quantities: D, Sr, tau, s (scalar lab-frame mass density)

  • Four primitive quantities: rho, ur, p, x (scalar concentration)

  • Gamma-law index of 4/3

The Python code assumes RK2 time stepping, although coefficients are written below for RK1 and low-storage RK3 as well. The C code hard-codes a PLM theta value of 2.0.

Functions

initial_condition(setup, mesh, i0, i1, time, xp)

Classes

Options(compute_wavespeed, rk_order)

Patch(setup, time, conserved, mesh, ...)

Buffers for the solution on a subset of the solution domain.

Physics()

Solver([setup, mesh, time, solution, ...])

Adapter class to drive the srhd_1d C extension module.