Sheath structure in the unignited mode.
Kinetic theory model of plasma sheath transition, analyzing sheath structure in unignited mode of thermionic converter
Sockol, P. M.↗
SEARCH · Search NASA
Results for “Plasma Sheath”
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.
Sheath structure in the unignited mode.
Kinetic theory model of plasma sheath transition applied to sheath structure in unignited mode of thermionic converter
Sockol, P. M.↗
Particle-in-cell simulations and passive bulk model of collisional capacitive discharge
Intermediate pressure (0.2–6.0 Torr), radio frequency capacitive discharges are of increasing importance in the thin film processing industry. However, existing low pressure (<0.1 Torr) discharge models are often used beyond their applicability in this pressure range. In this work, the authors performed one-dimensional particle-in-cell (PIC) simulations of a capacitive 2.5 cm gap argon discharge driven at 13.56 MHz with current density amplitude J 0 50 A/m 2 over a wide range of gas pressures p g = 0.04–20 Torr. For intermediate pressures, the PIC results showed that most of the ionization occurs near the plasma-sheath edges with very little occurring within the central bulk region. This led us to develop a “passive bulk” model of the discharge in which bulk electron heating does not contribute to ionization, and all the ionization is due solely to the sheath-heated electrons. The model assumes a constant ion-neutral mean free path λ i in the plasma bulk and either constant λ i or constant ion mobility μ i in the sheaths. The model is solved using an electron temperature iteration scheme. The PIC scalings of the discharge equilibrium parameters (e.g., electron density and temperature, sheath width, and voltage) with the input parameters p g and J 0 showed good agreement with those of the passive bulk model. A constant μ i sheath model showed the best agreement with the PIC results in the intermediate pressure regime of interest, but a constant λ i model is simpler to implement and showed reasonable agreement over a wide pressure range.
36 MATERIALS SCIENCE↗
Spacecraft-induced plasma energization and its role in flow phenomena
Plasma instabilities induced by orbiting vehicles can cause many important phenomena ranging from electron and ion heating and suprathermal electron tail energization, to enhanced ionization and optical emissions. We outline the basic collective processes leading to plasma energization near plasma sheaths and in regions of neutral gas streaming through plasma, and discuss the role of the induced collective effects in producing the optical emission spectra.
Papadopoulos, K.↗
Ion extraction from a plasma
An experimental investigation of the physical processes governing ion extraction from a plasma is presented. The screen hole plasma sheath of a multiaperture ion accelerator system is defined by equipotential plots for a variety of accelerator system geometries and operating conditions. A sheath thickness of at least fifteen Debye lengths is shown to be typical. The electron density variation within the sheath satisfies a Maxwell Boltzmann density distribution at an effective electron temperature dependent on the discharge plasma primary to Maxwellian electron density ratio. Plasma ion flow up to and through the sheath is predominately one dimensional and the ions enter the sheath with a modified Bohm velocity. Low values of the screen grid thickness to screen hole diameter ratio give good ion focusing and high extracted ion currents because of the effect of screen webbing on ion focusing.
Aston, G.↗
Simulations of Solar Wind Plasma Flow Around a Simple Solar Sail
In recent years, a number of solar sail missions of various designs and sizes have been proposed (e.g., Geostorm). Of importance to these missions is the interaction between the ambient solar wind plasma environment and the sail. Assuming a typical 1 AU solar wind environment of 400 km/s velocity, 3.5 cu cm density, ion temperature of approx.10 eV, electron temperature of 40 eV, and an ambient magnetic field strength of 10(exp -4) G, a first order estimate of the plasma interaction with square solar sails on the order of the sizes being considered for a Geostorm mission (50 m x 50 m and 75 m x 75 m corresponding to approx.2 and approx.3 times the Debye length in the plasma) is carried out. First, a crude current balance for the sail surface immersed in the plasma environment and in sunlight was used to estimate the surface potential of the model sails. This gave surface potentials of approx.10 V positive relative to the solar wind plasma. A 3-D, Electrostatic Particle-in-Cell (PIC) code was then used to simulate the solar wind flowing around the solar sail. It is assumed in the code that the solar wind protons can be treated as particles while the electrons follow a Boltzmann distribution. Next, the electric field and particle trajectories are solved self-consistently to give the proton flow field, the electrostatic field around the sail, and the plasma density in 3-D. The model sail was found to be surrounded by a plasma sheath within which the potential is positive compared to the ambient plasma and followed by a separate plasma wake which is negative relative to the plasma. This structure departs dramatically from a negatively charged plate such as might be found in the Earth s ionosphere on the night side where both the plate and its negative wake are contiguous. The implications of these findings are discussed as they apply to the proposed Geostorm solar sail mission.
Garrett, Henry B.↗
Collisionless sheath near an electrode in the presence of a magnetic field
Magnetic field effects on collisionless plasma sheath near planar electrode, solving electric and magnetic potentials nonlinear differential equations by numerical scheme
Hassan, H. A.↗
Computer models of the spacecraft wake
Until recently, computations of space plasma flow over a spacecraft have been unstable for ratios of spacecraft dimension to Debye length typical of the low Earth orbit environment. Calculations are presented of the spacecraft/environment interaction based on two computer codes, MACH and POLAR. MACH, an inside-out particle tracking code, was developed for the purpose of validating the physics of POLAR in regimes where these are no comprehensive theoretical or experimental results. While the spacecraft which can be treated by MACH are restricted to simple geometries, the methodology is more fundamental than POLAR. MACH generates self-consistent solutions within the context of quasisteady Vlasov plasma flow and achieves Debye ratios previously unobtainable. POLAR uses a three-dimensional finite-element representation of the vehicle in a staggered mesh. The plasma sheath is modeled by outside-in particle tracking. Solutions for the plasma flow, wake and vehicle charging are obtained by Vlasov-Poisson iteration; charge stabilization techniques make the results virtually insensitive to the Debye ratio. POLAR reproduces the Laframboise static plasma solutions for sperical probes and fits the Makita-Kuriki probe data for spheres in a flowing plasma in regions where comparisons are valid. POLAR and MACH solutions for the particle and electrostatic potential structure of the wake of a charged disk in a low-altitude flow are shown for Mach numbers 4, 5, and 8. New features of the solutions include ion focussing in the wake and a definitive determination of the sheath edge in the wake which shows that the sheath is not an equipotential.
Rubin, A. G.↗
Similarity for downscaled kinetic simulations of electrostatic plasmas: Reconciling the large system size with small Debye length
A simple similarity has been proposed for kinetic (e.g., particle-in-cell) simulations of plasma transport that can effectively address the long-standing challenge of reconciling the tiny Debye length with the vast system size. This applies to both transport in unmagnetized plasma and parallel transport in magnetized plasmas, where the characteristics length scales are given by the Debye length, collisional mean free paths, and the system or gradient lengths. The controlled scaled variables are the configuration space, x/L, and an artificial Coulomb Logarithm, L ln Λ, for collisions, while the scaled time, t/L, and electric field, LE, are automatic outcomes. The similarity properties are examined, demonstrating that the macroscopic transport physics is preserved through a similarity transformation while keeping the microscopic physics at its original scale of Debye length. To showcase the utility of this approach, two examples of 1D plasma transport problems were simulated using the VPIC code: the plasma thermal quench in tokamaks [Li et al., Nuclear Fusion 63, 066030 (2023)] and the plasma sheath in the high-recycling regime [Li et al., Physics of Plasmas 30, 063505 (2023)].
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Simulations of divertor plasmas with inverse sheaths
The effect of strong electron emission from material surfaces has been proposed to form an “inverse sheath”: a region with a positive potential relative to the near-wall plasma which prevents the flow of ions to the wall [M. D. Campanell, “Negative plasma potential relative to electronemitting surfaces,” Phys. Rev. E. 88, 033103 (2013); M. D. Campanell and M. V. Umansky, “Strongly emitting surfaces unable to float below plasma potential,” Phys. Rev. Lett. 116, 1–5 (2016); M. D. Campanell and G. R. Johnson, “Thermionic cooling of the target plasma to a sub-ev temperature,” Phys. Rev. Lett. 122, 1–5 (2019)]. We assess the viability of this regime in a tokamak device using the 2D edge plasma transport code UEDGE [T. Rognlien et al ., “A fully implicit, time dependent 2-D fluid code for modeling tokamak edge plasmas,” J. Nucl. Mater. 196–198, 347–351 (1992)]. Since the UEDGE code does not consider the sheath region directly, we apply boundary conditions at the divertor targets which emulate the physics of both “standard” and “inverse” sheath regimes [R. Masline et al ., “Influence of the inverse sheath on divertor plasma performance in tokamak edge plasma simulations,” Contrib. Plasma Phys. 60, e201900097 (2020)]. Using these boundary conditions, we perform scoping studies to assess plasma parameters near the target by varying the density at the core-edge interface. We observe a smooth transition in the resultant profiles of plasma parameters for the standard sheath, and a bifurcation across the simulation set for plasmas with an inverse sheath. The cause of this bifurcation is assessed by performing the parameter scan both with and without impurity radiation; we observe that the bifurcation persists in both cases, indicating that this bifurcation is caused by plasma recombination.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Surface voltage gradient role in high voltage solar array-plasma interaction: Center Director's discretionary fund
A large amount of experimental and analytical effort has been directed toward understanding the plasma sheath growth and discharge phenomena which lead to high voltage solar array-space plasma interactions. An important question which has not been addressed is how the surface voltage gradient on such an array may affect these interactions. The results of this study indicate that under certain conditions, the voltage gradient should be taken into account when evaluating the effect on a solar array operating in a plasma environment.
Carruth, M. R., Jr.↗
Rethinking Use of the OML Model in Electric Sail Development
In 1924, Irvin Langmuir and H. M. Mott-Smith published a theoretical model for the complex plasma sheath phenomenon in which they identified some very special cases which greatly simplified the sheath and allowed a closed solution to the problem. The most widely used application is for an electrostatic, or "Langmuir," probe in laboratory plasma. Although the Langmuir probe is physically simple (a biased wire) the theory describing its functional behavior and its current-voltage characteristic is extremely complex and, accordingly, a number of assumptions and approximations are used in the LMS model. These simplifications, correspondingly, place limits on the model's range of application. Adapting the LMS model to real-life conditions is the subject of numerous papers and dissertations. The Orbit-Motion Limited (OML) model that is widely used today is one of these adaptions that is a convenient means of calculating sheath effects. Since the Langmuir probe is a simple biased wire immersed in plasma, it is particularly tempting to use the OML equation in calculating the characteristics of the long, highly biased wires of an Electric Sail in the solar wind plasma. However, in order to arrive at the OML equation, a number of additional simplifying assumptions and approximations (beyond those made by Langmuir-Mott-Smith) are necessary. The OML equation is a good approximation when all conditions are met, but it would appear that the Electric Sail problem lies outside of the limits of applicability.
Stone, Nobie H.↗
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↗
Unique Results and Lessons Learned from the TSS Missions
In 1924, Irvin Langmuir and H. M. Mott-Smith published a theoretical model for the complex plasma sheath phenomenon in which they identified some very special cases which greatly simplified the sheath and allowed a closed solution to the problem. The most widely used application is for an electrostatic, or "Langmuir," probe in laboratory plasma. Although the Langmuir probe is physically simple (a biased wire) the theory describing its functional behavior and its current-voltage characteristic is extremely complex and, accordingly, a number of assumptions and approximations are used in the LMS model. These simplifications, correspondingly, place limits on the model's range of application. Adapting the LMS model to real-life conditions is the subject of numerous papers and dissertations. The Orbit-Motion Limited (OML) model that is widely used today is one of these adaptions that is a convenient means of calculating sheath effects. The OML equation for electron current collection by a positively biased body is simply: I is approximately equal to A x j(sub eo) x 2/√𝜋(phi)(exp ½) where A is the area of the body and phi is the electric potential on the body with respect to the plasma. Since the Langmuir probe is a simple biased wire immersed in plasma, it is particularly tempting to use the OML equation in calculating the characteristics of the long, highly biased wires of an Electric Sail in the solar wind plasma. However, in order to arrive at the OML equation, a number of additional simplifying assumptions and approximations (beyond those made by Langmuir-Mott-Smith) are necessary. The OML equation is a good approximation when all conditions are met, but it would appear that the Electric Sail problem lies outside of the limits of applicability.
Stone, Nobie H.↗
Aligning thermal and current quenches with a high density low-Z injection
The conventional approach for thermal quench (TQ) mitigation in a tokamak disruption is through a high- Z impurity injection that radiates away the plasma’s thermal energy before it reaches the wall. The downside is a robust Ohmic-to-runaway current conversion due to the radiatively clamped low post-thermal-quench electron temperature. An alternative approach is to deploy a low- Z (either deuterium or hydrogen) injection that aims to slow down the TQ, and ideally aligns it with the current quench (CQ). This approach has been investigated here via 3D MHD simulations using the PIXIE3D code. By boosting the hydrogen density, a fusion-grade plasma is dilutionally cooled at approximately the original pressure. Energy loss to the wall is controlled by a Bohm outflow condition at the boundary where the magnetic field intercepts a thin plasma sheath at the wall, in addition to Bremsstrahlung bulk losses. Robust MHD instabilities proceed as usual, while the collisionality of the plasma has been greatly increased and parallel transport is now in the Braginskii regime. The main conclusion of this study is that the decreased transport loss along open field lines due to a sufficient low- Z injection slows down the TQ rate to the order of 20 ms, aligned with the CQ timescale for a 15 MA ITER plasma.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Electric potential on solid spheres in a plasma
Derivation of the general expression for the potential on a solid sphere immersed in a plasma, showing the dependence of the potential on the radius (a) of the sphere and the radius (s) of the plasma sheath that develops around the sphere. In the limit where the radius a is much larger than the sheath thickness s-a, the well-known result for the potential on an infinite wall in contact with a plasma is recovered from this expression. At the other extreme, where s is much larger than a, the result derived by Spitzer (1941) for the potential on spherical grains in the interstellar plasma is obtained. Since the surface of the sphere forms a sink for the charged particles, there is a net drift of the plasma towards the surface. The effect of this drift on the potential is examined. Finally, for very small metallic spheres, an effect leading to a revision of the potential is discussed. This effect consists in a lowering of the potential barrier for the electrons due to the image force. The various effects limiting the potential on spheres are discussed.
De, B. R.↗
A computer model of solar panel-plasma interactions
High power solar arrays for satellite power systems are presently being planned with dimensions of kilometers, and with tens of kilovolts distributed over their surface. Such systems face many plasma interaction problems, such as power leakage to the plasma, particle focusing, and anomalous arcing. These effects cannot be adequately modeled without detailed knowledge of the plasma sheath structure and space charge effects. Laboratory studies of 1 by 10 meter solar array in a simulated low Earth orbit plasma are discussed. The plasma screening process is discussed, program theory is outlined, and a series of calibration models is presented. These models are designed to demonstrate that PANEL is capable of accurate self consistant space charge calculations. Such models include PANEL predictions for the Child-Langmuir diode problem.
Cooke, D. L.↗
3-dimensional current collection model
A three-dimensional, time dependent current collection model of a satellite has been developed for the TSS-1 system. The system has been simulated particularly for the Research of Plasma Electrodynamics (ROPE) experiment. The Maxwellian distributed particles with the geomagnetic field effects are applied in this numerical simulation. The preliminary results indicate that a ring current is observed surrounding the satellite in the equatorial plane. This ring current is found between the plasma sheath and the satellite surface and is oscillating with a time scale of approximately 1 microsec. This is equivalent to the electron plasma frequency. An hour glass shape of electron distribution was observed when the viewing direction is perpendicular to the equatorial plane. This result is consistent with previous findings from Linson (1969) and Antoniades et al. (1990). Electrons that are absorbed by the satellite are limited from the background ionosphere as indicated by Parker and Murphy (1967).
Hwang, Kai-Shen↗