Source author record

A. Pouquet

A. 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

26works
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

26 published item(s)

preprint2016arXiv

Depletion of Nonlinearity in Magnetohydrodynamic Turbulence: Insights from Analysis and Simulations

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

preprint2014arXiv

Evidence for Bolgiano-Obukhov scaling in rotating stratified turbulence using high-resolution direct numerical simulations

We report results on rotating stratified turbulence in the absence of forcing, with large-scale isotropic initial conditions, using direct numerical simulations computed on grids of up to 4096^3 points. The Reynolds and Froude numbers are respectively equal to Re=5.4 x 10^4 and Fr=0.0242. The ratio of the Brunt-Väisälä to the inertial wave frequency, N/f, is taken to be equal to 4.95, a choice appropriate to model the dynamics of the southern abyssal ocean at mid latitudes. This gives a global buoyancy Reynolds number R_B=ReFr^2=32, a value sufficient for some isotropy to be recovered in the small scales beyond the Ozmidov scale, but still moderate enough that the intermediate scales where waves are prevalent are well resolved. We concentrate on the large-scale dynamics, for which we find a spectrum compatible with the Bolgiano-Obukhov scaling, and confirm that the Froude number based on a typical vertical length scale is of order unity, with strong gradients in the vertical. Two characteristic scales emerge from this computation, and are identified from sharp variations in the spectral distribution of either total energy or helicity. A spectral break is also observed at a scale at which the partition of energy between the kinetic and potential modes changes abruptly, and beyond which a Kolmogorov-like spectrum recovers. Large slanted layers are ubiquitous in the flow in the velocity and temperature fields, with local overturning events indicated by small Richardson numbers, and a small large-scale enhancement of energy directly attributable to the effect of rotation is also observed.

preprint2014arXiv

Forced MHD turbulence in three dimensions using Taylor-Green symmetries

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

preprint2014arXiv

Large-scale anisotropy in stably stratified rotating flows

We present results from direct numerical simulations of the Boussinesq equations in the presence of rotation and/or stratification, both in the vertical direction. The runs are forced isotropically and randomly at small scales and have spatial resolutions of up to $1024^3$ grid points and Reynolds numbers of $\approx 1000$. We first show that solutions with negative energy flux and inverse cascades develop in rotating turbulence, whether or not stratification is present. However, the purely stratified case is characterized instead by an early-time, highly anisotropic transfer to large scales with almost zero net isotropic energy flux. This is consistent with previous studies that observed the development of vertically sheared horizontal winds, although only at substantially later times. However, and unlike previous works, when sufficient scale separation is allowed between the forcing scale and the domain size, the total energy displays a perpendicular (horizontal) spectrum with power law behavior compatible with $\sim k_\perp^{-5/3}$, including in the absence of rotation. In this latter purely stratified case, such a spectrum is the result of a direct cascade of the energy contained in the large-scale horizontal wind, as is evidenced by a strong positive flux of energy in the parallel direction at all scales including the largest resolved scales.

preprint2014arXiv

Magnetic field reversals and long-time memory in conducting flows

Employing a simple ideal magnetohydrodynamic model in spherical geometry,we show that the presence of either rotation or finite magnetic helicity is sufficient to induce dynamical reversals of the magnetic dipole moment. The statistical character of the model is similar to that of terrestrial magnetic field reversals, with the similarity being stronger when rotation is present.The connection between long time correlations, $1/f$ noise, and statistics of reversals is supported, consistent with earlier suggestions.

preprint2014arXiv

Stably stratified turbulence in the presence of large-scale forcing

