Stochastic collocation approach with adaptive mesh refinement for parametric uncertainty analysis

Abstract

Presence of a high-dimensional stochastic parameter space with discontinuities poses major computational challenges in analyzing and quantifying the effects of the uncertainties in a physical system. In this paper, we propose a stochastic collocation method with adaptive mesh refinement (SCAMR) to deal with high dimensional stochastic systems with discontinuities. Specifically, the proposed approach uses generalized polynomial chaos (gPC) expansion with Legendre polynomial basis and solves for the gPC coefficients using the least squares method. It also implements an adaptive mesh (element) refinement strategy which checks for abrupt variations in the output based on the second order gPC approximation error to track discontinuities or non-smoothness. In addition, the proposed method involves a criterion for checking possible dimensionality reduction and consequently, the decomposition of the full-dimensional problem to a number of lower-dimensional subproblems. Specifically, this criterion checks all the existing interactions between input dimensions of a specific problem based on the high-dimensional model representation (HDMR) method, and therefore automatically provides the subproblems which only involve interacting dimensions. The efficiency of the approach is demonstrated using both smooth and non-smooth function examples with input dimensions up to 300, and the 1 Corresponding author.approach is compared against other existing algorithms.

An adaptive high-order piecewise polynomial based sparse grid collocation method with applications

“…The error has been calculated using randomly generated 100000 sample points. We will use the following five functions which have been considered in previous work [18,19,3].…”

“…Theoretical justification will be provided, and applications in stochastic differential equations are considered. We note that many recent work in UQ has considered more efficient collocation schemes in higher dimensions and for functions with singularities [27,9,8,20,3]. Those techniques are not explored in this work.…”

An efficient optimization based microstructure reconstruction approach with multiple loss functions

“…The methods [14,15,16,17] are mostly analytical in nature and hence much faster than the optimization-based reconstruction approaches but their application is restricted to isotropic binary materials. Another class of microstructure reconstruction methods involves using deep learning [18,19], a machine learning approach that can be used amongst other tasks for surrogate modeling [20,21,22,23,24,25,26] and thus has been implemented successfully for a wide range of classification and regression tasks. Deep learning approaches, in particular convolutional neural networks, are particularly well-suited to handle image data, and have thus received attention lately from the materials research community to process microstructures image data for a variety of tasks [27,28].…”

Reinforcement Learning based Sequential Batch-sampling for Bayesian Optimal Experimental Design

“…In recent years, a plethora of research has focused on Bayesian design of experiments [8,9,10,11,12], both one-shot designs [13,14], myopic sequential design of experiments (SDOE) [15,16]. and adaptive refinement methods for surrogate modeling [17,18,19]. Batch designs for querying the expensive information source, using the so-called batch Bayesian optimization have been proposed in recent years [20,21,22,23,24], with promising results.…”