Analytical approximation and numerical studies of one-dimensional elliptic equation with random coefficients

Abstract

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, ω) with known covariance function. We derive the relation between the standard deviation of the solution T(x, ω) and the correlation length, η of K(x, ω). We observe that, the standard deviation, σT of the solution, T(x, ω), initially increases with the correlation length η up to a maximum value, σT, max at ηmax∼x(1−x)/3 and decreases beyond ηmax . We observe a scaling law between σT and η, that is, σT ∝ η1/2 for η → 0 and σT∝η−1/2 for η → ∞. We show that, for a small value of coefficient of variation (ɛK=σK/μK) of the random coefficient, the solution T(x, ω) can be approximated with a Gaussian random field regardless of the underlying probability distribution of K(x, ω). This approximation is valid for large value of εK, if the correlation length, η of input random field K(x, ω) 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, ω) 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.