high order geometric methods with exact conservation properties conservation laws in for example electromagnetism solid and fluid mechanics allow an exact discrete representation in terms of line surface and volume integrals in this paper we develop high order interpolants from any basis that constitutes a partition of unity which satisfy these integral relations exactly the resulting gradient curl and divergence conforming spaces have the property that the conservation laws become completely independent of the basis functions hence they are exactly satisfied at the coarsest level of discretization and on arbitrarily curved meshes as an illustration we apply our approach to b splines and compute a 2d stokes flow inside a lid driven cavity which displays amongst others a point wise divergence free velocity field

a locking free meshfree method for the simulation of shear deformable plates based on a mixed variational formulation abstract the problem of shear locking in the thin plate limit is a well known issue that must be overcome when discretising the reissner mindlin plate equations in this paper we present a shear locking free method utilising meshfree maximum entropy basis functions and rotated raviart thomas nedelec elements within a mixed variational formulation the formulation draws upon well known techniques in the finite element literature due to the inherent properties of the maximum entropy basis functions our method allows for the direct imposition of dirichlet essential boundary conditions in contrast to methods based on moving least squares basis functions we present benchmark problems that demonstrate the accuracy and performance of the proposed method

a hierarchical approach to adaptive local refinement in isogeometric analysis abstract adaptive local refinement is one of the key issues in isogeometric analysis in this article we present an adaptive local refinement technique for isogeometric analysis based on extensions of hierarchical b splines we investigate the theoretical properties of the spline space to ensure fundamental properties like linear independence and partition of unity furthermore we use concepts well established in finite element analysis to fully integrate hierarchical spline spaces into the isogeometric setting this also allows us to access a posteriori error estimation techniques numerical results for several different examples are given and they turn out to be very promising

some periodic point results in generalized metric spaces abstract we prove some fixed point and periodic point theorems for a map in generalized metric spaces an example is provided to support our result the results presented in this paper generalize several well known comparable results in the literature

local refinement of analysis suitable t splines we develop a local refinement algorithm for analysis suitable t splines which does not produce excessive propagation of control points we then demonstrate its use as an adaptive framework for isogeometric analysis analysis suitable t splines are a class of t splines which are linearly independent and form a partition of unity these properties coupled with local refinement make this class of t splines appealing as a basis for isogeometric analysis

two dimensional domain decomposition based on skeleton computation for parameterization and isogeometric analysis abstract this paper proposes a method for decomposing two dimensional domains into subdomains for parameterization and isogeometric analysis given a complex domain in a plane the skeleton of the domain is computed to guide the domain decomposition a continuous parameterization of the domain is then obtained by parameterizing each respective subdomain this parameterization method is applied with isogeometric analysis to solve numerical pdes over two dimensional domains examples are provided to demonstrate that the new parameterization method is superior to other state of the art parameterization techniques and that it performs better in isogeometric analysis

characterization of t splines with reduced continuity order on t meshes abstract the use of t splines 30 in isogeometric analysis 24 has been proposed in 5 as a tool to enhance the flexibility of isogeometric methods if t splines are a very general concept their success in isogeometric analysis relies upon some basic properties that needs to be true as e g i linear independence of blending functions and ii polynomial reproducibility at element level in this paper we study these properties for t splines of a reduced regularity order namely for t splines of degree p and regularity p 1 p 2 our results are both for odd and even degree under mild assumptions on the underlying t mesh t splines are shown to be linearly independent and the space they span is characterized in terms of piecewise polynomials on a topological extension of the t mesh also as p is odd we construct a new topological local refinement algorithm and demonstrate its locality properties through numerical examples

an iga discontinuous galerkin method on the union of overlapped patches abstract in isogeometric analysis iga the discontinuous galerkin method iga dg is often applied when the domain is composed of multiple non uniform rational b splines nurbs patches these patches are often required to be c 0 continuous to make the iga dg available and this requirement brings some problems to the cad process this paper extends the iga dg method to solve problems on the union of the overlapped patches with this method igaodg a problem on the overlapped domain is seen as multiple sub problems on the non trimmed patches and then the continuities of traces of the numerical solutions are penalized on the boundaries of non trimmed elements with the coefficients deduced from the classical dg method comparing with the iga methods which apply the cad boolean operations to deal with the overlapped domains the presented method carries out the integrations in the non trimmed patches thus no additional integration scheme needs to be considered to achieve the gaussian integration precision discussions on the presented discretization scheme and its convergence analysis show that this method is available on the domain created by the union of the matched mismatched and overlapped patches as predicted in the convergence analysis the presented method gives the optimal convergence rates in our numerical examples finally through one of its engineering applications this paper shows how to extend igaodg to solve problems that are defined on the trimmed patches

hierarchically refined and coarsened splines for moving interface problems with particular application to phase field models of prostate tumor growth abstract moving interface problems are ubiquitous in science and engineering to develop an accurate and efficient methodology for this class of problems we present algorithms for local h adaptivity of hierarchical b splines to be utilized in isogeometric analysis we extend bezier projection an efficient quadrature free local projection technique to the hierarchical setting in this case extraction operators may not be invertible to address this issue we develop a multi level reconstruction operator which maintains the locality properties of the projection we also introduce a balance parameter to control the overlap of hierarchical functions leading to improved numerical conditioning we apply our algorithms to the simulation of localized prostate cancer growth we model this disease using the phase field method and a set of diffusion reaction equations to account for the dynamics of nutrients and a key biomarker termed prostate specific antigen our results include examples on simple 2d and 3d domains and a more compelling tissue scale patient specific simulation which is run over a prostate anatomy extracted from medical images our methods for local h adaptivity efficiently capture the evolving interface between the tumor and the neighboring healthy tissue with remarkable accuracy in all cases

a fully locking free isogeometric approach for plane linear elasticity problems a stream function formulation in this paper we study plane incompressible elastic problems by means of a stream function formulation such that a divergence free displacement field can be computed from a scalar potential the numerical scheme is constructed within the framework of nurbs based isogeometric analysis and we take advantage of the high continuity guaranteed by nurbs basis functions in order to obtain the displacement field from the potential differentiation as a consequence the obtained numerical scheme is automatically locking free in the presence of the incompressibility constraint a discontinuous galerkin approach is proposed to deal with multiple mapped possibly multiply connected domains extensive numerical results are provided to show the method capabilities 2007 elsevier b v all rights reserved

