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.

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

NASA NTRS1997

Numerical Simulation of Rotation-Driven Plasma Transport In the Jovian Magnetosphere

A Jupiter version of the Rice Convection Model (RCM-J) was developed with support of an earlier NASA SR&T grant. The conversion from Earth to Jupiter included adding currents driven by centrifugal force, reversing the planetary magnetic field, and rescaling various parameters. A series of informative runs was carried out, all of them solving initial value problems. The simulations followed an initial plasma torus configuration as it fell apart by interchange instability. Some conclusions from the simulations were the following: 1. We confirmed that, for conventional values of the torus density and ionospheric conductance, the torus disintegrates by interchange instability on a time scale of approx. one day, which is 1-2 orders of magnitude shorter than the best estimates of the average residence time of plasma in the torus. 2. In the model, the instability could be slowed to an arbitrary degree by the addition of sufficient impounding energetic particles, as suggested earlier by Siscoe et al (1981). However, the observed energetic particles do not seem sufficient to guarantee impoundment (e.g., Mauk et al., 1996). 3. Whether inhibited by impoundment or not, the interchange was found to proceed by the formation of long fingers, which get thinner as they get longer. This picture differed dramatically from the conventional radial-diffusion picture (e.g., Siscoe and Summers (1981)), more superficially with the outward-moving-blob picture (Pontius and Hill, 1989). The obvious limitation of the original RCM-J was that it could not represent a plasma source. We could represent the decay of a pre-existing torus, but we could not represent the way ionization of material from Io continually replenishes the plasma. We consequently were precluded from studying a whole set of fundamental issues of torus theory, including whether the system can come to a steady state.

Wolf, Richard A.↗

NASA NTRS1993

Plasma dynamics driven by finite-width current filament and KV potential drops in ionosphere-magnetosphere coupling

We have simulated the plasma dynamics of a quasi-two-dimensional current filament in ionosphere-magnetosphere coupling. The simulation consists of a set of one-dimensional flux tube simulations with different imposed time-dependent, field-aligned currents. The dynamical interactions of the individual current filaments with the ionospheric and magnetospheric plasma generate kV field-aligned potential drops along the field lines, and the side-by-side display exhibits the evolution of the implied potential structure in the horizontal direction. The primary cross-field motion produced by the horizontal E-field (ExB drift) is perpendicular to both of the significant spatial directions, and is thus ignorable in this geometry. The effects of other cross-field drict processes are discussed. The simulation thus provides insight into the dynamical evolution of 2D potential structures driven by an imposed finite-width field-aligned current profile. Spatial and temporal variation of other plasma parameters (temperatures, density, etc.) that are determined primarily by parallel transport effects are also displayed.

Ganguli, Supriya B.↗

DOE OSTI2024

Effect of plasma initialization on 3D PIC simulation of Hall thruster azimuthal instability

Abstract The lack of understanding of the azimuthal instability and the resulting electron anomalous transport limits further improvement of Hall thrusters. Compared to theoretical and experimental approaches, the numerical particle-in-cell (PIC) simulation is a suitable and powerful tool, which has been widely applied to investigate the azimuthal instability, and great progress has been made in the past decades. However, PIC simulations are intrinsically computationally expensive, and it is realized that the Hall thruster azimuthal instability has a three dimensional nature. Therefore, massive 3D PIC simulation must be carried out to completely reveal the mechanism of the instability. In this paper, the effect of plasma initialization on 3D PIC simulation of Hall thruster azimuthal instability is studied as a starting point. It is found that by initializing with ion density and velocity fitting functions to the steady-state simulation results, a faster convergence can be obtained and the computational time can be reduced by about 1.5 times. Typical fitting functions of ion density, drifting velocity, and temperature are given, and the influence of different initialization profiles is presented.

Physics↗

NASA NTRS1999

A Simulation to Study Speed Distributions in a Solar Plasma

