Search NASA⌕ Search

SEARCH · Search NASA

Results for “neoclassical transport”

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 73 records · Page 4

Quasisymmetric Optimization of Nonaxisymmetry in Tokamaks

Predictive 3D optimization reveals a novel approach to modify a nonaxisymmetric magnetic perturbation to be entirely harmless for tokamaks, by essentially restoring quasisymmetry in perturbed particle orbits as much as possible. Such a quasisymmetric magnetic perturbation (QSMP) has been designed and successfully tested in the KSTAR and DIII-D tokamaks, demonstrating no performance degradation despite the large overall amplitudes of nonaxisymmetric fields and strong response otherwise expected in the tested plasmas. The results indicate that a quasisymmetric optimization is a robust path of error field correction across the resonant and nonresonant field spectrum in a tokamak, leveraging the prevailing concept of quasisymmetry for general 3D plasma confinement systems such as stellarators. Finally, the optimization becomes, in fact, a simple eigenvalue problem to the so-called torque response matrices if a perturbed equilibrium is calculated consistent with nonaxisymmetric neoclassical transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics Basis, Optimization, and Control for Integrated 3D Edge Long-pulse Tokamak Scenarios

The aim of this project with international collaborations (PI: J.-K. Park of Princeton Plasma Physics Laboratory) is to develop a unified physics basis and predictive capability for the control of edge-localized modes (ELMs) with optimized non-axisymmetric (3D) fields, by leveraging the unique research capabilities of international KSTAR tokamak in Korea. At the University of California, Irvine, we use a fusion simulation code GTC to study the effects of the 3D magnetic fields (including both magnetic islands and stochastic magnetic field lines) on microturbulence. We focus on the effects of 3D equilibrium on turbulent and neoclassical transport that leads to ELM suppression. During the project period, the UCI team of this collaboration has made significant progress in understanding turbulent transport in 3D geometry including both tokamaks with RMP and stellarators. Four UCI postdoctoral researchers, Drs. Javier Nicolau, Gyungjin Choi, Pengfei Liu, and Xishuo Wei were partially supported by this project and collaborating with the entire project team. We have published 4 journal papers and gave several invited talks. This project advance physics understand of 3D effects on microturbulence, specially effects of magnetic islands. It will contribute to the DOE SciDAC ISEP project in the area of GTC development for general 3D geometry capability, rigorous verification and validation, and readiness for the emerging exascale computing. Through collaborations with computational scientists, the project helps to advance high performance computing in the area of GPU acceleration, linear solver, and data management.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of the neoclassical ambipolar electric field in ion-root plasmas on W7-X

The role of the radial electric field in high-performance ion-root plasmas on Wendelstein 7-X (W7-X) is examined and compared with neoclassical predictions herein. The W7-X stellarator is the world’s first large-scale optimized stellarator. One of the important targets chosen for optimization during the W7-X design process was the reduction of core neoclassical heat transport. This optimization was targeted for reactor-relevant high-density plasmas with $T_e ≈ T_i$ in which the neoclassical ambipolar radial electric field is expected to be negative throughout the plasmas core. Measurements of the core radial electric field $(E_r)$ have confirmed that ion-root conditions (negative $(E_r)$ in the plasma core) have been achieved in W7-X with high-density plasmas and central ECRH. These measured $(E_r)$ profiles agree well with the neoclassical ambipolar $(E_r)$ predicted by the code SFINCS . This good agreement provides confidence in the validity of neoclassical calculations in high-density ion-root conditions, and enables initial studies on the role of neoclassical transport in the optimized high-density regime of W7-X. Profile measurements of electron temperature $(T_e)$, ion temperature $(T_i)$ and electron density $(n_e)$ along with approximations for the average value of $Z$ eff have been used as inputs to the sfincs code to calculate the ambipolar $E_r$ profile along with neoclassical ion and electron energy flux profiles $(Q$ NCi , $Q$ NCe $)$. Finally the total experimental energy input to the plasma from ECRH heating has been compared to the neoclassical energy fluxes to provide a first estimate for the fraction of transport that can be attributed to neoclassical processes in reactorrelevant high-density ion-root conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Charge-state independent anomalous transport for a wide range of different impurity species observed at Wendelstein 7-X

