3D inverse heat conduction analysis of a heterogeneous rock mass: application to Lascaux Cave
Abstract
The conservation of subterranean cultural heritage sites, particularly those located in karst environments, presents significant challenges due to the complexity of coupled thermal and geological processes. This study addresses the thermal behavior of Lascaux Cave (France), a UNESCO World Heritage site, by developing a 3D heterogeneous thermal model of the surrounding rock mass. The cave is approximately 235 m long, with up to 30 m of overlying rock mass. A geological partitioning was performed using electrical resistivity tomography (ERT) data and machine learning, allowing the rock mass volume to be divided into distinct clusters with homogeneous physical property. An inverse modeling framework, based on the Nelder–Mead simplex algorithm, was employed to estimate the optimal thermal diffusivity of each cluster. An objective function incorporating both phase shift and amplitude damping of periodic thermal signals was introduced to reduce sensitivity to uncertainties in initial temperature conditions. The results demonstrate that heterogeneity in thermal properties significantly influences the spatial distribution of temperatures and improves the accuracy of model predictions compared to conventional homogeneous approaches. This modeling approach provides valuable insight into the cave's thermal regime and establishes a foundation for predicting key conservation-related processes, such as natural convection, condensation, and localized degradation. Future work will focus on integrating the influence of water content and hydrological variations on heat transfer in the rock mass.
Related
- Multiscale Freestanding Vertically Aligned CNT Heat Sink Enabling Multimodal Heat Dissipation by Conduction, Convection, and Radiation
- 3D thermal modeling of Lascaux Cave: Dynamics of air circulation and implications for its microclimate
- The Application of Waste Coal Fly Ash as a Support for a Heterogeneous Catalyst for Biodiesel Production and a Comparison with the Engine Characteristics Using a Homogeneous Catalyst
- Multifunctional LiClO4 for Polymerization Regulation, Ionic Conduction, and Interphase Formation in In Situ Poly(tetrahydrofuran) Gel Polymer Electrolytes
- Correlated one-dimensional electron conduction in ferroelastic domain walls through nematic-like charge stripes
BibTeX
@article{lharti2026inverse,
title = {{3D inverse heat conduction analysis of a heterogeneous rock mass: application to Lascaux Cave}},
author = {Habiba Lharti and Fabien Salmon and Colette Sirieix and Joëlle Riss and Delphine Lacanette},
journal = {International Journal of Thermal Sciences},
year = {2026},
month = oct,
doi = {10.1016/j.ijthermalsci.2026.111445},
url = {https://imt-mines-ales.hal.science/hal-05840981v1},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - 3D inverse heat conduction analysis of a heterogeneous rock mass: application to Lascaux Cave
AU - Habiba Lharti
AU - Fabien Salmon
AU - Colette Sirieix
AU - Joëlle Riss
AU - Delphine Lacanette
PY - 2026
DA - 2026/10/20
JO - International Journal of Thermal Sciences
DO - 10.1016/j.ijthermalsci.2026.111445
UR - https://imt-mines-ales.hal.science/hal-05840981v1
AB - The conservation of subterranean cultural heritage sites, particularly those located in karst environments, presents significant challenges due to the complexity of coupled thermal and geological processes. This study addresses the thermal behavior of Lascaux Cave (France), a UNESCO World Heritage site, by developing a 3D heterogeneous thermal model of the surrounding rock mass. The cave is approximately 235 m long, with up to 30 m of overlying rock mass. A geological partitioning was performed using electrical resistivity tomography (ERT) data and machine learning, allowing the rock mass volume to be divided into distinct clusters with homogeneous physical property. An inverse modeling framework, based on the Nelder–Mead simplex algorithm, was employed to estimate the optimal thermal diffusivity of each cluster. An objective function incorporating both phase shift and amplitude damping of periodic thermal signals was introduced to reduce sensitivity to uncertainties in initial temperature conditions. The results demonstrate that heterogeneity in thermal properties significantly influences the spatial distribution of temperatures and improves the accuracy of model predictions compared to conventional homogeneous approaches. This modeling approach provides valuable insight into the cave's thermal regime and establishes a foundation for predicting key conservation-related processes, such as natural convection, condensation, and localized degradation. Future work will focus on integrating the influence of water content and hydrological variations on heat transfer in the rock mass.
SN - 1290-0729
ER -