Search

We investigate a shift in the critical temperature of rotating Bose-Einstein condensates mediated by the melting of the vortex lattice. Numerical simulations reveal that this temperature exhibits contrasting behavior depending on the system configuration: a negative shift occurs for fixed trap potentials due to the expansion of the condensate, while a positive shift is observed for fixed volumes, where vortex lattice rigidity suppresses thermal fluctuations. We introduce a vortex-energy model that captures the role of vortex interactions, the positional energy of the vortex lattice, as well as the phase transition and how the vortex lattice disappears. The findings provide insights into the thermodynamic properties of rotating condensates and the dynamics of vortex lattice melting, offering potential parallels with other quantum systems such as type-II superconductors.

How turbulent convective fluctuations organise to form large-scale structures in planetary atmospheres remains a question that eludes quantitative answers. The assumption that this process is the result of an inverse cascade was suggested half a century ago in two-dimensional fluids, but its applicability to atmospheric and oceanic flows remains heavily debated, hampering our understanding of the energy balance in planetary systems. We show with direct numerical simulations of spatial resolutions of 122882 $\times$ 384 points that rotating and stratified flows can support a bidirectional cascade of energy, in three dimensions, with a ratio of Rossby to Froude numbers comparable to that of the Earth's atmosphere. Our results establish that in dry atmospheres spontaneous order can arise via an inverse cascade to the largest spatial scales.

Active grids operated with random protocols are a standard way to generate large Reynolds number turbulence in wind and water tunnels. But anomalies in the decay and third-order scaling of active-grid turbulence have been reported. We combine Laser Doppler Velocimetry and hot-wire anemometry measurements in a wind tunnel, with machine learning techniques and numerical simulations, to gain further understanding on the reasons behind these anomalies. Numerical simulations that incorporate the statistical anomalies observed in the experimental velocity field near the active grid can reproduce the experimental anomalies observed later in the decay. The results indicate that anomalies in experiments near the active grid introduce correlations in the flow that persist for long times, and result in the flow being statistically different from homogeneous and isotropic turbulence.

We study density isolines in quantum turbulence under the Schramm-Loewner framework using direct numerical simulations of the truncated Gross-Pitaevskii equation, in both spherical and cylindrical traps with three-dimensional dynamics. Density isolines develop increasing complexity as turbulence matures. As the systems evolves towards a thermalized regime, it spontaneously develops conformal invariance. In contrast to other systems exhibiting conformal invariance, this system manifests it during the transition towards disorder rather than to self-organization. We discuss a link between this behavior in quantum turbulence and other 4-wave interacting systems.

The Basset-Boussinesq force is often perfunctorily neglected when studying small inertial particles in turbulence. This force arises from the diffusion of vorticity from the particles and, since it depends on the particles' past history, complicates the dynamics by transforming their equations of motion into integro-differential equations. However, this force is of the same order as other viscous forces acting on the particles, and beyond convenience, the reasons for neglecting it are unclear. This study addresses the following question: Under what conditions can the Basset-Boussinesq force be neglected in light particles in geophysical flows? We derive strict bounds for the magnitude of the Basset-Boussinesq force in stably stratified flows, in contexts of interest for geophysical turbulence. The bounds are validated by direct numerical simulations. The Basset-Boussinesq force is negligible when a buoyancy Stokes number $\textrm{Sb} = N \tau_p$ is small, where $N$ is the flow Brunt-V\"ais\"al\"a frequency and $\tau_p$ is the particle's Stokes time. Interestingly, for most oceanic particles this force may be negligible. Only for very strong stratification, or for particles with very large inertia, this force must be considered in the dynamics.

We show that the unsteadiness of turbulence has a drastic effect on turbulence parameters and in particle cluster formation. To this end we use direct numerical simulations of particle laden flows with a steady forcing that generates an unsteady large-scale flow. Particle clustering correlates with the instantaneous Taylor-based flow Reynolds number, and anti-correlates with its instantaneous turbulent energy dissipation constant. A dimensional argument for these correlations is presented. In natural flows, unsteadiness can result in extreme particle clustering, which is stronger than the clustering expected from averaged inertial turbulence effects.

We consider a quantum Otto cycle with an interacting Bose-Einstein condensate at finite temperature. We present a procedure to evolve this system in time in three spatial dimensions, in which closed (adiabatic) strokes are described by the Gross-Pitaevskii equation, and open (isochoric) strokes are modeled using a stochastic Ginzburg-Landau equation. We analyze the effect on the thermodynamic efficiency of the strength of interactions, the frequency of the harmonic trap, and the temperatures of the reservoirs. The efficiency has little sensitivity to changes in the temperatures, but decreases as interactions increase. However, stronger interactions allow for faster cycles and for substantial increases in power.

