Researcher profile

Mats G. Larson

Mats G. Larson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

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

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.

preprint2022arXiv

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

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

preprint2021arXiv

CutFEM Based on Extended Finite Element Spaces

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

preprint2021arXiv

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.

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.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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.