umu.sePublications
Change search

BETA
##### Publications (10 of 20) Show all publications
Elfverson, D., Larson, M. G. & Larsson, K. (2019). A new least squares stabilized Nitsche method for cut isogeometric analysis. Computer Methods in Applied Mechanics and Engineering, 349, 1-16
Open this publication in new window or tab >>A new least squares stabilized Nitsche method for cut isogeometric analysis
2019 (English)In: Computer Methods in Applied Mechanics and Engineering, ISSN 0045-7825, E-ISSN 1879-2138, Vol. 349, p. 1-16Article in journal (Refereed) Published
##### Abstract [en]

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 C1 splines and stabilize the standard Nitsche method by adding a 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 H2(Ω) and optimal order a priori error estimates in the energy and L2 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.

Elsevier, 2019
##### Keywords
Fictitious domain methods, Nitsche’s method, Least squares stabilization, Isogeometric analysis
##### National Category
Computational Mathematics
##### Identifiers
urn:nbn:se:umu:diva-156840 (URN)10.1016/j.cma.2019.02.011 (DOI)2-s2.0-85062154279 (Scopus ID)
##### Funder
Swedish Research Council, 2013-4708Swedish Research Council, 2017-03911Swedish Foundation for Strategic Research , AM13-0029eSSENCE - An eScience Collaboration Available from: 2019-03-01 Created: 2019-03-01 Last updated: 2019-06-13Bibliographically approved
Burman, E., Peter, H., Larson, M. G. & Larsson, K. (2019). Cut finite elements for convection in fractured domains. Computers & Fluids, 179, 726-734
Open this publication in new window or tab >>Cut finite elements for convection in fractured domains
2019 (English)In: Computers & Fluids, ISSN 0045-7930, E-ISSN 1879-0747, Vol. 179, p. 726-734Article in journal (Refereed) Published
##### Abstract [en]

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 is formulated in a compact form suitable for analysis using natural abstract directional derivative and divergence operators. The cut finite element method is posed on 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.

Elsevier, 2019
##### Keywords
Convection problems, Fractured domains, Mixed-dimensional domains, Galerkin least squares, A priori error estimates
##### National Category
Computational Mathematics
##### Identifiers
urn:nbn:se:umu:diva-158254 (URN)10.1016/j.compfluid.2018.07.022 (DOI)000467514000053 ()2-s2.0-85052134188 (Scopus ID)
##### Funder
Swedish Research Council, 2013-4708Swedish Research Council, 2017-03911Swedish Research Council, 2018-05262Swedish Foundation for Strategic Research , AM13-0029eSSENCE - An eScience Collaboration, - Available from: 2019-04-17 Created: 2019-04-17 Last updated: 2019-06-18Bibliographically approved
Burman, E., Elfverson, D., Hansbo, P., Larson, M. G. & Larsson, K. (2019). Cut topology optimization for linear elasticity with coupling to parametric nondesign domain regions. Computer Methods in Applied Mechanics and Engineering, 350, 462-479
Open this publication in new window or tab >>Cut topology optimization for linear elasticity with coupling to parametric nondesign domain regions
2019 (English)In: Computer Methods in Applied Mechanics and Engineering, ISSN 0045-7825, E-ISSN 1879-2138, Vol. 350, p. 462-479Article in journal (Refereed) Published
##### Abstract [en]

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.

Elsevier, 2019
##### Keywords
Material distribution topology optimization, Design and nondesign domain regions, Cut finite element method
##### National Category
Computational Mathematics
##### Identifiers
urn:nbn:se:umu:diva-157679 (URN)10.1016/j.cma.2019.03.016 (DOI)000468163500019 ()
##### Funder
Swedish Research Council, 2013-4708Swedish Research Council, 2017-03911Swedish Research Council, 2018-05262Swedish Foundation for Strategic Research , AM13-0029eSSENCE - An eScience Collaboration, - Available from: 2019-03-29 Created: 2019-03-29 Last updated: 2019-06-11Bibliographically approved
Jonsson, T., Larson, M. G. & Larsson, K. (2019). Graded Parametric CutFEM and CutIGA for Elliptic Boundary Value Problems in Domains with Corners. Computer Methods in Applied Mechanics and Engineering, 354, 331-350
Open this publication in new window or tab >>Graded Parametric CutFEM and CutIGA for Elliptic Boundary Value Problems in Domains with Corners
2019 (English)In: Computer Methods in Applied Mechanics and Engineering, ISSN 0045-7825, E-ISSN 1879-2138, Vol. 354, p. 331-350Article in journal (Refereed) Published
##### Abstract [en]

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 L2 norm. This is achieved by mapping to the reference domain where we employ a structured mesh.

