Numerical Modeling of Heat Transfer and Thermal Attenuation Lengths in Ventilated Caves

Atmospheric temperature variations are characterized by different frequency-modes including yearly and daily fluctuations as well as yearly average temperatures. In a ventilated cave, the thermal amplitude is attenuated with increasing distance from the upper and/or lower cave entrance. Convective heat flux for dry air is exchanged with cave wall and the heat transferred by conduction to the surrounding rock is attenuated within a certain distance. Here, we aim at determining the thermal attenuation length for the different modes along a cave and the surrounding rock. Distribution of amplitude of fluctuations along the cave is specified by using Fourier series expansion for dominant modes. In next steps, the effect of latent heat exchange at the cave wall due to evaporation or condensation will be studied. The main variables controlling thermal perturbations in ventilated caves are recognized and compared with field data by developing a numerical model based on the heat and mass transfer between the rock and the air. Moreover, our model aims to quantify the air and wall temperature profile along the conduit for a one-year. The model calculates water-vapor exchange rate along the cave showing the amount of consumed or produced water by evaporation or condensation along the cave and provides the rock temperature distribution in the surrounding hostrock. Our preliminary results contribute to a better understanding of the long-term cave dynamics and may support a quantitative interpretation of speleothem records.