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

Validation of the imaging neutral particle analyzer in nearly MHD quiescent plasmas using injected beam ions on DIII-D

The imaging neutral particle analyzer (INPA) is a scintillator-based diagnostic that provides energy and radially resolved measurements of confined fast ions. To verify operation, understanding, and modeling of the diagnostic performance, distributions of deuterium beam ions in low density, nearly MHD-quiescent plasmas are measured on the DIII-D tokamak using the INPA. Here, the distribution is altered by changing the pitch-angle scattering rate with electron cyclotron heating (ECH). Experimental images are mostly reproduced with simulation of core fast ions charge exchanging with injected beam and edge cold neutrals in classical cases. Image comparisons between the experiment and simulation show agreement with less than 25% discrepancies near the beam injection energy of 50 keV. Other strong passive signals, observed in a localized region on the image, are suggested to come from charge exchange between cold neutrals and fast ions around the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics basis for the Wisconsin HTS Axisymmetric Mirror (WHAM)

The Wisconsin high-temperature superconductor axisymmetric mirror experiment (WHAM) will be a high-field platform for prototyping technologies, validating interchange stabilization techniques and benchmarking numerical code performance, enabling the next step up to reactor parameters. A detailed overview of the experimental apparatus and its various subsystems is presented. WHAM will use electron cyclotron heating to ionize and build a dense target plasma for neutral beam injection of fast ions, stabilized by edge-biased sheared flow. At 25 keV injection energies, charge exchange dominates over impact ionization and limits the effectiveness of neutral beam injection fuelling. This paper outlines an iterative technique for self-consistently predicting the neutral beam driven anisotropic ion distribution and its role in the finite beta equilibrium. Beginning with recent work by Egedal et al. ( Nucl. Fusion , vol. 62, no. 12, 2022, p. 126053) on the WHAM geometry, we detail how the FIDASIM code is used to model the charge exchange sources and sinks in the distribution function, and both are combined with an anisotropic magnetohydrodynamic equilibrium solver method to self-consistently reach an equilibrium. We compare this with recent results using the CQL3D code adapted for the mirror geometry, which includes the high-harmonic fast wave heating of fast ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗
DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Adiabatic behavior of an elliptical vortex in a time-dependent external strain flow

The adiabatic behavior of two-dimensional (2D) fluid vortices subject to slowly varying, time-dependent external strain flows is studied theoretically, and experimentally using magnetized non-neutral electron plasmas. Here the E × B drift dynamics of the electrons perpendicular to the magnetic field are analogous to the motion of fluid vorticity under the 2D Euler equations describing ideal fluids. A low-dimensional elliptical vortex patch model is used along with a WKB approximation to derive a formula for the breaking of an adiabatic invariant due to the changing external flow. The invariant is interpreted as the amplitude of a perturbative oscillation about a stable fixed point corresponding to a steady elliptical vortex. Smooth, hyperbolic tangent and piecewise linear ramp functions are considered for the external strain time dependence. Standard exponential breaking is observed in the former case, whereas the latter exhibits a power-law breaking curve with periodic modulation. It is found that a driving term in the equations of motion contributes the majority of the breaking whereas the frequency variation plays a weaker, but significant role. The experimental data agree closely with the theoretical model. The most significant deviation is due to inviscid damping behavior associated with the smooth edges of the experimental vortices, which tends to reduce the amplitude of the oscillation. Finally, the results are compared and contrasted with other related experimental, numerical, and theoretical work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hermes-3

Hermes-3 is a plasma simulation model built on BOUT++, developed as part of research studies of the edge and divertor regions of magnetically confined fusion plasmas such as tokamaks. The main purpose of this model is multi-species simulation of plasmas containing a mixture of deuterium, tritium, helium and other plasma and neutral gas species.

Dudson, BenjaminD↗
DOE OSTI2023

Star-exoplanet interactions: A growing interdisciplinary field in heliophysics

