Source author record

Roland Grappin

Roland Grappin 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

13works
5topics
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

13 published item(s)

preprint2016arXiv

Alfvén-dynamo balance and magnetic excess in MHD turbulence

3D Magnetohydrodynamic (MHD) turbulent flows with initially magnetic and kinetic energies at equipartition spontaneously develop a magnetic excess (or residual energy), as well in numerical simulations and in the solar wind. Closure equations obtained in 1983 describe the residual spectrum as being produced by a dynamo source proportional to the total energy spectrum, balanced by a linear Alfvén damping term. A good agreement was found in 2005 with incompressible simulations; however, recent solar wind measurements disagree with these results. The previous dynamo-Alfvén theory is generalized to a family of models, leading to simple relations between residual and total energy spectra. We want to assess these models in detail against MHD simulations and solar wind data. The family of models is tested against compressible decaying MHD simulations with low Mach number, low cross-helicity, zero mean magnetic field, without or with expansion terms (EBM or expanding box model). A single dynamo-Alfvén model is found to describe correctly both solar wind scalings and compressible simulations without or with expansion. It is equivalent to the 1983-2005 closure equation but with critical balance of nonlinear turnover and linear Alfvén times, while the dynamo source term remains unchanged. The discrepancy with previous incompressible simulations is elucidated. The model predicts a linear relation between the spectral slopes of total and residual energies $m_R = -1/2 + 3/2 m_T$. Examining the solar wind data as in \cite{2013ApJ...770..125C}, our relation is found to be valid whatever the cross-helicity, even better so at high cross-helicity, with the total energy slope varying from $1.7$ to $1.55$.

preprint2016arXiv

Beyond the Maltese cross: geometry of turbulence between 0.2 and 1 AU

The spectral anisotropy of turbulent structures has been measured in the solar wind since 1990, relying on the assumption of axisymmetry about the mean magnetic field, B0. However, several works indicate that this hypothesis might be partially wrong, thus raising two questions: (i) is it correct to interpret measurements at 1 AU (the so-called Maltese cross) in term of a sum of slab and 2D turbulence? (ii) what information is really contained in the Maltese cross? We solve direct numerical simulations of the MHD equations including the transverse stretching exerted by the solar wind flow and study the genuine 3D anisotropy of turbulence as well as that one resulting from the assumption of axisymmetry about B0. We show that the evolution of the turbulent spectrum from 0.2 to 1 AU depends strongly on its initial anisotropy. An axisymmetric spectrum with respect to B0 keeps its axisymmetry, i.e., resists stretching perpendicular to radial, while an isotropic spectrum becomes essentially axisymmetric with respect to the radial direction. We conclude that close to the Sun, slow-wind turbulence has a spectrum that is axisymmetric around B0 and the measured 2D component at 1 AU describes the real shape of turbulent structures. On the contrary, fast-wind turbulence has a more isotropic spectrum at the source and becomes radially symmetric at 1 AU. Such structure is hidden by the symmetrization applied to the data that instead returns a slab geometry.

preprint2015arXiv

Anisotropy of third-order structure functions in MHD turbulence

The measure of the third-order structure function, Y, is employed in the solar wind to compute the cascade rate of turbulence. In the absence of a mean field B0=0, Y is expected to be isotropic (radial) and independent of the direction of increments, so its measure yields directly the cascade rate. For turbulence with mean field, as in the solar wind, Y is expected to become more two dimensional (2D), that is, to have larger perpendicular components, loosing the above simple symmetry. To get the cascade rate one should compute the flux of Y, which is not feasible with single-spacecraft data, thus measurements rely upon assumptions about the unknown symmetry. We use direct numerical simulations (DNS) of magneto-hydrodynamic (MHD) turbulence to characterize the anisotropy of Y. We find that for strong guide field B0=5 the degree of two-dimensionalization depends on the relative importance of shear and pseudo polarizations (the two components of an Alfvén mode in incompressible MHD). The anisotropy also shows up in the inertial range. The more Y is 2D, the more the inertial range extent differs along parallel and perpendicular directions. We finally test the two methods employed in observations and find that the so-obtained cascade rate may depend on the angle between B0 and the direction of increments. Both methods yield a vanishing cascade rate along the parallel direction, contrary to observations, suggesting a weaker anisotropy of solar wind turbulence compared to our DNS. This could be due to a weaker mean field and/or to solar wind expansion.

preprint2015arXiv

Imprints of expansion onto the local anisotropy of solar wind turbulence

We study the anisotropy of II-order structure functions defined in a frame attached to the local mean field in three-dimensional (3D) direct numerical simulations of magnetohydrodynamic turbulence, including or not the solar wind expansion. We simulate spacecraft flybys through the numerical domain by taking increments along the radial (wind) direction that forms an angle of $45^o$ with the ambient magnetic field. We find that only when expansion is taken into account, do the synthetic observations match the 3D anisotropy observed in the solar wind, including the change of anisotropy with scales. Our simulations also show that the anisotropy changes dramatically when considering increments oblique to the radial directions. Both results can be understood by noting that expansion reduces the radial component of the magnetic field at all scales, thus confining fluctuations in the plane perpendicular to the radial. Expansion is thus shown to affect not only the (global) spectral anisotropy, but also the local anisotropy of second-order structure functions by influencing the distribution of the local mean field, which enters this higher-order statistics.

