Cracks Learn from Their Mistakes ✛ ✎

Original: "A Phase-Field Study of Desiccation Crack Pattern Maturation under Drying-Wetting Cycles"
· Hiroki Fukushima, Satoshi Yukawa
arXiv:2609.11048v1 · 2026-09-10 · arXiv non-exclusive · independent analysis · 3 min
Repeated cycles of drying and rewetting make cracks in paste bend at 120-degree angles.

Context

This work sits somewhat off to the side: nearby in the archive are studies of collective behavior, stochastic modeling, and the properties of water, but there are almost no direct neighbors on the topic of drying pastes. The closest work concerns the fracture of materials such as plastic and steel, as well as the statistics of cracks in tissues. This paper fills a gap: it is the first to combine a phase-field fracture model with wetting and drying cycles to explain why crack intersection angles change over time.

Methods

The authors based their work on a phase-field fracture model, in which a crack is described not as a sharp line but as a diffuse damage field. To this model they added two new elements: crack healing during wetting and 'scars' — weakened zones that remain after previous cracks. The paste was represented as a thin elastic layer on a rigid substrate, and drying as a gradual reduction in volume. Material heterogeneity was set by random strength fields that persist from cycle to cycle. The model was solved using the finite element method, repeating drying and wetting cycles up to 50 times.

Results

Already after the first cycle, the distribution of crack intersection angles has a peak near 90 degrees, corresponding to T-shaped intersections. With each subsequent cycle, the fraction of angles near 120 degrees, that is, Y-shaped intersections, grows. The standard deviation of the angle from 120 degrees decreases exponentially with a characteristic time of 2.85 cycles, which agrees with experiments where the relaxation time is about 4 cycles. However, the peak at 90 degrees in the model barely decreases, indicating that some T-shaped intersections persist. Crack energy dominates the total energy of the system and also decreases exponentially with a similar relaxation time. The main contribution to the energy reduction comes from the shortening of the effective crack length, not from changes in their strength.

Implications

This work differs from its neighbors in that it is the first to reproduce the maturation of a crack pattern within a continuum model that explicitly accounts for wetting and drying cycles. Previous models either considered only a single drying cycle or used simplified spring networks that did not allow for energy analysis. Here, it was possible to show that geometric rearrangement and energy relaxation go hand in hand, and that the key parameter is crack length.

Future development

The model can be improved by describing scars from past cracks more realistically. Currently, weakening is introduced only from the second cycle onward, which creates an artificial jump between the first and second cycles. If a way is found to introduce scars smoothly, the model will come closer to experiment. It would also be worth tracking the fate of individual intersections to understand why some T-shaped nodes persist.

Impact

The results are important for understanding crack formation in drying soils, ceramics, and other materials subject to cyclic wetting and drying. They may also help in geology, in the study of columnar jointing in lava.

Next steps

The next step will be to identify and track individual crack intersections in the model, as well as to study the dependence on parameters that control material weakening. This will help clarify why some T-shaped intersections do not turn into Y-shaped ones.

Key open problems

This work directly addresses the unsolved problem of how memory of past deformations influences pattern formation in disordered media. In fracture mechanics, the criteria for choosing the crack intersection angle under cyclic loading remain an open question.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

You built a phase-field fracture model with healing and residual scars to reproduce how crack intersection angles change over repeated wetting and drying cycles. The model shows angles approaching 120° exponentially, with a characteristic time of about 3 cycles, and links this to a reduction in effective crack length.

Where its strength is

The strength of this work lies in its explicit treatment of cyclicity and its preservation of damage history, which makes it possible to trace the evolution of the pattern and connect the geometric rearrangement to energetic relaxation.

Where it could go next

Worth trying

  • Test how the characteristic relaxation time changes when the correlation length of the random strength field is varied.
  • Use topological invariants from the vertex model of tissue to quantify the maturity of the crack network.
  • Compare the energetic contributions under different drying regimes (for example, slow and fast) with the predictions of a plastic degradation model.

Why this may matter

Understanding how cracks rearrange under repeated cycles could help predict the durability of materials subject to periodic wetting and drying, and could also inform the control of crack patterns in industrial processes.

In experiments with coffee paste, after several cycles of wetting and drying the cracks line up at 120-degree angles, like a honeycomb. This is the angle at which three cracks can meet without creating excessive stress.

Key numbers

  • characteristic relaxation time of the angle: 2.85 cycles
  • characteristic relaxation time of the energy: 3.17 cycles
  • system: 100 × 100 mm
  • number of cycles in the simulation: up to 50
  • number of samples for statistics: 20
Concepts
charge cycle fracture Atlantic multidecadal oscillation scattering angle shear elasticity decay time damping rate