Engineering topics

Heidbrink, W. W.

Publications and source records attributed to Heidbrink, W. W..

18 records

DOE OSTI2024

A research program to measure the lifetime of spin polarized fuel

The use of spin polarized fuel could increase the deuterium-tritium (D-T) fusion cross section by a factor of 1.5 and, owing to alpha heating, increase the fusion power by an even larger factor. Issues associated with the use of polarized fuel in a reactor are identified. Theoretically, nuclei remain polarized in a hot fusion plasma. The similarity between the Lorentz force law and the Bloch equations suggests polarization can be preserved despite the rich electromagnetic spectrum present in a magnetic fusion device. The most important depolarization mechanisms can be tested in existing devices. The use of polarized deuterium and 3 He in an experiment avoids the complexities of handling tritium, while encompassing the same nuclear reaction spin-physics, making it a useful proxy to study issues associated with full D-T implementation. 3 He fuel with 65% polarization can be prepared by permeating optically-pumped 3 He into a shell pellet. Dynamically polarized 7 Li-D pellets can achieve 70% vector polarization for the deuterium. Cryogenically-frozen pellets can be injected into fusion facilities by special injectors that minimize depolarizing field gradients. Alternatively, polarized nuclei could be injected as a neutral beam. Once injected, the lifetime of the polarized fuel is monitored through measurements of escaping charged fusion products. Multiple experimental scenarios to measure the polarization lifetime in the DIII-D tokamak and other magnetic-confinement facilities are discussed, followed by outstanding issues that warrant further study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Time-resolved biphase signatures of quadratic nonlinearity observed in coupled Alfvén eigenmodes on the DIII-D tokamak

We report the detection of nonstationary quadratic coupling between toroidicity-induced Alfvén eigenmodes (TAEs) on sub-millisecond time scales. Identification of phase coherency between multiple TAEs and nonlinearly generated modes is facilitated by wavelet-based bicoherence analysis of time-series from inductive coils, taken from a DIII-D discharge heated by neutral beam injection (NBI). Characterization of nonlinear three-wave interaction is inferred by stationary local bispectrum phase (biphase) and confirmed via bandpass filtering. Biphase dynamics associated with prominent bispectral features are well-resolved in time and consistent with transient quadratic coupling. Onset and duration of nonlinearity are correlated with enhanced amplitude of participating TAEs; coincident changes in amplitude are observed for modes at difference frequency |f TAE,1 - f TAE,2 |.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Low-frequency shear Alfvén waves at DIII-D: Theoretical interpretation of experimental observations

The linear properties of the low-frequency shear Alfvén waves such as those associated with the beta-induced Alfvén eigenmodes (BAEs) and the low-frequency modes observed in reversed-magnetic-shear DIII-D discharges [W. Heidbrink et al., Nucl. Fusion 61, 066031 (2021)] are theoretically investigated and delineated based on the theoretical framework of the general fishbone-like dispersion relation (GFLDR). By adopting representative experimental equilibrium profiles, it is found that, even though both modes are predominantly of Alfvénic polarization, the low-frequency mode is a reactive unstable mode with weak coupling to the energetic particles, while the BAE involves a dissipative instability due to resonant excitation by the energetic ions. Thus, the low-frequency mode is more appropriately called a low-frequency Alfvén mode (LFAM). Moreover, the ascending frequency spectrum patterns of the experimentally observed BAEs and LFAMs can be theoretically reproduced by varying qmin and also be well interpreted based on the GFLDR. The present analysis illustrates the solid predictive capability of the GFLDR and its practical usefulness in enhancing the interpretative capability of both experimental and numerical simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Validation of energetic particle transport models for time-dependent integrated simulations of burning plasmas (Final ST Report)

The major goal of this project is to deliver a suite of validated energetic particle (EP) modeling tools for time-dependent integrated tokamak simulations (with the code TRANSP) that can be used for modeling plasma discharges featuring a substantial EP content, such as burning plasmas. During FY20, issues were resolved with an energetic particle transport model called TGLF-EP, with improvements made to the TGLF-EP workflow, generation of input files, and interpretation of runs on the NERSC supercomputing cluster. The recipient also gave input on simulations with the code TRANSP to reproduce the reference JET discharge, and assessed the accuracy and uncertainties of the simulated B field, current evolution, thermal plasma profiles, stored energy, etc. Input was also given on interpretation of fast-ion instability calculations and recommendations for modifications to JET afterglow discharge for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Effect of anisotropic fast ions on internal kink stability in DIII-D negative and positive triangularity plasmas

