Discrete and Continuous Dynamical Systems-Series B
An equation that arises in mathematical studies of the transport of pollutants in groundwater and of oil recovery processes is of the form: -del x . (kappa(x,.)del(x)u(x, omega)) = f(x), for x is an element of D, where kappa(x,.), the permeability tensor, is random and models the properties of the rocks, which are not know with certainty. Further, geostatistical models assume kappa(x, .) to be a log-normal random field. The use of Monte Carlo methods to approximate the expected value of u(x,.), higher moments, or other functionals of u(x,.), require solving similar system of equations many times as trajectories are considered, thus it becomes expensive and impractical. In this paper, we present and explain several advantages of using the White Noise probability space as a natural framework for this problem. Applying properly and timely the Wiener-Ito Chaos decomposition and an eigenspace decomposition, we obtain a symmetric positive definite linear system of equations whose solutions are the coefficients of a Galerkin-type approximation to the solution of the original equation. Moreover, this approach reduces the simulation of the approximation to u(x,omega) for a fixed omega, to the simulation of a finite number of independent normally distributed random variables.
First Page Number
Last Page Number
Roman, L. J., & Sarkis, M. (2006). Stochastic Galerkin Method for Elliptic Spdes: A White Noise Approach. Discrete and Continuous Dynamical Systems-Series B, 6(4), 941-955. Retrieved from https://digitalcommons.wpi.edu/mathematicalsciences-pubs/44