Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ion temperature”

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 253 records · Page 14

A suite of neutron time-of-flight detectors to measure hot-spot motion in direct-drive inertial confinement fusion experiments on OMEGA

A suite of six neutron time of flight detectors capable of measuring the primary DT neutron energy spectrum has recently been completed on the OMEGA Laser System. The detectors are positioned along multiple quasi-orthogonal lines of sight (LOS’s) in either a stand-alone, collinear, or antipodal configuration. The collinear detector configuration enables one to measure the neutron velocity along the detector LOS independent of the neutron-production history, while the antipodal detector configuration enables one to directly measure of the hot-spot velocity along the detector axis and Gamow velocity shift due to the plasma ion temperature. Using these six measurements of the primary DT neutron energy spectrum, the neutron-averaged hot-spot velocity and Gamow velocity shift of DT fusion neutrons have been measured in cryogenic experiments on OMEGA for the first time. Finally, the detector suite is operable between DT neutron yields of 1 × 10 13 to 2 × 10 14 with an uncertainty of ~14 km/s in the neutron-averaged hot-spot velocity measurement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Tungsten erosion during L-mode discharges in the DIII-D SAS-VW divertor

A series of L-mode plasma discharges was performed in the DIII-D tokamak to assess the impact of outer strike point (OSP) position and toroidal magnetic field direction on erosion and core contamination potential of the recently-installed, tungsten-coated Small Angle Slot (SAS-VW) divertor. In one discharge, in-slot emission spectroscopy measured an<48 % increase in the W gross erosion rate when the OSP was moved 3 cm outwards, away from the V-shaped vertex of the slot divertor. However, the effective W yield (erosion rate divided by the incident D flux) was, overall, insensitive to changes in OSP location. Consistently low estimates of the effective W yield based on measurements taken a few cm outwards from the vertex suggest potentially significant C surface contamination. No W emission signal was detected when orienting the toroidal magnetic field such that the ion B×∇B drift direction is pointed away from the X-point. However, measurements of W content in the plasma core for both toroidal magnetic field directions suggest the presence of additional, unmeasured sources of erosion. The difference in the measured core W density with OSP position is much greater than the difference in the measured erosion rates, which may suggest that the leakage of eroded impurities out of the divertor is governed primarily through the parallel ion temperature gradient and friction forces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advancing the physics basis for quasi-helically symmetric stellarators

A new optimized quasi-helically symmetric configuration is described that has the desirable properties of improved energetic particle confinement, reduced turbulent transport by three-dimensional shaping and non-resonant divertor capabilities. The configuration presented in this paper is explicitly optimized for quasi-helical symmetry, energetic particle confinement, neoclassical confinement and stability near the axis. Post optimization, the configuration was evaluated for its performance with regard to energetic particle transport, ideal magnetohydrodynamic stability at various values of plasma pressure and ion temperature gradient instability induced turbulent transport. The effects of discrete coils on various confinement figures of merit, including energetic particle confinement, are determined by generating single-filament coils for the configuration. Preliminary divertor analysis shows that coils can be created that do not interfere with expansion of the vessel volume near the regions of outgoing heat flux, thus demonstrating the possibility of operating a non-resonant divertor.

Physics↗

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ion temperature gradient (ITG), trapped electron mode (TEM) and electron temperature gradient (ETG) modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. Finally, it is also confirmed that thermodiffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic-ballooning-mode turbulence in low-average-magnetic-shear equilibria

In this work, kinetic-ballooning-mode (KBM) turbulence is studied via gyrokinetic flux-tube simulations in three magnetic equilibria that exhibit small average magnetic shear: the Helically Symmetric eXperiment (HSX), the helical-axis Heliotron-J and a circular tokamak geometry. For HSX, the onset of KBM being the dominant instability at low wavenumber occurs at a critical value of normalized plasma pressure $\beta ^{\rm KBM}_{\rm crit}$ that is an order of magnitude smaller than the magnetohydrodynamic (MHD) ballooning limit $\beta ^{\rm MHD}_{\rm crit}$ when a strong ion temperature gradient (ITG) is present. However, $\beta ^{\rm KBM}_{\rm crit}$ increases and approaches the MHD ballooning limit as the ITG tends to zero. For these configurations, $\beta ^{\rm KBM}_{\rm crit}$ also increases as the magnitude of the average magnetic shear increases, regardless of the sign of the normalized magnetic shear. Simulations of Heliotron-J and a circular axisymmetric geometry display behaviour similar to HSX with respect to $\beta ^{\rm KBM}_{\rm crit}$ . Despite large KBM growth rates at long wavelengths in HSX, saturation of KBM turbulence with $\beta > \beta _{\rm crit}^{\rm KBM}$ is achievable in HSX and results in lower heat transport relative to the electrostatic limit by a factor of roughly five. Nonlinear simulations also show that KBM transport dominates the dynamics when KBMs are destabilized linearly, even if KBM growth rates are subdominant to ITG growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phase contrast imaging measurements and numerical simulations of turbulent density fluctuations in gas-fuelled ECRH discharges in Wendelstein 7-X