Inertial particles in stably stratified flows play a fundamental role in geophysics, from the dynamics of nutrients in the ocean to the dispersion of pollutants in the atmosphere. We consider the Maxey-Riley equation for small neutrally buoyant inertial particles in the Boussinesq approximation, and discuss its limits of validity. We show that particles behave as forced damped oscillators, with different regimes depending on the particles Stokes number and the fluid Brunt-V\"ais\"al\"a frequency. Using direct numerical simulations we study the particles dynamics and we show that small neutrally buoyant particles in these flows tend to cluster in regions of low local vorticity. The particles, albeit small, behave fundamentally differently than tracers.

We study the long-time decay of rotating turbulence in Bose-Einstein condensates (BECs). We consider the Gross-Pitaevskii equation in a rotating frame of reference, and review different formulations for the Hamiltonian of a rotating BEC. We discuss how the energy can be decomposed, and present a method to generate out-of-equilibrium initial conditions. We also present a method to generate finite-temperature states of rotating BECs compatible with the Canonical or the Grand canonical ensembles. Finally, we integrate numerically rotating BECs in cigar-shaped traps. A transition is found in the system dynamics as the rotation rate is increased, with a final state of the decay of the turbulent flow compatible with an Abrikosov lattice in a finite-temperature thermalized state., Comment: 10 figures, 12 pages

We present a high-order numerical study of the dependence of the dynamo onset on different electromagnetic boundary conditions, in convecting Boussinesq flows forced by a temperature gradient. Perfectly conducting boundaries, vacuum, and mixed electromagnetic boundary conditions are considered, using a method that treats fields and boundary conditions with close to spectral accuracy. Having one or two conducting boundaries greatly facilitates dynamo action. For the mixed case it is shown that the critical magnetic Reynolds number becomes independent of the Rayleigh number, Ra, for sufficiently large Ra.

We study the three-dimensional clustering of velocity stagnation points, of nulls of the vorticity and of the Lagrangian acceleration, and of inertial particles in turbulent flows at fixed Reynolds numbers, but under different large-scale flow geometries. To this end, we combine direct numerical simulations of homogeneous and isotropic turbulence and of the Taylor-Green flow, with particle tracking velocimetry in a von K\'arm\'an experiment. While flows have different topologies (as nulls cluster differently), particles behave similarly in all cases, indicating that Taylor-scale neutrally buoyant particles cluster as inertial particles.

Since the idea of quantum turbulence was first proposed by Feynman, and later realized in experiments of superfluid helium and Bose-Einstein condensates, much emphasis has been put in finding signatures that distinguish quantum turbulence from its classical counterpart. Here we show that quantum turbulence in rotating condensates is fundamentally different from the classical case. While rotating quantum turbulence develops a negative temperature state with self-organization of the kinetic energy in quantized vortices, it also displays an anisotropic dissipation mechanism and a different, non-Kolmogorovian, scaling of the energy at small scales. This scaling is compatible with Vinen turbulence and is also found in recent simulations of condensates with multicharged vortices. An elementary explanation for the scaling is presented in terms of disorder in the vortices positions., Comment: 8 pages, 9 figures

We define and characterize multi-time Lagrangian structure functions using data stemming from two swirling flows with mean flow and turbulent fluctuations: A Taylor-Green numerical flow, and a von K\'arm\'an laboratory experiment. Data is obtained from numerical integration of tracers in the former case, and from three-dimensional particle tracking velocimetry measurements in the latter. Multi-time statistics are shown to decrease the contamination of large scales in the inertial range scaling. A time scale at which contamination from the mean flow becomes dominant is identified, with this scale separating two different Lagrangian scaling ranges. The results from the multi-time structure functions also indicate that Lagrangian intermittency is not a result of large-scale flow effects. The multi-time Lagrangian structure functions can be used without prior knowledge of the forcing mechanisms or boundary conditions, allowing their application in different flow geometries.

Topological properties of physical systems play a crucial role in our understanding of nature, yet their experimental determination remains elusive. We show that the mean helicity, a dynamical invariant in ideal flows, quantitatively affects trajectories of fluid elements: the linking number of Lagrangian trajectories depends on the mean helicity. Thus, a global topological invariant and a topological number of fluid trajectories become related, and we provide an empirical expression linking them. The relation shows the existence of long-term memory in the trajectories: the links can be made of the trajectory up to a given time, with particles positions in the past. This property also allows experimental measurements of mean helicity.