numerical solution to the time dependent maxwell equations in axisymmetric singular domains the singular complement method in this paper we present a method to solve numerically the axisymmetric time dependent maxwell equations in a singular domain in math methods appl sci 25 2002 49 math methods appl sci 26 2003 861 the mathematical tools and an in depth study of the problems posed in the meridian half plane were exposed the numerical method and experiments based on this theory are now described here it is also the generalization to axisymmetric problems of the singular complement method that we developed to solve maxwell equations in 2d singular domains see c r acad sci paris t 330 2000 391 it is based on a splitting of the space of solutions in a regular subspace and a singular one derived from the singular solutions of the laplace problem numerical examples are finally given to illustrate our purpose in particular they show how the singular complement method captures the singular part of the solution

three dimensional simulation of unstable gravity driven infiltration of water into a porous medium infiltration of water in dry porous media is subject to a powerful gravity driven instability although the phenomenon of unstable infiltration is well known its description using continuum mathematical models has posed a significant challenge for several decades the classical model of water flow in the unsaturated flow the richards equation is unable to reproduce the instability here we present a computational study of a model of unsaturated flow in porous media that extends the richards equation and is capable of predicting the instability and captures the key features of gravity fingering quantitatively the extended model is based on a phase field formulation and is fourth order in space the new model poses a set of challenges for numerical discretizations such as resolution of evolving interfaces stiffness in space and time treatment of singularly perturbed equations and discretization of higher order spatial partial differential operators we develop a numerical algorithm based on isogeometric analysis a generalization of the finite element method that permits the use of globally smooth basis functions leading to a simple and efficient discretization of higher order spatial operators in variational form we illustrate the accuracy efficiency and robustness of our method with several examples in two and three dimensions in both homogeneous and strongly heterogeneous media we simulate for the first time unstable gravity driven infiltration in three dimensions and confirm that the new theory reproduces the fundamental features of water infiltration into a porous medium our results are consistent with classical experimental observations that demonstrate a transition from stable to unstable fronts depending on the infiltration flux

duality and unified analysis of discrete approximations in structural dynamics and wave propagation comparison of p method finite elements with k method nurbs we study the discretization behavior of classical finite element and nurbs approximations on problems of structural vibrations and wave propagation we find that on the basis of equal numbers of degrees of freedom and bandwidth nurbs have superior approximation properties in fact we observe that the high mode behavior of classical finite elements is divergent with the order of approximation a surprisingly negative result on the other hand nurbs offer almost spectral approximation properties and all modes converge with increasing order of approximation

guaranteed and sharp a posteriori error estimates in isogeometric analysis we present functional type a posteriori error estimates in isogeometric analysis these estimates derived on functional grounds provide guaranteed and sharp upper bounds of the exact error in the energy norm moreover since these estimates do not contain any unknown generic constants they are fully computable and thus provide quantitative information on the error by exploiting the properties of non uniform rational b splines we present efficient computation of these error estimates the numerical realization and the quality of the computed error distribution are addressed the potential and the limitations of the proposed approach are illustrated using several computational examples

a nurbs based immersed methodology for fluid structure interaction abstract we introduce an isogeometric immersed and fully implicit formulation for fluid structure interaction fsi the method focuses on viscous incompressible flows and nonlinear hyperelastic incompressible solids which are a common case in various fields such as for example biomechanics in our fsi method we utilize an eulerian mesh on the whole domain and a lagrangian mesh on the solid domain the lagrangian mesh is arbitrarily located on top of the eulerian mesh in a non conforming fashion due to the formulation of our problem based on the immersed finite element method ifem we do not need mesh update or remeshing algorithms the fluid structure interface is the boundary of the lagrangian mesh but cuts arbitrarily the eulerian mesh the generalized method is used for time discretization and nurbs based isogeometric analysis is employed for the spatial discretization on both fluid and solid domains the information transfer between the two meshes is carried out using the nurbs functions which avoids the use of the so called discretized delta functions the higher order and especially the higher continuity of nurbs functions allow us to deal with severe mesh distortion in the lagrangian mesh in comparison with classical c 0 linear piecewise functions as we prove numerically our numerical solutions attain good agreement with theoretical data for free falling objects in two and three dimensions which confirms the feasibility of our methodology

explicit trace inequalities for isogeometric analysis and parametric hexahedral finite elements we derive new trace inequalities for nurbs mapped domains in addition to sobolev type inequalities we derive discrete trace inequalities for use in nurbs based isogeometric analysis all dependencies on shape size polynomial degree and the nurbs weighting function are precisely specified in our analysis and explicit values are provided for all bounding constants appearing in our estimates as hexahedral finite elements are special cases of nurbs our results specialize to parametric hexahedral finite elements and our analysis also generalizes to t spline based isogeometric analysis we compare the bounding constants appearing in our explicit trace inequalities with numerically computed optimal bounding constants and we discuss application of our results to a laplace problem we finish this paper with a brief exploration of so called patch wise trace inequalities for isogeometric analysis

algebraic multilevel preconditioning in isogeometric analysis construction and numerical studies abstract we present algebraic multilevel iteration amli methods for isogeometric discretization of scalar second order elliptic problems the construction of coarse grid operators and hierarchical complementary operators are given moreover for a uniform mesh on a unit interval the explicit representation of b spline basis functions for a fixed mesh size h is given for p 2 3 4 and for c 0 and c p 1 continuity the presented methods show h and almost p independent convergence rates supporting numerical results for convergence factor and iterations count for amli cycles v linear w nonlinear w are provided numerical tests are performed in two dimensions on square domain and quarter annulus and in three dimensions on quarter thick ring

