The SARS-CoV-2 is transmitted not only through coughing, but also through breathing, speaking or singing. We perform direct numerical simulations of the turbulent transport of potentially infectious aerosols in short conversations, involving repetitive phrases separated by quiescent intervals. We estimate that buoyancy effects due to droplet evaporation are small, and neglect them. A two-way conversation is shown to significantly reduce the aerosol exposure compared with a relative monologue by one person and relative silence of the other. This is because of the ‘cancelling’ effect produced by the two interacting speech jets. Unequal conversation is shown to significantly increase the infection risk to the person who talks less. Interestingly, a small height difference is worse for infection spread, due to reduced interference between the speech jets, than two faces at the same level. For small axial separation, speech jets show large oscillations and reach the other person intermittently. We suggest a range of lateral separations between two people to minimize transmission risk. A realistic estimate of the infection probability is provided by including exposure through the eyes and mouth, in addition to the more common method of using inhaled virions alone. We expect that our results will provide useful inputs to epidemiological models and to disease management.

We investigate statistically steady states of turbulent flows when molecular degrees of freedom, in particular vibration, are taken into account. Unlike laminar flows initially in thermal non-equilibrium which asymptotically relax towards thermal equilibrium, turbulent flows present persistent departures from thermal equilibrium. This is due to fluctuations in temperature and other thermodynamic variables, which are known to increase with turbulent Mach number. Analytical results are compared to direct numerical simulations at a range of Reynolds and Mach numbers as well as molecular parameters such as relaxation times. Turbulent fluctuations are also shown to disrupt the distribution of energy between translational–rotational–vibrational modes even if thermal equilibrium is attained instantaneously relative to turbulence time scales, an effect that increases with characteristic relaxation times. Because of the nonlinear relation between temperature and vibrational energy in equilibrium, the fluctuation of the latter could be strongly positively skewed with long tails in its probability density function. This effect is stronger in flows with strong temperature fluctuations and when vibrational modes are partially excited. Because of the finite-time relaxation of vibration, departures from equilibrium result in very strong transfers of energy from the translational–rotational mode to the vibrational mode. A simple spectral model can explain the stronger departures from thermal equilibrium observed at the small scales. The spectral behaviour of the instantaneous vibrational energy can be described by classical phenomenology for passive scalars.

The effects of rough surfaces on turbulent channel flow are modelled by an extra force term in the Navier–Stokes equations. This force term contains two parameters, related to the density and the height of the roughness elements, and a shape function, which regulates the influence of the force term with respect to the distance from the channel wall. This permits a more flexible specification of a rough surface than a single parameter such as the equivalent sand grain roughness. The effects of the roughness force term on turbulent channel flow have been investigated for a large number of parameter combinations and several shape functions by direct numerical simulations. It is possible to cover the full spectrum of rough flows ranging from hydraulically smooth through transitionally rough to fully rough cases. By using different parameter combinations and shape functions, it is possible to match the effects of different types of rough surfaces. Mean flow and standard turbulence statistics have been used to compare the results to recent experimental and numerical studies and a good qualitative agreement has been found. Outer scaling is preserved for the streamwise velocity for both the mean profile as well as its mean square fluctuations in all but extremely rough cases. The structure of the turbulent flow shows a trend towards more isotropic turbulent states within the roughness layer. In extremely rough cases, spanwise structures emerge near the wall and the turbulent state resembles a mixing layer. A direct comparison with the study of Ashrafian, Andersson & Manhart (Intl J. Heat Fluid Flow, vol. 25, 2004, pp. 373–383) shows a good quantitative agreement of the mean flow and Reynolds stresses everywhere except in the immediate vicinity of the rough wall. The proposed roughness force term may be of benefit as a wall model for direct and large-eddy numerical simulations in cases where the exact details of the flow over a rough wall can be neglected.

Well-resolved large-eddy simulations (LES) are performed in order to investigate flow phenomena and turbulence structure of the boundary layer along a supersonic compression ramp. The numerical simulations directly reproduce an available experimental result. The compression ramp has a deflection angle of $\beta\,{=}\,25^\circ$. The mean free-stream Mach number is $M_\infty\,{=}\,2.95$. The Reynolds number based on the incoming boundary-layer thickness is $Re_{\delta_0}\,{=}\,63\,560$ in accordance with the reference experiment. These simulations overcome deficiencies of earlier direct numerical simulations (DNS) and LES in terms of ramp-deflection angle, Reynolds number and spanwise size of the computational domain which is required for capturing the essential flow phenomena. The filtered conservation equations for mass, momentum and energy are solved with a high-order finite-difference scheme. The effect of subgrid scales is modelled by the approximate deconvolution model. About $18.5\,{\times}\,10^6$ grid points are used for discretizing the computational domain. To obtain mean flow and turbulence structure the flow is sampled 1272 times over 703 characteristic time scales of the incoming boundary layer. Statistical data are computed from these samples. An analysis of the data shows good agreement with the experiment in terms of mean quantities such as shock position, separation and reattachment location, skin-friction and surface-pressure distributions, and turbulence structure. The computational data confirm theoretical and experimental results on fluctuation amplification across the interaction region. In the wake of the main shock a shedding of shocklets is observed. The temporal behaviour of the coupled shock–separation system agrees well with experimental data. Unlike previous DNS the present simulation data provide indications of a large-scale shock motion. Also, evidence for the existence of three-dimensional large-scale streamwise structures, commonly referred to as Görtler-like vortices, is found.