Search

When a Bose-Einstein condensate (BEC) is driven out of equilibrium, density waves interact non-linearly and trigger turbulent cascades. In a turbulent BEC, energy is transferred towards small scales by a direct cascade, whereas the number of particles displays an inverse cascade toward large scales. In this work, we study analytically and numerically the direct and inverse cascades in wave-turbulent BECs. We analytically derive the Kolmogorov-Zakharov spectra, including the log-correction to the direct cascade scaling and the universal pre-factor constants for both cascades. We test and corroborate our predictions using high-resolution numerical simulations of the forced-dissipated Gross-Pitaevskii model in a periodic box and the corresponding wave-kinetic equation. Theoretical predictions and data are in excellent agreement, without adjustable parameters. Moreover, in order to connect with experiments, we test and validate our theoretical predictions using the Gross-Pitaevskii model with a confining cubic trap. Our results explain previous experimental observations and suggest new settings for future studies.

We study intermittency of circulation moments in turbulent superfluid helium by using experimental grid turbulence and numerical simulations of the Hall-Vinen-Bekarevich-Khalatnikov model. More precisely, we compute the velocity circulation $\Gamma_r$ in loops of size $r$ laying in the inertial range. For both, experimental and numerical data, the circulation variance shows a clear Kolmogorov scaling $\langle \Gamma_r^2 \rangle \sim r^{8/3}$ in the inertial range, independently of the temperature. Scaling exponents of high-order moments are comparable, within error bars, to previously reported anomalous circulation exponents in classical turbulence and low-temperature quantum turbulence numerical simulations., Comment: 12 pages, 10 figures

Turbulence is unique in its appeal across physics, mathematics and engineering. And yet a microscopic theory, starting from the basic equations of hydrodynamics, still eludes us. In the last decade or so, new directions at the interface of physics and mathematics have emerged, which strengthens the hope of ‘solving’ one of the oldest problems in the natural sciences. This two-part theme issue unites these new directions on a common platform emphasizing the underlying complementarity of the physicists’ and the mathematicians’ approaches to a remarkably challenging problem. This article is part of the theme issue ‘Scaling the turbulence edifice (part 1)’.

We study numerically the process of vortex nucleation at the wake of a moving object in superfluids using a generalized and non-local Gross-Pitaevskii model. The non-local potential is set to reproduce the roton minimum present in the excitation spectrum of superfluid helium. By applying numerically a Newton-Raphson method we determine the bifurcation diagram for different types of non-linearities and object sizes which allow for determining the corresponding critical velocities. In the case of a non-local potential, we observe that for small object sizes the critical velocity is simply determined by the Landau criterion for superfluidity whereas for large objects there is little difference between all models studied. Finally, we study dynamically in two and three dimensions how rotons and vortices are excited in the non-local model of superfluid., 9 pages, 8 figures

We test the predictions of the theory of weak wave turbulence by performing numerical simulations of the Gross-Pitaevskii equation (GPE) and the associated wave-kinetic equation (WKE). We consider an initial state localized in Fourier space, and we confront the solutions of the WKE obtained numerically with GPE data for both the wave-action spectrum and the probability density functions (PDFs) of the Fourier mode intensities. We find that the temporal evolution of the GPE data is accurately predicted by the WKE, with no adjustable parameters, for about two nonlinear kinetic times. Qualitative agreement between the GPE and the WKE persists also for longer times with some quantitative deviations that may be attributed to the combination of breakdown of the theoretical assumptions underlying the WKE as well as numerical issues. Furthermore, we study how the wave statistics evolves toward Gaussianity in a time scale of the order of the kinetic time.The excellent agreement between direct numerical simulations of the GPE and the WKE provides a new and solid ground to the theory of weak wave turbulence.

Experiments in the early 1980s have shown that a compact bundle of quantum vortex rings in superfluid helium remains coherent and travels a significant distance compared to its size. This is surprising because a single vortex ring, under the effect of friction with the background of thermal excitations (the so-called normal fluid), would quickly lose energy and decay. The observation of these long-lived vortex structures has remained an unsettled question since their experimental detection. In this work, by taking into account the fully coupled dynamics of superfluid vortices and the normal fluid, we show that the motion of the superfluid bundle generates a large-scale wake in the normal fluid which reduces the overall friction experienced by the bundle, enhancing its lifetime. This collective effect is similar to the drag reduction observed in systems of active, hydrodynamically cooperative agents such as bacteria in aqueous suspensions, fungal spores in the atmosphere and cyclists in pelotons., Comment: 10 pages, 4 figures

