Nonparametric Stochastic Subspaces via the Bootstrap for Characterizing Model Error
Akash Yadav, Ruda Zhang
TL;DR
This work tackles forward uncertainty quantification under epistemic and aleatory uncertainty by focusing on model-form error. It introduces SS-Bootstrap, a bootstrap-based, nonparametric stochastic subspace method within the stochastic reduced-order modeling framework, to characterize model error directly from empirical snapshots while preserving linear constraints. The approach requires only a single hyperparameter, beta, and yields sharper, yet reliable uncertainty estimates compared to Gaussian latent-variable alternatives like SS-PPCA, demonstrated on parametric linear static and linear dynamics problems with favorable computation times. The method is robust to non-Gaussian data and naturally extensible to other sampling-based model-reduction techniques and nonlinear problems, with future work aimed at model-error correction and broader MOR applications.
Abstract
Reliable forward uncertainty quantification in engineering requires methods that account for aleatory and epistemic uncertainties. In many applications, epistemic effects arising from uncertain parameters and model form dominate prediction error and strongly influence engineering decisions. Because distinguishing and representing each source separately is often infeasible, their combined effect is typically analyzed using a unified model-error framework. Model error directly affects model credibility and predictive reliability; yet its characterization remains challenging. To address this need, we introduce a bootstrap-based stochastic subspace model for characterizing model error in the stochastic reduced-order modeling framework. Given a snapshot matrix of state vectors, the method leverages the empirical data distribution to induce a sampling distribution over principal subspaces for reduced order modeling. The resulting stochastic model enables improved characterization of model error in computational mechanics compared with existing approaches. The method offers several advantages: (1) it is assumption-free and leverages the empirical data distribution; (2) it enforces linear constraints (such as boundary conditions) by construction; (3) it requires only one hyperparameter, significantly simplifying the training process; and (4) its algorithm is straightforward to implement. We evaluate the method's performance against existing approaches using numerical examples in computational mechanics and structural dynamics.
