An official website of the United States government
In this paper we construct a novel discretization of the Cahn-Hilliard equation coupled with the Navier-Stokes equations. The Cahn-Hilliard equation models the separation of a binary mixture. We construct a very simple time integration scheme for simulating the Cahn-Hilliard equation, which is based on splitting the fourth-order equation into two second-order Helmholtz equations. We combine the Cahn-Hilliard equation with the Navier-Stokes equations to simulate phase separation in a two-phase fluid flow in two dimensions. The scheme conserves mass and momentum and exhibits consistency between mass and momentum, allowing it to be used with large density ratios. We introduce a novel discretization of the surface tension force from the phase-field variable that has finite support around the transition region. The model has a parameter that allows it to transition from a smoothed continuum surface force to a fully sharp interface formulation. We show that our method achieves second-order accuracy, and we compare our method to previous work in a variety of experiments.
more »
« less
This content will become publicly available on January 25, 2026
Stability of cylindrical, multicomponent vesicles
Vesicles are important surrogate structures made up of multiple phospholipids and cholesterol distributed in the form of a lipid bilayer. Tubular vesicles can undergo pearling – i.e. formation of beads on the liquid thread akin to the Rayleigh–Plateau instability. Previous studies have inspected the effects of surface tension on the pearling instabilities of single-component vesicles. In this study, we perform a linear stability analysis on a multicomponent cylindrical vesicle. We solve the Stokes equations along with the Cahn–Hilliard equation to develop the linearized dynamic equations governing the vesicle shape and surface concentration fields. This helps us to show that multicomponent vesicles can undergo pearling, buckling and wrinkling even in the absence of surface tension, which is a significantly different result from studies on single-component vesicles. This behaviour arises due to the competition between the free energies of phase separation, line tension and bending for this multi-phospholipid system. We determine the conditions under which axisymmetric and non-axisymmetric modes are dominant, and supplement our results with an energy analysis that shows the sources for these instabilities. Lastly, we delve into a weakly nonlinear analysis where we solve the nonlinear Cahn–Hilliard equation in the weak deformation limit to understand how mode-mixing alters the late time dynamics of coarsening. We show that in many situations, the trends from our simulations qualitatively match recent experiments (Yanagisawaet al.,Phys. Rev. E, vol. 82, 2010, p. 051928).
more »
« less
- Award ID(s):
- 2147559
- PAR ID:
- 10569923
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 1003
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The strong functionalized Cahn–Hilliard equation models self assembly of amphiphilic polymers in solvent. It supports codimension one and two structures that each admit two classes of bifurcations: pearling, a short-wavelength in-plane modulation of interfacial width, and meandering, a long-wavelength instability that induces a transition to curve-lengthening flow. These two potential instabilities afford distinctive routes to changes in codimension and creation of non-codimensional defects such as end caps and Y-junctions. Prior work has characterized the onset of pearling, showing that it couples strongly to the spatially constant, temporally dynamic, bulk value of the chemical potential. We present a multiscale analysis of the competitive evolution of codimension one and two structures of amphiphilic polymers within the H−1 gradient flow of the strong Functionalized Cahn–Hilliard equation. Specifically we show that structures of each codimension transition from a curve lengthening to a curve shortening flow as the chemical potential falls through a corresponding critical value. The differences in these critical values quantify the competition between the morphologies of differing codimension for the amphiphilic polymer mass. We present a bifurcation diagram for the morphological competition and compare our results quantitatively to simulations of the full system and qualitatively to simulations of self-consistent mean field models and laboratory experiments. In particular we propose that the experimentally observed onset of morphological complexity arises from a transient passage through pearling instability while the associated flow is in the curve lengthening regime.more » « less
-
In this article, we consider a phase field model with different densities and viscosities for the coupled two-phase porous media flow and two-phase free flow, as well as the corresponding numerical simulation. This model consists of three parts: a Cahn–Hilliard–Darcy system with different densities/viscosities describing the porous media flow in matrix, a Cahn–Hilliard–Navier–Stokes system with different densities/viscosities describing the free fluid in conduit, and seven interface conditions coupling the flows in the matrix and the conduit. Based on the separate Cahn–Hilliard equations in the porous media region and the free flow region, a weak formulation is proposed to incorporate the two-phase systems of the two regions and the seven interface conditions between them, and the corresponding energy law is proved for the model. A fully decoupled numerical scheme, including the novel decoupling of the Cahn–Hilliard equations through the four phase interface conditions, is developed to solve this coupled nonlinear phase field model. An energy-law preservation is analyzed for the temporal semi-discretization scheme. Furthermore, a fully discretized Galerkin finite element method is proposed. Six numerical examples are provided to demonstrate the accuracy, discrete energy law, and applicability of the proposed fully decoupled scheme.more » « less
-
Abstract We propose two mass and heat energy conservative, unconditionally stable, decoupled numerical algorithms for solving the Cahn–Hilliard–Navier–Stokes–Darcy–Boussinesq system that models thermal convection of two‐phase flows in superposed free flow and porous media. The schemes totally decouple the computation of the Cahn–Hilliard equation, the Darcy equations, the heat equation, the Navier–Stokes equations at each time step, and thus significantly reducing the computational cost. We rigorously show that the schemes are conservative and energy‐law preserving. Numerical results are presented to demonstrate the accuracy and stability of the algorithms.more » « less
-
In this paper, we consider the numerical approximation for a phase field model of the coupled two-phase free flow and two-phase porous media flow. This model consists of Cahn– Hilliard–Navier–Stokes equations in the free flow region and Cahn–Hilliard–Darcy equations in the porous media region that are coupled by seven interface conditions. The coupled system is decoupled based on the interface conditions and the solution values on the interface from the previous time step. A fully discretized scheme with finite elements for the spatial discretization is developed to solve the decoupled system. In order to deal with the difficulties arising from the interface conditions, the decoupled scheme needs to be constructed appropriately for the interface terms, and a modified discrete energy is introduced with an interface component. Furthermore, the scheme is linearized and energy stable. Hence, at each time step one need only solve a linear elliptic system for each of the two decoupled equations. Stability of the model and the proposed method is rigorously proved. Numerical experiments are presented to illustrate the features of the proposed numerical method and verify the theoretical conclusions.more » « less