an isogeometric boundary element method for electromagnetic scattering with compatible b spline discretizations abstract we outline the construction of compatible b splines on 3d surfaces that satisfy the continuity requirements for electromagnetic scattering analysis with the boundary element method method of moments our approach makes use of non uniform rational b splines to represent model geometry and compatible b splines to approximate the surface current and adopts the isogeometric concept in which the basis for analysis is taken directly from cad geometry data the approach allows for high order approximations and crucially provides a direct link with cad data structures that allows for efficient design workflows after outlining the construction of div and curl conforming b splines defined over 3d surfaces we describe their use with the electric and magnetic field integral equations using a galerkin formulation we use bezier extraction to accelerate the computation of nurbs and b spline terms and employ h matrices to provide accelerated computations and memory reduction for the dense matrices that result from the boundary integral discretization the method is verified using the well known mie scattering problem posed over a perfectly electrically conducting sphere and the classic nasa almond problem finally we demonstrate the ability of the approach to handle models with complex geometry directly from cad without mesh generation

isogeometric spline forests in this paper we present isogeometric spline forests an isogeometric spline forest is a hierarchical spline representation capable of representing surfaces or volumes of arbitrarily complex geometry and topological genus spline forests can accomodate arbitrary degree and smoothness in the underlying hierarchical basis as well as non uniform knot interval configurations we describe adaptive h refinement and coarsening algorithms for isogeometric spline forests and develop a bezier extraction framework which provides a simple and efficient single level finite element description of the complex multi level unstructured hierarchical spline basis we then demonstrate the potential of spline forests as a basis for analysis in the context of transient advection diffusion problems where fully integrated adaptivity is demonstrated for the first time in an isogeometric simulation in all cases the adaptive process remains local even in the case of moving fronts and preserves exact geometry at the coarsest level of the discretization the accuracy and robustness of the approach is demonstrated in all cases

fast divergence conforming reduced basis methods for steady navier stokes flow abstract reduced basis methods rb methods or rbms form one of the most promising techniques to deliver numerical solutions of parametrized pdes in real time with reasonable accuracy for incompressible flow problems rbms based on lbb stable velocity pressure spaces do not generally inherit the stability of the underlying high fidelity model and instead additional stabilization techniques must be introduced one way of bypassing the loss of lbb stability in the rbm is to inflate the velocity space with supremizer modes this however deteriorates the performance of the rbm in the performance critical online stage as additional dofs must be introduced to retain stability while these dofs do not effectively contribute to accuracy of the rb approximation in this work we consider a velocity only rb approximation exploiting a solenoidal velocity basis the solenoidal reduced basis emerges directly from the high fidelity velocity solutions in the offline stage by means of piola transforms the solenoidality of the velocity space is retained under geometric transformations making the proposed rb method suitable also for the investigation of geometric parameters to ensure exact solenoidality of the high fidelity velocity solutions that constitute the rb we consider approximations based on divergence conforming compatible b splines we show that the velocity only rb method leads to a significant improvement in computational efficiency in the online stage and that the pressure solution can be recovered a posteriori at negligible extra cost we illustrate the solenoidal rb approach by modeling steady two dimensional navier stokes flow around a naca0015 airfoil at various angles of attack

three dimensional hierarchical mixed finite element approximations with enhanced primal variable accuracy abstract there are different possibilities of choosing balanced pairs of approximation spaces for dual flux and primal pressure variables to be used in discrete versions of the mixed finite element method for elliptic problems arising in fluid simulations three cases shall be studied and compared for discretized three dimensional formulations based on tetrahedral hexahedral and prismatic meshes the principle guiding the constructions is the property that the divergence of the dual space and the primal approximation space should coincide while keeping the same order of accuracy for the flux variable and varying the accuracy order of the primal variable there is the classic case of b d m k spaces based on tetrahedral meshes and polynomials of total degree k for the dual variable and k 1 for the primal variable showing stable simulations with optimal convergence rates of orders k 1 and k respectively another case is related to r t k and b d f m k 1 spaces for hexahedral and tetrahedral meshes respectively but holding for prismatic elements as well it gives identical approximation order k 1 for both primal and dual variables an improvement in accuracy obtained by increasing the degree of primal functions to k and by enriching the dual space with some properly chosen internal shape functions of degree k 1 while keeping degree k for the border fluxes a new type of approximation is proposed by further incrementing the order of some internal flux functions to k 2 and matching primal functions to k 1 higher than the border fluxes of degree k thus higher convergence rate of order k 2 is obtained for the primal variable using static condensation the global condensed system to be solved in all the cases has same dimension and structure which is proportional to the space dimension of the border fluxes for each element geometry illustrating results comparing the three different space configurations are presented for simulations based on hierarchical high order shape functions for h div conforming spaces which are specially constructed for affine tetrahedral hexahedral and prismatic meshes expected convergence rates are obtained for the flux pressure and divergence variables

planar domain parameterization for isogeometric analysis based on teichmuller mapping abstract given the boundary curve of a planar domain finding a parametric spline representation for the domain is called domain parameterization a good parameterization of the computational domain plays a key role in isogeometric analysis since it influences the accuracy of the subsequent analysis in this paper we propose a new approach for planar domain parameterization based on teichmuller map a special map in the class of quasi conformal map under given correspondence of four boundary curves a unique teichmuller map between a unit square and a computational domain can be obtained which guarantees a bijection map and minimizes the maximal conformality distortion we propose an efficient iterative algorithm to compute the teichmuller map based on alternating direction method of multipliers admm experimental results show that our method can produce more uniform parameterization and increase the accuracy as well as decrease the condition numbers of the stiffness matrices in isogeometric analysis than other state of the art approaches

shape optimization directly from cad an isogeometric boundary element approach using t splines abstract we develop a t spline isogeometric boundary element method igabem simpson et al 2012 scott et al 2013 simpson et al 2014 to shape sensitivity analysis and gradient based shape optimization in three dimensional linear elasticity contrary to finite element based isogeometric analysis iga approaches no parametrization of the volume is required hence the iterative optimization algorithm can be implemented directly from cad without any mesh generation or postprocessing step for returning the resulting structure to cad designers t splines also guarantee a water tight geometry without the manual geometrical repair work as with non uniform rational b splines nurbs we demonstrate the worth of the method by analysing problems with and without analytical solutions including engineering examples involving complex shapes additionally we provide all the derivations of the required sensitivities and the details pertaining to the geometries examined in the benchmarking to provide helpful reference problems for 3d shape optimization

