Source author record

Christian Merdon

Christian Merdon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Inf-sup stabilized Scott--Vogelius pairs on general simplicial grids by Raviart--Thomas enrichment

This paper considers the discretization of the Stokes equations with Scott--Vogelius pairs of finite element spaces on arbitrary shape-regular simplicial grids. A novel way of stabilizing these pairs with respect to the discrete inf-sup condition is proposed and analyzed. The key idea consists in enriching the continuous polynomials of order $k$ of the Scott--Vogelius velocity space with appropriately chosen and explicitly given Raviart--Thomas bubbles. This approach is inspired by [Li/Rui, IMA J. Numer. Anal, 2021], where the case $k=1$ was studied. The proposed method is pressure-robust, with optimally converging $\boldsymbol{H}^1$-conforming velocity and a small $\boldsymbol{H}(\mathrm{div})$-conforming correction rendering the full velocity divergence-free. For $k\ge d$, with $d$ being the dimension, the method is parameter-free. Furthermore, it is shown that the additional degrees of freedom for the Raviart--Thomas enrichment and also all non-constant pressure degrees of freedom can be condensated, effectively leading to a pressure-robust, inf-sup stable, optimally convergent $\boldsymbol{P}_k \times P_0$ scheme. Aspects of the implementation are discussed and numerical studies confirm the analytic results.

preprint2022arXiv

Pressure-robustness in the context of optimal control

This paper studies the benefits of pressure-robust discretizations in the scope of optimal control of incompressible flows. Gradient forces that may appear in the data can have a negative impact on the accuracy of state and control and can only be correctly balanced if their $L^2$-orthogonality onto discretely divergence-free test functions is restored. Perfectly orthogonal divergence-free discretizations or divergence-free reconstructions of these test functions do the trick and lead to much better analytic a priori estimates that are also validated in numerical examples.

preprint2020arXiv

A nonconforming pressure-robust finite element method for the Stokes equations on anisotropic meshes

Most classical finite element schemes for the (Navier-)Stokes equations are neither pressure-robust, nor are they inf-sup stable on general anisotropic triangulations. A lack of pressure-robustness may lead to large velocity errors, whenever the Stokes momentum balance is dominated by a strong and complicated pressure gradient. It is a consequence of a method, which does not exactly satisfy the divergence constraint. However, inf-sup stable schemes can often be made pressure-robust just by a recent, modified discretization of the exterior forcing term, using $\mathbf{H}(\operatorname{div})$-conforming velocity reconstruction operators. This approach has so far only been analyzed on shape-regular triangulations. The novelty of the present contribution is that the reconstruction approach for the Crouzeix-Raviart method, which has a stable Fortin operator on arbitrary meshes, is combined with results on the interpolation error on anisotropic elements for reconstruction operators of Raviart-Thomas and Brezzi-Douglas-Marini type, generalizing the method to a large class of anisotropic triangulations. Numerical examples confirm the theoretical results in a 2D and a 3D test case.

preprint2020arXiv

A pressure-robust discretization of Oseen's equation using stabilization in the vorticity equation

Discretization of Navier-Stokes' equations using pressure-robust finite element methods is considered for the high Reynolds number regime. To counter oscillations due to dominating convection we add a stabilization based on a bulk term in the form of a residual-based least squares stabilization of the vorticity equation supplemented by a penalty term on (certain components of) the gradient jump over the elements faces. Since the stabilization is based on the vorticity equation, it is independent of the pressure gradients, which makes it pressure-robust. Thus, we prove pressure-independent error estimates in the linearized case, known as Oseen's problem. In fact, we prove an $O(h^{k+\frac12})$ error estimate in the $L^2$-norm that is known to be the best that can be expected for this type of problem. Numerical examples are provided that, in addition to confirming the theoretical results, show that the present method compares favorably to the classical residual-based SUPG stabilization.

preprint2020arXiv

Divergence-preserving reconstructions on polygons and a really pressure-robust virtual element method for the Stokes problem

Non divergence-free discretisations for the incompressible Stokes problem may suffer from a lack of pressure-robustness characterised by large discretisations errors due to irrotational forces in the momentum balance. This paper argues that also divergence-free virtual element methods (VEM) on polygonal meshes are not really pressure-robust as long as the right-hand side is not discretised in a careful manner. To be able to evaluate the right-hand side for the test functions, some explicit interpolation of the virtual test functions is needed that can be evaluated pointwise everywhere. The standard discretisation via an $L^2$-best approximation does not preserve the divergence and so destroys the orthogonality between divergence-free test functions and possibly eminent gradient forces in the right-hand side. To repair this orthogonality and restore pressure-robustness another divergence-preserving reconstruction is suggested based on Raviart-Thomas approximations on local subtriangulations of the polygons. All findings are proven theoretically and are demonstrated numerically in two dimensions. The construction is also interesting for hybrid high-order methods on polygonal or polyhedral meshes.

preprint2019arXiv

A gradient-robust well-balanced scheme for the compressible isothermal Stokes problem

A novel notion for constructing a well-balanced scheme - a gradient-robust scheme - is introduced and a showcase application for a steady compressible, isothermal Stokes equations is presented. Gradient-robustness means that arbitrary gradient fields in the momentum balance are well-balanced by the discrete pressure gradient - if there is enough mass in the system to compensate the force. The scheme is asymptotic-preserving in the sense that it degenerates for low Mach numbers to a recent inf-sup stable and pressure-robust discretization for the incompressible Stokes equations. The convergence of the coupled FEM-FVM scheme for the nonlinear, isothermal Stokes equations is proved by compactness arguments. Numerical examples illustrate the numerical analysis, and show that the novel approach can lead to a dramatically increased accuracy in nearly-hydrostatic low Mach number flows. Numerical examples also suggest that a straight-forward extension to barotropic situations with nonlinear equations of state is feasible.

preprint2016arXiv

Divergence-free Reconstruction Operators for Pressure-Robust Stokes Discretizations With Continuous Pressure Finite Elements

Classical inf-sup stable mixed finite elements for the incompressible (Navier-)Stokes equations are not pressure-robust, i.e., their velocity errors depend on the continuous pressure. However, a modification only in the right hand side of a Stokes discretization is able to reestablish pressure-robustness, as shown recently for several inf-sup stable Stokes elements with discontinuous discrete pressures. In this contribution, this idea is extended to low and high order Taylor-Hood and mini elements, which have continuous discrete pressures. For the modification of the right hand side a velocity reconstruction operator is constructed that maps discretely divergence-free test functions to exactly divergence-free ones. The reconstruction is based on local $H(\mathrm{div})$-conforming flux equilibration on vertex patches, and fulfills certain orthogonality properties to provide consistency and optimal a-priori error estimates. Numerical examples for the incompressible Stokes and Navier-Stokes equations confirm that the new pressure-robust Taylor-Hood and mini elements converge with optimal order and outperform significantly the classical versions of those elements when the continuous pressure is comparably large.