Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma modeling”

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 487 records · Page 27

AI-Enhanced Computational Tools for Entry Systems Modeling

To advance the understanding of complex atmospheric entry phenomena, NASA’s Entry Systems Modeling (ESM) team [1] has developed high-fidelity computational tools addressing multiscale challenges, from material microstructures to full-scale heatshield response. This abstract highlights a subset of ESM tools, focusing on AI integration to enhance workflows and predictive modeling. - PuMA [2] computes effective material properties from high-resolution micro-CT scans, supporting TPS analysis for NASA missions. - TomoSAM [3] automates 3D tomography dataset segmentation for PuMA using the Segment Anything Model, reducing manual effort and improving accuracy. - PATO [4] models porous reactive materials under extreme conditions, with advancements such as unified solvers, mechanical erosion, and TPS coatings for NASA missions. - arcjetCV [5] employs deep learning to analyze arc jet test footage, measuring recession rates, shape changes, and shock standoff distances, bridging simulations, and experiments to reveal TPS ablation behavior. - ARCHeS [6] simulates arc heater plasma flows, modeling turbulence, radiation, and electromagnetic interactions to optimize arc heater performance, validate TPS under extreme conditions, and serve as a foundation for developing digital twins of arc heater facilities. - SPARTA [7] simulates rarefied hypersonic flows and gas-surface interactions for planetary entry missions, leveraging GPU architectures for scalable and efficient aerothermal and ablation analyses. AI-driven solutions, such as deep learning segmentation, have streamlined workflows in ESM tools and still hold significant potential to further accelerate processes and enhance automation in entry systems modeling. [1] Haskins, J.B. (2023), [2] Ferguson, J.C. (2018), [3] Meurisse, J.B.E. (2018), [4] Semeraro, F. (2023), [5] Quintart, A. (2024) [6] Meurisse, J.B.E. (2022), [7] Plimpton, S.J. (2019)

Predictive Modeling↗

Modeling Growth of Nanostructures in Plasmas

As semiconductor circuits shrink to CDs below 0.1 nm, it is becoming increasingly critical to replace and/or enhance existing technology with nanoscale structures, such as nanowires for interconnects. Nanowires grown in plasmas are strongly dependent on processing conditions, such as gas composition and substrate temperature. Growth occurs at specific sites, or step-edges, with the bulk growth rate of the nanowires determined from the equation of motion of the nucleating crystalline steps. Traditional front-tracking algorithms, such as string-based or level set methods, suffer either from numerical complications in higher spatial dimensions, or from difficulties in incorporating surface-intense physical and chemical phenomena. Phase field models have the robustness of the level set method, combined with the ability to implement surface-specific chemistry that is required to model crystal growth, although they do not necessarily directly solve for the advancing front location. We have adopted a phase field approach and will present results of the adatom density and step-growth location in time as a function of processing conditions, such as temperature and plasma gas composition.

Hwang, Helen H.↗

Multi-cathode unbalanced magnetron sputtering systems

Ion bombardment of a growing film during deposition is necessary in many instances to ensure a fully dense coating, particularly for hard coatings. Until the recent advent of unbalanced magnetron (UBM) cathodes, reactive sputtering had not been able to achieve the same degree of ion bombardment as other physical vapor deposition processes. The amount of ion bombardment of the substrate depends on the plasma density at the substrate, and in a UBM system the amount of bombardment will depend on the degree of unbalance of the cathode. In multi-cathode systems, the magnetic fields between the cathodes must be linked to confine the fast electrons that collide with the gas atoms. Any break in this linkage results in electrons being lost and a low plasma density. Modeling of the magnetic fields in a UBM cathode using a finite element analysis program has provided great insight into the interaction between the magnetic fields in multi-cathode systems. Large multi-cathode systems will require very strong magnets or many cathodes in order to maintain the magnetic field strength needed to achieve a high plasma density. Electromagnets offer the possibility of independent control of the plasma density. Such a system would be a large-scale version of an ion beam enhanced deposition (IBED) system, but, for the UBM system where the plasma would completely surround the substrate, the acronym IBED might now stand for Ion Blanket Enhanced Deposition.

Sproul, William D.↗

Plasma-gas interactions studies in a hybrid plume plasma rocket

