New Implementations of Surface Tension Forces for PLIC-VOF Methods

2010 ◽  
Author(s):  
Jesu´s Castro ◽  
Assensi Oliva ◽  
Eduardo Garci´a-Rivera
Keyword(s):  
Surface Tension ◽  
Surface Tension Force ◽  
Surface Force ◽  
Polar Coordinates ◽  
Staggered Grid ◽  
Reconstruction Algorithms ◽  
Interface Location ◽  
Parasitic Currents ◽  
Force Calculation ◽  
Surface Tension Forces

A comparative study has been carried out between different models of implementing surface tension for PLIC-VOF methods: i) the well known Continuum Surface Force (CSF); ii) Staggered Grid Interface Pressure (SGIP); iii) Modified Meier’s. As main difference with respect the CSF model, the two last use information of the interface location for the calculation of the surface tension forces, therefore, they need a model of interface reconstruction. The models will be tested under static and dynamic conditions, for the case of a drop and bubble between water and air. In order to reduce the ‘parasitic’ currents that is characteristic of this type of models, two different strategies are used: i) a new method for the calculation of the interface curvature based on the use of polar coordinates; ii) the use of kernels for smoothing the jump of the colour function and the interface. Additionally, the influence of the order of accuracy of the reconstruction algorithms employed that affect the curvature estimation and the surface tension force calculation is also checked. After testing the different cases, it can be concluded that SGIP and modified Meier perform in a comparable way, and much better than CSF model without kernel, that is only competitive with the use of the kernels. However, modified Meier is preferred respect SGIP due to its straitghforward implementation on unstructured meshes.

scholarly journals Comparison of Surface Tension Models for the Volume of Fluid Method

Processes ◽  
2019 ◽  
Vol 7 (8) ◽  
pp. 542 ◽  
Author(s):  
Kurian J. Vachaparambil ◽  
Kristian Etienne Einarsrud
Keyword(s):  
Surface Tension ◽  
Multiphase Flow ◽  
Capillary Rise ◽  
Surface Tension Force ◽  
Surface Force ◽  
Parallel Plates ◽  
Time Step ◽  
Mesh Resolution ◽  
Active Research ◽  
Spurious Currents

With the increasing use of Computational Fluid Dynamics to investigate multiphase flow scenarios, modelling surface tension effects has been a topic of active research. A well known associated problem is the generation of spurious velocities (or currents), arising due to inaccuracies in calculations of the surface tension force. These spurious currents cause nonphysical flows which can adversely affect the predictive capability of these simulations. In this paper, we implement the Continuum Surface Force (CSF), Smoothed CSF and Sharp Surface Force (SSF) models in OpenFOAM. The models were validated for various multiphase flow scenarios for Capillary numbers of 10 − 3 –10. All the surface tension models provide reasonable agreement with benchmarking data for rising bubble simulations. Both CSF and SSF models successfully predicted the capillary rise between two parallel plates, but Smoothed CSF could not provide reliable results. The evolution of spurious current were studied for millimetre-sized stationary bubbles. The results shows that SSF and CSF models generate the least and most spurious currents, respectively. We also show that maximum time step, mesh resolution and the under-relaxation factor used in the simulations affect the magnitude of spurious currents.

scholarly journals Relaxation and breakup of a cylindrical liquid column

1970 ◽  
Vol 39 (2) ◽  
pp. 57-64 ◽  
Author(s):  
Mohammad Ali ◽  
Akira Umemura
Keyword(s):  
Surface Tension ◽  
Stokes Equations ◽  
Capillary Wave ◽  
Three Dimensional ◽  
Surface Tension Force ◽  
Surface Force ◽  
Computer Code ◽  
Vof Method ◽  
Liquid Column ◽  
Navier Stokes

