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 559 records · Page 31

PLASMA INSTABILITIES IN THE CONTEXT OF CURRENT HELIUM SEDIMENTATION MODELS: DYNAMICAL IMPLICATIONS FOR THE ICM IN GALAXY CLUSTERS

Understanding whether Helium can sediment to the core of galaxy clusters is important for a number of problems in cosmology and astrophysics. All current models addressing this question are one-dimensional and do not account for the fact that magnetic fields can effectively channel ions and electrons, leading to anisotropic transport of momentum, heat, and particle diffusion in the weakly collisional intracluster medium (ICM). This anisotropy can lead to a wide variety of instabilities, which could be relevant for understanding the dynamics of heterogeneous media. In this paper, we consider the radial temperature and composition profiles as obtained from a state-of-the-art Helium sedimentation model and analyze its stability properties. We find that the associated radial profiles are unstable to different kinds of instabilities depending on the magnetic field orientation at all radii. The fastest growing modes are usually related to generalizations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability which operate in heterogeneous media. We find that the effect of sedimentation is to increase (decrease) the predicted growth rates in the inner (outer) cluster region. The unstable modes grow quickly compared to the sedimentation timescale. This suggests that the composition gradients as inferred from sedimentation models, which do not fully account for the anisotropic character of the weakly collisional environment, might not be very robust. Our results emphasize the subtleties involved in understanding the gas dynamics of the ICM and argue for the need of a comprehensive approach to address the issue of Helium sedimentation beyond current models.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental observation of the plasma potential with the screening temperature

An experimental investigation of the Boltzmann relation in the plasma bulk region and the presheath region was performed at low pressure inductively coupled plasma with bi-Maxwellian distribution. The spatial profiles of the electron energy distribution function (EEDF) were measured from the discharge center to the plasma boundary with an ion mean free path scale. The screening temperature was derived from the Einstein relation and the bi-Maxwellian EEDFs, and it has the same expression as the electron temperature in the ion directed velocity to the sheath edge in plasmas with the bi-Maxwellian EEDFs. Radial variations of the ambipolar potential and the presheath potential were calculated with the screening temperature and compared with the obtained value from the spatially measured EEDFs. It was shown that the ambipolar potential and the presheath potential are governed by the screening temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of neutral interactions on parallel transport and blob dynamics in gyrokinetic scrape-off layer simulations

The effect of neutral interactions on scrape-off layer (SOL) turbulence is investigated in a continuum gyrokinetic code that has been coupled to a continuum kinetic model of neutral transport. This extends the work of a previous paper, which compared two NSTX SOL simulations in simple helical geometry, one with neutrals and one without. The former included electron-impact ionization, charge exchange, and wall recycling. Here, the case with neutrals is compared to a gyrokinetic-only simulation that includes an effective ionization source to separate the effect of sourcing from charge exchange collisions. It is observed that sourcing accounts for many features of the simulated SOL with neutrals, including density and temperature magnitudes and reduced normalized density fluctuations, but differences persist. In particular, a flatter density profile results due to changes in parallel transport when neutral collisions are included, illustrating the importance of neutral drag on global plasma properties. An analysis of coherent turbulent structures, or blobs, in these simulations demonstrates the case with neutrals has slower and larger blobs. Here, a series of seeded blob simulations corroborates the blob velocity observation. In general, the blob motion does not contribute significantly to radial transport in these simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resonance broadened quasi-linear (RBQ) model for fast ion distribution relaxation due to Alfvénic eigenmodes

The burning plasma performance is limited by the confinement of the super-alfvénic fusion products such as alpha particles and the auxiliary heating ions capable of exciting the Alfvénic eigenmodes (AEs) (Gorelenkov et al 2014 Nucl. Fusion 54 125001). Here in this work the effect of AEs on fast ions is formulated within the quasi-linear (QL) theory generalized for this problem recently (Duarte 2017 PhD Thesis University of São Paulo, Brazil). The generalization involves the resonance line broadened interaction of energetic particles (EP) with AEs supplemented by the diffusion coefficients depending on EP position in the velocity space. A new resonance broadened QL code (or RBQ1D) based on this formulation allowing for EP diffusion in radial direction is built and presented in details. In RBQ1D applications we reduce the wave particle interaction (WPI) dynamics to 1D case when the particle kinetic energy is nearly constant. The diffusion equation for EP distribution evolution is then solved simultaneously for all particles along the angular momentum direction. We make initial applications of the RBQ1D to a DIII-D plasma with elevated q-profile where the beam ions show stiff transport properties (Collins et al (The DIII-D Team) 2016 Phys. Rev. Lett. 116 095001). In conclusion, AE driven fast ion profile relaxation is studied for validations of the QL approach in realistic conditions of beam ion driven instabilities in DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The intrinsic H I Lyman-alpha line profiles of late-type stars

