umu.sePublications
Change search
Link to record
Permanent link

Direct link
BETA
Publications (10 of 21) 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.

Place, publisher, year, edition, pages
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.

Place, publisher, year, edition, pages
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
Show others...
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.

Place, publisher, year, edition, pages
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.

Place, publisher, year, edition, pages
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
Burman, E., Elfverson, D., Hansbo, P., Larson, M. G. & Larsson, K. (2019). Hybridized CutFEM for Elliptic Interface Problems. SIAM Journal on Scientific Computing, 41(5), A3354-A3380
Open this publication in new window or tab >>Hybridized CutFEM for Elliptic Interface Problems
Show others...
2019 (English)In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 41, no 5, p. A3354-A3380Article in journal (Refereed) Published
Abstract [en]

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.

Place, publisher, year, edition, pages
Society for Industrial and Applied Mathematics, 2019
National Category
Computational Mathematics
Identifiers
urn:nbn:se:umu:diva-164715 (URN)10.1137/18M1223836 (DOI)000493897100026 ()
Funder
Swedish Foundation for Strategic Research , AM13-0029Swedish Research Council, 2013-4708Swedish Research Council, 2017-03911Swedish Research Council, 2018-05262eSSENCE - An eScience Collaboration, -
Available from: 2019-10-29 Created: 2019-10-29 Last updated: 2019-11-19Bibliographically 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
2018 (English)Report (Other academic)
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.

Publisher
p. 18
Keywords
electrodynamics, special relativity, coordinate-free calculus
National Category
Other Physics Topics
Research subject
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
2018 (English)Report (Other academic)
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  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  defines a product on the space of complex 4-vectors  so that  becomes an algebra isomorphic to  with  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 . For example, instead of the electromagnetic dyadic field  we now have the complex field variable . 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.

Publisher
p. 7
Keywords
electrodynamics, special relativity, coordinate-free calculus, spinors
National Category
Other Physics Topics
Research subject
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.

Place, publisher, year, edition, pages
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.

Place, publisher, year, edition, pages
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
Show others...
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
Organisations
Identifiers
ORCID iD: ORCID iD iconorcid.org/0000-0001-7838-1307

Search in DiVA

Show all publications