We perform two high resolution direct numerical simulations of stratified turbulence for Reynolds number equal to Re~25000 and Froude number respectively of Fr~0.1 and Fr~0.03. The flows are forced at large scale and discretized on an isotropic grid of 2048^3 points. Stratification makes the flow anisotropic and introduces two extra characteristic scales with respect to homogeneous isotropic turbulence: the buoyancy scale, L_B, and the Ozmidov scale, l_{oz}. The former is related to the number of layers that the flow develops in the direction of gravity, the latter is regarded as the scale at which isotropy is recovered. The values of L_B and l_{oz} depend on the Froude number and their absolute and relative size affect the repartition of energy among Fourier modes. By contrasting the behavior of the two simulated flows we identify some surprising similarities: after an initial transient the two flows evolve towards comparable values of the kinetic and potential enstrophy, and energy dissipation rate. Further similarities emerge at large scales: the same ratio between potential and total energy (~0.1) is spontaneously selected by the flows, and slow modes grow monotonically in both regimes causing a slow increase of the total energy in time. The axisymmetric total energy spectrum shows a wide variety of spectral slopes as a function of the angle between the imposed stratification and the wave vector. One-dimensional energy spectra computed in the direction parallel to gravity are flat from the forcing up to buoyancy scale. At intermediate scales a ~ k^{-3} parallel spectrum develops for the Fr ~ 0.03 run, whereas for weaker stratification, the saturation spectrum does not have enough scales to develop and instead one observes a power law compatible with Kolmogorov scaling. Finally, the spectrum of helicity is flat until L_B, as observed in the nocturnal planetary boundary layer.

preprint2013arXiv

Geophysical turbulence and the duality of the energy flow across scales

The ocean and the atmosphere, and hence the climate, are governed at large scale by interactions between pressure gradient, Coriolis and buoyancy forces. This leads to a quasi-geostrophic balance in which, in a two-dimensional-like fashion, the energy injected by solar radiation, winds or tides goes to large scales in what is known as an inverse cascade. Yet, except for Ekman friction, energy dissipation and turbulent mixing occur at small scale implying the formation of such scales associated with breaking of geostrophic dynamics through wave-eddy interactions \cite{ledwell_00, vanneste_13} or frontogenesis \cite{hoskins_72, mcwilliams_10}, in opposition to the inverse cascade. Can it be both at the same time? We exemplify here this dual behavior of energy with the help of three-dimensional direct numerical simulations of rotating stratified Boussinesq turbulence. We show that efficient small-scale mixing and large-scale coherence develop simultaneously in such geophysical and astrophysical flows, both with constant flux as required by theoretical arguments, thereby clearly resolving the aforementioned contradiction.

preprint2013arXiv

Ideal evolution of MHD turbulence when imposing Taylor-Green symmetries

We investigate the ideal and incompressible magnetohydrodynamic (MHD) equations in three space dimensions for the development of potentially singular structures. The methodology consists in implementing the four-fold symmetries of the Taylor-Green vortex generalized to MHD, leading to substantial computer time and memory savings at a given resolution; we also use a re-gridding method that allows for lower-resolution runs at early times, with no loss of spectral accuracy. One magnetic configuration is examined at an equivalent resolution of $6144^3$ points, and three different configurations on grids of $4096^3$ points. At the highest resolution, two different current and vorticity sheet systems are found to collide, producing two successive accelerations in the development of small scales. At the latest time, a convergence of magnetic field lines to the location of maximum current is probably leading locally to a strong bending and directional variability of such lines. A novel analytical method, based on sharp analysis inequalities, is used to assess the validity of the finite-time singularity scenario. This method allows one to rule out spurious singularities by evaluating the rate at which the logarithmic decrement of the analyticity-strip method goes to zero. The result is that the finite-time singularity scenario cannot be ruled out, and the singularity time could be somewhere between $t=2.33$ and $t=2.70.$ More robust conclusions will require higher resolution runs and grid-point interpolation measurements of maximum current and vorticity.

preprint2013arXiv

Inverse cascade behavior in freely decaying two-dimensional fluid turbulence

We present results from an ensemble of 50 runs of two-dimensional hydrodynamic turbulence with spatial resolution of 2048^2 grid points, and from an ensemble of 10 runs with 4096^2 grid points. All runs in each ensemble have random initial conditions with same initial integral scale, energy, enstrophy, and Reynolds number. When both ensemble- and time-averaged, inverse energy cascade behavior is observed, even in the absence of external mechanical forcing: the energy spectrum at scales larger than the characteristic scale of the flow follows a k^(-5/3) law, with negative flux, together with a k^(-3) law at smaller scales, and a positive flux of enstrophy. The source of energy for this behavior comes from the modal energy around the energy containing scale at t=0. The results shed some light into connections between decaying and forced turbulence, and into recent controversies in experimental studies of two-dimensional and magnetohydrodynamic turbulent flows.

