Showing 1–44 of 44 results for author: Verdini, A

Higher-Order Analysis of Three-Dimensional Anisotropy in Imbalanced Alfvénic Turbulence

Authors: Nikos Sioulas, Themistocles Zikopoulos, Chen Shi, Marco Velli, Trevor Bowen, Alfred Mallet, Luca Sorriso-Valvo, Andrea Verdini, B. D. G. Chandran, Mihailo M. Martinović, S. S. Cerri, Nooshin Davis, Corina Dunn

Abstract: We analyze in-situ observations of imbalanced solar wind turbulence to evaluate MHD turbulence models grounded in "Critical Balance" (CB) and "Scale-Dependent Dynamic Alignment" (SDDA). At energy injection scales, both outgoing and ingoing modes exhibit a weak cascade; a simultaneous tightening of SDDA is noted. Outgoing modes persist in a weak cascade across the inertial range, while ingoing mode… ▽ More We analyze in-situ observations of imbalanced solar wind turbulence to evaluate MHD turbulence models grounded in "Critical Balance" (CB) and "Scale-Dependent Dynamic Alignment" (SDDA). At energy injection scales, both outgoing and ingoing modes exhibit a weak cascade; a simultaneous tightening of SDDA is noted. Outgoing modes persist in a weak cascade across the inertial range, while ingoing modes shift to a strong cascade at $λ\approx 3 \times 10^{4} d_i$, with associated spectral scalings deviating from expected behavior due to "anomalous coherence" effects. The inertial range comprises two distinct sub-inertial segments. Beyond $λ\gtrsim 100 d_i$, eddies adopt a field-aligned tube topology, with SDDA signatures mainly evident in high amplitude fluctuations. The scaling exponents $ζ_{n}$ of the $n$-th order conditional structure functions, orthogonal to both the local mean field and fluctuation direction, align with the analytical models of Chandran et al. 2015 and Mallet et al. 2017, indicating "multifractal" statistics and strong intermittency; however, scaling in parallel and displacement components is more concave than predicted, possibly influenced by expansion effects. Below $λ\approx 100 d_i$, eddies become increasingly anisotropic, evolving into thin current sheet-like structures. Concurrently, $ζ_{n}$ scales linearly with order, marking a shift towards "monofractal" statistics. At $λ\approx 8 d_i$, the increase in aspect ratio halts, and the eddies become quasi-isotropic. This change may signal tearing instability, leading to reconnection, or result from energy redirection into the ion-cyclotron wave spectrum, aligning with the "helicity barrier". Our analysis utilizes 5-point structure functions, proving more effective than the traditional 2-point method in capturing steep scaling behaviors at smaller scales. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Temperature inversion in a confined plasma atmosphere: coarse-grained effect of temperature fluctuations at its base

Authors: Luca Barbieri, Emanuele Papini, Pierfrancesco Di Cintio, Simone Landi, Andrea Verdini, Lapo Casetti

Abstract: Prompted by the relevant problem of temperature inversion (i.e. gradient of density anti-correlated to the gradient of temperature) in astrophysics, we introduce a novel method to model a gravitationally confined multi-component collisionless plasma in contact with a fluctuating thermal boundary. We focus on systems with anti-correlated (inverted) density and temperature profiles, with application… ▽ More Prompted by the relevant problem of temperature inversion (i.e. gradient of density anti-correlated to the gradient of temperature) in astrophysics, we introduce a novel method to model a gravitationally confined multi-component collisionless plasma in contact with a fluctuating thermal boundary. We focus on systems with anti-correlated (inverted) density and temperature profiles, with applications to solar physics. The dynamics of the plasma is analytically described via the coupling of an appropriated coarse-grained distribution function and temporally coarse-grained Vlasov dynamics. We derive a stationary solution of the system and predict the inverted density and temperature profiles of the two-species for scenarios relevant for the corona. We validate our method by comparing the analytical results with kinetic numerical simulations of the plasma dynamics in the context of the two-species Hamiltonian mean-field model (HMF). Finally, we apply our theoretical framework to the problem of the temperature inversion in the solar corona obtaining density and temperature profiles in remarkably good agreement with the observations. △ Less

Submitted 18 June, 2024; v1 submitted 19 January, 2024; originally announced January 2024.

Comments: 27 pages, 12 Figs. Version accepted for publication in JPP

Temperature inversion in a gravitationally bound plasma: Case of the solar corona

Authors: Luca Barbieri, Lapo Casetti, Andrea Verdini, Simone Landi

Abstract: The temperature of the solar atmosphere increases from thousands to millions of degrees moving from the lower layer (chromosphere) to the outermost one (corona), while the density drops accordingly. The mechanism behind this phenomenon, known as a temperature inversion, is still unknown. In this work, we model a coronal loop as a collisionless plasma confined in a semicircular tube that is subject… ▽ More The temperature of the solar atmosphere increases from thousands to millions of degrees moving from the lower layer (chromosphere) to the outermost one (corona), while the density drops accordingly. The mechanism behind this phenomenon, known as a temperature inversion, is still unknown. In this work, we model a coronal loop as a collisionless plasma confined in a semicircular tube that is subject to the Sun's gravity and in thermal contact with a fully collisional chromosphere behaving as a thermostat at the loop's feet. By using kinetic $N$-particle simulations and analytical calculations, we show that rapid, intermittent, and short-lived heating events in the chromosphere drive the coronal plasma towards a non-equilibrium stationary state. The latter is characterized by suprathermal tails in the particles' velocity distribution functions, exhibiting temperature and density profiles strikingly similar to those observed in the atmosphere of the Sun. These results suggest that a million-Kelvin solar corona can be produced without the local deposition of heat in the upper layer of the atmosphere that is typically assumed by standard approaches. We find that suprathermal distribution functions in the corona are self-consistently produced instead of postulated a priori, in contrast to classical kinetic models based on a velocity filtration mechanism. △ Less

Submitted 20 December, 2023; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 6 pages, 3 figures, aa LATeX class. Accepted for publication in Astronomy and Astrophysics as a Letter

Journal ref: Astronomy & Astrophysics 681, L5 (2024)

The essential role of multi-point measurements in investigations of turbulence, three-dimensional structure, and dynamics: the solar wind beyond single scale and the Taylor Hypothesis

Authors: W. H. Matthaeus, S. Adhikari, R. Bandyopadhyay, M. R. Brown, R. Bruno, J. Borovsky, V. Carbone, D. Caprioli, A. Chasapis, R. Chhiber, S. Dasso, P. Dmitruk, L. Del Zanna, P. A. Dmitruk, Luca Franci, S. P. Gary, M. L. Goldstein, D. Gomez, A. Greco, T. S. Horbury, Hantao Ji, J. C. Kasper, K. G. Klein, S. Landi, Hui Li , et al. (27 additional authors not shown)

