A Numerical Method for Diffusion and Cahn-Hilliard Equations on Evolving Spherical Surfaces

Authors: Jyh-Yang Wu, Sheng-Gwo Chen

Abstract:

In this paper, we present a simple effective numerical geometric method to estimate the divergence of a vector field over a curved surface. The conservation law is an important principle in physics and mathematics. However, many well-known numerical methods for solving diffusion equations do not obey conservation laws. Our presented method in this paper combines the divergence theorem with a generalized finite difference method and obeys the conservation law on discrete closed surfaces. We use the similar method to solve the Cahn-Hilliard equations on evolving spherical surfaces and observe stability results in our numerical simulations.

Keywords: Conservation laws, diffusion equations, Cahn-Hilliard Equations, evolving surfaces.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1127164

References:

[1] M. Burger, L. He and C.-B. Schonlieb, “Cahn-Hilliard inpainting and a generalization for grayvalue images”, SIAM J. Imaging Sciences, vol. 2, 2009, pp. 1129-1167.
[2] J. Cahn, J. Hilliard, "Free energy of a non-uniform system. I. interfacial free energy", J. Chem. Phys vol. 28, 1958, pp. 258–267.
[3] S.-G. Chen and J.-Y. Wu, “Discrete conservation laws on curved surfaces”, SIAM J. Sci. Comput., vol. 35, 2013, pp. A719-A739.
[4] S.-G. Chen and J.-Y. Wu, “Discrete conservation laws on curved surfaces II: A dual approach”, SIAM J. Sci. Comput., vol. 36, 2014, pp. A1813-A1830.
[5] S.-G. Chen, M.-H. Chi, and J.-Y. Wu, “High-Order Algorithms for Laplace–Beltrami Operators and Geometric Invariants over Curved Surfaces”, vol 65, 2015, pp.839-865.
[6] S.-G. Chen and J.-Y. Wu, “Discrete conservation laws on evolving surfaces”, SIAM J. Sci. Comput., vol. 38, 2016, pp. A1725-A1742.
[7] M. DoCarmo, “Differential geometry of curves and surfaces, Prentice-Hall, Lodon, 1976.
[8] M. DoCarmo, “Riemannian geometry”, Birkhauser, Boston, 1992.
[9] I. Dolcetta, S. Vita, R. March, "Area-preserving curve-shortening flows: from phase separation to image processing", Interfaces Free Bound. vol. 4, 2002, pp. 325-343.
[10] Q. Dua, L. Ju and Li Tian, "Finite element approximation of the Cahn–Hilliard equation on surfaces", Comput. Methods Appl. Mech. Engrg., vol. 200, 2011, pp. 2458-2470.
[11] C. M. Elliott and T. Ranner, "Evolving surface finite element method for the Cahn–Hilliard equation", Numer. Math. vol.3, 2015, pp483-534.