Differential Growth 3

I recently squashed a few bugs in the differential growth simulation that prevented me from sharing the source earlier. It was satisfying. With that out of the way, here’s a link to the Grasshopper definition (note that it uses the current release of SpatialSlur.dll).

Growth density variation

As mentioned in the previous post, the simulation uses a half-edge data structure to handle refinement of the surface mesh as it grows. When an edge exceeds a specified maximum length, it splits (along with the triangles on either side). Edges are also spun to unify vertex valence throughout the mesh, resulting in well-proportioned triangles that are ideal for all sorts of things. The video above shows the effect of different maximum edge lengths on the growth process (from large to small).

The half-edge mesh is particularly well-suited to this process because it allows such local topological modifications to be executed in constant time. In other words, splitting an edge will always take the same amount of time regardless of how many other edges are in the mesh. This property is especially important here since we’re working with an algorithm that tends to dramatically increase its own input size from one iteration to the next. Without constant-time topological operators, things would quickly get out of hand.

Frame_0250
Growth 2 iteration 250
Frame_0350
Growth 2 iteration 350
Frame_0450
Growth 2 iteration 450

Growing Dendrites

Since stumbling upon this paper by Ryo Kobayashi, I’ve been spending a shameful amount of time growing dendrites on my computer.

Isotropic dendritic growth

The process shown in the video above makes use of a phase field model to track the evolution of the interface between a supercooled liquid and a solid. While this seemed like unfamiliar territory at first, a quick look under the hood revealed a set of reaction-diffusion-like equations. The system is represented by two co-dependent scalar fields – one representing temperature and one representing the phase. The latter is either 0 (liquid) or 1 (solid) everywhere in the system except at the interface where it grades smoothly between the two extremes. At each step, the current temperature contributes to the change in phase and the change in phase contributes to the change in temperature. Over time, this feedback leads to instabilities in the interface which result in all the fancy branching patterns.

The Grasshopper definition is available here. Note that it uses the current release of SpatialSlur.dll which can be downloaded here.

d_02
Isotropic dendritic growth 1
a_04
Isotropic dendritic growth 2
b_03_2
Isotropic dendritic growth 3
b_04_2
Isotropic dendritic growth 4