Source author record

K. Gustavsson

K. Gustavsson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

Fluid-inertia torques from particle-shape symmetry

Numerical simulation of particle motion in fluids at low particle Reynolds numbers is often based on empirical force and torque models obtained by fitting force and torque from ab-initio computations for simple particle shapes such as spheres, spheroids, or cylindrical disks and fibres. To do the same for more complex particles shapes, one needs to first know how particle shape constrains the dependence of force and torque on flow velocity, its gradient, and on particle orientation. Here we use symmetry analysis and perturbation theory to determine the form of the hydrodynamic torque on a particle settling in a quiescent fluid at low but non-zero particle Reynolds numbers, for particle shapes with different point-group symmetries. The symmetry conclusions are verified by comparing with explicit calculations for nearly spherical particles.

preprint2022arXiv

Alignment of elongated swimmers in a laminar and turbulent Kolmogorov flow

Many aquatic microorganisms are able to swim. In natural environments they typically do so in the presence of flows. In recent years it has been shown that the interplay of swimming and flows can give rise to interesting and biologically relevant phenomena, such as accumulation of microorganisms in specific flow regions and local alignment with the flow properties. Here, we consider a mechanical model for elongated microswimmers in a Kolmogorov flow, a prototypic shear flow, both in steady and in turbulent conditions. By means of direct numerical simulations, supported by analytical calculation in a simplified stochastic setting, we find that the alignment of the swimming direction with the local velocity is a general phenomenon. We also explore how the accumulation of microorganisms, typically observed in steady flows, is modified by the presence of unsteady fluctuations.

preprint2022arXiv

Bifurcations in droplet collisions

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

preprint2021arXiv

Inertial torque on a small spheroid in a stationary uniform flow

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.

preprint2019arXiv

Fractal catastrophes

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.

preprint2019arXiv

Theory for the effect of fluid inertia on the orientation of a small particle settling in turbulence

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.

preprint2016arXiv

Inertial-particle accelerations in turbulence: a Lagrangian closure

The distribution of particle accelerations in turbulence is intermittent, with non-Gaussian tails that are quite different for light and heavy particles. In this article we analyse a closure scheme for the acceleration fluctuations of light and heavy inertial particles in turbulence, formulated in terms of Lagrangian correlation functions of fluid tracers. We compute the variance and the flatness of inertial particle accelerations and we discuss their dependency on the Stokes number. The closure incorporates effects induced by the Lagrangian correlations along the trajectories of fluid tracers, and its predictions agree well with results of direct numerical simulations of inertial particles in turbulence, provided that the effects induced by the inertial preferential sampling of heavy/light particles outside/inside vortices are negligible. In particular, the scheme predicts the correct functional behaviour of the acceleration variance, as a function of Stokes, as well as the presence of a minimum/maximum for the flatness of the acceleration of heavy/light particles, in good qualitative agreement with numerical data. We also show that the closure works well when applied to the Lagrangian evolution of particles using a stochastic surrogate for the underlying Eulerian velocity field. Our results support the conclusion that there exist important contributions to the statistics of the acceleration of inertial particles independent of the preferential sampling. For heavy particles we observe deviations between the predictions of the closure scheme and direct numerical simulations, at Stokes numbers of order unity. For light particles the deviation occurs for larger Stokes numbers.

preprint2016arXiv

Preferential sampling and small-scale clustering of gyrotactic microswimmers in turbulence

Recent studies show that spherical motile micro-organisms in turbulence subject to gravitational torques gather in down-welling regions of the turbulent flow. By analysing a statistical model we analytically compute how shape affects the dynamics, preferential sampling, and small-scale spatial clustering. We find that oblong organisms may spend more time in up-welling regions of the flow, and that all organisms are biased to regions of positive fluid-velocity gradients in the upward direction. We analyse small-scale spatial clustering and find that oblong particles may either cluster more or less than spherical ones, depending on the strength of the gravitational torques.

preprint2016arXiv

Statistical model for collisions and recollisions of inertial particles in mixing flows

Finding a quantitative description of the rate of collisions between small particles suspended in mixing flows is a long-standing problem. Here we investigate the validity of a parameterisation of the collision rate for identical particles subject to Stokes force, based on results for relative velocities of heavy particles that were recently obtained within a statistical model for the dynamics of turbulent aerosols. This model represents the turbulent velocity fluctuations by Gaussian random functions. We find that the parameterisation gives quantitatively good results in the limit where the \lq ghost-particle approximation' applies. The collision rate is a sum of two contributions due to \lq caustics' and to \lq clustering'. Within the statistical model we compare the relative importance of these two collision mechanisms. The caustic formation rate is high when the particle inertia becomes large, and we find that caustics dominate the collision rate as soon as they form frequently. We compare the magnitude of the caustic contribution to the collision rate to the formation rate of caustics.

preprint2016arXiv

Statistical models for spatial patterns of heavy particles in turbulence

The dynamics of heavy particles suspended in turbulent flows is of fundamental importance for a wide range of questions in astrophysics, atmospheric physics, oceanography, and technology. Laboratory experiments and numerical simulations have demonstrated that heavy particles respond in intricate ways to turbulent fluctuations of the carrying fluid: non-interacting particles may cluster together and form spatial patterns even though the fluid is incompressible, and the relative speeds of nearby particles can fluctuate strongly. Both phenomena depend sensitively on the parameters of the system. This parameter dependence is difficult to model from first principles since turbulence plays an essential role. Laboratory experiments are also very difficult, precisely since they must refer to a turbulent environment. But in recent years it has become clear that important aspects of the dynamics of heavy particles in turbulence can be understood in terms of statistical models where the turbulent fluctuations are approximated by Gaussian random functions with appropriate correlation functions. In this review we summarise how such statistical-model calculations have led to a detailed understanding of the factors that determine heavy-particle dynamics in turbulence. We concentrate on spatial clustering of heavy particles in turbulence. This is an important question because spatial clustering affects the collision rate between the particles and thus the long-term fate of the system.

