SEARCH · Engineering Papers

Results for “drift wave turbulence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

142 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Scattering of radio frequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radio frequency (RF) electromagnetic waves that propagate through this region. It is important to know, whether used for diagnostics or for heating and current drive, the spectral properties of these scattered RF waves. The spectral changes influences the interpretation of the diagnostic-data obtained and the current and heating profiles. A full-wave, 3D electromagnetic code ScaRF (see Papadopoulos et al. 2019) has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma, anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can be for the study of scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U, DIII-D. For the case of medium-sized tokamaks, there's experimental evidence that drift waves and rippling modes are present in the edge region (see Ritz et al. 1984). Hence, we study the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes (see Papadopoulos et al. 2019). The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion (see Xiu & Karniadakis 2002) method in conjunction with the Smolyak sparse grid integration (see Papadopoulos et al. 2018) (PCE-SG). The PCE-SG method is proven accurate and much more efficient (roughly 2-orders of magnitude shorter execution time) compared to alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1991

Plasma convection and ion beam generation in the plasma sheet boundary layer

Because of the dawn-dusk electric field E(dd), plasma in the magnetotail convects from the lobe toward the central plasma sheet (CPS). In the absence of space or velocity diffusion due to plasma turbulence, convection would yield a steady state distribution function f = V exp (-2/3) g(v exp 2 V exp 2/3), where V is the flux tube volume. Starting with such a distribution function and a plasma beta which varies from beta greater than 1 in the CPS to beta much smaller than 1 in the lobe, the evolution of the ion distribution function was studied considering the combined effects of ion diffusion by kinetic Alfven waves (KAW) in the ULF frequency range (1-10 mHz) and convection due to E(dd) x B drift in the plasma sheet boundary layer (PSBL) and outer central plasma sheet (OCPS). The results show that, during the early stages after launching the KAWs, a beamlike ion distribution forms in the PSBL and at the same time the plasma density and temperature decrease in the OCPS. Following this stage, ions in the beams convect toward the CPS resulting in an increase of the plasma temperature in the OCPS.

Moghaddam-Taaheri, E.↗

NASA NTRS2010

Initial Results of DC Electric Fields, Associated Plasma Drifts, Magnetic Fields, and Plasma Waves Observed on the C/NOFS Satellite

Initial results are presented from the Vector Electric Field Investigation (VEFI) on the Air Force Communication/Navigation Outage Forecasting System (C/NOFS) satellite, a mission designed to understand, model, and forecast the presence of equatorial ionospheric irregularities. The VEFI instrument includes a vector DC electric field detector, a fixed-bias Langmuir probe operating in the ion saturation regime, a flux gate magnetometer, an optical lightning detector, and associated electronics including a burst memory. Compared to data obtained during more active solar conditions, the ambient DC electric fields and their associated E x B drifts are variable and somewhat weak, typically < 1 mV/m. Although average drift directions show similarities to those previously reported, eastward/outward during day and westward/downward at night, this pattern varies significantly with longitude and is not always present. Daytime vertical drifts near the magnetic equator are largest after sunrise, with smaller average velocities after noon. Little or no pre-reversal enhancement in the vertical drift near sunset is observed, attributable to the solar minimum conditions creating a much reduced neutral dynamo at the satellite altitude. The nighttime ionosphere is characterized by larger amplitude, structured electric fields, even where the plasma density appears nearly quiescent. Data from successive orbits reveal that the vertical drifts and plasma density are both clearly organized with longitude. The spread-F density depletions and corresponding electric fields that have been detected thus far have displayed a preponderance to appear between midnight and dawn. Associated with the narrow plasma depletions that are detected are broad spectra of electric field and plasma density irregularities for which a full vector set of measurements is available for detailed study. Finally, the data set includes a wide range of ELF/VLF/HF oscillations corresponding to a variety of plasma waves, in particular banded ELF hiss, whistlers, and lower hybrid wave turbulence triggered by lightning-induced sferics. The VEFI data represents a new set of measurements that are germane to numerous fundamental aspects of the electrodynamics and irregularities inherent to the Earth's low latitude ionosphere.

Pfaff, R.↗

DOE OSTI2021

