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

Peter Hansbo

Peter Hansbo 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 Science

Catalog footprint

What is connected

35works
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

35 published item(s)

preprint2022arXiv

A simple nonconforming tetrahedral element for the Stokes equations

In this paper we apply a nonconforming rotated bilinear tetrahedral element to the Stokes problem in $\mathbb{R}^3$. We show that the element is stable in combination with a piecewise linear, continuous, approximation of the pressure. This gives an approximation similar to the well known continuous $P^2-P^1$ Taylor$-$Hood element, but with fewer degrees of freedom. The element is a stable non-conforming low order element which fulfils Korn's inequality, leading to stability also in the case where the Stokes equations are written on stress form for use in the case of free surface flow.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Augmented Lagrangian approach to deriving discontinuous Galerkin methods for nonlinear elasticity problems

We use the augmented Lagrangian formalism to derive discontinuous Galerkin formulations for problems in nonlinear elasticity. In elasticity stress is typically a symmetric function of strain, leading to symmetric tangent stiffness matrices in Newtons method when conforming finite elements are used for discretization. By use of the augmented Lagrangian framework, we can also obtain symmetric tangent stiffness matrices in discontinuous Galerkin methods. We suggest two different approaches and give examples from plasticity and from large deformation hyperelasticity.

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

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

Nitsche's Finite Element Method for Model Coupling in Elasticity

We develop a Nitsche finite element method for a model of Euler--Bernoulli beams with axial stiffness embedded in a two--dimensional elastic bulk domain. The beams have their own displacement fields, and the elastic subdomains created by the beam network are triangulated independently and are coupled to the beams weakly by use of Nitsche's method in the framework of hybridization.

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

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

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
preprint2016arXiv

A Nitsche-type Method for Helmholtz Equation with an Embedded Acoustically Permeable Interface

We propose a new finite element method for Helmholtz equation in the situation where an acoustically permeable interface is embedded in the computational domain. A variant of Nitsche's method, different from the standard one, weakly enforces the impedance conditions for transmission through the interface. As opposed to a standard finite-element discretization of the problem, our method seamlessly handles a complex-valued impedance function $Z$ that is allowed to vanish. In the case of a vanishing impedance, the proposed method reduces to the classic Nitsche method to weakly enforce continuity over the interface. We show stability of the method, in terms of a discrete Gårding inequality, for a quite general class of surface impedance functions, provided that possible surface waves are sufficiently resolved by the mesh. Moreover, we prove an a priori error estimate under the assumption that the absolute value of the impedance is bounded away from zero almost everywhere. Numerical experiments illustrate the performance of the method for a number of test cases in 2D and 3D with different interface conditions.

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

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

Optimal Design of Fibre Reinforced Membrane Structures

A design problem of finding an optimally stiff membrane structure by selecting one-dimensional fiber reinforcements is formulated and solved. The membrane model is derived in a novel manner from a particular three-dimensional linear elastic orthotropic model by appropriate assumptions. The design problem is given in the form of two minimization statements, reminiscent of a Nash game. After finite element discretization, the separate treatment of each of the two minimization statements follows from classical results and methods of structural optimization: the stiffest orientation of reinforcing fibers coincides with principal stresses and the separate selection of density of fibers is a convex problem that can be solved by optimality criteria iterations. Numerical solutions are shown for two particular configurations. The first for a statically determined structure and the second for a statically undetermined one. The latter shows related but non-unique solutions.

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 for coupled bulk-surface problems on time-dependent domains

In this contribution we present a new computational method for coupled bulk-surface problems on time-dependent domains. The method is based on a space-time formulation using discontinuous piecewise linear elements in time and continuous piecewise linear elements in space on a fixed background mesh. The domain is represented using a piecewise linear level set function on the background mesh and a cut finite element method is used to discretize the bulk and surface problems. In the cut finite element method the bilinear forms associated with the weak formulation of the problem are directly evaluated on the bulk domain and the surface defined by the level set, essentially using the restrictions of the piecewise linear functions to the computational domain. In addition a stabilization term is added to stabilize convection as well as the resulting algebraic system that is solved in each time step. We show in numerical examples that the resulting method is accurate and stable and results in well conditioned algebraic systems independent of the position of the interface relative to the background mesh.

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 posteriori error estimates for continuous/discontinuous Galerkin approximations of the Kirchhoff-Love buckling problem

