Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA CONFINEMENT”

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 1,225 records · Page 68

Tritium Aspects of Fueling and Exhaust Pumping in Magnetic Fusion Energy

Magnetically confined fusion plasmas generate energy from deuterium-tritium (DT) fusion reactions that produce energetic 3.5 MeV alpha particles and 14 MeV neutrons. Since the DT fusion reaction rate is a strong function of plasma density, an efficient fueling source is needed to maintain high plasma density in such systems. Energetic ions in fusion plasmas are able to escape the confining magnetic fields at a much higher rate than the fusion reactions occur, thus dictating the fueling rate needed. These lost ions become neutralized and need to be pumped away as exhaust gas to be reinjected into the plasma as fuel atoms.The technology to fuel and pump fusion plasmas has to be inherently compatible with the tritium fuel. An ideal holistic solution would couple the pumping and fueling such that the pump exhaust is directly fed back into pellet formation without including impurity gases. This would greatly reduce the processing needs for the exhaust. Concepts to accomplish this are discussed along with the fueling and pumping needs for a DT fusion reactor.

Baylor, Larry↗

Direct observations of L-I-H and H-I-L transitions with the X-point reciprocating probe in ASDEX Upgrade

A reciprocating Langmuir probe was used to directly measure the behavior of turbulence and flows in the X-point region during transitions between low-(L) and high-confinement (H) mode in ASDEX Upgrade. The probe traverses the divertor horizontally in 140 ms, typically 2–5 cm below the X-point. Toroidal Mach number, density, floating potential (ϕ{sub f}), and electron temperature (T{sub e}) are measured. In the regime accessible to the probe (P{sub inj}<1.5 MW, line-integrated core density <4×10{sup 19} m{sup −2}), the L-H transition features an intermediate phase (I-phase), characterized by limit-cycle oscillations at 0.5–3 kHz [Conway et al., Phys. Rev. Lett. 106, 065001 (2011)]. The probe measurements reveal that this pulsing affects both the density and the toroidal Mach number. It is present in both the low-(LFS) and high-field sides (HFS) of the scrape-off layer, while high-amplitude broadband turbulence usually dominates the private-flux region. Profile comparisons between L-mode and I-phase show lower density in pulsing regions and small shifts in T{sub e}, directed oppositely on LFS and HFS, which are compensated by shifts in ϕ{sub f} to yield a surprisingly unchanged plasma potential profile. Directly observed L-I-phase transitions reveal that the onset of the pulsing is preceded by a fast 50% density drop in the HFS X-point region. Back transitions to L-mode occur essentially symmetrically, with the pulsing stopping first, followed by a fast recovery to L-mode density levels in the divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linking the micro and macro: L-H transition dynamics and threshold physics

The links between the microscopic dynamics and macroscopic threshold physics of the L → H transition are elucidated. Emphasis is placed on understanding the physics of power threshold scalings, and especially on understanding the minimum in the power threshold as a function of density P{sub thr} (n). By extending a numerical 1D model to evolve both electron and ion temperatures, including collisional coupling, we find that the decrease in P{sub thr} (n) along the low-density branch is due to the combination of an increase in collisional electron-to-ion energy transfer and an increase in the heating fraction coupled to the ions. Both processes strengthen the edge diamagnetic electric field needed to lock in the mean electric field shear for the L→H transition. The increase in P{sub thr} (n) along the high-density branch is due to the increase with ion collisionality of damping of turbulence-driven shear flows. Turbulence driven shear flows are needed to trigger the transition by extracting energy from the turbulence. Thus, we identify the critical transition physics components of the separatrix ion heat flux and the zonal flow excitation. The model reveals a power threshold minimum in density scans as a crossover between the threshold decrease supported by an increase in heat fraction received by ions (directly or indirectly, from electrons) and a threshold increase, supported by the rise in shear flow damping. The electron/ion heating mix emerges as important to the transition, in that it, together with electron-ion coupling, regulates the edge diamagnetic electric field shear. The importance of possible collisionless electron-ion heat transfer processes is explained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The role of toroidal rotation in the very high energy confinement quality observed in super H-mode experiments on DIII-D

In this paper, we report the key role that toroidal rotation and the related E × B shear physics played in the very high energy confinement quality (H 98y2 > 1.5) of super H-mode experiments on DIII-D. Experiments show that the energy confinement quality decreases when toroidal rotation decreases due to the decreased externally controlled torque per particle. Meanwhile, the total pedestal pressure in the experiments remains very high during the rotation and confinement quality change. TGYRO transport modeling suggests the contribution from rotation in the E × B shear is responsible for the confinement quality in excess of standard H-mode (H 98y2 ~ 1). CGYRO gyrokinetic simulations reveal the governing physics in the core plasma of super H-modes: significant up-shift of nonlinear the ITG critical gradient is observed when applying E × B shear physics in the modeling based on experimental data. Furthermore, the effects of other physical parameters and contribution from pedestal height, which may play minor roles in this study, are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of vertical displacement oscillatory modes and global Alfvén Eigenmodes in JET geometry

