We present results from a systematic numerical study of decaying turbulence in a dilute polymer solution by using a shell-model version of the FENE-P equations. Our study leads to an appealing definition of drag reduction for the case of decaying turbulence. We exhibit several new results, such as the potential-energy spectrum of the polymer, hitherto unobserved features in the temporal evolution of the kinetic-energy spectrum, and characterize intermittency in such systems. We compare our results with the GOY shell model for fluid turbulence.PACS numbers: 47.27.Gs, 83.60.YzThe phenomenon of drag reduction by polymer additives[1], whereby dilute solutions of linear, flexible, high-molecular-weight polymers exhibit frictional resistance to flow much lower than that of the pure solvent, has almost exclusively been studied within the context of statistically steady turbulent flows since the pioneering work of Toms [2]. By contrast, there is an extreme scarcity of results concerning the effects of polymer additives on decaying turbulence [3]. Experimental studies of decaying, homogeneous turbulence behind a grid indicate, for such dilute polymer solutions, a turbulent energy spectrum similar to that found without polymers [4,5]. However, flow visualization via die-injection tracers[5] and particle image velocimetry [6] show an inhibition of small-scale structures in the presence of polymer additives. To the best of our knowledge decaying turbulence in such polymer solutions has not been studied numerically. We initiate such a study here by using a shell model that is well suited to examining the effects of polymer additives in turbulent flows that are homogeneous and in which bounding walls have no direct role. We obtain several interesting results including a natural definition of the percentage dragreduction DR, which has been lacking for the case of decaying turbulence. We show that the dependence of DR on the polymer concentration c is in qualitative accord with experiments[1] as is the suppression of small-scale structures which we quantify by obtaining the filteredwavenumber-dependence of the flatness of the velocity field. We will use a shell-model version of the FENE-P (Finitely Extensible Nonlinear Elastic -Peterlin) [7,8] model for dilute polymer solutions that has often been used for studying viscoelastic effects since it contains the basic characteristics of molecular stretching, orientation and finite extensibility seen in polymer molecules. A direct numerical simulation of the FENE-P equations is computationally prohibitive. This motivates the use of a shell model that captures the essential features of the FENE-P equations. Recent studies [9] have exploited a formal analogy[10] of the FENE-P equations with those of magnetohydrodynamics (MHD) to construct such a shell model. We investigate decaying turbulence in a dilute polymer solution by developing a similar shell model for the FENE-P equations. The unforced FENE-P equations [7,8] arewhere p is the pressure, ν s the kinematic viscosity of the solven...