Polytopic symmetric interior-penalty approximation of elliptic SPDEs driven by spatial white noise

Abstract

We analyse symmetric interior-penalty discontinuous Galerkin approximations of elliptic stochastic partial differential equations driven by spatial Gaussian white noise on quasi-uniform polytopic meshes. The low regularity of the exact random field precludes the usual SIP consistency argument. We instead restrict the same continuum white-noise functional canonically to the discontinuous polynomial space and separate the error into stochastic forcing-restriction and deterministic discretisation components. For every fixed polynomial degree $k\geq1$, we prove the sharp two-sided estimate \[ \|Y-Y_h\|_{L^2(\Xi;L^2(\Omega))} \simeq h^{2-d/2}, \] giving rates $h$ in two dimensions and $h^{1/2}$ in three dimensions. A finite-dimensional lower bound establishes optimality. The analysis allows polytopic elements with arbitrarily many and arbitrarily small faces, including admissible non-convex elements. In two dimensions we also treat polygons with a single reentrant corner. Although the deterministic approximation exhibits an angle-dependent loss of regularity, the singular resolvent contribution is finite rank, and the sharp first-order stochastic rate is retained. The discontinuous structure yields an exact element-local sampler for the restricted white noise. Numerical experiments on irregular polygonal, non-convex polytopic and reentrant domains support the predicted convergence behaviour.