Marker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow around Irregular Geometries

Computational simulations of multiphase flow are challenging because many practical applications require adequate resolution of not only interfacial physics associated with moving boundaries with possible topological changes, but also around three-dimensional, irregular solid geometries. In this paper, we focus on the simulations of fluid/fluid dynamics around complex geometries, based on an Eulerian-Lagrangian framework. The approach uses two independent but related grid layouts to track the interfacial and solid boundary conditions, and is capable of capturing interfacial as well as multiphase dynamics. In particular, the stationary Cartesian grid with time dependent, local adaptive refinement is utilized to handle the computation of the transport equations, while the interface shape and movement are treated by marker-based triangulated surface meshes which freely move and interact with the Cartesian grid. The markers are also used to identify the location of solid boundaries and enforce the no-slip condition there. Issues related to the contact line treatment, topological changes of multiphase fronts during merger or breakup of objects, and necessary data structures and solution techniques are also highlighted. Selected test cases including spacecraft fuel tank flow management, and movement and rupture of interfaces associated with liquid plug flow are presented.

[1]  E. P. Symons,et al.  Draining characteristics of hemispherically bottomed cylinders in a low-gravity environment , 1978 .

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

[3]  J. Grotberg,et al.  Effects of inertia and gravity on liquid plug splitting at a bifurcation. , 2006, Journal of biomechanical engineering.

[4]  Roel Luppes,et al.  The numerical simulation of liquid sloshing on board spacecraft , 2007, J. Comput. Phys..

[5]  Thomas Erneux,et al.  Nonlinear rupture of free films , 1993 .

[6]  Mark Sussman,et al.  A sharp interface method for incompressible two-phase flows , 2007, J. Comput. Phys..

[7]  A. Chorin.  Numerical Solution of the Navier-Stokes Equations* , 1989 .

[8]  Zhu Wang,et al.  THE LEVEL SET METHOD ON ADAPTIVE CARTESIAN GRID FOR INTERFACE CAPTURING , 2004 .

[9]  Frank Losasso,et al.  A fast and accurate semi-Lagrangian particle level set method , 2005 .

[10]  Gretar Tryggvason,et al.  Computational Methods for Multiphase Flow: Frontmatter , 2007 .

[11]  Wei Shyy,et al.  Micro-scale drop dynamics for heat transfer enhancement , 2002 .

[12]  Wei Shyy,et al.  COMPUTATIONS OF DROP DYNAMICS WITH THE IMMERSED BOUNDARY METHOD, PART 1: NUMERICAL ALGORITHM AND BUOYANCY-INDUCED EFFECT , 2003 .

[13]  S. Osher,et al.  Spatially adaptive techniques for level set methods and incompressible flow , 2006 .

[14]  Zhilin Li,et al.  The immersed interface method for the Navier-Stokes equations with singular forces , 2001 .

[15]  Said I. Abdel-Khalik,et al.  Accurate representation of surface tension using the level contour reconstruction method , 2005 .

[16]  Wei Shyy,et al.  Multiphase computations using sharp and continuous interface techniques for micro-gravity applications , 2004 .

[17]  Ramesh Nallapati,et al.  A moving unstructured staggered mesh method for the simulation of incompressible free-surface flows , 2003 .

[18]  W. Shyy,et al.  A Unified Adaptive Cartesian Grid Method for Solid- Multiphase Fluid Dynamics with Moving Boundaries , 2007 .

[19]  W. Shyy,et al.  Regular Article: An Accurate Cartesian Grid Method for Viscous Incompressible Flows with Complex Immersed Boundaries , 1999 .

[20]  Matthew W. Williams,et al.  A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework , 2006, J. Comput. Phys..

[21]  Gretar Tryggvason,et al.  Numerical simulation of dendritic solidification with convection: two-dimensional geometry , 2002 .

[22]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[23]  B. Engquist,et al.  Discretization of Dirac delta functions in level set methods , 2005 .

[24]  S. G. Berenyi,et al.  Vapor ingestion phenomenon in hemispherically bottomed tanks in normal gravity and in weightlessness , 1970 .