isogeometric schwarz preconditioners for linear elasticity systems isogeometric schwarz preconditioners are constructed and analyzed for both compressible elasticity in primal formulation and almost incompressible elasticity in mixed formulation these preconditioners require the solution of local elasticity problems on overlapping subdomains forming a decomposition of the problem domain and the solution of a coarse elasticity problem associated with the subdomain coarse mesh an h analysis of the preconditioner for the primal formulation of compressible elasticity yields an optimal convergence rate bound that is scalable in the number of subdomains and is linear in the ratio between subdomain and overlap sizes extensive numerical experiments in 2d and 3d confirm this theoretical bound and show that an analogous bound holds for the mixed formulation of almost incompressible elasticity the numerical tests also show the good preconditioner performance with respect to the polynomial degree p and regularity k of the isogeometric basis functions as well as with respect to the presence of discontinuous elastic coefficients in composite materials and to domain deformation

isogeometric finite element data structures based on bezier extraction of t splines we develop finite element data structures for t splines based on bezier extraction generalizing our previous work for nurbs as in traditional finite element analysis the extracted bezier elements are defined in terms of a fixed set of polynomial basis functions the so called bernstein basis the bezier elements may be processed in the same way as in a standard finite element computer program utilizing exactly the same data processing arrays in fact only the shape function subroutine needs to be modified while all other aspects of a finite element program remain the same a byproduct of the extraction process is the element extraction operator this operator localizes the topological and global smoothness information to the element level and represents a canonical treatment of t junctions referred to as hanging nodes in finite element analysis and a fundamental feature of t splines a detailed example is presented to illustrate the ideas copyright 2011 john wiley sons ltd

an isogeometric collocation method using superconvergent points abstract we develop an iga collocation method modified by collocating at points other than the standard greville abscissae the method is related to orthogonal collocation used for solving differential equations and to the superconvergence theory therefore we refer to this method as super collocation iga sc by carefully choosing the collocation points it can be seen that the iga sc converges in the first derivative energy norms at rates similar to that of the galerkin solution this is different from the collocation at greville abscissae iga c where the convergence in energy norm for odd polynomial degrees is typically suboptimal the method is tested on 1d 2d and 3d numerical examples in which it is compared to iga c and galerkin s method iga g the comparison includes a detailed cost vs accuracy analysis which shows an improved efficiency of the proposed method in particular for odd polynomial degrees

a scalable block preconditioning strategy for divergence conforming b spline discretizations of the stokes problem abstract the recently introduced divergence conforming b spline discretizations allow the construction of smooth discrete velocity pressure pairs for viscous incompressible flows that are at the same time inf sup stable and pointwise divergence free when applied to the discretized stokes problem these spaces generate a symmetric and indefinite saddle point linear system the iterative method of choice to solve such system is the generalized minimum residual method this method lacks robustness and one remedy is to use preconditioners for linear systems of saddle point type a large family of preconditioners can be obtained by using a block factorization of the system in this paper we show how the nesting of black box solvers and preconditioners can be put together in a block triangular strategy to build a scalable block preconditioner for the stokes system discretized by divergence conforming b splines besides the known cavity flow problem we used for benchmark flows defined on complex geometries an eccentric annulus and hollow torus of an eccentric annular cross section

isogeometric analysis of the cahn hilliard phase field model abstract the cahn hilliard equation involves fourth order spatial derivatives finite element solutions are not common because primal variational formulations of fourth order operators are only well defined and integrable if the finite element basis functions are piecewise smooth and globally c 1 continuous there are a very limited number of two dimensional finite elements possessing c 1 continuity applicable to complex geometries but none in three dimensions we propose isogeometric analysis as a technology that possesses a unique combination of attributes for complex problems involving higher order differential operators namely higher order accuracy robustness two and three dimensional geometric flexibility compact support and most importantly the possibility of c 1 and higher order continuity a nurbs based variational formulation for the cahn hilliard equation was tested on two and three dimensional problems we present steady state solutions in two dimensions and for the first time in three dimensions to achieve these results an adaptive time stepping method is introduced we also present a technique for desensitizing calculations to dependence on mesh refinement this enables the calculation of topologically correct solutions on coarse meshes opening the way to practical engineering applications of phase field methodology

provably unconditionally stable second order time accurate mixed variational methods for phase field models we introduce provably unconditionally stable mixed variational methods for phase field models our formulation is based on a mixed finite element method for space discretization and a new second order accurate time integration algorithm the fully discrete formulation inherits the main characteristics of conserved phase dynamics namely mass conservation and nonlinear stability with respect to the free energy we illustrate the theory with the cahn hilliard equation but our method may be applied to other phase field models we also propose an adaptive time stepping version of the new time integration method we present some numerical examples that show the accuracy stability and robustness of the new method

pollution studies for high order isogeometric analysis and finite element for acoustic problems abstract it is well known that galerkin finite element methods suffer from pollution error when solving wave problems to reduce the pollution impact on the solution different approaches were proposed to enrich the finite element method with wave like functions so that the exact wavenumber is incorporated into the finite element approximation space solving wave problems with isogeometric analysis was also investigated in the literature where the superior behaviour of isogeometric analysis due to higher continuity in the underlying basis has been studied recently a plane wave enriched isogeometric analysis was introduced for acoustic problems however it remains unquantified the impact of these different approaches on the pollution or how they perform compared to each other in this work we show that isogeometric analysis outperforms finite element method in dealing with pollution we observe similar behaviour when both the methods are enriched with plane waves using higher order polynomials with fewer enrichment functions seems to improve the pollution compared to lower order polynomials with more functions however the latter still leads to smaller errors using similar number of degrees of freedom in conclusion we propose that partition of unity isogeometric analysis can be an efficient tool for wave problems as enrichment eliminates the need for domain re meshing at higher frequencies and also due to its ability to capture the exact geometry even on coarse meshes as well as its improved pollution behaviour

linear independence of the t spline blending functions associated with some particular t meshes based on the local refinement algorithm addressed in sederberg et al 2004 17 we analyze the linear independence of the bi cubic t spline blending functions corresponding to particular t meshes

