Search NASA⌕ Search

SEARCH · Search NASA

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 181 records · Page 10

Collaborative Research: Enabling multi-scale studies of magnetic reconnection with interpretable data-driven models

The development of accurate reduced descriptions and improved closures for magnetic reconnection is an important and a long‐standing challenge in plasma physics. The four‐fluid approach, and associated closures, that were investigated have the potential to improve the accuracy of plasma fluid models, capturing physical effects which would otherwise require a kinetic description. If successful, this approach could have an important impact for the modeling of laboratory and space plasmas. The major goals of this project were to develop new machine learning (ML) tools based on sparse and symbolic regression techniques, and to extract interpretable and generalizable reduced models (e.g., in the form of partial differential equations - PDEs) from data generated by first principles plasma simulations. Preserving interpretability of such data‐driven models is key to addressing the long‐standing theoretical and numerical challenges. Prior proof‐of‐principle studies have demonstrated the enormous potential of this approach, by recovering the well‐established hierarchy of plasma equations (from Vlasov to MHD) from data produced by particle‐in‐cell (PIC) simulations. Our goal in this project was to extend and apply these new tools to construct better kinetic closures for magnetic reconnection; to derive better models of particle injection and acceleration by this fundamental plasma process; and to use this understanding to accelerate the development of multi‐scale plasma algorithms. While our immediate focus was on the problem of magnetic reconnection, the tools that were will developed are general and applicable to other areas of plasma physics, and more broadly to many‐body phenomena. We anticipate that the development of these multi‐scale models will have a significant impact across different areas of plasma science, from fusion to space and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The dynamics of low-beta plasma clouds as simulated by a three-dimensional, electromagnetic particle code

A 3D electromagnetic and relativistic particle simulation code is used to investigate the dynamics of low-beta plasma clouds moving perpendicular to an ambient magnetic field in vacuum and in a background plasma. The simulations show the formation of the space charge sheaths at the sides of the cloud with the associated polarization electric field which facilitate the cross-field propagation, as well as the sheaths at the front and rear end of the cloud caused by the larger ion Larmor radius, which allows ions to move ahead and lag behind the electrons as they gyrate. It is found that, in a background plasma, electron and ion sheaths expand along the magnetic field at the same rate, whereas in vacuum the electron sheath expands much faster than the ion sheath. Sheath clouds are accelerated up to relativistic energies, which indicates that artificial plasma clouds released in the ionosphere or magnetosphere may generate optical emissions as energetic sheath electrons scatter in the upper atmosphere.

Neubert, T.↗

Injection of an overdense electron beam in space

A three-dimensional particle simulation model is used to study the injection and propagation of an overdense electron beam in the vicinity of a conductor in space. Beam electrons with a density of more than 100 times the ambient electron density are modeled using large-scale plasma simulations; in these simulations the surface area of the conductor is several thousands times that of the beam cross section at the injection point. The parameters of the simulations are chosen to allow the realistic simulation of active space experiments such as the Space Shuttle's Spacelab 2 electron beam mission. These simulations confirm space observations that an overdense electron beam can at least partially escape the near vicinity of the spacecraft, even in a fully ionized plasma. Once they have escaped from the vicinity of the spacecraft, these beam electrons should be able to propagate away freely until their energy is dissipated. It is suggested that such large-scale simulations could be useful in interpreting data from space experiments.

Okuda, Hideo↗

Informing ARC divertor design and plasma facing material selection through integrated modeling

This INFUSE 2023 project between UCLA and Commonwealth Fusion Systems used computer modeling to test whether tungsten materials can survive in CFS's ARC fusion reactor divertor. The team simulated plasma conditions and material responses, finding that tungsten-rhenium alloys resist grain growth better than pure tungsten, and that hydrogen buildup depends more on particle flux than temperature. The work helps CFS design durable plasma-facing components for their fusion power plant.

36 MATERIALS SCIENCE↗

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗

Plasma physics abstracts, 1 January - 31 December, 1986

