Two-dimensional model of phase segregation in liquid binary mixtures

Phase transitions affected by natural and forceful molecular interconversion

If a binary liquid mixture, composed of two alternative species with equal amounts, is quenched from a high temperature to a low temperature, below the critical point of demixing, then the mixture will phase separate through a process known as spinodal decomposition. However, if the two alternative species are allowed to interconvert, either naturally (e.g., the equilibrium interconversion of enantiomers) or forcefully (e.g., via an external source of energy or matter), then the process of phase separation may drastically change. In this case, depending on the nature of interconversion, two phenomena could be observed: either phase amplification, the growth of one phase at the expense of another stable phase, or microphase separation, the formation of nongrowing (steady-state) microphase domains. In this work, we phenomenologically generalize the Cahn–Hilliard theory of spinodal decomposition to include the molecular interconversion of species and describe the physical properties of systems undergoing either phase amplification or microphase separation. We apply the developed phenomenology to accurately describe the simulation results of three atomistic models that demonstrate phase amplification and/or microphase separation. We also discuss the application of our approach to phase transitions in polyamorphic liquids. Finally, we describe the effects of fluctuations of the order parameter in the critical region on phase amplification and microphase separation.

Triphase Separation of a Ternary Symmetric Highly Viscous Mixture

We discuss numerical results of diffusion-driven separation into three phases of a symmetric, three-component highly viscous liquid mixture after an instantaneous quench from the one-phase region into an unstable location within the tie triangle of its phase diagram. Our theoretical approach follows a diffuse-interface model of partially miscible ternary liquid mixtures that incorporates the one-parameter Margules correlation as a submodel for the enthalpic (so-called excess) component of the Gibbs energy of mixing, while its nonlocal part is represented based on a square-gradient (Cahn–Hilliard-type) modeling assumption. The governing equations for this phase-field ternary mixture model are simulated in 3D, showing the segregation kinetics in terms of basic segregation statistics, such as the integral scale of the pair-correlation function and the separation depth for each component. Based on the temporal evolution of the integral scales, phase separation takes place via the simultaneous growth of three phases up until a symmetry-breaking event after which one component continues to separate quickly, while phase separation for the other two seems to be delayed. However, inspection of the separation depths reveals that there can be no symmetry among the three components at any instant in time during a triphase segregation process.

Advancing micro-scale cooling by utilizing liquid-liquid phase separation

Achieving effective cooling within limited space is one of the key challenges for miniaturized product design. State-of-the-art micro-scale cooling enhancement techniques incorporate flow disturbances and boiling to reach high performance. However, these methods face the inherent issues of extra pressure drop, flow instability and dry-out that limits heat flux. Here we demonstrate that substantial cooling capability enhancement, up to 2.5 times, is realized by introducing the phase separation of a triethylamine (TEA)/water mixture at the micro-scale. Our experiments show that the enhancement behavior is closely related to the system’s initial composition, temperature, and flow conditions. Moreover, the mixture system exhibits reduced pressure drop after separation, which makes it more promising in serving practical applications. The results reveal new possibilities for liquid coolant selection and provide the experimental foundation for further research in this area.

Dissolution or Growth of a Liquid Drop via Phase-Field Ternary Mixture Model Based on the Non-Random, Two-Liquid Equation

We simulate the diffusion-driven dissolution or growth of a single-component liquid drop embedded in a continuous phase of a binary liquid. Our theoretical approach follows a diffuse-interface model of partially miscible ternary liquid mixtures that incorporates the non-random, two-liquid (NRTL) equation as a submodel for the enthalpic (so-called excess) component of the Gibbs energy of mixing, while its nonlocal part is represented based on a square-gradient (Cahn-Hilliard-type modeling) assumption. The governing equations for this phase-field ternary mixture model are simulated in 2D, showing that, for a single-component drop embedded in a continuous phase of a binary liquid (which is highly miscible with either one component of the continuous phase but is essentially immiscible with the other), the size of the drop can either shrink to zero or reach a stationary value, depending on whether the global composition of the mixture is within the one-phase region or the unstable range of the phase diagram.