The fundamental nature of turbulent density fluctuations in standard Wendelstein 7-X (W7-X) stellarator discharges is investigated experimentally via phase contrast imaging (PCI) in combination with gyrokinetic simulations with the code GENE. We find that density fluctuations are ion-temperature-gradient-driven and radially localised in the outer half of the plasma. It is shown that the line-integrated PCI measurements cover the right range of wavenumbers and a favourable toroidal and poloidal location to capture some of the strongest density fluctuations in W7-X. Due to the radial localisation of fluctuations, measured wavenumber–frequency spectra exhibit a dominant phase velocity, which can be related to the E × B rotation velocity at the radial position of a well in the neoclassical radial electric field. The match is robust against variations of heating power and line-integrated density, which is partly due to the localisation of fluctuations and partly due to effects of the radial gradient in the E × B velocity profile on the wavenumber–frequency spectrum. The latter effect is studied with a newly built synthetic PCI diagnostic and global gyrokinetic simulations with GENE-3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ultrafast multi-cycle terahertz measurements of the electrical conductivity in strongly excited solids

Key insights in materials at extreme temperatures and pressures can be gained by accurate measurements that determine the electrical conductivity. Free-electron laser pulses can ionize and excite matter out of equilibrium on femtosecond time scales, modifying the electronic and ionic structures and enhancing electronic scattering properties. The transient evolution of the conductivity manifests the energy coupling from high temperature electrons to low temperature ions. Here we combine accelerator-based, high-brightness multi-cycle terahertz radiation with a single-shot electro-optic sampling technique to probe the evolution of DC electrical conductivity using terahertz transmission measurements on sub-picosecond time scales with a multi-undulator free electron laser. Our results allow the direct determination of the electron-electron and electron-ion scattering frequencies that are the major contributors of the electrical resistivity.

36 MATERIALS SCIENCE↗

Nonlinear global gyrokinetic delta- f turbulence simulations in a quasi-axisymmetric stellarator

We use the global gyrokinetic stellarator code XGC-S to look at turbulence physics in geometry of the proposed quasi-axisymmetric stellarator, QUASAR. In this article, the XGC-S code is first verified for nonlinear simulations of ion temperature gradient-driven turbulence by a successful benchmark with the established nonlinear tokamak code XGC1 in a CYCLONE-like tokamak case. We then apply the XGC-S code to QUASAR for nonlinear turbulence physics. We find that the highly poloidally localized linear mode gives way to a more evenly distributed nonlinear turbulent structure. Furthermore, the radial extent of the mode increases and occupies a large fraction of the volume, despite the initially localized temperature profile, as the profiles relax.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dependence of the impurity transport on the dominant turbulent regime in ELM-y H-mode discharges on the DIII-D tokamak

