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

Erik Burman

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

ResearcherUnclaimed source record
math.NANumerical Analysismath.AP

Catalog footprint

What is connected

46works
3topics
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

46 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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Augmented Lagrangian finite element methods for contact problems

We propose two different Lagrange multiplier methods for contact problems derived from the augmented Lagrangian variational formulation. Both the obstacle problem, where a constraint on the solution is imposed in the bulk domain and the Signorini problem, where a lateral contact condition is imposed are considered. We consider both continuous and discontinuous approximation spaces for the Lagrange multiplier. In the latter case the method is unstable and a penalty on the jump of the multiplier must be applied for stability. We prove the existence and uniqueness of discrete solutions, best approximation estimates and convergence estimates that are optimal compared to the regularity of the solution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cut Finite Element Methods for Partial Differential Equations on Embedded Manifolds of Arbitrary Codimensions

We develop a theoretical framework for the analysis of stabilized cut finite element methods for the Laplace-Beltrami operator on a manifold embedded in $\mathbb{R}^d$ of arbitrary codimension. The method is based on using continuous piecewise polynomials on a background mesh in the embedding space for approximation together with a stabilizing form that ensures that the resulting problem is stable. The discrete manifold is represented using a triangulation which does not match the background mesh and does not need to be shape-regular, which includes level set descriptions of codimension one manifolds and the non-matching embedding of independently triangulated manifolds as special cases. We identify abstract key assumptions on the stabilizing form which allow us to prove a bound on the condition number of the stiffness matrix and optimal order a priori estimates. The key assumptions are verified for three different realizations of the stabilizing form including a novel stabilization approach based on penalizing the surface normal gradient on the background mesh. Finally, we present numerical results illustrating our results for a curve and a surface embedded in $\mathbb{R}^3$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Data assimilation for the heat equation using stabilized finite element methods

We consider data assimilation for the heat equation using a finite element space semi-discretization. The approach is optimization based, but the design of regularization operators and parameters rely on techniques from the theory of stabilized finite elements. The space semi-discretized system is shown to admit a unique solution. Combining sharp estimates of the numerical stability of the discrete scheme and conditional stability estimates of the ill-posed continuous pde-model we then derive error estimates that reflect the approximation order of the finite element space and the stability of the continuous model. Two different data assimilation situations with different stability properties are considered to illustrate the framework. Full detail on how to adapt known stability estimates for the continuous model to work with the numerical analysis framework is given in appendix.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fictitious domain method with boundary value correction using penalty-free Nitsche method

In this paper, we consider a fictitious domain approach based on a Nitsche type method without penalty. To allow for high order approximation using piecewise affine approximation of the geometry we use a boundary value correction technique based on Taylor expansion from the approximate to the physical boundary. To ensure stability of the method a ghost penalty stabilization is considered in the boundary zone. We prove optimal error estimates in the $H^1$-norm and estimates suboptimal by $\mathcal{O}(h^{\frac12})$ in the $L^2$-norm. The suboptimality is due to the lack of adjoint consistency of our formulation. Numerical results are provided to corroborate the theoretical study.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Full Gradient Stabilized Cut Finite Element Methods for Surface Partial Differential Equations

We propose and analyze a new stabilized cut finite element method for the Laplace-Beltrami operator on a closed surface. The new stabilization term provides control of the full $\mathbb{R}^3$ gradient on the active mesh consisting of the elements that intersect the surface. Compared to face stabilization, based on controlling the jumps in the normal gradient across faces between elements in the active mesh, the full gradient stabilization is easier to implement and does not significantly increase the number of nonzero elements in the mass and stiffness matrices. The full gradient stabilization term may be combined with a variational formulation of the Laplace-Beltrami operator based on tangential or full gradients and we present a simple and unified analysis that covers both cases. The full gradient stabilization term gives rise to a consistency error which, however, is of optimal order for piecewise linear elements, and we obtain optimal order a priori error estimates in the energy and $L^2$ norms as well as an optimal bound of the condition number. Finally, we present detailed numerical examples where we in particular study the sensitivity of the condition number and error on the stabilization parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galerkin least squares finite element method for the obstacle problem

We construct a consistent multiplier free method for the finite element solution of the obstacle problem. The method is based on an augmented Lagrangian formulation in which we eliminate the multiplier by use of its definition in a discrete setting. We prove existence and uniqueness of discrete solutions and optimal order a priori error estimates for smooth exact solutions. Using a saturation assumption we also prove an a posteriori error estimate. Numerical examples show the performance of the method and of an adaptive algorithm for the control of the discretization error.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Solving ill-posed control problems by stabilized finite element methods: an alternative to Tikhonov regularization