Diffusion-Driven Dissolution or Growth of a Liquid Drop Embedded in a Continuous Phase of Another Liquid via Phase-Field Ternary Mixture Model

We simulate the diffusion-driven dissolution or growth of a single-component (resp. two-component) drop embedded in a continuous phase of a binary (resp. single-component) liquid. Our theoretical approach follows a standard diffuse-interface model of partially miscible ternary liquid mixtures, which is based on a regular solution model assumption together with a Flory-Huggins and Cahn-Hilliard representation of the excess and nonlocal components of the Gibbs free energy of mixing. Based on 2D simulation results, we show that for a single-component drop embedded in a continuous phase of a binary liquid (which is highly miscible with either one component of the continuous phase but essentially immiscible with the other) the size of the drop can either shrink to zero or reach a stationary value, depending on whether the global composition of the mixture is within the one-phase region or the unstable range of the phase diagram. On the other hand, for an isolated two-component drop embedded in a continuous phase of a single-component liquid (which is essentially immiscible with either one component of the drop but miscible with the other) the size of the drop can either grow or shrink and, in particular, it will eventually go to zero if the global composition of the mixture is within the one-phase region; otherwise, for system locations in the unstable range the size of the drop tends to a constant value as the composition within the drop reaches its final equilibrium value.

Experimental investigations and phase-field simulations of triple-phase-separation kinetics within liquid ternary Co-Cu-Pb immiscible alloys

The phase-separation kinetics and microstructure evolution mechanisms of liquid ternary Co_{43}Cu_{40}Pb_{17} immiscible alloys are investigated by both the drop tube technique and phase-field method. Two successive phase separations take place during droplet falling and lead to the formation of a three-phase three-layer core-shell structure composed of a Co-rich core, a Cu-rich middle layer, and a Pb-rich shell. The Pb-rich shell becomes more and more conspicuous as droplet diameter decreases. Meanwhile, the Co-rich core center gradually moves away from the core-shell center. Theoretical analyses show that a larger temperature gradient inside a smaller alloy droplet induces the accelerated growth of the surface segregation shell during triple-phase separation. The residual Stokes motion and the asymmetric Marangoni convection result in the appearance of an eccentric Co-rich core and the core deviation degree is closely related to the droplet size and initial velocity. A three-dimensional phase-field model of ternary immiscible alloys, which considers the successive phase separations under the combined effects of Marangoni convection and surface segregation, is proposed to explore the formation mechanisms of three-phase core-shell structures. The simulated core-shell morphologies are consistent with the experimental observations, which verifies the model's validity in reproducing the core-shell dynamic evolution. Numerical results reveal that the development of three-phase three-layer core-shell structures can be attributed to the primary and then secondary phase separations dominated simultaneously by Marangoni convection and surface segregation. Furthermore, the effects of droplet temperature gradient on the growth kinetics of the surface segregation shell are analyzed in the light of phase-field theory.

Spinodal decomposition of chemically reactive binary mixtures

We simulate the influence of a reversible isomerization reaction on the phase segregation process occurring after spinodal decomposition of a deeply quenched regular binary mixture, restricting attention to systems wherein material transport occurs solely by diffusion. Our theoretical approach follows a diffuse-interface model of partially miscible binary mixtures wherein the coupling between reaction and diffusion is addressed within the frame of nonequilibrium thermodynamics, leading to a linear dependence of the reaction rate on the chemical affinity. Ultimately, the rate for an elementary reaction depends on the local part of the chemical potential difference since reaction is an inherently local phenomenon. Based on two-dimensional simulation results, we express the competition between segregation and reaction as a function of the Damköhler number. For a phase-separating mixture with components having different physical properties, a skewed phase diagram leads, at large times, to a system converging to a single-phase equilibrium state, corresponding to the absolute minimum of the Gibbs free energy. This conclusion continues to hold for the critical phase separation of an ideally perfectly symmetric binary mixture, where the choice of final equilibrium state at large times depends on the initial mean concentration being slightly larger or less than the critical concentration.

