Formulating bulk viscosity for neutron star simulations | Synapse