Next Article in Journal
Additive Manufacturing of High-Entropy Alloys: A Review
Next Article in Special Issue
Optimization and Entropy Production: Application to Carnot-Like Refrigeration Machines
Previous Article in Journal
A Novel Multi-Exposure Image Fusion Method Based on Adaptive Patch Structure
Previous Article in Special Issue
Magnetic Contribution to the Seebeck Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Triphase Separation of a Ternary Symmetric Highly Viscous Mixture

1
Departmento di Ingegneria dell’Energia, dei Sistemi, del Territorio e delle Costruzioni, Università of Pisa, Largo Lazzarino 1, 56122 Pisa, Italy
2
Departmento di Ingegneria Civile e Industriale, Università of Pisa, Largo Lazzarino 1, 56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(12), 936; https://doi.org/10.3390/e20120936
Submission received: 8 November 2018 / Revised: 30 November 2018 / Accepted: 3 December 2018 / Published: 6 December 2018
(This article belongs to the Special Issue Second Law: Survey and Application)

Abstract

:
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.

1. Introduction

A ternary mixture of partially miscible liquids, which is homogeneous at some high temperature in the one-phase region of its phase diagram, may separate into three phases after rapid quenching to an unstable location within the spinodal range of its phase diagram [1,2]. Presented herein is a numerical study of isothermal spinodal decomposition of a very viscous ternary liquid mixture in the limit of zero fluidity coefficient (further addressed below), following a thermal quench from a stable state in the one-phase region to an unstable location within the tie triangle of its phase diagram. Although our attention in this work will be restricted to nonreactive (symmetric) ternary liquid mixtures, it should be acknowledged that the behavior of reactive (not necessarily symmetric) ternary mixtures has received considerable attention over the past two decades [3,4,5,6,7,8,9,10]. In fact, even though spinodal decomposition in a number of reaction-diffusion systems is usually triggered by a thermal or compositional quench, actually the necessary displacement from a stable state in the one-phase region into the spinodal range can also occur by reaction, e.g., polymerization of a monomer when the polymer being formed is incompatible with one or more of the other components in the mixture. Such reaction-induced spinodal decomposition is particularly relevant to the design of new materials, such as, e.g., monolithic porous polymer or carbon materials, and consequently has been studied extensively over the past two decades, mostly by experiment [11,12,13,14,15,16,17,18]. However, in what follows, we will not address reaction-induced spinodal decomposition in any detail, since the primary objective of this work is to elucidate whether or not the assumption of a perfectly symmetric ternary system implies some kind of statistical symmetry during the triphase segregation process after spinodal decomposition triggered by a thermal quench. We will also briefly discuss segregation patterns (i.e., isosurfaces of concentration as a function of time) in two cases wherein phase separation for one component in a ternary mixture turns out to be faster or slower than that for the remaining components. Note that, since we focus on a ternary system with a vanishing ratio of capillary to viscous forces, the results reported herein may be considered as pertaining to ternary alloy systems or polymer blends. In fact, as is well known, the mechanical properties and rheology of ternary polymer blends are greatly influenced by their phase morphology, and, consequently, the most important application of phase separation dynamics to date has been to predict those phase morphologies as a function of controlling factors such as viscosity of components, composition, interfacial interaction between phases, and processing parameters. For the sake of the present discussion, we note that, while in a binary blend there are only two basic morphologies (particulate and bicontinuous), a large number of phase morphologies have been documented in ternary polymer blends, and, consequently, a taxonomy of phase morphologies has been the subject of a number of previous works; in particular, Nauman and He [19] showed 23 morphologies of which 18 were actually observed in their 2D simulation study but also included morphological types which reflected the topological limitations of 2D space. Once again, we emphasize that phase morphologies in a ternary polymer system are beyond the scope of the present work since the focus of this paper is on the triphase segregation kinetics following spinodal decomposition of a very viscous and symmetric ternary liquid mixture. Finally, we also acknowledge that significant progress in understanding the kinetics of ternary phase-separating mixtures in the absence of chemical reactions has been reported in a number of previous works [1,20,21,22,23,24,25].
An outline of the remainder of this paper is as follows. Below, we first present a brief review of the governing equations for our phase-field formulation, emphasizing that the ternary mixture model of interest herein is supposed to describe a nonideal liquid mixture with perfectly symmetric pure species; consequently, our phase-field formulation has been based on the one-parameter Margules correlation since that is the simplest excess Gibbs energy submodel one can write for a nonideal liquid mixture. In contrast, in some cases, more complex excess Gibbs energy submodels are required for modeling nonideal mixtures having a strong asymmetry among the component species (e.g., when two or more components are partially immiscible). That was specifically the case with the benzene-acetonitrile-water mixture considered in our previous work [26], where the phase-field ternary mixture model presented herein was generalized to incorporate the nonrandom, two-liquid (NRTL) equation [27,28,29] as an excess Gibbs energy submodel. Finally, in Section 3, using different initial concentration field perturbations, we discuss simulation results of a symmetric ternary mixture in terms of the integral scale of the pair-correlation function and the separation depth for each component. Conclusions are then presented.

