Search NASA⌕ Search

SEARCH · Search NASA

Results for “resistive instabilities”

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.

At least 235 records · Page 13

Transition from retrograde to prograde drift instabilities in a magnetron microdischarge

Rotating plasma structures, or “spokes,” in magnetized discharges characterized by perpendicular electric and magnetic fields have been seen in a growing number of studies and are believed to be the result of gradient-driven drift instabilities. Under certain conditions, we have shown the spoke’s rotation to be opposite to the E B direction, i.e., retrograde in its expected direction. Recently, we have linked such counter-intuitive rotation to a local inversion of the electric field. Here, we give further experimental evidence for this inversion and attempt to provide an explanation for a relatively distinct transition seen between retrograde and positive (prograde) drift. In the experiments, a partially magnetized plasma forms inside a magnetron device of mm radius operated with argon. Discharge current–voltage measurements are acquired for a range of argon fill pressure and inter-electrode spacing. We find two branches of operation—a low current branch of negative resistance, coinciding with the retrograde spoke rotation, and a higher current branch of positive resistance, coincident with prograde spoke rotation. We postulate that at low discharge currents, high magnetic field confinement leads to a large density gradient, causing more electron transport to the anode than that demanded by the external circuit. At higher currents, anomalous axial electron transport (across the magnetic field lines) becomes dominant, eliminating the conditions for field inversion. Here, the current thresholds for the field inversion are found to be sensitive to the argon pressure and inter-electrode spacing.

42 ENGINEERING↗

Simulations of saturated MHD activity in the HBT-EP tokamak

The NIMROD code is used to perform simulations of MHD activity in the HBT-EP tokamak, including the effect of a resistive wall. Linear simulations are conducted to scan current and pressure profiles to establish self-consistent baseline equilibria that result in MHD instabilities with low error between experimental diagnostics and synthetic equivalents. A family of such equilibria is found with relatively peaked current and pressure profiles. Nonlinear simulations are performed using the optimal equilibrium, which is linearly unstable to a Resistive Wall Mode—stable in the ideal wall limit and with an Alfvénic growth rate in the no-wall limit. Using resistive wall boundary conditions, nonlinear Hall MHD simulations yield a rotating, saturated n=1 mode similar to the type observed experimentally in HBT-EP discharges. An n=1 mode around the q=2 surface mediates a localized axisymmetric perturbation that flattens the plasma current from a linearly unstable state to a 3D, rotating, stable state. During the saturated period, periodic changes in inter-diagnostic phase lag allow for qualitative estimations of effective, macroscopic transport coefficients used in MHD closure models and could be extended to model post-disruption MHD structure rotation after thermal quenches seen in HBT-EP.

Arnold, David A. (ORCID:0009000059535761)↗

E region coupling effects on the Perkins spread F instability

The Perkins instability theory for mid-latitude spread F is generalized to include the effects of the E region on growth rates. A three-layer model, representing the E region, the F region, and a resistive layer of low Pedersen conductivity which separates them, is developed. The E region Pedersen current has a stabilizing effect on the perturbed F region, thereby reducing the growth rate of the instability. For spread F onset on September 17, 1974, the maximum growth rate is found to be only one third of the Perkins growth rate. The resistive layer reduces the electric field perturbation in the E region for short wavelengths but has little effect for long wavelengths.

Klevans, E. H.↗

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes

The Kelvin–Helmholtz (KH) instability of a shear layer with an initially uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with hydrodynamic KH, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited, thus transporting less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear because of the stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Local structure effects of carbon-doping on the phase change material Ge 2 Sb 2 Te 5

