Research on turbulent mixing processes is of great interest to those working on turbulent-reactive flows. In this paper, a detailed study has been performed for the evolution of scalar fields of different initial integral scales in decaying, homogeneous and isotropic turbulence using DNS technique. Passive scalar mixing in a cubical decaying, homogeneous, isotropic turbulence field is considered. The three-dimensional incompressible Navier-Stokes equations together with scalar equation are solved using Fractional Step Method. The convective and diffusive terms in governing equations are discretised by Compact Finite Difference Scheme. The 32 × 32 × 32 uniform staggered grids are used. The present simulation is performed at Taylor Reynolds number of 28.83. In this paper, the evolution of scalar RMS and scalar dissipation rate for different integral length scales has been presented. The initial velocity vector and Probability Density Function (PDF) of scalar at different eddy turn over time have also been presented.