Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELECTROSTATIC 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

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗

Non-linear Plasma Wake Growth of Electron Holes

An object’s wake in a plasma with small Debye length that drifts across the magnetic field is subject to electrostatic electron instabilities. Such situations include, for example, the moon in the solar wind wake and probes in magnetized laboratory plasmas. The instability drive mechanism can equivalently be considered drift down the potential-energy gradient or drift up the density gradient. The gradients arise because the plasma wake has a region of depressed density and electrostatic potential into which ions are attracted along the field. The non-linear consequences of the instability are analysed in this paper. At

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-organization of photoionized plasmas via kinetic instabilities

Abstract Self-organization in an unmagnetized collisionless plasma (in this paper) refers to formation of transient coherent structures such as collective oscillations (electrostatic waves) or magnetic fields resulting from so-called kinetic effects in the plasma. This topical review provides a comprehensive analysis of the self-organization of strong-field photoionized, non-equilibrium plasmas through kinetic instabilities. The authors propose and demonstrate a novel experimental platform that enables the formation of dense plasmas with known highly anisotropic and non-thermal electron velocity distribution functions on a timescale on the order of an inverse electron plasma frequency. We then show that such plasmas are highly susceptible to a hierarchy of kinetic instabilities, including two-stream, current filamentation and Weibel, that convert a fraction of the electron kinetic energy into electric and/or magnetic energy stored in self-organized structures. The electrostatic waves so produced are measured using a collective light (Thomson) scattering technique with femtosecond resolution as the kinetic instabilities aided by collisions eventually thermalize the plasma electrons. In addition, we describe a novel experimental technique that has made it possible to map the temporal evolution of the wavenumber spectrum of the thermal Weibel instability with picosecond resolution, which leads to the formation of quasi-static coherent magnetic fields with different topologies in photoionized plasmas. Finally, the paper summarizes the important results and discusses future directions on this topic.

Zhang, Chaojie (ORCID:0000000180353014)↗

Resonance cones in cold plasma: Origin, singularities, and power flow

In magnetized tenuous plasma, typical at the plasma edge of fusion devices, a nearly electrostatic wave mode with relatively enhanced electric field can propagate along a specific angle with the magnetic field. For this characteristic, it is known as a “resonance cone.” For instance, these waves can be excited by radio frequency antennas in the ion-cyclotron and lower-hybrid range of frequencies. We consider the resonance cones emitted by idealized spatially extended sources. In 2D, we use a novel geometric construction which generalizes the d'Alembert solution to curved boundaries/moving sources, and show, for the first time, that singular electric fields arise under these conditions, thereby bringing the resonance cones in line with the other resonances of the cold plasma theory. Still in 2D, we give an expression for the amount of power radiated by resonance cones in terms of surface quantities on the source, which is finite despite the singular electric field. In conclusion, we generalize the conclusions regarding the presence and location of singular electric fields to the 3D electromagnetic case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sheath constraints on turbulent magnetised plasmas

Abstract A solid target in contact with a plasma charges (negatively) to reflect the more mobile species (electrons) and thus keep the bulk plasma quasineutral. To shield the bulk plasma from the charged target, there is an oppositely (positively) charged sheath with a sharp electrostatic potential variation on the Debye length scale λ D . In magnetised plasmas where the magnetic field is inclined at an oblique angle α with the target, some of the sheath potential variation occurs also on the ion sound gyroradius length scale ρ S cos ⁡ α , caused by finite ion gyro-orbit distortion and losses. We consider a collisionless and steady-state magnetised plasma sheath whose thickness l ms ∼ max ( λ D , ρ S cos ⁡ α ) is smaller than the characteristic length scale L of spatial fluctuations in the bulk plasma, such that the limit l ms / L → 0 is appropriate. Spatial structures are assumed to be magnetic field-aligned. In the case of small magnetic field angle α ∼ δ ≡ ρ S / L ≪ 1 , electric fields tangential to the target transport ions towards the target via E × B drifts at a rate comparable to the one from parallel streaming. A generalised form of the kinetic Bohm–Chodura criterion at the sheath entrance is derived by requiring that the sheath electric field have a monotonic spatial decay far from the target. The criterion depends on tangential gradients of potential and ion distribution function, with additional nontrivial conditions.

kinetic Bohm–Chodura condition↗

Global gyrokinetic particle simulations of microturbulence in W7-X and LHD stellarators

