Researcher profile

Carsten Carstensen

Carsten Carstensen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
2topics
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

7 published item(s)

preprint2022arXiv

A priori and a posteriori error analysis of the lowest-order NCVEM for second-order linear indefinite elliptic problems

The nonconforming virtual element method (NCVEM) for the approximation of the weak solution to a general linear second-order non-selfadjoint indefinite elliptic PDE in a polygonal domain is analyzed under reduced elliptic regularity. The main tool in the a priori error analysis is the connection between the nonconforming virtual element space and the Sobolev space $H^1_0(Ω)$ by a right-inverse $J$ of the interpolation operator $I_h$. The stability of the discrete solution allows for the proof of existence of a unique discrete solution, of a discrete inf-sup estimate and, consequently, for optimal error estimates in the $H^1$ and $L^2$ norms. The explicit residual-based a posteriori error estimate for the NCVEM is reliable and efficient up to the stabilization and oscillation terms. Numerical experiments on different types of polygonal meshes illustrate the robustness of an error estimator and support the improved convergence rate of an adaptive mesh-refinement in comparison to the uniform mesh-refinement.

preprint2022arXiv

Adaptive guaranteed lower eigenvalue bounds with optimal convergence rates

Guaranteed lower Dirichlet eigenvalue bounds (GLB) can be computed for the $m$-th Laplace operator with a recently introduced extra-stabilized nonconforming Crouzeix-Raviart ($m=1$) or Morley ($m=2$) finite element eigensolver. Striking numerical evidence for the superiority of a new adaptive eigensolver motivates the convergence analysis in this paper with a proof of optimal convergence rates of the GLB towards a simple eigenvalue. The proof is based on (a generalization of) known abstract arguments entitled as the axioms of adaptivity. Beyond the known a priori convergence rates, a medius analysis is enfolded in this paper for the proof of best-approximation results. This and subordinated $L^2$ error estimates for locally refined triangulations appear of independent interest. The analysis of optimal convergence rates of an adaptive mesh-refining algorithm is performed in $3$D and highlights a new version of discrete reliability.

preprint2022arXiv

Direct guaranteed lower eigenvalue bounds with optimal a priori convergence rates for the bi-Laplacian

An extra-stabilised Morley finite element method (FEM) directly computes guaranteed lower eigenvalue bounds with optimal a priori convergence rates for the bi-Laplace Dirichlet eigenvalues. The smallness assumption $\min\{λ_h,λ\}h_{\max}^{4}$ $\le 184.9570$ in $2$D (resp. $\le 21.2912$ in $3$D) on the maximal mesh-size $h_{\max}$ makes the computed $k$-th discrete eigenvalue $λ_h\le λ$ a lower eigenvalue bound for the $k$-th Dirichlet eigenvalue $λ$. This holds for multiple and clusters of eigenvalues and serves for the localisation of the bi-Laplacian Dirichlet eigenvalues in particular for coarse meshes. The analysis requires interpolation error estimates for the Morley FEM with explicit constants in any space dimension $n\ge 2$, which are of independent interest. The convergence analysis in $3$D follows the Babuška-Osborn theory and relies on a companion operator for the Morley finite element method. This is based on the Worsey-Farin $3$D version of the Hsieh-Clough-Tocher macro element with a careful selection of center points in a further decomposition of each tetrahedron into $12$ sub-tetrahedra. Numerical experiments in $2$D support the optimal convergence rates of the extra-stabilised Morley FEM and suggest an adaptive algorithm with optimal empirical convergence rates.

preprint2022arXiv

Lowest-order equivalent nonstandard finite element methods for biharmonic plates