Recent DIII-D experiments show that sawtooth stability is strongly affected by anisotropic fast ions from neutral beam injection (NBI) in both negative and positive triangularity plasmas. Fast ions from co-current NBI are stabilizing for the sawtooth stability, resulting in longer sawtooth periods. On the other hand, fast ions from counter-current NBI are destabilizing, leading to small and frequent sawteeth. The relative change of sawtooth period and amplitude is more than a factor of two. These observations appear to hold in both plasma shapes. Non-perturbative toroidal modeling, utilizing the magnetohydrodynamic-kinetic hybrid stability code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503), reveals an asymmetric dependence of the stability of the n = 1 (n is the toroidal mode number) internal kink mode on the injection direction of NBI, being qualitatively consistent with the experimentally observed sawtooth behavior. The MARS-K modeling results suggest that anisotropic fast ions affect the mode growth rate and frequency through both adiabatic and non-adiabatic contributions. Here, the asymmetry of the internal kink mode instability relative to the NBI direction is mainly due to the non-adiabatic contribution of passing fast ions, which stabilize (destabilize) the internal kink with the co-(counter-) current NBI as compared to the fluid counterpart. However, finite orbit width (FOW) correction to passing particles partially cancels the asymmetry. Trapped particles are always stabilizing due to precessional drift resonance. Modeling also shows that fast ions affect the internal kink in a similar manner in both negative and positive triangularity plasmas, although being slightly more unstable in the negative triangularity configuration already in the fluid limit. The similarity is mainly attributed to the fact that the mode is localized in the plasma core region, with very similar eigenmode structures in both negative and positive configurations. Furthermore, MARS-K modeling indicates that other factors, such as the plasma rotation and the drift kinetic effects of thermal plasmas, weakly modify the mode stability as compared to the drift kinetic resonance effects and FOW correction of fast ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Modelling the Alfvén eigenmode induced fast-ion flow measured by an imaging neutral particle analyzer

An imaging neutral particle analyzer (INPA) provides energy and radially resolved measurements of the confined fast-ion population ranging from the high-field side to the edge on the midplane of the DIII-D tokamak. In recent experiments, it was used to diagnose fast-ion flow in the INPA-interrogated phase-space driven by multiple, marginally unstable Alfvén eigenmodes (AEs). Here, the key features of this measured fast-ion flow are: (I) a fast-ion flow from q min and the injection energy (81 keV) towards lower energies and plasma periphery.(II) A flow from the same location towards higher energies and the plasma core, (III) a phase-space 'hole' at the injected energy and plasma core and (IV) a pile-up at the plasma core at lower energies (~60 keV). Ad hoc energetic particle diffusivity modelling of TRANSP significantly deviates from the observation. Comparably, a reduced modelling, i.e. a combination of NOVA-K and ASCOT5 code with the measured mode structure and amplitude, generally reproduce some key features of the observed phase-space flow, but largely failed to interpret fast ion depletion near the plasma axis. At last, self-consistent, first-principle multi-phase hybrid simulations that include realistic neutral beam injection and collisions are able to reproduce most features of the time-resolved phase-space flow. During consecutive hybrid phases, an RSAE consistent with the experiment grows and saturates, redistributing the injected fast ions. The resulting synthetic INPA images are in good agreement with the measurement near the injection energy. The simulations track the fast-ion redistribution within the INPA range, confirming that the measured fast-ion flow follows streamlines defined by the intersection of phase-space surfaces of constant magnetic moment μ and constant E' = nE + ωP φ , where n and ω are the instability toroidal mode number and frequency, and E and P φ the ion energy and toroidal canonical momentum. Nonperturbative effects are required to reproduce the depletion of fast ions near the magnetic axis at the injection energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Energetic particle-induced geodesic acoustic modes on DIII-D

Various properties of the energetic particle-induced geodesic acoustic mode (EGAM) are explored in this large database analysis of DIII-D experimental data. EGAMs are n = 0 modes with m = 0 electrostatic potential fluctuations (where n/m = toroidal/poloidal mode number), m = 1 density fluctuations, and m = 2 magnetic fluctuations. The fundamental frequency (~20–40 kHz) of the mode is typically below that of the traditional geodesic acoustic mode frequency. EGAMs are most easily destabilized by beams in the counter plasma current (counter-I p ) direction as compared to co-Ip and off-axis beams. During counter beam injection, the mode frequency is found to have the strongest linear dependence (correlation coefficient r = –0.71) with the safety factor (q). Here, the stability of the mode in the space of q and poloidal beta (β p ) shows a clear boundary for the mode stability. The stability of the mode depends more strongly on damping rate than on fast-ion drive for a given injection geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Mode structure measurements of ion cyclotron emission and sub-cyclotron modes on DIII-D

