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

Dirk Praetorius

Dirk Praetorius 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-sciphysics.comp-phcond-mat.mes-hallmath-phmath.APmath.MPmath.OC

Catalog footprint

What is connected

50works
9topics
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

50 published item(s)

preprint2026arXiv

Multigoal-oriented adaptive finite element method with convergence rates

We formulate and analyze a goal-oriented adaptive finite element method for a symmetric linear elliptic partial differential equation (PDE) that can simultaneously deal with multiple linear goal functionals. In each step of the algorithm, only two linear finite element systems have to be solved. Moreover, all finite element solutions are computed with respect to the same discrete space, while the underlying triangulations are adapted to resolve all inherent singularities simultaneously. Unlike available results for such a setting in the literature, we give a thorough convergence analysis and verify that our algorithm guarantees, in an appropriate sense, even optimal convergence rates. Numerical experiments underline the derived theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Adaptive BEM for elliptic PDE systems, part II: Isogeometric analysis with hierarchical B-splines for weakly-singular integral equations

We formulate and analyze an adaptive algorithm for isogeometric analysis with hierarchical B-splines for weakly-singular boundary integral equations. We prove that the employed weighted-residual error estimator is reliable and converges at optimal algebraic rate. Numerical experiments with isogeometric boundary elements for the 3D Poisson problem confirm the theoretical results, which also cover general elliptic systems like linear elasticity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Adaptive FEM for parameter-errors in elliptic linear-quadratic parameter estimation problems

We consider an elliptic linear-quadratic parameter estimation problem with a finite number of parameters. A novel a priori bound for the parameter error is proved and, based on this bound, an adaptive finite element method driven by an a posteriori error estimator is presented. Unlike prior results in the literature, our estimator, which is composed of standard energy error residual estimators for the state equation and suitable co-state problems, reflects the faster convergence of the parameter error compared to the (co)-state variables. We show optimal convergence rates of our method; in particular and unlike prior works, we prove that the estimator decreases with a rate that is the sum of the best approximation rates of the state and co-state variables. Experiments confirm that our method matches the convergence rate of the parameter error.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The mass-lumped midpoint scheme for computational micromagnetics: Newton linearization and application to magnetic skyrmion dynamics

We discuss a mass-lumped midpoint scheme for the numerical approximation of the Landau-Lifshitz-Gilbert equation, which models the dynamics of the magnetization in ferromagnetic materials. In addition to the classical micromagnetic field contributions, our setting covers the non-standard Dzyaloshinskii-Moriya interaction, which is the essential ingredient for the enucleation and stabilization of magnetic skyrmions. Our analysis also includes the inexact solution of the arising nonlinear systems, for which we discuss both a constraint preserving fixed-point solver from the literature and a novel approach based on the Newton method. We numerically compare the two linearization techniques and show that the Newton solver leads to a considerably lower number of nonlinear iterations. Moreover, in a numerical study on magnetic skyrmions, we demonstrate that, for magnetization dynamics that are very sensitive to energy perturbations, the midpoint scheme, due to its conservation properties, is superior to the dissipative tangent plane schemes from the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Energy contraction and optimal convergence of adaptive iterative linearized finite element methods

We revisit a unified methodology for the iterative solution of nonlinear equations in Hilbert spaces. Our key observation is that the general approach from [Heid & Wihler, Math. Comp. 89 (2020), Calcolo 57 (2020)] satisfies an energy contraction property in the context of (abstract) strongly monotone problems. This property, in turn, is the crucial ingredient in the recent convergence analysis in [Gantner et al., arXiv:2003.10785]. In particular, we deduce that adaptive iterative linearized finite element methods (AILFEMs) lead to full linear convergence with optimal algebraic rates with respect to the degrees of freedom as well as the total computational time.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magnetostatics and micromagnetics with physics informed neural networks

