Hostname: page-component-78c5997874-j824f Total loading time: 0 Render date: 2024-11-10T15:50:04.076Z Has data issue: false hasContentIssue false
  • English
  • Français

Mass and Volume Conservation in Phase Field Models for Binary Fluids

Published online by Cambridge University Press:  03 June 2015

Jie Shen ,
Xiaofeng Yang  and
Qi Wang
Jie Shen*
Affiliation:
Department of Mathematics, Purdue University, West Lafayette, IN 46907, USA
Xiaofeng Yang*
Affiliation:
Department of Mathematics and NanoCenter at USC, University of South Carolina, Columbia, SC 29028, USA
Qi Wang*
Affiliation:
Department of Mathematics and NanoCenter at USC, University of South Carolina, Columbia, SC 29028, USA School of Mathematics, Nankai University, Tianjin 300071, PR. China Beijing Computational Science Research Center, Beijing 100084, PR. China
*
Email:shen@math.purdue.edu
Email:xfyang@math.sc.edu
*Corresponding author.Email:qwang@math.sc.edu
Get access

Abstract

The commonly used incompressible phase field models for non-reactive, binary fluids, in which the Cahn-Hilliard equation is used for the transport of phase variables (volume fractions), conserve the total volume of each phase as well as the material volume, but do not conserve the mass of the fluid mixture when densities of two components are different. In this paper, we formulate the phase field theory for mixtures of two incompressible fluids, consistent with the quasi-compressible theory [28], to ensure conservation of mass and momentum for the fluid mixture in addition to conservation of volume for each fluid phase. In this formulation, the mass-average velocity is no longer divergence-free (solenoidal) when densities of two components in the mixture are not equal, making it a compressible model subject to an internal con-straint. In one formulation of the compressible models with internal constraints (model 2), energy dissipation can be clearly established. An efficient numerical method is then devised to enforce this compressible internal constraint. Numerical simulations in confined geometries for both compressible and the incompressible models are carried out using spatially high order spectral methods to contrast the model predictions. Numerical comparisons show that (a) predictions by the two models agree qualitatively in the situation where the interfacial mixing layer is thin; and (b) predictions differ significantly in binary fluid mixtures undergoing mixing with a large mixing zone. The numerical study delineates the limitation of the commonly used incompressible phase field model using volume fractions and thereby cautions its predictive value in simulating well-mixed binary fluids.

Keywords

