Search NASA⌕ Search

Engineering topics

Hughes, J.

Publications and source records attributed to Hughes, J..

At least 19 records

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Divertor heat flux challenge and mitigation in SPARC

Owing to its high magnetic field, high power, and compact size, the SPARC experiment will operate with divertor conditions at or above those expected in reactor-class tokamaks. Power exhaust at this scale remains one of the key challenges for practical fusion energy. Based on empirical scalings, the peak unmitigated divertor parallel heat flux is projected to be greater than 10 GW m -2 . This is nearly an order of magnitude higher than has been demonstrated to date. Furthermore, the divertor parallel Edge-Localized Mode (ELM) energy fluence projections (~11–34 MJ m -2 ) are comparable with those for ITER. However, the relatively short pulse length (~25 s pulse, with a ~10 s flat top) provides the opportunity to consider mitigation schemes unsuited to long-pulse devices including ITER and reactors. The baseline scenario for SPARC employs a ~1 Hz strike point sweep to spread the heat flux over a large divertor target surface area to keep tile surface temperatures within tolerable levels without the use of active divertor cooling systems. In addition, SPARC operation presents a unique opportunity to study divertor heat exhaust mitigation at reactor-level plasma densities and power fluxes. Not only will SPARC test the limits of current experimental scalings and serve for benchmarking theoretical models in reactor regimes, it is also being designed to enable the assessment of long-legged and X-point target advanced divertor magnetic configurations. Experimental results from SPARC will be crucial to reducing risk for a fusion pilot plant divertor design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FY19 FES Theory Performance Target (Final Report)

Understanding the relevant turbulent transport mechanisms at the edge of a high-performance tokamak is essential for predicting and optimizing the H-mode pedestal structure in future burning plasma devices. Global electromagnetic gyrokinetic simulations will be performed based on representative experimental pedestal scenarios in order to clarify which instabilities are most important for each of the particle and heat transport channels. Edge transport modeling will be performed in order to estimate and bound the particle and heat sources—e.g., the ionization density source and the atomic energy loss channels due to ionization, charge exchange, and radiation. Comparisons will be made with data from the DIII-D, JET, C-Mod and NSTX or MAST experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surface heat flux feedback controlled impurity seeding experiments with Alcator C-Mod’s high-Z vertical target plate divertor: performance, limitations and implications for fusion power reactors

The Alcator C-Mod team has recently developed a feedback system to measure and control surface heat flux in real-time. The system uses real-time measurements of surface heat flux from surface thermocouples and a pulse-width modulated piezo valve to inject low-Z impurities (typically N2) into the private flux region. It has been used in C-Mod to mitigate peak surface heat fluxes >40 MW/m2 down to <10 MW/m2 while maintaining excellent core confinement, H98>1. While the system works quite well under relatively steady conditions, use of it during transients has revealed important limitations on feedback control of impurity seeding in conventional vertical target plate divertors. In some cases, the system is unable to avoid plasma reattachment to the divertor plate or the formation of a confinement-damaging x-point MARFE. This is due to the small operational window for mitigated heat flux in the parameters of incident plasma heat flux, plasma density, and impurity density as well as the relatively slow response of the impurity gas injection system compared to plasma transients. Given the severe consequences for failure of such a system to operate reliably in a reactor, there is substantial risk that the conventional vertical target plate divertor will not provide an adequately controllable system in reactor-class devices. These considerations motivate the need to develop passively stable, highly compliant divertor configurations and experimental facilities that can test such possible solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1

Here, a fast edge turbulence suppression event has been simulated in the electrostatic version of the gyrokinetic particle-in-cell code XGC1 in a realistic diverted tokamak edge geometry under neutral particle recycling. The results show that the sequence of turbulent Reynolds stress followed by neoclassical ion orbit-loss driven together conspire to form the sustaining radial electric field shear and to quench turbulent transport just inside the last closed magnetic flux surface. As a result, the main suppression action is located in a thin radial layer around ψ N ≃0.96–0.98, where ψ N is the normalized poloidal flux, with the time scale ~0.1 ms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surface heat flux feedback controlled impurity seeding experiments with Alcator C-Mod's high-Z vertical target plate divertor: performance, limitations and implications for fusion power reactors