Global gyrokinetic particle simulations of electrostatic ion temperature gradient (ITG) instability show that the most unstable eigenmode is localized to some magnetic fieldlines or discrete locations on the poloidal plane in the Wendelstein 7-X (W7-X) stellarator due to its mirror-like magnetic fields, which vary strongly in the toroidal direction and induce coupling of more toroidal harmonics (n) to form the linear eigenmode than in the Large Helical Device (LHD) stellarator. Nonlinear electrostatic simulation results show that self-generated zonal flows are the dominant saturation mechanism for the ITG instabilities in both the LHD and W7-X. Additionally, radial widths of the fluctuation intensity in both the LHD and W7-X are significantly broadened from the linear phase to the nonlinear phase due to turbulence spreading. Finally, nonlinear spectra in the W7-X are dominated by low-n harmonics, which can be generated both by nonlinear toroidal coupling of high-n harmonics and by linear toroidal coupling with large amplitude zonal flows due to the 3D equilibrium magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experiments and Simulations of the Dusty Properties of Hypervelocity Impact Plasmas

Hypervelocity impacts, referring to impacts occurring at several times the speed of sound in the target, result in inertial stresses greatly exceeding the material strength and produce behavior that is not fully understood. These impact phenomena are seen in space through impacts of orbital debris and meteoroids on spacecraft and airless natural bodies but can also be produced in laboratory environments using light-gas guns and electrostatic dust accelerators. When a hypervelocity impact occurs, plasma is generated both from thermal ionization and from pressure ionization, producing a wide range of plasma densities depending on the impactor’s velocity. This plasma is initially formed in the warm dense matter (WDM) regime and governed by high energy density (HED) hydrodynamics before it rapidly expands and dissipates, spanning many orders of magnitude in density and length scale. The plasma can also contain a dust component, which is particularly relevant for the lower-velocity range of hypervelocity impacts with parameters that are particularly relevant to Tokamak plasmas. Therefore, the properties and dynamics of these plasmas cannot be fully characterized without considering dusty plasma effects, which are ubiquitous but poorly understood.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Barkas effect in strongly magnetized plasmas

Strongly magnetized plasmas, which are characterized by the particle gyrofrequency exceeding the plasma frequency, exhibit novel transport properties. For example, recent work showed that the friction force on a test charge moving through a strongly magnetized plasma not only consists of the typical stopping power component but also includes components perpendicular to the test charge's velocity. However, these studies only considered test charges that have the same sign as the charge of the plasma particles. Here, we extend these calculations to the case of charges with opposite signs (such as an ion interacting with strongly magnetized electrons). This is done with both a novel generalized Boltzmann kinetic theory and molecular dynamics simulations. It is found that the friction force changes dramatically depending on the sign of the interacting charges. Likewise, the stopping power component for oppositely charged particles decreases in magnitude compared with like-charged particles, and the perpendicular components increase in magnitude. Moreover, the difference between the two cases increases as the gyrofrequency becomes larger compared with the plasma frequency. Finally, the electrical resistivity is calculated from the friction force, where it is found that strong magnetization in conjunction with oppositely charged interactions significantly decreases the parallel resistivity and increases the perpendicular resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron dynamics in radio frequency magnetron sputtering argon discharges with a dielectric target

Abstract We demonstrate a self-consistent and complete description of electron dynamics in a typical electropositive radio frequency magnetron sputtering (RFMS) argon discharge with a dielectric target. The electron dynamics, including the electron power absorption dynamics in one radio frequency (RF) period, is studied via a fully kinetic 2d3v particle-in-cell/Monte Carlo collision (PIC/MCC) electrostatic simulation. The interplay between the fundamental plasma parameters is analyzed through their spatiotemporal dynamics. Due to the influence of magnetic trap on the electron transport, a spatially dependent charging that perturbs the electric potential is observed on the dielectric target surface, resulting in a spatially dependent ion energy distribution along the target surface. The E × B drift-to-discharge current ratio is in approximate agreement with Bohm diffusion. The electron power absorption can be primarily decoupled into the positive Ohmic power absorption in the bulk plasma region and the negative pressure-induced power absorption near the target surface. Ohmic power absorption is the dominant electron power absorption mechanism, mostly contributed by the azimuthal electron current. The power absorption due to electron inertial effects is negligible on time-average. Both the maximum power absorption and dissipation of electrons appear in the bulk plasma region during the second half of the RF period, implying a strong electron trapping in magnetron discharges. The contribution of secondary electrons is negligible under typical RFMS discharge conditions.

Physics↗

Growth and Saturation of the Electron Drift Instability in a Crossed Field Plasma

