Researcher profile

Erik Burman

Erik Burman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

17 published item(s)

preprint2026arXiv

Variational data assimilation for the wave equation in heterogeneous media: Numerical investigation of stability

In recent years, several numerical methods for solving the unique continuation problem for the wave equation in a homogeneous medium with given data on the lateral boundary of the space-time cylinder have been proposed. This problem enjoys Lipschitz stability if the geometric control condition is fulfilled, which allows devising optimally convergent numerical methods. In this article, we investigate whether these results carry over to the case in which the medium exhibits a jump discontinuity. Our numerical experiments suggest a positive answer. However, we also observe that the presence of discontinuities in the medium renders the computations far more demanding than in the homogeneous case.

preprint2023arXiv

Data assimilation finite element method for the linearized Navier-Stokes equations with higher order polynomial approximation

In this article, we design and analyze an arbitrary-order stabilized finite element method to approximate the unique continuation problem for laminar steady flow described by the linearized incompressible Navier--Stokes equation. We derive quantitative local error estimates for the velocity, which account for noise level and polynomial degree, using the stability of the continuous problem in the form of a conditional stability estimate. Numerical examples illustrate the performances of the method with respect to the polynomial order and perturbations in the data. We observe that the higher order polynomials may be efficient for ill-posed problems, but are also more sensitive for problems with poor stability due to the ill-conditioning of the system.

preprint2022arXiv

A stabilized finite element method for inverse problems subject to the convection-diffusion equation. II: convection-dominated regime

We consider the numerical approximation of the ill-posed data assimilation problem for stationary convection-diffusion equations and extend our previous analysis in [Numer. Math. 144, 451--477, 2020] to the convection-dominated regime. Slightly adjusting the stabilized finite element method proposed for dominant diffusion, we draw upon a local error analysis to obtain quasi-optimal convergence along the characteristics of the convective field through the data set. The weight function multiplying the discrete solution is taken to be Lipschitz and a corresponding super approximation result (discrete commutator property) is proven. The effect of data perturbations is included in the analysis and we conclude the paper with some numerical experiments.

preprint2022arXiv

The augmented Lagrangian method as a framework for stabilised methods in computational mechanics

In this paper we will review recent advances in the application of the augmented Lagrange multiplier method as a general approach for generating multiplier--free stabilised methods. We first show how the method generates Galerkin/Least Squares type schemes for equality constraints and then how it can be extended to develop new stabilised methods for inequality constraints. Application to several different problems in computational mechanics is given.

preprint2021arXiv

A mechanically consistent model for fluid-structure interactions with contact including seepage

We present a new approach for the mechanically consistent modelling and simulation of fluid-structure interactions with contact. The fundamental idea consists of combining a relaxed contact formulation with the modelling of seepage through a porous layer of co-dimension 1 during contact. For the latter, a Darcy model is considered in a thin porous layer attached to a solid boundary in the limit of infinitesimal thickness. In combination with a relaxation of the contact conditions the computational model is both mechanically consistent and simple to implement. We analyse the approach in detailed numerical studies with both thick- and thin-walled solids, within a fully Eulerian and an immersed approach for the fluid-structure interaction and using fitted and unfitted finite element discretisations.

preprint2021arXiv

CutFEM Based on Extended Finite Element Spaces

We develop a general framework for construction and analysis of discrete extension operators with application to unfitted finite element approximation of partial differential equations. In unfitted methods so called cut elements intersected by the boundary occur and these elements must in general by stabilized in some way. Discrete extension operators provides such a stabilization by modification of the finite element space close to the boundary. More precisely, the finite element space is extended from the stable interior elements over the boundary in a stable way which also guarantees optimal approximation properties. Our framework is applicable to all standard nodal based finite elements of various order and regularity. We develop an abstract theory for elliptic problems and associated parabolic time dependent partial differential equations and derive a priori error estimates. We finally apply this to some examples of partial differential equations of different order including the interface problems, the biharmonic operator and the sixth order triharmonic operator.

preprint2021arXiv

Hybrid coupling of finite element and boundary element methods using Nitsche's method and the Calderon projection