Ge 2 Sb 2 Te 5 is used in phase change memory, a nonvolatile memory technology, due to its phase change properties. The primary advantage of phase change memory over the state-of-the-art (flash memory) is its simple and small device geometry, which allows for denser nodes and lower power consumption. In phase change memory, resistive heating induces fast switching between the high resistance amorphous and low resistance crystalline phases, corresponding to storage of low and high digital states, respectively. However, the instability of the amorphous phase of Ge 2 Sb 2 Te 5 presents issues with processing and long-term data storage; such issues can be resolved by C doping, which stabilizes the amorphous phase and raises the crystallization temperature. To better understand the local structural effects of C doping on Ge 2 Sb 2 Te 5 , in situ Ge K-edge X-ray absorption spectroscopy measurements were taken during heating of films with various C doping concentrations. Here, the range of structural transformation temperatures derived from X-ray absorption near-edge structure analysis across the C doping series proved narrower than crystallization temperatures reported in similar in situ X-ray diffraction experiments, which may reflect changes in local structure that precede long-range ordering during crystallization. In addition, rigorous extended X-ray absorption fine structure fitting across and between temperature series revealed effects of C doping on the rigidity of Ge–Te bonds at low (2 at% and 4 at%) C concentrations.

Langhout, John D.↗

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Instability and turbulent relaxation in a stochastic magnetic field

An analysis of instability dynamics in a stochastic magnetic field is presented for the tractable case of the resistive interchange. Externally prescribed static magnetic perturbations convert the eigenmode problem to a stochastic differential equation, which is solved by the method of averaging. The dynamics are rendered multi-scale, due to the size disparity between the test mode and magnetic perturbations. Maintaining quasi-neutrality at all orders requires that small-scale convective cell turbulence be driven by disparate scale interaction. Here, the cells in turn produce turbulent mixing of vorticity and pressure, which is calculated by fluctuation-dissipation type analyses, and are relevant to pump-out phenomena. The development of correlation between the ambient magnetic perturbations and the cells is demonstrated, showing that turbulence will 'lock on' to ambient stochasticity. Magnetic perturbations are shown to produce a magnetic braking effect on vorticity generation at large scale. Detailed testable predictions are presented. The relations of these findings to the results of available simulations and recent experiments are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lattice instabilities and phonon thermal transport in TlBr

Materials with simple crystal structure usually exhibit substantial lattice thermal conductivity ( κ latt ) due to a scarcity of strong phonon-phonon scattering channels. Yet here we find simple CsCl-structure TlBr to have an extremely low room temperature κ latt ~ 0.5 W / m-K , while at low temperature we find a high κ latt of 155 W / m-K . This extreme range of conductivity behavior results from lattice-instability-related anharmonicity, effective only at temperatures above ~ 30 K . Manipulating these instabilities (e.g., pressure, strain) may yield a pathway to ultralow thermal conductivity via intrinsic phonon resistance in similar materials.

36 MATERIALS SCIENCE↗

On the role of the lower hybrid drift instability in substorm dynamics

Recent studies of the lower hybrid drift instability have shed new light on the role of this mode in field-reversed plasmas. For substorm magnetotail conditions it is found that the lower hybrid drift instability can penetrate to the neutral line and can dissipate magnetic energy at a rate of approximately 4 x 10 to the 17th erg/s. Thus this instability is capable of playing a major role in the onset of substorms and providing resistivity for reconnection processes in the context of the neutral line substorm model.

Huba, J. D.↗

Characteristics of Si3N4-SiO2-Ce2O3 compositions sintered in high-pressure nitrogen

Full-density Si3N4-SiO2-Ce2O3 compositions were prepared by sintering with 2.5 MPa nitrogen pressure at temperatures of 1900 and 2090 C. Room-temperature flexural strengths near 700 MPa for sintered material compared favorably with the strength of hot-pressed material. At 1370 C, where flexural strengths as high as 363 MPa were obtained, it was observed that the coarsest structure was the strongest and the finest structure was the weakest. One of the compositions tested, Si3N4-8.7 wt pct SiO2-8.3 wt pct-Ce2O3, was found to have excellent 200-h oxidation resistance at 700, 1000, and 1370 C, without incidence of 700 to 1000 C phase instability and cracking.

