An accurate methodology for surface tension modeling in OpenFOAM

by   Abd Essamade Saufi, et al.

In this paper a numerical methodology for surface tension modeling is presented, with an emphasis on the implementation in the OpenFOAM framework. The methodology relies on a combination of (i) a well-balanced approach based on the Ghost Fluid Method (GFM), including the jump of density and pressure directly in the numerical discretization of the pressure equation, and (ii) Height Functions to evaluate the interface curvature, implemented, to the authors' knowledge, for the first time in OpenFOAM. The method is able to significantly reduce spurious currents (almost to machine accuracy) for a stationary droplet, showing second order convergence both for the curvature and the interface shape. Accurate results are also obtained for additional test cases such as translating droplets, capillary oscillations and rising bubbles, for which numerical results are comparable to what obtained by other numerical codes in the same conditions. Finally, the Height Functions method is extended to include the treatment of contact angles, both for sessile droplets and droplets suspended under the effect of gravity, showing a very good agreement with the theoretical prediction. The code works in parallel mode and details on the actual implementation in OpenFOAM are included to facilitate the reproducibility of the results.


page 5

page 16

page 18

page 19

page 21

page 40

page 41


A multigrid ghost-point level-set method for incompressible Navier-Stokes equations on moving domains with curved boundaries

In this paper we present a numerical approach to solve the Navier-Stokes...

A Tikhonov approach to level set curvature computation

In numerical simulations of two-phase flows, the computation of the curv...

Modeling blood flow in networks of viscoelastic vessels with the 1D augmented fluid-structure interaction system

A noteworthy aspect in blood flow modeling is the definition of the mech...

Degenerate area preserving surface Allen-Cahn equation and its sharp interface limit

We consider formal matched asymptotics to show the convergence of a dege...

Assessment of CFD capability for prediction of the Coandă effect

The tendency of a jet to stay attached to a flat or convex surface is ca...

Capillary Rise -- A Computational Benchmark for Wetting Processes

Four different numerical approaches are compared for the rise of liquid ...

Please sign up or login with your details

Forgot password? Click here to reset