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

Michael Feischl

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

ResearcherUnclaimed source record
math.NANumerical Analysismath.AP

Catalog footprint

What is connected

20works
3topics
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

20 published item(s)

preprint2022arXiv

Inf-sup stability implies quasi-orthogonality

We prove new optimality results for adaptive mesh refinement algorithms for non-symmetric, indefinite, and time-dependent problems by proposing a generalization of quasi-orthogonality which follows directly from the inf-sup stability of the underlying problem. This completely removes a central technical difficulty in modern proofs of optimal convergence of adaptive mesh refinement algorithms and leads to simple optimality proofs for the Taylor-Hood discretization of the stationary Stokes problem, a finite-element/boundary-element discretization of an unbounded transmission problem, and an adaptive time-stepping scheme for parabolic equations. The main technical tool are new stability bounds for the LU-factorization of matrices together with a recently established connection between quasi-orthogonality and matrix factorization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Higher-order linearly implicit full discretization of the Landau--Lifshitz--Gilbert equation

For the Landau--Lifshitz--Gilbert (LLG) equation of micromagnetics we study linearly implicit backward difference formula (BDF) time discretizations up to order $5$ combined with higher-order non-conforming finite element space discretizations, which are based on the weak formulation due to Alouges but use approximate tangent spaces that are defined by $L^2$-averaged instead of nodal orthogonality constraints. We prove stability and optimal-order error bounds in the situation of a sufficiently regular solution. For the BDF methods of orders $3$ to~$5$, this requires %a mild time step restriction $τ\leqslant ch$ and that the damping parameter in the LLG equations be above a positive threshold; this condition is not needed for the A-stable methods of orders $1$ and $2$, for which furthermore a discrete energy inequality irrespective of solution regularity is proved.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recurrent Neural Networks as Optimal Mesh Refinement Strategies

We show that an optimal finite element mesh refinement algorithm for a prototypical elliptic PDE can be learned by a recurrent neural network with a fixed number of trainable parameters independent of the desired accuracy and the input size, i.e., number of elements of the mesh. Moreover, for a general class of PDEs with solutions which are well-approximated by deep neural networks, we show that an optimal mesh refinement strategy can be learned by recurrent neural networks. This includes problems for which no optimal adaptive strategy is known yet.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sparse Compression of Expected Solution Operators

We show that the expected solution operator of prototypical linear elliptic partial differential operators with random coefficients is well approximated by a computable sparse matrix. This result is based on a random localized orthogonal multiresolution decomposition of the solution space that allows both the sparse approximate inversion of the random operator represented in this basis as well as its stochastic averaging. The approximate expected solution operator can be interpreted in terms of classical Haar wavelets. When combined with a suitable sampling approach for the expectation, this construction leads to an efficient method for computing a sparse representation of the expected solution operator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Eddy Current-LLG Equations-Part I: FEM-BEM Coupling

We analyse a numerical method for the coupled system of the eddy current equations in $\mathbb{R}^3$ with the Landau-Lifshitz-Gilbert equation in a bounded domain. The unbounded domain is discretised by means of finite-element/boundary-element coupling. Even though the considered problem is strongly nonlinear, the numerical approach is constructed such that only two linear systems per time step have to be solved. In this first part of the paper, we prove unconditional weak convergence (of a subsequence) of the finite-element solutions towards a weak solution. A priori error estimates will be presented in the second part.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Adaptive 2D IGA boundary element methods

We derive and discuss a posteriori error estimators for Galerkin and collocation IGA boundary element methods for weakly-singular integral equations of the first-kind in 2D. While recent own work considered the Faermann residual error estimator for Galerkin IGA boundary element methods, the present work focuses more on collocation and weighted- residual error estimators, which provide reliable upper bounds for the energy error. Our analysis allows piecewise smooth parametrizations of the boundary, local mesh-refinement, and related standard piecewise polynomials as well as NURBS. We formulate an adaptive algorithm which steers the local mesh-refinement and the multiplicity of the knots. Numerical experiments show that the proposed adaptive strategy leads to optimal convergence, and related IGA boundary element methods are superior to standard boundary element methods with piecewise polynomials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An abstract analysis of optimal goal-oriented adaptivity

