-
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.
-
CME Magnetic Structure and IMF Preconditioning Affecting SEP Transport
Authors:
Erika Palmerio,
Emilia K. J. Kilpua,
Olivier Witasse,
David Barnes,
Beatriz Sánchez-Cano,
Andreas J. Weiss,
Teresa Nieves-Chinchilla,
Christian Möstl,
Lan K. Jian,
Marilena Mierla,
Andrei N. Zhukov,
Jingnan Guo,
Luciano Rodriguez,
Patrick J. Lowrance,
Alexey Isavnin,
Lucile Turc,
Yoshifumi Futaana,
Mats Holmström
Abstract:
Coronal mass ejections (CMEs) and solar energetic particles (SEPs) are two phenomena that can cause severe space weather effects throughout the heliosphere. The evolution of CMEs, especially in terms of their magnetic structure, and the configuration of the interplanetary magnetic field (IMF) that influences the transport of SEPs are currently areas of active research. These two aspects are not ne…
▽ More
Coronal mass ejections (CMEs) and solar energetic particles (SEPs) are two phenomena that can cause severe space weather effects throughout the heliosphere. The evolution of CMEs, especially in terms of their magnetic structure, and the configuration of the interplanetary magnetic field (IMF) that influences the transport of SEPs are currently areas of active research. These two aspects are not necessarily independent of each other, especially during solar maximum when multiple eruptive events can occur close in time. Accordingly, we present the analysis of a CME that erupted on 2012 May 11 (SOL2012-05-11) and an SEP event following an eruption that took place on 2012 May 17 (SOL2012-05-17). After observing the May 11 CME using remote-sensing data from three viewpoints, we evaluate its propagation through interplanetary space using several models. Then, we analyse in-situ measurements from five predicted impact locations (Venus, Earth, the Spitzer Space Telescope, the Mars Science Laboratory en route to Mars, and Mars) in order to search for CME signatures. We find that all in-situ locations detect signatures of an SEP event, which we trace back to the May 17 eruption. These findings suggest that the May 11 CME provided a direct magnetic connectivity for the efficient transport of SEPs. We discuss the space weather implications of CME evolution, regarding in particular its magnetic structure, and CME-driven IMF preconditioning that facilitates SEP transport. Finally, this work remarks the importance of using data from multiple spacecraft, even those that do not include space weather research as their primary objective.
△ Less
Submitted 10 February, 2021;
originally announced February 2021.
-
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.