Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA SIMULATION”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Multi-channel transport experiments at Alcator C-Mod and comparison with gyrokinetic simulations

Multi-channel transport experiments have been conducted in auxiliary heated (Ion Cyclotron Range of Frequencies) L-mode plasmas at Alcator C-Mod [Marmar and Alcator C-Mod Group, Fusion Sci. Technol. 51(3), 3261 (2007)]. These plasmas provide good diagnostic coverage for measurements of kinetic profiles, impurity transport, and turbulence (electron temperature and density fluctuations). In the experiments, a steady sawtoothing L-mode plasma with 1.2 MW of on-axis RF heating is established and density is scanned by 20%. Measured rotation profiles change from peaked to hollow in shape as density is increased, but electron density and impurity profiles remain peaked. Ion or electron heat fluxes from the two plasmas are the same. The experimental results are compared directly to nonlinear gyrokinetic theory using synthetic diagnostics and the code GYRO [Candy and Waltz, J. Comput. Phys. 186, 545 (2003)]. We find good agreement with experimental ion heat flux, impurity particle transport, and trends in the fluctuation level ratio (T(tilde sign){sub e}/T{sub e})/(ñ{sub e}/n{sub e}), but underprediction of electron heat flux. We find that changes in momentum transport (rotation profiles changing from peaked to hollow) do not correlate with changes in particle transport, and also do not correlate with changes in linear mode dominance, e.g., Ion Temperature Gradient versus Trapped Electron Mode. The new C-Mod results suggest that the drives for momentum transport differ from drives for heat and particle transport. The experimental results are inconsistent with present quasilinear models, and the strong sensitivity of core rotation to density remains unexplained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Hollow cathode-based plasma contactor experiments for electrodynamic tether

The role plasma contactors play in effective electrodynamic tether operation is discussed. Hollow cathodes and hollow cathode-based plasma sources have been identified as leading candidates for the electrodynamic tether plasma contactor. Present experimental efforts to evaluate the suitability of these devices as plasma contactors are reviewed. This research includes the definition of preliminary plasma contactor designs, and the characterization of their operation as electron collectors from a simulated space plasma. The discovery of an 'ignited mode' regime of high contactor efficiency and low impedance is discussed, as well as is the application of recent models of the plasma coupling process to contactor operation. Results indicate that ampere-level electron currents can be exchanged between hollow cathode-based plasma contactors and a dilute plasma in this regime. A discussion of design considerations for plasma contactors is given which includes expressions defining the total mass flow rate and power requirements of plasma contactors operating in both the cathodic and anodic regimes, and correlation of this to the tether current. Finally, future ground and spaceflight experiments are proposed to resolve critical issues of plasma contactor operation.

Patterson, Michael J.↗
DOE OSTI2015

A one-dimensional study of the evolution of the microwave breakdown in air

The microwave breakdown in air is simulated numerically within a simple 1D model taking into account a perturbation of electromagnetic field by plasma. The simulations were performed using two qualitatively different codes. One of these codes is based on computation of Maxwell equations, whereas the other one utilizes an approximation of quasi-monochromatic electromagnetic field. There is a good agreement between simulation results obtained by using both codes. Calculations have been carried out in a wide range of air pressures and field frequencies; also varied were initial spatial distributions of plasma density. The results reveal strong dependence of the breakdown evolution on the relation between the field frequency and the gas pressure as well as on the presence of extended rarefied background plasma. At relatively low gas pressures (or high field frequencies), the breakdown process is accompanied by the stationary ionization wave propagating towards the incident electromagnetic wave. In the case of a high gas pressure (or a relatively low field frequency), the peculiarities of the breakdown are associated with a formation of plasma filament array. The extended background plasma can suppress formation of the plasma filament array completely even at high pressures (or low frequencies)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Energy-dependent implementation of secondary electron emission models in continuum kinetic sheath simulations

