A Bayesian numerical homogenization method for elliptic multiscale inverse problems

A new strategy based on numerical homogenization and Bayesian techniques for solving multiscale inverse problems is introduced. We consider a class of elliptic problems which vary at a microscopic scale, and we aim at recovering the highly oscillatory tensor from measurements of the fine scale solution at the boundary, using a coarse model based on numerical homogenization and model order reduction. We provide a rigorous Bayesian formulation of the problem, taking into account different possibilities for the choice of the prior measure. We prove well-posedness of the effective posterior measure and, by means of G-convergence, we establish a link between the effective posterior and the fine scale model. Several numerical experiments illustrate the efficiency of the proposed scheme and confirm the theoretical findings

STOCHASTIC MATCHED FILTERS FOR SIGNAL DETECTION APPLICATIONS

The stochastic matched filter (SMF) is a variation of the matched filter that can detect stochastic signals in noisy environments. Some earlier studies suggest that the SMF can be extended to the detection of frequency time-variant (nonstationary) signals, namely wideband modulated sonar in shallow water. This thesis considers the SMF algorithm first proposed by J.-F. Cavasillas in signal detection and estimation scenarios, and investigates its application to narrowband and chirp signals imbedded in white noise. In medium to high signal to noise ratio (SNR) values, results indicate that the SMF is a viable technique for signal detection and estimation, and could be employed in passive, real-time signal detection and estimation scenarios.Lieutenant Commander, United States NavyApproved for public release. Distribution is unlimited

Description Length Based Signal Detection in singular Spectrum Analysis

This paper provides an information theoretic analysis of the signal-noise separation problem in Singular Spectrum Analysis. We present a signal-plus-noise model based on the Karhunen-Loève expansion and use this model to motivate the construction of a minimum description length criterion that can be employed to select both the window length and the signal. We show that under very general regularity conditions the criterion will identify the true signal dimension with probability one as the sample size increases, and will choose the smallest window length consistent with the Whitney embedding theorem. Empirical results obtained using simulated and real world data sets indicate that the asymptotic theory is reflected in observed behaviour, even in relatively small samples.Karhunen-Loève expansion, minimum description length, signal-plus-noise model, Singular Spectrum Analysis, embedding

A spectral surrogate model for stochastic simulators computed from trajectory samples