Scattering of radiofrequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radiofrequency (RF) electromagnetic waves that propagate through this region. It is important to know the spectral properties of these scattered RF waves, whether used for diagnostics or for heating and current drive. The spectral changes influence the interpretation of the obtained diagnostic data, and the current and heating profiles. A full-wave, three-dimensional (3-D) electromagnetic code ScaRF has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method used for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma by the anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can also be used for the study of the scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U and DIII-D. For the case of medium-sized tokamaks, there is experimental evidence that drift waves and rippling modes are present in the edge region. Hence, we have studied the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes. The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion method in conjunction with the Smolyak sparse-grid integration, which is known as the PCE-SG method. The PCE-SG method is proven to be accurate and more efficient compared with alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Key Observations and Discoveries from the Vector Electric Field Investigation on the C/NOFS Satellite

Launched in 2008 and operating for 7.5 years, the Air Force Communication /Navigation Outage Forecasting System (C/NOFS) satellite included the Vector Electric Field Investigation (VEFI) designed and built at NASA’s Goddard Space Flight Center. VEFI successfully met its objectives: 1) investigate the role of ambient electric fields in initiating nighttime density depletions and turbulence; 2) determine the quasi-DC electric fields associated with abrupt density depletions, and 3) quantify the spectrum of the irregularities associated with density depletions, providing many key observations and discoveries including: -- Global (low latitude) vector DC electric fields at 16 s/sec revealing variations with local time, longitude, and season between extremely low and moderate solar activity-- Reversed zonal E x B drifts below the F-region ledge at sunset and simultaneous observations of large scale Kelvin-Helmholtz instabilities as seeds of spread-F -- Large expanses of quasi-coherent kilometer-scale vector wave observations (electric field and density) below 450 km and their discovery as a source of scintillations using the C/NOFS GPS -- Reversed zonal DC electric fields and simultaneous observations of afternoon counter electrojet-- Enhanced zonal DC electric fields at sunrise and full vector plasma flow continuity at the terminator-- Vector electric field and density irregularities extending to meter-scales within equatorial plasma depletions -- Intense electric field structure within the equatorial ionosphere at night in the absence of density depletions -- First observations of the westward equatorial electrojet in post-midnight ionosphere and possible association with downward meridional drifts-- Spaceborne Dst observations and implications for ring current asymmetries-- Vector electric and magnetic field Poynting flux within depletions and TIDs-- Measurements of ionospheric reflectance and Alfvenic waves associated with TIDs-- Observations of Alfven resonators -- ULF magnetic field structure within the nightime equatorial ionosphere-- Discovery of Schumman resonances in space and implications of ELF radiowave propagation-- Vector observations of 50-60 Hz powerline radiation without harmonics-- Ion cyclotron resonance absorption lines associated with ELF hiss and the identification of ambient ions-- Parallel electric fields associated with lower hybrid waves driven by thunderstorm lightning-related sferics-- Discovery of Z-mode radiation in the equatorial ionosphere including its associated with density depletions-- Optical lightning detector waveforms and electric field sferics observed up to altitudes of 800 km

Robert F Pfaff↗

NASA NTRS1984

Interferometer applications of VHF radars, part 3.6A

Using a spaced antenna setup of a VHF radar, the spatial distribution of amplitudes and phases of the radar echoes from the troposphere, stratosphere and mesosphere can be measured. Combining in a suitable analysis procedure the complex digital samples from the different receiving antennas is consistent with the radar interferometer method. In addition to the well-known parameters measured with the commonly applied Doppler and drifts methods, i.e., reflectivity and mean fluctuation velocity, the interferometer technique allows to measure the regular spectrum of the returns. This technique, which was first applied with the spaced antenna system of the SOUSY-VHF-Radar in W. Germany, as well as some first examples of results are described here. These comprise the measurements of the horizontal and vertical velocities of the mean flow as well as of turbulence structures, the aspect sensitivity and the tilt of layers from which the baroclinicity can be estimated. Particularly, results of interferometer measurements of the vertical and horizontal phase velocities and wavelengths of gravity waves in the stratosphere are displayed. The latter results are also discussed in terms of the generation and propagation of these waves.

Rottger, J.↗

NASA NTRS2016

Hall2De Simulations with a First-Principles Electron Transport Model Based on the Electron Cyclotron Drift Instability

