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

Journal Article · · Computer Methods in Applied Mechanics and Engineering
 [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)

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.

Research Organization:
Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
Sponsoring Organization:
USDOE
Grant/Contract Number:
AC05-00OR22725
OSTI ID:
1327763
Alternate ID(s):
OSTI ID: 1429198
Journal Information:
Computer Methods in Applied Mechanics and Engineering, Journal Name: Computer Methods in Applied Mechanics and Engineering; ISSN 0045-7825
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 23 works
Citation information provided by
Web of Science

References (45)

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
Isogeometric variational multiscale large-eddy simulation of fully-developed turbulent flow over a wavy wall journal September 2012
The role of continuity in residual-based variational multiscale modeling of turbulence journal June 2007
Large eddy simulation of turbulent Taylor–Couette flow using isogeometric analysis and the residual-based variational multiscale method journal May 2010
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

Cited By (1)

Dispersion-optimized quadrature rules for isogeometric analysis: Modified inner products, their dispersion properties, and optimally blended schemes journal June 2017