Countergradient turbulent transport in a plume with a crossflow

Direct numerical simulation of a turbulent forced buoyant plume in a crossflow is performed at a source Reynolds number Re0=1000\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{ Re }}_0=1000$$\end{document}, Richardson number Ri0=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{{Ri}}_0=1$$\end{document}, Prandtl number Pr=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{{Pr}}=1$$\end{document} and source-to-crossflow velocity ratio R0=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0=1$$\end{document}. The instantaneous and temporally averaged flow fields are assessed in detail, providing an overview of the flow dynamics. The velocity, temperature and pressure fields are used together with enstrophy fields to describe qualitatively the evolution of the plume as it is swept downstream by the crossflow, and the mechanisms involved in its evolution are outlined. The plume trajectory is determined quantitatively in a number of ways, and it is shown that the central streamline and the centre of buoyancy of the plume differ significantly-as with jets in crossflow, the central streamline is seen to follow the top of the plume, whereas the centre of buoyancy, by definition, describes the plume as a whole. We then investigate the turbulence properties inside the plume; in particular the eddy viscosity and diffusivity are presented, which are significant parameters in turbulence modelling. Assessment of turbulence production demonstrates the presence of regions where turbulence kinetic energy is redistributed to the kinetic energy of the mean flow, implying a negative eddy viscosity within certain regions of the domain. Similarly, the observation that the buoyancy flux and buoyancy gradient are anti-parallel in specific regions of the flow implies a negative eddy diffusivity in said regions, which must be realised in models of such flows in order to capture the countergradient transport of thermal properties. A characteristic eddy viscosity and diffusivity are presented, and shown to be approximately constant in the fully developed regime, resulting in a constant characteristic turbulent Prandtl number, in turn signifying self-similarity.