A high-order method to evolve in time electromagnetic and velocity fields in conducting fluids with non-periodic boundaries is presented. The method has a small overhead compared with fast FFT-based pseudospectral methods in periodic domains. It uses the magnetic vector potential formulation for accurately enforcing the null divergence of the magnetic field, and allowing for different boundary conditions including perfectly conducting walls or vacuum surroundings, two cases relevant for many astrophysical, geophysical, and industrial flows. A spectral Fourier continuation method is used to accurately represent all fields and their spatial derivatives, allowing also for efficient solution of Poisson equations with different boundaries. A study of conducting flows at different Reynolds and Hartmann numbers, and with different boundary conditions, is presented to study convergence of the method and the accuracy of the solenoidal and boundary conditions.

We observe the emergence of strong vertical drafts in direct numerical simulations of the Boussinesq equations in a range of parameters of geophysical interest. These structures, which appear intermittently in space and time, generate turbulence and enhance kinetic and potential energy dissipation, providing a possible explanation for the observed variability of the local energy dissipation in the bulk of oceanic flows, and the modulation of its probability distribution function. We show how, due to the extreme drafts, in runs with Froude numbers observable in geophysical scenarios, roughly 10% of the domain flow can account for up to 50% of the global volume dissipation, reminiscent of estimates based on oceanic models., Comment: 16 pages, 6 figures

The reduction of dimensionality of physical systems, specially in fluid dynamics, leads in many situations to nonlinear ordinary differential equations which have global invariant manifolds with algebraic expressions containing relevant physical information of the original system. We present a method to identify such manifolds, and we apply it to a reduced model for the Lagrangian evolution of field gradients in homogeneous and isotropic turbulence with a passive scalar., Comment: 6 pages, 3 figures

Based on direct numerical simulations with point-like inertial particles, with Stokes numbers, $\textrm{St}=0, 0.5$, $3$, and $6$, transported by homogeneous and isotropic turbulent flows, we present in this letter for the first time evidence for the existence of Markov property in Lagrangian turbulence. We show that the Markov property is valid for a finite step size larger than a Stokes number-dependent Einstein-Markov coherence time scale. This enables the description of multi-scale statistics of Lagrangian particles by Fokker-Planck equations, which can be embedded in an interdisciplinary approach linking the statistical description of turbulence with fluctuation theorems of non-equilibrium stochastic thermodynamics and local flow structures. The formalism allows estimation of the stochastic thermodynamics entropy exchange associated with the particles' Lagrangian trajectories. Entropy consuming trajectories of the particles are related to specific evolution of velocity increments through scales and may be seen as intermittent structures. Statistical features of Lagrangian paths and entropy values are thus fixed by the fluctuation theorems., Comment: Accepted in EPL, 6 pages, 4 figures

We present a numerical study of settling and clustering of small inertial particles in homogeneous and isotropic turbulence. Particles are denser than the fluid, but not in the limit of being much heavier than the displaced fluid. At fixed Reynolds and Stokes numbers we vary the fluid-to-particle mass ratio and the gravitational acceleration. The effect of varying one or the other is similar but not quite the same. We report non-monotonic behavior of the particles' velocity skewness and kurtosis with the second parameter, and an associated anomalous behavior of the settling velocity when compared to the free-fall Stokes velocity, including loitering cases. Clustering increases for increasing gravitational acceleration, and for decreasing fluid-to-particle mass ratio.

We analyze the vector nulls of velocity, Lagrangian acceleration, and vorticity, coming from direct numerical simulations of forced homogeneous isotropic turbulence at $Re_\lambda \in [40-610]$. We show that the clustering of velocity nulls is much stronger than those of acceleration and vorticity nulls. These acceleration and vorticity nulls, however, are denser than the velocity nulls. We study the scaling of clusters of these null points with $Re_\lambda$ and with characteristic turbulence lengthscales. We also analyze datasets of point inertial particles with Stokes numbers $St = 0.5$, 3, and 6, at $Re_\lambda = 240$. Inertial particles display preferential concentration with a degree of clustering that resembles some properties of the clustering of the Lagrangian acceleration nulls, in agreement with the proposed sweep-stick mechanism of clustering formation.