preprint2013arXiv

Turbulence comes in bursts in stably stratified flows

There is a clear distinction between simple laminar and complex turbulent fluids. But in some cases, as for the nocturnal planetary boundary layer, a stable and well-ordered flow can develop intense and sporadic bursts of turbulent activity which disappear slowly in time. This phenomenon is ill-understood and poorly modeled; and yet, it is central to our understanding of weather and climate dynamics. We present here a simple model which shows that in stably stratified turbulence, the stronger bursts can occur when the flow is expected to be more stable. The bursts are generated by a rapid non-linear amplification of energy stored in waves, and are associated with energetic interchanges between vertical velocity and temperature (or density) fluctuations. Direct numerical simulations on grids of 2048^3 points confirm this somewhat paradoxical result of measurably stronger events for more stable flows, displayed not only in the temperature and vertical velocity derivatives, but also in the amplitude of the fields themselves.

preprint2012arXiv

Helicity dynamics in stratified turbulence in the absence of forcing

A numerical study of decaying stably-stratified flows is performed. Relatively high stratification and moderate Reynolds numbers are considered, and a particular emphasis is placed on the role of helicity (velocity-vorticity correlations). The problem is tackled by integrating the Boussinesq equations in a periodic cubical domain using different initial conditions: a non-helical Taylor-Green (TG) flow, a fully helical Beltrami (ABC) flow, and random flows with a tunable helicity. We show that for stratified ABC flows helicity undergoes a substantially slower decay than for unstratified ABC flows. This fact is likely associated to the combined effect of stratification and large scale coherent structures. Indeed, when the latter are missing, as in random flows, helicity is rapidly destroyed by the onset of gravitational waves. A type of large-scale dissipative "cyclostrophic" balance can be invoked to explain this behavior. When helicity survives in the system it strongly affects the temporal energy decay and the energy distribution among Fourier modes. We discover in fact that the decay rate of energy for stratified helical flows is much slower than for stratified non-helical flows and can be considered with a phenomenological model in a way similar to what is done for unstratified rotating flows. We also show that helicity, when strong, has a measurable effect on the Fourier spectra, in particular at scales larger than the buoyancy scale for which it displays a rather flat scaling associated with vertical shear.

preprint2012arXiv

Inverse cascades in turbulence and the case of rotating flows

We first summarize briefly several properties concerning the dynamics of two-dimensional (2D) turbulence, with an emphasis on the inverse cascade of energy to the largest accessible scale of the system. In order to study a similar phenomenon in three-dimensional (3D) turbulence undergoing strong solid-body rotation, we test a previously developed Large Eddy Simulation (LES) model against a high-resolution direct numerical simulation of rotating turbulence on a grid of $3072^3$ points. We then describe new numerical results on the inverse energy cascade in rotating flows using this LES model and contrast the case of 2D versus 3D forcing, as well as non-helical forcing (i.e., with weak overall alignment between velocity and vorticity) versus the fully helical Beltrami case, both for deterministic and random forcing. The different scaling of the inverse energy cascade can be attributed to the dimensionality of the forcing, with, in general, either a $k_{\perp}^{-3}$ or a $k_{\perp}^{-5/3}$ energy spectrum of slow modes at large scales, perpendicular referring to the direction of rotation. We finally invoke the role of shear in the case of a strongly anisotropic deterministic forcing, using the so-called ABC flow.

preprint2012arXiv

Isotropisation at small scales of rotating helically-driven turbulence