The Lyman-alpha line of neutral hydrogen is probably the most important cooling channel for chromospheric plasma in late-type stars, yet it is also the least studied major line in the far ultraviolet. The scattering of much of the stellar Lyman-alpha flux by interstellar hydrogen, coupled with the geocoronal emission foreground, seriously complicates the analysis of the Lyman-alpha spectra. The influence of the local interstellar medium on the observed profiles was circumvented by observing stars with radial velocities sufficiently high to Doppler shift the center of the stellar emission line out of the interstellar absorption core. There are several stars that have high radial velocities by virtue of their presence in close binary systems. High resolution IUE (International Ultraviolet Explorer) spectra of Ly alpha line of two such eclipsing binary stars, AR Lac and TY Pyx, are obtained, at each orbital quadrature phase, when the projected orbital velocity is a maximum. By combining the spectra from opposite quadratures it is possible to piece together the entire stellar emission profiles. The third star in this study, delta Lep, is a single star with a high space velocity.

Neff, J. E.↗

Diagnosing the pedestal magnetic field and magnetohydrodynamics radial structure with pedestal–scrape of layer electron cyclotron emission radiation inversion in H-mode plasma (invited)

Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of plasma density/collisionality on divertor heat flux width

Both ASDEX-Upgrade (AUG) data and the generalized HD (GHD) model showed that the scrape-off width broadens as the density/collisionality increases [1, 2]. A series of BOUT++ transport simulations are performed to study the physics of the scaling characteristics of the divertor heat flux width vs density/collisionality via a plasma density scan with either fixed pressure profile or fixed temperature profile inside separatrix. Additionally, the simulations show that even in the drift dominated regime, the divertor heat flux width can be broadened due to the transition of the SOL residence time from the parallel particle flow time to the enhanced parallel conduction time as the collisionality/density increases as posited in the GHD model. In addition, the heat flux width is found to be proportional to the square root of ion mass for low collisionality while it has a weakly dependence on ion mass for high collisionality. Furthermore, our simulations show that as the density increases, the radial electric field (E r ) well shallows, which potentially weakens E r × B flow shear stabilization of turbulence at high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma edge and scrape-off layer turbulence in gyrokinetic simulations of negative triangularity plasmas

Gyrokinetic simulations in the long-wavelength or drift-kinetic limit are carried out of DIII-D inner-wall-limited (IWL) plasmas to investigate the effect of triangularity on edge and scrape-off layer (SOL) turbulence. The effect of neutral interactions and triangularity on plasma blobs is explored due to the impact blobs can have in setting the SOL width or introducing impurities through interactions with plasma-facing components. Seeded blob simulations with neutrals in shaped SOL scenarios demonstrate that increasing elongation, triangularity, or Shafranov shift decreases radial blob velocities, but neutral interactions have a minor effect. Fully turbulent simulations of DIII-D IWL plasmas include both open- and closed-field-line regions. The negative triangularity (NT) simulation has lower average core T e , lower normalized T e fluctuations, and lower fluxes, but a greater number of coherent structures (blobs) identified with increased size and velocity, on average. Density and electron temperature profiles are within a factor of 2 of experimental values. In conclusion, the increased trapped electron particle fraction in NT simulations is consistent with previous studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modelling enhanced confinement in drift-wave turbulence

Here, the results of modeling studies of an enhanced confinement in the drift wave turbulent plasma of the CSDX linear device are presented. The mechanism of enhanced confinement is investigated here using a reduced 1D, time-dependent model, which illustrates the exchange of enstrophy between two disparate scale structures: the mesoscale flow and profile, and the turbulence intensity fields. Mean density, mean vorticity, and turbulent potential enstrophy are the variables for this model. Total potential enstrophy is conserved in this model. Vorticity mixing occurs on a scale length related to an effective Rhines' scale of turbulence, and shrinks as both density and vorticity gradients steepen. Numerical results obtained from solution of the model agree well with the experimental data from CSDX showing: (i) a steepening of the mean density profile, indicating a radial transport barrier formation, (ii) the development of a radially sheared azimuthal flow velocity that coincides with the density steepening and initiates a turbulence quench, and (iii) negative Reynolds work values, indicating that fluctuations drive the shear flow. These observations as the magnitude of the magnetic field B increases are recovered using purely diffusive expressions for the vorticity and density fluxes. A new dimensionless turbulence parameter R DT -defined as the ratio of the integrated potential enstrophy transfer from turbulence to the flow, to the integrated potential enstrophy production due to relaxation of the density gradient is introduced as a turbulence collapse indicator that detects when the enhanced confinement state is triggered

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comprehensive comparisons of geodesic acoustic mode characteristics and dynamics between Tore Supra experiments and gyrokinetic simulations

