Source author record

M. Velli

M. Velli 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

33works
8topics
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

33 published item(s)

preprint2022arXiv

Features of magnetic field switchbacks in relation to the local-field geometry of large-amplitude Alfvénic oscillations: \emph{Wind} and \emph{PSP} observations

In this letter we report observations of magnetic switchback (SB) features near 1 au using data from the \emph{Wind} spacecraft. These features appear to be strikingly similar to the ones observed by the Parker Solar Probe mission (PSP) closer to the Sun: namely, one-sided spikes (or enhancements) in the solar-wind bulk speed $V$ that correlate/anti-correlate with the spikes seen in the radial-field component $B_R$. In the solar-wind streams that we analyzed, these specific SB features near 1 au are associated with large-amplitude Alfvénic oscillations that propagate outward from the sun along a local background (prevalent) magnetic field $\bf{B}_0$ that is nearly radial. We also show that, when $\bf{B}_0$ is nearly perpendicular to the radial direction, the large amplitude Alfvénic oscillations display variations in $V$ that are two-sided (i.e., $V$ alternately increases and decreases depending on the vector $Δ\bf{B}=\bf{B} - \bf{B}_0$). As a consequence, SBs may not appear always as one-sided spikes in $V$, especially at larger heliocentric distances where the local background field statistically departs from the radial direction. We suggest that SBs can be well described by large-amplitude Alfvénic fluctuations if the field rotation is computed with respect to a well-determined local background field that, in some cases, may deviate from the large-scale Parker field.

preprint2021arXiv

A solar source of Alfvénic magnetic field switchbacks: {\em in situ} remnants of magnetic funnels on supergranulation scales

One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfvénic magnetic field reversals termed 'switchbacks'. These $δB_R/B \sim \mathcal{O}(1$) fluctuations occur on a range of timescales and in {\em patches} separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate that patches of switchbacks are localized within the extensions of plasma structures originating at the base of the corona. These structures are characterized by an increase in alpha particle abundance, Mach number, plasma $β$ and pressure, and by depletions in the magnetic field magnitude and electron temperature. These intervals are in pressure-balance, implying stationary spatial structure, and the field depressions are consistent with overexpanded flux tubes. The structures are asymmetric in Carrington longitude with a steeper leading edge and a small ($\sim$1$^\circ$) edge of hotter plasma and enhanced magnetic field fluctuations. Some structures contain suprathermal ions to $\sim$85 keV that we argue are the energetic tail of the solar wind alpha population. The structures are separated in longitude by angular scales associated with supergranulation. This suggests that these switchbacks originate near the leading edge of the diverging magnetic field funnels associated with the network magnetic field - the primary wind sources. We propose an origin of the magnetic field switchbacks, hot plasma and suprathermals, alpha particles in interchange reconnection events just above the solar transition region and our measurements represent the extended regions of a turbulent outflow exhaust.

preprint2021arXiv

Flux ropes and dynamics of the heliospheric current sheet

Context. Solar Orbiter and PSP jointly observed the solar wind for the first time in June 2020, capturing data from very different solar wind streams, calm and Alfvénic wind as well as many dynamic structures. Aims. The aim here is to understand the origin and characteristics of the highly dynamic solar wind observed by the two probes, in particular in the vicinity of the heliospheric current sheet (HCS). Methods. We analyse the plasma data obtained by PSP and Solar Orbiter in situ during the month of June 2020. We use the Alfvén-wave turbulence MHD solar wind model WindPredict-AW, and perform two 3D simulations based on ADAPT solar magnetograms for this period. Results. We show that the dynamic regions measured by both spacecraft are pervaded with flux ropes close to the HCS. These flux ropes are also present in the simulations, forming at the tip of helmet streamers, i.e. at the base of the heliospheric current sheet. The formation mechanism involves a pressure driven instability followed by a fast tearing reconnection process, consistent with the picture of Réville et al. (2020a). We further characterize the 3D spatial structure of helmet streamer born flux ropes, which seems, in the simulations, to be related to the network of quasi-separatrices.