We present results of different approaches to model the evolution of the COVID-19 epidemic in Argentina, with a special focus on the megacity conformed by the city of Buenos Aires and its metropolitan area, including a total of 41 districts with over 13 million inhabitants. We first highlight the relevance of interpreting the early stage of the epidemic in light of incoming infectious travelers from abroad. Next, we critically evaluate certain proposed solutions to contain the epidemic based on instantaneous modifications of the reproductive number. Finally, we build increasingly complex and realistic models, ranging from simple homogeneous models used to estimate local reproduction numbers, to fully coupled inhomogeneous (deterministic or stochastic) models incorporating mobility estimates from cell phone location data. The models are capable of producing forecasts highly consistent with the official number of cases with minimal parameter fitting and fine-tuning. We discuss the strengths and limitations of the proposed models, focusing on the validity of different necessary first approximations, and caution future modeling efforts to exercise great care in the interpretation of long-term forecasts, and in the adoption of non-pharmaceutical interventions backed by numerical simulations., Comment: 16 pages, 14 figures

We present a sweep-stick mechanism for heavy particles transported by a turbulent flow under the action of gravity. Direct numerical simulations show that these particles preferentially explore regions of the flow with close to zero Lagrangian acceleration. However, the actual Lagrangian acceleration of the fluid elements where particles accumulate is not zero, and has a dependence on the Stokes number, the gravity acceleration, and the settling velocity of the particles., Comment: 8 pages, 9 figures

We present numerical evidence of a critical-like transition in an out-of-equilibrium mean-field description of a quantum system. By numerically solving the Gross-Pitaevskii equation we show that quantum turbulence displays an abrupt change between three-dimensional (3D) and two-dimensional (2D) behavior. The transition is observed both in quasi-2D flows in cubic domains (controlled by the amplitude of a 3D perturbation to the flow), as well as in flows in thin domains (controlled by the domain aspect ratio) in a configuration that mimics systems realized in laboratory experiments. In one regime the system displays a transfer of the energy towards smaller scales, while in the other the system displays a transfer of the energy towards larger scales and a coherent self-organization of the quantized vortices.

We present a Fourier Continuation-based parallel pseudospectral method for incompressible fluids in cuboid non-periodic domains. The method produces dispersionless and dissipationless derivatives with fast spectral convergence inside the domain, and with very high order convergence at the boundaries. Incompressibility is imposed by solving a Poisson equation for the pressure. Being Fourier-based, the method allows for fast computation of spectral transforms. It is compatible with uniform grids (although refined or nested meshes can also be implemented), which in turn allows for explicit time integration at sufficiently high Reynolds numbers. Using a new parallel code named SPECTER we illustrate the method with two problems: channel flow, and plane Rayleigh-B\'enard convection under the Boussinesq approximation. In both cases the method yields results compatible with previous studies using other high-order numerical methods, with mild requirements on the time step for stability.

We present a comparison of different particles' velocity and acceleration statistics in two paradigmatic turbulent swirling flows: the von K\'arm\'an flow in a laboratory experiment, and the Taylor-Green flow in direct numerical simulations. Tracers, as well as inertial particles, are considered. Results indicate that, in spite of the differences in boundary conditions and forcing mechanisms, scaling properties and statistical quantities reveal similarities between both flows, pointing to new methods to calibrate and compare models for particles dynamics in numerical simulations, as well as to characterize the dynamics of particles in simulations and experiments.

We show that the phase space of stratified turbulence mainly consists of two slow invariant manifolds with rich physics, embedded on a larger basin with fast evolution. A local invariant manifold in the vicinity of the fluid at equilibrium corresponds to waves, while a global invariant manifold corresponds to the onset of local convection. Using a reduced model derived from the Boussinesq equations, we propose that waves accumulate energy nonlinearly up to a point such that fluid elements escape from the local manifold and evolve fast to the global manifold, where kinetic energy can be more efficiently dissipated. After this, fluid elements return to the first manifold. As the stratification increases, the volume of the first manifold increases, and the second manifold becomes harder to access. This explains recent observations of enhanced intermittency and marginal instability in these flows. The reduced model also allows us to study structure formation, alignment of field gradients in the flow, and to identify balance relations that hold for each fluid element., Comment: 21 pages, 14 figures

