The linear eddy mixing model is used to predict the evolution of a decaying scalar field in statistically steady homogeneous turbulent flow over a wide range of Reynolds and Schmidt numbers. Model results at low Reynolds number and order unity Schmidt number are shown to be in good overall agreement with direct numerical simulations. Results at higher Schmidt and Reynolds numbers reproduce conventional scaling properties of the scalar statistics. Predictions of Schmidt number and Reynolds number sensitivity of the evolution of the scalar concentration probability density function are presented and interpreted.