Partial differential equations and variational problems can be solved with physics informed neural networks (PINNs). The unknown field is approximated with neural networks. Minimizing the residuals of the static Maxwell equation at collocation points or the magnetostatic energy, the weights of the neural network are adjusted so that the neural network solution approximates the magnetic vector potential. This way, the magnetic flux density for a given magnetization distribution can be estimated. With the magnetization as an additional unknown, inverse magnetostatic problems can be solved. Augmenting the magnetostatic energy with additional energy terms, micromagnetic problems can be solved. We demonstrate the use of physics informed neural networks for solving magnetostatic problems, computing the magnetization for inverse problems, and calculating the demagnetization curves for two-dimensional geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A short note on plain convergence of adaptive least-squares finite element methods

We show that adaptive least-squares finite element methods driven by the canonical least-squares functional converge under weak conditions on PDE operator, mesh-refinement, and marking strategy. Contrary to prior works, our plain convergence does neither rely on sufficiently fine initial meshes nor on severe restrictions on marking parameters. Finally, we prove that convergence is still valid if a contractive iterative solver is used to obtain the approximate solutions (e.g., the preconditioned conjugate gradient method with optimal preconditioner). The results apply within a fairly abstract framework which covers a variety of model problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive BEM for elliptic PDE systems, Part I: Abstract framework for weakly-singular integral equations

In the present work, we consider weakly-singular integral equations arising from linear second-order strongly-elliptic PDE systems with constant coefficients, including, e.g., linear elasticity. We introduce a general framework for optimal convergence of adaptive Galerkin BEM. We identify certain abstract properties for the underlying meshes, the corresponding mesh-refinement strategy, and the ansatz spaces that guarantee convergence at optimal algebraic rate of an adaptive algorithm driven by the weighted-residual error. These properties are satisfied, e.g., for discontinuous piecewise polynomials on simplicial meshes as well as certain ansatz spaces used for isogeometric analysis. Technical contributions include local inverse estimates for the (non-local) boundary integral operators associated to the PDE system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive IGAFEM with optimal convergence rates: T-splines

We consider an adaptive algorithm for finite element methods for the isogeometric analysis (IGAFEM) of elliptic (possibly non-symmetric) second-order partial differential equations. We employ analysis-suitable T-splines of arbitrary odd degree on T-meshes generated by the refinement strategy of [Morgenstern, Peterseim, Comput. Aided Geom. Design 34 (2015)] in 2D and [Morgenstern, SIAM J. Numer. Anal. 54 (2016)] in 3D. Adaptivity is driven by some weighted residual a posteriori error estimator. We prove linear convergence of the error estimator (which is equivalent to the sum of energy error plus data oscillations) with optimal algebraic rates with respect to the number of elements of the underlying mesh.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence and quasi-optimal cost of adaptive algorithms for nonlinear operators including iterative linearization and algebraic solver

We consider a second-order elliptic boundary value problem with strongly monotone and Lipschitz-continuous nonlinearity. We design and study its adaptive numerical approximation interconnecting a finite element discretization, the Banach-Picard linearization, and a contractive linear algebraic solver. We in particular identify stopping criteria for the algebraic solver that on the one hand do not request an overly tight tolerance but on the other hand are sufficient for the inexact (perturbed) Banach-Picard linearization to remain contractive. Similarly, we identify suitable stopping criteria for the Banach-Picard iteration that leave an amount of linearization error that is not harmful for the residual a-posteriori error estimate to steer reliably the adaptive mesh-refinement. For the resulting algorithm, we prove a contraction of the (doubly) inexact iterates after some amount of steps of mesh-refinement / linerization / algebraic solver, leading to its linear convergence. Moreover, for usual mesh-refinement rules, we also prove that the overall error decays at the optimal rate with respect to the number of elements (degrees of freedom) added with respect to the initial mesh. Finally, we prove that our fully adaptive algorithm drives the overall error down with the same optimal rate also with respect to the overall algorithmic cost expressed as the cumulated sum of the number of mesh elements over all mesh-refinement, linearization, and algebraic solver steps. Numerical experiments support these theoretical findings and illustrate the optimal overall algorithmic cost of the fully adaptive algorithm on several test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dörfler marking with minimal cardinality is a linear complexity problem

