Search

We solve for waves in a polytropic, stratified medium with a spatially varying background magnetic field that points along a horizontal $x$-direction, and with gravity that is directed along the vertical $z$-direction. Force balance determines the magnitude of the background magnetic field, $B_0^2 \sim z^{n+1}$, where $n$ is the polytropic index. Using numerical and asymptotic methods, we deduce an accurate and explicit dispersion relation for fast pressure-driven waves: $\Omega^2 \sim K\left(2m+n\right) \left[1 + (1/M_\mathrm{A})^2 (4-2\gamma+\cos^2\theta-3\cos^4\theta)/4 \right]$. Here, $\Omega$ is the frequency, $K$ the wavenumber, $\theta$ the angle the wave-vector makes with the background magnetic field, $M_\mathrm{A}$ the Alfv\'enic Mach number, and $m$ an integer representing the eigen state. Applications of our result are in magnetoseismology and nonlinear asteroseismology., Comment: 11 pages, 5 figures

Lagrangian pair dispersion provides insights into mixing in turbulent flows. By direct numerical simulations (DNS) we show that the statistics of pair dispersion in the randomly forced two-dimensional Burgers equation, which is a typical model of shock-dominated turbulence, is very different from its incompressible counterpart because Lagrangian particles get trapped in shocks. We develop a heuristic theoretical framework that accounts for this -- a generalization of the multifractal model -- whose prediction of the scaling of Lagrangian exit times agrees well with our DNS.

Elastic turbulence is the chaotic fluid motion resulting from elastic instabilities due to the addition of polymers in small concentrations at very small Reynolds ($\mbox{Re}$) numbers. Our direct numerical simulations show that elastic turbulence, though a low $\mbox{Re}$ phenomenon, has more in common with classical, Newtonian turbulence than previously thought. In particular, we find power-law spectra for kinetic energy $E(k) \sim k^{-4}$ and polymeric energy $E_{\rm p}(k) \sim k^{-3/2}$, independent of the Deborah ($\mbox{De}$) number. This is further supported by calculation of scale-by-scale energy budget which shows a balance between the viscous term and the polymeric term in the momentum equation. In real space, as expected, the velocity field is smooth, i.e., the velocity difference across a length scale $r$, $\delta u \sim r$ but, crucially, with a non-trivial sub-leading contribution $r^{3/2}$ which we extract by using the second difference of velocity. The structure functions of second difference of velocity up to order $6$ show clear evidence of intermittency/multifractality. We provide additional evidence in support of this intermittent nature by calculating moments of rate of dissipation of kinetic energy averaged over a ball of radius $r$, $\varepsilon_{r}$, from which we compute the multifractal spectrum.

Motivated by wind blowing over water, we use asymptotic methods to study the evolution of short wavelength interfacial waves driven by the combined action of these flows. We solve the Rayleigh equation for the stability of the shear flow, and construct a uniformly valid approximation for the perturbed streamfunction, or eigenfunction. We then expand the real part of the eigenvalue, the phase speed, in a power series of the inverse wavenumber and show that the imaginary part is exponentially small. We give expressions for the growth rates of the Miles (1957) and rippling (e.g., Young & Wolfe 2013) instabilities that are valid for an arbitrary shear flow. The accuracy of the results is demonstrated by a comparison with the exact solution of the eigenvalue problem in the case when both the wind and the current have an exponential profile.

We carry out a detailed study of dynamic multiscaling in the turbulent nonequilibrium, but statistically steady, state of the stochastically forced one-dimensional Burgers equation. We introduce the concept of $\textit{interval collapse times}$ $\tau_{\rm col}$, the time taken for an interval of length $\ell$, demarcated by a pair of Lagrangian tracers, to collapse at a shock. By calculating the dynamic scaling exponent of the order-$p$ moment of $\tau_{\rm col}$, we show that (a) there is $\textit{not one but an infinity of characteristic time scales}$ and (b) the probability distribution function of $\tau_{\rm col}$ is non-Gaussian and has a power-law tail. Our study is based on (a) a theoretical framework that allows us to obtain dynamic-multiscaling exponents analytically, (b) extensive direct numerical simulations, and (c) a careful comparison of the results of (a) and (b). We discuss possible generalizations of our work to dimensions $d >1 $, for the stochastically forced Burgers equation, and to other compressible flows that exhibit turbulence with shocks.

