Energy tunnels are typically comprised of tunnel linings that are equipped with heat transfer pipes, thereby converting them into tunnel ground heat exchangers (tunnel GHEs). This innovative technology offers renewable methods of space heating, cooling and hot water production. While still relatively uncommon, it is becoming increasingly popular due to the growing interest in sustainability and innovation, particularly in large infrastructure projects [8, 7, 10]. Additionally, the increasing need for clean and renewable sources of energy has made energy tunnels an attractive solution for addressing the global climate challenge [7, 11, 1].
 The geotechnical behaviour of energy tunnels that subject the soil to thermal loads is still not fully understood. Numerical analyses conducted in the past use simple soil constitutive models, such as the linear elastic model [2, 3, 9] and elastic-perfectly plastic Mohr-Coulomb model [6, 4, 5]. Using simple constitutive models for tunnel design may not accurately depict the complex behaviour of soil materials under different conditions, and the complex interplay among the deformation, flow and thermal models governing the soil behaviour. This can lead to calculated stresses and strains that differ from actual values, resulting in potentially unsafe designs [12, 13].
 Due to their multi-phase nature, saturated soils exhibit both reversible and irreversible thermal volumetric changes. The reversible volumetric change, resulting from thermal expansion/contraction, is typically assumed to occur within the elastic thermo-mechanical response of the soil. Conversely, the irreversible change is more appropriately treated within the elasto-plastic response regime.
 In addition to volumetric changes, thermal effects can impact the stiffness and strength of soil, ultimately influencing the load transfer mechanisms and stress distributions within the system. These changes, in turn, can lead to additional deformations in both the tunnel and the surrounding soil. It should be noted that simply applying proportional stress paths to any of the Mohr-Coulomb models will result in solely elastic responses and fail to reproduce the irrecoverability of the volumetric response. Therefore, the incorporation of an elasto-plastic framework for analysing thermo-mechanical coupling is required to obtain more accurate predictions with respect to changes in effective stress and ground displacement.
 In this study, the coupled thermo-mechanical behaviour of energy tunnels is investigated employing a thermo-elasto-plastic constitutive model developed in the framework of the critical state soil mechanics and considering incorporating thermally induced volumetric change in the model that involves adding a thermal strain component to the model. This is achieved by including a temperature-dependent parameter that accounts for the volumetric change due to change in temperature. The model is implemented in the finite element package COMSOL Multiphysics and used to study the thermo-mechanical behaviour of an energy tunnel in 2D plane strain condition (see Figure 1). The in-situ state of stress is reproduced, and the construction sequence is simulated (soil-structure convergence). The analysis focuses on the stresses and deformations that develop within and around energy tunnels due to the influence of geothermal operation.
 The proposed model is validated following several steps. First, the mechanical model is calibrated against field data from tunnels constructed in soft soils [10]. Then, thermo-mechanical effects are validated using existing monitoring data from real-scale tunnels available in literature [3,4,10,14]. For further validation, model estimations are compared with results from small-scale thermal-mechanical tests, such as thermal-triaxial tests. This comprehensive validation strategy thus ensures a rigorous and reliable assessment of our model’s performance in replicating soil behaviour near energy tunnels.