preprint2020arXiv

Localized magnetic field structures and their boundaries in the near-Sun solar wind from Parker Solar Probe measurements

One of the discoveries made by Parker Solar Probe during first encounters with the Sun is the ubiquitous presence of relatively small-scale structures standing out as sudden deflections of the magnetic field. They were called switchbacks as some of them show up the full reversal of the radial component of the magnetic field and then return to regular conditions. Analyzing the magnetic field and plasma perturbations associated with switchbacks we identify three types of structures with slightly different characteristics: 1. Alfvenic structures, where the variations of the magnetic field components take place while the magnitude of the field remains constant; 2. Compressional, the field magnitude varies together with changes of the components; 3. Structures manifesting full reversal of the magnetic field (extremal class of Alfvenic structures). Processing of structures boundaries and plasma bulk velocity perturbations lead to the conclusion that they represent localized magnetic field tubes with enhanced parallel plasma velocity and ion beta moving together with respect to surrounding plasma. The magnetic field deflections before and after the switchbacks reveal the existence of total axial current. The electric currents are concentrated on the relatively narrow boundary layers on the surface of the tubes and determine the magnetic field perturbations inside the tube. These currents are closed on the structure surface, and typically have comparable azimuthal and the axial components. The surface of the structure may also accommodate an electromagnetic wave, that assists to particles in carrying currents. We suggest that the two types of structures we analyzed here may represent the local manifestations of the tube deformations corresponding to a saturated stage of the Firehose instability development.

preprint2020arXiv

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.

preprint2020arXiv

Parker Solar Probe In-Situ Observations of Magnetic Reconnection Exhausts During Encounter 1

Magnetic reconnection in current sheets converts magnetic energy into particle energy. The process may play an important role in the acceleration and heating of the solar wind close to the Sun. Observations from Parker Solar Probe provide a new opportunity to study this problem, as it measures the solar wind at unprecedented close distances to the Sun. During the 1st orbit, PSP encountered a large number of current sheets in the solar wind through perihelion at 35.7 solar radii. We performed a comprehensive survey of these current sheets and found evidence for 21 reconnection exhausts. These exhausts were observed in heliospheric current sheets, coronal mass ejections, and regular solar wind. However, we find that the majority of current sheets encountered around perihelion, where the magnetic field was strongest and plasma beta was lowest, were Alfvénic structures associated with bursty radial jets and these current sheets did not appear to be undergoing local reconnection. We examined conditions around current sheets to address why some current sheets reconnected, while others did not. A key difference appears to be the degree of plasma velocity shear across the current sheets: The median velocity shear for the 21 reconnection exhausts was 24% of the Alfvén velocity shear, whereas the median shear across 43 Alfvénic current sheets examined was 71% of the Alfvén velocity shear. This finding could suggest that large, albeit sub-Alfvénic, velocity shears suppress reconnection. An alternative interpretation is that the Alfvénic current sheets are isolated rotational discontinuities which do not undergo local reconnection.

preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 1 -- CGL description and collisionless fluid hierarchy

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the CGL fluid model in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius (FLR) corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the 4th-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the "normal" closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the "normal" closure, and beyond the 4th-order moment, Landau fluid closures are required.

preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 2 -- Kinetic theory, Padé approximants and Landau fluid closures

In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(ζ)$, and the associated plasma response function $R(ζ)=1+ζZ(ζ)=-Z'(ζ)/2$. We then consider a 1D (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the 4th-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1D closures at higher-order moments than the 4th-order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the 4th-order moment level. In the Appendix A, we provide comprehensive tables with Padé approximants of $R(ζ)$ up to the 8th-pole order, with many given in an analytic form.

preprint2019arXiv

Clustering of Intermittent Magnetic and Flow Structures near Parker Solar Probe's First Perihelion -- A Partial-Variance-of-Increments Analysis

