PetIGA: A framework for high-performance isogeometric analysis

Lisandro Dalcin; N. Collier; Philippe Vignal; Adriano Mauricio Cortes; Victor M. Calo

doi:10.1016/j.cma.2016.05.011

PetIGA: A framework for high-performance isogeometric analysis

Lisandro Dalcin, N. Collier, Philippe Vignal, Adriano Mauricio Cortes, Victor M. Calo

Research output: Contribution to journal › Article › peer-review

123 Scopus citations

Abstract

We present PetIGA, a code framework to approximate the solution of partial differential equations using isogeometric analysis. PetIGA can be used to assemble matrices and vectors which come from a Galerkin weak form, discretized with Non-Uniform Rational B-spline basis functions. We base our framework on PETSc, a high-performance library for the scalable solution of partial differential equations, which simplifies the development of large-scale scientific codes, provides a rich environment for prototyping, and separates parallelism from algorithm choice. We describe the implementation of PetIGA, and exemplify its use by solving a model nonlinear problem. To illustrate the robustness and flexibility of PetIGA, we solve some challenging nonlinear partial differential equations that include problems in both solid and fluid mechanics. We show strong scaling results on up to 40964096 cores, which confirm the suitability of PetIGA for large scale simulations.

Original language	English (US)
Pages (from-to)	151-181
Number of pages	31
Journal	Computer Methods in Applied Mechanics and Engineering
Volume	308
DOIs	https://doi.org/10.1016/j.cma.2016.05.011
State	Published - May 26 2016

Bibliographical note

KAUST Repository Item: Exported on 2020-10-01
Acknowledgements: We would like to acknowledge the open source software packages that made this work possible: PETSc [24], NumPy [76], matplotlib [77], IPython [78]. We would like to thank Lina María Bernal Martinez, Gabriel Andres Espinosa Barrios, Federico Fuentes-Caycedo, Juan Camilo Mahecha Zambrano for their work on the hyper-elasticity implementation as a final project to the Non-linear Finite Element class taught by V.M. Calo and N. Collier for the Mechanical Engineering Department at Universidad de Los Andes in Bogotá, Colombia in July 2012. We would like to thank Adel Sarmiento Rodriguez for the visualization work on figures 12 and 13.

This work is part of the European Union’s Horizon 2020 research and innovation programme of the Marie Skłodowska-Curie grant agreement No. 644602. This work was also supported by the Center for Numerical Porous Media at King Abdullah University of Science and Technology and Agencia Nacional de Promoción Científica y Tecnológica grants PICT 0938–13, PICT 2660–14 and PICT-E 0191–14.

ASJC Scopus subject areas

General Physics and Astronomy
Mechanics of Materials
Mechanical Engineering
Computational Mechanics
Computer Science Applications

Access to Document

10.1016/j.cma.2016.05.011

https://repository.kaust.edu.sa/bitstream/10754/611201/1/1-s2.0-S0045782516303401-main.pdf

Cite this

@article{f4b3f87811a748028ea142945f785e80,

title = "PetIGA: A framework for high-performance isogeometric analysis",

abstract = "We present PetIGA, a code framework to approximate the solution of partial differential equations using isogeometric analysis. PetIGA can be used to assemble matrices and vectors which come from a Galerkin weak form, discretized with Non-Uniform Rational B-spline basis functions. We base our framework on PETSc, a high-performance library for the scalable solution of partial differential equations, which simplifies the development of large-scale scientific codes, provides a rich environment for prototyping, and separates parallelism from algorithm choice. We describe the implementation of PetIGA, and exemplify its use by solving a model nonlinear problem. To illustrate the robustness and flexibility of PetIGA, we solve some challenging nonlinear partial differential equations that include problems in both solid and fluid mechanics. We show strong scaling results on up to 40964096 cores, which confirm the suitability of PetIGA for large scale simulations.",

author = "Lisandro Dalcin and N. Collier and Philippe Vignal and Cortes, {Adriano Mauricio} and Calo, {Victor M.}",

note = "KAUST Repository Item: Exported on 2020-10-01 Acknowledgements: We would like to acknowledge the open source software packages that made this work possible: PETSc [24], NumPy [76], matplotlib [77], IPython [78]. We would like to thank Lina Mar{\'i}a Bernal Martinez, Gabriel Andres Espinosa Barrios, Federico Fuentes-Caycedo, Juan Camilo Mahecha Zambrano for their work on the hyper-elasticity implementation as a final project to the Non-linear Finite Element class taught by V.M. Calo and N. Collier for the Mechanical Engineering Department at Universidad de Los Andes in Bogot{\'a}, Colombia in July 2012. We would like to thank Adel Sarmiento Rodriguez for the visualization work on figures 12 and 13. This work is part of the European Union{\textquoteright}s Horizon 2020 research and innovation programme of the Marie Sk{\l}odowska-Curie grant agreement No. 644602. This work was also supported by the Center for Numerical Porous Media at King Abdullah University of Science and Technology and Agencia Nacional de Promoci{\'o}n Cient{\'i}fica y Tecnol{\'o}gica grants PICT 0938–13, PICT 2660–14 and PICT-E 0191–14.",

year = "2016",

month = may,

day = "26",

doi = "10.1016/j.cma.2016.05.011",

language = "English (US)",

volume = "308",

pages = "151--181",

journal = "Computer Methods in Applied Mechanics and Engineering",

issn = "0045-7825",

publisher = "Elsevier",

}