We study the spatio-temporal behavior of the Els\"asser variables describing magnetic and velocity field fluctuations, using direct numerical simulations of three-dimensional magnetohydrodynamic turbulence. We consider cases with relatively small, intermediate, and large values of a mean background magnetic field, and with null, small, and high cross-helicity (correlations between the velocity and the magnetic field). Wavenumber-dependent time correlation functions are computed for the different simulations. From these correlation functions, the decorrelation time is computed and compared with different theoretical characteristic times: the local non-linear time, the random-sweeping time, and the Alfv\'enic time. It is found that decorrelation times are dominated by sweeping effects for low values of the mean magnetic field and for low values of the cross-helicity, while for large values of the background field or of the cross-helicity and for wave vectors sufficiently aligned with the guide field, decorrelation times are controlled by Alfv\'enic effects. Finally, we observe counter-propagation of Alfv\'enic fluctuations due to reflections produced by inhomogeneities in the total magnetic field. This effect becomes more prominent in flows with large cross-helicity, strongly modifying the propagation of waves in turbulent magnetohydrodynamic flows., Comment: 16 pages, 40 figures

We present a reduced system of 7 ordinary differential equations that captures the time evolution of spatial gradients of the velocity and the temperature in fluid elements of stratified turbulent flows. We show the existence of invariant manifolds (further reducing the system dimensionality), and compare the results with data stemming from direct numerical simulations of the full incompressible Boussinesq equations in the stably stratified case. Numerical results accumulate over the invariant anifolds of the reduced system, indicating the system lives at the brink of an instability. Finally, we study the stability of the reduced system, and show that it is compatible with recent observations in stratified turbulence of non-monotonic dependence of intermittency with stratification., Comment: 7 pages, 4 figures

We present laboratory experiments of surface wave turbulence excited by paddles in the deep water regime. The free surface is seeded with buoyant particles that are advected and dispersed by the flow. Positions and velocities of the floaters are measured using particle tracking velocimetry. We study the statistics of velocity and acceleration of the particles, mean vertical displacements, single-particle horizontal dispersion, and the phenomenon of preferential concentration. Using a simple model together with the experimental data, we show that the time evolution of the particles has three characteristic processes that dominate the dynamics at different times: drag by surface waves at early times, trapping by horizontal eddies at intermediate times, and advection by a large-scale mean circulation at late times., Comment: 16 pages, 12 figures

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

We study the effect of different forcing functions and of the local gradient Richardson number Rig on the vertical mixing of Lagrangian tracers in stably stratified turbulence under the Boussinesq approximation, and present a wave and continuous-time random walk model for single- and two- particle vertical dispersion. The model consists of a random superposition of linear waves with their amplitude based on the observed Lagrangian spectrum of vertical velocity, and a random walk process to capture overturning that depends on the statistics of Rig among other Eulerian quantities. The model is in good agreement with direct numerical simulations of stratified turbulence, where single- and two-particle dispersion differs from the homogeneous and isotropic case. Moreover, the model gives insight into the mixture of linear and non-linear physics in the problem, as well as on the different processes responsible for vertical turbulent dispersion.

A parametric study of the magnetic dipole behavior in resistive incompressible MHD inside a rotating sphere is performed, using direct numerical simulations and considering Reynolds and Ekman numbers as controlling parameters. The tendency is to obtain geodynamo-like magnetic dipole reversal regimes for sufficiently small Ekman and large Reynolds numbers. The typical dipole latitude obtained in the reversal regime is around 40 degrees (with respect to the rotation axis of the sphere). A statistical analysis of waiting times between dipole reversals is also performed, obtaining a non-Poissonian distribution of waiting times, indicating long-term memory effects. We also report the presence of a $1/f$ frequency power spectrum in the magnetic dipole time-series, which also shows a tendency to grow toward lower frequencies as the Ekman number is decreased., Comment: 14 pages, 9 figures

An existing hybrid MPI-OpenMP scheme is augmented with a CUDA-based fine grain parallelization approach for multidimensional distributed Fourier transforms, in a well-characterized pseudospectral fluid turbulence code. Basics of the hybrid scheme are reviewed, and heuristics provided to show a potential benefit of the CUDA implementation. The method draws heavily on the CUDA runtime library to handle memory management, and on the cuFFT library for computing local FFTs. The manner in which the interfaces are constructed to these libraries, and ISO bindings utilized to facilitate platform portability, are discussed. CUDA streams are implemented to overlap data transfer with cuFFT computation. Testing with a baseline solver demonstrates significant aggregate speed-up over the hybrid MPI-OpenMP solver by offloading to GPUs on an NVLink-based test system. While the batch streamed approach provides little benefit with NVLink, we see a performance gain of 30% when tuned for the optimal number of streams on a PCIe-based system. It is found that strong GPU scaling is ideal, or slightly better than ideal, in all cases. In addition to speed-up measurements for the fiducial solver, we also consider several other solvers with different numbers of transform operations and find that aggregate speed-ups are nearly constant for all solvers., Comment: 17 pages, 8 figures; submitted to Parallel Computing