Most adaptive finite element strategies employ the Dörfler marking strategy to single out certain elements $\mathcal{M} \subseteq \mathcal{T}$ of a triangulation $\mathcal{T}$ for refinement. In the literature, different algorithms have been proposed to construct $\mathcal{M}$, where usually two goals compete: On the one hand, $\mathcal{M}$ should contain a minimal number of elements. On the other hand, one aims for linear costs with respect to the cardinality of $\mathcal{T}$. Unlike expected in the literature, we formulate and analyze an algorithm, which constructs a minimal set $\mathcal{M}$ at linear costs. Throughout, pseudocodes are given.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Instance-optimal goal-oriented adaptivity

We consider an adaptive finite element method with arbitrary but fixed polynomial degree $p \ge 1$, where adaptivity is driven by an edge-based residual error estimator. Based on the modified maximum criterion from [Diening et al, Found. Comput. Math. 16, 2016], we propose a goal-oriented adaptive algorithm and prove that it is instance optimal. More precisely, the goal-error is bounded by the product of the total errors (being the sum of energy error plus data oscillations) of the primal and the dual problem, and the proposed algorithm is instance optimal with respect to this upper bound. Numerical experiments underline our theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal convergence rates for goal-oriented FEM with quadratic goal functional

We consider a linear elliptic PDE and a quadratic goal functional. The goal-oriented adaptive FEM algorithm (GOAFEM) solves the primal as well as a dual problem, where the goal functional is always linearized around the discrete primal solution at hand. We show that the marking strategy proposed in [Feischl et al, SIAM J. Numer. Anal., 54 (2016)] for a linear goal functional is also optimal for quadratic goal functionals, i.e., GOAFEM leads to linear convergence with optimal convergence rates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Plain convergence of adaptive algorithms without exploiting reliability and efficiency

We consider h-adaptive algorithms in the context of the finite element method (FEM) and the boundary element method (BEM). Under quite general assumptions on the building blocks SOLVE, ESTIMATE, MARK, and REFINE of such algorithms, we prove plain convergence in the sense that the adaptive algorithm drives the underlying a posteriori error estimator to zero. Unlike available results in the literature, our analysis avoids the use of any reliability and efficiency estimate, but only relies on structural properties of the estimator, namely stability on non-refined elements and reduction on refined elements. In particular, the new framework thus covers also problems involving non-local operators like the fractional Laplacian or boundary integral equations, where (discrete) efficiency is (currently) not available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weak-strong uniqueness for the Landau-Lifshitz-Gilbert equation in micromagnetics

We consider the time-dependent Landau-Lifshitz-Gilbert equation. We prove that each weak solution coincides with the (unique) strong solution, as long as the latter exists in time. Unlike available results in the literature, our analysis also includes the physically relevant lower-order terms like Zeeman contribution, anisotropy, stray field, and the Dzyaloshinskii-Moriya interaction (which accounts for the emergence of magnetic Skyrmions). Moreover, our proof gives a template on how to approach weak-strong uniqueness for even more complicated problems, where LLG is (nonlinearly) coupled to other (nonlinear) PDE systems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Adaptive isogeometric boundary element methods with local smoothness control

In the frame of isogeometric analysis, we consider a Galerkin boundary element discretization of the hyper-singular integral equation associated with the 2D Laplacian. We propose and analyze an adaptive algorithm which locally refines the boundary partition and, moreover, steers the smoothness of the NURBS ansatz functions across elements. In particular and unlike prior work, the algorithm can increase and decrease the local smoothness properties and hence exploits the full potential of isogeometric analysis. We prove that the new adaptive strategy leads to linear convergence with optimal algebraic rates. Numerical experiments confirm the theoretical results. A short appendix comments on analogous results for the weakly-singular integral equation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Optimal convergence behavior of adaptive FEM driven by simple (h-h/2)-type error estimators

For some Poisson-type model problem, we prove that adaptive FEM driven by the (h-h/2)-type error estimators from [Ferraz-Leite, Ortner, Praetorius, Numer. Math. 116 (2010)] leads to convergence with optimal algebraic convergence rates. Besides the implementational simplicity, another striking feature of these estimators is that they can provide guaranteed lower bounds for the energy error with known efficiency constant 1.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The saturation assumption yields optimal convergence of two-level adaptive BEM

