A column pre-ordering strategy for the unsymmetric-pattern multifrontal method
ACM Transactions on Mathematical Software (TOMS)
Error analysis of a mixed finite element method for the Cahn-Hilliard equation
Numerische Mathematik
SIAM Journal on Numerical Analysis
Short note: Spontaneous shrinkage of drops and mass conservation in phase-field simulations
Journal of Computational Physics
A gradient stable scheme for a phase field model for the moving contact line problem
Journal of Computational Physics
Hi-index | 31.45 |
In this paper, a semi-implicit finite element method is presented for the coupled Cahn-Hilliard and Navier-Stokes equations with the generalized Navier boundary condition for the moving contact line problems. In our method, the system is solved in a decoupled way. For the Cahn-Hilliard equations, a convex splitting scheme is used along with a P1-P1 finite element discretization. The scheme is unconditionally stable. A linearized semi-implicit P2-P0 mixed finite element method is employed to solve the Navier-Stokes equations. With our method, the generalized Navier boundary condition is extended to handle the moving contact line problems with complex boundary in a very natural way. The efficiency and capacity of the present method are well demonstrated with several numerical examples.