Statistical finite elements for sequential data synthesis in solid dynamics
Pith reviewed 2026-05-10 14:49 UTC · model grok-4.3
The pith
A Bayesian filtering extension to statistical finite elements assimilates observational data sequentially into elastodynamic models by advancing the discretized state with a stochastic Newmark scheme.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Observational data are assimilated while the state of the spatially discretised finite element problem is advanced in time using the stochastic variant of the explicit Newmark scheme. The prior probability density of the state is obtained by solving an incremental probabilistic forward problem in which spatio-temporal Gaussian random fields for forcing, model misspecification and material parameters are specified via their stochastic PDE formulation. The resulting non-Gaussian prior is approximated by a perturbation approach, producing a Gaussian posterior with closed-form mean and covariance at every time step. The hyperparameters of the misspecification random field are calibrated by maxim
What carries the argument
The perturbation approximation that converts the non-Gaussian prior density arising from the stochastic PDEs into a Gaussian form, enabling closed-form Bayesian updates of the finite-element state within the statistical finite element framework.
If this is right
- At each time step the updated state estimate is Gaussian with an explicit mean vector and covariance matrix that incorporate both model uncertainty and fresh observations.
- Hyperparameters of the model-misspecification random field are determined automatically by maximising the marginal likelihood of the data observed up to the current time.
- The same sequential procedure applies without modification to one- and two-dimensional elastodynamic problems supplied with synthetic sensor data.
- Uncertainties from random forcing, parameter variability and structural model error propagate forward in time and are reduced at each measurement instant.
Where Pith is reading between the lines
- Because the posterior remains Gaussian after each update, the method could be embedded inside existing explicit time integrators with only modest extra linear-algebra cost per step.
- The same perturbation device might be reused for other time-dependent problems whose stochastic PDEs produce mildly non-Gaussian priors, such as diffusion or wave propagation outside solid mechanics.
- In engineering practice the calibrated misspecification field supplies a running diagnostic of where the finite-element model systematically deviates from reality.
- The closed-form posterior covariance could serve as the basis for adaptive mesh refinement or sensor placement decisions that minimise future uncertainty.
Load-bearing premise
The perturbation expansion remains accurate enough to represent the non-Gaussian prior induced by the stochastic PDEs for forcing, material parameters and model discrepancy.
What would settle it
A low-dimensional elastodynamic test problem in which the exact posterior can be computed by direct sampling or quadrature; the perturbation-based mean and covariance would then be compared directly to the sampled statistics for increasing levels of stochastic forcing or material variability.
Figures
read the original abstract
We present an approach for synthesising observational data with elastodynamic finite element models by extending the statistical finite element method (statFEM) framework. The proposed formulation adopts a Bayesian filtering approach to account for uncertainties in the data, the finite element model, and the discrepancies between the model and the physical system. Observational data are assimilated while the state of the spatially discretised finite element problem is advanced in time using the stochastic variant of the explicit Newmark scheme. The prior probability density of the state is obtained by solving an incremental probabilistic forward problem, and the corresponding posterior density is obtained by conditioning on the data available at each time step. In the probabilistic forward problem, spatio-temporal Gaussian random fields representing the forcing, model misspecification, and material parameters are specified via their stochastic PDE formulation. The resulting non-Gaussian prior probability density is approximated using a perturbation approach, yielding a Gaussian posterior with closed-form mean and covariance. The hyperparameters of the random field representing model misspecification are calibrated by maximising the marginal likelihood of the data. The proposed approach is illustrated on one- and two-dimensional elastodynamic examples with synthetic data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper extends the statistical finite element method (statFEM) to sequential data assimilation in elastodynamics. It advances the spatially discretized state using a stochastic explicit Newmark scheme while assimilating observations via Bayesian filtering. Spatio-temporal Gaussian random fields for forcing, model misspecification, and material parameters are defined through stochastic PDEs; the resulting non-Gaussian prior is approximated by perturbation to obtain a closed-form Gaussian posterior at each step. Hyperparameters of the misspecification field are calibrated by maximizing the marginal likelihood, and the method is demonstrated on 1D and 2D synthetic elastodynamic examples.
Significance. If the perturbation approximation remains accurate under time stepping, the approach would offer an efficient, closed-form route to uncertainty-aware data synthesis in finite-element solid dynamics, extending statFEM to sequential settings with explicit treatment of model discrepancy and parameter uncertainty. The combination of stochastic PDE priors, Newmark propagation, and marginal-likelihood calibration is technically novel and could support engineering applications where full Monte-Carlo propagation is prohibitive.
major comments (2)
- [§3 (probabilistic forward problem) and §5 (numerical examples)] The central claim that the perturbation expansion produces a sufficiently accurate Gaussian prior (and thus reliable closed-form posterior) at each time step is load-bearing, yet the manuscript provides no quantitative validation of the approximation error. In particular, material parameters enter the discrete stiffness and mass matrices nonlinearly, so even a first-order perturbation truncates higher-order moments that the explicit Newmark integrator can amplify over many steps. No a-posteriori error indicator, convergence study with respect to perturbation order, or comparison against Monte-Carlo propagation of the full non-Gaussian process appears in the numerical examples.
- [§5] The synthetic-data illustrations in §5 report only qualitative agreement between assimilated and reference solutions. Quantitative error norms (e.g., L2 or energy-norm discrepancies versus the true synthetic trajectory), calibration diagnostics for the marginal-likelihood hyperparameters, and baseline comparisons (deterministic Newmark, ensemble Kalman filter, or full Monte-Carlo statFEM) are absent, leaving the practical accuracy of the Gaussian posterior unquantified.
minor comments (2)
- [§3] Notation for the stochastic Newmark update and the perturbation expansion should be collected in a single table or appendix for readability.
- [§3.2] The manuscript should state explicitly whether the perturbation is first-order or higher and whether the resulting covariance is exact or approximate under the linearised map.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments that highlight important aspects of the perturbation approximation and the numerical validation. We respond to each major comment below and will revise the manuscript accordingly to strengthen the quantitative support for the claims.
read point-by-point responses
-
Referee: [§3 (probabilistic forward problem) and §5 (numerical examples)] The central claim that the perturbation expansion produces a sufficiently accurate Gaussian prior (and thus reliable closed-form posterior) at each time step is load-bearing, yet the manuscript provides no quantitative validation of the approximation error. In particular, material parameters enter the discrete stiffness and mass matrices nonlinearly, so even a first-order perturbation truncates higher-order moments that the explicit Newmark integrator can amplify over many steps. No a-posteriori error indicator, convergence study with respect to perturbation order, or comparison against Monte-Carlo propagation of the full non-Gaussian process appears in the numerical examples.
Authors: We agree that the accuracy of the first-order perturbation approximation for the non-Gaussian prior is a central assumption, particularly because material parameters enter the discrete operators nonlinearly and the explicit Newmark scheme can propagate higher-order effects over multiple steps. The manuscript justifies the approximation theoretically for small uncertainties and illustrates the overall procedure on synthetic data, but does not provide direct quantitative checks against the full non-Gaussian process. In the revised manuscript we will add, in §5, a comparison of the perturbation-derived mean and covariance against Monte-Carlo estimates obtained by sampling the stochastic forward problem for the 1D example; we will also include a simple a-posteriori indicator based on the size of the neglected second-order terms and examine its growth over time steps. revision: yes
-
Referee: [§5] The synthetic-data illustrations in §5 report only qualitative agreement between assimilated and reference solutions. Quantitative error norms (e.g., L2 or energy-norm discrepancies versus the true synthetic trajectory), calibration diagnostics for the marginal-likelihood hyperparameters, and baseline comparisons (deterministic Newmark, ensemble Kalman filter, or full Monte-Carlo statFEM) are absent, leaving the practical accuracy of the Gaussian posterior unquantified.
Authors: We concur that the current numerical examples emphasize visual agreement and would be strengthened by quantitative metrics. The manuscript demonstrates the assimilation procedure on 1D and 2D synthetic trajectories but reports only qualitative results. In the revision we will augment §5 with L2 and energy-norm error tables comparing the posterior mean and covariance to the known synthetic reference at selected time instants, together with the evolution of the marginal-likelihood-optimized hyperparameters. We will also add brief comparisons against a deterministic Newmark integrator and an ensemble Kalman filter on the same examples to quantify the benefit of the proposed closed-form Gaussian posterior. revision: yes
Circularity Check
No significant circularity; derivation is self-contained in standard Bayesian and FE techniques
full rationale
The paper advances the elastodynamic state via the stochastic explicit Newmark scheme, obtains the prior by solving an incremental probabilistic forward problem with Gaussian random fields for forcing/misspecification/parameters, approximates the resulting non-Gaussian density by perturbation to yield a closed-form Gaussian posterior, and calibrates misspecification hyperparameters by maximising the marginal likelihood. None of these steps reduces by construction to a fitted quantity renamed as a prediction, nor relies on a self-citation chain whose cited result is itself unverified; the perturbation and filtering steps are independent approximations grounded in external Bayesian and numerical methods rather than tautological redefinitions of the inputs.
Axiom & Free-Parameter Ledger
free parameters (1)
- hyperparameters of the model-misspecification random field
axioms (2)
- domain assumption Spatio-temporal Gaussian random fields for forcing, model misspecification, and material parameters can be specified via their stochastic PDE formulation.
- ad hoc to paper The perturbation approach yields a sufficiently accurate Gaussian approximation to the non-Gaussian prior density.
Reference graph
Works this paper leans on
-
[1]
K. Willcox, B. Segundo, The role of computational science in digital twins, Nature Computational Science 4 (2024) 147–149
work page 2024
-
[2]
D. Wagg, K. Worden, R. Barthorpe, P. Gardner, Digital twins: state-of-the-art and future directions for modeling and simulation in engineering dynamics applications, ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems, Part B: Mechanical Engineering 6 (2020) 030901
work page 2020
-
[3]
L. Sun, Z. Shang, Y . Xia, S. Bhowmick, S. Nagarajaiah, Review of bridge structural health monitoring aided by big data and artificial intelligence: from condition assessment to damage detection, Journal of Structural Engineering 146 (2020) 04020073:1–04020073:22
work page 2020
-
[4]
M. Girolami, E. Febrianto, G. Yin, F. Cirak, The statistical finite element method (statFEM) for coherent synthesis of observation data and model predictions, Computer Methods in Applied Mechanics and Engineering 375 (2021) 113533:1–113533:32
work page 2021
- [5]
-
[6]
E. Febrianto, L. Butler, M. Girolami, F. Cirak, Digital twinning of self-sensing structures using the statistical finite element method, Data- Centric Engineering 3 (2022) e31
work page 2022
-
[7]
¨O. D. Akyildiz, C. Duffin, S. Sabanis, M. Girolami, Statistical finite elements via Langevin dynamics, SIAM/ASA Journal on Uncertainty Quantification 10 (2022) 1560–1585
work page 2022
-
[8]
K.-J. Koh, F. Cirak, Stochastic PDE representation of random fields for large-scale Gaussian process regression and statistical finite element analysis, Computer Methods in Applied Mechanics and Engineering 417 (2023) 116358
work page 2023
-
[9]
A. Glyn-Davies, C. Duffin, I. Kazlauskaite, M. Girolami,¨O. D. Akyildiz, Statistical finite elements via interacting particle Langevin dynamics, SIAM/ASA Journal on Uncertainty Quantification 13 (2025) 1200–1227
work page 2025
-
[10]
V . Narouie, H. Wessels, F. Cirak, U. R ¨omer, Mechanical state estimation with a polynomial-chaos-based statistical finite element method, Computer Methods in Applied Mechanics and Engineering 441 (2025) 117970
work page 2025
-
[11]
T. Karvonen, F. Cirak, M. Girolami, Error analysis for a statistical finite element method, Journal of Multivariate Analysis 210 (2025) 105468
work page 2025
-
[12]
V . K. Narouie, J.-H. Urrea-Quintero, F. Cirak, H. Wessels, Unsupervised constitutive model discovery from sparse and noisy data, Computer Methods in Applied Mechanics and Engineering 452 (2026) 118722
work page 2026
- [13]
-
[14]
C. M. Bishop, Pattern Recognition and Machine Learning, Springer, 2006
work page 2006
- [15]
-
[16]
J. L. Beck, Bayesian system identification based on probability logic, Structural Control and Health Monitoring 17 (2010) 825–847
work page 2010
-
[17]
S. S ¨arkk¨a, L. Svensson, Bayesian filtering and smoothing, vol. 17, Cambridge University Press, 2023
work page 2023
-
[18]
R. E. Kalman, A new approach to linear filtering and prediction problems, Transactions of the ASME: Journal of Basic Engineering 82 (1960) 35–45
work page 1960
-
[19]
A. H. Jazwinski, Stochastic processes and filtering theory, Courier Corporation, 2007
work page 2007
-
[20]
E. A. Wan, R. Van Der Merwe, The unscented Kalman filter for nonlinear estimation, in: Proceedings of the IEEE 2000 adaptive systems for signal processing, communications, and control symposium (Cat. No. 00EX373), 153–158, 2000. 26
work page 2000
-
[21]
G. Evensen, Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics, Journal of Geophysical Research: Oceans 99 (1994) 10143–10162
work page 1994
- [22]
-
[23]
A. Corigliano, S. Mariani, Parameter identification in explicit structural dynamics: performance of the extended Kalman filter, Computer methods in applied mechanics and engineering 193 (2004) 3807–3835
work page 2004
-
[24]
E. Lourens, E. Reynders, G. De Roeck, G. Degrande, G. Lombaert, An augmented Kalman filter for force identification in structural dynamics, Mechanical systems and signal processing 27 (2012) 446–460
work page 2012
-
[25]
S. E. Azam, E. Chatzi, C. Papadimitriou, A dual Kalman filter approach for state estimation via output-only acceleration measurements, Mechanical systems and signal processing 60 (2015) 866–886
work page 2015
-
[26]
K. Maes, A. W. Smyth, G. De Roeck, G. Lombaert, Joint input-state estimation in structural dynamics, Mechanical Systems and Signal Processing 70 (2016) 445–466
work page 2016
-
[27]
H. Ebrahimian, R. Astroza, J. P. Conte, C. Papadimitriou, Bayesian optimal estimation for output-only nonlinear system and damage identi- fication of civil structures, Structural Control and Health Monitoring 25 (4) (2018) e2128
work page 2018
- [28]
-
[29]
S. S ¨arkk¨a, A. Solin, Applied stochastic differential equations, vol. 10, Cambridge University Press, 2019
work page 2019
-
[30]
K. Burrage, I. Lenane, G. Lythe, Numerical methods for second-order stochastic differential equations, SIAM Journal on Scientific Computing 29 (1) (2007) 245–264
work page 2007
-
[31]
R. Mannella, Quasisymplectic integrators for stochastic differential equations, Physical Review E 69 (4)
-
[32]
M. C. Kennedy, A. O’Hagan, Bayesian calibration of computer models, Journal of the Royal Statistical Society: Series B (Statistical Method- ology) 63 (2001) 425–464
work page 2001
- [33]
-
[34]
A. N ´ovoa, A. Racca, L. Magri, Inferring unknown unknowns: Regularized bias-aware ensemble Kalman filter, Computer Methods in Applied Mechanics and Engineering 418 (2024) 116502
work page 2024
-
[35]
F. Lindgren, H. Rue, J. Lindstr ¨om, An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach, Journal of the Royal Statistical Society Series B: Statistical Methodology 73 (4) (2011) 423–498
work page 2011
- [36]
-
[37]
F. Lindgren, D. Bolin, H. Rue, The SPDE approach for Gaussian and non-Gaussian fields: 10 years and still running, Spatial Statistics (2022) 100599:1–100599:28
work page 2022
-
[38]
R. G. Ghanem, P. D. Spanos, Stochastic finite elements: a spectral approach, Courier Corporation, 2003
work page 2003
-
[39]
D. Xiu, Numerical methods for stochastic computations: a spectral method approach, Princeton University Press, 2010
work page 2010
- [40]
-
[41]
H. G. Matthies, A. Keese, Galerkin methods for linear and nonlinear elliptic stochastic partial differential equations, Computer Methods in Applied Mechanics and Engineering 194 (2005) 1295–1331
work page 2005
-
[42]
W. K. Liu, T. Belytschko, A. Mani, Random field finite elements, International Journal for Numerical methods in Engineering 23 (10) (1986) 1831–1845
work page 1986
-
[43]
S. S ¨arkk¨a, A. Solin, On continuous-discrete cubature Kalman filtering, IFAC Proceedings V olumes 45 (16) (2012) 1221–1226
work page 2012
-
[44]
D. J. MacKay, Comparison of approximate methods for handling hyperparameters, Neural computation 11 (5) (1999) 1035–1068
work page 1999
-
[45]
Sudret, Polynomial chaos expansions and stochastic finite element methods, Taylor & Francis, 2015
B. Sudret, Polynomial chaos expansions and stochastic finite element methods, Taylor & Francis, 2015
work page 2015
- [46]
-
[47]
C. V . Mai, M. D. Spiridonakos, E. N. Chatzi, B. Sudret, Surrogate modelling for stochastic dynamical systems by combining NARX models and polynomial chaos expansions, 2016
work page 2016
- [48]
-
[49]
P. Bernard, G. Fleury, Stochastic Newmark scheme, Probabilistic Engineering Mechanics 17 (1) (2002) 45–61
work page 2002
-
[50]
D. Roy, M. K. Dash, Explorations of a family of stochastic Newmark methods in engineering dynamics, Computer Methods in Applied Mechanics and Engineering 194 (45) (2005) 4758–4796
work page 2005
- [51]
-
[52]
T. Weiland, M. Pf ¨ortner, P. Hennig, Flexible and efficient probabilistic PDE solvers through Gaussian Markov random fields, 2025
work page 2025
- [53]
-
[54]
S. Chu, J. Guilleminot, C. Kelly, B. Abar, K. Gall, Stochastic modeling and identification of material parameters on structures produced by additive manufacturing, Computer Methods in Applied Mechanics and Engineering 387 (2021) 114166
work page 2021
-
[55]
I. Ben-Yelun, A. O. Yuksel, F. Cirak, Robust topology optimisation of lattice structures with spatially correlated uncertainties, Structural and Multidisciplinary Optimization 67 (2) (2024) 16
work page 2024
-
[56]
Torquato, Random heterogeneous materials: microstructure and macroscopic properties, vol
S. Torquato, Random heterogeneous materials: microstructure and macroscopic properties, vol. 16, Springer, 2002
work page 2002
- [57]
-
[58]
Y . S. Jetti, M. Ostoja-Starzewski, Correlation structures of statistically isotropic stiffness and compliance TRFs through upscaling, Computer 27 Methods in Applied Mechanics and Engineering 432 (2024) 117356
work page 2024
- [59]
- [60]
- [61]
-
[62]
M. D. Hoffman, D. D. Blei, C. Wang, J. Paisley, Stochastic variational inference, Journal of Machine Learning Research 14 (2013) 1303–1347
work page 2013
-
[63]
K. B. Petersen, M. S. Pedersen, The matrix cookbook, Technical University of Denmark 7 (15) (2008) 510
work page 2008
-
[64]
N. Grønbech-Jensen, O. Farago, A simple and effective Verlet-type algorithm for simulating Langevin dynamics, Molecular Physics 111 (8) (2013) 983–991
work page 2013
-
[65]
S. Melchionna, Design of quasisymplectic propagators for Langevin dynamics, The Journal of Chemical Physics 127 (4) (2007) 044108
work page 2007
-
[66]
N. Grønbech-Jensen, Complete set of stochastic Verlet-type thermostats for correct Langevin simulations, Molecular Physics 118 (8) (2020) e1662506, ISSN 0026-8976. 28
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.