Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma sheaths & surfaces”

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 91 records · Page 5

Magnetic potential based formulation for linear and non-linear 3D RF sheath simulation

Abstract This paper reports a new numerical scheme to simulate the radio-frequency (RF) induced RF sheath, which is suitable for a large 3D simulation. In the RF sheath boundary model, the tangential component of the electric field ( E t ) is given by the gradient of a scalar electric field potential. We introduce two additional scalar potentials for the tangential components of the magnetic field, which effectively impose the normal electric displacement ( D n ) on the plasma sheath boundary condition via in-homogeneous Neumann boundary condition and constrain the tangential electric field on the surface as curl-free ( ∇ × E t = 0 ). In our approach, the non-linear sheath impedance is formulated as a natural extension of the large thickness (or asymptotic) sheath limit ( D n = 0 ), allowing for handling both asymptotic and non-linear regimes seamlessly. The new scheme is implemented using the Petra-M finite element method analysis framework and is verified with simulations in the literature. The significance of non-linearity is discussed in various plasma conditions. An application of this scheme to asymptotic RF sheath simulation on the WEST ICRF antenna side limiters is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of the inverse sheath on divertor plasma performance in tokamak edge plasma simulations

Recently, it was shown that strong electron thermionic emission from material walls could result in the formation of an “inverse sheath”, which prevents the flow of cold ions to the wall. Such regimes look very favorably from the point of view of plasma-material interactions at the edge of magnetic fusion devices, where the problem of the erosion of plasma facing components under ion irradiation is one of the key issues for developing of future magnetic fusion reactors. However, it is not clear whether such regimes are compatible with edge plasma parameters and heat removal requirements in fusion reactors.To address the issue of practicality of the “inverse sheath” regime for edge tokamak plasma conditions, we perform a set of numerical simulations with 2D edge plasma transport code UEDGE for a DIII-D-like geometry and magnetic configuration.To describe both “standard” and “inverse sheath” conditions within the framework of the UEDGE code (which does not consider the sheath region per se), at the material surfaces, we apply effective boundary conditions which emulate both “standard” and“inverse sheath” regimes. We demonstrate that for the same input parameters, spatial distributions of edge plasma parameters corresponding to detached divertor and“inverse sheath” regimes are similar, with only a few minor differences. We discuss the compatibility of “inverse sheath” regimes with core plasma parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional charging effects

The physics of three dimensional spacecraft charging is reviewed. The asymmetric plasma sheath created around a spacecraft can significantly alter the trajectories of incoming particles or repel secondary particles emitted from the surface back to the surface. Thus seemingly anomolous potential variations can be encountered. Depending on the sensitivities of the spacecraft and its instruments to such variations, charging can jeopardize mission performance or even mission life through enhanced arcing. Three dimensional interactions between the Space Station and the ambient plasma and magnetic field are increasingly important.

Garrett, Henry Berry↗

Effects of low-frequency voltage on nonlinear standing wave excitation, plasma uniformity, and ion dynamics in dual-frequency asymmetric capacitive discharges

Abstract It is known that in very-high-frequency (VHF) capacitively coupled plasmas, the higher harmonics generated by nonlinear sheath motion can enhance the standing wave effect (SWE), which can lead to center-peaked plasma density profiles. In this work, an improved nonlinear electromagnetic model incorporating a transmission line model, an electron momentum balance model, a bulk plasma model, a collisionless nonlinear numerical sheath model, and an ion Monte-Carlo collision (MCC) model is developed to study the effects of low-frequency (LF) voltage V L on the nonlinear standing wave excitation, plasma uniformity, and ion energy and angular distribution functions (IEDFs and IADFs) in dual-frequency (DF) asymmetric capacitive argon discharges at relatively low pressure of 3 Pa. The plasma diffusion in the radial direction and ion dynamics within the LF oscillating sheath are self-consistently considered. The LF voltage V L at 2 MHz varies from 0 to 700 V while the HF voltage V H at 60 MHz is fixed at 100 V. Simulation results indicate that without the addition of an LF source (i.e. V L = 0 V), there are a considerable number of high-order harmonics with short wavelengths, leading to significant SWE and central peak in the radial plasma density profile. Nevertheless, the high-order harmonic excitations tend to be weakened and merely occur around the phase of the full LF sheath collapse due to a shorter characteristic damping time of the surface waves as V L increases. This, combined with increased surface wavelengths of both the driving frequency and the higher harmonics at a higher V L , leads to suppressed standing waves and improved plasma uniformity. Meanwhile, the simulations show that both the low and the high energy peaks of IEDF move towards higher energies, and the energy peak separation width ΔEbecomes wider with the increase of V L . The IEDF at the radial center of the powered electrode exhibits a broader ΔEthan that at the edge. For the IADF, an increased V L results in more ions incident on the electrode with a smaller deflection angle. Because of a thinner sheath and a higher sheath voltage at the electrode center, the peak value of IADF at the electrode center is greater than that at the edge.

