Goal-Oriented Adaptivity using Unconventional Error Representations

In Goal-Oriented Adaptivity (GOA), the error in a Quantity of Interest (QoI) is represented using global error functions of the direct and adjoint problems. This error representation is subsequently bounded above by element-wise error indicators that are used to drive optimal refinements. In this wo...

Descripción completa

Detalles Bibliográficos
Autor: Darrigrand, V.
Tipo de recurso: tesis doctoral
Estado:Versión aceptada para publicación
Fecha de publicación:2017
País:España
Institución:Basque Center for Applied Mathematics (BCAM)
Repositorio:BIRD. BCAM's Institutional Repository Data
OAI Identifier:oai:bird.bcamath.org:20.500.11824/735
Acceso en línea:http://hdl.handle.net/20.500.11824/735
Access Level:acceso abierto
Palabra clave:Goal-Oriented Adaptivity
Error representation
id ES_552bcbb2a65dcb8fa3e71e57c8a9dca8
oai_identifier_str oai:bird.bcamath.org:20.500.11824/735
network_acronym_str ES
network_name_str España
repository_id_str
spelling Goal-Oriented Adaptivity using Unconventional Error RepresentationsDarrigrand, V.Goal-Oriented AdaptivityError representationIn Goal-Oriented Adaptivity (GOA), the error in a Quantity of Interest (QoI) is represented using global error functions of the direct and adjoint problems. This error representation is subsequently bounded above by element-wise error indicators that are used to drive optimal refinements. In this work, we propose to replace, in the error representation, the adjoint problem by an alternative operator. The main advantage of the proposed approach is that, when judiciously selecting such alternative operator, the corresponding upper bound of the error representation becomes sharper, leading to a more efficient GOA. These representations can be employed to design novel h, p, and hp energy-norm and goal-oriented adaptive algorithms. While the method can be applied to a variety of problems, in this Dissertation we first focus on one-dimensional (1D) problems, including Helmholtz and steady-state convection-dominated diffusion problems. Numerical results in 1D show that for the Helmholtz problem, it is advantageous to select the Laplace operator for the alternative error representation. Specifically, the upper bounds of the new error representation are sharper than the classical ones used in both energy-norm and goal-oriented adaptive methods, especially when the dispersion (pollution) error is significant. The 1D steady-state convection-dominated diffusion problem with homogeneous Dirichlet boundary conditions exhibits a boundary layer that produces a loss of numerical stability. The new error representation based on the Laplace operator delivers sharper error upper bounds. When applied to a p-GOA, the alternative error representation captures earlier the boundary layer, despite the existing spurious numerical oscillations. We then focus on the two- and three-dimensional (2D and 3D) Helmholtz equation. We show via extensive numerical experimentation that the upper bounds provided by the alternative error representations are sharper than the classical ones. When using the alternative error indicators, a naive p-adaptive process converges, whereas under the same conditions, the classical method fails and requires the use of the so-called Projection Based Interpolation (PBI) operator or some other technique to regain convergence. We also provide guidelines for finding operators delivering sharp error representation upper bounds.Basque Government Consolidated Research Group Grant IT649-13201720172017info:eu-repo/semantics/doctoralThesisinfo:eu-repo/semantics/acceptedVersionapplication/pdfhttp://hdl.handle.net/20.500.11824/735reponame:BIRD. BCAM's Institutional Repository Datainstname:Basque Center for Applied Mathematics (BCAM)Inglésinfo:eu-repo/grantAgreement/EC/H2020/644202info:eu-repo/grantAgreement/MINECO//SEV-2013-0323info:eu-repo/grantAgreement/MINECO//MTM2016-76329-Rinfo:eu-repo/grantAgreement/MINECO//MTM2016-81697-ERCinfo:eu-repo/grantAgreement/MINECO//MTM2013-40824-Pinfo:eu-repo/grantAgreement/Gobierno Vasco/BERC/BERC.2014-2017Reconocimiento-NoComercial-CompartirIgual 3.0 Españahttp://creativecommons.org/licenses/by-nc-sa/3.0/es/info:eu-repo/semantics/openAccessoai:bird.bcamath.org:20.500.11824/7352026-06-19T12:47:47Z
dc.title.none.fl_str_mv Goal-Oriented Adaptivity using Unconventional Error Representations
title Goal-Oriented Adaptivity using Unconventional Error Representations
spellingShingle Goal-Oriented Adaptivity using Unconventional Error Representations
Darrigrand, V.
Goal-Oriented Adaptivity
Error representation
title_short Goal-Oriented Adaptivity using Unconventional Error Representations
title_full Goal-Oriented Adaptivity using Unconventional Error Representations
title_fullStr Goal-Oriented Adaptivity using Unconventional Error Representations
title_full_unstemmed Goal-Oriented Adaptivity using Unconventional Error Representations
title_sort Goal-Oriented Adaptivity using Unconventional Error Representations
dc.creator.none.fl_str_mv Darrigrand, V.
author Darrigrand, V.
author_facet Darrigrand, V.
author_role author
dc.subject.none.fl_str_mv Goal-Oriented Adaptivity
Error representation
topic Goal-Oriented Adaptivity
Error representation
description In Goal-Oriented Adaptivity (GOA), the error in a Quantity of Interest (QoI) is represented using global error functions of the direct and adjoint problems. This error representation is subsequently bounded above by element-wise error indicators that are used to drive optimal refinements. In this work, we propose to replace, in the error representation, the adjoint problem by an alternative operator. The main advantage of the proposed approach is that, when judiciously selecting such alternative operator, the corresponding upper bound of the error representation becomes sharper, leading to a more efficient GOA. These representations can be employed to design novel h, p, and hp energy-norm and goal-oriented adaptive algorithms. While the method can be applied to a variety of problems, in this Dissertation we first focus on one-dimensional (1D) problems, including Helmholtz and steady-state convection-dominated diffusion problems. Numerical results in 1D show that for the Helmholtz problem, it is advantageous to select the Laplace operator for the alternative error representation. Specifically, the upper bounds of the new error representation are sharper than the classical ones used in both energy-norm and goal-oriented adaptive methods, especially when the dispersion (pollution) error is significant. The 1D steady-state convection-dominated diffusion problem with homogeneous Dirichlet boundary conditions exhibits a boundary layer that produces a loss of numerical stability. The new error representation based on the Laplace operator delivers sharper error upper bounds. When applied to a p-GOA, the alternative error representation captures earlier the boundary layer, despite the existing spurious numerical oscillations. We then focus on the two- and three-dimensional (2D and 3D) Helmholtz equation. We show via extensive numerical experimentation that the upper bounds provided by the alternative error representations are sharper than the classical ones. When using the alternative error indicators, a naive p-adaptive process converges, whereas under the same conditions, the classical method fails and requires the use of the so-called Projection Based Interpolation (PBI) operator or some other technique to regain convergence. We also provide guidelines for finding operators delivering sharp error representation upper bounds.
publishDate 2017
dc.date.none.fl_str_mv 2017
2017
2017
dc.type.none.fl_str_mv info:eu-repo/semantics/doctoralThesis
info:eu-repo/semantics/acceptedVersion
format doctoralThesis
status_str acceptedVersion
dc.identifier.none.fl_str_mv http://hdl.handle.net/20.500.11824/735
url http://hdl.handle.net/20.500.11824/735
dc.language.none.fl_str_mv Inglés
language_invalid_str_mv Inglés
dc.relation.none.fl_str_mv info:eu-repo/grantAgreement/EC/H2020/644202
info:eu-repo/grantAgreement/MINECO//SEV-2013-0323
info:eu-repo/grantAgreement/MINECO//MTM2016-76329-R
info:eu-repo/grantAgreement/MINECO//MTM2016-81697-ERC
info:eu-repo/grantAgreement/MINECO//MTM2013-40824-P
info:eu-repo/grantAgreement/Gobierno Vasco/BERC/BERC.2014-2017
dc.rights.none.fl_str_mv Reconocimiento-NoComercial-CompartirIgual 3.0 España
http://creativecommons.org/licenses/by-nc-sa/3.0/es/
info:eu-repo/semantics/openAccess
rights_invalid_str_mv Reconocimiento-NoComercial-CompartirIgual 3.0 España
http://creativecommons.org/licenses/by-nc-sa/3.0/es/
eu_rights_str_mv openAccess
dc.format.none.fl_str_mv application/pdf
dc.source.none.fl_str_mv reponame:BIRD. BCAM's Institutional Repository Data
instname:Basque Center for Applied Mathematics (BCAM)
instname_str Basque Center for Applied Mathematics (BCAM)
reponame_str BIRD. BCAM's Institutional Repository Data
collection BIRD. BCAM's Institutional Repository Data
repository.name.fl_str_mv
repository.mail.fl_str_mv
_version_ 1869408252035334144
score 15.300724