preprint2014arXiv

Evolution of turbulence in the expanding solar wind, a numerical study

We study the evolution of turbulence in the solar wind by solving numerically the full 3D magneto-hydrodynamic (MHD) equations embedded in a radial mean wind. The corresponding equations (expanding box model or EBM) have been considered earlier but never integrated in 3D simulations. Here, we follow the development of turbulence from 0.2 AU up to about 1.5 AU. Starting with isotropic spectra scaling as $k^{-1}$, we observe a steepening toward a $k^{-5/3}$ scaling in the middle of the wavenumber range and formation of spectral anisotropies. The advection of a plasma volume by the expanding solar wind causes a non-trivial stretching of the volume in directions transverse to radial and the selective decay of the components of velocity and magnetic fluctuations. These two effects combine to yield the following results. (i) Spectral anisotropy: gyrotropy is broken, and the radial wavevectors have most of the power. (ii) Coherent structures: radial streams emerge that resemble the observed microjets. (iii) Energy spectra per component: they show an ordering in good agreement with the one observed in the solar wind at 1 AU. The latter point includes a global dominance of the magnetic energy over kinetic energy in the inertial and $f^{-1}$ range and a dominance of the perpendicular-to-the-radial components over the radial components in the inertial range. We conclude that many of the above properties are the result of evolution during transport in the heliosphere, and not just the remnant of the initial turbulence close to the Sun.

preprint2014arXiv

Three-dimensional Iroshnikov-Kraichnan turbulence in a mean magnetic field

Forced, weak MHD turbulence with guide field is shown to adopt different regimes, depending on the magnetic excess of the large forced scales. When the magnetic excess is large enough, the classical perpendicular cascade with $5/3$ scaling is obtained, while when equipartition is imposed, an isotropic $3/2$ scaling appears in all directions with respect to the mean field (\cite{2010PhRvE..82b6406G} or GM10). We show here that the $3/2$ scaling of the GM10 regime is not ruled by a small-scale cross-helicity cascade, and propose that it is a 3D extension of a perpendicular weak Iroshnikov-Kraichnan (IK) cascade. We analyze in detail the structure functions in real space and show that they closely follow the critical balance relation both in the local frame and the global frame: we show that there is no contradiction between this and the isotropic $3/2$ scaling of the spectra. We propose a scenario explaining the spectral structure of the GM10 regime, that starts with a perpendicular weak IK cascade and extends to 3D by using quasi-resonant couplings. The quasi-resonance condition happens to reduce the energy flux in the same way as is done in the weak perpendicular cascade, so leading to a $3/2$ scaling in all directions. We discuss the possible applications of these findings to solar wind turbulence.

preprint2013arXiv

Rotational/ Compressional nature of the Magnetopause: application of the BV technique on a magnetopause case study

The magnetopause boundary implies two kinds of variations: a density/ temperature gradient and a magnetic field rotation. These two kinds are always observed in a close vicinity of each other, if not inseparably mixed. We present a case study from the Cluster data where the two are clearly separated and investigate the natures of both layers. We evidence that the first one is a slow shock while the second is a rotational discontinuity. The interaction between these two kinds of discontinuities is then studied with the help of 1,5-D magnetohydrodynamics simulations. The comparison with the data is quite positive and leads to think that most of the generic properties of the magnetopause may be interpreted in this sense.

preprint2012arXiv

Quasi-isotropic cascade in MHD turbulence with mean field

We propose a phenomenological theory of incompressible magnetohydrodynamic turbulence in the presence of a strong large-scale magnetic field, which establishes a link between the known anisotropic models of strong and weak MHD turbulence We argue that the Iroshnikov-Kraichnan isotropic cascade develops naturally within the plane perpendicular to the mean field, while oblique-parallel cascades with weaker amplitudes can develop, triggered by the perpendicular cascade, with a reduced flux resulting from a quasi-resonance condition. The resulting energy spectrum $E(k_\parallel,k_\bot)$ has the same slope in all directions. The ratio between the extents of the inertial range in the parallel and perpendicular directions is equal to $b_{rms}/B_0$. These properties match those found in recent 3D MHD simulations with isotropic forcing reported in [R. Grappin and W.-C. Müller, Phys. Rev. E \textbf{82}, 26406 (2010)].

preprint2012arXiv

Solar wind turbulent spectrum at plasma kinetic scales

The description of the turbulent spectrum of magnetic fluctuations in the solar wind in the kinetic range of scales is not yet completely established. Here, we perform a statistical study of 100 spectra measured by the STAFF instrument on the Cluster mission, which allows to resolve turbulent fluctuations from ion scales down to a fraction of electron scales, i.e. from $\sim 10^2$ km to $\sim 300$ m. We show that for $k_{\perp}ρ_e \in[0.03,3]$ (that corresponds approximately to the frequency in the spacecraft frame $f\in [3,300]$ Hz), all the observed spectra can be described by a general law $E(k_\perp)\propto k_\perp^{-8/3}\exp{(-k_\perp ρ_e)}$, where $k_{\perp}$ is the wave-vector component normal to the background magnetic field and $ρ_e$ the electron Larmor radius. This exponential tail found in the solar wind seems compatible with the Landau damping of magnetic fluctuations onto electrons.

