Conservative multigrid methods for Cahn–Hilliard fluids

Abstract

We develop a conservative, second-order accurate fully implicit discretization of the Navier-Stokes (NS) and CahnHilliard (CH) system that has an associated discrete energy functional. This system provides a diffuse-interface description of binary fluid flows with compressible or incompressible flow components [R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci. 454 (1998Sci. 454 ( ) 2617. In this work, we focus on the case of flows containing two immiscible, incompressible and density-matched components. The scheme, however, has a straightforward extension to multi-component systems. To efficiently solve the discrete system at the implicit time-level, we develop a nonlinear multigrid method to solve the CH equation which is then coupled to a projection method that is used to solve the NS equation. We demonstrate convergence of our scheme numerically in both the presence and absence of flow and perform simulations of phase separation via spinodal decomposition. We examine the separate effects of surface tension and external flow on the decomposition. We find surface tension driven flow alone increases coalescence rates through the retraction of interfaces. When there is an applied external shear, the evolution of the flow is nontrivial and the flow morphology repeats itself in time as multiple pinchoff and reconnection events occur. Eventually, the periodic motion ceases and the system relaxes to a global equilibrium. The equilibria we observe appears has a similar structure in all cases although the dynamics of the evolution is quite different. We view the work presented in this paper as preparatory for a detailed investigation of liquid-liquid interfaces with surface tension where the interfaces separate two immiscible fluids [On the pinchoff of liquid-liquid jets with surface tension, in preparation]. To this end, we also include a simulation of the pinchoff of a liquid thread under the Rayleigh instability at finite Reynolds number.

Nonlinear Modeling and Simulation of Tumor Growth

“…A feature of this approach is the incorporation of a novel continuum model of adhesion in this flux. Following the variational approach developed for diffuse interface models of multiphase flows and materials by Lowengrub and coworkers (e.g., [130,124,123,112,113,203,114]) and others (e.g. [84,101,11]) we introduce a continuum model of cell-cell and cellmatrix adhesion energy that can be written as an integral taken over the entire…”

Computational Phase‐Field Modeling

“…This theory has been widely studied from the physical, the mathematical and the computational points of view Jacqmin, 1999;Gal and Grasselli, 2010;Kay et al, 2008;Boyer et al, 2010;Gal and Grasselli, 2011;Lowengrub and Truskinovsky, 1998;Colli et al, 2012;Kim et al, 2004;Liu and Shen, 2003). To derive the theory, let us assume that we have two immiscible components with volume fractions ϕ 1 and ϕ 2 , respectively.…”

“…Two schemes which avoid constructing the secant, while still maintaining stability, are the implicit Taylor method by (Kim et al, 2004), given by (225)- (226) with…”

Three-dimensional coarsening dynamics of a conserved, nematic liquid crystal-isotropic fluid mixture

“…The same model, which stems from that considered in [35], has been used in the two-dimensional study of Mata et al [26] and the current work is a follow-up report on our findings for the corresponding three-dimensional system. Similar phase field models have been used extensively in phase separation [14,25,19,4,1,2,21,22,5,35,37,36,3,38,13,39,12,30,17,33].…”