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

Christoph Lehrenfeld

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

ResearcherUnclaimed source record
math.NANumerical AnalysisComputational Engineering, Finance, and Sciencephysics.comp-phSymbolic Computation

Catalog footprint

What is connected

13works
5topics
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

13 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
preprint2016arXiv

High order exactly divergence-free Hybrid Discontinuous Galerkin Methods for unsteady incompressible flows

In this paper we present an efficient discretization method for the solution of the unsteady incompressible Navier-Stokes equations based on a high order (Hybrid) Discontinuous Galerkin formulation. The crucial component for the efficiency of the discretization method is the disctinction between stiff linear parts and less stiff non-linear parts with respect to their temporal and spatial treatment. Exploiting the flexibility of operator-splitting time integration schemes we combine two spatial discretizations which are tailored for two simpler sub-problems: a corresponding hyperbolic transport problem and an unsteady Stokes problem. For the hyperbolic transport problem a spatial discretization with an Upwind Discontinuous Galerkin method and an explicit treatment in the time integration scheme is rather natural and allows for an efficient implementation. The treatment of the Stokes part involves the solution of linear systems. In this case a discretization with Hybrid Discontinuous Galerkin methods is better suited. We consider such a discretization for the Stokes part with two important features: H(div)-conforming finite elements to garantuee exactly divergence-free velocity solutions and a projection operator which reduces the number of globally coupled unknowns. We present the method, discuss implementational aspects and demonstrate the performance on two and three dimensional benchmark problems.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Removing the stabilization parameter in fitted and unfitted symmetric Nitsche formulations

In many situations with finite element discretizations it is desirable or necessary to impose boundary or interface conditions not as essential conditions -- i.e. through the finite element space -- but through the variational formulation. One popular way to do this is Nitsche's method. In Nitsche's method a stabilization parameter $λ$ has to be chosen "sufficiently large" to provide a stable formulation. Sometimes discretizations based on a Nitsche formulation are criticized because of the need to manually choose this parameter. While in the discontinuous Galerkin community variants of the Nitsche method -- known as "interior penalty" method in the DG context -- are known which do not require such a manually chosen stabilization parameter, this has not been considered for Nitsche formulations in other contexts. We introduce and analyse such a parameter-free variant for two applications of Nitsche's method. First, the classical Nitsche formulation for the imposition of boundary conditions with fitted meshes and secondly, an unfitted finite element discretizations for the imposition of interface conditions is considered. The introduced variants of corresponding Nitsche formulations do not change the sparsity pattern and can easily be implemented into existing finite element codes. The benefit of the new formulations is the removal of the Nitsche stabilization parameter $λ$ while keeping the stability properties of the original formulations for a "sufficiently large" stabilization parameter $λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High order unfitted finite element methods on level set domains using isoparametric mappings

We introduce a new class of unfitted finite element methods with high order accurate numerical integration over curved surfaces and volumes which are only implicitly defined by level set functions. An unfitted finite element method which is suitable for the case of piecewise planar interfaces is combined with a parametric mapping of the underlying mesh resulting in an isoparametric unfitted finite element method. The parametric mapping is constructed in a way such that the quality of the piecewise planar interface reconstruction is significantly improved allowing for high order accurate computations of (unfitted) domain and surface integrals. This approach is new. We present the method, discuss implementational aspects and present numerical examples which demonstrate the quality and potential of this method.

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
preprint2014arXiv

The Nitsche XFEM-DG space-time method and its implementation in three space dimensions

In the recent paper [C. Lehrenfeld, A. Reusken, SIAM J. Num. Anal., 51 (2013)] a new finite element discretization method for a class of two-phase mass transport problems is presented and analyzed. The transport problem describes mass transport in a domain with an evolving interface. Across the evolving interface a jump condition has to be satisfies. The discretization in that paper is a space-time approach which combines a discontinuous Galerkin (DG) technique (in time) with an extended finite element method (XFEM). Using the Nitsche method the jump condition is enforced in a weak sense. While the emphasis in that paper was on the analysis and one dimensional numerical experiments the main contribution of this paper is the discussion of implementation aspects for the spatially three dimensional case. As the space-time interface is typically given only implicitly as the zero-level of a level-set function, we construct a piecewise planar approximation of the space-time interface. This discrete interface is used to divide the space-time domain into its subdomains. An important component within this decomposition is a new method for dividing four-dimensional prisms intersected by a piecewise planar space-time interface into simplices. Such a subdivision algorithm is necessary for numerical integration on the subdomains as well as on the space-time interface. These numerical integrations are needed in the implementation of the Nitsche XFEM-DG method in three space dimensions. Corresponding numerical studies are presented and discussed.

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