2. Model Description

Below, we provide a brief summary of the governing equations for a diffuse-interface description of partially miscible ternary liquid mixtures, given that a detailed derivation of those equations has been presented elsewhere [2,26,30]. Consider a regular ternary mixture, whose component liquids have the same molar density, ρ . This mixture can be modeled within a diffuse-interface description by assuming that its free energy is the sum of a thermodynamic part and a nonlocal contribution [26,30,31,32,33,34], i.e.,
G = ρ R T V g ˜ d 3 r with g ˜ = g ( x 1 , x 2 ) a 2 2 x 1 · x 2 + x 1 · x 3 + x 2 · x 3 ,
where g is the (dimensionless) thermodynamic (i.e., coarse-grained) bulk free energy density, T the absolute temperature, V the volume, R the universal gas constant, while a is a characteristic length. Note that, since Equation (1) represents a coarse-grained expression of the free energy, the characteristic length is in no way equal to the actual interfacial thickness (still a is representative, on the mesoscale, of the typical interface thickness at local equilibrium). Assuming that the mixture has zero excess volume of mixing and zero excess entropy of mixing, the simplest expression [28,29,35] for the thermodynamic free energy density, g, corresponding to a perfectly symmetric, partially miscible ternary mixture, is the sum of an entropic, ideal part, and an enthalpic (so-called excess) part, with
g = g 0 + x 1 ln x 1 + x 2 ln x 2 + x 3 ln x 3 + Ψ ( x 1 x 2 + x 1 x 3 + x 2 x 3 ) .
Here, g 0 is the free energy of each pure component, while x i ( i = 1 , 2 , 3 ) denotes the molar (and mass) fraction of the ith species, and Ψ is the Margules parameter corresponding to any one component pair [35]. [Obviously, since we are considering an ideally, perfectly symmetric ternary mixture, the Margules parameter for the ( i , j ) component pair must be equal to those for the remaining pairs. Note that, for strongly nonideal systems (particularly for partially immiscible components when there is no symmetry among the pure species), we have shown a derivation of the species balance equations that incorporates NRTL as an excess Gibbs energy submodel [26], in place of the standard one-parameter Margules correlation employed in Equation (2).] When the mixture finds itself in a state of nonequilibrium, it evolves according to the equations governing dynamic processes, expressing standard conservation principles [36]. Herein, we assume isothermal conditions and therefore ignore the energy balance equation; consequently, the governing equations reduce to a pair of (nonreactive) species balance equations, that are normally coupled to the Navier–Stokes equation (in addition to the zero divergence constraint on the mass-averaged velocity that strictly applies to an isopycnic, regular mixture). Nondimensionalizing each species conservation equation (as well as the Navier–Stokes equation) based upon a diffusive scaling [37,38,39,40,41] gives rise to the so-called fluidity coefficient [42]
α R T a 2 M W ν D ,
with M W and ν denoting the molecular mass and kinematic viscosity of each pure component, while D denotes the (same) binary diffusivity for all component pairs. This dimensionless group can be introduced as an inverse capillary number [43], while it can also be interpreted as a Peclet number [33,34,38,39,40,41], i.e., the ratio of convective to diffusive mass fluxes in the species balance equations. In fact, in previous works [33,34,44,45,46,47,48,49], we have noted that, in low-viscosity systems, α is usually of order O ( 10 3 ÷ 10 6 ) , while highly viscous mixtures (e.g., polymer melts and alloys) correspond to a vanishing fluidity coefficient. In the latter case (which is the focus of the work reported herein), the diffuse-interface model describes a diffusive (or antidiffusive) separation process in the absence of flow, and the species balance equations assume the particularly simple form seen earlier [26,30]:
x 1 t = · x 1 x 2 μ ˜ 23 + x 1 ( 1 x 1 ) μ ˜ 13 ,
x 2 t = · x 1 x 2 μ ˜ 13 + x 2 ( 1 x 2 ) μ ˜ 23 .
These equations have been scaled based on the characteristic length a and the diffusive time a 2 / D . The generalized chemical potential differences in Equations (4) and (5) can be readily obtained from the Bakhuis relation [50], μ ˜ i 3 = δ g ˜ δ x i ( i = 1 , 2 ). In particular, we find
μ ˜ 13 = ln x 1 x 3 + Ψ ( x 3 x 1 ) a 2 2 2 x 1 x 3 ,
μ ˜ 23 = ln x 2 x 3 + Ψ ( x 3 x 2 ) a 2 2 2 x 2 x 3 .
These relations, together with the equalities μ ˜ i j = μ ˜ i k + μ ˜ k j and μ ˜ i j = μ ˜ j i , define all chemical potential differences for a ternary system. Finally, it is worth reiterating that Equations (4) and (5) constitute a system of fourth-order equations, which represents a generalization (specifically, a ternary version of Model B in the taxonomy of Hohenberg and Halperin [51]) of the classical Cahn–Hilliard equation to describe phase separation in binary mixtures [33,34,37,38,39,41,52,53,54].