We report mode structure measurements of coherent ion cyclotron emission (ICE) and sub-cyclotron modes on DIII-D. Through a dedicated experiment, we aimed to characterize a variety of modes via the upgraded ICE diagnostic in both L- and H-mode plasmas. In the L-mode plasmas, autopower spectrum peaks at harmonics of the ion cyclotron frequency fci were observed and are localized in the core of the plasma. Sub-cyclotron modes (f ~ 0.5f ci ) were also observed in L-mode plasmas when the toroidal magnetic field strength was lowered from B T = 2.17 T. In H-mode plasmas, many ICE harmonics localized to the plasma edge were observed, with some exceeding the nominal ICE diagnostic bandwidth of f ϵ [0, 100] MHz. Polarization estimates made using an orthogonal pair of vertical and horizontal magnetic pickup loops on the outer wall of the machine indicate that ICE and sub-cyclotron modes have compressional polarization at the plasma edge, the latter being consistent with simulation efforts for comparable DIII-D plasmas. For all modes, the same harmonics are observed on both the centerpost and outer wall loops, indicating that ICE is poloidally extended. Finally, toroidal mode numbers were calculated using three outer wall loops for both core ICE and sub-cyclotron modes in L-mode plasmas. Furthermore, the sub-cyclotron case served as a benchmark for our calculation method, with measured numbers of n ϵ [–24, –18] roughly agreeing with heuristic estimates of n ϵ [–20, –13]. Core ICE mode numbers were measured to be n ϵ [–10, 5], which is comparable to measurements made on other machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Experimental investigation of beam-ion losses induced by magnetic perturbations using the light ion beam probe technique in the ASDEX Upgrade tokamak

The impact of externally applied magnetic perturbations (MPs) on fast-ion losses has been investigated by means of the light ion beam probe (LIBP) technique in the ASDEX Upgrade tokamak. The LIBP technique allows to experimentally infer the fast-ion orbit displacement induced by MPs via first-orbit losses using scintillator based fast-ion loss detector (FILD) measurements. The fast-ion orbit displacement against different applied MP spectra has been studied. These shots were conducted in ELM mitigated H-mode plasmas. A rigid rotation of the MP coils was applied with a frequency of 1 Hz, with an n = 2 configuration and changing the differential phase between the upper and lower set of coils (ΔΦ ul ) on a shot-to-shot basis. Beam sources Q7 (tangential) and Q8 (radial) were used to probe different fast-ion orbits with FILD1. The measured fast-ion orbit displacement ranges from 3 to 20 mm approximately, and no qualitative difference is observed between ions from beam sources Q7 and Q8. The minimum is found for a ΔΦ ul ~ 50°, which is shifted with respect to the minimum of the plasma boundary displacement, found at ΔΦ ul ~ 0°. A first attempt to validate the orbit following code ASCOT—including the plasma response calculated with the MARS-F code—against these experimental measurements is performed. Finally, while the dependence of the first-orbit fast-ion displacement with ΔΦ ul does not match the experimental measurements, these simulations do capture other features such as the order of magnitude of the orbit displacement and the importance of the toroidal spectrum of the applied perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Regulation of Alfvén Eigenmodes by Microturbulence in Fusion Plasmas

Global gyrokinetic simulations of mesoscale reversed shear Alfven eigenmodes (RSAE) excited by energetic particles (EP) in fusion plasmas find that RSAE amplitude and EP transport are much higher than experimental levels at nonlinear saturation, but quickly diminish to very low levels after the saturation when background microturbulence is artificially suppressed. In contrast, in simulations coupling micro-meso scales, the RSAE amplitude and EP transport decrease drastically at the initial saturation but later increases to the experimental levels in the quasi-steady state with bursty dynamics due to regulation by thermal ion temperature gradient (ITG) microturbulence. The quasi-steady state EP transport is larger for a stronger microturbulence. Here, the RSAE amplitude in the quasi-steady state ITG-RSAE turbulence from gyrokinetic simulations, for the first time, agrees very well with experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

The radial phase variation of reversed-shear and toroidicity-induced Alfvén eigenmodes in DIII-D

