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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Simulation Tool for Dielectric Barrier Discharge Plasma Actuators at Atmospheric and Sub-Atmospheric Pressures: SBIR Phase I Final Report

This report is the final report of a SBIR Phase I project. It is identical to the final report submitted, after some proprietary information of administrative nature has been removed. The development of a numerical simulation tool for dielectric barrier discharge (DBD) plasma actuator is reported. The objectives of the project were to analyze and predict DBD operation at wide range of ambient gas pressures. It overcomes the limitations of traditional DBD codes which are limited to low-speed applications and have weak prediction capabilities. The software tool allows DBD actuator analysis and prediction for subsonic to hypersonic flow regime. The simulation tool is based on the VORPAL code developed by Tech-X Corporation. VORPAL's capability of modeling DBD plasma actuator at low pressures (0.1 to 10 torr) using kinetic plasma modeling approach, and at moderate to atmospheric pressures (1 to 10 atm) using hydrodynamic plasma modeling approach, were demonstrated. In addition, results of experiments with pulsed+bias DBD configuration that were performed for validation purposes are reported.

Likhanskii, Alexandre↗
NASA NTRS1983

Experimental measurements of the plasma sheath around pinhole defects in a simulated high-voltage solar array

An emissive Langmuir probe was used to measure the potentials within the plasma sheath developed around a hole in a simulated solar array at voltages between 50 and 450 V. The hole sizes were larger than actual pinhole defects; the plasma density was in the 10,000 per cu cm range, which is considerably lower than the density of 1,000,000 per cu cm found at low-earth-orbit altitudes. Despite these inadequacies in the simulation, the experiments indicate that this type of probe is a useful diagnostic technique for investigating the plasma sheaths developing around pinhole defects.

Gabriel, S. B.↗
DOE OSTI2024

Neoclassical transport of tungsten ion bundles in total-f neoclassical gyrokinetic simulations of a whole-volume JET-like plasma

Here, the application of a bundling technique to model the diverse charge states of tungsten impurity species in total-f gyrokinetic simulations is demonstrated. The gyrokinetic bundling method strategically groups tungsten ions of similar charge, optimizing computational efficiency. The initial radial configuration of these bundles and their respective charges are derived from a coronal approximation and the quasi-neutrality of the plasma. A low-density JET H-mode like plasma is simulated using the neoclassical version of XGC across the entire plasma volume, spanning from the magnetic axis to the divertor. An accumulation of tungsten is observed at the pedestal top, as a result of low-Z tungsten ions moving inward from the scrape-off-layer into the core region and high-Z tungsten ions moving outward from the core into the pedestal. This organization of the fluxes cannot be captured by a single tungsten-ion simulation. Large up-down poloidal asymmetries of tungsten form in the pedestal and strongly influence the direction of neoclassical fluxes. The temperature screening effect and its correlation with asymmetries are analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Studies on counterstreaming plasma expansion

Recent studies on counterstreaming plasma expansions are summarized. The basic phenomenon of plasma expansion is reviewed, and results from one-dimensional simulations of counterstreaming plasma expansion are discussed. Results from simulations based on an electrostatic particle-in-cell code, in which the dynamics of both the electrons and ions are exactly followed, are discussed. The formation of electrostatic shocks is addressed. Finally, results are presented on the ionospheric plasma expansion along the geomagnetic flux tubes by solving the hydrodynamic equations.

Singh, N.↗
DOE OSTI2024

Assessing the impact of alpha particles on thermal confinement in JET D-T plasmas through global GENE-Tango simulations

The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially global nonlinear electromagnetic GENE-Tango simulations, we compare results with/without alpha particles and alpha heating. Our findings reveal that alpha particles have a negligible impact on turbulent transport, with GENE-Tango converging to similar plasma profiles regardless of their inclusion as a kinetic species in GENE. On the other hand, alpha heating is found to contribute to the peaking of the electron temperature profiles, leading to a 1 keV drop on the on-axis electron temperature when alpha heating is neglected in Tango. The minimal impact of alpha particles on turbulent transport in this JET discharge–despite this being the shot with the highest fusion output–is attributed to the low content of fusion alpha in this discharge. To assess the potential impact of alpha particles on turbulent transport in regimes with higher alpha particle density, as expected in ITER and fusion reactors, we artificially increased the alpha particle concentration to levels expected for ITER. By performing global nonlinear GENE standalone simulations, we found that increasing the alpha particle density beyond five times the nominal value lead to significant overall turbulence destabilization. These results demonstrate that an increased alpha particle concentration can significantly impact transport properties under simulated JET experimental conditions. However, these findings cannot be directly extrapolated to ITER due to the substantial differences in parameters such as plasma size, magnetic field, plasma current, and thermal pressure.