The velocity circulation, a measure of the rotation of a fluid within a closed path, is a fundamental observable in classical and quantum flows. It is indeed a Lagrangian invariant in inviscid classical fluids. In quantum flows, circulation is quantized, taking discrete values that are directly related to the number and the orientation of thin vortex filaments enclosed by the path. By varying the size of such closed loop, the circulation provides a measure of the dependence of the flow structure on the considered scale. Here we consider the scale dependence of circulation statistics in quantum turbulence, using high resolution direct numerical simulations of a generalized Gross-Pitaevskii model. Results are compared to the circulation statistics obtained from simulations of the incompressible Navier-Stokes equations. When the integration path is smaller than the mean inter-vortex distance, the statistics of circulation in quantum turbulence displays extreme intermittent behavior due to the quantization of circulation, in stark contrast with the viscous scales of classical flows. In contrast, at larger scales, circulation moments display striking similarities with the statistics probed in the inertial range of classical turbulence. This includes the emergence of the power law scalings predicted from Kolmogorov's 1941 theory, as well as intermittency deviations that closely follow the recently proposed bifractal model for circulation moments in classical flows. To date, this is the most convincing evidence of intermittency in the large scales of quantum turbulence. Moreover, our results strongly reinforce the resemblance between classical and quantum turbulence, highlighting the universality of inertial range dynamics, including intermittency, across these two a priori very different systems., 12 pages, 7 figures; accepted in Phys. Rev. X

The dynamics of an active, finite-size and immiscible impurity in a dilute quantum fluid at finite temperature is characterized by means of numerical simulations of the projected Gross--Pitaevskii equation. The impurity is modeled as a localized repulsive potential and described with classical degrees of freedom. It is shown that impurities of different sizes thermalize with the fluid and undergo a stochastic dynamics compatible with an Ornstein--Uhlenbeck process at sufficiently large time-lags. The velocity correlation function and the displacement of the impurity are measured and an increment of the friction with temperature is observed. Such behavior is phenomenologically explained in a scenario where the impurity exchanges momentum with a dilute gas of thermal excitations, experiencing an Epstein drag., 10 pages, 7 figures

International audience; Particles are a widespread tool for obtaining information from fluid flows. When Eulerian data are unavailable, they may be employed to estimate flow fields or to identify coherent flow structures. Here we numerically examine the possibility of using particles to capture the dynamics of isolated vortex rings propagating in a quiescent fluid. The analysis is performed starting from numerical simulations of the Navier–Stokes and the Hall–Vinen–Bekarevich–Khalatnikov equations, respectively describing the dynamics of a Newtonian fluid and a finite-temperature superfluid. The flow-induced positions and velocities of particles suspended in the fluid are specifically used to compute the Lagrangian pseudovorticity field, a proxy for the Eulerian vorticity field recently employed in the context of superfluid $^{4}\textrm {He}$ experiments. We show that, when calculated from ideal Lagrangian tracers or from particles with low inertia, the pseudovorticity field can be accurately used to estimate the propagation velocity and the growth of isolated vortex rings, although the quantitative reconstruction of the corresponding vorticity fields remains challenging. On the other hand, particles with high inertia tend to preferentially sample specific flow regions, resulting in biased pseudovorticity fields that pollute the estimation of the vortex ring properties. Overall, this work neatly demonstrates that the Lagrangian pseudovorticity is a valuable tool for estimating the strength of macroscopic vortical structures in the absence of Eulerian data, which is, for example, the case for superfluid $^{4}\textrm {He}$ experiments.

The understanding of turbulent flows is one of the biggest current challenges in physics, as no first-principles theory exists to explain their observed spatio-temporal intermittency. Turbulent flows may be regarded as an intricate collection of mutually-interacting vortices. This picture becomes accurate in quantum turbulence, which is built on tangles of discrete vortex filaments. Here, we study the statistics of velocity circulation in quantum and classical turbulence. We show that, in quantum flows, Kolmogorov turbulence emerges from the correlation of vortex orientations, while deviations -- associated with intermittency -- originate from their non-trivial spatial arrangement. We then link the spatial distribution of vortices in quantum turbulence to the coarse-grained energy dissipation in classical turbulence, enabling the application of existent models of classical turbulence intermittency to the quantum case. Our results provide a connection between the intermittency of quantum and classical turbulence and initiate a promising path to a better understanding of the latter., Comment: 10 pages, 6 figures + supplementary information. Accepted for publication in Nature Communications

In a recent paper, Tanogami [Phys. Rev. E 103, 023106 (2021)2470-004510.1103/PhysRevE.103.023106] proposes a scenario for quantum turbulence where the energy spectrum at scales smaller than the intervortex distance is dominated by a quantum stress cascade, in opposition to Kelvin-wave cascade predictions. The purpose of the present Comment is to highlight some physical issues in the derivation of the quantum stress cascade, in particular to stress that quantization of circulation has been ignored.

We report an exact unique constant-flux power-law analytical solution of the wave kinetic equation for the turbulent energy spectrum, $E(k)=C_1 \sqrt{\varepsilon\, a c_{\rm s} }/k$, of acoustic waves in 2D with almost linear dispersion law, $\omega_k = c_{\rm s} k[1+(ak)^2]$, $ ak \ll 1$. Here $\varepsilon$ is the energy flux over scales, and $C_1$ is the universal constant which was found analytically. Our theory describes, for example, acoustic turbulence in 2D Bose-Einstein condensates (BECs). The corresponding 3D counterpart of turbulent acoustic spectrum was found over half a century ago, however, due to the singularity in 2D, no solution has been obtained until now. We show the spectrum $E(k)$ is realizable in direct numerical simulations of forced-dissipated Gross-Pitaevskii equation in the presence of strong condensate.

Finite-temperature quantum turbulence is often described in terms of two immiscible fluids that can flow with a non-zero mean relative velocity. Such out-of-equilibrium state is known as counterflow superfluid turbulence. We report here the emergence of a counterflow-induced inverse energy cascade in three-dimensional superfluid flows by performing extensive numerical simulations of the Hall-Vinen-Bekarevich-Khalatnikov model. As the intensity of the mean counterflow is increased, an abrupt transition, from a fully three-dimensional turbulent flow to a quasi two-dimensional system exhibiting a split cascade, is observed. The findings of this work could motivate new experimental settings to study quasi two-dimensional superfluid turbulence in the bulk of three-dimensional experiments. They might also find applications beyond superfluids in systems described by more than one fluid component., 5 pages, 4 figures

We statistically study vortex reconnections in quantum fluids by evolving different realizations of vortex Hopf links using the Gross--Pitaevskii model. Despite the time-reversibility of the model, we report a clear evidence that the dynamics of the reconnection process is time-irreversible, as reconnecting vortices tend to separate faster than they approach. Thanks to a matching theory devised concurrently in Proment and Krstulovic (arXiv:2005.02047), we quantitatively relate the origin of this asymmetry to the generation of a sound pulse after the reconnection event. Our results have the prospect of being tested in several quantum fluid experiments and, theoretically, may shed new light on the energy transfer mechanisms in both classical and quantum turbulent fluids., Comment: 5 pages, 4 figures

Reconnections between quantum vortex filaments in presence of trapped particles are investigated using numerical simulations of the Gross--Pitaevskii equation. Particles are described with classical degrees of freedom and modeled as highly repulsive potentials which deplete the superfluid. First, the case of a vortex dipole with a single particle trapped inside one of the vortices is studied. It is shown that the reconnection takes place at the position of the particle as a consequence of the symmetry breaking induced by it. The separation rate between the reconnecting points is compatible with the known dynamics of quantum vortex reconnections and it is independent of the particle mass and size. After the reconnection, the particle is pushed away with a constant velocity and its trajectory is deflected because of the transverse momentum exchange with the vortex filaments. The momentum exchanges between the particle, the vortex, and a density pulse are characterized. Finally, the reconnection of two linked rings, each of them with several initially randomly distributed particles is studied. It is observed that generically, reconnections take place at the location of trapped particles. It is shown that reconnection dynamics is unaffected for light particles., 9 pages, 7 figures

The evolution of a turbulent tangle of quantum vortices in presence of finite-size active particles is studied by means of numerical simulations of the Gross-Pitaevskii equation. Particles are modeled as potentials depleting the superfluid and described with classical degrees of freedom following a Newtonian dynamics. It is shown that particles do not modify the building-up and the decay of the superfluid Kolmogorov turbulent regime. It is observed that almost the totality of particles remains trapped inside quantum vortices, although they are occasionally detached and recaptured. The statistics of this process are presented and discussed. The particle Lagrangian dynamics is also studied. At large time scales, the velocity spectrum of particles is reminiscent of a classical Lagrangian turbulent behavior. At time-scales faster than the turnover time associated to the mean inter-vortex distance, the particle motion is dominated by oscillations due to Magnus effect. For light particles a non-classical scaling of the spectrum arises. The particle velocity and acceleration probability distribution functions are then studied. The decorrelation time of the particle acceleration is found to be shorter than in classical fluids, and related to the Magnus force experienced by the trapped particles., 17 pages, 11 figures

We statistically study of vortex reconnection in quantum fluids by evolving different realizations of vortex Hopf links using the Gross--Pitaevskii model. Despite the time-reversibility of the model, we report clear evidence that the dynamics of the reconnection process is time-irreversible, as reconnecting vortices tend to separate faster than they approach. Thanks to a matching theory devised concurrently in Proment and Krstulovic (2020), we quantitatively relate the origin of this asymmetry to the generation of a sound pulse after the reconnection event. Our results have the prospect of being tested in several quantum fluid experiments and, theoretically, may shed new light on the energy transfer mechanisms in both classical and quantum turbulent fluids.

We present the Fully cOUpled loCAl model of sUperfLuid Turbulence (FOUCAULT) that describes the dynamics of finite temperature superfluids. The superfluid component is described by the vortex filament method while the normal fluid is governed by a modified Navier-Stokes equation. The superfluid vortex lines and normal fluid components are fully coupled in a self-consistent manner by the friction force, which induces local disturbances in the normal fluid in the vicinity of vortex lines. The main focus of this work is the numerical scheme for distributing the friction force to the mesh points where the normal fluid is defined and for evaluating the velocity of the normalfluid on the Lagrangian discretization points along the vortex lines. In particular, we show that if this numerical scheme is not careful enough, spurious results may occur. The new scheme which we propose to overcome these difficulties is based on physical principles. Finally, we apply the new method to the problem of the motion of a superfluid vortex ring in a stationary normal fluid and in a turbulent normal fluid., 19 pages, 9 figures

In a concurrent work, Villois et al. 2020 reported the evidence that vortex reconnections in quantum fluids follow an irreversible dynamics, namely vortices separate faster than they approach; such time-asymmetry is explained by using simple conservation arguments. In this work we develop further these theoretical considerations and provide a detailed study of the vortex reconnection process for all the possible geometrical configurations of the order parameter (superfluid) wave function. By matching the theoretical description of incompressible vortex filaments and the linear theory describing locally vortex reconnections, we determine quantitatively the linear momentum and energy exchanges between the incompressible (vortices) and the compressible (density waves) degrees of freedom of the superfluid. We show theoretically and corroborate numerically, why a unidirectional density pulse must be generated after the reconnection process and why only certain reconnecting angles, related to the rates of approach and separations, are allowed. Finally, some aspects concerning the conservation of centre-line helicity during the reconnection process are discussed., Comment: 17 pages, 12 figures

Phase transitions of a finite-size two-dimensional quantum fluid of bosons in presence of impurities are studied by using the truncated Gross-Pitaevskii model. Impurities are described with classical degrees of freedom. A spontaneous clustering of impurities during the thermalization is observed. Depending on the interaction among impurities, such clusters can break due to thermal fluctuations at temperatures where the condensed fraction is still significant. Clustering is found to increase the condensation transition temperature. The condensation and the Berezinskii-Kosterlitz-Thouless transition temperatures are determined numerically and found to strongly depend on the volume occupied by the impurities. A relative increase up to a 20% of their respective values is observed, whereas their ratio remains approximately constant., Comment: 13 pages, 5 figures

International audience; In fluid turbulence, energy is transferred from one scale to another by an energy cascade that depends only on the energy-dissipation rate. It leads by dimensional arguments to the Kolmogorov 1941 (K41) spectrum. Here we show that the normal modes of vibrations in elastic plates also manifest an energy cascade with the same K41 spectrum in the fully nonlinear regime. In particular, we observe different patterns in the elastic deformations such as folds, developable cones, and even more complex stretching structures, in analogy with spots, swirls, vortices, and other structures in hydrodynamic turbulence. We show that the energy cascade is dominated by the kinetic contribution and that the stretching energy is at thermodynamical equilibrium. We characterize this energy cascade, the validity of the constant energy-dissipation rate over the scales. Finally, we discuss the role of intermittency using the correlation functions that exhibit anomalous exponents.

We study freely decaying quantum turbulence by performing high resolution numerical simulations of the Gross-Pitaevskii equation (GPE) in the Taylor-Green geometry. We use resolutions ranging from $1024^3$ to $4096^3$ grid points. The energy spectrum confirms the presence of both a Kolmogorov scaling range for scales larger than the intervortex scale $\ell$, and a second inertial range for scales smaller than $\ell$. Vortex line visualizations show the existence of substructures formed by a myriad of small-scale knotted vortices. Next, we study finite temperature effects in the decay of quantum turbulence by using the stochastic Ginzburg-Landau equation to generate thermal states, and then by evolving a combination of these thermal states with the Taylor-Green initial conditions using the GPE. We extract the mean free path out of these simulations by measuring the spectral broadening in the Bogoliubov dispersion relation obtained from spatio-temporal spectra, and use it to quantify the effective viscosity as a function of the temperature. Finally, in order to compare the decay of high temperature quantum and that of classical flows, and to further calibrate the estimations of viscosity from the mean free path in the GPE simulations, we perform low Reynolds number simulations of the Navier-Stokes equations., Comment: 15 figures, Supplemental material: Video M1 and M2

Motivated by systems in which droplets grow and shrink in a turbulence-driven supersaturation field, we investigate the problem of turbulent condensation in a general manner. Using direct numerical simulations we show that the turbulent fluctuations of the supersaturation field offer different conditions for the growth of droplets which evolve in time due to turbulent transport and mixing. Based on that, we propose a Lagrangian stochastic model for condensation and evaporation of small droplets in turbulent flows. It consists of a set of stochastic integro-differential equations for the joint evolution of the squared radius and the supersaturation along the droplet trajectories. The model has two parameters fixed by the total amount of water and the thermodynamic properties, as well as the Lagrangian integral timescale of the turbulent supersaturation. The model reproduces very well the droplet size distributions obtained from direct numerical simulations and their time evolution. A noticeable result is that, after a stage where the squared radius simply diffuses, the system converges exponentially fast to a statistical steady state independent of the initial conditions. The main mechanism involved in this convergence is a loss of memory induced by a significant number of droplets undergoing a complete evaporation before growing again. The statistical steady state is characterised by an exponential tail in the droplet mass distribution. These results reconcile those of earlier numerical studies, once these various regimes are considered., 24 pages, 12 figures

