On a marching level-set method for extended discontinuous Galerkin methods for incompressible two-phase flows

Martin Smuda,Florian Kummer

DOI: https://doi.org/10.48550/arXiv.2010.08417

2020-10-16

Abstract:In this work a solver for instationary two-phase flows on the basis of the extended Discontinuous Galerkin (extended DG/XDG) method is presented. The XDG method adapts the approximation space conformal to the position of the interface. This allows a sub-cell accurate representation of the incompressible Navier-Stokes equations in their sharp interface formulation. The interface is described as the zero set of a signed-distance level-set function and discretized by a standard DG method. For the interface, resp. level-set, evolution an extension velocity field is used and a two-staged algorithm is presented for its construction on a narrow-band. On the cut-cells a monolithic elliptic extension velocity method is adapted and a fast-marching procedure on the neighboring cells. The spatial discretization is based on a symmetric interior penalty method and for the temporal discretization a moving interface approach is adapted. A cell agglomeration technique is utilized for handling small cut-cells and topology changes during the interface motion. The method is validated against a wide range of typical two-phase surface tension driven flow phenomena including capillary waves, an oscillating droplet and the rising bubble benchmark.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to efficiently and accurately simulate interface motion and its impact on fluid dynamics in incompressible two - phase flow. Specifically, the paper focuses on developing a solver suitable for transient two - phase flow problems within the framework of the extended Discontinuous Galerkin (XDG) method. This method can adapt to changes in interface position, thereby achieving sub - cell - accurate representation of the sharp - interface form of the incompressible Navier - Stokes equations (NSE).

The method proposed in the paper mainly solves the following key problems:

1. **Interface Representation and Evolution**: The interface is described by a signed - distance level - set function defined by the zero set and discretized using the standard DG method. To maintain the signed - distance property and reduce the need for re - initialization, an extended velocity field is introduced, and a two - stage algorithm is proposed to construct this velocity field.

2. **Numerical Stability and Accuracy**: A single - block elliptic extended - velocity method is applied on cut cells, and a fast - marching procedure is adopted on adjacent cells. In addition, to handle small cut cells and topological changes during interface motion, a cell - aggregation technique is used.

3. **Time and Space Discretization**: For time discretization, a moving - interface method is adopted, which allows for high - order accuracy in time without full space - time element discretization. Space discretization is based on the symmetric interior - penalty method.

4. **Verification and Application**: This method has been verified through a series of typical two - phase surface - tension - driven flow phenomena, including capillary wave, oscillating droplet, and rising - bubble benchmark tests.

In summary, this paper aims to provide a new numerical method to simulate complex two - phase flow problems, especially those involving dynamic interface changes, in a more precise and efficient manner.