Heat flux along z-axis normalized to the 300 K value, which shows the presence of the xenon bubble and its effect on the local heat flux. Heat flux steadily decreases with increasing temperature as phonon transport becomes gradually more diffuse. This simulation was conducted using $ΛBates$⁠.