-
Error Estimates for Finite Element Approximations of Viscoelastic Dynamics: The Generalized Maxwell Model
Authors:
Martin Björklund,
Karl Larsson,
Mats G. Larson
Abstract:
We prove error estimates for a finite element approximation of viscoelastic dynamics based on continuous Galerkin in space and time, both in energy norm and in $L^2$ norm. The proof is based on an error representation formula using a discrete dual problem and a stability estimate involving the kinetic, elastic, and viscoelastic energies. To set up the dual error analysis and to prove the basic sta…
▽ More
We prove error estimates for a finite element approximation of viscoelastic dynamics based on continuous Galerkin in space and time, both in energy norm and in $L^2$ norm. The proof is based on an error representation formula using a discrete dual problem and a stability estimate involving the kinetic, elastic, and viscoelastic energies. To set up the dual error analysis and to prove the basic stability estimates, it is natural to formulate the problem as a system involving evolution equations for the viscoelastic stress, the displacements, and the velocities. The equations for the viscoelastic stress can, however, be solved analytically in terms of the deviatoric strain velocity, and therefore, the viscoelastic stress can be eliminated from the system, resulting in a system for displacements and velocities.
△ Less
Submitted 12 March, 2024; v1 submitted 21 November, 2023;
originally announced November 2023.
-
Critical time-step size analysis and mass scaling by ghost-penalty for immersogeometric explicit dynamics
Authors:
Stein K. F. Stoter,
Sai C. Divi,
E. Harald van Brummelen,
Mats G. Larson,
Frits de Prenter,
Clemens V. Verhoosel
Abstract:
In this article, we study the effect of small-cut elements on the critical time-step size in an immersogeometric context. We analyze different formulations for second-order (membrane) and fourth-order (shell-type) equations, and derive scaling relations between the critical time-step size and the cut-element size for various types of cuts. In particular, we focus on different approaches for the we…
▽ More
In this article, we study the effect of small-cut elements on the critical time-step size in an immersogeometric context. We analyze different formulations for second-order (membrane) and fourth-order (shell-type) equations, and derive scaling relations between the critical time-step size and the cut-element size for various types of cuts. In particular, we focus on different approaches for the weak imposition of Dirichlet conditions: by penalty enforcement and with Nitsche's method. The stability requirement for Nitsche's method necessitates either a cut-size dependent penalty parameter, or an additional ghost-penalty stabilization term is necessary. Our findings show that both techniques suffer from cut-size dependent critical time-step sizes, but the addition of a ghost-penalty term to the mass matrix serves to mitigate this issue. We confirm that this form of `mass-scaling' does not adversely affect error and convergence characteristics for a transient membrane example, and has the potential to increase the critical time-step size by orders of magnitude. Finally, for a prototypical simulation of a Kirchhoff-Love shell, our stabilized Nitsche formulation reduces the solution error by well over an order of magnitude compared to a penalty formulation at equal time-step size.
△ Less
Submitted 14 February, 2023;
originally announced February 2023.
-
Space-Time CutFEM on Overlapping Meshes: Simple Continuous Mesh Motion
Authors:
Mats G. Larson,
Anders Logg,
Carl Lundholm
Abstract:
We present a cut finite element method for the heat equation on two overlapping meshes: a stationary background mesh and an overlapping mesh that moves around inside/"on top" of it. Here the overlapping mesh is prescribed a simple continuous motion, meaning that its location as a function of time is continuous and piecewise linear. For the discrete function space, we use continuous Galerkin in spa…
▽ More
We present a cut finite element method for the heat equation on two overlapping meshes: a stationary background mesh and an overlapping mesh that moves around inside/"on top" of it. Here the overlapping mesh is prescribed a simple continuous motion, meaning that its location as a function of time is continuous and piecewise linear. For the discrete function space, we use continuous Galerkin in space and discontinuous Galerkin in time, with the addition of a discontinuity on the boundary between the two meshes. The finite element formulation is based on Nitsche's method and also includes an integral term over the space-time boundary between the two meshes that mimics the standard discontinuous Galerkin time-jump term. The simple continuous mesh motion results in a space-time discretization for which standard analysis methodologies either fail or are unsuitable. We therefore employ what seems to be a relatively new energy analysis framework that is general and robust enough to be applicable to the current setting. The energy analysis consists of a stability estimate that is slightly stronger than the standard basic one and an a priori error estimate that is of optimal order with respect to both time step and mesh size. We also present numerical results for a problem in one spatial dimension that verify the analytic error convergence orders.
△ Less
Submitted 1 March, 2023; v1 submitted 1 February, 2023;
originally announced February 2023.
-
Space-Time CutFEM on Overlapping Meshes: Simple Discontinuous Mesh Evolution
Authors:
Mats G. Larson,
Carl Lundholm
Abstract:
We present a cut finite element method for the heat equation on two overlapping meshes: a stationary background mesh and an overlapping mesh that evolves inside/"on top" of it. Here the overlapping mesh is prescribed a simple discontinuous evolution, meaning that its location, size, and shape as functions of time are discontinuous and piecewise constant. For the discrete function space, we use con…
▽ More
We present a cut finite element method for the heat equation on two overlapping meshes: a stationary background mesh and an overlapping mesh that evolves inside/"on top" of it. Here the overlapping mesh is prescribed a simple discontinuous evolution, meaning that its location, size, and shape as functions of time are discontinuous and piecewise constant. For the discrete function space, we use continuous Galerkin in space and discontinuous Galerkin in time, with the addition of a discontinuity on the boundary between the two meshes. The finite element formulation is based on Nitsche's method. The simple discontinuous mesh evolution results in a space-time discretization with a slabwise product structure between space and time which allows for existing analysis methodologies to be applied with only minor modifications. We follow the analysis methodology presented by Eriksson and Johnson in [1, 2]. The greatest modification is the introduction of a Ritzlike "shift operator" that is used to obtain the discrete strong stability needed for the error analysis. The shift operator generalizes the original analysis to some methods for which the discrete subspace at one time does not lie in the space of the stiffness form at the subsequent time. The error analysis consists of an a priori error estimate that is of optimal order with respect to both time step and mesh size. We also present numerical results for a problem in one spatial dimension that verify the analytic error convergence orders.
△ Less
Submitted 1 March, 2023; v1 submitted 31 January, 2023;
originally announced January 2023.
-
Hybridized Isogeometric Method for Elliptic Problems on CAD Surfaces with Gaps
Authors:
Tobias Jonsson,
Mats G. Larson,
Karl Larsson
Abstract:
We develop a method for solving elliptic partial differential equations on surfaces described by CAD patches that may have gaps/overlaps. The method is based on hybridization using a three-dimensional mesh that covers the gap/overlap between patches. Thus, the hybrid variable is defined on a three-dimensional mesh, and we need to add appropriate normal stabilization to obtain an accurate solution,…
▽ More
We develop a method for solving elliptic partial differential equations on surfaces described by CAD patches that may have gaps/overlaps. The method is based on hybridization using a three-dimensional mesh that covers the gap/overlap between patches. Thus, the hybrid variable is defined on a three-dimensional mesh, and we need to add appropriate normal stabilization to obtain an accurate solution, which we show can be done by adding a suitable term to the weak form. In practical applications, the hybrid mesh may be conveniently constructed using an octree to efficiently compute the necessary geometric information. We prove error estimates and present several numerical examples illustrating the application of the method to different problems, including a realistic CAD model.
△ Less
Submitted 7 March, 2023; v1 submitted 4 October, 2022;
originally announced October 2022.
-
Cut finite element method for divergence free approximation of incompressible flow: a Lagrange multiplier approach
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
In this note we design a cut finite element method for a low order divergence free element applied to a boundary value problem subject to Stokes' equations. For the imposition of Dirichlet boundary conditions we consider either Nitsche's method or a stabilized Lagrange multiplier method. In both cases the normal component of the velocity is constrained using a multiplier, different from the standa…
▽ More
In this note we design a cut finite element method for a low order divergence free element applied to a boundary value problem subject to Stokes' equations. For the imposition of Dirichlet boundary conditions we consider either Nitsche's method or a stabilized Lagrange multiplier method. In both cases the normal component of the velocity is constrained using a multiplier, different from the standard pressure approximation. The divergence of the approximate velocities is pointwise zero over the whole mesh domain, and we derive optimal error estimates for the velocity and pressures, where the error constant is independent of how the physical domain intersects the computational mesh, and of the regularity of the pressure multiplier imposing the divergence free condition.
△ Less
Submitted 13 November, 2023; v1 submitted 11 July, 2022;
originally announced July 2022.
-
The augmented Lagrangian method as a framework for stabilised methods in computational mechanics
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
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 d…
▽ More
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.
△ Less
Submitted 1 July, 2022;
originally announced July 2022.
-
Extension Operators for Trimmed Spline Spaces
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson,
Karl Larsson
Abstract:
We develop a discrete extension operator for trimmed spline spaces consisting of piecewise polynomial functions of degree $p$ with $k$ continuous derivatives. The construction is based on polynomial extension from neighboring elements together with projection back into the spline space. We prove stability and approximation results for the extension operator. Finally, we illustrate how we can use t…
▽ More
We develop a discrete extension operator for trimmed spline spaces consisting of piecewise polynomial functions of degree $p$ with $k$ continuous derivatives. The construction is based on polynomial extension from neighboring elements together with projection back into the spline space. We prove stability and approximation results for the extension operator. Finally, we illustrate how we can use the extension operator to construct a stable cut isogeometric method for an elliptic model problem.
△ Less
Submitted 13 May, 2022;
originally announced May 2022.
-
A simple nonconforming tetrahedral element for the Stokes equations
Authors:
Peter Hansbo,
Mats G. Larson
Abstract:
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-conf…
▽ More
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.
△ Less
Submitted 4 May, 2022;
originally announced May 2022.
-
On the Design of Locking Free Ghost Penalty Stabilization and the Relation to CutFEM with Discrete Extension
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
In this note, we develop a new stabilization mechanism for cut finite element methods that generalizes previous approaches of ghost penalty type in two ways: (1) The quantity that is stabilized and (2) The choice of elements that are connected in the stabilization. In particular, we can stabilize functionals of the discrete function such as finite element degrees of freedom. We subsequently show t…
▽ More
In this note, we develop a new stabilization mechanism for cut finite element methods that generalizes previous approaches of ghost penalty type in two ways: (1) The quantity that is stabilized and (2) The choice of elements that are connected in the stabilization. In particular, we can stabilize functionals of the discrete function such as finite element degrees of freedom. We subsequently show that the kernel of our ghost penalty operator defines a finite element space based on discrete extensions in the spirit of those introduced in Burman, E.; Hansbo, P. and Larson, M. G., CutFEM Based on Extended Finite Element Spaces, arXiv2101.10052, 2021.
△ Less
Submitted 3 May, 2022;
originally announced May 2022.
-
Augmented Lagrangian approach to deriving discontinuous Galerkin methods for nonlinear elasticity problems
Authors:
Peter Hansbo,
Mats G. Larson
Abstract:
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 tangen…
▽ More
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.
△ Less
Submitted 17 February, 2022;
originally announced February 2022.
-
The Finite Cell Method with Least Squares Stabilized Nitsche Boundary Conditions
Authors:
Karl Larsson,
Stefan Kollmannsberger,
Ernst Rank,
Mats G. Larson
Abstract:
We apply the recently developed least squares stabilized symmetric Nitsche method for enforcement of Dirichlet boundary conditions to the finite cell method. The least squares stabilized Nitsche method in combination with finite cell stabilization leads to a symmetric positive definite stiffness matrix and relies only on elementwise stabilization, which does not lead to additional fill in. We prov…
▽ More
We apply the recently developed least squares stabilized symmetric Nitsche method for enforcement of Dirichlet boundary conditions to the finite cell method. The least squares stabilized Nitsche method in combination with finite cell stabilization leads to a symmetric positive definite stiffness matrix and relies only on elementwise stabilization, which does not lead to additional fill in. We prove a priori error estimates and bounds on the condition numbers.
△ Less
Submitted 27 October, 2021;
originally announced October 2021.
-
Nitsche's Finite Element Method for Model Coupling in Elasticity
Authors:
Peter Hansbo,
Mats G. Larson
Abstract:
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.
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.
△ Less
Submitted 2 June, 2021;
originally announced June 2021.
-
Conservative Discontinuous Cut Finite Element Methods
Authors:
Mats G. Larson,
Sara Zahedi
Abstract:
We develop a conservative cut finite element method for an elliptic coupled bulk-interface problem. The method is based on a discontinuous Galerkin framework where stabilization is added in such a way that we retain conservation on macro elements containing one element with a large intersection with the domain and possibly a number of elements with small intersections. We derive error estimates an…
▽ More
We develop a conservative cut finite element method for an elliptic coupled bulk-interface problem. The method is based on a discontinuous Galerkin framework where stabilization is added in such a way that we retain conservation on macro elements containing one element with a large intersection with the domain and possibly a number of elements with small intersections. We derive error estimates and present confirming numerical results
△ Less
Submitted 5 May, 2021;
originally announced May 2021.
-
Error estimates for the Smagorinsky turbulence model: enhanced stability through scale separation and numerical stabilization
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
In the present work we show some results on the effect of the Smagorinsky model on the stability of the associated perturbation equation. We show that in the presence of a spectral gap, such that the flow can be decomposed in a large scale with moderate gradient and a small amplitude fine scale with arbitratry gradient, the Smagorinsky model admits stability estimates for perturbations, with expon…
▽ More
In the present work we show some results on the effect of the Smagorinsky model on the stability of the associated perturbation equation. We show that in the presence of a spectral gap, such that the flow can be decomposed in a large scale with moderate gradient and a small amplitude fine scale with arbitratry gradient, the Smagorinsky model admits stability estimates for perturbations, with exponential growth depending only on the large scale gradient. We then show in the context of stabilized finite element methods that the same result carries over to the approximation and that in this context, for suitably chosen finite element spaces the Smagorinsky model acts as a stabilizer yielding close to optimal error estimates in the $L^2$-norm for smooth flows in the pre-asymptotic high Reynolds number regime.
△ Less
Submitted 29 January, 2021;
originally announced February 2021.
-
CutFEM Based on Extended Finite Element Spaces
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
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…
▽ More
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.
△ Less
Submitted 25 January, 2021;
originally announced January 2021.
-
Explicit Time Stepping for the Wave Equation using CutFEM with Discrete Extension
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
In this note we develop a fully explicit cut finite element method for the wave equation. The method is based on using a standard leap frog scheme combined with an extension operator that defines the nodal values outside of the domain in terms of the nodal values inside the domain. We show that the mass matrix associated with the extended finite element space can be lumped leading to a fully expli…
▽ More
In this note we develop a fully explicit cut finite element method for the wave equation. The method is based on using a standard leap frog scheme combined with an extension operator that defines the nodal values outside of the domain in terms of the nodal values inside the domain. We show that the mass matrix associated with the extended finite element space can be lumped leading to a fully explicit scheme. We derive stability estimates for the method and provide optimal order a priori error estimates. Finally, we present some illustrating numerical examples.
△ Less
Submitted 10 November, 2020;
originally announced November 2020.
-
Comparison of Shape Derivatives using CutFEM for Ill-posed Bernoulli Free Boundary Problem
Authors:
Erik Burman,
Cuiyu He,
Mats G. Larson
Abstract:
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 optimiz…
▽ More
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.
△ Less
Submitted 21 August, 2020;
originally announced August 2020.
-
Low Regularity Estimates for CutFEM Approximations of an Elliptic Problem with Mixed Boundary Conditions
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
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 bounda…
▽ More
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.
△ Less
Submitted 6 July, 2020;
originally announced July 2020.
-
A cut finite element method for a model of pressure in fractured media
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
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 effici…
▽ More
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.
△ Less
Submitted 2 June, 2020;
originally announced June 2020.
-
Augmented Lagrangian Method for Thin Plates with Signorini Boundaries
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
We consider $C^1$-continuous approximations of the Kirchhoff plate problem in combination with a mesh dependent augmented Lagrangian method on a simply supported Signorini boundary.
We consider $C^1$-continuous approximations of the Kirchhoff plate problem in combination with a mesh dependent augmented Lagrangian method on a simply supported Signorini boundary.
△ Less
Submitted 15 December, 2019;
originally announced December 2019.
-
Cut Bogner-Fox-Schmit Elements for Plates
Authors:
Erik Burman,
Mats G. Larson,
Peter Hansbo
Abstract:
We present and analyze a method for thin plates based on cut Bogner-Fox-Schmit elements, which are $C^1$ elements obtained by taking tensor products of Hermite splines. The formulation is based on Nitsche's method for weak enforcement of essential boundary conditions together with addition of certain stabilization terms that enable us to establish coercivity and stability of the resulting system o…
▽ More
We present and analyze a method for thin plates based on cut Bogner-Fox-Schmit elements, which are $C^1$ elements obtained by taking tensor products of Hermite splines. The formulation is based on Nitsche's method for weak enforcement of essential boundary conditions together with addition of certain stabilization terms that enable us to establish coercivity and stability of the resulting system of linear equations. We also take geometric approximation of the boundary into account and we focus our presentation on the simply supported boundary conditions which is the most sensitive case for geometric approximation of the boundary.
△ Less
Submitted 1 November, 2019;
originally announced November 2019.
-
A Posteriori Error Estimates with Boundary Correction for a Cut Finite Element Method
Authors:
Erik Burman,
Cuiyu He,
Mats G. Larson
Abstract:
In this work we study a residual based a posteriori error estimation for the CutFEM method applied to an elliptic model problem. We consider the problem with non-polygonal boundary and the analysis takes into account the geometry and data approximation on the boundary. The reliability and efficiency are theoretically proved. Moreover, constants are robust with respect to how the domain boundary cu…
▽ More
In this work we study a residual based a posteriori error estimation for the CutFEM method applied to an elliptic model problem. We consider the problem with non-polygonal boundary and the analysis takes into account the geometry and data approximation on the boundary. The reliability and efficiency are theoretically proved. Moreover, constants are robust with respect to how the domain boundary cuts the mesh.
△ Less
Submitted 3 June, 2019;
originally announced June 2019.
-
Dirichlet Boundary Value Correction using Lagrange Multipliers
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
We propose a boundary value correction approach for cases when curved boundaries are approximated by straight lines (planes) and Lagrange multipliers are used to enforce Dirichlet boundary conditions. The approach allows for optimal order convergence for polynomial order up to 3. We show the relation to the Taylor series expansion approach used by Bramble, Dupont and Tomée [Math. Comp., 26:869--87…
▽ More
We propose a boundary value correction approach for cases when curved boundaries are approximated by straight lines (planes) and Lagrange multipliers are used to enforce Dirichlet boundary conditions. The approach allows for optimal order convergence for polynomial order up to 3. We show the relation to the Taylor series expansion approach used by Bramble, Dupont and Tomée [Math. Comp., 26:869--879, 1972] in the context of Nitsche's method and, in the case of inf-sup stable multiplier methods, prove a priori error estimates with explicit dependence on the meshsize and distance between the exact and approximate boundary.
△ Less
Submitted 17 March, 2019;
originally announced March 2019.
-
Graded Parametric CutFEM and CutIGA for Elliptic Boundary Value Problems in Domains with Corners
Authors:
Tobias Jonsson,
Mats G. Larson,
Karl Larsson
Abstract:
We develop a parametric cut finite element method for elliptic boundary value problems with corner singularities where we have weighted control of higher order derivatives of the solution to a neighborhood of a point at the boundary. Our approach is based on identification of a suitable mapping that grades the mesh towards the singularity. In particular, this mapping may be chosen without identify…
▽ More
We develop a parametric cut finite element method for elliptic boundary value problems with corner singularities where we have weighted control of higher order derivatives of the solution to a neighborhood of a point at the boundary. Our approach is based on identification of a suitable mapping that grades the mesh towards the singularity. In particular, this mapping may be chosen without identifying the opening angle at the corner. We employ cut finite elements together with Nitsche boundary conditions and stabilization in the vicinity of the boundary. We prove that the method is stable and convergent of optimal order in the energy norm and $L^2$ norm. This is achieved by mapping to the reference domain where we employ a structured mesh.
△ Less
Submitted 30 April, 2019; v1 submitted 20 December, 2018;
originally announced December 2018.
-
Hybridized CutFEM for Elliptic Interface Problems
Authors:
Erik Burman,
Daniel Elfverson,
Peter Hansbo,
Mats G. Larson,
Karl Larsson
Abstract:
We design and analyze a hybridized cut finite element method for elliptic interface problems. In this method very general meshes can be coupled over internal unfitted interfaces, through a skeletal variable, using a Nitsche type approach. We discuss how optimal error estimates for the method are obtained using the tools of cut finite element methods and prove a condition number estimate for the Sc…
▽ More
We design and analyze a hybridized cut finite element method for elliptic interface problems. In this method very general meshes can be coupled over internal unfitted interfaces, through a skeletal variable, using a Nitsche type approach. We discuss how optimal error estimates for the method are obtained using the tools of cut finite element methods and prove a condition number estimate for the Schur complement. Finally, we present illustrating numerical examples.
△ Less
Submitted 29 October, 2018;
originally announced October 2018.
-
A Stable Cut Finite Element Method for Partial Differential Equations on Surfaces: The Helmholtz-Beltrami Operator
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson,
Andre Massing
Abstract:
We consider solving the surface Helmholtz equation 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.…
▽ More
We consider solving the surface Helmholtz equation 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.Using a stabilized method combining Galerkin least squares stabilization and a penalty on the gradient jumps we obtain stability of the discrete formulation under the condition $h k < C$, where $h$ denotes the mesh size, $k$ the wave number and $C$ a constant depending mainly on the surface curvature $κ$, but not on the surface/mesh intersection. Optimal error estimates in the $H^1$ and $L^2$-norms follow.
△ Less
Submitted 9 October, 2018;
originally announced October 2018.
-
Cut Topology Optimization for Linear Elasticity with Coupling to Parametric Nondesign Domain Regions
Authors:
Erik Burman,
Daniel Elfverson,
Peter Hansbo,
Mats G. Larson,
Karl Larsson
Abstract:
We develop a density based topology optimization method for linear elasticity based on the cut finite element method. More precisely, the design domain is discretized using cut finite elements which allow complicated geometry to be represented on a structured fixed background mesh. The geometry of the design domain is allowed to cut through the background mesh in an arbitrary way and certain stabi…
▽ More
We develop a density based topology optimization method for linear elasticity based on the cut finite element method. More precisely, the design domain is discretized using cut finite elements which allow complicated geometry to be represented on a structured fixed background mesh. The geometry of the design domain is allowed to cut through the background mesh in an arbitrary way and certain stabilization terms are added in the vicinity of the cut boundary, which guarantee stability of the method. Furthermore, in addition to standard Dirichlet and Neumann conditions we consider interface conditions enabling coupling of the design domain to parts of the structure for which the design is already given. These given parts of the structure, called the nondesign domain regions, typically represents parts of the geometry provided by the designer. The nondesign domain regions may be discretized independently from the design domains using for example parametric meshed finite elements or isogeometric analysis. The interface and Dirichlet conditions are based on Nitsche's method and are stable for the full range of density parameters. In particular we obtain a traction-free Neumann condition in the limit when the density tends to zero.
△ Less
Submitted 6 March, 2019; v1 submitted 20 September, 2018;
originally announced September 2018.
-
A Stabilized Cut Streamline Diffusion Finite Element Method for Convection-Diffusion Problems on Surfaces
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson,
Andre Massing,
Sara Zahedi
Abstract:
We develop a stabilized cut finite element method for the stationary convection diffusion problem on a surface embedded in ${\mathbb{R}}^d$. The cut finite element method is based on using an embedding of the surface into a three dimensional mesh consisting of tetrahedra and then using the restriction of the standard piecewise linear continuous elements to a piecewise linear approximation of the s…
▽ More
We develop a stabilized cut finite element method for the stationary convection diffusion problem on a surface embedded in ${\mathbb{R}}^d$. The cut finite element method is based on using an embedding of the surface into a three dimensional mesh consisting of tetrahedra and then using the restriction of the standard piecewise linear continuous elements to a piecewise linear approximation of the surface. The stabilization consists of a standard streamline diffusion stabilization term on the discrete surface and a so called normal gradient stabilization term on the full tetrahedral elements in the active mesh. We prove optimal order a priori error estimates in the standard norm associated with the streamline diffusion method and bounds for the condition number of the resulting stiffness matrix. The condition number is of optimal order $O(h^{-1})$ for a specific choice of method parameters. Numerical example supporting our theoretical results are also included.
△ Less
Submitted 4 July, 2018;
originally announced July 2018.
-
A Cut Finite Element Method for Elliptic Bulk Problems with Embedded Surfaces
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson,
David Samvin
Abstract:
We propose an unfitted finite element method for flow in fractured porous media. The coupling across the fracture uses a Nitsche type mortaring, allowing for an accurate representation of the jump in the normal component of the gradient of the discrete solution across the fracture. The flow field in the fracture is modelled simultaneously, using the average of traces of the bulk variables on the f…
▽ More
We propose an unfitted finite element method for flow in fractured porous media. The coupling across the fracture uses a Nitsche type mortaring, allowing for an accurate representation of the jump in the normal component of the gradient of the discrete solution across the fracture. The flow field in the fracture is modelled simultaneously, using the average of traces of the bulk variables on the fractured. In particular the Laplace-Beltrami operator for the transport in the fracture is included using the average of the projection on the tangential plane of the fracture of the trace of the bulk gradient. Optimal order error estimates are proven under suitable regularity assumptions on the domain geometry. The extension to the case of bifurcating fractures is discussed. Finally the theory is illustrated by a series of numerical examples.
△ Less
Submitted 2 May, 2018;
originally announced May 2018.
-
A MultiMesh Finite Element Method for the Stokes Problem
Authors:
August Johansson,
Mats G. Larson,
Anders Logg
Abstract:
The multimesh finite element method enables the solution of partial differential equations on a computational mesh composed by multiple arbitrarily overlapping meshes. The discretization is based on a continuous--discontinuous function space with interface conditions enforced by means of Nitsche's method. In this contribution, we consider the Stokes problem as a first step towards flow application…
▽ More
The multimesh finite element method enables the solution of partial differential equations on a computational mesh composed by multiple arbitrarily overlapping meshes. The discretization is based on a continuous--discontinuous function space with interface conditions enforced by means of Nitsche's method. In this contribution, we consider the Stokes problem as a first step towards flow applications. The multimesh formulation leads to so called cut elements in the underlying meshes close to overlaps. These demand stabilization to ensure coercivity and stability of the stiffness matrix. We employ a consistent least-squares term on the overlap to ensure that the inf-sup condition holds. We here present the method for the Stokes problem, discuss the implementation, and verify that we have optimal convergence.
△ Less
Submitted 30 April, 2018;
originally announced May 2018.
-
MultiMesh Finite Elements with Flexible Mesh Sizes
Authors:
August Johansson,
Mats G. Larson,
Anders Logg
Abstract:
We analyze a new framework for expressing finite element methods on arbitrarily many intersecting meshes: multimesh finite element methods. The multimesh finite element method, first presented in [40], enables the use of separate meshes to discretize parts of a computational domain that are naturally separate; such as the components of an engine, the domains of a multiphysics problem, or solid bod…
▽ More
We analyze a new framework for expressing finite element methods on arbitrarily many intersecting meshes: multimesh finite element methods. The multimesh finite element method, first presented in [40], enables the use of separate meshes to discretize parts of a computational domain that are naturally separate; such as the components of an engine, the domains of a multiphysics problem, or solid bodies interacting under the influence of forces from surrounding fluids or other physical fields. Furthermore, each of these meshes may have its own mesh parameter.
In the present paper we study the Poisson equation and show that the proposed formulation is stable without assumptions on the relative sizes of the mesh parameters. In particular, we prove optimal order a priori error estimates as well as optimal order estimates of the condition number. Throughout the analysis, we trace the dependence of the number of intersecting meshes. Numerical examples are included to illustrate the stability of the method.
△ Less
Submitted 9 September, 2020; v1 submitted 17 April, 2018;
originally announced April 2018.
-
A New Least Squares Stabilized Nitsche Method for Cut Isogeometric Analysis
Authors:
Daniel Elfverson,
Mats G. Larson,
Karl Larsson
Abstract:
We derive a new stabilized symmetric Nitsche method for enforcement of Dirichlet boundary conditions for elliptic problems of second order in cut isogeometric analysis (CutIGA). We consider $C^1$ splines and stabilize the standard Nitsche method by adding certain elementwise least squares terms in the vicinity of the Dirichlet boundary and an additional term on the boundary which involves the tang…
▽ More
We derive a new stabilized symmetric Nitsche method for enforcement of Dirichlet boundary conditions for elliptic problems of second order in cut isogeometric analysis (CutIGA). We consider $C^1$ splines and stabilize the standard Nitsche method by adding certain elementwise least squares terms in the vicinity of the Dirichlet boundary and an additional term on the boundary which involves the tangential gradient. We show coercivity with respect to the energy norm for functions in $H^2(Ω)$ and optimal order a priori error estimates in the energy and $L^2$ norms. To obtain a well posed linear system of equations we combine our formulation with basis function removal which essentially eliminates basis functions with sufficiently small intersection with $Ω$. The upshot of the formulation is that only elementwise stabilization is added in contrast to standard procedures based on ghost penalty and related techniques and that the stabilization is consistent. In our numerical experiments we see that the method works remarkably well in even extreme cut situations using a Nitsche parameter of moderate size.
△ Less
Submitted 1 February, 2019; v1 submitted 16 April, 2018;
originally announced April 2018.
-
MultiMesh Finite Element Methods: Solving PDEs on Multiple Intersecting Meshes
Authors:
August Johansson,
Benjamin Kehlet,
Mats G. Larson,
Anders Logg
Abstract:
We present a new framework for expressing finite element methods on multiple intersecting meshes: multimesh finite element methods. The framework enables the use of separate meshes to discretize parts of a computational domain that are naturally separate; such as the components of an engine, the domains of a multiphysics problem, or solid bodies interacting under the influence of forces from surro…
▽ More
We present a new framework for expressing finite element methods on multiple intersecting meshes: multimesh finite element methods. The framework enables the use of separate meshes to discretize parts of a computational domain that are naturally separate; such as the components of an engine, the domains of a multiphysics problem, or solid bodies interacting under the influence of forces from surrounding fluids or other physical fields. Such multimesh finite element methods are particularly well suited to problems in which the computational domain undergoes large deformations as a result of the relative motion of the separate components of a multi-body system. In the present paper, we formulate the multimesh finite element method for the Poisson equation. Numerical examples demonstrate the optimal order convergence, the numerical robustness of the formulation and implementation in the face of thin intersections and rounding errors, as well as the applicability of the methodology.
In the accompanying paper~\cite{mmfem-2}, we analyze the proposed method and prove optimal order convergence and stability.
△ Less
Submitted 26 August, 2018; v1 submitted 23 March, 2018;
originally announced March 2018.
-
A simple embedded discrete fracture-matrix model for a coupled flow and transport problem in porous media
Authors:
Lars H. Odsæter,
Trond Kvamsdal,
Mats G. Larson
Abstract:
Accurate simulation of fluid flow and transport in fractured porous media is a key challenge in subsurface reservoir engineering. Due to the high ratio between its length and width, fractures can be modeled as lower dimensional interfaces embedded in the porous rock. We apply a recently developed embedded finite element method (EFEM) for the Darcy problem. This method allows for general fracture g…
▽ More
Accurate simulation of fluid flow and transport in fractured porous media is a key challenge in subsurface reservoir engineering. Due to the high ratio between its length and width, fractures can be modeled as lower dimensional interfaces embedded in the porous rock. We apply a recently developed embedded finite element method (EFEM) for the Darcy problem. This method allows for general fracture geometry, and the fractures may cut the finite element mesh arbitrarily. We present here a velocity model for EFEM and couple the Darcy problem to a transport problem for a passive solute. The main novelties of this work is a locally conservative velocity approximation derived from the EFEM solution, and the development of a lowest order upwind finite volume method for the transport problem. This numerical model is compatible with EFEM in the sense that the same computational mesh may be applied, so that we retain the same flexibility with respect to fracture geometry and meshing. Hence, our coupled solution strategy represents a simple approach in terms of formulation, implementation and meshing. We demonstrate our model by some numerical examples on both synthetic and realistic problems, including a benchmark study for single-phase flow. Despite the simplicity of the method, the results are promising.
△ Less
Submitted 9 March, 2018;
originally announced March 2018.
-
A cut finite element method with boundary value correction for the incompressible Stokes' equations
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
We design a cut finite element method for the incompressible Stokes equations on curved domains. The cut finite element method allows for the domain boundary to cut through the elements of the computational mesh in a very general fashion. To further facilitate the implementation we propose to use a piecewise affine discrete domain even if the physical domain has curved boundary. Dirichlet boundary…
▽ More
We design a cut finite element method for the incompressible Stokes equations on curved domains. The cut finite element method allows for the domain boundary to cut through the elements of the computational mesh in a very general fashion. To further facilitate the implementation we propose to use a piecewise affine discrete domain even if the physical domain has curved boundary. Dirichlet boundary conditions are imposed using Nitsche's method on the discrete boundary and the effect of the curved physical boundary is accounted for using the boundary value correction technique introduced for cut finite element methods in Burman, Hansbo, Larson, 'A cut finite element method with boundary value correction', Math. Comp. 87(310):633--657, 2018.
△ Less
Submitted 23 January, 2018;
originally announced January 2018.
-
Cut Finite Elements for Convection in Fractured Domains
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson,
Karl Larsson
Abstract:
We develop a cut finite element method (CutFEM) for the convection problem in a so called fractured domain which is a union of manifolds of different dimensions such that a $d$ dimensional component always resides on the boundary of a $d+1$ dimensional component. This type of domain can for instance be used to model porous media with embedded fractures that may intersect. The convection problem ca…
▽ More
We develop a cut finite element method (CutFEM) for the convection problem in a so called fractured domain which is a union of manifolds of different dimensions such that a $d$ dimensional component always resides on the boundary of a $d+1$ dimensional component. This type of domain can for instance be used to model porous media with embedded fractures that may intersect. The convection problem can be formulated in a compact form suitable for analysis using natural abstract directional derivative and divergence operators. The cut finite element method is based on using a fixed background mesh that covers the domain and the manifolds are allowed to cut through a fixed background mesh in an arbitrary way. We consider a simple method based on continuous piecewise linear elements together with weak enforcement of the coupling conditions and stabilization. We prove a priori error estimates and present illustrating numerical examples.
△ Less
Submitted 30 May, 2018; v1 submitted 18 January, 2018;
originally announced January 2018.
-
CutIGA with Basis Function Removal
Authors:
Daniel Elfverson,
Mats G. Larson,
Karl Larsson
Abstract:
We consider a cut isogeometric method, where the boundary of the domain is allowed to cut through the background mesh in an arbitrary fashion for a second order elliptic model problem. In order to stabilize the method on the cut boundary we remove basis functions which have small intersection with the computational domain. We determine criteria on the intersection which guarantee that the order of…
▽ More
We consider a cut isogeometric method, where the boundary of the domain is allowed to cut through the background mesh in an arbitrary fashion for a second order elliptic model problem. In order to stabilize the method on the cut boundary we remove basis functions which have small intersection with the computational domain. We determine criteria on the intersection which guarantee that the order of convergence in the energy norm is not affected by the removal. The higher order regularity of the B-spline basis functions leads to improved bounds compared to standard Lagrange elements.
△ Less
Submitted 2 January, 2018;
originally announced January 2018.
-
Primal dual mixed finite element methods for the elliptic Cauchy problem
Authors:
Erik Burman,
Mats. G. Larson,
Lauri Oksanen
Abstract:
We consider primal-dual mixed finite element methods for the solution of the elliptic Cauchy problem, or other related data assimilation problems. The method has a local conservation property. We derive a priori error estimates using known conditional stability estimates and determine the minimal amount of weakly consistent stabilization and Tikhonov regularization that yields optimal convergence…
▽ More
We consider primal-dual mixed finite element methods for the solution of the elliptic Cauchy problem, or other related data assimilation problems. The method has a local conservation property. We derive a priori error estimates using known conditional stability estimates and determine the minimal amount of weakly consistent stabilization and Tikhonov regularization that yields optimal convergence for smooth exact solutions. The effect of perturbations in data is also accounted for. A reduced version of the method, obtained by choosing a special stabilization of the dual variable, can be viewed as a variant of the least squares mixed finite element method introduced by Dardé, Hannukainen and Hyvönen in \emph{An {$H\sb {\sf{div}}$}-based mixed quasi-reversibility method for solving elliptic {C}auchy problems}, SIAM J. Numer. Anal., 51(4) 2013. The main difference is that our choice of regularization does not depend on auxiliary parameters, the mesh size being the only asymptotic parameter. Finally, we show that the reduced method can be used for defect correction iteration to determine the solution of the full method. The theory is illustrated by some numerical examples.
△ Less
Submitted 29 December, 2017;
originally announced December 2017.
-
Augmented Lagrangian and Galerkin least squares methods for membrane contact
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
In this paper, we propose a stabilised finite element method for the numerical solution of contact between a small deformation elastic membrane and a rigid obstacle. We limit ourselves to friction--free contact, but the formulation is readily extendable to more complex situations.
In this paper, we propose a stabilised finite element method for the numerical solution of contact between a small deformation elastic membrane and a rigid obstacle. We limit ourselves to friction--free contact, but the formulation is readily extendable to more complex situations.
△ Less
Submitted 13 November, 2017;
originally announced November 2017.
-
Stabilization of High Order Cut Finite Element Methods on Surfaces
Authors:
Mats G. Larson,
Sara Zahedi
Abstract:
We develop and analyze a stabilization term for cut finite element approximations of an elliptic second order partial differential equation on a surface embedded in $\mathbb{R}^d$. The new stabilization term combines properly scaled normal derivatives at the surface together with control of the jump in the normal derivatives across faces and provides control of the variation of the finite element…
▽ More
We develop and analyze a stabilization term for cut finite element approximations of an elliptic second order partial differential equation on a surface embedded in $\mathbb{R}^d$. The new stabilization term combines properly scaled normal derivatives at the surface together with control of the jump in the normal derivatives across faces and provides control of the variation of the finite element solution on the active three dimensional elements that intersect the surface. We show that the condition number of the stiffness matrix is $O(h^{-2})$, where $h$ is the mesh parameter. The stabilization term works for linear as well as for higher-order elements and the derivation of its stabilizing properties is quite straightforward, which we illustrate by discussing the extension of the analysis to general $n$-dimensional smooth manifolds embedded in $\mathbb{R}^d$, with codimension $d-n$. We also formulate properties of a general stabilization term that are sufficient to prove optimal scaling of the condition number and optimal error estimates in energy- and $L^2$-norm. We finally present numerical studies confirming our theoretical results.
△ Less
Submitted 22 August, 2018; v1 submitted 9 October, 2017;
originally announced October 2017.
-
A Simple Finite Element Method for Elliptic Bulk Problems with Embedded Surfaces
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
In this paper we develop a simple finite element method for simulation of embedded layers of high permeability in a matrix of lower permeability using a basic model of Darcy flow in embedded cracks. The cracks are allowed to cut through the mesh in arbitrary fashion and we take the flow in the crack into account by superposition. The fact that we use continuous elements leads to suboptimal converg…
▽ More
In this paper we develop a simple finite element method for simulation of embedded layers of high permeability in a matrix of lower permeability using a basic model of Darcy flow in embedded cracks. The cracks are allowed to cut through the mesh in arbitrary fashion and we take the flow in the crack into account by superposition. The fact that we use continuous elements leads to suboptimal convergence due to the loss of regularity across the crack. We therefore refine the mesh in the vicinity of the crack in order to recover optimal order convergence in terms of the global mesh parameter. The proper degree of refinement is determined based on an a priori error estimate and can thus be performed before the actual finite element computation is started. Numerical examples showing this effect and confirming the theoretical results are provided. The approach is easy to implement and beneficial for rapid assessment of the effect of crack orientation and may for example be used in an optimization loop.
△ Less
Submitted 4 September, 2017;
originally announced September 2017.
-
A simple approach for finite element simulation of reinforced plates
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
We present a new approach for adding Bernoulli beam reinforcements to Kirchhoff plates. The plate is discretised using a continuous/discontinuous finite element method based on standard continuous piecewise polynomial finite element spaces. The beams are discretised by the CutFEM technique of letting the basis functions of the plate represent also the beams which are allowed to pass through the pl…
▽ More
We present a new approach for adding Bernoulli beam reinforcements to Kirchhoff plates. The plate is discretised using a continuous/discontinuous finite element method based on standard continuous piecewise polynomial finite element spaces. The beams are discretised by the CutFEM technique of letting the basis functions of the plate represent also the beams which are allowed to pass through the plate elements. This allows for a fast and easy way of assessing where the plate should be supported, for instance, in an optimization loop.
△ Less
Submitted 5 June, 2017;
originally announced June 2017.
-
A Nitsche Method for Elliptic Problems on Composite Surfaces
Authors:
Peter Hansbo,
Tobias Jonsson,
Mats G. Larson,
Karl Larsson
Abstract:
We develop a finite element method for elliptic partial differential equations on so called composite surfaces that are built up out of a finite number of surfaces with boundaries that fit together nicely in the sense that the intersection between any two surfaces in the composite surface is either empty, a point, or a curve segment, called an interface curve. Note that several surfaces can inters…
▽ More
We develop a finite element method for elliptic partial differential equations on so called composite surfaces that are built up out of a finite number of surfaces with boundaries that fit together nicely in the sense that the intersection between any two surfaces in the composite surface is either empty, a point, or a curve segment, called an interface curve. Note that several surfaces can intersect along the same interface curve. On the composite surface we consider a broken finite element space which consists of a continuous finite element space at each subsurface without continuity requirements across the interface curves. We derive a Nitsche type formulation in this general setting and by assuming only that a certain inverse inequality and an approximation property hold we can derive stability and error estimates in the case when the geometry is exactly represented. We discuss several different realizations, including so called cut meshes, of the method. Finally, we present numerical examples.
△ Less
Submitted 25 May, 2017; v1 submitted 23 May, 2017;
originally announced May 2017.
-
Cut Finite Element Methods for Elliptic Problems on Multipatch Parametric Surfaces
Authors:
Tobias Jonsson,
Mats G. Larson,
Karl Larsson
Abstract:
We develop a finite element method for the Laplace--Beltrami operator on a surface described by a set of patchwise parametrizations. The patches provide a partition of the surface and each patch is the image by a diffeomorphism of a subdomain of the unit square which is bounded by a number of smooth trim curves. A patchwise tensor product mesh is constructed by using a structured mesh in the refer…
▽ More
We develop a finite element method for the Laplace--Beltrami operator on a surface described by a set of patchwise parametrizations. The patches provide a partition of the surface and each patch is the image by a diffeomorphism of a subdomain of the unit square which is bounded by a number of smooth trim curves. A patchwise tensor product mesh is constructed by using a structured mesh in the reference domain. Since the patches are trimmed we obtain cut elements in the vicinity of the interfaces. We discretize the Laplace--Beltrami operator using a cut finite element method that utilizes Nitsche's method to enforce continuity at the interfaces and a consistent stabilization term to handle the cut elements. Several quantities in the method are conveniently computed in the reference domain where the mappings impose a Riemannian metric. We derive a priori estimates in the energy and $L^2$ norm and also present several numerical examples confirming our theoretical results.
△ Less
Submitted 22 June, 2017; v1 submitted 21 March, 2017;
originally announced March 2017.
-
Finite element procedures for computing normals and mean curvature on triangulated surfaces and their use for mesh refinement
Authors:
Mirza Cenanovic,
Peter Hansbo,
Mats G. Larson
Abstract:
In this paper we consider finite element approaches to computing the mean curvature vector and normal at the vertices of piecewise linear triangulated surfaces. In particular, we adopt a stabilization technique which allows for first order $L^2$-convergence of the mean curvature vector and apply this stabilization technique also to the computation of continuous, recovered, normals using $L^2$-proj…
▽ More
In this paper we consider finite element approaches to computing the mean curvature vector and normal at the vertices of piecewise linear triangulated surfaces. In particular, we adopt a stabilization technique which allows for first order $L^2$-convergence of the mean curvature vector and apply this stabilization technique also to the computation of continuous, recovered, normals using $L^2$-projections of the piecewise constant face normals. Finally, we use our projected normals to define an adaptive mesh refinement approach to geometry resolution where we also employ spline techniques to reconstruct the surface before refinement. We compare or results to previously proposed approaches.
△ Less
Submitted 16 March, 2017;
originally announced March 2017.
-
Cut Finite Element Methods for Linear Elasticity Problems
Authors:
Peter Hansbo,
Mats G. Larson,
Karl Larsson
Abstract:
We formulate a cut finite element method for linear elasticity based on higher order elements on a fixed background mesh. Key to the method is a stabilization term which provides control of the jumps in the derivatives of the finite element functions across faces in the vicinity of the boundary. We then develop the basic theoretical results including error estimates and estimates of the condition…
▽ More
We formulate a cut finite element method for linear elasticity based on higher order elements on a fixed background mesh. Key to the method is a stabilization term which provides control of the jumps in the derivatives of the finite element functions across faces in the vicinity of the boundary. We then develop the basic theoretical results including error estimates and estimates of the condition number of the mass and stiffness matrices. We apply the method to the standard displacement problem, the frequency response problem, and the eigenvalue problem. We present several numerical examples including studies of thin bending dominated structures relevant for engineering applications. Finally, we develop a cut finite element method for fibre reinforced materials where the fibres are modeled as a superposition of a truss and a Euler-Bernoulli beam. The beam model leads to a fourth order problem which we discretize using the restriction of the bulk finite element space to the fibre together with a continuous/discontinuous finite element formulation. Here the bulk material stabilizes the problem and it is not necessary to add additional stabilization terms.
△ Less
Submitted 13 March, 2017;
originally announced March 2017.
-
Continuous/Discontinuous Finite Element Modelling of Kirchhoff Plate Structures in $\mathbb{R}^3$ Using Tangential Differential Calculus
Authors:
Peter Hansbo,
Mats G. Larson
Abstract:
We employ surface differential calculus to derive models for Kirchhoff plates including in-plane membrane deformations. We also extend our formulation to structures of plates. For solving the resulting set of partial differential equations, we employ a finite element method based on elements that are continuous for the displacements and discontinuous for the rotations, using $C^0$-elements for the…
▽ More
We employ surface differential calculus to derive models for Kirchhoff plates including in-plane membrane deformations. We also extend our formulation to structures of plates. For solving the resulting set of partial differential equations, we employ a finite element method based on elements that are continuous for the displacements and discontinuous for the rotations, using $C^0$-elements for the discretisation of the plate as well as for the membrane deformations. Key to the formulation of the method is a convenient definition of jumps and averages of forms that are $d$-linear in terms of the element edge normals.
△ Less
Submitted 14 February, 2017; v1 submitted 13 February, 2017;
originally announced February 2017.
-
A Stabilized Cut Finite Element Method for the Darcy Problem on Surfaces
Authors:
Peter Hansbo,
Mats G. Larson,
Andre Massing
Abstract:
We develop a cut finite element method for the Darcy problem on surfaces. The cut finite element method is based on embedding the surface in a three dimensional finite element mesh and using finite element spaces defined on the three dimensional mesh as trial and test functions. Since we consider a partial differential equation on a surface, the resulting discrete weak problem might be severely il…
▽ More
We develop a cut finite element method for the Darcy problem on surfaces. The cut finite element method is based on embedding the surface in a three dimensional finite element mesh and using finite element spaces defined on the three dimensional mesh as trial and test functions. Since we consider a partial differential equation on a surface, the resulting discrete weak problem might be severely ill conditioned. We propose a full gradient and a normal gradient based stabilization computed on the background mesh to render the proposed formulation stable and well conditioned irrespective of the surface positioning within the mesh. Our formulation extends and simplifies the Masud-Hughes stabilized primal mixed formulation of the Darcy surface problem proposed in [28] on fitted triangulated surfaces. The tangential condition on the velocity and the pressure gradient is enforced only weakly, avoiding the need for any tangential projection. The presented numerical analysis accounts for different polynomial orders for the velocity, pressure, and geometry approximation which are corroborated by numerical experiments. In particular, we demonstrate both theoretically and through numerical results that the normal gradient stabilized variant results in a high order scheme.
△ Less
Submitted 17 January, 2017;
originally announced January 2017.
-
Shape Optimization Using the Cut Finite Element Method
Authors:
Erik Burman,
Daniel Elfverson,
Peter Hansbo,
Mats G. Larson,
Karl Larsson
Abstract:
We present a cut finite element method for shape optimization in the case of linear elasticity. The elastic domain is defined by a level-set function, and the evolution of the domain is obtained by moving the level-set along a velocity field using a transport equation. The velocity field is the largest decreasing direction of the shape derivative that satisfies a certain regularity requirement and…
▽ More
We present a cut finite element method for shape optimization in the case of linear elasticity. The elastic domain is defined by a level-set function, and the evolution of the domain is obtained by moving the level-set along a velocity field using a transport equation. The velocity field is the largest decreasing direction of the shape derivative that satisfies a certain regularity requirement and the computation of the shape derivative is based on a volume formulation. Using the cut finite element method no re--meshing is required when updating the domain and we may also use higher order finite element approximations. To obtain a stable method, stabilization terms are added in the vicinity of the cut elements at the boundary, which provides control of the variation of the solution in the vicinity of the boundary. We implement and illustrate the performance of the method in the two--dimensional case, considering both triangular and quadrilateral meshes as well as finite element spaces of different order.
△ Less
Submitted 17 November, 2016;
originally announced November 2016.