SinoTechIntel Academic Portal
Open AccessDOI: 10.1016/j.ijmst.2026.01.008Original Research

Coupled TM-damage modeling and global sensitivity analysis of thermal spalling in heterogeneous rocks

Liyuan Liu¹,Mingshan Shi¹,Derek Elsworth¹,Tao Wang¹,Hongguang Ji¹,Le Zhang¹,Yaohui Li¹

University of Science and Technology Beijing

Read Executive PreviewQuick FAQ
Coupled TM-damage modeling and global sensitivity analysis of thermal spalling in heterogeneous rocks
Graphical Abstract / Figure
Published In
Academic Research Journal
Published:January 15, 2026Edition:Vol. 32, Issue 1 • pp. 100-112Citation:Liyuan Liu et al. (2026), Academic Research Journal
Impact FactorPeer-Reviewed Core
Sponsored Research Partner

Key Takeaways & Executive Findings

  • • Thermal stress dominates crack initiation and propagation during thermal spalling. • Lateral constraints significantly delay and suppress spalling. • Material heterogeneity markedly influences peak stress and damage modes within a certain range of thermal expansion coefficient. • DGSA effectively quantifies parameter influence and interactions, identifying dominant controls on spalling onset, severity, and damage morphology.
Sponsored Research Highlight

Abstract

Thermal spalling in heterogeneous rocks under rapid heating poses critical risks to deep mining and geothermal operations. In this study, we develop a coupled thermal–mechanical–damage (TM-D) model that explicitly incorporates Weibull distributed heterogeneity to a single fracture in rock, and validate it against ceramic quenching and granite acoustic emission experiments. Distance based generalized sensitivity analysis (DGSA) is applied to quantify the influence and interactions of key parameters, revealing the dominant controls on spalling onset, severity, and damage morphology. The results demonstrate that thermal stress dominates crack initiation and propagation, that lateral constraints can significantly delay and suppress spalling, and that material heterogeneity markedly influences peak stress and damage modes within a certain range of thermal expansion coefficient and has multiple effects on thermal spalling. This study provides a theoretical basis for quantitative assessment and parameter optimization of thermal spalling processes in rock masses.

1. Introduction

Thermal spalling in heterogeneous rocks under rapid heating is a critical hazard in high temperature subsurface engineering, including deep mining, geothermal drilling, and underground fire scenarios [1,2]. Sudden temperature gradients induce complex thermal and mechanical stresses, often leading to rapid crack initiation and surface rupture. While previous studies have advanced thermal–mechanical–damage (TM-D) modeling, most assume material homogeneity and overlook microstructural variability and parameter interactions [3], limiting predictive accuracy and applicability across diverse rock types. Moreover the absence of standardized testing protocols impedes inter study comparisons [4,5]. This study addresses these gaps by integrating a TM-D coupling framework with Weibull distributed heterogeneity and applying distance based generalized sensitivity analysis (DGSA) to quantify the influence and coupling of key parameters.

Modeling thermal spallation in rocks demands an appropriate framework for TM-D coupling. Prior research [4,6] on deep mining induced rock failure has quantified deformation and damage patterns at specific high temperatures and geostress levels, defining yield-damage criteria for homogeneous rock. More recent three-dimensional simulations incorporate material heterogeneity and local degradation via mesoscale renormalization, employing a maximum tensile-stress criterion for tension and a Drucker-Prager criterion for shear failure [7,8]. Nevertheless, analyses of heterogeneous thermal spalling in rocks subjected simultaneously to extreme temperatures and high confining pressure remain scarce.

SinoTechIntel Interactive Document Reader
Page 1–5 of Preview
100%
Download Full PDF

Loading authentic research manuscript (Pages 1–5)...

Sponsored Research Partner
Cite This Research Paper
Liyuan Liu, Mingshan Shi, Derek Elsworth, Tao Wang, Hongguang Ji, Le Zhang, Yaohui Li (2026). Coupled TM-damage modeling and global sensitivity analysis of thermal spalling in heterogeneous rocks. SinoTechIntel Verified Research. https://doi.org/10.1016/j.ijmst.2026.01.008
SinoTechIntel Academic & Legal Disclaimer

Research & Educational Purpose Only:The translations, structured abstracts, analytical annotations, and data reports provided by SinoTechIntel are intended exclusively for academic research, internal corporate R&D, and educational benchmarking. They do not constitute formal engineering, chemical safety, legal, or professional advice.

Copyright & Intellectual Property Notice: Original copyright of the underlying source articles and experimental data remains with the respective authors, institutions, and original publishing journals. SinoTechIntel claims intellectual property only over its proprietary translations, analytical syntheses, and AEO structured enhancements in accordance with international fair use and academic citation principles.

Frequently Asked Questions

What is thermal spalling in rocks?

Thermal spalling is a phenomenon where rapid heating induces thermal stresses in rocks, leading to crack initiation and surface rupture. It is a critical hazard in high-temperature subsurface engineering such as deep mining and geothermal drilling.