energetic particles↗
NASA NTRS1992

Interplanetary magnetic field changes and condensations in comet Halley's plasma tail

In a time-dependent three dimensional MHD simulation for cometary plasmas, Schmidt-Voigt (1989) could observe the formation of condensations in the plasma tail after a 90 degree change in the interplanetary magnetic field (IMF) sweeping over the comet. We investigated the IMF measurements of the Vega SC in the vicinity of the comet Halley for 90 degree changes in the clock angle and studied the relation between them and optical observations of condensations in the plasma tail. For the time interval 24 Feb. 86 to 14 Mar. 86, we could not find a correlation between such changes and the release of condensations from the cometary head.

Delva, Magda↗
NASA NTRS1981

Numerical simulation of the generation of electrostatic turbulence in the magnetotail

A two-dimensional plasma model is used to investigate the development of electrostatic turbulence in a magnetized plasma from plasma instabilities. The simulation consists of following the motion of 100,000 ions in their self-consistent electrostatic field. The electrons are treated as a constant neutralizing background. The instabilities modeled are driven by a ring-type velocity distribution and by interpenetrating ion beams in a time-variable magnetic field. Instability growth times are of the order of an ion gyroperiod in the case of the ring distribution and of the order of an ion plasma period in the case of the beam simulation. Maximum potential differences generated are of the order of the ion kinetic energies. These simulations demonstrate the cascade of wave energy to long wavelengths, thus showing the E x B turbulence can be generated from plasma microinstabilities. After the free energy feeding, the instabilities are exhausted, and wave energy at wavelengths less than an ion gyrodiameter decays quickly to equilibrium levels, while longer wavelength modes persist for much longer times. In one model with a time dependent, but spatially uniform, magnetic field the electric field energy at long wavelengths appeared to increase as a result of the increase of the magnetic field.

Swift, D. W.↗
NASA NTRS2018

Investigating the Response of Loop Plasma to Nanoflare Heating Using RADYN Simulations

We present the results of 1D hydrodynamic simulations of coronal loops that are subject to nanoflares, caused by either in situ thermal heating or nonthermal electron (NTE) beams. The synthesized intensity and Doppler shifts can be directly compared with Interface Region Imaging Spectrograph (IRIS) and Atmospheric Imaging Assembly (AIA) observations of rapid variability in the transition region (TR) of coronal loops, associated with transient coronal heating. We find that NTEs with high enough low-energy cutoff (EC) deposit energy in the lower TR and chromosphere, causing blueshifts (up to approximately 20 kilometers per second) in the IRIS Si IV lines, which thermal conduction cannot reproduce. The EC threshold value for the blueshifts depends on the total energy of the events (approximately 5 kiloelectronvolts for 1024 ergs, up to 15 kiloelectronvolts for 1025 ergs). The observed footpoint emission intensity and flows, combined with the simulations, can provide constraints on both the energy of the heating event and EC. The response of the loop plasma to nanoflares depends crucially on the electron density: significant Si IV intensity enhancements and flows are observed only for initially low-density loops (less than 10 (sup 9) per cubic centimeter). This provides a possible explanation of the relative scarcity of observations of significant moss variability. While the TR response to single heating episodes can be clearly observed, the predicted coronal emission (AIA 94 Angstroms) for single strands is below current detectability and can only be observed when several strands are heated closely in time. Finally, we show that the analysis of the IRIS Mg II chromospheric lines can help further constrain the properties of the heating mechanisms.

Sun: Corona↗
DOE OSTI2024

Development of spatially and temporally resolved electron density measurements for the assessment of radiation hydrodynamics simulations of laboratory X-ray photoionized plasmas