The plasma-material interactions present in multiple fusion and propulsion concepts between the flow of plasma through a channel and a material wall drive the emission of secondary electrons. This emission is capable of altering the fundamental structure of the sheath region, significantly changing the expected particle fluxes to the wall. The emission spectrum is separated into two major energy regimes, a peak of elastically backscattered primary electrons at the incoming energy, and cold secondary electrons inelastically emitted directly from the material. The ability of continuum kinetic simulations to accurately represent the secondary electron emission is limited by relevant models being formulated in terms of monoenergetic particle interactions which cannot be applied directly to the discrete distribution function. As a result, rigorous implementation of energy-dependent physics is often neglected in favor of simplified, constant models. We present here a novel implementation of semi-empirical models in the boundary of continuum kinetic simulations which allows the full range of this emission to be accurately captured in physically-relevant regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Experimental and simulation study of target biasing effects on plasma transport in linear plasma device MPS-LD

Linear plasma devices (LPDs) are important experimental platforms for investigating plasma–material interactions (PMI). In PMI experiments, it has been found that applying a target bias not only effectively modifies the incident ion energy, but also induces significant changes in the electron density and electron temperature, whereby the evolution of these plasma parameters is primarily governed by plasma transport processes. However, at present, the physical process and mechanism underlying such bias-induced variations remain unclear. In this work, biasing experiments under argon plasma discharge conditions were first carried out on the MPS-LD device. For the corresponding experiments, an electric potential model was newly developed based on the BOUT++ LPD module, enabling self-consistent simulations of plasma transport under biased conditions. Numerical simulations were then performed to reproduce the experimental results and to validate the accuracy of the proposed model. Finally, by combining experimental measurements with numerical simulations, a bias-voltage scan was performed to investigate how the electron density and electron temperature vary with the bias voltage (U bias ). The results show that applying negative bias decreases the target electron density (n e,T ) while increasing the target electron temperature (T e,T ). In contrast, positive bias increases both n e,T and T e,T ; however, at high positive bias, n e,T first reaches a maximum and subsequently decreases with further increases in U bias . The underlying physical mechanisms are analyzed using particle flux, momentum, and energy conservation. It indicates that the applied bias regulates the parallel electric field, thereby changing ion and electron velocities, and consequently affecting the electron density. At high positive bias, the ion velocity is further influenced by ion viscosity, leading to the reversal in n e,T . Meanwhile, the enhanced parallel electric field drives stronger currents, significantly increasing ion–electron frictional work and converting the input bias power into electron energy, which raises the electron temperature. In conclusion, these results contribute to a deeper understanding of the effects and mechanisms of biasing on plasma transport in the MPS-LD device.

BOUT++ simulation↗
DOE OSTI2025

Anisotropic Heating and Parallel Heat Flux in Electron-only Magnetic Reconnection with Intense Guide Fields

Electron-only reconnection (E-REC) is a process recently observed in the Earth’s magnetosheath, where magnetic reconnection occurs at electron kinetic scales, and ions do not couple to the reconnection process. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales in some settings. This paper investigates E-REC under different intensities of strong guide fields (the ratio between the guide field and the in-plane asymptotic field strength is 5, 10 and 20, respectively) via two-dimensional fully kinetic particle-in-cell simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. Similarly to previous experimental studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a hexapolar (six-polar) to an octopolar (eight-polar) structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and finite parallel electron heat flux density is observed. The three simulations exhibit features of the Chew–Goldberger–Low theory, with the level of consistency increasing as the guide field strength increases. This explains the electron parallel temperature variations and the shape of the EVDFs observed along the separatrices. This work may provide insights into the understanding of electron heating and parallel heat flux density in E-REC observed in the turbulent magnetosheath.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Simulation of the nonlinear evolution of electron plasma waves

Electrostatic waves driven by an electron beam in an ambient magnetized plasma were studied using a quasi-1D PIC simulation of electron plasma waves (i.e., Langmuir waves). The results disclose the presence of a process for moving wave energy from frequencies and wavenumbers predicted by linear theory to the Langmuir-like frequencies during saturation of the instability. A decay process for producing backward propagating Langmuir-like waves, along with low-frequency waves, is observed. The simulation results, however, indicate that the backscattering process is not the conventional Langmuir wave decay. Electrostatic waves near multiples of the electron plasma frequency are generated by wave-wave coupling during the nonlinear stage of the simulations, confirming the suggestion of Klimas (1983).

Nishikawa, K.-I.↗
NASA NTRS1981