Vertical Displacement Oscillatory Modes (VDOM), with frequency in the Alfvén range, are natural modes of oscillation of magnetically confined laboratory plasmas with elongated cross-section. These axisymmetric modes arise from the interaction between the plasma current, which is in equilibrium with currents flowing in external coils, and perturbed currents induced on a nearby conducting wall. The restoring force exerted by these perturbed currents on the vertical motion of the plasma column leads to its oscillatory behavior. An analytic model for VDOM was proposed based on an idealized 'straight tokamak' equilibrium with uniform equilibrium current density. This article introduces the first numerical simulations of VDOM in a realistic JET tokamak configuration, using the extended-MHD code NIMROD and drawing comparisons with Global Alfvén Eigenmodes (GAE). The results show qualitative agreement with analytic predictions regarding mode frequency and radial structure, supporting the identification of VDOM as a fundamental oscillation mode in tokamak plasmas. VDOM and GAE are modeled in a representative JET discharge, where axisymmetric perturbations with toroidal mode number n = 0 driven unstable by fast ions were observed. The two modes are examined separately using a forced oscillator within the NIMROD code, which enables a comparison of their characteristics and helps identify the experimentally observed mode possibly as a GAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temporal parallelization of edge plasma simulations using the parareal algorithm and the SOLPS code

We show that numerical modelling of edge plasma physics may be successfully parallelized in time. The parareal algorithm has been employed for this purpose and the SOLPS code package coupling the B2.5 finite-volume fluid plasma solver with the kinetic Monte-Carlo neutral code Eirene has been used as a test bed. The complex dynamics of the plasma and neutrals in the scrape-off layer (SOL) region makes this a unique application. It is demonstrated that a significant computational gain (more than an order of magnitude) may be obtained with this technique. The use of the IPS framework for event-based parareal implementation optimizes resource utilization and has been shown to significantly contribute to the computational gain.

97 MATHEMATICS AND COMPUTING↗

On the role of spectral resolution in velocity shear layer measurements by Doppler reflectometry

The signal quality of a Doppler reflectometer depends strongly on its spectral resolution, which is influenced by the microwave beam properties and the radius of curvature of the cutoff layer in the plasma. If measured close to a strong perpendicular velocity shear layer, the spectrum of the backscattered signal is influenced by different velocities. This can give rise to two Doppler shifted peaks in the spectrum as observed in TJ-II H-mode plasmas. It is shown by two-dimensional full wave simulations that the two peaks are separable provided the spectral resolution of the system is sufficient. However, if the spectral resolution is poor, the two peaks blend into one and yield an intermediate and incorrect velocity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The EPED pedestal model and edge localized mode-suppressed regimes: Studies of quiescent H-mode and development of a model for edge localized mode suppression via resonant magnetic perturbations

The EPED model predicts the H-mode pedestal height and width based upon two fundamental and calculable constraints: (1) onset of non-local peeling-ballooning modes at low to intermediate mode number, (2) onset of nearly local kinetic ballooning modes at high mode number. We present detailed tests of the EPED model in discharges with edge localized modes (ELMs), employing new high resolution measurements, and finding good quantitative agreement across a range of parameters. The EPED model is then applied for the first time to quiescent H-mode (QH), finding a similar level of agreement between predicted and observed pedestal height and width, and suggesting that the model can be used to predict the critical density for QH-mode operation. Finally, the model is applied toward understanding the suppression of ELMs with 3D resonant magnetic perturbations (RMP). Combining EPED with plasma response physics, a new working model for RMP ELM suppression is developed. We propose that ELMs are suppressed when a 'wall' associated with the RMP blocks the inward penetration of the edge transport barrier. A calculation of the required location of this 'wall' with EPED is consistent with observed profile changes during RMP ELM suppression and offers an explanation for the observed dependence on safety factor (q{sub 95}).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of magnetic fluctuations and rapid density decay of magnetospheric plasma in Ring Trap 1

The Ring Trap 1 device, a magnetospheric configuration generated by a levitated dipole field magnet, has created high-{beta} (local {beta} {approx} 70%) plasma by using electron cyclotron resonance heating (ECH). When a large population of energetic electrons is generated at low neutral gas pressure operation, high frequency magnetic fluctuations are observed. When the fluctuations are strongly excited, rapid loss of plasma was simultaneously observed especially in a quiet decay phase after the ECH microwave power is turned off. Although the plasma is confined in a strongly inhomogeneous dipole field configuration, the frequency spectra of the fluctuations have sharp frequency peaks, implying spatially localized sources of the fluctuations. The fluctuations are stabilized by decreasing the hot electron component below approximately 40%, realizing stable high-{beta} confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stationary density profiles in the Alcator C-mod tokamak

In the absence of an internal particle source, plasma turbulence will impose an intrinsic relationship between an inwards pinch and an outwards diffusion resulting in a stationary density profile. The Alcator C-mod tokamak utilizes RF heating and current drive so that fueling only occurs in the vicinity of the separatrix. Discharges that transition from L-mode to I-mode are seen to maintain a self-similar stationary density profile as measured by Thomson scattering. For discharges with negative magnetic shear, an observed rise of the safety factor in the vicinity of the magnetic axis appears to be accompanied by a decrease of electron density, qualitatively consistent with the theoretical expectations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of the transport shortfall in Alcator C-Mod L-mode plasmas

