Dendritic Solidification: How Metal Freezes Into Trees

Freezing is a race between heat and instability

When a pure liquid metal or alloy solidifies, the interface between solid and liquid rarely stays flat. A small random bump on the freezing front sticks slightly further out into the undercooled liquid, where the driving force for solidification is stronger, so it tends to grow even faster than its surroundings — an instability that, left unchecked, would run away immediately. What actually reins it in is surface tension, which penalises sharply curved interfaces and damps the smallest, tightest bumps. The competition between an unstable driving force at long wavelengths and a stabilising surface-tension penalty at short wavelengths selects a characteristic bump spacing, and repeated branching of that instability is what produces the tree-like dendrite — the classic snowflake-like shape seen in cast metal microstructure and, of course, in actual snowflakes.

The phase-field method: tracking an interface without tracking it

Directly tracking a solidification front that branches, merges and grows increasingly complex geometry is numerically brutal. The phase-field method sidesteps the problem by introducing a smooth order parameter φ that is 1 deep in the solid, 0 deep in the liquid, and transitions continuously across a thin diffuse interface — the sharp boundary is replaced by a steep but continuous gradient, and the interface's location and motion emerge automatically from solving a PDE for φ rather than being tracked explicitly:

∂φ/∂t = M · [ ∇²φ - (∂f/∂φ) ]  +  coupling to the temperature field
∂T/∂t = α∇²T + (L/cp) · (∂φ/∂t)

φ  = phase-field order parameter (0 = liquid, 1 = solid)
M  = interface mobility
f(φ) = double-well free-energy density with minima at φ=0 and φ=1
T  = temperature, α = thermal diffusivity
L  = latent heat of fusion, cp = specific heat

The φ equation drives the interface toward one of the two stable phases while the coupled temperature equation releases latent heat exactly where the interface advances, which is what makes the model self-consistently thermal rather than just geometric.

Undercooling: the fuel for growth

Undercooling — how far below the equilibrium freezing point the liquid actually sits — is the driving force for the whole process; a liquid sitting exactly at its melting point has no thermodynamic incentive to solidify at all. Larger undercooling means a bigger free-energy difference between liquid and solid, which means faster interface velocity and, because faster growth outruns the diffusion of latent heat away from the tip, generally finer, more branched dendrite structure. Modest undercooling produces coarse, slowly growing dendrites with widely spaced branches; deep undercooling produces fast, densely branched structures — a relationship metallurgists exploit deliberately, since dendrite arm spacing directly controls the mechanical properties of a cast part.

Anisotropy: why dendrites pick preferred directions

Real dendrites do not branch in arbitrary directions — in a cubic metal like nickel or iron they overwhelmingly grow along the crystal's ⟨100⟩ axes, producing the familiar cross-shaped snowflake symmetry at the microstructural level. This comes from anisotropic surface energy: the energy cost of the solid-liquid interface is not the same in every crystallographic direction, and the instability analysis that selects branch spacing also, in an anisotropic crystal, selects the specific directions where growth is energetically cheapest. Phase-field models capture this by making the interfacial energy and mobility terms depend on the local interface orientation relative to the crystal lattice, usually with a term like 1 + δ·cos(kθ) where k matches the crystal's rotational symmetry (k=4 for cubic crystals) and δ sets how strongly the anisotropy biases growth direction.

Why casters actually care

Dendrite arm spacing set during solidification is one of the strongest single predictors of a cast alloy's final mechanical properties — finer spacing generally means less microsegregation and better strength and ductility after heat treatment. Because that spacing is controlled by cooling rate and undercooling, foundries and continuous-casting operations tune solidification rate deliberately, and phase-field simulation is now a standard predictive tool for choosing casting parameters before committing to an expensive physical trial.

Frequently asked questions

Try it live

Everything above runs in your browser — open Crystal Growth 3D and change the parameters while it is running. Nothing is installed, nothing is uploaded, the whole model lives in one tab.