NAMay 3, 2016
Adaptive Mesh Refinement Strategies in Isogeometric Analysis - A Computational ComparisonPaul Hennig, Markus Kästner, Philipp Morgenstern et al.
We explain four variants of an adaptive finite element method with cubic splines and compare their performance in simple elliptic model problems. The methods in comparison are Truncated Hierarchical B-splines with two different refinement strategies, T-splines with the refinement strategy introduced by Scott et al. in 2012, and T-splines with an alternative refinement strategy introduced by some of the authors. In four examples, including singular and non-singular problems of linear elasticity and the Poisson problem, the H1-errors of the discrete solutions, the number of degrees of freedom as well as sparsity patterns and condition numbers of the discretized problem are compared.
NAMay 11, 2016
Globally structured 3D Analysis-suitable T-splines: definition, linear independence and m-graded local refinementPhilipp Morgenstern
This paper addresses the linear independence of T-splines that correspond to refinements of three-dimensional tensor-product meshes. We give an abstract definition of analysis-suitability, and prove that it is equivalent to dual-compatibility, wich guarantees linear independence of the T-spline blending functions. In addition, we present a local refinement algorithm that generates analysis-suitable meshes and has linear computational complexity in terms of the number of marked and generated mesh elements.
NASep 18, 2015
Complexity of hierarchical refinement for a class of admissible mesh configurationsAnnalisa Buffa, Carlotta Giannelli, Philipp Morgenstern et al.
An adaptive isogeometric method based on $d$-variate hierarchical spline constructions can be derived by considering a refine module that preserves a certain class of admissibility between two consecutive steps of the adaptive loop [6]. In this paper we provide a complexity estimate, i.e., an estimate on how the number of mesh elements grows with respect to the number of elements that are marked for refinement by the adaptive strategy. Our estimate is in the line of the similar ones proved in the finite element context, [3,24].
NADec 13, 2014
Analysis-suitable adaptive T-mesh refinement with linear complexityPhilipp Morgenstern, Daniel Peterseim
We present an efficient adaptive refinement procedure that preserves analysis-suitability of the T-mesh, this is, the linear independence of the T-spline blending functions. We prove analysis-suitability of the overlays and boundedness of their cardinalities, nestedness of the generated T-spline spaces, and linear computational complexity of the refinement procedure in terms of the number of marked and generated mesh elements.
NAOct 10, 2014
Multiscale Partition of UnityDaniel Peterseim, Patrick Henning, Philipp Morgenstern
We introduce a new Partition of Unity Method for the numerical homogenization of elliptic partial differential equations with arbitrarily rough coefficients. We do not restrict to a particular ansatz space or the existence of a finite element mesh. The method modifies a given partition of unity such that optimal convergence is achieved independent of oscillation or discontinuities of the diffusion coefficient. The modification is based on an orthogonal decomposition of the solution space while preserving the partition of unity property. This precomputation involves the solution of independent problems on local subdomains of selectable size. We deduce quantitative error estimates for the method that account for the chosen amount of localization. Numerical experiments illustrate the high approximation properties even for 'cheap' parameter choices.