Physics↗

A Coupled Plasma-Sheath Model for High Density Sources

High density, low pressure plasmas are used for etching and deposition in microelectronics fabrication processes. The process characteristics are strongly determined by the ion energy distribution (IED) and the ion flux arriving at the substrate that are responsible for desorption of etch products and neutral dissociation at the surface. The ion flux and energy are determined by a self- consistent modeling of the bulk plasma, where the ions and the neutral radicals are produced, and the sheath, where the ions are accelerated. Due to their widely different time scales, it is a formidable task to self-consistently resolve non-collisional sheath in a high density bulk plasma model. In this work, we first describe a coupled plasma-sheath model that attempts to resolve the non-collisional sheath in a reactor scale model. Second, we propose a semianalytical radio frequency (RF) sheath model to improve ion dynamics.

Bose, Deepak↗

NASA Talks: Harsh Lunar Environments

NASA Talks - Exploration drives technology, and as NASA works toward sustained operations on the lunar surface for a permanent moonbase, it is critical to understand how the location drives what’s needed to accomplish these goals. This is especially important as we land on the lunar South Pole, which will be different and riskier than the previous Apollo landing sites. The thermal conditions will be much colder, the lighting environment is challenged with low sun angles, the regolith holds a charge in the plasma sheath, and the radiation is significantly higher than in low Earth orbit (LEO). Stephanie Sipila, Deputy Manager of Extravehicular Activity and Human Surface Mobility Program (EHP) Technology, Integration and Partnerships Office, provided knowledge and insight into the game-changing technologies NASA is exploring to mitigate these risks.

Stephanie Sipila↗

Quantitative Measurement of Positive and Negative Ion Species Ejected from a Li–O–H Surface by Hydrogen and Noble Gas Ion Irradiation

We report sputtering yields of Li + , H - , O - , and OH x - ion species from an Li–O–H surface for H, D, He, Ne, and Ar ion irradiation at 45° incidence in the energy range of 30–2000 eV. A Li film was deposited on a stainless steel target using Li evaporators in the LTX-β vessel, using the LTX-β Sample Exposure Probe (SEP), which includes an ultrahigh vacuum suitcase for transferring targets without significant contamination from air exposure. The SEP was used to transfer the Li-coated target from LTX-β to a separate Sample Exposure Station (SES) to perform ion exposure measurements. The SEP was also used for characterization of the Li-coated target utilizing X-ray photoelectron spectroscopy in a different chamber, showing that the lithium film surface was oxidized. Ion exposures were performed using an electron cyclotron resonance plasma source in the SES. Sputtered/ejected species were sampled by a quadrupole mass spectrometer with capabilities for detecting positive and negative ions, and an energy filter for determining the mean kinetic energy of the ejected ion species. All ion irradiations caused Li + ions to be ejected, while causing impurity ions such as H + , H - , O - and OH - to be ejected. Measured ion energies of Li + ions from a Li–O–H surface suggested that the typical sheath potential on the divertor surface can trap sputtered Li + ions, which were previously reported as ~ 60% of total sputtered Li species from Li targets (Allain and Ruzic in Nucl Fusion 42:202, 2002). Hence, our results for the sputtering yields of ejected ion species and their associated ion energies from a Li–O–H surface indicates that lithium sputtering is suppressed and impurity removal is enhanced due to the sheath potential at the divertor surface for fusion reactor applications.

Plasma-material interaction↗

Multi-scale dynamics of atmospheric-pressure discharges ignited over liquid electrodes

A multiphase computational model of atmospheric-pressure gas discharges ignited over liquid electrodes is developed. The model takes into account both the liquid electrode motion under applied electric fields and plasma generation in the gas phase. The influence of the applied voltage and liquid properties (density, viscosity, and surface tension) on the liquid surface deformation is analyzed. It is shown that the shape of dynamic cones formed on the liquid surface differs from the static Taylor's cones. The influence of the liquid surface protrusions on gas breakdown dynamics is demonstrated. It is shown that the breakdown develops in two stages: first, a fast ionization wave propagates from the cathode to the anode, and then the cathode sheath collapses due to secondary electron emission from the liquid surface. It is shown that dynamics of the multiphase system containing gas, plasma, and liquid states is characterized by three disparate time scales: the fast electron time scale, the intermediate ion time scale, and the slow liquid dynamics time scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy-angle distribution of the ions in the RF sheath of ICRH antennas