The simulation of the geosynchronous Earth orbit plasma environment in Chamber A: An assessment of possible experimental investigations

The possible use of Chamber A for the replication or simulation of space plasma physics processes which occur in the geosynchronous Earth orbit (GEO) environment is considered. It is shown that replication is not possible and that scaling of the environmental conditions is required for study of the important instability processes. Rules for such experimental scaling are given. At the present time, it does not appear technologically feasible to satisfy these requirements in Chamber A. It is, however, possible to study and qualitatively evaluate the problem of vehicle charging at GEO. In particular, Chamber A is sufficiently large that a complete operational spacecraft could be irradiated by beams and charged to high potentials. Such testing would contribute to the assessment of the operational malfunctions expected at GEO and their possible correction. However, because of the many tabulated limitations in such a testing programs, its direct relevance to conditions expected in the geo environment remains questionable.

Bernstein, W.↗
DOE OSTI2025

Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing

The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

3-D Simulations of NSTAR Ion Thruster Plasma Interactions

Described is a Particle-in-Cell with Monte Carlo Collision code developed to perform detailed three-dimensional ion thruster simulations. To capture the full kinetic behavior of ion thruster plumes, both the electrons and ions are treated as test particles. Simulation results are given of the NSTAR ion thruster under ground test and in space conditions. Numerical results are compared.

NSTAR Ion Thruster Plasma Ionization Particles↗
DOE OSTI2021

Numerical heating in particle-in-cell simulations with Monte Carlo binary collisions

The binary Monte Carlo (MC) collision algorithm is a standard and robust method to include binary Coulomb collision effects in particle-in-cell (PIC) simulations of plasmas. Here we show that the coupling between PIC and MC algorithms can give rise to (nonphysical) numerical heating of the system that significantly exceeds that observed when these algorithms operate independently. We argue that this deleterious effect results from an inconsistency between the particle motion associated with MC collisions and the work performed by the collective electromagnetic field on the PIC grid. This inconsistency manifests as the (artificial) stochastic production of electromagnetic energy, which ultimately heats the plasma particles. Here, the MC-induced numerical heating can significantly impact the evolution of the simulated system for long simulation times (≳10 3 collision periods, for typical numerical parameters). We describe the source of the MC-induced numerical heating analytically and discuss strategies to minimize it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus

A particle-in-cell algorithm is derived with a canonical Poisson structure in the formalism of finite element exterior calculus. The resulting method belongs to the class of gauge-compatible splitting algorithms, which exactly preserve gauge symmetries and their associated conservation laws via the momentum map. We numerically demonstrate this time invariance of the momentum map and its usefulness in establishing precise initial conditions with a desired initial electric field and fixed background charge. The restriction of this canonical, finite element Poisson structure to the 1X2P $1\frac {1}{2}$ -dimensional phase space is also considered and simulated numerically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Multi-region approach to free-boundary three-dimensional tokamak equilibria and resistive wall instabilities

Free-boundary 3D tokamak equilibria and resistive wall instabilities are calculated using a new resistive wall model in the two-fluid M3D-C1 code. In this model, the resistive wall and surrounding vacuum region are included within the computational domain. This implementation contrasts with the method typically used in fluid codes in which the resistive wall is treated as a boundary condition on the computational domain boundary and has the advantage of maintaining purely local coupling of mesh elements. This new capability is used to simulate perturbed, free-boundary non-axisymmetric equilibria; the linear evolution of resistive wall modes; and the linear and nonlinear evolution of axisymmetric vertical displacement events (VDEs). Calculated growth rates for a resistive wall mode with arbitrary wall thickness are shown to agree well with the analytic theory. Equilibrium and VDE calculations are performed in diverted tokamak geometry, at physically realistic values of dissipation, and with resistive walls of finite width. Simulations of a VDE disruption extend into the current-quench phase, in which the plasma becomes limited by the first wall, and strong currents are observed to flow in the wall, in the SOL, and from the plasma to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Conceptual Design and Analysis of In-Vessel Components for the Materials Plasma Exposure eXperiment (MPEX)

