Using Wiener chaos expansion (WCE), we develop numerical algorithms for solving second-order linear parabolic stochastic partial differential equations (SPDEs). We propose a deterministic WCE-based algorithm for computing moments of the SPDE solutions without any use of the Monte Carlo technique. We also compare the proposed deterministic algorithm with two other numerical methods based on the Monte Carlo technique and demonstrate that the new method is more efficient for highly accurate solutions. Numerical tests verify that the scheme is of mean-square order) for diffusion and for diffusion-reaction SPDEs with constant or variable coefficients, where Δ is the time step, and N is the Wiener chaos order.
Introduction.In this paper we develop a new numerical method, based on nonlinear filtering ideas and spectral expansions, for advection-diffusion-reaction equations perturbed by random fluctuations, which form a broad class of second-order linear parabolic stochastic differential equations (SPDEs). The standard approach to constructing SPDE solvers starts with a space discretization of an SPDE, for which spectral methods (see, e.g., [4,10,14]), finite element methods (see, e.g., [1,8,32]) or spatial finite differences (see, e.g., [1,11,30,33]) can be used. The result of such a space discretization is a large system of ordinary stochastic differential equations (SDEs) which requires time discretization to complete a numerical algorithm. In [5,6] an SPDE is first discretized in time and then a finite element or finite difference method can be applied to this semidiscretization. Other numerical approaches include those making use of splitting techniques [2,17,12], quantization [9], or an approach based on the averaging-over-characteristic formula [26,27]. In [22,19] numerical algorithms based on the Wiener chaos expansion (WCE) were introduced for solving the nonlinear filtering problem for hidden Markov models. Since then the WCE-based numerical methods have been successfully developed in a number of directions (see, e.g., [13,31]).In computing moments of SPDE solutions, the existing approaches to solving SPDEs are usually complemented by the Monte Carlo technique. Consequently, in these approaches numerical approximations of SPDE moments have two errors: numerical integration error and Monte Carlo (statistical) error. To reach a high accuracy, we have to run a very large number of independent simulations of the SPDE to reduce the Monte Carlo error. Instead, here we exploit WCE numerical methods to construct a deterministic algorithm for computing moments of the SPDE solutions without any