Radio-frequency sheaths forming at the surface of ion cyclotron resonance heating (ICRH) actuators in fusion experiments are associated with enhanced impurity sputtering from ICRH plasma-facing components (PFCs). The minimization of impurity fluxes from the ICRH PFC is a vital task for the usage of ICRH systems. Capturing the ion kinetics is important for an accurate description of plasma–material interactions, because the ion dynamics plays a crucial role in RF sheaths. Here, we present a hybrid particle-in-cell (hPIC) model able to capture the kinetic behavior of the ions. We analyze the kinetic ion energy-angle distributions (IEADs) impacting the RF antenna and its dependence on different plasma and RF sheath parameters. In particular, the IEAD dependence on RF frequency and magnetic field alignment is investigated. Using hPIC, we simulated a case emulating the latest experimental campaign from JET. The simulation showed that under specific plasma and RF parameters, the kinetic motion of ions results in a cusp formation in the IEAD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Validation and applications of the POLAR code

The Potentials of Large Objects in the Auroral Region (POLAR) code was written to model the charging of large spacecraft by auroral electrons in low polar orbit. It has also proven to be a powerful tool for analysis of active and voltage-plasma interactions. POLAR combines 3-D object definition, direct solution to Poisson's equation, space charge transport, and surface charge computation. Physical effects treated by POLAR include the plasma wake and ram, the geomagnetic field, ionization, sheath and barrier formation, and photo and secondary electron generation and transport. Polar development is complete and it is now available to the aerospace community. The methodology of POLAR is reviewed but concentration is on the applications of POLAR to physical problems and observations with the intent of demonstrating both the code's capabilities, and its limitations. These applications include the Space Shuttle, the Defense Meteorological Satellite Program (DMSP) satellites, sounding rockets, plasma chamber experiments, and things not yet flown.

Cooke, David L.↗

Resonant absorption of Alfvenic surface waves and the heating of solar coronal loops

A first-order mass- and energy-balance model is developed for steady-state EUV 'coronal rain' loops that are not associated with postflare events and are often seen over sunspot umbrae. The model disregards variations in a loop's thermodynamic and magnetic properties along magnetic-field lines and yields average properties of the resonant absorption sheath, the boundary layer, and the loop's interior. Both irreversible heating by Alfvenic surface waves and energy transport via induced boundary-layer convection are taken into account in the analysis. Results that include predictions of the temperatures of the sheath and the boundary-layer plasma, the temperature of the interior plasma, the radiative output of the loop, and the filling factor associated with this radiation are given in terms of the period and velocity amplitude of the relevant surface wave as well as various parameters that characterize a typical 'coronal rain' loop. These results are shown to be consistent with the observed minimum and maximum temperatures as well as the radiative output of typical 'coronal rain' loops.

Ionson, J. A.↗

Computational investigation of incident ion angles and material erosion at rough graphite and silicon carbide divertor surfaces

We present a computational investigation of the dependence of material erosion on the incident ion angle at rough graphite and silicon carbide divertor surfaces. Ion angle distributions (IADs) for D plasmas at NSTX-U and DIII-D divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from NSTX-U and DIII-D experiments. The calculations found that the effective sputtering yields, the sputtering pattern, and the shadowed area are determined by the detailed surface topology rather than the root mean square roughness R RMS , which represents deviations from a flat surface. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA) (θ'). The mean surface inclination angle distribution (SIAD) (δ) was found to be a useful parameter to estimate the LIIA from the calculated IADs. We report global empirical formulas for the mean LIIA and fraction of the area shadowed from the main ions for D plasmas for rough surfaces with B-field incident angles α = 85°–89° as a function of the mean SIAD (δ). We propose the use of the mean LIIA (θ') to estimate the sputtering yield for rough surfaces from the angular dependence of the sputtering yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vidyut3d: A Non-Equilibrium Plasma Modeling Tool [SWR-24-101]

Vidyut3d is a massively-parallel plasma-fluid solver for low-temperature plasmas (LTPs) that supports both local field (LFA) and local mean energy (LMEA) approximations, as well as complex gas and surface-phase chemistry. The solver supports 2D and 3D domains, and uses AMReX's adaptive mesh refinement capabilities to increase the grid resolution around complex structures (e.g. streamer heads and sheaths) while maintaining a tractable problem size. Vidyut specializes in simulating various types of gas-phase discharges, as well as plasma/surface interactions and surface chemistry (e.g. for plasma-mediated catalysis applications). The solver also supports hybrid CPU/GPU parallelization strategies, and has demonstrated excellent scaling on various HPC architectures for problem sizes consisting of O(100 M) control volumes.

Sitaraman, Hariswaran↗

Electrostatic charging of spacecraft in response to electron beam injection