Here, the plasma volume averaged impurity confinement of selected charge states and impurity species has been characterized for the Stellarator Wendelstein 7-X (W7-X), covering a wide range of atomic charges (Z = 12–44) and atomic masses (M = 28–184). A comparison of the experimental findings to theoretical neoclassical and turbulent transport expectations suggests, aside from/in addition to the neoclassical transport, an additional significant anomalous transport mechanism, which is not inconsistent with the predictions of a turbulence dominated impurity transport and is in agreement with the experimental results from recent transport studies based on the direct measurements of impurity diffusion profiles, performed at W7-X.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear ion-scale microstability analysis of high and low-collisionality NSTX discharges and NSTX-U projections

Linear gyrokinetic simulations were conducted to investigate ion-gyroradius-scale micro-instability predictions for high-beta NSTX discharges and NSTX-U projections that span over an order of magnitude variation in collisionality. A complex mix of microtearing modes and hybrid trapped electron modes/kinetic ballooning modes (TEM/KBM) is predicted for all experimental or projected conditions. Ion temperature gradient (ITG) instabilities are typically stable in the NSTX discharges investigated, consistent with the observed neoclassical ion thermal transport. ITG thresholds inferred from the simulations are typically much higher than the experimental NSTX gradients, as well as the projected gradients in the NSTX-U scenario, which assumed ion temperatures limited by neoclassical transport only. The analysis suggests ITG instabilities are unlikely to contribute significant anomalous thermal losses in high-beta, lower collisionality NSTX-U scenarios. On the other hand, the NSTX experimental profiles and NSTX-U projections are predicted to be very close to the predicted onset of unstable KBM at most radii investigated. The proximity of the various discharges to the KBM instability threshold implies it may play an important role in setting profile shapes and limiting global energy confinement. Here, it remains to be understood and predicted how KBM contributes to multi-channel transport (thermal and particle transport, for both ions and electrons) in a way that is consistent with experimental inferences.

36 MATERIALS SCIENCE↗

MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients in stellarator plasmas

Abstract MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES/ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quasilinear theory and modelling of gyrokinetic turbulent transport in tokamaks

The theory, development, and validation of reduced quasilinear models of gyrokinetic turbulent transport in the closed flux surface core of tokamaks is reviewed. In combination with neoclassical collisional transport, these models are successful in accurately predicting core tokamak plasma temperature, density, rotation, and impurity profiles in a variety of confinement regimes. Refined experimental tests have been performed to validate the predictions of the quasilinear models, probing changes in the dominant gyrokinetic instabilities, as reflected in fluctuation measurements, cross-phases, and transport properties. These tests continue to produce a deeper understanding of the complex mix of instabilities at both electron and ion gyroradius scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impurity temperature screening in stellarators close to quasisymmetry

Impurity temperature screening is a favourable neoclassical phenomenon involving an outward radial flux of impurity ions from the core of fusion devices. Quasisymmetric magnetic fields lead to intrinsically ambipolar neoclassical fluxes that give rise to temperature screening for low enough $\unicode[STIX]{x1D702}^{-1}\equiv d\ln n/d\ln T$ . In contrast, neoclassical fluxes in generic stellarators will depend on the radial electric field, which is predicted to be inward for ion-root plasmas, potentially leading to impurity accumulation. Here, we examine the impurity particle flux in a number of approximately quasisymmetric stellarator configurations and parameter regimes while varying the amount of symmetry breaking in the magnetic field. For the majority of this work, neoclassical fluxes have been obtained using the SFINCS drift-kinetic equation solver with electrostatic potential $\unicode[STIX]{x1D6F7}=\unicode[STIX]{x1D6F7}(r)$ , where $r$ is a flux-surface label. Results indicate that achieving temperature screening is possible, but unlikely, at low-collisionality reactor-relevant conditions in the core. Thus, the small departures from symmetry in nominally quasisymmetric stellarators are large enough to significantly alter the neoclassical impurity transport. A further look at the magnitude of these fluxes when compared to a gyro-Bohm turbulence estimate suggests that neoclassical fluxes are small in configurations optimized for quasisymmetry when compared to turbulent fluxes.

Physics↗

Elements of H-mode pedestal structure

Abstract This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Soft X-ray measurements of transport and MHD activity in the core and edge NSTX plasma (Final Technical Report)

