Subgrid-scale modelling of surface tension within interface tracking-based Large Eddy and Interface Simulation of 3D interfacial flows

Abstract This paper introduces a novel modelling procedure for unresolved surface tension in interface tracking-based Large Eddy and Interface Simulation (LEIS) of interfacial flows. The approach involves modelling of the subgrid-scale (SGS) surface tension appearing in the filtered momentum equation, and explicitly the dominant SGS curvature part of that term. The procedure for modelling the SGS curvature is based on the analysis of variations in curvature resulting from using different filter widths, as well from using alternative discretization-based and direct-filtering techniques for generating coarse-filter curvature estimates. Unresolved surface tension modelling can act to smooth spatial curvature variations and dampen parasitic modes, and can alternatively restore interface wrinkling when super-grid scale resolution underpredicts local curvature. The paper outlines the principles involved in identifying local regions of unresolved surface tension, and in constructing a model that is robust across the range of interface topologies. The current implementation is based on the height-function curvature discretization scheme in Volume-of-Fluid (VOF) interfacial flow simulations. The procedure represents a SGS model in the true sense of the term – it is not a discretization scheme. The model should be used for Large Interface Simulation (LIS) of laminar interfacial flow modelling, too, complementing discretization schemes employed for resolved surface tension. Understandably the contribution of the model diminishes with mesh resolution as laminar flow does not continually introduce filter-scale interfacial deformations. The model may be more useful though for LEIS of turbulent interfacial flows, supporting discretization schemes for better capturing interface wrinkling, atomization and large bubble breakup.

[1]  M. Herrmann,et al.  An outline of an LES subgrid model for liquid / gas phase interface dynamics , 2008 .

[2]  Djamel Lakehal,et al.  Turbulence structure and interaction with steep breaking waves , 2011, Journal of Fluid Mechanics.

[3]  Test-case number 10 : Parasitic currents induced by surface tension ( PC ) March 2003 , 2004 .

[4]  Jean-François Remacle,et al.  A stabilized finite element method using a discontinuous level set approach for the computation of bubble dynamics , 2007, J. Comput. Phys..

[5]  P. Sagaut.  Large Eddy Simulation for Incompressible Flows , 2001 .

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

[7]  John R. Blake,et al.  Gas bubbles bursting at a free surface , 1993, Journal of Fluid Mechanics.

[8]  M. Rudman.  A Volume-Tracking Method for Incompressible Multifluid Flows with Large Density Variations , 1998 .

[9]  D. Kothe.  Perspective on Eulerian Finite Volume Methods for Incompressible Interfacial Flows , 1998 .

[10]  Nico de Jong,et al.  Nonspherical oscillations of ultrasound contrast agent microbubbles. , 2008, Ultrasound in medicine & biology.

[11]  O. Lebaigue,et al.  Jump conditions for filtered quantities at an under-resolved discontinuous interface. Part 1: Theoretical development , 2009 .

[12]  Djamel Lakehal,et al.  Large-eddy simulation of sheared interfacial flow , 2006 .

[13]  Djamel Lakehal,et al.  A Newton--Krylov Solver for Remapping-Based Volume-of-Fluid Methods , 2008, SIAM J. Sci. Comput..

[14]  M. Rudman,et al.  A 3D unsplit-advection volume tracking algorithm with planarity-preserving interface reconstruction , 2006 .

[15]  Djamel Lakehal.  LEIS for the prediction of turbulent multifluid flows applied to thermal-hydraulics applications , 2010 .

[16]  Olivier Simonin,et al.  DNS of the interaction between a deformable buoyant bubble and a spatially decaying turbulence: A priori tests for LES two-phase flow modelling , 2008 .

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

[18]  Djamel Lakehal,et al.  Multi-physics treatment in the vicinity of arbitrarily deformable gas-liquid interfaces , 2007, J. Comput. Phys..

[19]  Mark Sussman,et al.  High-order Techniques for Calculating Surface Tension Forces , 2022 .

[20]  F. Resch.  OCEANIC AIR BUBBLES AS GENERATORS OF MARINE AEROSOLS , 1986 .

[21]  D. Lakehal,et al.  Interface-turbulence interactions in large-scale bubbling processes , 2007 .

[22]  D. Jacqmin.  Regular Article: Calculation of Two-Phase Navier–Stokes Flows Using Phase-Field Modeling , 1999 .

[23]  Richard Manasseh,et al.  Cavitation microstreaming patterns in single and multiple bubble systems , 2007, Journal of Fluid Mechanics.

[24]  D. J. Torres,et al.  On the theory and computation of surface tension: the elimination of parasitic currents through energy conservation in the second-gradient method , 2002 .

[25]  D. B. Kothe,et al.  Convergence and accuracy of kernel-based continuum surface tension models , 1998 .

[26]  D. E. Spiel.  A hypothesis concerning the peak in film drop production as a function of bubble size , 1997 .

[27]  J. G. Hnat,et al.  Spherical cap bubbles and skirt formation , 1976 .

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

[29]  F. Resch,et al.  Distribution of the liquid aerosol produced from bursting bubbles in sea and distilled water , 1990 .

[30]  Markus Bussmann,et al.  Second-order accurate normals from height functions , 2008, J. Comput. Phys..

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

[32]  Axel Günther,et al.  Droplet production from disintegrating bubbles at water surfaces. Single vs. multiple bubbles , 2003 .

[33]  D. E. Spiel.  On the births of film drops from bubbles bursting on seawater , 1998 .

[34]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[35]  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..

[36]  F. Melo,et al.  Acoustic emission associated with the bursting of a gas bubble at the free surface of a non-Newtonian fluid. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[37]  Murray Rudman,et al.  Efficient simulation of surface tension-dominated flows through enhanced interface geometry interrogation , 2010, J. Comput. Phys..

[38]  S. Cummins,et al.  Estimating curvature from volume fractions , 2005 .

[39]  Pierre Sagaut,et al.  Towards large eddy simulation of isothermal two-phase flows: Governing equations and a priori tests , 2007 .

[40]  F. Faura,et al.  An improved height function technique for computing interface curvature from volume fractions , 2009 .

[41]  K. Hynynen,et al.  Ultrasound Enhanced Delivery of Molecular Imaging and Therapeutic Agents in Alzheimer's Disease Mouse Models , 2008, PloS one.

[42]  D. E. Spiel.  On the births of jet drops from bubbles bursting on water surfaces , 1995 .

[43]  Mark Sussman,et al.  A Stable and Efficient Method for Treating Surface Tension in Incompressible Two-Phase Flow , 2009, SIAM J. Sci. Comput..

[44]  Hiroshi Saeki,et al.  Three-dimensional vortex structures under breaking waves , 2005, Journal of Fluid Mechanics.

[45]  Jay P. Boris,et al.  On large eddy simulation using subgrid turbulence models Comment 1 , 1990 .

[46]  J. Kirby,et al.  Dynamics of surf-zone turbulence in a spilling breaker , 1996 .

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

[48]  Pierre Sagaut,et al.  Numerical simulation of phase separation and a priori two-phase LES filtering , 2008 .

[49]  A. H. Woodcock,et al.  Photographic Investigation of the Projection of Droplets by Bubbles Bursting at a Water Surface , 1954 .

[50]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .