Multi-physics treatment in the vicinity of arbitrarily deformable gas-liquid interfaces

A novel three-dimensional problem formulation is introduced for the simulation of turbulent interfacial multi-fluid flows. The strategy is built around the large eddy simulation (LES) concept, and can be employed for interfacial heat and mass transfer problems in which use can be made of either scalar transfer correlations, or exact mass/energy jump conditions. This multi-physics treatment capability at arbitrarily deformable interfaces translates into two main features: (i) a reconstructed distance function (RDF) is introduced to define a level-set interface-normal length scale, and (ii) an interfacial shear velocity is defined on the distance function support for further use in near-interface transport models. The solution algorithm uses VOF with piecewise planar interface reconstructions on a twice-as-fine mesh, and infers the convective mass fluxes from the interface solution to promote momentum conservation. The interfacial shear velocity defined on the distance function support is introduced to accommodate the asymptotic behaviour of turbulence approaching the interface in a proximity-dependent manner. Provided with highly accurate distance function data, the scheme generates near-interface damping functions that are second-order accurate and independent of interface orientation. The damping was found to be affected by errors introduced into shear velocity estimates by the high-frequency errors in the RDF scheme near the interface. The methodology has been applied for the simulation of a wave breaking scenario featuring multiple modes and interfacial length scales.

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

[2]  Pierre Lubin,et al.  Fully three-dimensional direct numerical simulation of a plunging breaker , 2003 .

[3]  Qun Zhao,et al.  Numerical simulation of breaking waves by a multi-scale turbulence model , 2004 .

[4]  Ernest R. Smith,et al.  Laboratory Study of Wave-Breaking Over Bars and Artificial Reefs , 1991 .

[5]  D. Lakehal.  Sub-grid scale modelling for the les of interfacial gas-liquid flows , 2005 .

[6]  Djamel Lakehal,et al.  Large-Eddy Simulation of Steep Water Waves , 2006 .

[7]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

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

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

[10]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[11]  Bernardus J. Geurts,et al.  A priori tests of large eddy simulation of the compressible plane mixing layer , 1995 .

[12]  Modélisation des écoulements avec interface gaz-liquide par simulation des grands tourbouillons (LES) des échelles de sous maille , 2005 .

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

[14]  A. Jenkins.  Do Strong Winds Blow Waves Flat , 2002 .

[15]  J. Kirby,et al.  Observation of undertow and turbulence in a laboratory surf zone , 1994 .

[16]  Djamel Lakehal,et al.  Direct numerical simulation of turbulence in a sheared air–water flow with a deformable interface , 2003, Journal of Fluid Mechanics.

[17]  D. Chopp.  Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[18]  D. B. Kothe,et al.  Modeling High Density Ratio Incompressible Interfacial Flows , 2002 .

[19]  Qiang Zhang,et al.  Three-Dimensional Front Tracking , 1998, SIAM J. Sci. Comput..

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

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

[22]  L. Rossi,et al.  Investigation of wall normal electromagnetic actuator for seawater flow control , 2002 .

[23]  Djamel Lakehal,et al.  Surface divergence models for scalar exchange between turbulent streams , 2004 .

[24]  Djamel Lakehal,et al.  Large-eddy simulation of bubbly turbulent shear flows , 2002 .

[25]  L YoungsD,et al.  Time-dependent multi-material flow with large fluid distortion. , 1982 .

[26]  Billy L. Edge,et al.  Ocean Wave Measurement and Analysis , 1994 .

[27]  P. Woodward,et al.  SLIC (Simple Line Interface Calculation) , 1976 .

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

[29]  V. S. Vaidhyanathan,et al.  Transport phenomena , 2005, Experientia.

[30]  Djamel Lakehal,et al.  LES of Turbulent Bubble Formation and Break-Up Based on Interface Tracking , 2004 .

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

[32]  John D. Fenton,et al.  A Fifth‐Order Stokes Theory for Steady Waves , 1985 .

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

[34]  Murray Rudman,et al.  A volume of fluid (VOF) method for the simulation of metallurgical flows , 2001 .

[35]  Erik Damgaard Christensen,et al.  Large eddy simulation of breaking waves , 2001 .

[36]  G. Yadigaroglu,et al.  Direct Numerical Simulation of turbulent heat transfer across a mobile, sheared gas-liquid interface , 2003 .

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