Tikhonov regularization is one of the most commonly used methods of regularization of ill-posed problems. In the setting of finite element solutions of elliptic partial differential control problems, Tikhonov regularization amounts to adding suitably weighted least squares terms of the control variable, or derivatives thereof, to the Lagrangian determining the optimality system. In this note we show that stabilization methods for discretely ill--posed problems developed in the setting of convection--dominated convection--diffusion problems, can be highly suitable for stabilizing optimal control problems, and that Tikhonov regularization will lead to less accurate discrete solutions. We consider data assimilation problems for Poisson's equation as illustration and derive new error estimates both for the the reconstruction of the solution from measured data and reconstruction of the source term from measured data. These estimates include both the effect of discretization error and error in measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stabilized nonconforming finite element methods for data assimilation in incompressible flows

We consider a stabilized nonconforming finite element method for data assimilation in incompressible flow subject to the Stokes' equations. The method uses a primal dual structure that allows for the inclusion of nonstandard data. Error estimates are obtained that are optimal compared to the conditional stability of the ill-posed data assimilation problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Penalty Free Nitsche Method and Nonconforming Finite Elements for the Signorini Problem

We design and analyse a Nitsche method for contact problems. Compared to the seminal work of Chouly and Hild (A Nitsche-based method for unilateral contact problems: numerical analysis. SIAM J. Numer. Anal. 51 (2013), no. 2) our method is constructed by expressing the contact conditions in a nonlinear function for the displacement variable instead of the lateral forces. The contact condition is then imposed using the nonsymmetric variant of Nitsche's method that does not require a penalty term for stability. Nonconforming piecewise affine elements are considered for the bulk discretization. We prove optimal error estimates in the energy norm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Cut Discontinuous Galerkin Method for the Laplace-Beltrami Operator

We develop a discontinuous cut finite element method (CutFEM) for the Laplace-Beltrami operator on a hypersurface embedded in $\mathbb{R}^d$. The method is constructed by using a discontinuous piecewise linear finite element space defined on a background mesh in $\mathbb{R}^d$. The surface is approximated by a continuous piecewise linear surface that cuts through the background mesh in an arbitrary fashion. Then a discontinuous Galerkin method is formulated on the discrete surface and in order to obtain coercivity, certain stabilization terms are added on the faces between neighboring elements that provide control of the discontinuity as well as the jump in the gradient. We derive optimal a priori error and condition number estimates which are independent of the positioning of the surface in the background mesh. Finally, we present numerical examples confirming our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Cut Finite Element Method with Boundary Value Correction

In this contribution we develop a cut finite element method with boundary value correction of the type originally proposed by Bramble, Dupont, and Thomee. The cut finite element method is a fictitious domain method with Nitsche type enforcement of Dirichlet conditions together with stabilization of the elements at the boundary which is stable and enjoy optimal order approximation properties. A computational difficulty is, however, the geometric computations related to quadrature on the cut elements which must be accurate enough to achieve higher order approximation. With boundary value correction we may use only a piecewise linear approximation of the boundary, which is very convenient in a cut finite element method, and still obtain optimal order convergence. The boundary value correction is a modified Nitsche formulation involving a Taylor expansion in the normal direction compensating for the approximation of the boundary. Key to the analysis is a consistent stabilization term which enables us to prove stability of the method and a priori error estimates with explicit dependence on the meshsize and distance between the exact and approximate boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A penalty free Nitsche type method for the weak imposition of boundary conditions in compressible and incompressible elasticity

In this paper, we study the stability of the nonsymmetric version of Nitsche's method without penalty for compressible and incompressible elasticity. For the compressible case we prove the convergence of the error in the $H^1$- and $L^2$-norms. In the incompressible case we use a Galerkin least squares pressure stabilization and we prove the convergence in the $H^1$-norm for the velocity and convergence of the pressure in the $L^2$-norm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Stabilized Cut Finite Element Method for the Three Field Stokes Problem

We propose a Nitsche-based fictitious domain method for the three field Stokes problem in which the boundary of the domain is allowed to cross through the elements of a fixed background mesh. The dependent variables of velocity, pressure and extra-stress tensor are discretised on the background mesh using linear finite elements. This equal order approximation is stabilized using a continuous interior penalty (CIP) method. On the unfitted domain boundary, Dirichlet boundary conditions are weakly enforced using Nitsche's method. We add CIP-like ghost penalties in the boundary region and prove that our scheme is inf-sup stable and that it has optimal convergence properties independent of how the domain boundary intersects the mesh. Additionally, we demonstrate that the condition number of the system matrix is bounded independently of the boundary location. We corroborate our theoretical findings with numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Edge-based nonlinear diffusion for finite element approximations of convection-diffusion equations and its relation to algebraic flux-correction schemes

