Simulation of two-fluid flows by the least-squares finite element method using a continuum surface tension model
✍ Scribed by Jie Wu; Sheng-Tao Yu; Bo-Nan Jiang
- Publisher
- John Wiley and Sons
- Year
- 1998
- Tongue
- English
- Weight
- 286 KB
- Volume
- 42
- Category
- Article
- ISSN
- 0029-5981
No coin nor oath required. For personal study only.
✦ Synopsis
In this paper a numerical procedure for simulating two-uid ows 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 uids of di erent properties are identiÿed through the use of a continuous ÿeld variable (colour function). The colour function assigns a unique constant (colour) to each uid. The interfaces between di erent uids are distinct due to sharp gradients of the colour function. The evolution of the interfaces is captured by solving the convective equation of the colour function. The CSF model is used as a means to treat surface tension e ect at the interfaces. Here a modiÿed version of the CSF model, proposed by Jacqmin, 3 is used to calculate the tension force. In the modiÿed version, the force term is obtained by calculating the divergence of a stress tensor deÿned by the gradient of the colour function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the di culty in approximating the curvatures of the interfaces.
The least-squares ÿnite element method (LSFEM) 4 is used to discretize the governing equation systems. The LSFEM has proven to be e ective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a uniÿed manner without any additional special treatment such as upwinding or artiÿcial dissipation.
Various bench mark tests have been carried out for both two-dimensional planar and axisymmetric ows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer. ?