The materials plasma exposure experiment (MPEX) is a linear plasma divertor simulator currently undergoing conceptual design. The facility will expose material samples to steady-state plasma fluxes to examine plasma–material interactions (PMIs) that are expected in the next generation of fusion devices. The plasmas will be generated by a helicon source, with electron and ion heating sources of up to 800 kW possible. The peak heat fluxes of the target are expected to be up to 10 MW/m 2 . The facility will be capable of handling low-activation neutron-irradiated samples in order to examine the multivariate effects of neutron damage and plasma fluence. Neutron-irradiated samples are planned to be roughly of 10-mm diameter; however, plasma-facing components up to 60×600 mm can be accommodated. The steady-state nature of the device will require the magnetic confinement of the plasma to be achieved with superconducting magnets, with a maximum on-axis field of 2.5 T. In addition, since MPEX will be a steady-state device, in-vessel components need to be water cooled. The primary in-vessel components will be the target, the dump plate, the limiter, the skimmers, and the microwave absorber. The conceptual design of these components is presented here, including analyses that confirm that the designs are adequate to meet the requirements of MPEX operation.

Lumsdaine, Arnold↗
NASA NTRS2010

Magnetic Shielding of the Acceleration Channel Walls in a Long-Life Hall Thruster

In a Qualification Life Test (QLT) of the BPT-4000 Hall thruster that recently accumulated greater than 10,000 h it was found that the erosion of the acceleration channel practically stopped after approximately 5,600 h. Numerical simulations of this thruster using a 2-D axisymmetric, magnetic field-aligned-mesh (MFAM) plasma solver reveal that the process that led to this significant reduction of the erosion was multifaceted. It is found that when the channel receded from its early-in-life geometry to its steady-state configuration several changes in the near-wall plasma and sheath were induced by the magnetic field that, collectively, constituted an effective shielding of the walls from any significant ion bombardment. Because all such changes in the behavior of the ionized gas near the eroding surfaces were caused by the topology of the magnetic field there, we term this process "magnetic shielding."

life modeling↗
DOE OSTI2015

Propagation of localized structures in relativistic magnetized electron-positron plasmas using particle-in-cell simulations

We use a particle-in-cell simulation to study the propagation of localized structures in a magnetized electron-positron plasma with relativistic finite temperature. We use as initial condition for the simulation an envelope soliton solution of the nonlinear Schrödinger equation, derived from the relativistic two fluid equations in the strongly magnetized limit. This envelope soliton turns out not to be a stable solution for the simulation and splits in two localized structures propagating in opposite directions. However, these two localized structures exhibit a soliton-like behavior, as they keep their profile after they collide with each other due to the periodic boundary conditions. We also observe the formation of localized structures in the evolution of a spatially uniform circularly polarized Alfvén wave. In both cases, the localized structures propagate with an amplitude independent velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Isolation and phase-space energization analysis of the instabilities in collisionless shocks

We analyse the generation of kinetic instabilities and their effect on the energization of ions in non-relativistic, oblique collisionless shocks using a 3D-3V (three spatial with three velocity components) simulation by dHybridR, a hybrid particle-in-cell code. At sufficiently high Mach number, quasi-perpendicular and oblique shocks can experience rippling of the shock surface caused by kinetic instabilities arising from free energy in the ion velocity distribution due to the combination of the incoming ion beam and the population of ions reflected at the shock front. To understand the role of the ripple on particle energization, we devise a new instability isolation method to identify the unstable modes underlying the ripple and interpret the results in terms of the governing kinetic instability. We generate velocity-space signatures using the field–particle correlation technique to look at energy transfer in phase space from the isolated instability driving the shock ripple, providing a viewpoint on the different dynamics of distinct populations of ions in phase space. Together, the field–particle correlation technique and our new instability isolation method provide a unique viewpoint on the different dynamics of distinct populations of ions in phase space and allow us to completely characterize the energetics of the collisionless shock under investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

New wave effects in nonstationary plasma

Through particle-in-cell simulations and analytics, a host of interesting and novel wave effects in nonstationary plasma are examined. In particular, Langmuir waves serve as a model system to explore wave dynamics in plasmas undergoing compression, expansion, and charge recombination. The entire wave life-cycle is explored, including wave excitation, adiabatic evolution and action conservation, nonadiabatic evolution and resonant wave-particle effects, collisional dissipation, and potential laboratory applications of the aforementioned phenomenology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