This paper introduces the series expansion for homogeneous, isotropic and mean square continuous random fields in the Euclidean space, which involves the Bessel function and the ultraspherical polynomial, but differs from the spectral representation in terms of the ordinary spherical harmonics that has more terms at each level.The series representation provides a simple and efficient approach for simulation of isotropic (non-Gaussian) random fields.