Source author record

B. Mehlig

B. Mehlig 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

45works
14topics
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

45 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

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

preprint2022arXiv

Unsteady and inertial dynamics of an active particle in a fluid

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.

preprint2021arXiv

Collisions of micron-sized, charged water droplets in still air

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.

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.

preprint2020arXiv

Alignment statistics of rods with the Lagrangian stretching direction in a channel flow

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).

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

Angular dynamics of a small particle in turbulence

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.

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

Nanoconfined circular and linear DNA - equilibrium conformations and unfolding kinetics

Studies of circular DNA confined to nanofluidic channels are relevant both from a fundamental polymer-physics perspective and due to the importance of circular DNA molecules in vivo. We here observe the unfolding of DNA from the circular to linear configuration as a light-induced double strand break occurs, characterize the dynamics, and compare the equilibrium conformational statistics of linear and circular configurations. This is important because it allows us to determine to which extent existing statistical theories describe the extension of confined circular DNA. We find that the ratio of the extensions of confined linear and circular DNA configurations increases as the buffer concentration decreases. The experimental results fall between theoretical predictions for the extended de Gennes regime at weaker confinement and the Odijk regime at stronger confinement. We show that it is possible to directly distinguish between circular and linear DNA molecules by measuring the emission intensity from the DNA. Finally, we determine the rate of unfolding and show that this rate is larger for more confined DNA, possibly reflecting the corresponding larger difference in entropy between the circular and linear configurations.

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.

preprint2015arXiv

A model for melting of confined DNA

When DNA molecules are heated they denature. This occurs locally so that loops of molten single DNA strands form, connected by intact double-stranded DNA pieces. The properties of this "melting" transition have been intensively investigated. Recently there has been a surge of interest in this question, caused by experiments determining the properties of partially bound DNA confined to nanochannels. But how does such confinement affect the melting transition? To answer this question we introduce, and solve a model predicting how confinement affects the melting transition for a simple model system by first disregarding the effect of self-avoidance. We find that the transition is smoother for narrower channels. By means of Monte-Carlo simulations we then show that a model incorporating self-avoidance shows qualitatively the same behaviour and that the effect of confinement is stronger than in the ideal case.

preprint2015arXiv

Effect of weak fluid inertia upon Jeffery orbits

We consider the rotation of small neutrally buoyant axisymmetric particles in a viscous steady shear flow. When inertial effects are negligible the problem exhibits infinitely many periodic solutions, the "Jeffery orbits". We compute how inertial effects lift their degeneracy by perturbatively solving the coupled particle-flow equations. We obtain an equation of motion valid at small shear Reynolds numbers, for spheroidal particles with arbitrary aspect ratios. We analyse how the linear stability of the \lq log-rolling\rq{} orbit depends on particle shape and find it to be unstable for prolate spheroids. This resolves a puzzle in the interpretation of direct numerical simulations of the problem. In general both unsteady and non-linear terms in the Navier-Stokes equations are important.

preprint2015arXiv

Entrainment of heterogeneous glycolytic oscillations in single cells

Cell signaling, gene expression, and metabolism are affected by cell-cell heterogeneity and random changes in the environment. The effects of such fluctuations on cell signaling and gene expression have recently been studied intensively using single-cell experiments. In metabolism heterogeneity may be particularly important because it may affect synchronisation of metabolic oscillations, an important example of cell-cell communication. This synchronisation is notoriously difficult to describe theoretically as the example of glycolytic oscillations shows: neither is the mechanism of glycolytic synchronisation understood nor the role of cell-cell heterogeneity. To pin down the mechanism and to assess its robustness and universality we have experimentally investigated the entrainment of glycolytic oscillations in individual yeast cells by periodic external perturbations. We find that oscillatory cells synchronise through phase shifts and that the mechanism is insensitive to cell heterogeneity (robustness) and similar for different types of external perturbations (universality).

preprint2015arXiv

Extension of nano-confined DNA: quantitative comparison between experiment and theory

