Analytical solutions of the reactant concentration inside porous spherical catalytic particles were obtained from unsteady reaction-diffusion equation by applying eigenfunction expansion method. Various surface concentrations as exponentially decaying or oscillating function were considered as boundary conditions to solve the unsteady partial differential equation as a function of radial distance and time. Dirac delta function was also used for the instantaneous injection of the reactant as the surface boundary condition to calculate average reactant concentration inside the particles as a function of time by Laplace transform. Besides spherical morphology, other geometries of particles, such as cylinder or slab, were considered to obtain the solution of the reaction-diffusion equation, and the results were compared with the solution in spherical coordinate. The concentration inside the particles based on calculation was compared with the bulk concentration of the reactant molecules measured by photocatalytic decomposition as a function of time.