The Alcator C-Mod team has recently developed a feedback system to measure and control surface heat flux in real-time. The system uses real-time measurements of surface heat flux from surface thermocouples and a pulse-width modulated piezo valve to inject low-Z impurities (typically N 2 ) into the private flux region. It has been used in C-Mod to mitigate peak surface heat fluxes >40 MW m -2 down to <10 MW m -2 while maintaining excellent core confinement, H 98 > 1. While the system works quite well under relatively steady conditions, use of it during transients has revealed important limitations on feedback control of impurity seeding in conventional vertical target plate divertors. In some cases, the system is unable to avoid plasma reattachment to the divertor plate or the formation of a confinement-damaging x-point MARFE. This is due to the small operational window for mitigated heat flux in the parameters of incident plasma heat flux, plasma density, and impurity density as well as the relatively slow response of the impurity gas injection system compared to plasma transients. Given the severe consequences for failure of such a system to operate reliably in a reactor, there is substantial risk that the conventional vertical target plate divertor will not provide an adequately controllable system in reactor-class devices. Finally, these considerations motivate the need to develop passively stable, highly compliant divertor configurations and experimental facilities that can test such possible solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effects of main-ion dilution on turbulence in low q 95 C-Mod ohmic plasmas, and comparisons with nonlinear GYRO

Recent experiments on C-mod seeding nitrogen into ohmic plasmas with q 95 = 3.4 found that the seeding greatly reduced long-wavelength (ITG-scale) turbulence. The long-wavelength turbulence that was reduced by the nitrogen seeding was localized to the region of r/a ≈ 0.85, where the turbulence is well above marginal stability (as evidenced by Q i /Q GB $\gg$ 1). The nonlinear gyrokinetic code GYRO was used to simulate the expected turbulence in these plasmas, and the simulated turbulent density fluctuations and turbulent energy fluxes quantitatively agreed with the experimental measurements both before and after the nitrogen seeding. Unexpectedly, the intrinsic rotation of the plasma was also found to be affected by the nitrogen seeding, in a manner apparently unrelated to a change in the electron-ion collisionality that was proposed by other experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and modeling uncertainties in the validation of lower hybrid current drive

Our work discusses sources of uncertainty in the validation of lower hybrid wave current drive simulations against experiments, by evolving self-consistently the magnetic equilibrium and the heating and current drive profiles, calculated with a combined toroidal ray tracing code and 3D Fokker–Planck solver. The simulations indicate a complex interplay of elements, where uncertainties in the input plasma parameters, in the models and in the transport solver combine and compensate each other, at times. It is concluded that ray-tracing calculations should include a realistic representation of the density and temperature in the region between the confined plasma and the wall, which is especially important in regimes where the LH waves are weakly damped and undergo multiple reflections from the plasma boundary. Uncertainties introduced in the processing of diagnostic data as well as uncertainties introduced by model approximations are assessed. We show that, by comparing the evolution of the plasma parameters in self-consistent simulations with available data, inconsistencies can be identified and limitations in the models or in the experimental data assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effects of dilution on turbulence and transport in C-Mod ohmic plasmas and comparisons with gyrokinetic simulations

