Variational Dynamics of Free Triple Junctions

We propose a level set framework for representing triple junctions either with or without free endpoints. For triple junctions without free endpoints, our method uses two level set functions to represent the three segments that constitute the structure. For free triple junctions, we extend our method using the free curve work of Schaeffer and Vese (J Math Imaging Vis, 1–17, 2013), Smereka (Phys D Nonlinear Phenom 138(3–4):282–301, 2000). For curves moving under length minimizing flows, it is well known that the endpoints either intersect perpendicularly to the boundary, do not intersect the boundary of the domain or the curve itself (free endpoints), or meet at triple junctions. Although many of these cases can be formulated within the level set framework, the case of free triple junctions does not appear in the literature. Therefore, the proposed free triple junction formulation completes the important curve structure representations within the level set framework. We derive an evolution equation for the dynamics of the triple junction under length and area minimizing flow. The resulting system of partial differential equations are both coupled and highly non-linear, so the system is solved numerically using the Sobolev preconditioned descent. Qualitative numerical experiments are presented on various triple junction and free triple junction configurations, as well as an example with a quadruple junction instability. Quantitative results show convergence of the preconditioned algorithm to the correct solutions.