The extension of DNA confined to nanochannels has been studied intensively and in detail. Yet quantitative comparisons between experiments and model calculations are difficult because most theoretical predictions involve undetermined prefactors, and because the model parameters (contour length, Kuhn length, effective width) are difficult to compute reliably, leading to substantial uncertainties. Here we use a recent asymptotically exact theory for the DNA extension in the "extended de Gennes regime" that allows us to compare experimental results with theory. For this purpose we performed new experiments, measuring the mean DNA extension and its standard deviation while varying the channel geometry, dye intercalation ratio, and ionic buffer strength. The experimental results agree very well with theory at high ionic strengths, indicating that the model parameters are reliable. At low ionic strengths the agreement is less good. We discuss possible reasons. Our approach allows, in principle, to measure the Kuhn length and effective width of a single DNA molecule and more generally of semiflexible polymers in solution.

preprint2015arXiv

Numerical analysis of the angular motion of a neutrally buoyant spheroid in shear flow at small Reynolds numbers

We numerically analyse the rotation of a neutrally buoyant spheroid in a shear flow at small shear Reynolds number. Using direct numerical stability analysis of the coupled nonlinear particle-flow problem we compute the linear stability of the log-rolling orbit at small shear Reynolds number, ${\rm Re}_a$. As ${\rm Re}_a \to 0$ and as the box size of the system tends to infinity we find good agreement between the numerical results and earlier analytical predictions valid to linear order in ${\rm Re}_a$ for the case of an unbounded shear. The numerical stability analysis indicates that there are substantial finite-size corrections to the analytical results obtained for the unbounded system. We also compare the analytical results to results of lattice-Boltzmann simulations to analyse the stability of the tumbling orbit at shear Reynolds numbers of order unity. Theory for an unbounded system at infinitesimal shear Reynolds number predicts a bifurcation of the tumbling orbit at aspect ratio $λ_{\rm c} \approx 0.137$ below which tumbling is stable (as well as log rolling). The simulation results show a bifurcation line in the $λ$-${\rm Re}_a$ plane that reaches $λ\approx0.1275$ at the smallest shear Reynolds number (${\rm Re}_a=1$) at which we could simulate with the lattice-Boltzmann code, in qualitative agreement with the analytical results.

preprint2015arXiv

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

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

preprint2015arXiv

Scaling regimes of a semi-flexible polymer in a rectangular channel

We derive scaling relations for the extension statistics and the confinement free energy for a semi-flexible polymer confined to a channel with a rectangular cross-section. Our motivation are recent numerical results [Gupta {\em et al.}, JCP {\bf 140} (2014) 214901] indicating that extensional fluctuations are quite different in rectangular channels compared to square channels. Our results are of direct relevance for interpreting current experiments on DNA molecules confined to nano-channels, as many experiments are performed for rectangular channels with large aspect ratios while theoretical and simulation results are usually obtained for square channels.

preprint2015arXiv

The role of inertia for the rotation of a nearly spherical particle in a general linear flow

We analyse the angular dynamics of a neutrally buoyant nearly spherical particle immersed in a steady general linear flow. The hydrodynamic torque acting on the particle is obtained by means of a reciprocal theorem, regular perturbation theory exploiting the small eccentricity of the nearly spherical particle, and assuming that inertial effects are small, but finite.

preprint2015arXiv

Tumbling of asymmetric microrods in a microchannel flow

We describe results of measurements of the orientational motion of glass microrods in a microchannel flow, following the orientational motion of particles with different shapes. We determine how the orientational dynamics depends on the shape of the particle and on its initial orientation. We find that the dynamics depends so sensitively on the degree to which particle axisymmetry is broken that it is difficult to find particles that are sufficiently axisymmetric so that they exhibit periodic tumbling ("Jeffery orbits"). The results of our measurements confirm earlier theoretical analysis predicting sensitive dependence on particle shape and its initial orientation. Our results illustrate the different types of orientational dynamics for asymmetric particles predicted by theory.

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

Confined polymers in the extended de Gennes regime