Stochastic simulators are non-deterministic computer models which provide a different response each time they are run, even when the input parameters are held at fixed values. They arise when additional sources of uncertainty are affecting the computer model, which are not explicitly modeled as input parameters. The uncertainty analysis of stochastic simulators requires their repeated evaluation for different values of the input variables, as well as for different realizations of the underlying latent stochasticity. The computational cost of such analyses can be considerable, which motivates the construction of surrogate models that can approximate the original model and its stochastic response, but can be evaluated at much lower cost. We propose a surrogate model for stochastic simulators based on spectral expansions. Considering a certain class of stochastic simulators that can be repeatedly evaluated for the same underlying random event, we view the simulator as a random field indexed by the input parameter space. For a fixed realization of the latent stochasticity, the response of the simulator is a deterministic function, called trajectory. Based on samples from several such trajectories, we approximate the latter by sparse polynomial chaos expansion and compute analytically an extended Karhunen-Lo\`eve expansion (KLE) to reduce its dimensionality. The uncorrelated but dependent random variables of the KLE are modeled by advanced statistical techniques such as parametric inference, vine copula modeling, and kernel density estimation. The resulting surrogate model approximates the marginals and the covariance function, and allows to obtain new realizations at low computational cost. We observe that in our numerical examples, the first mode of the KLE is by far the most important, and investigate this phenomenon and its implications

Ensemble Kalman filter for multiscale inverse problems

We present a novel algorithm based on the ensemble Kalman filter to solve inverse problems involving multiscale elliptic partial differential equations. Our method is based on numerical homogenization and finite element discretization and allows to recover a highly oscillatory tensor from measurements of the multiscale solution in a computationally inexpensive manner. The properties of the approximate solution are analysed with respect to the multiscale and discretization parameters, and a convergence result is shown to hold. A reinterpretation of the solution from a Bayesian perspective is provided, and convergence of the approximate conditional posterior distribution is proved with respect to the Wasserstein distance. A numerical experiment validates our methodology, with a particular emphasis on modelling error and computational cost

Numerical treatment of imprecise random fields in non-linear solid mechanics

The quantification and propagation of mixed uncertain material parameters in the context of solid mechanical finite element simulations is studied. While aleatory uncertainties appear in terms of spatial varying parameters, i.e. random fields, the epistemic character is induced by a lack of knowledge regarding the correlation length, which is therefore described by interval values. The concept and description of the resulting imprecise random fields is introduced in detail. The challenges occurring from interval valued correlation lengths are clarified. These include mainly the stochastic dimension, which can become very high under some circumstances, as well as the comparability of different correlation length scenarios with regard to the underlying truncation error of the applied Karhunen-Loève expansion. Additionally, the computation time can increase drastically, if the straightforward and robust double loop approach is applied. Sparse stochastic collocation method and sparse polynomial chaos expansion are studied to reduce the number of required sample evaluations, i.e. the computational cost. To keep the stochastic dimension as low as possible, the random fields are described by Karhunen-Loève expansion, using a modified exponential correlation kernel, which is advantageous in terms of a fast convergence while providing an analytic solution. Still, for small correlation lengths, the investigated approaches are limited by the curse of dimensionality. Furthermore, they turn out to be not suited for non-linear material models. As a straightforward alternative, a decoupled interpolation approach is proposed, offering a practical engineering estimate. For this purpose, the uncertain quantities only need to be propagated as a random variable and deterministically in terms of the mean values. From these results, the so-called absolutely no idea probability box (ani-p-box) can be obtained, bounding the results of the interval valued correlation length being between zero and infinity. The idea is, to interpolate the result of any arbitrary correlation length within this ani-p-box, exploiting prior knowledge about the statistical behaviour of the input random field corresponding to the correlation length. The new approach is studied for one- and two-dimensional random fields. Furthermore, linear and non-linear finite element models are used in terms of linear-elastic or elasto-plastic material laws, the latter including linear hardening. It appears that the approach only works satisfyingly for sufficiently smooth responses but an improvement by considering also higher order statistics is motivated for future research.DFG/SPP 1886/NA330/12-1/E

Modeling Spatially Varying Uncertainty in Composite Structures Using Lamination Parameters

An approach is presented for modeling spatially varying uncertainty in the ply orientations of composite structures. Lamination parameters are used with the aim of reducing the required number of random variables. Karhunen–Loève expansion is employed to decompose the uncertainty in each ply into a sum of random variables and spatially dependent functions. An intrusive polynomial chaos expansion is proposed to approximate the lamination parameters while preserving the separation of the random and spatial dependency. Closed-form expressions are derived for the expansion coefficients in two case studies; an initial example in which uncertainty is modeled using random variables, and a second random field example. The approach is compared against Monte Carlo simulation results for a variety of layups as well as closed-form expressions for the mean and covariance. By summing the polynomial chaos basis functions through the laminate thickness, the separation of the random and spatial dependency may be preserved at a laminate level and the number of random variables reduced for some minimum number of plies. The number of variables increases nonlinearly with the number of Karhunen–Loève expansion terms, and as such, the approach is only beneficial in low-order expansions using relatively few Karhunen–Loève expansion terms

Probabilistic micromechanical spatial variability quantification in laminated composites

SN and SS are grateful for the support provided through the Lloyd’s Register Foundation Centre. The Foundation helps to protect life and property by supporting engineering-related education, public engagement and the application of research.Peer reviewedPostprin

Model order reduction techniques in microelectromechanics

Ph.DDOCTOR OF PHILOSOPH