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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Threshold criterion for a space simulation beam-plasma discharge

An experimental and theoretical study of the threshold characteristics of a space simulation beam-plasma discharge with emphasis on density profiles and a density-dependent ignition criterion. The study included various beam-plasma conditions covering beam currents from 8 to 85 mA, beam energies from 0.8 to 2.0 keV, and magnetic fields at 0.9 and 1.5 G. The study included experimental determinations of radial profiles of electron density for each of the selected conditions extending from a low-density, pre-beam-plasma discharge state to a strong beam-plasma discharge condition. The experimental results are shown to agree with detailed model calculations, which consider the beam-plasma discharge to be produced by large-amplitude electron plasma waves resulting from the beam-plasma interaction.

Szuszczewicz, E. P.↗
NASA NTRS2010

Hybrid-PIC Computer Simulation of the Plasma and Erosion Processes in Hall Thrusters

HPHall software simulates and tracks the time-dependent evolution of the plasma and erosion processes in the discharge chamber and near-field plume of Hall thrusters. HPHall is an axisymmetric solver that employs a hybrid fluid/particle-in-cell (Hybrid-PIC) numerical approach. HPHall, originally developed by MIT in 1998, was upgraded to HPHall-2 by the Polytechnic University of Madrid in 2006. The Jet Propulsion Laboratory has continued the development of HPHall-2 through upgrades to the physical models employed in the code, and the addition of entirely new ones. Primary among these are the inclusion of a three-region electron mobility model that more accurately depicts the cross-field electron transport, and the development of an erosion sub-model that allows for the tracking of the erosion of the discharge chamber wall. The code is being developed to provide NASA science missions with a predictive tool of Hall thruster performance and lifetime that can be used to validate Hall thrusters for missions.

Hofer, Richard R.↗
DOE OSTI2024

Simulating energetic ions and enhanced fusion rates from ion-cyclotron resonance heating with a full-wave/Fokker–Planck model

Reproducing fast-ion enhanced fusion rates from ion-cyclotron resonance heating (ICRH) in tokamaks requires the self-consistent coupling of a full-wave solver and a Fokker–Planck solver, which evolves multiple simultaneously resonant ion species. We introduce a new self-consistent model that iterates the TORIC full-wave solver with the CQL3D Fokker–Planck solver using the integrated plasma simulator (IPS). This model evolves the bounce-averaged ion distribution functions in both parallel and perpendicular velocity-space with a quasilinear radio frequency (RF) diffusion operator valid in the ion finite Larmor radius (FLR) limit and the RF electric fields with the resultant non-Maxwellian FLR dielectric tensor. This produces non-Maxwellian ICRH simulations that are fully self-consistent, fast, and interoperable with integrated modeling frameworks, such as TRANSP/GACODE/IPS-FASTRAN. We demonstrate our model's capabilities by validating it against experimental data in Alcator C-Mod. We then perform the first RF heating simulations of SPARC using self-consistent non-Maxwellian ion distributions to investigate the potential to enhance fusion rates using ion cyclotron resonance heating generated fast ions.

Physics↗
DOE OSTI2025

WEST L-mode record long pulses guided by predictions using Integrated Modeling

A new record was set on the WEST Tokamak, designed to operate long duration plasmas in a tungsten (W) environment, with an injected energy of 1.15 GJ and a plasma duration 364s. Scenario development was supported by integrated modeling using the High Fidelity Plasma Simulator (HFPS), the European IMAS-coupled version of JETTO/JINTRAC, which integrates physics-driven modules into a unified framework. In particular, a reduced model for Lower-Hybrid heating and Current-Drive (LHCD) and the quasi-linear turbulent transport model TGLF are crucial for long pulses predictions up to the Last Closed Flux Surface (LCFS). Using this workflow, a 100 s reference discharge was modeled and plasma kinetic profiles and loop voltage were quantitatively well matched. In preparation for the recent long duration experiments, non-inductive current-drive actuators (I P , n e , P LHCD ) were varied to determine the operational domain going towards fully non-inductive discharges. In particular, decreasing the plasma current is shown to ease the access to such conditions, with a careful monitoring of (n e , P LHCD ) to avoid machine limitations. In addition, post-prediction experiments conducted within the investigated parameter range validated the predicted dependencies and were shown to be in quantitative agreement. Exploratory work on the use of ECCD for MHD stability purpose is also introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Principles of magnetohydrodynamic simulation in space plasmas

Attention is given to the philosophical as well as physical principles that are essential to the establishment of MHD simulation studies for solar plasma research, assuming the capabilities of state-of-the-art computers and emphasizing the importance of 'local' MHD simulation. Solar-terrestrial plasma space is divided into several elementary regions where a macroscopic elementary energy conversion process could conceivably occur; the local MHD simulation is defined as self-contained in each of the regions. The importance of, and the difficulties associated with, the boundary condition are discussed in detail. The roles of diagnostics and of the finite difference method are noted.

Sato, T.↗
NASA NTRS2012

Laboratory Measurements of the Dielectronic Recombination Satellite Transitions of He-Like FE XXV and H-Like FE XXVI

