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 379 records · Page 21

Turbulence spreading effects on the ELM size and SOL width

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width (λ q ) broadening in small ELM regimes. Here, this study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (E r ) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal E r well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width (λ q ) is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width (λ q ) broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux Γ ε at the separatrix can be enhanced by increasing either pedestal E r flow shear or local SOL pressure gradient. The λ q is broadened as the fluctuation energy intensity flux Γ ε at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detection of an electron beam in a high density plasma via an electrostatic probe

Here, an electron beam is detected by a 1D floating potential probe array in a relatively high density (10 12 –10 13 cm -3 ) and low temperature (~5 eV) plasma of the Magnetic Reconnection Experiment. Clear perturbations in the floating potential profile by the electron beam are observed. Based on the floating potential profile and a current balance equation to the probe array tips, the effective width of the electron beam is determined, from which we determine the radial and toroidal beam current density profiles. After the profile of the electron beam is specified from the measured beam current, we demonstrate the consistency of the current balance equation and the location of the perturbation is also in agreement with field line mapping. No significant broadening of the electron beam is observed after the beam propagates for tens of centimeters through the high density plasma. These results prove that the field line mapping is, in principle, possible in high density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of energetic particle phase space modifications by instabilities on integrated modeling

Tokamak plasmas can feature a large population of energetic particles (EP) from neutral beam injection or fusion reactions. In turn, energetic particles can drive instabilities, which affect the driving EP population leading to a distortion of the original EP distribution function and of quantities that depend on it. The latter include, for example, neutral beam (NB) current drive and plasma heating through EP thermalization. Those effects must be taken into account to enable reliable and quantitative simulations of discharges for present devices as well as predictions for future burning plasmas. Reduced models for EP transport are emerging as an effective tool for long time-scale integrated simulations of tokamak plasmas, possibly including the effects of instabilities on EP dynamics. Available models differ in how EP distribution properties are modified by instabilities, e.g. in terms of gradients in real or phase space. It is therefore crucial to assess to what extent different assumptions in the transport models affect predicted quantities such as EP profile, energy distribution, NB driven current and energy/momentum transfer to the thermal populations. A newly developed kick model, which includes modifications of the EP distribution by instabilities in both real and velocity space, is used in this work to investigate these issues. Coupled to TRANSP simulations, the kick model is used to analyze NB-heated NSTX and DIII-D discharges featuring unstable Alfvén eigenmodes (AEs). Results show that instabilities can strongly affect the EP distribution function, and modifications propagate to macroscopic quantities such as NB-driven current profile and NB power transferred to the thermal plasma species. Furthermore, those important aspects are only qualitatively captured by simpler fast ion transport models that are based on radial diffusion of energetic ions only.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CGYRO

CGYRO is a multi-species spectral gyrokinetic solver with sonic rotation capability. It is an Eulerian GK solver specifically designed and optimized for collisional, electromagnetic, multiscale simulation. A key algorithmic aspect of CGYRO is the radially spectral formulation used to reduce the complicated integral gyroaveraging kernel into a multiplication in wavenumber space, but retaining the ability to treat profile variation important for edge plasmas. This makes it suitable for next-generation computational systems that require high levels of parallel concurrency.

Belli, EmilyA↗

Strongly intermittent far scrape-off layer fluctuations in Alcator C-Mod plasmas close to the empirical discharge density limit

Abstract Intermittent plasma fluctuations in the boundary region of the Alcator C-Mod device were comprehensively investigated using data time-series from gas puff imaging and mirror Langmuir probe diagnostics. Fluctuations were sampled during stationary plasma conditions in ohmically heated, lower single null diverted configurations with scans in both line-averaged density and plasma current, with Greenwald density fractions up to 0.85. Utilizing a stochastic model, we describe the plasma fluctuations as a super-position of uncorrelated pulses, with large-amplitude events corresponding to blob-like filaments moving through the scrape-off layer (SOL). A deconvolution method is used to estimate the pulse arrival times and amplitudes. The analysis reveals a significant increase of pulse amplitudes and waiting times as the line-averaged density approaches the empirical discharge density limit. Broadened and flattened average radial profiles are thus accompanied by strongly intermittent and large-amplitude fluctuations. Although these filaments are arriving less frequently at high line-averaged densities, we show that there are significant increases in radial far-SOL particle and heat fluxes which will further enhance plasma–wall interactions. The stochastic model has been used as a framework for study of the scalings in the intermittency parameter, flux and mean amplitude and waiting times, and is being used to inform predictive capability for the effects of filamentary transport as a function of Greenwald fraction.

Physics↗

THE TOPOLOGY OF CANONICAL FLUX TUBES IN FLARED JET GEOMETRY

