Matrix-free weighted quadrature for a computationally efficient isogeometric k-method
From MaRDI portal
Publication:1986196
DOI10.1016/J.CMA.2018.04.029zbMATH Open1440.65230arXiv1712.08565OpenAlexW2796520447MaRDI QIDQ1986196FDOQ1986196
Authors: Giancarlo Sangalli, Mattia Tani 
Publication date: 8 April 2020
Published in: Computer Methods in Applied Mechanics and Engineering (Search for Journal in Brave)
Abstract: The -method is the isogeometric method based on splines (or NURBS, etc.) with maximum regularity. When implemented following the paradigms of classical finite element methods, the computational resources required by the method are prohibitive even for moderate degree. In order to address this issue, we propose a matrix-free strategy combined with weighted quadrature, which is an ad-hoc strategy to compute the integrals of the Galerkin system. Matrix-free weighted quadrature (MF-WQ) speeds up matrix operations, and, perhaps even more important, greatly reduces memory consumption. Our strategy also requires an efficient preconditioner for the linear system iterative solver. In this work we deal with an elliptic model problem, and adopt a preconditioner based on the Fast Diagonalization method, an old idea to solve Sylvester-like equations. Our numerical tests show that the isogeometric solver based on MF-WQ is faster than standard approaches (where the main cost is the matrix formation by standard Gaussian quadrature) even for low degree. But the main achievement is that, with MF-WQ, the -method gets orders of magnitude faster by increasing the degree, given a target accuracy. Therefore, we are able to show the superiority, in terms of computational efficiency, of the high-degree -method with respect to low-degree isogeometric discretizations. What we present here is applicable to more complex and realistic differential problems, but its effectiveness will depend on the preconditioner stage, which is as always problem-dependent. This situation is typical of modern high-order methods: the overall performance is mainly related to the quality of the preconditioner.
Full work available at URL: https://arxiv.org/abs/1712.08565
Recommendations
- Fast formation of isogeometric Galerkin matrices by weighted quadrature
- Fast formation and assembly of finite element matrices with application to isogeometric linear elasticity
- Quadrature rules in the isogeometric Galerkin method: state of the art and an introduction to weighted quadrature
- Fast formation of isogeometric Galerkin matrices via integration by interpolation and look-up
- Efficient matrix assembly in isogeometric analysis with hierarchical B-splines
Numerical computation using splines (65D07) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)
Cites Work
- The ubiquitous Kronecker product
- Duality and unified analysis of discrete approximations in structural dynamics and wave propagation: comparison of p-method finite elements with \(k\)-method NURBS
- Isogeometric analysis. Toward integration of CAD and FEA
- Some estimates for \(h\)-\(p\)-\(k\)-refinement in isogeometric analysis
- Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement
- Efficient quadrature for NURBS-based isogeometric analysis
- Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows
- Studies of refinement and continuity in isogeometric structural analysis
- Patient-specific isogeometric structural analysis of aortic valve closure
- \(n\)-widths, sup-infs, and optimality ratios for the \(k\)-version of the isogeometric finite element method
- A globally convergent matrix-free algorithm for implicit time-marching schemes arising in finite element analysis in fluids
- Robust and optimal multi-iterative techniques for Iga Galerkin linear systems
- Reduced Bézier element quadrature rules for quadratic and cubic splines in isogeometric analysis
- Direct solution of partial difference equations by tensor product methods
- The cost of continuity: a study of the performance of isogeometric finite elements using direct solvers
- Multigrid methods for isogeometric discretization
- An efficient implicit discontinuous spectral Galerkin method
- BPX-preconditioning for isogeometric analysis
- Fast simplicial finite element algorithms using Bernstein polynomials
- Mathematical analysis of variational isogeometric methods
- A new design for the implementation of isogeometric analysis in Octave and Matlab: GeoPDEs 3.0
- Isogeometric preconditioners based on fast solvers for the Sylvester equation
- Robust Multigrid for Isogeometric Analysis Based on Stable Splittings of Spline Spaces
- Optimal quadrature for univariate and tensor product splines
- Optimal quadrature rules for odd-degree spline spaces and their application to tensor-product-based isogeometric analysis
- Low rank tensor methods in Galerkin-based isogeometric analysis
- A preconditioning strategy for linear systems arising from nonsymmetric schemes in isogeometric analysis
- Fast formation of isogeometric Galerkin matrices by weighted quadrature
- Optimal and reduced quadrature rules for tensor product and hierarchically refined splines in isogeometric analysis
- Efficient matrix computation for tensor-product isogeometric analysis: the use of sum factorization
- Robust isogeometric preconditioners for the Stokes system based on the fast diagonalization method
- Bernstein-Bézier finite elements on tetrahedral-hexahedral-pyramidal partitions
- BPX preconditioners for isogeometric analysis using analysis-suitable T-splines
Cited In (42)
- A simplified Kirchhoff-Love large deformation model for elastic shells and its effective isogeometric formulation
- Space–time least–squares isogeometric method and efficient solver for parabolic problems
- Robust isogeometric preconditioners for the Stokes system based on the fast diagonalization method
- Isogeometric analysis on V-reps: first results
- Sum factorization techniques in isogeometric analysis
- A robust penalty coupling of non-matching isogeometric Kirchhoff-Love shell patches in large deformations
- An efficient solver for space-time isogeometric Galerkin methods for parabolic problems
- Fast immersed boundary method based on weighted quadrature
- A matrix-free isogeometric Galerkin method for Karhunen-Loève approximation of random fields using tensor product splines, tensor contraction and interpolation based quadrature
- A domain decomposition method for isogeometric multi-patch problems with inexact local solvers
- Explicit isogeometric collocation for the dynamics of three-dimensional beams undergoing finite motions
- Compressive isogeometric analysis
- The surrogate matrix methodology: accelerating isogeometric analysis of waves
- Isogeometric collocation for implicit dynamics of three-dimensional beams undergoing finite motions
- Phase-field modeling for polarization evolution in ferroelectric materials via an isogeometric collocation method
- The surrogate matrix methodology: low-cost assembly for isogeometric analysis
- Fast formation of isogeometric Galerkin matrices via integration by interpolation and look-up
- Overlapping additive Schwarz preconditioners for isogeometric collocation discretizations of linear elasticity
- Efficient extraction of hierarchical B-splines for local refinement and coarsening of isogeometric analysis
- A sparse-grid isogeometric solver
- A fully explicit isogeometric collocation formulation for the dynamics of geometrically exact beams
- Easy and efficient preconditioning of the isogeometric mass matrix
- A mixed integration point (MIP) formulation for hyperelastic Kirchhoff-Love shells for nonlinear static and dynamic analysis
- Robust approximation error estimates and multigrid solvers for isogeometric multi-patch discretizations
- Some remarks on spectral convergence and stability of iso-geometric analysis
- Numerical quadrature for Gregory quads
- Weighted isogeometric collocation based on spline projectors
- Effects of parameterization and knot placement techniques on primal and mixed isogeometric collocation formulations of spatial shear-deformable beams with varying curvature and torsion
- A minimal stabilization procedure for isogeometric methods on trimmed geometries
- On the fast assemblage of finite element matrices with application to nonlinear heat transfer problems
- A computational comparison between isogeometric analysis and spectral element methods: accuracy and spectral properties
- High‐fidelity tensor‐decomposition based matrix formation for isogeometric buckling analysis of laminated shells with solid‐shell formulation
- Mixed stress-displacement isogeometric collocation for nearly incompressible elasticity and elastoplasticity
- Patch-wise quadrature of trimmed surfaces in isogeometric analysis
- Efficient matrix assembly in isogeometric analysis with hierarchical B-splines
- Efficient matrix computation for isogeometric discretizations with hierarchical B-splines in any dimension
- Coupling a NURBS contact interface with a higher order finite element discretization for contact problems using the mortar method
- Explicit high-order generalized-\(\alpha\) methods for isogeometric analysis of structural dynamics
- An improved isogeometric collocation formulation for spatial multi-patch shear-deformable beams with arbitrary initial curvature
- Fast and multiscale formation of isogeometric matrices of microstructured geometric models
- Weighted quadrature for hierarchical B-splines
- A cheap preconditioner based on fast diagonalization method for matrix-free weighted-quadrature isogeometric analysis applied to nonlinear transient heat transfer problems
Uses Software
This page was built for publication: Matrix-free weighted quadrature for a computationally efficient isogeometric \(k\)-method
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1986196)
