Small heavy particles in a fluid flow respond to the flow on a timescale proportional to their inertia or Stokes number St. Their behavior has been thought to be gradually modified as St increases. We show, on the other hand, in the steady spatially periodic laminar Taylor Green vortex flow, that particle dynamics, and their effective diffusivity, actually change in an irregular, nonmonotonic, and sometimes discontinuous manner with increasing St. At St similar to 1, we show chaotic particle motion, contrasting with earlier conclusions for heavy particles in the same flow [Wang et al., Phys. Fluids 4, 1789 (1992)]. Particles may display trapped orbits, or unbounded diffusive or ballistic dispersion, with the vortices behaving like scatterers in a soft Lorentz gas [Klages et al., Phys. Rev. Lett. 122, 064102 (2019)]. The dynamics is nonergodic.

We use direct numerical simulations to study convection in rotating Rayleigh-Bénard convection in horizontally confined geometries of a given aspect ratio, with the walls held at fixed temperatures. We show that this arrangement is unconditionally unstable to flow that takes the form of wall-adjacent convection rolls. For wall temperatures close to the temperatures of the upper or lower boundaries, we show that the base state undergoes a Hopf bifurcation to a state comprised of spatiotemporal oscillations - 'wall modes' - precessing in a retrograde direction. We study the saturated nonlinear state of these modes, and show that the velocity boundary conditions at the upper and lower boundaries are crucial to the formation and propagation of the wall modes: asymmetric velocity boundary conditions at the upper and lower boundaries can lead to prograde wall modes, while stress-free boundary conditions at both walls can lead to wall modes that have no preferred direction of propagation.

We study the orientation dynamics of two-dimensional concavo-convex solid bodies that are denser than the fluid through which they fall under gravity. We show that the orientation dynamics of the body, quantified in terms of the angle 𝜙 relative to the horizontal, undergoes a transcritical bifurcation at a Reynolds number and a subcritical pitchfork bifurcation at a Reynolds number . For , the concave-downwards orientation of 𝜙=0 is unstable and bodies overturn into the 𝜙=𝜋 orientation. For , the falling body has two stable equilibria at 𝜙=0and𝜙=𝜋 for steady descent. For , the concave-downwards orientation of 𝜙=0 is again unstable and bodies that start concave-downwards exhibit overstable oscillations about the unstable fixed point, eventually tumbling into the stable 𝜙=𝜋 orientation. The at which the subcritical pitchfork bifurcation occurs is distinct from the 𝑅⁢𝑒 for the onset of vortex shedding, which causes the 𝜙=𝜋 equilibrium to also become unstable, with bodies fluttering about 𝜙=𝜋. The complex orientation dynamics of irregularly shaped bodies evidenced here are relevant in a wide range of settings, from the tumbling of hydrometeors to the settling of mollusk shells.

We study entrainment in dry thermals in neutrally and unstably stratified ambients, and moist thermals in dry-neutrally stratified ambients using direct numerical simulations. We find, in agreement with results of Lecoanet and Jeevanjee [J. Atmos. Sci. 76, 3785 (2019)], that turbulence plays a minor role in entrainment in dry thermals in a neutral ambient for Reynolds numbers Re≲10⁴. We then show that the net entrainment rate increases when the buoyancy of the thermals increases, either by condensation heating or because of an unstably stratified ambient. This is in contrast with the findings of Morrison et al. [J. Atmos. Sci. 78, 797 (2021)]. We also show that the role of turbulence is greater in these cases than in dry thermals and, significantly, that the combined action of condensation heating and turbulence creates intense small-scale vorticity, destroying the coherent vortex ring that is seen in dry and moist laminar thermals. These findings suggest that fully resolved simulations at Reynolds numbers significantly larger than the mixing transition Reynolds number Re=10⁴ are necessary to understand the role of turbulence in the entrainment in growing cumulus clouds, which consist of a series of thermals rising and decaying in succession.

