STRONG AND WEAK ERROR ESTIMATES FOR ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS WITH RANDOM COEFFICIENTS

We consider the problem of numerically approximating the solution of an elliptic partial differential equation with random coefficients and homogeneous Dirichlet boundary conditions. We focus on the case of a lognormal coefficient and deal with the lack of uniform coercivity and uniform boundedness with respect to the randomness. This model is frequently used in hydrogeology. We approximate this coefficient by a finite dimensional noise using a truncated Karhunen-Loeve expansion. We give estimates of the corresponding error on the solution, both a strong error estimate and a weak error estimate, that is, an estimate of the error commited on the law of the solution. We obtain a weak rate of convergence which is twice the strong one. In addition, we give a complete error estimate for the stochastic collocation method in this case, where neither coercivity nor boundedness is stochastically uniform. To conclude, we apply these results of strong and weak convergence to two classical cases of covariance kernel choices, the case of an exponential covariance kernel on a box and the case of an analytic covariance kernel, yielding explicit weak and strong convergence rates.