Hierarchical analysis-suitable T-splines: Formulation, Bézier extraction, and application as an adaptive basis for isogeometric analysis
In this paper hierarchical analysis-suitable T-splines (HASTS) are developed. The resulting spaces are a superset of both analysis-suitable T-splines and hierarchical B-splines. The additional flexibility provided by the hierarchy of T-spline spaces results in simple, highly localized refinement algorithms which can be utilized in a design or analysis context. A detailed theoretical formulation is presented including a proof of local linear independence for analysis-suitable T-splines, a requisite theoretical ingredient for HASTS. Bézier extraction is extended to HASTS simplifying the implementation of HASTS in existing finite element codes. The behavior of a simple HASTS refinement algorithm is compared to the local refinement algorithm for analysis-suitable T-splines demonstrating the superior efficiency and locality of the HASTS algorithm. Finally, HASTS are utilized as a basis for adaptive isogeometric analysis.
keywords:isogeometric analysis, hierarchical splines, adaptive mesh refinement, T-splines
In this work, a hierarchical extension of analysis-suitable T-splines is developed and utilized in the context of isogeometric design and analysis. We call this new spline description hierarchical analysis-suitable T-splines (HASTS). The class of HASTS is a strict superset of both analysis-suitable T-splines LiZhSeHuSc10 (); ScLiSeHu10 (); BeBuChSa12 (); BeBuSaVa12 (); LiScSe12 () and hierarchical B-splines FoBa88 (); SchDeScEvBoRaHu12 (); VuGiJuSi11 (); ScThEv13 (); GiJu13 (); GiJuSp12 ().
T-splines, introduced in the CAD community SeZhBaNa03 (), are a generalization of non-uniform rational B-splines (NURBS) which address fundamental limitations in NURBS-based design. For example, a T-spline can model a complicated design as a single, watertight geometry and are also locally refinable SeCaFiNoZhLy04 (); ScLiSeHu10 (). Since their advent they have emerged as an important technology across multiple disciplines and can be found in several major commercial CAD products TSManual12 (); Autodesk360 ().
Isogeometric analysis was introduced in HuCoBa04 () and described in detail in Cottrell:2009rp (). The isogeometric paradigm is simple: use the smooth spline basis that defines the geometry as the basis for analysis. As a result, exact geometry is introduced into the analysis, the smooth basis can be leveraged by the analysis EvBaBaHu09 (); HuEvRe13 (); CoHuRe07 (), and new innovative approaches to model design WaZhScHu11 (); LiZhHuScSe14 (), analysis SchDeScEvBoRaHu12 (); ScSiEvLiBoHuSe12 (); ScWuBl12 (); BeBaDeHsScHuBe09 (), optimization Wall08 (), and adaptivity Bazilevs2009 (); DoJuSi09 (); ScThEv13 (); ScThEv13 () are made possible. The use of T-splines as a basis for isogeometric analysis (IGA) has gained widespread attention across a number of application areas Bazilevs2009 (); ScBoHu10 (); ScLiSeHu10 (); Verhoosel:2010vn (); Verhoosel:2010ly (); BoScLaHuVe11 (); BeBaDeHsScHuBe09 (); SchDeScEvBoRaHu12 (); ScSiEvLiBoHuSe12 (); SiScTaThLi14 (); DiLoScWrTaZa13 (); HoReVeBo14 (); BaHsSc12 (); BuSaVa12 (); GiKoPoKaBeGeScHu14 (). Particular focus has been placed on the use of T-spline local refinement in an analysis context ScLiSeHu10 (); BoScLaHuVe11 (); Verhoosel:2010ly (); Verhoosel:2010vn ().
In the context of CAD, where a designer interacts directly with the geometry, T-spline local refinement is most useful if confined to a single level. In other words, all local refinement is done on one control mesh and all control points have similar influence on the shape of the surface. In this way, the geometric behavior of the surface is easily controlled through the manipulation of control points before and after refinement. In the context of analysis, however, where not all control points need to have a geometric interpretation, the single level restriction can be relaxed. This hierarchical point of view has important advantages:
Hierarchical local refinement remains completely local. Single level T-spline local refinement always entails a degree of nonlocal control point propagation ScLiSeHu10 ().
Hierarchical refinement and coarsening operations use a fixed control mesh which simplifies algorithmic developments, especially for parallel computations. Single level local refinement requires expensive mesh manipulation and modification operations.
Hierarchies of finite-dimensional subspaces are the natural setting for many optimized iterative solvers and preconditioning techniques for large-scale linear systems.
Initial investigations employing hierarchical B-spline refinement in the context of IGA have demonstrated the promise of the hierarchical approach KuVeZeBr14 (); SchDeScEvBoRaHu12 (); SchEvReScHu13 (); GrKrSch02 (); ScRa11 ().
HASTS inherit the design strengths of T-splines without the single level restriction. In this way, a complex T-spline design can be encapsulated in the first level of the hierarchy while higher levels can be leveraged to develop adaptive multiresolution schemes which are smooth, highly localized, geometrically exact, and appropriate for the analysis task at hand. We feel that this provides the appropriate mathematical foundation for the development of integrated isogeometric design and analysis methodologies for demanding applications in science and engineering. Note that, in this paper, we restrict our theoretical developments to HASTS defined over four-sided domains. However, extending HASTS to domains of arbitrary topological genus should be straightforward in the context of the recently introduced spline forest ScThEv13 ().
We note that in addition to T-splines, hierarchical B-splines, and NURBS a number of alternative spline technologies have been proposed as a basis for IGA with varying strengths and weaknesses. Truncated hierarchical B-splines (THB-splines) GiJuSp12 (); KiGiJu14 (); GiJuSp13 (); BeKiBrChMoOhKi14 () are a modification of hierarchical B-splines FoBa88 (); VuGiJuSi11 (); SchEvReScHu13 () which possess a partition of unity and enhanced numerical conditioning. B-spline forests ScThEv13 () are a generalization of hierarchical B-splines to surfaces and volumes of arbitrary topological genus. Polynomial splines over hierarchical T-meshes (PHT-splines) htspline2 (); htspline3 (); htspline4 (); htspline5 (), modified T-splines KaChDe13 (), and locally refined splines (LR-splines) DoLyPe13 (); Br13 () are closely related to T-splines with varying levels of smoothness and approaches to local refinement. Generalized B-splines MaPeSa11 (); CoMaPeSa10 () and T-splines BrBeChOhKi14 () enhance a piecewise polynomial spline basis by including non-polynomial functions, typically trigonometric or hyperbolic functions. Generalized splines permit the exact representation of conic sections without resorting to rational functions. Generalized splines can also be used to represent solution features with known non-polynomial characteristics exactly in certain circumstances.
1.1 Structure and content of the paper
In Section 2 the T-mesh is described and appropriate notational conventions are introduced. Analysis-suitable T-splines are then described in Section 3. The local linear independence of analysis-suitable T-splines is established in Section 4. Hierarchical analysis-suitable T-splines are then defined in Section 5. In preparation for their use in design and analysis a Bézier extraction framework is introduced in Section 6. HASTS are then utilized as a basis for isogeometric analysis in Section 7. In Section 8 we draw conclusions. We note that the paper has been written so the proof of local linear independence in Section 4 is self-contained and can be skipped if the reader is not interested in the detailed theory of analysis-suitable T-splines.
2 The T-mesh
The T-mesh is used to define the topological structure of the associated T-spline space. In other words, the T-mesh defines the basis functions and their relationship to one another. We closely follow the notational conventions introduced in LiScSe12 (); BeBuSaVa12 (); BeBuChSa12 ().
A T-mesh in two dimensions is a rectangular partition of such that all vertices have integer coordinates. All cells are rectangular, non-overlapping, and open. An edge is a horizontal or vertical line segment between vertices which does not intersect any cell. The valence of a vertex is the number of edges coincident to that vertex. Since all cells are assumed rectangular, only valence three (i.e., T-junction) or four is allowed for all vertices . The sets of horizontal and vertical coordinates in the T-mesh are denoted by and . The horizontal and vertical skeletons, and , of a T-mesh are the union of all horizontal and vertical edges, respectively, and associated vertices. The entire skeleton is denoted by .
We split into an active region and a frame region such that and , and where and are polynomial degrees. Note that both and are closed. Further, all T-meshes considered in this work are admissible as described in LiScSe12 (), a mild restriction always adopted in practice. The notation will indicate that can be created by adding vertices and edges to . Fig. 1 shows an example T-mesh.
2.1 Analysis-suitable T-meshes
Analysis-suitable T-splines (ASTS) were introduced in LiZhSeHuSc10 (). The analysis-suitability of a T-spline is dictated by the structure of the underlying T-mesh. We define face and edge extensions to be closed line segments that originate at T-junctions. For example, to define a horizontal face extension we trace out a horizontal line by moving in the direction of the missing edge until vertical edges or vertices are intersected. To define an edge extension we trace out a horizontal line by moving in the direction opposite the face extension until vertical edges or vertices are intersected. A T-junction extension includes both the face and edge extensions. Since extensions are defined as closed line segments they may intersect at their end points. An extended T-mesh, , is the T-mesh formed by adding the T-junction extensions to . The collection of rectangular cells in is denoted by . We say a T-mesh is analysis-suitable if no horizontal T-junction extension intersects a vertical T-junction extension. Face and edge extensions (along with analysis-suitability) are illustrated in Figure 2.
Anchors are used in the construction of T-spline blending functions. For an analysis-suitable T-mesh the anchors are located only in the active region and are defined as follows:
if and are odd the anchors are vertices. It is written as or equivalently .
if is even and is odd the anchors are horizontal edges. It is written as .
if is odd and is even the anchors are vertical edges. It is written as .
if and are even the anchors are cells. It is written as .
The set of all anchors is denoted by . The set of anchors for varying values of and are shown in Figure 3.
3 Analysis-suitable T-splines
Analysis-suitable T-splines form a useful subset of T-splines. ASTS maintain the important mathematical properties of the NURBS basis while providing an efficient and highly localized refinement capability. Several important properties of ASTS have been proven:
The blending functions are linearly independent for any choice of knots LiZhSeHuSc10 ().
The basis constitutes a partition of unity LiScSe12 ().
Each basis function is non-negative.
They can be generalized to arbitrary degree BeBuSaVa12 ().
An affine transformation of an analysis-suitable T-spline is obtained by applying the transformation to the control points. We refer to this as affine covariance. This implies that all “patch tests” (see Hug00 ()) are satisfied a priori.
They obey the convex hull property.
Optimal approximation LiScSe12 ().
The important properties of ASTS emanate directly from the topological properties of the underlying analysis-suitable T-mesh and resulting set of T-spline basis functions constructed from it.
3.1 T-spline basis functions, spaces, and geometry
Given a parametric domain we define global horizontal and vertical open knot vectors and , respectively. In other words,
As a result, every T-mesh vertex has the parametric representation . For reasons that will become apparent, we refer to a cell with positive parametric area as a Bézier element. The parametric domain of a Bézier element is denoted by . The set of all Bézier elements in a T-mesh is denoted by .
For each anchor we construct horizontal and vertical local index vectors and made up of increasing (but not necessarily consecutive) indices in and , respectively. Note that for odd, , and for even, . Similar relationships hold for . The procedure for determining local index vectors is shown in Fig. 4 for various polynomial degrees. To clarify this procedure we describe the anchors and associated local index vectors. In Figure 4a, and and thus the example anchor is the cell . The horizontal local index vector is and the vertical local index vector is We observe that subset of the vertical skeleton located at index does not span the entire height of the anchor cell, hence it is not included in the horizontal local index vector; similarly since subset of the horizontal skeleton located at index does not span the entire width of the cell it not included in the vertical local index vector. In Figure 4b, and thus the example anchor is the vertical edge . The horizontal local index vector is and the vertical local index vector is Similar to the prior example the subset of the vertical skeleton at index does not span the entire height of the anchor edge, hence it is not included in the horizontal local index vector. Figure 4c shows the case where and thus the example anchor is the horizontal edge . The horizontal local index vector is and the vertical local index vector is In the last case, shown in Figure 4d, and thus the example anchor is the vertex . The horizontal local index vector is and the vertical local index vector is
The T-spline blending function is given by
where and are B-spline basis functions associated with the local knot vectors and .
We define to be the set of all basis functions associated with a T-mesh. Given a weight for each a rational T-spline basis function can be written as
where is called a weight function. For clarity we will often suppress the dependence on the polynomial degrees and write the basis function as . Figure 5 shows several T-spline basis functions plotted in the parametric domain .
An ASTS space, denoted by , is the span of the blending functions in constructed from an analysis-suitable T-mesh. Given vector valued control points, , or , the geometry of a T-spline can be written as
4 Local linear independence of analysis-suitable T-splines
The local linear independence of analysis-suitable T-splines is an important theoretical result in its own right and is critical for our definition of hierarchical analysis-suitable T-splines in Section 5. The (global) linear independence of analysis-suitable T-splines was first shown in LiZhSeHuSc10 (). Local linear independence is a stronger result than global linear independence and is the notion of linear independence enjoyed by standard finite element bases. Local linear independence implies that the finite element basis is linearly independent over every element domain. For smooth locally-refined bases this element-level notion of linear independence may be lost.
Given a knot vector and degree we can recursively define B-spline basis functions as follows:
The de Boor algorithm Farin99 () provides a standard method for evaluating a B-spline, although other possibilities exist Se10 (); Ra89 (). A B-spline basis function can also be denoted by . Following Schu93 (), there exist dual functionals such that where the Kronecker delta is zero when and one otherwise. We have the following two results:
Suppose . Then
For any function , if , then for we have that
Given an analysis-suitable T-mesh and associated basis functions the set of functionals form a dual basis. Specifically, we have that
where and are dual basis functions corresponding to univariate B-splines Schu93 () with local knot vectors and , respectively.
4.2 Proof of local linear independence
Let be a cell from an analysis-suitable extended T-mesh with vertices and . If the basis function anchored at with local index vectors and is non-zero over then there must either exist an integer , , such that and or an integer , , such that and .
Suppose the lemma is false, then there exists at least one cell in that violates the condition. We denote this cell by be . Since violates the condition there exists at least one corner of that lies in where and . Without loss of generality we may assume that and . We have following three cases:
The corner is a vertex of the original T-mesh. This violates Lemma 3.2(a) in BeBuSaVa12 ().
The corner is the result of the intersection of two perpendicular T-junction extensions. This violates the assumption that the T-mesh is analysis-suitable.
The corner is the result of the intersection of one T-junction extension and a T-mesh edge. Without loss of generality, we assume the edge is a horizontal edge and the T-junction extension is vertical and is associated with a T-junction . As the edge cannot intersect the vertical line , it must terminate in a T-junction . Examining the extensions associated with and we find that they must intersect. This violates the assumption that the T-mesh is analysis-suitable.
Hence, such cell cannot exist. ∎
The basis for an analysis-suitable T-spline is locally linearly independent.
Let an arbitrary Bézier element, , from an analysis-suitable T-mesh be given. We denote the element vertices by and . Let be the set of anchors whose corresponding basis functions are non-zero over . Assume that for all . By Lemma 4.4, since is non-zero over , there exists either a , , such that and or an , , such that and By Lemma 4.3
Since for or for we have by Lemma 4.2 that
which completes the proof. ∎
5 Hierarchical analysis-suitable T-splines
A hierarchical T-spline space is constructed from a finite sequence of nested ASTS spaces, , , and bounded open index domains, , which define the nested domains for the hierarchy. Two important theoretical results for ASTS will be used in the construction of hierarchical analysis-suitable T-splines:
Given two analysis-suitable T-meshes with non-overlapping T-junction extensions, and , if , then .
Analysis-suitable T-splines are locally linear independent.
We note that to accommodate overlapping T-junction extensions requires a minor generalization of Theorem 5.1 which is not reproduced here to maintain clarity of exposition. For a complete description of the underlying theory we refer the interested reader to LiScSe12 (). The local linear independence of ASTS is proven in Section 4.
5.1 Sequences of analysis-suitable T-meshes
We construct a sequence of analysis-suitable T-meshes such that , , as follows:
Create from by subdividing each cell in into four congruent cells.
Extend T-junctions in until it is analysis-suitable and .
(b) Create from by subdividing Bézier elements.
(c) Extend T-junctions until is analysis-suitable and .
5.2 Hierarchical T-spline spaces
The hierarchical analysis-suitable T-spline basis can be constructed recursively in a manner analogous to that used for hierarchical B-splines VuGiJuSi11 ():
Recursively construct from by setting
We denote the number of functions in by . We call the space spanned by the functions in a hierarchical analysis-suitable T-spline space and denote it by . To make the ideas concrete a univariate hierarchical spline space is shown in Figure 7.
The linear independence of the functions in follows immediately from the definition of hierarchical T-splines and the local linear independence of ASTS (see Section 4).
The functions in the hierarchical basis are linearly independent.
See Lemma 2 in VuGiJuSi11 () ∎
Given , a sequence of hierarchical analysis-suitable T-spline bases, , .
See Lemma 3 in VuGiJuSi11 () ∎
6 Bézier extraction of hierarchical analysis-suitable T-splines
The Bézier extraction framework ScBoHu10 (); Borden:2010nx (); ScThEv13 () can be extended to HASTS in a straightforward fashion. Using Bézier extraction, the spline hierarchy is collapsed onto a single level finite element mesh which can then be processed by standard finite element codes without any explicit knowledge of HASTS algorithms or data structures.
6.1 Bernstein basis functions
The univariate Bernstein basis functions are written as
where and the binomial coefficient , . In CAGD, the Bernstein polynomials are usually defined over the unit interval , but in finite element analysis the biunit interval is preferred to take advantage of the usual domains for Gauss quadrature. The univariate Bernstein basis has the following properties:
Partition of unity.
Figure 8 shows the Bernstein basis for polynomial degrees .
We construct a bivariate Bernstein basis function of degree by where , , and , as the tensor product of univariate basis functions
6.2 The geometry of a hierarchical representation
In a single level T-spline, basis functions and control points have a one-to-one relationship and each control point influences the geometry in a similar manner. In a hierarchical context it is common to only associate control points with the functions in . This is the convention adopted in this paper. Note that by construction every blending function in can be written in terms of basis functions in (see Lemma 5.4). We call the functions in geometric blending functions. We use to denote the number of geometric blending functions.
Given vector valued control points, , or , and weights , the geometry of a hierarchical representation can be written as
where , is used to index the geometric blending functions, and is the weight function. The decoupling of geometry from the basis functions in is an additional complexity unique to hierarchical representations which is elegantly addressed via Bézier extraction.
6.3 Bézier Elements
The set of Bézier elements underlying a hierarchical T-spline are determined recursively in a manner similar to the basis. We denote the set of Bézier elements in a hierarchy by . We construct as follows:
Recursively construct from by setting
We denote the number of Bézier elements in by .
6.4 Element localization
Using standard techniques ScBoHu10 (); Borden:2010nx () it is possible to determine the set of functions in which are nonzero over any element in . This gives rise to a standard element connectivity map which, given an element index and local function index , returns a global function index . In other words . The reader is referred to Hug00 () for additional details on common approaches to finite element localization and the array. Note that can indicate an anchor or a global function index.
We write a rational hierarchical T-spline basis function, restricted to element , as
where and is the element weight function restricted to element . The element geometric map is the restriction of to element .
6.5 Bézier extraction
To present the basic ideas, Bézier extraction for a B-spline curve is shown graphically in Figure 9. Bézier extraction constructs a linear transformation defined by a matrix referred to as the extraction operator. The extraction operator maps a Bernstein polynomial basis defined on Bézier elements to the global spline basis. The transpose of the extraction operator maps the control points of the spline to the Bézier control points.
Each hierarchical basis function supported by element can be written in Bézier form as
where the dependence of the Bernstein polynomial on the polynomial degrees and has been suppressed for clarity. The overbar will be used to denote a quantity written in terms of the Bernstein basis defined over the element domain . The Bézier coefficients are computed using standard knot insertion techniques ScBoHu10 (). We denote the vector of hierarchical basis functions supported by element by and the vector of Bernstein basis functions by . We then have that
where is the element extraction operator (see ScThEv13 ()). In other words, the element extraction operator is composed of the Bézier coefficients .
We write the element weight function as
where is the number of geometric basis functions which are non-zero over element . We may also write the rational hierarchical basis functions as
Finally, the element geometric map can be written as
7 Computational Results
We illustrate the use of hierarchical T-splines in the context of isogeometric analysis. We consider problems that highlight the unique attributes of both hierarchical refinement and T-splines. The examples used are inspired by those found in ScLiSeHu10 (); VuGiJuSi11 (); ScThEv13 ().
7.1 A comparison between ASTS and HASTS local refinement
We compare local refinement of ASTS to local refinement of HASTS. When working with ASTS all refinement is performed on a single level whereas when working with HASTS this constraint is relaxed. For additional algorithmic details on local refinement of ASTS see ScLiSeHu10 (). We locally refine the T-spline ship hull design shown in Fig. 10 using both methods. The geometry is constructed using the Autodesk T-spline plugin for Rhino3d tspline_rhino (). T-splines are popular in ship hull design because an entire hull can be modeled by a single watertight surface with a minimal number of control points SeSe10 (). T-junctions can be used to efficiently model local features. Note that the initial T-spline of the hull contains just control points and Bézier elements.
We restrict the refinement region to the locations detailed in Fig. 11. It is assumed that the original design is too coarse to be used as a basis for analysis and additional resolution is required in the rectangular region followed by highly localized refinements along the region corresponding to the curve.
The HASTS refinement algorithm is based on the algorithm presented in ScThEv13 () for spline forests. The algorithm is element-based, meaning refinement is driven by the subdivision of Bézier elements. The hierarchical basis is then reextracted into the new hierarchical T-mesh topology to generate the new set of Bézier elements. A detailed description of the underlying algorithms, in the context of HASTS, will be postponed to a future publication. Figure 12 shows three HASTS local refinements along the curve shown in Figure 11. The elements are colored according to their level, , in the hierarchy. Note that no nonlocal propagation of local refinement occurs for HASTS. Only those elements specified for refinement are subdivided. This is possible due to the relaxation of the single level constraint inherent in ASTS. The refinements form a nested sequence of -continuous hierarchical analysis-suitable T-spline spaces. The geometry of the hull is exactly preserved during refinement. The final HASTS is composed of Bézier elements and basis functions. However, only geometric blending functions and control points are used to define the hull geometry.
As a comparison, Figure 13 shows the results of ASTS local refinement using the algorithm from ScLiSeHu10 (). The top figure shows the control points added during local refinement (black dots) along the curve. The region selected for refinement is shown in red. Observe the propagation of the control points away from the selected refinement region. The bottom figure shows the resulting Bézier elements after refinement. Superfluous control points and elements are added just to satisfy the single level constraint inherent in the definition of ASTS.
7.2 HASTS as an adaptive basis
We now consider HASTS as an adaptive basis for isogeometric analysis. We choose as a benchmark the advection skew to the mesh problem shown in Figure 14. This problem is advection dominated, with diffusivity of . Along the external boundary, the boundary conditions are selected such that sharp interior and boundary layers are present in the solution. In this case, degrees.
7.2.1 Problem Statement
Let be a bounded region in and assume has a piecewise smooth boundary . Let denote a general point in , and let the temperature at a point be denoted by . Given Dirichlet boundary data, , the steady-state advection-diffusion boundary value problem consists of finding the temperature such that
where and are the spatially varying solenoidal velocity vector and symmetric, positive-definite, diffusivity tensor, respectively. Note that in this paper we define where is a positive constant called the diffusivity coefficient. We employ SUPG BroHug82 () with a standard definition for the element stabilization parameter, .
7.2.2 A residual based error estimator
Note that this error estimator underestimates the error for diffusion-dominated flows but is adequate for the advection-dominated benchmark presented in this paper. Using standard techniques Hug00 () we use the element scaling
where are the mesh size distributions for and , respectively, and is the order of convergence of the method. The element size, , is the square root of the element area. We flag elements for refinement if . The adaptive process is repeated until a specified convergence tolerance is attained or a maximum number of hierarchical levels are introduced.
We solve the problem with biquadratic and bicubic hierarchical T-splines. The initial T-mesh for both cases is shown in Figure 15. Note that the initial T-mesh is locally refined to accommodate the presence of sharp boundary layers in the solution. Note that this refinement is not hierarchical. We have found that judiciously performing local refinement of the first level of a T-spline hierarchy to accommodate geometric features or boundary conditions leads to smaller hierarchies and more efficient solution procedures.
During each adaptive step the error is assessed as described in Section 7.2.2, elements are flagged for refinement and subdivided, and a new hierarchical basis is then extracted into the new hierarchical T-mesh topology. This generates a refined set of Bézier elements. The sequence of Bézier mesh refinements is shown for both the biquadratic and bicubic case in Figures 16–17. The sequence of biquadratic refinements form a nested sequence of -continuous HASTS spaces, whereas the sequence of bicubic refinements form a nested sequence of -continuous HASTS spaces. Note that fewer elements are required for convergence as the smoothness and order of the basis increases ScThEv13 ().
To illustrate the structure and distribution of the hierarchical basis the Greville abscissae ScSiEvLiBoHuSe12 () are plotted in Figures 18–20. Note that a linear parameterization was employed for all meshes and the level zero control points and blending functions define the geometry. The dots are scaled according to their level in the hierarchy; a larger dot denotes a lower level. The sequence of solutions are shown in Figures 21–22.
(a) Initial biquadratic Bézier mesh
(b) Initial bicubic Bézier mesh
(c) Second biquadratic Bézier mesh
(d) Second bicubic Bézier mesh
(e) Third biquadratic Bézier mesh
(f) Third bicubic Bézier mesh
(a) Fourth biquadratic Bézier mesh
(b) Fourth bicubic Bézier mesh
(c) Fifth biquadratic Bézier mesh
(d) Fifth bicubic Bézier mesh
(e) Sixth biquadratic Bézier mesh
(f) Sixth bicubic Bézier mesh
(a) Initial quadratic Greville abscissae
(b) Initial bicubic Greville abscissae
(c) Greville abscissae for second biquadratic mesh
(d) Greville abscissae for second bicubic mesh
(a) Greville abscissae for third biquadratic mesh
(b) Greville abscissae for third bicubic mesh
(c) Greville abscissae for fourth biquadratic mesh
(d) Greville abscissae for fourth bicubic mesh
(a) Greville abscissae for fifth biquadratic mesh
(b) Greville abscissae for fifth bicubic mesh
(c) Greville abscissae for sixth biquadratic mesh
(d) Greville abscissae for sixth bicubic mesh
(f) Solution 5
We have presented hierarchical analysis-suitable T-splines which is a superset of both analysis-suitable T-splines and hierarchical B-splines. We have also developed the necessary theoretical formulation of HASTS including a proof of the local linear independence of analysis-suitable T-splines. We presented a simple algorithm for the creation of nested T-spline spaces and also extended Bézier extraction to HASTS. We then demonstrated the potential of HASTS by comparing HASTS to a local refinement algorithm for T-splines which demonstrated the improved efficiency and locality of using a hierarchical approach. We also demonstrated the use of HASTS in the context of isogeometric analysis by solving the benchmark static skew advection problem.
In future work we will provide a detailed description of the underlying algorithms to perform hierarchical refinement in the context of HASTS. We will also consider hierarchical and refinements of T-splines. Finally, we intend to extend the definition of spline forests in ScThEv13 () to the T-spline regime. This will allow us to accommodate smooth interfaces and also interface directly with commercial T-spline products.
- (1) X. Li, J. Zheng, T. W. Sederberg, T. J. R. Hughes, M. A. Scott, On linear independence of T-spline blending functions, Computer Aided Geometric Design 29 (2012) 63 – 76.
- (2) M. A. Scott, X. Li, T. W. Sederberg, T. J. R. Hughes, Local refinement of analysis-suitable T-splines, Computer Methods in Applied Mechanics and Engineering 213 (2012) 206 – 222.
- (3) L. Beirão da Veiga, A. Buffa, D. Cho, G. Sangalli, Analysis-suitable T-splines are dual-compatible, Computer Methods in Applied Mechanics and Engineering 249 – 252 (2012) 42 – 51.
- (4) L. B. da Veiga, A. Buffa, G. Sangalli, R. Vázquez, Analysis-suitable T-splines of arbitrary degree: definition, linear independence, and approximation properties, Mathematical Models and Methods in Applied Sciences 23 (11) (2013) in press.
- (5) X. Li, M. A. Scott, Analysis-suitable T-splines: characterization, refineability, and approximation, Mathematical Models and Methods in Applied Science 24 (06) (2014) 1141–1164.
- (6) D. R. Forsey, R. H. Bartels, Hierarchical B-spline refinement, ACM SIGGRAPH Computer Graphics 22 (4) (1988) 205–212.
- (7) D. Schillinger, L. Dedé, M. A. Scott, J. A. Evans, M. J. Borden, E. Rank, T. J. R. Hughes, An isogeometric design-through-analysis methodology based on adaptive hierarchical refinement of NURBS, immersed boundary methods, and T-spline CAD surfaces, Computer Methods in Applied Mechanics and Engineering 249 – 252 (2012) 116 – 150.
- (8) A. Vuong, C. Giannelli, B. Jüttler, B. Simeon, A hierarchical approach to adaptive local refinement in isogeometric analysis, Computer Methods in Applied Mechanics and Engineering 200 (49 – 52) (2011) 3554 – 3567.
- (9) M. A. Scott, D. T. Thomas, E. J. Evans, Isogeometric spline forests, Computer Methods in Applied Mechanics and Engineering 269 (2014) 222 – 264.
- (10) C. Giannelli, B. Jüttler, Bases and dimensions of bivariate hierarchical tensor-product splines, Journal of Computational and Applied Mathematics 239 (0) (2013) 162 – 178.
- (11) C. Giannelli, B. Jüttler, H. Speleers, THB–splines: The truncated basis for hierarchical splines, Computer Aided Geometric Design 29 (7) (2012) 485 – 498.
- (12) T. W. Sederberg, J. Zheng, A. Bakenov, A. Nasri, T-splines and T-NURCCs, ACM Trans. Graph. 22 (2003) 477–484.
- (13) T. W. Sederberg, D. L. Cardon, G. T. Finnigan, N. S. North, J. Zheng, T. Lyche, T-spline simplification and local refinement, ACM Trans. Graph. 23 (2004) 276–283.
- (14) Autodesk, Autodesk T-Splines Plug-in for Rhino user manual, Autodesk (2012).
- (15) Autodesk, Inc., Autodesk Fusion 360, Autodesk, Inc. (2014).
- (16) T. J. R. Hughes, J. A. Cottrell, Y. Bazilevs, Isogeometric analysis: CAD, finite elements, NURBS, exact geometry, and mesh refinement, Computer Methods in Applied Mechanics and Engineering 194 (2005) 4135–4195.
- (17) J. A. Cottrell, T. J. R. Hughes, Y. Bazilevs, Isogeometric analysis: Toward Integration of CAD and FEA, Wiley, Chichester, 2009.
- (18) J. A. Evans, Y. Bazilevs, I. Babuška, T. J. R. Hughes, n-widths, sup-infs, and optimality ratios for the k-version of the isogeometric finite element method, Computer Methods in Applied Mechanics and Engineering 198 (21-26) (2009) 1726–1741.
- (19) T. J. R. Hughes, J. A. Evans, A. Reali, Finite element and NURBS approximations of eigenvalue, boundary-value, and initial-value problems, Computer Methods in Applied Mechanics and Engineering 272 (2014) 290 – 320.
- (20) J. A. Cottrell, T. J. R. Hughes, A. Reali, Studies of refinement and continuity in isogeometric analysis, Computer Methods in Applied Mechanics and Engineering 196 (2007) 4160–4183.
- (21) W. Wang, Y. Zhang, M. A. Scott, T. J. R. Hughes, Converting an unstructured quadrilateral mesh to a standard T-spline surface, Computational Mechanics, 48 (2011) 477 – 498.
- (22) L. Liu, Y. Zhang, T. J. R. Hughes, M. A. Scott, T. W. Sederberg, Volumetric T-spline construction using boolean operations, in: J. Sarrate, M. Staten (Eds.), Proceedings of the 22nd International Meshing Roundtable, Springer International Publishing, 2014, pp. 405–424.
- (23) M. A. Scott, R. N. Simpson, J. A. Evans, S. Lipton, S. P. A. Bordas, T. J. R. Hughes, T. W. Sederberg, Isogeometric boundary element analysis using unstructured T-splines, Computer Methods in Applied Mechanics and Engineering 254 (2013) 197 – 221.
- (24) R. Schmidt, R. Wüchner, K.-U. Bletzinger, Isogeometric analysis of trimmed NURBS geometries, Computer Methods in Applied Mechanics and Engineering 241–244 (2012) 93 – 111.
- (25) D. J. Benson, Y. Bazilevs, E. De Luycker, M. C. Hsu, M. A. Scott, T. J. R. Hughes, T. Belytschko, A generalized finite element formulation for arbitrary basis functions: From isogeometric analysis to XFEM, International Journal for Numerical Methods in Engineering 83 (2010) 765–785.
- (26) W. A. Wall, M. A. Frenzel, C. Cyron, Isogeometric structural shape optimization, Computer Methods in Applied Mechanics and Engineering 197 (2008) 2976–2988.
- (27) Y. Bazilevs, V. M. Calo, J. A. Cottrell, J. A. Evans, T. J. R. Hughes, S. Lipton, M. A. Scott, T. W. Sederberg, Isogeometric analysis using T-splines, Computer Methods in Applied Mechanics and Engineering 199 (5-8) (2010) 229–263.
- (28) M. Dörfel, B. Jüttler, B. Simeon, Adaptive isogeometric analysis by local h-refinement with T-splines, Computer Methods in Applied Mechanics and Engineering 199 (5–8) (2009) 264–275.
- (29) M. A. Scott, M. J. Borden, C. V. Verhoosel, T. W. Sederberg, T. J. R. Hughes, Isogeometric Finite Element Data Structures based on Bézier Extraction of T-splines, International Journal for Numerical Methods in Engineering, 88 (2011) 126 – 156.
- (30) C. V. Verhoosel, M. A. Scott, T. J. R. Hughes, R. de Borst, An isogeometric analysis approach to gradient damage models, International Journal for Numerical Methods in Engineering, 86 (2011) 115–134.
- (31) C. V. Verhoosel, M. A. Scott, R. de Borst, T. J. R. Hughes, An isogeometric approach to cohesive zone modeling, International Journal for Numerical Methods in Engineering, 87 (2011) 336 – 360.
- (32) M. J. Borden, M. A. Scott, C. V. Verhoosel, C. M. Landis, T. J. R. Hughes, A phase-field description of dynamic brittle fracture, Computer Methods in Applied Mechanics and Engineering 217 (2012) 77 – 95.
- (33) R. Simpson, M. Scott, M. Taus, D. Thomas, H. Lian, Acoustic isogeometric boundary element analysis, Computer Methods in Applied Mechanics and Engineering 269 (2014) 265–290.
- (34) R. Dimitri, L. D. Lorenzis, M. A. Scott, P. Wriggers, R. Taylor, G. Zavarise, Isogeometric large deformation frictionless contact using T-splines, Computer methods in applied mechanics and engineering 269 (2014) 394 – 414.
- (35) S. Hosseini, J. J. Remmers, C. V. Verhoosel, R. de Borst, An isogeometric continuum shell element for non-linear analysis, Computer Methods in Applied Mechanics and Engineering 271 (2014) 1 – 22.
- (36) Y. Bazilevs, M. C. Hsu, M. A. Scott, Isogeometric fluid-structure interaction analysis with emphasis on non-matching discretizations, and with application to wind turbines, Computer Methods in Applied Mechanics and Engineering 249 - 252 (2012) 28 – 41.
- (37) A. Buffa, G. Sangalli, R. Vázquez, Isogeometric methods for computational electromagnetics: B-spline and T-spline discretizations, Journal of Computational Physics 257, Part B (2014) 1291 – 1320.
- (38) A. I. Ginnis, K. V. Kostas, C. G. Politis, P. D. Kaklis, K. A. Belibassakis, T. P. Gerostathis, M. A. Scott, T. J. R. Hughes, Isogeometric boundary-element analysis for the wave-resistance problem using T-splines, Computer Methods in Applied Mechanics and Engineering submitted.
- (39) G. Kuru, C. Verhoosel, K. van der Zee, E. van Brummelen, Goal-adaptive isogeometric analysis with hierarchical splines, Computer Methods in Applied Mechanics and Engineering 270 (2014) 270 – 292.
- (40) D. Schillinger, J. A. Evans, A. Reali, M. A. Scott, T. J. R. Hughes, Isogeometric collocation: Cost comparison with Galerkin methods and extension to adaptive hierarchical NURBS discretizations, Computer Methods in Applied Mechanics and Engineering 267 (2013) 170 – 232.
- (41) E. Grinspun, P. Krysl, P. Schröder, CHARMS: a simple framework for adaptive simulation, ACM Transactions on Graphics 21 (3) (2002) 281–290.
- (42) D. Schillinger, E. Rank, An unfitted -adaptive finite element method based on hierarchical B-splines for interface problems of complex geometry, Computer Methods in Applied Mechanics and Engineering 200 (47–48) (2011) 3358 – 3380.
- (43) G. Kiss, C. Giannelli, B. Jüttler, Algorithms and data structures for truncated hierarchical B-splines, in: M. Floater, T. Lyche, M.-L. Mazure, K. Mørken, L. Schumaker (Eds.), Mathematical Methods for Curves and Surfaces, Vol. 8177 of Lecture Notes in Computer Science, Springer Berlin Heidelberg, 2014, pp. 304–323.
- (44) C. Giannelli, B. Jüttler, H. Speleers, Strongly stable bases for adaptively refined multilevel spline spaces, Advances in Computational Mathematics (2013) 1–32.
- (45) D. Berdinsky, T. wan Kim, C. Bracco, D. Cho, B. Mourrain, M. Oh, S. Kiatpanichgij, Dimensions and bases of hierarchical tensor-product splines, Journal of Computational and Applied Mathematics 257 (2014) 86 – 104.
- (46) J. Deng, F. Chen, X. Li, C. Hu, W. Tong, Z. Yang, Y. Feng, Polynomial splines over hierarchical T-meshes, Graphical Models 74 (2008) 76–86.
- (47) X. Li, J. Deng, F. Chen, Surface modeling with polynomial splines over hierarchical T-meshes, The Visual Computer 23 (2007) 1027–1033.
- (48) X. Li, J. Deng, F. Chen, Polynomial splines over general T-meshes, The Visual Computer 26 (2010) 277–286.
- (49) X. Li, J. Deng, F. Chen, The dimension of spline spaces over 3d hierarchical T-meshes, Journal of Information and Computational Science 3 (2006) 487–501.
- (50) H. Kang, F. Chen, J. Deng, Modified T-splines, Computer Aided Geometric Design 30 (9) (2013) 827 – 843.
- (51) T. Dokken, T. Lyche, K. F. Pettersen, Polynomial splines over locally refined box-partitions, Computer Aided Geometric Design 30 (3) (2013) 331–356.
- (52) A. Bressan, Some properties of LR-splines, Computer Aided Geometric Design 30 (8) (2013) 778 – 794.
- (53) C. Manni, F. Pelosi, M. L. Sampoli, Generalized B-splines as a tool in isogeometric analysis, Computer Methods in Applied Mechanics and Engineering 200 (5–8) (2011) 867 – 881.
- (54) P. Costantini, C. Manni, F. Pelosi, M. L. Sampoli, Quasi-interpolation in isogeometric analysis based on generalized B-splines, Computer Aided Geometric Design 27 (8) (2010) 656–668.
- (55) C. Bracco, D. Berdinsky, D. Cho, M. Oh, T. wan Kim, Trigonometric generalized T-splines, Computer Methods in Applied Mechanics and Engineering 268 (2014) 540 – 556.
- (56) T. J. R. Hughes, The Finite Element Method: Linear Static and Dynamic Finite Element Analysis, Dover Publications, Mineola, NY, 2000.
- (57) G. E. Farin, NURBS Curves and Surfaces: from Projective Geometry to Practical Use, A. K. Peters, Ltd., Natick, MA, 1999.
- (58) T. W. Sederberg, Computer Aided Geometric Design: Course Notes, http://www.tsplines.com/educationportal.html (2010).
- (59) L. Ramshaw, Blossoms are polar forms, Computer Aided Geometric Design 6 (4) (1989) 323–358.
- (60) L. Schumaker, Spline functions: basic theory, Krieger, 1993.
- (61) M. J. Borden, M. A. Scott, J. A. Evans, T. J. R. Hughes, Isogeometric finite element data structures based on Bézier extraction of NURBS, International Journal for Numerical Methods in Engineering, 87 (2011) 15 – 47.
- (62) Autodesk, Inc., http://www.tsplines.com/rhino/ (2011).
- (63) M. T. Sederberg, T. W. Sederberg, T-splines: A technology for marine design with minimal control points, http://www.tsplines.com/technicalpapers.html (2010).
- (64) A. N. Brooks, T. J. R. Hughes, Streamline upwind / Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations, Computer Methods in Applied Mechanics and Engineering 32 (1982) 199–259.
- (65) T. J. R. Hughes, Multiscale phenomena: Green’s functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods, Computer Methods in Applied Mechanics and Engineering 127 (1995) 387–401.
- (66) M. G. Larson, A. Målqvist, Adaptive variational multiscale methods based on a posteriori error estimation: Duality techniques for elliptic problems, in: B. Engquist, O. Runborg, P. Lötstedt (Eds.), Multiscale Methods in Science and Engineering, Vol. 44 of Lecture Notes in Computational Science and Engineering, Springer Berlin Heidelberg, 2005, pp. 181–193.
- (67) G. Hauke, M. H. Doweidar, M. Miana, The multiscale approach to error estimation and adaptivity, Computer Methods in Applied Mechanics and Engineering 195 (13–16) (2006) 1573 – 1593.
- (68) G. Hauke, M. H. Doweidar, S. Fuentes, Mesh adaptivity for the transport equation led by variational multiscale error estimators, International Journal for Numerical Methods in Fluids 69 (12) (2012) 1835–1850.