Search

We solve for the motion of charged particles in a helical time-periodic ABC (Arnold-Beltrami-Childress) magnetic field. The magnetic field lines of a stationary ABC field with coefficients A = B = C = 1 are chaotic, and we show that the motion of a charged particle in such a field is also chaotic at late times with positive Lyapunov exponent. We further show that in time-periodic ABC fields, the kinetic energy of a charged particle can increase indefinitely with time. At late times the mean kinetic energy grows as a power law in time with an exponent that approaches unity. For an initial distribution of particles, whose kinetic energy is uniformly distributed within some interval, the probability density function of kinetic energy is, at late times, close to a Gaussian but with steeper tails. [ABSTRACT FROM AUTHOR]

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 ≈ 1. 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 streamwise direction of the flow, but near the center line the statistics of orientation of the polymers is consistent with analogous results obtained recently in homogeneous and isotropic flows [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 different 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. [ABSTRACT FROM AUTHOR]

Evidence for the 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 nonhelical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of this magnetohydrodynamic instability appears to be the first of its kind and has properties similar to those in chiral symmetry breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are discussed. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We present direct numerical simulation of heavy inertial particles (dust) immersed in two-dimensional turbulent flow (gas). The dust particles are modeled as monodispersed heavy particles capable of modifying the flow through two-way coupling. By varying the Stokes number (St) and the mass-loading parameter (ϕ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 ϕm. In particular, clustering decreases as mass loading (ϕm) is increased. In the kinetic energy spectra of gas we show (i) the emergence of a different scaling regime and that (ii) the scaling exponent in this regime is not universal but a function of both St and ϕm . Using a scale-by-scale enstrophy budget analysis we show that in this emerged scaling regime, which we call the dust-dissipative range, viscous dissipation due to the gas balances the back-reaction from the dust. [ABSTRACT FROM AUTHOR]

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(τ), of a particle's energy over a time scale τ is non-Gaussian, and skewed toward 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 slow gain and fast loss. We find that the third moment of W(τ) scales as τ³ for small values of τ. We show that the PDF of power-input p is negatively skewed too; we use this skewness Ir as a measure of the time irreversibility and we demonstrate that it increases sharply with the Stokes number St for small St; this increase slows down at St≃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 slow gain and fast loss of the energy of the particles, because these PDFs possess exponential tails from which we infer the characteristic loss and gain times tloss and tgain, respectively, and we obtain tloss

We use direct numerical simulations to calculate the joint probability density function of the relative distance R and relative radial velocity component VR 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, D2. It was argued [K. Gustavsson and B. Mehlig, Phys. Rev. E 84, 045304 (2011); J. Turbul. 15, 34 (2014)] that the scale invariant part of the distribution has two asymptotic regimes: (1) |VR|«R, where the distribution depends solely on R, and (2) |VR|»R, where the distribution is a function of |VR| alone. The probability distributions in these two regimes are matched along a straight line: |VR|=z*R. Our simulations confirm that this is indeed correct. We further obtain D2 and Z* as a function of the Stokes number, St. The former depends nonmonotonically on St with a minimum at about St≈0.7 and the latter has only a weak dependence on St. [ABSTRACT FROM AUTHOR]

Small particles in suspension in a turbulent fluid have trajectories that do not follow the pathlines of the flow exactly. We investigate the statistics of the angle of deviation ɸ between the particle and fluid velocities. We show that, when the effects of particle inertia are small, the probability distribution function (PDF) Pɸ of this deviation angle shows a power-law region in which Pɸ~ɸ-4. We also find that the PDFs of the trajectory curvature κ and modulus θ of the torsion ϑ have power-law tails that scale, respectively, as Pκ~κ-5/2, as κ→∞, and Pθ~θ-3, as θ→∞: These exponents are in agreement with those previously observed for fluid pathlines. We propose a way to measure the complexity of heavy-particle trajectories by the number NI(t,St) of points (up until time t) at which the torsion changes sign. We present numerical evidence that nI(St)=limt→∞NI(t,St)/t ~St-Δ for large St, with Δ≃0.5. [ABSTRACT FROM AUTHOR]

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 simulations 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 tpers for which they remain in one of the four regions of the Q-R plane. We find that the PDFs of tpers display exponentially decaying tails for all four regions for tracers and heavy inertial particles. From these PDFs we extract characteristic time scales, which help us to quantify the time that such particles spend in vortical or strain-dominated regions of the flow. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We study the rheology of confined suspensions of neutrally buoyant rigid monodisperse spheres in plane-Couette flow using direct numerical simulations. We find that if the width of the channel is a (small) integer multiple of the sphere diameter, the spheres self-organize into two-dimensional layers that slide on each other and the effective viscosity of the suspension is significantly reduced. Each two-dimensional layer is found to be structurally liquidlike but its dynamics is frozen in time. [ABSTRACT FROM AUTHOR]

We investigate the 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 a high enough Galilei number (ratio of the buoyancy to viscous forces) the velocity power 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 (η). (b) For scales smaller than η, the physics of pseudo-turbulence is recovered.

We study the buckling of pressurized spherical shells by Monte Carlo simulations in which the detailed balance is explicitly broken-thereby driving the shell to be active, out of thermal equilibrium. Such a shell typically has either higher (active) or lower (sedate) fluctuations compared to one in thermal equilibrium depending on how the detailed balance is broken. We show that, for the same set of elastic parameters, a shell that is not buckled in thermal equilibrium can be buckled if turned active. Similarly a shell that is buckled in thermal equilibrium can unbuckle if sedated. Based on this result, we suggest that it is possible to experimentally design microscopic elastic shells whose buckling can be optically controlled.

The flow of Newtonian fluid at low Reynolds number is, in general, regular and time-reversible due to absence of nonlinear effects. For example, if the fluid is sheared by its boundary motion that is subsequently reversed, then all the fluid elements return to their initial positions. Consequently, mixing in microchannels happens solely due to molecular diffusion and is very slow. Here, we show, numerically, that the introduction of a single, freely floating, flexible filament in a time-periodic linear shear flow can break reversibility and give rise to chaos due to elastic nonlinearities, if the bending rigidity of the filament is within a carefully chosen range. Within this range, not only the shape of the filament is spatiotemporally chaotic, but also the flow is an efficient mixer. Overall, we find five dynamical phases: the shape of a stiff filament is time-invariant-either straight or buckled; it undergoes a period-two bifurcation as the filament is made softer; becomes spatiotemporally chaotic for even softer filaments but, surprisingly, the chaos is suppressed if bending rigidity is decreased further.

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.

We present a natural framework for studying the persistence problem in two-dimensional fluid turbulence by using the Okubo-Weiss parameter Λ to distinguish between vortical and extensional regions. We then use a direct numerical simulation 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 θ=2.9±0.2.

We carry out a direct numerical simulation (DNS) study that reveals the effects of polymers on statistically steady, forced, homogeneous, and 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; and signatures of the suppression of small-scale structures, including 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.

We carry out systematic and high-resolution studies of dynamo action in a shell model for magnetohydrodynamic (MHD) turbulence over wide ranges of the magnetic Prandtl number PrM and the magnetic Reynolds number ReM. Our study suggests that it is natural to think of dynamo onset as a nonequilibrium first-order phase transition between two different turbulent, but statistically steady, states. The ratio of the magnetic and kinetic energies is a convenient order parameter for this transition. By using this order parameter, we obtain the stability diagram (or nonequilibrium phase diagram) for dynamo formation in our MHD shell model in the (PrM-1,ReM) plane. The dynamo boundary, which separates dynamo and no-dynamo regions, appears to have a fractal character. We obtain a hysteretic behavior of the order parameter across this boundary and suggestions of nucleation-type phenomena.

Simulations of stochastically forced shear-flow turbulence in a shearing-periodic domain are used to study the spontaneous generation of large-scale flow patterns in the direction perpendicular to the plane of the shear. Based on an analysis of the resulting large-scale velocity correlations it is argued that the mechanism behind this phenomenon could be the mean-vorticity dynamo effect pioneered by Elperin, Kleeorin, and Rogachevskii [Phys. Rev. E 68, 016311 (2003)]. This effect is based on the anisotropy of the eddy viscosity tensor. One of its components may be able to replenish cross-stream mean flows by acting upon the streamwise component of the mean flow. Shear, in turn, closes the loop by acting upon the cross-stream mean flow to produce stronger streamwise mean flows. The diagonal component of the eddy viscosity is found to be of the order of the rms turbulent velocity divided by the wave number of the energy-carrying eddies.

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

We present the first study of the dynamic scaling or multiscaling of passive-scalar turbulence. For the Kraichnan version of passive-scalar turbulence we show analytically, in both Eulerian and quasi-Lagrangian frameworks, that simple dynamic scaling is obtained but with different dynamic exponents. By developing the multifractal model we show that dynamic multiscaling occurs in passive-scalar turbulence only if the advecting velocity field is itself multifractal. We substantiate our results by detailed numerical simulations in shell models of passive-scalar advection.

We study turbulence in the one-dimensional Burgers equation with a white-in-time, Gaussian random force that has a Fourier-space spectrum approximately 1/k, where k is the wave number. From very high-resolution numerical simulations, in the limit of vanishing viscosity, we find evidence for multiscaling of velocity structure functions which cannot be falsified by standard tests. We find a new artifact in which logarithmic corrections can appear disguised as anomalous scaling and conclude that bifractal scaling is likely.

We show that different ways of extracting time scales from time-dependent velocity structure functions lead to different dynamic-multiscaling exponents in fluid turbulence. These exponents are related to equal-time multiscaling exponents by different classes of bridge relations, which we derive. We check this explicitly by detailed numerical simulations of the Gledzer-Ohkitani-Yamada shell model for fluid turbulence. Our results can be generalized to any system in which both equal-time and time-dependent structure functions show multiscaling.

Active turbulent advection is considered in the context of magnetohydrodynamics. In this case, an auxiliary passive field bears no apparent connection to the active field. The scaling properties of the two fields are different. In the framework of a shell model, we show that the two-point structure function of the passive field has a unique zero mode, characterizing the scaling of this field only. In other words, the existence of statistical invariants for the decaying passive field carries no information on the scaling properties of the active field.