The Rayleigh-Taylor instability plays an important role in the dynamics of several astronomical objects, in particular, in supernovae (SN) evolution. In this paper we develop an analytical approach to study the stability analysis of spherical expansion of the SN ejecta by using a special transformation in the co-moving coordinate frame. We first study a non-stationary spherical expansion of a gas shell under the pressure of a central source. Then we analyze its stability with respect to a no radial, non spherically symmetric perturbation of the of the shell. We consider the case where the polytropic constant of the SN shell is $\gamma=5/3$ and we examine the evolution of a arbitrary shell perturbation. The dispersion relation is derived. The growth rate of the perturbation is found and its temporal and spatial evolution is discussed. The stability domain depends on the ejecta shell thickness, its acceleration, and the perturbation wavelength.