Effects of uncertainties in the domain on the solution of Neumann boundary value problems in two spatial dimensions