preprint2014arXiv

Clustering of particles falling in a turbulent flow

Spatial clustering of identical particles falling through a turbulent flow enhances the collision rate between the falling particles, an important problem in aerosol science. We analyse this problem using perturbation theory in a dimensionless parameter, the so-called Kubo number. This allows us to derive an analytical theory quantifying the spatial clustering. We find that clustering of small particles in incompressible random velocity fields may be reduced or enhanced by the effect of gravity (depending on the Stokes number of the particles) and may be strongly anisotropic.

preprint2014arXiv

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 in vortex tubes.

preprint2013arXiv

Inertial-particle dynamics in turbulent flows: caustics, concentration fluctuations, and random uncorrelated motion

We discuss the relation between three recent approaches of describing the dynamics and the spatial distribution of particles suspended in turbulent flows: phase-space singularities in the inertial particle dynamics (caustics), real-space singularities of the deformation tensor, and random uncorrelated motion. We discuss how the phase- and real-space singularities are related. Their formation is well understood in terms of a local theory. We discuss implications for random uncorrelated motion. Our results are supported by results of direct numerical simulations of inertial particles in model flows.

preprint2013arXiv

Lyapunov exponents for particles advected in compressible random velocity fields at small and large Kubo numbers

We calculate the Lyapunov exponents describing spatial clustering of particles advected in one- and two-dimensional random velocity fields at finite Kubo numbers Ku (a dimensionless parameter characterising the correlation time of the velocity field). In one dimension we obtain accurate results up to Ku ~ 1 by resummation of a perturbation expansion in Ku. At large Kubo numbers we compute the Lyapunov exponent by taking into account the fact that the particles follow the minima of the potential function corresponding to the velocity field. The Lyapunov exponent is always negative. In two spatial dimensions the sign of the maximal Lyapunov exponent λ_1 may change, depending upon the degree of compressibility of the flow and the Kubo number. For small Kubo numbers we compute the first four non-vanishing terms in the small-Ku expansion of the Lyapunov exponents. By resumming these expansions we obtain a precise estimate of the location of the path-coalescence transition (where λ_1 changes sign) for Kubo numbers up to approximately Ku = 0.5. For large Kubo numbers we estimate the Lyapunov exponents for a partially compressible velocity field by assuming that the particles sample those stagnation points of the velocity field that have a negative real part of the maximal eigenvalue of the matrix of flow-velocity gradients.

preprint2013arXiv

Relative velocities of inertial particles in turbulent aerosols

We compute the joint distribution of relative velocities and separations of identical inertial particles suspended in randomly mixing and turbulent flows. Our results are obtained by matching asymptotic forms of the distribution. The method takes into account spatial clustering of the suspended particles as well as singularities in their motion (so-called 'caustics'). It thus takes proper account of the fractal properties of phase space and the distribution is characterised in terms of the corresponding phase-space fractal dimension D_2. The method clearly exhibits universal aspects of the distribution (independent of the statistical properties of the flow): at small particle separations R and not too large radial relative speeds |V_R|, the distribution of radial relative velocities exhibits a universal power-law form ρ(V_R,R) \sim |V_R|^{D_2-d-1} provided that D_2 < d+1 (d is the spatial dimension) and that the Stokes number St is large enough for caustics to form. The range in V_R over which this power law is valid depends on R, on the Stokes number, and upon the nature of the flow. Our results are in good agreement with results of computer simulations of the dynamics of particles suspended in random velocity fields with finite correlation times. In the white-noise limit the results are consistent with those of [Gustavsson and Mehlig, Phys. Rev. E84 (2011) 045304].

preprint2013arXiv

Tumbling of small axisymmetric particles in random and turbulent flows

We analyse the tumbling of small non-spherical, 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 [Parsa et al., Phys. Rev. Lett. 109 (2012) 134501]. 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.

preprint2012arXiv

Distribution of velocity gradients and rate of caustic formation in turbulent aerosols at finite Kubo numbers

In a one-dimensional model for a turbulent aerosol (inertial particles suspended in a random flow) we compute the distributions of particle-velocity gradients and the rate of caustic formation at finite but small Kubo numbers Ku, for arbitrary Stokes numbers St. Our results are consistent with those obtained earlier in the limit of small Ku and and large St, such that Ku^2 St remains constant. We show how finite-time correlations and non-ergodic effects influence the inertial-particle dynamics at finite but small Kubo numbers.

preprint2010arXiv

Clustering of exponentially separating trajectories

It might be expected that trajectories for a dynamical system which has no negative Lyapunov exponent (implying exponential growth of small separations will not cluster together. However, clustering can occur such that the density $ρ(Δx)$ of trajectories within distance $Δx$ of a reference trajectory has a power-law divergence, so that $ρ(Δx)\sim Δx^{-β}$ when $Δx$ is sufficiently small, for some $0<β<1$. We demonstrate this effect using a random map in one dimension. We find no evidence for this effect in the chaotic logistic map, and argue that the effect is harder to observe in deterministic maps.