Plasma-gas interaction was investigated and the basic mechanisms for energy and particle transport. The solution approach assumes cylindrical geometry and includes a multiplicity of atomic reactions, and the presence of a strong magnetic field is described. The principal reactions are electron and ion impact ionization, and charge exchange between hot ions and cold neutrals. Radial particle and energy transport is mainly by diffusion. A modified Bohm diffusion model for plasma in the core of the plume, and classical neutral particle diffusion in the cooler regions of the flow are presented. Neutrals are allowed to free stream in the low density regions, where the collision mean-free-path becomes comparable, or even larger than the characteristic dimensions of the system.

Chang, F. R.↗

Calculation of the surface tension of liquid metals using a one-component-plasma reference system

The one-component-plasma (OCP) model is used as a reference system instead of the traditional hard-sphere fluid to calculate the liquid-vapor interfacial surface tension of liquid metals within the density functional formalism. The calculated surface tensions of the alkali metals are in excellent agreement with experiment. For the polyvalent metal Al, the result obtained is larger than experimental measurements. It is concluded that the OCP system is not suitable to describe the liquid-vapor phase transition in simple metals which have a nominal plasma parameter larger than the usual freezing value of about 178. The calculated interfacial widths in all cases are narrower than the expected experimental values.

Zeng, X. C.↗

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗

Temporal evolution of whistler growth in a cold plasma injection experiment

The evolution of whistler mode turbulence and particle participation in a cold plasma release are studied experimentally as part of the Active Magnetospheric Particle Tracers Explorer (AMPTE) program. The investigation used a simple time dependent cold plasma density model, and an atypical ambient radiation belt was assumed. It is shown that a cold lithium injection in the AMPTE parameter range can give rise to whistler mode turbulence with significant gain to the wave amplitude in a single pass through a flux tube. Whistler mode growth results in the pitch angle diffusion of energetic electrons, which are in turn precipitated. The rate of power input into the ionosphere from these precipitating electrons is roughly estimated at 2 ergs per sq cm/s. It is found that this value of the precipitated power is large enough to continue producing a visible aurora until the injected plasma in the flux tube is lost or destroyed by magnetospheric convection or other processes.

Ganguli, G.↗

Magnetohydrodynamic model of the Crab Nebula and its radiation

The interaction of the Crab pulsar with its surrounding nebula is modeled as a steady state, spherically symmetric relativistic wind which is confined by the outer layers of the supernova remnant. The pulsar wind is assumed to contain only positrons and electrons, and is terminated by a strong shock at a distance of 10 arcseconds from the pulsar, at the outer edge of the underluminous region. Down-stream of the shock, the flow decelerates and increases its pressure in order to match the boundary conditions imposed by outer nebula. Consistent flow solutions can be found if the magnetic luminosity of the pulsar wind is about 0.3 percent of its particle luminosity. The immediate postshock plasma is modeled as a power law distribution whose density and pressure moments satisfy the strong shock Rankine-Hugoniot relations. Using the downstream flow solution, the synchrotron continuum is calculated. The parameters of the pulsar wind and assumed power law spectral index are adjusted until the calculated total synchrotron power and the spectral distribution of the nebular continuum from optical to gamma-rays, are in reasonable accord with observations.

Coroniti, F. V.↗

Landau damping of magnetospherically reflected whistlers

Unducted VLF signals produced by lightning activity can form a population of magnetospherically reflected (MR) whistlers in the inner magnetosphere. It has been suggested recently that in the absence of significant attenuation such waves could merge into a broadband continuum with sufficient intensity to account for plasmaspheric hiss. To test this conjecture we have evaluated the path-integrated attenuation of MR whistlers along representative ray paths using the HOTRAY code. Using a realistic plasma distribution modeled on in-situ data, we find that the majority of MR waves experience significant damping after a few transits across the equator. This is primarily due to Landau resonance with suprathermal (0.1-1 keV) electrons. The attenuation is most pronounced for waves that propagate through the outer plasmasphere; this can readily account for the infrequent occurrence of multiple-hop MR waves for L greater than or equal to 3.5. Selected waves that originate at intermediate latitudes (15 deg is less than or equal to lambda is less than or equal to 35 deg) and whose ray paths are confined to the inner plasma- sphere may experience up to 10 magnetospheric reflections before substantial attentuation occurs. These waves should form the population of observed MR waves. Wave attenuation becomes more pronounced at higher frequencies; this can account for the absence of multiple-hop waves above 5 kHz. Weakly attenuated MR waves tend to migrate outward to the L shell, where their frequency is comparable to the equatorial lower hybrid frequency. The enhanced concentration of waves due to a merging of ray paths would produce a spectral feature that rises in frequency at lower L. This is quite distinct from the reported properties of plasmaspheric hiss, which maintains a constant frequency band throughout the entire plasmasphere. Furthermore, in the absence of mode conversion, waves below 500 Hz, which often form an important if not dominant part of the spectral properties of hiss, are unable to escape from the topside ionosphere in the whistler mode. Consequently, we conclude that unducted lightning signals cannot account for the origin of plasmaspheric hiss.

