Direct numerical simulations of passive scalars, with Prandtl numbers Pr=3, 5, and 7, advected by turbulence at three low Reynolds numbers were performed. The energy spectra are self-similar under the Kolmogorov scaling and exhibit behaviour consistent with many other investigations: a short inertial range for the highest Reynolds number and the universal exponential form of the spectrum for all Reynolds numbers in the dissipation range. In all cases the passive scalar spectra collapse to a single self-similar curve under the Batchelor scaling and exhibit the k−1 range followed by an exponential fall-off. We attribute the applicability of the Batchelor scaling to our low-Reynolds-number flows to the universality of the energy dissipation spectra. The Batchelor range is observed for wavenumbers in general agreement with experimental observations but smaller than predicted by the classical estimates. The discrepancy is caused by the fact that the velocity scales responsible for the generation of the Batchelor range are in the vicinity of the wavenumber of the maximum energy dissipation, which is one order of magnitude less than the Kolmogorov wavenumber used in the classical theory. Two different functional forms of passive scalar spectra proposed by Batchelor and Kraichnan were fitted to the simulation results and it was found that the Kraichnan model agrees very well with the data while the Batchelor formula displays systematic deviations from the data. Implications of these differences for the experimental procedures to measure the energy and passive scalar dissipation rates in oceanographic flows are discussed.