We show that the problem of describing the conformations of a semiflexible polymer confined to a channel can be mapped onto an exactly solvable model in the so-called extended de Gennes regime. This regime (where the polymer is neither weakly nor strongly confined) has recently been studied intensively experimentally and by means of computer simulations. The exact solution predicts precisely how the conformational fluctuations depend upon the channel width and upon the microscopic parameters characterising the physical properties of the polymer.

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.

preprint2014arXiv

The emergence of the rescue effect from the interaction of local stochastic dynamics within a metapopulation

Immigration can rescue local populations from extinction, helping to stabilise a metapopulation. Local population dynamics is important for determining the strength of this rescue effect, but the mechanistic link between local demographic parameters and the rescue effect at the metapopulation level has received very little attention by modellers. We develop an analytical framework that allows us to describe the emergence of the rescue effect from interacting local stochastic dynamics. We show this framework to be applicable to a wide range of spatial scales, providing a powerful and convenient alternative to individual-based models for making predictions concerning the fate of metapopulations. We show that the rescue effect plays an important role in minimising the increase in local extinction probability associated with high demographic stochasticity, but its role is more limited in the case of high local environmental stochasticity of recruitment or survival. While most models postulate the rescue effect, our framework provides an explicit mechanistic link between local dynamics and the emergence of the rescue effect, and more generally the stability of the whole metapopulation.

preprint2013arXiv

Aperiodic tumbling of microrods advected in a microchannel flow

We report on an experimental investigation of the tumbling of microrods in the shear flow of a microchannel (40 x 2.5 x 0.4 mm). The rods are 20 to 30 microns long and their diameters are of the order of 1 micron. Images of the centre-of-mass motion and the orientational dynamics of the rods are recorded using a microscope equipped with a CCD camera. A motorised microscope stage is used to track individual rods as they move along the channel. Automated image analysis determines the position and orientation of a tracked rods in each video frame. We find different behaviours, depending on the particle shape, its initial position, and orientation. First, we observe periodic as well as aperiodic tumbling. Second, the data show that different tumbling trajectories exhibit different sensitivities to external perturbations. These observations can be explained by slight asymmetries of the rods. Third we observe that after some time, initially periodic trajectories lose their phase. We attribute this to drift of the centre of mass of the rod from one to another stream line of the channel flow.

preprint2013arXiv

How Levins' dynamics emerges from a Ricker metapopulation model on the brink of extinction

Understanding the dynamics of metapopulations close to extinction is of vital importance for management. Levins-like models, in which local patches are treated as either occupied or empty, have been used extensively for this purpose, but they ignore the important role of local population dynamics. In this paper, we consider a stochastic metapopulation model where local populations follow a Ricker dynamics, and use this framework to investigate the behaviour of the metapopulation at the brink of extinction. As long as dispersal rates are not too large, the system is shown to have a time evolution consistent with Levins' dynamics. We derive analytical expressions for the colonisation and extinction rates (c and e) in Levins-type models in terms of reproduction, survival, and dispersal parameters of the local populations, providing an avenue to parameterising Levins-like models from the type of information on local demography that is available for a number of species. To facilitate applying our results, we provide a numerical implementation for computing c and e.

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

Monomer distributions and intra-chain collisions of a polymer confined to a channel

We study the conformations of a self-avoiding polymer confined to a channel by computing the cross-sectional distributions of the positions of its monomers. By means of Monte-Carlo simulations for a self-avoiding, freely-jointed chain we determine how the cross-sectional distribution for a given monomer depends on its location in the polymer, and how strongly this distribution is affected by self avoidance. To this end we analyze how the frequency of intra-chain collisions between monomers depends on their spatial position in the channel and on their location within the polymer. We show that most collisions occur between closely neighboring monomers. As a consequence the collision probability depends only weakly on the spatial position of the monomers. Our results explain why the effect of self-avoidance on the monomer distributions is weaker than predicted by mean-field theory. We discuss the relevance of our results for studies of DNA conformations in nanofluidic channels.