The photoionized plasma gas cell experiment is an established platform we use to make at-parameter (ξ > > 1 ergs cm s –1 ) measurements of plasma properties with application to high-energy astrophysical systems. We model the experiments with 1D radiation hydrodynamics simulations using the HELIOS-CR code to inform our understanding and assist in the interpretation of results. The simulations predict that the bulk of the plasma is in a quasi-uniform and hydrodynamically unperturbed state throughout the duration of the experiment. To evaluate this prediction, we introduced a photonic Doppler velocimetry (PDV) diagnostic to measure spatially and temporally resolved plasma electron density. The initial measurements were successful but had limitations that made model-data comparisons challenging. To address this, we re-designed the gas cell PDV diagnostic and doubled the number of measurement locations to sample across two thirds of the depth of the cell. We also present a comparison of the results from the upgraded PDV diagnostic to the HELIOS-CR simulations for the first time. As a result, the experimental data confirms the prediction of an unperturbed region in the bulk of the plasma but reveals discrepancies in the time evolution and spatial distribution of the simulated electron density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

STORM: Scrape-off layer turbulence in tokamak fusion reactors

The scrape-off layer of a tokamak fusion reactor carries the plasma exhaust from the hot core plasma to the material surfaces of the reactor vessel. The heat loads imposed by the exhaust are a critical limit on the performance of fusion power plants. Turbulent transport of the plasma regulates the width of the scrape-off layer plasma and must be modelled to understand the intensity of these heat loads. STORM is a plasma turbulence code capable of simulating three dimensional turbulence across the full scrape-off layer of a tokamak fusion reactor, using a drift reduced, collisional fluid model. STORM uses mostly finite difference schemes, with a staggered grid in the direction parallel to the magnetic field. We describe the model, geometry and initialisation options used by STORM, as well as the numerical methods, which are implemented using the BOUT++ plasma simulation framework. BOUT++ has been enhanced alongside the development of STORM, providing better support for staggered grid methods. We summarise these enhancements, including a detailed explanation of the parallel derivative methods, which underwent a major update for version 4 of BOUT++.

BOUT++↗
NASA NTRS1989

A model for plasma volume changes during short duration spaceflight

It is well established that plasma volume decreases during spaceflight and simulated weightlessness (bedrest). The decrement in plasma volume is thought to contribute to the orthostatic intolerance that has been observed in some crew members following spaceflight. To date, no studies have evaluated the effectiveness of fluid countermeasures of varying osmolality in the restoration of plasma volume and orthostatic tolerance in a controlled study. The overall objectives of this project were to: (1) provide a model that would rapidly and safely produce a fluid loss comparable to that which occurs during short duration spaceflight; and (2) design a study that would determine the optimal drink solution to restore orthostatic tolerance and describe the mechanism(s) whereby orthostatic tolerance is restored. In summary, Lasix can be used as a way of simulating the plasma volume changes that occur during short duration spaceflight. The total loss of plasma is comparable to spaceflight. Lasix is fast acting, and has relatively few side effects. The present design for evaluating the optimal fluid countermeasures will have important implications in restoring orthostatic tolerance and function in the latter stages of spaceflight when it is essential for safe operation of the spacecraft.

Davis, John E.↗
DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗
NASA NTRS1999

A Simulation to Study Speed Distributions in a Solar Plasma

