Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic plasmas”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗
NASA NTRS2017

Wave-Particle Energy Exchange Directly Observed in a Kinetic Alfven-Branch Wave

Alfven waves are fundamental plasma wave modes that permeate the universe. At small kinetic scales they provide a critical mechanism for the transfer of energy between electromagnetic fields and charged particles. These waves are important not only in planetary magnetospheres, heliospheres, and astrophysical systems, but also in laboratory plasma experiments and fusion reactors. Through measurement of charged particles and electromagnetic fields with NASAs Magnetospheric Multiscale (MMS) mission, we utilize Earths magnetosphere as a plasma physics laboratory. Here we confirm the conservative energy exchange between the electromagnetic field fluctuations and the charged particles that comprise an undamped kinetic Alfven wave. Electrons confined between adjacent wave peaks may have contributed to saturation of damping effects via non-linear particle trapping. The investigation of these detailed wave dynamics has been unexplored territory in experimental plasma physics and is only recently enabled by high-resolution MMS observations.

magnetohydrodynamic↗
DOE OSTI2023

Dynamic Structure of Magnetized High Energy Density Plasmas (Final Report)

This report summarizes results of project DE-SC0022202 ``Dynamic Structure of Magnetized High Energy Density Plasmas,'' which was originally funded for the period 08/01/2021 - 7/31/2022, but was no-cost extended to 7/31/2023. The primary objectives of the work, as stated in the original proposal, were ``develop a theory to describe a novel state of magnetized high energy density plasmas, benchmark the theory using first-principles molecular dynamics simulations, and apply the results to develop a means to interpret x-ray scattering measurements in future experiments.'' These were the goals of the original project proposal, which was for a 3 year project. The project was ultimately funded for 1 year, and so the work completed focused on the goals described for only the first year of the proposed work. This included further developing a kinetic theory for strongly magnetized plasmas, applying it to compute fluid-scale (MHD) transport coefficients, and to test the predictions using molecular dynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electron heating in kinetic-Alfvén-wave turbulence

We report analytical and numerical investigations of subion-scale turbulence in low-beta plasmas using a rigorous reduced kinetic model. We show that efficient electron heating occurs and is primarily due to Landau damping of kinetic Alfvén waves, as opposed to Ohmic dissipation. This collisionless damping is facilitated by the local weakening of advective nonlinearities and the ensuing unimpeded phase mixing near intermittent current sheets, where free energy concentrates. The linearly damped energy of electromagnetic fluctuations at each scale explains the steepening of their energy spectrum with respect to a fluid model where such damping is excluded (i.e., a model that imposes an isothermal electron closure). In conclusion, the use of a Hermite polynomial representation to express the velocity-space dependence of the electron distribution function enables us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Momentum conservation in current drive and alpha-channeling-mediated rotation drive

Alpha channeling uses waves to extract hot ash from a fusion plasma, transferring energy from the ash to the wave. It has been proposed that this process could create a radial electric field, efficiently driving E × B rotation. However, existing theories ignore the nonresonant particles, which play a critical role in enforcing momentum conservation in quasilinear theory. Because cross field charge transport and momentum conservation are fundamentally linked, this non-consistency throws the rotation drive into question. This paper has two main goals. First, we provide a pedantic and cohesive introduction to the recently developed simple, general, self-consistent quasilinear theory for electrostatic waves that explains the torques which allow for current drive parallel to the magnetic field, and charge extraction across it; a theory that has largely resolved the question of rotation drive by alpha channeling. We show how the theory reveals a fundamental difference between the reaction of nonresonant particles to plane waves that grow in time vs steady-state waves that have a nonuniform spatial structure, allowing rotation drive in the latter case while precluding it in the former, and we review the local and global conservation laws that lead to this result. Second, we provide two new results in support of the theory. First, we provide a novel two-particle Hamiltonian model that rigorously establishes the relationship between charge transport and momentum conservation. Second, we compare the new quasilinear theory to the oscillation-center theories of ponderomotive forces, showing how the latter often obscure the time-dependent nonresonant recoil, but ultimately lead to similar results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Electron-cyclotron-resonance (ECR) plasma acceleration