preprint2013arXiv

Orientational dynamics of weakly inertial axisymmetric particles in steady viscous flows

The orientational dynamics of weakly inertial axisymmetric particles in a steady flow is investigated. We derive an asymptotic equation of motion for the unit axial vector along the particle symmetry axis, valid for small Stokes number St, and for any axisymmetric particle in any steady linear viscous flow. This reduced dynamics is analysed in two ways, both pertain to the case of a simple shear flow. In this case inertia induces a coupling between precession and nutation. This coupling affects the dynamics of the particle, breaks the degeneracy of the Jeffery orbits, and creates two limiting periodic orbits. We calculate the leading-order Floquet exponents of the limiting periodic orbits and show analytically that prolate objects tend to a tumbling orbit, while oblate objects tend to a log-rolling orbit, in agreement with previous analytical and numerical results. Second, we analyse the role of the limiting orbits when rotational noise is present. We formulate the Fokker-Planck equation describing the orientational distribution of an axisymmetric particle, valid for small St and general Peclet number Pe. Numerical solutions of the Fokker-Planck equation, obtained by means of expansion in spherical harmonics, show that stationary orientational distributions are close to the inertia-free case when Pe St << 1, whereas they are determined by inertial effects, though small, when Pe >> 1/St >> 1.

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.

preprint2012arXiv

Metapopulation dynamics on the brink of extinction

We analyse metapopulation dynamics in terms of an individual-based, stochastic model of a finite metapopulation. We suggest a new approach, using the number of patches in the population as a large parameter. This approach does not require that the number of individuals per patch is large, neither is it necessary to assume a time-scale separation between local population dynamics and migration. Our approach makes it possible to accurately describe the dynamics of metapopulations consisting of many small patches. We focus on metapopulations on the brink of extinction. We estimate the time to extinction and describe the most likely path to extinction. We find that the logarithm of the time to extinction is proportional to the product of two vectors, a vector characterising the distribution of patch population sizes in the quasi-steady state, and a vector -- related to Fisher's reproduction vector -- that quantifies the sensitivity of the quasi-steady state distribution to demographic fluctuations. We compare our analytical results to stochastic simulations of the model, and discuss the range of validity of the analytical expressions. By identifying fast and slow degrees of freedom in the metapopulation dynamics, we show that the dynamics of large metapopulations close to extinction is approximately described by a deterministic equation originally proposed by Levins (1969). We were able to compute the rates in Levins' equation in terms of the parameters of our stochastic, individual-based model. It turns out, however, that the interpretation of the dynamical variable depends strongly on the intrinsic growth rate and carrying capacity of the patches. Only when the growth rate and the carrying capacity are large does the slow variable correspond to the number of patches, as envisaged by Levins. Last but not least, we discuss how our findings relate to other, widely used metapopulation models.

preprint2012arXiv

Orientational correlations in confined DNA

We study how the orientational correlations of DNA confined to nanochannels depend on the channel diameter D by means of Monte Carlo simulations and a mean-field theory. This theory describes DNA conformations in the experimentally relevant regime where the Flory-de Gennes theory does not apply. We show how local correlations determine the dependence of the end-to-end distance of the DNA molecule upon D. Tapered nanochannels provide the necessary resolution in D to study experimentally how the extension of confined DNA molecules depends upon D. Our experimental and theoretical results are in qualitative agreement.

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.

preprint2008arXiv

Shock fragmentation model for gravitational collapse

A cloud of gas collapsing under gravity will fragment. We present a new theory for this process, in which layers shocked gas fragment due to their gravitational instability. Our model explains why angular momentum does not inhibit the collapse process. The theory predicts that the fragmentation process produces objects which are significantly smaller than most stars, implying that accretion onto the fragments plays an essential role in determining the initial masses of stars. This prediction is also consistent with the hypothesis that planets can be produced by gravitational collapse.

preprint2003arXiv

Fano interference and cross-section fluctuations in molecular photodissociation