Elsevier, 2019
##### Keywords
Corner singularities, a priori error estimates, Cut finite element method, Cut isogeometric analysis
##### National Category
Computational Mathematics
##### Identifiers
urn:nbn:se:umu:diva-159747 (URN)10.1016/j.cma.2019.05.024 (DOI)
##### Funder
Swedish Foundation for Strategic Research , AM13-0029Swedish Research Council, 2013-4708Swedish Research Council, 2017-03911eSSENCE - An eScience Collaboration, - Available from: 2019-06-05 Created: 2019-06-05 Last updated: 2019-06-12Bibliographically approved
Larsson, J. & Larsson, K. (2018). An introduction to relativistic electrodynamics: Part I: Calculus with 4-vectors and 4-dyadics.
Open this publication in new window or tab >>An introduction to relativistic electrodynamics: Part I: Calculus with 4-vectors and 4-dyadics
##### Abstract [en]

The conventional way of introducing relativity when teaching electrodynamics is to leave Gibbs' vector calculus for a more general tensor calculus. This sudden change of formalism can be quite problematic for the students and we therefore in this two-part paper consider alternate approaches. In this Part I we use a simplified tensor formalism with 4-vectors and 4-dyadics (i.e., second order tensors built by 4-vectors) but with no tensors of higher order than two. This allows for notations in good contact with the coordinate-free Gibbs' vector calculus that the students already master. Thus we use boldface notations for 4-vectors and 4-dyadics without coordinates and index algebra to formulate Lorentz transformations, Maxwell's equations, the equation of the motion of charged particles and the stress-energy conservation law. By first working with this simplified tensor formalism the students will get better prepared to learn the standard tensor calculus needed in more advanced courses.

p. 18
##### Keywords
electrodynamics, special relativity, coordinate-free calculus
##### National Category
Other Physics Topics
Physics
##### Identifiers
urn:nbn:se:umu:diva-152747 (URN)
Available from: 2018-10-22 Created: 2018-10-22 Last updated: 2018-10-25
Larsson, J. & Larsson, K. (2018). An introduction to relativistic electrodynamics: Part II: Calculus with complex 4-vectors.
Open this publication in new window or tab >>An introduction to relativistic electrodynamics: Part II: Calculus with complex 4-vectors
##### Abstract [en]

The conventional way of introducing relativity when teaching electrodynamics is to leave Gibbs' vector calculus for a more general tensor calculus. This sudden change of formalism can be quite problematic for the students and we therefore in this two-part paper consider alternate approaches. The algebra $\mathbb{C}^{2 \times 2}$ of 2-by-2 complex matrices (sometimes presented in the form of Clifford algebra or complex quaternions) may be used for spinor related formulations of special relativity and electrodynamics. In this Part II we use this algebraic structure but with notations that fits in with the formalism of Part I. Each observer $\mathbf{e}_0$ defines a product on the space of complex 4-vectors $\mathbb{M}^{\mathbb{C}} = \mathbb{M} + i\mathbb{M}$ so that $\mathbb{M}^{\mathbb{C}}$ becomes an algebra isomorphic to $\mathbb{C}^{2\times 2}$ with $\mathbf{e}_0$ as algebra unit. The spacetime geometric equations of Part I become complex (spinor related) equations where the antisymmetric 4-dyadics have been replaced by complex 3-vectors, i.e., by elements in $\mathbf{e}_0^{\perp \mathbb{C}} = \mathbf{e}_0^\perp + i\mathbf{e}_0^\perp$. For example, instead of the electromagnetic dyadic field $\mathbf{F}$ we now have the complex field variable $\mathbf{f} = \mathbf{E}/c + i\mathbf{B}$. Some linear algebra together with the formalism of Gibbs' vector calculus (trivially allowing for complex 3-vectors) is sufficient for dealing with the equations in their complex form.

p. 7
##### Keywords
electrodynamics, special relativity, coordinate-free calculus, spinors
##### National Category
Other Physics Topics
Physics
##### Identifiers
urn:nbn:se:umu:diva-152751 (URN)
Available from: 2018-10-22 Created: 2018-10-22 Last updated: 2018-10-25
Hansbo, P., Larson, M. G. & Larsson, K. (2018). Cut Finite Element Methods for Linear Elasticity Problems. In: Stéphane P. A. Bordas; Erik Burman; Mats G. Larson; Maxim A. Olshanskii (Ed.), Geometrically Unfitted Finite Element Methods and Applications: Proceedings of the UCL Workshop 2016 (pp. 25-63). Springer
Open this publication in new window or tab >>Cut Finite Element Methods for Linear Elasticity Problems
2018 (English)In: Geometrically Unfitted Finite Element Methods and Applications: Proceedings of the UCL Workshop 2016 / [ed] Stéphane P. A. Bordas; Erik Burman; Mats G. Larson; Maxim A. Olshanskii, Springer, 2018, p. 25-63Chapter in book (Refereed)
##### Abstract [en]

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.