Magnetized plasma jets are generally modeled as magnetic flux tubes filled with flowing plasma governed by magnetohydrodynamics (MHD). We outline here a more fundamental approach based on flux tubes of canonical vorticity, where canonical vorticity is defined as the circulation of the species’ canonical momentum. This approach extends the concept of magnetic flux tube evolution to include the effects of finite particle momentum and enables visualization of the topology of plasma jets in regimes beyond MHD. A flared, current-carrying magnetic flux tube in an ion-electron plasma with finite ion momentum is thus equivalent to either a pair of electron and ion flow flux tubes, a pair of electron and ion canonical momentum flux tubes, or a pair of electron and ion canonical vorticity flux tubes. We examine the morphology of all these flux tubes for increasing electrical currents, different radial current profiles, different electron Mach numbers, and a fixed, flared, axisymmetric magnetic geometry. Calculations of gauge-invariant relative canonical helicities track the evolution of magnetic, cross, and kinetic helicities in the system, and show that ion flow fields can unwind to compensate for an increasing magnetic twist. The results demonstrate that including a species’ finite momentum can result in a very long collimated canonical vorticity flux tube even if the magnetic flux tube is flared. With finite momentum, particle density gradients must be normal to canonical vorticities, not to magnetic fields, so observations of collimated astrophysical jets could be images of canonical vorticity flux tubes instead of magnetic flux tubes.

79 ASTRONOMY AND ASTROPHYSICS↗

Status of the design of the ITER ECE diagnostic

In this study, the baseline design for the ITER electron cyclotron emission (ECE) diagnostic has entered the detailed preliminary design phase. Two plasma views are planned, a radial view and an oblique view that is sensitive to distortions in the electron momentum distribution near the average thermal momentum. Both views provide high spatial resolution electron temperature profiles when the momentum distribution remains Maxwellian. The ECE diagnostic system consists of the front-end optics, including two 1000 K calibration sources, in equatorial port plug EP9, the 70-1000 GHz transmission system from the front-end to the diagnostics hall, and the ECE instrumentation in the diagnostics hall. The baseline ECE instrumentation will include two Michelson interferometers that can simultaneously measure ordinary and extraordinary mode ECE from 70 to 1000 GHz, and two heterodyne radiometer systems, covering 122-230 GHz and 244-355 GHz. Significant design challenges include 1) developing highly-reliable 1000 K calibration sources and the associated shutters/mirrors, 2) providing compliant couplings between the front-end optics and the polarization splitter box that accommodate displacements of the vacuum vessel during plasma operations and bake out, 3) protecting components from damage due to stray ECH radiation and other intense millimeter wave emission and 4) providing the low-loss broadband transmission system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detection of an electron beam in a high density plasma via an electrostatic probe

An electron beam is detected by a 1D floating potential probe array in a relatively high density (10e12 - 10e13 cm-3) and low temperature (~ 5 eV) plasma of the Magnetic Reconnection Experiment (MRX). Clear perturbations in the floating potential profile by the electron beam are observed. Based on the floating potential profile and a current balance equation to the probe array tips, the effective width of the electron beam is determined, from which we determine the radial and toroidal beam current density profiles. After the profile of the electron beam is specified from the measured beam current, we demonstrate the consistency of the current balance equation and the location of the perturbation is also in agreement with field line mapping. No significant broadening of the electron beam is observed after the beam propagates for tens of centimeters through the high density plasma. These results prove that the field line mapping is, in principle, possible in high density plasmas.

field line mapping↗

Divertor detachment in the pre-fusion power operation phase in ITER during application of resonant magnetic perturbations

Detachment of the divertor plasma during application of resonant magnetic perturbation (RMP) fields is evaluated for hydrogen H-mode plasma during the first pre-fusion power operation (PFPO-1) phase in ITER by 3D plasma boundary modelling with EMC3–EIRENE. Plasma response effects from a linearized, resistive, single fluid MHD model are discussed, which includes partial screening of the externally applied field—but also field amplification near the separatrix. Furthermore, this field amplification is found to play a pivotal role for the magnetic footprint on the divertor targets, but is sensitive to model parameters. Extensions of the footprint beyond the straight portions of the ITER vertical divertor targets, optimized for high stationary heat flux handling, may be possible depending on the level of toroidal rotation in the plasma. Exhaust from the bulk plasma is guided by the helical corrugations (lobes) of the perturbed separatrix, and this results in an upstream heat flux that is distributed over these lobes with lower peak values than in the typical radial heat flux profiles seen in the absence of magnetic perturbations. As a consequence, an earlier onset (with respect to the upstream density) of detachment is found in the traditional strike zone when RMPs are applied, but secondary, non-axisymmetric strike locations appear—and those remain attached at temperatures above 10 eV. Neon seeding can mitigate these non-axisymmetric heat loads, but this becomes less efficient for large magnetic footprints.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inference of experimental radial impurity transport on Alcator C-Mod: Bayesian parameter estimation and model selection