Several years of work have combined plasma measurements in a laboratory Hall thruster and r-z numerical simulations with Hall2De to isolate the spatial variation of the anomalous collision frequency needed in Ohm’s law to produce the observed thruster behavior. This numerical solution is used here to test the validity of a first-principles model of the anomalous transport in these devices before such model is implemented self-consistently in r-z fluid codes like Hall2De. The first-principles model employs quasi-linear theory and is based on the hypothesis that the Electron Cyclotron Drift Instability (ECDI) excites ion acoustic turbulence that, in turn, enhances the effective collision frequency in these devices. We find that a model of the ECDI that assumes Maxwellian velocity distributions for electrons and singly-charged, main-beam, cold ions (Ti=0.07 eV) is insufficient to explain the expected variation of the anomalous collision frequency both in the interior and exterior of the acceleration channel. When warm ions (~0.5-3 eV) are accounted for, the ECDI model in the channel interior appears more promising but fails by orders of magnitude in the near plume region due to the much higher Landau damping of the ion acoustic waves there. This implies that either (a) some process allows the ECDI instability to remain uninhibited by classical Landau damping or, (b) that a different instability (or instabilities) altogether, also insusceptible to Landau damping, is/are active in this region. A previous hypothesis, that convection of wave energy generated by the ECDI in the channel plays a significant role in the near plume, is not supported by the results of the simulations.

Mikellides, Ioannis G.↗

DOE OSTI2021

Electromagnetic turbulence in increased β plasmas in the Large Plasma Device

The variation of pressure-gradient-driven turbulence with plasma $\beta$ (up to $\beta \approx 15\,\%$ ) is investigated in linear, magnetized plasma. The magnitude of magnetic fluctuations is observed to increase substantially with increasing $\beta$ . More importantly, parallel magnetic fluctuations are observed to dominate at higher $\beta$ values, with $\delta B_\parallel / \delta B_\perp \approx 2$ and $\delta B / B_0 \approx 1\,\%$ . Parallel magnetic fluctuations are strongly correlated with density fluctuations and the two are observed to be out of phase. The relative magnitude of and cross-phase between density and parallel magnetic field fluctuations are consistent with the dynamic pressure balance ( $P+{B_{0}^2}/{2\mu _0} = \textrm {constant}$ ). A local slab model theory for electromagnetic, modified drift Alfvén waves, including parallel magnetic fluctuations, shows partial agreement with experimental observations.

Physics↗

DOE OSTI2020

Simulation of divertor heat flux width on EAST by BOUT++ transport code

The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Zonal flow excitation in electron-scale tokamak turbulence

Abstract The derivation of an intermediate-scale gyrokinetic-electron theory in nonuniform tokamak plasmas (Chen et al 2021 Nucl. Fusion 61 066017 ) has shown that a Navier–Stokes type nonlinearity couples electron-temperature-gradient (ETG) modes and zonal flow (ZF) modes with wavelengths much shorter than the ion gyroradius but much longer than the electron gyroradius. This intermediate-scale ETG-ZF coupling is typically stronger than the Hasegawa–Mima type nonlinearity characteristic of the fluid approximation and is predicted to lead to relevant ZF generation and ETG mode regulation. Electron-scale, continuum, gyrokinetic simulation results are presented here which include both single-mode ETG and full-spectrum ETG turbulence. The ZF generation due to single ETG modes is investigated and the single-mode intermediate-scale results are found to be in agreement with theory. The full-spectrum results are then presented and explained qualitatively in terms of the single-mode results. It is found that the ETG-driven ZFs regulate intermediate-scale electron heat flux transport to levels in the predicted range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The Impact of Anomalous Resistivity in Vacuum Contaminant Plasmas on the Electrothermal Instability

This manuscript presents an assessment of the electrothermal instability (ETI) in the presence of anomalous resistivity (AR) in vacuum contaminant plasmas (VCP) when applied to a magnetized liner inertial fusion (MagLIF)-like load.Pulsed-power driven dielectrically coated metallic liners, like in MagLIF, experience the current-driven electrothermal instability which occurs when a material’s resistivity changes with temperature and is subject to ohmic heating. Large scale pulsed-power facilities that use magnetically insulated transmission lines (MITL) have been shown to generate low-density plasma which enters the target chamber and coalesces around the load. The low-density high-temperature vacuum contaminant plasmas (VCP) can parasitically divert current from the load through causing a short in the anode-cathode gap inside the target chamber. Resistive magnetohydrodynamic (MHD) simulations of these VCP experience unphysical runaway ohmic heating due to under predicting the resistivity by using a purely collisional resistivity model.AR provides a physics-based way to address this runaway heating through increasing the resistivity in a proportional way with the drift speed. In this work, 1D simulations probe the effect that AR in VCP has on the magnetic diffusion rate, and 2D simulations show how this effect manifests in the nonlinear striation form of the ETI for a MagLIF-like load. Beryllium and aluminum dielectrically coated liners are used for the 1D and 2D simulations in this work. The1D simulations show that a VCP causes a delay in the current delivery to the load by upwards of 8 ns at 60 ns into atypical current pulse for an 1 × 10 7 A scale pulsed-power accelerator. The 2D simulations show the delay observed in the 1D simulations is reduced substantially (4 ns) by 2D turbulence that disrupts the VCP layer early in time. Here, the 2Dsimulated ETI growth varies across AR models, more so for the beryllium liner than the aluminum liner because the beryllium liner shows an enhanced rate of penetration for the magnetic diffusion wave in comparison to the aluminum liner. The 2D simulations show the bulk dielectric thickness varies across AR models with the Davidson AR model beingthe largest and the Buneman AR model being the smallest, and in connection with the thickness the Rayleigh-Taylor bubble-spike distances varies correspondingly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2014

