Abstract:The efficient discretization of a multidimensional random field with high definition and large geometric size remains a significant challenge. Compared with the simulation of one-dimensional or two-dimensional random fields, the generation of three-dimensional (3-D) random fields using the traditional Karhunen–Loève (K–L) expansion method tends to require a relatively long computational time and a larger amount of physical memory. In the present study, a decomposed K–L expansion scheme is proposed, which is applicable when a separable autocorrelation function (ACF) is used. The separability in the ACF is a precondition for the implementation of the proposed method. The proposed method decomposes the discretization of a 3-D random field into that of separate one-dimensional random fields and computes the eigenpairs in each dimension respectively. The accuracy and efficiency of the proposed method are demonstrated, and the numerical solutions of the eigenpairs used for the random field discretization are validated by comparing them with the theoretical solutions. Comparisons with the traditional K–L expansion method showed that the proposed scheme significantly reduced the computational time and memory requirements. This makes it potentially useful for the discretization of a multidimensional random field with a large geometric size and high definition, which can be used in stochastic finite element analysis. The proposed random field discretization method was applied to analyze a 3-D saturated clay slope under undrained conditions. The results of the stochastic finite element analysis suggest that the 3-D model is rather advantageous for estimating slope stability when the spatial variability of soil properties is considered. The 3-D model can capture more spatial details in soil distribution and its effects on failure mechanisms and safety factors from a theoretical point of view. Therefore, the proposed decomposed K–L expansion method for 3-D random field discretization has strong application potential in geotechnical problems.

An Interpolation Method for Karhunen–Loève Expansion of Random Field Discretization

A Decomposed Karhunen–Loève Expansion Scheme for the Discretization of Multidimensional Random Fields in Geotechnical Variability Analysis

A Chebyshev Polynomial-Based Galerkin Method for the Discretization of Spatially Varying Random Properties

Numerical Simulation of Random Fields with a High-Order Polynomial Based Ritz–Galerkin Approach

A Galerkin Isogeometric Method for Karhunen-Loeve Approximation of Random Fields

Simulation of Multi-Dimensional Random Fields by Karhunen–Loève Expansion

On the Numerical Approximation of the Karhunen-Loève Expansion for Random Fields with Random Discrete Data

Research on Application of Karhunen-Loeve Expansion to Simulating Anisotropic Random Field of Soil Property

A convergent adaptive finite element stochastic Galerkin method based on multilevel expansions of random fields

Karhunen-Loeve Series Based on Monte-Carlo Finite Element Method in Stochastic Analysis

Karhunen-Loéve Expansion for Random Earthquake Excitations

Local Error Estimates for Radial Basis Function Interpolation of Scattered Data

Stochastic Finite Element Methods with the Euclidean Degree for Partial Differential Equations with Random Inputs

Error Estimates for the Interpolating Moving Least-Squares Method in N -Dimensional Space

An alternating low-rank projection approach for partial differential equations with random inputs

Local Multilevel Scattered Data Interpolation

A Note on Richardson Extrapolation of Galerkin Methods for Eigenvalue Problems of Fredholm Integral Equations

Interpolating Element-Free Galerkin Method for the Regularized Long Wave Equation and Its Error Analysis

Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods

Adaptive Fourier decomposition-type sparse representations versus the Karhunen-Loeve expansion for decomposing stochastic processes

A Non Grid-Based Interpolation Scheme for the Eigenvalue Problem