3. Results and Discussion

Numerical methods employed herein for integrating the species balance equations above in a periodic box are exactly the same as discussed in a previous work [26] and therefore will not be repeated below. In fact, due to our assumption of periodic boundaries, initially we attempted simulations based on a pseudospectral spatial discretization [30,38,39,40,41,49], in conjunction with an adaptive temporal scheme (which had been successfully employed in simulations of dissolution or growth of a multicomponent drop into the continuous phase of another liquid [26,30]). Subsequently, even though our pseudospectral code achieved a stable integration, we realized (based on numerical tests) that a pseudospectral discretization does not preserve the scalar boundedness of the concentration fields. Consequently, the temporal evolution of the separation depth for some component usually went out of bounds, since at some locations in the computational domain the mass fraction for that component became slightly larger than one, or less than zero (or both) for a significant fraction of simulation time. However, assuming the same explicit adaptive temporal scheme as in our pseudospectral code, we also found that issues of scalar boundedness could be circumvented by changing the spatial discretization from pseudospectral to finite volume. (Incidentally, we chose a standard finite volume discretization based on a cell-centered variable arrangement, and using central differencing for both interpolation and differentiation; consequently, our scheme is second-order accurate in space since it relies on a uniform Cartesian grid.) As to the time marching, we also noticed that a stable integration could be achieved using the standard fourth-order Runge–Kutta scheme. Finally, we chose this last temporal scheme (in conjunction with a standard finite-volume discretization) in our production runs for calculating the results presented below.
In our simulations we chose Ψ = 4 , since that was the value for the Margules parameter in the triphase separation simulation of an ideally perfectly symmetric ternary mixture reported by Park et al. [2]. In fact, for Ψ = 4 > 8 3 the phase diagram for the mixture is a tie triangle [1], whose vertices represent the compositions of three coexisting phases (see Figure 1) at local equilibrium. If we let x ξ = ( x , x , 1 2 x ) , x η = ( 1 2 x , x , x ) , and x ζ = ( x , 1 2 x , x ) , the coordinates for the vertices are readily found by observing that all (thermodynamic) chemical potential differences reduce to either a trivial identity or the following relation
μ 21 t h = ln x 1 2 x + Ψ ( 1 3 x ) = 0 x * = 0.023 .
Hence we find x ξ = ( 0.023 , 0.023 , 0.954 ) , x η = ( 0.954 , 0.023 , 0.023 ) , and x ζ = ( 0.023 , 0.954 , 0.023 ) .
We investigate (isothermal) triphase separation in an ideally perfectly symmetric and highly viscous (zero fluidity) ternary mixture, which is instantaneously quenched from a stable state having the initial composition x A = 1 3 , 1 3 , 1 3 in the one-phase region to an unstable state (at a smaller temperature) corresponding to point A in the phase diagram in Figure 1. Assuming an instantaneous quench to a uniform temperature, the initial field for each mass fraction is specified as random (delocalized) concentration fluctuations superimposed on a uniform x i , 0 = 1 3 composition. From this simulation (denoted as case I, more precisely defined below), we show isosurfaces of phase η vs. those for phase ξ at one particular instant in time in Figure 2, suggesting that phase separation for component 3 is faster than that for component 1. For a more quantitative characterization of the the phase-separation kinetics in such a highly viscous system, we looked at the temporal evolution of three characteristic length scales of single-phase microdomains, defined as
L i ( t ) = 1 x i , r m s 2 k | x ^ i ( k ) | 2 | k | ( i = 1 , 2 , 3 ) ,
with x ^ i ( k ) denoting the Fourier coefficient for the ith mass fraction at wavevector k . In other words, L i is the same as the the integral scale of the pair-correlation function associated with the ith mass fraction. As can be seen (Figure 3), at first the integral scales for different components are essentially coincident as they undergo a rapid transient with an undershoot behavior that delays phase separation, after which the three phases seem to grow together, with a coarsening rate that could be expressed approximately as L i τ 0.07 (where τ 10 5 t ); note, however, that there does not seem to be enough of a scaling range for this power law to be physically meaningful. Finally, at about τ 4 × 10 4 , a symmetry-breaking event occurs as component 3 continues to separate quickly, while phase separation for the other two components seems to be delayed. From our simulation results, this behavior seems to be robust, in that very similar temporal histories are obtained for the reciprocal first and second moments of the structure factor associated with the ith mass fraction, i.e.,
λ i ( 1 ) = k | x ^ i ( k ) | 2 k | x ^ i ( k ) | 2 | k | , and λ i ( 2 ) = k | x ^ i ( k ) | 2 k | x ^ i ( k ) | 2 | k | 2 1 / 2 .
We also looked at the temporal evolution of the separation depth for each component, defined as
s i = x i ( r ) x i , 0 x i , e q ( r ) x i , 0 ,
where
x i , e q ( r ) = 1 2 x * if x i ( r ) > x i , 0 , x * if x i ( r ) < x i , 0 .
Clearly, this definition is the same as that employed for studying phase separation in a binary system [38,39,41,49,52,54,55,56]; in fact, the reason that the same definition can be brought to bear on a ternary system is that, based on the mixture phase diagram in Figure 1, both phases in any one pair of coexisting phases at equilibrium possess the same x i , e q (with i denoting the component at the opposite vertex in the phase diagram). As can be seen (Figure 4), phase separation for component 3 is faster than that for the remaining two components, thus confirming the (qualitative) conclusions that were reached by observing Figure 2. In Figure 4 we can see that s 3 at any instant in time is larger than s 1 and s 2 , but for a very small subrange of simulated times wherein all three separation depths show a very rapid increase to their steady state values. In fact, in this case (I) the initial fields for the first two components were taken to be numerically identical in order to have s 1 = s 2 at the initial time; consequently, conservation of mass necessarily implies s 3 ( 0 ) > s 1 ( 0 ) [in fact, it is easily seen that s 3 ( 0 ) 2 s 1 ( 0 ) ].
In another simulation (denoted as case II below), the initial conditions were generated such that the random noise for component 2 exactly canceled that for 1 at the initial time; this necessarily implies s 3 ( 0 ) = 0 , along with s 1 ( 0 ) = s 2 ( 0 ) . In this case (see Figure 5), phase separation for component 3 is slower than that for the remaining two components, given that the separation depth for component 3 remains smaller (and in fact close to zero) than that for 1 and 2 for some time after the initial time, after which all separation depths show a very rapid increase to their steady state values. Once again, even though s 1 and s 2 remain identical for all time, the steady state value for 3 turns out to be slightly larger than that for the first two components. Also note that in this simulation the temporal evolution of the integral scales is very different from that corresponding to case I; in fact, even though the integral scales for 1 and 2 are essentially identical for all time (see Figure 6), the integral scale for 3 is systematically smaller than that for 1 (or 2) while showing a similar decrease–increase behavior. Once again, this temporal evolution seems to be robust, in that very similar temporal histories are also displayed by the reciprocal first and second moments of the structure factor.
Summarizing, in simulation case I [such that s 3 ( 0 ) > s 1 ( 0 ) = s 2 ( 0 ) ] based on the temporal evolution of the integral scales, it seemed as though a symmetry-breaking event took place after which the integral scale vs. time dependences showed a different phase-separating behavior for components 1 and 2 as opposed to 3. However, looking at the separation depth in this case reveals that phase separation for component 3 is always faster than that for the first two components, implying that there can be no symmetry among the three components at any instant in time. A similar conclusion is arrived at by looking at simulation case II [such that 0 = s 3 ( 0 ) < s 1 ( 0 ) = s 2 ( 0 ) ], wherein phase separation for component 3 was always slower than that for the remaining components. Furthermore, in this case the integral scale for 3 was always less than that for 1 (or 2), providing even stronger evidence against any presumption of symmetry among the three components during a triphase segregation process.
As a final note, it should be acknowledged that, even though a perfectly symmetric ternary mixture does not exist in practice, our phase-field ternary mixture model has been successfully employed to simulate the so-called continuous-phase/dispersed-circular/dispersed-circular (a.k.a. dual discrete particle) phase morphology for a phase-separating symmetric polymer-polymer-polymer system [57]. Although improved agreement between model predictions and experimental results of phase morphology for a polystyrene/PMMA/polybutadiene blend was obtained in simulations that allowed for different diffusivities and Flory parameters [58], simulation results that instead assumed equality between all diffusivities, Flory parameters, and chain lengths showed that a phase-field ternary mixture model [59] is able to capture the main features of the previously noted phase morphology [58,60].