Visualizing Flutter Mechanism as Traveling Wave Through Animation of Simulation Results for the Semi-Span Super-Sonic Transport Wind-Tunnel Model

It has long been recognized that frequency and phasing of structural modes in the presence of airflow play a fundamental role in the occurrence of flutter. Animation of simulation results for the long, slender Semi-Span Super-Sonic Transport (S4T) wind-tunnel model demonstrates that, for the case of mass-ballasted nacelles, the flutter mode can be described as a traveling wave propagating downstream. Such a characterization provides certain insights, such as (1) describing the means by which energy is transferred from the airflow to the structure, (2) identifying airspeed as an upper limit for speed of wave propagation, (3) providing an interpretation for a companion mode that coalesces in frequency with the flutter mode but becomes very well damped, (4) providing an explanation for bursts of response to uniform turbulence, and (5) providing an explanation for loss of low frequency (lead) phase margin with increases in dynamic pressure (at constant Mach number) for feedback systems that use sensors located upstream from active control surfaces. Results from simulation animation, simplified modeling, and wind-tunnel testing are presented for comparison. The simulation animation was generated using double time-integration in Simulink of vertical accelerometer signals distributed over wing and fuselage, along with time histories for actuated control surfaces. Crossing points for a zero-elevation reference plane were tracked along a network of lines connecting the accelerometer locations. Accelerometer signals were used in preference to modal displacement state variables in anticipation that the technique could be used to animate motion of the actual wind-tunnel model using data acquired during testing. Double integration of wind-tunnel accelerometer signals introduced severe drift even with removal of both position and rate biases such that the technique does not currently work. Using wind-tunnel data to drive a Kalman filter based upon fitting coefficients to analytical mode shapes might provide a better means to animate the wind tunnel data.

Christhilf, David M.↗

NASA NTRS2009

Initial Results from the Vector Electric Field Investigation on the C/NOFS Satellite

Initial results are presented from the Vector Electric Field Investigation (VEFI) on the Air Force Communication/Navigation Outage Forecasting System (C/NOFS) satellite, a mission designed to understand, model, and forecast the presence of equatorial ionospheric irregularities. The VEFI instrument includes a vector DC electric field detector, a fixed-bias Langmuir probe operating in the ion saturation regime, a flux gate magnetometer, an optical lightning detector, and associated electronics including a burst memory. The DC electric field detector has revealed zonal and meridional electric fields that undergo a diurnal variation, typically displaying eastward and outward-directed fields during the day and westward and downward-directed fields at night. In general, the measured DC electric field amplitudes are in the 0.5-2 mV/m range, corresponding to I3 x B drifts of the order of 30-150 m/s. What is surprising is the high degree of large-scale (10's to 100's of km) structure in the DC electric field, particularly at night, regardless of whether well-defined spread-F plasma density depletions are present. The spread-F density depletions and corresponding electric fields that have been detected thus far have displayed a preponderance to appear between midnight and dawn. Associated with the narrow plasma depletions that are detected are broad spectra of electric field and plasma density irregularities for which a full vector set of measurements is available for detailed study. On some occasions, localized regions of low frequency (< 8 Hz) magnetic field broadband irregularities have been detected, suggestive of filamentary currents, although there is no one-to-one correspondence of these waves with the observed plasma density depletions, at least within the data examined thus far. Finally, the data set includes a wide range of ELF/VLF/HF waves corresponding to a variety of plasma waves, in particular banded ELF hiss, whistlers, and lower hybrid wave turbulence triggered by lightning-induced sferics. The VEFI data set represents a treasure trove of measurements that are germane to numerous fundamental aspects of the electrodynamics and irregularities inherent to the Earth's low latitude ionosphere.

