On two-phase flow solvers in irregular domains with contact line

Mathieu Lepilliez (1, 2, 3) , Elena Roxana Popescu (1) , Frédéric Gibou (4) , Sébastien Tanguy (1)

1 IMFT - Institut de mécanique des fluides de Toulouse
2 CNES - Centre National d'Études Spatiales [Toulouse]
3 Airbus [France]
4 UC Santa Barbara - University of California [Santa Barbara]

We present numerical methods that enable the direct numerical simulation of two-phase flows in irregular domains. A method is presented to account for surface tension effects in a mesh cell containing a triple line between the liquid, gas and solid phases. Our numerical method is based on the level-set method to capture the liquid–gas interface and on the single-phase Navier–Stokes solver in irregular domain proposed in [35]to impose the solid boundary in an Eulerian framework. We also present a strategy for the implicit treatment of the viscous term and how to impose both a Neumann boundary condition and a jump condition when solving for the pressure field. Special care is given on how to take into account the contact angle, the no-slip boundary condition for the velocity field and the volume forces. Finally, we present numerical results in two and three spatial dimensions evaluating our simulations with several benchmarks.

Mots clés: Sharp interface methods; Irregular domains; Implicit viscosity; Contact lines; Level set method; Ghost-Fluid Method

Domaines: Mécanique des fluides [physics.class-ph]; Physique Numérique [physics.comp-ph]

Mathieu Lepilliez, Elena Roxana Popescu, Frédéric Gibou, Sébastien Tanguy. On two-phase flow solvers in irregular domains with contact line. Journal of Computational Physics, 2016, 321, pp.1217-1251. ⟨10.1016/j.jcp.2016.06.013⟩. ⟨hal-01349346⟩