Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA RADIAL PROFILES”

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 397 records · Page 22

Simulation of triggering and evolution of ELM by pellet injection in EAST under BOUT++ framework

A BOUT ++ three-field magnetohydrodynamic model is employed to study the triggering and evolution of edge localized mode (ELM) by Li pellets injected along the outer mid-plane in the EAST configuration. The linear simulation shows that compared with a large deposition on the pedestal top (scenario I), a smaller deposition within the steep-gradient pedestal region (scenario II) can stimulate much larger linear growth rates of all-n peeling-ballooning modes (PBMs). The nonlinear simulation shows that there exists a pellet size threshold for ELM triggering for two deposition locations; the threshold for scenario I predicted in the present study matches the EAST observation well. Comparison of the two scenarios reveals that a smaller deposition is sufficient to trigger an ELM in a much shorter time in scenario II, whose ELM size is comparable to that in scenario I. This conclusion confirms previous DIII-D and ASDEX-Upgrade observations, suggesting that the steep-gradient pedestal region is a favorable deposition location for ELM triggering with minimum pellet size. Simulation analyses also find that the positive radial gradient of the hump-like pressure profile in the outer mid-plane induced by the pellet deposition plays a different role in the two scenarios. In scenario I, the force resulting from the gradient hinders the outflow of core plasmas and in return, the perturbation is suppressed from spreading inwards after ELM crashes. In scenario II, with a sizable deposition, the gradient results in another competitive perturbation growth region during the linear phase, thus dispersing the free energy and reducing the efficiency of destabilizing PBMs by pellet injection. The suppressing effect of saturated zonal flow on other modes, the short ELM fast crash phase, and the restricting transport effect of the positive radial pressure gradient work together to constrain the pedestal energy loss, especially when the pellet deposition amount is high.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluid and gyrokinetic turbulence in open field-line, helical plasmas

Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently, it has become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas [T. N. Bernard et al., Phys. Plasmas 26, 042301 (2019)]. Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs that can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, and improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Langmuir and Spatially Resolved OES Diagnostics for Arc Jet Shock and Boundary Layer Measurements Relevant to Atmospheric Re-Entry Environments

The second-generation miniature Arc Jet Research Chamber (mARC) ground testing facility at NASA Ames requires physical characterization of its supersonic plasma flow to verify numerical models and material sample testing. Two types of plasma diagnostics will be utilized: Langmuir probes and optical emission spectroscopy (OES). Several single and triple Langmuir probes have been constructed and will be inserted parallel to the flow via a sweep arm to measure post-shock electron temperature and number density along radial profiles at different axial distances. Due to the high heat flux (200-3000 W/cm2) sweep durations are necessarily short (< 0.1 s). Measuring these properties near the probe surface provide characterization of the post-shock environment. This is an important metric for validating coupled CFD-material codes, and will provide insight into longstanding thermocouple anomalies attributed to charging of heatshield surfaces. An OES system has also been assembled to obtain spatially resolved density and temperature for various flow species. The system consists of a convex lens that forms an image of the plasma onto a 16-channel linear fiber array, which will simultaneously capture the spectra from sixteen different locations. The light emission from the plasma is then focused by a spherical mirror into a spectrometer and recorded with a camera. These measurements non-intrusively provide a new spatially resolved flow characterization capability. This will initially be used to resolve issues with shock radiation interference in pyrometer measurements of material surface temperature. The design of each diagnostic and any initial results will be discussed.

Ethan Main-Loam Leong↗

Observation of weakly damped modes using high resolution measurement of turbulence in a dipole confined plasma

When plasma is confined by a levitated superconducting current ring, the turbulent radial particle flux changes the direction depending upon the gradient of the particle number within tubes of equal magnetic flux [Garnier et al., Phys. Plasmas 24, 012506 (2017)]. When the density profile is highly peaked, the particle flux is outward. When the density profile is less peaked, turbulence drives an inward “particle pinch.” This change in the direction of particle flux coincides with a reversal of the mean toroidal propagation direction of turbulent fluctuations. Previously, only the mean wavenumber of the fluctuations was reported. The mean wavenumber was computed from the ensemble correlation between two probes. Here, the full frequency-wavenumber power spectrum of the saturated turbulence is calculated by applying Capon’s “maximum likelihood method” to floating potential measurements from an array of probes. The power spectrum shows that the highest intensity modes are due to unstable entropy modes, and the less intense modes coincide with weakly damped convective cells, which rotate toroidally in the direction opposite to the unstable modes. Because the weakly damped modes make up a significant fraction of the fluctuating power, these measurements imply that saturated turbulence in a dipole-confined plasma involves the excitation of weakly damped rotating convective cells.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculation of the radial electric field from a modified Ohm's law

