: Quantum impurity models (QIMs) are ubiquitous throughout physics. As simplified toy models they provide crucial insights for understanding more complicated strongly correlated systems, while in their own right are accurate descriptions of many experimental platforms. In equilibrium, their physics is well understood and have proven a testing ground for many powerful theoretical tools, both numerical and analytical, in use today. Their nonequilibrium physics is much less studied and understood. However, the recent advancements in nonequilibrium integrable quantum systems through the development of generalized hydrodynamics (GHD) coupled with the fact that many archetypal QIMs are in fact integrable presents an enticing opportunity to enhance our understanding of these systems. We take a step towards this by expanding the framework of GHD to incorporate integrable interacting QIMs. We present a set of Bethe-Boltzmann type equations which incorporate the effects of impurity scattering and discuss the new aspects which include entropy production. These impurity GHD equations are then used to study a bipartioning quench wherein a relevant backscattering impurity is included at the location of the bipartition. The density and current profiles are studied as a function of the impurity strength and expressions for the entanglement entropy and full counting statistics are derived.
