Source author record

Pablo D. Mininni

Pablo D. Mininni appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

Multi-time structure functions and the Lagrangian scaling of turbulence

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ármán 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.

preprint2022arXiv

Turbulence generation by large-scale extreme vertical drafts and the modulation of local energy dissipation in stably stratified geophysical flows

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.

preprint2022arXiv

Vector potential-based MHD solver for non-periodic flows using Fourier continuation expansions

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.

preprint2021arXiv

Broken mirror symmetry of tracer's trajectories in turbulence

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.

preprint2021arXiv

Settling and clustering of particles of moderate mass density in turbulence

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.

preprint2020arXiv

Fourier continuation method for incompressible fluids with boundaries

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

preprint2020arXiv

Preferential concentration of free-falling heavy particles in turbulence

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.

preprint2020arXiv

Velocity and acceleration statistics in particle-laden turbulent swirling flows

We present a comparison of different particles' velocity and acceleration statistics in two paradigmatic turbulent swirling flows: the von Kármán 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.

preprint2019arXiv

Spatio-temporal behavior of magnetohydrodynamic fluctuations with cross-helicity and background magnetic field

We study the spatio-temporal behavior of the Elsässer 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énic 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énic effects. Finally, we observe counter-propagation of Alfvénic 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.

preprint2016arXiv

Simulations of the Kelvin-Helmholtz instability driven by coronal mass ejections in the turbulent corona

Recent high resolution AIA/SDO images show evidence of the development of the Kelvin-Helmholtz instability, as coronal mass ejections (CMEs) expand in the ambient corona. A large-scale magnetic field mostly tangential to the interface is inferred, both on the CME and on the background sides. However, the magnetic field component along the shear flow is not strong enough to quench the instability. There is also observational evidence that the ambient corona is in a turbulent regime, and therefore the criteria for the development of the instability are a-priori expected to differ from the laminar case. To study the evolution of the Kelvin-Helmholtz instability with a turbulent background, we perform three-dimensional simulations of the incompressible magnetohydrodynamic equations. The instability is driven by a velocity profile tangential to the CME-corona interface, which we simulate through a hyperbolic tangent profile. The turbulent background is generated by the application of a stationary stirring force. We compute the instability growth-rate for different values of the turbulence intensity, and find that the role of turbulence is to attenuate the growth. The fact that the Kelvin-Helmholtz instability is observed, sets an upper limit to the correlation length of the coronal background turbulence.

preprint2015arXiv

Spatio-temporal detection of Kelvin waves in quantum turbulence simulations

We present evidence of Kelvin excitations in space-time resolved spectra of numerical simulations of quantum turbulence. Kelvin waves are transverse and circularly polarized waves that propagate along quantized vortices, for which the restitutive force is the tension of the vortex line, and which play an important role in theories of superfluid turbulence. We use the Gross-Pitaevskii equation to model quantum flows, letting an initial array of well-organized vortices develop into a turbulent bundle of intertwined vortex filaments. By achieving high spatial and temporal resolution we are able to calculate space-time resolved mass density and kinetic energy spectra. Evidence of Kelvin and sound waves is clear in both spectra. Identification of the waves allows us to extract the spatial spectrum of Kelvin waves, clarifying their role in the transfer of energy

preprint2014arXiv

Kelvin-Helmholtz versus Hall Magneto-shear instability in astrophysical flows

We study the stability of shear flows in a fully ionized plasma. Kelvin-Helmholtz is a well known, macroscopic and ideal shear-driven instability. In sufficiently low density plasmas, also the microscopic Hall magneto-shear instability can take place. We performed three-dimensional simulations of the Hall-MHD equations where these two instabilities are present, and carried out a comparative study. We find that when the shear flow is so intense that its vorticity surpasses the ion-cyclotron frequency of the plasma, the Hall magneto-shear instability is not only non-negligible, but it actually displays growth rates larger than those of the Kelvin-Helmholtz instability.

preprint2013arXiv

On the emergence of helicity in rotating stratified turbulence

We perform numerical simulations of decaying rotating stratified turbulence and show, in the Boussinesq framework, that helicity (velocity-vorticity correlation), as observed in super-cell storms and hurricanes, is spontaneously created due to an interplay between buoyancy and rotation common to large-scale atmospheric and oceanic flows. Helicity emerges from the joint action of eddies and of inertia-gravity waves (with inertia and gravity with respective associated frequencies $f$ and $N$), and it occurs when the waves are sufficiently strong. For $N/f < 3$ the amount of helicity produced is correctly predicted by a quasi-linear balance equation. Outside this regime, and up to the highest Reynolds number obtained in this study, namely $Re\approx 10000$, helicity production is found to be persistent for $N/f$ as large as $\approx 17$, and for $ReFr^2$ and $ReRo^2 $ respectively as large as $\approx 100$ and $\approx 24000$.