[25]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[26]  Wei Shyy,et al.  Computations of Multiphase Fluid Flows Using Marker-Based Adaptive, Multilevel Cartesian Grid Method , 2007 .

[27]  X. Jiang,et al.  Numerical simulation of the head-on collision of two equal-sized drops with van der Waals forces , 2007 .

[28]  L. J. Hastings,et al.  An experimental study of the behavior of a sloshing liquid subjected to a sudden reduction in acceleration , 1968 .

[29]  Xiaolin Li,et al.  Robust Computational Algorithms for Dynamic Interface Tracking in Three Dimensions , 1999, SIAM J. Sci. Comput..

[30]  Wei Shyy,et al.  Computational Modeling for Fluid Flow and Interfacial Transport (Dover Books on Engineering) , 1993 .

[31]  Luis Gustavo Nonato,et al.  A front-tracking/front-capturing method for the simulation of 3D multi-fluid flows with free surfaces , 2004 .

[32]  Brian Wetton,et al.  Computation of a Moving Drop/Bubble on a Solid Surface using a Front-Tracking Method , 2004 .

[33]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[34]  Qiang Zhang,et al.  Front tracking in two and three dimensions , 1998 .

[35]  Gretar Tryggvason,et al.  Computations of Boiling Flows , 2004 .

[36]  Wei Shyy,et al.  Computational modeling for multiphase flows with spacecraft application , 2007 .

[37]  H. Udaykumar,et al.  Sharp interface Cartesian grid method I: An easily implemented technique for 3D moving boundary computations , 2005 .

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

[39]  Cheng-Feng Tai,et al.  Multiphase/Multidomain Computations Using Continuum and Lattice–Boltzmann Methods , 2006 .

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

[41]  Meng-Sing Liou,et al.  A robust and accurate approach to computing compressible multiphase flow: Stratified flow model and AUSM+-up scheme , 2007, J. Comput. Phys..

[42]  Gretar Tryggvason,et al.  Numerical simulation of dendritic solidification with convection: three-dimensional flow , 2004 .

[43]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[44]  Hideki Fujioka,et al.  Steady propagation of a liquid plug in a two-dimensional channel. , 2004, Journal of biomechanical engineering.

[45]  Clinton P. T. Groth,et al.  A Mesh Adjustment Scheme for Embedded Boundaries , 2006 .

[46]  Marc Garbey,et al.  Three-dimensional adaptive, Cartesian grid method for multiphase flow computations , 2005 .

[47]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[48]  Seungwon Shin,et al.  Modeling three-dimensional multiphase flow using a level contour reconstruction method for front tracking without connectivity , 2002 .

[49]  S. Abdel-Khalik,et al.  Direct three-dimensional numerical simulation of nucleate boiling using the level contour reconstruction method , 2005 .

[50]  F. Sotiropoulos,et al.  A hybrid Cartesian/immersed boundary method for simulating flows with 3D, geometrically complex, moving bodies , 2005 .

[51]  H. S. Udaykumar,et al.  A Sharp Interface Cartesian Grid Methodfor Simulating Flows with ComplexMoving Boundaries , 2001 .

[52]  Rajkeshar Singh,et al.  Three-dimensional marker-based multiphase flow computation using adaptive Cartesian grid techniques , 2006 .

[53]  B. Wetton,et al.  Analysis and computation of immersed boundaries, with application to pulp fibres , 1997 .

[54]  Wei Shyy,et al.  Three-dimensional adaptive Cartesian grid method with conservative interface restructuring and reconstruction , 2007, J. Comput. Phys..

[55]  Raj M. Manglik,et al.  Short-Time-Transient Surfactant Dynamics and Marangoni Convection Around Boiling Nuclei , 2003 .

[56]  L. M. Hocking.  A moving fluid interface on a rough surface , 1976, Journal of Fluid Mechanics.

[57]  M. Aftosmis.  Solution adaptive cartesian grid methods for aerodynamic flows with complex geometries , 1997 .

[58]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[59]  Elias Balaras,et al.  An embedded-boundary formulation for large-eddy simulation of turbulent flows interacting with moving boundaries , 2006, J. Comput. Phys..

[60]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .