Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MHD Instability”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Mixing of the Interstellar and Solar Plasmas at the Heliospheric Interface

From the ideal MHD perspective, the heliopause is a tangential discontinuity that separates the solar wind plasma from the local interstellar medium plasma. There are physical processes, however, that make the heliopause permeable. They can be subdivided into kinetic and MHD categories. Kinetic processes occur on small length and time scales, and cannot be resolved with MHD equations. On the other hand, MHD instabilities of the heliopause have much larger scales and can be easily observed by spacecraft. The heliopause may also be a subject of magnetic reconnection. In this paper, we discuss mechanisms of plasma mixing at the heliopause in the context of Voyager 1 observations. Numerical results are obtained with a Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS), which is a package of numerical codes capable of performing adaptive mesh refinement simulations of complex plasma flows in the presence of discontinuities and charge exchange between ions and neutral atoms. The flow of the ionized component is described with the ideal MHD equations, while the transport of atoms is governed either by the Boltzmann equation or multiple Euler gas dynamics equations. The code can also treat nonthermal ions and turbulence produced by them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
DOE OSTI2016

Discovery of stationary operation of quiescent H-mode plasmas with net-zero neutral beam injection torque and high energy confinement on DIII-D

Recent experiments in DIII-D [J. L. Luxon et al., in Plasma Physics and Controlled Nuclear Fusion Research 1996 (International Atomic Energy Agency, Vienna, 1987), Vol. I, p. 159] have led to the discovery of a means of modifying edge turbulence to achieve stationary, high confinement operation without Edge Localized Mode (ELM) instabilities and with no net external torque input. Eliminating the ELM-induced heat bursts and controlling plasma stability at low rotation represent two of the great challenges for fusion energy. By exploiting edge turbulence in a novel manner, we achieved excellent tokamak performance, well above the H{sub 98y2} international tokamak energy confinement scaling (H{sub 98y2} = 1.25), thus meeting an additional confinement challenge that is usually difficult at low torque. The new regime is triggered in double null plasmas by ramping the injected torque to zero and then maintaining it there. This lowers E × B rotation shear in the plasma edge, allowing low-k, broadband, electromagnetic turbulence to increase. In the H-mode edge, a narrow transport barrier usually grows until MHD instability (a peeling ballooning mode) leads to the ELM heat burst. However, the increased turbulence reduces the pressure gradient, allowing the development of a broader and thus higher transport barrier. A 60% increase in pedestal pressure and 40% increase in energy confinement result. An increase in the E × B shearing rate inside of the edge pedestal is a key factor in the confinement increase. Strong double-null plasma shaping raises the threshold for the ELM instability, allowing the plasma to reach a transport-limited state near but below the explosive ELM stability boundary. The resulting plasmas have burning-plasma-relevant β{sub N} = 1.6–1.8 and run without the need for extra torque from 3D magnetic fields. To date, stationary conditions have been produced for 2 s or 12 energy confinement times, limited only by external hardware constraints. Stationary operation with improved pedestal conditions is highly significant for future burning plasma devices, since operation without ELMs at low rotation and good confinement is key for fusion energy production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Magnetohydrodynamic thermal instabilities in cool inhomogeneous atmospheres

The stability of magnetic loops to current-driven filamentation instabilities is investigated. The unperturbed atmosphere is assumed to be composed of an (upper) isothermal optically thin low-density portion and a (lower) higher-density portion which is in radiative equilibrium; in both cases, the atmosphere is in hydrostatic equilibrium, so that gravitational stratification is taken into account. In order to provide specific equilibrium conditions for evaluation of the dispersion relation, conditions appropriate for the surface of a solar-type star are adopted; i.e., a fairly low temperature (T = 5000 K) appropriate for a 'precoronal' state associated, for example, with magnetic flux emerging from photospheric levels under the action of magnetic buoyancy. A linear stability analysis is performed, and numerical results show that physically plausible current densities, which would be generated by typical loop-footpoint motions, are effective in driving MHD instabilities in such a plasma. The instability growth rates are strongly dependent on the assumed current density distribution and on the density scale height.

Bodo, G.↗
DOE OSTI2023

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Heating and current drive requirements towards steady state operation in ITER

Steady state scenarios envisaged for ITER aim at optimizing the bootstrap current, while maintaining sufficient confinement and stability to provide the necessary fusion yield. Non-inductive scenarios will need to operate with Internal Transport Barriers (ITBs) in order to reach adequate fusion gain at typical currents of 9 MA. However, the large pressure gradients associated with ITBs in regions of weak or negative magnetic shear can be conducive to ideal MHD instabilities, reducing the no-wall limit. The E × B flow shear from toroidal plasma rotation is expected to be low in ITER, with a major role in the ITB dynamics being played by magnetic geometry. Combinations of H/CD sources that maintain weakly reversed magnetic shear profiles throughout the discharge are the focus of this work. Time-dependent transport simulations indicate that, with a trade-off of the EC equatorial and upper launcher, the formation and sustainment of quasi-steady state ITBs could be demonstrated in ITER with the baseline heating configuration. However, with proper constraints from peeling-ballooning theory on the pedestal width and height, the fusion gain and the maximum non-inductive current are below the ITER target. Upgrades of the heating and current drive system in ITER, like the use of Lower Hybrid current drive, could overcome these limitations, sustaining higher non-inductive current and confinement, more expanded ITBs which are ideal MHD stable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Progress toward fully noninductive discharge operation in DIII-D using off-axis neutral beam injection

The initial experiments on off-axis neutral beam injection into high noninductive current fraction (f{sub NI}), high normalized pressure (β{sub N}) discharges in DIII-D [J. L. Luxon, Fusion Sci. Technol. 48, 828 (2005)] have demonstrated changes in the plasma profiles that increase the limits to plasma pressure from ideal low-n instabilities. The current profile is broadened and the minimum value of the safety factor (q{sub min}) can be maintained above 2 where the profile of the thermal component of the plasma pressure is found to be broader. The off-axis neutral beam injection results in a broadening of the fast-ion pressure profile. Confinement of the thermal component of the plasma is consistent with the IPB98(y,2) scaling, but global confinement with q{sub min}>2 is below the ITER-89P scaling, apparently as a result of enhanced transport of fast ions. A 0-D model is used to examine the parameter space for f{sub NI}=1 operation and project the requirements for high performance steady-state discharges. Fully noninductive solutions are found with 4<β{sub N}<5 and bootstrap current fraction near 0.5 for a weak shear safety factor profile. A 1-D model is used to show that a f{sub NI}=1 discharge at the top of this range of β{sub N} that is predicted stable to n=1, 2, and 3 ideal MHD instabilities is accessible through further broadening of the current and pressure profiles with off-axis neutral beam injection and electron cyclotron current drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2023

Nonlinear control of safety factor gradient in tokamaks using spatially variable electron cyclotron current drives

Active control of plasma properties may be necessary to achieve stable operation of next-generation tokamaks over large time scales. Such control algorithms can regulate the plasma properties to avoid the onset of magnetohydrodynamic (MHD) instabilities. For instance, the global and local properties of the safety factor profile are linked to the onset of neoclassical tearing modes (NTMs). This work proposes a model-based control approach for deterring/delaying NTMs through active modulation of one of the safety factor properties - the gradient of the safety factor profile at a particular rational safety factor surface. In particular, a novel control-oriented model for the local safety factor gradient is developed. The nonlinear control model is governed by a nonautonomous ordinary differential equation that accounts for a given rational safety factor surface’s spatial variation over time. Further, to improve the controllability of the spatially evolving parameter, the control model treats ECH&CD positions, along with noninductive powers, as controllable variables. A nonlinear control algorithm based on feedback linearization with optimization is synthesized to achieve the objective of regulating the safety factor gradient around a given target. The proposed algorithm allocates optimal ECH&CD positions, in addition to auxiliary powers, at each time instant as the rational safety factor surface drifts to locations with a low control authority. Stability guarantees of the proposed control law are also discussed in this work. The proposed algorithm is tested for a DIII-D tokamak scenario in nonlinear simulations carried out using the Control Oriented Transport SIMulator (COTSIM). Both fixed and moving ECH&CD cases are studied, and their outcomes are compared. Simulation results demonstrate that enthusiastic regulation of the safety factor gradient can be achieved during the ramp-up and flat-top phases of tokamak operation in both fixed and moving ECH&CD cases. However, real-time updates of ECH&CD positions can prevent the saturation of auxiliary powers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Comparative H-mode density limit studies in JET and AUG