The eigenfunction of an instability contains information about energy flow in the wave. Here, the amplitude and phase of electron cyclotron emission radiometer data from hundreds of DIII-D reversed shear Alfvén eigenmodes (RSAE) and toroidicity-induced Alfvén eigenmodes (TAE) are analyzed along the outboard horizontal midplane. The radial phase profile can be flat, linearly rising or falling, convex or concave; in other words, a wide variety of shapes is observed. For a particular mode, often the radial phase profile remains approximately constant as the mode evolves in time but sometimes it changes rapidly. Many TAEs and some RSAEs have phase profiles that are rather flat where the mode amplitude is largest but rise steadily by ~2π at large major radius. Rapid phase changes are observed when the frequencies of an RSAE and TAE overlap and the modes couple. The phase profile depends weakly on the fast-ion gradient that would appear in the absence of wave-induced transport. Linear and quadratic fits to the phase profiles, together with many plasma parameters, are assembled into RSAE and TAE databases. In both cases, large variability is observed. For RSAEs, the strongest phase dependencies are on electron temperature T e , RSAE mode frequency, and the density of carbon impurities. For TAEs, the strongest dependencies are on beam power and major radius of the mode. In general, the average RSAE radial phase profile is essentially flat but the TAE profile has non-zero slope and curvature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Prediction of the energetic particle redistribution by an improved critical gradient model and analysis of the transport threshold