The linear growth and nonlinear energy transfer of the electron drift instability (EDI) are experimentally measured in the plume of a low-temperature, Hall effect discharge. A frequency-based bispectral analysis technique applied to fast ion density fluctuation measurements shows a growth rate function that is qualitatively similar to predictions from the linear instability dispersion relation, but an order of magnitude smaller. Calculation of the nonlinear transfer function indicates multiple three-wave interactions between high-frequency resonances of the instability in addition to an inverse energy cascade toward lower-frequency modes. Furthermore, these results are discussed in the context of recent theoretical, numerical, and experimental efforts on the EDI in Hall effect discharges and how the EDI may impact anomalous cross field transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reconciling calculations and measurements of inverse bremsstrahlung absorption

It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

XGC-6D

XGC-6D is a geometric full-f Particle-In-Cell code on unstructured meshes including XGC-like unstructured mesh, tetrahedral mesh, etc. The code is particularly suitable to study boundary plasma with steep gradient and frequency higher than ion gyro-frequency. References: Wang Z, Qin H, Sturdevant B, Chang CS. Geometric electrostatic particle-in-cell algorithm on unstructured meshes. Journal of Plasma Physics. 2021;87(4):905870406.

Wang, Zhenyu [Princeton Plasma Physics Laboratory ↗

XGC-6D

XGC-6D is a geometric full-f Particle-In-Cell code on unstructured meshes including XGC-like unstructured mesh, tetrahedral mesh, etc. The code is particularly suitable to study boundary plasma with steep gradient and frequency higher than ion gyro-frequency. References: Wang Z, Qin H, Sturdevant B, Chang CS. Geometric electrostatic particle-in-cell algorithm on unstructured meshes. Journal of Plasma Physics. 2021;87(4):905870406.

Wang, Zhenyu [Princeton Plasma Physics Laboratory ↗

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electromagnetic formalism

This is a companion paper to the previous work [P. H. Yoon, Phys. Plasmas 31, 032309 (2024)] in which the nonlinear susceptibilities of weakly turbulent magnetized plasma are derived under a simplifying assumption of electrostatic interaction. The present paper extends the analysis to a general situation of electromagnetic interaction. The main novelty of the previous and present papers is that by employing the Bessel function addition theorem, the mathematical definitions for the susceptibilities are substantially simplified, a procedure that has not been discussed in the existing literature. In the present paper, a full set of Maxwell’s equations are considered in conjunction with the nonlinear Vlasov equation, which is solved by a perturbative method. The result is a fully general nonlinear susceptibility, given in tensorial form, which is applicable for weakly turbulent magnetized plasmas.

Materials Science↗

Effects of Coulomb coupling on friction in strongly magnetized plasmas

The friction force on a test particle traveling through a plasma that is both strongly coupled and strongly magnetized is studied using molecular dynamics simulations. In addition to the usual stopping power component aligned antiparallel to the velocity, a transverse component that is perpendicular to both the velocity and Lorentz force is observed. This component, which was recently discovered in weakly coupled plasmas, is found to increase in both absolute and relative magnitude in the strongly coupled regime. Strong coupling is also observed to induce a third component of the friction force in the direction of the Lorentz force. These first-principles simulations reveal novel physics associated with collisions in strongly coupled, strongly magnetized, plasmas that are not predicted by existing kinetic theories. The effect is expected to influence macroscopic transport in a number of laboratory experiments and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shivalik Plasma Device-I, a glow discharge device to study the collective dynamics of dusty plasma

Using a uniquely configured glow discharge-based Shivalik Plasma Device-I, we demonstrate a variety of collective phenomena in dusty plasma away from the glow discharge region. The cylindrical glass device produces plasma using parallel disc-shaped electrodes with a smaller anode size than the cathode. The dust microparticles are initially sprinkled over the grounded cathode. These particles acquire a significant negative charge upon plasma formation, resulting in their levitation due to the balance between the Coulomb force and gravity. The new device supports the levitation of a big-sized (10 × 8 × 5 cm 3 ) three-dimensional dust cloud over the glass surface. It contrasts the dusty plasma formations in-between electrodes reported earlier. As the discharge voltage varies from high to low, the dust cloud travels from over the glass surface to between the electrodes. A complex interplay of dust void over the cathode, a sharp density gradient, and gravity lead to self-excitation of collective dust phenomena. It includes dust density waves (phase velocity, v ph ~4 cm/s), dust cloud oscillation (frequency, f = 5 Hz), sheared flow (flow velocity, v f ~ 1 cm/s), and multiple-sized dust vortices. These dust vortices provided an excellent platform for studying turbulent mixing phenomena. The power spectrum analysis agreed with two-dimensional Kolmogorov power-law scaling. This is an ideal dusty plasma apparatus where we can create or move the dust cloud to a location of choice from the glass surface to in-between the electrodes and excite one among many collective dust dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