During the Parker Solar Probe's (PSP) first perihelion pass, the spacecraft reached within a heliocentric distance of \(\sim 37~R_\odot\) and observed numerous magnetic and flow structures characterized by sharp gradients. To better understand these intermittent structures in the young solar wind, an important property to examine is their degree of correlation in time and space. To this end, we use the well-tested Partial Variance of Increments (PVI) technique to identify intermittent events in FIELDS and SWEAP observations of magnetic and proton-velocity fields (respectively) during PSP's first solar encounter, when the spacecraft was within 0.25 au from the Sun. We then examine distributions of waiting times between events with varying separation and PVI thresholds. We find power-law distributions for waiting times shorter than a characteristic scale comparable to the correlation time, suggesting a high degree of correlation that may originate in a clustering process. Waiting times longer than this characteristic time are better described by an exponential, suggesting a random memory-less Poisson process at play. These findings are consistent with near-Earth observations of solar wind turbulence. The present study complements the one by Dudok de Wit et al. (2020, present volume), which focuses on waiting times between observed "switchbacks" in the radial magnetic field.

preprint2019arXiv

Enhanced Energy Transfer Rate in Solar Wind Turbulence Observed near the Sun from Parker Solar Probe

Direct evidence of an inertial-range turbulent energy cascade has been provided by spacecraft observations in heliospheric plasmas. In the solar wind, the average value of the derived heating rate near 1 au is $\sim 10^{3}\, \mathrm{J\,kg^{-1}\,s^{-1}}$, an amount sufficient to account for observed departures from adiabatic expansion. Parker Solar Probe (PSP), even during its first solar encounter, offers the first opportunity to compute, in a similar fashion, a fluid-scale energy decay rate, much closer to the solar corona than any prior in-situ observations. Using the Politano-Pouquet third-order law and the von Kármán decay law, we estimate the fluid-range energy transfer rate in the inner heliosphere, at heliocentric distance $R$ ranging from $54\,R_{\odot}$ (0.25 au) to $36\,R_{\odot}$ (0.17 au). The energy transfer rate obtained near the first perihelion is about 100 times higher than the average value at 1 au. This dramatic increase in the heating rate is unprecedented in previous solar wind observations, including those from Helios, and the values are close to those obtained in the shocked plasma inside the terrestrial magnetosheath.

preprint2019arXiv

Observations of Energetic-Particle Population Enhancements along Intermittent Structures near the Sun from Parker Solar Probe

Observations at 1 au have confirmed that enhancements in measured energetic particle fluxes are statistically associated with "rough" magnetic fields, i.e., fields having atypically large spatial derivatives or increments, as measured by the Partial Variance of Increments (PVI) method. One way to interpret this observation is as an association of the energetic particles with trapping or channeling within magnetic flux tubes, possibly near their boundaries. However, it remains unclear whether this association is a transport or local effect; i.e., the particles might have been energized at a distant location, perhaps by shocks or reconnection, or they might experience local energization or re-acceleration. The Parker Solar Probe (PSP), even in its first two orbits, offers a unique opportunity to study this statistical correlation closer to the corona. As a first step, we analyze the separate correlation properties of the energetic particles measured by the \isois instruments during the first solar encounter. The distribution of time intervals between a specific type of event, i.e., the waiting time, can indicate the nature of the underlying process. We find that the \isois observations show a power-law distribution of waiting times, indicating a correlated (non-Poisson) distribution. Analysis of low-energy \isois data suggests that the results are consistent with the 1 au studies, although we find hints of some unexpected behavior. A more complete understanding of these statistical distributions will provide valuable insights into the origin and propagation of solar energetic particles, a picture that should become clear with future PSP orbits.

preprint2019arXiv

Observations of heating along intermittent structures in the inner heliosphere from PSP data

