Search

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 (J. Fluid Mech., vol. 3, 1957, pp. 185-204) and rippling (e.g. Young & Wolfe, J. Fluid Mech., vol. 739,2014, pp. 276-307) 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. [ABSTRACT FROM AUTHOR]

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, ρa/ρw ≡ ∈ ≪ 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 c0 of surface waves. At order ∈, the eigenvalue has a finite imaginary part, which implies growth. Miles (J. Fluid Mech., vol. 3, 1957, pp. 185-204) showed that the growth rate is proportional to the square modulus of the leading-order eigenfunction evaluated at the so-called critical level z = zc, where the wind speed is equal to c0 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. [ABSTRACT FROM AUTHOR]

Aim: We present evidence for finite magnetic helicity density in the heliosphere and numerical models thereof, and relate it to the magnetic field properties of the dynamo in the solar convection zone. Methods: We use simulations and solar wind data to compute magnetic helicity either directly from the simulations or indirectly using time series of the skew-symmetric components of the magnetic correlation tensor. Results: We find that the solar dynamo produces negative magnetic helicity at small scales and positive at large scales. However, in the heliosphere these properties are reversed and the magnetic helicity is now positive at small scales and negative at large scales. We explain this by the fact that a negative diffusive magnetic helicity flux corresponds to a positive gradient of magnetic helicity, which leads to a change of sign from negative to positive values at some radius in the northern hemisphere.

The rate of collision and the relative velocities of the colliding particles in turbulent flows are a crucial part of several natural phenomena, e.g. rain formation in warm clouds and planetesimal formation in protoplanetary discs. The particles are often modelled as passive, but heavy and inertial. Within this model, large relative velocities emerge due to formation of singularities (caustics) of 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 (St). The best approximation to our data is J ~exp(-C/St), in the limit St→0 where C is a non-universal constant. [ABSTRACT FROM AUTHOR]

As expected, we observe Kolmogorov scaling Graph $E(k)\sim k^{-5/3}$ for Graph $k and a pseudo-turbulence scaling Graph $E(k) \sim k^{-3}$ for Graph $k>k {d}$. Hence for both of these cases we obtain Graph $E(k) \sim k^{-3}$ for Graph $k > k {d}$ and Graph $E(k)\sim k^{-5/3}$ for Graph $k. In figure 13 we plot the p.d.f. of Graph ${\rm \pi} {K}(\boldsymbol {x})$ with Graph $K=k {d}$ for Graph $b=0$, Graph $b=0.13$ and Graph $b=\infty$. For all the runs, Graph $L=2{\rm \pi}$ and Graph $d=0.46$, and the dimensionless numbers Graph $Ga= 302$, Graph $Bo =1.8$, Graph $At = 0.04$ and Graph $\varPhi = 1.64\,\%$ are kept fixed. \end{equation} In figure 3( I a i , I b i ) we plot the velocity streamlines of the average velocity field Graph ${\boldsymbol {u}}^{CM}({\boldsymbol {\xi }})$ for Graph $b=\infty$ (Graph ${\texttt {R0}}$) and Graph $b=0.13$ (Graph ${\texttt {R3}}$). [Extracted from the article]

The effect of turbulence on the mass transfer between a fluid and embedded small heavy inertial particles that experience surface reactions is studied. For simplicity, the surface reaction, which takes place when a gas phase reactant is converted to a gas phase product at the external surface of the particles, is unimolar and isothermal. 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 Damköhler numbers, where the Damköhler number is defined as the ratio of the turbulent and chemical time scales, and it may totally control the mass transfer rate for Damköhler numbers larger than unity. A model that describes how this effect should be incorporated into existing simulation tools that utilize the Reynolds averaged Navier-Stokes approach is presented. [ABSTRACT FROM AUTHOR]

We present an overview of the statistical properties of turbulence in two-dimensional (2D) fluids. After a brief recapitulation of well-known results for statistically homogeneous and isotropic 2D fluid turbulence, we give an overview of recent progress in this field for such 2D turbulence in conducting fluids, fluids with polymer additives, binary-fluid mixtures, and superfluids; we also discuss the statistical properties of particles advected by 2D turbulent fluids. [ABSTRACT FROM AUTHOR]

We study the role of turbulence in magnetic reconnection, within the framework of magnetohydrodynamics, using three-dimensional direct numerical simulations. For small turbulent intensity, we find that the reconnection rate obeys Sweet–Parker scaling. For large enough turbulent intensity, reconnection rate departs significantly from Sweet–Parker behaviour, becomes almost a constant as a function of the Lundquist number.We further study energization of test-particles in the same set-up. We find that the speed of the energized particles obeys a Maxwellian distribution, whose variance also obeys Sweet–Parker scaling for small turbulent intensity but depends weakly on the Lundquist number for large turbulent intensity. Furthermore, the variance is found to increase with the strength of the reconnecting magnetic field. [ABSTRACT FROM AUTHOR]

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öhler numbers. The particle clustering significantly slows down the reaction rate for increasing Damköhler 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. [ABSTRACT FROM AUTHOR]

We consider strongly stratified forced turbulence in a plane-parallel layer with helicity and corresponding large-scale dynamo action in the lower part and non-helical turbulence in the upper. The magnetic field is found to develop strongly concentrated bipolar structures near the surface. They form elongated bands with a sharp interface between opposite polarities. Unlike earlier experiments with imposed magnetic field, the inclusion of rotation does not strongly suppress the formation of these structures. We perform a systematic numerical study of this phenomenon by varying magnetic Reynolds number, scale-separation ratio, and Coriolis number. We focus on the formation of a current sheet between bipolar regions where reconnection of oppositely oriented field lines occurs. We determine the reconnection rate by measuring either the inflow velocity in the vicinity of the current sheet or by measuring the electric field in the reconnection region. We demonstrate that for large Lundquist numbers, S > 10³, the reconnection rate is nearly independent of S in agreement with results of recent numerical simulations performed by other groups in simpler settings. [ABSTRACT FROM AUTHOR]