Topics addressed include: ion-cyclotron waves; plasma waves; solar wind lithium releases; bow shock; Pi2 wave bursts; auroral kilometric radiation; ion energization; magnetic field corrections; electric fields; magnetospheric processes; electron acceleration; inner heliosphere; nightside auroral zone; computerized simulation; plasma wave turbulence; and magnetohydrodynamic waves in plasma sheets.

Gurnett, D. A.↗

Experimental investigations of a uranium plasma pertinent to a self-sustaining plasma source

The research is pertinent to the realization of a self-sustained fissioning plasma for applications such as nuclear propulsion, closed cycle MHD power generation using a plasma core reactor, and heat engines such as the nuclear piston engine, as well as the direct conversion of fission energy into optical radiation (nuclear pumped lasers). Diagnostic measurement methods and experimental devices simulating plasma core reactor conditions are discussed. Studies on the following topics are considered: (1) ballistic piston compressor (U-235); (2) high pressure uranium plasma (natural uranium); (3) sliding spark discharge (natural uranium); (4) fission fragment interaction (He-3 and U-235); and (5) nuclear pumped lasers (He-3 and U-235).

Schneider, R. T.↗

A Structure-Preserving Decorated Particle Method for the Vlasov-Poisson System

We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by macro-particles with position and velocity. In contrast, Scovel-Weinstein decorated particles involve additional shape degrees of freedom, while maintaining a finite-dimensional reduction with Hamiltonian structure inherited from the continuum model. Although the original work established this structure three decades ago, its computational potential has remained largely unexplored. We present a practical implementation of the Scovel-Weinstein model and compare it with a standard PIC algorithm. Numerical experiments demonstrate that macro-particles in standard PIC can be replaced by far fewer decorated particles while retaining comparable accuracy. This decorated particle approach offers a new structure-preserving paradigm for kinetic plasma simulation.

65M75, 70H05, 70G65↗

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gyrophase-bunched Electrons: Cluster Observations and Simulations

The interaction of the solar wind with the Earth's bow shock is responsible for a number of phenomena both upstream and downstream. One of the least understood and studied of these is the generation of gyrophase-bunched electrons. We will describe initial work searching for and analyzing such events in (Plasma Electron And Current Experiment (PEACE) data from CLUSTER. Plasma simulations will complement the data analysis. Recent analyses suggest that gyrophase-bunched electrons are present rather frequently in the upstream region, due either to energization in the reflection off the shock front or, more likely, to phase trapping in locally produced whistler wave fields.

Gurgiolo, C.↗

Modeling of Sheath Ion-Molecule Reactions in Plasma Enhanced Chemical Vapor Deposition of Carbon Nanotubes

In many plasma simulations, ion-molecule reactions are modeled using ion energy independent reaction rate coefficients that are taken from low temperature selected-ion flow tube experiments. Only exothermic or nearly thermoneutral reactions are considered. This is appropriate for plasma applications such as high-density plasma sources in which sheaths are collisionless and ion temperatures 111 the bulk p!asma do not deviate significantly from the gas temperature. However, for applications at high pressure and large sheath voltages, this assumption does not hold as the sheaths are collisional and ions gain significant energy in the sheaths from Joule heating. Ion temperatures and thus reaction rates vary significantly across the discharge, and endothermic reactions become important in the sheaths. One such application is plasma enhanced chemical vapor deposition of carbon nanotubes in which dc discharges are struck at pressures between 1-20 Torr with applied voltages in the range of 500-700 V. The present work investigates The importance of the inclusion of ion energy dependent ion-molecule reaction rates and the role of collision induced dissociation in generating radicals from the feedstock used in carbon nanotube growth.

Hash, David B.↗

Numerical simulation of turbulent flows with a three-dimensional vortex-in-cell method

A three-dimensional vortex-in-cell method has been developed for the evaluation of local flow fields due to a family of vortex filaments which employs the principles and architecture of a code developed for magnetic field evaluation in plasma simulations. The computational effort in the new method, as compared to 'vortex pushing' by direct Biot-Savarat interaction, increases directly, rather than quadratically, with the number of vortex elements. The method is well suited for studying large number of vortex filaments or rings and can be used for simulating continuous vorticity.