The solar wind proton temperature at 1-au has been found to be correlated with small-scale intermittent magnetic structures, i.e., regions with enhanced temperature are associated with coherent structures such as current sheets. Using Parker Solar Probe data from the first encounter, we study this association using measurements of radial proton temperature, employing the Partial Variance of Increments (PVI) technique to identify intermittent magnetic structures. We observe that the probability density functions of high-PVI events have higher median temperatures than those with lower PVI, The regions in space where PVI peaks were also locations that had enhanced temperatures when compared with similar regions suggesting a heating mechanism in the young solar wind that is associated with intermittency developed by a nonlinear turbulent cascade.n the immediate vicinity.

preprint2016arXiv

The "ideal" tearing mode: theory and resistive MHD simulations

Classical MHD reconnection theories, both the stationary Sweet-Parker model and the tearing instability, are known to provide rates which are too slow to explain the observations. However, a recent analysis has shown that there exists a critical threshold on current sheet's thickness, namely a/L~S^(-1/3), beyond which the tearing modes evolve on fast macroscopic Alfvenic timescales, provided the Lunquist number S is high enough, as invariably found in solar and astrophysical plasmas. Therefore, the classical Sweet-Parker scenario, for which the diffusive region scales as a/L~S^(-1/2) and thus can be up to ~100 times thinner than the critical value, is likely to be never realized in nature, as the current sheet itself disrupts in the elongation process. We present here two-dimensional, compressible, resistive MHD simulations, with S ranging from 10^5 to 10^7, that fully confirm the linear analysis. Moreover, we show that a secondary plasmoid instability always occurs when the same critical scaling is reached on the local, smaller scale, leading to a cascading explosive process, reminiscent of the flaring activity.

preprint2015arXiv

Ion kinetic energy conservation and magnetic field strength constancy in multi-fluid solar wind Alfvénic turbulence

We investigate properties of the plasma fluid motion in the large amplitude low frequency fluctuations of highly Alfvénic fast solar wind. We show that protons locally conserve total kinetic energy when observed from an effective frame of reference comoving with the fluctuations. For typical properties of the fast wind, this frame can be reasonably identified by alpha particles, which, owing to their drift with respect to protons at about the Alfvén speed along the magnetic field, do not partake in the fluid low frequency fluctuations. Using their velocity to transform proton velocity into the frame of Alfvénic turbulence, we demonstrate that the resulting plasma motion is characterized by a constant absolute value of the velocity, zero electric fields, and aligned velocity and magnetic field vectors as expected for unidirectional Alfvénic fluctuations in equilibrium. We propose that this constraint, via the correlation between velocity and magnetic field in Alfvénic turbulence, is at the origin of the observed constancy of the magnetic field: while the constant velocity corresponding to constant energy can be only observed in the frame of the fluctuations, the correspondingly constant total magnetic field, invariant for Galilean transformations, remains the observational signature, in the spacecraft frame, of the constant total energy in the Alfvén turbulence frame.

preprint2015arXiv

Observational Signatures of Coronal Loop Heating and Cooling Driven by Footpoint Shuffling

The evolution of a coronal loop is studied by means of numerical simulations of the fully compressible three-dimensional magnetohydrodynamic equations using the HYPERION code. The footpoints of the loop magnetic field are advected by random motions. As a consequence the magnetic field in the loop is energized and develops turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets: energy is deposited at small scales where heating occurs. Dissipation is non-uniformly distributed so that only a fraction of the coronal mass and volume gets heated at any time. Temperature and density are highly structured at scales which, in the solar corona, remain observationally unresolved: the plasma of our simulated loop is multi-thermal, where highly dynamical hotter and cooler plasma strands are scattered throughout the loop at sub-observational scales. Numerical simulations of coronal loops of 50000 km length and axial magnetic field intensities ranging from 0.01 to 0.04 Tesla are presented. To connect these simulations to observations we use the computed number densities and temperatures to synthesize the intensities expected in emission lines typically observed with the Extreme ultraviolet Imaging Spectrometer (EIS) on Hinode. These intensities are used to compute differential emission measure distributions using the Monte Carlo Markov Chain code, which are very similar to those derived from observations of solar active regions. We conclude that coronal heating is found to be strongly intermittent in space and time, with only small portions of the coronal loop being heated: in fact, at any given time, most of the corona is cooling down.