We consider the convergence of adaptive BEM for weakly-singular and hypersingular integral equations associated with the Laplacian and the Helmholtz operator in 2D and 3D. The local mesh-refinement is driven by some two-level error estimator. We show that the adaptive algorithm drives the underlying error estimates to zero. Moreover, we prove that the saturation assumption already implies linear convergence of the error with optimal algebraic rates.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Computational micromagnetics with Commics

We present our open-source Python module Commics for the study of the magnetization dynamics in ferromagnetic materials via micromagnetic simulations. It implements state-of-the-art unconditionally convergent finite element methods for the numerical integration of the Landau-Lifshitz-Gilbert equation. The implementation is based on the multiphysics finite element software Netgen/NGSolve. The simulation scripts are written in Python, which leads to very readable code and direct access to extensive post-processing. Together with documentation and example scripts, the code is freely available on GitLab.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Iterative solution and preconditioning for the tangent plane scheme in computational micromagnetics

The tangent plane scheme is a time-marching scheme for the numerical solution of the nonlinear parabolic Landau-Lifshitz-Gilbert equation (LLG), which describes the time evolution of ferromagnetic configurations. Exploiting the geometric structure of LLG, the tangent plane scheme requires only the solution of one linear variational form per time-step, which is posed in the discrete tangent space determined by the nodal values of the current magnetization. We develop an effective solution strategy for the arising constrained linear systems, which is based on appropriate Householder reflections. We derive possible preconditioners, which are (essentially) independent of the time-step, and prove that the preconditioned GMRES algorithm leads to linear convergence. Numerical experiments underpin the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A self-consistent spin-diffusion model for micromagnetics

We propose a three-dimensional micromagnetic model that dynamically solves the Landau-Lifshitz-Gilbert equation coupled to the full spin-diffusion equation. In contrast to previous methods, we solve for the magnetization dynamics and the electric potential in a self-consistent fashion. This treatment allows for an accurate description of magnetization dependent resistance changes. Moreover, the presented algorithm describes both spin accumulation due to smooth magnetization transitions and due to material interfaces as in multilayer structures. The model and its finite-element implementation are validated by current driven motion of a magnetic vortex structure. In a second experiment, the resistivity of a magnetic multilayer structure in dependence of the tilting angle of the magnetization in the different layers is investigated. Both examples show good agreement with reference simulations and experiments respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive vertex-centered finite volume methods with convergence rates

We consider the vertex-centered finite volume method with first-order conforming ansatz functions. The adaptive mesh-refinement is driven by the local contributions of the weighted-residual error estimator. We prove that the adaptive algorithm leads to linear convergence with generically optimal algebraic rates for the error estimator and the sum of energy error plus data oscillations. While similar results have been derived for finite element methods and boundary element methods, the present work appears to be the first for adaptive finite volume methods, where the lack of the classical Galerkin orthogonality leads to new challenges.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Basic noise mechanisms of heat-assisted-magnetic recording

Heat-assisted magnetic recording (HAMR) is expected to be a key technology to significantly increase the areal storage density of magnetic recording devices. At high temperatures thermally induced noise becomes a major problem, which must be overcome in order to reliably write magnetic bits with narrow transitions. We propose an elementary model based on the effective recording time window (ERTW) to compute the switching probability of bits during HARM of bit-patterned media. With few assumptions this analytical model allows to gain deeper insights into the basic noise mechanisms like AC and DC noise. Finally, we discuss strategies to reduce noise and to increase the areal storage density of both bit-patterned as well as granular media.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Influence of grain size and exchange interaction on the LLB modeling procedure

Reliably predicting bit-error rates in realistic heat-assisted magnetic recording simulations is a challenging task. Integrating the Landau-Lifshitz-Bloch (LLB) equation can reduce the computational effort to determine the magnetization dynamics in the vicinity of the Curie temperature. If one aims that these dynamics coincide with trajectories calculated from the atomistic Landau-Lifshitz-Gilbert equation, one has to carefully model required temperature dependent material functions such as the zero-field equilibrium magnetization as well as the parallel and normal susceptibilities. We present an extensive study on how these functions depend on grain size and exchange interactions. We show that, if the size or the exchange constant of a reference grain is modified, the material functions can be scaled, according to the changed Curie temperature, yielding negligible errors. This is shown to be valid for volume changes of up to $\pm 40$ % and variations of the exchange constant of up to $\pm10$ %. Besides the temperature dependent material curves, computed switching probabilities also agree well with probabilities separately determined for each system. Our study suggest that there is no need to recalculate the required LLB input functions for each particle. Within the presented limits it is sufficient to scale them to the Curie temperature of the altered system.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A three-dimensional spin-diffusion model for micromagnetics

