Equal higher order analysis of an unfitted discontinuous Galerkin method for Stokes flow systems

Abstract

In this work, we analyze an unfitted discontinuous Galerkin discretization for the numerical solution of the Stokes system based on equal higher-order discontinuous velocities and pressures. This approach combines the best from both worlds, firstly the advantages of a piece-wise discontinuous high-order accurate approximation and secondly the advantages of an unfitted to the true geometry grid around possibly complex objects and/or geometrical deformations. Utilizing a fictitious domain framework, the physical domain of interest is embedded in an unfitted background mesh and the geometrically unfitted discretization is built upon symmetric interior penalty discontinuous Galerkin formulation. To enhance stability we enrich the discrete variational formulation with a pressure stabilization term. Moreover, the present contribution adopts high order ghost penalty strategies to address the ill conditioning of the system matrix caused by small truncated elements with respect to the unfitted boundary. Motivated by continuous unfitted FEM [21,74,75] along with other unfitted mesh surveys grounded on discontinuous spaces [10,44,45,73], we use proper velocity and pressure ghost penalties defined on faces of cut cells to establish a robust high-order method, in spite of the cell agglomeration technique usually applied on dG methods. The current presentation should prove valuable in engineering applications where special emphasis is placed on the optimal effective approximation attaining much smaller relative errors in coarser meshes. Inf-sup stability, the optimal order of convergence, and the condition number sensitivity with respect to cut configuration are investigated. Numerical examples verify the theoretical results.

arXiv comment31 pages

References in corpus (13)

• High order unfitted finite element methods on level set domains using isoparametric mappings
• A cut-cell finite volume - finite element coupling approach for fluid-structure interaction in compressible flow
• An Immersed Interface Method for Discrete Surfaces
• A Nitsche cut finite element method for the Oseen problem with general Navier boundary conditions
• A stabilized cut discontinuous Galerkin framework: I. Elliptic boundary value and interface problems
• A Reduced Basis approach for PDEs on parametrized geometries based on the Shifted Boundary Finite Element Method and application to a Stokes Flow
• A Reduced-Order Shifted Boundary Method for Parametrized incompressible Navier-Stokes equations
• CutFEM without cutting the mesh cells: a new way to impose Dirichlet and Neumann boundary conditions on unfitted meshes
• A finite element method by patch reconstruction for the Stokes problem using mixed formulations
• A Reduced Order Approach for the Embedded Shifted Boundary FEM and a Heat Exchange System on Parametrized Geometries
• An unfitted discontinuous Galerkin scheme for conservation laws on evolving surfaces
• Random geometries for optimal control PDE problems based on fictitious domain FEMS and cut elements
• A stabilized DG cut cell method for discretizing the linear transport equation