preprint2013arXiv

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

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.

preprint2013arXiv

Field lines twisting in a noisy corona: implications for energy storage and release, and initiation of solar eruptions

We present simulations modeling closed regions of the solar corona threaded by a strong magnetic field where localized photospheric vortical motions twist the coronal field lines. The linear and nonlinear dynamics are investigated in the reduced magnetohydrodynamic regime in Cartesian geometry. Initially the magnetic field lines get twisted and the system becomes unstable to the internal kink mode, confirming and extending previous results. As typical in this kind of investigations, where initial conditions implement smooth fields and flux-tubes, we have neglected fluctuations and the fields are laminar until the instability sets in. But previous investigations indicate that fluctuations, excited by photospheric motions and coronal dynamics, are naturally present at all scales in the coronal fields. Thus, in order to understand the effect of a photospheric vortex on a more realistic corona, we continue the simulations after kink instability sets in, when turbulent fluctuations have already developed in the corona. In the nonlinear stage the system never returns to the simple initial state with ordered twisted field lines, and kink instability does not occur again. Nevertheless field lines get twisted, but in a disordered way, and energy accumulates at large scales through an inverse cascade. This energy can subsequently be released in micro-flares or larger flares, when interaction with neighboring structures occurs or via other mechanisms. The impact on coronal dynamics and CMEs initiation is discussed.

preprint2012arXiv

Interchange reconnection in a turbulent Corona

Magnetic reconnection at the interface between coronal holes and loops, so-called interchange reconnection, can release the hotter, denser plasma from magnetically confined regions into the heliosphere, contributing to the formation of the highly variable slow solar wind. The interchange process is often thought to develop at the apex of streamers or pseudo-streamers, near Y and X-type neutral points, but slow streams with loop composition have been recently observed along fanlike open field lines adjacent to closed regions, far from the apex. However, coronal heating models, with magnetic field lines shuffled by convective motions, show that reconnection can occur continuously in unipolar magnetic field regions with no neutral points: photospheric motions induce a magnetohydrodynamic turbulent cascade in the coronal field that creates the necessary small scales, where a sheared magnetic field component orthogonal to the strong axial field is created locally and can reconnect. We propose that a similar mechanism operates near and around boundaries between open and closed regions inducing a continual stochastic rearrangement of connectivity. We examine a reduced magnetohydrodynamic model of a simplified interface region between open and closed corona threaded by a strong unipolar magnetic field. This boundary is not stationary, becomes fractal, and field lines change connectivity continuously, becoming alternatively open and closed. This model suggests that slow wind may originate everywhere along loop-coronal hole boundary regions, and can account naturally and simply for outflows at and adjacent to such boundaries and for the observed diffusion of slow wind around the heliospheric current sheet.

preprint2012arXiv

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

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

preprint2012arXiv

Reflection and dissipation of Alfvén waves in interstellar clouds

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.

preprint2012arXiv

Turbulent Coronal Heating Mechanisms: Coupling of Dynamics and Thermodynamics

Context. Photospheric motions shuffle the footpoints of the strong axial magnetic field that threads coronal loops giving rise to turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets where energy is deposited at small-scales and the heating occurs. Previous studies show that current sheets thickness is orders of magnitude smaller than current state of the art observational resolution (~700 km). Aim. In order to understand coronal heating and interpret correctly observations it is crucial to study the thermodynamics of such a system where energy is deposited at unresolved small-scales. Methods. Fully compressible three-dimensional magnetohydrodynamic simulations are carried out to understand the thermodynamics of coronal heating in the magnetically confined solar corona. Results. We show that temperature is highly structured at scales below observational resolution and nonhomogeneously distributed so that only a fraction of the coronal mass and volume gets heated at each time. Conclusions. This is a multi-thermal system where hotter and cooler plasma strands are found one next to the other also at sub-resolution scales and exhibit a temporal dynamics.