overlapping schwarz methods for isogeometric analysis we construct and analyze an overlapping schwarz preconditioner for elliptic problems discretized with isogeometric analysis the preconditioner is based on partitioning the domain of the problem into overlapping subdomains solving local isogeometric problems on these subdomains and solving an additional coarse isogeometric problem associated with the subdomain mesh we develop an h analysis of the preconditioner showing in particular that the resulting algorithm is scalable and its convergence rate depends linearly on the ratio between subdomain and overlap sizes for fixed polynomial degree p and regularity k of the basis functions numerical results in two and three dimensional tests show the good convergence properties of the preconditioner with respect to the isogeometric discretization parameters h p k number of subdomains n overlap size and also jumps in the coefficients of the elliptic operator

a large deformation frictional contact formulation using nurbs based isogeometric analysis this paper focuses on the application of nurbs based isogeometric analysis to coulomb frictional contact problems between deformable bodies in the context of large deformations a mortar based approach is presented to treat the contact constraints whereby the discretization of the continuum is performed with arbitrary order nurbs as well as c0 continuous lagrange polynomial elements for comparison purposes the numerical examples show that the proposed contact formulation in conjunction with the nurbs discretization delivers accurate and robust predictions results of lower quality are obtained from the lagrange discretization as well as from a different contact formulation based on the enforcement of the contact constraints at every integration point on the contact surface copyright 2011 john wiley sons ltd

trigonometric generalized t splines abstract the paper s main aim consists of extending the t spline approach to trigonometric generalized b splines a particularly relevant case of non polynomial splines such goal can be achieved by a careful revision of some results concerning the basic properties of the trigonometric generalized b splines and by a formalization of the concept of t splines in the trigonometric setting moreover fundamental for the use of this new tool is the study of the noteworthy case with constant frequencies and of the linear independence of the corresponding blending functions which can be proved to be strongly linked to the linear independence of the polynomial blending functions associated to the same t mesh

weak coupling for isogeometric analysis of non matching and trimmed multi patch geometries nitsche s method can be used as a coupling tool for non matching discretizations by weakly enforcing interface constraints we explore the use of weak coupling based on nitsche s method in the context of higher order and higher continuity b splines and nurbs we demonstrate that weakly coupled spline discretizations do not compromise the accuracy of isogeometric analysis we show that the combination of weak coupling with the finite cell method opens the door for a truly isogeometric treatment of trimmed b spline and nurbs geometries that eliminates the need for costly reparameterization procedures we test our methodology for several relevant technical problems in two and three dimensions such as gluing together trimmed multi patches and connecting non matching meshes that contain b spline basis functions and standard triangular finite elements the results demonstrate that the concept of nitsche based weak coupling in conjunction with the finite cell method has the potential to considerably increase the flexibility of the design through analysis process in isogeometric analysis

weakening the tight coupling between geometry and simulation in isogeometric analysis from sub and super geometric analysis to geometry independent field approximation gift this paper presents an approach to generalize the concept of isogeometric analysis by allowing different spaces for the parameterization of the computational domain and for the approximation of the solution field the method inherits the main advantage of isogeometric analysis ie preserves the original exact computer aided design geometry for example given by nonuniform rational b splines but allows pairing it with an approximation space which is more suitable flexible for analysis for example t splines lr splines truncated hierarchical b splines and pht splines this generalization offers the advantage of adaptive local refinement without the need to reparameterize the domain and therefore without weakening the link with the computer aided design model we demonstrate the use of the method with different choices of geometry and field spaces and show that despite the failure of the standard patch test the optimum convergence rate is achieved for nonnested spaces

new rectangular plate elements based on twist kirchhoff theory abstract we introduce a new framework for the development of thin plate finite elements the twist kirchhoff theory a family of rectangular plate elements is derived that takes advantage of the special structure of this new theory particular attention is focused on the lowest order member of the family an eight degree of freedom four node element with mid side rotations whose stiffness matrix is exactly computed with one point gaussian quadrature we prove a convergence theorem for it and various error estimates these are also generalized to the higher order elements in the family numerical tests corroborate the theoretical results

isogeometric analysis and harmonic stator rotor coupling for simulating electric machines abstract this work proposes isogeometric analysis as an alternative to classical finite elements for simulating electric machines through the spline based isogeometric discretization it is possible to parametrize the circular arcs exactly thereby avoiding any geometrical error in the representation of the air gap where a high accuracy is mandatory to increase the generality of the method and to allow rotation the rotor and the stator computational domains are constructed independently as multipatch entities the two subdomains are then coupled using harmonic basis functions at the interface which gives rise to a saddle point problem the properties of isogeometric analysis combined with harmonic stator rotor coupling are presented the results and performance of the new approach are compared to the ones for a classical finite element method using a permanent magnet synchronous machine as an example

integrating cad with abaqus isogeometric analysis iga has been developed for more than a decade however the usage of iga is by far limited mostly within academic community the lack of automatic or semi automatic software platform of iga is one of the main bottlenecks that prevent iga from wide applications in industry in this paper we present a comprehensive iga software platform that allows iga to be incorporated into existing commercial software such as abaqus heading one step further to bridge the gap between design and analysis the proposed iga software framework takes advantage of user defined elements in abaqus linking with general iges files from commercial computer aided design packages rhino specific files and mesh data the platform includes all the necessary modules of the design through analysis pipeline pre processing surface and volumetric t spline construction analysis and post processing several practical application problems are studied to demonstrate the capability of the proposed software platform

isogeometric analysis and proper orthogonal decomposition for parabolic problems we investigate the combination of isogeometric analysis iga and proper orthogonal decomposition pod based on the galerkin method for model order reduction of linear parabolic partial differential equations for the proposed fully discrete scheme the associated numerical error features three components due to spatial discretization by iga time discretization with the theta scheme and eigenvalue truncation by pod first we prove a priori error estimates of the spatial iga semi discrete scheme then we show stability and prove a priori error estimates of the space time discrete scheme and the fully discrete iga theta pod galerkin scheme numerical tests are provided to show efficiency and accuracy of nurbs based iga for model order reduction in comparison with standard finite element based pod techniques