65M7076M2276T99Phase field modelcompressibilitymultiphase fluid flowsspectral methods
Type
Research Article
Information
Communications in Computational Physics , Volume 13 , Issue 4 , April 2013 , pp. 1045 - 1065
Copyright
Copyright © Global Science Press Limited 2013

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Bird, , Stewart, , and Lightfoot, , Transport Phenomena, John Wiley and Sons, 2002.Google Scholar
[2]Probstein, , Physicochemical Hydrodynamics, John Wiley and Sons, 1994.CrossRefGoogle Scholar
[3]Bird, B., Armstrong, R., Hassager, O., Dynamics of Polymeric Liquids, 2nd Ed., Vol. 2, John Wiley and Sons, New York, 1987.Google Scholar
[4]Beris, A. N. and Edwards, B., Thermodynamics of Flowing Systems, Oxford University Press, Oxford, 1994.Google Scholar
[5]Lindley, B., Wang, Q. and Zhang, T., Multicomponent models for biofilm flows, Discrete and Continuous Dynamic Systems- Series B, 15(2) (2011), 417–456.Google Scholar
[6]Cahn, J. W. and Hilliard, J. E.Free energy of a nonuniform system. i: interfacial free energy. J. Chem. Phys., 28 (1959), 258–267.Google Scholar
[7]Cahn, J. W. and Hilliard, J. E.Free energy of a nonuniform system-iii: Nucleation in a 2- component incompressible fluid. J. Chem. Phys., 31(3) (1959), 688–699.CrossRefGoogle Scholar
[8] L. Chen, Q. and Yang, W., Computer simulation of the dynamics of a quenched system with large number of non-conserved order parameters, Phys. Rev. B, 50 (1994), 15752–15756.CrossRefGoogle Scholar
[9]Chen, L. Q., Phase-field modeling for microstructure evolution, Annu. Rev. Mater. Res., 32 (2002), 113–140.CrossRefGoogle Scholar
[10]Chen, L. Q. and Wang, Y., The Continuum Field Approach to Modeling Microstructural Evolution, J. Miner Met. Mater. Soc. 48 (12) (1996), 13–18.CrossRefGoogle Scholar
[11]Doi, M., Introduction to Polymer Physics, Clarendon Press, Oxford, 1996.Google Scholar
[12]Du, Q., Liu, C., Ryham, R. and Wang, X., Phase field modeling of the spontaneous curvature effect in cell membranes, Comm. Pur. Applied. Anal., 4 (2005), 537–548.CrossRefGoogle Scholar
[13]Du, Q., Liu, C., Ryham, R. and Wang, X., A phase field formulation of the Willmore problem, Nonlinearity, 18 (2005), 1249–1267.CrossRefGoogle Scholar
[14]Du, Q., Liu, C. and Wang, X., A Phase Field Approach in the Numerical Study of the Elastic Bending Energy for Vesicle Membranes, J. Comp. Phy., 198 (2004), 450–468.Google Scholar
[15]Du, Q., Liu, C. and Wang, X., Retrieving topological information for phase field models, SIAM Journal on Applied Mathematics, 65 (2005), 1913–1932.Google Scholar
[16]Du, Q., Liu, C. and Wang, X., Simulating the Deformation of Vesicle Membranes under Elastic Bending Energy in Three Dimensions, J. Comp. Phys., 212 (2005), 757–777.Google Scholar
[17]Forest, M. G. and Wang, Q., Hydrodynamic theories for blends of flexible polymer and nematic polymers, Physical Review E, 72 (2005), 041805.Google Scholar
[18]Forest, M. G., Liao, Q. and Wang, Qi, 2-D Kinetic Theory for Polymer Particulate Nanocomposites, Communication in Computational Physics, 7 (2) (2010), 250–282.CrossRefGoogle Scholar
[19]Feng, J. J., Liu, C., Shen, J. and Yue, P., Transient Drop Deformation upon Startup of Shear in Viscoelastic Fluids, Fluids. Phys. Fluids, 17 (2005), 123101.Google Scholar
[20]Flory, P. J., Principles of Polymer Chemistry, Cornell University Press, Ithaca, NY, 1953.Google Scholar
[21]Hobayashi, R., Modeling and numerical simulations of dendritic crystal growth, Physica D, 63 (1993), 410–423.CrossRefGoogle Scholar
[22]Karma, A. and Rappel, W., Phase-Field Model of Dendritic Sidebranching with Thermal Noise, Phys. Rev. E, 60 (1999), 3614–3625.Google Scholar
[23]Hua, l. Jinsong, Lin, Ping, Chun, , Liu, , Wang, Qi, Energy Law Preserving C° Finite Element Schemes for Phase Field Models in Two-phase Flow Computations, J. Comp. Phys., 230 (19) (2011), 7115–7131.Google Scholar
[24]Liu, C. and Walkington, N. J., An Eulerian description of fluids containing visco-hyperelastic particles, Arch. Rat. Mech. Ana., 159 (2001), 229–252.Google Scholar
[25]Liu, C. and Shen, J., A phase field model for the mixture of two incompressible fluids and its approximation by a fourier-spectral method, Physica D, 179 (2003), 211–228.Google Scholar
[26]Li, Y., Hu, S., Liu, Z., and Chen, L., Phase-field model of domain structures in ferroelectric thin films, Appl. Phys. Lett., 78 (2001), 3878–3880.CrossRefGoogle Scholar
[27]Lu, W. and Suo, Z., Dynamics of nanoscale pattern formation of an epitaxial monolayer, J. Mech. Phys. Solids, 49 (2001), 1937–1950.CrossRefGoogle Scholar
[28]Lowengrub, J. and Truskinovsky, L., Quasi-incompressible Cahn-Hilliard fluids and topological transitions, R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci., 454 (1998), 2617–2654.CrossRefGoogle Scholar
[29]McFadden, G., Wheeler, A., Braun, R., Coriell, S., and Sekerka, R., Phys. Rev. E, 48 (1998), 20162024.Google Scholar
[30]Seol, D. J., Hu, S. Y., Li, Y. L., Shen, J., Oh, K. H. and Chen, L. Q., Three-dimensional Phase-Field Modeling of Spinodal Decomposition in Constrained Films, Acta Materialia 51 (2003), 5173–5185.CrossRefGoogle Scholar
[31]Shen, J. and Yang, X., An efficient moving mesh spectral method for the phase-field model of two phase flows, J. Comput. Phys., 228 (2009), 2978–2992.Google Scholar
[32]Shen, J. and Yang, X., Energy Stable Schemes for Cahn-Hilliard phase-field model of two- phase incompressible flows, Chinese Ann. Math. series B, 31 (2010), 743–758.Google Scholar
[33]Shen, J. and Yang, X., A phase-field model and its numerical approximation for two-phase incompressible flows with different densities and viscositites, SIAM J. Sci. Comput., 32(3) (2010), 1159–1179.Google Scholar
[34]Tadmor, E., Phillips, R., and Ortiz, M., Mixed Atomistic and Continuum Models of Deformation in Solids, Langmuir, 12 (1996), 4529–4534.Google Scholar
[35]Wang, X., Asymptotic Analysis of Phase Field Formulations of Bending Elasticity Models, submitted to SIAM Mathematical Analysis (2006).Google Scholar
[36]Wang, Q., E, W., Liu, C., and Zhang, P., Kinetic theories for flows of nonhomogeneous rodlike liquid crystalline polymers with a nonlocal intermolecular potential, Physical Review E, 65(5) (2002), 0515041–0515047.Google Scholar
[37]Wang, Q., A hydrodynamic theory of nematic liquid crystalline polymers of different configurations, Journal of Chemical Physics, 116 (2002), 9120–9136.Google Scholar
[38]Wang, Q., Forest, M. G. and Zhou, R., A hydrodynamic theory for solutions of nonhomogeneous nematic liquid crystalline polymers with density variations, J. of Fluid Engineering, 126 (2004), 180–188.Google Scholar
[39]Wang, Y. and Chen, C. L., Simulation of microstructure evolution. In Methods in Materials Research, N. Ksufmann, Ed. E., Abbaschian, R., Bocarsly, A., Chien, C. L., Dollimore, D., et al., (1999), 2a3.1–2a3.23.Google Scholar
[40]Wheeler, A., McFadden, G., and Boettinger, W., Proc. R. Soc. London Ser. A, 452 (1996), 495525.Google Scholar
[41]Wise, S. M., Lowengrub, J. S., Kim, J. S. and Johnson, W. C., Efficient phase-field simulation of quantum dot formation in a strained heteroepitaxial film, Superlattices and Microstructures, 36 (2004) 293–304.CrossRefGoogle Scholar
[42]Yang, X., Feng, J., Liu, C. and Shen, J., Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method, J. Comput. Phys., 218 (2006), 417–428.Google Scholar
[43]Yue, P., Feng, J. J., Liu, C., and Shen, J., A diffuse-interface method for simulating two-phase flows of complex fluids, J. Fluid Mech., 515 (2004), 293–317.Google Scholar
[44]Yue, P., Feng, J. J., Liu, C., and Shen, J., Diffuse-interface simulations of drop coalescence and retraction in viscoelastic fluids, J. Non-Newtonian Fluid Mech., 129 (2005), 163–176.Google Scholar
[45]Zhang, T. Y., Cogan, N., and Wang, Q., Phase Field Models for Biofilms. II. 2-D Numerical Simulations of Biofilm-Flow Interaction, Communications in Computational Physics, 4 (2008), 72–101.Google Scholar
[46]Zhang, T. Y. and Wang, Q., Cahn-Hilliard vs Singular Cahn-Hilliard Equations in Phase Field Modeling, Communications in Computational Physics, 7(2) (2010), 362–382.Google Scholar