We present numerical evidence of how three-dimensionalization occurs at small scale in rotating turbulence with Beltrami (ABC) forcing, creating helical flow. The Zeman scale $\ell_Ω$ at which the inertial and eddy turn-over times are equal is more than one order of magnitude larger than the dissipation scale, with the relevant domains (large-scale inverse cascade of energy, dual regime in the direct cascade of energy $E$ and helicity $H$, and dissipation) each moderately resolved. These results stem from the analysis of a large direct numerical simulation on a grid of $3072^3$ points, with Rossby and Reynolds numbers respectively equal to 0.07 and $2.7\times 10^4$. At scales smaller than the forcing, a helical wave-modulated inertial law for the energy and helicity spectra is followed beyond $\ell_Ω$ by Kolmogorov spectra for $E$ and $H$. Looking at the two-dimensional slow manifold, we also show that the helicity spectrum breaks down at $\ell_Ω$, a clear sign of recovery of three-dimensionality in the small scales.

preprint2011arXiv

Conformal invariance in three-dimensional rotating turbulence

We examine three--dimensional turbulent flows in the presence of solid-body rotation and helical forcing in the framework of stochastic Schramm-Löwner evolution curves (SLE). The data stems from a run on a grid of $1536^3$ points, with Reynolds and Rossby numbers of respectively 5100 and 0.06. We average the parallel component of the vorticity in the direction parallel to that of rotation, and examine the resulting $<ω_\textrm{z}>_\textrm{z}$ field for scaling properties of its zero-value contours. We find for the first time for three-dimensional fluid turbulence evidence of nodal curves being conformal invariant, belonging to a SLE class with associated Brownian diffusivity $κ=3.6\pm 0.1$. SLE behavior is related to the self-similarity of the direct cascade of energy to small scales in this flow, and to the partial bi-dimensionalization of the flow because of rotation. We recover the value of $κ$ with a heuristic argument and show that this value is consistent with several non-trivial SLE predictions.

preprint2011arXiv

High Reynolds number magnetohydrodynamic turbulence using a Lagrangian model

With the help of a model of magnetohydrodynamic (MHD) turbulence tested previously, we explore high Reynolds number regimes up to equivalent resolutions of 6000^3 grid points in the absence of forcing and with no imposed uniform magnetic field. For the given initial condition chosen here, with equal kinetic and magnetic energy, the flow ends up being dominated by the magnetic field, and the dynamics leads to an isotropic Iroshnikov-Kraichnan energy spectrum. However, the locally anisotropic magnetic field fluctuations perpendicular to the local mean field follow a Kolmogorov law. We find that the ratio of the eddy turnover time to the Alfven time increases with wavenumber, contrary to the so-called critical balance hypothesis. Residual energy and helicity spectra are also considered; the role played by the conservation of magnetic helicity is studied, and scaling laws are found for the magnetic helicity and residual helicity spectra. We put these results in the context of the dynamics of a globally isotropic MHD flow which is locally anisotropic because of the influence of the strong large-scale magnetic field, leading to a partial equilibration between kinetic and magnetic modes for the energy and the helicity.

preprint2011arXiv

On the emergence of very long time fluctuations and 1/f noise in ideal flows

This study shows the connection between three previously observed but seemingly unrelated phenomena in hydrodynamic (HD) and magnetohydrodynamic (MHD) turbulent flows, involving the emergence of fluctuations occurring on very long time scales: the low-frequency 1/f noise in the power frequency spectrum, the delayed ergodicity of complex valued amplitude fluctuations in wavenumber space, and the spontaneous flippings or reversals of large scale fields. Direct numerical simulations of ideal MHD and HD are employed in three space dimensions, at low resolution, for long periods of time, and with high accuracy to study several cases: Different geometries, presence of rotation and/or a uniform magnetic field, and different values of the associated conserved global quantities. It is conjectured that the origin of all these long-time phenomena is rooted in the interaction of the longest wavelength fluctuations available to the system with fluctuations at much smaller scales. The strength of this non-local interaction is controlled either by the existence of conserved global quantities with a back-transfer in Fourier space, or by the presence of a slow manifold in the dynamics.

preprint2010arXiv

Experimental evidence of conformal invariance in soap film turbulent flows

