A structure-preserving method of fundamental solutions for the multi-phase Mullins–Sekerka flow

A charge simulation method, a variant of the method of fundamental solutions, is applied to approximate the multi-phase Mullins–Sekerka flow in R 2 and in a half-plane H bounded by a Neumann wall W . In the underlying mathematical model, interfaces driven by their curvature are coupled through a harmonic chemical-potential field. Each chemical potential is represented by fundamental solutions centered at charge points off the curve, so no bulk mesh is required. The scheme treats curve networks separating several phases at triple junctions, including phases that occupy more than one region; on the half-plane boundary, the homogeneous Neumann condition is imposed exactly by image charges, and mobile contacts stay orthogonal to the wall. The discretization is structure-preserving in the sense that every bounded phase area is conserved to machine precision at the velocity level by a null-space projection of the discrete area constraints. The proposed scheme is assessed through a convergence test against an exact three-concentric-circle solution.