Main ion dilution has been predicted by gyrokinetic simulations to have a significant effect on ion thermal transport in C-Mod ohmic plasmas. This effect was verified experimentally with a specific set of experiments on C-Mod in which ohmic deuterium plasmas across the linear ohmic confinement (LOC) through the saturated ohmic confinement (SOC) regimes were diluted by seeding with nitrogen gas (Z=7) injection. The seeding was observed to increase the normalized ion temperature gradients (ITGs) by up to 30% without a corresponding increase in the gyrobohm normalized ion energy flux, indicating a change in either the stiffness or the critical ion temperature gradient associated with ITG turbulence. The seeding also reversed the direction of the intrinsic toroidal rotation in plasmas slightly above the normal intrinsic rotation reversal critical density. GYRO simulations of the seeded and unseeded plasmas show that the seeding affected both the critical gradient and the stiffness. For plasmas in the LOC regime, the dilution primarily increased the critical gradient, while for plasmas in the SOC regime the dilution primarily decreased the stiffness. At r/a = 0.8, where the experimental fluxes were above marginal stability, local GYRO predicted and experimental energy fluxes agreed, except for Q i in the SOC regime where GYRO under-predicted the experimental energy flux. At r/a = 0.6, where the experimental fluxes were close to marginally stable, local GYRO predicted ITG modes to be strongly unstable and are responsible for both Q i and Q e (with Q i > Q e ), as opposed to the experiment where Q i < Q e . In contrast, global GYRO in this region predicted the ITG modes to be closer to marginal stability, and accurately predict the experimental Q i when the T i profile is modified within experimental uncertainties. Further, the fact that Q e is always less than Q i in the r/a = 0.6 simulations with k θ ρ s ≤1 indicates that high-k electron temperature gradient driven (ETG) modes must be included in future simulations and may be responsible for the electron energy transport in this case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-resonant destabilization of (1/1) internal kink mode by suprathermal electron pressure

New experimental observations are reported on the structure and dynamics of short-lived periodic (1, 1) “fishbone”-like oscillations that appear during radio frequency heating and current-drive experiments in tokamak plasmas. For the first time, measurements can directly relate changes in the high energy electrons to the mode onset, saturation, and damping. In the relatively high collisionality of Alcator C-Mod with lower hybrid current drive, the instability appears to be destabilized by the non-resonant suprathermal electron pressure—rather than by wave-particle resonance, rotates toroidally with the plasma and grows independently of the (1, 1) sawtooth crash driven by the thermal plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Correlation ECE diagnostic in Alcator C-Mod

Correlation ECE (CECE) is a diagnostic technique that allows measurement of small amplitude electron temperature, Te, fluctuations through standard cross-correlation analysis methods. In Alcator C-Mod, a new CECE diagnostic has been installed[Sung RSI 2012], and interesting phenomena have been observed in various plasma conditions. We find that local Te fluctuations near the edge (ρ ~ 0:8) decrease across the linearto- saturated ohmic confinement transition, with fluctuations decreasing with increasing plasma density[Sung NF 2013], which occurs simultaneously with rotation reversals[Rice NF 2011]. Te fluctuations are also reduced across core rotation reversals with an increase of plasma density in RF heated L-mode plasmas, which implies that the same physics related to the reduction of Te fluctuations may be applied to both ohmic and RF heated L-mode plasmas. In I-mode plasmas, we observe the reduction of core Te fluctuations, which indicates changes of turbulence occur not only in the pedestal region but also in the core across the L/I transition[White NF 2014]. The present CECE diagnostic system in C-Mod and these experimental results are described in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Planetary Data System - A Case Study in the Development and Management of Meta-Data for a Scientific Digital Library

The Planetary Data System (PDS) is an active science data archive managed by scientists for NASA's planetary science community. With the advent of the World Wide Web the majority of the archive has been placed on-line as a science digital libraty for access by scientists, the educational community, and the general public.

Data System Meta-data scientific digital library↗

Analysis of LAC Observations of Clusters of Galaxies and Supernova Remnants

The following publications are included and serve as the final report: The X-ray Spectrum of Abell 665; Clusters of Galaxies; Ginga Observation of an Oxygen-rich Supernova Remnant; Ginga Observations of the Coma Cluster and Studies of the Spatial Distribution of Iron; A Measurement of the Hubble Constant from the X-ray Properties and the Sunyaev-Zel'dovich Effect of Abell 2218; Non-polytropic Model for the Coma Cluster; and Abundance Gradients in Cooling Flow Clusters: Ginga LAC (Large Area Counter) and Einstein SSS (Solid State Spectrometer) Spectra of A496, A1795, A2142, and A2199.

