An efficient numerical method for simulating multiphase flows using a diffuse interface model

This paper presents a new diffuse interface model for multiphase incompressible immiscible fluid flows with surface tension and buoyancy effects. In the new model, we employ a new chemical potential that can eliminate spurious phases at binary interfaces, and consider a phase-dependent variable mobility to investigate the effect of the mobility on the fluid dynamics. We also significantly improve the computational efficiency of the numerical algorithm by adapting the recently developed scheme for the multiphase-field equation. To illustrate the robustness and accuracy of the diffuse interface model for surface tension- and buoyancy-dominant multi-component fluid flows, we perform numerical experiments, such as equilibrium phase-field profiles, the deformation of drops in shear flow, a pressure field distribution, the efficiency of the proposed scheme, a buoyancy-driven bubble in ambient fluids, and the mixing of a six-component mixture in a gravitational field. The numerical result obtained by the present model and solution algorithm is in good agreement with the analytical solution and, furthermore, we not only remove the spurious phase-field profiles, but also improve the computational efficiency of the numerical solver.

[1]  Ingo Steinbach,et al.  A generalized field method for multiphase transformations using interface fields , 1999 .

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

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

[4]  Junseok Kim,et al.  Accurate contact angle boundary conditions for the Cahn–Hilliard equations , 2011 .

[5]  Patrick Patrick Anderson,et al.  Capillary spreading of a droplet in the partially wetting regime using a diffuse-interface model , 2007, Journal of Fluid Mechanics.

[6]  A. Tornberg,et al.  A numerical method for two phase flows with insoluble surfactants , 2011 .

[7]  Patrick Tabeling,et al.  Droplet breakup in microfluidic junctions of arbitrary angles. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  Junseok Kim,et al.  On the long time simulation of the Rayleigh–Taylor instability , 2011 .

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

[10]  Patrick Patrick Anderson,et al.  Cahn–Hilliard modeling of particles suspended in two‐phase flows , 2012 .

[11]  M. Renardy,et al.  PROST: a parabolic reconstruction of surface tension for the volume-of-fluid method , 2002 .

[12]  Vassili S. Sochnikov,et al.  Level set calculations of the evolution of boundaries on a dynamically adaptive grid , 2003 .

[13]  Xiaofeng Yang,et al.  Mass and Volume Conservation in Phase Field Models for Binary Fluids , 2013 .

[14]  C. L. Tucker,et al.  Microstructural evolution in polymer blends , 2003 .

[15]  Junseok Kim,et al.  Phase field computations for ternary fluid flows , 2007 .

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

[17]  Lingling Shui,et al.  Multiphase flow in microfluidic systems --control and applications of droplets and interfaces. , 2007, Advances in colloid and interface science.

[18]  Junseok Kim,et al.  Numerical simulation of the three-dimensional Rayleigh-Taylor instability , 2013, Comput. Math. Appl..

[19]  Jie Shen,et al.  A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method , 2003 .

[20]  Mario De Menech.  Modeling of droplet breakup in a microfluidic T-shaped junction with a phase-field model , 2006 .

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

[22]  U. Ozen,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

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

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

[25]  E. Puckett,et al.  A High-Order Projection Method for Tracking Fluid Interfaces in Variable Density Incompressible Flows , 1997 .

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

[27]  M. Sommerfeld,et al.  Multiphase Flows with Droplets and Particles , 2011 .

[28]  T. N. Stevenson,et al.  Fluid Mechanics , 2021, Nature.

[29]  Harald Garcke,et al.  On anisotropic order parameter models for multi-phase system and their sharp interface limits , 1998 .

[30]  I. Steinbach,et al.  A phase field concept for multiphase systems , 1996 .

[31]  Hans-Jörg Bart,et al.  Investigating the process of liquid-liquid extraction by means of computational fluid dynamics , 1996 .

[32]  Patrick Patrick Anderson,et al.  Diffuse Interface Modeling of Droplet Impact on a Pre‐Patterned Solid Surface , 2005 .

[33]  Sankaran Sundaresan,et al.  Modeling the hydrodynamics of multiphase flow reactors: Current status and challenges , 2000 .

[34]  Arash Karimipour,et al.  Mixed convection of copper-water nanofluid in a shallow inclined lid driven cavity using the lattice Boltzmann method , 2014 .

[35]  Andrea Mazzino,et al.  Phase-field model for the Rayleigh–Taylor instability of immiscible fluids , 2008, Journal of Fluid Mechanics.

[36]  J. Taylor,et al.  Overview no. 113 surface motion by surface diffusion , 1994 .

[37]  Wei Shyy,et al.  Multiphase Dynamics in Arbitrary Geometries on Fixed Cartesian Grids , 1997 .

[38]  S. Foroughi,et al.  Lattice Boltzmann method on quadtree grids for simulating fluid flow through porous media: A new automatic algorithm , 2013 .

[39]  John-Chang Chen,et al.  Onset of entrainment between immiscible liquid layers due to rising gas bubbles , 1988 .

[40]  G. Taylor.  The instability of liquid surfaces when accelerated in a direction perpendicular to their planes. I , 1950, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[41]  John W. Cahn,et al.  On Spinodal Decomposition , 1961 .

[42]  S. Zaleski,et al.  Volume-of-Fluid Interface Tracking with Smoothed Surface Stress Methods for Three-Dimensional Flows , 1999 .

[43]  J. Strutt.  Scientific Papers: Investigation of the Character of the Equilibrium of an Incompressible Heavy Fluid of Variable Density , 2009 .

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

[45]  J. Dukowicz.  A particle-fluid numerical model for liquid sprays , 1980 .

[46]  Jang-Min Park,et al.  A ternary model for double-emulsion formation in a capillary microfluidic device. , 2012, Lab on a chip.

[47]  D. Weitz,et al.  Monodisperse Double Emulsions Generated from a Microcapillary Device , 2005, Science.

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

[49]  R. Chella,et al.  Mixing of a two-phase fluid by cavity flow. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[50]  Britta Nestler,et al.  Phase-field model for solidification of a monotectic alloy with convection , 2000 .

[51]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy and Free Energy of a Nonuniform System. III. Nucleation in a Two‐Component Incompressible Fluid , 2013 .

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

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

[54]  S. Osher,et al.  A Level Set Formulation of Eulerian Interface Capturing Methods for Incompressible Fluid Flows , 1996 .

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

[56]  Junseok Kim,et al.  Regularized Dirac delta functions for phase field models , 2012 .

[57]  Britta Nestler,et al.  A multi-phase-field model of eutectic and peritectic alloys: numerical simulation of growth structures , 2000 .

[58]  Junseok Kim,et al.  A comparison study of the Boussinesq and the variable density models on buoyancy-driven flows , 2012 .

[59]  Jun Li,et al.  A Class of Conservative Phase Field Models for Multiphase Fluid Flows , 2014 .

[60]  Vittorio Cristini,et al.  A new volume-of-fluid formulation for surfactants and simulations of drop deformation under shear at a low viscosity ratio , 2002 .

[61]  J. López,et al.  On the reinitialization procedure in a narrow‐band locally refined level set method for interfacial flows , 2005 .

[62]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[63]  Britta Nestler,et al.  A multigrid solver for phase field simulation of microstructure evolution , 2008, Math. Comput. Simul..

[64]  Junseok Kim,et al.  Buoyancy-driven mixing of multi-component fluids in two-dimensional tilted channels , 2013 .

[65]  Franck Boyer,et al.  Numerical schemes for a three component Cahn-Hilliard model , 2011 .