Source author record

Annick Pouquet

Annick Pouquet 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
12topics
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

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.

preprint2021arXiv

Helical Fluid and (Hall)-MHD Turbulence: a Brief Review

Helicity, a measure of the breakage of reflectional symmetry representing the topology of turbulent flows, contributes in a crucial way to their dynamics and to their fundamental statistical properties. We review several of their main features, both new and old, as the discovery of bi-directional cascades or the role of helical vortices in the enhancement of large-scale magnetic fields in the dynamo problem. The dynamical contribution in magnetohydrodynamic (MHD) of the cross-correlation between velocity and induction is discussed as well. We consider next how turbulent transport is affected by helical constraints, in particular in the context of magnetic reconnection and fusion plasmas {under one- and two-fluid approximations}. Central issues on how to construct turbulence models for non-reflectionally symmetric helical flows are reviewed, including in the presence of shear, and we finally briefly mention the possible role of helicity in the development of strongly localized quasi-singular structures {at small scale).

preprint2020arXiv

Coupling large eddies and waves in turbulence: Case study of magnetic helicity at the ion inertial scale

In turbulence, for neutral or conducting fluids, a large ratio of scales is excited because of the possible occurrence of inverse cascades to large, global scales together with direct cascades to small, dissipative scales, as observed in the atmosphere and oceans, or in the solar environment. In this context, using direct numerical simulations with forcing, we analyze scale dynamics in the presence of magnetic fields with a generalized Ohm's law including a Hall current. The ion inertial length epsilon_H serves as the control parameter at fixed Reynolds number. Both the magnetic and generalized helicity -- invariants in the ideal case -- grow linearly with time, as expected from classical arguments. The cross-correlation between the velocity and magnetic field grows as well, more so in relative terms for a stronger Hall current. We find that the helical growth rates vary exponentially with epsilon_H, provided the ion inertial scale resides within the inverse cascade range. These exponential variations are recovered phenomenologically using simple scaling arguments. They are directly linked to the wavenumber power-law dependence of generalized and magnetic helicity, k^(-2), in their inverse ranges. This illustrates and confirms the important role of the interplay between large and small scales in the dynamics of turbulent flows.

preprint2020arXiv

Single-particle Lagrangian statistics from direct numerical simulations of rotating-stratified turbulence

Geophysical fluid flows are predominantly turbulent and often strongly affected by the Earth's rotation, as well as by stable density stratification. Using direct numerical simulations of forced Boussinesq equations, we study the influence of these effects on the motion of fluid particles, focusing on cases where the frequencies associated with rotation and stratification (RaS), $N$ and $f$ respectively, are held at a fixed ratio $N/f=5$. As the intensity of RaS increases, a sharp transition is observed between a regime dominated by eddies to a regime dominated by waves, which can also be seemingly described by simply comparing the time scale $1/N$ and $τ_η$ (the Kolmogorov time scale). We perform a detailed study of Lagrangian statistics of acceleration, velocity and related quantities in the two regimes. The flow anisotropy induces a clear difference between particle motion in the horizontal and vertical directions. In the regime $Nτ_η<1$, acceleration statistics in both horizontal and vertical directions, exhibit well known characteristics of isotropic turbulence. In contrast for $Nτ_η>1$, they are directly influenced by imposed RaS. The Lagrangian velocity statistics exhibit visible anisotropy for all runs; nevertheless the degree of anisotropy becomes very strong in the regime $Nτ_η>1$. We find that in the regime $Nτ_η<1$, rotation enhances the mean displacement of particles in horizontal planes at short times, but inhibits them at longer times. This inhibition of horizontal displacement becomes stronger for $Nτ_η>1$, with no clear diffusive behavior. Displacements in the vertical direction are always inhibited. The inhibition becomes extremely strong when $Nτ_η>1$, with the particles almost being trapped horizontally.

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

Rotating turbulence under "precession-like" perturbation

The effects of changing the orientation of the rotation axis on homogeneous turbulence is considered. We perform direct numerical simulations on a periodic box of $1024^3$ grid points, where the orientation of the rotation axis is changed (a) at a fixed time instant (b) regularly at time intervals commensurate with the rotation time scale. The former is characterized by a dominant inverse energy cascade whereas in the latter, the inverse cascade is stymied due to the recurrent changes in the rotation axis resulting in a strong forward energy transfer and large scale structures that resemble those of isotropic turbulence.

preprint2015arXiv

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

preprint2014arXiv

Restricted Equilibrium and the Energy Cascade in Rotating and Stratified Flows

Most of the turbulent flows appearing in nature (e.g. geophysical and astrophysical flows) are subjected to strong rotation and stratification. These effects break the symmetries of classical, homogenous isotropic turbulence. In doing so, they introduce a natural decomposition of phase space in terms of wave modes and potential vorticity modes. The appearance of a new time scale associated to the propagation of waves, in addition to the eddy turnover time, increases the complexity of the energy transfers between the various scales; nonlinearly interacting waves may dominate at some scales while balanced motion may prevail at others. In the end, it is difficult to predict \emph{a priori} if the energy cascades downscale as in homogeneous isotropic turbulence, upscale as expected from balanced dynamics, or follows yet another phenomenology. In this paper, we suggest a theoretical approach based on equilibrium statistical mechanics for the ideal system, inspired from the restricted partition function formalism introduced in metastability studies. In this framework, we show analytically that in the presence of rotation, when the dynamics is restricted to the slow modes, the equilibrium energy spectrum features an infrared divergence characteristic of an inverse cascade regime, whereas this is not the case for purely stratified flows.

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

A review of the possible role of constraints in MHD turbulence

