The a posteriori finite element method (APFEM), a framework for efficient parametric study and Bayesian inferences

Stochastic methods have recently been the subject of increased attention in Computational Mechanics for their ability to account for the stochasticity of both material parameters and geometrical features in their predictions. Among them, the Galerkin Stochastic Finite Element Method (GSFEM) was show...

Full description

Bibliographic Details
Main Authors: Ammouche, Y, Jérusalem, A
Format: Journal article
Language:English
Published: Elsevier 2023
_version_ 1826310463200165888
author Ammouche, Y
Jérusalem, A
author_facet Ammouche, Y
Jérusalem, A
author_sort Ammouche, Y
collection OXFORD
description Stochastic methods have recently been the subject of increased attention in Computational Mechanics for their ability to account for the stochasticity of both material parameters and geometrical features in their predictions. Among them, the Galerkin Stochastic Finite Element Method (GSFEM) was shown to be particularly efficient and able to provide accurate output statistics, although at the cost of intrusive coding and additional theoretical algebraic efforts. In this method, distributions of the stochastic parameters are used as inputs for the solver, which in turn outputs nodal displacement distributions in one simulation. Here, we propose an extension of the GSFEM—termed the A posteriori Finite Element Method or APFEM—where uniform distributions are taken by default to allow for parametric studies of the inputs of interest as a postprocessing step after the simulation. Doing so, APFEM only requires the knowledge of the vertices of the parameter space. In particular, one key advantage of APFEM is its use in the context of Bayesian inferences, where the random evaluations required by the Bayesian setting (usually done through Monte Carlo) can be done exactly without the need for further simulations. Finally, we demonstrate the potential of APFEM by solving forward models with parametric boundary conditions in the context of (i) metamaterial design and (ii) pitchfork bifurcation of the buckling of a slender structure; and demonstrate the flexibility of its use for Bayesian inference by (iii) inferring friction coefficient of a half plane in a contact mechanics problem and (iv) inferring the stiffness of a brain region in the context of cancer surgical planning.
first_indexed 2024-03-07T07:50:48Z
format Journal article
id oxford-uuid:f90f5372-007f-4cdc-99b0-c1f1168d7791
institution University of Oxford
language English
last_indexed 2024-03-07T07:50:48Z
publishDate 2023
publisher Elsevier
record_format dspace
spelling oxford-uuid:f90f5372-007f-4cdc-99b0-c1f1168d77912023-07-12T08:25:54ZThe a posteriori finite element method (APFEM), a framework for efficient parametric study and Bayesian inferencesJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:f90f5372-007f-4cdc-99b0-c1f1168d7791EnglishSymplectic ElementsElsevier2023Ammouche, YJérusalem, AStochastic methods have recently been the subject of increased attention in Computational Mechanics for their ability to account for the stochasticity of both material parameters and geometrical features in their predictions. Among them, the Galerkin Stochastic Finite Element Method (GSFEM) was shown to be particularly efficient and able to provide accurate output statistics, although at the cost of intrusive coding and additional theoretical algebraic efforts. In this method, distributions of the stochastic parameters are used as inputs for the solver, which in turn outputs nodal displacement distributions in one simulation. Here, we propose an extension of the GSFEM—termed the A posteriori Finite Element Method or APFEM—where uniform distributions are taken by default to allow for parametric studies of the inputs of interest as a postprocessing step after the simulation. Doing so, APFEM only requires the knowledge of the vertices of the parameter space. In particular, one key advantage of APFEM is its use in the context of Bayesian inferences, where the random evaluations required by the Bayesian setting (usually done through Monte Carlo) can be done exactly without the need for further simulations. Finally, we demonstrate the potential of APFEM by solving forward models with parametric boundary conditions in the context of (i) metamaterial design and (ii) pitchfork bifurcation of the buckling of a slender structure; and demonstrate the flexibility of its use for Bayesian inference by (iii) inferring friction coefficient of a half plane in a contact mechanics problem and (iv) inferring the stiffness of a brain region in the context of cancer surgical planning.
spellingShingle Ammouche, Y
Jérusalem, A
The a posteriori finite element method (APFEM), a framework for efficient parametric study and Bayesian inferences
title The a posteriori finite element method (APFEM), a framework for efficient parametric study and Bayesian inferences
title_full The a posteriori finite element method (APFEM), a framework for efficient parametric study and Bayesian inferences
title_fullStr The a posteriori finite element method (APFEM), a framework for efficient parametric study and Bayesian inferences
title_full_unstemmed The a posteriori finite element method (APFEM), a framework for efficient parametric study and Bayesian inferences
title_short The a posteriori finite element method (APFEM), a framework for efficient parametric study and Bayesian inferences
title_sort a posteriori finite element method apfem a framework for efficient parametric study and bayesian inferences
work_keys_str_mv AT ammouchey theaposteriorifiniteelementmethodapfemaframeworkforefficientparametricstudyandbayesianinferences
AT jerusalema theaposteriorifiniteelementmethodapfemaframeworkforefficientparametricstudyandbayesianinferences
AT ammouchey aposteriorifiniteelementmethodapfemaframeworkforefficientparametricstudyandbayesianinferences
AT jerusalema aposteriorifiniteelementmethodapfemaframeworkforefficientparametricstudyandbayesianinferences