Traditionally, heliophysics is characterized as the study of the near-Earth space environment, where plasmas and neutral gases originating from the Earth, the Sun, and other solar system bodies interact in ways that are detectable only through in-situ or close-range (usually within ~10 AU) remote sensing. As a result, heliophysics has data from the space environment around a handful of solar system objects, in particular the Sun and Earth. Comparatively, astrophysics has data from an extensive array of objects, but is more limited in temporal, spatial, and wavelength information from any individual object. Thus, our understanding of planetary space environments as a complex, multi-dimensional network of specific interacting systems may in the past have seemed to have little to do with the highly diverse space environments detected through astrophysical methods. Recent technological advances have begun to bridge this divide. Exoplanetary studies are opening up avenues to study planetary environments beyond our solar system, with missions like Kepler, TESS, and JWST, along with increasing capabilities of ground-based observations. At the same time, heliophysics studies are pushing beyond the boundaries of our heliosphere with Voyager, IBEX, and the future IMAP mission. The interdisciplinary field of star-exoplanet interactions is a critical, growing area of study that enriches heliophysics. A multidisciplinary approach to heliophysics enables us to better understand universal processes that operate in diverse environments, as well as the evolution of our solar system and extreme space weather. The expertise, data, theory, and modeling tools developed by heliophysicists are crucial in understanding the space environments of exoplanets, their host stars, and their potential habitability. The mutual benefit that heliophysics and exoplanetary studies offer each other depends on strong, continuing solar system-focused and Earth-focused heliophysics studies. The heliophysics discipline requires new targeted funding to support inter-divisional opportunities, including small multi-disciplinary research projects, large collaborative research teams, and observations targeting the heliophysics of planetary and exoplanet systems. Here we discuss areas of heliophysics-relevant exoplanetary research, observational opportunities and challenges, and ways to promote the inclusion of heliophysics within the wider exoplanetary community.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A

Divertor detachment is a promising method to solve the power exhaust problem in tokamak devices or even in future magnetic fusion reactors. In this work, a detachment experiment (HL-2A shot #38008) with mixed gas seeding (60% nitrogen and 40% deuterium) is simulated using the SD1D module in BOUT++. In the process from attachment to detachment, the target electron temperature and the target ion saturation current in simulations are found to be consistent with the experimental results measured by Langmuir probes on the target plate. In order to understand the underlying detachment mechanism on HL-2A, this work analyses the role of different particle species in the cases with different seeding rates. It shows that the plasma density varies little and the density of neutrals (D and D 2 ) slightly decreases with the increase of nitrogen seeding rate, such that plasma–neutral interactions cannot effectively reduce plasma energy and plasma momentum in the divertor. The case with a high seeding rate predicts that increasing seeding rate cannot reach a target temperature lower than 2.5 eV, which is the required temperature for strong plasma–molecule interactions. Thus, the plasma–molecule interactions may not be important in the divertor during nitrogen seeding. Here this work also studies the parallel forces including (1) force due to the parallel electric field, (2) friction forces, (3) ion- and electron-thermal forces, and (4) collisional reactions (e.g. charge exchange recombination and ionisation). It is found that the friction force between nitrogen ions and other particle species (primarily D + ) is the dominant force towards the target, while ion- and electron-thermal forces are the dominant force pushing nitrogen ions back to upstream. Parallel forces determine the parallel distribution of nitrogen impurities, and therefore decide the region of nitrogen radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of particle transport, neutrals and radiation in magnetically-confined plasmas with Aurora

In this work, we present Aurora, an open-source package for particle transport, neutrals and radiation modeling in magnetic confinement fusion plasmas. Aurora's modern multi-language interface enables simulations of 1.5D impurity transport within high-performance computing frameworks, particularly for the inference of particle transport coefficients. A user-friendly Python library allows simple interaction with atomic rates from the Atomic Data and Atomic Structure database as well as other sources. This enables a range of radiation predictions, both for power balance and spectroscopic analysis. We discuss here the superstaging approximation for complex ions, as a way to group charge states and reduce computational cost, demonstrating its wide applicability within the Aurora forward model and beyond. Aurora also facilitates neutral particle analysis, both from experimental spectroscopic data and other simulation codes. Leveraging Aurora's capabilities to interface SOLPS-ITER results, we demonstrate that charge exchange is unlikely to affect the total radiated power from the ITER core during high performance operation. Finally, we describe the ImpRad module in the one modeling framework for integrated task framework, developed to enable experimental analysis and transport inferences on multiple devices using Aurora.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Disorder-induced heating as a mechanism for fast neutral gas heating in atmospheric pressure plasmas

Abstract Recent findings suggest that ions are strongly correlated in atmospheric pressure plasmas if the ionization fraction is sufficiently high ( ≳ 10 − 5 ). A consequence is that ionization causes disorder-induced heating (DIH), which triggers a significant rise in ion temperature on a picosecond timescale. This is followed by a rise in the neutral gas temperature on a longer timescale of up to nanoseconds due to ion–neutral temperature relaxation. The sequence of DIH and ion–neutral temperature relaxation suggests a new mechanism for ultrafast neutral gas heating. Previous work considered only the case of an instantaneous ionization pulse, whereas the ionization pulse extends over nanoseconds in many experiments. Here, molecular dynamics simulations are used to analyze the evolution of ion and neutral gas temperatures for a gradual ionization over several nanoseconds. The results are compared with published experimental results from a nanosecond pulsed discharge, showing good agreement with a measurement of fast neutral gas heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Isotope effects on transport in LHD

Abstract Isotope effects are one of the most important issues for predicting future reactor operations. Large helical device (LHD) is the presently working largest stellarator/helical device using super conducting helical coils. In LHD, deuterium experiments started in 2017. Extensive studies regarding isotope effects on transport have been carried out. In this paper, the results of isotope effect studies in LHD are reported. The systematic studies were performed adjusting operational parameters and nondimensional parameters. In L mode like normal confinement plasma, where internal and edge transport barriers are not formed, the scaling of global energy confinement time ( τ E ) with operational parameters shows positive mass dependence ( M 0.27 ; where M is effective ion mass) in electron cyclotron heating plasma and no mass dependence ( M 0.0 ) in neutral beam injection heating plasma. The non-negative ion mass dependence is anti-gyro-Bohm scaling. The role of the turbulence in isotope effects was also found by turbulence measurements and gyrokinetic simulation. Better accessibility to electron and ion internal transport barrier (ITB) plasma is found in deuterium (D) plasma than in hydrogen (H). Gyro kinetic non-linear simulation shows reduced ion heat flux due to the larger generation of zonal flow in deuterium plasma. Peaked carbon density profile plays a prominent role in reducing ion energy transport in ITB plasma. This is evident only in plasma with deuterium ions. New findings on the mixing and non-mixing states of D and H particle transports are reported. In the mixing state, ion particle diffusivities are higher than electron particle diffusivities and D and H ion density profiles are almost identical. In the non-mixing state, ion particle diffusivity is much lower than electron diffusivity. Deuterium and hydrogen ion profiles are clearly different. Different turbulence structures were found in the mixing and non-mixing states suggesting different turbulence modes play a role.

Physics↗
DOE OSTI2021

Testing the DIII-D co/counter off-axis neutral beam injected power and ability to balance injected torque

DIII-D has undergone a major upgrade and successfully injected high power off-axis neutral beams (~4 MW) in both co-current and counter-current directions. This capability of high power co/counter steerable off-axis neutral beams on a major tokamak opens a unique parameter space of broad pressure and current profiles for high beta steady-state advanced tokamak (AT) scenarios, while retaining the ability to balance the injected torque for low rotation studies. This co/counter off-axis neutral beam capability is being used to validate physics-based energetic particle and thermal transport models for designing next-step facilities based on the steady-state AT approach. This paper reports on the critical evaluation of the transmitted power and energetic ion population produced by this heating and current drive system, which is assessed through visible imaging, neutron measurements and rotation profile measurements at balanced torque. Minimal losses of neutral beam power have been achieved by optimizing the strongly focused ion sources required to pass through the aperture. Tilting of the ion source has been guided by fast visible imaging and resulted in neutral beam injection along the design centerline with empirical characterization of each beam's divergence derived from the imaging data and used in the NUBEAM description of the beam injection. Through exclusive power injection of each neutral beam into MHD quiescent plasmas across a range of neutral beam voltage, perveance and plasma current we conclude that a modest reduction (~10%–15%) of transmitted power compared to on-axis, standard focus has been incurred. We report corrections that more accurately represent the injected power. Good ability to balance the neutral beam torque has been demonstrated by injecting the new off-axis counter injecting beam against the existing off-axis co-injected beam in 2.0 T, 1.0 MA, MHD quiescent L-mode plasmas. Furthermore, the torque balance studies verify the ability to operate with balanced injection, which is critical for achieving low torque and low rotation operation for physics studies and in ITER demonstration discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparing the kinetics of ionized and neutral atoms from single and multi-element laser-produced plasmas

In this work, kinetics of ion and neutral atom emission features were compared for nanosecond laser-produced plasmas generated from several metal targets (i.e., Al, Ti, Zr, Nb, Ta) and an alloy containing all of these as principal alloying elements. Plasmas were produced by focusing 6 ns, 1064 nm pulses from an Nd:YAG laser on the targets of interest in a vacuum. A Faraday cup was used for collecting ion temporal features while spatially and temporally resolved emission spectroscopy was used for measuring the optical time of flight of various neutral atomic transitions. Our results highlight that most probable ion and atom velocities decay with increasing atomic mass. Trends for ions from the alloy target represent a weighted average where all ions contribute. For both ions and atoms, velocities decrease with increasing heat of vaporization and melting temperature, consistent with the thermal mechanisms that contribute to nanosecond laser ablation. Kinetic energies for neutral atoms from pure metal targets have some variability with atomic mass, whereas kinetic energies for atoms from the alloy target are more similar. These more similar kinetic energies observed for neutral atoms in the multi-element plasma may be attributed to collisions between species from all elements in the Knudsen layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A diagnostic to measure neutral-atom density in fusion-research plasmas

For this study, a femtosecond two-photon-absorption laser-induced-fluorescence (TALIF) diagnostic was designed, installed, and operated on the Princeton-Field-Reversed Configuration-2 device to provide non-invasive measurements of the time and spatially resolved neutral-atom densities in its plasmas. Calibration of the H o density was accomplished by comparison with Kr TALIF. Measurements on plasmas formed of either H 2 or Kr fill gases allowed examination of nominally long and short ionization mean-free-path regimes. With multi-kW plasma heating and H 2 fill gas, a spatially uniform H o density of order 10 17 m –3 was measured with better than ±2 mm and 10 µs resolution. Under similar plasma conditions but with Kr fill gas, a 3-fold decrease in the in-plasma Kr density was observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2022

How sheath properties change with gas pressure: modeling and simulation

Abstract Particle-in-cell simulations are used to study how neutral pressure influences plasma properties at the sheath edge. The high rate of ion–neutral collisions at pressures above several mTorr are found to cause a decrease in the ion velocity at the sheath edge (collisional Bohm criterion), a decrease in the edge-to-center density ratio ( h l factor), and an increase in the sheath width and sheath potential drop. A comparison with existing analytic models generally indicates favorable agreement, but with some distinctions. One is that models for the h l factor need to be made consistent with the collisional Bohm criterion. With this and similar corrections, a comprehensive fluid-based model of the plasma boundary transition is constructed that compares well with the simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