We investigate spectral properties of buoyancy driven bubbly flows. Using high-resolution numerical simulations and phenomenology of homogeneous turbulence, we identify the relevant energy transfer mechanisms. We find: (a) At high enough Galilei number (ratio of the buoyancy to viscous forces) the kinetic energy spectrum shows the Kolmogorov scaling with a power law exponent $-5/3$ for the range of scales between the bubble diameter and the dissipation scale ($\eta$). (b) For scales smaller than $\eta$, the physics of pseudo-turbulence is recovered.

We use direct numerical simulations to study homogeneous, and isotropic turbulent flows of dilute polymer solutions at high Reynolds and Deborah numbers. We find that for small wavenumbers $k$, the kinetic energy spectrum shows Kolmogorov--like behavior which crosses over at a larger $k$ to a novel, elastic scaling regime, $E(k) \sim k^{-\xi}$, with $\xi \approx 2.3$. We study the contribution of the polymers to the flux of kinetic energy through scales, and find that it can be decomposed into two parts: one increase in effective viscous dissipation, and a purely elastic contribution that dominates over the nonlinear flux in the range of $k$ over which the elastic scaling is observed. The multiscale balance between the two fluxes determines the crossover wavenumber which depends non-monotically on the Deborah number. Consistently, structure functions also show two scaling ranges, with intermittency present in both of them in equal measure.

The rate of collision and the relative velocities of the colliding particles in turbulent flows is a crucial part of several natural phenomena, e.g., rain formation in warm clouds and planetesimal formation in a protoplanetary disks. The particles are often modeled as passive, but heavy and inertial. Within this model, large relative velocities emerge due to formation of singularities (caustics) of in the gradient matrix of the velocities of the particles. Using extensive direct numerical simulations of heavy particles in both two (direct and inverse cascade) and three dimensional turbulent flows we calculate the rate of formation of caustics, $J$ as a function of the Stokes number (${\rm St}$).The best approximation to our data is $J \sim \exp(-C/{\rm St})$, in the limit ${\rm St} \to 0 $ where $C$ is a non-universal constant., Comment: Consistent with the version accepted for publication in Philosophical Transactions of the Royal Society A

When wind blows over water, ripples are generated on the water surface. These ripples can be regarded as perturbations of the wind field, which is modelled as a parallel inviscid flow. For a given wavenumber $k$, the perturbed streamfunction of the wind field and the complex phase speed are the eigenfunction and the eigenvalue of the so-called Rayleigh equation in a semi-infinite domain. Because of the small air-water density ratio, $\rho_{\rm{a}}/\rho_{\rm{w}}\equiv\epsilon<<1$, the wind and the ripples are weakly coupled, and the eigenvalue problem can be solved perturbatively. At the leading order, the eigenvalue is equal to the phase speed $c_0$ of surface waves. At order $\epsilon$, the eigenvalue has a finite imaginary part, which implies growth. Miles (1957) showed that the growth rate is proportional to the square modulus of the leading-order eigenfunction evaluated at the so-called critical level $z=z_c$, where the wind speed is equal to $c_0$ and the waves extract energy from the wind. Here, we construct uniform asymptotic approximations of the leading-order eigenfunction for long waves, which we use to calculate the growth rate as a function of $k$. In the strong wind limit, we find that the fastest growing wave is such that the aerodynamic pressure is in phase with the wave slope. The results are confirmed numerically.

We present a Direct Numerical Simulation (DNS) study of buoyancy-driven bubbly flows in the presence of large scale driving that generates turbulence. On increasing the turbulence intensity: (a) the bubble trajectories become more curved, and (b) the average rise velocity of the bubbles decreases. We find that the energy spectrum of the flow shows a pseudo-turbulence scaling for length scales smaller than the bubble diameter and a Kolmogorov scaling for scales larger than the bubble diameter. We conduct a scale-by-scale energy budget analysis to understand the scaling behaviour observed in the spectrum. Although our bubbles are weakly buoyant, the statistical properties of our DNS are consistent with the experiments that investigate turbulence modulation by air bubbles in water.

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., Comment: 6 pages, 3 figures, supplementary material

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., Comment: Accepted for publication in Astrophysical J

