Source author record

Martin Vohralík

Martin Vohralík 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

4works
2topics
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

4 published item(s)

preprint2022arXiv

$p$-robust equilibrated flux reconstruction in ${\boldsymbol H}(\mathrm{curl})$ based on local minimizations. Application to a posteriori analysis of the curl-curl problem

We present a local construction of H(curl)-conforming piecewise polynomials satisfying a prescribed curl constraint. We start from a piecewise polynomial not contained in the H(curl) space but satisfying a suitable orthogonality property. The procedure employs minimizations in vertex patches and the outcome is, up to a generic constant independent of the underlying polynomial degree, as accurate as the best-approximations over the entire local versions of H(curl). This allows to design guaranteed, fully computable, constant-free, and polynomial-degree-robust a posteriori error estimates of Prager-Synge type for Nédélec finite element approximations of the curl-curl problem. A divergence-free decomposition of a divergence-free H(div)-conforming piecewise polynomial, relying on over-constrained minimizations in Raviart-Thomas spaces, is the key ingredient. Numerical results illustrate the theoretical developments.

preprint2020arXiv

Equivalence of local-and global-best approximations, a simple stable local commuting projector, and optimal $hp$ approximation estimates in $H(\mathrm{div})$

Given an arbitrary function in H(div), we show that the error attained by the global-best approximation by H(div)-conforming piecewise polynomial Raviart-Thomas-Nédélec elements under additional constraints on the divergence and normal flux on the boundary, is, up to a generic constant, equivalent to the sum of independent local-best approximation errors over individual mesh elements, without constraints on the divergence or normal fluxes. The generic constant only depends on the shape-regularity of the underlying simplicial mesh, the space dimension, and the polynomial degree of the approximations. The analysis also gives rise to a stable, local, commuting projector in H(div), delivering an approximation error that is equivalent to the local-best approximation. We next present a variant of the equivalence result, where robustness of the constant with respect to the polynomial degree is attained for unbalanced approximations. These two results together further enable us to derive rates of convergence of global-best approximations that are fully optimal in both the mesh size h and the polynomial degree p, for vector fields that only feature elementwise the minimal necessary Sobolev regularity. We finally show how to apply our findings to derive optimal a priori hp-error estimates for mixed and least-squares finite element methods applied to a model diffusion problem.

preprint2020arXiv

Guaranteed a posteriori bounds for eigenvalues and eigenvectors: multiplicities and clusters

This paper presents a posteriori error estimates for conforming numerical approximations of eigenvalue clusters of second-order self-adjoint elliptic linear operators with compact resolvent. Given a cluster of eigenvalues, we estimate the error in the sum of the eigenvalues, as well as the error in the eigenvectors represented through the density matrix, i.e., the orthogonal projector on the associated eigenspace. This allows us to deal with degenerate (multiple) eigenvalues within the framework. All the bounds are valid under the only assumption that the cluster is separated from the surrounding smaller and larger eigenvalues; we show how this assumption can be numerically checked. Our bounds are guaranteed and converge with the same speed as the exact errors. They can be turned into fully computable bounds as soon as an estimate on the dual norm of the residual is available, which is presented in two particular cases: the Laplace eigenvalue problem discretized with conforming finite elements, and a Schr{ö}dinger operator with periodic boundary conditions of the form $--$Δ$ + V$ discretized with planewaves. For these two cases, numerical illustrations are provided on a set of test problems.

preprint2020arXiv

Polynomial-degree-robust H(curl)-stability of discrete minimization in a tetrahedron

We prove that the minimizer in the Nédélec polynomial space of some degree p > 0 of a discrete minimization problem performs as well as the continuous minimizer in H(curl), up to a constant that is independent of the polynomial degree p. The minimization problems are posed for fields defined on a single non-degenerate tetrahedron in R^3 with polynomial constraints enforced on the curl of the field and its tangential trace on some faces of the tetrahedron. This result builds upon [L. Demkowicz, J. Gopalakrishnan, J. Schöberl SIAM J. Numer. Anal. 47 (2009), 3293--3324] and [M. Costabel, A. McIntosh, Math. Z. 265 (2010), 297--320] and is a fundamental ingredient to build polynomial-degree-robust a posteriori error estimators when approximating the Maxwell equations in several regimes leading to a curl-curl problem.