Pfaff, R.↗

DOE OSTI2021

Simulations of heat fluxes in an ELMy H-mode discharge on HL-2A

In order to study the distribution and evolution of the transient heat flux on HL-2A during edge-localized-mode (ELM) bursts, the BOUT++ electromagnetic six-field two-fluid model is used to simulate the pedestal collapse under the lower single-null divertor geometry. The equilibrium profiles of HL-2A ELMy H-mode discharge No. 24 953 are adopted as the initial condition in the original case. In this instance, linear analysis shows that the resistive ballooning mode (RBM) and drift-Alfven wave are unstable to this equilibrium, and RBM is the dominant instability. The evolutions of the radial heat fluxes at the outer mid-plane and heat fluxes to the inner and outer targets during the ELM event are presented. Six more equilibria are constructed based on the original case to find out the influence of the pedestal profiles on the peak electron heat flux. The results indicate that the heat flux increases with temperature and/or density, and the theoretical analysis and simulation results consistently show that the heat flux q ∥e is proportional to n e0,SEP $T^{\frac32}_{e0,SEP}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Heat pulse propagation and anomalous electron heat transport measurements on the optimized stellarator W7-X

The optimized stellarator Wendelstein 7-X (W7-X) is designed to have an approximately quasi-isodynamic magnetic configuration with reduced neoclassical transport in comparison to a classical stellarator, and turbulent transport is expected to be a significant source of anomalous heat transport across the plasma minor radius. The ion temperature gradient driven mode and the trapped electron mode (TEM) are thought to be responsible for the ion-scale turbulence in W7-X plasmas with volume averaged pressure below 1%. In this work, the electron temperature gradient driven turbulence is shown to be a good candidate for the explanation of the observed electron heat flux, in the inner plasma region where the density gradient is weak (in the outer region, a relatively stronger density gradient would drive additional TEM turbulence). Further, the experimental electron heat transport measured during electron cyclotron resonant heating power and plasma density scans is compared to neoclassical predictions, and the stiffness in the electron heat transport measured during transient transport experiments is presented in three common magnetic configurations of W7-X. In low-(β) plasma discharges, the stiffness in the electron heat flux, quantified by the ratio of the heat pulse to power balance diffusivity, ${\chi }_{\mathrm{e}}^{\mathrm{H}\mathrm{P}}/{\chi }_{\mathrm{e}}^{\mathrm{P}\mathrm{B}}$, is measured to be less than 2, and trend downwards with increasing collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Acoustic instability driven by cosmic-ray streaming

We study the linear stability of compressional waves in a medium through which cosmic rays stream at the Alfven speed due to strong coupling with Alfven waves. Acoustic waves can be driven unstable by the cosmic-ray drift, provided that the streaming speed is sufficiently large compared to the thermal sound speed. Two effects can cause instability: (1) the heating of the thermal gas due to the damping of Alfven waves driven unstable by cosmic-ray streaming; and (2) phase shifts in the cosmic-ray pressure perturbation caused by the combination of cosmic-ray streaming and diffusion. The instability does not depend on the magnitude of the background cosmic-ray pressure gradient, and occurs whether or not cosmic-ray diffusion is important relative to streaming. When the cosmic-ray pressure is small compared to the gas pressure, or cosmic-ray diffusion is strong, the instability manifests itself as a weak overstability of slow magnetosonic waves. Larger cosmic-ray pressure gives rise to new hybrid modes, which can be strongly unstable in the limits of both weak and strong cosmic-ray diffusion and in the presence of thermal conduction. Parts of our analysis parallel earlier work by McKenzie & Webb (which were brought to our attention after this paper was accepted for publication), but our treatment of diffusive effects, thermal conduction, and nonlinearities represent significant extensions. Although the linear growth rate of instability is independent of the background cosmic-ray pressure gradient, the onset of nonlinear eff ects does depend on absolute value of DEL (vector differential operator) P(sub c). At the onset of nonlinearity the fractional amplitude of cosmic-ray pressure perturbations is delta P(sub C)/P(sub C) approximately (kL) (exp -1) much less than 1, where k is the wavenumber and L is the pressure scale height of the unperturbed cosmic rays. We speculate that the instability may lead to a mode of cosmic-ray transport in which plateaus of uniform cosmic-ray pressure are separated by either laminar or turbulent jumps in which the thermal gas is subject to intense heating.

Begelman, Mitchell C.↗