Electron beam injections from spacecraft now constitute a major activity in space research. Here, the charging level of a conducting surface when an electron beam is injected from it is investigated. Injections into both vacuum and an ambient plasma are considered. When a Maxwellian beam is injected into vacuum, the surface changes to a potential much greater than the average beam energy. The dependence of this excess is examined by considering beams with water-bag types of velocity distribution functions in which no electron has a velocity V(max) above a certain value. The electric field distribution in the electron sheath near the surface is determined by the pressure distribution. Thus, the surface potential is determined not only by V(max) but by all the beam parameters. The ambient plasma reduces the charging level and causes an oscillation in the surface potential. The oscillation frequency is the electron-plasma frequency associated with the ambient plasma.

Singh, Nagendra↗

Anode power deposition in quasi-steady MPD arcs

The power deposited in the anode of a quasi-steady MPD accelerator has been measured directly by thermocouples attached to the inside surface of a shell anode which provide a local measurement of anode heat flux. The results over a range of arc currents from 5.5 to 44 kiloamperes and argon mass flows from 1 g/sec to 48 g/sec show that the fraction of the total input power deposited in the anode decreases drastically from 50% at an arc power of 200 kW to 10% at 20 MW, and that anode power is not uniformly deposited in the anode. A theoretical model of the anode heat transfer, including effects of anode work function, electron thermal energy, and anode sheath, can be brought into reasonable agreement with the measurements, provided the effective range of the conduction electrons from within the discharge plasma to the anode surface is properly acknowledged.

Saber, A. J.↗

Current-Produced Magnetic Field Effects on Current Collection

Current collection by an infinitely long, conducting cylinder in a magnetized plasma, taking into account the magnetic field of the collected current, is discussed. A region of closed magnetic surfaces disconnects the cylinder from infinity. Due to this, the collected current depends on the ratio between this region and the plasma sheath region and, under some conditions, current reduction arises. It is found that the upper-bound limit of current collection is reduced due to this change of magnetic field topology. The effect can be essential even if the orbit-limited model of current collection is valid. This model is used to find the reduction of the total current collected by a cylinder (e.g., a bare tether). Such effect strongly depends on plasma density. The results are applied to a tether system in the ionosphere. For this case, it is found that current reduction can be significant for long tethers in typical dayside ionospheric conditions.

Khazanov, G. V.↗

Current-Produced Magnetic Field Effects on Current Collection

Current collection by an infinitely long, conducting cylinder in a magnetized plasma, taking into account the magnetic field of the collected current, is discussed. A region of closed magnetic surfaces disconnects the cylinder from infinity. Owing to this, the collected current depends on the ratio between this region and the plasma sheath region and, under some conditions, current reduction arises. It is found that the upper bound limit of current collection is reduced due to this change of magnetic field topology. The effect can be substantial even if the orbit-limited model of current collection is valid. This model is used to find the reduction of the total current collected by a cylinder (e.g., a bare tether). It is shown that this effect strongly depends on plasma density. The results are applied to a tether system in the ionosphere, and it is found that current reduction can be significant for long tethers in typical dayside ionospheric conditions.

Khazanov, G. V.↗

Current-Produced Magnetic Field Effects on Current Collection

Current collection by an infinitely long, conducting cylinder in a magnetized plasma, taking into account the magnetic field of the collected current, is discussed. A region of closed magnetic surfaces disconnects the cylinder from infinity. Due to this, the collected current depends on the ratio between this region and the plasma sheath region and, under some conditions, current reduction arises. The current collection along a realistic "bare wire" space tether is considered. A number of factors are taken into account, including the resistance of the wire and shielding resulting from the current-induced magnetic field produced by current flow in the tether. The plasma density, tether length and radius, the geomagnetic field strength and angle to the orbital velocity vector were all used as parameters in the study. It is shown that magnetic shielding for certain tether system configurations, when combined with particular values of the governing parameters, significantly reduces the collected current. Specifically, it is shown that an electrodynamic tether in the "thruster" mode suffers greater reduction from magnetic shielding than a tether with the same characteristics deployed in the "generator" mode. We find that, for both modes, current-induced magnetic shielding becomes more significant as plasma density and wire radius increase. The same is true for the dependence on the angle of the geomagnetic field to the orbital velocity vector and the motion-induced electric field for the generator mode For the thruster mode, the effect is larger for smaller angles. In both operating modes, the shielding is more important for smaller angles between the tether and magnetic field. In addition to the above dependencies, the effect for the thruster mode essentially depends on the tether length. In general, any parametric change that increases tether current, relative to the strength of the electric field between the tether and the ambient plasma, will increase the shielding effect (thereby reducing overall tether current compared to the case without considering these effects). The effect of the induced magnetic shielding can be significant in some cases.

Khazanov, George V.↗