A least-squares/finite element method for the numerical solution of the Navier-Stokes-Cahn-Hilliard system modeling the motion of the contact line

  • Authors:
  • Qiaolin He;Roland Glowinski;Xiao-Ping Wang

  • Affiliations:
  • Department of Mathematics, Sichuan University, Chengdu 610064, PR China and Department of Mathematics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong;Department of Mathematics, University of Houston, Houston, TX 77204, USA and Institute of Advanced Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong;Department of Mathematics and KAUST-HKUST Micro/Nanofludic Joint Laboratory, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2011

Quantified Score

Hi-index 31.46

Visualization

Abstract

In this article we discuss the numerical solution of the Navier-Stokes-Cahn-Hilliard system modeling the motion of the contact line separating two immiscible incompressible viscous fluids near a solid wall. The method we employ combines a finite element space approximation with a time discretization by operator-splitting. To solve the Cahn-Hilliard part of the problem, we use a least-squares/conjugate gradient method. We also show that the scheme has the total energy decaying in time property under certain conditions. Our numerical experiments indicate that the method discussed here is accurate, stable and efficient.