We present a fully Bayesian approach for the inference of radial profiles of impurity transport coefficients and compare its results to neoclassical, gyrofluid and gyrokinetic modeling. Using nested sampling, the Bayesian Impurity Transport InferencE (BITE) framework can handle complex parameter spaces with multiple possible solutions, offering great advantages in interpretative power and reliability with respect to previously demonstrated methods. BITE employs a forward model based on the pySTRAHL package, built on the success of the well-known STRAHL code [Dux, IPP Report, 2004], to simulate impurity transport in magnetically-confined plasmas. In this paper, we focus on calcium (Ca, Z=20) Laser Blow-Off injections into Alcator C-Mod plasmas. Multiple Ca atomic lines are diagnosed via high-resolution X-ray Imaging Crystal Spectroscopy and Vacuum Ultra-Violet measurements. We analyze a sawtoothing I-mode discharge for which neoclassical and turbulent (quasilinear and nonlinear) predictions are also obtained. Here, we find good agreement in diffusion across the entire radial extent, while turbulent convection and density profile peaking are estimated to be larger in experiment than suggested by theory. Efforts and challenges associated with the inference of experimental pedestal impurity transport are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An optical-input Maximum Likelihood Estimation feedback system demonstrated on tokamak horizontal equilibrium control

A readily parallelized Maximum Likelihood Estimation (MLE) algorithm with linear computational complexity is demonstrated in real time using only measurements from an extreme ultraviolet (EUV) diagnostic to control the horizontal position of a tokamak plasma. A set of trial emissivity profiles are parameterized by the control quantity of interest (R m ), and the MLE is identified from the profile which minimizes the signal reconstruction residual. The algorithm depends on an empirically determined likelihood function with exponential form. EUV emission (λ ≈ 15eV-1keV) is captured in a poloidal plane by four 16-channel AXUV diodes mounted at different poloidal angles with radial and angular resolution sufficient to discern plasma equilibrium evolution in HBT-EP. Calculations of the plasma major radius by the system are consistent within diagnostic uncertainty for the majority of the discharge with those of: a weighted average of vertical soft X-ray or EUV chords, magnetic sensors, and an equilibrium reconstruction. The feedback system corrects for a horizontal displacement of the major radius equal to 20% of the plasma minor radius by adjusting the vertical field produced from 40 in-vessel control coils in real time. The MLE calculation is performed on a GPU in a 15 μs cycle, with similar performance in this application to a simple weighted average of vertical chords. Finally, results demonstrate horizontal position control using magnetic actuators and an optical observer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotopic mass effects of tritium-fueled high-performance TFTR supershots

An increase in total stored energy correlated with the addition of tritium fuel was observed in supershots during the TFTR DT campaign. This supershot regime had strikingly high, centrally peaked ion and electron temperatures, and the largest neutron emission rates observed in TFTR. This paper presents a study of the causes of this increase in stored energy in supershots. Twenty-six supershots have been recently reanalyzed with the TRANSP plasma analysis code. Early TRANSP simulations did not accurately match the measured magnitude and time evolutions of the neutron emission rates. This mismatch is attributed to neglecting apparent increases of trace amounts of heavy impurities during neutral beam injection. The new TRANSP runs were tuned to accommodate this and match the measured global neutron emission rates. These new runs also had improved fidelity in predicting the time histories and radial dependencies of measured DT neutron emission rates. That in turn adds confidence in the simulated thermal deuterium and tritium density profiles that are needed for calculating the average hydrogenic atomic mass profiles. Six subsets of these supershots had well matched toroidal field B tor , plasma current I p , flux geometry, and total injected neutral beam power. The mix of D and T beam ions was varied for different discharges. The magnitude of the increase of the thermal ion energy W i with added tritium was relatively small, and the total W tot increase is dominated by the increase in fast beam ions with T. Analyses at times before the occurrences of deleterious MHD instabilities yielded scaling of W tot with the volume-average isotopic mass consistent with previous publications. The relative fraction of fast energy ions is expected to be small in practical tokamak reactors. Thus the increase in stored energy W tot observed in TFTR supershots does not appear likely to be significantly helpful for producing useful fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global gyrokinetic stability of collisionless microtearing modes in large aspect ratio tokamaks