Modeling dispersed solid phases in fluids still represents a computational challenge when considering a small-scale coupling in wide systems, such as the atmosphere or industrial processes at high Reynolds numbers. A numerical method is here introduced for simulating the dynamics of diffusive heavy inertial particles in turbulent flows. The approach is based on the position/velocity phase–space particle distribution. The discretization of velocities is inspired from lattice Boltzmann methods and is chosen to match discrete displacements between two time steps. For each spatial position, the time evolution of particles momentum is approximated by a finite-volume approach. The proposed method is tested for particles experiencing a Stokes viscous drag with a prescribed fluid velocity field in one dimension using a random flow, and in two dimensions with the solution to the forced incompressible Navier–Stokes equations. Results show good agreements between Lagrangian and Eulerian dynamics for both spatial clustering and the dispersion in particle velocities. In particular, the proposed method, in contrast to hydrodynamical Eulerian descriptions of the dispersed phase, is able to reproduce fine particle kinetic phenomena, such as caustic formation or trajectory crossings. This indicates the suitability of this approach at large Stokes numbers for situations where details of collision processes are important.

Particles are today the main tool to study superfluid turbulence and visualize quantum vortices. In this work, we study the dynamics and the spatial distribution of particles in co-flow and counterflow superfluid helium turbulence in the framework of the two-fluid Hall-Vinen-Bekarevich-Khalatnikov (HVBK) model. We perform three-dimensional numerical simulations of the HVBK equations along with the particle dynamics that depends on the motion of both fluid components. We find that, at low temperatures, where the superfluid mass fraction dominates, particles strongly cluster in vortex filaments regardless of their physical properties. At higher temperatures, as viscous drag becomes important and the two components become tightly coupled, the clustering dynamics in the coflowing case approach those found in classical turbulence, while under strong counterflow, the particle distribution is dominated by the quasi-two-dimensionalization of the flow., Comment: 8 pages, 4 figures, accepted in Physical Review Fluids

Weak wave turbulence has been observed on a thin elastic plates since the work by During et al. [Phys. Rev. Lett. 97, 025503 (2006)PRLTAO0031-900710.1103/PhysRevLett.97.025503]. Here we report theoretical, experimental, and numerical studies of wave turbulence in a forced thin elastic plate submitted to increasing tension. When increasing the tension (or decreasing the bending stiffness of the plate) the plate evolves progressively from a plate into an elastic membrane as in drums. We first consider a thin plate and increase the tension in experiments and numerical simulations. We observe that the system remains in a state of weak turbulence of weakly dispersive waves. This observation is in contrast with what has been observed in water waves when decreasing the water depth, which also changes the waves from dispersive to weakly dispersive. The weak turbulence observed in the deep water case evolves into a solitonic regime. Here no such transition is observed for the stretched plate. We then apply the weak turbulence theory to the membrane case and show with numerical simulations that indeed the weak turbulence framework remains valid for the membrane and no formation of singular structures (shocks) should be expected in contrast with acoustic wave turbulence.

This paper investigates the effect of inertia on the dynamics of elongated chains to go beyond the overdamped case that is often used to study such systems. For that purpose, numerical simulations are performed considering the motion of freely jointed bead-rod chains in an extensional flow in the presence of thermal noise. The coil-stretch transition and the tumbling instability are characterized as a function of three parameters: the Peclet number, the Stokes number, and the chain length. Numerical results show that the coil-stretch transition remains when inertia is present and that it depends nonlinearly on the Stokes and Peclet numbers. Theoretical and numerical analyses also highlight the role of intermediate stable configurations in the dynamics of elongated chains: chains can indeed remain trapped for a certain time in these configurations, especially while undergoing a tumbling event.

One of the main features of superfluids is the presence of topological defects with quantised circulation. These objects are known as quantum vortices and exhibit a hydrodynamic behaviour. Nowadays, particles are the main experimental tool used to visualise quantum vortices and to study their dynamics. We use a self-consistent model based on the three-dimensional Gross-Pitaevskii (GP) equation to explore theoretically and numerically the attractive interaction between particles and quantised vortices at very low temperature. Particles are described as localised potentials depleting the superfluid and following Newtonian dynamics. We are able to derive analytically a reduced central-force model that only depends on the classical degrees of freedom of the particle. Such model is found to be consistent with the GP simulations. We then generalised the model to include deformations of the vortex filament. The resulting long-range mutual interaction qualitatively reproduces the observed generation of a cusp on the vortex filament during the particle approach. Moreover, we show that particles can excite Kelvin waves on the vortex filament through a resonance mechanism even if they are still far from it., Comment: 13 pages, 5 figures