Based on the theory of critical gradient model (CGM) and following the simulation method proposed by Waltz et al., a combination of TGLFEP and EPtran code is employed to predict the energetic particle (EP) transport induced by Alfvén eigenmodes (AEs). To be consistent with the experiment, recent improvements to the simulation method include consideration of threshold evolution and orbit loss due to finite orbit width. The revised CGM is applied to simulate two DIII-D experimental discharges (#142111 and #153071). It well reproduces the experimental profiles with multiple unstable AEs and large-scale EP transport. Discharge #142111 had previously been simulated using a nonlinear MHD-kinetic code MEGA with a transport mechanism based on stochasticity induced by overlapping AE. By comparing the simulated EP profiles, we find that the AE transport threshold is approximated by both the MEGA nonlinear stability threshold and the proposed CGM threshold (error <5% for single n and <17% for multiple n simulation). Both of them are larger than the linear stability threshold of the most unstable AE mode by a quantity of the order of the flux needed to sustain EP transport by the background turbulence. Additionally, we have also applied the improved CGM to simulate the α particle redistribution for a China Fusion Engineering Test Reactor steady state scenario. Because of the clear separation between the AE unstable region and the loss cone, only a moderate α particle loss of ~9.6% is predicted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Orbit tomography of energetic particle distribution functions

Both fast ions and runaway electrons are described by distribution functions, the understanding of which are of critical importance for the success of future fusion devices such as ITER. Typically, energetic particle diagnostics are only sensitive to a limited subsection of the energetic particle phase-space which is often insufficient for model validation. However, previous publications show that multiple measurements of a single spatially localized volume can be used to reconstruct a distribution function of the energetic particle velocity-space by using the diagnostics' velocity-space weight functions, i.e. velocity-space Tomography. In this work we use the recently formulated orbit weight functions to remove the restriction of spatially localized measurements and present orbit tomography, which is used to reconstruct the 3D phase-space distribution of all energetic particle orbits in the plasma. Through a transformation of the orbit distribution, the full energetic particle distribution function can be determined in the standard {energy, pitch, r, z}-space. We benchmark the technique by reconstructing the fast-ion distribution function of an MHD-quiescent DIII-D discharge using synthetic and experimental FIDA measurements. We also use the method to study the redistribution of fast ions during a sawtooth crash at ASDEX upgrade using FIDA measurements. Lastly, a comparison between the orbit tomography and velocity-space tomography is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Fast wave interferometer for ion density measurement on DIII-D

=A fast wave interferometer (FWI), which can measure ion mass density, has been developed on DIII-D for its use on future fusion reactors, as well as for the study of ion behavior in current plasma devices. The frequency of the fast waves used for the FWI is around 60 MHz, and require antennas and coaxial cables or waveguides, which, unlike traditional mirror-based optical interferometers, are less susceptible to neutron/gamma-ray radiation and are relatively immune to impurity deposition and erosion as well as alignment issues. The bulk ion density evaluated using FWI show good agreement with that derived from CO 2 interferometry within about 15%. When the ion mass density measurement by FWI is combined with an electron density measurement from CO 2 interferometry, Z eff measurements are also enabled and are in agreement with those from visible Bremsstrahlung measurements. Additionally, large-bandwidth FWI measurements clearly resolve 10–100 kHz coherent modes and demonstrate its potential as a core fluctuation diagnostic, sensitive to both magnetic and ion density perturbations.

Instruments & Instrumentation↗

DOE OSTI2021

Visualization of Fast Ion Phase-Space Flow Driven by Alfvén Instabilities

Fast ion phase-space flow, driven by Alfven eigenmodes (AEs), is measured by an imaging neutral particle analyzer in the DIII-D tokamak. The flow firstly appears near the minimum safety factor at the injection energy of neutral beams, and then moves radially inward and outward by gaining and losing energy, respectively. The flow trajectories in phase space align well with the intersection lines of the constant magnetic moment surfaces and constant E – (ω/n)P ζ surfaces, where E, P ζ are energy and toroidal canonical momentum of ions; ω and n are angular frequencies and toroidal mode numbers of AEs. It is found that the flow is so destructive that the thermalization of fast ions is no longer observed in regions of strong interaction. Here, the measured phase-space flow is consistent with nonlinear hybrid kinetic-magnetohydrodynamics simulation. Calculations of the relatively narrow phase-space islands reveal that fast ions must transition between different flow trajectories to experience large-scale phase-space transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Testing the DIII-D co/counter off-axis neutral beam injected power and ability to balance injected torque

DIII-D has undergone a major upgrade and successfully injected high power off-axis neutral beams (~4 MW) in both co-current and counter-current directions. This capability of high power co/counter steerable off-axis neutral beams on a major tokamak opens a unique parameter space of broad pressure and current profiles for high beta steady-state advanced tokamak (AT) scenarios, while retaining the ability to balance the injected torque for low rotation studies. This co/counter off-axis neutral beam capability is being used to validate physics-based energetic particle and thermal transport models for designing next-step facilities based on the steady-state AT approach. This paper reports on the critical evaluation of the transmitted power and energetic ion population produced by this heating and current drive system, which is assessed through visible imaging, neutron measurements and rotation profile measurements at balanced torque. Minimal losses of neutral beam power have been achieved by optimizing the strongly focused ion sources required to pass through the aperture. Tilting of the ion source has been guided by fast visible imaging and resulted in neutral beam injection along the design centerline with empirical characterization of each beam's divergence derived from the imaging data and used in the NUBEAM description of the beam injection. Through exclusive power injection of each neutral beam into MHD quiescent plasmas across a range of neutral beam voltage, perveance and plasma current we conclude that a modest reduction (~10%–15%) of transmitted power compared to on-axis, standard focus has been incurred. We report corrections that more accurately represent the injected power. Good ability to balance the neutral beam torque has been demonstrated by injecting the new off-axis counter injecting beam against the existing off-axis co-injected beam in 2.0 T, 1.0 MA, MHD quiescent L-mode plasmas. Furthermore, the torque balance studies verify the ability to operate with balanced injection, which is critical for achieving low torque and low rotation operation for physics studies and in ITER demonstration discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Isotope dependence of beta-induced Alfvén eigenmode (BAE) and low frequency mode (LFM) stability in DIII-D

The stability of beta-induced Alfvén eigenmodes (BAE) and the low frequency modes (LFMs) that were formerly called beta-induced Alfvén-acoustic eigenmodes is discussed. Here, after a brief summary of previous publications on the stability in DIII-D beam-heated, reversed-shear, deuterium plasmas with deuterium neutral beam injection (NBI), new observations in mixed hydrogen and deuterium plasmas are reported. With deuterium NBI, BAEs are at least as unstable in mixed-species plasmas as in deuterium plasmas; however, with hydrogen NBI, the BAEs are stable. In contrast, the LFMs are unaffected by changes in beam species, consistent with the previous observation that LFMs are not driven by high-energy beam ions. As predicted by theory, the LFMs appear more unstable in mixed species plasmas than in pure deuterium discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Propagation of a narrow plasma beam in an oblique magnetic field

The propagation of an intense neutralized ion beam (v is about 5 x 10 exp 8 cm/sec, n is about 10 exp 10/cu cm) through a large insulated vacuum chamber is measured as a function of magnetic field strength and direction. When the beam propagates parallel to the applied field, beam divergence is reduced. When the beam propagates perpendicular to the applied fields, the downstream beam density decreases with increasing field strength. When the beam velocity vector intersects the magnetic field at an oblique angle, beam propagation is determined primarily by the perpendicular component of the field.

Heidbrink, W. W.↗