Springer, 2018
##### Series
Lecture Notes in Computational Science and Engineering, ISSN 1439-7358, E-ISSN 2197-7100 ; 121
##### National Category
Computational Mathematics
##### Identifiers
urn:nbn:se:umu:diva-145934 (URN)10.1007/978-3-319-71431-8 (DOI)978-3-319-71430-1 (ISBN)978-3-319-71431-8 (ISBN)
Available from: 2018-03-22 Created: 2018-03-22 Last updated: 2018-06-09
Elfverson, D., Larson, M. G. & Larsson, K. (2018). CutIGA with basis function removal. Advanced Modeling and Simulation in Engineering Sciences, 5(6)
Open this publication in new window or tab >>CutIGA with basis function removal
2018 (English)In: Advanced Modeling and Simulation in Engineering Sciences, ISSN 2213-7467, Vol. 5, no 6Article in journal (Refereed) Published
##### Abstract [en]

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.

Springer, 2018
##### National Category
Computational Mathematics
##### Identifiers
urn:nbn:se:umu:diva-145669 (URN)10.1186/s40323-018-0099-2 (DOI)
Available from: 2018-03-13 Created: 2018-03-13 Last updated: 2018-06-09Bibliographically approved
Burman, E., Elfverson, D., Hansbo, P., Larson, M. & Larsson, K. (2018). Shape optimization using the cut finite element method. Computer Methods in Applied Mechanics and Engineering, 328, 242-261
Open this publication in new window or tab >>Shape optimization using the cut finite element method
2018 (English)In: Computer Methods in Applied Mechanics and Engineering, ISSN 0045-7825, E-ISSN 1879-2138, Vol. 328, p. 242-261Article in journal (Refereed) Published
##### Abstract [en]

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.

##### Place, publisher, year, edition, pages
Lausanne: Elsevier, 2018
##### Keywords
CutFEM, Shape optimization, Level-set, Fictitious domain method, Linear elasticity
##### National Category
Computational Mathematics
##### Identifiers
urn:nbn:se:umu:diva-140281 (URN)10.1016/j.cma.2017.09.005 (DOI)000416218500011 ()
##### Funder
Swedish Foundation for Strategic Research , AM13-0029Swedish Research Council, 2011-4992Swedish Research Council, 2013-4708eSSENCE - An eScience Collaboration Available from: 2017-10-04 Created: 2017-10-04 Last updated: 2018-06-09Bibliographically approved
Larsson, K. & Larson, M. G. (2017). A continuous/discontinuous Galerkin method and a priori error estimates for the biharmonic problem on surfaces. Mathematics of Computation, 86(308), 2613-2649
Open this publication in new window or tab >>A continuous/discontinuous Galerkin method and a priori error estimates for the biharmonic problem on surfaces
2017 (English)In: Mathematics of Computation, ISSN 0025-5718, E-ISSN 1088-6842, Vol. 86, no 308, p. 2613-2649Article in journal (Refereed) Published
##### Abstract [en]

We present a continuous/discontinuous Galerkin method for approximating solutions to a fourth order elliptic PDE on a surface embedded in R-3. A priori error estimates, taking both the approximation of the surface and the approximation of surface differential operators into account, are proven in a discrete energy norm and in L-2 norm. This can be seen as an extension of the formalism and method originally used by Dziuk ( 1988) for approximating solutions to the Laplace-Beltrami problem, and within this setting this is the first analysis of a surface finite element method formulated using higher order surface differential operators. Using a polygonal approximation inverted right perpendicular(h) of an implicitly defined surface inverted right perpendicular we employ continuous piecewise quadratic finite elements to approximate solutions to the biharmonic equation on inverted right perpendicular. Numerical examples on the sphere and on the torus confirm the convergence rate implied by our estimates.

##### National Category
Computational Mathematics
##### Identifiers
urn:nbn:se:umu:diva-79207 (URN)10.1090/mcom/3179 (DOI)000404567600003 ()
##### Note

Originally published in manuscript form with title [A continuous/discontinuous Galerkin method for the biharmonic problem on surfaces]

Available from: 2013-08-13 Created: 2013-08-13 Last updated: 2018-06-08Bibliographically approved
##### Identifiers
ORCID iD: orcid.org/0000-0001-7838-1307

#### Search in DiVA

Show all publications