preprint2012arXiv

Anisotropy and non-universality in scaling laws of the large scale energy spectrum in rotating turbulence

Rapidly rotating turbulent flow is characterized by the emergence of columnar structures that are representative of quasi-two dimensional behavior of the flow. It is known that when energy is injected into the fluid at an intermediate scale $L_f$, it cascades towards smaller as well as larger scales. In this paper we analyze the flow in the \textit{inverse cascade} range at a small but fixed Rossby number, {$\mathcal{R}o_f \approx 0.05$}. Several {numerical simulations with} helical and non-helical forcing functions are considered in periodic boxes with unit aspect ratio. In order to resolve the inverse cascade range with {reasonably} large Reynolds number, the analysis is based on large eddy simulations which include the effect of helicity on eddy viscosity and eddy noise. Thus, we model the small scales and resolve explicitly the large scales. We show that the large-scale energy spectrum has at least two solutions: one that is consistent with Kolmogorov-Kraichnan-Batchelor-Leith phenomenology for the inverse cascade of energy in two-dimensional (2D) turbulence with a {$\sim k_{\perp}^{-5/3}$} scaling, and the other that corresponds to a steeper {$\sim k_{\perp}^{-3}$} spectrum in which the three-dimensional (3D) modes release a substantial fraction of their energy per unit time to 2D modes. {The spectrum that} emerges {depends on} the anisotropy of the forcing function{,} the former solution prevailing for forcings in which more energy is injected into 2D modes while the latter prevails for isotropic forcing. {In the case of anisotropic forcing, whence the energy} goes from the 2D to the 3D modes at low wavenumbers, large-scale shear is created resulting in another time scale $τ_{sh}$, associated with shear, {thereby producing} a $\sim k^{-1}$ spectrum for the {total energy} with the 2D modes still following a {$\sim k_{\perp}^{-5/3}$} scaling.

preprint2012arXiv

Not Much Helicity is Needed to Drive Large Scale Dynamos

Understanding the in situ amplification of large scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity f_h,C has not been previously well quantified. Using direct numerical simulations for a simple helical dynamo, we show that f_h,C decreases as the ratio of forcing to large scale wave numbers k_F/k_min increases. From the condition that a large scale helical dynamo must overcome the backreaction from any non-helical field on the large scales, we develop a theory that can explain the simulations. For k_F/k_min>8 we find f_h,C< 3%, implying that very small helicity fractions strongly influence magnetic spectra for even moderate scale separation.

preprint2012arXiv

Thermalization and free decay in Surface Quasi-Geostrophic flows

We derive statistical equilibrium solutions of the truncated inviscid surface quasi-geostrophic (SQG) equations, and verify the validity of these solutions at late times in numerical simulations of the truncated SQG equations. The results indicate enstrophy thermalizes while energy can condense at the gravest modes, in agreement with previous indications of a direct cascade of enstrophy and an inverse cascade of energy in forced-dissipative SQG systems. At early times, the truncated inviscid SQG simulations show a behavior reminiscent of forced-dissipative SQG turbulence, and we identify spectral scaling laws for the energy and enstrophy spectra. Finally, a comparison between viscous and inviscid simulations allows us to identify free-decay similarity laws for the enstrophy in SQG turbulence at very large Reynolds number.

preprint2010arXiv

A hybrid MPI-OpenMP scheme for scalable parallel pseudospectral computations for fluid turbulence

A hybrid scheme that utilizes MPI for distributed memory parallelism and OpenMP for shared memory parallelism is presented. The work is motivated by the desire to achieve exceptionally high Reynolds numbers in pseudospectral computations of fluid turbulence on emerging petascale, high core-count, massively parallel processing systems. The hybrid implementation derives from and augments a well-tested scalable MPI-parallelized pseudospectral code. The hybrid paradigm leads to a new picture for the domain decomposition of the pseudospectral grids, which is helpful in understanding, among other things, the 3D transpose of the global data that is necessary for the parallel fast Fourier transforms that are the central component of the numerical discretizations. Details of the hybrid implementation are provided, and performance tests illustrate the utility of the method. It is shown that the hybrid scheme achieves near ideal scalability up to ~20000 compute cores with a maximum mean efficiency of 83%. Data are presented that demonstrate how to choose the optimal number of MPI processes and OpenMP threads in order to optimize code performance on two different platforms.