We investigate the large-scale intermittency of vertical velocity and temperature, and the mixing properties of stably stratified turbulent flows using both Lagrangian and Eulerian fields from direct numerical simulations, in a parameter space relevant for the atmosphere and the oceans. Over a range of Froude numbers of geophysical interest ($\approx 0.05-0.3$) we observe very large fluctuations of the vertical components of the velocity and the potential temperature, localized in space and time, with a sharp transition leading to non-Gaussian wings of the probability distribution functions. This behavior is captured by a simple model representing the competition between gravity waves on a fast time-scale and nonlinear steepening on a slower time-scale. The existence of a resonant regime characterized by enhanced large-scale intermittency, as understood within the framework of the proposed model, is then linked to the emergence of structures in the velocity and potential temperature fields, localized overturning and mixing. Finally, in the same regime we observe a linear scaling of the mixing efficiency with the Froude number and an increase of its value of roughly one order of magnitude.

We perform a study on the evolution of helical quantum turbulence at different temperatures by solving numerically the Gross-Pitaevskii and the Stochastic Ginzburg-Landau equations, using up to $4096^3$ grid points with a pseudospectral method. We show that for temperatures close to the critical the fluid described by these equations can act as a classical viscous flow, with the decay of the incompressible kinetic energy and the helicity becoming exponential. The transition from this behavior to the one observed at zero temperature is smooth as a function of temperature. Moreover, the presence of strong thermal effects can inhibit the development of a proper turbulent cascade. We provide anzats for the effective viscosity and friction as a function of the temperature

The spectrally truncated, or finite dimensional, versions of several equations of inviscid flows display transient solutions which match their viscous counterparts, but which eventually lead to thermalized states in which energy is in equipartition between all modes. Recent advances in the study of the Burgers equation show that the thermalization process is triggered after the formation of sharp localized structures within the flow called "tygers". We show that the process of thermalization first takes place in well defined subdomains, before engulfing the whole space. Using spatio-temporal analysis on data from numerical simulations, we study propagation of tygers and find that they move at a well defined mean speed that can be obtained from energy conservation arguments.

In this paper, we present an overview of concepts and data concerning inverse cascades of excitation towards scales larger than the forcing scale in a variety of contexts, from two-dimensional fluids and wave turbulence, to geophysical flows in the presence of rotation and stratification. We briefly discuss the role of anisotropy in the occurrence and properties of such cascades. We then show that the cascade of some invariant, for example the total energy, may be transferred through nonlinear interactions to both the small scales and the large scales, with in each case a constant flux. This is in contrast to the classical picture, and we illustrate such a dual cascade in the context of atmospheric and oceanic observations, direct numerical simulations and modeling. We also show that this dual cascade of total energy can in fact be decomposed in some cases into separate cascades of the kinetic and potential energies, provided the Froude and Rossby numbers are small enough. In all cases, the potential energy flux remains small, of the order of 10% or less relative to the kinetic energy flux. Finally, we demonstrate that, in the small-scale inertial range, approximate equipartition between potential and kinetic modes is obtained, leading to an energy ratio close to one, with strong departure at large scales due to the dominant kinetic energy inverse cascade and piling-up at the lowest spatial frequency, and at small scales due to unbalanced dissipation processes, even though the Prandtl number is equal to one.

Using spatio-temporal spectra we show direct evidence of excitation of magnetosonic and Alfv\'en waves in three-dimensional compressible magnetohydrodynamic turbulence at small Mach numbers. For the plasma pressure dominated regime, or high $\beta$ regime (with $\beta$ the ratio between fluid and magnetic pressure), and for the magnetic pressure dominated regime, or low $\beta$ regime, we study magnetic field fluctuations parallel and perpendicular to a guide magnetic field $\textbf{B}_0$. In the low $\beta$ case we find excitation of compressible and incompressible fluctuations, with a transfer of energy towards Alfv\'enic modes and to a lesser extent towards magnetosonic modes. In particular, we find signatures of the presence of fast magnetosonic waves in a scenario compatible with that of weak turbulence. In the high $\beta$ case, fast and slow magnetosonic waves are present, with no clear trace of Alfv\'en waves, and a significant part of the energy is carried by two-dimensional turbulent eddies.