generalized b splines as a tool in isogeometric analysis the concept of isogeometric analysis has been proposed in 13 where nurbs are considered as basis of the analysis thanks to their ability to construct an exact geometric model in several practical applications and to their popularity in commercial cad systems in this paper we propose an alternative to the rational model presenting an isogeometric analysis approach based on generalized b splines geometric models exactly represented by generalized b splines include those generated by nurbs moreover generalized b splines possess all fundamental properties of algebraic b splines and nurbs including classical refinement processes as h p k refinements finally since generalized b splines are not confined to rational functions they behave completely similar to algebraic b splines with respect to differentiation and integration this seems to be of interest in the treatment of some relevant problems

overlapping schwarz methods for isogeometric analysis based on generalized b splines abstract generalized b splines gb splines are a special class of tchebycheff b splines that are smooth piecewise function with sections in more general spaces gb splines allow for an exact representation of conic sections as well as transcendental curves and thus they become very attractive for geometrical modeling and numerical simulation in this paper we present overlapping schwarz preconditioners for elliptic problems discretized with isogeometric analysis based on gb splines an h analysis of the proposed preconditioners shows an optimal convergence rate bound that is scalable in the number of subdomains and that is linear in the ratio between subdomain and overlap sizes numerical results in two and three dimensional tests confirm this analysis and also illustrate the good convergence properties of the preconditioner with respect to the discretization parameters the domain deformation and the jumping coefficients

a new procedure for the construction of hierarchical high order hdiv and hcurl finite element spaces this paper considers a systematic procedure for the construction of a hierarchy of high order finite element approximation for hdiv and hcurl spaces based on triangular and quadrilateral partitions of bidimensional domains the principle is to choose an appropriate set of vectors based on the geometry of each element which are multiplied by an available set of h 1 hierarchical scalar basic functions this strategy produces vector basis functions with continuous normal or tangent components on the elements interfaces properties that characterise functions in hdiv or hcurl respectively we also present a numerical study to evaluate the correct balancedness of the resulting hdiv spaces of degree k and l 2 spaces of degree k 1 on the resolution of the mixed formulation for a steklov eigenvalue problem

enhancing isogeometric analysis by a finite element based local refinement strategy abstract while isogeometric analysis has the potential to close the gap between computer aided design and finite element methods the underlying structure of nurbs non uniform rational b splines is a weakness when it comes to local refinement we propose a hybrid method that combines a globally c 1 continuous piecewise polynomial finite element basis with rational nurbs mappings in such a way that an isoparametric setting and exact geometry representation are preserved we define this basis over t meshes with a hierarchical structure that allows locally restricted refinement combined with a state of the art a posteriori error estimator we present an adaptive refinement procedure this concept is successfully demonstrated with the laplace equation advection diffusion problems and linear elasticity problems

accurate efficient and iso geometrically flexible collocation methods for phase field models we propose new collocation methods for phase field models our algorithms are based on isogeometric analysis a new technology that makes use of functions from computational geometry such as for example non uniform rational b splines nurbs nurbs exhibit excellent approximability and controllable global smoothness and can represent exactly most geometries encapsulated in computer aided design cad models these attributes permitted us to derive accurate efficient and geometrically flexible collocation methods for phase field models the performance of our method is demonstrated by several numerical examples of phase separation modeled by the cahn hilliard equation we feel that our method successfully combines the geometrical flexibility of finite elements with the accuracy and simplicity of pseudo spectral collocation methods and is a viable alternative to classical collocation methods

a hybrid igafem igabem formulation for two dimensional stationary magnetic and magneto mechanical field problems abstract isogeometric analysis iga can bridge the gap between geometrical and numerical modelling to this end the same functions used in computer aided design are applied to represent geometry and approximate field variables in the numerical model the concept has already been implemented to solve field problems using finite element fem and boundary element methods bem but coupling of both methods has not been applied in the current work an isogeometric fem bem coupling is proposed and applied to two dimensional stationary magnetic field problems while fem is used to model magnetisable bodies allowing for heterogeneous structures and non linear constitutive behaviour the bem domain accounts for the surrounding free space both methods are coupled on the surface of the magnetisable body due to this hybrid igafem bem approach no meshing of free space is necessary and truncation errors are avoided for problems to be solved on open infinite or semi infinite domains once the solution for the magnetic problem is obtained using the hybrid method igafem is used to solve a magneto mechanical field problem with one sided coupling in a subdomain of the magnetic problem this one sided coupling is realised by a magnetic stress tensor computed from the solution of the stationary magnetic field problem from the comparison of error norms and convergence rates for nurbs and discretisations based on lagrangian polynomials smaller errors and similar convergence rates are found for the proposed method for the same polynomial order of the basis functions and a comparable mesh size

isogeometric methods for computational electromagnetics b spline and t spline discretizations in this paper we introduce methods for electromagnetic wave propagation based on splines and on t splines we define spline spaces which form a de rham complex and following the isogeometric paradigm we map them on domains which are piecewise spline or nurbs geometries we analyze their geometric and topological structure as related to the connectivity of the underlying mesh and we present degrees of freedom together with their physical interpretation the theory is then extended to the case of meshes with t junctions leveraging on the recent theory of t splines the use of t splines enhance our spline methods with local refinement capability and numerical tests show the efficiency and the accuracy of the techniques we propose

the value of continuity refined isogeometric analysis and fast direct solvers abstract we propose the use of highly continuous finite element spaces interconnected with low continuity hyperplanes to maximize the performance of direct solvers starting from a highly continuous isogeometric analysis iga discretization we introduce c 0 separators to reduce the interconnection between degrees of freedom in the mesh by doing so both the solution time and best approximation errors are simultaneously improved we call the resulting method refined isogeometric analysis riga to illustrate the impact of the continuity reduction we analyze the number of floating point operations flops computational times and memory required to solve the linear system obtained by discretizing the laplace problem with structured meshes and uniform polynomial orders theoretical estimates demonstrate that an optimal continuity reduction may decrease the total computational time by a factor between p 2 and p 3 with p being the polynomial order of the discretization numerical results indicate that our proposed refined isogeometric analysis delivers a speed up factor proportional to p 2 in a 2 d mesh with four million elements and p 5 the linear system resulting from riga is solved 22 times faster than the one from highly continuous iga in a 3 d mesh with one million elements and p 3 the linear system is solved 15 times faster for the refined than the maximum continuity isogeometric analysis

