Conservative finite volume scheme for BDNK relativistic dissipative hydrodynamics*

We present the first conservative finite volume numerical scheme for the equations of BDNK theory, a promising framework for incorporating dissipative effects such as viscosity and heat conduction into relativistic fluid models. Such effects have been shown to be important in modeling the quark-gluon plasma produced in heavy-ion collisions, and recent studies have suggested that they may also be needed to model binary neutron star mergers at the level of accuracy required for next-generation gravitational wave detectors.  Our scheme serves as an early step toward the application of BDNK theory to systems such as these, and has been designed to possess the same enhanced accuracy and stability properties typically required of relativistic ideal fluid solvers. In this talk, I briefly review BDNK theory, outline the construction of the scheme, highlight its behavior in a suite of multidimensional test problems, and comment on directions for future work.