摘要

In this work, we study a one-dimensional elliptic equation with a random coefficient and derive an explicit analytical approximation. We model the random coefficient with a spatially varying random field, K(x, co) with known covariance function. We derive the relation between the standard deviation of the solution T(x, co) and the correlation length, 77 of IC(x, co). We observe that, the standard deviation, UT of the solution, T(x, co), initially increases with the correlation length eta up to a maximum value, QT, max at eta(max) root x(1 x)/3 and decreases beyond nmax. We observe a scaling law between criand 77, that is, QT alpha eta(1/2) for eta -> 0 and sigma(T) proportional to n(-1/2) for eta -> infinity. We show that, for a small value of coefficient of variation (sK = 1.4K) of the random coefficient, the solution T(x, co) can be approximated with a Gaussian random field regardless of the underlying probability distribution of K(x, co). This approximation is valid for large value of 6K, if the correlation length, s of input random field (epsilon(K) = sigma(K)/mu(k)) is small. We compare the analytical results with numerical ones obtained from Monte-Carlo method and polynomial chaos based stochastic collocation method. Under aforementioned conditions, we observe a good agreement between the numerical simulations and the analytical results. For a given random coefficient K(x, co) with known mean and variance we can quickly estimate the variance of the solution at any location for a given correlation length. If the correlation length is not available which is the case in most practical situations, we can still use this analytical solution to estimate the maximum variance of the solution at any location.

  • 出版日期2016-5