Recovery-based error estimation and adaptivity using high-order splines over hierarchical T-meshes

Cosmin Anitescu, Md Naim Hossain, Timon Rabczuk

Research output: Contribution to journalArticle

33 Citations (Scopus)

Abstract

An adaptive higher-order method based on a generalization of polynomial/rational splines over hierarchical T-meshes (PHT/RHT-splines) is introduced. While most problems considered in isogeometric analysis can be solved efficiently when the solution is smooth, many non-trivial simulations have rough solutions. This can be caused, for example, by the presence of re-entrant corners in the domain. For such problems, a tensor-product basis is less suitable for resolving the singularities that appear, as refinement propagates throughout the computational domain. Hierarchical bases and adaptivity allow for a more efficient way of dealing with singularities, by adding more degrees of freedom only where they are necessary to improve the approximation. In order to drive the adaptive refinement, an efficient recovery-based error estimator is proposed in this work. The estimator produces a “recovered solution” which is a more accurate approximation than the computed numerical solution. Several 2D and 3D numerical investigations with PHT-splines of higher order and greater continuity show good performance compared to uniform refinement in terms of degrees of freedom and computational cost.

Original languageEnglish
Pages (from-to)638-662
Number of pages25
JournalComputer Methods in Applied Mechanics and Engineering
Volume328
DOIs
Publication statusPublished - 2018 Jan 1
Externally publishedYes

Fingerprint

splines
Splines
Error analysis
mesh
recovery
estimators
Recovery
degrees of freedom
approximation
continuity
Tensors
polynomials
Polynomials
tensors
costs
products
Costs
simulation

Keywords

  • Error estimation
  • Hierarchical splines
  • Isogeometric analysis
  • Linear elasticity
  • Superconvergent patch recovery

ASJC Scopus subject areas

  • Computational Mechanics
  • Mechanics of Materials
  • Mechanical Engineering
  • Physics and Astronomy(all)
  • Computer Science Applications

Cite this

Recovery-based error estimation and adaptivity using high-order splines over hierarchical T-meshes. / Anitescu, Cosmin; Hossain, Md Naim; Rabczuk, Timon.

In: Computer Methods in Applied Mechanics and Engineering, Vol. 328, 01.01.2018, p. 638-662.

Research output: Contribution to journalArticle

@article{6181e91ab29b4594b7f2bf0512c72715,
title = "Recovery-based error estimation and adaptivity using high-order splines over hierarchical T-meshes",
abstract = "An adaptive higher-order method based on a generalization of polynomial/rational splines over hierarchical T-meshes (PHT/RHT-splines) is introduced. While most problems considered in isogeometric analysis can be solved efficiently when the solution is smooth, many non-trivial simulations have rough solutions. This can be caused, for example, by the presence of re-entrant corners in the domain. For such problems, a tensor-product basis is less suitable for resolving the singularities that appear, as refinement propagates throughout the computational domain. Hierarchical bases and adaptivity allow for a more efficient way of dealing with singularities, by adding more degrees of freedom only where they are necessary to improve the approximation. In order to drive the adaptive refinement, an efficient recovery-based error estimator is proposed in this work. The estimator produces a “recovered solution” which is a more accurate approximation than the computed numerical solution. Several 2D and 3D numerical investigations with PHT-splines of higher order and greater continuity show good performance compared to uniform refinement in terms of degrees of freedom and computational cost.",
keywords = "Error estimation, Hierarchical splines, Isogeometric analysis, Linear elasticity, Superconvergent patch recovery",
author = "Cosmin Anitescu and Hossain, {Md Naim} and Timon Rabczuk",
year = "2018",
month = "1",
day = "1",
doi = "10.1016/j.cma.2017.08.032",
language = "English",
volume = "328",
pages = "638--662",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0045-7825",
publisher = "Elsevier",

}

TY - JOUR

T1 - Recovery-based error estimation and adaptivity using high-order splines over hierarchical T-meshes

AU - Anitescu, Cosmin

AU - Hossain, Md Naim

AU - Rabczuk, Timon

PY - 2018/1/1

Y1 - 2018/1/1

N2 - An adaptive higher-order method based on a generalization of polynomial/rational splines over hierarchical T-meshes (PHT/RHT-splines) is introduced. While most problems considered in isogeometric analysis can be solved efficiently when the solution is smooth, many non-trivial simulations have rough solutions. This can be caused, for example, by the presence of re-entrant corners in the domain. For such problems, a tensor-product basis is less suitable for resolving the singularities that appear, as refinement propagates throughout the computational domain. Hierarchical bases and adaptivity allow for a more efficient way of dealing with singularities, by adding more degrees of freedom only where they are necessary to improve the approximation. In order to drive the adaptive refinement, an efficient recovery-based error estimator is proposed in this work. The estimator produces a “recovered solution” which is a more accurate approximation than the computed numerical solution. Several 2D and 3D numerical investigations with PHT-splines of higher order and greater continuity show good performance compared to uniform refinement in terms of degrees of freedom and computational cost.

AB - An adaptive higher-order method based on a generalization of polynomial/rational splines over hierarchical T-meshes (PHT/RHT-splines) is introduced. While most problems considered in isogeometric analysis can be solved efficiently when the solution is smooth, many non-trivial simulations have rough solutions. This can be caused, for example, by the presence of re-entrant corners in the domain. For such problems, a tensor-product basis is less suitable for resolving the singularities that appear, as refinement propagates throughout the computational domain. Hierarchical bases and adaptivity allow for a more efficient way of dealing with singularities, by adding more degrees of freedom only where they are necessary to improve the approximation. In order to drive the adaptive refinement, an efficient recovery-based error estimator is proposed in this work. The estimator produces a “recovered solution” which is a more accurate approximation than the computed numerical solution. Several 2D and 3D numerical investigations with PHT-splines of higher order and greater continuity show good performance compared to uniform refinement in terms of degrees of freedom and computational cost.

KW - Error estimation

KW - Hierarchical splines

KW - Isogeometric analysis

KW - Linear elasticity

KW - Superconvergent patch recovery

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

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

U2 - 10.1016/j.cma.2017.08.032

DO - 10.1016/j.cma.2017.08.032

M3 - Article

AN - SCOPUS:85032879540

VL - 328

SP - 638

EP - 662

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0045-7825

ER -