Dual weighted residual based error control for nonstationary convection-dominated equations: potential or ballast?

Even though substantial progress has been made in the numerical approximation of convection-dominated problems, its major challenges remain in the scope of current research. In particular, parameter robust a posteriori error estimates for quantities of physical interest and adaptive mesh refinement strategies with proved convergence are still missing. Here, we study numerically the potential of the Dual Weighted Residual (DWR) approach applied to stabilized finite element methods to further enhance the quality of approximations. The impact of a strict application of the DWR methodology is particularly focused rather than the reduction of computational costs for solving the dual problem by interpolation or localization.Comment: arXiv admin note: text overlap with arXiv:1803.1064

CutFEM and ghost stabilization techniques for higher order space-time discretizations of the Navier-Stokes equations

We propose and analyze computationally a new fictitious domain method, based on higher order space-time finite element discretizations, for the simulation of the nonstationary, incompressible Navier-Stokes equations on evolving domains. The physical domain is embedded into a fixed computational mesh such that arbitrary intersections of the moving domain's boundaries with the background mesh occur. The potential of such cut finite element techniques for higher order space-time finite element methods has rarely been studied in the literature so far and deserves further elucidation. The key ingredients of the approach are the weak formulation of Dirichlet boundary conditions by Nitsche's method, the flexible and efficient integration over all types of intersections of cells by moving boundaries and the spatial extension of the discrete physical quantities to the entire computational background mesh including fictitious (ghost) subdomains of fluid flow. Thereby, an expensive remeshing and adaptation of the sparse matrix data structure are avoided and the computations are accelerated. To prevent spurious oscillations caused by irregular intersections of mesh cells, a penalization, defining also implicitly the extension to ghost domains, is added. These techniques are embedded in an arbitrary order, discontinuous Galerkin discretization of the time variable and an inf-sup stable discretization of the spatial variables. The parallel implementation of the matrix assembly is described. The optimal order convergence properties of the algorithm are illustrated in a numerical experiment for an evolving domain. The well-known 2d benchmark of flow around a cylinder as well as flow around moving obstacles with arising cut cells and fictitious domains are considered further

Structure preserving discontinuous Galerkin approximation of a hyperbolic-parabolic system

We study the numerical approximation of a coupled hyperbolic-parabolic system by a family of discontinuous Galerkin space-time finite element methods. The model is rewritten as a first-order evolutionary problem that is treated by the unified abstract solution theory of R.\ Picard. To preserve the mathematical structure of the evolutionary equation on the fully discrete level, suitable generalizations of the distribution gradient and divergence operators on broken polynomial spaces on which the discontinuous Galerkin approach is built on are defined. Well-posedness of the fully discrete problem and error estimates for the discontinuous Galerkin approximation in space and time are proved

A geometric multigrid method for space-time finite element discretizations of the Navier-Stokes equations and its application to 3d flow simulation

We present a parallelized geometric multigrid (GMG) method, based on the cell-based Vanka smoother, for higher order space-time finite element methods (STFEM) to the incompressible Navier--Stokes equations. The STFEM is implemented as a time marching scheme. The GMG solver is applied as a preconditioner for GMRES iterations. Its performance properties are demonstrated for 2d and 3d benchmarks of flow around a cylinder. The key ingredients of the GMG approach are the construction of the local Vanka smoother over all degrees of freedom in time of the respective subinterval and its efficient application. For this, data structures that store pre-computed cell inverses of the Jacobian for all hierarchical levels and require only a reasonable amount of memory overhead are generated. The GMG method is built for the \emph{deal.II} finite element library. The concepts are flexible and can be transferred to similar software platforms.Comment: Key updates of this revision: - Added Subsection 5.2 "Parallel scaling", in which a strong scaling benchmark is performed - Added Subsection 5.3 "Parameter robustness regarding v", where the robustness of the proposed numerical scheme, regarding changes in the viscosity, is computationally analyze