extended isogeometric boundary element method xibem for two dimensional helmholtz problems isogeometric analysis is a topic of considerable interest in the field of numerical analysis the boundary element method bem requires only the bounding surface of geometries to be described this makes non uniform rational b splines nurbs which commonly describe the bounding curve or surface of geometries in cad software appear to be a natural tool for the approach this isogeometric analysis bem igabem provides accuracy benefits over conventional bem schemes due to the analytical geometry provided by nurbs when applied to wave problems it has been shown that enriching bem approximations with a partition of unity basis in what has become known as the pu bem allows highly accurate solutions to be obtained with a much reduced number of degrees of freedom this paper combines these approaches and presents an extended isogeometric bem xibem which uses partition of unity enriched nurbs functions this new approach provides benefits which surpass those of both the igabem and the pu bem two numerical examples are given a single scattering cylinder and a multiple scatterer made up of two capsules and a cylinder

divergence conforming discretization for stokes problem on locally refined meshes using lr b splines abstract to solve the incompressible flow problems using isogeometric analysis the div compatible spline spaces were originally introduced by buffa et al 2011 and later developed by evans 2011 in this paper we extend the div compatible spline spaces with local refinement capability using locally refined lr b splines over rectangular domains we argue that the spline spaces generated on locally refined meshes will satisfy compatibility provided they span the entire function spaces as governed by mourrain 2014 dimension formula we will in this work use the structured refined lr b splines as introduced by johannessen et al 2014 further we consider these div compatible lr b spline spaces to approximate the velocity and pressure fields in mixed discretization for stokes problem and a set of standard benchmark tests are performed to show the stability efficiency and the conservation properties of the discrete velocity fields in adaptive isogeometric analysis

isogeometric analysis in electromagnetics b splines approximation abstract we introduce a new discretization scheme for maxwell equations in two space dimension inspired by the new paradigm of isogeometric analysis introduced in hughes et al 2005 16 we propose an algorithm based on the use of bivariate b splines spaces suitably adapted to electromagnetics we construct b splines spaces of variable interelement regularity on the parametric domain these spaces and their push forwards on the physical domain form a de rham diagram and we use them to solve the maxwell source and eigen problem our scheme has the following features i is adapted to treat complex geometries ii is spectral correct iii provides regular e g globally c 0 discrete solutions of maxwell equations

weighted regularization of maxwell equations in polyhedral domains we present a new method of regularizing time harmonic maxwell equations by a bf grad div term adapted to the geometry of the domain this method applies to polygonal domains in two dimensions as well as to polyhedral domains in three dimensions in the presence of reentrant corners or edges the usual regularization is known to produce wrong solutions due the non density of smooth fields in the variational space we get rid of this undesirable effect by the introduction of special weights inside the divergence integral standard finite elements can then be used for the approximation of the solution this method proves to be numerically efficient

adjoint based airfoil optimization with adaptive isogeometric discontinuous galerkin method abstract in this work an adjoint based airfoil shape optimization algorithm is developed based on the adaptive isogeometric discontinuous galerkin method for compressible euler equations to investigate the significance of each design variable of airfoil b spline parameterization we first parameterize the airfoil by b spline curve approximation with some control points viewed as design variables and build the b spline representation of the flow field with the curve to apply the goal oriented h adaptive isogeometric dg method for flow solution then we compute and employ the discrete adjoint solutions for both multi target error estimation in adaptive mesh refinement with the isogeometric nature not only all the geometrical cells but also the numerical basis functions can be analytically expressed by the design variables indicating that the numerical solutions and objective could be differentiable with respect to those variables consequently the gradient is totally computed in an accurate approach and the sensitivity analysis is thus improved by reducing the spatial discretization error and introducing the analytical expression of derivative to reveal the key parameters for optimization in an intuitive and efficient manner although the sqp optimization algorithm is adopted in the paper the given accurate gradient can be applied to any gradient based optimization algorithm the proposed algorithm is demonstrated on rae2822 airfoil with inviscid transonic flow where the shape is optimized to minimize the drag coefficient at a constrained lift and airfoil area the numerical results show that the drag is much more sensitive to the design variables near the tailing edge at the beginning but sensitivity is reduced when optimal

contact treatment in isogeometric analysis with nurbs abstract we study nurbs based isogeometric analysis of contact problems and compare with standard c 0 continuous lagrange finite elements a knot to surface kts algorithm is developed to treat the contact constraints with nurbs contact surface discretizations qualitative studies deliver satisfactory results for various finite deformation frictionless thermoelastic contact problems quantitative studies based on the hertz problem suggest the need for a relaxation of the mechanical contact constraints that appear in the standard kts approach the improved mortar based kts algorithm delivers robust and accurate results for nurbs discretizations based on numerical examples we conclude that nurbs based isogeometric analysis is a viable technology for contact problems and offers potential accuracy as well as convergence improvements over c 0 continuous finite elements

multi patch discontinuous galerkin isogeometric analysis for wave propagation explicit time stepping and efficient mass matrix inversion abstract we present a class of spline finite element methods for time domain wave propagation which are particularly amenable to explicit time stepping the proposed methods utilize a discontinuous galerkin discretization to enforce continuity of the solution field across geometric patches in a multi patch setting which yields a mass matrix with convenient block diagonal structure over each patch we show how to accurately and efficiently invert mass matrices in the presence of curved geometries by using a weight adjusted approximation of the mass matrix inverse this approximation restores a tensor product structure while retaining provable high order accuracy and semi discrete energy stability we also estimate the maximum stable timestep for spline based finite elements and show that the use of spline spaces results in less stringent cfl restrictions than equivalent c 0 or discontinuous finite element spaces finally we explore the use of optimal knot vectors based on l 2 n widths we show how the use of optimal knot vectors can improve both approximation properties and the maximum stable timestep and present a simple heuristic method for approximating optimal knot positions numerical experiments confirm the accuracy and stability of the proposed methods