A modified Ohm's Law, derived from the conservation of deuterium and carbon ions and electron momentum and the requirement for charge neutrality, yields an expression for the radial electric field, E r , in the edge pedestal region in terms of the motional electric field due to the carbon and deuterium ion rotation velocities as well as pressure gradients and the radial plasma current. This analytical Ohm's Law model for E r is first shown to be consistent with the conventional "experimental" electric field calculated from the carbon radial momentum balance using experimental carbon rotation and pressure gradient measurements when experimental profiles are used to evaluate the Ohm's Law in three DIII-D [Luxon, Nucl. Fusion 42, 614 (2002)] representative discharges (for L-mode, H-mode, and Resonant Magnetic Perturbation operating regimes). In order to test the practical predictive ability of the modified Ohm's Law, the calculations were repeated using rotation velocities calculated with neoclassical rotation models instead of measured rotation velocities. The Ohm's Law predicted E r using theoretical rotation velocities did not agree with the "experimental" E r as well as the Ohm's Law prediction using experimental rotation velocities, indicating that more accurate models for predicting edge rotation velocity are thus needed in order to have a validated predictive model of E r in the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of gas-puff Z-pinches on COBRA

Gas-puff Z-pinch experiments were conducted on the 1 MA, 200 ns pulse duration Cornell Beam Research Accelerator (COBRA) pulsed power generator in order to achieve an understanding of the dynamics and instability development in the imploding and stagnating plasma. The triple-nozzle gas-puff valve, pre-ionizer, and load hardware are described. Specific diagnostics for the gas-puff experiments, including a Planar Laser Induced Fluorescence system for measuring the radial neutral density profiles along with a Laser Shearing Interferometer and Laser Wavefront Analyzer for electron density measurements, are also described. The results of a series of experiments using two annular argon (Ar) and/or neon (Ne) gas shells (puff-on-puff) with or without an on- (or near-) axis wire are presented. For all of these experiments, plenum pressures were adjusted to hold the radial mass density profile as similar as possible. Initial implosion stability studies were performed using various combinations of the heavier (Ar) and lighter (Ne) gasses. Implosions with Ne in the outer shell and Ar in the inner were more stable than the opposite arrangement. Current waveforms can be adjusted on COBRA and it was found that the particular shape of the 200 ns current pulse affected on the duration and diameter of the stagnated pinched column and the x-ray yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ENERGY CONTENT AND PROPAGATION IN TRANSVERSE SOLAR ATMOSPHERIC WAVES

Recently, a significant amount of transverse wave energy has been estimated propagating along solar atmospheric magnetic fields. However, these estimates have been made with the classic bulk Alfven wave model which assumes a homogeneous plasma. In this paper, the kinetic, magnetic, and total energy densities and the flux of energy are computed for transverse MHD waves in one-dimensional cylindrical flux tube models with a piecewise constant or continuous radial density profile. There are fundamental deviations from the properties for classic bulk Alfven waves. (1) There is no local equipartition between kinetic and magnetic energy. (2) The flux of energy and the velocity of energy transfer have, in addition to a component parallel to the magnetic field, components in the planes normal to the magnetic field. (3) The energy densities and the flux of energy vary spatially, contrary to the case of classic bulk Alfven waves. This last property has the important consequence that the energy flux computed with the well known expression for bulk Alfven waves could overestimate the real flux by a factor in the range 10-50, depending on the flux tube equilibrium properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Rocket spectroscopy of zeta Orionis.

Analysis of a spectrum of zeta Ori extending from 922 to 1453 A with approximately 0.8 A resolution recorded at rocket altitudes. All lines used in existing models of stellar atmospheres appear in the recorded spectrum with the exception of those masked by telluric N2 or strong P Cygni-type profiles and by an O V line at 1371.29 A. Fifteen multiplets of subordinate lines have been reliably identified, indicating an approximate range of excitation from 0 to 50 eV. Transitions in C III (1176 A), N III (991 A), N V (1239, 1243 A), O VI (1032, 1038 A), Si IV (1394, 1403 A), S IV (1063, 1074 A), and S VI (933, 944 A) have been observed as P Cygni-type profiles presumably arising in a circumstellar envelope. The degree of ionization, transitions present, and mean radial velocities are all consistent with viewing the envelope as a hot (about 100,000 K), rarefied plasma in which collisional ionization is important. Interstellar lines in C I (1277, 1280 A), C II (1036, 1334 A), N I (1134-1135 A), N I (1200-1201 A), N II (1084-1086 A), O I (1302, 1305 A), Si II (1190, 1193 A), Si II (1260 A), and Si II (1304 A) have been definitely identified. Other transitions in Ar II, S I, C I, and Fe II are tentatively identified. The equivalent width of the L alpha line is found to be 10.4 plus or minus 1.6 A, corresponding to a columnar density of 2.0 plus or minus 0.7 x 10 to the 20th per cu cm.

