Thickness-based adaptive mesh refinement methods for multi-phase flow simulations with thin regions

In numerical simulations of multi-scale, multi-phase flows, grid refinement is required to resolve regions with small scales. A notable example is liquid-jet atomization and subsequent droplet dynamics. It is essential to characterize the detailed flow physics with variable length scales with high fidelity, in order to elucidate the underlying mechanisms. In this paper, two thickness-based mesh refinement schemes are developed based on distance- and topology-oriented criteria for thin regions with confining wall/plane of symmetry and in any situation, respectively. Both techniques are implemented in a general framework with a volume-of-fluid formulation and an adaptive-mesh-refinement capability. The distance-oriented technique compares against a critical value, the ratio of an interfacial cell size to the distance between the mass center of the cell and a reference plane. The topology-oriented technique is developed from digital topology theories to handle more general conditions. The requirement for interfacial mesh refinement can be detected swiftly, without the need of thickness information, equation solving, variable averaging or mesh repairing. The mesh refinement level increases smoothly on demand in thin regions. The schemes have been verified and validated against several benchmark cases to demonstrate their effectiveness and robustness. These include the dynamics of colliding droplets, droplet motions in a microchannel, and atomization of liquid impinging jets. Overall, the thickness-based refinement technique provides highly adaptive meshes for problems with thin regions in an efficient and fully automatic manner. (C) 2014 Elsevier Inc. All rights reserved.

[1]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[2]  V. Cristini,et al.  Adaptive unstructured volume remeshing - II: Application to two- and three-dimensional level-set simulations of multiphase flow , 2005 .

[3]  Stéphane Popinet,et al.  An accurate adaptive solver for surface-tension-driven interfacial flows , 2009, J. Comput. Phys..

[4]  François Guibault,et al.  Three-dimensional anisotropic geometric metrics based on local domain curvature and thickness , 2005, Comput. Aided Des..

[5]  Vittorio Cristini,et al.  An adaptive mesh algorithm for evolving surfaces: simulation of drop breakup and coalescence , 2001 .

[6]  J. M. Bush,et al.  On the collision of laminar jets: fluid chains and fishbones , 2004, Journal of Fluid Mechanics.

[7]  Biao Zhou,et al.  Experimental and mechanistic description of merging and bouncing in head-on binary droplet collision , 2008 .

[8]  Vigor Yang,et al.  Impinging Jet Dynamics , 2012 .

[9]  S. Popinet,et al.  HIGH-FIDELITY SIMULATIONS OF IMPINGING JET ATOMIZATION , 2013 .

[10]  Héctor D. Ceniceros,et al.  A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation , 2010 .

[11]  R. LeVeque.  High-resolution conservative algorithms for advection in incompressible flow , 1996 .

[12]  M. Dupin,et al.  Simulation of a microfluidic flow-focusing device. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  David P. Schmidt,et al.  Adaptive tetrahedral meshing in free-surface flow , 2005 .

[14]  Shaoping Quan,et al.  Simulations of multiphase flows with multiple length scales using moving mesh interface tracking with adaptive meshing , 2011, J. Comput. Phys..

[15]  Ann S. Almgren,et al.  An adaptive level set approach for incompressible two-phase flows , 1997 .

[16]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[17]  Gihun Son,et al.  Numerical study of obstacle configuration for droplet splitting in a microchannel , 2013 .

[18]  Guoqun Zhao,et al.  Adaptive hexahedral mesh generation based on local domain curvature and thickness using a modified grid-based method , 2007 .

[19]  A. M. Khokhlov,et al.  Fully Threaded Tree for Adaptive Refinement Fluid Dynamics Simulations , 1998 .

[20]  Chung King Law,et al.  Regimes of coalescence and separation in droplet collision , 1997, Journal of Fluid Mechanics.

[21]  Michihisa Tsutahara,et al.  Three-dimensional lattice Boltzmann simulations of droplet formation in a cross-junction microchannel , 2008 .

[22]  Alexei M. Khokhlov,et al.  Fully Threaded Tree Algorithms for Adaptive Refinement Fluid Dynamics Simulations , 1997, astro-ph/9701194.

[23]  N. Ashgriz,et al.  Coalescence and separation in binary collisions of liquid drops , 1990, Journal of Fluid Mechanics.

[24]  Ricardo H. Nochetto,et al.  A hybrid variational front tracking-level set mesh generator for problems exhibiting large deformations and topological changes , 2010, J. Comput. Phys..

[25]  Nicolas Bremond,et al.  Decompressing emulsion droplets favors coalescence. , 2008, Physical review letters.

[26]  Markus Bussmann,et al.  Adaptive VOF with curvature‐based refinement , 2007 .

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

[28]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

[29]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[30]  Stéphane Popinet,et al.  A front-tracking algorithm for accurate representation of surface tension , 1999 .