Laser blow-off injections of aluminum and tungsten have been performed on the DIII-D tokamak to investigate the variation of impurity transport in a set of dedicated ion and electron heating scans with a fixed value of the external torque. The particle transport is quantified via the Bayesian inference method, which, constrained by a combination of a charge exchange recombination spectroscopy, soft X-ray measurements, and VUV spectroscopy provides a detailed uncertainty quantification of transport coefficients. Contrasting discharge phases with a dominant electron and ion heating reveal 3-fold drop in the impurity confinement time and order of magnitude increase in midradius impurity diffusion, when additional electron heating is applied. Further, the calculated stationary aluminum density profiles reverse from peaked in electron heated to hollow in the ion heated case, following a similar trend as electron and carbon density. Comparable values of a core diffusion have been observed for W and Al ions, while differences in the propagation dynamics of these impurities are attributed to pedestal and edge transport. Modeling of the core transport with non-linear gyrokinetics code CGYRO, significantly underpredicts the magnitude of the variation in Al transport. Diffusion increases 3-times steeper with additional electron heat flux and 10-times lower diffusion is observed in ion heated case than predicted by the modeling. Here, the CGYRO model quantitatively matches the increase of the Al diffusion when approaching the linear threshold for the transition from the ion temperature gradient (ITG) to trapped electron mode (TEM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced pedestal H-mode at low edge ion collisionality on NSTX

The Enhanced Pedestal (EP) H-mode regime is an attractive wide-pedestal high-βp scenario for the National Spherical Torus Experiment Upgrade (NSTX-U) and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (βN > 5), and significant bootstrap fraction (fBS > 0.6) at Ip/BT = 2 MA/T. This regime is realized when the edge ion collisionality becomes sufficiently small that a positive feedback interaction occurs between a reduction in the ion neoclassical energy transport and an increase in the particle transport from pressure-driven edge instabilities. The EP H-mode was most often observed as a transition following a large edge-localized mode in conditions with low edge neutral recycling. It is hypothesized that the onset of pressure-driven instabilities prior to the full recovery of the neutral density leads to a temporary period with an elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy and magnetic spectroscopy in order to describe the evolution of the edge particle transport mechanisms during the ELM recovery and the saturated EP H-mode state. The observations are consistent with the hypothesis that the onset of pressure-driven edge instabilities, such as the kinetic ballooning mode and kink-peeling mode, can be responsible for the increased particle transport in EP H-mode.

Battaglia, D. J. (ORCID:0000000188979740)↗

Global gyrokinetic particle simulations of microturbulence in W7-X and LHD stellarators

Global gyrokinetic particle simulations of electrostatic ion temperature gradient (ITG) instability show that the most unstable eigenmode is localized to some magnetic fieldlines or discrete locations on the poloidal plane in the Wendelstein 7-X (W7-X) stellarator due to its mirror-like magnetic fields, which vary strongly in the toroidal direction and induce coupling of more toroidal harmonics (n) to form the linear eigenmode than in the Large Helical Device (LHD) stellarator. Nonlinear electrostatic simulation results show that self-generated zonal flows are the dominant saturation mechanism for the ITG instabilities in both the LHD and W7-X. Additionally, radial widths of the fluctuation intensity in both the LHD and W7-X are significantly broadened from the linear phase to the nonlinear phase due to turbulence spreading. Finally, nonlinear spectra in the W7-X are dominated by low-n harmonics, which can be generated both by nonlinear toroidal coupling of high-n harmonics and by linear toroidal coupling with large amplitude zonal flows due to the 3D equilibrium magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of local electrostatic gyrokinetic simulation of driftwave instability in field-reversed configuration

Local linear simulations of ion temperature gradient (ITG) instabilities using a gyrokinetic particle code GKD1D have been verified by comparing with the analytic dispersion relation in the slab geometry and by benchmarking with a global gyrokinetic code GTC-X in a realistic field-reversed configuration (FRC). Results of ITG simulations using either adiabatic or kinetic electrons exhibit excellent agreements between the two codes for linear frequencies, growth rates, and mode structures. The ITG modes in the FRC scrape-off layer region are mostly dominated by the poloidal m = 1 harmonic. Kinetic electrons enhance the subdominant m = 2 harmonic and the even parity dominates the ITG mode structure. Kinetic electrons increase the growth rates for the low-n modes, but decrease the growth rates for the high-n modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deficiencies in compression and yield in x-ray-driven implosions

This paper analyzes x-ray–driven implosions that are designed to be less sensitive to 2-D and 3-D effects in hohlraum and capsule physics. Key performance metrics including the burn-averaged ion temperature, hot-spot areal density, and fusion yield are found to agree with simulations where the design adiabat (internal pressure) is multiplied by a factor of 1.4. Furthermore, these results motivate the development of a simple model for interpreting experimental data, which is then used to quantify how improvements in compression could help achieve ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First coupled GENE–XGC microturbulence simulations

Covering the core and the edge region of a tokamak, respectively, the two gyrokinetic turbulence codes Gyrokinetic Electromagnetic Numerical Experiment (GENE) and X-point Gyrokinetic Code (XGC) have been successfully coupled by exchanging three-dimensional charge density data needed to solve the gyrokinetic Poisson equation over the entire spatial domain. Certain challenges for the coupling procedure arise from the fact that the two codes employ completely different numerical methods. This includes, in particular, the necessity to introduce mapping procedures for the transfer of data between the unstructured triangular mesh of XGC and the logically rectangular grid (in a combination of real and Fourier space) used by GENE. Constraints on the coupling scheme are also imposed by the use of different time integrators. First, coupled simulations are presented. We have considered collisionless ion temperature gradient turbulence, in both circular and fully shaped plasmas. Coupled simulations successfully reproduce both GENE and XGC reference results, confirming the validity of the code coupling approach toward a whole device model. Here, many lessons learned in the present context, in particular, the need for a coupling procedure as flexible as possible, should be valuable to our and other efforts to couple different kinds of codes in pursuit of a more comprehensive description of complex real-world systems and will drive our further developments of a whole device model for fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Suppression of ITG turbulence due to spectral shift during biasing induced H-mode on HBT-EP

Investigations of biasing induced H-mode discharges on HBT-EP show that the edge turbulence is consistent with the ion temperature gradient) mode and have allowed for the controlled observation of the effect of applied flow shear on the turbulence. Measurements of the radial wavenumber spectrum of floating potentials at the edge show that the turbulence intensity decreases with increasing shift in the spectrum average &#x27E8; k r &#x27E9; when increasing amounts of bias probe voltage (and increasing amounts of flow shear) is applied. This is in agreement with the spectral shift model [Staebler et al., Phys. Rev. Lett. 110, 055003 2013] for turbulence suppression via sheared flow. A shift in the wavenumber spectrum occurs at applied electrode voltages and currents below the threshold needed for an L–H transition, and a dithering transition is obtained when biasing near the threshold. Suppression of blob-filament turbulence in the scrape-off layer (SOL) precedes the L–H transition, with the SOL turbulence remaining low throughout the dithering phase, despite the modulation of turbulence levels in the nearby edge. This demonstrates that the SOL turbulence “decouples” from the edge turbulence. The spectral shift in the measured radial wavenumber is corroborated by the direct measurement of eddy tilt angle using a novel time delay analysis technique first developed for Doppler reflectometry [Pinzon et al., Plasma Phys. Controlled Fusion 61, 105009 (2019)] but adapted here for floating potential measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