Sanders, W. A.↗

The magnetohydrodynamics of current sheets

Examples of current sheets are summarized and their formation is described. A universal phenomenon in cosmic plasmas is the creation of sheets off intense current near X-type neutral points (where the magnetic field vanishes). These sheets are important as sites where the magnetic-field energy is converted efficiently into heat and bulk kinetic energy and where particles can be accelerated to high energies. Examples include disruptions in laboratory tokamaks, substorms in the earth's magnetosphere, and flares on the sun. The basic behavior of a one-dimensional sheet is presented, together with an account of the linear tearing-mode instability that can cause the field lines in such a sheet to reconnect. Such reconnection may develop in different ways: it may arise from a spontaneous instability or it may be driven, either from outside by motions or locally by a resistivity enhancement. Various processes are described that may occur during the nonlinear development of tearing, along with the many numerical and laboratory experiments that are aiding our understanding of this intriguing cosmical process.

Priest, E. R.↗

Numerical simulation of MPD thruster flows with anomalous transport

Anomalous transport effects in an Ar self-field coaxial MPD thruster are presently studied by means of a fully 2D two-fluid numerical code; its calculations are extended to a range of typical operating conditions. An effort is made to compare the spatial distribution of the steady state flow and field properties and thruster power-dissipation values for simulation runs with and without anomalous transport. A conductivity law based on the nonlinear saturation of lower hybrid current-driven instability is used for the calculations. Anomalous-transport simulation runs have indicated that the resistivity in specific areas of the discharge is significantly higher than that calculated in classical runs.

Caldo, Giuliano↗

Stability of the Halley cometosheath with resistivity and plasma motion

The MHD stability of the cometary inner sheath determined by the balance between the inward Lorentz body force and the outward ion-neutral drag force is investigated by numerically solving the wave equations which include resistivity, plasma motion and plasma pressure with the help of two-point boundary value method. The eigenvalues and the eigenfunctions are obtained numerically by treating the cometary inner sheath as a layer of finite thickness, bounded by the contact surface, that is, the diamagnetic cavity boundary. To gain insight into the problem, certain limiting cases of the wave equations are also discussed. The diamagnetic cavity boundary and the adjacent layer of about 100-km thickness of Comet Halley is found to be unstable. The effects of finite plasma pressure, dissociative recombination, mass loading due to photoionization, resistivity, and plasma motion are found to be stabilizing but are unable to quench the instability completely. Motion of the Halley ionopause has been confirmed by observations.

Srivastava, Krishna M.↗

Composite Fuselage Technology

Work was conducted over a ten-year period to address the central issue of damage in primary load-bearing aircraft composite structure, specifically fuselage structure. This included the three facets of damage resistance, damage tolerance, and damage arrest. Experimental, analytical, and numerical work was conducted in order to identify and better understand the mechanisms that control the structural behavior of fuselage structures in their response to the three aspects of damage. Furthermore, work was done to develop straightforward design methodologies that can be employed by structural designers in preliminary design stages to make intelligent choices concerning the material, layup, and structural configurations so that a more efficient structure with structural integrity can be designed and built. Considerable progress was made towards achieving these goals via this work. In regard to damage tolerance considerations, the following were identified as important effects: composite layup and associated orthotropy/structural anisotropy, specifics of initial local damage mechanisms, role of longitudinal versus hoop stress, and large deformation and associated geometric nonlinearity. Means were established to account for effects of radius and for the nonlinear response. In particular, nondimensional parameters were identified to characterize the importance of nonlinearity in the response of pressurized cylinders. This led to the establishment of a iso-nonlinear-error plot for reference in structural design. Finally, in the case of damage tolerance, the general approach of the original methodology to predict the failure pressure involving extending basic plate failure data by accounting for the local stress intensification was accomplished for the general case by accounting for the mechanisms noted by utilizing the capability of the STAGS finite element code and numerically calculating the local stress intensification for the particular configuration to be considered. For the issue of damage arrest, placement of and configuration of stiffeners (including stiffener curvature), and magnitude and orientation of principal strains due to local bending were found to be key considerations. Means were established to account for stiffener effectiveness quantitatively based on radius, slit size, stiffener curvature' and relative bending stifffiesses involved. Geometric nonlinearity was also found to play an - 24 - important role here. Furthermore, it was determined that damage propagation is controlled by different mechanisms (hoop stress versus flapping stress and the associated factors involved in each) depending upon the direction of damage propagation. This latter item results in an inability to scale these phenomena in one test due to the different factors involved. Finally, the importance of shell curvature and associated instability in response to transverse loading including impact were found to be important considerations in damage resistance. A technique, involving asymmetric meshing of a finite element mesh, was developed to predict this behavior and showed excellent correlation with experimental results. Further details of these ten years of work are presented herein with references made to the fourteen documents produced during this work where full details can be found. Implications of this work are discussed and recommendations made. Although it is clear that there is more work to be done to fully understand composite fuselage technology and specifically the overall issue of damage in primary load-bearing composite structures, important understanding and capability has been extended via this work.

