An XFEM Implementation of a Projection Method for 3D Incompressible Two-Fluid Flows with Arbitrary High Contrasts in Material Properties

An XFEM Implementation of a Projection Method for 3D Incompressible Two-Fluid Flows with Arbitrary High Contrasts in Material Properties

Year:    2018

Communications in Computational Physics, Vol. 24 (2018), Iss. 3 : pp. 593–622

Abstract

This paper presents an XFEM implementation of a projection algorithm to compute in an Eulerian framework 3D incompressible two-fluid flows with arbitrary high contrasts in material properties. It is designed to deal with both strong and weak discontinuities across the interface for pressure and velocity fields, respectively. A classical enrichment function accounts for velocity gradient discontinuities across the interface and a new quadratic enrichment function accounts for pressure discontinuities across the interface. A splitting of two-fluid elements is performed to achieve accurate numerical integrations, meanwhile a scaling coefficient accounting for both physical and geometrical considerations alleviates ill-conditioning. Various validations have been carried and very good solution accuracy is achieved even on coarse meshes, as from the minimal mesh not conforming to the interface. This implementation enables to compute accurate solutions regardless of discontinuity magnitude (arbitrary high contrast in material properties) and mesh size of two-fluid elements, which can constitute a decisive advantage for large size 3D computations.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/cicp.OA-2017-0152

Communications in Computational Physics, Vol. 24 (2018), Iss. 3 : pp. 593–622

Published online:    2018-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    30

Keywords:    Incompressible two-fluid flows strong and weak field discontinuities Extended Finite Element Method projection algorithms.