The statistical finite element method (statFEM) for coherent synthesis of observation data and model predictions
arXiv:1905.06391 · doi:10.1016/j.cma.2020.113533
Abstract
The increased availability of observation data from engineering systems in operation poses the question of how to incorporate this data into finite element models. To this end, we propose a novel statistical construction of the finite element method that provides the means of synthesising measurement data and finite element models. The Bayesian statistical framework is adopted to treat all the uncertainties present in the data, the mathematical model and its finite element discretisation. From the outset, we postulate a data-generating model which additively decomposes data into a finite element, a model misspecification and a noise component. Each of the components may be uncertain and is considered as a random variable with a respective prior probability density. The prior of the finite element component is given by a conventional stochastic forward problem. The prior probabilities of the model misspecification and measurement noise, without loss of generality, are assumed to have zero-mean and known covariance structure. Our proposed statistical model is hierarchical in the sense that each of the three random components may depend on non-observable random hyperparameters. Because of the hierarchical structure of the statistical model, Bayes rule is applied on three different levels in turn to infer the posterior densities of the three random components and hyperparameters. On level one, we determine the posterior densities of the finite element component and the true system response using the prior finite element density given by the forward problem and the data likelihood. On the next level, we infer the hyperparameter posterior densities from their respective priors and the marginal likelihood of the first inference problem. Finally, on level three we use Bayes rule to choose the most suitable finite element model in light of the observed data by computing the model posteriors.
References in corpus (3)
- Geometric MCMC for Infinite-Dimensional Inverse Problems
- A computational framework for infinite-dimensional Bayesian inverse problems. Part I: The linearized case, with application to global seismic inversion
- The synthesis of data from instrumented structures and physics-based models via Gaussian processes
Cited by in corpus (19)
- Reduced and All-at-Once Approaches for Model Calibration and Discovery in Computational Solid Mechanics
- On generative models as the basis for digital twins
- Fully probabilistic deep models for forward and inverse problems in parametric PDEs
- Variational Bayesian Approximation of Inverse Problems using Sparse Precision Matrices
- A probabilistic finite element method based on random meshes: Error estimators and Bayesian inverse problems
- Low-rank statistical finite elements for scalable model-data synthesis
- A Bayesian constitutive model selection framework for biaxial mechanical testing of planar soft tissues: application to porcine aortic valves
- Inferring Displacement Fields from Sparse Measurements Using the Statistical Finite Element Method
- Variational Inference for Nonlinear Inverse Problems via Neural Net Kernels: Comparison to Bayesian Neural Networks, Application to Topology Optimization
- Variational Bayesian surrogate modelling with application to robust design optimisation
- Radiative Transfer as a Bayesian Linear Regression problem
- Stochastic Inference of Plate Bending from Heterogeneous Data: Physics-informed Gaussian Processes via Kirchhoff-Love Theory
- Probabilistic Linear Solvers for Machine Learning
- Probabilistic Numerical Method of Lines for Time-Dependent Partial Differential Equations
- Unsupervised Constitutive Model Discovery from Sparse and Noisy Data
- Statistical reduced order modelling for the parametric Helmholtz equation
- Active Learning with Adaptive Non-Stationary Kernel for Continuous-Fidelity Surrogate Models
- A Regularized Ensemble Kalman Filter for Stochastic Phase Field Models of Brittle Fracture
- Sampling Methods for Bayesian Inference Involving Convergent Noisy Approximations of Forward Maps