preprint2012arXiv

Transition from weak to strong cascade in MHD turbulence

The transition from weak to strong turbulence when passing from large to small scales in magnetohydrodynamic (MHD) turbulence with guide field is a cornerstone of anisotropic turbulence theory. We present the first check of this transition, using the Shell-RMHD which combines a shell model of perpendicular nonlinear coupling and linear propagation along the guide field. This model allows us to reach Reynolds numbers around $10^6$. We obtain surprisingly good agreement with the theoretical predictions, with a reduced perpendicular energy spectrum scaling as $k^{-2}$ at large scales and as $k^{-5/3}$ at small scales, where critical balance between nonlinear and propagation time is reached. However, even in the strong regime, a high level of excitation is found in the weak coupling region of Fourier space, which is due to the rich frequency spectrum of large eddies. A corollary is that the reduced parallel spectral slope is not a definite test of the spectral anisotropy, contrary to standard belief.

preprint2011arXiv

Coronal heating in coupled photosphere-chromosphere-coronal systems: turbulence and leakage

Coronal loops act as resonant cavities for low frequency fluctuations that are transmitted from the deeper layers of the solar atmosphere and are amplified in the corona, triggering nonlinear interactions. However trapping is not perfect, some energy leaks down to the chromosphere, thus limiting the turbulence development and the associated heating. We consider the combined effects of turbulence and leakage in determining the energy level and associated heating rate in models of coronal loops which include the chromosphere and transition region. We use a piece-wise constant model for the Alfven speed and a Reduced MHD - Shell model to describe the interplay between turbulent dynamics in the direction perpendicular to the mean field and propagation along the field. Turbulence is sustained by incoming fluctuations which are equivalent, in the line-tied case, to forcing by the photospheric shear flows. While varying the turbulence strength, we compare systematically the average coronal energy level (E) and dissipation rate (D) in three models with increasing complexity: the classical closed model, the semi-open corona model, and the corona-chromosphere (or 3-layer) model, the latter two models allowing energy leakage. We find that: (i) Leakage always plays a role (even for strong turbulence), E and D are systematically lower than in the line-tied model. (ii) E is close to the resonant prediction, i.e., assuming effective turbulent correlation time longer than the Alfven coronal crossing time (Ta). (iii) D is close to the value given by the ratio of photospheric energy divided by Ta (iv) The coronal spectra exibits an inertial range with 5/3 spectral slope, and a large scale peak of trapped resonant modes that inhibit nonlinear couplings. (v) In the realistic 3-layer model, the two-component spectrum leads to a damping time equal to the Kolmogorov time reduced by a factor u_rms/Va_corona

preprint2011arXiv

Coupling the solar dynamo and the corona: wind properties, mass and momentum losses during an activity cycle

We study the connections between the sun's convection zone and the evolution of the solar wind and corona. We let the magnetic fields generated by a 2.5D axisymmetric kinematic dynamo code (STELEM) evolve in a 2.5D axisymmetric coronal isothermal MHD code (DIP). The computations cover an 11 year activity cycle. The solar wind's asymptotic velocity varies in latitude and in time in good agreement with the available observations. The magnetic polarity reversal happens at different paces at different coronal heights. Overall sun's mass loss rate, momentum flux and magnetic braking torque vary considerably throughout the cycle. This cyclic modulation is determined by the latitudinal distribution of the sources of open flux and solar wind and the geometry of the Alfvén surface. Wind sources and braking torque application zones also vary accordingly.

preprint2010arXiv

Scaling and anisotropy of magnetohydrodynamic turbulence in a strong mean magnetic field

We present a new analysis of the anisotropic spectral energy distribution in incompressible magnetohydrodynamic (MHD) turbulence permeated by a strong mean magnetic field. The turbulent flow is generated by high-resolution pseudo-spectral direct numerical simulations with large-scale isotropic forcing. Examining the radial energy distribution for various angles $θ$ with respect to $\mathbf{B}_0$ reveals a specific structure which remains hidden when not taking axial symmetry with respect to $B_0$ into account. For each direction, starting at the forced large-scales, the spectrum first exhibits an amplitude drop around a wavenumber $k_0$ which marks the start of a scaling range and goes on up to a dissipative wavenumber $k_d(θ)$. The 3D spectrum for $k \ge k_0$ is described by a single $θ$-independent functional form $F(k/k_d)$, the scaling law being the same in every direction. The previous properties still hold when increasing the mean field from $B_0=5$ up to $B_0=10 \ b_{rms}$, as well as when passing from resistive to ideal flows. We conjecture that at fixed $B_0$ the direction-independent scaling regime is reached when increasing the Reynolds number above a threshold which raises with increasing $B_0$. Below that threshold critically balanced turbulence is expected.