Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma properties”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A generalized Boltzmann kinetic theory for strongly magnetized plasmas with application to friction

Coulomb collisions in plasmas are typically modeled using the Boltzmann collision operator, or its variants, which apply to weakly magnetized plasmas in which the typical gyroradius of particles significantly exceeds the Debye length. Conversely, O’Neil has developed a kinetic theory to treat plasmas that are so strongly magnetized that the typical gyroradius of particles is much smaller than the distance of the closest approach in a binary collision. Here, we develop a generalized collision operator that applies across the full range of magnetization strength. Since there is no closed-form solution for the scattering cross section when plasma is strongly magnetized, the input to the collision operator is obtained by numerically calculating particle trajectories. To demonstrate novel physics associated with strong magnetization, it is used to compute the friction force on a massive test charge. In addition to the traditional stopping power component, this is found to exhibit a transverse component that is perpendicular to both the velocity and Lorentz force vectors in the strongly magnetized regime, as was predicted recently using linear response theory. As a result, good agreement is found between the collision theory and linear response theory in the regime in which both apply, but the new collision theory also applies to stronger magnetization strength regimes than the linear response theory is expected to apply in.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Asymmetry between deuterium and tritium turbulent particle flows

The asymmetry between deuterium (D) and tritium (T) turbulent particle fluxes in mixed D–T plasmas is studied with numerical simulations of nonlinear gyrokinetic turbulence in ion temperature gradient-dominated and trapped electron mode-dominated regimes. At 50–50 D–T concentration, the asymmetry, or flow separation, between D and T fluxes is such that the tritium is better confined than the deuterium in both regimes. To supplement the nonlinear simulations, an analytic quasilinear theory of the particle flux symmetry breaking is developed and is valid for general electron dynamics. Furthermore, this theory correctly predicts the numerically computed deviation in the ion density fraction (from 50–50), or ion density gradient (from the electron gradient), required to restore equal deuterium and tritium fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

RF sheath induced sputtering on Proto-MPEX part 2: Impurity transport modeling and experimental comparison

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a pulsed linear plasma device that uses a radio frequency (RF) helicon antenna and a helicon ceramic window as the plasma ionization source. Modeling performed within this work on the helicon operations predicts hot spots on the helicon ceramic window due to an RF-induced sheath. The RF sheath potential leads to physical sputtering of the ceramic window predominantly where the helicon antenna sits around the helicon window. The role of oxygen in sputtering within Proto-MPEX is also investigated and is needed to match the experimental data. Experiments on Proto-MPEX show reasonable agreement with the modeled erosion/redeposition pattern on the helicon window, and the modeled impurity flux to the target shows radial transport within the device as expected from classical collisional impurity transport. Temperature screening of impurities is shown to be needed for Proto-MPEX to reduce impurity flux to the target. The modeling approach used in this work can be applied to other plasma RF applications for the determination of impurity production and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fueling limits in a cylindrical viscosity-limited reactor

Recently, a method to achieve a “natural hot-ion mode” was suggested by utilizing ion viscous heating in a rotating plasma with a fixed boundary. We explore the steady-state solution to the Braginskii equations and find the parameter regime in which a significant temperature difference between ions and electrons can be sustained in a driven steady state. Here, the threshold for this effect occurs at ρ i ≳ 0.1R. An analytic, leading order low flow solution is obtained, and a numerical, moderate Mach number M ≲ 2 is investigated. The limitation is found to be at moderate Mach numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Use of intrinsic hysteresis for the active control of internal transport barriers in magnetically confined fusion plasmas

In magnetically confined fusion devices, control of internal transport barriers (ITBs) is important both to enhance and suppress the turbulent transport to improve confinement control. Barrier control should allow for the improvement of confinement to aid in achieving the needed fusion criteria while also permitting the degradation of confinement to control profiles and clean the device by moving out the impurities accumulated near the core. In this work, we present a novel control scenario that takes advantage of the hysteresis intrinsic to transport barriers to easily cycle through enhanced and degraded confinement regimes. The control scenario is illustrated using a five-field simplified transport model for an ITB using typical parameters of a neutral beam injection-heated DIII-D tokamak discharge. Pellets and ion cyclotron resonance frequency power are used as control knobs for this active control scenario. These knobs adequately modify at will the local gradients and, therefore, the growth rates and shearing rates, allowing for easy and efficient control of the barrier by taking advantage of the barrier hysteresis. The result is a control cycle that could be operated with a relatively small amount of power in high performance regimes which, nowadays, typically require large power to control. It may also have advantages to avoid, or at least ameliorate, the appearance of magnetohydrodynamic instabilities in the barrier region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Preface for the Sherwood Fusion Theory 2021–2022 special collection