We implement a finite-element scheme that solves the Landau-Lifshitz-Gilbert equation coupled to a diffusion equation accounting for spin-polarized currents. The latter solves for the spin accumulation not only in magnetic materials but also in nonmagnetic conductors. The presented method incorporates the model by Slonczewski for the description of spin torque in magnetic multilayers as well as the model of Zhang and Li for the description of current driven domain-wall motion. Furthermore it is able to do both resolve the time evolution of the spin accumulation or treat it in an adiabatic fashion by the choice of sufficiently large time steps.

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
preprint2015arXiv

Areal density optimizations for heat-assisted-magnetic recording of high density bit-patterned media

Heat-assisted-magnetic recording (HAMR) is hoped to be the future recording technique for high density storage devices. Nevertheless, there exist several realizations strategies. With a coarse-grained Landau-Lifshitz-Bloch (LLB) model we investigate in detail benefits and disadvantages of continuous and pulsed laser spot recording of shingled and conventional bit-patterned media. Additionally we compare single phase grains and bits having a bilayer structure with graded Curie temperature, consisting of a hard magnetic layer with high $T_{\mathrm{C}}$ and a soft magnetic one with low $T_{\mathrm{C}}$, respectively. To describe the whole write process as realistic as possible a distribution of the grain sizes and Curie temperatures, a displacement jitter of the head and the bit positions are considered. For all these cases we calculate bit error rates of various grain patterns, temperatures and write head positions to optimize the achievable areal storage density. Within our analysis shingled HAMR with a continuous laser pulse moving over the medium reaches the best results, and thus having the highest potential to become the next generation storage device.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Heat assisted magnetic recording of bit patterned media beyond 10 Tb/in$^2$

The limits of the areal storage density as can be achieved with heat assisted magnetic recording (HAMR) are still an open issue. We want to address this central question and present the design of a possible bit patterned medium with an areal storage density above 10 Tb/in$^2$. The model uses hard magnetic recording grains with 5 nm diameter and 10 nm height. It assumes a realistic distribution of the Curie temperature of the underlying material as well as a realistic distribution of the grain size and the grain position. In order to compute the areal density we analyze the detailed switching behavior of a recording bit under different external conditions, which allows to compute the bit error rate of a recording process (shingled and conventional) for different grain spacings and write head positions. Hence, we are able to optimize the areal density of the presented medium.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

L2-orthogonal projections onto finite elements on locally refined meshes are H1-stable

We merge and extend recent results which prove the H1-stability of the L2-orthogonal projection onto standard finite element spaces, provided that the underlying simplicial triangulation is appropriately graded. For lowest-order Courant finite elements S1(T) in Rd with d>=2, we prove that such a grading is always ensured for adaptive meshes generated by newest vertex bisection. For higher-order finite elements Sp(T) with p>=1, we extend existing bounds on the polynomial degree with a computer-assisted proof. We also consider L2-orthogonal projections onto certain subspaces of Sp(T) which incorporate zero Dirichlet boundary conditions resp. an integral mean zero property.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal additive Schwarz methods for the $hp$-BEM: the hypersingular integral operator in 3D on locally refined meshes

We propose and analyze an overlapping Schwarz preconditioner for the $p$ and $hp$ boundary element method for the hypersingular integral equation in 3D. We consider surface triangulations consisting of triangles. The condition number is bounded uniformly in the mesh size $h$ and the polynomial order $p$. The preconditioner handles adaptively refined meshes and is based on a local multilevel preconditioner for the lowest order space. Numerical experiments on different geometries illustrate its robustness.

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

Existence of $\mathcal{H}$-matrix approximants to the inverses of BEM matrices: the simple-layer operator