preprint2011arXiv

A model for magnetically coupled sympathetic eruptions

Sympathetic eruptions on the Sun have been observed for several decades, but the mechanisms by which one eruption can trigger another one remain poorly understood. We present a 3D MHD simulation that suggests two possible magnetic trigger mechanisms for sympathetic eruptions. We consider a configuration that contains two coronal flux ropes located within a pseudo-streamer and one rope located next to it. A sequence of eruptions is initiated by triggering the eruption of the flux rope next to the streamer. The expansion of the rope leads to two consecutive reconnection events, each of which triggers the eruption of a flux rope by removing a sufficient amount of overlying flux. The simulation qualitatively reproduces important aspects of the global sympathetic event on 2010 August 1 and provides a scenario for so-called twin filament eruptions. The suggested mechanisms are applicable also for sympathetic eruptions occurring in other magnetic configurations.

preprint2011arXiv

Magnetohydrodynamic Turbulent Cascade of Coronal Loop Magnetic Fields

The Parker model for coronal heating is investigated through a high resolution simulation. An inertial range is resolved where fluctuating magnetic energy E_M (k_perp) \propto k_\perp^{-2.7} exceeds kinetic energy E_K (k_\perp) \propto k_\perp^{-0.6}. Increments scale as δb_\ell \simeq \ell^{-0.85} and δu_\ell \simeq \ell^{+0.2} with velocity increasing at small scales, indicating that magnetic reconnection plays a prime role in this turbulent system. We show that spectral energy transport is akin to standard magnetohydrodynamic (MHD) turbulence even for a system of reconnecting current sheets sustained by the boundary. In this new MHD turbulent cascade, kinetic energy flows are negligible while cross-field flows are enhanced, and through a series of "reflections" between the two fields, cascade more than half of the total spectral energy flow.

preprint2010arXiv

Diamagnetic and Expansion Effects on the Observable Properties of the Slow Solar Wind in a Coronal Streamer

The plasma density enhancements recently observed by the Large-Angle Spectrometric Coronagraph (LASCO) instrument onboard the Solar and Heliospheric Observatory (SOHO) spacecraft have sparked considerable interest. In our previous theoretical study of the formation and initial motion of these density enhancements it is found that beyond the helmet cusp of a coronal streamer the magnetized wake configuration is resistively unstable, that a traveling magnetic island develops at the center of the streamer, and that density enhancements occur within the magnetic islands. As the massive magnetic island travels outward, both its speed and width increase. The island passively traces the acceleration of the inner part of the wake. In the present paper a few spherical geometry effects are included, taking into account either the radial divergence of the magnetic field lines and the average expansion suffered by a parcel of plasma propagating outward, using the Expanding Box Model (EBM), and the diamagnetic force due to the overall magnetic field radial gradients, the so-called melon-seed force. It is found that the values of the acceleration and density contrasts can be in good agreement with LASCO observations, provided the spherical divergence of the magnetic lines starts beyond a critical distance from the Sun and the initial stage of the formation and acceleration of the plasmoid is due to the cartesian evolution of MHD instabilities. This result provides a constraint on the topology of the magnetic field in the coronal streamer.

preprint2010arXiv

Heating of coronal loops: weak MHD turbulence and scaling laws

To understand the nonlinear dynamics of the Parker scenario for coronal heating, long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry are carried out. A loop is modeled as a box extended along the direction of the strong magnetic field $B_0$ in which the system is embedded. At the top and bottom plates, which represent the photosphere, velocity fields mimicking photospheric motions are imposed. We show that the nonlinear dynamics is described by different regimes of MHD anisotropic turbulence, with spectra characterized by intertial range power laws whose indexes range from Kolmogorov-like values ($\sim 5/3$) up to $\sim 3$. We briefly describe the bearing for coronal heating rates.

preprint2010arXiv

Shear Photospheric Forcing and the Origin of Turbulence in Coronal Loops