A so-called 'transport shortfall,' where ion and electron heat fluxes and turbulence are underpredicted by gyrokinetic codes, has been robustly identified in DIII-D L-mode plasmas for {rho}>0.55[T. L. Rhodes et al., Nucl. Fusion 51(6), 063022 (2011); and C. Holland et al., Phys. Plasmas 16(5), 052301 (2009)]. To probe the existence of a transport shortfall across different tokamaks, a dedicated scan of auxiliary heated L-mode discharges in Alcator C-Mod are studied in detail with nonlinear gyrokinetic simulations for the first time. Two discharges, only differing by the amount of auxiliary heating are investigated using both linear and nonlinear simulation of the GYRO code [J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)]. Nonlinear gyrokinetic simulation of the low and high input power discharges reveals a discrepancy between simulation and experiment in only the electron heat flux channel of the low input power discharge. However, both discharges demonstrate excellent agreement in the ion heat flux channel, and the high input power discharge demonstrates simultaneous agreement with experiment in both the electron and ion heat flux channels. A summary of linear and nonlinear gyrokinetic results and a discussion of possible explanations for the agreement/disagreement in each heat flux channel is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Heating of thermal non-equilibrium ions by Alfvén wave via nonresonant interaction

Pickup of thermal non-equilibrium ions by Alfvén wave via nonresonant wave-particle interaction is investigated by means of analytical test-particle theory. Some interesting and new results are found. No matter what the initial velocity distribution is, if the background magnetic field, the Alfvén speed, and the Alfvén magnetic field are fixed, the average parallel velocity never changes when t→∞. Heating effects in the perpendicular and parallel direction just depend on the initial temperature, and the perpendicular temperature increase is more prominent. It is noted that the heating effect of thermal non-equilibrium ions (Kappa ions) is weaker than that of the Maxwellian. This phenomenon may be relative to the heating of ions in the solar corona as well as in some toroidal confinement fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast wave direct electron heating in advanced inductive and ITER baseline scenario discharges in DIII-D

Fast Wave (FW) heating and electron cyclotron heating (ECH) are used in the DIII-D tokamak to study plasmas with low applied torque and dominant electron heating characteristic of burning plasmas. FW heating via direct electron damping has reached the 2.5 MW level in high performance ELMy H-mode plasmas. In Advanced Inductive (AI) plasmas, core FW heating was found to be comparable to that of ECH, consistent with the excellent first-pass absorption of FWs predicted by ray-tracing models at high electron beta. FW heating at the ∼2 MW level to ELMy H-mode discharges in the ITER Baseline Scenario (IBS) showed unexpectedly strong absorption of FW power by injected neutral beam (NB) ions, indicated by significant enhancement of the D-D neutron rate, while the intended absorption on core electrons appeared rather weak. The AI and IBS discharges are compared in an effort to identify the causes of the different response to FWs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Short wavelength turbulence generated by shear in the quiescent H-mode edge on DIII–D

A region of turbulence with large radial wavenumber (k{sub r}ρ{sub s}>1) is found in the high-shear portion of the plasma edge in Quiescent H-mode (QH-mode) on DIII–D using the Phase Contrast Imaging (PCI) diagnostic. At its peak outside the minimum of the E{sub r} well, the turbulence exhibits large amplitude n{sup ~}/n∼40%, with large radial wavenumber |k{sup ¯}{sub r}/k{sup ¯}{sub θ}|∼11 and short radial correlation length L{sub r}/ρ{sub i}∼0.2. The turbulence inside the E{sub r} well minimum is characterized by the opposite sign in radial wavenumber from that of turbulence outside the minimum, consistent with the expected effects of velocity shear. The PCI diagnostic provides a line-integrated measurement of density fluctuations, so data are taken during a scan of plasma position at constant parameters to allow the PCI to sample a range in k{sub r}/k{sub θ}. Analysis of the Doppler shift and plasma geometry allows the turbulence to be localized to a narrow region 3 mm inside the last closed flux surface, outside the minimum of the E{sub r} well. The turbulence amplitude and radial wavenumber and correlation length are determined by fitting the PCI results with a simple non-isotropic turbulence model with two regions of turbulence. These PCI observations, made in QH-mode, are qualitatively similar to those made in standard edge localized modes (ELM)-free H-mode and between ELMs, suggesting a similar role for large k{sub r} turbulence there.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sustained spheromaks with ideal n = 1 kink stability and pressure confinement

Increasing the helicity injector drive frequency up to 68.5 kHz on the Helicity Injected Torus-Steady Inductive (HIT-SI) experiment has produced spheromaks with current amplifications of 3.8, ideal n = 1 kink stability, improved toroidal symmetry and pressure confinement. Current centroid calculations from surface magnetic probes show an outward shift in the magnetic field at frequencies above 50 kHz. Grad-Shafranov equilibria indicate pressure confinement at higher injector operating frequencies. The minimum characteristic frequency needed to achieve this confining effect on HIT-SI plasmas is found to be approximately 30 kHz by analysis of the density fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