Dynamic evolution process of multilayer core-shell microstructures within containerlessly solidifying Fe(50)Sn(50) immiscible alloy

Multilayer core-shell structures are frequently formed in polymers and alloys when temperature and concentration fields are well symmetrical spatially. Here we report that two- to five-layer core-shell microstructures were the dominant structural morphology of a binary Fe(50)Sn(50) immiscible alloy solidified under the containerless and microgravity states within a drop tube. Three dimensional phase field simulation reveals that both the uniformly dispersive structure and the multilayer core-shells are the various metastable and transitional states of the liquid phase separation process. Only the two-layer core-shell is the most stable microstructure with the lowest chemical potential. Because of the suppression of Stokes motion, solutal Marangoni migration becomes important to drive the evolution of core-shell structures.

Droplet motion in one-component fluids on solid substrates with wettability gradients

Droplet motion on solid substrates has been widely studied not only because of its importance in fundamental research but also because of its promising potentials in droplet-based devices developed for various applications in chemistry, biology, and industry. In this paper, we investigate the motion of an evaporating droplet in one-component fluids on a solid substrate with a wettability gradient. As is well known, there are two major difficulties in the continuum description of fluid flows and heat fluxes near the contact line of droplets on solid substrates, namely, the hydrodynamic (stress) singularity and thermal singularity. To model the droplet motion, we use the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)] for the hydrodynamic equations in the bulk region, supplemented with the boundary conditions at the fluid-solid interface. In this continuum hydrodynamic model, various physical processes involved in the droplet motion can be taken into account simultaneously, e.g., phase transitions (evaporation or condensation), capillary flows, fluid velocity slip, and substrate cooling or heating. Due to the use of the phase field method (diffuse interface method), the hydrodynamic and thermal singularities are resolved automatically. Furthermore, in the dynamic van der Waals theory, the evaporation or condensation rate at the liquid-gas interface is an outcome of the calculation rather than a prerequisite as in most of the other models proposed for evaporating droplets. Numerical results show that the droplet migrates in the direction of increasing wettability on the solid substrates. The migration velocity of the droplet is found to be proportional to the wettability gradients as predicted by Brochard [Langmuir 5, 432 (1989)]. The proportionality coefficient is found to be linearly dependent on the ratio of slip length to initial droplet radius. These results indicate that the steady migration of the droplets results from the balance between the (conservative) driving force due to the wettability gradient and the (dissipative) viscous drag force. In addition, we study the motion of droplets on cooled or heated solid substrates with wettability gradients. The fast temperature variations from the solid to the fluid can be accurately described in the present approach. It is observed that accompanying the droplet migration, the contact lines move through phase transition and boundary velocity slip with their relative contributions mostly determined by the slip length. The results presented in this paper may lead to a more complete understanding of the droplet motion driven by wettability gradients with a detailed picture of the fluid flows and phase transitions in the vicinity of the moving contact line.

Dynamic scaling in phase separation kinetics for quasi-two-dimensional membranes

We consider the dynamics of phase separation in lipid bilayer membranes, modeled as flat two-dimensional liquid sheets within a bulk fluid, both in the creeping flow approximation. We present scaling arguments that suggest asymptotic coarsening in these systems is characterized by a length scale R(t) ~ t(1/2) for critical (bicontinuous) phase separation and R(t) ~t(1/3) for off-critical concentrations (droplet morphology). In this limit, the bulk fluid is the primary source of dissipation. We also address these questions with continuum stochastic hydrodynamic simulations. We see evidence of scaling violation in critical phase separation, where isolated circular domains coarsen slower than elongated ones. However, we also find a region of apparent scaling where R(t) ~ t(1/2) is observed. This appears to be due to the competition of thermal and hydrodynamic effects. We argue that the diversity of scaling exponents measured in experiment and prior simulations can in part be attributed to certain measurements lying outside the asymptotic long-length-scale regime, and provide a framework to help understand these results. We also discuss a few simple generalizations to confined membranes and membranes in which inertia is relevant.