Instability of capillary wave and breakup of a square cylindrical liquid column during its relaxation have been investigated numerically by simulating three-dimensional Navier-Stokes equations. For this investigation a computer code based on volume-of-fluid (VOF) method has been developed and validated with published experimental results. The result shows that the agreement of numerical simulation is quite well with the experimental data. The code is then used to study the capillary wave and breakup phenomena of the liquid column. The investigation shows the underlying physics during relaxation of the square cylindrical liquid column, illustrates the formation and propagation of capillary wave, and breakup processes. The breakup behavior for the present configuration of the liquid column shows some significant differences from those predicted by conventional jet atomization theories. The formation of capillary wave is initiated by the surface tension on the sharp edge of the square end of the cylinder and the propagation of the wave occurs due to the effect of surface tension force on the motion of the fluid. The propagation of capillary wave to the end of liquid column causes a disturbance in the system and makes the waves unstable which initiates the breakup of the liquid column. The characteristics of the capillary wave show that the amplitude of the swell grows faster than the neck of the wave and that of the tip wave grows much faster than the other waves. The velocity of the liquid particle is dominated by the pressure in the liquid column. Keywords: Instability; Continuum surface force; Liquid disintegration; Capillary wave; Surface tension; VOF method doi:10.3329/jme.v39i2.1847 Journal of Mechanical Engineering, Vol. ME39, No. 2, Dec. 2008 57-64

Numerical Study on Impingement of a Droplet upon a Liquid Film

2010 ◽  
Vol 44-47 ◽  
pp. 2499-2503
Author(s):  
Hong Liu ◽  
Mao Zhao Xie ◽  
Su Chun Wang ◽  
Ming Jia
Keyword(s):  
Surface Tension ◽  
Liquid Film ◽  
Stokes Equations ◽  
Numerical Study ◽  
Surface Tension Force ◽  
Surface Force ◽  
Navier Stokes ◽  
Initial Kinetic Energy ◽  
Droplet Impingement ◽  
Fine Grid

This paper reports progress in the numerical simulations of a droplet impingement upon the wall film of the same liquid. The full Navier-Stokes equations are solved in axisymmetric formulation. The surface tension force is modeled by a continuum surface force (CSF) model. An adapting local refinement technique is used to provide the fine grid coupled by the volume-of fluid (VOF) method for tracking the interface between the gas and the droplet and liquid film. Results indicate that the motion behavior of droplet impingement upon the liquid film is dominantly influenced by the initial kinetic energy and the thickness of the film as well as the surface tension and the liquid viscosity.

A Numerical Method for Capillarity-Driven Free Surface Flows

2005 ◽  
Author(s):  
Albert Y. Tong ◽  
Zhaoyuan Wang
Keyword(s):  
Surface Tension ◽  
Free Surface ◽  
Level Set ◽  
Surface Tension Force ◽  
Surface Force ◽  
Volume Of Fluid ◽  
Free Surface Flows ◽  
Tension Force ◽  
Navier Stokes Equation ◽  
Spurious Currents

The continuum surface force (CSF) method has been extensively employed in the volume-of-fluid (VOF), level set (LS) and front tracking methods to model surface tension force. It is a robust method requiring relatively easy implementation. However, it is known to generate spurious currents near the interface that may lead to disastrous interface instabilities and failures of grid convergence. A different surface tension implementation algorithm, referred to as the pressure boundary method (PBM), is introduced in this study. The surface tension force is incorporated into the Navier-Stokes equation via a capillary pressure gradient while the free surface is tracked by a coupled level set and volume-of-fluid (CLSVOF) method. It has been shown that the spurious currents are greatly reduced by the present method with the sharp pressure boundary condition preserved. The numerical results of several cases have been compared with data reported in the literature and are found to be in a close agreement.

A Level Set Method for Simulation of Coalescence of Droplets

2005 ◽  
Author(s):  
A. Salih ◽  
S. Ghosh Moulic
Keyword(s):  
Surface Tension ◽  
Level Set Method ◽  
Level Set ◽  
Stokes Equations ◽  
Surface Tension Force ◽  
Tension Force ◽  
Staggered Grid ◽  
Navier Stokes ◽  
Level Set Function ◽  
Set Function