We have carried out a numerical simulation of a plasma with characteristics similar to those found in the core of the Sun. Particular emphasis is placed on the Coulomb interaction between the ions and electrons, which could result in a relative velocity distribution different from the Maxwell-Boltzmann (MB) distribution generally assumed for a plasma. The fact that the distribution may not exactly follow the MB distribution could have very important consequences for a variety of problems in solar physics, especially the neutrino problem. Very briefly, the neutrino problem is that the observed neutrino detections from the Sun are smaller than what the standard solar theory predicts. In Section I we introduce the problem and in section II we discuss the approach to try to solve the problem: i.e., a molecular dynamics approach. In section III we provide details about the integration method, and any simplifications that can be applied to the problem. In section IV (the core of this report) we state our results. First for the specific case of 1000 particles and then for other cases with different number of particles. In section V we summarize our findings and state our conclusions. Sections VI VII and VIII provide the list of figures, reference material and acknowledgements respectively.

Cheeseman, Peter↗

NASA NTRS1965

Plasma chambers and scaling

Laboratory simulation and scaling of plasmas found in ionosphere, magnetosphere, and interplanetary space

MAGNETOSPHERE↗

NASA NTRS1990

Cold plasma heating in the plasma sheet boundary layer - Theory and simulations

Satellite observations in recent years have confirmed that the plasma sheet boundary layer is a permanent feature of the earth's magnetotail located between the lobe and central plasma sheet during both quiet and active magnetic periods. Distinct features of the boundary layer include field aligned ion beams and intense electrostatic emissions known as broadband electrostatic noise. Since the plasma sheet boundary layer is a spatial feature of the magnetotail, within it will occur thermal mixing of the resident warm boundary layer plasma with inflowing (convecting) cold ionospheric plasma. A theoretical study involving linear theory and nonlinear numerical particle simulations is presented which examines ion beam instabilities in the presence of a thermally mixed hot and cold background plasma. It is found that the free energy in the ion beams can heat the cool ionospheric plasma to ambient plasma sheet boundary layer temperatures via broadband electrostatic noise. These results, along with recent observational reports that ionospheric outflow can account for measured plasma sheet densities, suggest that the ionospheric role in plasma sheet dynamics and content may be as large as the solar wind.

Schriver, David↗

DOE OSTI2025

Visualization techniques for the gyrokinetic tokamak simulation code

Gyrokinetic simulations of plasma microturbulence in tokamaks are challenging to visualize because the compute grid follows the magnetic field lines that spiral around the torus. We have overcome this challenge by developing three new approaches that improve visualization of gyrokinetics. Our techniques work directly with the topology of magnetic flux surfaces where the simulation stores variables in concentric rings on poloidal planes (vertical cross sections of the torus). Our visualization preview step triangulates each consecutive pair of rings to display the data on a poloidal plane. The second visualization technique follows spiral field lines around the torus and constructs polygons to visualize a flux surface. Third, the poloidal triangles are connected between planes to form prisms that compose a 3-D model of the entire torus. The visualization workflow produces detailed geometry that matches the high resolution, irregular compute grid for every time step. The surface and solid models are displayed in scientific visualization programs to effectively explore and communicate the results, including fluctuation of electron density, ion temperature, and electrostatic potential. Highly detailed renderings verify plasma behavior along magnetic field lines over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Understanding extraction limits of plasma cathodes with experiment and simulation

The project is focused on computations for enhanced ionization near the exit orifice of a plasma cathode and how the cathode extraction electrode geometry and stray magnetic field structure influence available extractable current. The computational effort will take place at the Princeton Collaborative Research Laboratory (PCRF). The specific modeling platform to be utilized from PCRF is the electrostatic direct implicit particle in cell (EDIPIC) code. EDIPIC will be used to study ionization processes in the electron extraction region of the cathode. The control of electron energy is a key aspect to minimizing power losses in the plasma. Exiting electrons ionize gas leaving the orifice and in the process produce the dense plasma from which large amounts of current can be extracted. Project activities include studying processes that determine the energy distribution of the exiting electrons thus providing insight into how to optimize the source. The model will be validated and developed using experimental measurements as warranted. The familiarity and expertise of the team at PPPL with this type of plasma and the challenges posed greatly enhances the likelihood of gaining insight into the extraction sheath under the duration of the proposed activity. The effort is expected to shed a great deal of insight into electron extraction processes and pathways for optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Current collection through the transport of electrons across magnetic field lines