We present a series of numerical simulations aimed at understanding the nature and origin of turbulence in coronal loops in the framework of the Parker model for coronal heating. A coronal loop is studied via reduced magnetohydrodynamics simulations in Cartesian geometry. A uniform and strong magnetic field threads the volume between the two photospheric planes, where a velocity field in the form of a 1D shear flow pattern is present. Initially the magnetic field which developes in the coronal loop is a simple map of the photospheric velocity field. This initial configuration is unstable to a multiple tearing instability which develops islands with X and O points in the plane orthogonal to the axial field. Once the nonlinear stage sets in the system evolution is characterized by a regime of MHD turbulence dominated by magnetic energy. A well developed power law in energy spectra is observed and the magnetic field never returns to the simple initial state mapping the photospheric flow. The formation of X and O points in the planes orthogonal to the axial field allows the continued and repeated formation and dissipation of small scale current sheets where the plasma is heated. We conclude that the observed turbulent dynamics are not induced by the complexity of the pattern that the magnetic field-lines footpoints follow but they rather stem from the inherent nonlinear nature of the system.

preprint2010arXiv

The Parker Scenario for Coronal Heating as an MHD Turbulence Problem

The Parker or field line tangling model of coronal heating is investigated through long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry within the framework of reduced magnetohydrodynamics (RMHD). Slow photospheric motions induce a Poynting flux which saturates by driving an anisotropic turbulent cascade dominated by magnetic energy and characterized by current sheets elongated along the axial magnetic field. Increasing the value of the axial magnetic field different regimes of MHD turbulence develop with a bearing on coronal heating rates. In physical space magnetic field lines at the scale of convection cells appear only slightly bended in agreement with observations of large loops of current (E)UV and X-ray imagers.

preprint2009arXiv

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

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.

preprint2009arXiv

Coupling Photosphere and Corona: Linear and Turbulent Regimes

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.

preprint2009arXiv

Turbulence, Energy Transfers and Reconnection in Compressible Coronal Heating Field-line Tangling Models

MHD turbulence has long been proposed as a mechanism for the heating of coronal loops in the framework of the Parker scenario for coronal heating. So far most of the studies have focused on its dynamical properties without considering its thermodynamical and radiative features, because of the very demanding computational requirements. In this paper we extend this previous research to the compressible regime, including an energy equation, by using HYPERION, a new parallelized, viscoresistive, three-dimensional compressible MHD code. HYPERION employs a Fourier collocation -- finite difference spatial discretization, and uses a third-order Runge-Kutta temporal discretization. We show that the implementation of a thermal conduction parallel to the DC magnetic field induces a radiative emission concentrated at the boundaries, with properties similar to the chromosphere--transition region--corona system.

preprint2007arXiv

Alfven Waves and Turbulence in the Solar Atmosphere and Solar Wind

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.

preprint2007arXiv

Shell-models of RMHD turbulence and the heating of solar coronal loops

A simplified non-linear numerical model for the development of incompressible magnetohydrodynamics (MHD) in the presence of a strong magnetic field B0 and stratification, nicknamed Shell-Atm, is presented. In planes orthogonal to the mean field, the non-linear incompressible dynamics is replaced by 2D shell-models for the complex variables u and b, allowing one to reach large Reynolds numbers while at the same time carrying out sufficiently long time integrations to obtain a good statistics at moderate computational cost. The shell-models of different planes are coupled by Alfven waves propagating along B0. The model may be applied to open or closed magnetic field configurations where the axial field dominates and the plasma pressure is low; here we apply it to the specific case of a magnetic loop of the solar corona heated via turbulence driven by photospheric motions, and we use statistics for its analysis. The Alfven waves interact non-linearly and form turbulent spectra in the directions perpendicular and, via propagation, also parallel to the mean field. A heating function is obtained, and is shown to be intermittent; the average heating is consistent with values required for sustaining a hot corona, and is proportional to the aspect ratio of the loop to the power -1.5; characteristic properties of heating events are distributed as power-laws. Cross-correlations show a delay of dissipation compared to energy content.