An official website of the United States government
The properties of collisionless shocks, like the density jump, are usually derived from magnetohydrodynamics (MHD), where isotropic pressures are assumed. Yet, in a collisionless plasma, an external magnetic field can sustain a stable anisotropy. We have already devised a model for the kinetic history of the plasma through the shock front ( J. Plasma Phys. , vol. 84, issue 6, 2018, 905840604), allowing to self-consistently compute the downstream anisotropy, and hence the density jump, in terms of the upstream parameters. This model deals with the case of a parallel shock, where the magnetic field is normal to the front both in the upstream and the downstream. Yet, MHD also allows for shock solutions, the so-called switch-on solutions, where the field is normal to the front only in the upstream. This article consists in applying our model to these switch-on shocks. While MHD offers only one switch-on solution within a limited range of Alfvén Mach numbers, our model offers two kinds of solutions within a slightly different range of Alfvén Mach numbers. These two solutions are most likely the outcome of the intermediate and fast MHD shocks under our model. While the intermediate and fast shocks merge in MHD for the parallel case, they do not within our model. For simplicity, the formalism is restricted to non-relativistic shocks in pair plasmas where the upstream is cold.
more »
« less
Kinetic simulations of strongly magnetized parallel shocks: deviations from MHD jump conditions
ABSTRACT Shocks waves are a ubiquitous feature of many astrophysical plasma systems, and an important process for energy dissipation and transfer. The physics of these shock waves are frequently treated/modelled as a collisional, fluid magnetohydrodynamic (MHD) discontinuity, despite the fact that many shocks occur in the collisionless regime. In light of this, using fully kinetic, 3D simulations of non-relativistic, parallel propagating collisionless shocks comprised of electron-positron plasma, we detail the deviation of collisionless shocks form MHD predictions for varying magnetization/Alfvénic Mach numbers, with particular focus on systems with Alfénic Mach numbers much smaller than sonic Mach numbers. We show that the shock compression ratio decreases for sufficiently large upstream magnetic fields, in agreement with theoretical predictions from previous works. Additionally, we examine the role of magnetic field strength on the shock front width. This work reinforces a growing body of work that suggest that modelling many astrophysical systems with only a fluid plasma description omits potentially important physics.
more »
« less
- Award ID(s):
- 1936393
- PAR ID:
- 10348637
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 509
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 2084 to 2090
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The study of collisionless shocks and their role in cosmic ray acceleration has gained importance through observations and simulations, driving interest in reproducing these conditions in laboratory experiments using high-power lasers. In this work, we examine the role of three-dimensional (3D) effects in ion acceleration in quasi-perpendicular shocks under laboratory-relevant conditions. Using hybrid particle-in-cell (PIC) simulations (kinetic ions and fluid electrons), we explore how the Alfvénic and sonic Mach numbers, along with plasma beta, influence ion energization, unlocked only in 3D, and establish scaling criteria for when conducting 3D simulations is necessary. Our results show that efficient ion acceleration requires Alfvénic Mach numbers ≥25 and sonic Mach numbers ≥13, with plasma-β≤5. We theoretically found that, while two-dimensional (2D) simulations suffice for current laboratory-accessible shock conditions, 3D effects become crucial for shock velocities exceeding 1000 km/s and experiments sustaining the shock for at least 10 ns. We surveyed previous laboratory experiments on collisionless shocks and found that 3D effects are unimportant under those conditions, implying that one-dimensional and 2D simulations should be enough to model the accelerated ion spectra. However, we do find that the same experiments are realistically close to accessing the regime relevant to 3D effects, an exciting prospect for future laboratory efforts. We propose modifications to past experimental configurations to optimize and control 3D effects on ion acceleration. These proposed experiments could be used to benchmark plasma astrophysics kinetic codes and/or employed as controllable sources of energetic particles.more » « less
-
Abstract Collisionless low-Mach-number shocks are abundant in astrophysical and space plasma environments, exhibiting complex wave activity and wave–particle interactions. In this paper, we present 2D Particle-in-Cell (PIC) simulations of quasi-perpendicular nonrelativistic (vsh≈ (5500–22000) km s−1) low-Mach-number shocks, with a specific focus on studying electrostatic waves in the shock ramp and precursor regions. In these shocks, an ion-scale oblique whistler wave creates a configuration with two hot counterstreaming electron beams, which drive unstable electron acoustic waves (EAWs) that can turn into electrostatic solitary waves (ESWs) at the late stage of their evolution. By conducting simulations with periodic boundaries, we show that the EAW properties agree with linear dispersion analysis. The characteristics of ESWs in shock simulations, including their wavelength and amplitude, depend on the shock velocity. When extrapolated to shocks with realistic velocities (vsh≈ 300 km s−1), the ESW wavelength is reduced to one-tenth of the electron skin depth and the ESW amplitude is anticipated to surpass that of the quasi-static electric field by more than a factor of 100. These theoretical predictions may explain a discrepancy, between PIC and satellite measurements, in the relative amplitude of high- and low-frequency electric field fluctuations.more » « less
-
Shock waves in plasma are usually dealt with using magnetohydrodynamics (MHD). Yet, MHD entails the assumption of a short mean free path, which is not fulfilled in a collisionless plasma. Recently, for pair plasmas, we devised a model allowing one to account for kinetic effects within a MHD-like formalism. Its relies on an estimate of the anisotropy generated when crossing the front, with a subsequent assessment of the stability of this anisotropy in the downstream. We solved our model for parallel, perpendicular and switch-on shocks. Here we bridge between all these cases by treating the problem of an arbitrarily, but coplanar, oriented magnetic field. Even though the formalism presented is valid for anisotropic upstream temperatures, only the case of a cold upstream is solved. We find extra solutions which are not part of the MHD catalogue, and a density jump that is notably less in the quasi-parallel, highly magnetized, regime. Given the complexity of the calculations, this work is mainly devoted to the presentation of the mathematical aspect of our model. A forthcoming article will be devoted to the physics of the shocks here defined.more » « less
-
We investigate three-dimensional (3-D) bow shocks in a highly collisional magnetized aluminium plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ( $$M_S \sim 7$$ ), super-Alfvénic ( $$M_A > 1$$ ) magnetized flows with frozen-in magnetic flux ( $$R_M \gg 1$$ ). These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. The 3-D simulations performed using the resistive magnetohydrodynamic (MHD) code Gorgon confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/stab3110
