-
Simulation Models for Exploring Magnetic Reconnection
Authors:
Michael Shay,
Subash Adhikari,
Naoki Beesho,
Joachim Birn,
Jorg Buechner,
Paul Cassak,
Li-Jen Chen,
Yuxi Chen,
Giulia Cozzani,
Jim Drake,
Fan Guo,
Michael Hesse,
Neeraj Jain,
Yann Pfau-Kempf,
Yu Lin,
Yi-Hsin Liu,
Mitsuo Oka,
Yuri A. Omelchenko,
Minna Palmroth,
Oreste Pezzi,
Patricia H. Reiff,
Marc Swisdak,
Frank Toffoletto,
Gabor Toth,
Richard A. Wolf
Abstract:
Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This chapter reviews many of these…
▽ More
Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This chapter reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).
△ Less
Submitted 9 June, 2024;
originally announced June 2024.
-
Flux-tube-dependent propagation of Alfvén waves in the solar corona
Authors:
Chaitanya Prasad Sishtla,
Jens Pomoell,
Emilia Kilpua,
Simon Good,
Farhad Daei,
Minna Palmroth
Abstract:
Alfvén-wave turbulence has emerged as an important heating mechanism to accelerate the solar wind. The generation of this turbulent heating is dependent on the presence and subsequent interaction of counter-propagating alfvén waves. This requires us to understand the propagation and evolution of alfvén waves in the solar wind in order to develop an understanding of the relationship between turbule…
▽ More
Alfvén-wave turbulence has emerged as an important heating mechanism to accelerate the solar wind. The generation of this turbulent heating is dependent on the presence and subsequent interaction of counter-propagating alfvén waves. This requires us to understand the propagation and evolution of alfvén waves in the solar wind in order to develop an understanding of the relationship between turbulent heating and solar-wind parameters. We aim to study the response of the solar wind upon injecting monochromatic single-frequency alfvén waves at the base of the corona for various magnetic flux-tube geometries. We used an ideal magnetohydrodynamic (MHD) model using 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. Alfvén waves were found to be reflected due to the development of the parametric decay instability (PDI). Further investigation revealed that the PDI was suppressed both by efficient reflections at low frequencies as well as magnetic flux-tube geometries.
△ Less
Submitted 20 April, 2022;
originally announced April 2022.
-
ULF wave transmission across collisionless shocks: 2.5D local hybrid simulations
Authors:
Primoz Kajdic,
Yann Pfau-Kempf,
Lucile Turc,
Andrew P Dimmock,
Minna Palmroth,
Kazue Takahashi,
Eemilia Kilpua,
Jan Soucek,
Naoko Takahashi,
Luis Preisser,
Xochitl Blanco-Cano,
Domenico Trotta,
David Burgess
Abstract:
We study the interaction of upstream ultra-low frequency (ULF) waves with collisionless shocks by analyzing the outputs of eleven 2D local hybrid simulation runs. Our simulated shocks have Alfvénic Mach numbers between 4.29-7.42 and their $θ_{BN}$ angles are 15$^\circ$, 30$^\circ$, 45$^\circ$ and 50$^\circ$. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude…
▽ More
We study the interaction of upstream ultra-low frequency (ULF) waves with collisionless shocks by analyzing the outputs of eleven 2D local hybrid simulation runs. Our simulated shocks have Alfvénic Mach numbers between 4.29-7.42 and their $θ_{BN}$ angles are 15$^\circ$, 30$^\circ$, 45$^\circ$ and 50$^\circ$. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude of the upstream waves exhibit a complex dependence on the shock's M$_A$ and $θ_{BN}$. The wavelength positively correlates with both parameters, with the dependence on $θ_{BN}$ being much stronger. The amplitude of the ULF waves is proportional to the product of the reflected beam velocity and density, which also depend on M$_A$ and $θ_{BN}$. The interaction of the ULF waves with the shock causes large-scale (several tens of upstream ion inertial lengths) shock rippling. The properties of the shock ripples are related to the ULF wave properties, namely thier wavelength and amplitude. In turn, the ripples have a large impact on the ULF wave transmission across the shock because they change local shock properties ($θ_{BN}$, strength), so that different sections of the same ULF wave front encounter shock with different characteristics. Downstream fluctuations do not resemble the upstream waves in terms the wavefront extension, orientation or their wavelength. However some features are conserved in the Fourier spectra of downstream compressive waves that present a bump or flattening at wavelengths approximately corresponding to those of the upstream ULF waves. In the transverse downstream spectra these features are weaker.
△ Less
Submitted 25 January, 2022;
originally announced January 2022.
-
Electron acceleration and radio emission following the early interaction of two coronal mass ejections
Authors:
D. E. Morosan,
E. Palmerio,
J. E. Räsänen,
E. K. J. Kilpua,
J. Magdalenić,
B. J. Lynch,
A. Kumari,
J. Pomoell,
M. Palmroth
Abstract:
Context. Coronal mass ejections (CMEs) are large eruptions of magnetised plasma from the Sun that are often accompanied by solar radio bursts produced by accelerated electrons. Aims. A powerful source for accelerating electron beams are CME-driven shocks, however, there are other mechanisms capable of accelerating electrons during a CME eruption. So far, studies have relied on the traditional clas…
▽ More
Context. Coronal mass ejections (CMEs) are large eruptions of magnetised plasma from the Sun that are often accompanied by solar radio bursts produced by accelerated electrons. Aims. A powerful source for accelerating electron beams are CME-driven shocks, however, there are other mechanisms capable of accelerating electrons during a CME eruption. So far, studies have relied on the traditional classification of solar radio bursts into five groups (Type I-V) based mainly on their shapes and characteristics in dynamic spectra. Here, we aim to determine the origin of moving radio bursts associated with a CME that do not fit into the present classification of the solar radio emission. Methods. By using radio imaging from the Nançay Radioheliograph, combined with observations from the Solar Dynamics Observatory, Solar and Heliospheric Observatory, and Solar Terrestrial Relations Observatory spacecraft, we investigate the moving radio bursts accompanying two subsequent CMEs on 22 May 2013. We use three-dimensional reconstructions of the two associated CME eruptions to show the possible origin of the observed radio emission. Results. We identified three moving radio bursts at unusually high altitudes in the corona that are located at the northern CME flank and move outwards synchronously with the CME. The radio bursts correspond to fine-structured emission in dynamic spectra with durations of ~1 s, and they may show forward or reverse frequency drifts. Since the CME expands closely following an earlier CME, a low coronal CME-CME interaction is likely responsible for the observed radio emission.
△ Less
Submitted 1 September, 2020; v1 submitted 24 August, 2020;
originally announced August 2020.
-
Three-dimensional reconstruction of multiple particle acceleration regions during a coronal mass ejection
Authors:
D. E. Morosan,
E. Palmerio,
J. Pomoell,
R. Vainio,
M. Palmroth,
E. K. J. Kilpua
Abstract:
Context. Some of the most prominent sources for particle acceleration in our Solar System are large eruptions of magnetised plasma from the Sun called coronal mass ejections (CMEs). These accelerated particles can generate radio emission through various mechanisms. Aims. CMEs are often accompanied by a variety of solar radio bursts with different shapes and characteristics in dynamic spectra. Radi…
▽ More
Context. Some of the most prominent sources for particle acceleration in our Solar System are large eruptions of magnetised plasma from the Sun called coronal mass ejections (CMEs). These accelerated particles can generate radio emission through various mechanisms. Aims. CMEs are often accompanied by a variety of solar radio bursts with different shapes and characteristics in dynamic spectra. Radio bursts directly associated with CMEs often show movement in the direction of CME expansion. Here, we aim to determine the emission mechanism of multiple moving radio bursts that accompanied a flare and CME that took place on 14 June 2012. Methods. We used radio imaging from the Nancay Radioheliograph, combined with observations from the Solar Dynamics Observatory and Solar Terrestrial Relations Observatory spacecraft, to analyse these moving radio bursts in order to determine their emission mechanism and three-dimensional (3D) location with respect to the expanding CME. Results. In using a 3D representation of the particle acceleration locations in relation to the overlying coronal magnetic field and the CME propagation, for the first time, we provide evidence that these moving radio bursts originate near the CME flanks and some that are possible signatures of shock-accelerated electrons following the fast CME expansion in the low corona. Conclusions. The moving radio bursts, as well as other stationary bursts observed during the eruption, occur simultaneously with a type IV continuum in dynamic spectra, which is not usually associated with emission at the CME flanks. Our results show that moving radio bursts that could traditionally be classified as moving type IVs can represent shock signatures associated with CME flanks or plasma emission inside the CME behind its flanks, which are closely related to the lateral expansion of the CME in the low corona.
△ Less
Submitted 23 January, 2020;
originally announced January 2020.
-
Particle Telescope aboard FORESAIL-1: simulated performance
Authors:
Philipp Oleynik,
Rami Vainio,
Hannu-Pekka Hedman,
Arttu Punkkinen,
Risto Punkkinen,
Lassi Salomaa,
Tero Säntti,
Jarno Tuominen,
Pasi Virtanen,
Alexandre Bosser,
Pekka Janhunen,
Emilia Kilpua,
Minna Palmroth,
Jaan Praks,
Andris Slavinskis,
Syed R. U. Kakakhel,
Juhani Peltonen,
Juha Plosila,
Jani Tammi,
Hannu Tenhunen,
Tomi Westerlund
Abstract:
The Particle Telescope (PATE) of FORESAIL-1 mission is described. FORESAIL-1 is a CubeSat mission to polar Low Earth Orbit. Its scientific objectives are to characterize electron precipitation from the radiation belts and to observe energetic neutral atoms (ENAs) originating from the Sun during the strongest solar flares. For that purpose, the 3-unit CubeSat carries a particle telescope that measu…
▽ More
The Particle Telescope (PATE) of FORESAIL-1 mission is described. FORESAIL-1 is a CubeSat mission to polar Low Earth Orbit. Its scientific objectives are to characterize electron precipitation from the radiation belts and to observe energetic neutral atoms (ENAs) originating from the Sun during the strongest solar flares. For that purpose, the 3-unit CubeSat carries a particle telescope that measures energetic electrons in the nominal energy range of 80--800 keV in seven energy channels and energetic protons at 0.3--10 MeV in ten channels. In addition, particles penetrating the whole telescope at higher energies will be measured in three channels: one $>$800 keV electron channel, two integral proton channels at $>$10 MeV energies. The instrument contains two telescopes at right angles to each other, one measuring along the spin axis of the spacecraft and one perpendicular to it. During a spin period (nominally 15 s), the rotating telescope will, thus, deliver angular distributions of protons and electrons, at 11.25-degree clock-angle resolution, which enables one to accurately determine the pitch-angle distribution and separate the trapped and precipitating particles. During the last part of the mission, the rotation axis will be accurately pointed toward the Sun, enabling the measurement of the energetic hydrogen from that direction. Using the geomagnetic field as a filter and comparing the rates observed by the two telescopes, the instrument can observe the solar ENA flux for events similar to the only one so far observed in December 2006. We present the Geant4-simulated energy and angular response functions of the telescope and assess its sensitivity showing that they are adequate to address the scientific objectives of the mission.
△ Less
Submitted 14 November, 2019;
originally announced November 2019.
-
Vlasov methods in space physics and astrophysics
Authors:
Minna Palmroth,
Urs Ganse,
Yann Pfau-Kempf,
Markus Battarbee,
Lucile Turc,
Thiago Brito,
Maxime Grandin,
Sanni Hoilijoki,
Arto Sandroos,
Sebastian von Alfthan
Abstract:
This paper reviews Vlasov-based numerical methods used to model plasma in space physics and astrophysics. Plasma consists of collectively behaving charged particles that form the major part of baryonic matter in the Universe. Many concepts ranging from our own planetary environment to the Solar system and beyond can be understood in terms of kinetic plasma physics, represented by the Vlasov equati…
▽ More
This paper reviews Vlasov-based numerical methods used to model plasma in space physics and astrophysics. Plasma consists of collectively behaving charged particles that form the major part of baryonic matter in the Universe. Many concepts ranging from our own planetary environment to the Solar system and beyond can be understood in terms of kinetic plasma physics, represented by the Vlasov equation. We introduce the physical basis for the Vlasov system, and then outline the associated numerical methods that are typically used. A particular application of the Vlasov system is Vlasiator, the world's first global hybrid-Vlasov simulation for the Earth's magnetic domain, the magnetosphere. We introduce the design strategies for Vlasiator and outline its numerical concepts ranging from solvers to coupling schemes. We review Vlasiator's parallelisation methods and introduce the used high-performance computing (HPC) techniques. A short review of verification, validation and physical results is included. The purpose of the paper is to present the Vlasov system and introduce an example implementation, and to illustrate that even with massive computational challenges, an accurate description of physics can be rewarding in itself and significantly advance our understanding. Upcoming supercomputing resources are making similar efforts feasible in other fields as well, making our design options relevant for others facing similar challenges.
△ Less
Submitted 17 August, 2018;
originally announced August 2018.