Thorne, Richard M.↗

First Results from ARTEMIS, a New Two-Spacecraft Lunar Mission: Counter-Streaming Plasma Populations in the Lunar Wake

We present observations from the first passage through the lunar plasma wake by one of two spacecraft comprising ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun), a new lunar mission that re-tasks two of five probes from the THEMIS magnetospheric mission. On Feb 13, 2010, ARTEMIS probe P1 passed through the wake at 3.5 lunar radii downstream from the Moon, in a region between those explored by Wind and the Lunar Prospector, Kaguya, Chandrayaan, and Chang'E missions. ARTEMIS observed interpenetrating proton, alpha particle, and electron populations refilling the wake along magnetic field lines from both flanks. The characteristics of these distributions match expectations from self-similar models of plasma expansion into vacuum, with an asymmetric character likely driven by a combination of a tilted interplanetary magnetic field and an anisotropic incident solar wind electron population. On this flyby, ARTEMIS provided unprecedented measurements of the interpenetrating beams of both electrons and ions naturally produced by the filtration and acceleration effects of electric fields set up during the refilling process. ARTEMIS also measured electrostatic oscillations closely correlated with counter-streaming electron beams in the wake, as previously hypothesized but never before directly measured. These observations demonstrate the capability of the comprehensively instrumented ARTEMIS spacecraft and the potential for new lunar science from this unique two spacecraft constellation.

Halekas, J. S.↗

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↗

Sunlit Io Atmospheric [O I] 6300 A and the Plasma Torus

A large database of sunlit Io [O I] 6300A emission, acquired over the period 1990-1999, with extensive coverage of Io orbital phase angle phi and System III longitude lambda(sub III), exhibits significant long-term and short-term variations in [O I] 6300A emission intensities. The long-term average intensity shows a clear dependence on lambda(sub III), which establishes conclusively that the emission is produced by the interaction between Io's atmosphere and the plasma torus. Two prominent average intensity maxima, 70 deg to 90 deg wide, are centered at lambda(sub III) about 130deg. and about 295 deg. A comparison of data from October 1998 with a three-dimensional plasma torus model, based upon electron impact excitation of atomic oxygen, suggests a basis for study of the torus interaction with Io's atmosphere. The observed short-term, erratic [O I] 6300A intensity variations fluctuate approximately 20% to 50% on time scale of tens of minutes with less frequent fluctuations of a factor of about 2. The most likely candidate to produce these fluctuations is a time-variable energy flux of field-aligned nonthermal electrons identified recently in Galileo PLS data. If true, the short-term [O I] intensity fluctuations may be related to variable field-aligned currents driven by inward and outward torus plasma transport and/or transient high-latitude, field-aligned potential drops. A correlation between the intensity and emission line width indicates molecular dissociation may contribute significantly to the [O I] 6300A emission. The nonthermal electron energy flux may produce O(1-D) by electron impact dissociation of SO2 and SO, with the excess energy going into excitation of O and its kinetic energy. The [O I] 6300A emission database establishes Io as a valuable probe of the torus, responding to local conditions at Io's position.

Oliversen, Ronald J.↗

First Results from ARTEMIS, A New Two-Spacecraft Lunar Mission: Counter-Streaming Plasma Populations in the Lunar Wake

We present observations from the first passage through the lunar plasma wake by one of two spacecraft comprising ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun), a new lunar mission that re-tasks two of five probes from the THEMIS magnetospheric mission. On Feb 13, 2010, ARTEMIS probe P1 passed through the wake at approximately 3.5 lunar radii downstream from the Moon, in a region between those explored by Wind and the Lunar Prospector, Kaguya, Chandrayaan, and Chang'E missions. ARTEMIS observed interpenetrating proton, alpha particle, and electron populations refilling the wake along magnetic field lines from both flanks. The characteristics of these distributions match expectations from self-similar models of plasma expansion into vacuum, with an asymmetric character likely driven by a combination of a tilted interplanetary magnetic field and an anisotropic incident solar wind electron population. On this flyby, ARTEMIS provided unprecedented measurements of the interpenetrating beams of both electrons and ions naturally produced by the filtration and acceleration effects of electric fields set up during the refilling process. ARTEMIS also measured electrostatic oscillations closely correlated with counter-streaming electron beams in the wake, as previously hypothesized but never before directly measured. These observations demonstrate the capability of the comprehensively instrumented ARTEMIS spacecraft and the potential for new lunar science from this unique two spacecraft constellation.