We study the spreading of viruses, such as SARS-CoV-2, by airborne aerosols, via a new first-passage-time problem for Lagrangian tracers that are advected by a turbulent flow: By direct numerical simulations of the three-dimensional (3D) incompressible, Navier-Stokes equation, we obtain the time $t_R$ at which a tracer, initially at the origin of a sphere of radius $R$, crosses the surface of the sphere \textit{for the first time}. We obtain the probability distribution function $\mathcal{P}(R,t_R)$ and show that it displays two qualitatively different behaviors: (a) for $R \ll L_{\rm I}$, $\mathcal{P}(R,t_R)$ has a power-law tail $\sim t_R^{-\alpha}$, with the exponent $\alpha = 4$ and $L_{\rm I}$ the integral scale of the turbulent flow; (b) for $l_{\rm I} \lesssim R $, the tail of $\mathcal{P}(R,t_R)$ decays exponentially. We develop models that allow us to obtain these asymptotic behaviors analytically. We show how to use $\mathcal{P}(R,t_R)$ to develop social-distancing guidelines for the mitigation of the spreading of airborne aerosols with viruses such as SARS-CoV-2.

We present Direct Numerical Simulation (DNS) of heavy inertial particles (dust) immersed in two-dimensional turbulent flow (gas). The dust are modeled as mono-dispersed heavy particles capable of modifying the flow through two-way coupling. By varying the Stokes number (St) and the mass-loading parameter $({\phi}_{\rm m})$, we study the clustering phenomenon and the gas phase kinetic energy spectra. We find that the dust-dust correlation dimension $(d_2)$ also depends on ${\phi}_{\rm m}$. In particular, clustering decreases as mass-loading $({\phi}_{\rm m})$, is increased. In the kinetic energy spectra of gas we show: (i) emergence of a new scaling regime, (ii) the scaling exponent in this regime is not unique but rather a function of both St and $({\phi}_{\rm m} )$. Using a scale-by-scale enstrophy budget analysis we show in the new scaling regime, viscous dissipation due to the gas balances back-reaction from the dust.

We perform direct numerical simulations of the flow through a model of a deformable porous medium. Our model is a two-dimensional hexagonal lattice, with defects, of soft elastic cylindrical pillars, with elastic shear modulus $G$, immersed in a liquid. We use a two-phase approach: the liquid phase is a viscous fluid and the solid phase is modeled as an incompressible viscoelastic material, whose complete nonlinear structural response is considered. We observe that the Darcy flux ($q$) is a nonlinear function -- steeper than linear -- of the pressure-difference ($\Delta P$) across the medium. Furthermore, the flux is larger for a softer medium (smaller $G$). We construct a theory of this super-linear behavior by modelling the channels between the solid cylinders as elastic channels whose walls are made of a material with a linear constitutive relation but can undergo large deformation. Our theory further predicts that the flow permeability is a universal function of $\Delta P/G$, which is confirmed by the present simulations., Comment: 6 pages, 3 figures, Some minor changes (including the title) from the previous submission

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., Comment: 12 pages, 3 figures

We consider the kinematic fluctuation dynamo problem in a flow that is random, white-in-time, with both solenoidal and potential components. This model is a generalization of the well-studied Kazantsev model. If both the solenoidal and potential parts have the same scaling exponent, then, as the compressibility of the flow increases, the growth rate decreases but remains positive. If the scaling exponents for the solenoidal and potential parts differ, in particular if they correspond to typical Kolmogorov and Burgers values, we again find that an increase in compressibility slows down the growth rate but does not turn it off. The slow down is, however, weaker and the critical magnetic Reynolds number is lower than when both the solenoidal and potential components display the Kolmogorov scaling. Intriguingly, we find that there exist cases, when the potential part is smoother than the solenoidal part, for which an increase in compressibility increases the growth rate. We also find that the critical value of the scaling exponent above which a dynamo is seen is unity irrespective of the compressibility. Finally, we realize that the dimension $d = 3$ is special, since for all other values of $d$ the critical exponent is higher and depends on the compressibility., Comment: 12 pages, 6 figures