singular higher order complete vector bases for finite methods this paper presents new singular curl and divergence conforming vector bases that incorporate the edge conditions singular bases complete to arbitrarily high order are described in a unified and consistent manner for curved triangular and quadrilateral elements the higher order basis functions are obtained as the product of lowest order functions and silvester lagrange interpolatory polynomials with specially arranged arrays of interpolation points the completeness properties are discussed and these bases are proved to be fully compatible with the standard high order regular vector bases used in adjacent elements the curl divergence conforming singular bases guarantee tangential normal continuity along the edges of the elements allowing for the discontinuity of normal tangential components adequate modeling of the curl divergence and removal of spurious modes solutions these singular high order bases should provide more accurate and efficient numerical solutions of both surface integral and differential problems sample numerical results confirm the faster convergence of these bases on wedge problems

anisotropic nurbs approximation in isogeometric analysis abstract we investigate the approximation error of mapped nurbs spaces such as those used in isogeometric analysis in accordance to the isoparametric paradigm we extend the results of 6 to the case of anisotropic meshes such as those needed in the presence of boundary layers the theoretical results cover anisotropic h refinement and are supported by numerical tests involving different geometries

enriched isogeometric collocation for two dimensional time harmonic acoustics abstract in this work the isogeometric collocation iga c is paired with two types of enrichment plane wave pw and generalized harmonic polynomial ghp functions to solve 2d problems for the helmholtz equation a parametric study is conducted and a detailed assessment of the performance of the method in a number of benchmark problems is provided three different collocation methods are tested for the non enriched formulation namely greville abscissae ga approximated cauchy galerkin acg and superconvergent sc points showing that the acg scheme is the best choice in terms of overall error convergence rate and ease of solving the linear system then the influence of the number of shape functions in the original and enriched basis location and number of collocation points and wave number over both the convergence rate as well as the condition number of the stiffness matrix are studied the numerical results show that 1 there is an improvement over the non enriched formulation 2 the improvement depends on the choice of the number and type of enrichment and 3 the pollution error is not completely alleviated with the enriched formulations

isogeometric shape optimization of vibrating membranes we consider a model problem of isogeometric shape optimization of vibrating membranes whose shapes are allowed to vary freely the main obstacle we face is the need for robust and inexpensive extension of a b spline parametrization from the boundary of a domain onto its interior a task which has to be performed in every optimization iteration we experiment with two numerical methods one is based on the idea of constructing a quasi conformal mapping whereas the other is based on a spring based mesh model for carrying out this task which turn out to work sufficiently well in the present situation we perform a number of numerical experiments with our isogeometric shape optimization algorithm and present smooth optimized membrane shapes our conclusion is that isogeometric analysis fits well with shape optimization

proper orthogonal decomposition with supg stabilized isogeometric analysis for reduced order modelling of unsteady convection dominated convection diffusion reaction problems abstract we consider reduced order modelling of unsteady convection dominated convection diffusion reaction problems with proper orthogonal decomposition pod in combination with isogeometric analysis isogeometric analysis has potential advantages in exact geometry representations efficient mesh generation different h p and k refinements and smooth bspline nurbs basis functions in order to compensate the oscillations caused by the convection dominated effect the streamline upwind petrov galerkin supg stabilization method is used both in generation of snapshots and pod galerkin method based on the recent novel and promising discretization method isogeometric analysis we propose a new fully discrete supg stabilized scheme the associated numerical error features three components due to spatial discretization by isogeometric analysis with supg stabilization time discretization with the crank nicolson scheme and modes truncation by pod we show a priori error estimates of the fully discrete scheme and give suitable stabilization parameters numerically a variety of two and three dimensional benchmark tests and numerical examples are provided to show the effectiveness accuracy and efficiency of the reduced order modelling methods by virtue of potential advantages of isogeometric analysis

a simple algorithm for obtaining nearly optimal quadrature rules for nurbs based isogeometric analysis we develop new quadrature rules for isogeometric analysis based on the solution of a local nonlinear problem a simple and robust algorithm is developed to determine the rules which are exact for important b spline spaces of uniform and geometrically stretched knot spacings we consider both periodic and open knot vector configurations and illustrate the efficiency of the rules on selected boundary value problems we find that the rules are almost optimally efficient but much easier to obtain than optimal rules which require the solution of global nonlinear problems that are often ill posed

combined instrumental variable and subspace fitting approach to parameter estimation of noisy input output systems the paper considers the problem of estimating the parameters of linear discrete time systems from noise corrupted input output measurements under fairly general conditions the output and input noises may be auto correlated and they may be cross correlated as well by using the instrumental variable iv principle a covariance matrix is obtained the singular vectors of which bear complete information on the parameters of the system under study a weighted subspace fitting wsf procedure is then employed on the sample singular vectors to derive estimates of the parameters of the system the combined iv wsf method proposed in the present paper is noniterative and simple to use its large sample statistical performance is analyzed in detail and the theoretical results so obtained are used to predict the behavior of the method in samples with practical lengths several numerical examples are included to show the agreement between the theoretically predicted and the empirically observed performances

isogeometric divergence conforming b splines for the unsteady navier stokes equations divergence conforming b splines are developed for application to the incompressible navier stokes equations on geometrically mapped domains these enable smooth pointwise divergence free solutions and thus satisfy mass conservation in the strongest possible sense semi discrete methods based on divergence conforming b splines are shown to conserve linear and angular momentum and satisfy balance laws for energy vorticity enstrophy and helicity these are geometric structure preserving quantities and numerical simulations that are sensitive to them are shown to be qualitatively correct and quantitatively accurate the methods developed are anticipated to open new doors to the practical calculation of complex flows and to studies of their physical behavior

numerical solution to the time dependent maxwell equations in two dimensional singular domains in this paper we present a method to solve numerically the time dependent maxwell equations in nonsmooth and nonconvex domains indeed the solution is not of regularity h1 in space in general moreover the space of h1 regular fields is not dense in the space of solutions thus an h1 conforming finite element method can fail even with mesh refinement the situation is different than in the case of the laplace problem or of the lame system for which mesh refinement or the addition of conforming singular functions work to cope with this difficulty the singular complement method is introduced this method consists of adding some well chosen test functions these functions are derived from the singular solutions of the laplace problem also the scm preserves the interesting features of the original method easiness of implementation low memory requirements small cost in terms of the cpu time to ascertain its validity some concrete problems are solved numerically