Convergence analysis of a numerical scheme for the Cahn-Hilliard-Navier-Stokes system with dynamical boundary condition and its application to moving contact line problem
Abstract
A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Navier-Stokes system, combined with a dynamical boundary condition. Such a physical system has potential applications in the moving contact line problem. The boundary profile is governed by a lower-dimensional energy potential, coupled with a non-homogeneous boundary condition for the phase variable. In the numerical design, a convex-splitting approach is applied to the chemical potential in both bulk and surface levels, which leads to a highly coupled nonlinear system. A semi-implicit discretization is taken in the nonlinear fluid convection, as well as the coupled terms between the fluid motion and phase variable evolution. A careful finite difference approximation and convexity analysis reveals that such a numerical system could be represented as a non-symmetric and monotone mapping associated with the fluid convection. In turn, the unique solvability is valid based on the monotonicity argument. The total energy stability analysis is obtained through a careful summation-by-part calculation. In particular, an optimal rate convergence analysis is theoretically established in this work. The discrete mass conservation of the exact solution is required to preserve the mean-zero property of the error function, so that the associated discrete H_h^{-1} norm is well-defined. A combination of the Fourier projection and an auxiliary function is applied to overcome this difficulty. Furthermore, an approach of rough and refined error estimates concludes the desired convergence result. Some numerical results are presented in this article, which demonstrate the robustness of the proposed numerical scheme. In our knowledge, this work provides a theoretical proof of convergence analysis and error estimate for a numerical scheme to the moving contact line problem, for the first time in the literature.
Related
- Numerical approximation of nonlinear fractional dynamical systems for coastal wave modeling: rigorous analysis of convergence and computational efficiency
- Accelerated Convergence of a Second-Order Dynamical System and its Application to Splitting Algorithms for Comonotone Inclusions
- CONVERGENCE ANALYSIS OF A DISCRETE DUALITY FINITE-VOLUME METHOD FOR CAHN-HILLIARD EQUATIONS WITH SURFACTANTS
- Convergence rate of stochastic gradient descent under Łojasiewicz type condition and polynomially controlled noise
- Mathematical and Numerical Analysis of Tsunami Generation and Propagation: The Initial Tsunami Condition
BibTeX
@misc{guo2026convergence,
title = {{Convergence analysis of a numerical scheme for the Cahn-Hilliard-Navier-Stokes system with dynamical boundary condition and its application to moving contact line problem}},
author = {Yunzhuo Guo and Cheng Wang and Xianmin Xu and Zhengru Zhang},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.28954},
eprint = {2609.28954},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.28954},
}
RIS
TY - UNPB
TI - Convergence analysis of a numerical scheme for the Cahn-Hilliard-Navier-Stokes system with dynamical boundary condition and its application to moving contact line problem
AU - Yunzhuo Guo
AU - Cheng Wang
AU - Xianmin Xu
AU - Zhengru Zhang
PY - 2026
DA - 2026/09/24
PB - arXiv
DO - 10.48550/arxiv.2609.28954
UR - https://arxiv.org/abs/2609.28954
AB - A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Navier-Stokes system, combined with a dynamical boundary condition. Such a physical system has potential applications in the moving contact line problem. The boundary profile is governed by a lower-dimensional energy potential, coupled with a non-homogeneous boundary condition for the phase variable. In the numerical design, a convex-splitting approach is applied to the chemical potential in both bulk and surface levels, which leads to a highly coupled nonlinear system. A semi-implicit discretization is taken in the nonlinear fluid convection, as well as the coupled terms between the fluid motion and phase variable evolution. A careful finite difference approximation and convexity analysis reveals that such a numerical system could be represented as a non-symmetric and monotone mapping associated with the fluid convection. In turn, the unique solvability is valid based on the monotonicity argument. The total energy stability analysis is obtained through a careful summation-by-part calculation. In particular, an optimal rate convergence analysis is theoretically established in this work. The discrete mass conservation of the exact solution is required to preserve the mean-zero property of the error function, so that the associated discrete H_h^{-1} norm is well-defined. A combination of the Fourier projection and an auxiliary function is applied to overcome this difficulty. Furthermore, an approach of rough and refined error estimates concludes the desired convergence result. Some numerical results are presented in this article, which demonstrate the robustness of the proposed numerical scheme. In our knowledge, this work provides a theoretical proof of convergence analysis and error estimate for a numerical scheme to the moving contact line problem, for the first time in the literature.
ER -