We present an extensive numerical study of the time irreversibility of the dynamics of heavy inertial particles in three-dimensional, statistically homogeneous and isotropic turbulent flows. We show that the probability density function (PDF) of the increment, $W(\tau)$, of a particle's energy over a time-scale $\tau$ is non-Gaussian, and skewed towards negative values. This implies that, on average, particles gain energy over a period of time that is longer than the duration over which they lose energy. We call this $\textit{slow gain}$ and $\textit{fast loss}$. We find that the third moment of $W(\tau)$ scales as $\tau^3$, for small values of $\tau$. We show that the PDF of power-input $p$ is negatively skewed too; we use this skewness ${\rm Ir}$ as a measure of the time-irreversibility and we demonstrate that it increases sharply with the Stokes number ${\rm St}$, for small ${\rm St}$; this increase slows down at ${\rm St} \simeq 1$. Furthermore, we obtain the PDFs of $t^+$ and $t^-$, the times over which $p$ has, respectively, positive or negative signs, i.e., the particle gains or loses energy. We obtain from these PDFs a direct and natural quantification of the the slow-gain and fast-loss of the particles, because these PDFs possess exponential tails, whence we infer the characteristic loss and gain times $t_{\rm loss}$ and $t_{\rm gain}$, respectively; and we obtain $t_{\rm loss} < t_{\rm gain}$, for all the cases we have considered. Finally, we show that the slow-gain in energy of the particles is equally likely in vortical or strain-dominated regions of the flow; in contrast, the fast-loss of energy occurs with greater probability in the latter than in the former.