In the present paper, we discuss a numerical method based on the level set algorithm to simulate two-phase fluid flow systems. Surface tension force at the fluid interface is implemented through the CSF model of Brackbill et al. [1]. The incompressible Navier-Stokes equations were solved on a staggered grid using an explicit projection method. A fifth-order WENO [2] scheme was used for advancing the level set function. We improved the implementation of WENO scheme by staggering the level set function. The Navier-Stokes part of the code was validated by computing the standard lid-driven cavity flow and the free surface part of the code was validated by advecting the interface in a prescribed velocity field. The Young-Laplace law for a static drop has been verified to validate the implementation of surface tension force. We simulated the coalescence of two drops under zero-gravity condition and evaluated the mass conservation property of the level set method.

scholarly journals On sharp surface force model: Effect of sharpening coefficient

2020 ◽  
Vol 3 (3) ◽  
pp. 226-232 ◽  
Author(s):  
Kurian J. Vachaparambil ◽  
Kristian Etienne Einarsrud
Keyword(s):  
Surface Tension ◽  
Capillary Rise ◽  
Surface Tension Force ◽  
Surface Force ◽  
Volume Of Fluid ◽  
Tension Force ◽  
Force Model ◽  
Relative Ease

Abstract Amongst the multitude of approaches available in literature to reduce spurious velocities in Volume of Fluid approach, the Sharp Surface Force (SSF) model is increasingly being used due to its relative ease to implement. The SSF approach relies on a user-defined parameter, the sharpening coefficient, which determines the extent of the smeared nature of interface used to determine the surface tension force. In this paper, we use the SSF model implemented in OpenFOAM® to investigate the effect of this sharpening coefficient on spurious velocities and accuracy of dynamic, i.e., capillary rise, and static bubble simulations. Results show that increasing the sharpening coefficient generally reduces the spurious velocities in both static and dynamic cases. Although static millimeter sized bubbles were simulated with the whole range of sharpening coefficients, sub-millimeter sized bubbles show nonphysical behavior for values larger than 0.3. The accuracy of the capillary rise simulations has been observed to change non-linearly with the sharpening coefficient. This work illustrates the importance of using an optimized value of the sharpening coefficient with respect to spurious velocities and accuracy of the simulation.

Numerical Analysis of Heat Transfer of a Flow Confined by Wire Screen in Lithium Bromide Absorption Process

2014 ◽  
Vol 348 ◽  
pp. 40-50 ◽  
Author(s):  
Herbert Obame Mve ◽  
Romuald Rullière ◽  
Rémi Goulet ◽  
Phillippe Haberschill
Keyword(s):  
Heat Transfer ◽  
Surface Tension ◽  
Numerical Analysis ◽  
Lithium Bromide ◽  
Surface Tension Force ◽  
Surface Force ◽  
Geometrical Parameters ◽  
Force Model ◽  
Vof Model ◽  
Static Contact

The present study deals with the numerical analysis of heat transfer inside a lithium bromide solution flowing down between finely meshed plastic wire screens. These screens confine the flow through capillary action while allowing the water vapor transfer inside an innovative absorber technology. The complex menisci shape formed on the confinement grid level, where the surface tension forces are of first importance, are reconstructed by a volume of fluid (VOF) model. A continuum surface force model is used to account for the surface tension force. A static contact angle is used to define the wall adhesion. A new algorithm, consisting to set an unique constant temperature at the liquid/vapour interface and to determine the evolution of heat transfer characteristics over the simulation domain, has been implemented and validated by analytical solutions. A parametric study has been conducted to determine the effect of the inlet velocity and the geometrical parameters (wire diameter and the number of divisions).

INFLUENCE OF SURFACE TENSION FORCES ON FILMWISE CONDENSATION INTENSITY

2016 ◽  
Vol 7 (1-2) ◽  
pp. 93-102
Author(s):  
Vasily Buz ◽  
Konstantin Goncharov ◽  
Henry F. Smirnov
Keyword(s):  
Surface Tension ◽  
Surface Tension Forces
Sign in / Sign up

Export Citation Format

Share Document