4. Conclusions

We have presented 3D numerical results of an isothermal phase-field ternary mixture model in application to a symmetric and highly viscous (zero fluidity parameter) three-component liquid mixture as it separates into three phases after a thermal quench from a stable state in the one-phase region. In particular, we have shown the temporal evolution of basic segregation statistics such as the integral scale of the pair-correlation function and the separation depth for each component. In fact, for a ternary system the latter quantity can be introduced using the same definition of separation depth as employed for measuring the phase separation rate in a binary system. In addition, we have shown that the triphase segregation process considered herein is inherently asymmetrical since it can occur in either one of two ways: depending on the details of the noise used for building the initial conditions, phase separation for the third component will be faster or slower than that for the remaining components.

Author Contributions

Roberto Mauri suggested the problem. Andrea Lamorgese wrote computer codes and performed the simulations, and also wrote the paper. Finally, Roberto Mauri provided useful comments on an early draft of the manuscript.

Acknowledgments

This work was supported by the PROT-ONE Project, FAR-FAS 2014, CUP 4421.02102014.07200099, funded by Regione Toscana.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
2DTwo-Dimensional
3DThree-Dimensional
NRTLNon-Random, Two-Liquid

References

  1. Huang, C.; Olvera de la Cruz, M.; Swift, B.W. Phase separation of ternary mixtures: Symmetric polymer blends. Macromolecules 1995, 28, 7996–8005. [Google Scholar] [CrossRef]
  2. Park, J.M.; Mauri, R.; Anderson, P.D. Phase separation of viscous ternary liquid mixtures. Chem. Eng. Sci. 2012, 80, 270–278. [Google Scholar] [CrossRef]
  3. Liu, B.; Tong, C.; Yang, Y. The kinetics and phase patterns in a ternary mixture coupled with chemical reaction of A + B ⇌ C. J. Phys. Chem. B 2001, 105, 10091–10100. [Google Scholar] [CrossRef]
  4. Tong, C.; Zhang, H.; Yang, Y. Phase separation dynamics and reaction kinetics of ternary mixture coupled with interfacial chemical reaction. J. Phys. Chem. B 2002, 106, 7869–7877. [Google Scholar] [CrossRef]
  5. Tong, C.; Yang, Y. Phase-separation dynamics of a ternary mixture coupled with reversible chemical reaction. J. Chem. Phys. 2002, 116, 1519–1529. [Google Scholar] [CrossRef]
  6. Good, K.; Kuksenok, O.; Buxton, G.A.; Ginzburg, V.V.; Balazs, A.C. Effect of hydrodynamic interactions on the evolution of chemically reactive ternary mixtures. J. Chem. Phys. 2004, 121, 6052–6063. [Google Scholar] [CrossRef] [PubMed]
  7. Travasso, R.D.M.; Buxton, G.A.; Kuksenok, O.; Good, K.; Balazs, A.C. Modeling the morphology and mechanical properties of sheared ternary mixtures. J. Chem. Phys. 2005, 122, 194906. [Google Scholar] [CrossRef] [PubMed]
  8. Kuksenok, O.; Travasso, R.D.M.; Balazs, A.C. Dynamics of ternary mixtures with photosensitive chemical reactions: Creating three-dimensionally ordered blends. Phys. Rev. E 2006, 74, 011502. [Google Scholar] [CrossRef]
  9. Shukutani, T.; Myojo, T.; Nakanishi, H.; Norisuye, T.; Tran-Cong-Miyata, Q. Tricontinuous morphology of ternary polymer blends driven by photopolymerization: Reaction and phase separation kinetics. Macromolecules 2014, 47, 4380–4386. [Google Scholar] [CrossRef]
  10. Tran-Cong-Miyata, Q.; Nakanishi, H. Phase separation of polymer mixtures driven by photochemical reactions: Current status and perspectives. Polym. Int. 2017, 66, 213–222. [Google Scholar] [CrossRef]
  11. Zucchi, I.A.; Galante, M.J.; Borrajo, J.; Williams, R.J. A model system for the thermodynamic analysis of reaction-induced phase separation: Solutions of polystyrene in bifunctionalepoxy/amine monomers. Macromol. Chem. Phys. 2004, 205, 676–683. [Google Scholar] [CrossRef]
  12. Chan, P.K.; Rey, A.D. Polymerization-induced phase separation. 1. Droplet size selection mechanism. Macromolecules 1996, 29, 8934–8941. [Google Scholar] [CrossRef]
  13. Chan, P.K.; Rey, A.D. Polymerization-induced phase separation. 2. Morphological analysis. Macromolecules 1997, 30, 2135–2143. [Google Scholar] [CrossRef]
  14. Lee, K.W.D.; Chan, P.K.; Feng, X. A computational study of the polymerization-induced phase separation phenomenon in polymer solutions under a temperature gradient. Macromol. Theory Simul. 2003, 12, 413–424. [Google Scholar] [CrossRef]
  15. Oh, J.; Rey, A. Computational simulation of polymerization-induced phase separation under a temperature gradient. Comput. Theor. Polymer Sci. 2001, 11, 205–217. [Google Scholar] [CrossRef]
  16. Stubbs, J.M.; Durant, Y.G.; Sundberg, D.C. Polymer phase separation in composite latex particles. 1. Considerations for the nucleation and growth mechanism. C.R. Chimie 2003, 6, 1217–1232. [Google Scholar] [CrossRef]
  17. Sutton, D.; Stanford, J.; Ryan, A. Reaction-induced phase-separation in polyoxyethylene/polystyrene blends. II. Structure development. J. Macromol. Sci. B 2004, 43, 233–251. [Google Scholar] [CrossRef]
  18. Girard-Reydet, E.; Sautereau, H.; Pascault, J.P.; Keates, P.; Navard, P.; Thollet, G.; Vigier, G. Reaction-induced phase separation mechanisms in modified thermosets. Polymer 1998, 39, 2269–2279. [Google Scholar] [CrossRef]
  19. Nauman, E.B.; He, D.Q. Morphology predictions for ternary polymer blends undergoing spinodal decomposition. Polymer 1994, 35, 2243–2255. [Google Scholar] [CrossRef]
  20. Huang, C.; Olvera de la Cruz, M. Adsorption of a minority component in polymer blend interfaces. Phys. Rev. E 1996, 53, 812–819. [Google Scholar] [CrossRef]
  21. Tafa, K.; Puri, S.; Kumar, D. Kinetics of phase separation in ternary mixtures. Phys. Rev. E 2001, 64, 056139. [Google Scholar] [CrossRef] [PubMed]
  22. Ma, Y.Q. Domain patterns in ternary mixtures with different interfacial properties. J. Chem. Phys. 2001, 114, 3734–3738. [Google Scholar] [CrossRef]
  23. Lakshmi, K.C.; Kumar, P.B.S. Scale invariance in coarsening of binary and ternary fluids. Phys. Rev. E 2003, 67, 011507. [Google Scholar] [CrossRef] [PubMed]
  24. Puri, S.; Kumar, D. Autocorrelation functions for phase separation in ternary mixtures. Phys. Rev. E 2004, 70, 051501. [Google Scholar] [CrossRef] [PubMed]
  25. Singh, A.; Puri, S. Phase separation in ternary fluid mixtures: a molecular dynamics study. Soft Matter 2015, 11, 2213–2219. [Google Scholar] [CrossRef] [PubMed]
  26. Lamorgese, A.; Mauri, R. Dissolution or growth of a liquid drop via phase-field ternary mixture model based on the non-random, two-liquid equation. Entropy 2018, 20, 125. [Google Scholar] [CrossRef]
  27. Renon, H.; Prausnitz, J.M. Local compositions in thermodynamic excess functions for liquid mixtures. AIChE. J. 1968, 14, 135–144. [Google Scholar] [CrossRef]
  28. Prausnitz, J.M.; Lichtenthaler, R.N.; Gomes de Azevedo, E. Molecular Thermodynamics of Fluid-Phase Equilibria, 3rd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 1999; Chapter 6. [Google Scholar]
  29. O’Connell, J.P.; Haile, J.M. Thermodynamics: Fundamentals for Applications; Cambridge University Press: New York, NY, USA, 2005; Chapter 11. [Google Scholar]
  30. Lamorgese, A.; Mauri, R. Diffusion-driven dissolution or growth of a liquid drop embedded in a continuous phase of another liquid via phase-field ternary mixture model. Langmuir 2017, 33, 13125–13132. [Google Scholar] [CrossRef]
  31. Cahn, J.W.; Hilliard, J.E. Free energy of a nonuniform system. I. Interfacial free energy. J. Chem. Phys. 1958, 28, 258–267. [Google Scholar] [CrossRef]
  32. Cahn, J.W.; Hilliard, J.E. Free energy of a nonuniform system. III. Nucleation in a two-component incompressible fluid. J. Chem. Phys. 1959, 31, 688–699. [Google Scholar] [CrossRef]
  33. Lamorgese, A.; Molin, D.; Mauri, R. Phase-field approach to multiphase flow modeling. Milan J. Math. 2011, 79, 597–642. [Google Scholar] [CrossRef]
  34. Lamorgese, A.; Mauri, R.; Sagis, L.M.C. Modeling soft interface dominated systems: A comparison of phase field and Gibbs dividing surface models. Phys. Rep. 2017, 675, 1–54. [Google Scholar] [CrossRef]
  35. Sandler, S.I. Chemical, Biochemical, and Engineering Thermodynamics, 4th ed.; Wiley: New York, NY, USA, 2006; Chapter 9. [Google Scholar]
  36. Serrin, J. Mathematical principles of classical fluid mechanics. In Encyclopedia of Physics, Vol. VIII/1; Springer: Berlin, Germany, 1959; pp. 125–263. [Google Scholar]
  37. Vladimirova, N.; Malagoli, A.; Mauri, R. Two-dimensional model of phase segregation in liquid binary mixtures with an initial concentration gradient. Chem. Eng. Sci. 2000, 55, 6109–6128. [Google Scholar] [CrossRef]
  38. Lamorgese, A.; Mauri, R. Phase separation of liquid mixtures. In Nonlinear Dynamics and Control in Process Engineering; Continillo, G., Crescitelli, S., Giona, M., Eds.; Springer: Rome, Italy, 2002; pp. 139–152. [Google Scholar]
  39. Lamorgese, A.; Mauri, R. Nucleation and spinodal decomposition of liquid mixtures. Phys. Fluids 2005, 17, 034107. [Google Scholar] [CrossRef]
  40. Lamorgese, A.; Mauri, R. Mixing of macroscopically quiescent liquid mixtures. Phys. Fluids 2006, 18, 044107. [Google Scholar] [CrossRef]
  41. Lamorgese, A.; Mauri, R. Diffuse-Interface modeling of phase segregation in liquid mixtures. Int. J. Multiphase Flow 2008, 34, 987–995. [Google Scholar] [CrossRef]
  42. Tanaka, H.; Araki, T. Spontaneous double phase separation induced by rapid hydrodynamic coarsening in two-dimensional fluid mixtures. Phys. Rev. Lett. 1998, 81, 389–392. [Google Scholar] [CrossRef]
  43. Vladimirova, N.; Malagoli, A.; Mauri, R. Diffusiophoresis of two-dimensional liquid droplets in a phase-separating system. Phys. Rev. E 1999, 60, 2037–2044. [Google Scholar] [CrossRef]
  44. Gupta, R.; Mauri, R.; Shinnar, R. Liquid-liquid extraction using the composition-induced phase separation process. Ind. Eng. Chem. Res. 1996, 35, 2360–2368. [Google Scholar] [CrossRef]
  45. Gupta, R.; Mauri, R.; Shinnar, R. Phase separation of liquid mixtures in the presence of surfactants. Ind. Eng. Chem. Res. 1999, 38, 2418–2424. [Google Scholar] [CrossRef]
  46. Mauri, R.; Califano, F.; Calvi, E.; Gupta, R.; Shinnar, R. Convection-driven phase segregation of deeply quenched liquid mixtures. J. Chem. Phys. 2003, 118, 8841–8846. [Google Scholar] [CrossRef]
  47. Santonicola, G.; Mauri, R.; Shinnar, R. Phase separation of initially inhomogeneous liquid mixtures. Ind. Eng. Chem. Res. 2001, 40, 2004–2010. [Google Scholar] [CrossRef]
  48. Lamorgese, A.; Mauri, R. Phase-field modeling of interfacial dynamics in emulsion flows: Nonequilibrium surface tension. Int. J. Multiphase Flow 2016, 85, 164–172. [Google Scholar] [CrossRef]
  49. Lamorgese, A.; Mauri, R. Phase-field modeling of mixing/demixing of regular binary mixtures with a composition-dependent viscosity. J. Appl. Phys. 2017, 121, 134302. [Google Scholar] [CrossRef]
  50. Prigogine, I.; Defay, R. Chemical Thermodynamics, 4th ed.; Longmans: London, UK, 1967. [Google Scholar]
  51. Hohenberg, P.C.; Halperin, B.I. Theory of dynamic critical phenomena. Rev. Mod. Phys. 1977, 49, 435–479. [Google Scholar] [CrossRef]
  52. Vladimirova, N.; Malagoli, A.; Mauri, R. Two-dimensional model of phase segregation in liquid binary mixtures. Phys. Rev. E 1999, 60, 6968–6977. [Google Scholar] [CrossRef]
  53. Lamorgese, A.; Mauri, R. Liquid mixture convection during phase separation in a temperature gradient. Phys. Fluids 2011, 23, 034102. [Google Scholar] [CrossRef]
  54. Lamorgese, A.; Mauri, R. Spinodal decomposition of chemically reactive binary mixtures. Phys. Rev. E 2016, 94, 022605. [Google Scholar] [CrossRef]
  55. Vladimirova, N.; Malagoli, A.; Mauri, R. Diffusion-driven phase separation of deeply quenched mixtures. Phys. Rev. E 1998, 58, 7691–7699. [Google Scholar] [CrossRef]
  56. Lamorgese, A.G.; Mauri, R. Diffuse-interface modeling of liquid-vapor phase separation in a van der Waals fluid. Phys. Fluids 2009, 21, 044107. [Google Scholar] [CrossRef] [Green Version]
  57. Cavanaugh, T.J.; Buttle, K.; Turner, J.N.; Nauman, E.B. Confirmation of predicted polymer blend morphologies from bulk specimens. Polymer 1998, 39, 3611–3621. [Google Scholar] [CrossRef]
  58. Alfarraj, A.A.; Nauman, E.B. Spinodal decomposition in ternary systems with significantly different component diffusivities. Macromol. Theory Simul. 2007, 16, 627–631. [Google Scholar] [CrossRef]
  59. Nauman, E.B.; He, D.Q. Nonlinear diffusion and phase separation. Chem. Eng. Sci. 2001, 56, 1999–2018. [Google Scholar] [CrossRef]
  60. Alfarraj, A.A. Muticomponent Diffusion and Phase Separation in Asymmetric Ternary Polymer Systems. Ph.D. Thesis, Rensselaer Polytechnic Institute, New York, NY, USA, 2007. [Google Scholar]
