SIMULATION OF 2-FLUID FLOWS BY THE LEAST-SQUARES FINITE-ELEMENT METHOD USING A CONTINUUM SURFACE-TENSION MODEL

Authors
Citation
J. Wu et al., SIMULATION OF 2-FLUID FLOWS BY THE LEAST-SQUARES FINITE-ELEMENT METHOD USING A CONTINUUM SURFACE-TENSION MODEL, International journal for numerical methods in engineering, 42(4), 1998, pp. 583-600
Citations number
31
Categorie Soggetti
Mathematics,Engineering,Mathematics
ISSN journal
00295981
Volume
42
Issue
4
Year of publication
1998
Pages
583 - 600
Database
ISI
SICI code
0029-5981(1998)42:4<583:SO2FBT>2.0.ZU;2-1
Abstract
In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols(1) and the Continuum Surface Force (CSF) model developed by Brackbill et al.(2) In the VOF method fluids of di fferent properties are identified through the use of a continuous held variable (colour function). The colour function assigns a unique cons tant (colour) to each fluid. The interfaces between different fluids a re distinct due to sharp gradients of the colour function. The evoluti on of the interfaces is captured by solving the convective equation of the colour function. The CSF model is used as a means to treat surfac e tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin,(3) is used to calculate the tension force . In the modified version, the force term is obtained by calculating t he divergence of a stress tensor defined by the gradient of the colour function. In its analytical form, this stress formulation is equivale nt to the original CSF model.(2) Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the int erfaces. The least-squares finite element method (LSFEM)(4) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present app lications. The LSFEM handles all the equations in a unified manner wit hout any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both t wo-dimensional planar and axisymmetric flows, including a dam breaking , oscillating and stationary bubbles and a conical liquid sheet in a p ressure swirl atomizer. (C) 1998 John Wiley & Sons, Ltd.