Polynomial chaos representation of spatio-temporal random fields from experimental measurements

Two numerical techniques are proposed to construct a polynomial chaos (PC) representation of an arbitrary second-order random vector. In the first approach, a PC representation is constructed by matching a target joint probability density function (pdf) based on sequential conditioning (a sequence of conditional probability relations) in conjunction with the Rosenblatt transformation. In the second approach, the PC representation is obtained by having recourse to the Rosenblatt transformation and simultaneously matching a set of target marginal pdfs and target Spearman's rank correlation coefficient (SRCC) matrix. Both techniques are applied to model an experimental spatio-temporal data set, exhibiting strong non-stationary and non-Gaussian features. The data consists of a set of oceanographic temperature records obtained from a shallow-water acoustics transmission experiment [1]. The measurement data, observed over a finite denumerable subset of the indexing set of the random process, is treated as a collection of observed samples of a second-order random vector that can be treated as a finite-dimensional approximation of the original random field. A set of properly ordered conditional pdfs, that uniquely characterizes the target joint pdf, in the first approach and a set of target marginal pdfs and a target SRCC matrix, in the second approach, are estimated from available experimental data. Digital realizations sampled from the constructed PC representations based on both schemes capture the observed statistical characteristics of the experimental data with sufficient accuracy. The relative advantages and disadvantages of the two proposed techniques are also highlighted.

[1]  P. A. G. van der Geest,et al.  An algorithm to generate samples of multi-variate distributions with correlated marginals , 1998 .

[2]  R. A. Silverman,et al.  Special functions and their applications , 1966 .

[3]  D. B. Preston.  Spectral Analysis and Time Series , 1983 .

[4]  Shane G. Henderson,et al.  Behavior of the NORTA method for correlated random vector generation as the dimension increases , 2003, TOMC.

[5]  P. Embrechts,et al.  Risk Management: Correlation and Dependence in Risk Management: Properties and Pitfalls , 2002 .

[6]  Robert H. Headrick,et al.  An overview of the 1995 SWARM shallow-water internal wave acoustic scattering experiment , 1997 .

[7]  Roger Ghanem,et al.  Ingredients for a general purpose stochastic finite elements implementation , 1999 .

[8]  Roger G. Ghanem,et al.  On the construction and analysis of stochastic models: Characterization and propagation of the errors associated with limited data , 2006, J. Comput. Phys..

[9]  Ronald L. Wasserstein,et al.  Monte Carlo: Concepts, Algorithms, and Applications , 1997 .

[10]  M. Rosenblatt.  Remarks on a Multivariate Transformation , 1952 .

[11]  Thorsten Rheinländer.  Risk Management: Value at Risk and Beyond , 2003 .

[12]  Stephen P. Boyd,et al.  Semidefinite Programming , 1996, SIAM Rev..

[13]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[14]  W. T. Martin,et al.  The Orthogonal Development of Non-Linear Functionals in Series of Fourier-Hermite Functionals , 1947 .

[15]  G. Phillips.  Interpolation and Approximation by Polynomials , 2003 .

[16]  Kok-Kwang Phoon,et al.  Convergence study of the truncated Karhunen–Loeve expansion for simulation of stochastic processes , 2001 .

[17]  Wolfgang Hörmann,et al.  Automatic Nonuniform Random Variate Generation , 2011 .

[18]  N. Wiener.  The Homogeneous Chaos , 1938 .

[19]  R. Ghanem,et al.  Multi-resolution analysis of wiener-type uncertainty propagation schemes , 2004 .

[20]  E. Lehmann,et al.  Nonparametrics: Statistical Methods Based on Ranks , 1976 .

[21]  R. Nelsen.  An Introduction to Copulas , 1998 .

[22]  A. Izenman.  Recent Developments in Nonparametric Density Estimation , 1991 .

[23]  D. W. Scott,et al.  Multivariate Density Estimation, Theory, Practice and Visualization , 1992 .

[24]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[25]  Roger G. Ghanem,et al.  Asymptotic Sampling Distribution for Polynomial Chaos Representation of Data: A Maximum Entropy and Fisher information approach , 2006, Proceedings of the 45th IEEE Conference on Decision and Control.

[26]  Sonjoy Das,et al.  Model, identification & analysis of complex stochastic systems: Applications in stochastic partial differential equations and multiscale mechanics , 2008 .

[27]  Roger Ghanem,et al.  A Method for the Validation of Predictive Computations Using a Stochastic Approach , 2004 .

[28]  Christoph Schwab,et al.  Karhunen-Loève approximation of random fields by generalized fast multipole methods , 2006, J. Comput. Phys..

[29]  Christian Soize,et al.  Maximum likelihood estimation of stochastic chaos representations from experimental data , 2006 .

[30]  G. Karniadakis,et al.  An adaptive multi-element generalized polynomial chaos method for stochastic differential equations , 2005 .

[31]  G. Tolstov.  Fourier Series , 1962 .

[32]  Roger Ghanem,et al.  Analysis of Eigenvalues and Modal Interaction of Stochastic Systems , 2005 .

[33]  Roger G. Ghanem,et al.  Physical Systems with Random Uncertainties: Chaos Representations with Arbitrary Probability Measure , 2005, SIAM J. Sci. Comput..

[34]  Roger M. Cooke,et al.  Uncertainty Analysis with High Dimensional Dependence Modelling , 2006 .

[35]  H. Najm,et al.  Uncertainty quantification in reacting-flow simulations through non-intrusive spectral projection , 2003 .

[36]  R. Nelsen.  An Introduction to Copulas (Springer Series in Statistics) , 2006 .

[37]  R. Ghanem,et al.  Iterative solution of systems of linear equations arising in the context of stochastic finite elements , 2000 .

[38]  Roger Ghanem,et al.  Simulation of multi-dimensional non-gaussian non-stationary random fields , 2002 .

[39]  P. Billingsley,et al.  Probability and Measure , 1980 .

[40]  Lynch,et al.  Acoustic field variability induced by time evolving internal wave fields , 2000, The Journal of the Acoustical Society of America.

[41]  George Deodatis,et al.  Simulation of Highly Skewed Non-Gaussian Stochastic Processes , 2001 .

[42]  M. Loève.  Probability Theory II , 1978 .

[43]  Jeroen A. S. Witteveen,et al.  Modeling physical uncertainties in dynamic stall induced fluid-structure interaction of turbine blades using arbitrary polynomial chaos , 2007 .

[44]  Shane G. Henderson,et al.  Chessboard Distributions and Random Vectors with Specified Marginals and Covariance Matrix , 2002, Oper. Res..

[45]  A. Izenman.  Review Papers: Recent Developments in Nonparametric Density Estimation , 1991 .

[46]  H. O. Lancaster.  Some properties of the bivariate normal distribution considered in the form of a contingency table , 1957 .

[47]  Fawwaz T. Ulaby,et al.  A computationally efficient multivariate maximum-entropy density estimation (MEDE) technique , 2004, IEEE Transactions on Geoscience and Remote Sensing.