BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Dirk Praetorius

Dirk Praetorius contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysiscond-mat.mtrl-scimath-phmath.APmath.MPmath.OCphysics.comp-ph

Trust snapshot

Quick read

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

21 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