In this paper we discuss a hybridised method for FEM-BEM coupling. The coupling from both sides use a Nitsche type approach to couple to the trace variable. This leads to a formulation that is robust and flexible with respect to approximation spaces and can easily be combined as a building block with other hybridised methods. Energy error norm estimates and the convergence of Jacobi iterations are proved and the performance of the method is illustrated on some computational examples.

preprint2020arXiv

A cut finite element method for a model of pressure in fractured media

We develop a robust cut finite element method for a model of diffusion in fractured media consisting of a bulk domain with embedded cracks. The crack has its own pressure field and can cut through the bulk mesh in a very general fashion. Starting from a common background bulk mesh, that covers the domain, finite element spaces are constructed for the interface and bulk subdomains leading to efficient computations of the coupling terms. The crack pressure field also uses the bulk mesh for its representation. The interface conditions are a generalized form of conditions of Robin type previously considered in the literature which allows the modeling of a range of flow regimes across the fracture. The method is robust in the following way: 1. Stability of the formulation in the full range of parameter choices; and 2. Not sensitive to the location of the interface in the background mesh. We derive an optimal order a priori error estimate and present illustrating numerical examples.

preprint2020arXiv

A pressure-robust discretization of Oseen's equation using stabilization in the vorticity equation

Discretization of Navier-Stokes' equations using pressure-robust finite element methods is considered for the high Reynolds number regime. To counter oscillations due to dominating convection we add a stabilization based on a bulk term in the form of a residual-based least squares stabilization of the vorticity equation supplemented by a penalty term on (certain components of) the gradient jump over the elements faces. Since the stabilization is based on the vorticity equation, it is independent of the pressure gradients, which makes it pressure-robust. Thus, we prove pressure-independent error estimates in the linearized case, known as Oseen's problem. In fact, we prove an $O(h^{k+\frac12})$ error estimate in the $L^2$-norm that is known to be the best that can be expected for this type of problem. Numerical examples are provided that, in addition to confirming the theoretical results, show that the present method compares favorably to the classical residual-based SUPG stabilization.

preprint2020arXiv

Application of a minimal compatible element to incompressible and nearly incompressible continuum mechanics

