Picture of neon light reading 'Open'

Discover open research at Strathprints as part of International Open Access Week!

23-29 October 2017 is International Open Access Week. The Strathprints institutional repository is a digital archive of Open Access research outputs, all produced by University of Strathclyde researchers.

Explore recent world leading Open Access research content this Open Access Week from across Strathclyde's many research active faculties: Engineering, Science, Humanities, Arts & Social Sciences and Strathclyde Business School.

Explore all Strathclyde Open Access research outputs...

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

[img]
Preview
Text (strathprints002129)
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.