An unfitted interior penalty discontinuous Galerkin method for incompressible Navier–Stokes two‐phase flow

Abstract

SUMMARY A discontinuous Galerkin method for the solution of the immiscible and incompressible two‐phase flow problem based on the nonsymmetric interior penalty method is presented. Therefore, the incompressible Navier–Stokes equation is solved for a domain decomposed into two subdomains with different values of viscosity and density as well as a singular surface tension force. On the basis of a piecewise linear approximation of the interface, meshes for both phases are cut out of a structured mesh. The discontinuous finite elements are defined on the resulting Cartesian cut‐cell mesh and may therefore approximate the discontinuities of the pressure and the velocity derivatives across the interface with high accuracy. As the mesh resolves the interface, regularization of the density and viscosity jumps across the interface is not required. This preserves the local conservation property of the velocity field even in the vicinity of the interface and constitutes a significant advantage compared with standard methods that require regularization of these discontinuities and cannot represent the jumps and kinks in pressure and velocity. A powerful subtessellation algorithm is incorporated to allow the usage of standard time integrators (such as Crank–Nicholson) on the time‐dependent mesh. The presented discretization is applicable to both the two‐dimensional and three‐dimensional cases. The performance of our approach is demonstrated by application to a two‐dimensional benchmark problem, allowing for a thorough comparison with other numerical methods. Copyright © 2012 John Wiley & Sons, Ltd.

The non-symmetric Nitsche method for the parameter-free imposition of weak boundary and coupling conditions in immersed finite elements

“…On the other hand, the non-symmetric Nitsche method leads to unsymmetric system matrices, convergence is currently only proved for basis functions of degree p ≥ 2, and convergence rates of the L 2 error are suboptimal [3,4,46,58,59]. A stabilization term can mitigate the reduced L 2 accuracy [5,34,43,75], but again leads to the question of estimating appropriate parameters.…”

A parameter-free variational coupling approach for trimmed isogeometric thin shells

“…In line with [77], we could remove all oscillations by re-introducing moderate stabilization. Although we did not see similar oscillations in any other example we computed, the potential of increasing robustness by moderate stabilization should be further examined in the future.…”

“…In the context of the non-symmetric interior penalty discontinous Galerkin method [55,57], penalty-free nonsymmetric Nitsche formulations have been reported to lead to oscillations near interfaces [77]. One way to effectively reduce these oscillations is to introduce a stabilization term, that has the same form as in the symmetric Nitsche method.…”

“…The two patches are coupled weakly with the non-symmetric Nitsche method, where we add the stabilization terms (76). Figure 32 plots the bending moment computed with the non-symmetric Nitsche method at different values C=C S =C N , which determines the level of stabilization via relations (77). We observe that the parameter-free variant leads to local oscillations at the coupling interface, which can be mitigated by increasing the level of stabilization.…”

The Unfitted Discontinuous Galerkin Method for Solving the EEG Forward Problem

“…The method was first introduced in the context of micro-scale simulations in porous media [29] and has since then found different applications [30]- [32]. It combines the advantages of the DG-FEM approach on a structured hexahedral mesh with implicit surface representations given by level set functions and thereby the avoidance of skull leakages and inappropriate modeling of smooth tissue surfaces.…”