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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2022

3D full wave fast wave modeling with realistic HHFW antenna geometry and SOL plasma in NSTX-U

This paper reports the significant advancement of our ability to model and to understand how RF waves interact with the SOL plasma, by developing for the full torus 3D scrape-off layer (SOL) plasma simulation together with the antenna and core plasma in the NSTX-U device in HHFW frequency regime for a cold plasma model. The present paper extends the previous modeling (Bertelli et al 2020 AIP Conf. Proc. 2254 030001) to a full 3D device geometry including realistic antenna geometry in order to capture a more realistic experimental situation including the fully 3D wave field effects and the antenna plasma interaction in the SOL plasma, and, at the same time, resulting core wave propagation. Here, the central tool of the present work is the Petra-M code, which is a newly developed state-of-the-art generic electromagnetic simulation tool for modeling RF wave propagation based on MFEM (http://mfem.org), an open source scalable C++ finite element method library.

3D RF modeling↗

NASA NTRS1991

Procedures for Including Secondary Electron Emission in Numerical Simulations of Plasma-Insulator Interactions

Previous Monte Carlo simulations provide a data base for properties of secondary electron emission (SEE) from insulators and metals. Incident primary electrons are considered at energies up to 1200 eV. The behavior of secondary electrons is characterized by (1) yield vs. primary energy E(sub p), (2) distribution vs. secondary energy E(sub s), and (3) distribution vs. angle of emission theta. Special attention is paid to the low energy range E(sub p) up to 50 eV, where the number and energy of secondary electrons is limited by the finite band gap of the insulator. For primary energies above 50 eV the SEE yield curve can be conveniently parameterized by a Haffner formula. The energy distribution of secondary electrons is described by an empirical formula with average energy about 8.0 eV. The angular distribution of secondaries is slightly more peaked in the forward direction than the customary cos theta distribution. Empirical formulas and parameters are given for all yield and distribution curves. Procedures and algorithms are described for using these results to find the SEE yield, and then to choose the energy and angle of emergence of each secondary electron. These procedures can readily be incorporated into numerical simulations of plasma-solid surface interactions in low earth orbit.

Beyst, Brian↗

DOE OSTI2016

Effect of a pulsating electric field on ECR heating in the CERA-RX(C) X-ray generator

3D particle-in-cell plasma simulations for the field configurations implemented in the CERA-RX(C) ECR X-ray generator (2.45 GHz) have been conducted. Dependences of the energy spectra of electrons incident on the target electrode on the amplitude and frequency of pulsations of the electric field in a megahertz range are derived. The simulation data are compared with the results of the full-scale experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Neutral Vlasov kinetic theory of magnetized plasmas

The low-frequency limit of Maxwell equations is considered in the Maxwell-Vlasov system. This limit produces a neutral Vlasov system that captures essential features of plasma dynamics, while neglecting radiation effects. Euler-Poincaré reduction theory is used to show that the neutral Vlasov kinetic theory possesses a variational formulation in both Lagrangian and Eulerian coordinates. By construction, the new model recovers all collisionless neutral models employed in plasma simulations. Then, comparisons between the neutral Vlasov system and hybrid kinetic-fluid models are presented in the linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Verification of particle simulation of radio frequency waves in fusion plasmas

Radio frequency (RF) waves can offer heating, current and flow drive, as well as instability control for steady state operations of fusion experiments. A particle simulation model has been designed in this work to provide a first-principles tool for studying the RF nonlinear interactions with plasmas. In this model, ions are considered as fully kinetic particles using the Vlasov equation and electrons are treated as guiding centers using the drift kinetic equation. This model has been implemented in a global gyrokinetic toroidal code using real electron-to-ion mass ratio. To confirm the model, linear simulations of ion plasma oscillation, ion Bernstein wave, and lower hybrid wave are carried out in cylindrical geometry and found to agree well with analytic predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP

Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

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↗

DOE OSTI2016

A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma

In order to enable kinetic simulation of non-thermal edge plasmas at a reduced computational cost, a new hybrid-Lagrangian δf scheme has been developed that utilizes the phase space grid in addition to the usual marker particles, taking advantage of the computational strengths from both sides. The new scheme splits the particle distribution function of a kinetic equation into two parts. Marker particles contain the fast space-time varying, δf, part of the distribution function and the coarse-grained phase-space grid contains the slow space-time varying part. The coarse-grained phase-space grid reduces the memory-requirement and the computing cost, while the marker particles provide scalable computing ability for the fine-grained physics. Weights of the marker particles are determined by a direct weight evolution equation instead of the differential form weight evolution equations that the conventional delta-f schemes use. The particle weight can be slowly transferred to the phase space grid, thereby reducing the growth of the particle weights. The non-Lagrangian part of the kinetic equation – e.g., collision operation, ionization, charge exchange, heat-source, radiative cooling, and others – can be operated directly on the phase space grid. Deviation of the particle distribution function on the velocity grid from a Maxwellian distribution function – driven by ionization, charge exchange and wall loss – is allowed to be arbitrarily large. The numerical scheme is implemented in the gyrokinetic particle code XGC1, which specializes in simulating the tokamak edge plasma that crosses the magnetic separatrix and is in contact with the material wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Non-thermal electron acceleration in low Mach number collisionless shocks. II. Firehose-mediated Fermi acceleration and its dependence on pre-shock conditions

Electron acceleration to non-thermal energies is known to occur in low Mach number (M{sub s} ≲ 5) shocks in galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Using two-dimensional (2D) particle-in-cell (PIC) plasma simulations, we showed in Paper I that electrons are efficiently accelerated in low Mach number (M{sub s} = 3) quasi-perpendicular shocks via a Fermi-like process. The electrons bounce between the upstream region and the shock front, with each reflection at the shock resulting in energy gain via shock drift acceleration. The upstream scattering is provided by oblique magnetic waves that are self-generated by the electrons escaping ahead of the shock. In the present work, we employ additional 2D PIC simulations to address the nature of the upstream oblique waves. We find that the waves are generated by the shock-reflected electrons via the firehose instability, which is driven by an anisotropy in the electron velocity distribution. We systematically explore how the efficiency of wave generation and of electron acceleration depend on the magnetic field obliquity, the flow magnetization (or equivalently, the plasma beta), and the upstream electron temperature. We find that the mechanism works for shocks with high plasma beta (≳ 20) at nearly all magnetic field obliquities, and for electron temperatures in the range relevant for galaxy clusters. Our findings offer a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1994

Auroral plasma transport processes in the presence of kV potential structures

We have simulated plasma transport processes in the presence of a quasi-two-dimensional current filament, that generated kV potential structure in the auroral region. The simulation consists of a set of one-dimensional flux tube simulations with different imposed time-dependent, field-aligned currents. The model uses the 16 moment system of equations and simultaneously solves coupled continuity and momentum equations and equations describing the transport along the magnetic field lines of parallel and perpendicular thermal energy and heat flows for each species. The lower end of the simulation is at an altitude of 800 km, in the collisional topside ionosphere, while the upper end is at 10 R(sub E) in the magnetosphere. The plasma consists of hot electrons and protons of magnetospheric origin and low-energy electrons, protons, and oxygen ions of ionospheric origin. The dynamical interaction of the individual current filaments with ionospheric and magnetospheric plasma generates a potential structure in the horizontal direction and kilovolt field-aligned potential drops along the field lines. The side-by-side display exhibits the evolution of the implied potential structure in the horizontial direction. In the presence of this potential structure and parallel electric field ionospheric plasma density is depleted and velocity is reduced, while density enhancement and increased velocity is observed in magnetospheric plasma. The ionospheric and magnetospheric electron temperatures increase below 2 R(sub E) due to magnetic mirror force on converging geomagnetic field lines. The primary cross-field motion produced by the horizontal E field (E x B drift) is perpendicular to both of the significant spatial directions and is thus ignorable in this geometry. The effects of other cross-field drift processes are discussed. The simulation thus provides insight into the dynamical evolution of two-dimensional potential structures driven by an imposed finite width, field-aligned current profile.

Ganguli, Supriya B.↗

NASA NTRS2002

Test and Analysis Capabilities of the Space Environment Effects Team at Marshall Space Flight Center

Marshall Space Flight Center has developed world-class space environmental effects testing facilities to simulate the space environment. The combined environmental effects test system exposes temperature-controlled samples to simultaneous protons, high- and low-energy electrons, vacuum ultraviolet (VUV) radiation, and near-ultraviolet (NUV) radiation. Separate chambers for studying the effects of NUV and VUV at elevated temperatures are also available. The Atomic Oxygen Beam Facility exposes samples to atomic oxygen of 5 eV energy to simulate low-Earth orbit (LEO). The LEO space plasma simulators are used to study current collection to biased spacecraft surfaces, arcing from insulators and electrical conductivity of materials. Plasma propulsion techniques are analyzed using the Marshall magnetic mirror system. The micro light gas gun simulates micrometeoroid and space debris impacts. Candidate materials and hardware for spacecraft can be evaluated for durability in the space environment with a variety of analytical techniques. Mass, solar absorptance, infrared emittance, transmission, reflectance, bidirectional reflectance distribution function, and surface morphology characterization can be performed. The data from the space environmental effects testing facilities, combined with analytical results from flight experiments, enable the Environmental Effects Group to determine optimum materials for use on spacecraft.

Finckenor, M. M.↗

NASA NTRS1987

Simulations of auroral plasma processes - Electric fields, waves and particles

Plasma processes driven by current sheets of finite thicknesses in an ambient magnetized plasma are studied using a 2 1/2 dimensional particle-in-cell code, and similarities are found between simulated plasma processes and those observed in the auroral plasma. Current sheets are shown to be bounded by large perpendicular electric fields occurring near their edges above the conducting boundary. Shaped potential structures form when the current sheets are narrow, and when the current sheets are wide, potential structures develop a significant parallel potential drop such that the electrons are accelerated upwards. Downward parallel electric fields of variable strength are noted in the downward current region, and double layer formation is seen in both narrow and wide current sheets. High frequency oscillations near the electron plasma frequency and its harmonic are seen, and low frequency waves are observed.

Singh, Nagendra↗

DOE OSTI2016

THE IBEX RIBBON AND THE PICKUP ION RING STABILITY IN THE OUTER HELIOSHEATH. I. THEORY AND HYBRID SIMULATIONS

The nearly circular band of energetic neutral atom emission dominating the field of view of the Interplanetary Boundary Explorer ( IBEX ) satellite, is most commonly attributed to the effect of charge exchange of secondary pickup ions (PUIs) gyrating about the magnetic field in the outer heliosheath and the interstellar space beyond. Several models for the PUI dynamics of this mechanism have been proposed, each requiring either strong or weak scattering of the initial pitch angle. Conventional wisdom states that ring distributions tend to generate waves and scatter onto a shell on timescales too short for charge exchange to occur. We performed a careful study of ring and thin shell proton distribution stability using theoretical tools and hybrid plasma simulations. We show that the kinetic behavior of a freshly injected proton ring is a far more complicated process than previously thought. In the presence of a warm Maxwellian core, narrower rings could be more stable than broader toroidal distributions. The scattered rings possess a fine structure that can only be revealed using very large numbers of macroparticles in a simulation. It is demonstrated that a “stability gap” in ring temperature exists where the protons could retain large gyrating anisotropies for years, and the wave activity could remain below the level of the ambient magnetic fluctuations in interstellar space. In the directions away from the ribbon, however, a partial shell distribution is more likely to be unstable, leading to significant scattering into one hemisphere in velocity space. The process is accompanied by turbulence production, which is puzzling given the very low level of magnetic fluctuations measured in the outer heliosheath by Voyager 1 .

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1982

Plasma flow measurements in a simulated low earth orbit plasma

The employment of large, higher power solar arrays for space operation has been considered, taking into account a utilization of high operating voltages. In connection with the consideration of such arrays, attention must be given to the fact that the ambient environment of space contains a tenuous low energy plasma which can interact with the high voltage array causing power 'leakage' and arcing. An investigation has been conducted with the aim to simulate the behavior of such an array in low-earth-orbit (LEO). During the experiments, local concentrations of the 'leakage' current were observed when the panel was at a high voltage. These concentrations could overload or damage a small area of cells in a large string. It was hypothesized that this effect was produced by electrostatic focusing of the particles by the sheath fields. To verify this experimentally, an end-effect Langmuir probe was employed. The obtained results are discussed.

Gabriel, S. B.↗

DOE OSTI2018

FOCUS: A full-orbit CUDA solver for particle simulations in magnetized plasmas

A Full-Orbit CUDA Solver for particle simulations in plasmas, FOCUS, has been developed. The code follows exact particle trajectories by solving Newton’s equations with the Lorentz force. The code can use fields calculated analytically or numerically and it is also possible to use the information provided by equilibrium reconstruction and transport codes. FOCUS has an elastic collisions module which covers the whole particle energy range in magnetic fusion devices. Moreover, an atomic collision module was also included to simulate the interaction of the test particles with neutral or partially ionized species. Regarding performance, the main feature is that the code runs on GPUs, allowing the simulation of a large number of particles using moderate computational resources. Different versions of the code have been used in several papers. Here we present a complete description of its capabilities, the basis physics included and detailed information about the numerical algorithms employed. Finally, some examples are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Progress with the COGENT Edge Kinetic Code: Implementing the Fokker-Plank Collision Operator

Here, COGENT is a continuum gyrokinetic code for edge plasma simulations being developed by the Edge Simulation Laboratory collaboration. The code is distinguished by application of a fourth-order finite-volume (conservative) discretization, and mapped multiblock grid technology to handle the geometric complexity of the tokamak edge. The distribution function F is discretized in v∥ – μ (parallel velocity – magnetic moment) velocity coordinates, and the code presently solves an axisymmetric full-f gyro-kinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. COGENT capabilities are extended by implementing the fully nonlinear Fokker-Plank operator to model Coulomb collisions in magnetized edge plasmas. The corresponding Rosenbluth potentials are computed by making use of a finite-difference scheme and multipole-expansion boundary conditions. Details of the numerical algorithms and results of the initial verification studies are discussed. (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Full-wave feasibility study of magnetic diagnostic based on O-X mode conversion and oblique reflectometry imaging

An innovative millimeter wave diagnostic is proposed to measure the local magnetic field and the edge current as a function of the minor radius in the pedestal region. The idea behind such diagnostic is to localize and characterize a direction of reduced reflectivity at the O-mode cutoff layer. We modeled the wave scattering and mode-conversion processes by means of the finite-element COMSOL Multiphysics code in two dimensions (2D). Sensitivity studies were performed for parameters mocking up DIII-D plasmas. Simulations confirmed the presence of a minimum in reflectivity of an externally injected O-mode beam, and confirmed that this minimum depends on the magnetic field at the cutoff, as expected from the OX mode conversion physics. This study gives confidence in the feasibility of the diagnostic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Preliminary design of a tangentially viewing imaging bolometer for NSTX-U

The infrared imaging video bolometer (IRVB) measures plasma radiated power images using a thin metal foil. Two different designs with a tangential view of NSTX-U are made assuming a 640 × 480 (1280 × 1024) pixel, 30 (105) fps, 50 (20) mK, IR camera imaging the 9 cm × 9 cm × 2 μm Pt foil. The foil is divided into 40 × 40 (64 × 64) IRVB channels. This gives a spatial resolution of 3.4 (2.2) cm on the machine mid-plane. The noise equivalent power density of the IRVB is given as 113 (46) μW/cm{sup 2} for a time resolution of 33 (20) ms. Synthetic images derived from Scrape Off Layer Plasma Simulation data using the IRVB geometry show peak signal levels ranging from ∼0.8 to ∼80 (∼0.36 to ∼26) mW/cm{sup 2}.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