We consider the question of approximating the inverse $\mathbf W = \mathbf V^{-1}$ of the Galerkin stiffness matrix $\mathbf V$ obtained by discretizing the simple-layer operator $V$ with piecewise constant functions. The block partitioning of $\mathbf W$ is assumed to satisfy any of the standard admissibility criteria that are employed in connection with clustering algorithms to approximate the discrete BEM operator $\mathbf V$. We show that $\mathbf W$ can be approximated by blockwise low-rank matrices such that the error decays exponentially in the block rank employed. Similar exponential approximability results are shown for the Cholesky factorization of $\mathbf V$.

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
preprint2014arXiv

Spin-polarized transport in ferromagnetic multilayers: An unconditionally convergent FEM integrator

We propose and analyze a decoupled time-marching scheme for the coupling of the Landau-Lifshitz-Gilbert equation with a quasilinear diffusion equation for the spin accumulation. This model describes the interplay of magnetization and electron spin accumulation in magnetic and non-magnetic multilayer structures. Despite the strong nonlinearity of the overall PDE system, the proposed integrator requires only the solution of two linear systems per time-step. Unconditional convergence of the integrator towards weak solutions is proved.

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

H-matrix approximability of the inverses of FEM matrices

We study the question of approximability for the inverse of the FEM stiffness matrix for (scalar) second order elliptic boundary value problems by blockwise low rank matrices such as those given by the H-matrix format. We show that exponential convergence in the local block rank can be achieved. We also show that exponentially accurate LU-decompositions in the H-matrix format are possible for the stiffness matrices arising in the FEM. Unlike prior works, our analysis avoids any coupling of the block rank r and the mesh width h and also covers mixed Dirichlet-Neumann-Robin boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On a decoupled linear FEM integrator for Eddy-current-LLG

We propose a numerical integrator for the coupled system of the eddy-current equation with the nonlinear Landau-Lifshitz-Gilbert equation. The considered effective field contains a general field contribution, and we particularly cover exchange, anisotropy, applied field, and magnetic field (stemming from the eddy-current equation). Even though the considered problem is nonlinear, our scheme requires only the solution of two linear systems per time-step. Moreover, our algorithm decouples both equations so that in each time-step, one linear system is solved for the magnetization, and afterwards one linear system is solved for the magnetic field. Unconditional convergence -- at least of a subsequence -- towards a weak solution is proved, and our analysis even provides existence of such weak solutions. Numerical experiments with a micromagnetic benchmark problem underline the performance of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Landau-Lifshitz-Gilbert equation with magnetostriction

To describe and simulate dynamic micromagnetic phenomena, we consider a coupled system of the nonlinear Landau-Lifshitz-Gilbert equation and the conservation of momentum equation. This coupling allows to include magnetostrictive effects into the simulations. Existence of weak solutions has recently been shown in [Carbout et al. 2011]. In our contribution, we give an alternate proof which additionally provides an effective numerical integrator. The latter is based on lowest-order finite elements in space and a linear-implicit Euler time-stepping. Despite the nonlinearity, only two linear systems have to be solved per timestep, and the integrator fully decouples both equations. Finally, we prove unconditional convergence---at least of a subsequence---towards, and hence existence of, a weak solution of the coupled system, as timestep size and spatial mesh-size tend to zero. Numerical experiments conclude the work and shed new light on the existence of blow-up in micromagnetic simulations.

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
preprint2012arXiv

On 2D Newest Vertex Bisection: Optimality of Mesh-Closure and H1-Stability of L2-Projection

Newest vertex bisection (NVB) is a popular local mesh-refinement strategy for regular triangulations which consist of simplices. For the 2D case, we prove that the mesh-closure step of NVB, which preserves regularity of the triangulation, is quasi-optimal and that the corresponding L2-projection onto lowest-order Courant finite elements (P1-FEM) is always H1-stable. Throughout, no additional assumptions on the initial triangulation are imposed. Our analysis thus improves results of Binev, Dahmen & DeVore (Numer. Math. 97, 2004), Carstensen (Constr. Approx. 20, 2004), and Stevenson (Math. Comp. 77, 2008) in the sense that all assumptions of their theorems are removed. Consequently, our results relax the requirements under which adaptive finite element schemes can be mathematically guaranteed to convergence with quasi-optimal rates.

0 citations0 reviewsNo code linkedOpen access