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

Christoph Lehrenfeld

Christoph Lehrenfeld contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phSymbolic Computation

Trust snapshot

Quick read

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

8 published item(s)

preprint2022arXiv

Isoparametric unfitted BDF -- Finite element method for PDEs on evolving domains

We propose a new discretization method for PDEs on moving domains in the setting of unfitted finite element methods, which is provably higher-order accurate in space and time. In the considered setting, the physical domain that evolves essentially arbitrarily through a time-independent computational background domain, is represented by a level set function. For the time discretization, the application of standard time stepping schemes that are based on finite difference approximations of the time derivative is not directly possible, as the degrees of freedom may get active or inactive across such a finite difference stencil in time. In [Lehrenfeld, Olshanskii. An Eulerian finite element method for PDEs in time-dependent domains. ESAIM: M2AN, 53:585--614, 2019] this problem is overcome by extending the discrete solution at every timestep to a sufficiently large neighborhood so that all the degrees of freedom that are relevant at the next time step stay active. But that paper focuses on low-order methods. We advance these results with introducing and analyzing realizable techniques for the extension to higher order. To obtain higher-order convergence in space and time, we combine the BDF time stepping with the isoparametric unfitted FEM. The latter has been used and analyzed for several stationary problems before. However, for moving domains the key ingredient in the method, the transformation of the underlying mesh, becomes time-dependent which gives rise to some technical issues. We treat these with special care, carry out an a priori error analysis and two numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Robust finite element discretizations for a simplified Galbrun's equation

Driven by the challenging task of finding robust discretization methods for Galbrun's equation, we investigate conditions for stability and different aspects of robustness for different finite element schemes on a simplified version of the equations. The considered PDE is a second order indefinite vector-PDE which remains if only the highest order terms of Galbrun's equation are taken into account. A key property for stability is a Helmholtz-type decomposition which results in a strong connection between stable discretizations for Galbrun's equation and Stokes and nearly incompressible linear elasticity problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An unfitted Eulerian finite element method for the time-dependent Stokes problem on moving domains

We analyse a Eulerian Finite Element method, combining a Eulerian time-stepping scheme applied to the time-dependent Stokes equations using the CutFEM approach with inf-sup stable Taylor-Hood elements for the spatial discretisation. This is based on the method introduced by Lehrenfeld \& Olshanskii [ESAIM: M2AN 53(2):585--614] in the context of a scalar convection-diffusion problems on moving domains, and extended to the non-stationary Stokes problem on moving domains by Burman, Frei \& Massing [arXiv:1910.03054 [math.NA]] using stabilised equal-order elements. The analysis includes the geometrical error made by integrating over approximated levelset domains in the discrete CutFEM setting. The method is implemented and the theoretical results are illustrated using numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Locking free and gradient robust H(div)-conforming HDG methods for linear elasticity

Robust discretization methods for (nearly-incompressible) linear elasticity are free of volume-locking and gradient-robust. While volume-locking is a well-known problem that can be dealt with in many different discretization approaches, the concept of gradient-robustness for linear elasticity is new. We discuss both aspects and propose novel Hybrid Discontinuous Galerkin (HDG) methods for linear elasticity. The starting point for these methods is a divergence-conforming discretization. As a consequence of its well-behaved Stokes limit the method is gradient-robust and free of volume-locking. To improve computational efficiency, we additionally consider discretizations with relaxed divergence-conformity and a modification which re-enables gradient-robustness, yielding a robust and quasi-optimal discretization also in the sense of HDG superconvergence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sweeping preconditioners for stratified media in the presence of reflections

In this paper we consider sweeping preconditioners for time harmonic wave propagation in stratified media, especially in the presence of reflections. In the most famous class of sweeping preconditioners Dirichlet-to-Neumann operators for half-space problems are approximated through absorbing boundary conditions. In the presence of reflections absorbing boundary conditions are not accurate resulting in an unsatisfactory performance of these sweeping preconditioners. We explore the potential of using more accurate Dirichlet-to-Neumann operators within the sweep. To this end, we make use of the separability of the equation for the background model. While this improves the accuracy of the Dirichlet-to-Neumann operator, we find both from numerical tests and analytical arguments that it is very sensitive to perturbations in the presence of reflections. This implies that even if accurate approximations to Dirichlet-to-Neumann operators can be devised for a stratified medium, sweeping preconditioners are limited to very small perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Divergence-free tangential finite element methods for incompressible flows on surfaces

In this work we consider the numerical solution of incompressible flows on two-dimensional manifolds. Whereas the compatibility demands of the velocity and the pressure spaces are known from the flat case one further has to deal with the approximation of a velocity field that lies only in the tangential space of the given geometry. Abandoning $H^1$-conformity allows us to construct finite elements which are -- due to an application of the Piola transformation -- exactly tangential. To reintroduce continuity (in a weak sense) we make use of (hybrid) discontinuous Galerkin techniques. To further improve this approach, $H(\operatorname{div}_Γ)$-conforming finite elements can be used to obtain exactly divergence-free velocity solutions. We present several new finite element discretizations. On a number of numerical examples we examine and compare their qualitative properties and accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimal preconditioners for Nitsche-XFEM discretizations of interface problems

In the past decade, a combination of unfitted finite elements (or XFEM) with the Nitsche method has become a popular discretization method for elliptic interface problems. This development started with the introduction and analysis of this Nitsche-XFEM technique in the paper [A. Hansbo, P. Hansbo, Comput. Methods Appl. Mech. Engrg. 191 (2002)]. In general, the resulting linear systems have very large condition numbers, which depend not only on the mesh size $h$, but also on how the interface intersects the mesh. This paper is concerned with the design and analysis of optimal preconditioners for such linear systems. We propose an additive subspace preconditioner which is optimal in the sense that the resulting condition number is independent of the mesh size $h$ and the interface position. We further show that already the simple diagonal scaling of the stifness matrix results in a condition number that is bounded by $ch^{-2}$, with a constant $c$ that does not depend on the location of the interface. Both results are proven for the two-dimensional case. Results of numerical experiments in two and three dimensions are presented, which illustrate the quality of the preconditioner.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Computer Algebra meets Finite Elements: an Efficient Implementation for Maxwell's Equations

We consider the numerical discretization of the time-domain Maxwell's equations with an energy-conserving discontinuous Galerkin finite element formulation. This particular formulation allows for higher order approximations of the electric and magnetic field. Special emphasis is placed on an efficient implementation which is achieved by taking advantage of recurrence properties and the tensor-product structure of the chosen shape functions. These recurrences have been derived symbolically with computer algebra methods reminiscent of the holonomic systems approach.

0 citations0 reviewsNo code linkedOpen access