preprint2010arXiv

Hall-MHD small-scale dynamos

Much of the progress in our understanding of dynamo mechanisms has been made within the theoretical framework of magnetohydrodynamics (MHD). However, for sufficiently diffuse media, the Hall effect eventually becomes non-negligible. We present results from three dimensional simulations of the Hall-MHD equations subjected to random non-helical forcing. We study the role of the Hall effect in the dynamo efficiency for different values of the Hall parameter, using a pseudospectral code to achieve exponentially fast convergence. We also study energy transfer rates among spatial scales to determine the relative importance of the various nonlinear effects in the dynamo process and in the energy cascade. The Hall effect produces a reduction of the direct energy cascade at scales larger than the Hall scale, and therefore leads to smaller energy dissipation rates. Finally, we present results stemming from simulations at large magnetic Prandtl numbers, which is the relevant regime in hot and diffuse media such a the interstellar medium.

preprint2010arXiv

Structures in magnetohydrodynamic turbulence: detection and scaling

We present a systematic analysis of statistical properties of turbulent current and vorticity structures at a given time using cluster analysis. The data stems from numerical simulations of decaying three-dimensional (3D) magnetohydrodynamic turbulence in the absence of an imposed uniform magnetic field; the magnetic Prandtl number is taken equal to unity, and we use a periodic box with grids of up to 1536^3 points, and with Taylor Reynolds numbers up to 1100. The initial conditions are either an X-point configuration embedded in 3D, the so-called Orszag-Tang vortex, or an Arn'old-Beltrami-Childress configuration with a fully helical velocity and magnetic field. In each case two snapshots are analyzed, separated by one turn-over time, starting just after the peak of dissipation. We show that the algorithm is able to select a large number of structures (in excess of 8,000) for each snapshot and that the statistical properties of these clusters are remarkably similar for the two snapshots as well as for the two flows under study in terms of scaling laws for the cluster characteristics, with the structures in the vorticity and in the current behaving in the same way. We also study the effect of Reynolds number on cluster statistics, and we finally analyze the properties of these clusters in terms of their velocity-magnetic field correlation. Self-organized criticality features have been identified in the dissipative range of scales. A different scaling arises in the inertial range, which cannot be identified for the moment with a known self-organized criticality class consistent with MHD. We suggest that this range can be governed by turbulence dynamics as opposed to criticality, and propose an interpretation of intermittency in terms of propagation of local instabilities.

preprint2005arXiv

Shell to shell energy transfer in MHD, Part I: steady state turbulence

We investigate the transfer of energy from large scales to small scales in fully developed forced three-dimensional MHD-turbulence by analyzing the results of direct numerical simulations in the absence of an externally imposed uniform magnetic field. Our results show that the transfer of kinetic energy from the large scales to kinetic energy at smaller scales, and the transfer of magnetic energy from the large scales to magnetic energy at smaller scales, are local, as is also found in the case of neutral fluids, and in a way that is compatible with Kolmogorov (1941) theory of turbulence. However, the transfer of energy from the velocity field to the magnetic field is a highly non-local process in Fourier space. Energy from the velocity field at large scales can be transfered directly into small scale magnetic fields without the participation of intermediate scales. Some implications of our results to MHD turbulence modeling are also discussed.

preprint2005arXiv

Shell to shell energy transfer in MHD, Part II: Kinematic dynamo

We study the transfer of energy between different scales for forced three-dimensional MHD turbulent flows in the kinematic dynamo regime. Two different forces are examined: a non-helical Taylor Green flow with magnetic Prandtl number P_M=0.4, and a helical ABC flow with P_M=1. This analysis allows us to examine which scales of the velocity flow are responsible for dynamo action, and identify which scales of the magnetic field receive energy directly from the velocity field and which scales receive magnetic energy through the cascade of the magnetic field from large to small scales. Our results show that the turbulent velocity fluctuations are responsible for the magnetic field amplification in the small scales (small scale dynamo) while the large scale field is amplified mostly due to the large scale flow. A direct cascade of the magnetic field energy from large to small scales is also present and is a complementary mechanism for the increase of the magnetic field in the small scales. Input of energy from the velocity field in the small magnetic scales dominates over the energy that is cascaded down from the large scales until the large-scale peak of the magnetic energy spectrum is reached. At even smaller scales, most of the magnetic energy input is from the cascading process.