Source author record

Carl-Martin Pfeiler

Carl-Martin Pfeiler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis physics.comp-ph

Catalog footprint

What is connected

6works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

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.

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.

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.

preprint2016arXiv

Higher order unfitted FEM for Stokes interface problems

We consider the discretization of a stationary Stokes interface problem in a velocity-pressure formulation. The interface is described implicitly as the zero level of a scalar function as it is common in level set based methods. Hence, the interface is not aligned with the mesh. An unfitted finite element discretization based on a Taylor-Hood velocity-pressure pair and an XFEM (or CutFEM) modification is used for the approximation of the solution. This allows for the accurate approximation of solutions which have strong or weak discontinuities across interfaces which are not aligned with the mesh. To arrive at a consistent, stable and accurate formulation we require several additional techniques. First, a Nitsche-type formulation is used to implement interface conditions in a weak sense. Secondly, we use the ghost penalty stabilization to obtain an inf-sup stable variational formulation. Finally, for the highly accurate approximation of the implicitly described geometry, we use a combination of a piecewise linear interface reconstruction and a parametric mapping of the underlying mesh. We introduce the method and discuss results of numerical examples.

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.