Numerical simulations are used to investigate the collection of electrons by a long conducting cylinder in a flowing plasma. The plasma flow simulates the relative motion between a spacecraft and the plasma. The sheath structures and the levels of electron current collections for the cases with and without an ambient magnetic field are examined. It is found that for the flow perpendicular to the magnetic field, the current is considerably enhanced depending on the relative drift velocity. In the case of a nonzero magnetic field perpendicular to the cylinder axis, the potential structure is a 2D double layer with dimensions L-perpendicular much less than L-parallel. L-perpendicular is found to be the current limiting radius given by the Parker-Murphy model (1967). For the flow along the ambient magnetic field, the electron current is found to be smaller than that for the flow perpendicular to this field. This is explained in terms of the potential structures.

Singh, N.↗

DOE OSTI2025

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

Transient effects in beam-plasma interactions in a space simulation chamber stimulated by a fast pulse electron gun

Interest in the interaction of electron beams with plasma generated by ionization caused by the primary electron beam was stimulated by the need to develop special vacuum tubes to operate in the kMHz frequency region. The experiments of Getty and Smullin (1963) indicated that the interaction of an energetic electron beam with its self-produced plasma resulted in the emission of wave energy over a wide range of frequencies associated with cyclotron and longitudinal plasma instabilities. This enhanced the thermal plasma density in the vicinity of the beam, and the term Beam-Plasma Discharge (BPD) was employed to described this phenomenon. The present investigation is concerned with some of the transient phenomena associated with wave emission during the beam switch-on and switch-off periods. Results are presented on the changes in electron energy spectra on a time scale of tens of milliseconds following beam switch-on. The results are discussed in terms of the beam plasma discharge phenomenon.

Raitt, W. J.↗

NASA NTRS2012

Three-Dimensional Numerical Simulations of Equatorial Spread F: Results and Observations in the Pacific Sector

A three-dimensional numerical simulation of plasma density irregularities in the postsunset equatorial F region ionosphere leading to equatorial spread F (ESF) is described. The simulation evolves under realistic background conditions including bottomside plasma shear flow and vertical current. It also incorporates C/NOFS satellite data which partially specify the forcing. A combination of generalized Rayleigh-Taylor instability (GRT) and collisional shear instability (CSI) produces growing waveforms with key features that agree with C/NOFS satellite and ALTAIR radar observations in the Pacific sector, including features such as gross morphology and rates of development. The transient response of CSI is consistent with the observation of bottomside waves with wavelengths close to 30 km, whereas the steady state behavior of the combined instability can account for the 100+ km wavelength waves that predominate in the F region.

ALTAIR↗

NASA NTRS1987

Core plasma in the magnetosphere

The paper presents a review of new findings related to core plasma (which includes ions and electrons with energies less than 50 eV) identified in studies of magnetosphere during 1983-1986. These findings include the cleft ion fountain; the plasmaspheric heavy ion torus; core molecular ions; high core-plasma densities in the plasma-sheet boundary layer; intense transverse heating of core plasmas in the equatorial regions; the supersonic polar wind (H/+/ and He/+/); toroidal or ring heavy-ion distribution functions in the auroral region; conic-to-field-aligned ion signatures of inverted-V auroral events; and the large-scale electron density structure of the magnetosphere. The progress made on modeling core plasma includes simulation of thermal helium heating at the equator; kinetic models of the plasma transport from localized ionospheric regions; time-dependent hydrodynamic models of the polar wind and plasmasphere-ionosphere coupling; and kinetic models of plasmasphere refilling.