For the case of approximation of convection--diffusion equations using piecewise affine continuous finite elements a new edge-based nonlinear diffusion operator is proposed that makes the scheme satisfy a discrete maximum principle. The diffusion operator is shown to be Lipschitz continuous and linearity preserving. Using these properties we provide a full stability and error analysis, which, in the diffusion dominated regime, shows existence, uniqueness and optimal convergence. Then the algebraic flux correction method is recalled and we show that the present method can be interpreted as an algebraic flux correction method for a particular definition of the flux limiters. The performance of the method is illustrated on some numerical test cases in two space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stabilised finite element methods for ill-posed problems with conditional stability

In this paper we discuss the adjoint stabilised finite element method introduced in, E. Burman, Stabilized finite element methods for nonsymmetric, noncoercive and ill-posed problems. Part I: elliptic equations, SIAM Journal on Scientific Computing, and how it may be used for the computation of solutions to problems for which the standard stability theory given by the Lax-Milgram Lemma or the Babuska-Brezzi Theorem fails. We pay particular attention to ill-posed problems that have some conditional stability property and prove (conditional) error estimates in an abstract framework. As a model problem we consider the elliptic Cauchy problem and provide a complete numerical analysis for this case. Some numerical examples are given to illustrate the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stabilized CutFEM for the Convection Problem on Surfaces

We develop a stabilized cut finite element method for the convection problem on a surface based on continuous piecewise linear approximation and gradient jump stabilization terms. The discrete piecewise linear surface cuts through a background mesh consisting of tetrahedra in an arbitrary way and the finite element space consists of piecewise linear continuous functions defined on the background mesh. The variational form involves integrals on the surface and the gradient jump stabilization term is defined on the full faces of the tetrahedra. The stabilization term serves two purposes: first the method is stabilized and secondly the resulting linear system of equations is algebraically stable. We establish stability results that are analogous to the standard meshed flat case and prove $h^{3/2}$ order convergence in the natural norm associated with the method and that the full gradient enjoys $h^{3/4}$ order of convergence in $L^2$. We also show that the condition number of the stiffness matrix is bounded by $h^{-2}$. Finally, our results are verified by numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A stabilized cut finite element method for partial differential equations on surfaces: The Laplace-Beltrami operator

We consider solving the Laplace-Beltrami problem on a smooth two dimensional surface embedded into a three dimensional space meshed with tetrahedra. The mesh does not respect the surface and thus the surface cuts through the elements. We consider a Galerkin method based on using the restrictions of continuous piecewise linears defined on the tetrahedra to the surface as trial and test functions. The resulting discrete method may be severely ill-conditioned, and the main purpose of this paper is to suggest a remedy for this problem based on adding a consistent stabilization term to the original bilinear form. We show optimal estimates for the condition number of the stabilized method independent of the location of the surface. We also prove optimal a priori error estimates for the stabilized method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A stabilized nonconforming finite element method for the elliptic Cauchy problem

In this paper we propose a nonconforming finite element method for the solution of the ill-posed elliptic Cauchy problem. We prove error estimates using continuous dependence estimates in the $L^2$-norm. The effect of perturbations in data on the estimates is investigated. The recently derived framework from \cite{Bu13,Bu14} is extended to include the case of nonconforming approximation spaces and we show that the use of such spaces allows us to reduce the amount of stabilization necessary for convergence, even in the case of ill-posed problems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cut Finite Element Methods for Coupled Bulk-Surface Problems

We develop a cut finite element method for a second order elliptic coupled bulk-surface model problem. We prove a priori estimates for the energy and $L^2$ norms of the error. Using stabilization terms we show that the resulting algebraic system of equations has a similar condition number as a standard fitted finite element method. Finally, we present a numerical example illustrating the accuracy and the robustness of our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Error estimates for shock capturing finite element approximations of the one dimensional Burgers' equation

We propose an error analysis in weak norms of a shock capturing finite element method for the Burgers' equation. The estimates can be related to estimates of certain filtered quantities and are robust in the inviscid limit. Using a total variation apriori bound on the discrete solution and an interpolation inequality error estimates in $L^p$-norms can are obtained using interpolation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Error estimates for stabilized finite element methods applied to ill-posed problems

We propose an analysis for the stabilized finite element methods proposed in, E. Burman, Stabilized finite element methods for nonsymmetric, noncoercive, and ill-posed problems. Part I: Elliptic equations. SIAM J. Sci. Comput., 35(6) 2013, valid in the case of ill-posed problems for which only weak continuous dependence can be assumed. A priori and a posteriori error estimates are obtained without assuming coercivity or inf-sup stability of the continuous problem. A numerical example illustrates the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Robust error estimates in weak norms for advection dominated transport problems with rough data

