Efficient extraction of hierarchical B-Splines for local refinement and coarsening of Isogeometric Analysis

Abstract The main motivation for hierarchical B-Splines in Isogeometric Analysis is to perform a local refinement that is globally not a tensor product. Compared to standard knot insertion, this allows to increase refinement locality. Moreover, the implementation of hierarchical refinement on existing codes can be facilitated by Bezier extraction. This approach was initially developed for unrefined patches and then extended to local refinement. In case of B-Spline patches that are not locally refined, extraction and projection algorithms can be optimized by leveraging the tensor structure. However, these optimizations are not straightforwardly applicable to hierarchical B-Splines. In this contribution, we show an approach where all extraction operations are formulated to work directly on the univariate extraction operators, without explicitly computing and storing the full tensor product. In particular, we consider the extraction of functions, and projection of degrees of freedom for refinement and coarsening. A slightly modified Bezier projection is proposed to fully exploit the developed algorithms. We finally compare our approach to an explicit computation of the extraction operator in both a linear transient and a nonlinear example.

[1]  Andreas Keller,et al.  Finite Element Approximation with Hierarchical B-Splines , 2014, Curves and Surfaces.

[2]  C. Loan.  The ubiquitous Kronecker product , 2000 .

[3]  Tom Lyche,et al.  Polynomial splines over locally refined box-partitions , 2013, Comput. Aided Geom. Des..

[4]  F. Cirak,et al.  A subdivision-based implementation of the hierarchical b-spline finite element method , 2013 .

[5]  F. Brezzi,et al.  An "immersed" finite element method based on a locally anisotropic remeshing for the incompressible Stokes problem , 2015 .

[6]  Rafael Vázquez,et al.  Algorithms for the implementation of adaptive isogeometric methods using hierarchical splines , 2016 .

[7]  T. Hughes,et al.  Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows , 2007 .

[8]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

[9]  D. J. Benson,et al.  Patient-specific isogeometric structural analysis of aortic valve closure , 2015 .

[10]  Trond Kvamsdal,et al.  On the similarities and differences between Classical Hierarchical, Truncated Hierarchical and LR B-splines , 2015 .

[11]  Alexander Düster,et al.  Numerical integration for nonlinear problems of the finite cell method using an adaptive scheme based on moment fitting , 2019, Comput. Math. Appl..

[12]  Yuri Bazilevs,et al.  3D simulation of wind turbine rotors at full scale. Part II: Fluid–structure interaction modeling with composite blades , 2011 .

[13]  Markus Kästner,et al.  Bézier extraction and adaptive refinement of truncated hierarchical NURBS , 2016 .

[14]  Alessandro Reali,et al.  Suitably graded THB-spline refinement and coarsening: Towards an adaptive isogeometric analysis of additive manufacturing processes , 2018, Computer Methods in Applied Mechanics and Engineering.

[15]  Wolfgang Böhm.  On the efficiency of knot insertion algorithms , 1985, Comput. Aided Geom. Des..

[16]  Paulo Fernandes,et al.  Efficient descriptor-vector multiplications in stochastic automata networks , 1998, JACM.

[17]  John A. Evans,et al.  Bézier projection: A unified approach for local projection and quadrature-free refinement and coarsening of NURBS and T-splines with particular application to isogeometric design and analysis , 2014, 1404.7155.

[18]  Carlotta Giannelli,et al.  Adaptive isogeometric methods with hierarchical splines: error estimator and convergence , 2015, 1502.00565.

[19]  D. Owen,et al.  Computational methods for plasticity : theory and applications , 2008 .

[20]  Giancarlo Sangalli,et al.  Fast formation of isogeometric Galerkin matrices by weighted quadrature , 2016, 1605.01238.

[21]  Xin Li,et al.  Hierarchical T-splines: Analysis-suitability, Bézier extraction, and application as an adaptive basis for isogeometric analysis , 2014, 1404.4346.

[22]  Hendrik Speleers,et al.  Strongly stable bases for adaptively refined multilevel spline spaces , 2014, Adv. Comput. Math..

[23]  T. Hughes,et al.  Local refinement of analysis-suitable T-splines , 2012 .

[24]  Víctor D. Fachinotti,et al.  Analytical solutions of the thermal field induced by moving double‐ellipsoidal and double‐elliptical heat sources in a semi‐infinite body , 2011 .

[25]  Ernst Rank,et al.  A hierarchical computational model for moving thermal loads and phase changes with applications to selective laser melting , 2017, Comput. Math. Appl..

[26]  Giancarlo Sangalli,et al.  Matrix-free weighted quadrature for a computationally efficient isogeometric k-method , 2017, Computer Methods in Applied Mechanics and Engineering.

[27]  Victor M. Calo,et al.  The role of continuity in residual-based variational multiscale modeling of turbulence , 2007 .

[28]  Hendrik Speleers,et al.  THB-splines: The truncated basis for hierarchical splines , 2012, Comput. Aided Geom. Des..

[29]  David R. Forsey,et al.  Hierarchical B-spline refinement , 1988, SIGGRAPH.

[30]  Rida T. Farouki,et al.  On the numerical condition of Bernstein-Bézier subdivision processes , 1990 .

[31]  T. Hughes,et al.  B¯ and F¯ projection methods for nearly incompressible linear and non-linear elasticity and plasticity using higher-order NURBS elements , 2008 .

[32]  Cv Clemens Verhoosel,et al.  A phase-field description of dynamic brittle fracture , 2012 .

[33]  Michael A. Scott,et al.  Isogeometric spline forests , 2014 .

[34]  Andrea Bressan,et al.  A Versatile Strategy for the Implementation of Adaptive Splines , 2016, MMCS.

[35]  Alessandro Reali,et al.  Multi-level Bézier extraction for hierarchical local refinement of Isogeometric Analysis , 2017 .

[36]  Michael C. H. Wu,et al.  Isogeometric Kirchhoff–Love shell formulations for general hyperelastic materials , 2015 .

[37]  Alessandro Reali,et al.  Isogeometric Analysis of Structural Vibrations , 2006 .

[38]  Thomas J. R. Hughes,et al.  Hierarchically refined and coarsened splines for moving interface problems, with particular application to phase-field models of prostate tumor growth , 2017 .

[39]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[40]  Yuri Bazilevs,et al.  Dynamic and fluid–structure interaction simulations of bioprosthetic heart valves using parametric design with T-splines and Fung-type material models , 2015, Computational mechanics.

[41]  Eid H. Doha,et al.  Integrals of Bernstein polynomials: An application for the solution of high even-order differential equations , 2011, Appl. Math. Lett..

[42]  John A. Evans,et al.  Isogeometric finite element data structures based on Bézier extraction of NURBS , 2011 .

[43]  Anh-Vu Vuong,et al.  Finite Element Concepts and Bezier Extraction in Hierarchical Isogeometric Analysis , 2014 .

[44]  Bert Jüttler,et al.  THB-splines: An effective mathematical technology for adaptive refinement in geometric design and isogeometric analysis , 2016 .

[45]  Markus Kästner,et al.  Projection and transfer operators in adaptive isogeometric analysis with hierarchical B-splines , 2018, Computer Methods in Applied Mechanics and Engineering.

[46]  Günther Greiner,et al.  Interpolating and approximating scattered 3D-data with hierarchical tensor product B-splines , 2010 .