TY - JOUR
T1 - Direct Discretization Method for the Cahn–Hilliard Equation on an Evolving Surface
AU - Li, Yibao
AU - Qi, Xuelin
AU - Kim, Junseok
N1 - Funding Information:
Acknowledgements Y.B. Li is supported by National Natural Science Foundation of China(Nos. 11601416, 11631012). The corresponding author (J.S. Kim) was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2016R1D1A1B03933243). The authors thank the reviewers for the constructive and helpful comments concerning the revision of this article.
PY - 2018/11/1
Y1 - 2018/11/1
N2 - We propose a simple and efficient direct discretization scheme for solving the Cahn–Hilliard (CH) equation on an evolving surface. By using a conservation law and transport formulae, we derive the CH equation on evolving surfaces. An evolving surface is discretized using an unstructured triangular mesh. The discrete CH equation is defined on the surface mesh and its dual surface polygonal tessellation. The evolving triangular surfaces are then realized by moving the surface nodes according to a given velocity field. The proposed scheme is based on the Crank–Nicolson scheme and a linearly stabilized splitting scheme. The scheme is second-order accurate, with respect to both space and time. The resulting system of discrete equations is easy to implement, and is solved by using an efficient biconjugate gradient stabilized method. Several numerical experiments are presented to demonstrate the performance and effectiveness of the proposed numerical scheme.
AB - We propose a simple and efficient direct discretization scheme for solving the Cahn–Hilliard (CH) equation on an evolving surface. By using a conservation law and transport formulae, we derive the CH equation on evolving surfaces. An evolving surface is discretized using an unstructured triangular mesh. The discrete CH equation is defined on the surface mesh and its dual surface polygonal tessellation. The evolving triangular surfaces are then realized by moving the surface nodes according to a given velocity field. The proposed scheme is based on the Crank–Nicolson scheme and a linearly stabilized splitting scheme. The scheme is second-order accurate, with respect to both space and time. The resulting system of discrete equations is easy to implement, and is solved by using an efficient biconjugate gradient stabilized method. Several numerical experiments are presented to demonstrate the performance and effectiveness of the proposed numerical scheme.
KW - Cahn–Hilliard equation
KW - Evolving surface
KW - Laplace–Beltrami operator
KW - Triangular surface mesh
UR - http://www.scopus.com/inward/record.url?scp=85047663085&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85047663085&partnerID=8YFLogxK
U2 - 10.1007/s10915-018-0742-6
DO - 10.1007/s10915-018-0742-6
M3 - Article
AN - SCOPUS:85047663085
VL - 77
SP - 1147
EP - 1163
JO - Journal of Scientific Computing
JF - Journal of Scientific Computing
SN - 0885-7474
IS - 2
ER -