We present experimental evidence of statistical conformal invariance in isocontours of fluid thickness in experiments of two-dimensional turbulence using soap films. A Schlieren technique is used to visualize regions of the flow with constant film thickness, and association of isocontours with Schramm-Löwner evolution (SLE) is used to identify conformal invariance. In experiments where an inverse energy cascade develops, statistical evidence is consistent with such an association. The diffusivity of the associated one-dimensional Brownian process is close to 8/3, a value previously identified in isocontours of vorticity in high-resolution numerical simulations of two-dimensional turbulence (D. Bernard et al., Nature Phys. 2, 124, 2006). In experiments where the inverse energy cascade is not sufficiently developed, no statistical evidence of conformal invariance is found.

preprint2010arXiv

Large scale behavior and statistical equilibria in rotating flows

We examine long-time properties of the ideal dynamics of three--dimensional flows, in the presence or not of an imposed solid-body rotation and with or without helicity (velocity-vorticity correlation). In all cases the results agree with the isotropic predictions stemming from statistical mechanics. No accumulation of excitation occurs in the large scales, even though in the dissipative rotating case anisotropy and accumulation, in the form of an inverse cascade of energy, are known to occur. We attribute this latter discrepancy to the linearity of the term responsible for the emergence of inertial waves. At intermediate times, inertial energy spectra emerge that differ somewhat from classical wave-turbulence expectations, and with a trace of large-scale excitation that goes away for long times. These results are discussed in the context of partial two-dimensionalization of the flow undergoing strong rotation as advocated by several authors.

preprint2009arXiv

Finite dissipation and intermittency in magnetohydrodynamics

We present an analysis of data stemming from numerical simulations of decaying magnetohydrodynamic (MHD) turbulence up to grid resolution of 1536^3 points and up to Taylor Reynolds number of 1200. The initial conditions are such that the initial velocity and magnetic fields are helical and in equipartition, while their correlation is negligible. Analyzing the data at the peak of dissipation, we show that the dissipation in MHD seems to asymptote to a constant as the Reynolds number increases, thereby strengthening the possibility of fast reconnection events in the solar environment for very large Reynolds numbers. Furthermore, intermittency of MHD flows, as determined by the spectrum of anomalous exponents of structure functions of the velocity and the magnetic field, is stronger than for fluids, confirming earlier results; however, we also find that there is a measurable difference between the exponents of the velocity and those of the magnetic field, as observed recently in the solar wind. Finally, we discuss the spectral scaling laws that arise in this flow.

preprint2009arXiv

Helical rotating turbulence. Part II. Intermittency, scale invariance and structures

We study the intermittency properties of the energy and helicity cascades in two 1536^3 direct numerical simulations of helical rotating turbulence. Symmetric and anti-symmetric velocity increments are examined, as well as probability density functions of the velocity field and of the helicity density. It is found that the direct cascade of energy to small scales is scale invariant and non-intermittent, whereas the direct cascade of helicity is highly intermittent. Furthermore, the study of structure functions of different orders allows us to identify a recovery of isotropy of strong events at very small scales in the flow. Finally, we observe the juxtaposition in space of strong laminar and persistent helical columns next to time-varying vortex tangles, the former being associated with the self-similarity of energy and the latter with the intermittency of helicity.

preprint2009arXiv

On the lack of universality in decaying magnetohydrodynamic turbulence

Using computations of three-dimensional magnetohydrodynamic (MHD) turbulence with a Taylor-Green flow, whose inherent time-independent symmetries are implemented numerically, and in the absence of either a forcing function or an imposed uniform magnetic field, we show that three different inertial ranges for the energy spectrum may emerge for three different initial magnetic fields, the selecting parameter being the ratio of the Alfven to the eddy turnover time. Equivalent computational grids range from 128^3 to 2048^3 points with a unit magnetic Prandtl number and a Taylor Reynolds number of up to 1500 at the peak of dissipation. We also show convergence of our results with Reynolds number. Our study is consistent with previous findings of a variety of energy spectra in MHD turbulence by studies performed in the presence of both a forcing term with a given correlation time and a strong, uniform magnetic field. In contrast to the previous studies, however, the ratio of characteristic time scales here can only be ascribed to the intrinsic nonlinear dynamics of the flows under study.

