Generalized conservative phase field model and its lattice Boltzmann scheme for multicomponent multiphase flows

Abstract In this paper a generalized conservative phase field-lattice Boltzmann model is proposed for immiscible incompressible flows with any number of components. In this phase field model, the each phase field variable satisfies a generalized eikonal equation and the volume fraction constraint. By introducing a quadratic functional, the dynamics of phase field variables can be treated as a gradient flow method. The phase field equations are reformulated based on the mass conservative law and the reduction-consistency conditions. Moreover, based on the generalized free energy functional, a family of generalized continuous surface tension force formulations is deduced by using the virtual work method. A lattice Boltzmann model (LBM) is developed to solve the generalized conservative phase field equations and the hydrodynamics equations. In the numerical simulation section, the motion of circular interfaces and Poiseuille flow with four phases are simulated first. The convergence rate and reduction-consistent of present model for the phase field equations are validated. The computational accuracy of present model for the hydrodynamics equations is also validated. Then serval two-dimensional (2D) examples with pairwise surface tension effect, including the droplets immersed in another fluid with four fluid phases, spreading of liquid lenses and spinodal decomposition with four fluid phases, are simulated. The obtained numerical results are in good agreement with the analytical solutions. Four continuous surface tension force formulations are also comparative studied. Finally, as a practical application, a three-dimensional (3D) bubble rising in a three-layer liquid is investigated.

[1]  Hongkai Zhao,et al.  A fast sweeping method for Eikonal equations , 2004, Math. Comput..

[2]  Abbas Fakhari,et al.  Conservative phase-field lattice-Boltzmann model for ternary fluids , 2018, J. Comput. Phys..

[3]  Elie Bretin,et al.  A modified phase field approximation for mean curvature flow with conservation of the volume , 2009, 0904.0098.

[4]  Shuonan Wu,et al.  Multiphase Allen-Cahn and Cahn-Hilliard models and their discretizations with the effect of pairwise surface tensions , 2016, J. Comput. Phys..

[5]  Samuel M. Allen,et al.  Mechanisms of phase transformations within the miscibility gap of Fe-rich Fe-Al alloys , 1976 .

[6]  Zhenhua Chai,et al.  A comparative study of local and nonlocal Allen-Cahn equations with mass conservation , 2018, International Journal of Heat and Mass Transfer.

[7]  Li-Shi Luo,et al.  A weighted multiple-relaxation-time lattice Boltzmann method for multiphase flows and its application to partial coalescence cascades , 2017, J. Comput. Phys..

[8]  Ross T. Whitaker,et al.  A Fast Iterative Method for Eikonal Equations , 2008, SIAM J. Sci. Comput..

[9]  S. He,et al.  Phase-field-based lattice Boltzmann model for incompressible binary fluid systems with density and viscosity contrasts. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[11]  Junseok Kim,et al.  A generalized continuous surface tension force formulation for phase-field models for multi-component immiscible fluid flows , 2009 .

[12]  Chang Shu,et al.  Diffuse interface model for incompressible two-phase flows with large density ratios , 2007, J. Comput. Phys..

[13]  Héctor D. Ceniceros,et al.  Computation of multiphase systems with phase field models , 2002 .

[14]  Wenjing Ye,et al.  A Flux-Corrected Phase-Field Method for Surface Diffusion , 2017 .

[15]  Song Zheng,et al.  Phase-field-theory-based lattice Boltzmann equation method for N immiscible incompressible fluids. , 2019, Physical review. E.

[16]  Geoffrey Ingram Taylor,et al.  The formation of emulsions in definable fields of flow , 1934 .

[17]  Junseok Kim,et al.  An efficient and accurate numerical algorithm for the vector-valued Allen-Cahn equations , 2012, Comput. Phys. Commun..

[18]  Jie Shen,et al.  Modeling and numerical approximation of two-phase incompressible flows by a phase-field approach , 2011 .

[19]  Abbas Fakhari,et al.  Phase-field modeling by the method of lattice Boltzmann equations. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  Junseok Kim,et al.  A second-order accurate non-linear difference scheme for the N -component Cahn–Hilliard system , 2008 .

[21]  Junseok Kim,et al.  A new conservative vector-valued Allen-Cahn equation and its fast numerical method , 2017, Comput. Phys. Commun..

[22]  J. Rubinstein,et al.  Nonlocal reaction−diffusion equations and nucleation , 1992 .

