An official website of the United States government
Abstract We study the linear stability of a planar interface separating two fluids in relative motion, focusing on conditions appropriate for the boundaries of relativistic jets. The jet is magnetically dominated, whereas the ambient wind is gas-pressure-dominated. We derive the most general form of the dispersion relation and provide an analytical approximation of its solution for an ambient sound speed much smaller than the jet Alfvén speedvA, as appropriate for realistic systems. The stability properties are chiefly determined by the angleψbetween the wavevector and the jet magnetic field. Forψ=π/2, magnetic tension plays no role, and our solution resembles the one of a gas-pressure-dominated jet. Here, only sub-Alfvénic jets are unstable ( , wherevis the shear velocity andθthe angle between the velocity and the wavevector). Forψ= 0, the free energy in the velocity shear needs to overcome the magnetic tension, and only super-Alfvénic jets are unstable ( , with Γwthe wind adiabatic index). Our results have important implications for the propagation and emission of relativistic magnetized jets.
more »
« less
Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows
ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $$M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$$, where $$v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.
more »
« less
- Award ID(s):
- 2108201
- PAR ID:
- 10429325
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 524
- Issue:
- 1
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 90-99
- Size(s):
- p. 90-99
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We investigate the validity of Taylor’s hypothesis (TH) in the analysis of velocity and magnetic field fluctuations in Alfvénic solar wind streams measured by Parker Solar Probe (PSP) during the first four encounters. The analysis is based on a recent model of the spacetime correlation of magnetohydrodynamic (MHD) turbulence, which has been validated in high-resolution numerical simulations of strong reduced MHD turbulence. We use PSP velocity and magnetic field measurements from 24 h intervals selected from each of the first four encounters. The applicability of TH is investigated by measuring the parameter ϵ = δu 0 /√2 V ⊥ , which quantifies the ratio between the typical speed of large-scale fluctuations, δu 0 , and the local perpendicular PSP speed in the solar wind frame, V ⊥ . TH is expected to be applicable for ϵ ≲ 0.5 when PSP is moving nearly perpendicular to the local magnetic field in the plasma frame, irrespective of the Alfvén Mach number M A = V SW ∕ V A , where V SW and V A are the local solar wind and Alfvén speed, respectively. For the four selected solar wind intervals, we find that between 10 and 60% of the time, the parameter ϵ is below 0.2 and the sampling angle (between the spacecraft velocity in the plasma frame and the local magnetic field) is greater than 30°. For angles above 30°, the sampling direction is sufficiently oblique to allow one to reconstruct the reduced energy spectrum E ( k ⊥ ) of magnetic fluctuations from its measured frequency spectra. The spectral indices determined from power-law fits of the measured frequency spectrum accurately represent the spectral indices associated with the underlying spatial spectrum of turbulent fluctuations in the plasma frame. Aside from a frequency broadening due to large-scale sweeping that requires careful consideration, the spatial spectrum can be recovered to obtain the distribution of fluctuation’s energy across scales in the plasma frame.more » « less
-
Alfvén waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfvén waves, showing that an $$E_{B_{\perp }}(k_{\perp }) \propto k_{\perp }^{-2}$$ energy spectrum develops as a result of the weak turbulence cascade in relativistic magnetohydrodynamics and its infinitely magnetized (force-free) limit. The plasma turbulence ubiquitously generates current sheets, which act as locations where magnetic energy dissipates. We show that current sheets form as a natural result of nonlinear interactions between counter-propagating Alfvén waves. These current sheets form owing to the compression of elongated eddies, driven by the shear induced by growing higher-order modes, and undergo a thinning process until they break-up into small-scale turbulent structures. We explore the formation of current sheets both in overlapping waves and in localized wave packet collisions. The relativistic interaction of localized Alfvén waves induces both Alfvén waves and fast waves, and efficiently mediates the conversion and dissipation of electromagnetic energy in astrophysical systems. Plasma energization through reconnection in current sheets emerging during the interaction of Alfvén waves can potentially explain X-ray emission in black hole accretion coronae and neutron star magnetospheres.more » « less
-
A growing body of evidence suggests that the solar wind is powered to a large extent by an Alfvén-wave (AW) energy flux. AWs energize the solar wind via two mechanisms: heating and work. We use high-resolution direct numerical simulations of reflection-driven AW turbulence (RDAWT) in a fast-solar-wind stream emanating from a coronal hole to investigate both mechanisms. In particular, we compute the fraction of the AW power at the coronal base ( $$P_\textrm {AWb}$$ ) that is transferred to solar-wind particles via heating between the coronal base and heliocentric distance $$r$$ , which we denote by $$\chi _{H}(r)$$ , and the fraction that is transferred via work, which we denote by $$\chi _{W}(r)$$ . We find that $$\chi _{W}(r_{A})$$ ranges from 0.15 to 0.3, where $$r_{A}$$ is the Alfvén critical point. This value is small compared with one because the Alfvén speed $$v_{A}$$ exceeds the outflow velocity $$U$$ at $$r < r_{A}$$ , so the AWs race through the plasma without doing much work. At $$r>r_{A}$$ , where $$v_{A} < U$$ , the AWs are in an approximate sense ‘stuck to the plasma’, which helps them do pressure work as the plasma expands. However, much of the AW power has dissipated by the time the AWs reach $$r=r_{A}$$ , so the total rate at which AWs do work on the plasma at $$r>r_{A}$$ is a modest fraction of $$P_\textrm {AWb}$$ . We find that heating is more effective than work at $$r < r_{A}$$ , with $$\chi _{H}(r_{A})$$ ranging from 0.5 to 0.7. The reason that $$\chi _{H} \geq 0.5$$ in our simulations is that an appreciable fraction of the local AW power dissipates within each Alfvén-speed scale height in RDAWT, and there are a few Alfvén-speed scale heights between the coronal base and $$r_{A}$$ . A given amount of heating produces more magnetic moment in regions of weaker magnetic field. Thus, paradoxically, the average proton magnetic moment increases robustly with increasing $$r$$ at $$r>r_{A}$$ , even though the total rate at which AW energy is transferred to particles at $$r>r_{A}$$ is a small fraction of $$P_\textrm {AWb}$$ .more » « less
-
Abstract Particle acceleration behind a shock wave due to interactions between magnetic islands in the heliosphere has attracted attention in recent years. The downstream acceleration may yield a continuous increase of particle flux downstream of the shock wave. Although it is not obvious how the downstream magnetic islands are produced, it has been suggested that current sheets are involved in the generation of magnetic islands due to their interaction with a shock wave. We perform 2D hybrid kinetic simulations to investigate the interaction between multiple current sheets and a shock wave. In the simulation, current sheets are compressed by the shock wave and a tearing instability develops at the compressed current sheets downstream of the shock. As the result of this instability, the electromagnetic fields become turbulent and magnetic islands form well downstream of the shock wave. We find a “post-cursor” region in which the downstream flow speed normal to the shock wave in the downstream rest frame is decelerated to ∼ 1 V A immediately behind the shock wave, where V A is the upstream Alfvén speed. The flow speed then gradually decelerates to 0 accompanied by the development of the tearing instability. We also observe an efficient production of energetic particles above 100 E 0 during the development of the instability some distance downstream of the shock wave, where E 0 = m p V A 2 and m p is the proton mass. This feature corresponds to Voyager observations showing that the anomalous cosmic-ray intensity increase begins some distance downstream of the heliospheric termination shock.more » « less