Plasma↗

Modeling of Solid Waste Processing Options in BIO-Plex

BIO-Plex is a ground-based test bed currently under development by NASA for testing technologies and practices that may be utilized in future long-term life support missions. All aspects of such an Advanced Life Support (ALS) System must be considered to confidently construct a reliable system, which will not only allow the crew to survive in harsh environments, but allow the crew time to perform meaningful research. Effective handling of solid wastes is a critical aspect of the system, especially when recovery of resources contained in the waste is required. This is particularly important for ALS Systems configurations that include a Biomass Production Chamber. In these cases, significant amounts of inedible biomass waste may be produced, which can ultimately serve as a repository of necessary resources for sustaining life, notably carbon, water, and plant nutrients. Numerous biological and physicochemical solid waste processing options have been considered. Biological options include composting, aerobic digestion, and anaerobic digestion. Physicochemical options include pyrolysis, SCWO (supercritical water oxidation), various incineration configurations, microwave incineration, magnetically assisted gasification, and low temperature plasma reaction. Modeling of these options is a necessary step to assist in the design process. A previously developed top-level model of BIO-Plex implemented in MATLAB Simulink (r) for the use of systems analysis and design has been adopted for this analysis. Presently, this model only considered incineration for solid waste processing. Present work, reported here, includes the expansion of this model to include a wider array of solid waste processing options selected from the above options, bearing in mind potential, near term solid waste treatment systems. Furthermore, a trade study has also been performed among these solid waste processing technologies in an effort to determine the ideal technology for long-term life support missions.

Rodriguez, Luis F.↗

Simulation of space plasma phenomena

The importance of computer modeling in investigating physical phenomena is discussed, noting that modern computers can follow the time evolution of systems containing many millions of degrees of freedom, all of which are simultaneously interacting with each other. The two types of models used in simulating plasma, that is, particle and fluid, are described. Attention is given to particle models and some recent results obtained concerning the diffuse aurora. The simulations are shown to give clear evidence that a type of nonlinear cyclotron resonance process for the cold particles is responsible for the saturation observed in these calculations. It is also responsible for the cold electron heating and it is found that at saturation the rate of heating of the cold electrons essentially balances the rate at which energy is fed into unstable waves.

Dawson, J. M.↗

Low altitude signature of the plasma sheet boundary layer: Observations and model

Low-altitude spacecraft on magnetotail field lines often detect a distinctive signature in the precipitating ion flux. A velocity-dispersed ion structure is often observed near the poleward boundary of the auroral oval. At the low-latitude edge of this structure an absence of precipitating ions is seen, previously referred to as 'the gap,' separating the velocity -- dispersed ions at the higher latitudes from the more diffuse, plasma sheet-like ions at lower latitudes. We present a model of low-altitude particle precipitation that reproduces these observed features in the ion spectra and provides a quantitative estimate of the downtail plasma sheet properties. The model calculations are compared with observations from the Akebono spacecraft. In this model, the dispersed ion velocity signature maps to a region in the distant plasma sheet where the plasma has a field-aligned bulk flow. The gap maps to a region in the distant magnetotail where the ion fluxes are below the detection threshold of the instrument, due to the low plasma sheet density and temperature in that region.

Onsager, T. G.↗

A nonequilibrium model for a moderate pressure hydrogen microwave discharge plasma

This document describes a simple nonequilibrium energy exchange and chemical reaction model to be used in a computational fluid dynamics calculation for a hydrogen plasma excited by microwaves. The model takes into account the exchange between the electrons and excited states of molecular and atomic hydrogen. Specifically, electron-translation, electron-vibration, translation-vibration, ionization, and dissociation are included. The model assumes three temperatures, translational/rotational, vibrational, and electron, each describing a Boltzmann distribution for its respective energy mode. The energy from the microwave source is coupled to the energy equation via a source term that depends on an effective electric field which must be calculated outside the present model. This electric field must be found by coupling the results of the fluid dynamics and kinetics solution with a solution to Maxwell's equations that includes the effects of the plasma permittivity. The solution to Maxwell's equations is not within the scope of this present paper.

Scott, Carl D.↗