In this note we will explore some applications of the recently constructed piecewise affine, $H^1$-conforming element that fits in a discrete de Rham complex [Christiansen and Hu, Generalized finite element systems for smooth differential forms and Stokes' problem. Numer. Math. 140 (2018)]. In particular we show how the element leads to locking free methods for incompressible elasticity and viscosity robust methods for the Brinkman model.

preprint2020arXiv

Comparison of Shape Derivatives using CutFEM for Ill-posed Bernoulli Free Boundary Problem

In this paper we discuss a level set approach for the identification of an unknown boundary in a computational domain. The problem takes the form of a Bernoulli problem where only the Dirichlet datum is known on the boundary that is to be identified, but additional information on the Neumann condition is available on the known part of the boundary. The approach uses a classical constrained optimization problem, where a cost functional is minimized with respect to the unknown boundary, the position of which is defined implicitly by a level set function. To solve the optimization problem a steepest descent algorithm using shape derivatives is applied. In each iteration the cut finite element method is used to obtain high accuracy approximations of the pde-model constraint for a given level set configuration without re-meshing. We consider three different shape derivatives. First the classical one, derived using the continuous optimization problem (optimize then discretize). Then the functional is first discretized using the CutFEM method and the shape derivative is evaluated on the finite element functional (discretize then optimize). Finally we consider a third approach, also using a discretized functional. In this case we do not perturb the domain, but consider a so-called boundary value correction method, where a small correction to the boundary position may be included in the weak boundary condition. Using this correction the shape derivative may be obtained by perturbing a distance parameter in the discrete variational formulation. The theoretical discussion is illustrated with a series of numerical examples showing that all three approaches produce similar result on the proposed Bernoulli problem.

preprint2020arXiv

Fully discrete loosely coupled Robin-Robin scheme for incompressible fluid-structure interaction: stability and error analysis

We consider a fully discrete loosely coupled scheme for incompressible fluid-structure interaction based on the time semi-discrete splitting method introduced in {\emph{[Burman, Durst \& Guzmán, arXiv:1911.06760]}}. The splittling method uses a Robin-Robin type coupling that allows for a segregated solution of the solid and the fluid systems, without inner iterations. For the discretisation in space we consider piecewise affine continuous finite elements for all the fields and ensure the inf-sup condition by using a Brezzi-Pitkäranta type pressure stabilization. The interfacial fluid-stresses are evaluated in a variationally consistent fashion, that is shown to admit an equivalent Lagrange multiplier formulation. We prove that the method is unconditionally stable and robust with respect to the amount of added-mass in the system. Furthermore, we provide an error estimate that shows the error in the natural energy norm for the system is $\mathcal O\big(\sqrt{T}(\sqrt{Δt} + h)\big)$ where $T$ is the final time, $Δt$ the time-step length and $h$ the space discretization parameter.

preprint2020arXiv

Low Regularity Estimates for CutFEM Approximations of an Elliptic Problem with Mixed Boundary Conditions

We show error estimates for a cut finite element approximation of a second order elliptic problem with mixed boundary conditions. The error estimates are of low regularity type where we consider the case when the exact solution $u \in H^s$ with $s\in (1,3/2]$. For Nitsche type methods this case requires special handling of the terms involving the normal flux of the exact solution at the the boundary. For Dirichlet boundary conditions the estimates are optimal, whereas in the case of mixed Dirichlet-Neumann boundary conditions they are suboptimal by a logarithmic factor.

preprint2020arXiv

Stability and error analysis of a splitting method using Robin-Robin coupling applied to a fluid-structure interaction problem

We analyze a splitting method for a canonical fluid structure interaction problem. The splittling method uses a Robin-Robin boundary condition, explicit strategy. We prove the method is stable and, furthermore, we provide an error estimate that shows the error at the final time $T$ is $O(\sqrt{TΔt})$ where $Δt$ is the time step.

preprint2019arXiv

A stabilized finite element method for inverse problems subject to the convection-diffusion equation. I: diffusion-dominated regime

The numerical approximation of an inverse problem subject to the convection--diffusion equation when diffusion dominates is studied. We derive Carleman estimates that are on a form suitable for use in numerical analysis and with explicit dependence on the Péclet number. A stabilized finite element method is then proposed and analysed. An upper bound on the condition number is first derived. Combining the stability estimates on the continuous problem with the numerical stability of the method, we then obtain error estimates in local $H^1$- or $L^2$-norms that are optimal with respect to the approximation order, the problem's stability and perturbations in data. The convergence order is the same for both norms, but the $H^1$-estimate requires an additional divergence assumption for the convective field. The theory is illustrated in some computational examples.

preprint2018arXiv

Boundary element methods with weakly imposed boundary conditions

We consider boundary element methods where the Calderón projector is used for the system matrix and boundary conditions are weakly imposed using a particular variational boundary operator designed using techniques from augmented Lagrangian methods. Regardless of the boundary conditions, both the primal trace variable and the flux are approximated. We focus on the imposition of Dirichlet, mixed Dirichlet--Neumann, and Robin conditions. A salient feature of the Robin condition is that the conditioning of the system is robust also for stiff boundary conditions. The theory is illustrated by a series of numerical examples.

preprint2017arXiv

Software frameworks for integral equations in electromagnetic scattering based on Calderón identities

In recent years there have been tremendous advances in the theoretical understanding of boundary integral equations for Maxwell problems. In particular, stable dual pairing of discretisation spaces have been developed that allow robust formulations of the preconditioned electric field, magnetic field and combined field integral equations. Within the BEM++ boundary element library we have developed implementations of these frameworks that allow an intuitive formulation of the typical Maxwell boundary integral formulations within a few lines of code. The basis of these developments is an efficient and robust implementations of Calderón identities together with a product algebra that hides and automates most technicalities involved in assembling Galerkin boundary integral equations. In this paper we demonstrate this framework and use it to derive very simple and robust software formulations of the standard preconditioned electric field, magnetic field and regularised combined field integral equations for Maxwell.