The present book develops the mathematical and numerical analysis of linear, elliptic and parabolic partial differential equations (PDEs) with coefficients whose logarithms are modelled as Gaussian random fields (GRFs), in polygonal and polyhedral physical domains.