Second order buckling theory involves a one-way coupled coupled problem where the stress tensor from a plane stress problem appears in an eigenvalue problem for the fourth order Kirchhoff plate. In this paper we present an a posteriori error estimate for the critical buckling load and mode corresponding to the smallest eigenvalue and associated eigenvector. A particular feature of the analysis is that we take the effect of approximate computation of the stress tensor and also provide an error indicator for the plane stress problem. The Kirchhoff plate is discretized using a continuous/discontinuous finite element method which uses standard continuous piecewise polynomial finite element spaces which can also be used to solve the plane stress problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Stabilized Finite Element Method for the Darcy Problem on Surfaces

We consider a stabilized finite element method for the Darcy problem on a surface based on the Masud-Hughes formulation. A special feature of the method is that the tangential condition of the velocity field is weakly enforced through the bilinear form and that standard parametric continuous polynomial spaces on triangulations can be used. We prove optimal order a priori estimates that take the approximation of the geometry and the solution into account.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cut finite element modeling of linear membranes

We construct a cut finite element method for the membrane elasticity problem on an embedded mesh using tangential differential calculus. Both free membranes and membranes coupled to 3D elasticity are considered. The discretization comes from a Galerkin method using the restriction of 3D basis funtions (linear or trilinear) to the surface representing the membrane. In the case of coupling to 3D elasticity, we view the membrane as giving additional stiffness contributions to the standard stiffness matrix resulting from the discretization of the three-dimensional continuum.

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 cut finite element method for a Stokes interface problem

We present a finite element method for the Stokes equations involving two immiscible incompressible fluids with different viscosities and with surface tension. The interface separating the two fluids does not need to align with the mesh. We propose a Nitsche formulation which allows for discontinuities along the interface with optimal a priori error estimates. A stabilization procedure is included which ensures that the method produces a well conditioned stiffness matrix independent of the location of the interface.

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

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

Minimal surface computation using a finite element method on an embedded surface

We suggest a finite element method for computing minimal surfaces based on computing a discrete Laplace-Beltrami operator operating on the coordinates of the surface. The surface is a discrete representation of the zero level set of a distance function using linear tetrahedral finite elements, and the finite element discretization is done on the piecewise planar isosurface using the shape functions from the background three dimensional mesh used to represent the distance function. A recently suggested stabilization scheme is a crucial component in the method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stabilized Finite Element Approximation of the Mean Curvature Vector on Closed Surfaces

We develop a stabilized discrete Laplace-Beltrami operator that is used to compute an approximate mean curvature vector which enjoys convergence of order one in L2. The stabilization is of gradient jump type and we consider both standard meshed surfaces and so called cut surfaces that are level sets of piecewise linear distance functions. We prove a priori error estimates and verify the theoretical results numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Locking Free Quadrilateral Continuous/Discontinuous Finite Element Methods for the Reissner-Mindlin Plate

We develop a finite element method with continuous displacements and discontinuous rotations for the Mindlin-Reissner plate model on quadrilateral elements. To avoid shear locking, the rotations must have the same polynomial degree in the parametric reference plane as the parametric derivatives of the displacements, and obey the same transformation law to the physical plane as the gradient of displacements. We prove optimal convergence, uniformly in the plate thickness, and provide numerical results that confirm our estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variational Formulation of Curved Beams in Global Coordinates

In this paper we derive a variational formulation for a linear curved beam which is natively expressed in global Cartesian coordinates. During derivation the beam midline is assumed to be implicitly described by a vector distance function which eliminates the need for local coordinates. The only geometrical information appearing in the final expressions for the governing equations is the tangential direction, and thus there is no need to introduce normal directions along the curve. As a consequence zero or discontinuous curvature, for example at inflection points, pose no difficulty in this formulation. Kinematic assumptions encompassing both Timoshenko and Euler--Bernoulli beam theories are considered. With the exception of truly three dimensional formulations, models for curved beams found in literature are typically derived in the Frenet frame defined by the geometry of the beam midline. While it is intuitive to formulate curved beam models in these local coordinates, the Frenet frame suffers from ambiguity and sudden changes of orientation in straight sections of the beam. Based on the variational formulation we implement finite element models using global Cartesian degrees of freedom and discuss curvature coupling effects and locking. Numerical comparisons with classical solutions for both straight and curved cantilever beams under a tip load are given, as well as numerical examples illustrating curvature coupling effects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Intrinsic finite element modeling of a linear membrane shell problem

A Galerkin finite element method for the membrane elasticity problem on a meshed surface is constructed by using two-dimensional elements extended into three dimensions. The membrane finite element model is established using the intrinsic approach suggested by [Delfour and Zolésio, A boundary differential equation for thin shells. J. Differential Equations, 119(2):426--449, 1995].

0 citations0 reviewsNo code linkedOpen access