Surface tension implementation for Gensmac 2D

Filho, Antonio Castelo and Mangiavacchi, Norberto and Tomé, Murilo F. and Cuminato, José A. and Fortuna, Armando O. and de Oliveira, Juliana and Ferreira, Valdemir G. and McKee, S. (2001) Surface tension implementation for Gensmac 2D. Journal of the Brazilian Society of Mechanical Sciences, 23 (4). pp. 523-532. ISSN 0100-7386 (https://doi.org/10.1590/S0100-73862001000400013)

[thumbnail of strathprints002129]
Preview
Text. Filename: strathprints002129.pdf
Accepted Author Manuscript
License: Creative Commons Attribution-NonCommercial 4.0 logo

Download (296kB)| Preview

Abstract

In the present work we describe a method which allows the incorporation of surface tension into the GENSMAC2D code. This is achieved on two scales. First on the scale of a cell, the surface tension effects are incorporated into the free surface boundary conditions through the computation of the capillary pressure. The required curvature is estimated by fitting a least square circle to the free surface using the tracking particles in the cell and in its close neighbors. On a sub-cell scale, short wavelength perturbations are filtered out using a local 4-point stencil which is mass conservative. An efficient implementation is obtained through a dual representation of the cell data, using both a matrix representation, for ease at identifying neighbouring cells, and also a tree data structure, which permits the representation of specific groups of cells with additional information pertaining to that group. The resulting code is shown to be robust, and to produce accurate results when compared with exact solutions of selected fluid dynamic problems involving surface tension.