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) and the discrete chemical potential bounded in ℓ(0 , T; L2) , 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.


Mathematics and Statistics


National Science Foundation, Grant 1312701

Keywords and Phrases

35K35; 35K55; 65M12; 65M60

International Standard Serial Number (ISSN)


Document Type

Article - Journal

Document Version


File Type





© 2023 Springer, All rights reserved.

Publication Date

01 Nov 2017