sailfish.solvers.srhd_2d

Two-dimensional relativistic hydro solver supporting homologous mesh motion.

The solver configuration is:

  • Spherical polar coordinates in 2d

  • Possible homologous radial expansion

  • Four conserved quantities: D, Sr, Sq, tau

  • Four primitive quantities: rho, ur, uq, p

  • 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, j0, ...)

Classes

Options(compute_wavespeed, rk_order, ...)

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

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

Physics(jet_mdot, jet_gamma_beta, jet_theta, ...)

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

Adapter class to drive the srhd_1d C extension module.