Researcher profile

Christian Merdon

Christian Merdon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 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.