Buoyancy-driven mixing of multi-component fluids in two-dimensional tilted channels

Abstract Buoyancy-driven mixing of multi-component incompressible immiscible fluids in two-dimensional tilted channels is studied numerically using a phase-field model. This paper extends the previous work [K.C. Sahu, S.P. Vanka, A multiphase lattice Boltzmann study of buoyancy-induced mixing in a tilted channel, Comput. Fluids 50 (2011) 199–215] to the multi-component (more than two) fluid case. The mixing dynamics are governed by the modified Navier–Stokes equations and the multi-component convective Cahn–Hilliard equations. A finite difference method is used to discretize the governing system. To solve the equations efficiently and accurately, we employ Chorin’s projection method for the modified Navier–Stokes equations, and the recently developed practically unconditionally stable method for the multi-component Cahn–Hilliard equations. We numerically investigate the effects of various density ratios, tilt angles, Reynolds numbers, and Weber numbers on the interface structures and front velocities. The trends observed in simulations with multi-component fluids are consistent with previous numerical results for two-component fluids.

[1]  E. Hinch,et al.  Transient buoyancy-driven front dynamics in nearly horizontal tubes , 2007 .

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

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

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

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

[6]  E. Mello,et al.  Numerical study of the Cahn–Hilliard equation in one, two and three dimensions , 2004, cond-mat/0410772.

[7]  Richard Welford,et al.  A multigrid finite element solver for the Cahn-Hilliard equation , 2006, J. Comput. Phys..

[8]  Héctor D. Ceniceros,et al.  A nonstiff, adaptive mesh refinement-based method for the Cahn-Hilliard equation , 2007, J. Comput. Phys..

[9]  E. Hinch,et al.  Front dynamics and macroscopic diffusion in buoyant mixing in a tilted tube , 2007 .

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

[11]  E. Hinch,et al.  Buoyancy driven miscible front dynamics in tilted tubes , 2004 .

[12]  A. Peirce,et al.  Unsteady axial mixing by natural convection in a vertical column , 1992 .

[13]  D. J. Eyre,et al.  An Unconditionally Stable One-Step Scheme for Gradient Systems , 1997 .

[14]  H. Ding,et al.  Pressure-driven miscible two-fluid channel flow with density gradients , 2008 .

[15]  I. Frigaard,et al.  Influence of an imposed flow on the stability of a gravity current in a near horizontal duct , 2010 .

[16]  E. Hinch,et al.  Laser-induced fluorescence measurements of buoyancy driven mixing in tilted tubes , 2006 .

[17]  Junseok Kim,et al.  A numerical method for the Cahn–Hilliard equation with a variable mobility , 2007 .

[18]  A. Prasad,et al.  Instability due to Viscosity Stratification Downstream of a Centerline Injector , 2008 .

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

[20]  Qiang Du,et al.  A Fourier Spectral Moving Mesh Method for the Cahn-Hilliard Equation with Elasticity , 2009 .

[21]  Bahar Firoozabadi,et al.  2D numerical simulation of density currents using the SPH projection method , 2013 .

[22]  Andreas Prohl,et al.  Numerical analysis of the Cahn-Hilliard equation and approximation for the Hele-Shaw problem, Part II: Error analysis and convergence of the interface , 2001 .

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

[24]  T. Benjamin.  Gravity currents and related phenomena , 1968, Journal of Fluid Mechanics.

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

[26]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

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

[28]  C. M. Elliott,et al.  On the Cahn-Hilliard equation with degenerate mobility , 1996 .

[29]  Stuart B. Dalziel,et al.  Gravity currents produced by lock exchange , 2004, Journal of Fluid Mechanics.

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

[31]  Daniel D. Joseph,et al.  CORE-ANNULAR FLOWS , 1997 .

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

[33]  Franck Boyer,et al.  Numerical study of viscoelastic mixtures through a Cahn–Hilliard flow model , 2004 .

[34]  Chakrabarti,et al.  Numerical study of the Cahn-Hilliard equation in three dimensions. , 1988, Physical review letters.

[35]  W. Rockwell Geyer,et al.  Gravity currents: In the environment and the laboratory , 1989 .

[36]  Jie Shen,et al.  Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[37]  S. Vanka,et al.  A multiphase lattice Boltzmann study of buoyancy-induced mixing in a tilted channel , 2011 .

[38]  Y. Renardy.  Stability of the interface in two-layer couette flow of upper convected maxwell liquids , 1988 .

[39]  Charles M. Elliott,et al.  The Cahn–Hilliard equation with a concentration dependent mobility: motion by minus the Laplacian of the mean curvature , 1996, European Journal of Applied Mathematics.

[40]  Surya Pratap Vanka,et al.  Multiphase lattice Boltzmann simulations of buoyancy-induced flow of two immiscible fluids with different viscosities , 2012 .

[41]  J. Magnaudet,et al.  Effects of channel geometry on buoyancy-driven mixing , 2008 .

[42]  E. Hinch,et al.  Buoyant mixing of miscible fluids in tilted tubes , 2004 .

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