Identification of the mechanisms for the H-mode density limit in machines with fully metallic walls, and their scaling to future devices is essential to find for these machines the optimal operational boundaries with the highest attainable density and confinement. Systematic investigations of H-mode density limit plasmas in experiments with deuterium external gas fuelling have been performed on machines with fully metallic walls, JET and AUG and results have been compared with one another. Basically, the operation phases are identical for both tokamaks: the stable H-mode phase, degrading H-mode phase, breakdown of the H-mode with energy confinement deterioration usually accompanied by a dithering cycling phase, followed by the l-mode phase. The observed H-mode density limit on both machines is found close to the Greenwald limit (n/n GW = 0.8–1.1 in the observed magnetic configurations). The similar behavior of the radiation on both tokamaks demonstrates that the density limit (DL) is neither related to additional energy losses from the confined region by radiation, nor to an inward collapse of the hot discharge core induced by overcooling of the plasma periphery by radiation. It was observed on both machines that detachment, as well as the X-point MARFE itself, does not trigger a transition in the confinement regime and thus does not present a limit on the plasma density. It is the plasma confinement, most likely determined by edge parameters, which is ultimately responsible for the transition from H- to l-mode. The measured Greenwald fractions are found to be consistent with the predictions from different theoretical models, based on MHD instability theory in the near-SOL.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The Explosion Mechanism of Core-Collapse Supernovae: Progress in Supernova Theory and Experiments

The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention to its asymmetric character. We highlight a series of successful studies connecting observations of supernova remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well as the need to gain a better understanding of hydrodynamical and MHD instabilities such as SASI and neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative system where buoyancy effects are absent. This dynamical system can be studied numerically and also experimentally with a water fountain. Lastly, we analyse the potential of this complementary research tool for supernova theory. We also review its potential for public outreach in science museums.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Fast-ion physics in SPARC

Potential loss of energetic ions including alphas and radio-frequency tail ions due to classical orbit effects and magnetohydrodynamic instabilities (MHD) are central physics issues in the design and experimental physics programme of the SPARC tokamak. The expected loss of fusion alpha power due to ripple-induced transport is computed for the SPARC tokamak design by the ASCOT and SPIRAL orbit-simulation codes, to assess the expected surface heating of plasma-facing components. We find good agreement between the ASCOT and SPIRAL simulation results not only in integrated quantities (fraction of alpha power loss) but also in the spatial, temporal and pitch-angle dependence of the losses. If the toroidal field (TF) coils are well-aligned, the SPARC edge ripple is small (0.15–0.30 %), the computed ripple-induced alpha power loss is small ( ~0.25% ) and the corresponding peak surface power density is acceptable ( 244 kW m -2 ). However, the ripple and ripple-induced losses increase strongly if the TF coils are assumed to suffer increasing magnitudes of misalignment. Surface heat loads may become problematic if the TF coil misalignment approaches the centimetre level. Ripple-induced losses of the energetic ion tail driven by ion cyclotron range of frequency (ICRF) heating are not expected to generate significant wall or limiter heating in the nominal SPARC plasma scenario. Because the expected classical fast-ion losses are small, SPARC will be able to observe and study fast-ion redistribution due to MHD including sawteeth and Alfvén eigenmodes (AEs). SPARC's parameter space for AE physics even at moderate Q is shown to reasonably overlap that of the demonstration power plant ARC (Sorbom et al., Fusion Engng Des., vol. 100, 2015, p. 378), and thus measurements of AE mode amplitude, spectrum and associated fast-ion transport in SPARC would provide relevant guidance about AE behaviour expected in ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Bootstrap current modeling in M3D-C1