The Sherwood Fusion Theory conference has a 60-year history and has been held nearly annually for at least five decades. With theoretical and computational research directed toward fusion energy as an organizing principle, the typical Sherwood program encompasses a wide range of plasma physics topics pertaining to a variety of confinement concepts along with advancements in algorithms and numerical methods for fusion plasma simulations. Additionally, the typical annual meeting comprises roughly a dozen invited talks plus a small number of plenary speakers and several poster sessions. The program committee presents several awards each year to graduate (or undergraduate) students with outstanding poster presentations. Beginning in 2021, both invited speakers and student poster award winners were asked to submit manuscripts for this special collection, which combines contributions from the 2021 and 2022 Sherwood meetings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of a central barrier potential on asymmetry-induced transport in a coaxial Malmberg–Penning trap

Experiments are presented measuring the radial particle flux produced when a central barrier potential is applied to the central ring of a coaxial Malmberg-Penning trap at the same time as the voltages producing the main asymmetry potential $\phi_1(r)\cos{(kz)}\cos{(l\theta-\omega t)}$. When a negative DC voltage is applied to the entire central ring, the peak flux produced by the main asymmetry is reduced by a factor $e^{(V_{eff}/1.2)}$ where $V_{eff}$ is the effective applied potential in volts. This barrier potential also increases the frequency of this peak flux. When asymmetric $\pm$ voltages are applied to the two azimuthally-divided halves of the central ring, DC voltages increase the radial flux while AC voltages decrease it. It is shown that these asymmetric barrier voltages produce their own transport and that the resulting modification of the plasma may be a factor in changing the observed flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical implications of including drifts in SOLPS-ITER simulations of EAST

The inclusion of drifts in plasma edge codes like SOLPS-ITER is required to match simulation data with experimental profiles. However, this remains numerically challenging. In this paper, the effect of some numerical factors on the final plasma solution is investigated. This study is performed on three EAST simulations in the upper single null configuration: an attached purely deuterium case, an attached case with limited Ne-seeding, and a detached Ne-seeded case. The effects of the anomalous conductivity and anomalous thermo-electric coefficient on the plasma potential are investigated. Next, the effect of the employed grids is shown. In order to investigate these effects, accurate drift simulations are needed. Therefore, the employed time step and numerical parameters are discussed for the three studied simulations. For all presented simulations, it is verified that the restriction of the grid to the first flux surface tangent to the main chamber wall is sufficient within the context of non-extended simulations. This means that the main power dissipation takes place inside the simulated domain, and only a small fraction of the power is leaving the B2.5 grid through the grid boundary closest to the first wall. Finally, the effect of drifts on the asymmetry between the inner and outer target for EAST simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quantifying fission yields at the National Ignition Facility using depleted uranium foil experiments

There are programs for high-Z shell experiments at the National Ignition Facility (NIF). For shells made of actinide material, a quantitative fission diagnostic is needed in order to determine how much fission took place and whether the fission was sufficient to produce a non-negligible heat source in the burning capsule. Here, we present a viable coupled experimental and theoretical technique for making quantitative fission measurements possible. The proposed scheme involves using a well-characterized set of depleted uranium foils outside an NIF capsule to verify the conversion of xenon and krypton fission fragments collected at the Radiochemical Analysis of Gaseous Samples (RAGS) facility into total fission yield. We present the calculations needed for this conversion, including the decays in and out of fission fragment chains during the RAGS pump-down of the NIF chamber.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗
DOE OSTI2025

On transport of electrons in a magnetized plasma sheath

Electron trajectories are analyzed in the plasma sheath region in a magnetic field that intersects the surface at a shallow angle. Depending on the size of the Debye width relative to the electron gyroradius, λ D /ρ e , an electron either drifts in the sheath electric field or follows a sequence of reflections (“bouncing”) on the surface before reversing its parallel velocity. For thin sheath, λ D /ρ e ≪ 1, the conservation of the electron magnetic moment can be violated, depending on the ratio of the electron perpendicular and parallel velocities V ⊥ /V || and the magnetic field pitch angle with respect to the surface. Small-scale waviness of the surface can strongly amplify this effect, leading to violation of the magnetic moment conservation even for λ D /ρ e ≳ 1, which opens a channel of anomalous radial transport for electrons in plasma devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗
DOE OSTI2020

Cold-hot coupled waves in a flowing magnetized plasma

Abstract Nonlinear coupling of cold and hot waves in a flowing magnetized plasma is analyzed with the Vlasov equation. An analytical solution is obtained for cold waves of a small amplitude (weak flow) and a long wavelength. The distribution function is obtained by integrating the kinetic equation along a perturbed phase-space trajectory for a time-varying plasma flow. The kinetic description presents a generalized dispersion relation that involves resonances depending on cold and hot wave dispersions. Coherent fluid motion leads to radiation peaks in addition to the cyclotron harmonics, where the wavenumber of the cold wave determines the peak frequencies. The peaks appear narrow when the wave propagates perpendicular to the time-averaged flow while they become broad due to the Doppler effect when the wave propagates parallel to the flow. Fully kinetic particle-in-cell simulations corroborate the theoretical predictions. The dispersion relation and resulting wave spectra provide information about plasma parameters and flow properties.

