An official website of the United States government
ABSTRACT Recently, cosmic rays (CRs) have emerged as a leading candidate for driving galactic winds. Small-scale processes can dramatically affect global wind properties. We run two-moment simulations of CR streaming to study how sound waves are driven unstable by phase-shifted CR forces and CR heating. We verify linear theory growth rates. As the sound waves grow non-linear, they steepen into a quasi-periodic series of propagating shocks; the density jumps at shocks create CR bottlenecks. The depth of a propagating bottleneck depends on both the density jump and its velocity; ΔPc is smaller for rapidly moving bottlenecks. A series of bottlenecks creates a CR staircase structure, which can be understood from a convex hull construction. The system reaches a steady state between growth of new perturbations, and stair mergers. CRs are decoupled at plateaus, but exert intense forces and heating at stair jumps. The absence of CR heating at plateaus leads to cooling, strong gas pressure gradients and further shocks. If bottlenecks are stationary, they can drastically modify global flows; if their propagation times are comparable to dynamical times, their effects on global momentum and energy transfer are modest. The CR acoustic instability is likely relevant in thermal interfaces between cold and hot gas, as well as galactic winds. Similar to increased opacity in radiative flows, the build-up of CR pressure due to bottlenecks can significantly increase mass outflow rates, by up to an order of magnitude. It seeds unusual forms of thermal instability, and the shocks could have distinct observational signatures, on ∼kpc scales.
more »
« less
Building a weak shockwave from linear modes
In shockwave theory, the density, velocity and pressure jumps are derived from the conservation equations. Here, we address the physics of a weak shock the other way around. We first show that the density profile of a weak shockwave in a fluid can be expressed as a sum of linear acoustic modes. The shock so built propagates at the speed of sound and matter is exactly conserved at the front crossing. Yet, momentum and energy are only conserved up to order 0 in powers of the shock amplitude. The density, velocity and pressure jumps are similar to those of a fluid shock, and an equivalent Mach number can be defined. A similar process is possible in magnetohydrodynamics. Yet, such a decomposition is found impossible for collisionless shocks due to the dispersive nature of ion acoustic waves. Weakly nonlinear corrections to their frequency do not solve the problem. Weak collisionless shocks could be inherently nonlinear, non-amenable to any linear superposition. Or they could be non-existent, as hinted by recent works.
more »
« less
- Award ID(s):
- 1816420
- PAR ID:
- 10337489
- Date Published:
- Journal Name:
- Journal of Plasma Physics
- Volume:
- 88
- Issue:
- 1
- ISSN:
- 0022-3778
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this paper, nonlinearity associated with intense ultrasound is studied by using the one-dimensional motion of nonlinear shock wave in an ideal fluid. In nonlinear acoustics, the wave speed of different segments of a waveform is different, which causes distortion in the waveform and can result in the formation of a shock (discontinuity). Acoustic pressure of high-intensity waves causes particles in the ideal fluid to vibrate forward and backward, and this disturbance is of relatively large magnitude due to high-intensities, which leads to nonlinearity in the waveform. In this research, this vibration of fluid due to the intense ultrasonic wave is modeled as a fluid pushed by one complete cycle of piston. In a piston cycle, as it moves forward, it causes fluid particles to compress, which may lead to the formation of a shock (discontinuity). Then as the piston retracts, a forward-moving rarefaction, a smooth fan zone of continuously changing pressure, density, and velocity is generated. When the piston stops at the end of the cycle, another shock is sent forward into the medium. The variation in wave speed over the entire waveform is calculated by solving a Riemann problem. This study examined the interaction of shocks with a rarefaction. The flow field resulting from these interactions shows that the shock waves are attenuated to a Mach wave, and the pressure distribution within the flow field shows the initial wave is dissipated. The developed theory is applied to waves generated by 20 KHz, 500 KHz, and 2 MHz transducers with 50, 150, 500, and 1500 W power levels to explore the effect of frequency and power on the generation and decay of shock waves. This work enhances the understanding of the interactions of high-intensity ultrasonic waves with fluids.more » « less
-
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
-
We have performed two-dimensional hybrid simulations of non-relativistic collisionless shocks in the presence of pre-existing energetic particles (‘seeds’); such a study applies, for instance, to the re-acceleration of galactic cosmic rays (CRs) in supernova remnant (SNR) shocks and solar wind energetic particles in heliospheric shocks. Energetic particles can be effectively reflected and accelerated regardless of shock inclination via a process that we call diffusive shock re-acceleration. We find that re-accelerated seeds can drive the streaming instability in the shock upstream and produce effective magnetic field amplification. This can eventually trigger the injection of thermal protons even at oblique shocks that ordinarily cannot inject thermal particles. We characterize the current in reflected seeds, finding that it tends to a universal value $$J\simeq en_{\text{CR}}v_{\text{sh}}$$ , where $$en_{\text{CR}}$$ is the seed charge density and $$v_{\text{sh}}$$ is the shock velocity. When applying our results to SNRs, we find that the re-acceleration of galactic CRs can excite the Bell instability to nonlinear levels in less than $${\sim}10~\text{yr}$$ , thereby providing a minimum level of magnetic field amplification for any SNR shock. Finally, we discuss the relevance of diffusive shock re-acceleration also for other environments, such as heliospheric shocks, galactic superbubbles and clusters of galaxies.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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/S0022377821001252
