Abstract:In this paper, we present a novel second order in time mixed finite element scheme for the Cahn–Hilliard–Navier–Stokes equations with matched densities. The scheme combines a standard second order Crank–Nicolson method for the Navier–Stokes equations and a modification to the Crank–Nicolson method for the Cahn–Hilliard equation. In particular, a second order Adams-Bashforth extrapolation and a trapezoidal rule are included to help preserve the energy stability natural to the Cahn–Hilliard equation. We show that our scheme is unconditionally energy stable with respect to a modification of the continuous free energy of the PDE system. Specifically, the discrete phase variable is shown to be bounded in ℓ∞0,T;L∞\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\ell ^\infty \left( 0,T;L^\infty \right) $$\end{document} and the discrete chemical potential bounded in ℓ∞0,T;L2\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\ell ^\infty \left( 0,T;L^2\right) $$\end{document}, for any time and space step sizes, in two and three dimensions, and for any finite final time T. We subsequently prove that these variables along with the fluid velocity converge with optimal rates in the appropriate energy norms in both two and three dimensions.

A Finite Element Method for the Numerical Solution of the Coupled Cahn-Hilliard and Navier-Stokes System for Moving Contact Line Problems

A quasi-Lagrangian finite element method for the Navier-Stokes equations in a time-dependent domain

3D Adaptive Finite Element Method for a Phase Field Model for the Moving Contact Line Problems

An Energy-stable Finite Element Method for the Simulation of Moving Contact Lines in Two-phase Flows

Finite element methods for a class of continuum models for immiscible flows with moving contact lines

A Numerical Scheme for Coupled CFD and CSD Simulation of 3-D Fluid-Membrane Interactions

A Conservative Finite Volume-Particle Hybrid Method for Simulation of Incompressible Interfacial Flow

Convergence analysis and error estimates for a second order accurate finite element method for the Cahn–Hilliard–Navier–Stokes system

A high-order accurate unconditionally stable bound-preserving numerical scheme for the Cahn-Hilliard-Navier-Stokes equations

The subdivision-based IGA-EIEQ numerical scheme for the Navier–Stokes equations coupled with Cahn–Hilliard phase-field model of two-phase incompressible flow on complex curved surfaces

Numerical Solution of the Navier‐Stokes Equations by the Finite Element Method

Moving Mesh Finite Element Method for Unsteady Navier-Stokes Flow

Allen–Cahn–Navier–Stokes–Voigt Systems with Moving Contact Lines

Well-posedness of Navier-Stokes/Cahn-Hilliard equations modeling the dynamics of contact line in a channel

A simple and efficient finite difference scheme to the Cahn–Hilliard–Navier–Stokes system equations

Finite Element Iterative Methods for the 3D Steady Navier--Stokes Equations

A second order accuracy preserving method for moving contact lines with Stokes flow

Decoupled Modified Characteristic Finite Element Method with Different Subdomain Time Steps for Nonstationary Dual-Porosity-Navier-Stokes Model

Moving Mesh Finite Element Methods for the Incompressible Navier-Stokes Equations.

Convergence analysis of a decoupled pressure-correction SAV-FEM for the Cahn–Hilliard–Navier–Stokes model

A new decoupled unconditionally stable scheme and its optimal error analysis for the Cahn-Hilliard-Navier-Stokes equations