Porous geomaterials—from deep clay formations and granitic fracture networks to water-bearing aquifers—exhibit complex, non-linear feedback loops between heat conduction and advection, Darcy fluid seepage, and poroelastic matrix deformation.
Click any vertex or directional coupling pathway in the interactive network below to discover how microscopic pore mechanisms translate into macroscopic continuum field equations and geological phenomena.
The coupled THM behavior is formulated by combining linear momentum balance for the solid skeleton, mass conservation for the compressible pore fluid, and energy conservation for the multi-phase mixture.
Solves the coupled 1D Biot-Terzaghi-Thermal system in real time. Observe how thermal expansion initially triggers rapid surface heave and severe pore overpressure, followed by hydraulic drainage and consolidation settlement.
Interact directly with the porous continuum! Click or drag across the grid to inject heat, pump fluid, apply compressive mechanical loads, or draw impermeable barriers. Watch pore tracer particles respond to non-isothermal Darcy velocity fields.
Fluid injection increases pore pressure ($p$), shifting effective stresses left. Simultaneously, cold water injection induces thermal contraction ($\Delta \sigma_{\text{th}} < 0$), reducing confining stress. When the Mohr circle touches the failure envelope, catastrophic shear rupture occurs!
Tracks the multi-decade coupled evolution of a high-level nuclear waste deposition hole. Canister radioactive heat decays according to isotope half-lives, while water from the host rock infiltrates the compacted bentonite buffer, generating high swelling pressure ($\sim 7\,\text{MPa}$) to seal technological voids.
How THM processes govern modern energy, geo-environmental, and climate resilience projects across the globe.
Repositories in Switzerland (Opalinus Clay), France (Callovo-Oxfordian), and Sweden/Finland (KBS-3 Granite) must verify that canister decay heat does not exceed the illitization threshold ($100^\circ\text{C}$) or induce fracturing by thermal pressurization before bentonite swelling fully seals the excavation damage zone (EDZ).
High-pressure cold water injection into deep crystalline basements (e.g. Utah FORGE, Basel, Soultz-sous-Forêts) causes coupled hydraulic shear stimulation and thermoelastic contraction of fracture walls, enhancing permeability by orders of magnitude while managing induced seismicity risk.
Injecting millions of tons of supercritical $CO_2$ creates localized overpressure plumes and intense thermal cooling near the wellhead (Joule-Thomson cooling). Geomechanical modeling ensures caprock integrity is preserved against hydraulic fracturing or fault reactivation.
Excessive groundwater extraction reduces pore water pressure, transferring overburden stress onto the compressible clayey aquitard matrix ($\Delta \sigma' = -\Delta p$). Inelastic pore collapse has caused over 9 meters of land subsidence in California's San Joaquin Valley and Mexico City.
Compare the characteristic time scales of thermal diffusion ($t_{\text{th}} = L^2 / c_{\text{th}}$) and hydraulic dissipation ($t_h = L^2 / c_v$), and calculate Thermal Péclet numbers across experimental rock types.