Volume 13, Issue 2
A Multiple Scalar Auxiliary Variables Approach to the Energy Stable Scheme of the Moving Contact Line Problem

Numer. Math. Theor. Meth. Appl., 13 (2020), pp. 539-568.

Published online: 2020-03

Cited by

Export citation
• Abstract

In this article, we explore the numerical approximation for a phase-field model to the moving contact line (MCL) problem, governed by the Cahn-Hilliard equation with a dynamic contact angle condition. A difficult and challenging issue comes from the high order, nonlinear and coupled time marching system with a small dimensionless parameter describing the interface thickness. The scalar auxiliary variable (SAV) approach has been recently developed for a large class of gradient flows in principle. Due to the nonlinear gradient terms in both the bulk and boundary energies, we introduce multiple scalar auxiliary variables (MSAV) approach rather than single SAV scheme to deal with phase field contact line model. The MSAV scheme numerically gives a better description of the contact line dynamic  in terms of the difference between the modified and original energies. Moreover, the resulting numerical schemes enjoy the advantages of second order accuracy, unconditional energy stability, and only require solving a linear decoupled system with constant coefficients at each time step. Numerical experiments are presented to verify the effectiveness and efficiency of the proposed schemes for a wide range of mobility parameters and phenomenological boundary relaxation parameters.

35K35, 65M12, 65Z05, 65P40

kangfd@sustech.edu.cn (Fengdai Kang)

zhangz@sustech.edu.cn (Zhen Zhang)

• BibTex
• RIS
• TXT
@Article{NMTMA-13-539, author = {Kang , Fengdai and Zhang , Zhen}, title = {A Multiple Scalar Auxiliary Variables Approach to the Energy Stable Scheme of the Moving Contact Line Problem}, journal = {Numerical Mathematics: Theory, Methods and Applications}, year = {2020}, volume = {13}, number = {2}, pages = {539--568}, abstract = {

In this article, we explore the numerical approximation for a phase-field model to the moving contact line (MCL) problem, governed by the Cahn-Hilliard equation with a dynamic contact angle condition. A difficult and challenging issue comes from the high order, nonlinear and coupled time marching system with a small dimensionless parameter describing the interface thickness. The scalar auxiliary variable (SAV) approach has been recently developed for a large class of gradient flows in principle. Due to the nonlinear gradient terms in both the bulk and boundary energies, we introduce multiple scalar auxiliary variables (MSAV) approach rather than single SAV scheme to deal with phase field contact line model. The MSAV scheme numerically gives a better description of the contact line dynamic  in terms of the difference between the modified and original energies. Moreover, the resulting numerical schemes enjoy the advantages of second order accuracy, unconditional energy stability, and only require solving a linear decoupled system with constant coefficients at each time step. Numerical experiments are presented to verify the effectiveness and efficiency of the proposed schemes for a wide range of mobility parameters and phenomenological boundary relaxation parameters.

}, issn = {2079-7338}, doi = {https://doi.org/10.4208/nmtma.OA-2019-0154}, url = {http://global-sci.org/intro/article_detail/nmtma/15491.html} }
TY - JOUR T1 - A Multiple Scalar Auxiliary Variables Approach to the Energy Stable Scheme of the Moving Contact Line Problem AU - Kang , Fengdai AU - Zhang , Zhen JO - Numerical Mathematics: Theory, Methods and Applications VL - 2 SP - 539 EP - 568 PY - 2020 DA - 2020/03 SN - 13 DO - http://doi.org/10.4208/nmtma.OA-2019-0154 UR - https://global-sci.org/intro/article_detail/nmtma/15491.html KW - Moving contact line, unconditionally energy stable, Cahn-Hilliard equation, multiple scalar auxiliary variables. AB -

In this article, we explore the numerical approximation for a phase-field model to the moving contact line (MCL) problem, governed by the Cahn-Hilliard equation with a dynamic contact angle condition. A difficult and challenging issue comes from the high order, nonlinear and coupled time marching system with a small dimensionless parameter describing the interface thickness. The scalar auxiliary variable (SAV) approach has been recently developed for a large class of gradient flows in principle. Due to the nonlinear gradient terms in both the bulk and boundary energies, we introduce multiple scalar auxiliary variables (MSAV) approach rather than single SAV scheme to deal with phase field contact line model. The MSAV scheme numerically gives a better description of the contact line dynamic  in terms of the difference between the modified and original energies. Moreover, the resulting numerical schemes enjoy the advantages of second order accuracy, unconditional energy stability, and only require solving a linear decoupled system with constant coefficients at each time step. Numerical experiments are presented to verify the effectiveness and efficiency of the proposed schemes for a wide range of mobility parameters and phenomenological boundary relaxation parameters.

Fengdai Kang & Zhen Zhang. (2020). A Multiple Scalar Auxiliary Variables Approach to the Energy Stable Scheme of the Moving Contact Line Problem. Numerical Mathematics: Theory, Methods and Applications. 13 (2). 539-568. doi:10.4208/nmtma.OA-2019-0154
Copy to clipboard
The citation has been copied to your clipboard