Show simple item record

dc.contributor.authorGarcia, D.
dc.contributor.authorPardo, D.
dc.contributor.authorDalcin, L.
dc.contributor.authorCalo, Victor
dc.date.accessioned2018-05-18T07:57:49Z
dc.date.available2018-05-18T07:57:49Z
dc.date.created2018-05-18T00:23:16Z
dc.date.issued2018
dc.identifier.citationGarcia, D. and Pardo, D. and Dalcin, L. and Calo, V. 2018. Refined Isogeometric Analysis for a preconditioned conjugate gradient solver. Computer Methods in Applied Mechanics and Engineering. 335: pp. 490-509.
dc.identifier.urihttp://hdl.handle.net/20.500.11937/67251
dc.identifier.doi10.1016/j.cma.2018.02.006
dc.description.abstract

Starting from a highly continuous Isogeometric Analysis (IGA) discretization, refined Isogeometric Analysis (rIGA) introduces C 0 hyperplanes that act as separators for the direct LU factorization solver. As a result, the total computational cost required to solve the corresponding system of equations using a direct LU factorization solver dramatically reduces (up to a factor of 55) (Garcia et al., 2017). At the same time, rIGA enriches the IGA spaces, thus improving the best approximation error. In this work, we extend the complexity analysis of rIGA to the case of iterative solvers. We build an iterative solver as follows: we first construct the Schur complements using a direct solver over small subdomains (macro-elements). We then assemble those Schur complements into a global skeleton system. Subsequently, we solve this system iteratively using Conjugate Gradients (CG) with an incomplete LU (ILU) preconditioner. For a 2D Poisson model problem with a structured mesh and a uniform polynomial degree of approximation, rIGA achieves moderate savings with respect to IGA in terms of the number of Floating Point Operations (FLOPs) and computational time (in seconds) required to solve the resulting system of linear equations. For instance, for a mesh with four million elements and polynomial degree p=3, the iterative solver is approximately 2.6 times faster (in time) when applied to the rIGA system than to the IGA one. These savings occur because the skeleton rIGA system contains fewer non-zero entries than the IGA one. The opposite situation occurs for 3D problems, and as a result, 3D rIGA discretizations provide no gains with respect to their IGA counterparts when considering iterative solvers.

dc.publisherElsevier BV
dc.titleRefined Isogeometric Analysis for a preconditioned conjugate gradient solver
dc.typeJournal Article
dcterms.source.volume335
dcterms.source.startPage490
dcterms.source.endPage509
dcterms.source.issn0045-7825
dcterms.source.titleComputer Methods in Applied Mechanics and Engineering
curtin.departmentSchool of Earth and Planetary Sciences (EPS)
curtin.accessStatusOpen access


Files in this item

Thumbnail

This item appears in the following Collection(s)

Show simple item record