We perform direct numerical simulations (DNS) of passive heavy inertial particles (dust) in homogeneous and isotropic two-dimensional turbulent flows (gas) for a range of Stokes number, ${\rm St} < 1$, using both Lagrangian and Eulerian approach (with a shock-capturing scheme). We find that: The dust-density field in our Eulerian simulations have the same correlation dimension $d_2$ as obtained from the clustering of particles in the Lagrangian simulations for ${\rm St} < 1$; The cumulative probability distribution function of the dust-density coarse-grained over a scale $r$ in the inertial range has a left-tail with a power-law fall-off indicating presence of voids; The energy spectrum of the dust-velocity has a power-law range with an exponent that is same as the gas-velocity spectrum except at very high Fourier modes; The compressibility of the dust-velocity field is proportional to ${\rm St}^2$. We quantify the topological properties of the dust-velocity and the gas-velocity through their gradient matrices, called $\mathcal{A}$ and $\mathcal{B}$, respectively. The topological properties of $\mathcal{B}$ are the same in Eulerian and Lagrangian frames only if the Eulerian data are weighed by the dust-density -- a correspondence that we use to study Lagrangian properties of $\mathcal{A}$. In the Lagrangian frame, the mean value of the trace of $\mathcal{A} \sim - \exp(-C/{\rm St}$, with a constant $C\approx 0.1$. The topology of the dust-velocity fields shows that as ${\rm St} increases the contribution to negative divergence comes mostly from saddles and the contribution to positive divergence comes from both vortices and saddles. Compared to the Eulerian case, the density-weighed Eulerian case has less inward spirals and more converging saddles. Outward spirals are the least probable topological structures in both cases., Comment: abstract condensed to fit arxiv format

We use direct numerical simulations to calculate the joint probability density function of the relative distance $R$ and relative radial velocity component $V_R$ for a pair of heavy inertial particles suspended in homogeneous and isotropic turbulent flows. At small scales the distribution is scale invariant, with a scaling exponent that is related to the particle-particle correlation dimension in phase space, $D_2$. It was argued [1, 2] that the scale invariant part of the distribution has two asymptotic regimes: (1) $|V_R| \ll R$ where the distribution depends solely on $R$; and (2) $|V_R| \gg R$ where the distribution is a function of $|V_R|$ alone. The probability distributions in these two regimes are matched along a straight line $|V_R| = z^\ast R$. Our simulations confirm that this is indeed correct. We further obtain $D_2$ and $z^\ast$ as a function of the Stokes number, ${\rm St}$. The former depends non-monotonically on ${\rm St}$ with a minimum at about ${\rm St} \approx 0.7$ and the latter has only a weak dependence on ${\rm St}$.

We study suspensions of deformable (viscoelastic) spheres in a Newtonian solvent in plane Couette geometry, by means of direct numerical simulations. We find that in the limit of vanishing inertia the effective viscosity $\mu$ of the suspension increases as the volume-fraction occupied by the spheres $\Phi$ increases and decreases as the elastic modulus of the spheres $G$ decreases; the function $\mu(\Phi,G)$ collapses to an universal function, $\mu(\Phi_e)$, with a reduced effective volume fraction $\Phi_e(\Phi,G)$. Remarkably, the function $\mu(\Phi_e)$ is the well-known Eilers fit that describes the rheology of suspension of rigid spheres at all $\Phi$. Our results suggest new ways to interpret macro-rheology of blood.

The effect of turbulence on the mass and heat transfer between small heavy inertial particles (HIP) and an embedding fluid is studied. Two effects are identified. The first effect is due to the relative velocity between the fluid and the particles, and a model for the relative velocity is presented. The second effect is due to the clustering of particles, where the mass transfer rate is inhibited due to the rapid depletion of the consumed species inside the dense particle clusters. This last effect is relevant for large Damkohler numbers and it may totally control the mass transfer rate for Damkohler numbers larger than unity. A model that describes how this effect should be incorporated into existing particle simulation tools is presented., Comment: 14 pages, 7 figures

We obtain the probability distribution functions (PDFs) of the time that a Lagrangian tracer or a heavy inertial particle spends in vortical or strain-dominated regions of a turbulent flow, by carrying out direct numerical simulation (DNS) of such particles advected by statistically steady, homogeneous and isotropic turbulence in the forced, three-dimensional, incompressible Navier-Stokes equation. We use the two invariants, $Q$ and $R$, of the velocity-gradient tensor to distinguish between vortical and strain-dominated regions of the flow and partition the $Q-R$ plane into four different regions depending on the topology of the flow; out of these four regions two correspond to vorticity-dominated regions of the flow and two correspond to strain-dominated ones. We obtain $Q$ and $R$ along the trajectories of tracers and heavy inertial particles and find out the time $\mathrm{t_{pers}}$ for which they remain in one of the four regions of the $Q-R$ plane. We find that the PDFs of $\mathrm{t_{pers}}$ display exponentially decaying tails for all four regions for tracers and heavy inertial particles. From these PDFs we extract characteristic times scales, which help us to quantify the time that such particles spend in vortical or strain-dominated regions of the flow.

The effect of turbulence on the heterogeneous (solid-fluid) reactions of solid particles is studied numerically with Direct Numerical Simulations (DNS). A simplified reaction system is used, where the solid-fluid reaction is represented by a single isothermal reaction step. It is found that, due to the clustering of particles by the isotropic turbulence, the overall reaction rate is entirely controlled by the turbulence for large Damk\"ohler numbers. The particle clustering significantly slows down the reaction rate for increasing Damk\"ohler numbers which reaches an asymptotic limit that can be analytically derived. This implies that the effect of turbulence on heterogeneously reacting particles should be included in models that are used in CFD simulations of e.g. char burnout in combustors or gasifiers. Such a model, based on the chemical and turbulent time scales, is here proposed for the heterogeneous reaction rate in the presence of turbulence., Comment: Published online in the Proceedings of the Combustion Institute

We show, by direct numerical simulations, that heavy inertial particles (characterized by Stokes number $\St$) in inhomogeneously forced statistically stationary isothermal turbulent flows cluster at the minima of mean-square turbulent velocity. Two turbulent transport processes, turbophoresis and turbulent diffusion together determine the spatial distribution of the particles. If the turbulent diffusivity is assumed to scale with turbulent root mean square velocity, as is the case for homogeneous turbulence, the turbophoretic coefficient can be calculated. Indeed, for the above assumption, the non-dimensional product of the turbophoretic coefficient and the rms velocity is shown to increase with $St$ for small $St$, reach a maxima for $St\approx 10$ and decrease as $\sim St^{-0.33}$ for large $St$., Comment: 7 pages, 3 Figures

We uncover universal statistical properties of the trajectories of heavy inertial particles in three-dimensional, statistically steady, homogeneous, and isotropic turbulent flows by extensive direct numerical simulations. We show that the probability distribution functions (PDFs) $P(\phi)$, of the angle $\phi$ between the Eulerian velocity ${\bf u}$ and the particle velocity ${\bf v}$, at this point and time, shows a power-law region in which $P(\phi) \sim \phi^{-\gamma}$, with a new universal exponent $\gamma \simeq 4$. Furthermore, the PDFs of the trajectory curvature $\kappa$ and modulus $\theta$ of the torsion $\vartheta$ have power-law tails that scale, respectively, as $P(\kappa) \sim \kappa^{-h_\kappa}$, as $\kappa \to \infty$, and $P(\theta) \sim \theta^{-h_\theta}$, as $\theta \to \infty$, with exponents $h_\kappa \simeq 2.5$ and $h_\theta \simeq 3$ that are universal to the extent that they do not depend on the Stokes number ${\rm St}$ (given our error bars). We also show that $\gamma$, $h_\kappa$ and $h_\theta$ can be obtained by using simple stochastic models. We characterize the complexity of heavy-particle trajectories by the number $N_{\rm I}(t,{\rm St})$ of points (up until time $t$) at which $\vartheta$ changes sign. We show that $n_{\rm I}({\rm St}) \equiv \lim_{t\to\infty} \frac{N_{\rm I}(t,{\rm St})}{t} \sim {\rm St}^{-\Delta}$, with $\Delta \simeq 0.4$ a universal exponent.

We obtain, by extensive direct numerical simulations, trajectories of heavy inertial particles in two-dimensional, statistically steady, homogeneous, and isotropic turbulent flows, with friction. We show that the probability distribution function $\mathcal{P}(\kappa)$, of the trajectory curvature $\kappa$, is such that, as $\kappa \to \infty$, $\mathcal{P}(\kappa) \sim \kappa^{-h_{\rm r}}$, with $h_{\rm r} = 2.07 \pm 0.09$. The exponent $h_{\rm r}$ is universal, insofar as it is independent of the Stokes number ($\rm{St}$) and the energy-injection wave number. We show that this exponent lies within error bars of their counterparts for trajectories of Lagrangian tracers. We demonstrate that the complexity of heavy-particle trajectories can be characterized by the number $N_{\rm I}(t,{\rm St})$ of inflection points (up until time $t$) in the trajectory and $n_{\rm I} ({\rm St}) \equiv \lim_{t\to\infty} \frac{N_{\rm I} (t,{\rm St})}{t} \sim {\rm St}^{-\Delta}$, where the exponent $\Delta = 0.33 \pm0.02$ is also universal., Comment: submitted to Phys. Rev. Lett. The supplementary information is added as an appendix

Guided by extensive numerical simulations, we propose a microfluidic device that can sort elastic capsules by their deformability. The device consists of a duct embedded with a semi-cylindrical obstacle, and a diffuser which further enhances the sorting capability. We demonstrate that the device can operate reasonably well under changes in the initial position of the the capsule. The efficiency of the device remains essentially unaltered under small changes of the obstacle shape (from semi-circular to semi-elliptic cross-section). Confinement along the direction perpendicular to the plane of the device increases its efficiency. This work is the first numerical study of cell sorting by a realistic microfluidic device., Comment: preprint, 7 pages 6 figures

As a test bed for the growth of protoplanetary bodies in a turbulent circumstellar disk we examine the fate of a boulder using direct numerical simulations of particle seeded gas flowing around it. We provide an accurate description of the flow by imposing no-slip and non-penetrating boundary conditions on the boulder surface using the immersed boundary method pioneered by Peskin (2002). Advected by the turbulent disk flow, the dust grains collide with the boulder and we compute the probability density function (PDF) of the normal component of the collisional velocity. Through this examination of the statistics of collisional velocities we test the recently developed concept of collisional fusion which provides a physical basis for a range of collisional velocities exhibiting perfect sticking. A boulder can then grow sufficiently rapidly to settle into a Keplerian orbit on disk evolution time scales., Comment: Astrophysical Journal, in press

Using direct numerical simulations, we verify that "flow IV" of Roberts (1972) exhibits dynamo action dominated by horizontally averaged large-scale magnetic field. With the test-field method we compute the turbulent magnetic diffusivity and find that it is negative and overcomes the molecular diffusivity, thus explaining quantitatively the large-scale dynamo for magnetic Reynolds numbers above $\approx8$. As expected for a dynamo of this type, but contrary to $\alpha$-effect dynamos, the two horizontal field components grow independently of each other and have arbitrary amplitude ratios and phase differences. Small length scales of the mean magnetic field are shown to be stabilized by the turbulent magnetic diffusivity becoming positive at larger wavenumbers. Oscillatory decaying or growing solutions have also been found in certain wavenumber intervals and sufficiently large values of the magnetic Reynolds number. For magnetic Reynolds numbers below $\approx0.5$ the turbulent magnetic diffusivity is confirmed to be positive, as expected for all incompressible flows. Earlier claims of a dynamo driven by a modified Taylor-Green flow through negative eddy diffusivity could not be confirmed., Comment: 7 pages, 9 figures, accepted to MNRAS

Shear thickening appears as an increase of the viscosity of a dense suspension with the shear rate, sometimes sudden and violent at high volume fraction. Its origin for noncolloidal suspension with non-negligible inertial effects is still debated. Here we consider a simple shear flow and demonstrate that fluid inertia causes a strong microstructure anisotropy that results in the formation of a shadow region with no relative flux of particles. We show that shear thickening at finite inertia can be explained as an increase of the effective volume fraction when considering the dynamically excluded volume due to these shadow regions

We study spontaneous breakdown of chiral symmetry during the nonlinear evolution of the Tayler instability. We start with an initial steady state of zero helicity. Within linearized perturbation calculations, helical perturbations of this initial state have the same growth rate for either sign of helicity. Direct numerical simulations (DNS) of the fully nonlinear equations, however, show that an infinitesimal excess of one sign of helicity in the initial perturbation gives rise to a saturated helical state. We further show that this symmetry breaking can be described by weakly nonlinear finite--amplitude equations with undetermined coefficients which can be deduced solely from symmetry consideration. By fitting solutions of the amplitude equations to data from DNS we further determine the coefficients of the amplitude equations., Comment: 8 pages, 5 figures, Published on Phys. Rev. E

In two-dimensional forced Navier-Stokes turbulence, energy cascades to the largest scales in the system to form a pair of coherent vortices known as the Bose condensate. We show, both numerically and analytically, that the energy condensation saturates and the system reaches a statistically stationary state. The time scale of saturation is inversely proportional to the viscosity and the saturation energy level is determined by both the viscosity and the force. We further show that, without sufficient resolution to resolve the small-scale enstrophy spectrum, numerical simulations can give a spurious result for the saturation energy level. We also find that the movement of the condensate is similar to the motion of an inertial particle with an effective drag force. Furthermore, we show that the profile of the saturated coherent vortices can be described by a Gaussian core with exponential wings., Comment: 9 pages, 8 figures, updated to match the published version

We present the first demonstration of the negative effective magnetic pressure instability in direct numerical simulations of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta. By the action of this instability, an initially uniform horizontal magnetic field forms flux concentrations whose scale is large compared to the turbulent scale. We further show that the magnetic energy of these large-scale structures is only weakly dependent on the magnetic Reynolds number, provided its value is large enough for the instability to be excited. Our results support earlier mean-field calculations and analytic work which identified this instability. Applications to the formation of active regions in the Sun are discussed., Comment: 4 pages, 5 figures, published in ApJL

We obtain, by extensive direct numerical simulations, time-dependent and equal-time structure functions for the vorticity, in both quasi-Lagrangian and Eulerian frames, for the direct-cascade regime in two-dimensional fluid turbulence with air-drag-induced friction. We show that differ- ent ways of extracting time scales from these time-dependent structure functions lead to different dynamic-multiscaling exponents, which are related to equal-time multiscaling exponents by different classes of bridge relations; for a representative value of the friction we verify that, given our error bars, these bridge relations hold., Comment: Published in Phys. Rev. Lett, Supplemental materials added as appendix

We present direct numerical simulations of turbulent channel flow with passive Lagrangian polymers. To understand the polymer behavior we investigate the behavior of infinitesimal line elements and calculate the probability distribution function (PDF) of finite-time Lyapunov exponents and from them the corresponding Cramer's function for the channel flow. We study the statistics of polymer elongation for both the Oldroyd-B model (for Weissenberg number $\Wi <1$) and the FENE model. We use the location of the minima of the Cramer's function to define the Weissenberg number precisely such that we observe coil-stretch transition at $\Wi\approx1$. We find agreement with earlier analytical predictions for PDF of polymer extensions made by Balkovsky, Fouxon and Lebedev [Phys. Rev. Lett., 84, 4765 (2000).] for linear polymers (Oldroyd-B model) with $\Wi<1$ and by Chertkov [Phys. Rev. Lett., 84, 4761 (2000).] for nonlinear FENE-P model of polymers. For $\Wi>1$ (FENE model) the polymer are significantly more stretched near the wall than at the center of the channel where the flow is closer to homogenous isotropic turbulence. Furthermore near the wall the polymers show a strong tendency to orient along the stream-wise direction of the flow but near the centerline the statistics of orientation of the polymers is consistent with analogous results obtained recently in homogeneous and isotropic flows., Comment: consistent with the published version

Evidence for a parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of an instability in magnetohydrodynamics appears to be the first of its kind and is similar to chiral symmetry breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are being discussed., Comment: consistent with the published version

We present a natural framework for studying the persistence problem in two-dimensional fluid turbulence by using the Okubo-Weiss parameter $\Lambda$ to distinguish between vortical and extensional regions. We then use a direct numerical simulation (DNS) of the two-dimensional, incompressible Navier--Stokes equation with Ekman friction to study probability distribution functions (PDFs) of the persistence times of vortical and extensional regions by employing both Eulerian and Lagrangian measurements. We find that, in the Eulerian case, the persistence-time PDFs have exponential tails; by contrast, this PDF for Lagrangian particles, in vortical regions, has a power-law tail with an exponent $\theta=2.9\pm0.2$., Comment: consistent with the published version

We carry out a direct numerical simulation (DNS) study that reveals the effects of polymers on statistically steady, forced, homogeneous, isotropic fluid turbulence. We find clear manifestations of dissipation-reduction phenomena: On the addition of polymers to the turbulent fluid, we obtain a reduction in the energy dissipation rate, a significant modification of the fluid energy spectrum, especially in the deep-dissipation range, a suppression of small-scale intermittency, and a decrease in small-scale vorticity filaments. We also compare our results with recent experiments and earlier DNS studies of decaying fluid turbulence with polymer additives., Comment: consistent with the published version

We systematise the study of dynamic multiscaling of time-dependent structure functions in different models of passive-scalar and fluid turbulence. We show that, by suitably normalising these structure functions, we can eliminate their dependence on the origin of time at which we start our measurements and that these normalised structure functions yield the same linear bridge relations that relate the dynamic-multiscaling and equal-time exponents for statistically steady turbulence. We show analytically, for both the Kraichnan Model of passive-scalar turbulence and its shell model analogue, and numerically, for the GOY shell model of fluid turbulence and a shell model for passive-scalar turbulence, that these exponents and bridge relations are the same for statistically steady and decaying turbulence. Thus we provide strong evidence for dynamic universality, i.e., dynamic-multiscaling exponents do not depend on whether the turbulence decays or is statistically steady., Comment: 29 pages, 18 figures. For figures with better resolution, please download the preprint available at the website of the International Collaboration of Turbulence Research

We give an overview of the progress that has been made in recent years in understanding the dynamic multiscaling of homogeneous, isotropic turbulence and related problems. We emphasise the similarity of this problem with the dynamic scaling of time-dependent correlation functions in the vicinity of a critical point in, e.g., a spin system. The universality of dynamic-multiscaling exponents in fluid turbulence is explored by detailed simulations of the GOY shell model for fluid turbulence., Comment: 13 pages, 10 figures, Submitted for the Proceedings of the Statphys 23 Conference

The existence of drag reduction by polymer additives, well established for wall-bounded turbulent flows, is controversial in homogeneous, isotropic turbulence. To settle this controversy we carry out a high-resolution direct numerical simulation (DNS) of decaying, homogeneous, isotropic turbulence with polymer additives. Our study reveals clear manifestations of drag-reduction-type phenomena: On the addition of polymers to the turbulent fluid we obtain a reduction in the energy dissipation rate, a significant modification of the fluid energy spectrum especially in the deep-dissipation range, a suppression of small-scale intermittency, and a decrease in small-scale vorticity filaments., Comment: 4 pages, 3 figures