The large-scale structures in the ocean and the atmosphere are in geostrophic balance, and a conduit must be found to channel the energy to the small scales where it can be dissipated. In turbulence this takes the form of an energy cascade, whereas one possible mechanism in a balanced flow at large scales is through the formation of fronts, a common occurrence in geophysical dynamics. We show in this paper that an iconic configuration in laboratory and numerical experiments for the study of turbulence, that of the Taylor-Green or von K\'arm\'an swirling flow, can be suitably adapted to the case of fluids with large aspect ratios, leading to the creation of an imposed large-scale vertical shear. To this effect we use direct numerical simulations of the Boussinesq equations without net rotation and with no small-scale modeling, and with this idealized Taylor-Green set-up. Various grid spacings are used, up to $2048^2\times 256$ spatial points. The grids are always isotropic, with box aspect ratios of either $1:4$ or $1:8$. We find that when shear and stratification are comparable, the imposed shear layer resulting from the forcing leads to the formation of multiple fronts and filaments which destabilize and further evolve into a turbulent flow in the bulk, with a sizable amount of dissipation and mixing, and with a cycle of front creation, instability, and development of turbulence. The results depend on the vertical length scales for shear and for stratification, with stronger large-scale gradients being generated when the two length scales are comparable., Comment: 19 pages, 15 figures, several simulations added in this new version

We present models for single-particle dispersion in vertical and horizontal directions of stably stratified flows. The model in the vertical direction is based on the observed Lagrangian spectrum of the vertical velocity, while the model in the horizontal direction is a combination of a continuous-time eddy-constrained random walk process with a contribution to transport from horizontal winds. Transport at times larger than the Lagrangian turnover time is not universal and dependent on these winds. The models yield results in good agreement with direct numerical simulations of stratified turbulence, for which single-particle dispersion differs from the well studied case of homogeneous and isotropic turbulence.

While in classical turbulence helicity depletes nonlinearity and can alter the evolution of turbulent flows, in quantum turbulence its role is not fully understood. We present numerical simulations of the free decay of a helical quantum turbulent flow using the Gross-Pitaevskii equation at high spatial resolution. The evolution has remarkable similarities with classical flows, which go as far as displaying a dual transfer of incompressible kinetic energy and helicity to small scales. Spatio-temporal analysis indicates that both quantities are dissipated at small scales through non-linear excitation of Kelvin waves and the subsequent emission of phonons. At the onset of the decay, the resulting turbulent flow displays polarized large scale structures and unpolarized patches of quiescense reminiscent of those observed in simulations of classical turbulence at very large Reynolds numbers., Comment: Fixed typos

