Fast Isogeometric Boundary Element Method based on Independent Field Approximation

An isogeometric boundary element method for problems in elasticity is presented, which is based on an independent approximation for the geometry, traction and displacement field. This enables a flexible choice of refinement strategies, permits an efficient evaluation of geometry related information, a mixed collocation scheme which deals with discontinuous tractions along non-smooth boundaries and a significant reduction of the right hand side of the system of equations for common boundary conditions. All these benefits are achieved without any loss of accuracy compared to conventional isogeometric formulations. The system matrices are approximated by means of hierarchical matrices to reduce the computational complexity for large scale analysis. For the required geometrical bisection of the domain, a strategy for the evaluation of bounding boxes containing the supports of NURBS basis functions is presented. The versatility and accuracy of the proposed methodology is demonstrated by convergence studies showing optimal rates and real world examples in two and three dimensions.Comment: 32 pages, 27 figure

The INTERNODES method for the treatment of non-conforming multipatch geometries in Isogeometric Analysis

In this paper we apply the INTERNODES method to solve second order elliptic problems discretized by Isogeometric Analysis methods on non-conforming multiple patches in 2D and 3D geometries. INTERNODES is an interpolation-based method that, on each interface of the configuration, exploits two independent interpolation operators to enforce the continuity of the traces and of the normal derivatives. INTERNODES supports non-conformity on NURBS spaces as well as on geometries. We specify how to set up the interpolation matrices on non-conforming interfaces, how to enforce the continuity of the normal derivatives and we give special attention to implementation aspects. The numerical results show that INTERNODES exhibits optimal convergence rate with respect to the mesh size of the NURBS spaces an that it is robust with respect to jumping coefficients.Comment: Accepted for publication in Computer Methods in Applied Mechanics and Engineerin

Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion

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 result in less stringent CFL restrictions than equivalent piecewise continuous or discontinuous finite element spaces. Finally, we explore the use of optimal knot vectors based on L2 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

Certified Reduced Basis Method for Affinely Parametric Isogeometric Analysis NURBS Approximation

In this work we apply reduced basis methods for parametric PDEs to an isogeometric formulation based on NURBS. We propose an integrated and complete work pipeline from CAD to parametrization of domain geometry, then from full order to certified reduced basis solution. IsoGeometric Analysis (IGA), as well as reduced basis methods for parametric PDEs growing research themes in scientific computing and computational mechanics. Their combination enhances the solution of some class of problems, especially the ones characterized by parametrized geometries. This work shows that it is also possible for some class of problems to deal with affine geometrical parametrization combined with a NURBS IGA formulation. In this work we show a certification of accuracy and a complete integration between IGA formulation and parametric certified greedy RB formulation by introducing two numerical examples in heat transfer with different parametrization