BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Carsten Carstensen

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

ResearcherUnclaimed source record
math.NANumerical Analysiscond-mat.mtrl-scimath.AP

Catalog footprint

What is connected

13works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Axioms of adaptivity for separate marking

Mixed finite element methods with flux errors in $H(div)$-norms and div-least-squares finite element methods require a separate marking strategy in obligatory adaptive mesh-refining. The refinement indicator $σ^2(\mathcal T,K)=η^2(\mathcal T,K)+μ^2(K)$ of a finite element domain $K$ in an admissible triangulation $\mathcal T$ consists of some residual-based error estimator $η(\mathcal T,K)$ with some reduction property under local mesh-refining and some data approximation error $μ(K)$. Separate marking means either Dörfler marking if $μ^2(\mathcal T) \leq κη^2(\mathcal T)$ or otherwise an optimal data approximation algorithm runs with controlled accuracy as established in [Carstensen, Rabus, Math.Comp. 2011; Rabus, J.Numer.Math. 2015]. The axioms are abstract and sufficient conditions on the estimators $η(\mathcal T,K)$ and data approximation errors $μ(K)$ for optimal asymptotic convergence rates. The enfolded set of axioms simplifies \cite{CFP14} for collective marking, treats separate marking established for the first time in an abstract framework, generalizes [Carstensen, Park, SIAM J.Numer.Anal. 2015] for least-squares schemes, and extends [Carstensen, Rabus, Math.Comp. 2011] to the mixed FEM with flux error control in $H(div)$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Geometry of polycrystals and microstructure

We investigate the geometry of polycrystals, showing that for polycrystals formed of convex grains the interior grains are polyhedral, while for polycrystals with general grain geometry the set of triple points is small. Then we investigate possible martensitic morphologies resulting from intergrain contact. For cubic-to-tetragonal transformations we show that homogeneous zero-energy microstructures matching a pure dilatation on a grain boundary necessarily involve more than four deformation gradients. We discuss the relevance of this result for observations of microstructures involving second and third-order laminates in various materials. Finally we consider the more specialized situation of bicrystals formed from materials having two martensitic energy wells (such as for orthorhombic to monoclinic transformations), but without any restrictions on the possible microstructure, showing how a generalization of the Hadamard jump condition can be applied at the intergrain boundary to show that a pure phase in either grain is impossible at minimum energy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparison results for the Stokes equations

This paper enfolds a medius analysis for the Stokes equations and compares different finite element methods (FEMs). A first result is a best approximation result for a P1 non-conforming FEM. The main comparison result is that the error of the P2-P0-FEM is a lower bound to the error of the Bernardi-Raugel (or reduced P2-P0) FEM, which is a lower bound to the error of the P1 non-conforming FEM, and this is a lower bound to the error of the MINI-FEM. The paper discusses the converse direction, as well as other methods such as the discontinuous Galerkin and pseudostress FEMs. Furthermore this paper provides counterexamples for equivalent convergence when different pressure approximations are considered. The mathematical arguments are various conforming companions as well as the discrete inf-sup condition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Error analysis of nonconforming and mixed FEMs for second-order linear non-selfadjoint and indefinite elliptic problems

The state-of-the art proof of a global inf-sup condition on mixed finite element schemes does not allow for an analysis of truly indefinite, second-order linear elliptic PDEs. This paper, therefore, first analyses a nonconforming finite element discretization which converges owing to some a priori $L^2$ error estimates even for reduced regularity on non-convex polygonal domains. An equivalence result of that nonconforming finite element scheme to the mixed finite element method (MFEM) leads to the well-posedness of the discrete solution and to a priori error estimates for the MFEM. The explicit residual-based a posteriori error analysis allows some reliable and efficient error control and motivates some adaptive discretization which improves the empirical convergence rates in three computational benchmarks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Axioms of Adaptivity

This paper aims first at a simultaneous axiomatic presentation of the proof of optimal convergence rates for adaptive finite element methods and second at some refinements of particular questions like the avoidance of (discrete) lower bounds, inexact solvers, inhomogeneous boundary data, or the use of equivalent error estimators. Solely four axioms guarantee the optimality in terms of the error estimators. Compared to the state of the art in the temporary literature, the improvements of this article can be summarized as follows: First, a general framework is presented which covers the existing literature on optimality of adaptive schemes. The abstract analysis covers linear as well as nonlinear problems and is independent of the underlying finite element or boundary element method. Second, efficiency of the error estimator is neither needed to prove convergence nor quasi-optimal convergence behavior of the error estimator. In this paper, efficiency exclusively characterizes the approximation classes involved in terms of the best-approximation error and data resolution and so the upper bound on the optimal marking parameters does not depend on the efficiency constant. Third, some general quasi-Galerkin orthogonality is not only sufficient, but also necessary for the $R$-linear convergence of the error estimator, which is a fundamental ingredient in the current quasi-optimality analysis due to Stevenson 2007. Finally, the general analysis allows for equivalent error estimators and inexact solvers as well as different non-homogeneous and mixed boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Uniform convergence and a posteriori error estimation for assumed stress hybrid finite element methods

Assumed stress hybrid methods are known to improve the performance of standard displacement-based finite elements and are widely used in computational mechanics. The methods are based on the Hellinger-Reissner variational principle for the displacement and stress variables. This work analyzes two existing 4-node hybrid stress quadrilateral elements due to Pian and Sumihara [Int. J. Numer. Meth. Engng, 1984] and due to Xie and Zhou [Int. J. Numer. Meth. Engng, 2004], which behave robustly in numerical benchmark tests. For the finite elements, the isoparametric bilinear interpolation is used for the displacement approximation, while different piecewise-independent 5-parameter modes are employed for the stress approximation. We show that the two schemes are free from Poisson-locking, in the sense that the error bound in the a priori estimate is independent of the relevant Lame constant $λ$. We also establish the equivalence of the methods to two assumed enhanced strain schemes. Finally, we derive reliable and efficient residual-based a posteriori error estimators for the stress in $L^{2}$-norm and the displacement in $H^{1}$-norm, and verify the theoretical results by some numerical experiments.

0 citations0 reviewsNo code linkedOpen access