A structure-preserving finite element method for the multi-phase Mullins-Sekerka problem with triple junctions
Tokuhiro Eto, Harald Garcke, Robert Nürnberg
TL;DR
The paper addresses numerical approximation of sharp-interface multi-phase Mullins-Sekerka flow with triple junctions by introducing a variational, fully discrete parametric finite element method (BGN-type). The scheme is unconditionally stable and preserves two structural properties critical to the model: the weighted interfacial length dissipates as energy decreases, and the enclosed phase areas are exactly conserved, with a nonlinear variant achieving exact area preservation at the fully discrete level. A key feature is the tangential motion of vertices, which yields asymptotically equidistributed meshes and eliminates the need for remeshing. The authors generalize the approach to arbitrary multi-phase systems, provide convergence tests, and showcase extensive 3- and 4-phase evolutions under equal and unequal surface energies, including topology-changing scenarios, demonstrating robustness and practical applicability for complex curve networks.
Abstract
We consider a sharp interface formulation for the multi-phase Mullins-Sekerka flow. The flow is characterized by a network of curves evolving such that the total surface energy of the curves is reduced, while the areas of the enclosed phases are conserved. Making use of a variational formulation, we introduce a fully discrete finite element method. Our discretization features a parametric approximation of the moving interfaces that is independent of the discretization used for the equations in the bulk. The scheme can be shown to be unconditionally stable and to satisfy an exact volume conservation property. Moreover, an inherent tangential velocity for the vertices on the discrete curves leads to asymptotically equidistributed vertices, meaning no remeshing is necessary in practice. Several numerical examples, including a convergence experiment for the three-phase Mullins-Sekerka flow, demonstrate the capabilities of the introduced method.
