The value of continuity: Refined isogeometric analysis and fast direct solvers
Abstract
We propose the use of highly continuous finite element spaces interconnected with low continuity hyperplanes to maximize the performance of direct solvers. Starting from a highly continuous Isogeometric Analysis (IGA) discretization, we introduce -separators to reduce the interconnection between degrees of freedom in the mesh. By doing so, both the solution time and best approximation errors are simultaneously improved. We call the resulting method “refined Isogeometric Analysis (rIGA)”. To illustrate the impact of the continuity reduction, we analyze the number of Floating Point Operations (FLOPs), computational times, and memory required to solve the linear system obtained by discretizing the Laplace problem with structured meshes and uniform polynomial orders. Theoretical estimates demonstrate that an optimal continuity reduction may decrease the total computational time by a factor between and , with being the polynomial order of the discretization. Numerical results indicate that our proposed refined isogeometric analysis delivers a speed-up factor proportional to . In a mesh with four million elements and , the linear system resulting from rIGA is solved 22 times faster than the one from highly continuous IGA. In a mesh with one million elements and , the linear system is solved 15 times faster for the refined thanmore »
- Authors:
-
- Basque Center for Applied Mathematics (BCAM), Bilbao (Spain)
- Basque Center for Applied Mathematics (BCAM), Bilbao (Spain); Univ. of the Basque Country UPV.EHU, Leioa (Spain); Ikerbasque (Basque Foundation for Sciences), Bilbao (Spain)
- King Abdullah Univ. of Science and Technology (KAUST), Thuwal (Saudi Arabia); Consejo Nacional de Investigaciones Cientificas y Tecnicas, Santa Fe (Argentina); Univ. Nacional del Litoral, Santa Fe (Argentina)
- AGH Univ. of Science and Technology, Krakow (Poland)
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
- Curtin Univ., Bentley, WA (Australia); Commonwealth Scientific and Industrial Research Organisation (CSIRO), Kensington, WA (Australia)
- Publication Date:
- Research Org.:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 1327763
- Alternate Identifier(s):
- OSTI ID: 1429198
- Grant/Contract Number:
- AC05-00OR22725
- Resource Type:
- Journal Article: Accepted Manuscript
- Journal Name:
- Computer Methods in Applied Mechanics and Engineering
- Additional Journal Information:
- Journal Name: Computer Methods in Applied Mechanics and Engineering; Journal ID: ISSN 0045-7825
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; isogeometric analysis (IGA); finite element analysis (FEA); refined isogeometric analysis (rIGA); direct solvers; multi-frontal solvers; k-refinement
Citation Formats
Garcia, Daniel, Pardo, David, Dalcin, Lisandro, Paszynski, Maciej, Collier, Nathan, and Calo, Victor M. The value of continuity: Refined isogeometric analysis and fast direct solvers. United States: N. p., 2016.
Web. doi:10.1016/j.cma.2016.08.017.
Garcia, Daniel, Pardo, David, Dalcin, Lisandro, Paszynski, Maciej, Collier, Nathan, & Calo, Victor M. The value of continuity: Refined isogeometric analysis and fast direct solvers. United States. https://doi.org/10.1016/j.cma.2016.08.017
Garcia, Daniel, Pardo, David, Dalcin, Lisandro, Paszynski, Maciej, Collier, Nathan, and Calo, Victor M. 2016.
"The value of continuity: Refined isogeometric analysis and fast direct solvers". United States. https://doi.org/10.1016/j.cma.2016.08.017. https://www.osti.gov/servlets/purl/1327763.
@article{osti_1327763,
title = {The value of continuity: Refined isogeometric analysis and fast direct solvers},
author = {Garcia, Daniel and Pardo, David and Dalcin, Lisandro and Paszynski, Maciej and Collier, Nathan and Calo, Victor M.},
abstractNote = {We propose the use of highly continuous finite element spaces interconnected with low continuity hyperplanes to maximize the performance of direct solvers. Starting from a highly continuous Isogeometric Analysis (IGA) discretization, we introduce -separators to reduce the interconnection between degrees of freedom in the mesh. By doing so, both the solution time and best approximation errors are simultaneously improved. We call the resulting method “refined Isogeometric Analysis (rIGA)”. To illustrate the impact of the continuity reduction, we analyze the number of Floating Point Operations (FLOPs), computational times, and memory required to solve the linear system obtained by discretizing the Laplace problem with structured meshes and uniform polynomial orders. Theoretical estimates demonstrate that an optimal continuity reduction may decrease the total computational time by a factor between and , with being the polynomial order of the discretization. Numerical results indicate that our proposed refined isogeometric analysis delivers a speed-up factor proportional to . In a mesh with four million elements and , the linear system resulting from rIGA is solved 22 times faster than the one from highly continuous IGA. In a mesh with one million elements and , the linear system is solved 15 times faster for the refined than the maximum continuity isogeometric analysis.},
doi = {10.1016/j.cma.2016.08.017},
url = {https://www.osti.gov/biblio/1327763},
journal = {Computer Methods in Applied Mechanics and Engineering},
issn = {0045-7825},
number = ,
volume = ,
place = {United States},
year = {Wed Aug 24 00:00:00 EDT 2016},
month = {Wed Aug 24 00:00:00 EDT 2016}
}
Web of Science
Works referenced in this record:
Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement
journal, October 2005
- Hughes, T. J. R.; Cottrell, J. A.; Bazilevs, Y.
- Computer Methods in Applied Mechanics and Engineering, Vol. 194, Issue 39-41
The cost of continuity: A study of the performance of isogeometric finite elements using direct solvers
journal, March 2012
- Collier, Nathan; Pardo, David; Dalcin, Lisandro
- Computer Methods in Applied Mechanics and Engineering, Vol. 213-216
The Multifrontal Solution of Indefinite Sparse Symmetric Linear
journal, September 1983
- Duff, I. S.; Reid, J. K.
- ACM Transactions on Mathematical Software, Vol. 9, Issue 3
On the computational efficiency of isogeometric methods for smooth elliptic problems using direct solvers: ON THE COMPUTATIONAL EFFICIENCY OF ISOGEOMETRIC METHODS
journal, September 2014
- Collier, N.; Dalcin, L.; Calo, V. M.
- International Journal for Numerical Methods in Engineering, Vol. 100, Issue 8
Computational complexity and memory usage for multi-frontal direct solvers used in p finite element analysis
journal, January 2011
- Calo, Victor M.; Collier, Nathaniel O.; Pardo, David
- Procedia Computer Science, Vol. 4
Isogeometric analysis of the Cahn–Hilliard phase-field model
journal, September 2008
- Gómez, Héctor; Calo, Victor M.; Bazilevs, Yuri
- Computer Methods in Applied Mechanics and Engineering, Vol. 197, Issue 49-50
Phase Field Modeling Using PetIGA
journal, January 2013
- Vignal, Philippe A.; Collier, Nathan; Calo, V. M.
- Procedia Computer Science, Vol. 18
An energy-stable convex splitting for the phase-field crystal equation
journal, October 2015
- Vignal, P.; Dalcin, L.; Brown, D. L.
- Computers & Structures, Vol. 158
Isogeometric Fluid–structure Interaction Analysis with Applications to Arterial Blood Flow
journal, June 2006
- Bazilevs, Y.; Calo, V. M.; Zhang, Y.
- Computational Mechanics, Vol. 38, Issue 4-5
Isogeometric fluid-structure interaction: theory, algorithms, and computations
journal, August 2008
- Bazilevs, Y.; Calo, V. M.; Hughes, T. J. R.
- Computational Mechanics, Vol. 43, Issue 1
An immersogeometric variational framework for fluid–structure interaction: Application to bioprosthetic heart valves
journal, February 2015
- Kamensky, David; Hsu, Ming-Chen; Schillinger, Dominik
- Computer Methods in Applied Mechanics and Engineering, Vol. 284
Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows
journal, December 2007
- Bazilevs, Y.; Calo, V. M.; Cottrell, J. A.
- Computer Methods in Applied Mechanics and Engineering, Vol. 197, Issue 1-4
Isogeometric variational multiscale modeling of wall-bounded turbulent flows with weakly enforced boundary conditions on unstretched meshes
journal, February 2010
- Bazilevs, Y.; Michler, C.; Calo, V. M.
- Computer Methods in Applied Mechanics and Engineering, Vol. 199, Issue 13-16
Isogeometric variational multiscale large-eddy simulation of fully-developed turbulent flow over a wavy wall
journal, September 2012
- Chang, K.; Hughes, T. J. R.; Calo, V. M.
- Computers & Fluids, Vol. 68
The role of continuity in residual-based variational multiscale modeling of turbulence
journal, June 2007
- Akkerman, I.; Bazilevs, Y.; Calo, V. M.
- Computational Mechanics, Vol. 41, Issue 3
Large eddy simulation of turbulent Taylor–Couette flow using isogeometric analysis and the residual-based variational multiscale method
journal, May 2010
- Bazilevs, Y.; Akkerman, I.
- Journal of Computational Physics, Vol. 229, Issue 9
Simulation of laminar and turbulent concentric pipe flows with the isogeometric variational multiscale method
journal, January 2013
- Ghaffari Motlagh, Yousef; Ahn, Hyung Taek; Hughes, Thomas J. R.
- Computers & Fluids, Vol. 71
Discretizations in isogeometric analysis of Navier–Stokes flow
journal, October 2011
- Nielsen, Peter Nørtoft; Gersborg, Allan Roulund; Gravesen, Jens
- Computer Methods in Applied Mechanics and Engineering, Vol. 200, Issue 45-46
Isogeometric Analysis and error estimates for high order partial differential equations in fluid dynamics
journal, October 2014
- Tagliabue, Anna; Dedè, Luca; Quarteroni, Alfio
- Computers & Fluids, Vol. 102
Coupling Navier-stokes and Cahn-hilliard Equations in a Two-dimensional Annular flow Configuration
journal, January 2015
- Vignal, Philippe; Sarmiento, Adel; Côrtes, Adriano M. A.
- Procedia Computer Science, Vol. 51
Modeling Phase-transitions Using a High-performance, Isogeometric Analysis Framework
journal, January 2014
- Vignal, Philippe; Dalcin, Lisandro; Collier, Nathan. O.
- Procedia Computer Science, Vol. 29
Studies of refinement and continuity in isogeometric structural analysis
journal, September 2007
- Cottrell, J. A.; Hughes, T. J. R.; Reali, A.
- Computer Methods in Applied Mechanics and Engineering, Vol. 196, Issue 41-44
Isogeometric analysis of structural vibrations
journal, August 2006
- Cottrell, J. A.; Reali, A.; Bazilevs, Y.
- Computer Methods in Applied Mechanics and Engineering, Vol. 195, Issue 41-43
A fully “locking-free” isogeometric approach for plane linear elasticity problems: A stream function formulation
journal, December 2007
- Auricchio, F.; da Veiga, L. Beirão; Buffa, A.
- Computer Methods in Applied Mechanics and Engineering, Vol. 197, Issue 1-4
Robustness of isogeometric structural discretizations under severe mesh distortion
journal, January 2010
- Lipton, S.; Evans, J. A.; Bazilevs, Y.
- Computer Methods in Applied Mechanics and Engineering, Vol. 199, Issue 5-8
Isogeometric Analysis of Hyperelastic Materials Using PetIGA
journal, January 2013
- Bernal, L. M.; Calo, V. M.; Collier, N.
- Procedia Computer Science, Vol. 18
Isogeometric analysis of the isothermal Navier–Stokes–Korteweg equations
journal, May 2010
- Gomez, Hector; Hughes, Thomas J. R.; Nogueira, Xesús
- Computer Methods in Applied Mechanics and Engineering, Vol. 199, Issue 25-28
IsoGeometric Analysis: Stable elements for the 2D Stokes equation
journal, April 2010
- Buffa, A.; de Falco, C.; Sangalli, G.
- International Journal for Numerical Methods in Fluids, Vol. 65, Issue 11-12
Isogeometric divergence-conforming B-splines for the unsteady Navier–Stokes equations
journal, May 2013
- Evans, John A.; Hughes, Thomas J. R.
- Journal of Computational Physics, Vol. 241
Improving stability of stabilized and multiscale formulations in flow simulations at small time steps
journal, February 2010
- Hsu, M. -C.; Bazilevs, Y.; Calo, V. M.
- Computer Methods in Applied Mechanics and Engineering, Vol. 199, Issue 13-16
Weak Dirichlet boundary conditions for wall-bounded turbulent flows
journal, November 2007
- Bazilevs, Y.; Michler, C.; Calo, V. M.
- Computer Methods in Applied Mechanics and Engineering, Vol. 196, Issue 49-52
Patient-specific vascular NURBS modeling for isogeometric analysis of blood flow
journal, May 2007
- Zhang, Yongjie; Bazilevs, Yuri; Goswami, Samrat
- Computer Methods in Applied Mechanics and Engineering, Vol. 196, Issue 29-30
Patient-specific isogeometric fluid–structure interaction analysis of thoracic aortic blood flow due to implantation of the Jarvik 2000 left ventricular assist device
journal, September 2009
- Bazilevs, Y.; Gohean, J. R.; Hughes, T. J. R.
- Computer Methods in Applied Mechanics and Engineering, Vol. 198, Issue 45-46
Mathematical modeling of coupled drug and drug-encapsulated nanoparticle transport in patient-specific coronary artery walls
journal, August 2011
- Hossain, Shaolie S.; Hossainy, Syed F. A.; Bazilevs, Yuri
- Computational Mechanics, Vol. 49, Issue 2
Multiphysics model for blood flow and drug transport with application to patient-specific coronary artery flow
journal, August 2008
- Calo, V. M.; Brasher, N. F.; Bazilevs, Y.
- Computational Mechanics, Vol. 43, Issue 1
Isogeometric analysis in electromagnetics: B-splines approximation
journal, March 2010
- Buffa, A.; Sangalli, G.; Vázquez, R.
- Computer Methods in Applied Mechanics and Engineering, Vol. 199, Issue 17-20
ISOGEOMETRIC ANALYSIS: APPROXIMATION, STABILITY AND ERROR ESTIMATES FOR h-REFINED MESHES
journal, July 2006
- Bazilevs, Y.; BeirÃO Da Veiga, L.; Cottrell, J. A.
- Mathematical Models and Methods in Applied Sciences, Vol. 16, Issue 07
Nested Dissection of a Regular Finite Element Mesh
journal, April 1973
- George, Alan
- SIAM Journal on Numerical Analysis, Vol. 10, Issue 2
A frontal solution program for finite element analysis
journal, January 1970
- Irons, Bruce M.
- International Journal for Numerical Methods in Engineering, Vol. 2, Issue 1
Performance evaluation of block-diagonal preconditioners for the divergence-conforming B-spline discretization of the Stokes system
journal, November 2015
- Côrtes, A. M. A.; Coutinho, A. L. G. A.; Dalcin, L.
- Journal of Computational Science, Vol. 11
Energy exchange analysis in droplet dynamics via the Navier–Stokes–Cahn–Hilliard model
journal, May 2016
- Espath, L. F. R.; Sarmiento, A. F.; Vignal, P.
- Journal of Fluid Mechanics, Vol. 797
Micropolar Fluids Using B-spline Divergence Conforming Spaces
journal, January 2014
- Sarmiento, Adel; Garcia, Daniel; Dalcin, Lisandro
- Procedia Computer Science, Vol. 29
PetIGA-MF: A multi-field high-performance toolbox for structure-preserving B-splines spaces
journal, January 2017
- Sarmiento, A. F.; Côrtes, A. M. A.; Garcia, D. A.
- Journal of Computational Science, Vol. 18
A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling
journal, January 2001
- Amestoy, Patrick R.; Duff, Iain S.; L'Excellent, Jean-Yves
- SIAM Journal on Matrix Analysis and Applications, Vol. 23, Issue 1
Hybrid scheduling for the parallel solution of linear systems
journal, February 2006
- Amestoy, Patrick R.; Guermouche, Abdou; L’Excellent, Jean-Yves
- Parallel Computing, Vol. 32, Issue 2
Works referencing / citing this record:
Dispersion-optimized quadrature rules for isogeometric analysis: Modified inner products, their dispersion properties, and optimally blended schemes
journal, June 2017
- Puzyrev, Vladimir; Deng, Quanling; Calo, Victor
- Computer Methods in Applied Mechanics and Engineering, Vol. 320