Lagace, Paul A.↗

Successful Surface Treatments for Reducing Instabilities in Advanced Nickel-base Superalloys for Turbine Blades

An optimized carburization treatment has been developed to mitigate instabilities that form in the microstructures of advanced turbine airfoil materials. Current turbine airfoils consist of a single crystal superalloy base that provides the mechanical performance of the airfoil, a thermal barrier coating (TBC) that reduces the temperature of the base superalloy, and a bondcoat between the superalloy and the TBC, that improves the oxidation and corrosion resistance of the base superalloy and the spallation resistance of the TBC. Advanced nickel-base superalloys containing high levels of refractory metals have been observed to develop an instability called secondary reaction zone (SRZ), which can form beneath diffusion aluminide bondcoats. This instability between the superalloy and the bondcoat has the potential of reducing the mechanical properties of thin-wall turbine airfoils. Controlled gas carburization treatments combined with a prior stress relief heat treatment and adequate surface preparation have been utilized effectively to minimize the formation of SRZ. These additional processing steps are employed before the aluminide bondcoat is deposited and are believed to change the local chemistry and local stresses of the surface of the superalloy. This paper presents the detailed processing steps used to reduce SRZ between platinum aluminide bondcoats and advanced single crystal superalloys.

Locci, Ivan E.↗

Generation of lower-hybrid noise by superthermal cross-field ion currents

A kinetic model for generating lower-hybrid waves by a low-density, superthermal ion beam directed across a magnetic field is investigated. Growth rate calculations demonstrate the salient features of the process, in particular, the small threshold drift of several background ion thermal speeds required for instability. The mechanism is suitable for collisionless dissipation in systems which evince the characteristics of anomalous resistivity in the presence of relatively small electrical currents.

Barbosa, D. D.↗

Alfven shock trains

The Cohen-Kulsrud-Burgers equation (CKB) is used to consider the nonlinear evolution of resistive, quasi-parallel Alfven waves subject to a long-wavelength, plane-polarized, monochromatic instability. The instability saturates by nonlinear steepening, which proceeds until the periodic waveform develops an interior scale length comparable to the dissipation length; a fast or an intermediate shock then forms. The result is a periodic train of Alfven shocks of one or the other type. For propagation strictly parallel to the magnetic field, there will be two shocks per instability wavelength. Numerical integration of the time-dependent CKB equation shows that an initial, small-amplitude growing wave asymptotes to a stable, periodic stationary wave whose analytic solution specifies how the type of shock embedded in the shock train, and the amplitude and speed of the shock train, depend on the strength and phase of the instability. Waveforms observed upstream of the earth's bowshock and cometary shocks resemble those calculated here.

Malkov, M. A.↗