-
Short term vs. long term: optimization of microswimmer navigation on different time horizons
Authors:
Navid Mousavi,
Jingran Qiu,
Lihao Zhao,
Bernhard Mehlig,
Kristian Gustavsson
Abstract:
We use reinforcement learning to find strategies that allow microswimmers in turbulence to avoid regions of large strain. This question is motivated by the hypothesis that swimming microorganisms tend to avoid such regions to minimise the risk of predation. We ask which local cues a microswimmer must measure to efficiently avoid such straining regions. We find that it can succeed without direction…
▽ More
We use reinforcement learning to find strategies that allow microswimmers in turbulence to avoid regions of large strain. This question is motivated by the hypothesis that swimming microorganisms tend to avoid such regions to minimise the risk of predation. We ask which local cues a microswimmer must measure to efficiently avoid such straining regions. We find that it can succeed without directional information, merely by measuring the magnitude of the local strain. However, the swimmer avoids straining regions more efficiently if it can measure the sign of local strain gradients. We compare our results with those of an earlier study [Mousavi et al. arxiv:2309.09641] where a short-time expansion was used to find optimal strategies. We find that the short-time strategies work well in some cases but not in others. We derive a new theory that explains when the time-horizon matters for our optimisation problem, and when it does not. We find the strategy with best performance when the time-horizon coincides with the correlation time of the turbulent fluctuations. We also explain how the update frequency (the frequency at which the swimmer updates its state) affects the found strategies. We find that higher update frequencies yield better performance, as long as the time between updates is smaller than the correlation time of the flow.
△ Less
Submitted 30 April, 2024;
originally announced April 2024.
-
Efficient survival strategy for zooplankton in turbulence
Authors:
Navid Mousavi,
Jingran Qiu,
Bernhard Mehlig,
Lihao Zhao,
Kristian Gustavsson
Abstract:
Zooplankton in a quiescent environment can detect predators by hydrodynamic sensing, triggering powerful escape responses. Since turbulent strain tends to mask the hydrodynamic signal, the organisms should avoid such regions, but it is not known how they accomplish this. We found a simple, robust, and highly efficient strategy, that relies on measuring the sign of gradients of squared strain. Plan…
▽ More
Zooplankton in a quiescent environment can detect predators by hydrodynamic sensing, triggering powerful escape responses. Since turbulent strain tends to mask the hydrodynamic signal, the organisms should avoid such regions, but it is not known how they accomplish this. We found a simple, robust, and highly efficient strategy, that relies on measuring the sign of gradients of squared strain. Plankton following this strategy show very strong spatial clustering, and align against the local flow velocity, facilitating mate finding and feeding. The strategy has the potential to reconcile competing fitness pressures.
△ Less
Submitted 8 April, 2024; v1 submitted 18 September, 2023;
originally announced September 2023.
-
Lagrangian supersaturation fluctuations at the cloud edge
Authors:
J. Fries,
G. Sardina,
G. Svensson,
A. Pumir,
B. Mehlig
Abstract:
Evaporation of cloud droplets accelerates when turbulence mixes dry air into the cloud, affecting droplet-size distributions in atmospheric clouds, combustion sprays, and jets of exhaled droplets. The challenge is to model local correlations between droplet numbers, sizes, and supersaturation, which determine supersaturation fluctuations along droplet paths (Lagrangian fluctuations). We derived a…
▽ More
Evaporation of cloud droplets accelerates when turbulence mixes dry air into the cloud, affecting droplet-size distributions in atmospheric clouds, combustion sprays, and jets of exhaled droplets. The challenge is to model local correlations between droplet numbers, sizes, and supersaturation, which determine supersaturation fluctuations along droplet paths (Lagrangian fluctuations). We derived a statistical model that accounts for these correlations. Its predictions are in quantitative agreement with results of direct numerical simulations, and it explains the key mechanisms at play.
△ Less
Submitted 6 September, 2023;
originally announced September 2023.
-
Caustic formation in a non-Gaussian model for turbulent aerosols
Authors:
J. Meibohm,
L. Sundberg,
B. Mehlig,
K. Gustavsson
Abstract:
Caustics in the dynamics of heavy particles in turbulence accelerate particle collisions. The rate $\mathscr{J}$ at which these singularities form depends sensitively on the Stokes number St, the non-dimensional inertia parameter. Exact results for this sensitive dependence have been obtained using Gaussian statistical models for turbulent aerosols. However, direct numerical simulations of heavy p…
▽ More
Caustics in the dynamics of heavy particles in turbulence accelerate particle collisions. The rate $\mathscr{J}$ at which these singularities form depends sensitively on the Stokes number St, the non-dimensional inertia parameter. Exact results for this sensitive dependence have been obtained using Gaussian statistical models for turbulent aerosols. However, direct numerical simulations of heavy particles in turbulence yield much larger caustic-formation rates than predicted by the Gaussian theory. In order to understand possible mechanisms explaining this difference, we analyse a non-Gaussian statistical model for caustic formation in the limit of small St. We show that at small St, $\mathscr{J}$ depends sensitively on the tails of the distribution of Lagrangian fluid-velocity gradients. This explains why different authors obtained different St-dependencies of $\mathscr{J}$ in numerical-simulation studies. The most-likely gradient fluctuation that induces caustics at small St, by contrast, is the same in the non-Gaussian and Gaussian models. Direct-numerical simulation results for particles in turbulence show that the optimal fluctuation is similar, but not identical, to that obtained by the model calculations.
△ Less
Submitted 15 May, 2024; v1 submitted 20 July, 2023;
originally announced July 2023.
-
Statistical models for the dynamics of heavy particles in turbulence
Authors:
J. Bec,
K. Gustavsson,
B. Mehlig
Abstract:
When very small particles are suspended in a fluid in motion, they tend to follow the flow. How such tracer particles are mixed, transported, and dispersed by turbulent flow has been successfully described by statistical models. Heavy particles, with mass densities larger than that of the carrying fluid, can detach from the flow. This results in preferential sampling, small-scale fractal clusterin…
▽ More
When very small particles are suspended in a fluid in motion, they tend to follow the flow. How such tracer particles are mixed, transported, and dispersed by turbulent flow has been successfully described by statistical models. Heavy particles, with mass densities larger than that of the carrying fluid, can detach from the flow. This results in preferential sampling, small-scale fractal clustering, and large collision velocities. To describe these effects of particle inertia, it is necessary to consider both particle positions and velocities in phase space. In recent years, statistical phase-space models have significantly contributed to our understanding of inertial-particle dynamics in turbulence. These models help to identify the key mechanisms and non-dimensional parameters governing the particle dynamics, and have made qualitative, and in some cases quantitative predictions. This article reviews statistical phase-space models for the dynamics of small, yet heavy, spherical particles in turbulence. We evaluate their effectiveness by comparing their predictions with results from numerical simulations and laboratory experiments, and summarise their successes and failures. Annu. Rev. Fluid Mech. 56: In press. DOI: 10.1146/annurev-fluid-032822-014140.
△ Less
Submitted 12 September, 2023; v1 submitted 3 April, 2023;
originally announced April 2023.
-
Inertial angular dynamics of non-spherical atmospheric particles
Authors:
Taraprasad Bhowmick,
Johannes Seesing,
Kristian Gustavsson,
Johannes Guettler,
Alain Pumir,
Bernhard Mehlig,
Yong Wang,
Gholamhossein Bagheri
Abstract:
Cloud-ice crystals, volcanic ash, and microplastic are ubiquitous in the atmosphere. These non-spherical particles are small, but their mass density is much greater than that of air. Little is known about their inertial dynamics, mainly because experiments with such heavy, sub-millimetre particles in air are difficult. We tracked the inertial dynamics of heavy sub-millimetre spheroids through stil…
▽ More
Cloud-ice crystals, volcanic ash, and microplastic are ubiquitous in the atmosphere. These non-spherical particles are small, but their mass density is much greater than that of air. Little is known about their inertial dynamics, mainly because experiments with such heavy, sub-millimetre particles in air are difficult. We tracked the inertial dynamics of heavy sub-millimetre spheroids through still air and observed that their orientations fluctuate considerably, in stark contrast to the rapid alignment seen in high-density fluids such as water. A model, that quantitatively describes the resulting transient oscillations of the particle orientation, shows that the oscillations are due to particle inertia, and allows us to study the effect of particle shape and volume, beyond the parameters of the experiment. We discuss implications for the angular dynamics of such particles in turbulent air. We conclude that the particle inertia can significantly delay the alignment and increase angular fluctuations. This has significant implications for the statistics of particle orientation, affecting settling velocities and atmospheric residence times, collision/aggregation mechanism and how the particles scatter and absorb solar radiation.
△ Less
Submitted 7 March, 2023;
originally announced March 2023.
-
Hydrodynamic force on a small squirmer moving with a time-dependent velocity at small Reynolds numbers
Authors:
T. Redaelli,
F. Candelier,
R. Mehaddi,
C. Eloy,
B. Mehlig
Abstract:
We calculate the hydrodynamic force on a small spherical, unsteady squirmer moving with a time-dependent velocity in a fluid at rest, taking into account convective and unsteady fluid inertia effects in perturbation theory. Our results generalise those of Lovalenti and Brady (1993) from passive to active spherical particles. We find that convective inertia changes the history-contribution to the h…
▽ More
We calculate the hydrodynamic force on a small spherical, unsteady squirmer moving with a time-dependent velocity in a fluid at rest, taking into account convective and unsteady fluid inertia effects in perturbation theory. Our results generalise those of Lovalenti and Brady (1993) from passive to active spherical particles. We find that convective inertia changes the history-contribution to the hydrodynamic force, as it does for passive particles. We determine how the hydrodynamic force depends on the swimming gait of the unsteady squirmer. Since swimming breaks the spherical symmetry of the problem, the force is not completely determined by the outer solution of the asymptotic matching problem, as it is for passive spheres. There are additional contributions due to the inhomogeneous solution of the inner problem. We also compute the disturbance flow, illustrating convective and unsteady effects when the particle experiences a sudden start followed by a sudden stop.
△ Less
Submitted 19 December, 2023; v1 submitted 16 September, 2022;
originally announced September 2022.
-
Critical charges for droplet collisions
Authors:
A. Dubey,
G. P. Bewley,
K. Gustavsson,
B. Mehlig
Abstract:
The collision efficiency of uncharged micron-sized water droplets in air is determined by the breakdown of hydrodynamics at droplet separations of the order of the mean-free path, by van-der-Waals forces, or a combination of the two. In contrast, electrostatic forces determine the collision efficiency of charged droplets if the charge is large enough. To find the charge for which the transition to…
▽ More
The collision efficiency of uncharged micron-sized water droplets in air is determined by the breakdown of hydrodynamics at droplet separations of the order of the mean-free path, by van-der-Waals forces, or a combination of the two. In contrast, electrostatic forces determine the collision efficiency of charged droplets if the charge is large enough. To find the charge for which the transition to charge-dominated collisions occurs, we computed the collision efficiency of charged, hydrodynamically-interacting droplets settling in quiescent air, including the breakdown of hydrodynamics at small interfacial distances. For oppositely charged droplets, the transition occurs when a saddle point of the relative droplet-dynamics exits the region where the hydrodynamics breaks down. For droplets with radii $16\,μ$m and $20\,μ$m, this occurs at $\sim 10^3$ elementary charges $e$. For smaller charges, the collision efficiency depends upon the Kn number (defined as the ratio of the mean-free-path of air to the mean droplet radius), whereas for larger charges it does not. For droplets charged with the same polarity, the critical charge is $\sim 10^4\,e$ for the above radii.
△ Less
Submitted 12 September, 2022;
originally announced September 2022.
-
Inertial torque on a squirmer
Authors:
F. Candelier,
J. Qiu,
L. Zhao,
G. Voth,
B. Mehlig
Abstract:
A small spheroid settling in a quiescent fluid experiences an inertial torque that aligns it so that it settles with its broad side first. Here we show that an active particle experiences such a torque too, as it settles in a fluid at rest. For a spherical squirmer, the torque is $\boldsymbol{T}^\prime = -{\tfrac{9}{8}} m_f (\boldsymbol{v}_s^{(0)} \wedge \boldsymbol{v}_g^{(0)})$ where…
▽ More
A small spheroid settling in a quiescent fluid experiences an inertial torque that aligns it so that it settles with its broad side first. Here we show that an active particle experiences such a torque too, as it settles in a fluid at rest. For a spherical squirmer, the torque is $\boldsymbol{T}^\prime = -{\tfrac{9}{8}} m_f (\boldsymbol{v}_s^{(0)} \wedge \boldsymbol{v}_g^{(0)})$ where $\boldsymbol{v}_s^{(0)}$ is the swimming velocity, $\boldsymbol{v}_g^{(0)}$ is the settling velocity in the Stokes approximation, and $m_f$ is the equivalent fluid mass. This torque aligns the swimming direction against gravity: swimming up is stable, swimming down is unstable.
△ Less
Submitted 7 September, 2022;
originally announced September 2022.
-
Second-order inertial forces and torques on a sphere in a viscous steady linear flow
Authors:
Fabien Candelier,
Rabah Mehaddi,
Bernhard Mehlig,
Jacques Magnaudet
Abstract:
We compute the full set of second-order inertial corrections to the instantaneous force and torque acting on a small spherical rigid particle moving unsteadily in a general steady linear flow. This is achieved by using matched asymptotic expansions and formulating the problem in a coordinate system co-moving with the background flow. Effects of the fluid-velocity gradients are assumed to be small,…
▽ More
We compute the full set of second-order inertial corrections to the instantaneous force and torque acting on a small spherical rigid particle moving unsteadily in a general steady linear flow. This is achieved by using matched asymptotic expansions and formulating the problem in a coordinate system co-moving with the background flow. Effects of the fluid-velocity gradients are assumed to be small, but to dominate over those of the velocity difference between the body and fluid, which makes the results essentially relevant to nearly neutrally buoyant particles. The outer solution (which at first order is responsible for the Basset-Boussinesq history force at short time and for shear-induced forces such as the Saffman lift force at long time) is expressed via a flow-dependent tensorial kernel. The second-order inner solution brings a number of different contributions to the force and torque. Some are proportional to the relative translational or angular acceleration between the particle and fluid, while others take the form of products of the rotation/strain rate of the background flow and the relative translational or angular velocity between the particle and fluid. Adding the outer and inner contributions, the known added-mass force or the spin-induced lift force are recovered, and new effects involving the rotation/strain rate of the background flow are revealed. The resulting force and torque equations provide a rational extension of the classical Basset-Boussinesq-Oseen equation incorporating all first- and second-order fluid inertia effects resulting from both unsteadiness and velocity gradients of the carrying flow.
△ Less
Submitted 17 November, 2022; v1 submitted 24 August, 2022;
originally announced August 2022.
-
Caustics in turbulent aerosols form along the Vieillefosse line at weak particle inertia
Authors:
Jan Meibohm,
Kristian Gustavsson,
Bernhard Mehlig
Abstract:
Caustic singularities of the spatial distribution of particles in turbulent aerosols enhance collision rates and accelerate coagulation. Here we investigate how and where caustics form at weak particle inertia, by analysing a three-dimensional Gaussian statistical model for turbulent aerosols in the persistent limit, where the flow varies slowly compared with the particle relaxation time. In this…
▽ More
Caustic singularities of the spatial distribution of particles in turbulent aerosols enhance collision rates and accelerate coagulation. Here we investigate how and where caustics form at weak particle inertia, by analysing a three-dimensional Gaussian statistical model for turbulent aerosols in the persistent limit, where the flow varies slowly compared with the particle relaxation time. In this case, correlations between particle- and fluid-velocity gradients are strong, and caustics are induced by large, strain-dominated excursions of the fluid-velocity gradients. These excursions must cross a characteristic threshold in the plane spanned by the invariants $Q$ and $R$ of the fluid-velocity gradients. Our method predicts that the most likely way to reach this threshold is by a unique ``optimal fluctuation'' that propagates along the Vieillefosse line, $27R^2/4 +Q^3=0$. We determine the shape of the optimal fluctuation as a function of time and show that it is dominant in numerical statistical-model simulations even for moderate particle inertia.
△ Less
Submitted 16 February, 2023; v1 submitted 28 July, 2022;
originally announced July 2022.
-
Bifurcations in droplet collisions
Authors:
A. Dubey,
K. Gustavsson,
G. Bewley,
B. Mehlig
Abstract:
Saffman and Turner (1957) argued that the collision rate for droplets in turbulence increases as the turbulent strain rate increases. But the numerical simulations of Dhanasekaran et al. (2021) in a steady straining flow show that the Saffman-Turner model is oversimplified because it neglects droplet-droplet interactions. These result in a complex dependence of the collision rate on the strain rat…
▽ More
Saffman and Turner (1957) argued that the collision rate for droplets in turbulence increases as the turbulent strain rate increases. But the numerical simulations of Dhanasekaran et al. (2021) in a steady straining flow show that the Saffman-Turner model is oversimplified because it neglects droplet-droplet interactions. These result in a complex dependence of the collision rate on the strain rate and on the differential settling speed. Here we show that this dependence is explained by a sequence of bifurcations in the collision dynamics. We compute the bifurcation diagram when strain is aligned with gravity, and show that it yields important insights into the collision dynamics. First, the steady-state collision rate remains non-zero in the limit Kn $\to0$, contrary to the common assumption that the collision rate tends to zero in this limit (Kn is a non-dimensional measure of the mean free path of air). Second, the non-monotonic dependence of the collision rate on the differential settling speed is explained by a grazing bifurcation. Third, the bifurcation analysis explains why so-called "closed trajectories" appear and disappear. Fourth, our analysis predicts strong spatial clustering near certain saddle points, where the effects of strain and differential settling cancel
△ Less
Submitted 16 January, 2022;
originally announced January 2022.
-
Experimental validation of fluid inertia models for a cylinder settling in a quiescent flow
Authors:
F. Cabrera,
M. Z. Sheikh,
B. Mehlig,
N. Plihon,
M. Bourgoin,
A. Pumir,
A. Naso
Abstract:
The precise description of the motion of anisotropic particles in a flow rests on the understanding of the force and torque acting on them. Here, we study experimentally small, very elongated particles settling in a fluid at small Reynolds number. In our experiments, we can, to a very good approximation, relate the rate of rotation of cylindrical tungsten rods, of aspect ratios β= 8 and β= 16, set…
▽ More
The precise description of the motion of anisotropic particles in a flow rests on the understanding of the force and torque acting on them. Here, we study experimentally small, very elongated particles settling in a fluid at small Reynolds number. In our experiments, we can, to a very good approximation, relate the rate of rotation of cylindrical tungsten rods, of aspect ratios β= 8 and β= 16, settling in pure glycerol to the torque they are experiencing. This allows us to compare the measured torque with expressions obtained either in the slender-rod limit, or in the case of spheroids. Both theories predict a simple angle dependence for the torque, which is found to capture very well the experimental results. Surprisingly, the slender-rod approximation predicts much better the results for β= 8, than for β= 16. In the latter case, the expression obtained for a spheroid provides a better approximation. The translational dynamics is shown to be in qualitative agreement with the slender-rod and spheroid models, the former one being found to represent better the experimental data.
△ Less
Submitted 13 July, 2021;
originally announced July 2021.
-
Collisions of micron-sized, charged water droplets in still air
Authors:
G. Magnusson,
A. Dubey,
R. Kearney,
G. P. Bewley,
B. Mehlig
Abstract:
We investigate the effect of electrical charge on collisions of hydrodynamically interacting, micron-sized water droplets settling through quiescent air. The relative dynamics of charged droplets is determined by hydrodynamic interactions, particle and fluid inertia, and electrostatic forces. We analyse the resulting relative dynamics of oppositely charged droplets by determining its fixed points…
▽ More
We investigate the effect of electrical charge on collisions of hydrodynamically interacting, micron-sized water droplets settling through quiescent air. The relative dynamics of charged droplets is determined by hydrodynamic interactions, particle and fluid inertia, and electrostatic forces. We analyse the resulting relative dynamics of oppositely charged droplets by determining its fixed points and their stable and unstable manifolds. The stable manifold of a saddle point forms a separatrix that separates colliding trajectories from those that do not collide. The qualitative conclusions from this theory are in excellent agreement with experiments.
△ Less
Submitted 22 June, 2021;
originally announced June 2021.
-
Unsteady and inertial dynamics of an active particle in a fluid
Authors:
T. Redaelli,
F. Candelier,
R. Mehaddi,
B. Mehlig
Abstract:
It is well known that the reversibility of Stokes flow makes it difficult for small microorganisms to swim. Inertial effects break this reversibility, allowing new mechanisms of propulsion and feeding. Therefore it is important to understand the effects of unsteady and fluid inertia on the dynamics of microorganisms in flow. In this work, we show how to translate known inertial effects for non-mot…
▽ More
It is well known that the reversibility of Stokes flow makes it difficult for small microorganisms to swim. Inertial effects break this reversibility, allowing new mechanisms of propulsion and feeding. Therefore it is important to understand the effects of unsteady and fluid inertia on the dynamics of microorganisms in flow. In this work, we show how to translate known inertial effects for non-motile organisms to motile ones, from passive to active particles. The method relies on a principle used earlier by Legendre and Magnaudet (1997) to deduce inertial corrections to the lift force on a bubble from the inertial drag on a solid sphere, using the fact that small inertial effects are determined by the far field of the disturbance flow. The method allows for example to compute the inertial effect of unsteady fluid accelerations on motile organisms, and the inertial forces such organisms experience in steady shear flow. We explain why the method fails to describe the effect of convective fluid inertia.
△ Less
Submitted 20 June, 2022; v1 submitted 4 May, 2021;
originally announced May 2021.
-
Navigation of micro-swimmers in steady flow: the importance of symmetries
Authors:
J. Qiu,
N. Mousavi,
K. Gustavsson,
C. Xu,
B. Mehlig,
L. Zhao
Abstract:
Marine microorganisms must cope with complex flow patterns and even turbulence as they navigate the ocean. To survive they must avoid predation and find efficient energy sources. A major difficulty in analysing possible survival strategies is that the time series of environmental cues in non-linear flow is complex, and that it depends on the decisions taken by the organism. One way of determining…
▽ More
Marine microorganisms must cope with complex flow patterns and even turbulence as they navigate the ocean. To survive they must avoid predation and find efficient energy sources. A major difficulty in analysing possible survival strategies is that the time series of environmental cues in non-linear flow is complex, and that it depends on the decisions taken by the organism. One way of determining and evaluating optimal strategies is reinforcement learning. In a proof-of-principle study, Colabrese et al. [Phys. Rev. Lett. (2017)] used this method to find out how a micro-swimmer in a vortex flow can navigate towards the surface as quickly as possible, given a fixed swimming speed. The swimmer measured its instantaneous swimming direction and the local flow vorticity in the laboratory frame, and reacted to these cues by swimming either left, right, up, or down. However, usually a motile microorganism measures the local flow rather than global information, and it can only react in relation to the local flow, because in general it cannot access global information (such as up or down in the laboratory frame). Here we analyse optimal strategies with local signals and actions that do not refer to the laboratory frame. We demonstrate that symmetry-breaking is required in order to learn vertical migration in a meaningful way. Using reinforcement learning we analyse the emerging strategies for different sets of environmental cues that microorganisms are known to measure.
△ Less
Submitted 22 April, 2021;
originally announced April 2021.
-
Key parameters for droplet evaporation and mixing at the cloud edge
Authors:
J. Fries,
G. Sardina,
G. Svensson,
B. Mehlig
Abstract:
The distribution of liquid water in ice-free clouds determines their radiative properties, a significant source of uncertainty in weather and climate models. Evaporation and turbulent mixing cause a cloud to display large variations in droplet-number density, but quite small variations in droplet size [Beals et al. (2015)]. Yet direct numerical simulations of the joint effect of evaporation and mi…
▽ More
The distribution of liquid water in ice-free clouds determines their radiative properties, a significant source of uncertainty in weather and climate models. Evaporation and turbulent mixing cause a cloud to display large variations in droplet-number density, but quite small variations in droplet size [Beals et al. (2015)]. Yet direct numerical simulations of the joint effect of evaporation and mixing near the cloud edge predict quite different behaviors, and it remains an open question how to reconcile these results with the experimental findings. To infer the history of mixing and evaporation from observational snapshots of droplets in clouds is challenging because clouds are transient systems. We formulated a statistical model that provides a reliable description of the evaporation-mixing process as seen in direct numerical simulations, and allows to infer important aspects of the history of observed droplet populations, highlighting the key mechanisms at work, and explaining the differences between observations and simulations.
△ Less
Submitted 4 February, 2021;
originally announced February 2021.
-
Paths to caustic formation in turbulent aerosols
Authors:
Jan Meibohm,
Vikash Pandey,
Akshay Bhatnagar,
Kristian Gustavsson,
Dhrubaditya Mitra,
Prasad Perlekar,
B. Mehlig
Abstract:
The dynamics of small, yet heavy, identical particles in turbulence exhibits singularities, called caustics, that lead to large fluctuations in the spatial particle-number density, and in collision velocities. For large particle, inertia the fluid velocity at the particle position is essentially a white-noise signal and caustic formation is analogous to Kramers escape. Here we show that caustic fo…
▽ More
The dynamics of small, yet heavy, identical particles in turbulence exhibits singularities, called caustics, that lead to large fluctuations in the spatial particle-number density, and in collision velocities. For large particle, inertia the fluid velocity at the particle position is essentially a white-noise signal and caustic formation is analogous to Kramers escape. Here we show that caustic formation at small particle inertia is different. Caustics tend to form in the vicinity of particle trajectories that experience a specific history of fluid-velocity gradients, characterised by low vorticity and a violent strain exceeding a large threshold. We develop a theory that explains our findings in terms of an optimal path to caustic formation that is approached in the small inertia limit.
△ Less
Submitted 15 June, 2021; v1 submitted 15 December, 2020;
originally announced December 2020.
-
Planetesimals on eccentric orbits erode rapidly
Authors:
Lukas Cedenblad,
Noemi Schaffer,
Anders Johansen,
B. Mehlig,
Dhrubaditya Mitra
Abstract:
We investigate the possibility of erosion of planetesimals in a protoplanetary disk. We use theory and direct numerical simulations (Lattice Boltzmann Method) to calculate the erosion of large -- much larger than the mean-free-path of gas molecules -- bodies of different shapes in flows. We find that erosion follows a universal power-law in time, at intermediate times, independent of the Reynolds…
▽ More
We investigate the possibility of erosion of planetesimals in a protoplanetary disk. We use theory and direct numerical simulations (Lattice Boltzmann Method) to calculate the erosion of large -- much larger than the mean-free-path of gas molecules -- bodies of different shapes in flows. We find that erosion follows a universal power-law in time, at intermediate times, independent of the Reynolds number of the flow and the initial shape of the body. Consequently, we estimate that planetesimals in eccentric orbits, of even very small eccentricity, rapidly (in about hundred years) erodes away if the semi-major axis of their orbit lies in the inner disk -- less than about $10$ au. Even planetesimals in circular orbits erode away in approximately ten thousand years if the semi-major axis of their orbits are $\lessapprox 0.6$au.
△ Less
Submitted 30 September, 2021; v1 submitted 29 November, 2020;
originally announced November 2020.
-
Effect of particle inertia on the alignment of small ice crystals in turbulent clouds
Authors:
K. Gustavsson,
M. Z. Sheikh,
A. Naso,
A. Pumir,
B. Mehlig
Abstract:
Small non-spherical particles settling in a quiescent fluid tend to orient so that their broad side faces down, because this is a stable fixed point of their angular dynamics at small particle Reynolds number. Turbulence randomises the orientations to some extent, and this affects the reflection patterns of polarised light from turbulent clouds containing ice crystals. An overdamped theory predict…
▽ More
Small non-spherical particles settling in a quiescent fluid tend to orient so that their broad side faces down, because this is a stable fixed point of their angular dynamics at small particle Reynolds number. Turbulence randomises the orientations to some extent, and this affects the reflection patterns of polarised light from turbulent clouds containing ice crystals. An overdamped theory predicts that turbulence-induced fluctuations of the orientation are very small when the settling number Sv (a dimensionless measure of the settling speed) is large. At small Sv, by contrast, the overdamped theory predicts that turbulence randomises the orientations. This overdamped theory neglects the effect of particle inertia. Therefore we consider here how particle inertia affects the orientation of small crystals settling in turbulent air. We find that it can significantly increase the orientation variance, even when the Stokes number St (a dimensionless measure of particle inertia) is quite small. We identify different asymptotic parameter regimes where the tilt-angle variance is proportional to different inverse powers of Sv. We estimate parameter values for ice crystals in turbulent clouds and show that they cover several of the identified regimes. The theory predicts how the degree of alignment depends on particle size, shape and turbulence intensity, and that the strong horizontal alignment of small crystals is only possible when the turbulent energy dissipation is weak, of the order of $1\,$cm$^2$/s$^3$ or less.
△ Less
Submitted 4 February, 2021; v1 submitted 22 July, 2020;
originally announced July 2020.
-
Lord Kelvin's isotropic helicoid
Authors:
Darci Collins,
Rami J. Hamati,
Fabien Candelier,
Kristian Gustavsson,
Bernhard Mehlig,
Greg A. Voth
Abstract:
Nearly 150 years ago, Lord Kelvin proposed the isotropic helicoid, a particle with isotropic yet chiral interactions with a fluid, so that translation couples to rotation. An implementation of his design fabricated with a three-dimensional printer is found experimentally to have no detectable translation-rotation coupling, although the particle point-group symmetry allows this coupling. We explain…
▽ More
Nearly 150 years ago, Lord Kelvin proposed the isotropic helicoid, a particle with isotropic yet chiral interactions with a fluid, so that translation couples to rotation. An implementation of his design fabricated with a three-dimensional printer is found experimentally to have no detectable translation-rotation coupling, although the particle point-group symmetry allows this coupling. We explain these results by demonstrating that in Stokes flow, the chiral coupling of such isotropic helicoids made out of non-chiral vanes is due only to hydrodynamic interactions between these vanes. Therefore it is small. In summary, Kelvin's predicted isotropic helicoid exists, but only as a weak breaking of a symmetry of non-interacting vanes in Stokes flow.
△ Less
Submitted 21 July, 2021; v1 submitted 15 June, 2020;
originally announced June 2020.
-
Inertial torque on a small spheroid in a stationary uniform flow
Authors:
F. Jiang,
L. Zhao,
H. I. Andersson,
K. Gustavsson,
A. Pumir,
B. Mehlig
Abstract:
How anisotropic particles rotate and orient in a flow depends on the hydrodynamic torque they experience. Here we compute the torque acting on a small spheroid in a uniform flow by numerically solving the Navier-Stokes equations. Particle shape is varied from oblate (aspect ratio $λ= 1/6$) to prolate ($λ= 6$), and we consider low and moderate particle Reynolds numbers (${\rm Re} \le 50$). We demon…
▽ More
How anisotropic particles rotate and orient in a flow depends on the hydrodynamic torque they experience. Here we compute the torque acting on a small spheroid in a uniform flow by numerically solving the Navier-Stokes equations. Particle shape is varied from oblate (aspect ratio $λ= 1/6$) to prolate ($λ= 6$), and we consider low and moderate particle Reynolds numbers (${\rm Re} \le 50$). We demonstrate that the angular dependence of the torque, predicted theoretically for small particle Reynolds numbers remains qualitatively correct for Reynolds numbers up to ${\rm Re} \sim 10$. The amplitude of the torque, however, is smaller than the theoretical prediction, the more so as ${\rm Re}$ increases. For Re larger than $10$, the flow past oblate spheroids acquires a more complicated structure, resulting in systematic deviations from the theoretical predictions. Overall, our numerical results provide a justification of recent theories for the orientation statistics of ice-crystals settling in a turbulent flow.
△ Less
Submitted 6 February, 2021; v1 submitted 12 May, 2020;
originally announced May 2020.
-
Alignment statistics of rods with the Lagrangian stretching direction in a channel flow
Authors:
Z. Cui,
A. Dubey,
L. Zhao,
B. Mehlig
Abstract:
In homogeneous isotropic turbulence, slender rods are known to align with the Lagrangian stretching direction. However, how the degree of alignment depends on the aspect ratio of the rod is not understood. Moreover, many flows of practical interest are anisotropic and inhomogeneous. Here we study the alignment of rods with the Lagrangian stretching direction in a channel flow, which is approximate…
▽ More
In homogeneous isotropic turbulence, slender rods are known to align with the Lagrangian stretching direction. However, how the degree of alignment depends on the aspect ratio of the rod is not understood. Moreover, many flows of practical interest are anisotropic and inhomogeneous. Here we study the alignment of rods with the Lagrangian stretching direction in a channel flow, which is approximately homogeneous and isotropic near the center but inhomogeneous and anisotropic near the walls. Our main question is how the distribution of relative angles between a rod and the Lagrangian stretching direction depends on the aspect ratio of the rod and upon the distance of the rod from the channel wall. We find that the distribution exhibits two regimes: a plateau at small angles that corresponds to random uncorrelated motion, and power-law tails that describe large excursions. The variance of the relative angle is described by the width of the plateau. We find that slender rods near the channel center align better with the Lagrangian stretching direction, compared to those near the channel wall. These observations are explained in terms of simple statistical models based on Jeffery's equation, qualitatively near the channel center and quantitatively near the channel wall. Lastly we discuss the consequences of our results for the distribution of relative angles between the orientations of nearby rods (Zhao et al., Phys. Rev. Fluids, vol. 4, 2019, 054602).
△ Less
Submitted 3 February, 2020;
originally announced February 2020.
-
Intercellular communication induces glycolytic synchronization waves between individually oscillating cells
Authors:
Martin Mojica-Benavides,
David D. van Niekerk,
Mite Mijalkov,
Jacky L. Snoep,
Bernhard Mehlig,
Giovanni Volpe,
Mattias Goksör,
Caroline B. Adiels
Abstract:
Metabolic oscillations in single cells underlie the mechanisms behind cell synchronization and cell-cell communication. For example, glycolytic oscillations mediated by biochemical communication between cells may synchronize the pulsatile insulin secretion by pancreatic tissue, and a link between glycolytic synchronization anomalies and type-2 diabetes has been hypotesized. Cultures of yeast cells…
▽ More
Metabolic oscillations in single cells underlie the mechanisms behind cell synchronization and cell-cell communication. For example, glycolytic oscillations mediated by biochemical communication between cells may synchronize the pulsatile insulin secretion by pancreatic tissue, and a link between glycolytic synchronization anomalies and type-2 diabetes has been hypotesized. Cultures of yeast cells have provided an ideal model system to study synchronization and propagation waves of glycolytic oscillations in large populations. However, the mechanism by which synchronization occurs at individual cell-cell level and overcome local chemical concentrations and heterogenic biological clocks, is still an open question because of experimental limitations in sensitive and specific handling of single cells. Here, we show how the coupling of intercellular diffusion with the phase regulation of individual oscillating cells induce glycolytic synchronization waves. We directly measure the single-cell metabolic responses from yeast cells in a microfluidic environment and characterize a discretized cell-cell communication using graph theory. We corroborate our findings with simulations based on a kinetic detailed model for individual yeast cells. These findings can provide insight into the roles cellular synchronization play in biomedical applications, such as insulin secretion regulation at the cellular level.
△ Less
Submitted 26 March, 2020; v1 submitted 11 September, 2019;
originally announced September 2019.
-
Alignment of non-spherical active particles in chaotic flows
Authors:
M. Borgnino,
K. Gustavsson,
F. De Lillo,
G. Boffetta,
M. Cencini,
B. Mehlig
Abstract:
We study the orientation statistics of spheroidal, axisymmetric microswimmers, with shapes ranging from disks to rods, swimming in chaotic, moderately turbulent flows. Numerical simulations show that rod-like active particles preferentially align with the flow velocity. To explain the underlying mechanism we solve a statistical model via perturbation theory. We show that such alignment is caused b…
▽ More
We study the orientation statistics of spheroidal, axisymmetric microswimmers, with shapes ranging from disks to rods, swimming in chaotic, moderately turbulent flows. Numerical simulations show that rod-like active particles preferentially align with the flow velocity. To explain the underlying mechanism we solve a statistical model via perturbation theory. We show that such alignment is caused by correlations of fluid velocity and its gradients along particle paths combined with fore-aft symmetry breaking due to both swimming and particle nonsphericity. Remarkably, the discovered alignment is found to be a robust kinematical effect, independent of the underlying flow evolution. We discuss its possible relevance for aquatic ecology.
△ Less
Submitted 5 September, 2019;
originally announced September 2019.
-
Fractal catastrophes
Authors:
J. Meibohm,
K. Gustavsson,
J. Bec,
B. Mehlig
Abstract:
We analyse the spatial inhomogeneities ('spatial clustering') in the distribution of particles accelerated by a force that changes randomly in space and time. To quantify spatial clustering, the phase-space dynamics of the particles must be projected to configuration space. Folds of a smooth phase-space manifold give rise to catastrophes ('caustics') in this projection. When the inertial particle…
▽ More
We analyse the spatial inhomogeneities ('spatial clustering') in the distribution of particles accelerated by a force that changes randomly in space and time. To quantify spatial clustering, the phase-space dynamics of the particles must be projected to configuration space. Folds of a smooth phase-space manifold give rise to catastrophes ('caustics') in this projection. When the inertial particle dynamics is damped by friction, however, the phase-space manifold converges towards a fractal attractor. It is believed that caustics increase spatial clustering also in this case, but a quantitative theory is missing. We solve this problem by determining how projection affects the distribution of finite-time Lyapunov exponents. Applying our method in one spatial dimension we find that caustics arising from the projection of a dynamical fractal attractor ('fractal catastrophes') make a distinct and universal contribution to the distribution of spatial finite-time Lyapunov exponents. Our results explain a projection formula for the spatial fractal correlation dimension, and how a fluctuation relation for the distribution of finite-time Lyapunov exponents for white-in-time Gaussian force fields breaks upon projection. We explore the implications of our results for heavy particles in turbulence, and for wave propagation in random media.
△ Less
Submitted 8 October, 2019; v1 submitted 21 May, 2019;
originally announced May 2019.
-
Importance of fluid inertia for the orientation of spheroids settling in turbulent flow
Authors:
Muhammad Zubair Sheikh,
Kristian Gustavsson,
Diego Lopez,
Emmanuel Leveque,
Bernhard Mehlig,
Alain Pumir,
Aurore Naso
Abstract:
How non-spherical particles orient as they settle in a flow has important practical implications in a number of scientific and engineering problems. In a quiescent fluid, a slowly settling particle orients so that it settles with its broad side first. This is an effect of the torque due to convective inertia of the fluid set in motion by the settling particle, which maximises the drag experienced…
▽ More
How non-spherical particles orient as they settle in a flow has important practical implications in a number of scientific and engineering problems. In a quiescent fluid, a slowly settling particle orients so that it settles with its broad side first. This is an effect of the torque due to convective inertia of the fluid set in motion by the settling particle, which maximises the drag experienced by the particle. Turbulent flows tend to randomise the particle orientation. Recently the settling of non-spherical particles in turbulence was analysed neglecting the effect of convective fluid inertia, but taking into account the effect of the turbulent fluid-velocity gradients on the particle orientation. These studies reached the opposite conclusion, namely that a rod settles preferentially with its tip first, wheras a disk settles with its edge first, therefore minimizing the drag on the particle. Here, we consider both effects, the convective inertial torque as well as the torque due to fluctuating velocity gradients, and ask under which circumstances either one or the other dominate. To this end we estimate the ratio of the magnitudes of the two torques. Our estimates suggest that the fluid-inertia torque prevails in high-Reynolds number flows. In this case non-spherical particles are expected to settle with a maximal drag. But when the Reynolds number is small then the torque due to fluid-velocity gradients may dominate, causing the particle to settle with its broad side first.
△ Less
Submitted 27 August, 2019; v1 submitted 16 April, 2019;
originally announced April 2019.
-
Theory for the effect of fluid inertia on the orientation of a small particle settling in turbulence
Authors:
K. Gustavsson,
M. Z. Sheikh,
D. Lopez,
A. Naso,
A. Pumir,
B. Mehlig
Abstract:
Ice crystals settling through a turbulent cloud are rotated by turbulent velocity gradients. In the same way, turbulence affects the orientation of aggregates of organic matter settling in the ocean. In fact most solid particles encountered in Nature are not spherical, and their orientation affects their settling speed, as well as collision rates between particles. Therefore it is important to und…
▽ More
Ice crystals settling through a turbulent cloud are rotated by turbulent velocity gradients. In the same way, turbulence affects the orientation of aggregates of organic matter settling in the ocean. In fact most solid particles encountered in Nature are not spherical, and their orientation affects their settling speed, as well as collision rates between particles. Therefore it is important to understand the distribution of orientations of non-spherical particles settling in turbulence. Here we study the angular dynamics of small prolate spheroids settling in homogeneous isotropic turbulence. We consider a limit of the problem where the fluid torque due to convective inertia dominates, so that rods settle essentially horizontally. Turbulence causes the orientation of the settling particles to fluctuate, and we calculate their orientation distribution for prolate spheroids with arbitrary aspect ratios for large settling number Sv (a dimensionless measure of the settling speed), assuming small Stokes number St (a dimensionless measure of particle inertia). This overdamped theory predicts that the orientation distribution is very narrow at large Sv, with a variance proportional to ${\rm Sv}^{-4}$. By considering the role of particle inertia, we analyse the limitations of the overdamped theory, and determine its range of applicability. Our predictions are in excellent agreement with numerical simulations of simplified models of turbulent flows. Finally we contrast our results with those of an alternative theory predicting that the orientation variance scales as ${\rm Sv}^{-2}$ at large Sv.
△ Less
Submitted 31 March, 2019;
originally announced April 2019.
-
Heavy particles in a persistent random flow with traps
Authors:
J. Meibohm,
B. Mehlig
Abstract:
We study a one-dimensional model for heavy particles in a compressible fluid. The fluid-velocity field is modelled by a persistent Gaussian random function, and the particles are assumed to be weakly inertial. Since one-dimensional fluid-velocity fields are always compressible, the model exhibits spatial trapping regions where particles tend to accumulate. We determine the statistics of fluid-velo…
▽ More
We study a one-dimensional model for heavy particles in a compressible fluid. The fluid-velocity field is modelled by a persistent Gaussian random function, and the particles are assumed to be weakly inertial. Since one-dimensional fluid-velocity fields are always compressible, the model exhibits spatial trapping regions where particles tend to accumulate. We determine the statistics of fluid-velocity gradients in the vicinity of these traps and show how this allows to determine the spatial Lyapunov exponent and the rate of caustic formation. We compare our analytical results with numerical simulations of the model and explore the limits of validity of the theory. Finally, we discuss implications for higher-dimensional systems.
△ Less
Submitted 5 August, 2019; v1 submitted 12 February, 2019;
originally announced February 2019.
-
Collision fluctuations of lucky droplets with superdroplets
Authors:
Xiang-Yu Li,
Bernhard Mehlig,
Gunilla Svensson,
Axel Brandenburg,
Nils E. L. Haugen
Abstract:
It was previously shown that the superdroplet algorithm for modeling the collision-coalescence process can faithfully represent mean droplet growth in turbulent clouds. But an open question is how accurately the superdroplet algorithm accounts for fluctuations in the collisional aggregation process. Such fluctuations are particularly important in dilute suspensions. Even in the absence of turbulen…
▽ More
It was previously shown that the superdroplet algorithm for modeling the collision-coalescence process can faithfully represent mean droplet growth in turbulent clouds. But an open question is how accurately the superdroplet algorithm accounts for fluctuations in the collisional aggregation process. Such fluctuations are particularly important in dilute suspensions. Even in the absence of turbulence, Poisson fluctuations of collision times in dilute suspensions may result in substantial variations in the growth process, resulting in a broad distribution of growth times to reach a certain droplet size. We quantify the accuracy of the superdroplet algorithm in describing the fluctuating growth history of a larger droplet that settles under the effect of gravity in a quiescent fluid and collides with a dilute suspension of smaller droplets that were initially randomly distributed in space ('lucky droplet model'). We assess the effect of fluctuations upon the growth history of the lucky droplet and compute the distribution of cumulative collision times. The latter is shown to be sensitive enough to detect the subtle increase of fluctuations associated with collisions between multiple lucky droplets. The superdroplet algorithm incorporates fluctuations in two distinct ways: through the random spatial distribution of superdroplets and through the Monte Carlo collision algorithm involved. Using specifically designed numerical experiments, we show that both on their own give an accurate representation of fluctuations. We conclude that the superdroplet algorithm can faithfully represent fluctuations in the coagulation of droplets driven by gravity.
△ Less
Submitted 3 July, 2022; v1 submitted 17 October, 2018;
originally announced October 2018.
-
Relative velocities in bi-disperse turbulent aerosols: simulations and theory
Authors:
Akshay Bhatnagar,
K. Gustavsson,
B. Mehlig,
Dhrubaditya Mitra
Abstract:
We perform direct numerical simulations of a bi-disperse suspension of heavy spherical particles in forced, homogeneous, and isotropic three-dimensional turbulence. We compute the joint distribution of relative particle distances and longitudinal relative velocities between particles of different sizes, and compare the results with recent theoretical predictions [Meibohm et al. Phys. Rev. E 96 (20…
▽ More
We perform direct numerical simulations of a bi-disperse suspension of heavy spherical particles in forced, homogeneous, and isotropic three-dimensional turbulence. We compute the joint distribution of relative particle distances and longitudinal relative velocities between particles of different sizes, and compare the results with recent theoretical predictions [Meibohm et al. Phys. Rev. E 96 (2017) 061102] for the shape of this distribution. We also compute the moments of relative velocities as a function of particle separation, and compare with the theoretical predictions. We observe good agreement.
△ Less
Submitted 20 December, 2018; v1 submitted 27 September, 2018;
originally announced September 2018.
-
Stochastic unfolding of nanoconfined DNA: experiments, model and Bayesian analysis
Authors:
Jens Krog,
Mohammadreza Alizadehheidari,
Erik Werner,
Santosh Kumar Bikkarolla,
Jonas O. Tegenfeldt,
Bernhard Mehlig,
Michael A. Lomholt,
Fredrik Westerlund,
Tobias Ambjornsson
Abstract:
Nanochannels provide means for detailed experiments on the effect of confinement on biomacromolecules, such as DNA. We here introduce a model for the complete unfolding of DNA from the circular to linear configuration. Two main ingredients are the entropic unfolding force as well as the friction coefficient for the unfolding process, and we describe the associated dynamics by a non-linear Langevin…
▽ More
Nanochannels provide means for detailed experiments on the effect of confinement on biomacromolecules, such as DNA. We here introduce a model for the complete unfolding of DNA from the circular to linear configuration. Two main ingredients are the entropic unfolding force as well as the friction coefficient for the unfolding process, and we describe the associated dynamics by a non-linear Langevin equation. By analyzing experimental data where DNA molecules are photo-cut and unfolded inside a nanochannel, our model allows us to extract values for the unfolding force as well as the friction coefficient for the first time. In order to extract numerical values for these physical quantities, we employ a recently introduced Bayesian inference framework. We find that the determined unfolding force is in agreement with estimates from a simple Flory type argument. The estimated friction coefficient is in agreement with theoretical estimates for motion of a cylinder in a channel. We further validate the estimated friction constant by extracting this parameter from DNA's center-of-mass motion before and after unfolding, yielding decent agreement.
△ Less
Submitted 30 January, 2020; v1 submitted 8 August, 2018;
originally announced August 2018.
-
Condensational and collisional growth of cloud droplets in a turbulent environment
Authors:
Xiang-Yu Li,
Axel Brandenburg,
Gunilla Svensson,
Nils Haugen,
Bernhard Mehlig,
Igor Rogachevskii
Abstract:
We investigate the effect of turbulence on the combined condensational and collisional growth of cloud droplets by means of high resolution direct numerical simulations of turbulence and a superparticle approximation for droplet dynamics and collisions. The droplets are subject to turbulence as well as gravity, and their collision and coalescence efficiencies are taken to be unity. We solve the th…
▽ More
We investigate the effect of turbulence on the combined condensational and collisional growth of cloud droplets by means of high resolution direct numerical simulations of turbulence and a superparticle approximation for droplet dynamics and collisions. The droplets are subject to turbulence as well as gravity, and their collision and coalescence efficiencies are taken to be unity. We solve the thermodynamic equations governing temperature, water-vapor mixing ratio, and the resulting supersaturation fields together with the Navier-Stokes equation. We find that the droplet-size distribution broadens with increasing Reynolds number and/or mean energy dissipation rate. Turbulence affects the condensational growth directly through supersaturation fluctuations, and it influences collisional growth indirectly through condensation. Our simulations show for the first time that, in the absence of the mean updraft cooling, supersaturation fluctuation-induced broadening of droplet-size distributions enhances the collisional growth. This is contrary to classical (non-turbulent) condensational growth, which leads to a growing mean droplet size, but a narrower droplet-size distribution. Our findings, instead, show that condensational growth facilitates collisional growth by broadening the size distribution in the tails at an early stage of rain formation. With increasing Reynolds numbers, evaporation becomes stronger. This counteracts the broadening effect due to condensation at late stages of rain formation. Our conclusions are consistent with results of laboratory experiments and field observations, and show that supersaturation fluctuations are important for precipitation.
△ Less
Submitted 26 January, 2020; v1 submitted 31 July, 2018;
originally announced July 2018.
-
Fractal dimensions and trajectory crossings in correlated random walks
Authors:
A. Dubey,
J. Meibohm,
K. Gustavsson,
B. Mehlig
Abstract:
We study spatial clustering in a discrete, one-dimensional, stochastic, toy model of heavy particles in turbulence and calculate the spectrum of multifractal dimensions $D_q$ as functions of a dimensionless parameter, $α$, that plays the role of an inertia parameter. Using the fact that it suffices to consider the linearized dynamics of the model at small separations, we find that…
▽ More
We study spatial clustering in a discrete, one-dimensional, stochastic, toy model of heavy particles in turbulence and calculate the spectrum of multifractal dimensions $D_q$ as functions of a dimensionless parameter, $α$, that plays the role of an inertia parameter. Using the fact that it suffices to consider the linearized dynamics of the model at small separations, we find that $D_q =D_2/(q-1)$ for $q=2,3,\ldots$. The correlation dimension $D_2$ turns out to be a non-analytic function of the inertia parameter in this model. We calculate $D_2$ for small $α$ up to the next-to-leading order in the non-analytic term.
△ Less
Submitted 21 June, 2018;
originally announced June 2018.
-
Time-dependent lift and drag on a rigid body in a viscous steady linear flow
Authors:
Fabien Candelier,
Bernhard Mehlig,
Jacques Magnaudet
Abstract:
We compute the leading-order inertial corrections to the instantaneous force acting on a rigid body moving with a time-dependent slip velocity in a linear flow field, assuming that the variation of the undisturbed flow at the body scale is much larger than the slip velocity between the body and the fluid. Motivated by applications to turbulent particle-laden flows, we seek a formulation allowing t…
▽ More
We compute the leading-order inertial corrections to the instantaneous force acting on a rigid body moving with a time-dependent slip velocity in a linear flow field, assuming that the variation of the undisturbed flow at the body scale is much larger than the slip velocity between the body and the fluid. Motivated by applications to turbulent particle-laden flows, we seek a formulation allowing this force to be determined for an arbitrarily-shaped body moving in a general linear flow. We express the equations governing the flow disturbance in a non-orthogonal coordinate system moving with the undisturbed flow and solve the problem using matched asymptotic expansions. The use of the co-moving coordinates enables the leading-order inertial corrections to the force to be obtained at any time in an arbitrary linear flow field. We then specialize this approach to compute the time-dependent force components for a sphere moving in three canonical flows: solid body rotation, planar elongation, and uniform shear. We discuss the behaviour and physical origin of the different force components in the short-time and quasi-steady limits. Last, we illustrate the influence of time-dependent and quasi-steady inertial effects by examining the sedimentation of prolate and oblate spheroids in a pure shear flow.
△ Less
Submitted 14 June, 2018;
originally announced June 2018.
-
Distribution of label spacings for genome mapping in nanochannels
Authors:
D. Ödman,
E. Werner,
K. D. Dorfman,
C. R. Doering,
B. Mehlig
Abstract:
In genome mapping experiments, long DNA molecules are stretched by confining them to very narrow channels, so that the locations of sequence-specific fluorescent labels along the channel axis provide large-scale genomic information. It is difficult, however, to make the channels narrow enough so that the DNA molecule is fully stretched. In practice its conformations may form hairpins that change t…
▽ More
In genome mapping experiments, long DNA molecules are stretched by confining them to very narrow channels, so that the locations of sequence-specific fluorescent labels along the channel axis provide large-scale genomic information. It is difficult, however, to make the channels narrow enough so that the DNA molecule is fully stretched. In practice its conformations may form hairpins that change the spacings between internal segments of the DNA molecule, and thus the label locations along the channel axis. Here we describe a theory for the distribution of label spacings that explains the heavy tails observed in distributions of label spacings in genome mapping experiments.
△ Less
Submitted 13 June, 2018; v1 submitted 30 March, 2018;
originally announced March 2018.
-
Inertial drag on a sphere settling in a stratified fluid
Authors:
R. Mehaddi,
F. Candelier,
B. Mehlig
Abstract:
We compute the drag force on a sphere settling slowly in a quiescent, linearly stratified fluid. Stratification can significantly enhance the drag experienced by the settling particle. The magnitude of this effect depends on whether fluid-density transport around the settling particle is due to diffusion, to advection by the disturbance flow caused by the particle, or due to both. It therefore mat…
▽ More
We compute the drag force on a sphere settling slowly in a quiescent, linearly stratified fluid. Stratification can significantly enhance the drag experienced by the settling particle. The magnitude of this effect depends on whether fluid-density transport around the settling particle is due to diffusion, to advection by the disturbance flow caused by the particle, or due to both. It therefore matters how efficiently the fluid disturbance is convected away from the particle by fluid-inertial terms. When these terms dominate, the Oseen drag force must be recovered. We compute by perturbation theory how the Oseen drag is modified by diffusion and stratification. Our results are in good agreement with recent direct-numerical simulation studies of the problem at small Reynolds numbers and large (but not too large) Froude numbers.
△ Less
Submitted 28 February, 2018;
originally announced February 2018.
-
Spinning and tumbling of micron-sized triangles in a micro-channel shear flow
Authors:
J. Fries,
M. Vijay Kumar,
B. Mekonnen Mihiretie,
D. Hanstorp,
B. Mehlig
Abstract:
We report on measurements of the angular dynamics of micron-sized equilaterally triangular platelets suspended in a micro-channel shear flow. Our measurements confirm that such particles spin and tumble like a spheroid in a simple shear. Since the triangle has corners we can observe the spinning directly. In general the spinning frequency is different from the tumbling frequency, and the spinning…
▽ More
We report on measurements of the angular dynamics of micron-sized equilaterally triangular platelets suspended in a micro-channel shear flow. Our measurements confirm that such particles spin and tumble like a spheroid in a simple shear. Since the triangle has corners we can observe the spinning directly. In general the spinning frequency is different from the tumbling frequency, and the spinning is affected by tumbling. This gives rise to doubly-periodic angular dynamics.
△ Less
Submitted 22 December, 2017;
originally announced December 2017.
-
Effect of turbulence on collisional growth of cloud droplets
Authors:
Xiang-Yu Li,
A. Brandenburg,
G. Svensson,
N. E. L. Haugen,
B. Mehlig,
I. Rogachevskii
Abstract:
We investigate the effect of turbulence on the collisional growth of um-sized droplets through high- resolution numerical simulations with well resolved Kolmogorov scales, assuming a collision and coalescence efficiency of unity. The droplet dynamics and collisions are approximated using a superparticle approach. In the absence of gravity, we show that the time evolution of the shape of the drople…
▽ More
We investigate the effect of turbulence on the collisional growth of um-sized droplets through high- resolution numerical simulations with well resolved Kolmogorov scales, assuming a collision and coalescence efficiency of unity. The droplet dynamics and collisions are approximated using a superparticle approach. In the absence of gravity, we show that the time evolution of the shape of the droplet-size distribution due to turbulence-induced collisions depends strongly on the turbulent energy-dissipation rate, but only weakly on the Reynolds number. This can be explained through the energy dissipation rate dependence of the mean collision rate described by the Saffman-Turner collision model. Consistent with the Saffman-Turner collision model and its extensions, the collision rate increases as the square root of the energy dissipation rate even when coalescence is invoked. The size distribution exhibits power law behavior with a slope of -3.7 between a maximum at approximately 10 um up to about 40 um. When gravity is invoked, turbulence is found to dominate the time evolution of an initially monodisperse droplet distribution at early times. At later times, however, gravity takes over and dominates the collisional growth. We find that the formation of large droplets is very sensitive to the turbulent energy dissipation rate. This is due to the fact that turbulence enhances the collisional growth between similar sized droplets at the early stage of raindrop formation. The mean collision rate grows exponentially, which is consistent with the theoretical prediction of the continuous collisional growth even when turbulence-generated collisions are invoked. This consistency only reflects the mean effect of turbulence on collisional growth.
△ Less
Submitted 15 August, 2018; v1 submitted 27 November, 2017;
originally announced November 2017.
-
Passive directors in turbulence
Authors:
L. Zhao,
K. Gustavsson,
R. Ni,
S. Kramel,
G. Voth,
H. I. Andersson,
B. Mehlig
Abstract:
In experiments and numerical simulations we measured angles between the symmetry axes of small spheroids advected in turbulence ("passive directors"). Since turbulent strains tend to align nearby spheroids, one might think that their relative angles are quite small. We show that this intuition fails in general because angles between the symmetry axes of nearby particles are anomalously large. We i…
▽ More
In experiments and numerical simulations we measured angles between the symmetry axes of small spheroids advected in turbulence ("passive directors"). Since turbulent strains tend to align nearby spheroids, one might think that their relative angles are quite small. We show that this intuition fails in general because angles between the symmetry axes of nearby particles are anomalously large. We identify two mechanisms that cause this phenomenon. First, the dynamics evolves to a fractal attractor despite the fact that the fluid velocity is spatially smooth at small scales. Second, this fractal forms steps akin to scar lines observed in the director patterns for random or chaotic two-dimensional maps.
△ Less
Submitted 12 February, 2019; v1 submitted 19 July, 2017;
originally announced July 2017.
-
Emergent scar lines in chaotic advection of passive directors
Authors:
Bardia Hejazi,
Bernhard Mehlig,
Greg A. Voth
Abstract:
We examine the spatial field of orientations of slender fibers that are advected by a two-dimensional fluid flow. The orientation field of these passive directors are important in a wide range of industrial and geophysical flows. We introduce emergent scar lines as the dominant coherent structures in the orientation field of passive directors in chaotic flows. Previous work has identified the exis…
▽ More
We examine the spatial field of orientations of slender fibers that are advected by a two-dimensional fluid flow. The orientation field of these passive directors are important in a wide range of industrial and geophysical flows. We introduce emergent scar lines as the dominant coherent structures in the orientation field of passive directors in chaotic flows. Previous work has identified the existence of scar lines where the orientation rotates by π over short distances, but the lines that were identified disappeared as time progressed. As a result, earlier work focused on topological singularities in the orientation field which we find to play a negligible role at long times. We use the standard map as a simple time-periodic two-dimensional (2D) flow that produces Lagrangian chaos. This class of flows produce persistent patterns in passive scalar advection, and we find that a different kind of persistent pattern develops in the passive director orientation field. We identify the mechanism by which emergent scar lines grow to dominate these patterns at long times in complex flows. Emergent scar lines form where the recent stretching of the fluid element is perpendicular to earlier stretching. Thus these scar lines can be labeled by their age, defined as the time since their stretching reached a maximum.
△ Less
Submitted 4 December, 2017; v1 submitted 22 June, 2017;
originally announced June 2017.
-
Intrinsic viscosity of a suspension of weakly Brownian ellipsoids in shear
Authors:
G. Almondo,
J. Einarsson,
J. R. Angilella,
B. Mehlig
Abstract:
We analyze the angular dynamics of triaxial ellipsoids in a shear flow subject to weak thermal noise. By numerically integrating an overdamped angular Langevin equation, we find the steady angular probability distribution for a range of triaxial particle shapes. From this distribution we compute the intrinsic viscosity of a dilute suspension of triaxial particles. We determine how the viscosity de…
▽ More
We analyze the angular dynamics of triaxial ellipsoids in a shear flow subject to weak thermal noise. By numerically integrating an overdamped angular Langevin equation, we find the steady angular probability distribution for a range of triaxial particle shapes. From this distribution we compute the intrinsic viscosity of a dilute suspension of triaxial particles. We determine how the viscosity depends on particle shape in the limit of weak thermal noise. While the deterministic angular dynamics depends very sensitively on particle shape, we find that the shape dependence of the intrinsic viscosity is weaker, in general, and that suspensions of rod-like particles are the most sensitive to breaking of axisymmetry. The intrinsic viscosity of a dilute suspension of triaxial particles is smaller than that of a suspension of axisymmetric particles with the same volume, and the same ratio of major to minor axis lengths.
△ Less
Submitted 13 June, 2018; v1 submitted 19 May, 2017;
originally announced May 2017.
-
One-parameter scaling theory for DNA extension in a nanochannel
Authors:
E. Werner,
G. K. Cheong,
D. Gupta,
K. D. Dorfman,
B. Mehlig
Abstract:
Experiments measuring DNA extension in nanochannels are at odds with even the most basic predictions of current scaling arguments for the conformations of confined semiflexible polymers such as DNA. We show that a theory based on a weakly self-avoiding, one-dimensional "telegraph" process collapses experimental data and simulation results onto a single master curve throughout the experimentally re…
▽ More
Experiments measuring DNA extension in nanochannels are at odds with even the most basic predictions of current scaling arguments for the conformations of confined semiflexible polymers such as DNA. We show that a theory based on a weakly self-avoiding, one-dimensional "telegraph" process collapses experimental data and simulation results onto a single master curve throughout the experimentally relevant region of parameter space and explains the mechanisms at play.
△ Less
Submitted 13 December, 2017; v1 submitted 12 May, 2017;
originally announced May 2017.
-
Optical manipulation for studies of collisional dynamics of micron-sized droplets under gravity
Authors:
M. Ivanov,
K. Chang,
I. Galinskiy,
B. Mehlig,
D. Hanstorp
Abstract:
A new experimental technique for creating and imaging collisions of micron-sized droplets settling under gravity is presented. A pair of glycerol droplets is suspended in air by means of two optical traps. The droplet relative velocities are determined by the droplet sizes. The impact parameter is precisely controlled by positioning the droplets using the two optical traps. The droplets are releas…
▽ More
A new experimental technique for creating and imaging collisions of micron-sized droplets settling under gravity is presented. A pair of glycerol droplets is suspended in air by means of two optical traps. The droplet relative velocities are determined by the droplet sizes. The impact parameter is precisely controlled by positioning the droplets using the two optical traps. The droplets are released by turning off the trapping light using electro-optical modulators. The motion of the sedimenting droplets is then captured by two synchronized high-speed cameras, at a frame rate of up to 63 kHz. The method allows the direct imaging of the collision of droplets without the influence of the optical confinement imposed by the trapping force. The method will facilitate efficient studies of the microphysics of neutral, as well as charged, liquid droplets and their interactions with light, electric field and thermodynamic environment, such as temperature or vapor concentration.
△ Less
Submitted 31 March, 2017;
originally announced March 2017.
-
Relative velocities in bidisperse turbulent suspensions
Authors:
J. Meibohm,
L. Pistone,
K. Gustavsson,
B. Mehlig
Abstract:
We investigate the distribution of relative velocities between small heavy particles of different sizes in turbulence by analysing a statistical model for bidisperse turbulent suspensions, containing particles with two different Stokes numbers. This number, ${\rm St}$, is a measure of particle inertia which in turn depends on particle size. When the Stokes numbers are similar, the distribution exh…
▽ More
We investigate the distribution of relative velocities between small heavy particles of different sizes in turbulence by analysing a statistical model for bidisperse turbulent suspensions, containing particles with two different Stokes numbers. This number, ${\rm St}$, is a measure of particle inertia which in turn depends on particle size. When the Stokes numbers are similar, the distribution exhibits power-law tails, just as in the case of equal ${\rm St}$. The power-law exponent is a non-analytic function of the mean Stokes number $\overline{\rm St}$, so that the exponent cannot be calculated in perturbation theory around the advective limit. When the Stokes-number difference is larger, the power law disappears, but the tails of the distribution still dominate the relative-velocity moments, if $\overline{\rm St}$ is large enough.
△ Less
Submitted 18 July, 2017; v1 submitted 5 March, 2017;
originally announced March 2017.
-
Statistical model for the orientation of non-spherical particles settling in turbulence
Authors:
K. Gustavsson,
J. Jucha,
A. Naso,
E. Lévêque,
A. Pumir,
B. Mehlig
Abstract:
The orientation of small anisotropic particles settling in a turbulent fluid determines some essential properties of the suspension. We show that the orientation distribution of small heavy spheroids settling through turbulence can be accurately predicted by a simple Gaussian statistical model that takes into account particle inertia and provides a quantitative understanding of the orientation dis…
▽ More
The orientation of small anisotropic particles settling in a turbulent fluid determines some essential properties of the suspension. We show that the orientation distribution of small heavy spheroids settling through turbulence can be accurately predicted by a simple Gaussian statistical model that takes into account particle inertia and provides a quantitative understanding of the orientation distribution on the problem parameters when fluid inertia is negligible. Our results open the way to a parameterisation of the distribution of ice-crystals in clouds, and potentially leads to an improved understanding of radiation reflection, or particle aggregation through collisions in clouds.
△ Less
Submitted 27 July, 2017; v1 submitted 24 February, 2017;
originally announced February 2017.
-
Hairpins in the conformations of a confined polymer
Authors:
E. Werner,
A. Jain,
A. Muralidhar,
K. Frykholm,
T. St Clere Smithe,
J. Fritzsche,
F. Westerlund,
K. D. Dorfman,
B. Mehlig
Abstract:
If a semiflexible polymer confined to a narrow channel bends around by 180 degrees, the polymer is said to exhibit a hairpin. The equilibrium extension statistics of the confined polymer are well understood when hairpins are vanishingly rare or when they are plentiful. Here we analyze the extension statistics in the intermediate situation via experiments with DNA coated by the protein RecA, which…
▽ More
If a semiflexible polymer confined to a narrow channel bends around by 180 degrees, the polymer is said to exhibit a hairpin. The equilibrium extension statistics of the confined polymer are well understood when hairpins are vanishingly rare or when they are plentiful. Here we analyze the extension statistics in the intermediate situation via experiments with DNA coated by the protein RecA, which enhances the stiffness of the DNA molecule by approximately one order of magnitude. We find that the extension distribution is highly non-Gaussian, in good agreement with Monte Carlo simulations of confined discrete wormlike chains. We develop a simple model that qualitatively explains the form of the extension distribution. The model shows that the tail of the distribution at short extensions is determined by conformations with one hairpin.
△ Less
Submitted 30 January, 2018; v1 submitted 17 November, 2016;
originally announced November 2016.
-
Spherical particle sedimenting in weakly viscoelastic shear flow
Authors:
J. Einarsson,
B. Mehlig
Abstract:
We consider the dynamics of a small spherical particle driven through an unbounded viscoelastic shear flow by an external force. We give analytical solutions to both the mobility problem (velocity of forced particle) and the resistance problem (force on fixed particle), valid to second order in the dimensionless Deborah and Weissenberg numbers, which represent the elastic relaxation time of the fl…
▽ More
We consider the dynamics of a small spherical particle driven through an unbounded viscoelastic shear flow by an external force. We give analytical solutions to both the mobility problem (velocity of forced particle) and the resistance problem (force on fixed particle), valid to second order in the dimensionless Deborah and Weissenberg numbers, which represent the elastic relaxation time of the fluid relative to the rate of translation and the imposed shear rate. We find a shear-induced lift at $O({\rm Wi})$, a modified drag at $O({\rm De}^2)$ and $O({\rm Wi}^2)$, and a second lift that is orthogonal to the first, at $O({\rm Wi}^2)$. The relative importance of these effects depends strongly on the orientation of the forcing relative to the shear. We discuss how these forces affect the terminal settling velocity in an inclined shear flow. We also describe a new basis set of symmetric Cartesian tensors, and demonstrate how they enable general tensorial perturbation calculations such as the present theory. In particular this scheme allows us to write down a solution to the inhomogenous Stokes equations, required by the perturbation expansion, by a sequence of algebraic manipulations well suited to computer implementation.
△ Less
Submitted 11 May, 2017; v1 submitted 18 September, 2016;
originally announced September 2016.
-
Angular dynamics of small crystals in viscous flow
Authors:
J. Fries,
J. Einarsson,
B. Mehlig
Abstract:
The angular dynamics of a very small ellipsoidal particle in a viscous flow decouples from its translational dynamics, and the particle angular velocity is given by Jeffery's theory. It is known that cuboid particles share these properties. In the literature a special case is most frequently discussed, namely that of axisymmetric particles with a continuous rotation symmetry. Here we compute the a…
▽ More
The angular dynamics of a very small ellipsoidal particle in a viscous flow decouples from its translational dynamics, and the particle angular velocity is given by Jeffery's theory. It is known that cuboid particles share these properties. In the literature a special case is most frequently discussed, namely that of axisymmetric particles with a continuous rotation symmetry. Here we compute the angular dynamics of crystals that possess a discrete rotation symmetry and certain mirror symmetries, but that do not have a continuous rotation symmetry. We give examples of such particles that nevertheless obey Jeffery's theory. But there are other examples where the angular dynamics is determined by a more general equation of motion.
△ Less
Submitted 23 January, 2017; v1 submitted 25 July, 2016;
originally announced July 2016.
-
Angular dynamics of a small particle in turbulence
Authors:
F. Candelier,
J. Einarsson,
B. Mehlig
Abstract:
We compute the angular dynamics of a neutrally buoyant nearly spherical particle immersed in an unsteady fluid. We assume that the particle is small, that its translational slip velocity is negligible, and that unsteady and convective inertia are small perturbations. We derive an approximation for the torque on the particle that determines the first inertial corrections to Jeffery's equation. Thes…
▽ More
We compute the angular dynamics of a neutrally buoyant nearly spherical particle immersed in an unsteady fluid. We assume that the particle is small, that its translational slip velocity is negligible, and that unsteady and convective inertia are small perturbations. We derive an approximation for the torque on the particle that determines the first inertial corrections to Jeffery's equation. These corrections arise as a consequence of local vortex stretching, and can be substantial in turbulence where local vortex stretching is strong and closely linked to the irreversibility of turbulence.
△ Less
Submitted 2 October, 2016; v1 submitted 23 July, 2016;
originally announced July 2016.