Insight into vortex reconnections in superfluids is presented, making use of analytical results and numerical simulations of the Gross-Pitaevskii model. Universal aspects of the reconnection process are investigated by considering different initial vortex configurations and making use of a recently developed tracking algorithm to reconstruct the vortex filaments. We show that during a reconnection event the vortex lines approach and separate always according to the time scaling $ \delta \sim t^{1/2} $ with prefactors that depend on the vortex configuration. We also investigate the behavior of curvature and torsion close to the reconnection point, demonstrating analytically that the curvature can exhibit a self-similar behavior that might be broken by the development of shocklike structures in the torsion.

An exact result concerning the energy transfers between non-linear waves of thin elastic plate is derived. Following Kolmogorov's original ideas in hydrodynamical turbulence, but applied to the F\"oppl-von K\'arm\'an equation for thin plates, the corresponding K\'arm\'an-Howarth-Monin relation and an equivalent of the $\frac{4}{5}$-Kolmogorov's law is derived. A third-order structure function involving increments of the amplitude, velocity and the Airy stress function of a plate, is proven to be equal to $-\varepsilon\, \ell$, where $\ell$ is a length scale in the inertial range at which the increments are evaluated and $\varepsilon$ the energy dissipation rate. Numerical data confirm this law. In addition, a useful definition of the energy fluxes in Fourier space is introduced and proven numerically to be flat in the inertial range. The exact results derived in this Letter are valid for both, weak and strong wave-turbulence. They could be used as a theoretical benchmark of new wave-turbulence theories and to develop further analogies with hydrodynamical turbulence.

A dilute system of reacting particles transported by fluid flows is considered. The particles react as A+A→∅ with a given rate when they are within a finite radius of interaction. The system is described in terms of the joint n-point number spatial density that it is shown to obey a hierarchy of transport equations. An analytic solution is obtained in the dilute or, which is equivalent, the long-time limit by using a Lagrangian approach where statistical averages are performed along non-reacting trajectories. In this limit, it is shown that the moments of the number of particles have an exponential decay rather than the algebraic prediction of standard mean-field approaches. The effective reaction rate is then related to Lagrangian pair statistics by a large-deviation principle. A phenomenological model is introduced to study the qualitative behavior of the effective rate as a function of the interaction length, the degree of chaoticity of the dynamics and the compressibility of the carrier flow. Exact computations, obtained via a Feynman–Kac approach, in a smooth, compressible, random delta-correlated-in-time Gaussian velocity field support the proposed heuristic approach.

Low-temperature grid generated turbulence is investigated by using numerical simulations of the Gross-Pitaevskii equation. The statistics of regularized velocity increments are studied. Increments of the incompressible velocity are found to be skewed for turbulent states. Results are later confronted with the (quasi) homogeneous and isotropic Taylor-Green flow, revealing the universality of the statistics. For this flow, the statistics are found to be intermittent and a Kolmogorov constant close to the one of classical fluid is found for the second order structure function.

The development and decay of a turbulent vortex tangle driven by the Gross-Pitaevskii equation is studied. Using a recently-developed accurate and robust tracking algorithm, all quantised vortices are extracted from the fields. The Vinen's decay law for the total vortex length with a coefficient that is in quantitative agreement with the values measured in Helium II is observed. The topology of the tangle is then studied showing that linked rings may appear during the decay. The tracking also allows for determining the statistics of small-scales quantities of vortex lines, exhibiting large fluctuations of curvature and torsion. Finally, the temporal evolution of the Kelvin wave spectrum is obtained providing evidence of the development of a weak-wave turbulence cascade., Comment: 8 pages, 8 figures

A short review is given of recent papers on the relaxation to (incompressible) absolute equilibrium. A new algorithm to construct absolute equilibrium of spectrally truncated compressible flows is described. The algorithm uses stochastic processes based on the Clebsch representation of the velocity field to generate density and velocity fields that follow by construction the absolute equilibrium stationary probability. The new method is shown to reproduce the well-known Gaussian results in the incompressible limit. The irrotational compressible absolute equilibrium case is characterized and the distribution is shown to be non-Gaussian. The high-temperature compressible spectra are found not to obey k2 scaling. Finally, oscillating behavior in constant-pressure variable-temperature relaxation is obtained, suggesting the presence of second sound.

