European Journal of Mechanics, B/Fluids, vol.115, 2026 (SCI-Expanded, Scopus)
Building upon a modified Karman–Pohlhausen technique, a recent study by Panfilov (2021) employed spherical coordinates to solve the heat transport problem in a heterogeneous domain surrounding a cavity storing cryogenic fluids underground. This analysis revealed the formation of an ice ring around the cavity, acting as a protective barrier against flooding from the stored material. This present work expands on that research by introducing heat generation and absorption into the media, aiming to analyze the temporal evolution of temperature and its impact on ice ring formation. Such heat exchange could be caused by seasonal fluctuations or geothermal activity. Motivated by these real-world influences, we extend the theoretical framework presented in Panfilov (2021) to investigate the universal evolution of the temperature field in the cavity, insulation, and rock regions. This study will track the emergence, persistence (dependent on heat balance), and eventual disappearance of the ice zone while determining its maximum thickness as a function of various parameters. We anticipate that heat generation will accelerate heat transfer between zones, reducing the perturbation length and consequently shortening the lifespan of the ice ring. Conversely, heat absorption will slow down thermal wave propagation by increasing the perturbation time length, thereby prolonging the freezing front of the ice ring and extending the life of both the ice crust and the cryogenic liquid within the underground cavity.