We consider mixing problems in the form of transient convection--diffusion equations with a velocity vector field with multiscale character and rough data. We assume that the velocity field has two scales, a coarse scale with slow spatial variation, which is responsible for advective transport and a fine scale with small amplitude that contributes to the mixing. For this problem we consider the estimation of filtered error quantities for solutions computed using a finite element method with symmetric stabilization. A posteriori error estimates and a priori error estimates are derived using the multiscale decomposition of the advective velocity to improve stability. All estimates are independent both of the Péclet number and of the regularity of the exact solution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stabilized finite element methods for nonsymmetric, noncoercive and ill-posed problem. Part II: hyperbolic equations

In this paper we consider stabilised finite element methods for hyperbolic transport equations without coercivity. Abstract conditions for the convergence of the methods are introduced and these conditions are shown to hold for three different stabilised methods: the Galerkin least squares method, the continuous interior penalty method and the discontinuous Galerkin method. We consider both the standard stabilisation methods and the optimisation based method introduced in \cite{part1}. The main idea of the latter is to write the stabilised method in an optimisation framework and select the discrete function for which a certain cost functional, in our case stabilisation term, is minimised. Some numerical examples illustrate the theoretical investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Projection stabilisation of Lagrange multipliers for the imposition of constraints on interfaces and boundaries

Projection stabilisation applied to general Lagrange multiplier finite element methods is introduced and analysed in an abstract framework. We then consider some applications of the stabilised methods: (i) the weak imposition of boundary conditions, (ii) multiphysics coupling on unfitted meshes, (iii) a new interpretation of the classical residual stabilised Lagrange multiplier method introduced in H. J. C. Barbosa and T. J. R. Hughes, The finite element method with Lagrange multipliers on the boundary: circumventing the Babuška-Brezzi condition, Comput. Methods Appl. Mech. Engrg., 85(1):109--128, 1991 .

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Robust error estimates for stabilized finite element approximations of the two dimensional Navier-Stokes equations with application to implicit large eddy simulation

We consider error estimates in weak parametrised norms for stabilized finite element approximations of the two-dimensional Navier-Stokes' equations. These weak norms can be related to the norms of certain filtered quantities, where the parameter of the norm, relates to the filter width. Under the assumption of the existence of a certain decomposition of the solution, into large eddies and fine scale fluctuations, the constants of the estimates are proven to be independent of both the Reynolds number and the Sobolev norm of the exact solution. Instead they exhibit exponential growth with a coefficient proportional to the maximum gradient of the large eddies. The error estimates are on a posteriori form, but using Sobolev injections valid on finite element spaces and the properties of the stabilization operators the residuals may be upper bounded uniformly, leading to robust a priori error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stabilised finite element methods for non-symmetric, non-coercive and ill-posed problems. Part I: elliptic equations

In this paper we propose a new method to stabilise non-symmetric indefinite problems. The idea is to solve a forward and an adjoint problem simultaneously using a suitable stabilised finite element method. Both stabilisation of the element residual and jumps of certain derivatives of the discrete solution over element faces may be used. Under the assumption of well posedness of the partial differential equation and its associated adjoint problem we prove optimal error estimates in $H^1$ and $L^2$ norms in an abstract framework. Some examples of problems that are neither symmetric nor coercive, but that enter the abstract framework are given. First we treat indefinite convection-diffusion equations, with non-solenoidal transport velocity and either pure Dirichlet conditions or pure Neumann conditions and then a Cauchy problem for the Helmholtz operator. Some numerical illustrations are given.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Continuous Interior Penalty Finite Element Method for Helmholtz Equation with High Wave Number: One Dimensional Analysis

This paper addresses the properties of Continuous Interior Penalty (CIP) finite element solutions for the Helmholtz equation. The $h$-version of the CIP finite element method with piecewise linear approximation is applied to a one-dimensional model problem. We first show discrete well posedness and convergence results, using the imaginary part of the stabilization operator, for the complex Helmholtz equation. Then we consider a method with real valued penalty parameter and prove an error estimate of the discrete solution in the $H^1$-norm, as the sum of best approximation plus a pollution term that is the order of the phase difference. It is proved that the pollution can be eliminated by selecting the penalty parameter appropriately. As a result of this analysis, thorough and rigorous understanding of the error behavior throughout the range of convergence is gained. Numerical results are presented that show sharpness of the error estimates and highlight some phenomena of the discrete solution behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A penalty free non-symmetric Nitsche type method for the weak imposition of boundary conditions

In this note we show that the non-symmetric version of the classical Nitsche's method for the weak imposition of boundary conditions is stable without penalty term. We prove optimal $H^1$-error estimates and $L^2$-estimates that are suboptimal with half an order in $h$. Both the pure diffusion and the convection--diffusion problems are discussed.

0 citations0 reviewsNo code linkedOpen access