An official website of the United States government
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
Criteria for ion acceleration in laboratory magnetized quasi-perpendicular collisionless shocks: When are 2D simulations enough?
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
- Award ID(s):
- 2206607
- PAR ID:
- 10591974
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- Physics of Plasmas
- Volume:
- 32
- Issue:
- 5
- ISSN:
- 1070-664X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.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
-
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
-
Abstract Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions that reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase-space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique’s utility for simulation, laboratory, and spacecraft analysis.more » « less
