Search

Start Over Domain Natural Sciences

91. A software package for plasma facing component analysis and design: the Heat flux Engineering Analysis Toolkit (HEAT)

Author(s):
Looby, Tom; Reinke, Matthew; Wingen, Andreas; Menard, Jonathan; Gerhardt, Stefan; Gray, Travis; Donovan, David; Unterberg, Ezekial; Klabacha, Jonathan; Messineo, Mike
Abstract:
The engineering limits of plasma facing components (PFCs) constrain the allowable operational space of tokamaks. Poorly managed heat fluxes that push the PFCs beyond their limits not only degrade core plasma performance via elevated impurities, but can also result in PFC failure due to thermal stresses or melting. Simple axisymmetric assumptions fail to capture the complex interaction between 3D PFC geometry and 2D or 3D plasmas. This results in fusion systems that must either operate with increased risk or reduce PFC loads, potentially through lower core plasma performance, to maintain a nominal safety factor. High precision 3D heat flux predictions are necessary to accurately ascertain the state of a PFC given the evolution of the magnetic equilibrium. A new code, the Heat flux Engineering Analysis Toolkit (HEAT), has been developed to provide high precision 3D predictions and analysis for PFCs. HEAT couples many otherwise disparate computational tools together into a single open source python package. Magnetic equilibrium, engineering CAD, finite volume solvers, scrape off layer plasma physics, visualization, high performace computing, and more, are connected in a single web-based user interface. Linux users may use HEAT without any software prerequisites via an appImage. This manuscript introduces HEAT, discusses the software architecture, presents first HEAT results, and outlines physics modules in development.
Type:
Dataset
Issue Date:
March 2021

92. Hybrid simulations of sub-cyclotron compressional and global Alfven Eigenmode stability in spherical tokamaks

Author(s):
Lestz, J.B.; Belova, E.V.; Gorelenkov, N.N
Abstract:
A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron frequency compressional (CAE) and global (GAE) Alfven Eigenmodes in low aspect ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam power in NSTX experiments, prompting further examination of the conditions for their excitation. Linear simulations are performed with the hybrid MHD-kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Three distinct types of modes are found in simulations -- co-CAEs, cntr-GAEs, and co-GAEs -- with differing spectral and stability properties. The simulations reveal that unstable GAEs are more ubiquitous than unstable CAEs, consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from simulations and estimated analytically for relevant sources not present in the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.
Type:
Dataset
Issue Date:
March 2021

93. Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

Author(s):
Mollen Albert; Adams Mark F.; Knepley Matthew G.; Hager Robert; Chang C. S.
Abstract:
The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a 5-dimensional (5D) continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular 2D velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.
Type:
Dataset
Issue Date:
March 2021

94. The updated ITPA global H-mode confinement database: description and analysis

Author(s):
Verdoolaege, G.; Kaye, S.M.; Angioni, C.; Kardaunn, O.W.J.F.; Maslov, M.; Romanelli, M.; Ryter, F.; Thomsen, K.
Abstract:
The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.
Type:
Dataset
Issue Date:
March 2021

95. Wave-kinetic approach to zonal-flow dynamics: recent advances

Author(s):
Zhu, Hongxuan; Dodin, I. Y.
Abstract:
Basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of wave-kinetic theory. Recently, this framework has been re-examined from first principles, which has led to more accurate yet still tractable "improved" wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow "curvature" (the second radial derivative of the flow velocity) on dynamics and stability of drift waves and zonal flows. We overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.
Type:
Dataset
Issue Date:
March 2021

96. Data for: Three-Dimensional Morphometry of Ooids in Oolites: a new tool for more accurate and precise paleoenvironmental interpretation

Author(s):
Howes, Bolton; Mehra, Akshay; Maloof, Adam
Abstract:
The prevalence of ooids in the stratigraphic record, and their association with shallow-water carbonate environments, make ooids an important paleoenvironmental indicator. Recent advances in the theoretical understanding of ooid morphology, along with empirical studies from Turks and Caicos, Great Salt Lake, and The Bahamas, have demonstrated that the morphology of ooids is indicative of depositional environment and hydraulic conditions. To apply this knowledge from modern environments to the stratigraphic record of Earth history, researchers measure the size and shape of lithified ooids on two-dimensional surfaces (i.e., thin sections or polished slabs), often assuming that random 2D slices intersect the nuclei and that the orientation of the ooids is known. Here we demonstrate that these assumptions rarely are true, resulting in errors of up to 35% on metrics like major axis length. We present a method for making 3D reconstructions by serial grinding and imaging, which enables accurate measurement of the morphology of individual ooids within an oolite, as well as the sorting and porosity of a sample. We also provide three case studies that use the morphology of ooids in oolites to extract environmental information. Each case study demonstrates that 2D measurements can be useful if the environmental signal is large relative to the error from 2D measurements. However, 3D measurements substantially improve the accuracy and precision of environmental interpretations. This study focuses on oolites, but errors from 2D measurements are not unique to oolites; this method can be used to extract accurate grain and porosity measurements from any lithified granular sample.
Type:
Dataset
Issue Date:
22 February 2021