Buneman, O.↗

Modeling and verification of dynamic field ionization for laser-target interactions

Integrating field ionization models into kinetic plasma simulations is required for a variety of applications, especially when field strengths vary from low to high regimes, such as in laser-target interactions. The introduction of new physics models into kinetic codes requires a rigorous verification of their accuracy through well-defined verification problems. In this work, the field ionization model that has been included in the kinetic plasma code VPIC is presented, along with the detailed approach adopted for its integration. This model includes a comprehensive range of field ionization mechanisms: multiphoton ionization, tunneling ionization, and barrier suppression ionization. New verification problems employed to evaluate the ionization model's fidelity are outlined, and the simulation parameters that affect the accuracy of simulation results are explored. Additionally, this work addresses the impact of field ionization on computational performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Description of FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

The urgency to deliver fusion power is growing now more than ever, with increasing pressure for both public programs and private companies to meet milestones timelines and overcome significant remaining technical challenges to ensure growth of a nascent fusion industry in time to meet rapidly growing clean energy demands. With incredible advancements in computation and years of investment in fusion model development and validation, integrated modeling is poised to fill a key role in accelerating the timeline to a fusion pilot plant (FPP). Future fusion pilot plants will operate in regimes far beyond current experience, and device design will rely on physics-based prediction and extrapolation. Many concepts will also rely on simulation to assess safety (shielding, tritium management, materials activation and lifetimes), economics and scalability before the decision to build. Importantly, integrated simulation can be used to reveal and solve the complexities of system integration that may otherwise not be apparent in physical components or models developed in isolation. New experimental test facilities that produce relevant conditions to validate and resolve key technical challenges for various subsystems (materials, blankets, fuel cycle, etc.) have been repeatedly called for by the fusion community but are not yet realized. Integrated modeling has an important role in identifying realistic load conditions (thermal, electromagnetic, plasma, neutron and photon loads, etc.) and defining the components and experiments for these test facilities in order to ensure meaningful validation that sufficiently reduces modeling uncertainties and technical risk for the full integrated reactor. The Fusion REactor Design and Assessment (FREDA) SciDAC project is building a component-based integrated modeling framework & data structure to enable self-consistent, multi-fidelity, iterative optimization workflows for the fusion reactor design process. FREDA aims to shorten the time to viable designs by providing a set of flexible workflows to support the various stages of the design process using an integrated model hierarchy, ranging from the simple analytic descriptions to the highest fidelity, theory-based plasma and engineering modeling developed by the fusion and fission communities. These tools are expected to be needed for timely support of FPP design in the milestone program and in the FIRE collaboratives. The plasma simulation backbone of FREDA is IPS-FASTRAN with newly developed coupled Core-Edge Pedestal-SOL (CESOL) workflows, which is being extended to the far-SOL region up to the plasma facing components. FREDA incorporates the FERMI engineering modeling suite and will enable self-consistent evaluation of the thermal shields, limiters, blanket, magnets, and other surrounding structures with predictions of temperatures, erosion, dpa, activation, tritium generation and transport, creep, corrosion, material degradation, etc. Parametric generation of 3D CAD enables rapid iteration of component geometry in response to plasma and loading specifications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation studies of plasma waves in the electron foreshock - The transition from reactive to kinetic instability

Particle simulation experiments were used to analyze the electron beam-plasma instability. It is shown that there is a transition from the reactive state of the electron beam-plasma instability to the kinetic instability of Langmuir waves. Quantitative tests, which include an evaluation of the dispersion relation for the evolving non-Maxwellian beam distribution, show that a quasi-linear theory describes the onset of this transition and applies again fully to the kinetic stage. This stage is practically identical to the late stage seen in simulations of plasma waves in the electron foreshock described by Dum (1990).

Dum, C. T.↗

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