Shear stabilization of a collision-dominated plasma.
Magnetic shear effect on stability of electron- neutral collision dominated plasma, noting helical perturbations onset
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.
Magnetic shear effect on stability of electron- neutral collision dominated plasma, noting helical perturbations onset
A general stability analysis is given of the Kevin-Helmholtz instability, for the case of sheared MHD flow of finite thickness in a compressible plasma which allows for the arbitrary orientation of the magnetic field, velocity flow, and wave vector in the plane perpendicular to the velocity gradient. The stability problem is reduced to the solution of a single second-order differential equation including a gravitational term to represent the coupling between the Kelvin-Helmholtz mode and the interchange mode. Compressibility and a magnetic field component parallel to the flow are found to be stabilizing effects, with destabilization of only the fast magnetosonic mode in the transverse case, and the presence of both Alfven and slow magnetosonic components in the parallel case. Analysis results are used in a discussion of the stability of sheared plasma flow at the magnetopause boundary and in the solar wind.
Magnetic field nonuniformity and its effect on stability of theta pinch plasma
The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.
A search for causes of intermittent mid-term (about two hours) instability in emission signals from an inductively coupled plasma led to adoption of a tandem spray-chamber arrangement and subsequently to use of a surfactant (Triton X-100) to mitigate the remaining and newly found instabilities. Through a series of investigations, abrupt signal excursions in the tandem setup were traced to droplet coagulation and drainage inside the glass tube that connected the two spray chambers. However, the signal shifts were not the result of sample-solution release directly but rather to the influence of the underlying factors on plasma behavior. Experiments tailored to the study included not only examination of temporal signal behavior but also collection of long-term videos and measurement of radiofrequency characteristics of the plasma. The addition of a surfactant, Triton X-100, for signal-stability improvement is applicable not only to systems that employ tandem spray chambers but also to conventional single Scott-type chamber arrangements. Further, use of the surfactant was unsuccessful in overcoming “acid effects”, either of the steady-state or transient nature, and did not alter plasma background or analyte signals significantly.
The temperature dependence of the scattering and absorption coefficients for a set of freestanding plasma-sprayed 8 wt% yttria-stabilized zirconia (8YSZ) thermal barrier coatings (TBCs) was determined at temperatures up to 1360 C in a wavelength range from 1.2 micrometers up to the 8YSZ absorption edge. The scattering and absorption coefficients were determined by fitting the directional-hemispherical reflectance and transmittance values calculated by a four-flux Kubelka Munk method to the experimentally measured hemispherical-directional reflectance and transmittance values obtained for five 8YSZ thicknesses. The scattering coefficient exhibited a continuous decrease with increasing wavelength and showed no significant temperature dependence. The scattering is primarily attributed to the relatively temperature-insensitive refractive index mismatch between the 8YSZ and its internal voids. The absorption coefficient was very low (less than 1 per centimeter) at wavelengths between 2 micrometers and the absorption edge and showed a definite temperature dependence that consisted of a shift of the absorption edge to shorter wavelengths and an increase in the weak absorption below the absorption edge with increasing temperature. The shift in the absorption edge with temperature is attributed to strongly temperature-dependent multiphonon absorption. While TBC hemispherical transmittance beyond the absorption edge can be predicted by a simple exponential decrease with thickness, below the absorption edge, typical TBC thicknesses are well below the thickness range where a simple exponential decrease in hemispherical transmittance with TBC thickness is expected. [Correction added after online publication August 11, 2009: "edge to a shorter wavelengths" has been updated as edge to shorter wavelengths."]
A theoretical treatment of the linear stability of sub-Alfvenic plasma expansion is developed. The theory is fully kinetic and includes finite-beta effects, collisional effects, and neutral gas flow. A variety of results are obtained, and are applied to the the AMPTE magnetotail release, the NRL laser experiment, and the upcoming CRRES GTO releases.
GE Steam Power, Inc. (GE) proposed to improve reliability, flexibility, and economics of an existing coal-fired power plant by applying a new advanced technology developed by GE, a plasma-assisted pulverized fuel firing system. The objective of this program is to demonstrate the achievement of lower load by improved flame stabilization and therefore lower operating costs in a full-scale field installation at coal-fired electric utility. GE’s Plasma Ignition and Combustion Stabilization System is designed to operate continuously to support low load operation. With the plasma on, the flame will be attached and stable, removing the firing system as a limitation to low load operation. In addition, GE’s exclusively from ABENZ company licensed AC based technology has a 90+% system efficiency compared to all other systems at which are DC and operate with ~75% efficiency. Plant operating costs are lowered by eliminating use of expensive support fuel as well as the ability to operate at lower loads. The utilities’ ability to better match the demand curve will result in significant savings. Maintenance is lower for an AC system than a DC system as it operates at lower current. This eliminates the need for a demineralized cooling water system and provides longer electrode life which translates into both material and labor savings. It is the objective of GE to not only demonstrate the additional low load achievable with a plasma system after best achievable tuning, sensor and software approach has been exhausted, but also the increased stability of the flame at all loads with plasma assistance as well as cost savings at all low loads using plasma instead of oil. Upon successful completion of this project, GE will have sufficient field experience to rapidly deploy the Plasma Technology. The project objectives were achieved through the implementation of plasma ignitor technology at PacifiCorp Hunter Station Unit 3. A plasma ignitor system was retrofitted on ten wall-fired burners, Mill 3-4 combustion system. The Plasma Ignitors installed at Hunter proved that this GE technology is a direct and complete replacement for the original oil ignitors. The Hunter Unit 3 burner management system allows the plasma system to be used in all applications that originally required oil to be burned. This includes any time the Mill 3-4 is started or stopped for any reason including boiler starts, load changes, and low load support.
Not provided.
The aim of the present study is to analyze the effect of the electron cyclotron heating (ECH) on the linear stability of Alfvén eigenmodes (AEs) and energetic particle modes (EPMs) triggered by energetic ions in Heliotron J plasma. The analysis is performed using the FAR3d code that solves a reduced MHD model to describe the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The simulations reproduce the AE/EPM stability trends observed in the experiments as the electron temperature (T e ) increases, modifying the thermal plasma β, EP β and EP slowing-down time. Particularly, the $n/m = 1/2$ EPM and $2/4$ Global AE are stabilized in the low-bumpiness (LB) configuration due to an enhancement of the continuum, finite Larmor radius and e-i Landau damping effects as the thermal β increases. On the other hand, a larger ECH injection power cannot stabilize the AE/EPM in medium-bumpiness and high-bumpiness (HB) configurations because the damping effects are weaker compared to the LB case, unable to balance the further destabilization induced by an enhanced EP resonance as the EP slowing-down time and EP β increases with T e .
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.
Polycarbonate films are subjected to plasma treatment in a number of applications such as improving adhesion between polycarbonate and silicon alloy in protective and optical coatings. The changes in surface chemistry due to plasma treatment have tendency to revert back. Thus stability of the plasma induced changes on polymer surfaces over desired time period is very important. The objective of this study was to examine the effect of ageing on atmospheric pressure helium-plasma treated polycarbonate (PC) sample as a function of treatment time. The ageing effects were studied over a period of 10 days. The samples were plasma treated for 0.5, 2, 5 and 10 minutes. Contact angle measurements were made to study surface energy changes. Modification of surface chemical structure was examined using, X-ray Photoelectron Spectroscopy (XPS). Contact angle measurements on untreated and plasma treated surfaces were made immediately, 24, 48, 72 and 96 hrs after treatment. Contact angle decreased from 93 deg for untreated sample to 30 deg for sample plasma treated for 10 minutes. After 10 days the contact angles for the 10 minute plasma treated sample increased to 67 deg, but it never reverted back to that of untreated surface. Similarly the O/C ratio increased from 0.136 for untreated sample to 0.321 for 10 minute plasma treated sample indication increase in surface energy.
An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.
The linearized equation of motion is given in a Lagrangian representation for a rotating plasma with anisotropic pressure. A WKB theory is developed for large-n ballooning modes in an axisymmetric configuration with field-aligned and rigid toroidal flows. In the presence of field-aligned flows, it is shown that a resonance occurs which is strongly suggestive of a generalized mirror instability. In the presence of toroidal rotation, a possible stabilizing effect is identified for P(normal) greater than P(parallel). Finally, as a special case of the theory, the necessary and sufficient conditions for stability in a static, anisotropic plasma are obtained.
This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.
Abstract The influence of energetic particles (EPs) on the ideal internal kink mode, in rotating tokamak plasmas, is numerically investigated by simultaneously solving MHD-kinetic hybrid equations together with a toroidal momentum balance equation utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). The neoclassical toroidal viscous (NTV) torque, induced by precessional drift resonances of trapped energetic particles, acts as the momentum sink term to damp the plasma flow. Quasi-linear initial value simulations show local reduction of the flow amplitude and enhancement of the flow shear near the q = 1 rational surface ( q is the safety factor) due to EP induced NTV. Both effects in turn destabilize the internal kink mode. These numerical findings are robust against the initial linear stability of internal kink, the initial plasma flow profile, as well as the equilibrium distribution model for EPs.
The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.
Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.