Physics↗
DOE OSTI2022

Predicting resistive wall mode stability in NSTX through balanced random forests and counterfactual explanations

Abstract Recent progress in the disruption event characterization and forecasting framework has shown that machine learning guided by physics theory can be easily implemented as a supporting tool for fast computations of ideal stability properties of spherical tokamak plasmas. In order to extend that idea, a customized random forest (RF) classifier that takes into account imbalances in the training data is hereby employed to predict resistive wall mode (RWM) stability for a set of high beta discharges from the NSTX spherical tokamak. More specifically, with this approach each tree in the forest is trained on samples that are balanced via a user-defined over/under-sampler. The proposed approach outperforms classical cost-sensitive methods for the problem at hand, in particular when used in conjunction with a random under-sampler, while also resulting in a threefold reduction in the training time. In order to further understand the model’s decisions, a diverse set of counterfactual explanations based on determinantal point processes (DPP) is generated and evaluated. Via the use of DPP, the underlying RF model infers that the presence of hypothetical magnetohydrodynamic activity would have prevented the RWM from concurrently going unstable, which is a counterfactual that is indeed expected by prior physics knowledge. Given that this result emerges from the data-driven RF classifier and the use of counterfactuals without hand-crafted embedding of prior physics intuition, it motivates the usage of counterfactuals to simulate real-time control by generating the β N levels that would have kept the RWM stable for a set of unstable discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Direct-drive laser fusion: status, plans and future

Laser-direct drive (LDD), along with laser indirect (X-ray) drive (LID) and magnetic drive with pulsed power, is one of the three viable inertial confinement fusion approaches to achieving fusion ignition and gain in the laboratory. The LDD programme is primarily being executed at both the Omega Laser Facility at the Laboratory for Laser Energetics and at the National Ignition Facility (NIF) at Lawrence Livermore National Laboratory. LDD research at Omega includes cryogenic implosions, fundamental physics including material properties, hydrodynamics and laser–plasma interaction physics. LDD research on the NIF is focused on energy coupling and laser–plasma interactions physics at ignition-scale plasmas. Limited implosions on the NIF in the ‘polar-drive’ configuration, where the irradiation geometry is configured for LID, are also a feature of LDD research. The ability to conduct research over a large range of energy, power and scale size using both Omega and the NIF is a major positive aspect of LDD research that reduces the risk in scaling from OMEGA to megajoule-class lasers. Furthermore, the paper will summarize the present status of LDD research and plans for the future with the goal of ultimately achieving a burning plasma in the laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Comparison of the static structure factor at long wavelengths for a dusty plasma liquid and other liquids

Especially small values of the static structure factor S(k) at long wavelengths, i.e., small k, were obtained in an analysis of experimental data, for a two-dimensional dusty plasma in its liquid state. For comparison, an analysis of S(k) data was carried out for many previously published experiments with other liquids. Here, the latter analysis indicates that the magnitude of S(k) at small k is typically in a range 0.02–0.13. In contrast, the corresponding value for a dusty plasma liquid was found to be as small as 0.0139. Another basic finding for the dusty plasma liquid is that S(k) at small k generally increases with temperature, with its lowest value, noted above, occurring near the melting point. Simulations were carried out for the dusty plasma liquid, and their results are generally consistent with the experiment. Since a dusty plasma has a soft interparticle interaction, our findings support earlier theoretical suggestions that a useful design strategy for creating materials having exceptionally low values of S(0), so-called hyperuniform materials, is the use of a condensed material composed of particles that interact softly at their periphery.

2-dimensional systems↗
DOE OSTI2023

Finite-size effects in the static structure factor S ( k ) and S ( 0 ) for a two-dimensional Yukawa liquid

Finite-size effects in the static structure factor S(k) are analyzed for an amorphous substance. As the number of particles is reduced, S(0) increases greatly, up to an order of magnitude. Meanwhile, there is a decrease in the height of the first peak S peak. These finite-size effects are modeled accurately by the Binder formula for S(0) and our empirical formula for S peak. Procedures are suggested to correct for finite-size effects in S(k) data and in the hyperuniformity index H ≡ S(0)/S peak. These principles generally apply to S(k) obtained from particle positions in noncrystalline substances. As a result, the amorphous substance we simulate is a two-dimensional liquid, with a soft Yukawa interaction modeling a dusty plasma experiment.

2-dimensional systems↗