Abstract: Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer h… ▽ More Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer heliosphere. This shortcoming leaves numerous central questions incompletely answered. Deficiencies remain in at least two important subjects, Space Weather and fundamental plasma turbulence theory, due to a lack of a more complete understanding of the space-time structure of dynamic plasmas. Only with multispacecraft measurements over suitable spans of spatial separation and temporal duration can these ambiguities be resolved. We note that these characterizations apply to turbulence across a wide range of scales, and also equally well to shocks, flux ropes, magnetic clouds, current sheets, stream interactions, etc. In the following, we will describe the basic requirements for resolving space-time structure in general, using turbulence' as both an example and a principal target or study. Several types of missions are suggested to resolve space-time structure throughout the Heliosphere. △ Less

Submitted 26 November, 2022; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: White Paper submitted to: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. arXiv admin note: substantial text overlap with arXiv:1903.06890

Ion-scale transition of plasma turbulence: Pressure-strain effect

Authors: Petr Hellinger, Victor Montagud-Camps, Luca Franci, Lorenzo Matteini, Emanuele Papini, Andrea Verdini, Simone Landi

Abstract: We investigate properties of solar wind-like plasma turbulence using direct numerical simulations. We analyze the transition from large, magnetohydrodynamic (MHD) scales to the ion characteristic ones using two-dimensional hybrid (fluid electrons, kinetic ions) simulations. To capture and quantify turbulence properties, we apply the Karman-Howarth-Monin (KHM) equation for compressible Hall MHD (ex… ▽ More We investigate properties of solar wind-like plasma turbulence using direct numerical simulations. We analyze the transition from large, magnetohydrodynamic (MHD) scales to the ion characteristic ones using two-dimensional hybrid (fluid electrons, kinetic ions) simulations. To capture and quantify turbulence properties, we apply the Karman-Howarth-Monin (KHM) equation for compressible Hall MHD (extended by considering the plasma pressure as a tensor quantity) to the numerical results. The KHM analysis indicates that the transition from MHD to ion scales (the so called ion break in the power spectrum) results from a combination of an onset of Hall physics and of an effective dissipation owing to the pressure-strain energy-exchange channel and resistivity. We discuss the simulation results in the context of the solar wind. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 7 pages, 3 figures

Spectral transfer and Kármán-Howarth-Monin equations for compressible Hall magnetohydrodynamics

Authors: Petr Hellinger, Emanuele Papini, Andrea Verdini, Simone Landi, Luca Franci, Lorenzo Matteini, Victor Montagud-Camps

Abstract: We derive two new forms of the Kármán-Howarth-Monin equation for decaying compressible Hall magnetohydrodynamic (MHD) turbulence. We test them on results of a weakly-compressible, two-dimensional, moderate-Reynolds-number Hall MHD simulation and compare them with an isotropic spectral transfer (ST) equation. The KHM and ST equations are automatically satisfied during the whole simulation owing to… ▽ More We derive two new forms of the Kármán-Howarth-Monin equation for decaying compressible Hall magnetohydrodynamic (MHD) turbulence. We test them on results of a weakly-compressible, two-dimensional, moderate-Reynolds-number Hall MHD simulation and compare them with an isotropic spectral transfer (ST) equation. The KHM and ST equations are automatically satisfied during the whole simulation owing to the periodic boundary conditions and have complementary cumulative behavior. They are used here to analyze the onset of turbulence and its properties when it is fully developed. These approaches give equivalent results characterizing: the decay of the kinetic + magnetic energy at large scales, the MHD and Hall cross-scale energy transfer/cascade, the pressure dilatation, and the dissipation. The Hall cascade appears when the MHD one brings the energy close to the ion inertial range and is related to the formation of reconnecting current sheets. At later times, the pressure-dilation energy-exchange rate oscillates around zero with no net effect on the cross-scale energy transfer when averaged over a period of its oscillations. A reduced one-dimensional analysis suggests that all three methods may be useful to estimate the energy cascade rate from in situ observations. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 11 pages, 7 figues

Investigating the origin of the FIP effect with a shell turbulence model

Authors: Victor Réville, Alexis P. Rouillard, Marco Velli, Andrea Verdini, Éric Buchlin, Michaël Lavarra, Nicolas Poirier

Abstract: The enrichment of coronal loops and the slow solar wind with elements that have low First Ionisation Potential, known as the FIP effect, has often been interpreted as the tracer of a common origin. A current explanation for this FIP fractionation rests on the influence of ponderomotive forces and turbulent mixing acting at the top of the chromosphere. The implied wave transport and turbulence mech… ▽ More The enrichment of coronal loops and the slow solar wind with elements that have low First Ionisation Potential, known as the FIP effect, has often been interpreted as the tracer of a common origin. A current explanation for this FIP fractionation rests on the influence of ponderomotive forces and turbulent mixing acting at the top of the chromosphere. The implied wave transport and turbulence mechanisms are also key to wave-driven coronal heating and solar wind acceleration models. This work makes use of a shell turbulence model run on open and closed magnetic field lines of the solar corona to investigate with a unified approach the influence of magnetic topology, turbulence amplitude and dissipation on the FIP fractionation. We try in particular to assess whether there is a clear distinction between the FIP effect on closed and open field regions. △ Less

Submitted 5 January, 2021; originally announced January 2021.

Comments: 10 pages, 6 figures, accepted for publication in Frontiers in Astronomy and Space Sciences

The Solar Orbiter Science Activity Plan: translating solar and heliospheric physics questions into action

Authors: I. Zouganelis, A. De Groof, A. P. Walsh, D. R. Williams, D. Mueller, O. C. St Cyr, F. Auchere, D. Berghmans, A. Fludra, T. S. Horbury, R. A. Howard, S. Krucker, M. Maksimovic, C. J. Owen, J. Rodriiguez-Pacheco, M. Romoli, S. K. Solanki, C. Watson, L. Sanchez, J. Lefort, P. Osuna, H. R. Gilbert, T. Nieves-Chinchilla, L. Abbo, O. Alexandrova , et al. (160 additional authors not shown)

Abstract: Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operat… ▽ More Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operations are essential to address the following four top-level science questions: (1) What drives the solar wind and where does the coronal magnetic field originate? (2) How do solar transients drive heliospheric variability? (3) How do solar eruptions produce energetic particle radiation that fills the heliosphere? (4) How does the solar dynamo work and drive connections between the Sun and the heliosphere? Maximising the mission's science return requires considering the characteristics of each orbit, including the relative position of the spacecraft to Earth (affecting downlink rates), trajectory events (such as gravitational assist manoeuvres), and the phase of the solar activity cycle. Furthermore, since each orbit's science telemetry will be downloaded over the course of the following orbit, science operations must be planned at mission level, rather than at the level of individual orbits. It is important to explore the way in which those science questions are translated into an actual plan of observations that fits into the mission, thus ensuring that no opportunities are missed. First, the overarching goals are broken down into specific, answerable questions along with the required observations and the so-called Science Activity Plan (SAP) is developed to achieve this. The SAP groups objectives that require similar observations into Solar Orbiter Observing Plans (SOOPs), resulting in a strategic, top-level view of the optimal opportunities for science observations during the mission lifetime. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 20 pages, 1 figure, accepted by Astronomy & Astrophysics

Journal ref: A&A 642, A3 (2020)

Comparing turbulent cascades and heating vs spectral anisotropy in solar wind via direct simulations

Authors: Victor Montagud-Camps, Roland Grappin, Andrea Verdini

Abstract: In a previous work (MGV18), we showed numerically that the turbulent cascade generated by quasi-2D structures (with wave vectors mostly-perpendicular to the mean magnetic field) is able to generate a temperature profile close to the one observed in solar wind ($\simeq 1/R$) in the range 0.2 $\le R \le$ 1 au. Theory, observations and numerical simulations point to another robust structure, the radi… ▽ More In a previous work (MGV18), we showed numerically that the turbulent cascade generated by quasi-2D structures (with wave vectors mostly-perpendicular to the mean magnetic field) is able to generate a temperature profile close to the one observed in solar wind ($\simeq 1/R$) in the range 0.2 $\le R \le$ 1 au. Theory, observations and numerical simulations point to another robust structure, the radial-slab, with dominant wave vectors along the radial: we study here the efficiency of the radial-slab cascade in building the $1/R$ temperature profile. As in MGV18, we solve the three-dimensional MHD equations including expansion to simulate the turbulent evolution. We find that an isotropic distribution of wave vectors with large cross helicity at 0.2 au, along with a large wind expansion rate, lead again to a temperature decay rate close to $1/R$ but with a radial-slab anisotropy at 1 au. Surprisingly, the turbulent cascade concentrates in the plane transverse to the radial direction, displaying 1D spectra with scalings close to $k^{-5/3}$ in this plane. This supports both the idea of turbulent heating of the solar wind, and the existence of two different turbulent cascades, quasi-2D and radial slab, at the origin of the heating. We conclude that sampling the radial spectrum in the solar wind may give but a poor information on the real cascade regime and rate when the radial slab is a non-negligible part of turbulence. △ Less

Submitted 31 August, 2020; originally announced August 2020.

Comments: 13 pages, 15 figures. Submitted to ApJ

Modeling Kelvin-Helmholtz instability-driven turbulence with hybrid simulations of Alfvénic turbulence

Authors: Luca Franci, Julia E. Stawarz, Emanuele Papini, Petr Hellinger, Takuma Nakamura, David Burgess, Simone Landi, Andrea Verdini, Lorenzo Matteini, Robert Ergun, Olivier Le Contel, Per-Arne Lindqvist

Abstract: Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfvénic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and… ▽ More Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfvénic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfvén wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath. △ Less

Submitted 17 November, 2019; originally announced November 2019.

Comments: 11 pages, 6 figures, submitted to The Astrophysical Journal

Turbulence vs. fire hose instabilities: 3-D hybrid expanding box simulations

Authors: P. Hellinger, L. Matteini, S. Landi, L. Franci, A. Verdini, E. Papini

Abstract: The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-hel… ▽ More The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-helicity. A turbulent cascade rapidly develops and leads to a weak proton heating that is not sufficient to overcome the expansion-driven perpendicular cooling. The plasma system eventually drives the parallel and oblique fire hose instabilities that generate quasi-monochromatic wave packets that reduce the proton temperature anisotropy. The fire hose wave activity has a low amplitude with wave vectors quasi-parallel/oblique with respect to the ambient magnetic field outside of the region dominated by the turbulent cascade and is discernible in one-dimensional power spectra taken only in the direction quasi-parallel/oblique with respect to the ambient magnetic field; at quasi-perpendicular angles the wave activity is hidden by the turbulent background. These waves are partly reabsorbed by protons and partly couple to and participate in the turbulent cascade. Their presence reduces kurtosis, a measure of intermittency, and the Shannon entropy but increases the Jensen-Shannon complexity of magnetic fluctuations; these changes are weak and anisotropic with respect to the ambient magnetic field and it's not clear if they can be used to indirectly discern the presence of instability-driven waves. △ Less

Submitted 21 August, 2019; originally announced August 2019.

Comments: 12 pages, 10 figures

Three-dimensional local anisotropy of velocity fluctuations in the solar wind

Authors: Andrea Verdini, Roland Grappin, Olga Alexandrova, Luca Franci, Simone Landi, Lorenzo Matteini, Emanuele Papini

Abstract: We analyse velocity fluctuations in the solar wind at magneto-fluid scales in two datasets, extracted from Wind data in the period 2005-2015, that are characterised by strong or weak expansion. Expansion affects measurements of anisotropy because it breaks axisymmetry around the mean magnetic field. Indeed, the small-scale three-dimensional local anisotropy of magnetic fluctuations (δB) as measure… ▽ More We analyse velocity fluctuations in the solar wind at magneto-fluid scales in two datasets, extracted from Wind data in the period 2005-2015, that are characterised by strong or weak expansion. Expansion affects measurements of anisotropy because it breaks axisymmetry around the mean magnetic field. Indeed, the small-scale three-dimensional local anisotropy of magnetic fluctuations (δB) as measured by structure functions (SF_B) is consistent with tube-like structures for strong expansion. When passing to weak expansion, structures become ribbon-like because of the flattening of SFB along one of the two perpendicular directions. The power-law index that is consistent with a spectral slope -5/3 for strong expansion now becomes closer to -3/2. This index is also characteristic of velocity fluctuations in the solar wind. We study velocity fluctuations (δV) to understand if the anisotropy of their structure functions (SF_V ) also changes with the strength of expansion and if the difference with the magnetic spectral index is washed out once anisotropy is accounted for. We find that SF_V is generally flatter than SF_B. When expansion passes from strong to weak, a further flattening of the perpendicular SF_V occurs and the small-scale anisotropy switches from tube-like to ribbon-like structures. These two types of anisotropy, common to SF_V and SF_B, are associated to distinct large-scale variance anisotropies of δB in the strong- and weak-expansion datasets. We conclude that SF_V shows anisotropic three-dimensional scaling similar to SF_B, with however systematic flatter scalings, reflecting the difference between global spectral slopes. △ Less

Submitted 8 April, 2019; originally announced April 2019.

Comments: accepted in MNRAS

Spectral anisotropies and intermittency of plasma turbulence at ion kinetic scales

Authors: Simone Landi, Luca Franci, Emanuele Papini, Andrea Verdini, Lorenzo Matteini, Petr Hellinger

Abstract: By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories a… ▽ More By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories although the non linear energy transfer is still in the perpendicular direction, only advected in the parallel direction as expected balancing the non-linear energy transfer time and the decorrelation time. Such result can be explained by a phenomenological model based on the formation of strong intermittent two-dimensional structures in the plane perpendicular to the local mean field that have some prescribed aspect ratio eventually depending on the scale. This model support the idea that small scales structures, such as reconnecting current sheets, contribute significantly to the formation of the turbulent cascade at kinetic scales. △ Less

Submitted 17 April, 2019; v1 submitted 8 April, 2019; originally announced April 2019.

Comments: 9 pages, 4 figures

[Plasma 2020 Decadal] The essential role of multi-point measurements in turbulence investigations: the solar wind beyond single scale and beyond the Taylor Hypothesis

Authors: W. H. Matthaeus, R. Bandyopadhyay, M. R. Brown, J. Borovsky, V. Carbone, D. Caprioli, A. Chasapis, R. Chhiber, S. Dasso, P. Dmitruk, L. Del Zanna, P. A. Dmitruk, Luca Franci, S. P. Gary, M. L. Goldstein, D. Gomez, A. Greco, T. S. Horbury, Hantao Ji, J. C. Kasper, K. G. Klein, S. Landi, Hui Li, F. Malara, B. A. Maruca , et al. (24 additional authors not shown)

Abstract: This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for su… ▽ More This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence. △ Less

Submitted 16 March, 2019; originally announced March 2019.

Comments: White paper submitted to the PLASMA 2020 Decadal Survey Committee

Von Karman-Howarth equation for Hall magnetohydrodynamics: Hybrid simulations

Authors: Petr Hellinger, Andrea Verdini, Simone Landi, Luca Franci, Lorenzo Matteini

Abstract: A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is… ▽ More A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is valid. In the sub-ion range the cascade continues via the Hall term but the dynamic law derived in the framework of incompressible Hall MHD equations is obtained only in a low plasma beta simulation. For a higher beta plasma the cascade rate decreases in the sub-ion range and the change becomes more pronounced as the plasma beta increases. This break in the cascade flux can be ascribed to non thermal (kinetic) features or to others terms in the dynamical equation that are not included in the Hall-MHD incompressible approximation. △ Less

Submitted 26 March, 2018; originally announced March 2018.

Comments: 3 Figures, 5 pages, ApJL, submitted

3D Anisotropy of Solar Wind Turbulence, Tubes, or Ribbons?

Authors: Andrea Verdini, Roland Grappin, Olga Alexandrova, Sonny Lion

Abstract: We study the anisotropy with respect to the local magnetic field of turbulent magnetic fluctuations at magnetofluid scales in the solar wind. Previous measurements in the fast solar wind obtained axisymmetric anisotropy, despite that the analysis method allows nonaxisymmetric structures. These results are probably contaminated by the wind expansion that introduces another symmetry axis, namely, th… ▽ More We study the anisotropy with respect to the local magnetic field of turbulent magnetic fluctuations at magnetofluid scales in the solar wind. Previous measurements in the fast solar wind obtained axisymmetric anisotropy, despite that the analysis method allows nonaxisymmetric structures. These results are probably contaminated by the wind expansion that introduces another symmetry axis, namely, the radial direction, as indicated by recent numerical simulations. These simulations also show that while the expansion is strong, the principal fluctuations are in the plane perpendicular to the radial direction. Using this property, we separate 11 yr of Wind spacecraft data into two subsets characterized by strong and weak expansion and determine the corresponding turbulence anisotropy. Under strong expansion, the small-scale anisotropy is consistent with the Goldreich & Sridhar critical balance. As in previous works, when the radial symmetry axis is not eliminated, the turbulent structures are field-aligned tubes. Under weak expansion, we find 3D anisotropy predicted by the Boldyrev model, that is, turbulent structures are ribbons and not tubes. However, the very basis of the Boldyrev phenomenology, namely, a cross-helicity increasing at small scales, is not observed in the solar wind: the origin of the ribbon formation is unknown. △ Less

Submitted 30 January, 2019; v1 submitted 27 February, 2018; originally announced February 2018.

Comments: Published version with ERRATUM (fig7) attached

Journal ref: 2018ApJ...853...85V

Three-dimensional simulations of solar wind turbulence with the hybrid code CAMELIA

Authors: L. Franci, P. Hellinger, M. Guarrasi, C. H. K. Chen, E. Papini, A. Verdini, L. Matteini, S. Landi

Abstract: We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particle-in-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good ag… ▽ More We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particle-in-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good agreement with previous two-dimensional (2D) HPIC simulations, especially in the kinetic range of scales, and reproduce several features observed in solar wind spectra. By providing scaling tests on many different architectures and convergence studies, we prove CAMELIA to represent a very efficient, accurate and reliable tool for investigating the develpoment of the turbulent cascade in the solar wind, being able to cover simultaneously several decades in wavenumber, also in 3D. △ Less

Submitted 11 December, 2017; originally announced December 2017.

Comments: 13 pages, 10 figures, submitted to the Proceedings of the 12th International Conference on Numerical Modeling of Space Plasma Flows (ASTRONUM-2017)

Solar wind turbulent cascade from MHD to sub-ion scales: large-size 3D hybrid particle-in-cell simulations

Authors: Luca Franci, Simone Landi, Andrea Verdini, Lorenzo Matteini, Petr Hellinger

Abstract: Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfvènic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic fiel… ▽ More Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfvènic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic field. The omnidirectional magnetic field spectrum shows a double power-law behavior over almost two decades in wavenumber, with a Kolmogorov-like index at large scales, a spectral break around ion scales, and a steepening at sub-ion scales. Power laws are also observed in the spectra of the ion bulk velocity, density, and electric field, both at magnetohydrodynamic (MHD) and at kinetic scales. Despite the complex structure, the omnidirectional spectra of all fields at ion and sub-ion scales are in remarkable quantitative agreement with those of a two-dimensional (2D) simulation with similar physical parameters. This provides a partial, a-posteriori validation of the 2D approximation at kinetic scales. Conversely, at MHD scales, the spectra of the density and of the velocity (and, consequently, of the electric field) exhibit differences between the 2D and 3D cases. Although they can be partly ascribed to the lower spatial resolution, the main reason is likely the larger importance of compressible effects in a full geometry. Our findings are also in remarkable quantitative agreement with solar wind observations. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 11 pages, 8 figures, submitted to The Astrophysical Journal

Sub-structure formation in starless cores

Authors: C. Toci, D. Galli, A. Verdini, L. Del Zanna, S. Landi

Abstract: Motivated by recent observational searches of sub-structure in starless molecular cloud cores, we investigate the evolution of density perturbations on scales smaller than the Jeans length embedded in contracting isothermal clouds, adopting the same formalism developed for the expanding Universe and the solar wind. We find that initially small amplitude, Jeans-stable perturbations (propagating as… ▽ More Motivated by recent observational searches of sub-structure in starless molecular cloud cores, we investigate the evolution of density perturbations on scales smaller than the Jeans length embedded in contracting isothermal clouds, adopting the same formalism developed for the expanding Universe and the solar wind. We find that initially small amplitude, Jeans-stable perturbations (propagating as sound waves in the absence of a magnetic field), are amplified adiabatically during the contraction, approximately conserving the wave action density, until they either become nonlinear and steepen into shocks at a time $t_{\rm nl}$, or become gravitationally unstable when the Jeans length decreases below the scale of the perturbations at a time $t_{\rm gr}$. We evaluate analytically the time $t_{\rm nl}$ at which the perturbations enter the non-linear stage using a Burgers' equation approach, and we verify numerically that this time marks the beginning of the phase of rapid dissipation of the kinetic energy of the perturbations. We then show that for typical values of the rms Mach number in molecular cloud cores, $t_{\rm nl}$ is smaller than $t_{\rm gr}$, and therefore density perturbations likely dissipate before becoming gravitational unstable. Solenoidal modes grow at a faster rate than compressible modes, and may eventually promote fragmentation through the formation of vortical structures. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: 8 pages, 4 figures

Magnetic reconnection as a driver for a sub-ion scale cascade in plasma turbulence

Authors: L. Franci, S. S. Cerri, F. Califano, S. Landi, E. Papini, A. Verdini, L. Matteini, F. Jenko, P. Hellinger

Abstract: A new path for the generation of a sub-ion scale cascade in collisionless plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of power-l… ▽ More A new path for the generation of a sub-ion scale cascade in collisionless plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of power-law energy spectra below the so-called ion break occurs as soon as the first magnetic reconnection events take place, regardless of the actual state of the turbulent cascade at MHD scales. In both simulations, the reconnection-mediated small-scale energy spectrum of parallel magnetic fluctuations exhibits a very stable spectral slope of ~-2.8, whether or not a large-scale turbulent cascade has already fully developed. Once a quasi-stationary turbulent state is achieved, the spectrum of the total magnetic fluctuations settles towards a spectral index of -5/3 in the MHD range and of ~-3 at sub-ion scales. △ Less

Submitted 19 July, 2017; originally announced July 2017.

Comments: 6 pages, 5 figures

Mirror instability in the turbulent solar wind

Authors: P. Hellinger, S. Landi, L. Matteini, A. Verdini, L. Franci

Abstract: The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition betwe… ▽ More The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at sub-ion scales. The imposed expansion (taking a strictly transverse ambient magnetic field) leads to generation of an important perpendicular proton temperature anisotropy that eventually drives the mirror instability. This instability generates large-amplitude, nonpropagating, compressible, pressure-balanced magnetic structures in a form of magnetic enhancements/humps that reduce the perpendicular temperature anisotropy. △ Less

Submitted 17 March, 2017; originally announced March 2017.

Comments: ApJ, 7 pages, 9 figures. arXiv admin note: text overlap with arXiv:1508.03159

Plasma beta dependence of the ion-scale spectral break of solar wind turbulence: high-resolution 2D hybrid simulations

Authors: Luca Franci, Simone Landi, Lorenzo Matteini, Andrea Verdini, Petr Hellinger

Abstract: We investigate properties of the ion-scale spectral break of solar wind turbulence by means of two-dimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, di… ▽ More We investigate properties of the ion-scale spectral break of solar wind turbulence by means of two-dimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, distributed over three orders of magnitude, from 0.01 to 10. In all the cases, once turbulence is fully developed, we observe a power-law spectrum of the fluctuating magnetic field on large scales (in the inertial range) with a spectral index close to -5/3, while in the sub-ion range we observe another power-law spectrum with a spectral index systematically varying with $β$ (from around -3.6 for small values to around -2.9 for large ones). The two ranges are separated by a spectral break around ion scales. The length scale at which this transition occurs is found to be proportional to the ion inertial length, $d_i$, for $β\ll 1$ and to the ion gyroradius, $ρ_i = d_i\sqrtβ$, for $β\gg 1$, i.e., to the larger between the two scales in both the extreme regimes. For intermediate cases, i.e., $β\sim 1$, a combination of the two scales is involved. We infer an empiric relation for the dependency of the spectral break on $β$ that provides a good fit over the whole range of values. We compare our results with in situ observations in the solar wind and suggest possible explanations for such a behavior. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: 8 pages, 4 figures, accepted for publication on The Astrophysical Journal

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

Authors: Andrea Verdini, Roland Grappin

Abstract: 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 turbule… ▽ More 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. △ Less

Submitted 28 September, 2016; originally announced September 2016.

Comments: 8 pages, 8 figures, accepted for publication in ApJ

Two-dimensional Hybrid Simulations of Kinetic Plasma Turbulence: Current and Vorticity vs Proton Temperature

Authors: Luca Franci, Petr Hellinger, Lorenzo Matteini, Andrea Verdini, Simone Landi

Abstract: Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate… ▽ More Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate that, although regions where the proton temperature is enhanced and temperature anisotropies are larger correspond approximately to regions where many thin current sheets form, no firm quantitative evidence supports the idea of a direct causality between the two phenomena. On the other hand, quite a clear correlation between the behavior of the proton temperature and the out-of-plane vorticity is obtained. △ Less

Submitted 4 April, 2016; originally announced April 2016.

Comments: 4 pages, 2 figures, Proceedings of the Fourteenth International Solar Wind Conference

Journal ref: AIP Conf. Proc. 1720, 040003 (2016)

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

Authors: Roland Grappin, Wolf-Christian Müller, Andrea Verdini

Abstract: 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 d… ▽ More 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$. △ Less

Submitted 11 March, 2016; originally announced March 2016.

Comments: 7 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 589, A131 (2016)

Plasma turbulence and kinetic instabilities at ion scales in the expanding solar wind

Authors: Petr Hellinger, Lorenzo Matteini, Simone Landi, Andrea Verdini, Luca Franci, Pavel M. Travnicek

Abstract: The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvénic fluctuations which have energy equipartition… ▽ More The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvénic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at ion scales. The turbulent cascade leads to an overall anisotropic proton heating, protons are heated in the perpendicular direction, and, initially, also in the parallel direction. The imposed expansion leads to generation of a large parallel proton temperature anisotropy which is at later stages partly reduced by turbulence. The turbulent heating is not sufficient to overcome the expansion-driven perpendicular cooling and the system eventually drives the oblique firehose instability in a form of localized nonlinear wave packets which efficiently reduce the parallel temperature anisotropy. This work demonstrates that kinetic instabilities may coexist with strong plasma turbulence even in a constrained 2-D regime. △ Less

Submitted 29 September, 2015; v1 submitted 13 August, 2015; originally announced August 2015.

Comments: 6 pages, 5 figures, in press in ApJL

Journal ref: Astrophys. J. Lett. (2015), 811, L32

High-resolution hybrid simulations of kinetic plasma turbulence at proton scales

Authors: Luca Franci, Simone Landi, Lorenzo Matteini, Andrea Verdini, Petr Hellinger

Abstract: We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbu… ▽ More We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe a MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of -5/3 and -3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ~-3. A -2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell and the resolution employed, the perpendicular proton temperature strongly depends on these parameters. △ Less

Submitted 19 June, 2015; originally announced June 2015.

Comments: 15 pages, 13 figures, submitted to ApJ

Imprints of expansion onto the local anisotropy of solar wind turbulence

Authors: Andrea Verdini, Roland Grappin

Abstract: 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 ambien… ▽ More 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. △ Less

Submitted 10 June, 2015; originally announced June 2015.

Comments: 5 pages, 5 figures, accepted in ApJL

Solar wind turbulence from MHD to sub-ion scales: high-resolution hybrid simulations

Authors: Luca Franci, Andrea Verdini, Lorenzo Matteini, Simone Landi, Petr Hellinger

Abstract: We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several propert… ▽ More We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm's law. △ Less

Submitted 18 April, 2015; v1 submitted 18 March, 2015; originally announced March 2015.

Comments: 5 pages, 5 figures; introduction and conclusions changed, references updated, accepted for publication in ApJL

Anisotropy of third-order structure functions in MHD turbulence

Authors: Andrea Verdini, Roland Grappin, Petr Hellinger, Simone Landi, Wolf-Christian Müller

Abstract: 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 (… ▽ More 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. △ Less

Submitted 2 April, 2015; v1 submitted 16 February, 2015; originally announced February 2015.

Comments: Some text editing and typos corrected, 13 pages, 6 figures, to be published in ApJ

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

Authors: Yue Dong, Andrea Verdini, Roland Grappin

Abstract: 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 scali… ▽ More 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. △ Less

Submitted 29 August, 2014; originally announced September 2014.

Comments: 14 pages, 12 figures, accepted in ApJ

Parametric decay of parallel and oblique Alfven waves in the expanding solar wind

Authors: L. Del Zanna, L. Matteini, S. Landi, A. Verdini, M. Velli

Abstract: The long-term evolution of large-amplitude Alfven waves propagating in the solar wind is investigated by performing two-dimensional MHD simulations within the expanding box model. The linear and nonlinear phases of the parametric decay instability are studied for both circularly polarized waves in parallel propagation and for arc-polarized waves in oblique propagation. The non-monochromatic case i… ▽ More The long-term evolution of large-amplitude Alfven waves propagating in the solar wind is investigated by performing two-dimensional MHD simulations within the expanding box model. The linear and nonlinear phases of the parametric decay instability are studied for both circularly polarized waves in parallel propagation and for arc-polarized waves in oblique propagation. The non-monochromatic case is also considered. In the oblique case, the direct excitation of daughter modes transverse to the local background field is found for the first time in an expanding environment, and this transverse cascade seems to be favored for monochromatic mother waves. The expansion effect reduces the instability growth rate, and it can even suppress its onset for the lowest frequency modes considered here, possibly explaining the persistence of these outgoing waves in the solar wind. △ Less

Submitted 22 July, 2014; originally announced July 2014.

Comments: 20 pages, 9 figures, accepted for publication in Joural of Plasma Physics

Journal ref: J. Plasma Phys. 81 (2015) 325810102

Validating a Time-Dependent Wave-Turbulence-Driven Model of the Solar Wind

Authors: Roberto Lionello, Marco Velli, Cooper Downs, Jon A. Linker, Zoran Mikić, Andrea Verdini

Abstract: Although the mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still being actively investigated, it is largely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process. \citet{2010ApJ...708L.116V} presented a model for heating and accelerating the solar wind based on the turbulent dissipation… ▽ More Although the mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still being actively investigated, it is largely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process. \citet{2010ApJ...708L.116V} presented a model for heating and accelerating the solar wind based on the turbulent dissipation of Alfvén waves. We first use a time-dependent model of the solar wind to reproduce one of \citeauthor{2010ApJ...708L.116V}'s solutions; then we extend its application to the case when the energy equation includes thermal conduction and radiation losses, and the upper chromosphere is part of the computational domain. Using this model, we explore parameter space and describe the characteristics of a fast-solar-wind solution. We discuss how this formulation may be applied to a 3D MHD model of the coron a and solar wind \citep{2009ApJ...690..902L}. △ Less

Submitted 17 February, 2014; originally announced February 2014.

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

Authors: Roland Grappin, Wolf-Christian Müller, Andrea Verdini, Özgür Gürcan

Abstract: 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… ▽ More 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. △ Less

Submitted 13 December, 2014; v1 submitted 12 December, 2013; originally announced December 2013.

Comments: Major re-write of manuscript

Coupling the solar surface and the corona: coronal rotation, Alfvén wave-driven polar plumes

Authors: R. F. Pinto, R. Grappin, M. Velli, A. Verdini

Abstract: The dynamical response of the solar corona to surface and sub-surface perturbations depends on the chromospheric stratification, and specifically on how efficiently these layers reflect or transmit incoming Alfvén waves. While it would be desirable to include the chromospheric layers in the numerical simulations used to study such phenomena, that is most often not feasible. We defined and tested a… ▽ More The dynamical response of the solar corona to surface and sub-surface perturbations depends on the chromospheric stratification, and specifically on how efficiently these layers reflect or transmit incoming Alfvén waves. While it would be desirable to include the chromospheric layers in the numerical simulations used to study such phenomena, that is most often not feasible. We defined and tested a simple approximation allowing the study of coronal phenomena while taking into account a parametrised chromospheric reflectivity. We addressed the problems of the transmission of the surface rotation to the corona and that of the generation of polar plumes by Alfvén waves (Pinto et al., 2010, 2011). We found that a high (yet partial) effective chromospheric reflectivity is required to properly describe the angular momentum balance in the corona and the way the surface differential rotation is transmitted upwards. Alfvén wave-driven polar plumes maintain their properties for a wide range of values for the reflectivity, but they become bursty (and eventually disrupt) when the limit of total reflection is attained. △ Less

Submitted 2 January, 2013; originally announced January 2013.

Comments: Solar Wind 13: Proceedings of the Thirteenth International Solar Wind Conference

Reflection and dissipation of Alfvén waves in interstellar clouds

Authors: C. Pinto, A. Verdini, D. Galli, M. Velli

Abstract: Context: Supersonic nonthermal motions in molecular clouds are often interpreted as long-lived magnetohydrodynamic (MHD) waves. The propagation and amplitude of these waves is affected by local physical characteristics, most importantly the gas density and the ionization fraction. Aims: We study the propagation, reflection and dissipation of Alfvén waves in molecular clouds deriving the behavior… ▽ More Context: Supersonic nonthermal motions in molecular clouds are often interpreted as long-lived magnetohydrodynamic (MHD) waves. The propagation and amplitude of these waves is affected by local physical characteristics, most importantly the gas density and the ionization fraction. Aims: We study the propagation, reflection and dissipation of Alfvén waves in molecular clouds deriving the behavior of observable quantities such as the amplitudes of velocity fluctuations and the rate of energy dissipation. Methods: We formulated the problem in terms of Elsässer variables for transverse MHD waves propagating in a one-dimensional inhomogeneous medium, including the dissipation due to collisions between ions and neutrals and to a nonlinear turbulent cascade treated in a phenomenological way. We considered both steady-state and time-dependent situations and solved the equations of the problem numerically with an iterative method and a Lax-Wendroff scheme, respectively. Results: Alfvén waves incident on overdense regions with density profiles typical of cloud cores embedded in a diffuse gas suffer enhanced reflection in the regions of the steepest density gradient, and strong dissipation in the core's interior. These effects are especially significant when the wavelength is intermediate between the critical wavelength for propagation and the typical scale of the density gradient. For larger wave amplitudes and/or steeper input spectra, the effects of the perpendicular turbulent cascade result in a stronger energy dissipation in the regions immediately surrounding the dense core. Conclusions: The results may help to interpret the sharp decrease of line width observed in the environments of low-mass cloud cores in several molecular transitions. △ Less

Submitted 18 July, 2012; originally announced July 2012.

Comments: 9 pages, 6 figures

Transition from weak to strong cascade in MHD turbulence

Authors: Andrea Verdini, Roland Grappin

Abstract: 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 rea… ▽ More 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. △ Less

Submitted 29 June, 2012; v1 submitted 10 April, 2012; originally announced April 2012.

Comments: 4 pages, 5 figures, accepted in PRL

Time delays in quasi-periodic pulsations observed during the X2.2 solar flare on 2011 February 15

Authors: L. Dolla, C. Marqué, D. B. Seaton, T. Van Doorsselaere, M. Dominique, D. Berghmans, C. Cabanas, A. De Groof, W. Schmutz, A. Verdini, M. J. West, J. Zender, A. N. Zhukov

Abstract: We report observations of quasi-periodic pulsations (QPPs) during the X2.2 flare of 2011 February 15, observed simultaneously in several wavebands. We focus on fluctuations on time scale 1-30 s and find different time lags between different wavebands. During the impulsive phase, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) channels in the range 25-100 keV lead all the other ch… ▽ More We report observations of quasi-periodic pulsations (QPPs) during the X2.2 flare of 2011 February 15, observed simultaneously in several wavebands. We focus on fluctuations on time scale 1-30 s and find different time lags between different wavebands. During the impulsive phase, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) channels in the range 25-100 keV lead all the other channels. They are followed by the Nobeyama RadioPolarimeters at 9 and 17 GHz and the Extreme Ultra-Violet (EUV) channels of the Euv SpectroPhotometer (ESP) onboard the Solar Dynamic Observatory (SDO). The Zirconium and Aluminum filter channels of the Large Yield Radiometer (LYRA) onboard the Project for On-Board Autonomy (PROBA2) satellite and the SXR channel of ESP follow. The largest lags occur in observations from the Geostationary Operational Environmental Satellite (GOES), where the channel at 1-8 Å leads the 0.5-4 Å channel by several seconds. The time lags between the first and last channels is up to 9 s. We identified at least two distinct time intervals during the flare impulsive phase, during which the QPPs were associated with two different sources in the Nobeyama RadioHeliograph at 17 GHz. The radio as well as the hard X-ray channels showed different lags during these two intervals. To our knowledge, this is the first time that time lags are reported between EUV and SXR fluctuations on these time scales. We discuss possible emission mechanisms and interpretations, including flare electron trapping. △ Less

Submitted 28 March, 2012; originally announced March 2012.

Journal ref: ApJ 2012, 749, L16

On the origin of the 1/f spectrum in the solar wind magnetic field

Authors: A. Verdini, R. Grappin, R. Pinto, M. Velli

Abstract: We present a mechanism for the formation of the low frequency 1/f magnetic spectrum based on numerical solutions of a shell reduced-MHD model of the turbulent dynamics inside the sub-Alfvénic solar wind. We assign reasonably realistic profiles to the wind speed and the density along the radial direction, and a radial magnetic field. Alfvén waves of short periodicity (600 s) are injected at the bas… ▽ More We present a mechanism for the formation of the low frequency 1/f magnetic spectrum based on numerical solutions of a shell reduced-MHD model of the turbulent dynamics inside the sub-Alfvénic solar wind. We assign reasonably realistic profiles to the wind speed and the density along the radial direction, and a radial magnetic field. Alfvén waves of short periodicity (600 s) are injected at the base of the chromosphere, penetrate into the corona and are partially reflected, thus triggering a turbulent cascade. The cascade is strong for the reflected wave while it is weak for the outward propagating waves. Reflection at the transition region recycles the strong turbulent spectrum into the outward weak spectrum, which is advected beyond the Alfvénic critical point without substantial evolution. There, the magnetic field has a perpendicular power-law spectrum with slope close to the Kolmogorov -5/3. The parallel spectrum is inherited from the frequency spectrum of large (perpendicular) eddies. The shape is a double power-law with slopes of -1 and -2 at low and high frequencies respectively, the position of the break depending on the injected spectrum. We suggest that the double power-law spectrum measured by Helios at 0.3 AU, where the average magnetic field is not aligned with the radial (contrary to our assumptions) results from the combination of such different spectral slopes. At low frequency the parallel spectrum dominates with its characteristic 1/f shape, while at higher frequencies its steep spectral slope (-2) is masked by the more energetic perpendicular spectrum (slope -5/3). △ Less

Submitted 9 April, 2012; v1 submitted 28 March, 2012; originally announced March 2012.

Comments: 5 pages, 4 figures, accepted for publication in ApJL, V2: typo corrected in eq.1, color figures

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

Authors: Andrea Verdini, Roland Grappin, Marco Velli

Abstract: 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 a… ▽ More 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 △ Less

Submitted 22 November, 2011; originally announced November 2011.

Comments: 15 pages, 15 figures, Accepted for publication in A&A

Coupling Photosphere and Corona: Linear and Turbulent Regimes

Authors: A. Verdini, R. Grappin, M. Velli

Abstract: In a recent work Grappin et al. [1] have shown that low- frequency movements can be transmitted from one footpoint to the other along a magnetic loop, thus mimicking a friction effect of the corona on the photosphere, and invalidating the line-tying approximation. We consider here successively the effect of high frequencies and turbulent damping on the process. We use a very simple atmospheric m… ▽ More In a recent work Grappin et al. [1] have shown that low- frequency movements can be transmitted from one footpoint to the other along a magnetic loop, thus mimicking a friction effect of the corona on the photosphere, and invalidating the line-tying approximation. We consider here successively the effect of high frequencies and turbulent damping on the process. We use a very simple atmospheric model which allows to study analytically the laminar case, and to study the turbulent case both using simple phenomenological arguments and a more sophisticated turbulence model [2]. We find that, except when turbulent damping is such that all turbulence is damped during loop traversal, coupling still occurs between distant footpoints, and moreover the coronal field induced by photospheric movements saturates at finite values. △ Less

Submitted 8 December, 2009; originally announced December 2009.

Comments: proceedings of the Solar Wind 12 conference (St. Malo)

Journal ref: AIP Conf.Proc.1216:28-31,2010

A turbulence-driven model for heating and acceleration of the fast wind in coronal holes

Authors: A. Verdini, M. Velli, W. H. Matthaeus, S. Oughton, P. Dmitruk

Abstract: A model is presented for generation of fast solar wind in coronal holes, relying on heating that is dominated by turbulent dissipation of MHD fluctuations transported upwards in the solar atmosphere. Scale-separated transport equations include large-scale fields, transverse Alfvenic fluctuations, and a small compressive dissipation due to parallel shears near the transition region. The model acc… ▽ More A model is presented for generation of fast solar wind in coronal holes, relying on heating that is dominated by turbulent dissipation of MHD fluctuations transported upwards in the solar atmosphere. Scale-separated transport equations include large-scale fields, transverse Alfvenic fluctuations, and a small compressive dissipation due to parallel shears near the transition region. The model accounts for proton temperature, density, wind speed, and fluctuation amplitude as observed in remote sensing and in situ satellite data. △ Less

Submitted 27 November, 2009; originally announced November 2009.

Comments: accepted for publication in ApJL

Turbulence in the sub-Alfvénic solar wind driven by reflection of low-frequency Alfvén waves

Authors: A. Verdini, M. Velli, E. Buchlin

Abstract: We study the formation and evolution of a turbulent spectrum of Alfvén waves driven by reflection off the solar wind density gradients, starting from the coronal base up to 17 solar radii, well beyond the Alfvénic critical point. The background solar wind is assigned and 2D shell models are used to describe nonlinear interactions. We find that the turbulent spectra are influenced by the nature o… ▽ More We study the formation and evolution of a turbulent spectrum of Alfvén waves driven by reflection off the solar wind density gradients, starting from the coronal base up to 17 solar radii, well beyond the Alfvénic critical point. The background solar wind is assigned and 2D shell models are used to describe nonlinear interactions. We find that the turbulent spectra are influenced by the nature of reflected waves. Close to the base, these give rise to a flatter and steeper spectrum for the outgoing and reflected waves respectively. At higher heliocentric distance both spectra evolve toward an asymptotic Kolmogorov spectrum. The turbulent dissipation is found to account for at least half of the heating required to sustain the background imposed solar wind and its shape is found to be determined by the reflection-determined turbulent heating below 1.5 solar radii. Therefore reflection and reflection-driven turbulence are shown to play a key role in the accelerationof the fast solar wind and origin of the turbulent spectrum found at 0.3 AU in the heliosphere. △ Less

Submitted 15 May, 2009; originally announced May 2009.

Comments: Accepted for publication in ApJL

Journal ref: Astrophys.J.700:L39-L42,2009

Alfven Waves and Turbulence in the Solar Atmosphere and Solar Wind

Authors: A. Verdini, M. Velli

Abstract: We solve the problem of propagation and dissipation of Alfvenic turbulence in a model solar atmosphere consisting of a static photosphere and chromosphere, transition region, and open corona and solar wind, using a phenomenological model for the turbulent dissipation based on wave reflection. We show that most of the dissipation for a given wave-frequency spectrum occurs in the lower corona, and… ▽ More We solve the problem of propagation and dissipation of Alfvenic turbulence in a model solar atmosphere consisting of a static photosphere and chromosphere, transition region, and open corona and solar wind, using a phenomenological model for the turbulent dissipation based on wave reflection. We show that most of the dissipation for a given wave-frequency spectrum occurs in the lower corona, and the overall rms amplitude of the fluctuations evolves in a way consistent with observations. The frequency spectrum, for a Kolmogorov-like slope, is not found to change dramatically from the photosphere to the solar wind, however it does preserve signatures of transmission throughout the lower atmospheric layers, namely oscillations in the spectrum at high frequencies reminiscent of the resonances found in the linear case. These may disappear once more realistic couplings for the non-linear terms are introduced, or if time-dependent variability of the lower atmospheric layer is introduced. △ Less

Submitted 7 February, 2007; originally announced February 2007.

Comments: 27 pages, 9 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.662:669-676,2007