• 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.