skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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 » the maximum continuity isogeometric analysis.« less

Authors:
 [1];  [2]; ORCiD logo [3];  [4];  [5];  [6]
  1. Basque Center for Applied Mathematics (BCAM), Bilbao (Spain)
  2. Basque Center for Applied Mathematics (BCAM), Bilbao (Spain); Univ. of the Basque Country UPV.EHU, Leioa (Spain); Ikerbasque (Basque Foundation for Sciences), Bilbao (Spain)
  3. 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)
  4. AGH Univ. of Science and Technology, Krakow (Poland)
  5. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  6. 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}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 23 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement
journal, October 2005


The cost of continuity: A study of the performance of isogeometric finite elements using direct solvers
journal, March 2012


The Multifrontal Solution of Indefinite Sparse Symmetric Linear
journal, September 1983


On the computational efficiency of isogeometric methods for smooth elliptic problems using direct solvers: ON THE COMPUTATIONAL EFFICIENCY OF ISOGEOMETRIC METHODS
journal, September 2014


Computational complexity and memory usage for multi-frontal direct solvers used in p finite element analysis
journal, January 2011


Isogeometric analysis of the Cahn–Hilliard phase-field model
journal, September 2008


Phase Field Modeling Using PetIGA
journal, January 2013


An energy-stable convex splitting for the phase-field crystal equation
journal, October 2015


Isogeometric Fluid–structure Interaction Analysis with Applications to Arterial Blood Flow
journal, June 2006


Isogeometric fluid-structure interaction: theory, algorithms, and computations
journal, August 2008


An immersogeometric variational framework for fluid–structure interaction: Application to bioprosthetic heart valves
journal, February 2015


Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows
journal, December 2007


Isogeometric variational multiscale modeling of wall-bounded turbulent flows with weakly enforced boundary conditions on unstretched meshes
journal, February 2010


The role of continuity in residual-based variational multiscale modeling of turbulence
journal, June 2007


Simulation of laminar and turbulent concentric pipe flows with the isogeometric variational multiscale method
journal, January 2013


Discretizations in isogeometric analysis of Navier–Stokes flow
journal, October 2011


Isogeometric Analysis and error estimates for high order partial differential equations in fluid dynamics
journal, October 2014


Coupling Navier-stokes and Cahn-hilliard Equations in a Two-dimensional Annular flow Configuration
journal, January 2015


Modeling Phase-transitions Using a High-performance, Isogeometric Analysis Framework
journal, January 2014


Studies of refinement and continuity in isogeometric structural analysis
journal, September 2007


Isogeometric analysis of structural vibrations
journal, August 2006


A fully “locking-free” isogeometric approach for plane linear elasticity problems: A stream function formulation
journal, December 2007


Robustness of isogeometric structural discretizations under severe mesh distortion
journal, January 2010


Isogeometric Analysis of Hyperelastic Materials Using PetIGA
journal, January 2013


Isogeometric analysis of the isothermal Navier–Stokes–Korteweg equations
journal, May 2010


IsoGeometric Analysis: Stable elements for the 2D Stokes equation
journal, April 2010


Isogeometric divergence-conforming B-splines for the unsteady Navier–Stokes equations
journal, May 2013


Improving stability of stabilized and multiscale formulations in flow simulations at small time steps
journal, February 2010


Weak Dirichlet boundary conditions for wall-bounded turbulent flows
journal, November 2007


Patient-specific vascular NURBS modeling for isogeometric analysis of blood flow
journal, May 2007


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


Mathematical modeling of coupled drug and drug-encapsulated nanoparticle transport in patient-specific coronary artery walls
journal, August 2011


Multiphysics model for blood flow and drug transport with application to patient-specific coronary artery flow
journal, August 2008


Isogeometric analysis in electromagnetics: B-splines approximation
journal, March 2010


ISOGEOMETRIC ANALYSIS: APPROXIMATION, STABILITY AND ERROR ESTIMATES FOR h-REFINED MESHES
journal, July 2006


Nested Dissection of a Regular Finite Element Mesh
journal, April 1973


A frontal solution program for finite element analysis
journal, January 1970


Performance evaluation of block-diagonal preconditioners for the divergence-conforming B-spline discretization of the Stokes system
journal, November 2015


Energy exchange analysis in droplet dynamics via the Navier–Stokes–Cahn–Hilliard model
journal, May 2016


Micropolar Fluids Using B-spline Divergence Conforming Spaces
journal, January 2014


PetIGA-MF: A multi-field high-performance toolbox for structure-preserving B-splines spaces
journal, January 2017


A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling
journal, January 2001


Hybrid scheduling for the parallel solution of linear systems
journal, February 2006


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