In a dedicated collisionality scan in Tore Supra, the geodesic acoustic mode (GAM) is detected and identified with the Doppler backscattering technique. Observations are compared to the results of a simulation with the gyrokinetic code GYSELA. We found that the GAM frequency in experiments is lower than predicted by simulation and theory. Moreover, the disagreement is higher in the low collisionality scenario. Bursts of non harmonic GAM oscillations have been characterized with filtering techniques, such as the Hilbert-Huang transform. When comparing this dynamical behaviour between experiments and simulation, the probability density function of GAM amplitude and the burst autocorrelation time are found to be remarkably similar. In the simulation, where the radial profile of GAM frequency is continuous, we observed a phenomenon of radial phase mixing of the GAM oscillations, which could influence the burst autocorrelation time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗

On the toroidal plasma rotations induced by lower hybrid waves

A theoretical model is developed to explain the plasma rotations induced by lower hybrid waves in Alcator C-Mod. In this model, torodial rotations are driven by the Lorentz force on the bulk-electron flow across flux surfaces, which is a response of the plasma to the resonant-electron flow across flux surfaces induced by the lower hybrid waves. The flow across flux surfaces of the resonant electrons and the bulk electrons are coupled through the radial electric field initiated by the resonant electrons, and the friction between ions and electrons transfers the toroidal momentum to ions from electrons. An improved quasilinear theory with gyrophase dependent distribution function is developed to calculate the perpendicular resonant-electron flow. Toroidal rotations are determined using a set of fluid equations for bulk electrons and ions, which are solved numerically by a finite-difference method. Numerical results agree well with the experimental observations in terms of flow profile and amplitude. The model explains the strong correlation between torodial flow and internal inductance observed experimentally, and predicts both counter-current and co-current flows, depending on the perpendicular wave vectors of the lower hybrid waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory and simulation of quasilinear transport from external magnetic field perturbations in a DIII-D plasma

The linear response profiles for the 3D perturbed magnetic fields, currents, ion velocities, plasma density, pressures, and electric potential from low-n external resonant magnetic field perturbations (RMPs) are obtained from the collisional two-fluid M3D-C1 code [N. M. Ferraro and S. C. Jardin, J. Comput. Phys. 228, 7742 (2009)]. A newly developed post-processing RMPtran code computes the resulting quasilinear E×B and magnetic (J×B) radial transport flows with respect to the unperturbed flux surfaces in all channels. RMPtran simulations focus on ion (center of mass) particle and transient non-ambipolar current flows, as well as the toroidal angular momentum flow. The paper attempts to delineate the RMP transport mechanisms that might be responsible for the RMP density pump-out seen in DIII-D [M. A. Mahdavi and J. L. Luxon, Fusion Sci. Technol. 48, 2 (2005)]. Experimentally, the starting high toroidal rotation does not brake to a significantly lower rotation after the pump-out suggesting that convective and E×B transport mechanisms dominate. The direct J×B torque from the transient non-ambipolar radial current expected to accelerate plasma rotation is shown to cancel much of the Maxwell stress J×B torque expected to brake the plasma rotation. The dominant E×B Reynolds stress accelerates rotation at the top of the pedestal while braking rotation further down the pedestal.

Waltz, R. E. (ORCID:0000000252034830)↗

Assessment of the strength of kinetic effects of parallel electron transport in the SOL and divertor of JET high radiative H-mode plasmas using EDGE2D-EIRENE and KIPP codes

The kinetic code for plasma periphery (KIPP) was used here to assess the importance of the kinetic effects of parallel electron transport in the scrape-off layer (SOL) and divertor of JET high radiative H-mode inter-ELM plasma conditions with the ITER-like wall and strong nitrogen (N 2 ) injection. Plasma parameter profiles along a magnetic field from one of the EDGE2D-EIRENE simulation cases were used as an input for KIPP runs. Profiles were maintained by particle and power sources. KIPP generated electron distribution functions, f e , parallel power fluxes, electron–ion thermoforces, Debye sheath potential drops and electron sheath transmission factors at divertor targets. For heat fluxes in the main SOL, KIPP results showed deviations from classical (e.g. Braginskii) fluxes by factors typically of ~1.5, sometimes up to 2, with the flux limiting for more upstream positions and flux enhancement near entrances to the divertor. In the divertor, at the same time, for radial positions closer to the separatrix, very large heat flux enhancement factors of up to ten or even higher, indicative of a strong nonlocal heat transport, were found at the outer target, with heat power flux density exhibiting bump-on-tail features at high energies. Under such extreme conditions, however, contributions of conductive power fluxes to total power fluxes were strongly reduced, with convective power fluxes becoming comparable, or sometimes exceeding, conductive power fluxes. Electron–ion thermoforce, on the other hand, which is known to be determined mostly by thermal and subthermal electrons, was found to be in good agreement with Braginskii formulas, including the Z eff dependence. Overall, KIPP results indicate, at least for the plasma conditions used in this modelling, a sizable, but not dominant, effect of kinetics on parallel electron transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of ion and electron heat transport of high- T e ECH heated discharges in the large helical device