Bootstrap current plays a crucial role in the equilibrium of magnetically confined plasmas, particularly in quasi-symmetric stellarators and in tokamaks, where it can represent bulk of the electric current density. Accurate modeling of this current is essential for understanding the magnetohydrodynamic (MHD) equilibrium and stability of these configurations. This study expands the modeling capabilities of M3D-C1, an extended-MHD code, by implementing self-consistent physics models for bootstrap current. It employs two analytical frameworks: a generalized Sauter model (Sauter et al. 1999 Phys. Plasmas vol. 6, no. 7, pp. 2834–2839), and a revised Sauter-like model (Redl et al. 2021 Phys. Plasmas vol. 28, no. 2, pp. 022502). The isomorphism described by Landreman et al. (2022 Phys. Rev. Lett. vol. 128, pp. 035001) is employed to apply these models to quasi-symmetric stellarators. The implementation in M3D-C1 is benchmarked against neoclassical codes, including NEO, XGCa and SFINCS, showing excellent agreement. These improvements allow M3D-C1 to self-consistently calculate the neoclassical contributions to plasma current in axisymmetric and quasi-symmetric configurations, providing a more accurate representation of the plasma behavior in these configurations. A workflow for evaluating the neoclassical transport using SFINCS with arbitrary toroidal equilibria calculated using M3D-C1 is also presented. This workflow enables a quantitative evaluation of the error in the Sauter-like model in cases that deviate from axi- or quasi-symmetry (e.g. through the development of an MHD instability).

fusion plasma↗
DOE OSTI2012

Lithium-based surfaces controlling fusion plasma behavior at the plasma-material interface

The plasma-material interface and its impact on the performance of magnetically confined thermonuclear fusion plasmas are considered to be one of the key scientific gaps in the realization of nuclear fusion power. At this interface, high particle and heat flux from the fusion plasma can limit the material's lifetime and reliability and therefore hinder operation of the fusion device. Lithium-based surfaces are now being used in major magnetic confinement fusion devices and have observed profound effects on plasma performance including enhanced confinement, suppression and control of edge localized modes (ELM), lower hydrogen recycling and impurity suppression. The critical spatial scale length of deuterium and helium particle interactions in lithium ranges between 5-100 nm depending on the incident particle energies at the edge and magnetic configuration. Lithium-based surfaces also range from liquid state to solid lithium coatings on a variety of substrates (e.g., graphite, stainless steel, refractory metal W/Mo/etc., or porous metal structures). Temperature-dependent effects from lithium-based surfaces as plasma facing components (PFC) include magnetohydrodynamic (MHD) instability issues related to liquid lithium, surface impurity, and deuterium retention issues, and anomalous physical sputtering increase at temperatures above lithium's melting point. The paper discusses the viability of lithium-based surfaces in future burning-plasma environments such as those found in ITER and DEMO-like fusion reactor devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Compact solid-state neutral particle analyzer in current mode

Solid state neutral particle analyzer (ssNPA) arrays are operated in current mode on the DIII-D tokamak and the National Spherical Torus Experiment (NSTX). Compared with conventional pulse-counting NPAs, current-mode operation sacrifices energy resolution to obtain economical, high-bandwidth, pitch-angle resolved measurements. With the success from a new three-channel near-vertical-view current mode ssNPA on DIII-D, the apertures on an existing array on NSTX were expanded to increase the particle influx. The sightlines of both arrays intersect heating beams, enabling both active and passive charge exchange measurements. The spatial resolution at beam intersection is typically 5 cm on both devices. Directly deposited ultra-thin foils on the detector surface block stray photons below the energy of 1 keV and also set low energy threshold about 25 keV for deuterium particle detection. Oscillations in neutral flux produced by high frequency magnetohydrodynamics (MHD) instabilities are readily detected.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