A review of some of the issues that have arisen over the years concerning the energy distribution among scales for magnetohydrodynamics (MHD) turbulence is given here. A variety of tools are employed to that effect, and a central role is played by taking into consideration the ideal (non-dissipative) invariants, namely the total energy, the magnetic helicity and the cross-correlations between the velocity and the magnetic field (concentrating on the three-dimensional case). These concepts, based mostly on theory, models and direct numerical simulations, are briefly put in the context of observations, in particular the solar wind, and some of the remaining open questions are delineated as well. New results on ideal MHD dynamics in three dimensions on equivalent grids of up to $6144^3$ points using the Taylor-Green flow generalized to MHD are also mentioned.

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

Long-time properties of MHD turbulence and the role of symmetries

We investigate long-time properties of three-dimensional MHD turbulence in the absence of forcing and examine in particular the role played by the quadratic invariants of the system and by the symmetries of the initial configurations. We observe that, when sufficient accuracy is used, initial conditions with a high degree of symmetries, as in the absence of helicity, do not travel through parameter space over time whereas by perturbing these solutions either explicitly or implicitly using for example single precision for long times, the flows depart from their original behavior and can become either strongly helical, or have a strong alignment between the velocity and the magnetic field. When the symmetries are broken, the flows evolve towards different end states, as predicted by statistical arguments for non-dissipative systems with the addition of an energy minimization principle, as already analyzed in \cite{stribling_90} for random initial conditions using a moderate number of Fourier modes. Furthermore, the alignment properties of these flows, between velocity, vorticity, magnetic potential, induction and current, correspond to the dominance of two main regimes, one helically dominated and one in quasi-equipartition of kinetic and magnetic energy. We also contrast the scaling of the ratio of magnetic energy to kinetic energy as a function of wavenumber to the ratio of eddy turn-over time to Alfvén time as a function of wavenumber. We find that the former ratio is constant with an approximate equipartition for scales smaller than the largest scale of the flow whereas the ratio of time scales increases with increasing wavenumber.

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.

preprint2011arXiv

Alfvén waves and ideal two-dimensional Galerkin truncated magnetohydrodynamics

We investigate numerically the dynamics of two-dimensional Euler and ideal magnetohydrodynamics (MHD) flows in systems with a finite number of modes, up to $4096^2$, for which several quadratic invariants are preserved by the truncation and the statistical equilibria are known. Initial conditions are the Orszag-Tang vortex with a neutral X-point centered on a stagnation point of the velocity field in the large scales. In MHD, we observe that the total energy spectra at intermediate times and intermediate scales correspond to the interactions of eddies and waves, $E_T(k)\sim k^{-3/2}$. Moreover, no dissipative range is visible neither for Euler nor for MHD in two dimensions; in the former case, this may be linked to the existence of a vanishing turbulent viscosity whereas in MHD, the numerical resolution employed may be insufficient. When imposing a uniform magnetic field to the flow, we observe a lack of saturation of the formation of small scales together with a significant slowing-down of their equilibration, with however a cut-off independent partial thermalization being reached at intermediate scales.

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

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.

preprint2010arXiv

The effect of subfilter-scale physics on regularization models

The subfilter-scale (SFS) physics of regularization models are investigated to understand the regularizations' performance as SFS models. The strong suppression of spectrally local SFS interactions and the conservation of small-scale circulation in the Lagrangian-averaged Navier-Stokes alpha-model (LANS-alpha) is found to lead to the formation of rigid bodies. These contaminate the superfilter-scale energy spectrum with a scaling that approaches k^1 as the SFS spectra is resolved. The Clark-alpha and Leray-alpha models, truncations of LANS-alpha, do not conserve small-scale circulation and do not develop rigid bodies. LANS-alpha, however, is closest to Navier-Stokes in intermittency properties. All three models are found to be stable at high Reynolds number. Differences between L^2 and H^1 norm models are clarified. For magnetohydrodynamics (MHD), the presence of the Lorentz force as a source (or sink) for circulation and as a facilitator of both spectrally nonlocal large to small scale interactions as well as local SFS interactions prevents the formation of rigid bodies in Lagrangian-averaged MHD (LAMHD-alpha). We find LAMHD-alpha performs well as a predictor of superfilter-scale energy spectra and of intermittent current sheets at high Reynolds numbers. We expect it may prove to be a generally applicable MHD-LES.

preprint2008arXiv

Generalized Eulerian-Lagrangian description of Navier-Stokes and resistive MHD dynamics

New generalized equations of motion for the Weber-Clebsch potentials that describe both the Navier-Stokes and MHD dynamics are derived. These depend on a new parameter, which has dimensions of time for Navier-Stokes and inverse velocity for MHD. Direct numerical simulations are performed. For Navier-Stokes, the generalized formalism captures the intense reconnection of vortices of the Boratav, Pelz and Zabusky flow, in agreement with the previous study by Ohkitani and Constantin. For MHD, the new formalism is used to detect magnetic reconnection in several flows: the 3D Arnold, Beltrami and Childress (ABC) flow and the (2D and 3D) Orszag-Tang vortex. It is concluded that periods of intense activity in the magnetic enstrophy are correlated with periods of increasingly frequent resettings. Finally, the positive correlation between the sharpness of the increase in resetting frequency and the spatial localization of the reconnection region is discussed.

preprint2006arXiv

Dynamo action at low magnetic Prandtl numbers: mean flow vs. fully turbulent motion

We compute numerically the threshold for dynamo action in Taylor-Green swirling flows. Kinematic calculations, for which the flow field is fixed to its time averaged profile, are compared to dynamical runs for which both the Navier-Stokes and the induction equations are jointly solved. The kinematic instability is found to have two branches, for all explored Reynolds numbers. The dynamical dynamo threshold follows these branches: at low Reynolds number it lies within the low branch while at high kinetic Reynolds number it is close to the high branch.

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.