Smith, A. M.↗

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport dynamics of self-consistent, near-marginal drift-wave turbulence. I. Investigation of the ability of external flows to tune the non-diffusive dynamics

The reduction of turbulent transport across sheared flow regions has been known for a long time in magnetically confined toroidal plasmas. However, details of the dynamics are still unclear, in particular, in what refers to the changes caused by the flow on the nature of radial transport itself. In Paper II, we have shown in a simplified model of drift wave turbulence that, when the background profile is allowed to evolve self-consistently with fluctuations, a variety of transport regimes ranging from superdiffusive to subdiffusive open up depending on the properties of the underlying turbulence. In this study, we show that externally applied sheared flows can, under the proper conditions, cause the transport dynamics to be diffusive or subdiffusive

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory of a steady-state nonisothermal positive column in a magnetic field.

The steady state of a nonisothermal cylindrical discharge in a longitudinal magnetic field is described by three moment equations in the plasma approximation for electrons and ions. Volume generation of electron-ion pairs by single-stage ionization, the presence of axial current, and collisions with neutrals are considered. The theory yields radial profiles of the charged particle velocities, density, potential, electron and ion temperatures, and demonstrates similarity laws for the positive column. The results are compared with two-moment theories and with experimental data on helium, argon, and mercury vapor found in the literature for a wide range of pressures.

Ilic, D. B.↗

A summary of whistlers observed by Voyager 1 at Jupiter

We summarize the Voyager 1 observations of whistlers at Jupiter in order to provide a basis for further analyses of the density profile of the Io plasma torus as well as to support studies of atmospheric lightning at Jupiter. All the whistlers detected by Voyager 1 fell into three general regions in the torus at radial distances ranging beteen 5 and 6 R sub J. An analysis of the broadband wave amplitudes measured by the Voyayer 1 plasma wave instrument and estimates of the peak whistler amplitudes imply the grouping of whistlers was due to variations in the sensitivity of the receiver to whistlers and not to variations in the source or propagation paths of the whistlers. The whistler dispersions are presented in statistical form for each of the three groups of events and analyzed in view of the structure of the Io plasma torus as determined by plasma probe measurements. The results of these analyses give source locations for the whistlers at the foot of the magnetic field lines threading the torus in both hemispheres and over a range of longitudes.

Kurth, W. S.↗

Behavior of lithium ions in the turbulent near-wall tokamak plasma under heating of ions and electrons of the main plasma

Turbulent dynamics of the near-wall tokamak plasma is simulated by numerically solving the nonlinear reduced Braginskii magnetohydrodynamic equations with allowance for a lithium ion admixture. The effects of turbulence and radiation of the admixture are analyzed in the framework of a self-consistent approach. The radial distributions of the radiative loss power and the density of Li{sup 0} atoms and Li{sup +1} ions are obtained as functions of the electron and ion temperatures of the main plasma in the near-wall layer. The results of numerical simulations show that supply of lithium ions into the low-temperature near-wall plasma substantially depends on whether the additional power is deposited into the electron or ion component of the main plasma. If the electron temperature in the layer increases (ECR heating), then the ion density drops. At the same time, an increase in the temperature of the main ions (ICR heating) leads to an increase in the density of Li{sup +1} ions. The results of numerical simulations are explained by the different influence of the electron and ion temperatures on the atomic processes governing the accumulation and loss of particles in the balance equations for neutral Li{sup 0} atoms and Li{sup +1} ions in the admixture. The radial profile of the electron temperature and the corresponding distribution of the radiative loss power for different densities of neutral Li{sup 0} atoms on the wall are obtained. The calculations show that the presence of Li{sup +1} ions affects turbulent transport of the main ions. In this case, the electron heat flux increases by 20–30% with increasing Li{sup +1} density, whereas the flux of the main ions drops by nearly the same amount. The radial profile of the turbulent flux of lithium ions is obtained. It is demonstrated that the appearance of the pinch effect is related to the positive density gradient of lithium ions across the calculation layer. For the parameters of the T-10 tokamak, the effect of radiative cooling of the near-wall plasma layer becomes appreciable when the near-wall density of neutral lithium atoms exceeds 7 × 10{sup 11} cm{sup −3}. In this case, the density of radiative loss power in the center of the layer is estimated to be about 500–600 kW/m{sup 3}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-resolution tangential absolute extreme ultraviolet arrays for radiated power density measurements on NSTX-U