We (Peter Cheeseman of NASA Ames/Caelum Research) & Jose Luis Alvarellos of the SJSU Physics Department/SJSU Foundation.) 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 1 we introduce the problem and in section 2 we discuss the approach to try to solve the problem: i.e., a molecular dynamics approach. In section 3 we provide details about the integration method, and any simplifications that can be applied to the problem. In section 4 (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 5 we summarize our findings and state our conclusions. Sections 6 and 7 provide the list of figures, reference material and acknowledgments respectively.

Cheeseman, Peter↗
DOE OSTI2026

Semiclassical treatment of bottomonium suppression and regeneration in 𝑝 + Pb collisions

Here, we study bottomonium suppression in 𝑝 + Pb relative to 𝑝 + 𝑝 collisions at center-of-mass energies of $\sqrt{s_{NN}}$ = 5.02 and 8.16 TeV. Specifically, we combine cold nuclear matter effects (nuclear modifications of the parton densities, energy loss, and momentum broadening) with those from hot nuclear matter (suppression and regeneration) by implementing the formation of a quark-gluon plasma in hydrodynamic simulations. Bottomonium transport in the quark-gluon plasma is evaluated semiclassically, employing two different reaction rates. The first includes quasifree inelastic scattering and gluodissociation employing a perturbative coupling to the medium. The second is based on in-medium 𝑇-matrix calculations where the input potential is constrained by lattice quantum chromodynamics to extract the bottomonium masses and dissociation rates. These semiclassical results are compared to previous calculations in an open quantum system approach and to the experimental data. Predictions for 𝜒 𝑏 suppression at $\sqrt{s_{NN}}$ = 8.16 TeV are also presented.

Physics - Physics of elementary particles and fiel↗
NASA NTRS1987

Solar transition region and coronal response to heating rate perturbations

The solar transition region is in a dynamic state characterized by impulsively upflowing plasma and continually downflowing plasma. Using numerical simulations, the conjecture that the areas of downflowing plasma are simply the base regions of coronal loops in which the heating rate is gradually decreasing and the areas of upflowing plasma are the base regions of coronal loops in which the heating rate is gradually increasing is examined. The calculations suggest that gradually reducing or increasing the heating in a magnetic flux tube will not result in plasma motions that are similar to those that are observed at high spatial resolution in the UV.

Mariska, John T.↗
NASA NTRS1986

The plasma dynamics of hypersonic spacecraft: Applications of laboratory simulations and active in situ experiments

Attempts to gain an understanding of spacecraft plasma dynamics via experimental investigation of the interaction between artificially synthesized, collisionless, flowing plasmas and laboratory test bodies date back to the early 1960's. In the past 25 years, a number of researchers have succeeded in simulating certain limited aspects of the complex spacecraft-space plasma interaction reasonably well. Theoretical treatments have also provided limited models of the phenomena. Several active experiments were recently conducted from the space shuttle that specifically attempted to observe the Orbiter-ionospheric interaction. These experiments have contributed greatly to an appreciation for the complexity of spacecraft-space plasma interaction but, so far, have answered few questions. Therefore, even though the plasma dynamics of hypersonic spacecraft is fundamental to space technology, it remains largely an open issue. A brief overview is provided of the primary results from previous ground-based experimental investigations and the preliminary results of investigations conducted on the STS-3 and Spacelab 2 missions. In addition, several, as yet unexplained, aspects of the spacecraft-space plasma interaction are suggested for future research.

Stone, N. H.↗
NASA NTRS1982

Numerical simulation of the interaction of the plasma sheet with the lobes of the Earth's magnetotail

Codes involving one and two spatial dimensions and three velocity dimensions were used to model the Earth's magnetotail. It was shown that the magnetotail can become inflated as a consequence of low energy plasma convection toward the neutral plane. The computer study exhibits a conversion of both magnetic field energy and of energy supplied by the convection electric field into particle energy. The numerical simulations suggest that much of the magnetotail substorm morphology may be a simple consequence of an increase, followed by a decrease, in the convection electric field, without the requirement of any magnetospheric size scale plasma instability or other disruptive processes. It is also concluded that the presence of the convection electric field and a continuing replenishment of low energy particles in the magnetotail are both necessary for maintenance of the magnetotail.

Swift, D. W.↗
NASA NTRS2007

Numerical Simulations of Flow Separation Control in Low-Pressure Turbines using Plasma Actuators

A recently introduced phenomenological model to simulate flow control applications using plasma actuators has been further developed and improved in order to expand its use to complicated actuator geometries. The new modeling approach eliminates the requirement of an empirical charge density distribution shape by using the embedded electrode as a source for the charge density. The resulting model is validated against a flat plate experiment with quiescent environment. The modeling approach incorporates the effect of the plasma actuators on the external flow into Navier Stokes computations as a body force vector which is obtained as a product of the net charge density and the electric field. The model solves the Maxwell equation to obtain the electric field due to the applied AC voltage at the electrodes and an additional equation for the charge density distribution representing the plasma density. The new modeling approach solves the charge density equation in the computational domain assuming the embedded electrode as a source therefore automatically generating a charge density distribution on the surface exposed to the flow similar to that observed in the experiments without explicitly specifying an empirical distribution. The model is validated against a flat plate experiment with quiescent environment.

Suzen, Y. B.↗