Hughes, J.↗

Pinatubo Global- to Micro-Scale Evolution: A Unified Picture from Space, Air, and Ground Measurements

We combine space, air, and ground measurements to develop a composite picture of the post-Pinatubo aerosol, and assess the consistency and uncertainties of various measurement and retrieval techniques. impactor and optical counter measurements, as well as retrievals from optical depth spectra, paint a generally consistent picture of the evolution of particle effective radii, R(sub eff). In the first month after the eruption, although particle numbers increased by orders of magnitude, R(sub eff) was similar to the preeruption value of 4.2 micrometers, because both small (r less than 0.25 micrometers) and large (r greater than 0.6 micrometers) particles increased in number, Over the next 3-6 months, R(sub eff) increased rapidly to about 0.5 micrometers. In general, R(sub eff) continued to increase for about a year after the eruption. The peak wavelength of optical depth spectra increased from initial values of less than 0.42 micrometers to values between 0.78 and 1 micrometer. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between the global maps of 0.5 and 1.0-micrometer optical depth derived from the AVHRR and SAGE satellite measurements. It also sets a context for evaluating remaining uncertainties in each of these satellite data products. We also make consensus recommendations for particle composition, shape, and temperature- and wavelength-dependent refractive index, and show how the latter effect on backscatter spectra can influence particle sizes retrieved from multiwavelength lidar measurements.

Russell, Philip B.↗

Pinatubo Aerosol Global-to-Micro-Scale Evolution: A Unified Picture From Space, Air, and Ground Measurements

We combine a variety of measurements to develop a composite picture of the post-Pinatubo aerosol and assess the consistency and uncertainties of the measurement and retrieval techniques Satellite infrared spectroscopy, particle morphology, and evaporation temperature measurements are in accord with theoretical calculations In showing a dominant particle composition of H2SO4-H2O mature with H2SO4 weight traction of 65-80% for most stratospheric temperatures and humidities. Important exceptions are: (1) the presence of volcanic ash at all attitudes initially and in a layer just above the tropopause until at least March 1992, and (2) much smaller H2SO4 weight fractions at the low temperatures attained In high latitude winters and at the tropical tropopause. Laboratory spectroscopy and theoretical calculations yield wavelength- and temperature-dependent refractive indices for the dominant H2SO4-H2O droplets. These in turn permit derivation of particle size spectra from measured optical depth spectra for comparison to direct measurements by impactors and optical counters, All three techniques paint a generally consistent picture of the evolution of R(sub eff) the effective, or area-weighted, particle radius. In the first month after the eruption, although particle numbers increased by orders of magnitude, R(sub eff) was similar to the pre-eruption value of 0.1 to 0.2 microns because both small (r less than 0.2 microns) and large (r greater than 0.6 microns) particles increased in number. Over the next 3-6 months, R(sub eff) increased to about 0.5 microns, reflecting particle growth through condensation and coagulation. In general, R(sub eff) continued to increase for about a year after the eruption. Extinction spectra computed from in situ size distribution measurements are consistent with optical depth measurements, which show spectra with maxima initially at wavelengths of 0.42 microns or less, and thereafter progressively increasing to between 0.78 and 1 micron. Not until 1993 does optical depth spectra begin to show a clear return to the preeruption signature of maximizing at the shortest visible wavelengths or in the near UV. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between the global maps of 0.5- and 1.0-kilometer optical depth derived from the AVHRR and SAGE satellite measurements. It also sets a context for evaluating remaining uncertainties in each of these satellite data products. We also show how the effects of wavelength-dependent refractive index on backscatter spectra can influence particle sizes retrieved from multiwavelength lidar measurements.

Russell, Philip B.↗

A Unified Picture of Pinatubo Aerosol Global-to Micro-Scale Evolution, From Space, Air, and Ground Measurements

