A new contact line treatment for a conservative level set method

As is well known, for PDEs that enjoy a conservation or dissipation property, numerical schemes that inherit this property are often advantageous in that the schemes are fairly stable and give qualitatively better numerical solutions in practice. Lately, Furihata and Matsuo have developed the so-called “discrete variational derivative method” that automatically constructs energy preserving or dissipative finite difference schemes. Although this method was originally developed on uniform meshes, the use of non-uniform meshes is of importance for multi-dimensional problems. On the other hand, the theories of discrete differential forms have received much attention recently. These theories provide a discrete analogue of the vector calculus on general meshes. In this paper, we show that the discrete variational derivative method and the discrete differential forms by Bochev and Hyman can be combined. Applications to the Cahn–Hilliard equation and the Klein–Gordon equation on triangular meshes are provided as demonstrations. We also show that the schemes for these equations are H1-stable under some assumptions. In particular, one for the nonlinear Klein–Gordon equation is obtained by combination of the energy conservation property and the discrete Poincaré inequality, which are the temporal and spacial structures that are preserved by the above methods.