Figure 1. Phase diagram for an ideally perfectly symmetric ternary mixture with Ψ = 4 . Point A denotes an (unstable) equilibrium state [with x A = ( 1 3 , 1 3 , 1 3 ) ] corresponding to the initial mixture composition.
Figure 1. Phase diagram for an ideally perfectly symmetric ternary mixture with Ψ = 4 . Point A denotes an (unstable) equilibrium state [with x A = ( 1 3 , 1 3 , 1 3 ) ] corresponding to the initial mixture composition.
Entropy 20 00936 g001
Figure 2. Isosurfaces of phase η ( x 1 = 0.95 , top) vs. those for phase ξ ( x 3 = 0.95 , bottom) at the same (non-dimensional) time τ = 1.43 × 10 3 from 3D simulation (case I) of a ternary mixture with Ψ = 4 on a 128 3 grid.
Figure 2. Isosurfaces of phase η ( x 1 = 0.95 , top) vs. those for phase ξ ( x 3 = 0.95 , bottom) at the same (non-dimensional) time τ = 1.43 × 10 3 from 3D simulation (case I) of a ternary mixture with Ψ = 4 on a 128 3 grid.
Entropy 20 00936 g002
Figure 3. Integral scales vs. time (diffusive time units) from 3D phase separation simulation (case I) of a perfectly symmetric ternary mixture with Ψ = 4 and global composition x A = 1 3 , 1 3 , 1 3 on a 128 3 grid.
Figure 3. Integral scales vs. time (diffusive time units) from 3D phase separation simulation (case I) of a perfectly symmetric ternary mixture with Ψ = 4 and global composition x A = 1 3 , 1 3 , 1 3 on a 128 3 grid.
Entropy 20 00936 g003
Figure 4. Separation depths vs. time (diffusive time units) from 3D phase-field simulation (case I) of a perfectly symmetric ternary mixture with Ψ = 4 and global composition x A = 1 3 , 1 3 , 1 3 on a 128 3 grid with s 1 ( 0 ) = s 2 ( 0 ) and s 3 ( 0 ) s 1 ( 0 ) + s 2 ( 0 ) .
Figure 4. Separation depths vs. time (diffusive time units) from 3D phase-field simulation (case I) of a perfectly symmetric ternary mixture with Ψ = 4 and global composition x A = 1 3 , 1 3 , 1 3 on a 128 3 grid with s 1 ( 0 ) = s 2 ( 0 ) and s 3 ( 0 ) s 1 ( 0 ) + s 2 ( 0 ) .
Entropy 20 00936 g004
Figure 5. Separation depths vs. time (diffusive time units) from 3D phase-field simulation (case II) of a perfectly symmetric ternary mixture with Ψ = 4 and global composition x A = 1 3 , 1 3 , 1 3 on a 128 3 grid with s 1 ( 0 ) = s 2 ( 0 ) and s 3 ( 0 ) = 0 .
Figure 5. Separation depths vs. time (diffusive time units) from 3D phase-field simulation (case II) of a perfectly symmetric ternary mixture with Ψ = 4 and global composition x A = 1 3 , 1 3 , 1 3 on a 128 3 grid with s 1 ( 0 ) = s 2 ( 0 ) and s 3 ( 0 ) = 0 .
Entropy 20 00936 g005
Figure 6. Integral scales vs. time (diffusive time units) from 3D phase separation simulation (case II) of a perfectly symmetric ternary mixture with Ψ = 4 and global composition x A = 1 3 , 1 3 , 1 3 on a 128 3 grid with s 1 ( 0 ) = s 2 ( 0 ) and s 3 ( 0 ) = 0 .
Figure 6. Integral scales vs. time (diffusive time units) from 3D phase separation simulation (case II) of a perfectly symmetric ternary mixture with Ψ = 4 and global composition x A = 1 3 , 1 3 , 1 3 on a 128 3 grid with s 1 ( 0 ) = s 2 ( 0 ) and s 3 ( 0 ) = 0 .
Entropy 20 00936 g006

Share and Cite

MDPI and ACS Style

Lamorgese, A.; Mauri, R. Triphase Separation of a Ternary Symmetric Highly Viscous Mixture. Entropy 2018, 20, 936. https://doi.org/10.3390/e20120936

AMA Style

Lamorgese A, Mauri R. Triphase Separation of a Ternary Symmetric Highly Viscous Mixture. Entropy. 2018; 20(12):936. https://doi.org/10.3390/e20120936

Chicago/Turabian Style

Lamorgese, Andrea, and Roberto Mauri. 2018. "Triphase Separation of a Ternary Symmetric Highly Viscous Mixture" Entropy 20, no. 12: 936. https://doi.org/10.3390/e20120936

APA Style

Lamorgese, A., & Mauri, R. (2018). Triphase Separation of a Ternary Symmetric Highly Viscous Mixture. Entropy, 20(12), 936. https://doi.org/10.3390/e20120936

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop