Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutrals in plasma”

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

3D Exact Analytical Solutions of Two-fluid Plasma, Magnetohydrodynamics, and Neutral Fluid Equations for the Creation of Ordered Structures as well as Jet-like Flows

The 3D exact analytical solutions of ideal two-fluid plasma, single-fluid plasma, and neutral fluid equations have been found using physically justifiable assumptions. Surprisingly these solutions satisfy all nonlinearities in the systems. It is pointed out that these solutions explain the fundamental mechanism behind the creation of a vast variety of ordered structures in plasmas and fluids. In the limiting case of 2D dependence of fields, the theoretical model for plasma is applied to explain the formation of spicules in the solar chromosphere. It is pointed out that the main contribution of electron (ion) baroclinic vectors is to produce vorticity in the plasma, and that magnetic field generation is coupled with the flow of both electrons and ions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Time-dependent modeling of coupled plasma-wall dynamics

Effects of dynamic wall outgassing on divertor plasma during an ELM-like heat pulse are modelled using newly coupled edge plasma code UEDGE and all reaction-diffusion code FACE. Different divertor regimes are simulated depending on the degree of the plasma detachment. It is shown that hydrogen outgassing from the divertor target plate induced by ELM pulse in the semi-detached regime can trigger deep plasma detachment. The induced detachment does not spontaneously reverse after the pulse due to ineffective implantation of low energy plasma hydrogen in the target material. Furthermore, we demonstrate that plasma neutrals play essential role in heat transport in a detached divertor. The modelled hydrogen dynamics in the wall material demonstrates that the heat pulse leads to enhanced hydrogen transport toward the wall bulk.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Investigation of divertor detachment induced through neon seeding and density ramp on HL-3

A new self-consistent 1D scrape-off layer model has been recently developed in BOUT++ framework, named SD1D, which includes equations for various particle species (e.g. main plasma, neutrals and impurities) and couples open databases like ADAS and AMJUEL. It is able to quickly and effectively simulate divertor detachment experiments. In this work, a typical detachment experiment (shot #6270) on HL-3 with neon seeding is simulated using the SD1D code. It is found that the target electron temperature and the target ion saturation current in the simulations are consistent with experimental results measured by Langmuir probes on the target plate. The variation of D α radiation intensity in the divertor is qualitatively similar to the measured D α signal. Following the experimental validations, different upstream densities are set in the simulations to study the impurity distribution under different plasma density conditions. It is found that increasing upstream density can be helpful for the control of the neon radiation front (closer to the target). In this work we also compare two detachment regimes in simulations. Based on the same initial experimental parameters (shot #6270) on HL-3, a scan of upstream density and a scan of neon seeding rate are carried out respectively. It is found that the role of atomic and molecular processes is different in the two detachment regimes. The current density roll-over is ascribed to a drop in the divertor ion source, and the variation of D α radiation intensity via different excitation channels is associated with the relevant collisional reaction sources.

BOUT++↗
DOE OSTI2026

Post-puff SOL broadening on MAST-U under high-recycling conditions: evidence consistent with cross-field transport changes

Transient broadening of the scrape-off layer (SOL) density profile can modify main-chamber first-wall particle fluxes and divertor loading, yet its control parameters remain debated between divertor-regime transitions, neutral dynamics and changes in cross-field transport. We investigate fueling-driven SOL density-profile evolution and post-fueling relaxation on MAST-Upgrade (MAST-U) in ohmic L-mode, using two otherwise similar double-null Conventional Divertor discharges (I p = 450 kA, B T = 0.33 T) with identical 50 ms low-field side gas puffs; the only intentional difference is the puff start time. Upstream Thomson scattering shows that both discharges develop a transient far-SOL density profile modification, expressed as an increased SOL-width metric λ n e and a far-SOL enhancement consistent with a shoulder-like signature. In the earlier-puff case, the SOL-width metric remains elevated after puff termination and the post-puff decay time scales are systematically longer across the analysed radii. Outer-target Langmuir probes indicate high-recycling conditions during the analysed window (few-eV T e with radially peaked j sat,∥ and q ∥ profiles, without signatures of deep detachment). The outer-divertor collisionality proxy Λ div is elevated in both cases and does not discriminate between the different post-puff persistence. The target-integrated ion flux proxy ∫ J sat dA evolves nearly identically in both discharges when aligned to the puff start. Taken together, these observations suggest that the late/post-puff upstream SOL evolution is not set by parallel exhaust to the outer targetalone and is more consistent with upstream cross-field redistribution, with a possible role for plasma–neutral coupling and fueling geometry.

MAST-U↗
DOE OSTI2024

Density profiles in stellarators: an overview of particle transport, fuelling and profile shaping studies at TJ-II

We provide an overview of activities carried out at the TJ-II stellarator for improving our understanding of- and developing plasma physics models for particle density profiles in stellarators. Namely, we report on recent progress in turbulent particle transport simulation, validation of pellet deposition models, density profile shaping for performance control and new experimental techniques for edge turbulence and plasma-neutral interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A gyrokinetic simulation model for 2D equilibrium potential in the scrape-off layer of a field-reversed configuration

The equilibrium potential structure in the scrape-off layer (SOL) of the field-reversed configuration (FRC) can be affected by the penetration of edge biasing applied at the divertor ends. The primary focus of the paper is to establish a formulation that accurately captures both parallel and radial variations of the two-dimensional (2D) potential in SOL. The formulation mainly describes a quasi-neutral plasma with a logical sheath boundary. A full-f gyrokinetic ion model and a massless electron model are implemented in the GTC-X code to solve for the self-consistent equilibrium potential, given fixed radial potential profiles at the boundaries. The first essential point of this 2D model lies in its ability to couple radial and parallel dynamics stemming from resistive currents and drag force on ions. The model successfully recovers the fluid force balance and continuity equations. These collisional effects on 2D potential mainly appear through the density profile changes, modifying the potential through electron pressure gradient. This means an accurate prescription of electron density and temperature profiles is important in predicting the potential structure in the FRC SOL. The Debye sheath potential and the potential profiles applied at the boundaries can be additional factors contributing to the 2D variations in SOL. This comprehensive full-f scheme holds promise for future investigations into turbulent transport in the presence of the self-consistent 2D potential together with the non-Maxwellian distributions and open boundary conditions in the FRC SOL.

Physics↗
DOE OSTI2024

Impact of gas injection location and divertor surface material on ITER fusion power operation phase divertor performance assessed with SOLPS-ITER *

Abstract The ITER divertor design and performance assessment, primarily based on the SOLPS-4.3 burning plasma database (Pitts R. et al 2019 Nucl. Mater. Energy 20 100696), assumes the use of beryllium (Be) as the divertor surface material and the injection of gas from the main chamber top. However, the current ITER baseline favors gas injection from the more toroidally symmetric sub-divertor region. This paper evaluates the implications of these assumptions for divertor performance in the ITER fusion power operation phase. The impact of the divertor surface material and the gas injection location on the main ions mirrors the hydrogen only low power phase scenario shown in Park J.-S. et al (2020 Nucl. Fusion 61 016021). However, during burning plasma operation, extrinsic impurity seeding will be required. In the case of neon (Ne), studied here, impurity retention is influenced by both the divertor surface material and the fueling location. Neon leakage increases due to more energetic reflection from tungsten than beryllium, but equivalent divertor performance can be achieved by adjusting the neon seeding rate. While the impurity seeding location does not affect the distributions of impurity or radiation, the fueling location does. Top fueling provides local ionization sources mainly in the mid-SOL under detached conditions, enhancing divergences of the flux there (source-driven flow), bringing stagnation points close to the fueling location, and equilibrating flows towards both targets. In contrast, the global flow pattern (in the absence of fluid drifts) in the case of sub-divertor fueling is biased towards the inner target. Impurity flows, driven by force balance, largely mirror those of the main ion flow, including the stagnation point. The case with top fueling enhances Ne retention and corresponding radiation in the outer divertor, effectively reducing the total and peak target heat fluxes by 20%–40%, compared to the case with divertor fueling. Meanwhile, the case with outer target fueling also achieves similar reductions by enhancing plasma-neutral interactions. These results suggest the possibility that the selection of the fueling location and throughput can be used as an actuator to control impurity divertor retention and divertor radiation asymmetry.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Full-orbit simulations of fast-ion charge-exchange losses induced by neutral particles outside the plasma last closed flux surface

The interaction between neutrals outside tokamak plasmas and fast ions that are generated by neutral beam injection and pass outside the plasma is studied. It is shown that at beam energies (up to ~100~keV) fast-ion losses are significant, up to 20% in the studied case, because of charge exchange reactions between the neutrals and the fast ions. About half of the fast ions that neutralize outside the plasma get lost to the wall while the other half re-enter the plasma and ionize again, often on better confined orbits which is visible as an enhanced beam power deposition in the plasma core (rho ~< 0.7). As a result, the the plasma performance as measured from the beam-thermal neutrons is hardly changed despite the significant fast-ion losses toward the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Diocotron modes in pure electron plasmas in the APEX levitating dipole trap

Abstract Observations of persistent toroidal modes in dipole-trapped pure electron plasmas are presented. Non-neutral electron plasmas were confined by the poloidal magnetic field of the APEX (A Positron Electron eXperiment) levitating dipole trap. The negative bias of the electron emitter was gated to control the filling time. Large amplitude, narrow-band oscillations were consistently observed in wall probe measurements after the fill. The oscillatory image-charge currents induced by the trapped electron plasma typically persisted for many seconds with fundamental frequencies, f 0 = 50 to 500 kHz, that scale with the E → × B → drift velocity. We attribute this feature to the dipole-equivalent of the diocotron mode that is commonly observed in non-neutral plasmas in linear traps.

Deller, A. (ORCID:0000000234301501)↗
DOE OSTI2026

Machine Learning Approach to Modeling of Neutral Particles Transport in Plasma

A propagator‐based approach is investigated for Monte‐Carlo (MC) modeling of neutral particle transport in fusion boundary plasmas. The propagator is based on a Green's function for the neutral kinetic equation, which depends on the plasma profiles. A neural network (NN)‐based model for the propagator provides a fast and accurate solution for the neutral distribution function in plasma. Preliminary results from a small 1D test problem look encouraging. The proposed approach, a propagator‐based NN model for neutral transport in plasma, has potential for generalization to higher dimensions and efficient coupling with plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Collisional-energy-cascade model for nonthermal velocity distributions of neutral atoms in plasmas

Nonthermal velocity distributions with much greater tails than a Maxwellian have been observed for radical atoms in plasmas for a long time. Historically, such velocity distributions have been modeled by a two-temperature Maxwell distribution. In this paper, I propose a model based on collisional energy cascade, which has been studied in the field of granular materials. In the collisional energy cascade, a particle ensemble undergoes energy input at the high-energy region, entropy production by elastic collisions among particles, and energy dissipation. For radical atoms, energy input may be caused by the Franck-Condon energy of molecular dissociation or charge-exchange collision with hot ions, and the input energy is eventually dissipated by collisions with the walls. Here, I show that the steady-state velocity distribution in the collisional energy cascade is approximated by the generalized Mittag-Leffler distribution, which is a one-parameter extension of the Maxwell distribution. This parameter indicates the degree of the nonthermality and is related to the relative importance of energy dissipation over entropy production. This model is compared with a direct molecular dynamics simulation for simplified gaseous systems with energy input, as well as some experimentally observed velocity distributions of light radicals in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Spectroscopic camera analysis of the roles of molecularly assisted reaction chains during detachment in JET L-mode plasmas

The roles of the molecularly assisted ionization (MAI), recombination (MAR) and dissociation (MAD) reaction chains with respect to the purely atomic ionization and recombination processes were studied experimentally during detachment in low-confinement mode (L-mode) plasmas in JET with the help of experimentally inferred divertor plasma and neutral conditions, extracted previously from filtered camera observations of deuterium Balmer emission, and the reaction coefficients provided by the ADAS, AMJUEL and H2VIBR atomic and molecular databases. The direct contribution of MAI and MAR in the outer divertor particle balance was found to be inferior to the electron–atom ionization (EAI) and electron–ion recombination (EIR). Near the outer strike point, a strong atom source due to the D$^+_2$ -driven MAD was, however, observed to correlate with the onset of detachment at outer strike point temperatures of T e,osp = 0.9 –2.0 eV via increased plasma-neutral interactions before the increasing dominance of EIR at T e,osp < 0.9 eV, followed by increasing degree of detachment. The analysis was supported by predictions from EDGE2D-EIRENE simulations which were in qualitative agreement with the experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neutral Atomic Beam Techniques for Diagnosis of Edge and Scrape-Off Layer Plasmas (2017-2019) [and subsequently titled] Neutral Atomic Beam Techniques for Diagnosis of Interior and Scrape-Off Layer Plasmas (2019-2022) (Final Report)

This is the Final Technical Report for the DOE Measurement Innovation grant DE-SC0017998 titled Neutral Atomic Beam Techniques for Diagnosis of Edge and Scrape-Off Layer Plasmas (2017-2019) and subsequently titled Neutral Atomic Beam Techniques for Diagnosis of Interior and Scrape-Off Layer Plasmas (2019-2022). The grant was funded by the Office of Fusion Energy Sciences (FES) for the period 1 September 2017 – 31 May 2022. Work performed through the grant has resulted in diagnostic innovations and new measurement capabilities that could improve understanding of magnetic equilibrium, edge localized modes in tokamaks, and bootstrap and Pfirsch-Schlüter currents in stellarators. We have performed simulations of a novel non-hydrogen beam-based diagnostic that would determine a local value (in the plasma) of the poloidal flux function by measuring the toroidal angle of the beam particles' velocity at a detector external to the plasma. The approach is similar to that of a Heavy Ion Beam Probe (HIBP), but would be able to operate with a smaller, more economical detection systems and, in doing so, may enable substantial extension of scenarios in which beam-based diagnostics are deployed. Key results include: Identification of the neutral and singly charged atoms that are most appropriate for use in the beam based diagnostic and the ionization cross-sections that are important for studying the plasma edge of a large device; Development of a promising geometry for operation in which the beam is injected horizontally near the torus midplane; Implementation of simulations using banana orbit trajectories of secondary beam ions to access plasma regions that are commonly in-accessible to an HIBP diagnostic, including sample volumes near the X-point of a diverted tokamak; Assessment of the ability of the diagnostic concept (originally developed under the assumption of, and reliant on, symmetry) to work in stellarators by performing simulations with both quasi-symmetric fields and non-quasi-symmetric fields; Measurement of noise levels with a prototype detector on a stellarator and exploration of techniques to minimize effects of noise on a diagnostic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

1D model of tokamak scrape-off layer and divertor plasmas: Dimensionless parameters and the results of numerical simulations

Solutions of a 1D model of the scrape-off layer and a divertor plasma in a tokamak are analyzed on the basis of dimensionless parameters governing the physics of both plasma and neutral transport and atomic physics processes. Numerical simulations based on fluid plasma and Monte Carlo neutral modeling confirm earlier results obtained with simplified analytic considerations and numerical simulations performed with fluid neutral description and demonstrate bifurcation of plasma parameters. Further, these simulations also confirm that the major physical reasons for the rollover of the plasma flux to the target, which is the key feature of the phenomenon of divertor plasma detachment, are the impurity radiation loss and plasma recombination. It is suggested that simulations of the edge plasma based on dimensionless parameters provide a powerful tool for checking the consistency of the physics incorporated in existing complex 2D codes used for edge plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