The radiated-power-density diagnostic on the equatorial midplane for the NSTX-U tokamak will be upgraded to measure the radial structure of the photon emissivity profile with an improved radial resolution. This diagnostic will enhance the characterization and studies of power balance, impurity transport, and MHD. The layout and response expected of the new system is shown for different plasma conditions and impurity concentrations. The effect of toroidal rotation driving poloidal asymmetries in the core radiation from high-Z impurities is also addressed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Verification and validation of linear gyrokinetic and kinetic-MHD simulations for internal kink instability in DIII-D tokamak

Here, verification and linear validation of the internal kink instability in tokamak have been performed for both gyrokinetic (GTC) and kinetic-MHD codes (GAM-solver, M3D-C1-K, NOVA, XTOR-K). Using realistic magnetic geometry and plasma profiles from the same equilibrium reconstruction of the DIII-D shot #141216, these codes exhibit excellent agreement for the growth rate and mode structure of the internal kink mode when all kinetic effects are suppressed. The simulated radial mode structures, obtained from linear simulations, are in reasonable agreement with the normalised electron cyclotron emission measurement after adjusting, within the experimental uncertainty, the safety factor q = 1 flux-surface location in the equilibrium reconstruction. Compressible magnetic perturbations strongly destabilize the kink, while poloidal variations of the equilibrium current density reduce the growth rate of the kink. Furthermore, kinetic effects of thermal ions are found to decrease the kink growth rate in kinetic-MHD simulations, but increase the kink growth rate in gyrokinetic simulations, due to the additional drive of the ion temperature gradient and parallel electric field. Kinetic thermal electrons are found to have negligible effects on the internal kink instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High energy density Z-pinch plasmas using flow stabilization

The ZaP Flow Z-Pinch research project[1] at the University of Washington investigates the effect of sheared flows on MHD instabilities. Axially flowing Z-pinch plasmas are produced that are 100 cm long with a 1 cm radius. The plasma remains quiescent for many radial Alfvén times and axial flow times. The quiescent periods are characterized by low magnetic mode activity measured at several locations along the plasma column and by stationary visible plasma emission. Plasma evolution is modeled with high-resolution simulation codes – Mach2, WARPX, NIMROD, and HiFi. Plasma flow profiles are experimentally measured with a multi-chord ion Doppler spectrometer. A sheared flow profile is observed to be coincident with the quiescent period, and is consistent with classical plasma viscosity. Equilibrium is determined by diagnostic measurements: interferometry for density; spectroscopy for ion temperature, plasma flow, and density[2]; Thomson scattering for electron temperature; Zeeman splitting for internal magnetic field measurements[3]; and fast framing photography for global structure. Wall stabilization has been investigated computationally and experimentally by removing 70% of the surrounding conducting wall to demonstrate no change in stability behavior.[4] Experimental evidence suggests that the plasma lifetime is only limited by plasma supply and current waveform. The flow Z-pinch concept provides an approach to achieve high energy density plasmas,[5] which are large, easy to diagnose, and persist for extended durations. A new experiment, ZaP-HD, has been built to investigate this approach by separating the flow Z-pinch formation from the radial compression using a triaxial-electrode configuration. This innovation allows more detailed investigations of the sheared flow stabilizing effect, and it allows compression to much higher densities than previously achieved on ZaP by reducing the linear density and increasing the pinch current. Experimental results and scaling analyses will be presented. In addition to studying fundamental plasma science and high energy density physics, the ZaP and ZaP-HD experiments can be applied to laboratory astrophysics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport-driven toroidal rotation with general viscosity profile

Abstract Using the assumption of a weak normalized turbulent viscosity, usually valid in practice, the modulated-transport model (Stoltzfus-Dueck 2012 Phys. Plasmas 19 055908) is generalized to allow the turbulent transport coefficient to vary in an arbitrary way on radial and poloidal position. The new approach clarifies the physical interpretation of the earlier results and significantly simplifies the calculation, via a boundary-layer asymptotic method. Rigorous detailed appendices verify the result of the simple boundary-layer calculation, also demonstrating that it achieves the claimed order of accuracy and providing a concrete prediction for the strong plasma flows in the immediate vicinity of the last closed flux surface. The new formulas are used to predict plasma rotation at the core-edge boundary, in cases with and without externally applied torque. Dimensional formulas and extensive discussion are provided, to support experimental application of the new model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