Bernhard Mehlig

Caustics in turbulent aerosols

Networks of caustics can occur in the distribution of particles suspended in a randomly moving gas. These can facilitate coagulation of particles by bringing them into close proximity, even in cases where the trajectories do not coalesce. The evolution of these caustic patterns depends upon the Lyapunov exponents λ1, λ2 of the suspended particles, as well as the rate J at which particles encounter caustics. We develop a theory determining the quantities J, λ1, λ2 from the statistical properties of the gas flow, in the limit of short correlation times.

Caustic activation of rain showers

We show quantitatively how the collision rate of droplets of visible moisture in turbulent air increases very abruptly as the intensity of the turbulence passes a threshold, due to the formation of fold caustics in their velocity field. The formation of caustics is an activated process, in which a measure of the intensity of the turbulence, termed the Stokes number St, is analogous to temperature in a chemical reaction: the rate of collision contains a factor exp(-C/St). Our results are relevant to the long-standing problem of explaining the rapid onset of rainfall from convecting clouds. Our theory does not involve spatial clustering of particles.

Extreme Female Promiscuity in a Non-Social Invertebrate Species

Background While males usually benefit from as many matings as possible, females often evolve various methods of resistance to matings. The prevalent explanation for this is that the cost of additional matings exceeds the benefits of receiving sperm from a large number of males. Here we demonstrate, however, a strongly deviating pattern of polyandry. Methodology/Principal Findings We analysed paternity in the marine snail Littorina saxatilis by genotyping large clutches (53–79) of offspring from four females sampled in their natural habitats. We found evidence of extreme promiscuity with 15–23 males having sired the offspring of each female within the same mating period. Conclusions/Significance Such a high level of promiscuity has previously only been observed in a few species of social insects. We argue that genetic bet-hedging (as has been suggested earlier) is unlikely to explain such extreme polyandry. Instead we propose that these high levels are examples of convenience polyandry: females accept high numbers of matings if costs of refusing males are higher than costs of accepting superfluous matings.

Path coalescence transition and its applications

We analyze the motion of a system of particles subjected to a random force fluctuating in both space and time, and experiencing viscous damping. When the damping exceeds a certain threshold, the system undergoes a phase transition; the particle trajectories coalesce. We analyze this transition by mapping it to a Kramers problem which we solve exactly. In the limit of weak random force we characterize the dynamics by computing the rate at which caustics are crossed, and the statistics of the particle density in the coalescing phase. Last but not least we describe possible realizations of the effect, ranging from trajectories of raindrops on perspex surfaces to animal migration patterns.

Distribution of relative velocities in turbulent aerosols.

We compute the distribution of relative velocities for a one-dimensional model of heavy particles suspended in a turbulent flow, quantifying the caustic contribution to the moments of relative velocities. The same principles determine the corresponding caustic contribution in d spatial dimensions. The distribution of relative velocities Δv at small separations R acquires the universal form ρ(Δv,R)∼R(d-1)|Δv|(D(2)-2d) for large (but not too large) values of |Δv|. Here D(2) is the phase-space correlation dimension. Our conclusions are in excellent agreement with numerical simulations of particles suspended in a randomly mixing flow in two dimensions, and in quantitative agreement with published data on direct numerical simulations of particles in turbulent flows.

Clustering by Mixing Flows

We calculate the Lyapunov exponents for particles suspended in a random three-dimensional flow, concentrating on the limit where the viscous damping rate is small compared to the inverse correlation time. In this limit Lyapunov exponents are obtained as a power series in epsilon, a dimensionless measure of the particle inertia. Although the perturbation generates an asymptotic series, we obtain accurate results from a Padé-Borel summation. Our results prove that particles suspended in an incompressible random mixing flow can show pronounced clustering when the Stokes number is large and we characterize two distinct clustering effects which occur in that limit.

Shape-dependence of particle rotation in isotropic turbulence

We consider the rotation of neutrally buoyant axisymmetric particles suspended in isotropic turbulence. Using laboratory experiments as well as numerical and analytical calculations, we explore how particle rotation depends upon particle shape. We find that shape strongly affects orientational trajectories, but that it has negligible effect on the variance of the particle angular velocity. Previous work has shown that shape significantly affects the variance of the tumbling rate of axisymmetric particles. It follows that shape affects the spinning rate in a way that is, on average, complementary to the shape-dependence of the tumbling rate. We confirm this relationship using direct numerical simulations, showing how tumbling rate and spinning rate variances show complementary trends for rod-shaped and disk-shaped particles. We also consider a random but non-turbulent flow. This allows us to explore which of the features observed for rotation in turbulent flow are due to the effects of particle alignment i...

Tumbling of small axisymmetric particles in random and turbulent flows.

We analyze the tumbling of small nonspherical, axisymmetric particles in random and turbulent flows. We compute the orientational dynamics in terms of a perturbation expansion in the Kubo number, and obtain the tumbling rate in terms of Lagrangian correlation functions. These capture preferential sampling of the fluid gradients, which in turn can give rise to differences in the tumbling rates of disks and rods. We show that this is a weak effect in Gaussian random flows. But in turbulent flows persistent regions of high vorticity cause disks to tumble much faster than rods, as observed in direct numerical simulations [S. Parsa, E. Calzavarini, F. Toschi, and G. A. Voth, Phys. Rev. Lett. 109, 134501 (2012)]. For larger particles (at finite Stokes numbers), rotational and translational inertia affects the tumbling rate and the angle at which particles collide, due to the formation of rotational caustics.

Coagulation by random velocity fields as a Kramers problem.

We analyze the motion of a system of particles suspended in a fluid which has a random velocity field. There are coagulating and noncoagulating phases. We show that the phase transition is related to a Kramers problem, and we use this to determine the phase diagram in two dimensions, as a function of the dimensionless inertia of the particles, epsilon, and a measure of the relative intensities of potential and solenoidal components of the velocity field, Gamma. We find that the phase line is described by a function which is nonanalytic at epsilon=0, and which is related to escape over a barrier in the Kramers problem. We discuss the physical realizations of this phase transition.

Rotation of a spheroid in a simple shear at small Reynolds number

We derive an effective equation of motion for the orientational dynamics of a neutrally buoyant spheroid suspended in a simple shear flow, valid for arbitrary particle aspect ratios and to linear order in the shear Reynolds number. We show how inertial effects lift the degeneracy of the Jeffery orbits and determine the stabilities of the log-rolling and tumbling orbits at infinitesimal shear Reynolds numbers. For prolate spheroids we find stable tumbling in the shear plane, log-rolling is unstable. For oblate particles, by contrast, log-rolling is stable and tumbling is unstable provided that the aspect ratio is larger than a critical value. When the aspect ratio is smaller than this value tumbling turns stable, and an unstable limit cycle is born.