98. Data for: 'How is sea level change encoded in carbonate stratigraphy?'

Author(s):
Geyman, Emily C.; Maloof, Adam C.; Dyer, Blake
Abstract:
The history of organismal evolution, seawater chemistry, and paleoclimate is recorded in layers of carbonate sedimentary rock. Meter-scale cyclic stacking patterns in these carbonates often are interpreted as representing sea level change. A reliable sedimentary proxy for eustasy would be profoundly useful for reconstructing paleoclimate, since sea level responds to changes in temperature and ice volume. However, the translation from water depth to carbonate layering has proven difficult, with recent surveys of modern shallow water platforms revealing little correlation between carbonate facies (i.e., grain size, sedimentary bed forms, ecology) and water depth. We train a convolutional neural network with satellite imagery and new field observations from a 3,000 km2 region northwest of Andros Island (Bahamas) to generate a facies map with 5 m resolution. Leveraging a newly-published bathymetry for the same region, we test the hypothesis that one can extract a signal of water depth change, not simply from individual facies, but from sequences of facies transitions analogous to vertically stacked carbonate strata. Our Hidden Markov Model (HMM) can distinguish relative sea level fall from random variability with ∼90% accuracy. Finally, since shallowing-upward patterns can result from local (autogenic) processes in addition to forced mechanisms such as eustasy, we search for statistical tools to diagnose the presence or absence of external forcings on relative sea level. With a new data-driven forward model that simulates how modern facies mosaics evolve to stack strata, we show how different sea level forcings generate characteristic patterns of cycle thicknesses in shallow carbonates, providing a new tool for quantitative reconstruction of ancient sea level conditions from the geologic record.
Type:
Dataset
Issue Date:
1 February 2021

99. Inversion technique to obtain local ion temperature profiles for an axisymmetric plasma with toroidal and radial velocities

Author(s):
Bell, Ronald E.
Abstract:
A matrix inversion technique is derived to calculate local ion temperature from line-integrated measurements of an extended emission source in an axisymmetric plasma which exactly corrects for both toroidal velocity and radial velocity components. Local emissivity and toroidal velocity can be directly recovered from line-integrated spectroscopic measurements, but an independent measurement of the radial velocity is necessary to complete the temperature inversion. The extension of this technique to handle the radial velocity is relevant for magnetic reconnection and merging compression devices where temperature inversion from spectroscopic measurements is desired. A simulation demonstrates the effects of radial velocity on the determination of ion temperature.
Type:
Dataset
Issue Date:
February 2021

100. Microtearing Instabilities and Electron Thermal Transport in Low and High Collisionality NSTX Discharges

Author(s):
Rafiq T; Kaye S; Guttenfelder W; Weiland J; Schuster E; Anderson J; Luo L;
Abstract:
Microtearing mode (MTM) real frequency, growth rate, magnetic fluctuation amplitude and resulting electron thermal transport are studied in systematic NSTX scans of relevant plasma parameters. The dependency of the MTM real frequency and growth rate on plasma parameters, suitable for low and high collision NSTX discharges, is obtained by using the reduced MTM transport model [T. Rafiq, et al., Phys. Plasmas 23, 062507 (2016)]. The plasma parameter dependencies are compared and found to be consistent with the results obtained from MTM using the Gyrokinetic GYRO code. The scaling trend of collision frequency and plasma beta is found to be consistent with the global energy confinement trend observed in the NSTX experiment. The strength of the magnetic fluctuation is found to be consistent with the gyrokinetic estimate.In earlier studies, it was found that the version of the Multi-Mode (MM) anomalous transport model, which did not contain the effect of MTMs, provided an appropriate description of the electron temperature profiles in standard tokamak discharges and not in spherical tokamaks. When the MM model, which involves transport associated with MTMs, is incorporated in the TRANSP code and is used in the study of electron thermal transport in NSTX discharges, it is observed that the agreement with the experimental electron temperature profile is substantially improved.
Type:
Dataset
Issue Date:
February 2021