We build on recent developments in the study of fluid turbulence [Gibbon \textit{et al.} Nonlinearity 27, 2605 (2014)] to define suitably scaled, order-$m$ moments, $D_m^{\pm}$, of $\omega^\pm= \omega \pm j$, where $\omega$ and $j$ are, respectively, the vorticity and current density in three-dimensional magnetohydrodynamics (MHD). We show by mathematical analysis, for unit magnetic Prandtl number $P_M$, how these moments can be used to identify three possible regimes for solutions of the MHD equations; these regimes are specified by inequalities for $D_m^{\pm}$ and $D_1^{\pm}$. We then compare our mathematical results with those from our direct numerical simulations (DNSs) and thus demonstrate that 3D MHD turbulence is like its fluid-turbulence counterpart insofar as all solutions, which we have investigated, remain in \textit{only one of these regimes}; this regime has depleted nonlinearity. We examine the implications of our results for the exponents $q^{\pm}$ that characterize the power-law dependences of the energy spectra $\mathcal{E}^{\pm}(k)$ on the wave number $k$, in the inertial range of scales. We also comment on (a) the generalization of our results to the case $P_M \neq 1$ and (b) the relation between $D_m^{\pm}$ and the order-$m$ moments of gradients of hydrodynamic fields, which are used in characterizing intermittency in turbulent flows., Comment: 14 pages, 3 figures

We study the front shape of avalanches produced in a rotating two-dimensional drum partially filled with small glass beads. We focus our work on the study of the length and shape of granular fronts, in particular how they do scale in terms of the physical parameters involved. A single scaling law for the length is found. This scaling law is also relevant for the behavior of the full shape of the fronts. More than 300 different fronts shape, for different values of the parameters collapse into an universal curve.

