Numerical modeling of bulk viscosity in neutron stars | Synapse