An efficient threshold dynamics method for wetting on rough surfaces
The threshold dynamics method developed by Merriman, Bence and Osher(MBO) is an efficent method for simulating the motion by mean curvature flow when the interface is away from the solid boundary. Direct generalization of the MBO type method to the wetting problems with interface intersecting the solid boundary is not easy because solving heat equation on general domain with wetting boundary condition is not as efficient as that for the original MBO method. The dynamics of the contact point also follows a different dynamic law compared to interface dynamics away from the boundary. In this paper, we develop an efficient volume preserving threshold dynamics method for wetting on rough surfaces, which is based on minimization of the weighted surface area functional over a extended domain that includes the solid phase. The method is simple, stable with the complexity $O(N \log N)$ per time step and it is not sensitive to the inhomogeneity or roughness of the solid boundary.