Kraichnan seminal ideas on inverse cascades yielded new tools to study common phenomena in geophysical turbulent flows. In the atmosphere and the oceans, rotation and stratification result in a flow that can be approximated as two-dimensional at very large scales, but which requires considering three-dimensional effects to fully describe turbulent transport processes and non-linear phenomena. Motions can thus be classified into two classes: fast modes consisting of inertia-gravity waves, and slow quasi-geostrophic modes for which the Coriolis force and horizontal pressure gradients are close to balance. In this paper we review previous results on the strength of the inverse cascade in rotating and stratified flows, and then present new results on the effect of varying the strength of rotation and stratification (measured by the ratio $N/f$ of the Brunt-V\"ais\"ala frequency to the Coriolis frequency) on the amplitude of the waves and on the flow quasi-geostrophic behavior. We show that the inverse cascade is more efficient in the range of $N/f$ for which resonant triads do not exist, $1/2 \le N/f \le 2$. We then use the spatio-temporal spectrum, and characterization of the flow temporal and spatial scales, to show that in this range slow modes dominate the dynamics, while the strength of the waves (and their relevance in the flow dynamics) is weaker., Comment: 17 pages, 13 figures, new version has all changes in the paper published in the Two Dimensional Turbulence Focus Issue of Physics of Fluids

In this paper we analyze the effect of dynamical three-dimensional MHD turbulence on test particle acceleration, and compare how this evolving system affects particle energization by current sheets interaction, against frozen-in-time fields. To do this we analize the ensamble particle acceleration for static electromagnetic fields extracted from direct numerical simulations of the MHD equations, and compare with the dynamical fields. We show that a reduction in particle acceleration in the dynamical model results from the particle trapping in the field lines, which forces the particles to remain in a moving current sheet that suppress the longer exposure at the strong electric field gradients located between structures, which is an efficient particle acceleration mechanism. In addition, we analize the effect of anisotropy caused by the mean magnetic field. It is well known that for sufficiently strong external fields, the system suffers a transition towards a two-dimensional flow. This causes an increment in the size of the coherent structures, resulting in a magnetized state of the particles and the reduction of the particle energization., Comment: 13 pages,7 figures

Using direct numerical simulations of three-dimensional magnetohydrodynamic (MHD) turbulence the spatio-temporal behavior of magnetic field fluctuations is analyzed. Cases with relatively small, medium and large values of a mean background magnetic field are considered. The (wavenumber) scale dependent time correlation function is directly computed for different simulations, varying the mean magnetic field value. From this correlation function the time decorrelation is computed and compared with different theoretical times, namely, the local non-linear time, the random sweeping time, and the Alfv\'enic time, the latter being a wave effect. It is observed that time decorrelations are dominated by sweeping effects, and only at large values of the mean magnetic field and for wave vectors mainly aligned with this field time decorrelations are controlled by Alfv\'enic effects., Comment: 12 pages, 13 figures

We study the transition in dimensionality of a three-dimensional magnetohydrodynamic flow forced only mechanically, when the strength of a magnetic guiding field is gradually increased. We use numerical simulations to consider cases in which the mechanical forcing injects (or not) helicity in the flow. As the guiding field is increased, the strength of the magnetic field fluctuations decrease as a power law of the guiding field intensity. We show that for strong enough guiding fields, the helical magnetohydrodynamic flow can become almost two-dimensional. In this case, the mechanical energy can undergo a process compatible with an inverse cascade, being transferred preferentially towards scales larger than the forcing scale. The presence of helicity changes the spectral scaling of the small magnetic field fluctuations, and affects the statistics of the velocity field and of the velocity gradients. Moreover, at small scales the dynamics of the flow becomes dominated by a direct cascade of helicity, which can be used to derive scaling laws for the velocity field., Comment: 11 pages, 11 figures

Nonlinear triadic interactions are at the heart of our understanding of turbulence. In flows where waves are present modes must not only be in a triad to interact, but their frequencies must also satisfy an extra condition: the interactions that dominate the energy transfer are expected to be resonant. We derive equations that allow direct measurement of the actual degree of resonance of each triad in a turbulent flow. We then apply the method to the case of rotating turbulence, where eddies coexist with inertial waves. We show that for a range of wave numbers, resonant and near-resonant triads are dominant, the latter allowing a transfer of net energy towards two-dimensional modes that would be inaccessible otherwise. The results are in good agreement with approximations often done in theories of rotating turbulence, and with the mechanism of parametric instability proposed to explain the development of anisotropy in such flows. We also observe that, at least for the moderate Rossby numbers studied here, marginally near-resonant and non-resonant triads play a non-negligible role in the coupling of modes.

The effect of compressibility in charged particle energization by magnetohydrodynamic (MHD) fields is studiedin the context of test particle simulations. This problem is relevant to the solar wind and the solar corona due to the compressible nature of the flow in those astrophysical scenarios. We consider turbulent electromagnetic fields obtained from direct numerical simulations of the MHD equations with a strong background magnetic field. In order to explore the flow compressibilty effect over the particle dynamics we performed different numerical experiments: an incompressible case, and two weak compressible cases with Mach number M = 0.1 and M = 0.25. We analyze the behavior of protons and electrons in those turbulent fields, which are well known to form aligned current sheets in the direction of the guide magnetic field. What we call protons and electrons are test particles with scales comparable to (for protons) and much smaller than (for electrons) the dissipative scale of MHD turbulence, maintaining the correct mass ratio me /mi. For these test particles we show that compressibility enhances the efficiency of proton acceleration, and that the energization is caused by perpendicular electric fields generated between currents sheets. On the other hand, electrons remain magnetized and display an almost adiabatic motion, with no effect of compressibility observed. Another set of numerical experiments takes into account two fluid modifications, namely electric field due to Hall effect and electron pressure gradient. We show that the electron pressure has an important contribution to electron acceleration allowing highly parallel energization. In contrast, no significant effect of these additional terms is observed for the protons., Comment: 10 pages

Helicity is a topological invariant that measures the linkage and knottedness of lines, tubes and ribbons. As such, it has found myriads of applications in astrophysics and solar physics, in fluid dynamics, in atmospheric sciences, and in biology. In quantum flows, where topology-changing reconnection events are a staple, helicity appears as a key quantity to study. However, the usual definition of helicity is not well posed in quantum vortices, and its computation based on counting links and crossings of vortex lines can be downright impossible to apply in complex and turbulent scenarios. We present a new definition of helicity which overcomes these problems. With it, we show that only certain reconnection events conserve helicity. In other cases helicity can change abruptly during reconnection. Furthermore, we show that these events can also excite Kelvin waves, which slowly deplete helicity as they interact nonlinearly, thus linking the theory of vortex knots with observations of quantum turbulence.