A comprehensive suite of innovative diagnostic instrumentation and techniques has been developed by the Johns Hopkins collaboration for measurements of transport and MHD activity in the core and edge of the National Spherical Torus Experiment (NSTX) at the Princeton Plasma Physics Laboratory (PPPL). These include ultrafast multi-energy poloidal soft X-ray arrays, high spatial resolution multi-energy tangential soft X-ray arrays, high spectral resolution transmission grating imaging VUV spectrometers and ultrafast bolometric arrays, which measure over a broad range of spatial, spectral, and temporal scales the emission from ohmic, RF, and beam driven NSTX plasmas. The instrumentation was designed to fulfill the specific measurement needs of the Spherical Torus (ST) magnetic fusion concept and is accompanied by a comprehensive set of numerical tools for data analysis and modeling, aimed at determining as many as possible plasma parameters and at extracting as much as possible physics information from NSTX experiments. The Johns Hopkins diagnostics have been and continue to be instrumental for essentially all the NSTX experiments. In addition, our own research and experiments using these diagnostics has led over the years to fundamental conclusions for the ST fusion physics, such as the establishment of near neoclassical impurity transport in the NSTX core and edge plasma, and the critical role of the Global Alfven Eigenmode (GAE) MHD activity for the ST core electron thermal transport. The Johns Hopkins Diagnostics continue to be operational on the upgraded NSTX-U and will provide key plasma information in its upcoming high field, current and power experimental campaigns. The present report summarizes the main technical characteristics and measurement capabilities of the diagnostic instrumentation and methods developed within our research and describes a key physics result achieved during this research. Success in our research was made possible also by the excellent support and cooperation of the NSTX research team and of PPPL researchers, who are warmly acknowledged here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic fields with general omnigenity

Omnigenity is a desirable property of toroidal magnetic fields that ensures confinement of trapped particles. Confining charged particles is a basic requirement for any fusion power plant design, but it can be difficult to satisfy with the non-axisymmetric magnetic fields used by the stellarator approach. Every ideal magnetohydrodynamic equilibrium previously found to approximate omnigenity has been either axisymmetric, quasi-symmetric or has poloidally closed contours of magnetic field strength $B$ . However, general omnigenous equilibria are a much larger design space than these subsets. A new model is presented and employed in the DESC stellarator optimization suite to represent and discover the full parameter space of omnigenous equilibria. Although exact omnigenity aside from quasi-symmetry is impossible, these results reveal that excellent particle confinement can be achieved in practice. Examples far from quasi-symmetry with poloidally, helically and toroidally closed $B$ contours are attained with DESC and shown to have low neoclassical collisional transport and fast particle losses.

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)↗

Confinement in electron heated plasmas in Wendelstein 7-X and ASDEX Upgrade; the necessity to control turbulent transport

In electron (cyclotron) heated plasmas, in both ASDEX Upgrade (L-mode) and Wendelstein 7-X, clamping of the ion temperature occurs at $T_\mathrm{i}$ ~ 1.5 keV independent of magnetic configuration. The ions in such plasmas are heated through the energy exchange power as ${n}_{\mathrm{e}}^{2}({T}_{\mathrm{e}}-{T}_{\mathrm{i}})/{T}_{\mathrm{e}}^{3/2}$, which offers a broad ion heating profile, similar to that offered by alpha heating in future thermonuclear fusion reactors. However, the predominant electron heating may put an additional constraint on the ion heat transport, as the ratio $T_\mathrm{e}T_\mathrm{i}$ > 1 can exacerbates ITG/TEM core turbulence. Therefore, in practical terms the strongly 'stiff' core transport translates into Ti-clamping in electron heated plasmas. Due to this clamping, electron heated L-mode scenarios, with standard gas fueling, in either tokamaks or stellarators may struggle to reach high normalized ion temperature gradients required in a compact fusion reactor. The comparison shows that core heat transport in neoclassically optimized stellarators is driven by the same mechanisms as in tokamaks. The absence of a strong H-mode temperature edge pedestal in stellarators, sofar (which, like in tokamaks, could lift the clamped temperature-gradients in the core), puts a strong requirement on reliable and sustainable core turbulence suppression techniques in stellarators.

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 ELM-free high-betap scenario for NSTX-U and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (betaN > 5) and significant bootstrap fraction (f_BS > 0.6) at I_p/B_T = 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. EP H-mode was most often observed as a transition following a large ELM 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 elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy (BES) 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 KBM and kink-peeling, can be responsible for the increased particle transport in EP H-mode.

edge and boundary physics↗

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement.The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking into account the linear two-fluid plasma response.With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate, and suppression of the electron heat flux in the steepest part of the edge pedestal.In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.

DIII-D↗