A research effort directed at analytically and experimentally investigating electron-cyclotron-resonance (ECR) plasma acceleration is outlined. In addition, relevant past research is reviewed. Also, the prospects for application of ECR plasma acceleration to spacecraft propulsion are described. It is shown that previously unexplained losses in converting microwave power to directed kinetic power via ECR plasma acceleration can be understood in terms of diffusion of energized plasma to the physical walls of the accelerator. It is also argued that line radiation losses due to electron-ion and electron-atom inelastic collisions should be less than estimated in past research. Based on this new understanding, the expectation now exists that very efficient ECR plasma accelerators can be designed for application to high-specific-impulse spacecraft propulsion.

Sercel, Joel C.↗
DOE OSTI2022

Drift kinetic theory of the NTM magnetic islands in a finite beta general geometry tokamak plasma

In (Imada et al 2019 Nucl. Fusion 59 046016 and references therein) a new 4D drift kinetic nonlinear theory, valid in the limit of a low beta, small inverse aspect ratio, circular cross section, toroidal geometry, to describe the plasma response to the neoclassical tearing mode (NTM) magnetic perturbation is derived. In (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) this theory is reduced in a low collisionality limit, which allows a dimensionality reduction to a 3D problem to efficiently resolve the collisional dissipation layer in the vicinity of the trapped-passing boundary. (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) adopts an improved model for the magnetic drift frequency, which reduces the threshold magnetic island half-width from $8.73\rho_{b i}$, where $\rho_{b i}$ is the trapped ion banana orbit width, to $1.46 \rho_{b i}$, making it in closer agreement with experimental observations for the large aspect ratio tokamak equilibrium. In the present paper, the theory is extended to a high beta, arbitrary tokamak geometry to capture the plasma shaping effects on the NTM threshold physics with the focus on the non-zero triangularity discharges that are known to have a strong impact on the plasma MHD stability. First, it is found that the higher triangularity plasma is more prone to NTMs which is in agreement with the 2/1 tearing mode onset relative frequency measurements in DIII-D. Second, the NTM threshold dependence on the tokamak inverse aspect ratio obtained in (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) is refined and extended to a finite aspect ratio limit. Third, the NTM threshold dependence on poloidal beta is obtained and benchmarked against the EAST threshold island width measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Processes in relativistic plasmas

The establishment and maintenance of a Boltzmann distribution in particle kinetic energies is investigated for a plasma with theta = KTe/mc-squared much greater than unity, where m is the electron mass. It is shown that thermalization of the electron gas by binary collisions is not sufficiently effective to maintain the equilibrium distribution when other processes that perturb the equilibrium are taken into account. Electron-positron pair production in electron-electron and electron-ion collisions, and perturbations of a Boltzmann distribution by nonthermal processes are evaluated. Thermalization by means of other mechanisms, such as interaction with plasma waves is discussed, and the opacity of a relativistic plasma is computed for Compton scattering, pair production in the fields of electrons and ions, inverse bremsstrahlung, and synchrotron self-absorption.

Gould, R. J.↗
DOE OSTI2020

Expanded capabilities of the CarMa code in modeling resistive wall mode dynamics with 3-D conductors

In this work, an improved version of the CarMa code is presented, called CarMa-D, for the analysis of resistive wall modes (RWMs) in fusion devices, simultaneously considering the effects of volumetric three-dimensional conducting structures, plasma dynamics, toroidal rotation or drift-kinetic damping. Unlike static CarMa, CarMa-D does not rely on the simplifying assumptions such as neglecting the plasma mass, toroidal rotation and kinetic damping. The new coupling strategy is based on matrix-based Padé rational functions approximation of plasma a response. The arising mathematical model is formally equivalent to the original CarMa model, but with a higher number of degrees of freedom to model the dynamics of the plasma. CarMa-D overcomes the main limitations of the original CarMa, in particular: (i) the massless assumption for the plasma is removed, allowing modeling of global modes growing both on ideal kink time scales and in the typical RWM growth rate regime, with a suitable treatment of the model; (ii) the effects of toroidal plasma flow and drift kinetic damping can be included into the new model, providing a powerful tool to study macroscopic phenomena where both plasma dynamics and 3-D conducting structures play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Factors affecting ion kinetic temperature, number density, and containment time in the NASA Lewis bumpy-torus plasma

