Braginskii Viscosity in Cosmological Simulations of Galaxy Clusters: Implementation, Validation, and First Application

Abstract

We present the implementation of an anisotropic viscosity solver within the magnetohydrodynamics (MHD) framework of the TreeSPH code OpenGadget3. The solver models anisotropic viscous transport along magnetic field lines following the Braginskii formulation and includes physically motivated limiters based on the mirror and firehose instability thresholds, which constrain the viscous stress in weakly collisional plasmas. To validate the implementation, we performed a suite of standard test problems—including two variants of the sound wave test, circularly and linearly polarized Alfvén waves, fast magnetosonic wave, and the Kelvin–Helmholtz instability—both with and without the plasma-instability limiters. The results show excellent agreement with the Arepo implementation of a similar anisotropic viscosity model, confirming the accuracy and robustness of our method. Our formulation integrates seamlessly within the individual adaptive timestepping framework of OpenGadget3, avoiding the need for subcycling. This provides efficient and stable time integration while maintaining physical consistency. Finally, we applied the new solver to a cosmological zoom-in simulation of a galaxy cluster as a proof-of-concept application, demonstrating its capability to model anisotropic transport and plasma microphysics in realistic large-scale environments. Our implementation offers a versatile and computationally efficient tool for studying anisotropic viscosity in magnetized astrophysical systems.