Hydrodynamic effects on spinodal decomposition kinetics in planar lipid bilayer membranes

The formation and dynamics of spatially extended compositional domains in multicomponent lipid membranes lie at the heart of many important biological and biophysical phenomena. While the thermodynamic basis for domain formation has been explored extensively in the past, domain growth in the presence of hydrodynamic interactions both within the (effectively) two-dimensional membrane and in the three-dimensional solvent in which the membrane is immersed has received little attention. In this work, we explore the role of hydrodynamic effects on spinodal decomposition kinetics via continuum simulations of a convective Cahn-Hilliard equation for membrane composition coupled to the Stokes equation. Our approach explicitly includes hydrodynamics both within the planar membrane and in the three-dimensional solvent in the viscously dominated flow regime. Numerical simulations reveal that dynamical scaling breaks down for critical lipid mixtures due to distinct coarsening mechanisms for elongated versus more isotropic compositional lipid domains. The breakdown in scaling should be readily observable in experiments on model membrane systems.

Evidence of convective heat transfer enhancement induced by spinodal decomposition

Spinodal decomposition can be driven by either diffusion or self-induced convection; the importance of convection relative to diffusion depends on the Péclet number, defined as the ratio between convective and diffusive mass fluxes. Diffusion is the dominating mechanism of phase segregation when the Péclet number is small - i.e., when viscosity and diffusivity are large - or when the domain characteristic size is small. For low-viscosity mixtures, convection is the dominating process and the segregation is very rapid as it takes a few seconds compared to the hours needed in the case of pure diffusion. In such cases, strong convective motion of the phase segregating domains is generated even in small-size systems and is almost independent of the temperature difference as long as it is below the transition value. We study experimentally the enhancement of heat transfer in a 1-mm -thick cell. A water-acetonitrile-toulene mixture is quenched into a two-phase region so as to induce convection-driven spinodal decomposition. The heat transfer rate is measured and compared to that obtained in the absence of convective motion. A substantial reduction in the cooling time obtains in the case of spinodal decomposition. The heat transfer enhancement induced by this self-induced, disordered but effectively convective effect may be exploited in the cooling or heating of small-scale systems whereby forced convection cannot be achieved because of the small sizes involved. A scaling analysis of the data based on the diffuse interface H model for a symmetric mixture near the equilibrium point yields very encouraging agreement and insights.

Separation of liquid phases in giant vesicles of ternary mixtures of phospholipids and cholesterol

We use fluorescence microscopy to directly observe liquid phases in giant unilamellar vesicles. We find that a long list of ternary mixtures of high melting temperature (saturated) lipids, low melting temperature (usually unsaturated) lipids, and cholesterol produce liquid domains. For one model mixture in particular, DPPC/DOPC/Chol, we have mapped phase boundaries for the full ternary system. For this mixture we observe two coexisting liquid phases over a wide range of lipid composition and temperature, with one phase rich in the unsaturated lipid and the other rich in the saturated lipid and cholesterol. We find a simple relationship between chain melting temperature and miscibility transition temperature that holds for both phosphatidylcholine and sphingomyelin lipids. We experimentally cross miscibility boundaries both by changing temperature and by the depletion of cholesterol with beta-cyclodextrin. Liquid domains in vesicles exhibit interesting behavior: they collide and coalesce, can finger into stripes, and can bulge out of the vesicle. To date, we have not observed macroscopic separation of liquid phases in only binary lipid mixtures.