The degree of toroidal symmetry of the plasma, the number of midplane electrode rings, the configuration of electrode rings, and the location of the diagnostic instruments with respect to the electrode rings used to generate the plasma are discussed. Impurities were deliberately introduced into the plasma, and the effects of the impurity fraction on ion kinetic temperature and electron number density were observed. It is concluded that, if necessary precautions are taken, the plasma communicates extremely well along the magnetic field lines and displays a high degree of symmetry from sector to sector for a wide range of electrode ring configurations and operating conditions. Finally, some characteristic data taken under nonoptimized conditions are presented, which include the highest electron number density and the longest particle containment time (1.9 msec) observed. Also, evidence from a paired comparison test is presented which shows that the electric field acting along the minor radius of the toroidal plasma improves the plasma density and the calculated containment time more than an order of magnitude if the electric field points inward, relative to the values observed when it points (and pushes ions) radially outward.

Roth, J. R.↗
DOE OSTI2020

Generalized slab universal instability and its appearance in pair plasma

In this work, a generalized linear dispersion relation of electromagnetic slab universal modes is derived, taking into account arbitrary ion charge state, electron finite Larmor radius (FLR) effects, and Debye shielding effects. As a consequence, it provides more accurate predictions and is applicable to a wider range of plasmas. We find that electron FLR effects have a weakly stabilizing effect on the slab universal instability in electron–ion plasma, while Debye shielding strongly stabilizes this instability when λ D approaches ρ i (λ D is the Debye length and ρi is the ion gyroradius). In particular, we examine the stability criterion for this instability in electron–positron pair plasmas and find that the instability persists in this simplest plasma system as long as the pair plasma number density exceeds the critical value n c = B 2 /(8πm e c 2 ).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wavelength modulation laser-induced fluorescence for plasma characterization

Laser-Induced Fluorescence (LIF) spectroscopy is an essential tool for probing ion and atom velocity distribution functions (VDFs) in complex plasmas. VDFs carry information about the kinetic properties of species that is critical for plasma characterization. Accurate interpretation of these functions is challenging due to factors such as multicomponent distributions, broadening effects, and background emissions. Our research investigates the use of Wavelength Modulation (WM) LIF to enhance the sensitivity of VDF measurements. Unlike standard Amplitude Modulation (AM) methods, WM–LIF measures the derivative of the LIF signal. This approach makes variations in VDF shape more pronounced. VDF measurements with WM–LIF were investigated with both numerical modeling and experimental measurements. The developed model enables the generation of both WM and AM signals, facilitating comparative analysis of fitting outcomes. Experiments were conducted in a weakly collisional argon plasma with magnetized electrons and non-magnetized ions. Measurements of the argon ion VDFs employed a narrow-band tunable diode laser, which scanned the 4p 4 D 7/2 –3d 4 F 9/2 transition centered at 664.553 nm in vacuum. A lock-in amplifier detected the second harmonic WM signal, which was generated by modulating the laser wavelength with an externally controlled piezo-driven mirror of the diode laser. Finally, our findings indicate that the WM–LIF signal is more sensitive to fitting parameters, allowing for better identification of VDF parameters such as the number of distribution components, their temperatures, and velocities. In addition, WM–LIF can serve as an independent method to verify AM measurements and is particularly beneficial in environments with substantial light noise or background emissions, such as those involving thermionic cathodes and reflective surfaces.

47 OTHER INSTRUMENTATION↗
NASA NTRS2014

The Link Between Shocks, Turbulence, and Magnetic Reconnection in Collisionless Plasmas

Global hybrid (electron fluid, kinetic ions) and fully kinetic simulations of the magnetosphere have been used to show surprising interconnection between shocks, turbulence and magnetic reconnection. In particular collisionless shocks with their reflected ions that can get upstream before retransmission can generate previously unforeseen phenomena in the post shocked flows: (i) formation of reconnecting current sheets and magnetic islands with sizes up to tens of ion inertial length. (ii) Generation of large scale low frequency electromagnetic waves that are compressed and amplified as they cross the shock. These 'wavefronts' maintain their integrity for tens of ion cyclotron times but eventually disrupt and dissipate their energy. (iii) Rippling of the shock front, which can in turn lead to formation of fast collimated jets extending to hundreds of ion inertial lengths downstream of the shock. The jets, which have high dynamical pressure, 'stir' the downstream region, creating large scale disturbances such as vortices, sunward flows, and can trigger flux ropes along the magnetopause. This phenomenology closes the loop between shocks, turbulence and magnetic reconnection in ways previously unrealized. These interconnections appear generic for the collisionless plasmas typical of space, and are expected even at planar shocks, although they will also occur at curved shocks as occur at planets or around ejecta.

MAGNETOSPHERE↗
DOE OSTI2022

Verification of neoclassical toroidal viscosity induced by energetic particles

