Computing curvature for volume of fluid methods using machine learning

Abstract In spite of considerable progress, computing curvature in Volume of Fluid (VOF) methods continues to be a challenge. The goal is to develop a function or a subroutine that returns the curvature in computational cells containing an interface separating two immiscible fluids, given the volume fraction in the cell and the adjacent cells. Currently, the most accurate approach is to fit a curve (2D), or a surface (3D), matching the volume fractions and finding the curvature by differentiation. Here, a different approach is examined. A synthetic data set, relating curvature to volume fractions, is generated using well-defined shapes where the curvature and volume fractions are easily found and then machine learning is used to fit the data (training). The resulting function is used to find the curvature for shapes not used for the training and implemented into a code to track moving interfaces. The results suggest that using machine learning to generate the relationship is a viable approach that results in reasonably accurate predictions.

[1]  Gretar Tryggvason,et al.  Direct Numerical Simulations of Gas–Liquid Multiphase Flows: Preface , 2011 .

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

[3]  R. Scardovelli,et al.  A mixed markers and volume-of-fluid method for the reconstruction and advection of interfaces in two-phase and free-boundary flows , 2003 .

[4]  Sandro Manservisi,et al.  On the properties and limitations of the height function method in two-dimensional Cartesian geometry , 2011, J. Comput. Phys..

[5]  R. Scardovelli,et al.  A surface marker algorithm coupled to an area-preserving marker redistribution method for three-dimensional interface tracking , 2004 .

[6]  S. Zaleski,et al.  Modelling Merging and Fragmentation in Multiphase Flows with SURFER , 1994 .

[7]  Brian L. Smith,et al.  A novel technique for including surface tension in PLIC-VOF methods , 2002 .

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

[9]  Olivier Desjardins,et al.  A mesh-decoupled height function method for computing interface curvature , 2015, J. Comput. Phys..

[10]  Sandro Manservisi,et al.  Vofi - A library to initialize the volume fraction scalar field , 2016, Comput. Phys. Commun..

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

[12]  S. Popinet.  Numerical Models of Surface Tension , 2018 .

[13]  Sandro Manservisi,et al.  Numerical integration of implicit functions for the initialization of the VOF function , 2015 .

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

[15]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

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

[17]  M. Renardy,et al.  PROST: a parabolic reconstruction of surface tension for the volume-of-fluid method , 2002 .