Horwitz, J. L.↗

NASA NTRS1988

A three-dimensional MHD simulation of the interaction of the solar wind with Comet Halley

The interaction between the solar wind and cometary plasmas is simulated using a three-dimensional time-dependent MHD simulation model, and the results are compared with the recent satellite observations of Comet Halley. The model, which includes cometary mass loading, reproduces many of the features observed by the Suisei probe and the Giottot, including the weak bow shock, the enhancement of the magnetic field in front of the contact surface, and the plasma temperature increase across the bow shock (while it decreased near the comet).

Ogino, Tatsuki↗

DOE OSTI2024

Theoretical and global simulation analysis of collisional microtearing modes

Microtearing modes (MTMs) are suggested as a candidate for anomalous thermal transport in tokamak H-mode discharges. This study investigates MTMs in tokamak plasmas, employing simulations in the BOUT++ framework. It simplifies and linearizes the governing equations in detailed linear simulations. The study meticulously evaluates various conductivity models under diverse plasma conditions and collision regimes. The research thoroughly assesses different conductivity models across a range of plasma conditions and collision regimes. A unified dispersion relation that includes both MTM and Drift-Alfvén Wave (DAW) instabilities is derived, showing that DAW and MTM instabilities occur at varying distances from the rational surface. Specifically, MTMs become unstable near the rational surface but stabilize farther away, while drift-Alfvén instability appears farther from the rational surface. The study also re-derives MTM dispersion relations using Ohm's law and the vorticity equation, providing a thorough analysis of electromagnetic and electrostatic interactions in tokamaks. Global simulations demonstrate an inverse correlation between MTM growth rates and collisionality, and a direct correlation with temperature gradients. The nonalignment of the rational surface with the peak ω*e stabilizes the MTMs. Nonlinear simulations highlight electron temperature relaxation as the primary saturation mechanism for MTMs, with magnetic flutter identified as the dominant mode of electron thermal transport.

Fan, K. (ORCID:0000000227518809)↗

DOE OSTI2025

Data from "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions"

Data and input files related to the paper "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions" (https://doi.org/10.1116/6.0004027). This includes the final DP model used in all simulations, training data set, example input files to run DeepMD (with LAMMPS), and data tables summarizing the results obtained from the simulations.

machine learning models↗

DOE OSTI2025

Confinement performance predictions for a high field axisymmetric tandem mirror

This paper presents a Hammir tandem mirror confinement performance analysis based on Realta Fusion’s first-of-a-kind model for axisymmetric magnetic mirror fusion performance. This model uses an integrated end plug simulation model including, heating, equilibrium and transport combined with a new formulation of the plasma operation contours (POPCONs) technique for the tandem mirror central cell. Using this model in concert with machine learning optimization techniques, it is shown that an end plug utilizing high temperature superconducting magnets and modern neutral beams enables a classical tandem mirror pilot plant producing a fusion gain Q > 5. The approach here represents an important advance in tandem mirror design. The high-fidelity end plug model enables calculations of heating and transport in the highly non-Maxwellian end plug to be made more accurately. The detailed end plug modelling performed in this work has highlighted the importance of classical radial transport and neutral beam absorption efficiency on end plug viability. The central cell POPCON technique allows consideration of a wide range of parameters in the relatively simple near-Maxwellian central cell, facilitating the selection of more optimal central cell plasmas. These advances make it possible to find more conservative classical tandem mirror fusion pilot plant operating points with lower temperatures, neutral beam energies and end plug performance requirements than designs in the literature. Despite being more conservative, it is shown that these operating points have sufficient confinement performance to serve as the basis of a viable fusion pilot plant provided that they can be stabilized against magnetohydrodynamic and trapped particle modes.

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.↗