Weak convergence of finite element approximations of linear stochastic evolution equations with additive noise II. Fully discrete schemes
We present an abstract framework for analyzing the weak error of fully discrete approximation schemes for linear evolution equations driven by additive Gaussian noise. First, an abstract representation formula is derived for sufficiently smooth test functions. The formula is then applied to the wave equation, where the spatial approximation is done via the standard continuous finite element method and the time discretization via an I-stable rational approximation to the exponential function. It is found that the rate of weak convergence is twice that of strong convergence. Furthermore, in contrast to the parabolic case, higher order schemes in time, such as the Crank-Nicolson scheme, are worthwhile to use if the solution is not very regular. Finally we apply the theory to parabolic equations and detail a weak error estimate for the linearized Cahn-Hilliard-Cook equation as well as comment on the stochastic heat equation.
weak convergence finite element stochastic additive noise wave equation time discretization rational approximation parabolic equation error estimate cahn-hilliard-cook equation heat equation hyperbolic equation wiener process 65m60 60h15 60h35 65c30