This paper presents the implementation of the eXtended Finite Element Method (XFEM) in the general-purpose commercial software package COMSOL Multiphysics for multi-field thermo_x0002_hydromechanical problems in discontinuous porous media. To this end, an exclusive enrichment strategy is proposed in compliance with the COMSOL modelling structure. COMSOL modules and physics interfaces are adopted to take account of the relevant physical processes involved in thermo_x0002_hydromechanical coupling analysis, namely: the mechanical deformation, fluid flow in porous media and heat transfer. The model preprocessing, level-set-updates, coupling of the relevant physics and postprocessing procedures are performed adopting a coherent utilisation of the COMSOL’s built-in features along with the COMSOL’s LiveLink for MATLAB functions. The capabilities and performance of the proposed approach are investigated by examining several multi-field thermo_x0002_hydro-mechanical simulations with reference to the block caving process in mines.