[23]  Raoyang Zhang,et al.  A Lattice Boltzmann Scheme for Incompressible Multiphase Flow and Its Application in Simulation of Rayleigh-Taylor Instability , 1998 .

[24]  Chunfeng Zhou,et al.  Spontaneous shrinkage of drops and mass conservation in phase-field simulations , 2007, J. Comput. Phys..

[25]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[26]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[27]  P. B. Mumford,et al.  Emulsions , 1944 .

[28]  Ying Sun,et al.  Sharp interface tracking using the phase-field equation , 2007, J. Comput. Phys..

[29]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[30]  D. Chopp,et al.  A projection method for motion of triple junctions by level sets , 2002 .

[31]  Jung-Il Choi,et al.  Multi-component Cahn-Hilliard system with different boundary conditions in complex domains , 2016, J. Comput. Phys..

[32]  Junseok Kim,et al.  Comparison study of the conservative Allen-Cahn and the Cahn-Hilliard equations , 2016, Math. Comput. Simul..

[33]  Alfredo Soldati,et al.  Mass-conservation-improved phase field methods for turbulent multiphase flow simulation , 2018, Acta Mechanica.

[34]  Frederick Stern,et al.  A highly scalable massively parallel fast marching method for the Eikonal equation , 2015, J. Comput. Phys..

[35]  Suchuan Dong,et al.  Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters , 2014, J. Comput. Phys..

[36]  Z. Chai,et al.  Lattice Boltzmann modeling of three-phase incompressible flows. , 2016, Physical review. E.

[37]  James J. Feng,et al.  A diffuse-interface method for simulating two-phase flows of complex fluids , 2004, Journal of Fluid Mechanics.

[38]  Abbas Fakhari,et al.  Improved locality of the phase-field lattice-Boltzmann model for immiscible fluids at high density ratios. , 2017, Physical review. E.

[39]  D. M. Anderson,et al.  DIFFUSE-INTERFACE METHODS IN FLUID MECHANICS , 1997 .

[40]  Franck Boyer,et al.  Study of a three component Cahn-Hilliard flow model , 2006 .

[41]  Junseok Kim,et al.  An efficient numerical method for simulating multiphase flows using a diffuse interface model , 2015 .

[42]  Junseok Kim,et al.  A practically unconditionally gradient stable scheme for the N-component Cahn-Hilliard system , 2012 .

[43]  Franck Boyer,et al.  Hierarchy of consistent n-component Cahn–Hilliard systems , 2014 .

[44]  Jung-Il Choi,et al.  A phase-field fluid modeling and computation with interfacial profile correction term , 2016, Commun. Nonlinear Sci. Numer. Simul..

[45]  Junseok Kim,et al.  Phase field modeling and simulation of three-phase flows , 2005 .

[46]  Yang Hu,et al.  Hybrid Allen-Cahn-based lattice Boltzmann model for incompressible two-phase flows: The reduction of numerical dispersion. , 2019, Physical review. E.

[47]  J. Lowengrub,et al.  Conservative multigrid methods for Cahn-Hilliard fluids , 2004 .

[48]  Manuel Quezada de Luna,et al.  A monolithic conservative level set method with built-in redistancing , 2019, J. Comput. Phys..

[49]  Suchuan Dong,et al.  Wall-bounded multiphase flows of N immiscible incompressible fluids: Consistency and contact-angle boundary condition , 2016, J. Comput. Phys..

[50]  Chunming Li,et al.  Level set evolution without re-initialization: a new variational formulation , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[51]  Tomohiro Takaki,et al.  Multi-phase-field modeling using a conservative Allen–Cahn equation for multiphase flow , 2019, Computers & Fluids.

[52]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[53]  Pao-Hsiung Chiu,et al.  A conservative phase field method for solving incompressible two-phase flows , 2011, J. Comput. Phys..

[54]  Suchuan Dong,et al.  Multiphase flows of N immiscible incompressible fluids: A reduction-consistent and thermodynamically-consistent formulation and associated algorithm , 2017, J. Comput. Phys..

[55]  Suchuan Dong,et al.  An efficient algorithm for incompressible N-phase flows , 2014, J. Comput. Phys..

[56]  Harald Garcke,et al.  A MultiPhase Field Concept: Numerical Simulations of Moving Phase Boundaries and Multiple Junctions , 1999, SIAM J. Appl. Math..

[57]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .