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

Mats G. Larson

Mats G. Larson 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

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

42 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
preprint2021arXiv

The Finite Cell Method with Least Squares Stabilized Nitsche Boundary Conditions

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.

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

MultiMesh Finite Elements with Flexible Mesh Sizes

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.

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

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galerkin least squares finite element method for the obstacle problem

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiscale methods for problems with complex geometry

We propose a multiscale method for elliptic problems on complex domains, e.g. domains with cracks or complicated boundary. For local singularities this paper also offers a discrete alternative to enrichment techniques such as XFEM. We construct corrected coarse test and trail spaces which takes the fine scale features of the computational domain into account. The corrections only need to be computed in regions surrounding fine scale geometric features. We achieve linear convergence rate in energy norm for the multiscale solution. Moreover, the conditioning of the resulting matrices is not affected by the way the domain boundary cuts the coarse elements in the background mesh. The analytical findings are verified in a series of numerical experiments.

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

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 Nitsche-based cut finite element method for a fluid--structure interaction problem

We present a new composite mesh finite element method for fluid--structure interaction problems. The method is based on surrounding the structure by a boundary-fitted fluid mesh which is embedded into a fixed background fluid mesh. The embedding allows for an arbitrary overlap of the fluid meshes. The coupling between the embedded and background fluid meshes is enforced using a stabilized Nitsche formulation which allows us to establish stability and optimal order \emph{a priori} error estimates, see~\cite{MassingLarsonLoggEtAl2013}. We consider here a steady state fluid--structure interaction problem where a hyperelastic structure interacts with a viscous fluid modeled by the Stokes equations. We evaluate an iterative solution procedure based on splitting and present three-dimensional numerical examples.

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

Continuous Piecewise Linear Finite Elements for the Kirchhoff-Love Plate Equation

A family of continuous piecewise linear finite elements for thin plate problems is presented. We use standard linear interpolation of the deflection field to reconstruct a discontinuous piecewise quadratic deflection field. This allows us to use discontinuous Galerkin methods for the Kirchhoff-Love plate equation. Three example reconstructions of quadratic functions from linear interpolation triangles are presented: a reconstruction using Morley basis functions, a fully quadratic reconstruction, and a more general least squares approach to a fully quadratic reconstruction. The Morley reconstruction is shown to be equivalent to the Basic Plate Triangle. Given a condition on the reconstruction operator, a priori error estimates are proved in energy norm and $L^2$ norm. Numerical results indicate that the Morley reconstruction/Basic Plate Triangle does not converge on unstructured meshes while the fully quadratic reconstruction show optimal convergence.

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

High Order Cut Finite Element Methods for the Stokes Problem

We develop a high order cut finite element method for the Stokes problem based on general inf-sup stable finite element spaces. We focus in particular on composite meshes consisting of one mesh that overlaps another. The method is based on a Nitsche formulation of the interface condition together with a stabilization term. Starting from inf-sup stable spaces on the two meshes, we prove that the resulting composite method is indeed inf-sup stable and as a consequence optimal \emph{a~priori} error estimates hold.

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

$L^2$-error estimates for finite element approximations of boundary fluxes

We prove quasi-optimal a priori error estimates for finite element approximations of boundary normal fluxes in the $L^2$-norm. Our results are valid for a variety of different schemes for weakly enforcing Dirichlet boundary conditions including Nitsche's method, and Lagrange multiplier methods. The proof is based on an error representation formula that is derived by using a discrete dual problem with $L^2$-Dirichlet boundary data and combines a weighted discrete stability estimate for the dual problem with anisotropic interpolation estimates in the boundary zone.

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

A Posteriori Error Analysis of Component Mode Synthesis for the Frequency Response Problem

We consider the frequency response problem and derive a posteriori error estimates for the discrete error in a reduced finite element model obtained using the component mode synthesis (CMS) method. We provide estimates in a linear quantity of interest and the energy norm. The estimates reflect to what degree each CMS subspace influence the overall error in the reduced solution. This enables automatic error control through adaptive algorithms that determine suitable dimensions of each subspace. We illustrate the theoretical results by including several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A stabilized Nitsche overlapping mesh method for the Stokes problem

We develop a Nitsche-based formulation for a general class of stabilized finite element methods for the Stokes problem posed on a pair of overlapping, non-matching meshes. By ex- tending the least-squares stabilization to the overlap region, we prove that the method is stable, consistent, and optimally convergent. To avoid an ill-conditioned linear algebra system, the scheme is augmented by a least-squares term measuring the discontinuity of the solution in the overlap region of the two meshes. As a consequence, we may prove an estimate for the condition number of the resulting stiffness matrix that is independent of the location of the interface. Finally, we present numerical examples in three spatial dimensions illustrating and confirming the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Adaptive Finite Element Splitting Method for the Incompressible Navier-Stokes Equations

We present an adaptive finite element method for the incompressible Navier--Stokes equations based on a standard splitting scheme (the incremental pressure correction scheme). The presented method combines the efficiency and simplicity of a splitting method with the powerful framework offered by the finite element method for error analysis and adaptivity. An a posteriori error estimate is derived which expresses the error in a goal functional of interest as a sum of contributions from spatial discretization, time discretization and a term that measures the deviation of the splitting scheme from a pure Galerkin scheme (the computational error). Numerical examples are presented which demonstrate the performance of the adaptive algorithm and high quality efficiency indices. It is further demonstrated that the computational error of the Navier--Stokes momentum equation is linear in the size of the time step while the computational error of the continuity equation is quadratic in the size of the time step.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient implementation of finite element methods on non-matching and overlapping meshes in 3D

In recent years, a number of finite element methods have been formulated for the solution of partial differential equations on complex geometries based on non-matching or overlapping meshes. Examples of such methods include the fictitious domain method, the extended finite element method, and Nitsche's method. In all of these methods, integrals must be computed over cut cells or subsimplices which is challenging to implement, especially in three space dimensions. In this note, we address the main challenges of such an implementation and demonstrate good performance of a fully general code for automatic detection of mesh intersections and integration over cut cells and subsimplices. As a canonical example of an overlapping mesh method, we consider Nitsche's method which we apply to Poisson's equation and a linear elastic problem.

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