Spectrum of passive scalars of high molecular diffusivity in turbulent mixing

We consider the mixing of passive scalars transported in turbulent flow, with a molecular diffusivity that is large compared to the kinematic viscosity of the fluid. This particular case of mixing has not received much attention in experiment or simulation even though the first putative theory, due to Batchelor, Howells & Townsend (J. Fluid Mech., vol. 5, 1959, pp. 134-139), is now more than 50 years old. We study the problem using direct numerical simulation of decaying scalar fields in steadily sustained homogeneous turbulence as the Schmidt number (the ratio of the kinematic viscosity of the fluid to the molecular diffusivity of the scalar) is allowed to vary from 1/8 to 1/2048 for two values of the microscale Reynolds number, R_λ\ ≈ 140 and ≈ 240. The simulations show that the passive scalar spectrum assumes a slope of-17/3 in a range of scales, as predicted by the theory, when the Schmidt number is small and the Reynolds number is simultaneously large. The observed agreement between theory and simulation in the prefactor in the spectrum is not perfect. We assess the reasons for this discrepancy by a careful examination of the scalar evolution equation in the light of the assumptions of the theory, and conclude that the finite range of scales resolved in simulations is the main reason. Numerical issues specific to the regime of very low Schmidt numbers are also addressed briefly.