We perform direct numerical simulations of the equations of magnetohydrodynamics with external random forcing and in the presence of gravity. The domain is divided into two parts: a lower layer where the forcing is helical and an upper layer where the helicity of the forcing is zero with a smooth transition in between. At early times, a large-scale helical dynamo develops in the bottom layer. At later times the dynamo saturates, but the vertical magnetic field continues to develop and rises to form dynamic bipolar structures at the top, which later disappear and reappear. Some of the structures look similar to δ spots observed in the Sun. This is the first example of magnetic flux concentrations, owing to strong density stratification, from self-consistent dynamo simulations that generate bipolar, super-equipartition strength, magnetic structures whose energy density can exceeds the turbulent kinetic energy by even a factor of 10. [ABSTRACT FROM PUBLISHER]

The cosmic magnetic fields in regions of low plasma pressure and large currents, such as in interstellar space and gaseous nebulae, are force-free in the sense that the Lorentz force vanishes. The three-dimensional Arnold-Beltrami-Childress (ABC) field is an example of a force-free, helical magnetic field. In fluid dynamics, ABC flows are steady state solutions of the Euler equation. The ABC magnetic field lines exhibit a complex and varied structure that is a mix of regular and chaotic trajectories in phase space. The characteristic features of field line trajectories are illustrated through the phase space distribution of finite-distance and asymptotic-distance Lyapunov exponents. In regions of chaotic trajectories, an ensemble-averaged variance of the distance between field lines reveals anomalous diffusion--in fact, superdiffusion--of the field lines. The motion of charged particles in the force-free ABC magnetic fields is different from the flow of passive scalars in ABC flows. The particles do not necessarily follow the field lines and display a variety of dynamical behavior depending on their energy, and their initial pitch-angle. There is an overlap, in space, of the regions in which the field lines and the particle orbits are chaotic. The time evolution of an ensemble of particles, in such regions, can be divided into three categories. For short times, the motion of the particles is essentially ballistic; the ensemble-averaged, mean square displacement is approximately proportional to t², where t is the time of evolution. The intermediate time region is defined by a decay of the velocity autocorrelation function--this being a measure of the time after which the collective dynamics is independent of the initial conditions. For longer times, the particles undergo superdiffusion--the mean square displacement is proportional to tα, where α > 1, and is weakly dependent on the energy of the particles. These super-diffusive characteristics, both of magnetic field lines and of particles moving in these fields, strongly suggest that theories of transport in three-dimensional chaotic magnetic fields need a shift from the usual paradigm of quasilinear diffusion. [ABSTRACT FROM AUTHOR]

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]

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. Advected by the turbulent disk flow, the dust grains collide with the boulder and we compute the probability density function 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 timescales. [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]

The article presents a study which explores the nonmonotonic behavior of elastic range scaling in polymeric turbulence at large Reynolds and Deborah numbers. It mentions the use of direct numerical simulations to study homogeneous and isotropic turbulent flows of dilute polymer solutions at high Reynolds and Deborah numbers.

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]

ABSTRACT We consider mean-field dynamo models with fluctuating α effect, both with and without large-scale shear. The α effect is chosen to be Gaussian white noise with zero mean and a given covariance. In the presence of shear, we show analytically that (in infinitely large domains) the mean-squared magnetic field shows exponential growth. The growth rate of the fastest growing mode is proportional to the shear rate. This result agrees with earlier numerical results of Yousef et al. and the recent analytical treatment by Heinemann, McWilliams & Schekochihin who use a method different from ours. In the absence of shear, an incoherent α2 dynamo may also be possible. We further show by explicit calculation of the growth rate of third- and fourth-order moments of the magnetic field that the probability density function of the mean magnetic field generated by this dynamo is non-Gaussian. [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]

The chiral symmetry breaking properties of the Tayler instability are discussed. Effective amplitude equations are determined in one case. This model has three free parameters that are determined numerically. Comparison with chiral symmetry breaking in biochemistry is made. [ABSTRACT FROM PUBLISHER]

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]

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e., when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential. [ABSTRACT FROM AUTHOR]

We present a unified three-dimensional model of the convection zone and upper atmosphere of the Sun in spherical geometry. In this model, magnetic fields, generated by a helically forced dynamo in the convection zone, emerge without the assistance of magnetic buoyancy. We use an isothermal equation of state with gravity and density stratification. Recurrent plasmoid ejections, which rise through the outer atmosphere, is observed. In addition, the current helicity of the small–scale field is transported outwards and form large structures like magnetic clouds. [ABSTRACT FROM PUBLISHER]

We summarise recent results form direct numerical simulations of both non-rotating helically forced and rotating convection driven MHD equations in spherical wedge-shape domains. In the former, using perfect-conductor boundary conditions along the latitudinal boundaries we observe oscillations, polarity reversals and equatorward migration of the large-scale magnetic fields. In the latter we obtain angular velocity with cylindrical contours and large-scale magnetic field which shows oscillations, polarity reversals but poleward migration. The occurrence of these behviours in direct numerical simulations is clearly of interest. However the present models as they stand are not directly applicable to the solar dynamo problem. Nevertheless, they provide general insights into the operation of turbulent dynamos. [ABSTRACT FROM PUBLISHER]

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]