-
Solar wind data analysis aided by synthetic modeling: a better understanding of plasma-frame variations from temporal data
Authors:
Norbert Magyar,
Jaye Verniero,
Adam Szabo,
Jiyuan Zhang,
Tom Van Doorsselaere
Abstract:
In-situ measurements of the solar wind, a turbulent and anisotropic plasma flow originating at the Sun, are mostly carried out by single spacecraft, resulting in one-dimensional time series. The conversion of these measurements to the spatial frame of the plasma is a great challenge, but required for direct comparison of the measurements with MHD turbulence theories. Here we present a toolkit, bas…
▽ More
In-situ measurements of the solar wind, a turbulent and anisotropic plasma flow originating at the Sun, are mostly carried out by single spacecraft, resulting in one-dimensional time series. The conversion of these measurements to the spatial frame of the plasma is a great challenge, but required for direct comparison of the measurements with MHD turbulence theories. Here we present a toolkit, based on the synthetic modeling of solar wind fluctuations as two-dimensional noise maps with adjustable spectral and power anisotropy, that can help with the temporal-spatial conversion of real data. Specifically, by following the spacecraft trajectory through a noise map (relative velocity and angle relative to some mean magnetic field) with properties tuned to mimic those of the solar wind, the likelihood that the temporal data fluctuations represent parallel or perpendicular fluctuations in the plasma frame can be quantified by correlating structure functions of the noise map. Synthetic temporal data can also be generated, which can provide a testing ground for analysis applied to the solar wind data. We demonstrate this tool by investigating Parker Solar Probe's E7 trajectory and data, and showcase several possible ways in which it can be used. We find that whether temporal variations in the spacecraft frame come from parallel or perpendicular variations in the plasma frame strongly depends on the spectral and power anisotropy of the measured wind. Data analysis assisted by such underlying synthetic models as presented here could open up new ways to interpret measurements in the future, specifically in the more reliable determination of plasma frame quantities from temporal measurements.
△ Less
Submitted 21 May, 2024;
originally announced May 2024.
-
Validity of using Elsässer variables to study the interaction of compressible solar wind fluctuations with a coronal mass ejection
Authors:
Chaitanya Prasad Sishtla,
Jens Pomoell,
Norbert Magyar,
Emilia Kilpua,
Simon Good
Abstract:
Alfvénic fluctuations, as modelled by the non-linear interactions of Alfvén waves of various scales, are seen to dominate solar wind turbulence. However, there is also a non-negligible component of non-Alfvénic fluctuations. The Elsässer formalism, which is central to the study of Alfvénic turbulence due to its ability to differentiate between parallel and anti-parallel Alfvén waves, cannot strict…
▽ More
Alfvénic fluctuations, as modelled by the non-linear interactions of Alfvén waves of various scales, are seen to dominate solar wind turbulence. However, there is also a non-negligible component of non-Alfvénic fluctuations. The Elsässer formalism, which is central to the study of Alfvénic turbulence due to its ability to differentiate between parallel and anti-parallel Alfvén waves, cannot strictly separate wavemodes in the presence of compressive magnetoacoustic waves. In this study, we analyse the deviations generated in the Elsässer formalism as density fluctuations are naturally generated through the propagation of a linearly polarised Alfvén wave. The study was performed in the context of a coronal mass ejection (CME) propagating through the solar wind, which enables the creation of two solar wind regimes, pristine wind and a shocked CME sheath, where the Elsässer formalism can be evaluated. In these two regimes we studied the deviations of the Elsässer formalism in separating parallel and anti-parallel components of Alfvénic solar wind perturbations generated by small-amplitude density fluctuations. We used an ideal 2.5D magnetohydrodynamic (MHD) model with an adiabatic equation of state. An Alfvén pump wave was injected into the quiet solar wind by perturbing the transverse magnetic field and velocity components. This wave subsequently generates density fluctuations through the ponderomotive force. A CME was injected by inserting a flux-rope modelled as a magnetic island into the quasi-steady solar wind. The presence of density perturbations creates an approximately 10% deviation in the Elsässer variables and reflection coefficient for the Alfvén waves as well as a deviation of approximately 0.1 in the cross helicity in regions containing both parallel and anti-parallel fluctuations.
△ Less
Submitted 14 February, 2024;
originally announced February 2024.
-
The magnetohydrodynamic equations in terms of waveframe variables
Authors:
Tom Van Doorsselaere,
Norbert Magyar,
M. Valeria Sieyra,
Marcel Goossens
Abstract:
Generalising the Elsasser variables, we introduce the Q-variables. These are more flexible than the Elsasser variables, because they also allow to track waves with phase speeds different than the Alfven speed. We rewrite the MHD equations with these Q-variables. We consider also the linearised version of the resulting MHD equations in a uniform plasma, and recover the classical Alfven waves, but a…
▽ More
Generalising the Elsasser variables, we introduce the Q-variables. These are more flexible than the Elsasser variables, because they also allow to track waves with phase speeds different than the Alfven speed. We rewrite the MHD equations with these Q-variables. We consider also the linearised version of the resulting MHD equations in a uniform plasma, and recover the classical Alfven waves, but also separate the fast and slow magnetosonic waves in upward and downward propagating waves. Moreover, we show that the Q-variables may also track the upward and downward propagating surface Alfven waves in a non-uniform plasma, displaying the power of our generalisation. In the end, we lay the mathematical framework for driving solar wind models with a multitude of wave drivers.
△ Less
Submitted 15 January, 2024;
originally announced January 2024.
-
The Effect of the Parametric Decay Instability on the Morphology of Coronal Type III Radio Bursts
Authors:
Chaitanya Prasad Sishtla,
Immanuel Christopher Jebaraj,
Jens Pomoell,
Norbert Magyar,
Marc Pulupa,
Emilia Kilpua,
Stuart D. Bale
Abstract:
The nonlinear evolution of Alfvén waves in the solar corona leads to the generation of Alfvénic turbulence. This description of the Alfvén waves involves parametric instabilities where the parent wave decays into slow mode waves giving rise to density fluctuations. These density fluctuations, in turn, play a crucial role in the modulation of the dynamic spectrum of type III radio bursts, which are…
▽ More
The nonlinear evolution of Alfvén waves in the solar corona leads to the generation of Alfvénic turbulence. This description of the Alfvén waves involves parametric instabilities where the parent wave decays into slow mode waves giving rise to density fluctuations. These density fluctuations, in turn, play a crucial role in the modulation of the dynamic spectrum of type III radio bursts, which are observed at the fundamental of local plasma frequency and are sensitive to the local density. During observations of such radio bursts, fine structures are detected across different temporal ranges. In this study, we examine density fluctuations generated through the parametric decay instability (PDI) of Alfvén waves as a mechanism to generate striations in the dynamic spectrum of type III radio bursts using magnetohydrodynamic simulations of the solar corona. An Alfvén wave is injected into the quiet solar wind by perturbing the transverse magnetic field and velocity components which subsequently undergo the PDI instability. The type III burst is modelled as a fast-moving radiation source that samples the background solar wind as it propagates to emit radio waves. We find the simulated dynamic spectrum to contain striations directly affected by the multi-scale density fluctuations in the wind.
△ Less
Submitted 16 December, 2023;
originally announced December 2023.
-
On the Role of Alfvenic Fluctuations as Mediators of Coherence within Interplanetary Coronal Mass Ejections: Investigation of Multi-Spacecraft Measurements at 1 au
Authors:
Camilla Scolini,
Noe Lugaz,
Reka M. Winslow,
Charles J. Farrugia,
Norbert Magyar,
Fabio Bacchini
Abstract:
Interplanetary coronal mass ejections (ICMEs) are defined as ''coherent'' if they are capable of responding to external perturbations in a collective manner. This implies that information must be able to propagate across ICME structures, and if this is not the case, single-point in-situ measurements cannot be considered as indicative of global ICME properties. Here, we investigate the role of Alfv…
▽ More
Interplanetary coronal mass ejections (ICMEs) are defined as ''coherent'' if they are capable of responding to external perturbations in a collective manner. This implies that information must be able to propagate across ICME structures, and if this is not the case, single-point in-situ measurements cannot be considered as indicative of global ICME properties. Here, we investigate the role of Alfvenic fluctuations (AFs) as mediators of ICME coherence. We consider multi-point magnetic field and plasma measurements of 10 ICMEs observed by the ACE and Wind spacecraft at 1 au at longitudinal separations of 0.5°-0.7°. For each event, we analyze the Alfvenicity in terms of the residual energy and cross helicity of fluctuations, and the coherence in terms of the magnetic correlation between Wind and ACE. We find that ~65% and 90% of ICME sheaths and magnetic ejecta (MEs), respectively, present extended AFs covering at least 20% of the structure. Cross helicity suggests AFs of solar and interplanetary origin may co-exist in the ICME population at 1 au. AFs are mainly concentrated downstream of shocks and in the back of MEs. The magnetic field is poorly correlated within sheaths, while the correlation decreases from the front to the back of the MEs for most magnetic field components. AFs are also associated with lower magnetic field correlations. This suggests either that ICME coherence is not mediated by Alfven waves, implying that the coherence scale may be smaller than previously predicted, or that the magnetic field correlation is not a measure of coherence.
△ Less
Submitted 7 December, 2023;
originally announced December 2023.
-
Phase mixing and the 1/f spectrum in the solar wind
Authors:
Norbert Magyar,
Tom Van Doorsselaere
Abstract:
The origin and evolution of the 1/f power law observed in the energy spectrum of solar coronal and solar wind fluctuations at scales of around an hour is not entirely understood. Several existing theories aim at explaining it, involving both linear and nonlinear mechanisms. An often overlooked property of the solar corona and solar wind is their highly inhomogeneous nature. In this paper we invest…
▽ More
The origin and evolution of the 1/f power law observed in the energy spectrum of solar coronal and solar wind fluctuations at scales of around an hour is not entirely understood. Several existing theories aim at explaining it, involving both linear and nonlinear mechanisms. An often overlooked property of the solar corona and solar wind is their highly inhomogeneous nature. In this paper we investigate the linear evolution of pure Alfvén and surface Alfvén waves propagating through a plasma which is inhomogeneous across the magnetic field. The inhomogeneity is given by density, which we model to be two-dimensional colored noise, with power spectral slopes ranging from -2 to -1. Alfvén waves propagate independently on individual magnetic field lines, and eventually get completely out of phase through the process of phase mixing, leading to unrealistic spectra. When the coupling between the inhomogeneous background and the propagating waves is fully accounted for, transverse waves such as surface Alfvén waves (also referred to as kink or Alfvénic) appear, showing collective wave behavior of neighboring magnetic field lines with different Alfvén speeds. We show that the linear cascade of surface Alfvén wave energy, induced by phase mixing and resonant absorption, leads to a perpendicular wave energy spectrum which tends to the perpendicular power spectrum of the background density. Based on our model, we propose that a perpendicular density power spectrum of 1/f in the solar corona can induce, through linear processes, the 1/f spectrum of the fluctuations that is observed at the largest scales.
△ Less
Submitted 18 August, 2022;
originally announced August 2022.
-
High frequency decayless waves with significant energy in Solar Orbiter/EUI observations
Authors:
Elena Petrova,
Norbert Magyar,
Tom Van Doorsselaere,
David Berghmans
Abstract:
High-frequency wave phenomena present a great deal of interest as one of the possible candidates to contribute to the energy input required to heat the corona as a part of the AC heating theory. However, the resolution of imaging instruments up until the Solar Orbiter have made it impossible to resolve the necessary time and spatial scales. The present paper reports on high-frequency transverse mo…
▽ More
High-frequency wave phenomena present a great deal of interest as one of the possible candidates to contribute to the energy input required to heat the corona as a part of the AC heating theory. However, the resolution of imaging instruments up until the Solar Orbiter have made it impossible to resolve the necessary time and spatial scales. The present paper reports on high-frequency transverse motions in a small loop located in a quiet Sun region of the corona. The oscillations were observed with the HRIEUV telescope (17.4 nm) of the EUI instrument onboard the Solar Orbiter. We detect two transverse oscillations in short loops with lengths of 4.5 Mm and 11 Mm. The shorter loop displays an oscillation with a 14 s period and the longer a 30 s period. Despite the high resolution, no definitive identification as propagating or standing waves is possible. The velocity amplitudes are found to be equal to 72 km/s and 125 km/s, respectively, for the shorter and longer loop. Based on that, we also estimated the values of the energy flux contained in the loops - the energy flux of the 14 s oscillation is 1.9 kW m^-2 and of the 30 s oscillation it is 6.5 kW m^-2 . While these oscillations have been observed in the Quiet Sun, their energy fluxes are of the same order as the energy input required to heat the active solar corona. Numerical simulations were performed in order to reproduce the observed oscillations. The correspondence of the numerical results to the observations provide support to the energy content estimates for the observations. Such high energy densities have not yet been observed in decayless coronal waves, and this is promising for coronal heating models based on wave damping.
△ Less
Submitted 11 May, 2022;
originally announced May 2022.
-
The dynamics and observability of circularly polarized kink waves
Authors:
N. Magyar,
T. Duckenfield,
T. Van Doorsselaere,
V. M. Nakariakov
Abstract:
Context. Kink waves are routinely observed in coronal loops. Resonant absorption is a well-accepted mechanism that extracts energy from kink waves. Nonlinear kink waves are know to be affected by the Kelvin-Helmholtz instability. However, all previous numerical studies consider linearly polarized kink waves. Aims. We study the properties of circularly polarized kink waves on straight plasma cylind…
▽ More
Context. Kink waves are routinely observed in coronal loops. Resonant absorption is a well-accepted mechanism that extracts energy from kink waves. Nonlinear kink waves are know to be affected by the Kelvin-Helmholtz instability. However, all previous numerical studies consider linearly polarized kink waves. Aims. We study the properties of circularly polarized kink waves on straight plasma cylinders, for both standing and propagating waves, and compare them to the properties of linearly polarized kink waves. Methods. We use the code MPI-AMRVAC to solve the full 3D Magnetohydrodynamic (MHD) equations for a straight magnetic cylinder, excited by both standing and propagating circularly polarized kink (m = 1) modes. Results. The damping due to resonant absorption is independent of the polarization state. The morphology or appearance of the induced resonant flow is different for the two polarizations, however, there are essentially no differences in the forward-modeled Doppler signals. For nonlinear oscillations, the growth rate of small scales is determined by the total energy of the oscillation rather than the perturbation amplitude. We discuss possible implications and seismological relevance.
△ Less
Submitted 30 December, 2021;
originally announced December 2021.
-
Kink oscillations of coronal loops
Authors:
V. M. Nakariakov,
S. A. Anfinogentov,
P. Antolin,
R. Jain,
D. Y. Kolotkov,
E. G. Kupriyanova,
D. Li,
N. Magyar,
G. Nistico,
D. J. Pascoe,
A. K. Srivastava,
J. Terradas,
S. Vashegani Farahani,
G. Verth,
D. Yuan,
I. V. Zimovets
Abstract:
Kink oscillations of coronal loops, i.e., standing kink waves, is one of the most studied dynamic phenomena in the solar corona. The oscillations are excited by impulsive energy releases, such as low coronal eruptions. Typical periods of the oscillations are from a few to several minutes, and are found to increase linearly with the increase in the major radius of the oscillating loops. It clearly…
▽ More
Kink oscillations of coronal loops, i.e., standing kink waves, is one of the most studied dynamic phenomena in the solar corona. The oscillations are excited by impulsive energy releases, such as low coronal eruptions. Typical periods of the oscillations are from a few to several minutes, and are found to increase linearly with the increase in the major radius of the oscillating loops. It clearly demonstrates that kink oscillations are natural modes of the loops, and can be described as standing fast magnetoacoustic waves with the wavelength determined by the length of the loop. Kink oscillations are observed in two different regimes. In the rapidly decaying regime, the apparent displacement amplitude reaches several minor radii of the loop. The damping time which is about several oscillation periods decreases with the increase in the oscillation amplitude, suggesting a nonlinear nature of the damping. In the decayless regime, the amplitudes are smaller than a minor radius, and the driver is still debated. The review summarises major findings obtained during the last decade, and covers both observational and theoretical results. Observational results include creation and analysis of comprehensive catalogues of the oscillation events, and detection of kink oscillations with imaging and spectral instruments in the EUV and microwave bands. Theoretical results include various approaches to modelling in terms of the magnetohydrodynamic wave theory. Properties of kink oscillations are found to depend on parameters of the oscillating loop, such as the magnetic twist, stratification, steady flows, temperature variations and so on, which make kink oscillations a natural probe of these parameters by the method of magnetohydrodynamic seismology.
△ Less
Submitted 23 September, 2021;
originally announced September 2021.
-
Non-linear damping of standing kink waves computed with Elsasser variables
Authors:
Tom Van Doorsselaere,
Marcel Goossens,
Norbert Magyar,
Michael S. Ruderman,
Rajab Ismayilli
Abstract:
In a previous paper, we computed the energy density and the non-linear energy cascade rate for transverse kink waves using Elsasser variables. In this paper, we focus on the standing kink waves, which are impulsively excited in coronal loops by external perturbations. We present an analytical calculation to compute the damping time due to the non-linear development of the Kelvin-Helmholtz instabil…
▽ More
In a previous paper, we computed the energy density and the non-linear energy cascade rate for transverse kink waves using Elsasser variables. In this paper, we focus on the standing kink waves, which are impulsively excited in coronal loops by external perturbations. We present an analytical calculation to compute the damping time due to the non-linear development of the Kelvin-Helmholtz instability. The main result is that the damping time is inversely proportional to the oscillation amplitude. We compare the damping times from our formula with the results of numerical simulations and observations. In both cases we find a reasonably good match. The comparison with the simulations show that the non-linear damping dominates in the high amplitude regime, while the low amplitude regime shows damping by resonant absorption. In the comparison with the observations, we find a power law inversely proportional to the amplitude $η^{-1}$ as an outer envelope for our Monte Carlo data points.
△ Less
Submitted 29 April, 2021;
originally announced April 2021.
-
Could switchbacks originate in the lower solar atmosphere? II. Propagation of switchbacks in the solar corona
Authors:
Norbert Magyar,
Dominik Utz,
Robertus Erdélyi,
Valery M. Nakariakov
Abstract:
The magnetic switchbacks observed recently by the Parker Solar Probe have raised the question about their nature and origin. One of the competing theories of their origin is the interchange reconnection in the solar corona. In this scenario, switchbacks are generated at the reconnection site between open and closed magnetic fields, and are either advected by an upflow or propagate as waves into th…
▽ More
The magnetic switchbacks observed recently by the Parker Solar Probe have raised the question about their nature and origin. One of the competing theories of their origin is the interchange reconnection in the solar corona. In this scenario, switchbacks are generated at the reconnection site between open and closed magnetic fields, and are either advected by an upflow or propagate as waves into the solar wind. In this paper we test the wave hypothesis, numerically modelling the propagation of a switchback, modelled as an embedded Alfvén wave packet of constant magnetic field magnitude, through the gravitationally stratified solar corona with different degrees of background magnetic field expansion. While switchbacks propagating in a uniform medium with no gravity are relatively stable, as reported previously, we find that gravitational stratification together with the expansion of the magnetic field act in multiple ways to deform the switchbacks. These include WKB effects, which depend on the degree of magnetic field expansion, and also finite-amplitude effects, such as the symmetry breaking between nonlinear advection and the Lorentz force. In a straight or radially expanding magnetic field the propagating switchbacks unfold into waves that cause minimal magnetic field deflections, while a super-radially expanding magnetic field aids in maintaining strong deflections. Other important effects are the mass uplift the propagating switchbacks induce and the reconnection and drainage of plasmoids contained within the switchbacks. In the Appendix, we examine a series of setups with different switchback configurations and parameters, which broaden the scope of our study.
△ Less
Submitted 20 April, 2021;
originally announced April 2021.
-
Could switchbacks originate in the lower solar atmosphere? I. Formation mechanisms of switchbacks
Authors:
Norbert Magyar,
Dominik Utz,
Robertus Erdélyi,
Valery M. Nakariakov
Abstract:
The recent rediscovery of magnetic field switchbacks or deflections embedded in the solar wind flow by the Parker Solar Probe mission lead to a huge interest in the modelling of the formation mechanisms and origin of these switchbacks. Several scenarios for their generation were put forth, ranging from lower solar atmospheric origins by reconnection, to being a manifestation of turbulence in the s…
▽ More
The recent rediscovery of magnetic field switchbacks or deflections embedded in the solar wind flow by the Parker Solar Probe mission lead to a huge interest in the modelling of the formation mechanisms and origin of these switchbacks. Several scenarios for their generation were put forth, ranging from lower solar atmospheric origins by reconnection, to being a manifestation of turbulence in the solar wind, and so on. Here we study some potential formation mechanisms of magnetic switchbacks in the lower solar atmosphere, using three-dimensional magneto-hydrodynamic (MHD) numerical simulations. The model is that of an intense flux tube in an open magnetic field region, aiming to represent a magnetic bright point opening up to an open coronal magnetic field structure, e.g. a coronal hole. The model is driven with different plasma flows in the photosphere, such as a fast up-shooting jet, as well as shearing flows generated by vortex motions or torsional oscillations. In all scenarios considered, we witness the formation of magnetic switchbacks in regions corresponding to chromospheric heights. Therefore, photospheric plasma flows around the foot-points of intense flux tubes appear to be suitable drivers for the formation of magnetic switchbacks in the lower solar atmosphere. Nevertheless, these switchbacks do not appear to be able to enter the coronal heights of the simulation in the present model. In conclusion, based on the presented simulations, switchbacks measured in the solar wind are unlikely to originate from photospheric or chromospheric dynamics.
△ Less
Submitted 5 March, 2021;
originally announced March 2021.
-
MHD Waves in open coronal structures
Authors:
D. Banerjee,
S. Krishna Prasad,
V. Pant,
J. A. McLaughlin,
P. Antolin,
N. Magyar,
L. Ofman,
H. Tian,
T. Van Doorsselaere,
I. De Moortel,
T. Wang
Abstract:
Modern observatories have revealed the ubiquitous presence of magnetohydrodynamic waves in the solar corona. The propagating waves (in contrast to the standing waves) are usually originated in the lower solar atmosphere which makes them particularly relevant to coronal heating. Furthermore, open coronal structures are believed to be the source regions of solar wind, therefore, the detection of MHD…
▽ More
Modern observatories have revealed the ubiquitous presence of magnetohydrodynamic waves in the solar corona. The propagating waves (in contrast to the standing waves) are usually originated in the lower solar atmosphere which makes them particularly relevant to coronal heating. Furthermore, open coronal structures are believed to be the source regions of solar wind, therefore, the detection of MHD waves in these structures is also pertinent to the acceleration of solar wind. Besides, the advanced capabilities of the current generation telescopes have allowed us to extract important coronal properties through MHD seismology. The recent progress made in the detection, origin, and damping of both propagating slow mangetoacoustic waves and kink (Alfvénic) waves is presented in this review article especially in the context of open coronal structures. Where appropriate, we give an overview on associated theoretical modelling studies. A few of the important seismological applications of these waves are discussed. The possible role of Aflvénic waves in the acceleration of solar wind is also touched upon.
△ Less
Submitted 21 September, 2021; v1 submitted 16 December, 2020;
originally announced December 2020.
-
Coronal heating by MHD waves
Authors:
Tom Van Doorsselaere,
Abhishek K. Srivastava,
Patrick Antolin,
Norbert Magyar,
Soheil Vasheghani Farahani,
Hui Tian,
Dmitrii Y. Kolotkov,
Leon Ofman,
Mingzhe Guo,
Iñigo Arregui,
Ineke De Moortel,
David Pascoe
Abstract:
The heating of the solar chromosphere and corona to the observed high temperatures, imply the presence of ongoing heating that balances the strong radiative and thermal conduction losses expected in the solar atmosphere. It has been theorized for decades that the required heating mechanisms of the chromospheric and coronal parts of the active regions, quiet-Sun, and coronal holes are associated wi…
▽ More
The heating of the solar chromosphere and corona to the observed high temperatures, imply the presence of ongoing heating that balances the strong radiative and thermal conduction losses expected in the solar atmosphere. It has been theorized for decades that the required heating mechanisms of the chromospheric and coronal parts of the active regions, quiet-Sun, and coronal holes are associated with the solar magnetic fields. However, the exact physical process that transport and dissipate the magnetic energy which ultimately leads to the solar plasma heating are not yet fully understood. The current understanding of coronal heating relies on two main mechanism: reconnection and MHD waves that may have various degrees of importance in different coronal regions. In this review we focus on recent advances in our understanding of MHD wave heating mechanisms. First, we focus on giving an overview of observational results, where we show that different wave modes have been discovered in the corona in the last decade, many of which are associated with a significant energy flux, either generated in situ or pumped from the lower solar atmosphere. Afterwards, we summarise the recent findings of numerical modelling of waves, motivated by the observational results. Despite the advances, only 3D MHD models with Alfvén wave heating in an unstructured corona can explain the observed coronal temperatures compatible with the quiet Sun, while 3D MHD wave heating models including cross-field density structuring are not yet able to account for the heating of coronal loops in active regions to their observed temperature.
△ Less
Submitted 2 December, 2020;
originally announced December 2020.
-
Three-dimensional simulations of the inhomogeneous Low Solar Wind
Authors:
N. Magyar,
V. M. Nakariakov
Abstract:
In the near future, Parker Solar Probe will put theories about the dynamics and nature of the transition between the solar corona and the solar wind to stringent tests. The most popular mechanism aimed to explain the dynamics of the nascent solar wind, including its heating and acceleration is magnetohydrodynamic (MHD) turbulence. Most of the previous models focus on nonlinear cascade induced by i…
▽ More
In the near future, Parker Solar Probe will put theories about the dynamics and nature of the transition between the solar corona and the solar wind to stringent tests. The most popular mechanism aimed to explain the dynamics of the nascent solar wind, including its heating and acceleration is magnetohydrodynamic (MHD) turbulence. Most of the previous models focus on nonlinear cascade induced by interactions of outgoing Alfvén waves and their reflections, ignoring effects that might be related to perpendicular structuring of the solar coronal plasma, despite overwhelming evidence for it. In this paper, for the first time, we analyse through 3D MHD numerical simulations the dynamics of the perpendicularly structured solar corona and solar wind, from the low corona to 15 R_Sun. We find that background structuring has a strong effect on the evolution of MHD turbulence, on much faster time scales than in the perpendicularly homogeneous case. On time scales shorter than nonlinear times, linear effects related to phase mixing result in a 1/f perpendicular energy spectrum. As the turbulent cascade develops, we observe a perpendicular (parallel) energy spectrum with the power law index of -3/2 or -5/3 (-2), a steeper perpendicular magnetic field than velocity spectrum, and a strong build-up of negative residual energy. We conclude that the turbulence is most probably generated by the self-cascade of the driven transverse kink waves, referred to previously as `uniturbulence', which might represent the dominant nonlinear energy cascade channel in the pristine solar wind.
△ Less
Submitted 1 December, 2020;
originally announced December 2020.
-
Wave pressure and energy cascade rate of kink waves computed with Elsasser variables
Authors:
Tom Van Doorsselaere,
Bo Li,
Marcel Goossens,
Bogdan Hnat,
Norbert Magyar
Abstract:
Numerical simulations have revealed a new type of turbulence of unidirectional waves in a plasma that is perpendicularly structured (Magyar et al. 2017), named uniturbulence. For this new type of turbulence, the transverse structuring modifies the upward propagating wave to have both Elsasser variables, leading to the well-known perpendicular cascade. In this paper, we study an analytical descript…
▽ More
Numerical simulations have revealed a new type of turbulence of unidirectional waves in a plasma that is perpendicularly structured (Magyar et al. 2017), named uniturbulence. For this new type of turbulence, the transverse structuring modifies the upward propagating wave to have both Elsasser variables, leading to the well-known perpendicular cascade. In this paper, we study an analytical description of the non-linear evolution of kink waves in a cylindrical flux tube, which are prone to uniturbulence. We show that they lead to a non-linear cascade for both propagating and standing waves. We calculate explicit expressions for the wave pressure and energy cascade rate. The computed damping rate τ/P depends on the density contrast of the flux tube and the background plasma and is inversely proportional to the amplitude of the kink wave. The dependence on the density contrast shows that it plays a role especially in the lower solar corona. These expressions may be added in Alfven wave driven models of the solar atmosphere (such as AWSOM, van der Holst et al. 2014), modifying it to UAWSOM (Uniturbulence and Alfven Wave Solar Model).
△ Less
Submitted 30 July, 2020;
originally announced July 2020.
-
Standing kink waves in sigmoid solar coronal loops: implications for coronal seismology
Authors:
N. Magyar,
V. M. Nakariakov
Abstract:
Using full three-dimensional magnetohydrodynamic numerical simulations, we study the effects of magnetic field sigmoidity or helicity on the properties of the fundamental kink oscillation of solar coronal loops. Our model consists of a single denser coronal loop, embedded in a plasma with dipolar force-free magnetic field with a constant alpha-parameter. For the loop with no sigmoidity, we find th…
▽ More
Using full three-dimensional magnetohydrodynamic numerical simulations, we study the effects of magnetic field sigmoidity or helicity on the properties of the fundamental kink oscillation of solar coronal loops. Our model consists of a single denser coronal loop, embedded in a plasma with dipolar force-free magnetic field with a constant alpha-parameter. For the loop with no sigmoidity, we find that the numerically determined oscillation period of the fundamental kink mode matches the theoretical period calculated using WKB theory. In contrast, with increasing sigmoidity of the loop, the actual period is increasingly smaller than the one estimated by WKB theory. Translated through coronal seismology, increasing sigmoidity results in magnetic field estimates which are increasingly shifting towards higher values, and even surpassing the average value for the highest alpha value considered. Nevertheless, the estimated range of the coronal magnetic field value lies within the mimimal/maximal limits, proving the robustness coronal seismology. We propose that the discrepancy in the estimations of the absolute value of the force-free magnetic field could be exploited seismologically to determine the free energy of coronal loops, if averages of the internal magnetic field and density can be reliably estimated by other methods.
△ Less
Submitted 29 April, 2020;
originally announced April 2020.
-
Understanding uniturbulence: self-cascade of MHD waves in the presence of inhomogeneities
Authors:
N. Magyar,
T. Van Doorsselaere,
M. Goossens
Abstract:
It is widely accepted in the MHD turbulence community that the nonlinear cascade of wave energy requires counter-propagating Alfvénic wave-packets, along some mean magnetic field. This fact is an obvious outcome of the MHD equations under the assumptions of incompressibility and homogeneity. Despite attempts to relax these assumptions in the context of MHD turbulence, the central idea of turbulenc…
▽ More
It is widely accepted in the MHD turbulence community that the nonlinear cascade of wave energy requires counter-propagating Alfvénic wave-packets, along some mean magnetic field. This fact is an obvious outcome of the MHD equations under the assumptions of incompressibility and homogeneity. Despite attempts to relax these assumptions in the context of MHD turbulence, the central idea of turbulence generation persists. However, once the assumptions of incompressiblity and homogeneity break down, the generally accepted picture of turbulent cascade generation is not universal. In this paper, we show that perpendicular inhomogeneities (across the mean magnetic field) lead to propagating wave solutions which are necessarily described by co-propagating Elsässer fields, already in the incompressible case. One simple example of these wave solutions is the surface Alfvén wave on a planar discontinuity across the magnetic field. We show through numerical simulations how the nonlinear self-deformation of these unidirectionally propagating waves leads to a cascade of wave energy across the magnetic field. The existence of this type of unidirectional cascade might have an additional strong effect on the turbulent dissipation rate of dominantly outward propagating Alfvénic waves in structured plasma, as in the solar corona and solar wind.
△ Less
Submitted 24 July, 2019;
originally announced July 2019.
-
Investigating `dark' energy in the solar corona using forward modeling of MHD waves
Authors:
Vaibhav Pant,
Norbert Magyar,
Tom Van Doorsselaere,
Richard J. Morton
Abstract:
It is now well established that the Alfvénic waves are ubiquitous in the solar corona. However, the Alfvénic wave energy estimated from the Doppler velocity measurements in the corona was found to be four orders of magnitude less than that estimated from non-thermal line widths. McIntosh & De Pontieu (2012) suggested that this discrepancy in energy might be due to the line-of-sight (LOS) superposi…
▽ More
It is now well established that the Alfvénic waves are ubiquitous in the solar corona. However, the Alfvénic wave energy estimated from the Doppler velocity measurements in the corona was found to be four orders of magnitude less than that estimated from non-thermal line widths. McIntosh & De Pontieu (2012) suggested that this discrepancy in energy might be due to the line-of-sight (LOS) superposition of the several oscillating structures, which can lead to an underestimation of the Alfvénic wave amplitudes and energies. McIntosh & De Pontieu (2012) termed this coronal `dark' or `hidden' energy. However, their simulations required the use of an additional, unknown source of Alfvénic wave energy to provide agreement with measurements of the coronal non-thermal line widths. In this study, we investigate the requirement of this unknown source of additional `dark' energy in the solar corona using gravitationally stratified 3D magnetohydrodynamic (MHD) simulations of propagating waves. We excite the transverse MHD waves and generate synthetic observations for the Fe XIII emission line. We establish that the LOS superposition greatly reduces the Doppler velocity amplitudes and increases the non-thermal line widths. Importantly, our model generates the observed wedge-shaped correlation between Doppler velocities and non-thermal line widths. We find that the observed wave energy is only 0.2-1\% of the true wave energy which explains 2-3 orders of magnitude of the energy discrepancy. We conclusively establish that the true wave energies are hidden in the non-thermal line widths. Hence, our results rule out the requirement for an additional `dark' energy in the solar corona.
△ Less
Submitted 26 June, 2019;
originally announced June 2019.
-
The nature of Elsässer variables in compressible MHD
Authors:
N. Magyar,
T. Van Doorsselaere,
M. Goossens
Abstract:
The Elsässer variables are often used in studies of plasma turbulence, in helping differentiate between MHD waves propagating parallel or anti-parallel to the main magnetic field. While for pure Alfvén waves in a homogeneous plasma the method is strictly valid, we show that compressible, magnetoacoustic waves are in general described by both Elsässer variables. Furthermore, in a compressible and i…
▽ More
The Elsässer variables are often used in studies of plasma turbulence, in helping differentiate between MHD waves propagating parallel or anti-parallel to the main magnetic field. While for pure Alfvén waves in a homogeneous plasma the method is strictly valid, we show that compressible, magnetoacoustic waves are in general described by both Elsässer variables. Furthermore, in a compressible and inhomogeneous plasma, the pure MHD waves (Alfvén, fast and slow) are no longer normal modes, but waves become linearly coupled or display mixed properties of Alfvén and magnetoacoustic nature. These waves are necessarily described by both Elsässer variables and therefore the Elsässer formalism cannot be used to strictly separate parallel and anti-parallel propagating waves. Nevertheless, even in an inhomogeneous plasma, for a highly Alfvénic wave the Elsässer variable corresponding to the propagation direction appears still dominating. We suggest that for Alfvénic waves, the relative amplitude of Elsässer variables depends on the local degree of inhomogeneity and other plasma and wave properties. This finding has implications for turbulence studies in inhomogeneous and compressible plasmas, such as the solar corona and solar wind.
△ Less
Submitted 5 February, 2019;
originally announced February 2019.
-
The polytropic index of solar coronal plasma in sunspot fan loops and its temperature dependence
Authors:
S. Krishna Prasad,
J. O. Raes,
T. Van Doorsselaere,
N. Magyar,
D. B. Jess
Abstract:
Observations of slow magneto-acoustic waves have been demonstrated to possess a number of applications in coronal seismology. Determination of the polytropic index ($γ$) is one such important application. Analysing the amplitudes of oscillations in temperature and density corresponding to a slow magneto-acoustic wave, the polytropic index in the solar corona has been calculated and based on the ob…
▽ More
Observations of slow magneto-acoustic waves have been demonstrated to possess a number of applications in coronal seismology. Determination of the polytropic index ($γ$) is one such important application. Analysing the amplitudes of oscillations in temperature and density corresponding to a slow magneto-acoustic wave, the polytropic index in the solar corona has been calculated and based on the obtained value it has been inferred that thermal conduction is highly suppressed in a very hot loop in contrast to an earlier report of high thermal conduction in a relatively colder loop. In this study, using SDO/AIA data, we analysed slow magneto-acoustic waves propagating along sunspot fan loops from 30 different active regions and computed polytropic indices for several loops at multiple spatial positions. The obtained $γ$ values vary from 1.04$\pm$0.01 to 1.58$\pm$0.12 and most importantly display a temperature dependence indicating higher $γ$ at hotter temperatures. This behaviour brings both the previous studies to agreement and perhaps implies a gradual suppression of thermal conduction with increase in temperature of the loop. The observed phase shifts between temperature and density oscillations, however, are substantially larger than that expected from a classical thermal conduction and appear to be influenced by a line-of-sight integration effect on the emission measure.
△ Less
Submitted 19 October, 2018;
originally announced October 2018.
-
Assessing the Capabilities of Dynamic Coronal Seismology of Alfvénic Waves through Forward Modeling
Authors:
N. Magyar,
T. Van Doorsselaere
Abstract:
Coronal seismology is a diagnostic tool used in solar physics for measuring parameters that are otherwise hard to measure; of these parameters, magnetic field values are arguably the most important. The parameters are inferred by combining observations of waves with magnetohydrodynamic (MHD) wave theory. To date, coronal seismology has successfully been applied to various single-oscillation events…
▽ More
Coronal seismology is a diagnostic tool used in solar physics for measuring parameters that are otherwise hard to measure; of these parameters, magnetic field values are arguably the most important. The parameters are inferred by combining observations of waves with magnetohydrodynamic (MHD) wave theory. To date, coronal seismology has successfully been applied to various single-oscillation events. Such events are relatively rare, resulting in rare occasions to use diagnostics. Ubiquitous waves in the solar atmosphere might, however, allow for the possibility of dynamic coronal seismology, which involves the continuous inversions of coronal parameters and would constitute a huge leap forward in many areas of solar physics. In this paper, we investigate the robustness and accuracy of magnetic field diagnostics applied to forward-modeled 3D MHD simulations of propagating Alfvénic waves. We find that the seismologically measured magnetic field values are reassuringly close to the input value (within 20%) for a range of setups studied, providing encouragement and confidence for the further development of dynamic coronal seismology.
△ Less
Submitted 16 April, 2018; v1 submitted 6 April, 2018;
originally announced April 2018.
-
Effect of magnetic twist on nonlinear transverse kink oscillations
Authors:
J. Terradas,
N. Magyar,
T. Van Doorsselaere
Abstract:
Magnetic twist is thought to play an important role in many structures of the solar atmosphere. One of the effects of twist is to modify the properties of the eigenmodes of magnetic tubes. In the linear regime standing kink solutions are characterized by a change in polarization of the transverse displacement along the twisted tube. In the nonlinear regime, magnetic twist affects the development o…
▽ More
Magnetic twist is thought to play an important role in many structures of the solar atmosphere. One of the effects of twist is to modify the properties of the eigenmodes of magnetic tubes. In the linear regime standing kink solutions are characterized by a change in polarization of the transverse displacement along the twisted tube. In the nonlinear regime, magnetic twist affects the development of shear instabilities that appear at the tube boundary when it is oscillating laterally. These Kelvin-Helmholtz instabilities (KHI) are produced either by the jump in the azimuthal component of the velocity at the edge of the sharp boundary between the internal and external part of the tube, or either by the continuous small length-scales produced by phase-mixing when there is a smooth inhomogeneous layer. In this work the effect of twist is consistently investigated by solving the time-dependent problem including the process of energy transfer to the inhomogeneous layer. It is found that twist always delays the appearance of the shear instability but for tubes with thin inhomogeneous layers the effect is relatively small for moderate values of twist. On the contrary, for tubes with thick layers, the effect of twist is much stronger. This can have some important implications regarding observations of transverse kink modes and the KHI itself.
△ Less
Submitted 19 December, 2017;
originally announced December 2017.
-
Generalized phase mixing: Turbulence-like behaviour from unidirectionally propagating MHD waves
Authors:
Norbert Magyar,
Tom Van Doorsselaere,
Marcel Goossens
Abstract:
We present the results of three-dimensional (3D) ideal magnetohydrodynamics (MHD) simulations on the dynamics of a perpendicularly inhomogeneous plasma disturbed by propagating Alfvénic waves. Simpler versions of this scenario have been extensively studied as the phenomenon of phase mixing. We show that, by generalizing the textbook version of phase mixing, interesting phenomena are obtained, such…
▽ More
We present the results of three-dimensional (3D) ideal magnetohydrodynamics (MHD) simulations on the dynamics of a perpendicularly inhomogeneous plasma disturbed by propagating Alfvénic waves. Simpler versions of this scenario have been extensively studied as the phenomenon of phase mixing. We show that, by generalizing the textbook version of phase mixing, interesting phenomena are obtained, such as turbulence-like behavior and complex current-sheet structure, a novelty in longitudinally homogeneous plasma excited by unidirectionally propagating waves. This constitutes an important finding for turbulence-related phenomena in astrophysics in general, relaxing the conditions that have to be fulfilled in order to generate turbulent behavior.
△ Less
Submitted 2 November, 2017; v1 submitted 8 February, 2017;
originally announced February 2017.
-
Damping of nonlinear standing kink oscillations: a numerical study
Authors:
N. Magyar,
T. Van Doorsselaere
Abstract:
We aim to study the standing fundamental kink mode of coronal loops in the nonlinear regime, investigating the changes in energy evolution in the cross-section and oscillation amplitude of the loop which are related to nonlinear effects, in particular to the development of the Kelvin-Helmholtz instability (KHI). We run idea, high-resolution three-dimensional (3D) magnetohydrodynamics (MHD) simulat…
▽ More
We aim to study the standing fundamental kink mode of coronal loops in the nonlinear regime, investigating the changes in energy evolution in the cross-section and oscillation amplitude of the loop which are related to nonlinear effects, in particular to the development of the Kelvin-Helmholtz instability (KHI). We run idea, high-resolution three-dimensional (3D) magnetohydrodynamics (MHD) simulations, studying the influence of the initial velocity amplitude and the inhomogeneous layer thickness. We model the coronal loop as a straight, homogeneous magnetic flux tube with an outer inhomogeneous layer, embedded in a straight, homogeneous magnetic field. We find that, for low amplitudes which do not allow for the KHI to develop during the simulated time, the damping time agrees with the theory of resonant absorption. However, for higher amplitudes, the presence of KHI around the oscillating loop can alter the loop's evolution, resulting in a significantly faster damping than predicted by the linear theory in some cases. This questions the accuracy of seismological methods applied to observed damping profiles, based on linear theory.
△ Less
Submitted 22 September, 2016;
originally announced September 2016.
-
The instability and non-existence of multi-stranded loops, when driven by transverse waves
Authors:
N. Magyar,
T. Van Doorsselaere
Abstract:
In recent years, omni-present transverse waves have been observed in all layers of the solar atmosphere. Coronal loops are often modeled as a collection of individual strands, in order to explain their thermal behaviour and appearance. We perform 3D ideal MHD simulations to study the effect of a continuous small amplitude transverse footpoint driving on the internal structure of a coronal loop com…
▽ More
In recent years, omni-present transverse waves have been observed in all layers of the solar atmosphere. Coronal loops are often modeled as a collection of individual strands, in order to explain their thermal behaviour and appearance. We perform 3D ideal MHD simulations to study the effect of a continuous small amplitude transverse footpoint driving on the internal structure of a coronal loop composed of strands. The output is also converted to synthetic images, corresponding to the AIA 171 A and 193 A passbands, using FoMo. We show that the multi-stranded loop ceases to exist in the traditional sense of the word, because the plasma is efficiently mixed perpendicularly to the magnetic field, with the Kelvin-Helmholtz instability acting as the main mechanism. The final product of our simulation is mixed loop with density structures on a large range of scales, resembling a power-law. Thus, multi-stranded loops are unstable to driving by transverse waves, and this raises a strong doubt on the usability and applicability of coronal loop models consisting of independent strands.
△ Less
Submitted 14 April, 2016;
originally announced April 2016.
-
Forward Modelling of Propagating Slow Waves in Coronal Loops and Their Frequency-Dependent Damping
Authors:
Sudip Mandal,
Norbert Magyar,
Ding Yuan,
Tom Van Doorsselaere,
Dipankar Banerjee
Abstract:
Propagating slow waves in coronal loops exhibit a damping which depends upon the frequency of the waves. In this study we aim to investigate the relationship of the damping length (L$_d$) with the frequency of the propagating wave. We present a 3-D coronal loop model with uniform density and temperature and investigate the frequency dependent damping mechanism for the four chosen wave periods. We…
▽ More
Propagating slow waves in coronal loops exhibit a damping which depends upon the frequency of the waves. In this study we aim to investigate the relationship of the damping length (L$_d$) with the frequency of the propagating wave. We present a 3-D coronal loop model with uniform density and temperature and investigate the frequency dependent damping mechanism for the four chosen wave periods. We include the thermal conduction to damp the waves as they propagate through the loop. The numerical model output has been forward modelled to generate synthetic images of SDO/AIA 171 Å and 193 Å channels. The use of forward modelling, which incorporates the atomic emission properties into the intensity images, allows us to directly compare our results with the real observations. The results show that the damping lengths vary linearly with the periods. We also measure the contributions of the emission properties on the damping lengths by using density values from the simulation. In addition to that} we have also calculated the theoretical dependence of L$_d$ with wave periods and showed that it is consistent with the results we obtained from the numerical modelling and earlier observations.
△ Less
Submitted 1 February, 2016;
originally announced February 2016.
-
Numerical simulations of transverse oscillations in radiatively cooling coronal loops
Authors:
N. Magyar,
T. Van Doorsselaere,
A. Marcu
Abstract:
We aim to study the influence of radiative cooling on the standing kink oscillations of a coronal loop. Using the FLASH code, we solved the 3D ideal magnetohydrodynamic equations. Our model consists of a straight, density enhanced and gravitationally stratified magnetic flux tube. We perturbed the system initially, leading to a transverse oscillation of the structure, and followed its evolution fo…
▽ More
We aim to study the influence of radiative cooling on the standing kink oscillations of a coronal loop. Using the FLASH code, we solved the 3D ideal magnetohydrodynamic equations. Our model consists of a straight, density enhanced and gravitationally stratified magnetic flux tube. We perturbed the system initially, leading to a transverse oscillation of the structure, and followed its evolution for a number of periods. A realistic radiative cooling is implemented. Results are compared to available analytical theory. We find that in the linear regime (i.e. low amplitude perturbation and slow cooling) the obtained period and damping time are in good agreement with theory. The cooling leads to an amplification of the oscillation amplitude. However, the difference between the cooling and non-cooling cases is small (around 6% after 6 oscillations). In high amplitude runs with realistic cooling, instabilities deform the loop, leading to increased damping. In this case, the difference between cooling and non-cooling is still negligible at around 12%. A set of simulations with higher density loops are also performed, to explore what happens when the cooling takes place in a very short time (tcool = 100 s). We strengthen the results of previous analytical studies that state that the amplification due to cooling is ineffective, and its influence on the oscillation characteristics is small, at least for the cases shown here. Furthermore, the presence of a relatively strong damping in the high amplitude runs even in the fast cooling case indicates that it is unlikely that cooling could alone account for the observed, flare-related undamped oscillations of coronal loops. These results may be significant in the field of coronal seismology, allowing its application to coronal loop oscillations with observed fading-out or cooling behaviour.
△ Less
Submitted 29 October, 2015;
originally announced October 2015.
-
Appearance of an accretion disk perturbed by fractional Brownian Motion density
Authors:
Gabriela Mocanu,
Norbert Magyar,
Anabele Pardi,
Alexandru Marcu
Abstract:
This paper aims to investigate/map the effects that perturbations applied to an accretion disk might produce on the registered Light Curves (LC). The case of accretion disks around supermassive active black holes (AGNs) is studied with the goal to explain some of the statistical properties of the observed IntraDay Variability (IDV). The region producing optical IDV is perturbed by allowing it to d…
▽ More
This paper aims to investigate/map the effects that perturbations applied to an accretion disk might produce on the registered Light Curves (LC). The case of accretion disks around supermassive active black holes (AGNs) is studied with the goal to explain some of the statistical properties of the observed IntraDay Variability (IDV). The region producing optical IDV is perturbed by allowing it to develop a mass density of a fractional Brownian Motion-like type. The light curves and spectral slopes are calculated and compared to observational data for different Hurst parameters. The spectral slopes of the simulated light curves vary in the range $(0.4,2.5)$. The agreement with observational data shows that a magnetized disk subjected to stochastic perturbations can produce some of the features observed in the light curves.
△ Less
Submitted 1 February, 2014;
originally announced February 2014.