An analysis of the radial electric field and heat transport, both for ions and electrons, is presented for a high-${{T}_{\text{e}}}$ electron cyclotron heated (ECH) discharge on the large helical device (LHD). Transport analysis is done using the task3d transport suite utilizing experimentally measured profiles for both ions and electrons. Ion temperature and perpendicular flow profiles are measured using the recently installed x-ray imaging crystal spectrometer diagnostic (XICS), while electron temperature and density profiles are measured using Thomson scattering. The analysis also includes calculated ECH power deposition profiles as determined through the travis ray-tracing code. This is the first time on LHD that this type of integrated transport analysis with measured ion temperature profiles has been performed without NBI, allowing the heat transport properties of plasmas with only ECH heating to be more clearly examined. For this study, a plasma discharge is chosen which develops a high central electron temperature (${{T}_{\text{eo}}}=9$ keV) at moderately low densities (${{n}_{\text{eo}}}=1.5\times {{10}^{19}}$ m-3). The experimentally determined transport properties from task3d are compared to neoclassical predictions as calculated by the gsrake and fortec-3d codes. The predicted electron fluxes are seen to be an order of magnitude less than the measured fluxes, indicating that electron transport is largely anomalous, while the neoclassical and measured ion heat fluxes are of the same magnitude. Neoclassical predictions of a strong positive ambipolar electric field (${{E}_{\text{r}}}$ ) in the plasma core are validated through comparisons to perpendicular flow measurements from the XICS diagnostic. Furthermore, this provides confidence that the predictions are producing physically meaningful results for the particle fluxes and radial electric field, which are a key component in correctly predicting plasma confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Full-f gyrokinetic simulation over a confinement time

A long time ion temperature gradient driven turbulence simulation over a confinement time is performed using the full-f gyrokinetic Eulerian code GT5D. The convergence of steady temperature and rotation profiles is examined, and it is shown that the profile relaxation can be significantly accelerated when the simulation is initialized with linearly unstable temperature profiles. In the steady state, the temperature profile and the ion heat diffusivity are self-consistently determined by the power balance condition, while the intrinsic rotation profile is sustained by complicated momentum transport processes without momentum input. The steady turbulent momentum transport is characterized by bursty non-diffusive fluxes, and the resulting turbulent residual stress is consistent with the profile shear stress theory [Y. Camenen et al., “Consequences of profile shearing on toroidal momentum transport,” Nucl. Fusion 51, 073039 (2011)] in which the residual stress depends not only on the profile shear and the radial electric field shear but also on the radial electric field itself. Based on the toroidal angular momentum conservation, it is found that in the steady null momentum transport state, the turbulent residual stress is cancelled by the neoclassical counterpart, which is greatly enhanced in the presence of turbulent fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interplay between fluctuation driven toroidal axisymmetric flows and resistive ballooning mode turbulence

An interplay between fluctuation driven toroidal axisymmetric flows (convective cell modes) and resistive ballooning mode turbulence after the pedestal collapse is numerically studied by a four-field reduced MHD model in the BOUT++ framework. The strong flow shear suppresses the radial transport of pressure filaments, and the pressure profile in the pedestal region is partially recovered. As a result, a secondary instability is quasilinearly excited, which yields a secondary collapse. Here, the subsequent damped oscillation is also analyzed by phase diagram analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel convection and E × B drifts in the TCV snowflake divertor and their effects on target heat-fluxes

Parallel convection and E × B drifts act together to redistribute heat between the strike-points in the low field side snowflake minus (LFS SF–). The cumulative heat convection from both mechanisms is enhanced near the secondary X-point and is shown to dominate over heat conduction, partly explaining why the LFS SF– distributes power more evenly than the single null (SN) or other snowflake (SF) configurations. Pressure profiles at the entrance of the divertor are strongly affected by the position of the secondary X-point and magnetic field direction indicating the importance of E × B drifts. Pressure drops of up to 50% appear between the outer-midplane (OMP) and the divertor entrance enhancing the role of parallel heat convection. The electron temperature and density profiles and the radial turbulent fluxes measured at the OMP are largely unaffected by the changes in divertor geometry, even on flux surfaces where the connection length is infinite.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