We present laboratory spectra of dielectronic recombination (DR) satellite transitions attached to the He-like and H-like iron resonance lines obtained with the NASA Goddard Space Flight Center X-ray calorimeter and produced by a thermal plasma simu1ation technique on the EBIT-I electron beam ion trap at the Lawrence Livermore National Laboratory. We demonstrate that the calorimeter has sufficient spectral resolution in the 6-9 keV range to provide reliable measurements not only of standard DR satellite to resonance line intensities but also of DR satellite to DR satellite ratios that can be used to diagnose nonthermal electron distributions. Electron temperatures derived from the measured line intensities are consistent with the temperature of the simulated plasma. Temperature measurements based on DR satellite transitions have significant advantages over those based on collisional ionization equilibrium or continuum shape. Thus, successful demonstration of this method with the X-ray calorimeter is an important step fur its application in X-ray astronomy.

Gu, M. F.↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗
NASA NTRS1990

Particle simulations in magnetospheric plasmas

In view of the recent remarkable advancement of computer technology and simulation software, simulation studies are one of the most powerful academic tools for establishment of quantitative space physics and modelling of our space environment. The complex nature encountered in space plasma physics has motivated considerable development in computer simulations, which have played an essential role in the development of space plasma theory. This report describes research undertaken to understand physical processes involved in plasma waves observed in the magnetospheric plasmas, and associated nonlinear phenomena such as heating, diffusion, and acceleration of particles due to excited waves. The research explains and clarifies the observational data both qualitatively and quantitatively.

Nishikawa, Ken-Ichi↗
DOE OSTI2024

Encoding of linear kinetic plasma problems in quantum circuits via data compression

We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation Aψ = b to be solved by using the quantum signal processing algorithm. The latter requires encoding of matrix A in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode A in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗
DOE OSTI2025

Minimizing CGYRO HPC Communication Costs in Ensembles with XGYRO by Sharing the Collisional Constant Tensor Structure

First-principles fusion plasma simulations are both compute and memory intensive, and CGYRO is no exception. The use of many HPC nodes to fit the problem in the available memory thus results in significant communication overhead, which is hard to avoid for any single simulation. That said, most fusion studies are composed of ensembles of simulations, so we developed a new tool, named XGYRO, that executes a whole ensemble of CGYRO simulations as a single HPC job. By treating the ensemble as a unit, XGYRO can alter the global buffer distribution logic and apply optimizations that are not feasible on any single simulation, but only on the ensemble as a whole. The main saving comes from the sharing of the collisional constant tensor structure, since its values are typically identical between parameter-sweep simulations. This data structure dominates the memory consumption of CGYRO simulations, so distributing it among the whole ensemble results in drastic memory savings for each simulation, which in turn results in overall lower communication overhead.

CGYRO↗
NASA NTRS1973

Development of a closed pore insulation material

A closed pore ceramic foam insulation material (CPI) has been developed that offers possibilities for use as a reusable external heat shield for the NASA manned space shuttle. The outstanding characteristics of CPI are: (1) negligible water absorption due to a noninterconnecting network of cells; (2) high emittance at room and elevated temperature; (3) ability to survive at least 10 simulated reentry cycles to 1500 K using radiant heat lamps to simulate the reentry heat fluxes; (4) ability to survive, with no change in properties or appearance, at least 10 simulated plasma arc jet cycles to 1500 K (with the exception of some stress cracks induced either by the unduly severe nature of the initial arc splash heating pulse or by improper mechanical holding of the specimen in the test fixture); (5) strength (flexure); and (6) a low thermal conductivity throughout the temperature range of interest for the space shuttle.

Tobin, A.↗
NASA NTRS1983

Experimental modelling of satellite wakes in auroral arcs

Preliminary measurements have been made in a large laboratory discharge device configured to simulate plasma wake phenomena. A large relative electron drift impinging on a nonconducting disc gives rise to a perturbed density struture downstream, as well as a reflected beam of electrons upstream. Velocity distribution measurements were made using a novel energy analyzer with angular resolution.

Wild, N.↗
NASA NTRS1984

Collisionless dissipation in quasi-perpendicular shocks

Microscopic dissipation processes in quasi-perpendicular shocks are studied by two-dimensional plasma simulations in which electrons and ions are treated as particles moving in self-consistent electric and magnetic fields. Cross-field currents induce substantial turbulence at the shock front reducing the reflected ion fraction, increasing the bulk ion temperature behind the shock, doubling the average magnetic ramp thickness, and enhancing the upstream field aligned electron heat flow. The short scale length magnetic fluctuations observed in the bow shock are probably associated with this turbulence.

Forslund, D. W.↗
NASA NTRS1987

Plasma contactors for electrodynamic tethers

Plasma contactors could be used to ground satellites to space plasma to acquire a flow of electrons to propel or power the satellites. A tether would cut across geomagnetic field lines, producing a potential difference between the ends of the tether. Closing the connection between the ends would form a circuit in which an electrical load could be inserted. Design constraints of the circuit are low impedance and a fully reversible high current. The contactor would generate a neutral plasma to connect to the ionospheric plasma. The surface area of the connection would have to be kept large enough for the current density to be equal to the random electron current density in the unperturbed space plasmas. The other contactor would feed electrons and draw ions from the space plasma. Experimental results from spaceborne and ground-based space plasma simulator tests of hollow cathodes that have shown that multiampere currents can be collected are described.

Patterson, Michael J.↗