We propose a novel second order in time numerical scheme for Cahn-Hilliard-NavierStokes phase field model with matched density. The scheme is based on second order convex-splitting for the Cahn-Hilliard equation and pressure-projection for the Navier-Stokes equation. We show that the scheme is mass-conservative, satisfies a modified energy law and is therefore unconditionally stable. Moreover, we prove that the scheme is unconditionally uniquely solvable at each time step by exploring the monotonicity associated with the scheme. Thanks to the weak coupling of the scheme, we design an efficient Picard iteration procedure to further decouple the computation of Cahn-Hilliard equation and Navier-Stokes equation. We implement the scheme by the mixed finite element method. Ample numerical experiments are performed to validate the accuracy and efficiency of the numerical scheme.
In this paper, we develop a series of linear, unconditionally energy stable numerical schemes for solving the classical phase field crystal model. The temporal discretizations are based on the first order Euler method, the second order backward differentiation formulas (BDF2) and the second order Crank-Nicolson method, respectively. The schemes lead to linear elliptic equations to be solved at each time step, and the induced linear systems are symmetric positive definite. We prove that all three schemes are unconditionally energy stable rigorously. Various classical numerical experiments in 2D and 3D are performed to validate the accuracy and efficiency of the proposed schemes.
We study the well-posedness of a coupled Cahn-Hilliard-Stokes-Darcy system which is a diffuse-interface model for essentially immiscible two phase incompressible flows with matched density in a karstic geometry. Existence of finite energy weak solution that is global in time is established in both 2D and 3D. Weak-strong uniqueness property of the weak solutions is provided as well.
Multiphase flow phenomena are ubiquitous. Common examples include coupled atmosphere and ocean system (air and water), oil reservoir (water, oil and gas), cloud and fog (water vapor, water and air). Multiphase flows also play an important role in many engineering and environmental science applications.In some applications such as flows in unconfined karst aquifers, karst oil reservoir, proton membrane exchange fuel cell, multiphase flows in conduits and in porous media must be considered together. Geometric configurations that contain both conduit (or vug) and porous media are termed karstic geometry. Despite the importance of the subject, little work has been done on multi-phase flows in karstic geometry.In this paper we present a family of phase field (diffusive interface) models for two phase flow in karstic geometry. These models together with the associated interface boundary conditions are derived utilizing Onsager's extremum principle. The models derived enjoy physically important energy laws. Uniquely solvable first and second order in time numerical schemes that preserve the associated energy law are presented as well.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.