TY - JOUR

T1 - PetIGA: A framework for high-performance isogeometric analysis

AU - Dalcin, Lisandro

AU - Collier, N.

AU - Vignal, Philippe

AU - Cortes, Adriano Mauricio

AU - Calo, Victor M.

N1 - KAUST Repository Item: Exported on 2020-10-01 Acknowledgements: We would like to acknowledge the open source software packages that made this work possible: PETSc [24], NumPy [76], matplotlib [77], IPython [78]. We would like to thank Lina María Bernal Martinez, Gabriel Andres Espinosa Barrios, Federico Fuentes-Caycedo, Juan Camilo Mahecha Zambrano for their work on the hyper-elasticity implementation as a final project to the Non-linear Finite Element class taught by V.M. Calo and N. Collier for the Mechanical Engineering Department at Universidad de Los Andes in Bogotá, Colombia in July 2012. We would like to thank Adel Sarmiento Rodriguez for the visualization work on figures 12 and 13. This work is part of the European Union’s Horizon 2020 research and innovation programme of the Marie Skłodowska-Curie grant agreement No. 644602. This work was also supported by the Center for Numerical Porous Media at King Abdullah University of Science and Technology and Agencia Nacional de Promoción Científica y Tecnológica grants PICT 0938–13, PICT 2660–14 and PICT-E 0191–14.

PY - 2016/5/26

Y1 - 2016/5/26

N2 - We present PetIGA, a code framework to approximate the solution of partial differential equations using isogeometric analysis. PetIGA can be used to assemble matrices and vectors which come from a Galerkin weak form, discretized with Non-Uniform Rational B-spline basis functions. We base our framework on PETSc, a high-performance library for the scalable solution of partial differential equations, which simplifies the development of large-scale scientific codes, provides a rich environment for prototyping, and separates parallelism from algorithm choice. We describe the implementation of PetIGA, and exemplify its use by solving a model nonlinear problem. To illustrate the robustness and flexibility of PetIGA, we solve some challenging nonlinear partial differential equations that include problems in both solid and fluid mechanics. We show strong scaling results on up to 40964096 cores, which confirm the suitability of PetIGA for large scale simulations.

AB - We present PetIGA, a code framework to approximate the solution of partial differential equations using isogeometric analysis. PetIGA can be used to assemble matrices and vectors which come from a Galerkin weak form, discretized with Non-Uniform Rational B-spline basis functions. We base our framework on PETSc, a high-performance library for the scalable solution of partial differential equations, which simplifies the development of large-scale scientific codes, provides a rich environment for prototyping, and separates parallelism from algorithm choice. We describe the implementation of PetIGA, and exemplify its use by solving a model nonlinear problem. To illustrate the robustness and flexibility of PetIGA, we solve some challenging nonlinear partial differential equations that include problems in both solid and fluid mechanics. We show strong scaling results on up to 40964096 cores, which confirm the suitability of PetIGA for large scale simulations.

UR - http://hdl.handle.net/10754/611201

UR - http://linkinghub.elsevier.com/retrieve/pii/S0045782516303401

UR - http://www.scopus.com/inward/record.url?scp=84973352693&partnerID=8YFLogxK

U2 - 10.1016/j.cma.2016.05.011

DO - 10.1016/j.cma.2016.05.011

M3 - Article

SN - 0045-7825

VL - 308

SP - 151

EP - 181

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

ER -

PetIGA: A framework for high-performance isogeometric analysis

Abstract

Bibliographical note

ASJC Scopus subject areas

Access to Document

Other files and links

Fingerprint

Cite this