Linear full radius gyrokinetic calculations show the existence of unstable microtearing modes (MTMs) in purely collisionless, high temperature, large aspect ratio tokamak plasmas. The present study takes into account fully gyrokinetic highly passing ions and electrons. The global 2-D structures of the collisionless mode with full radius coupling of the poloidal modes is obtained and compared with another electromagnetic mode, namely, the Alfvén Ion Temperature Gradient (AITG) mode (or Kinetic Ballooning Mode, KBM) for the same equilibrium profile. Several important characteristics of the modes are brought out and compared, such as a clear signature in the symmetry properties of the two modes, the plasma–β dependence, and radial and poloidal length scales of the electrostatic and magnetic vector potential fluctuations. Extensive parameter scans for this collisionless microtearing mode reveal the scaling of the growth rate with β and the electron temperature gradient η{sub e}. Scans at different β values show an inverse relationship between the η{sub e} threshold and β, leading to a stability diagram, and implying that the mode might exist at moderate to strong temperature gradients for finite β plasmas in large aspect ratio tokamaks. In contrast to small aspect ratio tokamaks where the trapped electron magnetic drift resonance is found to be important, in large aspect ratio tokamaks, a strong destabilization due to the magnetic drift resonance of passing electrons is observed and is identified as a possible collisionless drive mechanism for the collisionless MTM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Discharge Characterization of 40 cm-Microwave ECR Ion Source and Neutralizer

Discharge characteristics of a 40 cm, 2.45 GHz Electron Cyclotron Resonance (ECR) ion thruster discharge chamber and neutralizer were acquired. Thruster bulk discharge plasma characteristics were assessed using a single Langmuir probe. Total extractable ion current was measured as a function of input microwave power and flow rate. Additionally, radial ion current density profiles at the thruster.s exit plane were characterized using five equally spaced Faraday probes. Distinct low and high density operating modes were observed as discharge input power was varied from 0 to 200 W. In the high mode, extractable ion currents as high as 0.82 A were measured. Neutralizer emission current was characterized as a function of flow rate and microwave power. Neutralizer extraction currents as high as 0.6 A were measured.

Foster, John E.↗

The motional Stark effect diagnostic for ITER

An overview of the plans for the motional Stark effect (MSE) diagnostic installation on the International Thermonuclear Experimental Reactor (ITER) is presented. The MSE diagnostic uniquely provides spatially localized magnetic field measurements inside the plasma. These are used to constrain equilibrium reconstructions to determine q(r), the safety factor as a function of minor radius. Meeting the system requirements to deliver q-profiles and related quantities with the specified radial resolution of 20 points over the minor radius, 10 ms time resolution, and better than 10% accuracy is challenging. MSE systems observe the D/H-α emission near 656.3 nm from neutral beams. As the beam atoms traverse the magnetic field, B⃗, at high velocity, v⃗, they experience a Lorentz electric field, v⃗×B⃗, which causes the spectral emission to be split and polarized due to the Stark effect. Traditional MSE-LP (line polarization) measurements determine the direction of the magnetic field in the observation volume using polarimetric analysis of the detected light. The harsh conditions of ITER are expected to deposit thin films of contaminants on the first mirror, which would alter the polarization state of reflected light significantly. On ITER, the combination of high magnetic field strength and high energy beams makes the Stark spectrum resolution suitable for the determination of the magnetic field magnitude from the line shift, so this approach has been selected. Every aspect of the measurement system must be planned for the burning plasma environment and carefully analyzed ahead of time. Current status and plans for the system are presented.

Instruments & Instrumentation↗

Finite-difference multiple fluid solution for source-driven rotation in highly magnetized linear plasma device

The rotation profile of a magnetized plasma cylinder composed of multiple fluids is investigated analytically, expanding on previous results. The analytic steady-state solution is used as a benchmark for a time-dependent multiple-fluid finite-difference code, MITNS: Multiple-Ion Transport Numerical Solver. Magnetic field evolution is taken into account, both analytically and numerically. Its details are shown to be of importance when particles are allowed out of the domain. MITNS reproduces the asymptotic expansion results for a small parameter $δ<<<1$. For $\sqrt{m_elm_i} ~ δ << 1$ , a slightly different regime, dominated by viscosity-induced transport of ions, is found numerically and analytically. Finally, this verification supports the use of this code for more complex time-dependent calculations in the future. Additionally, we derive the angular velocity profile of each species due to radial particle and charge fluxes of various strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bohm-criterion approximation versus optimal matched solution for a cylindrical probe in radial-motion theory

The theory of positive-ion collection by a probe immersed in a low-pressure plasma was reviewed and extended by Allen et al. [Proc. Phys. Soc. 70, 297 (1957)]. The numerical computations for cylindrical and spherical probes in a sheath region were presented by F. F. Chen [J. Nucl. Energy C 7, 41 (1965)]. Here, in this paper, the sheath and presheath solutions for a cylindrical probe are matched through a numerical matching procedure to yield “matched” potential profile or “M solution.” The solution based on the Bohm criterion approach “B solution” is discussed for this particular problem. The comparison of cylindrical probe characteristics obtained from the correct potential profile (M solution) and the approximated Bohm-criterion approach are different. This raises questions about the correctness of cylindrical probe theories relying only on the Bohm-criterion approach. Also the comparison between theoretical and experimental ion current characteristics shows that in an argon plasma the ions motion towards the probe is almost radial.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