TY - JOUR
T1 - A numerical method for the quasi-incompressible Cahn-Hilliard-Navier-Stokes equations for variable density flows with a discrete energy law
AU - Guo, Z.
AU - Lin, P.
AU - Lowengrub, J. S.
PY - 2014/11/1
Y1 - 2014/11/1
N2 - In this paper, we investigate numerically a diffuse interface model for the Navier-Stokes equation with fluid-fluid interface when the fluids have different densities [48]. Under minor reformulation of the system, we show that there is a continuous energy law underlying the system, assuming that all variables have reasonable regularities. It is shown in the literature that an energy law preserving method will perform better for multiphase problems. Thus for the reformulated system, we design a finite element method and a special temporal scheme where the energy law is preserved at the discrete level. Such a discrete energy law (almost the same as the continuous energy law) for this variable density two-phase flow model has never been established before with finite element. A Newton method is introduced to linearise the highly non-linear system of our discretization scheme. Some numerical experiments are carried out using the adaptive mesh to investigate the scenario of coalescing and rising drops with differing density ratio. The snapshots for the evolution of the interface together with the adaptive mesh at different times are presented to show that the evolution, including the break-up/pinch-off of the drop, can be handled smoothly by our numerical scheme. The discrete energy functional for the system is examined to show that the energy law at the discrete level is preserved by our scheme.
AB - In this paper, we investigate numerically a diffuse interface model for the Navier-Stokes equation with fluid-fluid interface when the fluids have different densities [48]. Under minor reformulation of the system, we show that there is a continuous energy law underlying the system, assuming that all variables have reasonable regularities. It is shown in the literature that an energy law preserving method will perform better for multiphase problems. Thus for the reformulated system, we design a finite element method and a special temporal scheme where the energy law is preserved at the discrete level. Such a discrete energy law (almost the same as the continuous energy law) for this variable density two-phase flow model has never been established before with finite element. A Newton method is introduced to linearise the highly non-linear system of our discretization scheme. Some numerical experiments are carried out using the adaptive mesh to investigate the scenario of coalescing and rising drops with differing density ratio. The snapshots for the evolution of the interface together with the adaptive mesh at different times are presented to show that the evolution, including the break-up/pinch-off of the drop, can be handled smoothly by our numerical scheme. The discrete energy functional for the system is examined to show that the energy law at the discrete level is preserved by our scheme.
UR - http://www.scopus.com/inward/record.url?scp=84905986864&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2014.07.038
DO - 10.1016/j.jcp.2014.07.038
M3 - Article
AN - SCOPUS:84905986864
SN - 0021-9991
VL - 276
SP - 486
EP - 507
JO - Journal of Computational Physics
JF - Journal of Computational Physics
ER -