We derive an expression for the total photodissociation cross section of a molecule incorporating both indirect processes that proceed through excited resonances, and direct processes. We show that this cross section exhibits generalized Beutler-Fano line shapes in the limit of isolated resonances. Assuming that the closed system can be modeled by random matrix theory, we derive the statistical properties of the photodissociation cross section and find that they are significantly affected by the direct processes. We identify a unique signature of the direct processes in the cross-section distribution in the limit of isolated resonances.

preprint2002arXiv

Clustering of SNPs along a chromosome: can the neutral model be rejected?

Single nucleotide polymorphisms (SNPs) often appear in clusters along the length of a chromosome. This is due to variation in local coalescent times caused by,for example, selection or recombination. Here we investigate whether recombination alone (within a neutral model) can cause statistically significant SNP clustering. We measure the extent of SNP clustering as the ratio between the variance of SNPs found in bins of length $l$, and the mean number of SNPs in such bins, $σ^2_l/μ_l$. For a uniform SNP distribution $σ^2_l/μ_l=1$, for clustered SNPs $σ^2_l/μ_l > 1$. Apart from the bin length, three length scales are important when accounting for SNP clustering: The mean distance between neighboring SNPs, $Δ$, the mean length of chromosome segments with constant time to the most recent common ancestor, $\el$, and the total length of the chromosome, $L$. We show that SNP clustering is observed if $Δ< \el \ll L$. Moreover, if $l\ll \el \ll L$, clustering becomes independent of the rate of recombination. We apply our results to the analysis of SNP data sets from mice, and human chromosomes 6 and X. Of the three data sets investigated, the human X chromosome displays the most significant deviation from neutrality.

preprint2000arXiv

An exact-diagonalization study of rare events in disordered conductors

We determine the statistical properties of wave functions in disordered quantum systems by exact diagonalization of one-, two- and quasi-one dimensional tight-binding Hamiltonians. In the quasi-one dimensional case we find that the tails of the distribution of wave-function amplitudes are described by the non-linear sigma-model. In two dimensions, the tails of the distribution function are consistent with a recent prediction based on a direct optimal fluctuation method.

preprint1999arXiv

Statistical properties of eigenvectors in non-Hermitian Gaussian random matrix ensembles

Statistical properties of eigenvectors in non-Hermitian random matrix ensembles are discussed, with an emphasis on correlations between left and right eigenvectors. Two approaches are described. One is an exact calculation for Ginibre's ensemble, in which each matrix element is an independent, identically distributed Gaussian complex random variable. The other is a simpler calculation using $N^{-1}$ as an expansion parameter, where $N$ is the rank of the random matrix: this is applied to Girko's ensemble. Consequences of eigenvector correlations which may be of physical importance in applications are also discussed. It is shown that eigenvalues are much more sensitive to perturbations than in the corresponding Hermitian random matrix ensembles. It is also shown that, in problems with time-evolution governed by a non- Hermitian random matrix, transients are controlled by eigenvector correlations.

preprint1997arXiv

Universal spectral properties of spatially periodic quantum systems with chaotic classical dynamics

We consider a quasi one-dimensional chain of N chaotic scattering elements with periodic boundary conditions. The classical dynamics of this system is dominated by diffusion. The quantum theory, on the other hand, depends crucially on whether the chain is disordered or invariant under lattice translations. In the disordered case, the spectrum is dominated by Anderson localization whereas in the periodic case, the spectrum is arranged in bands. We investigate the special features in the spectral statistics for a periodic chain. For finite N, we define spectral form factors involving correlations both for identical and non-identical Bloch numbers. The short-time regime is treated within the semiclassical approximation, where the spectral form factor can be expressed in terms of a coarse-grained classical propagator which obeys a diffusion equation with periodic boundary conditions. In the long-time regime, the form factor decays algebraically towards an asymptotic constant. In the limit $N\to\infty$, we derive a universal scaling function for the form factor. The theory is supported by numerical results for quasi one-dimensional periodic chains of coupled Sinai billiards.