preprint2009arXiv

Rotating helical turbulence. Part I. Global evolution and spectral behavior

We present results from two 1536^3 direct numerical simulations of rotating turbulence where both energy and helicity are injected into the flow by an external forcing. The dual cascade of energy and helicity towards smaller scales observed in isotropic and homogeneous turbulence is broken in the presence of rotation, with the development of an inverse cascade of energy now coexisting with direct cascades of energy and helicity. In the direct cascade range, the flux of helicity dominates over that of energy at low Rossby number. These cascades have several consequences for the statistics of the flow. The evolution of global quantities and of the energy and helicity spectra is studied, and comparisons with simulations at different Reynolds and Rossby numbers at lower resolution are done to identify scaling laws.

preprint2009arXiv

The dynamics of unforced turbulence at high Reynolds number for Taylor-Green vortices generalized to MHD

We study decaying magnetohydrodynamics (MHD) turbulence stemming from the evolution of the Taylor-Green (TG) flow generalized recently to MHD, with equal viscosity and magnetic resistivity and up to equivalent grid resolutions of 2048^3 points. A pseudo-spectral code is used in which the symmetries of the velocity and magnetic fields have been implemented, allowing for sizable savings in both computer time and usage of memory at a given Reynolds number. The flow is non-helical, and at initial time the kinetic and magnetic energies are taken to be equal and concentrated in the large scales. After testing the validity of the method on grids of 512^3 points, we analyze the data on the large grids up to Taylor Reynolds numbers of ~2200. We find that the global temporal evolution is accelerated in MHD, compared to the corresponding neutral fluid case. We also observe an interval of time when such configurations have quasi-constant total dissipation, time during which statistical properties are determined after averaging over of the order of two turn-over times. A weak turbulence spectrum obtains which is also given in terms of its anisotropic components. Finally, we contrast the development of small-scale eddies with two other initial conditions for the magnetic field and briefly discuss the structures that develop, and which display a complex array of current and vorticity sheets with clear rolling-up and folding.

preprint2009arXiv

The interplay between helicity and rotation in turbulence: implications for scaling laws and small-scale dynamics

Invariance properties of physical systems govern their behavior: energy conservation in turbulence drives a wide distribution of energy among modes, observed in geophysical or astrophysical flows. In ideal hydrodynamics, the role of helicity conservation (correlation between velocity and its curl, measuring departures from mirror symmetry) remains unclear since it does not alter the energy spectrum. However, with solid body rotation, significant differences emerge between helical and non-helical flows. We first outline several results, like the energy and helicity spectral distribution and the breaking of strict universality for the individual spectra. Using massive numerical simulations, we then show that small-scale structures and their intermittency properties differ according to whether helicity is present or not, in particular with respect to the emergence of Beltrami-core vortices (BCV) that are laminar helical vertical updrafts. These results point to the discovery of a small parameter besides the Rossby number; this could relate the problem of rotating helical turbulence to that of critical phenomena, through renormalization group and weak turbulence theory. This parameter can be associated with the adimensionalized ratio of the energy to helicity flux to small scales, the three-dimensional energy cascade being weak and self-similar.

preprint2008arXiv

Cascades, thermalization and eddy viscosity in helical Galerkin truncated Euler flows

The dynamics of the truncated Euler equations with helical initial conditions are studied. Transient energy and helicity cascades leading to Kraichnan helical absolute equilibrium at small scales are obtained for the first time. The results of [Cichowlas et al. Phys. Rev. Lett. 95, 264502 (2005)] are extended to helical flows. Similarities between the turbulent transient evolution of the ideal (time-reversible) system and viscous helical flows are found. The observed differences in the behavior of truncated Euler and (constant viscosity) Navier-Stokes are qualitatively understood using the concept of eddy viscosity. The large scales of truncated Euler equations are then shown to follow quantitatively an effective Navier-Stokes dynamics based on a variable (scale dependent) eddy viscosity.

preprint2007arXiv

A comparison of spectral element and finite difference methods using statically refined nonconforming grids for the MHD island coalescence instability problem

A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.