Research Papers

A Class of Conservative Phase Field Models for Multiphase Fluid Flows

[+] Author and Article Information
Jun Li

School of Mathematical Sciences and LPMC,
Nankai University,
Tianjin 300071, China
e-mail: nkjunli@gmail.com

Qi Wang

Department of Mathematics,
Interdisciplinary Mathematics Institute and
NanoCenter at USC,
University of South Carolina,
Columbia, SC 29208
School of Mathematics,
Nankai University,
Tianjin 300071, China
e-mail: qwang@math.sc.edu

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received August 12, 2012; final manuscript received March 28, 2013; accepted manuscript posted May 7, 2013; published online September 16, 2013. Assoc. Editor: Nesreen Ghaddar.

J. Appl. Mech 81(2), 021004 (Sep 16, 2013) (8 pages) Paper No: JAM-12-1385; doi: 10.1115/1.4024404 History: Received August 12, 2012; Revised March 28, 2013; Accepted May 07, 2013

Hydrodynamic phase field models for multiphase fluids formulated using volume fractions of incompressible fluid components do not normally conserve mass. In this paper, we formulate phase field theories for mixtures of multiple incompressible fluids, using volume fractions, to ensure conservation of mass and momentum for the fluid mixture as well as the total volume for each fluid phase. In this formulation, the mass-average velocity is nonsolenoidal when the densities of incompressible fluid components in the mixture are not equal, making it a bona fide compressible model subject to an internal constraint. Derivation of mass conservation and energy dissipation in phase field models based on both Allen–Cahn dynamics and Cahn–Hilliard dynamics are presented. One salient feature of the phase field models is that the hydrostatic pressure is coupled with the transport of the volume fractions making the momentum transport and the volume fraction transport fully coupled in light of the mass conservation. Near equilibrium dynamics are explored using a linear analysis. In the case of binary fluid mixtures, one potential growth mode is identified in all the models for a class of free energy, which has been adopted for multiphase fluids. The growth is either absent for all waves or of a longwave feature.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Du, Q., Liu, C., Ryham, R., and Wang, X., 2005, “Phase Field Modeling of the Spontaneous Curvature Effect in Cell Membranes,” Comm. Pure Appl. Anal., 4, pp. 537–548. [CrossRef]
Du, Q., Liu, C., Ryham, R., and Wang, X., 2005, “A Phase Field Formulation of the Willmore Problem,” Nonlinearity, 18, pp. 1249–1267. [CrossRef]
Du, Q., Liu, C., and Wang, X., 2004, “A Phase Field Approach in the Numerical Study of the Elastic Bending Energy for Vesicle Membranes,” J. Comput. Phys., 198, pp. 450–468. [CrossRef]
Du, Q., Liu, C., and Wang, X., 2005, “Retrieving Topological Information for Phase Field Models,” SIAM J. Appl. Math., 65, pp. 1913–1932. [CrossRef]
Du, Q., Liu, C., and Wang, X., 2005, “Simulating the Deformation of Vesicle Membranes Under Elastic Bending Energy in Three Dimensions,” J. Comput. Phys., 212, pp. 757–777. [CrossRef]
Shen, J., and Yang, X., 2009, “An Efficient Moving Mesh Spectral Method for the Phase-Field Model of Two Phase Flows,” J. Comput. Phys., 228, pp. 2978–2992. [CrossRef]
Shen, J., and Yang, X., 2010, “Energy Stable Schemes for Cahn-Hilliard Phase-Field Model of Two-Phase Incompressible Flows,” Chin. Ann. Math., Ser. B, 31, pp. 743–758. [CrossRef]
Shen, J., and Yang, X., 2010, “A Phase-Field Model and Its Numerical Approximation for Two-Phase Incompressible Flows With Different Densities and Viscositites,” SIAM J. Sci. Comput. (USA), 32, pp. 1159–1179. [CrossRef]
Feng, J. J., Liu, C., Shen, J., and Yue, P., 2005, “Transient Drop Deformation Upon Startup of Shear in Viscoelastic Fluids,” Phys. Fluids, 17, p. 123101. [CrossRef]
Hua, J., Lin, P., Liu, C., and Wang, Q., 2011, “Energy Law Preserving C0 Finite Element Schemes for Phase Field Models in Two-Phase Flow Computations,” J. Comput. Phys., 230, pp. 7115–7131. [CrossRef]
Liu, C., and Walkington, N. J., 2001, “An Eulerian Description of Fluids Containing Visco-Hyperelastic Particles,” Arch. Ration. Mech. Anal., 159, pp. 229–252. [CrossRef]
Lowengrub, J., and Truskinovsky, L., 1998, “Quasi-Incompressible Cahn–Hilliard Fluids and Topological Transitions,” Proc. R. Soc. London, Ser. A, 454, pp. 2617–2654. [CrossRef]
Liu, C., and Shen, J., 2003, “A Phase Field Model for the Mixture of Two Incompressible Fluids and Its Approximation by a Fourier-Spectral Method,” Physica D, 179, pp. 211–228. [CrossRef]
Wise, S. M., Lowengrub, J. S., Kim, J. S., and Johnson, W. C., 2004, “Efficient Phase-Field Simulation of Quantum Dot Formation in a Strained Heteroepitaxial Film,” Superlattices Microstruct., 36, pp. 293–304. [CrossRef]
Wheeler, A., McFadden, G., and Boettinger, W., 1996, “Phase-Field Model for Solidification of a Eutectic Alloy,” Proc. R. Soc. London, Ser. A, 452, pp. 495–525. [CrossRef]
Yang, X., Feng, J., Liu, C., and Shen, J., 2006, “Numerical Simulations of Jet Pinching-Off and Drop Formation Using an Energetic Variational Phase-Field Method,” J. Comput. Phys., 218, pp. 417–428. [CrossRef]
Yue, P., Feng, J. J., Liu, C., and Shen, J., 2004, “A Diffuse-Interface Method for Simulating Two-Phase Flows of Complex Fluids,” ASME J. Fluid Mech., 515, pp. 293–317. [CrossRef]
Yue, P., Feng, J. J., Liu, C., and Shen, J., 2005, “Diffuse-Interface Simulations of Drop Coalescence and Retraction in Viscoelastic Fluids,” J. Non-Newtonian Fluid Mech., 129, pp. 163–176. [CrossRef]
Zhang, T. Y., Cogan, N., and Wang, Q., 2008, “Phase Field Models for Biofilms. II. 2-D Numerical Simulations of Biofilm-Flow Interaction,” Comm. Comp. Phys., 4(1), pp. 72–101.
Kim, J., and Lowengrub, J., 2005, “Phase Field Modeling and Simulation of Three-Phase Flows,” Interface Free Boundaries, 7, pp. 435–466. [CrossRef]
Hobayashi, R., 1993, “Modeling and Numerical Simulations of Dendritic Crystal Growth,” Physica D, 63, pp. 410–423. [CrossRef]
McFadden, G., Wheeler, A., Braun, R., Coriell, S., and Sekerka, R., 1998, “Phase-Field Models for Anisotropic Interfaces,” Phys. Rev. E, 48, pp. 2016–2024. [CrossRef]
Karma, A., and Rappel, W., 1999, “Phase-Field Model of Dendritic Sidebranching With Thermal Noise,” Phys. Rev. E, 60, pp. 3614–3625. [CrossRef]
Chen, L. Q., and Yang, W., 1994, “Computer Simulation of the Dynamics of a Quenched System With Large Number of Non-Conserved Order Parameters,” Phys. Rev. B, 50, pp. 15752–15756. [CrossRef]
Chen, L. Q., 2002, “Phase-Field Modeling for Microstructure Evolution,” Annu. Rev. Mater. Res., 32, pp. 113–140. [CrossRef]
Li, Y., Hu, S., Liu, Z., and Chen, L., 2001, “Phase-Field Model of Domain Structures in Ferroelectric Thin Films,” Appl. Phys. Lett., 78, pp. 3878–3880. [CrossRef]
Seol, D. J., Hu, S. Y., Li, Y. L., Shen, J., Oh, K. H., and Chen, L. Q., 2003, “Three-Dimensional Phase-Field Modeling of Spinodal Decomposition in Constrained Films,” Acta Mater., 51, pp. 5173–5185. [CrossRef]
Lu, W., and Suo, Z., 2001, “Dynamics of Nanoscale Pattern Formation of an Epitaxial Monolayer,” J. Mech. Phys. Solids, 49, pp. 1937–1950. [CrossRef]
Chen, L. Q., and Wang, Y., 1996, “The Continuum Field Approach to Modeling Microstructural Evolution,” JOM, 48, pp. 13–18. [CrossRef]
Wang, Y., and Chen, C. L., 1999, “Simulation of Microstructure Evolution,” Methods in Materials Research, John Wiley and Sons, New York, pp. 2a3.1–2a3.23.
Tadmor, E., Phillips, R., and Ortiz, M., 1996, “Mixed Atomistic and Continuum Models of Deformation in Solids,” Langmuir, 12, pp. 4529–4534. [CrossRef]
Forest, M. G., Liao, Q., and Wang, Q., 2010, “2-D Kinetic Theory for Polymer Particulate Nanocomposites,” Comm. Comp. Phys., 7, pp. 250–282. [CrossRef]
Zhang, T. Y., and Wang, Q., 2010, “Cahn–Hilliard vs Singular Cahn–Hilliard Equations in Phase Field Modeling,” Comm. Comp. Phys., 7, pp. 362–382. [CrossRef]
Lindley, B., Wang, Q., and Zhang, T., 2011, “Multicomponent Models for Biofilm Flows,” Discrete Contin. Dyn. Syst., Ser. B, 15, pp. 417–456. [CrossRef]
Cahn, J. W., and Hilliard, J. E., 1959, “Free Energy of a Nonuniform System. I. Interfacial Free Energy,” J. Chem. Phys., 28, pp. 258–267. [CrossRef]
Cahn, J. W., and Hilliard, J. E., 1959, “Free Energy of a Nonuniform System-III: Nucleation in a 2-Component Incompressible Fluid,” J. Chem. Phys., 31, pp. 688–699. [CrossRef]
Flory, P. J., 1953, Principles of Polymer Chemistry, Cornell University, Ithaca, NY.
Doi, M., 1996, Introduction to Polymer Physics, Clarendon Press, Oxford, UK.
Coates, D., 1995, “Polymer-Dispersed Liquid Crystals,” J. Mater. Chem., 5, pp. 2063–2072. [CrossRef]
Larson, R. G., 1999, The Structure and Rheology of Complex Fluids, Oxford University, New York.
Shen, J., Yang, X., and Wang, Q., 2013, “Mass Conserved Phase Field Model for Binary Fluids,” Comm. Comp. Phys., 13, pp. 1045–1065. [CrossRef]
Bird, R. B., Stewart, W. E., and Lightfoot, E. N., 2002, Transport Phenomena, Wiley, New York.
Beris, A. N., and Edwards, B., 1994, Thermodynamics of Flowing Systems, Oxford University, Oxford, UK.
Probstein, R. F., 1994, Physicochemical Hydrodynamics, Wiley, New York.
Zhang, T. Y., Cogan, N., and Wang, Q., 2008, “Phase Field Models for Biofilms. I. Theory and 1-D Simulations,” SIAM J. Appl. Math., 69(3), pp. 641–669. [CrossRef]