The popular (piecewise) quadratic schemes for the biharmonic equation based on triangles are the nonconforming Morley finite element, the discontinuous Galerkin, the $C^0$ interior penalty, and the WOPSIP schemes. Those methods are modified in their right-hand side $F\in H^{-2}(Ω)$ replaced by $F\circ (JI_{\rm M}) $ and then are quasi-optimal in their respective discrete norms. The smoother $JI_{\rm M}$ is defined for a piecewise smooth input function by a (generalized) Morley interpolation $I_{\rm M}$ followed by a companion operator $J$. An abstract framework for the error analysis in the energy, weaker and piecewise Sobolev norms for the schemes is outlined and applied to the biharmonic equation. Three errors are also equivalent in some particular discrete norm from [Carstensen, Gallistl, Nataraj: Comparison results of nonstandard $P_2$ finite element methods for the biharmonic problem, ESAIM Math. Model. Numer. Anal. (2015)] without data oscillations. This paper extends the work [Veeser, Zanotti: Quasi-optimal nonconforming methods for symmetric elliptic problems, SIAM J. Numer. Anal. 56 (2018)] to the discontinuous Galerkin scheme and adds error estimates in weaker and piecewise Sobolev norms.

preprint2022arXiv

Nonconforming virtual elements for the biharmonic equation with Morley degrees of freedom on polygonal meshes

The lowest-order nonconforming virtual element extends the Morley triangular element to polygons for the approximation of the weak solution $u\in V:=H^2_0(Ω)$ to the biharmonic equation. The abstract framework allows (even a mixture of) two examples of the local discrete spaces $V_h(P)$ and a smoother allows rough source terms $F\in V^*=H^{-2}(Ω)$. The a priori and a posteriori error analysis in this paper circumvents any trace of second derivatives by some computable conforming companion operator $J:V_h\to V$ from the nonconforming virtual element space $V_h$. The operator $J$ is a right-inverse of the interpolation operator and leads to optimal error estimates in piecewise Sobolev norms without any additional regularity assumptions on $u\in V$. As a smoother the companion operator modifies the discrete right-hand side and then allows a quasi-best approximation. An explicit residual-based a posteriori error estimator is reliable and efficient up to data oscillations. Numerical examples display the predicted empirical convergence rates for uniform and optimal convergence rates for adaptive mesh-refinement.

preprint2021arXiv

A priori and a posteriori error analysis of the Crouzeix-Raviart and Morley FEM with original and modified righthand sides

This article on nonconforming schemes for $m$ harmonic problems simultaneously treats the Crouzeix-Raviart ($m=1$) and the Morley finite elements ($m=2$) for the original and for modified right-hand side $F$ in the dual space $V^*:=H^{-m}(Ω)$ to the energy space $V:=H^{m}_0(Ω)$. The smoother $J:V_{\rm nc} \to V$ in this paper is a companion operator, that is a linear and bounded right-inverse to the nonconforming interpolation operator $I_{\rm nc} :V\to V_{\rm nc}$, and modifies the discrete right-hand side $F_h:=F\circ J \in V_{\rm nc}^*$. The best-approximation property of the modified scheme from Veeser et al. (2018) is recovered and complemented with an analysis of the convergence rates in weaker Sobolev norms. Examples with oscillating data show that the original method may fail to enjoy the best-approximation property but can also be better than the modified scheme. The a~posteriori analysis of this paper concerns data oscillations of various types in a class of right-hand sides $F\in V^*$. The reliable error estimates involve explicit constants and can be recommended for explicit error control of the piecewise energy norm. The efficiency follows solely up to data oscillations and examples illustrate this can be problematic.

preprint2020arXiv

Morley Finite Element Method for the von Kármán Obstacle Problem

This paper focusses on the von Kármán equations for the moderately large deformation of a very thin plate with the convex obstacle constraint leading to a coupled system of semilinear fourth-order obstacle problem and motivates its nonconforming Morley finite element approximation. The first part establishes the well-posedness of the von Kármán obstacle problem and also discusses the uniqueness of the solution under an a priori and an a posteriori smallness condition on the data. The second part of the article discusses the regularity result of Frehse from 1971 and combines it with the regularity of the solution on a polygonal domain. The third part of the article shows an a priori error estimate for optimal convergence rates for the Morley finite element approximation to the von Kármán obstacle problem for small data. The article concludes with numerical results that illustrates the requirement of smallness assumption on the data for optimal convergence rate.