The Effect of Anomalous Resistivity on Electrothermal Instability Growth
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Beams of energetic runaway electrons are generated during disruptions in tokamaks, and fluid models are used to study their effects on macroscale dynamics. Linear computations of a massless, runaway electron beam coupled to MHD plasma show that resistive hose instabilities grow faster than tearing modes at large resistivity. Eigenvalue results with reduced models of the resistive hose instability are compared with results from the full MHD and beam system, showing that the resistive hose decouples from any plasma response. An estimate of plasma temperature at which growth of the resistive hose dominates tearing for post-disruption DIII-D plasma parameters is in a physically relevant regime.
Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.
A spectral numerical-simulation code and ancillary diagnostics are used here to analyze the behavior of resistive magnetic tearing in the presence of short-wavelength turbulence. It is found that, in general, the anomalous resistivity defined by Biskamp and Welter (1983) due to short-wavelength turbulence tends to return quickly towards zero even when artificially supported away from zero, and that its effect on tearing-mode behavior is not consistent with its interpretation as a resistivity. The behavior reported by Biskamp and Welter and the behavior observed here are analyzed. It is argued that, while not meaningful as a true resistivity, the anomalous resistivity parameter is a useful diagnostic showing the energy balance of the system.
The toroidal distribution of runaway electrons (REs) striking the centre post in DIII-D post-disruption RE final loss events is simulated utilizing the MARS-F code combined with the REORBIT module, by tracing the guiding-centre drift orbits of test REs in the presence of external (applied) plus internal (intrinsic instability) 3-D fields. To better recover the experimental results, three different equilibria with different safety factor profiles are adopted. Dominant resistive magnetohydrodynamic instabilities are found to be m/n=2/1, 3/2 and 1/1 modes (where m and n are the poloidal and toroidal mode numbers). Externally applied n=1 resonant magnetic perturbations (RMPs), with both even and odd parity configuration, are considered in each equilibrium. Locking of MHD instability to the applied RMPs is simulated by minimizing the total perturbed magnetic energy in the plasma. Here, the simulated toroidal locations, where most REs strike in the presence of an edge m/n= 2/1 instability and the n=1 RMPs in even parity, are found to agree well with experiments. Modeling also captures the measured n = 2 RE loss pattern on the wall for experiments where an edge m/n=3/2 instability is dominant instead of the m/n=2/1 mode. The key overall findings include (a) the measured global RE toroidal impact pattern in experiments can be explained by a locked resistive kink instability, (b) the toroidal phase of the locked mode is well predicted by minimizing the total perturbed magnetic energy, and (c) this edge locked mode determines the toroidal phase of the peak RE wall impact.
Previous stability analysis of NSTX equilibrium with lithium-conditioning demonstrates that enhanced resistivity due to the increased effective charge number Z eff (i.e. increased impurity level) can provide a stabilizing effect on low-n edge localized modes. This paper extends the resistivity stabilizing effect to the intermediate-n peeling-ballooning (PB) instabilities with the linear stability analysis of EAST high-confinement mode equilibria in NIMROD two-fluid calculations. However, in this work, the resistivity stabilizing effect on PB instabilities in the EAST tokamak appears weaker than that found in NSTX. This work may provide better insight into the physical mechanism behind the beneficial effects of impurities on pedestal stability.
The multi-pulse test line (MPTL) presently under T development at the Dual Axis Radiography Hydrodynamic Test (DARHT) facility may have several meters of drifting beam transport without external magnetic focusing field. The absence of external magnetic focusing is the worst case for growth of the resistive-wall instability. This instability is usually considered to be strictly a long-pulse (low-frequency) problem, but has been theoretically shown to couple the multiple pulses in a long pulse train, such as are expected to be the subject of experimentation on the MPTL. In this note, we explore the resistive wall instability in the MPTL parameter range using our LAMDA beam dynamics code. (LAMDA stands for Linear Accelerator Model for DARHT)
Flute, whistler, loss-cone, and resistive plasma instabilities in magnetosphere
Explore the source record for details and available documents.
The effects of velocity shear on the resistive drift wave instability in the non-adiabatic limit and the toroidal ion temperature gradient instability are investigated for a plasma of inhomogeneous density/ion temperature, respectively. For the resistive drift wave, we find that the instability growth rate decreases monotonically with increasing magnitude of shear, but we find that complete stabilization is impossible. For the ion temperature gradient instability, we find that the standard WKB approximation is insufficient to describe the full behavior of the instability, and that an analysis of the localized eigenmode problem reveals two separate unstable solutions, which the WKB approximation does not predict. The impact of flow shear on these two new unstable solutions is discussed. In both resistive drift wave and ion temperature gradient instabilities, the sheared flow causes a shifting, tilting, and sharpening of the electrostatic potential eddies.
The paper analyzes the 3D resistive tearing instability in a magnetotail configuration by means of a resistive MHD code for various forms of the energy equation, including ohmic and compressional heating. A complete derivation of the set of MHD equations from a kinetic formalism is presented, with all simplifying assumptions noted. Emphasis is placed on the changes in magnetotail evolution caused by a variation of the assumed value of the polytropic exponent gamma, representing the ratio of specific heat in the energy equation which closes the set of MHD equations. Within the gamma range of 2/3-4, a tearing mode is found to develop in all cases. Certain properties, such as the formation and ejection of a plasmoid, associated with the formation of the reversal region of the north-south component of the magnetic field, and strong earthward and tailward flow fields prevail as qualitative features independent of the value of gamma chosen. Ohmic dissipation is found to be the major source of the heating, while adiabatic or nonadiabatic compression may contribute primarily in the closed field region earthward of the reconnection site.
Abstract Recent MHD simulations find that the reversal of the local resistive interchange parameter D R from negative to positive due to impurity radiation cooling is able to linearly destabilize the resistive tearing mode (TM) growth in a tokamak above a threshold in impurity level. A layer of perturbed Pfirsch–Schlüter current density and resistivity are also induced by the impurity radiation, which further govern the TM growth and saturation in the nonlinear stage. The impurity threshold and the TM growth strongly depend on the parallel thermal conductivity, and such a dependence derives from the impact on D R of the fast parallel thermal equilibration along the helical magnetic field lines.
We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.
The stability of the geomagnetic tail is investigated on the basis of three dimensional resistive magnetohydrodynamic simulations, using different dynamic constraints and different initial equilibria. Different forms of the energy equation for isotropic pressure are found to have no significant effect on the dynamic growth of a resistive tearing instability, which is responsible for near Earth reconnection, plasmoid formation and ejection, and the generation of fast plasma flows. The constraints of a modified double adiabatic approach, however, can quench the tearing instability through the development of large, mirror type, anisotropies in the boundary regions of the plasma sheet, unless isotropization occurs on fast, nearly Alfvenic, time scales. The presence of a net cross tail magnetic field component B(sub yN) can reduce the growth of the instability without complete stabilization. An increase of B(sub z) from midnight toward the tail flanks, however, by more than a factor of about 3, apparently completely stabilizes the tearing mode. Stabilization and destabilization thus may depend on properties and constraints (and their release) in regions other than the neutral sheet where reconnection is initiated.
Combustion instability, where unsteady heat release couples with acoustic modes, has long been an area of concern in liquid rocket engines. Accurate modeling of the acoustic normal modes of the combustion chamber is important to understanding and preventing combustion instability. The injector resistance can have a significant influence on the chamber normal mode shape, and hence on the system stability.
The new Scorpius linear induction electron accelerator is under development for multi-pulse flash radiography of large, explosively-driven hydrodynamic experiments. Beam physics from the cathode to the target was examined with computer simulations, including envelope, ray-trace, and particle-in cell (PIC) codes. Beam instabilities investigated included beam breakup (BBU), image displacement, diocotron, parametric envelope, ion hose, and the resistive wall instability. Beam corkscrew motion and emittance growth from beam mismatch were also studied. The results of these simulations is documented in a series of reports. In this report the computer codes and physical models used for these simulations are described. The conclusion of this study is that Scorpius will produce and accelerate a beam with radiographic quality equivalent to the present accelerators at Los Alamos National Laboratory if the same engineering standards and construction details are upheld.
The linear theory of the generation of electrostatic bursts of noise by electrons trapped in chorus wave packets is developed for a finite temperature electron beam and a Maxwellian elecron and ion background. The growth rates determined qualitatively in good agreement with those obtained by previous authors from a more idealized model. Two connected instability mechanisms seem to be occurring: a beam plasma (electron-ion two-stream) instability commonly associated with intensification of the chorus power levels, and a transitional or borderline resistive medium instability commonly associated with chorus hooks. The physical reasons for the two mechanisms is discussed. In the second case electron beams are difficult to identify in the particle data. An expression is obtained for the maximum growth rate in terms of the ratios of the beam and electron thermal velocities to the beam velocity, and of the beam density to plasma density. It is anticipated that this may allow the observed peak in the electrostatic noise spectrum to be used as a diagnostic for the beam characteristics. Previously announced in STAR as N84-12832
A one-dimensional non-Maxwellian model for the steady neutral sheet at the center of the earth's plasma sheet is constructed, and its instability properties are examined with respect to the lower-hybrid-drift wave. It is shown that in the limit of a thin neutral sheet, the distribution of modified Alper and Boltzmann distributions is a more powerful source of free energy than the usual drifting Maxwellian. It is also shown that most regions within the neutral sheet are unstable, and the frequency spectrum of the unstable waves is nearly the same as in the drifting Maxwellian case. The growth rate is found to be enhanced, and the saturation level of the fluctuating field and the anomalous resistivity calculated in the entire sheet are found to be 3 to 4 times larger than those in a Maxwellian sheet.