Elliptic equations of higher stochastic order

¹ Department of Mathematics, USC, Los Angeles, CA 90089, USA. lototsky@math.usc.edu; http://www-rcf.usc.edu/~lototsky
² Division of Applied Mathematics, Brown University, Providence, RI 02912, USA. rozovsky@dam.brown.edu
³ Department of Mathematics, Center for Computation and Technology, Louisiana State University, Baton Rouge, LA 70803, USA. xlwan@math.lsu.edu

Received: 3 July 2009

Abstract

This paper discusses analytical and numerical issues related to elliptic equations with random coefficients which are generally nonlinear functions of white noise. Singularity issues are avoided by using the Itô-Skorohod calculus to interpret the interactions between the coefficients and the solution. The solution is constructed by means of the Wiener Chaos (Cameron-Martin) expansions. The existence and uniqueness of the solutions are established under rather weak assumptions, the main of which requires only that the expectation of the highest order (differential) operator is a non-degenerate elliptic operator. The deterministic coefficients of the Wiener Chaos expansion of the solution solve a lower-triangular system of linear elliptic equations (the propagator). This structure of the propagator insures linear complexity of the related numerical algorithms. Using the lower triangular structure and linearity of the propagator, the rate of convergence is derived for a spectral/hp finite element approximation. The results of related numerical experiments are presented.