The phenomenology of turbulent relative dispersion is revisited. A heuristic scenario is proposed, in which pairs of tracers undergo a succession of independent ballistic separations during time intervals whose lengths fluctuate. This approach suggests that the logarithm of the distance between tracers self-averages and performs a continuous-time random walk. This leads to specific predictions for the probability distribution of separations, that differ from those obtained using scale-dependent eddy-diffusivity models (e.g. in the framework of Richardson's approach). Such predictions are tested against high-resolution simulations and shed new lights on the explosive separation between tracers., 10 pages, 5 figures

We examine the scaling laws of magnetohydrodynamic (MHD) turbulence for three different types of forcing functions and imposing at all times the fourfold symmetries of the Taylor-Green (TG) vortex generalized to MHD; no uniform magnetic field is present and the magnetic Prandtl number is equal to unity. We also include pumping in the induction equation, and we take the three configurations studied in the decaying case in Lee et al. [Phys. Rev. E 81, 016318 (2010)]. To that effect, we employ direct numerical simulations up to an equivalent resolution of ${2048}^{3}$ grid points. We find that, similarly to the case when the forcing is absent, different spectral indices for the total energy spectrum emerge, corresponding to either a Kolmogorov law, an Iroshnikov-Kraichnan law that arises from the interactions of turbulent eddies and Alfv\'en waves, or to weak turbulence when the large-scale magnetic field is strong. We also examine the inertial range dynamics in terms of the ratios of kinetic to magnetic energy, and of the turnover time to the Alfv\'en time, and analyze the temporal variations of these quasiequilibria.

Heavy inertial particles transported by a turbulent flow are shown to concentrate in the regions where an advected passive scalar, such as temperature, displays very strong front-like discontinuities. This novel effect is responsible for extremely high levels of fluctuations for the passive field sampled by the particles that impacts the heat fluxes exchanged between the particles and the surrounding fluid. Instantaneous and averaged heat fluxes are shown to follow strongly intermittent statistics and anomalous scaling laws., Comment: 4 pages, 5 figures

The evolution of a Taylor-Green forced magnetohydrodynamic (MHD) system showing dynamo activity is analyzed via direct numerical simulations. The statistical properties of the velocity and magnetic field in Eulerian coordinates and along trajectories of fluid elements (Lagrangian coordinates) are found to change between the kinematic, non-linear and saturated regime. Fluid element (tracer) trajectories change from chaotic quasi-isotropic (kinematic phase) to mean magnetic field aligned (saturated phase). The probability density functions (PDFs) of the magnetic field change from strongly non-Gaussian in the kinematic to quasi-Gaussian PDFs in the saturated regime so that their flatness give a precise handle on the definition of the limiting points of the three regimes. Also the statistics of the fluctuations of the kinetic and magnetic energy along fluid trajectories change. All this goes along with a dramatic increase of the correlation time of velocity and magnetic field fluctuations experienced by tracers significantly exceeding one turbulent large-eddy turn-over time. A remarkable consequence is an intermittent scaling regime of the Lagrangian magnetic field structure functions at unusually long time scales.

We investigate the ideal and incompressible magnetohydrodynamic (MHD) equations in three space dimensions for the development of potentially singular structures. The methodology consists in implementing the fourfold symmetries of the Taylor-Green vortex generalized to MHD, leading to substantial computer time and memory savings at a given resolution; we also use a regridding method that allows for lower-resolution runs at early times, with no loss of spectral accuracy. One magnetic configuration is examined at an equivalent resolution of 6144 3 points and three different configurations on grids of 4096 3 points. At the highest resolution, two different current and vorticity sheet systems are found to collide, producing two successive accelerations in the development of small scales. At the latest time, a convergence of magnetic field lines to the location of maximum current is probably leading locally to a strong bending and directional variability of such lines. A novel analytical method, based on sharp analysis inequalities, is used to assess the validity of the finite-time singularity scenario. This method allows one to rule out spurious singularities by evaluating the rate at which the logarithmic decrement of the analyticity-strip method goes to zero. The result is that the finite-time singularity scenario cannot be ruled out, and the singularity time could be somewhere between t = 2.33 and t = 2.70 . More robust conclusions will require higher resolution runs and grid-point interpolation measurements of maximum current and vorticity. Fil: Brachet, M. E.. Centre National de la Recherche Scientifique. Ecole Normale Superieur; Francia Fil: Bustamante, M. D.. University College Dublin; Irlanda Fil: Krstulovic, G.. Observatoire de la Cote D; Francia Fil: Mininni, Pablo Daniel. National Center for Atmospheric Research; Estados Unidos. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; Argentina Fil: Pouquet, A.. National Center for Atmospheric Research; Estados Unidos Fil: Rosenberg, D.. National Center for Atmospheric Research; Estados Unidos

We present an accurate and robust numerical method to track quantized vortex lines in a superfluid described by the Gross-Pitaevskii equation. By utilizing the pseudo-vorticity field of the associated complex scalar order parameter of the superfluid, we are able to track the topological defects of the superfluid and reconstruct the vortex lines which correspond to zeros of the field. Throughout, we assume our field is periodic to allow us to make extensive use of the Fourier representation of the field and its derivatives in order to retain spectral accuracy. We present several case studies to test the precision of the method which include the evaluation of the curvature and torsion of a torus vortex knot, and the measurement of the Kelvin wave spectrum of a vortex line and a vortex ring. The method we present makes no a-priori assumptions on the geometry of the vortices and is therefore applicable to a wide range of systems such as a superfluid in a turbulent state that is characterised by many vortex rings coexisting with sound waves. This allows us to track the positions of the vortex filaments in a dense turbulent vortex tangle and extract statistical information about the distribution of the size of the vortex rings and the inter-vortex separations. In principle, the method can be extended to track similar topological defects arising in other physical systems., 22 pages, 8 figures

We study the statistical properties of the Kelvin waves propagating along quantized superfluid vortices driven by the Gross-Pitaevskii equation. No artificial forcing or dissipation is added. Vortex positions are accurately tracked. This procedure directly allows us to obtain the Kevin-waves occupation-number spectrum. Numerical data obtained from long time integration and ensemble-average over initial conditions supports the spectrum proposed in [L'vov and Nazarenko, JETP Lett 91, 428 (2010)]. Kelvin wave modes in the inertial range are found to be Gaussian as expected by weak-turbulence predictions. Finally the dissipative range of the Kelvin-wave spectrum is studied. Strong non-Gaussian fluctuations are observed in this range.

A mass-ejection model in a time-dependent random environment with both temporal and spatial correlations is introduced. When the environment has a finite correlation length, individual particle trajectories are found to diffuse at large times with a displacement distribution that approaches a Gaussian. The collective dynamics of diffusing particles reaches a statistically stationary state, which is characterized in terms of a fluctuating mass density field. The probability distribution of density is studied numerically for both smooth and non-smooth scale-invariant random environments. Competition between trapping in the regions where the ejection rate of the environment vanishes and mixing due to its temporal dependence leads to large fluctuations of mass. These mechanisms are found to result in the presence of intermediate power-law tails in the probability distribution of the mass density. For spatially differentiable environments, the exponent of the right tail is shown to be universal and equal to −3/2. However, at small values, it is found to depend on the environment. Finally, spatial scaling properties of the mass distribution are investigated. The distribution of the coarse-grained density is shown to possess some rescaling properties that depend on the scale, the amplitude of the ejection rate and the Holder exponent of the environment.

We present a numerical study of the magnetic field generated by the Taylor-Green vortex. We show that periodic boundary conditions can be used to mimic realistic boundary conditions by prescribing the symmetries of the velocity and magnetic fields. This gives insight into some problems of central interest for dynamos: the possible effect of velocity fluctuations on the dynamo threshold, and the role of boundary conditions on the threshold and on the geometry of the magnetic field generated by dynamo action. In particular, we show that an axial dipolar dynamo similar to the one observed in a recent experiment can be obtained with an appropriate choice of the symmetries of the magnetic field. The nonlinear saturation is studied and a simple model explaining the magnetic Prandtl number dependence of the super- and subcritical nature of the dynamo transition is given.