Computing bounds for linear functionals of exact weak solutions to Poisson's equation

We present a method for Poisson’s equation that computes guaranteed upper and lower bounds for the values of piecewise-polynomial linear functional outputs of the exact weak solution of the infinite-dimensional continuum problem with piecewise-polynomial forcing. The method results from exploiting t...

Descripción completa

Detalles Bibliográficos
Autores: Sauer-Budge, A. M., Bonet Carbonell, Javier|||0000-0002-0430-5181, Huerta, Antonio|||0000-0003-4198-3798, Peraire Guitart, Jaume
Tipo de recurso: artículo
Fecha de publicación:2004
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:upcommons.upc.edu:2117/7997
Acceso en línea:https://hdl.handle.net/2117/7997
https://dx.doi.org/10.1137/S0036142903425045
Access Level:acceso abierto
Palabra clave:Poisson algebras
Poisson, Equació de
Àrees temàtiques de la UPC::Matemàtiques i estadística::Àlgebra
id ES_52e9f6e61395ada36e6b4c8e99a7c526
oai_identifier_str oai:upcommons.upc.edu:2117/7997
network_acronym_str ES
network_name_str España
repository_id_str
spelling Computing bounds for linear functionals of exact weak solutions to Poisson's equationSauer-Budge, A. M.Bonet Carbonell, Javier|||0000-0002-0430-5181Huerta, Antonio|||0000-0003-4198-3798Peraire Guitart, JaumePoisson algebrasPoisson, Equació deÀrees temàtiques de la UPC::Matemàtiques i estadística::ÀlgebraWe present a method for Poisson’s equation that computes guaranteed upper and lower bounds for the values of piecewise-polynomial linear functional outputs of the exact weak solution of the infinite-dimensional continuum problem with piecewise-polynomial forcing. The method results from exploiting the Lagrangian saddle point property engendered by recasting the output problem as a constrained minimization problem. Localization is achieved by Lagrangian relaxation and the bounds are computed by appeal to a local dual problem. The proposed method computes approximate Lagrange multipliers using traditional finite element approximations to calculate a primal and an adjoint solution along with well known hybridization techniques to calculate interelement continuity multipliers. The computed bounds hold uniformly for any level of refinement, and in the asymptotic convergence regime of the finite element method, the bound gap decreases at twice the rate of the energy norm measure of the error in the finite element solution. Given a finite element solution and its output adjoint solution, the method can be used to provide a certificate of precision for the output with an asymptotic complexity that is linear in the number of elements in the finite element discretization. The elemental contributions to the bound gap are always positive and hence lend themselves to be used as adaptive indicators, as we demonstrate with a numerical example.Peer Reviewed20042004-01-0120102010-07-02journal articlehttp://purl.org/coar/resource_type/c_6501VoRhttp://purl.org/coar/version/c_970fb48d4fbd8a85info:eu-repo/semantics/articleapplication/pdfhttps://hdl.handle.net/2117/7997https://dx.doi.org/10.1137/S0036142903425045reponame:UPCommons. Portal del coneixement obert de la UPCinstname:Universitat Politècnica de Catalunya (UPC)Inglésengopen accesshttp://purl.org/coar/access_right/c_abf2info:eu-repo/semantics/openAccessoai:upcommons.upc.edu:2117/79972026-05-27T15:37:01Z
dc.title.none.fl_str_mv Computing bounds for linear functionals of exact weak solutions to Poisson's equation
title Computing bounds for linear functionals of exact weak solutions to Poisson's equation
spellingShingle Computing bounds for linear functionals of exact weak solutions to Poisson's equation
Sauer-Budge, A. M.
Poisson algebras
Poisson, Equació de
Àrees temàtiques de la UPC::Matemàtiques i estadística::Àlgebra
title_short Computing bounds for linear functionals of exact weak solutions to Poisson's equation
title_full Computing bounds for linear functionals of exact weak solutions to Poisson's equation
title_fullStr Computing bounds for linear functionals of exact weak solutions to Poisson's equation
title_full_unstemmed Computing bounds for linear functionals of exact weak solutions to Poisson's equation
title_sort Computing bounds for linear functionals of exact weak solutions to Poisson's equation
dc.creator.none.fl_str_mv Sauer-Budge, A. M.
Bonet Carbonell, Javier|||0000-0002-0430-5181
Huerta, Antonio|||0000-0003-4198-3798
Peraire Guitart, Jaume
author Sauer-Budge, A. M.
author_facet Sauer-Budge, A. M.
Bonet Carbonell, Javier|||0000-0002-0430-5181
Huerta, Antonio|||0000-0003-4198-3798
Peraire Guitart, Jaume
author_role author
author2 Bonet Carbonell, Javier|||0000-0002-0430-5181
Huerta, Antonio|||0000-0003-4198-3798
Peraire Guitart, Jaume
author2_role author
author
author
dc.subject.none.fl_str_mv Poisson algebras
Poisson, Equació de
Àrees temàtiques de la UPC::Matemàtiques i estadística::Àlgebra
topic Poisson algebras
Poisson, Equació de
Àrees temàtiques de la UPC::Matemàtiques i estadística::Àlgebra
description We present a method for Poisson’s equation that computes guaranteed upper and lower bounds for the values of piecewise-polynomial linear functional outputs of the exact weak solution of the infinite-dimensional continuum problem with piecewise-polynomial forcing. The method results from exploiting the Lagrangian saddle point property engendered by recasting the output problem as a constrained minimization problem. Localization is achieved by Lagrangian relaxation and the bounds are computed by appeal to a local dual problem. The proposed method computes approximate Lagrange multipliers using traditional finite element approximations to calculate a primal and an adjoint solution along with well known hybridization techniques to calculate interelement continuity multipliers. The computed bounds hold uniformly for any level of refinement, and in the asymptotic convergence regime of the finite element method, the bound gap decreases at twice the rate of the energy norm measure of the error in the finite element solution. Given a finite element solution and its output adjoint solution, the method can be used to provide a certificate of precision for the output with an asymptotic complexity that is linear in the number of elements in the finite element discretization. The elemental contributions to the bound gap are always positive and hence lend themselves to be used as adaptive indicators, as we demonstrate with a numerical example.
publishDate 2004
dc.date.none.fl_str_mv 2004
2004-01-01
2010
2010-07-02
dc.type.none.fl_str_mv journal article
http://purl.org/coar/resource_type/c_6501
VoR
http://purl.org/coar/version/c_970fb48d4fbd8a85
dc.type.openaire.fl_str_mv info:eu-repo/semantics/article
format article
dc.identifier.none.fl_str_mv https://hdl.handle.net/2117/7997
https://dx.doi.org/10.1137/S0036142903425045
url https://hdl.handle.net/2117/7997
https://dx.doi.org/10.1137/S0036142903425045
dc.language.none.fl_str_mv Inglés
eng
language_invalid_str_mv Inglés
language eng
dc.rights.none.fl_str_mv open access
http://purl.org/coar/access_right/c_abf2
dc.rights.openaire.fl_str_mv info:eu-repo/semantics/openAccess
rights_invalid_str_mv open access
http://purl.org/coar/access_right/c_abf2
eu_rights_str_mv openAccess
dc.format.none.fl_str_mv application/pdf
dc.source.none.fl_str_mv reponame:UPCommons. Portal del coneixement obert de la UPC
instname:Universitat Politècnica de Catalunya (UPC)
instname_str Universitat Politècnica de Catalunya (UPC)
reponame_str UPCommons. Portal del coneixement obert de la UPC
collection UPCommons. Portal del coneixement obert de la UPC
repository.name.fl_str_mv
repository.mail.fl_str_mv
_version_ 1869408081399513088
score 15,301603