Resulting density distributions of the numerical simulations for different Milky Way models and TDE rates. The snapshots are taken at the end of the simulations when the forward shocks reach the size of the EBs and show the detailed bubble structure, including the elapsed time since the first TDE. The different subfigures correspond to the results of the exponential halo model with energy and mass injections every 100–1000 yr (top left), the β-model halo with TDEs every 10 kyr (top right), the β-model halo with TDEs every 100 kyr (bottom left), and the multicomponent model with TDEs every 100 kyr (bottom right).