We combine a variety of measurements to develop a composite picture of the post-Pinatubo aerosol and assess the consistency and uncertainties of the measurement and retrieval techniques. Satellite infrared spectroscopy, particle morphology, and evaporation temperature measurements are in accord with theoretical calculations in showing a dominant particle composition of H2SO4-H2O mixture, with H2SO4 weight fraction of 65-80% for most stratospheric temperatures and humidities. Important exceptions are: (1) the presence of volcanic ash at all altitudes initially and in a layer just above the tropopause until at least March 1992, and (2) much smaller H2SO4 weight fractions at the low temperatures attained in high latitude winters and at the tropical tropopause, Laboratory spectroscopy and theoretical calculations yield wavelength- and temperature-dependent refractive indices for the dominant H2SO4-H2O droplets. These in turn permit derivation of particle size spectra from measured optical depth spectra, for comparison to direct measurements by impactors and optical counters. All three techniques paint a generally consistent picture of the evolution of R(sub eff), the effective, or area-weighted, particle radius. In the first month after the eruption, although particle numbers increased by orders of magnitude, R(sub eff) was similar to the preemption value of 0.1 to 0.2 microns, because both small (r less than 0.2 microns) and large (r greater than 0.6 micron particles increased in number. Over the next 3-6 months, R(sub eff) increased to about 0.5 microns reflecting particle growth through condensation and coagulation. In general, R(sub eff) continued to increase for about a year after the eruption. Extinction spectra computed from in situ size distribution measurements are consistent with optical depth measurements, which show spectra with maxima initially at wavelengths of 0.42 microns or less, and thereafter progressively increasing to between 0.78 and 1 micron. Not until 1993 do optical depth spectra begin to show a clear return to the preemption signature of maximizing at the shortest visible wavelengths or in the near UV. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between the global maps of 0.5- 1.0- micron optical depth derived from the AVHRR and SAGE satellite measurements.

Russell, Philip B.↗

HRMA calibration handbook: EKC gravity compensated XRCF models

This document, consisting of hardcopy printout of explanatory text, figures, and tables, represents one incarnation of the AXAF high resolution mirror assembly (HRMA) Calibration Handbook. However, as we have envisioned it, the handbook also consists of electronic versions of this hardcopy printout (in the form of postscript files), the individual scripts which produced the various figures and the associated input data, the model raytrace files, and all scripts, parameter files, and input data necessary to generate the raytraces. These data are all available electronically as either ASCII or FITS files. The handbook is intended to be a living document and will be updated as new information and/or fabrication data on the HRMA are obtained, or when the need for additional results are indicated. The SAO Mission Support Team (MST) is developing a high fidelity HRMA model, consisting of analytical and numerical calculations, computer software, and databases of fundamental physical constants, laboratory measurements, configuration data, finite element models, AXAF assembly data, and so on. This model serves as the basis for the simulations presented in the handbook. The 'core' of the model is the raytrace package OSAC, which we have substantially modified and now refer to as SAOsac. One major structural modification to the software has been to utilize the UNIX binary pipe data transport mechanism for passing rays between program modules. This change has made it possible to simulate rays which are distributed randomly over the entrance aperture of the telescope. It has also resulted in a highly efficient system for tracing large numbers of rays. In one application to date (the analysis of VETA-I ring focus data) we have employed 2 x 10(exp 7) rays, a substantial improvement over the limit of 1 x 10(exp 4) rays in the original OSAC module. A second major modification is the manner in which SAOsac incorporates low spatial frequency surface errors into the geometric raytrace. The original OSAC included the ability to use Legendre-Fourier polynomials to describe deviations from the basic optical prescription. To this we have added bicubic splines to address a deficiency in the handling of the sharper deformations in the areas of mirror support pads. SAO has developed software (TRANS-FIT) to translate the most common finite element analysis models into these forms for incorporation into the raytrace program.

Tananbaum, H. D.↗