This is a video of the zero-temperature Glauber dynamics of the Ising model in the 3-dimensional cubic lattice, with “hexagonal” boundary conditions (which is equivalent to the usual flip dynamics on rhombus tilings of the hexagon). It converges to an asymptotic shape exhibiting an arctic circle.

The first version is of size 50 in each direction:

Here is a larger version (side length 100):