Self-similarity and scaling of thermal shock fractures
  The problem of crack pattern formation due to thermal shock loading at the surface of half-space is solved numerically using two-dimensional boundary element method. The results of numerical simulations with 100-200 random simultaneously growing and interacting cracks are used to obtain scaling relations for crack length and spacing. The numerical results predict that such process of pattern formation with quasi-static crack growth is not stable and at some point the excess energy leads to unstable propagation of one of the longest crack. The onset of instability has also been determined from numerical results. 
Sergejs Tarasovs and Ahmad Ghassemi
cond-mat.mtrl-sci cond-mat.soft