How does material heterogeneity affect thermal spalling?

Material heterogeneity, modeled using Weibull distribution, significantly influences peak stress and damage modes. Within a certain range of thermal expansion coefficient, heterogeneity can alter the spalling behavior, showing multiple effects on the process.

What is DGSA and how is it used in this study?

DGSA (Distance-based Generalized Sensitivity Analysis) is a method to quantify the influence and interactions of model parameters. In this study, it is applied to identify dominant controls on spalling onset, severity, and damage morphology, revealing key parameter effects.

What are the practical implications of this research?

The findings provide a theoretical basis for quantitative assessment and parameter optimization of thermal spalling processes, aiding in the design of safer and more efficient deep mining and geothermal operations.

How was the coupled TM-D model validated?

The model was validated against ceramic quenching and granite acoustic emission experiments, ensuring its accuracy in predicting thermal spalling behavior.

Recommended Scientific Literature & Research Partners

Related Technical Papers & Translations

Research Paper
A Novel Approach for Enhanced Brain Tumor Segmentation Using Multimodal MRI and Deep Learning

A Novel Approach for Enhanced Brain Tumor Segmentation Using Multimodal MRI and Deep Learning

Brain tumor segmentation from multimodal MRI is crucial for diagnosis and treatment planning. In this study, we propose a novel deep learning framework that integrates structural and functional imaging modalities to improve segmentation accuracy. Our method employs a multi-scale attention mechanism and a hybrid loss function to handle class imbalance and boundary ambiguity. Evaluated on the BraTS benchmark, our approach achieves state-of-the-art performance, with Dice scores of 0.91, 0.87, and 0.84 for whole tumor, core, and enhancing tumor, respectively. Furthermore, we demonstrate the generalizability of our model across different scanners and protocols. Our findings suggest that the proposed method can significantly aid clinical decision-making and surgical planning.

Read Abstract & PDF
Research Paper
Investigation of coupled acoustic and electrical responses and early warning approaches during re-loading of damaged coal

Investigation of coupled acoustic and electrical responses and early warning approaches during re-loading of damaged coal

Initial damage from engineering disturbances in deep coal mining degrades mechanical properties and heightens dynamic-hazard risks, challenging conventional monitoring. This study probes the coupled acoustic-electrical responses of initially damaged coal under reloading and develops a multi-parameter, multi-level dynamic integrated early-warning model. Using a true-triaxial Split Hopkinson Pressure Bar (SHPB) system, we prepared specimens with graded damage by varying static deviatoric stresses and dynamic impacts. Uniaxial compression reloading was conducted with synchronous acoustic emission (AE) and resistivity monitoring. Joint time-domain responses of force, acoustics, and electricity delineated distinct loading stages. Time-frequency features were extracted via Fourier and wavelet transforms; crack architecture was quantified by 3D AE localization and fractal-dimension analysis. Initial damage markedly reduced load-bearing capacity. Resistivity decreased sharply with increasing deviatoric stress, while cumulative AE counts increased strongly. The AE spectrum evolved from bimodal to broadband with low- and high-frequency enhancement. The resistivity spectrum showed progressive bandwidth broadening, energy amplification, and high-frequency advancement. The AE spatial fractal dimension rose significantly during compaction. An integrated warning system combining multiscale entropy fusion, Temporal Convolutional Network (TCN)-Transformer forecasting, recurrence-network analysis, and a Bayesian framework yielded a 28.4 s lead time, offering a theoretical basis and technical pathway for intelligent prevention of dynamic hazards.

Read Abstract & PDF
Research Paper
Influence of aggregate particle size on fracture behavior and energy evolution of cemented rockfill in the post-peak stage

Influence of aggregate particle size on fracture behavior and energy evolution of cemented rockfill in the post-peak stage

Cemented rockfill (CRF) combines structural support with sustainable reuse of coal-derived solid waste. This study integrates digital image correlation, acoustic emission monitoring, and finite–discrete element simulations to investigate mechanical behavior, fracture development, and energy evolution of CRF containing 54% aggregate content with three grain-size distributions (5–10, 10–20, and 20–30 mm). Results indicate finer aggregates raise compressive strength and elastic modulus, and increase post-peak softening and residual stiffness. Fracture patterns transition from dominantly unidirectional failure in coarse specimens to pronounced X-shaped conjugate shear in fine specimens, with cracks initiating at boundaries and propagating inward. The proportion of failed joints at comparable strains decreases markedly with finer gradation, reflecting a more homogeneous crack network that enhances post-peak load retention and produces frequent minor stress fluctuations. Energy analyses reveal a coarse > medium > fine ordering in cumulative dissipation; however, finer aggregates delay rapid kinetic and dissipative energy release, promoting slower energy redistribution and improved load resistance. These findings quantify how aggregate gradation controls deformational mechanisms, crack topology, and energy partitioning, and provide design guidance for optimizing aggregate size and cementitious composition to enhance ductility, energy absorption, and structural reliability of CRF in underground engineering.

Read Abstract & PDF