Grahic Jump Location
Fig. 1

Growth rates for phase field models for water and glycerin mixtures. The growth rate in the Cahn–Hilliard I model is higher than that in the Cahn–Hilliard II just like the situation for the Allen–Cahn models. The parameter values used in the figure: water density ρ1 = 1000 kg/m3, glycerin density ρ2 = 1261 kg/m3, water volume fraction ϕ0 = 0.7, ρ0 = ϕ0ρ1 + ρ2 (1 − ϕ0), η = 0.4503 Pa · s, ν = 0.4517 Pa · s, Γ1 = 10−5 kg/s2, Γ2 = 100 kg/s2m, water viscosity is 8.9 × 10−4 Pa · s, and glycerin viscosity is 1.499 Pa · s. (a) Two Allen–Cahn models, where λ1 = 10−4 kg/sm, λ2 = 10−3 kg/sm. The fastest growth rate is 446.2457104 at k = 1059.064423 in AC-I model and 446.2398776 at k = 1055.023079 in AC-II model. (b). Two Allen–Cahn models, where λ1 = 104 kg/sm, λ2 = 103 kg/sm. The fastest growth rate is 0.003655986108 at k = 3.031398953 in the AC-I model and 0.003439753578 at k = 2.940352436 in the AC-II model. The two Allen–Cahn models yield very similar growth behavior. (c) Two Cahn–Hilliard models, where λ1 = 10−4 m3s/kg, λ2 = 10−3 m3s/kg. The fastest growth rate is 403.3730366 at k = 1121.072073 in the CH-I model and 400.9478476 at k = 1124.400744 in the CH-II model. In general, we conclude that the phase field models predict quantitatively comparable results in near equilibrium dynamics in any parameter regimes.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In