DNS of passive scalars in turbulent pipe flow

We study the statistics of passive scalars at Pr = 1, for turbulent flow within a smooth straight pipe of circular cross section up to Re-tau approximate to 6000 using direct numerical simulation (DNS) of the Navier-Stokes equations. While featuring a general organisation similar to the axial velocity field, passive scalar fields show additional energy at small wavenumbers, resulting in a higher degree of mixing and in a k(-4/3) spectral inertial range. The DNS results highlight logarithmic growth of the inner-scaled bulk and mean centreline scalar values with the friction Reynolds number, implying an estimated scalar von Karman constant k theta approximate to 0.459, which also nicely fits the mean scalar profiles. The DNS data are used to synthesise a modified form of the classical predictive formula of Kader & Yaglom (Intl J. Heat Mass Transfer, vol. 15 (12), 1972, pp. 2329-2351), which points to some shortcomings of the original formulation. Universality of the mean core scalar profile in defect form is recovered, with very nearly parabolic shape. Logarithmic growth of the buffer-layer peak of the scalar variance is found in the Reynolds number range under scrutiny, which well conforms with Townsend's attached-eddy hypothesis, whose validity is also supported by the spectral maps. The behaviour of the turbulent Prandtl number shows good universality in the outer wall layer, with values Pr-t approximate to 0.84, as also found in previous studies, but closer to unity near the wall, where existing correlations do not reproduce the observed trends.