In this work, a diffusion equation model (DEM) is applied to a room acoustics case for in-depth sound field analysis. Background of the theory, the governing and boundary equations specifically applicable to this study are presented. A three-dimensional geometric model of a monumental worship space is composed. The DEM is solved over this model in a finite element framework to obtain sound energy densities. The sound field within the monument is numerically assessed; spatial sound energy distributions and flow vector analysis are conducted through the time-dependent DEM solutions.