The thermal particles contributed neoclassical toroidal viscosity (NTV) have been successfully developed and explored by many impressive works such as the study by Shaing et al. [Phys. Plasmas 10, 1443 (2003)] and Zhu et al. [Phys. Rev. Lett. 96, 225002 (2006)]. In this work, the scope of the NTV study is extended to explore the contribution of energetic particles (EPs) through both theory and experiments. In theory, the existence of the NTV torque due to the precessional drift resonance of trapped EPs is identified based on the equivalence between the NTV torque and the perturbed drift kinetic energy [J. Park, Phys. Plasmas 18, 110702 (2011)]. Toroidal modeling with the Magneto Resistive Spectrum - drift Kinetic code [Y. Liu, Phys. Plasmas 15, 112503 (2008)], based on this equivalence, indicates that trapped EPs can contribute a significant amount of the NTV torque. Meanwhile, this work also focuses on developing the dedicated DIII-D experiments in the presence of the n = 2 external magnetic perturbation to verify the EP induced NTV (EP-NTV) by measuring the change of the NTV torque while varying the angle and the voltage of the neutral beam injection. However, the developed experiments have been unable to create conditions necessary to clearly demonstrate the presence of EP-NTV. The main challenge is separating the resonant and non-resonant momentum transport responses in the plasma. Finally, the experience, gained from this study, can help the further exploration of EP-NTV in the future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MARS-F/K modeling of plasma response and fast ion losses due to RMP in KSTAR

Abstract The toroidal single-fluid magnetohydrodynamic (MHD) code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and the MHD-kinetic hybrid code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) are utilized to study the plasma response to the n = 1 ( n is the toroidal mode number) resonant magnetic perturbation (RMP), applied to suppress the type-I edge localized mode (ELM) in a KSTAR discharge. Both the resistive-rotating and ideal-static plasma models identify strong screening of the resonant radial field harmonics of the applied RMP due to the plasma response, and predict a strong edge-peeling response of the plasma which is consistent with the optimal ELM control coil current configuration adopted in experiment. The RMP-induced radial displacement of the plasma, computed by the resistive-rotating plasma model, agrees reasonably well with that reconstructed from the measured data in the plasma core. Taking into account the drift kinetic response of fast ions, MARS-K hybrid modeling also finds quantitative agreement of the plasma core fluid pressure perturbation with experiment. Based on the MARS-F computed plasma response, a guiding-center orbit-tracing simulation finds about 0.3% of fast ion losses due the n = 1 RMP in the KSTAR ELM control experiment considered. Most losses are associated with counter-current fast ions located near the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Gyrokinetic profile prediction and validation of a negative triangularity plasma in ASDEX Upgrade

In this work, gyrokinetic simulations are performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73–93) for a negative triangularity H-mode plasma in ASDEX Upgrade, and compared with experimental measurements. The PORTALS framework (Rodriguez-Fernandez et al 2024 Nucl. Fusion 64 076034) is used to accelerate the prediction of kinetic profiles for this plasma, using surrogate modeling and Bayesian optimization. Ion heat flux, electron heat flux, and electron particle flux are simultaneously matched across the simulated radial regime of the plasma (normalized radius $r/a = 0.35-0.90$), and the resulting ion temperature, electron temperature, and electron density profiles match well with the experimental profile data within this radial range. A synthetic Correlation Electron Cyclotron Emission diagnostic is applied to find well-matched electron temperature fluctuation properties between simulation and experiment. The flux-matched profiles provide a basis for investigation of the turbulence nature across the plasma radius, revealing the dominance of Trapped Electron Mode turbulence at $r/a = 0.35$, the dominance of Ion Temperature Gradient turbulence at $r/a = 0.55$, 0.75, and 0.83, and an instability boundary at $r/a = 0.90$.

gyrokinetic simulation↗
NASA NTRS1979

Theory of the large-amplitude plane magnetoacoustic wave propagating transverse to the magnetic field in a hot collisionless plasma

An exact solution of the kinetic and electromagnetic equations for a large-amplitude plane magnetoacoustic wave propagating transverse to the magnetic field in a hot collisionless plasma is presented. The solution gives simple relations among the magnetic-field strength, density, stress tensor, and plasma velocity, all of which are measurable in the interplanetary plasma. These relations are independent of the electron and ion velocity distributions, subject to certain restrictions on 'high-velocity tails.' The magnetic field of the wave is linearly polarized. The wave steepens to form a shock much as the analogous waves of MHD theory do.

Barnes, A.↗