We provide an abstract framework for optimal goal-oriented adaptivity for finite element methods and boundary element methods in the spirit of [Carstensen et al., Comput. Math. Appl. 67 (2014)]. We prove that this framework covers standard discretizations of general second-order linear elliptic PDEs and hence generalizes available results [Mommer & Stevenson, SIAM J. Numer. Anal. 47 (2009); Becker et al., SIAM J. Numer. Anal. 49 (2011)] beyond the Poisson equation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence of adaptive BEM and adaptive FEM-BEM coupling for estimators without h-weighting factor

We analyze adaptive mesh-refining algorithms in the frame of boundary element methods (BEM) and the coupling of finite elements and boundary elements (FEM-BEM). Adaptivity is driven by the two-level error estimator proposed by Ernst P. Stephan, Norbert Heuer, and coworkers in the frame of BEM and FEM-BEM or by the residual error estimator introduced by Birgit Faermann for BEM for weakly-singular integral equations. We prove that in either case the usual adaptive algorithm drives the associated error estimator to zero. Emphasis is put on the fact that the error estimators considered are {not even globally equivalent to weighted-residual error estimators for which recently convergence with quasi-optimal algebraic rates has been derived.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multiscale modeling in micromagnetics: existence of solutions and numerical integration

Various applications ranging from spintronic devices, giant magnetoresistance sensors, and magnetic storage devices, include magnetic parts on very different length scales. Since the consideration of the Landau-Lifshitz-Gilbert equation (LLG) constrains the maximum element size to the exchange length within the media, it is numerically not attractive to simulate macroscopic parts with this approach. On the other hand, the magnetostatic Maxwell equations do not constrain the element size, but cannot describe the short-range exchange interaction accurately. A combination of both methods allows to describe magnetic domains within the micromagnetic regime by use of LLG and also considers the macroscopic parts by a non-linear material law using the Maxwell equations. In our work, we prove that under certain assumptions on the non-linear material law, this multiscale version of LLG admits weak solutions. Our proof is constructive in the sense that we provide a linear-implicit numerical integrator for the multiscale model such that the numerically computable finite element solutions admit weak $H^1$-convergence (at least for a subsequence) towards a weak solution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reliable and efficient a posteriori error estimation for adaptive IGA boundary element methods for weakly-singular integral equations

We consider the Galerkin boundary element method (BEM) for weakly-singular integral equations of the first-kind in 2D. We analyze some residual-type a posteriori error estimator which provides a lower as well as an upper bound for the unknown Galerkin BEM error. The required assumptions are weak and allow for piecewise smooth parametrizations of the boundary, local mesh-refinement, and related standard piecewise polynomials as well as NURBS. In particular, our analysis gives a first contribution to adaptive BEM in the frame of isogeometric analysis (IGABEM), for which we formulate an adaptive algorithm which steers the local mesh-refinement and the multiplicity of the knots. Numerical experiments underline the theoretical findings and show that the proposed adaptive strategy leads to optimal convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adaptive FEM with optimal convergence rates for a certain class of non-symmetric and possibly non-linear problems

We analyze adaptive mesh-refining algorithms for conforming finite element discretizations of certain non-linear second-order partial differential equations. We allow continuous polynomials of arbitrary, but fixed polynomial order. The adaptivity is driven by the residual error estimator. We prove convergence even with optimal algebraic convergence rates. In particular, our analysis covers general linear second-order elliptic operators. Unlike prior works for linear non-symmetric operators, our analysis avoids the interior node property for the refinement, and the differential operator has to satisfy a Gårding inequality only. If the differential operator is uniformly elliptic, no additional assumption on the initial mesh is posed.

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
preprint2013arXiv

Convergence and Quasi-Optimality of Adaptive FEM with Inhomogeneous Dirichlet Data

We consider the solution of a second order elliptic PDE with inhomogeneous Dirichlet data by means of adaptive lowest-order FEM. As is usually done in practice, the given Dirichlet data are discretized by nodal interpolation. As model example serves the Poisson equation with mixed Dirichlet-Neumann boundary conditions. For error estimation, we use an edge-based residual error estimator which replaces the volume residual contributions by edge oscillations. For 2D, we prove convergence of the adaptive algorithm even with quasi-optimal convergence rate. For 2D and 3D, we show convergence if the nodal interpolation operator is replaced by the L^2-projection or the Scott-Zhang quasi-interpolation operator. As a byproduct of the proof, we show that the Scott-Zhang operator converges pointwise to a limiting operator as the mesh is locally refined. This property might be of independent interest besides the current application. Finally, numerical experiments conclude the work.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Each H^{1/2}-stable projection yields convergence and quasi-optimality of adaptive FEM with inhomogeneous Dirichlet data in R^d

We consider the solution of second order elliptic PDEs in $\R^d$ with inhomogeneous Dirichlet data by means of an $h$-adaptive FEM with fixed polynomial order $p\in\N$. As model example serves the Poisson equation with mixed Dirichlet-Neumann boundary conditions, where the inhomogeneous Dirichlet data are discretized by use of an $H^{1/2}$-stable projection, for instance, the $L^2$-projection for $p=1$ or the Scott-Zhang projection for general $p\ge1$. For error estimation, we use a residual error estimator which includes the Dirichlet data oscillations. We prove that each $H^{1/2}$-stable projection yields convergence of the adaptive algorithm even with quasi-optimal convergence rate. Numerical experiments with the $L^2$- and Scott-Zhang projection conclude the work.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

ZZ-type aposteriori error estimators for adaptive boundary element methods on a curve

In the context of the adaptive finite element method (FEM), ZZ-error estimators named after Zienkiewicz and Zhu are mathematically well-established and widely used in practice. In this work, we propose and analyze ZZ-type error estimators for the adaptive boundary element method (BEM). We consider weakly-singular and hyper-singular integral equations and prove, in particular, convergence of the related adaptive mesh-refining algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Classical FEM-BEM coupling methods: nonlinearities, well-posedness, and adaptivity

We consider a (possibly) nonlinear interface problem in 2D and 3D, which is solved by use of various adaptive FEM-BEM coupling strategies, namely the Johnson-Nédélec coupling, the Bielak-MacCamy coupling, and Costabel's symmetric coupling. We provide a framework to prove that the continuous as well as the discrete Galerkin solutions of these coupling methods additionally solve an appropriate operator equation with a Lipschitz continuous and strongly monotone operator. Therefore, the coupling formulations are well-defined, and the Galerkin solutions are quasi-optimal in the sense of a Céa-type lemma. For the respective Galerkin discretizations with lowest-order polynomials, we provide reliable residual-based error estimators. Together with an estimator reduction property, we prove convergence of the adaptive FEM-BEM coupling methods. A key point for the proof of the estimator reduction are novel inverse-type estimates for the involved boundary integral operators which are advertized. Numerical experiments conclude the work and compare performance and effectivity of the three adaptive coupling procedures in the presence of generic singularities.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence of adaptive FEM for some elliptic obstacle problem with inhomogeneous Dirichlet data

In this work, we show the convergence of adaptive lowest-order FEM (AFEM) for an elliptic obstacle problem with non-homogeneous Dirichlet data, where the obstacle χ is restricted only by χ in H^2(Ω). The adaptive loop is steered by some residual based error estimator introduced in Braess, Carstensen & Hoppe (2007) that is extended to control oscillations of the Dirichlet data, as well. In the spirit of Cascon et al. (2008), we show that a weighted sum of energy error, estimator, and Dirichlet oscillations satisfies a contraction property up to certain vanishing energy contributions. This result extends the analysis of Braess, Carstensen & Hoppe (2007) and Page & Praetorius (2010) to the case of non-homogeneous Dirichlet data as well as certain non-affine obstacles and introduces some energy estimates to overcome the lack of nestedness of the discrete spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Inverse estimates for elliptic boundary integral operators and their application to the adaptive coupling of FEM and BEM

We prove inverse-type estimates for the four classical boundary integral operators associated with the Laplace operator. These estimates are used to show convergence of an h-adaptive algorithm for the coupling of a finite element method with a boundary element method which is driven by a weighted residual error estimator.

0 citations0 reviewsNo code linkedOpen access