Multiphase Dynamics in Arbitrary Geometries on Fixed Cartesian Grids

In this work, a mixed Eulerian?Lagrangian algorithm, calledELAFINT(Eulerian Lagrangian algorithm for interface tracking) is developed further and applied to compute flows with solid?fluid and fluid?fluid interfaces. The method is capable of handling fluid flows in the presence of both irregularly shaped solid boundaries and moving boundaries on a fixed Cartesian grid. The field equations are solved on the underlying fixed grid using a collocated variable, pressure-based formulation. The moving boundary is tracked explicitly by the Lagrangian translation of marker particles. The moving boundary passes through the grid and the immersed boundary technique is used to handle its interaction with the underlying grid. The internal solid boundaries are dealt with by using a cut-cell technique. Particular attention is directed toward conservation and consistency in the vicinity of both phase boundaries. The complex geometry feature has been tested for a variety of flow problems. The performance of the immersed boundary representation is demonstrated in the simulation of Newtonian liquid drops. The combination of the two features is then employed in the simulation of motion of drops through constricted tubes. The capabilities developed here can be useful for solving flow problems involving moving and stationary complex boundaries.

[1]  D. B. Kothe,et al.  RIPPLE: A NEW MODEL FOR INCOMPRESSIBLE FLOWS WITH FREE SURFACES , 1991 .

[2]  Bram van Leer,et al.  A simulation technique for 2-D unsteady inviscid flows around arbitrarily moving and deforming bodies of arbitrary geometry , 1993 .

[3]  John B. Bell,et al.  Cartesian grid method for unsteady compressible flow in irregular regions , 1995 .

[4]  Wei Shyy,et al.  Development of pressure-based composite multigrid methods for complex fluid flows , 1996 .

[5]  L. Sirovich,et al.  Modeling a no-slip flow boundary with an external force field , 1993 .

[6]  W. Olbricht,et al.  The deformation and breakup of liquid drops in low Reynolds number flow through a capillary , 1992 .

[7]  In Seok Kang,et al.  Numerical solution of axisymmetric, unsteady free‐boundary problems at finite Reynolds number. I. Finite‐difference scheme and its application to the deformation of a bubble in a uniaxial straining flow , 1987 .

[8]  N. Ashgriz,et al.  A spine-flux method for simulating free surface flows , 1995 .

[9]  Lawrence Sirovich,et al.  Direct numerical simulation of turbulent flow over a modeled riblet covered surface , 1995, Journal of Fluid Mechanics.

[10]  Dynamics of drops in cavity flows: Aggregation of high viscosity ratio drops , 1996 .

[11]  E. Bender.  Numerical heat transfer and fluid flow. Von S. V. Patankar. Hemisphere Publishing Corporation, Washington – New York – London. McGraw Hill Book Company, New York 1980. 1. Aufl., 197 S., 76 Abb., geb., DM 71,90 , 1981 .

[12]  Howard A. Stone,et al.  Relaxation and breakup of an initially extended drop in an otherwise quiescent fluid , 1989, Journal of Fluid Mechanics.

[13]  Gretar Tryggvason,et al.  Bifurcation of tracked scalar waves , 1986 .

[14]  Stanley Osher,et al.  A Hybrid Method for Moving Interface Problems with Application to the Hele-Shaw Flow , 1997 .

[15]  Hideaki Miyata,et al.  Finite-difference simulation of breaking waves , 1986 .

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

[17]  Wei Shyy,et al.  A GRID-SUPPORTED MARKER PARTICLE SCHEME FOR INTERFACE TRACKING , 1995 .

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

[19]  L. G. Leal,et al.  Numerical solution of free-boundary problems in fluid mechanics. Part 1. The finite-difference technique , 1984, Journal of Fluid Mechanics.

[20]  W. Shyy,et al.  Simulation of interfacial instabilities during solidification. I: Conduction and capillarity effects , 1995 .

[21]  S. Majumdar.  Role of underrelaxation in momentum interpolation for calculation of flow with nonstaggered grids , 1988 .

[22]  D. P. Young,et al.  A locally refined rectangular grid finite element method: application to computational fluid dynamics and computational physics , 1990 .

[23]  Fue-Sang Lien,et al.  A general non-orthogonal collocated finite volume algorithm for turbulent flow at all speeds incorporating second-moment turbulence-transport closure, Part 1: Computational implementation , 1994 .

[24]  J. Sethian.  Level set methods : evolving interfaces in geometry, fluid mechanics, computer vision, and materials science , 1996 .

[25]  L. Fauci,et al.  A computational model of aquatic animal locomotion , 1988 .

[26]  Wei Shyy,et al.  Computational Fluid Dynamics with Moving Boundaries , 1995 .

[27]  W. Shyy,et al.  ELAFINT - A mixed Eulerian-Lagrangian method for fluid flows with complex and moving boundaries , 1994 .

[28]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[29]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[30]  D. Juric,et al.  A Front-Tracking Method for Dendritic Solidification , 1996 .

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

[32]  Howard A. Stone,et al.  The influence of initial deformation on drop breakup in subcritical time-dependent flows at low Reynolds numbers , 1989, Journal of Fluid Mechanics.

[33]  D. D. Zeeuw,et al.  An adaptively refined Cartesian mesh solver for the Euler equations , 1993 .

[34]  Michael J. Miksis,et al.  Dynamics of a drop in a constricted capillary tube , 1994, Journal of Fluid Mechanics.

[35]  Charles S. Peskin,et al.  Improved Volume Conservation in the Computation of Flows with Immersed Elastic Boundaries , 1993 .

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

[37]  Wei Shyy,et al.  Calculation of meniscus shapes and transport processes in float zone , 1995 .

[38]  J. Quirk.  An alternative to unstructured grids for computing gas dynamic flows around arbitrarily complex two-dimensional bodies , 1994 .

[39]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[40]  Gretar Tryggvason,et al.  Head-on collision of drops—A numerical investigation , 1993 .

[41]  P. Roe.  Book Review Computational Modelling for Fluid Flow and Interfacial Transport , 1994 .

[42]  M. Manga.  Dynamics of drops in branched tubes , 1996, Journal of Fluid Mechanics.

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

[44]  C. Peskin.  Numerical analysis of blood flow in the heart , 1977 .

[45]  Vladimir Hlavacek,et al.  Numerical Solution of Free-Boundary Problems , 1996 .

[46]  L. G. Leal,et al.  The creeping motion of immiscible drops through a converging/diverging tube , 1983, Journal of Fluid Mechanics.