Hager, Robert; Ku, Seung-Hoe; Sharma, Amil Y.; Churchill, Randy Michael; Chang, C. S.; Scheinberg, Aaron
The simplified delta-f mixed-variable/pull-back electromagnetic simulation algorithm implemented in XGC for core plasma simulations by Cole et al. [Phys. Plasmas 28, 034501 (2021)] has been generalized to a total-f electromagnetic algorithm that can include, for the first time, the boundary plasma in diverted magnetic geometry with neutral particle recycling, turbulence and neoclassical physics.
The delta-f mixed-variable/pull-back electromagnetic implementation is based on the pioneering work by Kleiber and Mischenko et al. [Kleiber et al., Phys. Plasmas 23, 032501 (2016); Mishchenko et al., Comput. Phys. Commun. 238, 194 (2019)].
An electromagnetic demonstration simulation is performed in a DIII-D-like, H-mode boundary plasma, including a corresponding comparative electrostatic simulation, which confirms that the electromagnetic simulation is necessary for a higher fidelity understanding of the electron particle and heat transport even at the low-beta pedestal foot in the vicinity of the magnetic separatrix.
The dynamic interplay between the core and the edge plasma has important consequences in the confinement and heating of fusion plasma. The transport of the Scrape-Off-Layer (SOL) plasma imposes boundary conditions on the core plasma, and neutral transport through the SOL influences the core plasma sourcing. In order to better study these effects in a self-consistent, time-dependent fashion with reasonable turn-around time, a reduced model is needed. In this paper we introduce the SOL Box Model, a reduced SOL model that calculates the plasma temperature and density in the SOL given the core-to-edge particle and power fluxes and recycling coefficients. The analytic nature of the Box Model allows one to readily incorporate SOL physics in time-dependent transport solvers for pulse design applications in the control room. Here we demonstrate such a coupling with the core transport solver TRANSP and compare the results with density and temperature measurements, obtained through Thomson scattering and Langmuir probes, of an NSTX discharge. Implications for future interpretive and predictive simulations are discussed.
Notterman, Daniel A; Schneper, Lisa M; Drake, Amanda; Piyasena, Chinthika
This entry contains the data used in the PLOS ONE publication entitled, "Characteristics of salivary telomere length shortening in preterm infants" by Schneper et al. The objective of the study was to examine the association between gestational age, telomere length (TL) and rate of shortening in newborns. Genomic DNA was isolated from buccal samples of 39 term infants at birth and one year and 32 preterm infants at birth, term-adjusted age (40 weeks post-conception) and age one-year corrected for gestational duration. Telomere length was measured by quantitative real-time PCR. Demographic and clinical data were collected during clinic or research visits and from hospital records. Socioeconomic status was estimated using the deprivation category (DEPCAT) scores derived from the Carstairs score of the subject's postal code.
Schwartz, Jacob A.; Ricks, Wilson; Kolemen, Egemen; Jenkins, Jesse D.
Fusion could be a part of future decarbonized electricity systems, but it will need to compete with other technologies.
In particular, pulsed tokamaks plants have a unique operational mode, and evaluating
which characteristics make them economically competitive can help select between design pathways.
Using a capacity expansion and operations model,
we determined cost thresholds for pulsed tokamaks to reach a range of penetration levels in a future decarbonized US Eastern Interconnection.
The required capital cost to reach a fusion capacity of 100 GW varied from $3000 to $7200/kW,
and the equilibrium penetration increases rapidly with decreasing cost.
The value per unit power capacity depends on the variable operational cost and on cost of its competition, particularly fission, much more than on the pulse cycle parameters.
These findings can therefore provide initial cost targets for fusion more generally in the United States.
Geyman, Emily C.; Wu, Ziman; Nadeau, Matthew D.; Edmonsond, Stacey; Turner, Andrew; Purkis, Sam J.; Howes, Bolton; Dyer, Blake; Ahm, Anne-Sofie C.; Yao, Nan; Deutsch, Curtis A.; Higgins, John A.; Stolper, Daniel A.; Maloof, Adam C.
Carbonate mud represents one of the most important geochemical archives for reconstructing ancient climatic, environmental, and evolutionary change from the rock record. Mud also represents a major sink in the global carbon cycle. Yet, there remains no consensus about how and where carbonate mud is formed. In this contribution, we present new geochemical data that bear on this problem, including stable isotope and minor and trace element data from carbonate sources in the modern Bahamas such as ooids, corals, foraminifera, and green algae.
Data set corresponding to "NAPS: Integrating pose estimation and tag-based tracking." This dataset contains the corresponding videos, tracking scripts, and SLEAP models along with SLEAP, NAPS, and ArUco tracking results.
The usage of permanent magnets to shape the confining field of a stellarator has the potential to reduce or eliminate the need for non-planar coils. As a proof-of-concept for this idea, we have developed a procedure for designing an array of cubic permanent magnets that works in tandem with a set of toroidal-field coils to confine a stellarator plasma. All of the magnets in the design are constrained to have identical geometry and one of three polarization types in order to simplify fabrication while still producing sufficient field accuracy. We present some of the key steps leading to the design, including the geometric arrangement of the magnets around the device, the procedure for optimizing the polarizations according to three allowable magnet types, and the choice of magnet types to be used. We apply these methods to design an array of rare-Earth permanent magnets that can be paired with a set of planar toroidal-field coils to confine a quasi-axisymmetric plasma with a toroidal magnetic field strength of about 0.5 T on axis.
Zhu, Hongxuan; Stoltzfus-Dueck, T; Hager, R; Ku, S; Chang, C. S.
Ion orbit loss is considered important for generating the radially inward electric field Er in a tokamak edge plasma. In particular, this effect is emphasized in diverted tokamaks with a magnetic X point. In neoclassical equilibria, Coulomb collisions can scatter ions onto loss orbits and generate a radially outward current, which in steady state is balanced by the radially inward current from viscosity. To quantitatively measure this loss-orbit current in an edge pedestal, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. As the first application of this diagnostic, a neoclassical DIII-D H-mode plasma is studied using gyrokinetic ions and adiabatic electrons. The validity of the diagnostic is demonstrated by studying the collisional relaxation of Er in the core. After this demonstration, the loss-orbit current is numerically measured in the edge pedestal in quasisteady state. In this plasma, it is found that the radial electric force on ions from Er approximately balances the ion radial pressure gradient in the edge pedestal, with the radial force from the plasma flow term being a minor component. The effect of orbit loss on Er is found to be only mild.
An important goal of stellarator optimization is to achieve good confinement of
energetic particles such as, in the case of a reactor, alphas created by Deuterium-Tritium
(D-T) fusion. In this work, a fixed-boundary stellarator equilibrium was re-optimized for
energetic particle confinement via a two-step process: first, by minimizing deviations from quasi-axisymmetry (QA) on a single flux surface near the mid-radius, and secondly by maintaining
this improved quasi-axisymmetry while minimizing the analytical quantity ΓC , which represents
the angle between magnetic flux surfaces and contours of J||, the second adiabatic invariant.
This was performed multiple times, resulting in a group of equilibria with significantly reduced
energetic particle losses, as evaluated by Monte Carlo simulations of alpha particles in scaled-up
versions of the equilibria. This is the first time that energetic particle losses in a QA stellarator
have successfully been reduced by optimizing ΓC . The relationship between energetic particle
losses and metrics such as QA error (Eqa) and ΓC in this set of equilibria were examined via
statistical methods and a nearly linear relationship between volume-averaged ΓC and prompt
particle losses was found.
Recent U.S. fusion development strategy reports all recommend that the U.S. should pursue innovative science and technology to enable construction of a Fusion Pilot Plant (FPP) that produces net electricity from fusion at low capital cost. Compact tokamaks have been proposed as a means of potentially reducing the capital cost of a fusion pilot plant. However, compact steady-state tokamak FPPs face the challenge of integrating a high fraction of self-driven current with high core confinement, plasma pressure, and high divertor parallel heat flux. This integration is sufficiently challenging that a dedicated sustained-high-power-density (SHPD) tokamak facility is proposed by the U.S. community as the optimal way to close this integration gap. Performance projections for the steady-state tokamak FPP regime are presented and a preliminary SHPD device with substantial flexibility in lower aspect ratio (A=2-2.5), shaping, and divertor configuration to narrow gaps to a FPP is described.
Sharma, A. Y.; Cole, M. D. J.; Görler, T.; Chen, Y.; Hatch, D. R.; Guttenfelder, W.; Hager, R.; Sturdevant, B. J.; Ku, S.; Chang, C. S.
Plasma shaping may have a stronger effect on global turbulence in tight-aspect-ratio tokamaks than in conventional-aspect-ratio tokamaks due to the higher toroidicity and more acute poloidal asymmetry in the magnetic field. In addition, previous local gyrokinetic studies have shown that it is necessary to include parallel magnetic field perturbations in order to accurately compute growth rates of electromagnetic modes in tight-aspect-ratio tokamaks. In this work, the effects of elongation and triangularity on global, ion-scale, linear electromagnetic modes are studied at NSTX aspect ratio and high plasma beta using the global gyrokinetic particle-in-cell code XGC. The effects of compressional magnetic perturbations are approximated via a well-known modification to the particle drifts that was developed for flux-tube simulations [N. Joiner et al., Phys. Plasmas 17, 072104 (2010)], without proof of its validity in a global simulation. Magnetic equilibria are re-constructed for each distinct plasma profile that is used. Coulomb collision effects are not considered. Within the limitations imposed by the present study, it is found that linear growth rates of electromagnetic modes (collisionless microtearing modes and kinetic ballooning modes) are significantly reduced by NSTX-like shaping. For example, growth rates of kinetic ballooning modes at high beta are reduced to the level of that of collisionless trapped electron modes.
Stoltzfus-Dueck, T; Hornsby, W A; Grosshauser, S R
Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.
Wang, Yin; Gilson, Erik P.; Ebrahimi, Fatima; Goodman, Jeremy; Caspary, Kyle J.; Winarto, Himawan W.; Ji, Hantao
This dataset provides the source data of figures in the main text of the paper "Identification of a non-axisymmetric mode in laboratory experiments searching for standard magnetorotational instability" accepted by Nature Communications.