We study, by direct numerical simulations in two and three dimensions, the instability caused by the settling and evaporation of water droplets out of a cloudy layer saturated with vapor into a dry subcloud ambient under conditions where mammatus clouds were shown to form but with the addition of background shear. We show that shear changes the type of cloud formation qualitatively from mammatus-like to a newly identified cloud type called asperitas. Intermediate levels of shear are shown to be needed. Shear suppresses the growth of small-scale perturbations, giving rise to smooth, long-lasting structures and smaller rates of mixing. Three-dimensionality is shown to make a qualitative difference, unlike in mammatus clouds. We also show that under non-cloud-like conditions, the instability can be very different.

We use well-resolved numerical simulations to study the combined effects of buoyancy, pressure-driven shear and rotation on the melt rate and morphology of a layer of pure solid overlying its liquid phase in three dimensions at a Rayleigh number Ra = 1.25 x 10(5) . During thermal convection, we find that the rate of melting of the solid phase varies non-monotonically with the strength of the imposed shear flow. In the absence of rotation, depending on whether buoyancy or shear dominates the flow, we observe either domes or ridges aligned in the direction of the shear flow, respectively. Furthermore, we show that the geometry of the phase boundary has important effects on the magnitude and evolution of the heat flux in the liquid layer. In the presence of rotation, the strength of which is characterized by the Rossby number, Ro, we observe that for Ro = O(1), the mean flow in the interior is perpendicular to the direction of the constant horizontal applied pressure gradient. As the magnitude of this pressure gradient increases, the geometry of solid-liquid interface evolves from the voids characteristic of melting by rotating convection, to grooves oriented perpendicular or obliquely to the direction of the pressure gradient.

We study the role of phase change and thermal noise in particle transport in turbulent flows. We employ a toy model to extract the main physics: Condensing droplets are modelled as heavy particles which grow in size, the ambient flow is modelled as a two-dimensional Taylor-Green flow consisting of an array of vortices delineated by separatrices, and thermal noise are modelled as uncorrelated Gaussian white noise. In general, heavy inertial particles are centrifuged out of regions of high vorticity and into regions of high strain. In cellular flows, we find, in agreement with earlier results, that droplets with Stokes numbers smaller than a critical value, St<Stcr, remain trapped in the vortices in which they are initialized, while larger droplets move ballistically away from their initial positions by crossing separatrices. We independently vary the Péclet number Pe characterizing the amplitude of thermal noise and the condensation rate Π to study their effects on the critical Stokes number for droplet trapping, as well as on the final states of motion of the droplets. We find that the imposition of thermal noise, or of a finite condensation rate, allows droplets of St<Stcr to leave their initial vortices. We find that the effects of thermal noise become negligible for growing droplets and that growing droplets achieve ballistic motion when their Stokes numbers become O(1). We also find an intermediate regime prior to attaining the ballistic state, in which droplets move diffusively away from their initial vortices in the presence of thermal noise.

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.

This article describes the dynamics of small inertial particles centrifuging out of a single vortex. It shows the importance of caustics formation in the vicinity of a single vortex: both for particle collisions and void formation. From these single-vortex studies we provide estimates of the role of caustics in high Reynolds number turbulence, and in the case of clouds, estimate how they may help in rain initiation by bridging the droplet-growth bottleneck. We briefly describe how the Basset-Boussinesq history force may be calculated by a method which does not involve huge memory costs, and we provide arguments for its possible importance for droplets in turbulence. We discuss how phase change could render cloud turbulence fundamentally different from turbulence in other situations.

The COVID-19 pandemic has inspired several studies on the fluid dynamics of respiratory events. Here, we propose a computational approach in which respiratory droplets are coarse-grained into an Eulerian liquid field advected by the fluid streamlines. A direct numerical simulation is carried out for a moist cough using a closure model for space-time dependence of the evaporation time scale. Stokes-number estimates are provided, for the initial droplet size of 10 mu m, which are found to be ≪1, thereby justifying the neglect of droplet inertia, over the duration of the simulation. Several important features of the moist-cough flow reported in the literature using Lagrangian tracking methods have been accurately captured using our scheme. Some new results are presented, including the evaporation time for a 'mild' cough, a saturation-temperature diagram and a favourable correlation between the vorticity and liquid fields. The present approach can be extended for studying the long-range transmission of virus-laden droplets.