Curtin University Homepage
  • Library
  • Help
    • Admin

    espace - Curtin’s institutional repository

    JavaScript is disabled for your browser. Some features of this site may not work without it.
    View Item 
    • espace Home
    • espace
    • Curtin Research Publications
    • View Item
    • espace Home
    • espace
    • Curtin Research Publications
    • View Item

    Refined Isogeometric Analysis for a preconditioned conjugate gradient solver

    265700.pdf (349.0Kb)
    Access Status
    Open access
    Authors
    Garcia, D.
    Pardo, D.
    Dalcin, L.
    Calo, Victor
    Date
    2018
    Type
    Journal Article
    
    Metadata
    Show full item record
    Citation
    Garcia, 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.
    Source Title
    Computer Methods in Applied Mechanics and Engineering
    DOI
    10.1016/j.cma.2018.02.006
    ISSN
    0045-7825
    School
    School of Earth and Planetary Sciences (EPS)
    URI
    http://hdl.handle.net/20.500.11937/67251
    Collection
    • Curtin Research Publications
    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.

    Related items

    Showing items related by title, author, creator and subject.

    • A sequential computational approach to optimal control problems for differential-algebraic systems based on efficient implicit Runge–Kutta integration
      Jiang, C.; Xie, K.; Yu, C.; Yu, M.; Wang, H.; He, Y.; Teo, Kok Lay (2018)
      Efficient and reliable integrators are indispensable for the design of sequential solvers for optimal control problems involving continuous dynamics, especially for real-time applications. In this paper, optimal control ...
    • A review of block Krylov subspace methods for multisource electromagnetic modelling
      Puzyrev, Volodymyr; Cela, J. (2015)
      Practical applications of controlled-source electromagnetic (EM) modelling require solutions for multiple sources at several frequencies, thus leading to a dramatic increase of the computational cost. In this paper, we ...
    • Exploring the use of SPIKE-based solvers on large electromagnetic modeling
      Rodriguez Bernabeu, S.; Puzyrev, Vladimir; Hanzich, M.; Fernández, S. (2016)
      Frequency-domain seismic and electromagnetic modeling requires solving the linear systems resulting from the discretization of the corresponding time-harmonic equations. Geophysical inversion is typically performed using ...
    Advanced search

    Browse

    Communities & CollectionsIssue DateAuthorTitleSubjectDocument TypeThis CollectionIssue DateAuthorTitleSubjectDocument Type

    My Account

    Admin

    Statistics

    Most Popular ItemsStatistics by CountryMost Popular Authors

    Follow Curtin

    • 
    • 
    • 
    • 
    • 

    CRICOS Provider Code: 00301JABN: 99 143 842 569TEQSA: PRV12158

    Copyright | Disclaimer | Privacy statement | Accessibility

    Curtin would like to pay respect to the Aboriginal and Torres Strait Islander members of our community by acknowledging the traditional owners of the land on which the Perth campus is located, the Whadjuk people of the Nyungar Nation; and on our Kalgoorlie campus, the Wongutha people of the North-Eastern Goldfields.