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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impurity transport in PISCES-RF

Linear plasma devices (LPD) utilizing a helicon plasma source, a high density light ion source, can generate impurities due to progressive erosion of the radio frequency (RF) transmission window caused by rectified sheath voltage. These source-born impurities can entrain and be transported by the plasma toward a target, affecting plasma-material interaction studies. Earlier work on material testing in Prototype-Materials Plasma Exposure eXperiment at ORNL revealed significant source impurity deposition on downstream targets. However, using a similar RF source, no target impurity deposition is observed in Plasma Interaction Surface Component Experimental Station (PISCES)-RF despite evidence of RF window erosion in the source region, thereby motivating the present work. Experimentally, using various magnetic field configurations upstream of the PISCES-RF plasma source and seeding titanium (Ti) impurities at various axial locations, impurity transport and deposition along the machine axis were investigated. It was found that Ti deposition was localized to the side of the plasma source where the Ti impurity was seeded. In contrast, aluminum (Al) deposition, originating from the sputtering of the helicon window, occurred predominantly upstream of the plasma source, suggesting an asymmetry in the axial transport of eroded RF window material. These observations suggest a stagnation of the parallel plasma flow immediately downstream of the plasma source, with impurity ions remaining unmagnetized near the source upstream. Al deposition in magnetic field-free regions in PISCES-RF indicates that sputtered Al impurities likely remained neutral due to their large ionization mean-free path under PISCES-RF conditions. Plasma modeling and simulation supported this, indicating that Al-neutrals transport toward the helicon source upstream for low electron density cases. It was found that the Larmor radius of the Al ions was greater than the plasma radius towards the source upstream and remained weakly magnetized in PISCES-RF, meaning that plasma source-born impurities are not efficiently entrained in the plasma flow. These findings provide critical insights into impurity transport in helicon plasma-based LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Observations of non-linear plasmon damping in dense plasmas

Here, we present simulations using finite-temperature density-functional-theory molecular-dynamics to calculate dynamic dielectric properties in warm dense aluminum. The comparison between exchange-correlation functionals in the Perdew, Burke, Ernzerhof approximation, Strongly Constrained and Appropriately Normed Semilocal Density Functional, and Heyd, Scuseria, Ernzerhof (HSE) approximation indicates evident differences in the electron transition energies, dc conductivity, and Lorenz number. The HSE calculations show excellent agreement with x-ray scattering data as well as dc conductivity and absorption measurements. These findings demonstrate non-Drude behavior of the dynamic conductivity above the Cooper minimum that needs to be taken into account to determine optical properties in the warm dense matter regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

iVPIC: A low-­dispersion, energy-­conserving relativistic PIC solver for LPI simulations

We have developed a novel low-­dispersion, exactly energy-­conserving PIC algorithm for the relativistic Vlasov-­Maxwell system. The approach features an exact energy conservation theorem while preserving the favorable performance and numerical dispersion properties of explicit PIC. The new algorithm has the potential to enable much longer laser-­plasma-­interaction (LPI) simulations than are currently possible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

TRANSP-based optimization towards tokamak scenario development

An optimization approach that incorporates the predictive transport code TRANSP is proposed for tokamak scenario development. Optimization methods are often employed to develop open-loop control strategies to aid access to high performance tokamak scenarios. Generally, the optimization approaches use control-oriented models, i.e. models that are reduced in complexity and prediction accuracy as compared to physics-oriented transport codes such as TRANSP. In the presented approach, an optimization procedure using the TRANSP code to simulate the tokamak plasma is considered for improved predictive capabilities. As a test case, the neutral beam injection (NBI) power is optimized to develop a control strategy that maximizes the noninductive current fraction during the ramp-up phase for NSTX-U. Simulation studies towards the achievement of noninductive ramp-up in NSTX-U have already been carried out with the TRANSP code. The optimization-based approach presented in this work is used to maximize the noninductive current fraction during ramp-up in NSTX-U, demonstrating that the scenario development task can be automated. An additional test case considers optimization of the current ramp rate in DIII-D for obtaining a stationary plasma characterized by a flat loop voltage profile in the flattop phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Full- f gyrokinetic simulation of turbulence in a helical open-field-line plasma

Curvature-driven turbulence in a helical open-field-line plasma is investigated using electrostatic five-dimensional gyrokinetic continuum simulations in an all-bad-curvature helical-slab geometry. Parameters for a National Spherical Torus Experiment scrape-off-layer plasma are used in the model. Here, the formation and convective radial transport of plasma blobs is observed, and it is shown that the radial particle-transport levels are several times higher than diffusive Bohm-transport estimates. By reducing the strength of the poloidal magnetic field, the profile of the heat flux to the divertor plate is observed to broaden.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Effects of plasmoid formation on sawtooth process in a tokamak

For realistic values of the Lundquist number in tokamak plasmas, 1/1 magnetic islands lead to the formation of a secondary thin current sheet, which breaks up into a chain of small magnetic islands, called plasmoids. The role of plasmoid dynamics during the sawtooth reconnection process in fusion plasmas remains an unresolved issue. In this study, systematic simulations are performed to investigate the resistive internal kink mode using the full resistive magnetohydrodynamics equations implemented in the NIMROD code in a simplified tokamak geometry. For Lundquist number S≥1.6×10 7 , the secondary current sheet is found to be unstable to plasmoids during the nonlinear resistive kink mode evolution with a critical aspect ratio of the current sheet of ~70. The merging of small plasmoids leads to the formation of a monster plasmoid that can significantly affect the primary island evolution. This may provide an explanation for the partial reconnection observed in sawtooth experiments.

Computer simulation↗
NASA NTRS1993

A Remote Laser-mass Spectrometer for Determination of Elemental Composition

Determination of the elemental composition of lunar, asteroid, and planetary surfaces is a major concern for science and resource utilization of space. The science associated with the development of a satellite or lunar rover laser-mass spectrometer instrument is presented here. The instrument would include a pulsed laser with sufficient energy to create a plasma on a remote surface. Ions ejected from this plasma travel back to the spacecraft or rover, where they are analyzed by a time-of-flight mass spectrometer, giving the elemental and isotope composition. This concept is based on the LIMA-D instrument on board the former Soviet Union Phobos-88 spacecraft sent to Mars. A laser-mass spectrometer placed on a rover or satellite would substantially improve the data return over alternative techniques. The spatial resolution would be centimeters, and a complete mass spectrum could be achieved in one laser shot. An experiment is described that demonstrates these features. A 400 mj Nd:YAG laser is focused, to an intensity of 10(exp 11) w/sq cm, onto a Al, Ag, Cu, Ge, or lunar simulant target. A plasma forms from which ions are ejected. Some of these ions travel down an 18-m evacuated flight tube to a microchannel plate detector. Alternatively, the ions are captured by an ion trap where they are stored until pulsed into a 1-m time-of-flight mass spectrometer, giving the elemental composition of the remote surface. A television camera monitors the plasma plume shape, and a photodiode monitors the temporal plasma emission . With this system, ions of Al, Ag, Cu, Ge, and lunar simulant have been detected at 18 m. The mass spectrum from the ion trap and 1-m time-of-flight tube will be presented.

Deyoung, R. J.↗
NASA NTRS1992

The scaling of relativistic double-year widths - Poisson-Vlasov solutions and particle-in-cell simulations

The study of relativistic plasma double layers is described through the solution of the one-dimensional, unmagnetized, steady-state Poisson-Vlasov equations and by means of one-dimensional, unmagnetized, particle-in-cell simulations. The thickness vs potential-drop scaling law is extended to relativistic potential drops and relativistic plasma temperatures. The transition in the scaling law for 'strong' double layers suggested by analytical two-beam models by Carlqvist (1982) is confirmed, and causality problems of standard double-layer simulation techniques applied to relativistic plasma systems are discussed.

Sulkanen, Martin E.↗
NASA NTRS2010

Jovian Plasma Torus Interaction with Europa: 3D Hybrid Kinetic Simulation. First results

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa-moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements, (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy etal.,2007;Shematovichetal.,2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyro radius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions).Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization,electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma:(a) with O(++) ions; (b) with O(++) and S(++) ions. The majority of O2 atmosphere is thermal with an extended cold population (Cassidyetal.,2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.

Lipatov, A. S.↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. II. General Relativistic Module

Black hole (BH) environments often host plasmas that are fully collisionless or contain intrinsically collisionless regions, including relativistic jets and coronae, where particle energization is ubiquitous. Capturing the physics of these systems requires numerical methods capable of modeling relativistic, magnetized, collisionless plasmas in strong gravitational fields. In this work, we introduce the general relativistic module for Entity—the first open-source, coordinate-agnostic, performance-portable, particle-in-cell code. The code enables fast axisymmetric simulations of collisionless plasmas around BHs on any modern high-performance computing architecture (both GPUs and CPUs).

Galishnikova, Alisa [Flatiron Institute, New York,↗
DOE OSTI2019

TRANSP-based Optimization Towards Tokamak Scenario Development

An optimization approach that incorporates the predictive transport code TRANSP is proposed for tokamak scenario development. Optimization methods are often employed to develop open-loop control strategies to aid access to high performance tokamak scenarios. In general, the optimization approaches use control-oriented models, i.e. models that are reduced in complexity and prediction accuracy as compared to physics-oriented transport codes such as TRANSP. In the presented approach, an optimization procedure using the TRANSP code to simulate the tokamak plasma is considered for improved predictive capabilities. As a test case, the neutral beam injection (NBI) power is optimized to develop a control strategy that maximizes the non-inductive current fraction during the ramp-up phase for NSTX-U. Simulation studies towards the achievement of non-inductive ramp up in NSTX-U have already been carried out with the TRANSP code. The optimization-based approach proposed in this work is used to maximize the non-inductive current fraction during ramp-up in NSTX-U, demonstrating that the scenario development task can be automated. An additional test case considers optimization of the current ramp rate in DIII-D for obtaining a stationary plasma characterized bya flat loop voltage profile in the flattop phase.

DIII-D↗
DOE OSTI2025

Verification of a Fluid-Based Plasma-Edge Model Within the Multiphysics Object-Oriented Simulation Environment (MOOSE) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. These models must incorporate the interconnected phenomena of these devices, including plasma physics, neutronics, first wall interactions, and tritium transport. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification studies have been conducted using newly developed magnetic plasma capabilities. These involved convergence analyses utilizing the method of manufactured solutions to verify new operators and case studies. A modular approach was taken here to demonstrate increasingly complicated simulation scenarios, which included a singular fluid with uniform magnetic field case, a singular fluid with spatially varying magnetic field case, and a coupled multifluid case.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

UEDGE modelling of detached divertor operation for long-leg divertor geometries in ARC

In this work, simulations of the ARC scrape-off layer and divertor configuration are carried out using UEDGE to perform the first power-handling assessment for the divertor design. A convective–diffusive anomalous transport model is implemented in the code to reproduce the radial plasma profiles anticipated for the device at the outer mid-plane. Levels of exhaust power from the core are varied, and a 0.5% fixed fraction neon impurity seeding is included in investigations to find steady-state detached solutions. Initial studies employing a super-X divertor configuration have shown that a stable detached divertor-operational window exists for reactor exhaust powers in the range of 72–92 MW, which is close to projected levels of exhaust power of 105 MW across the separatrix anticipated for ARC. Under these conditions, the divertor heat flux is fully dissipated by radiation and radial losses to the sidewalls of the divertor channel, leading to acceptable levels of power flux density on plasma-facing components. Simulations are being extended to study the performance of the X-point target geometry and to explore the sensitivity to the physics model and input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Cold Plasma Measurements

We have continued the simulation campaign in support of our ongoing magnetospheric cold plasma research project. This project aims to develop the next-generation particle instruments to measure the properties of the cold particle populations in the Earth’s magnetosphere. For this purpose, simulations have been performed with a Particle-In-Cell (PIC) code called the Curvilinear PIC (CPIC). The code is formulated in curvilinear geometry and couples the standard PIC algorithm with algorithms for the generation and adaptation of the underlaying computational mesh. It conforms to complex objects like spacecraft and it can place more grid points in regions where higher resolution is needed. The code also features a scalable solver based on the multigrid algorithm and it is fully parallelized via domain decomposition and MPI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Relativistic Eulerian Vlasov simulations of the amplification of seed pulses by Brillouin backscattering in plasmas

We apply an Eulerian Vlasov code to study the amplification by Brillouin scattering of a short seed laser pulse by a long pump laser pulse in an underdense plasma. The stimulated Brillouin backscattering interaction is the coupling of the pump and seed electromagnetic waves propagating in opposite directions, and the ion plasma wave. The code solves the one-dimensional relativistic Vlasov-Maxwell set of equations. Large amplitude ion waves are generated. In the simulations we present, the density plateau of the plasma is n{sub e}=0.3 n{sub c} (n{sub c} is the critical density), which excludes spurious stimulated Raman scattering amplification (which can occur only if n{sub e}<n{sub c}/4). We also varied the duration and/or amplitude of the short input seed pulse to study how these influence its subsequent behaviour. An initially broad pulse grows more rapidly than an initially narrow pulse. Furthermore, for an initially broader seed pulse, towards the end of the simulation, it is seen to become narrower and to gradually detach from the trailing signal. On the contrary, initially very narrow seed pulses are seen to broaden. The absence of noise in the Vlasov simulations allows to simulate long plasma amplifier lengths, and to follow the evolution of the system with a fully kinetic description and with an accurate representation of the phase-space structures of distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

DBD Plasma Actuators for Flow Control in Air Vehicles and Jet Engines - Simulation of Flight Conditions in Test Chambers by Density Matching

Dielectric Barrier Discharge (DBD) Plasma actuators for active flow control in aircraft and jet engines need to be tested in the laboratory to characterize their performance at flight operating conditions. DBD plasma actuators generate a wall-jet electronically by creating weakly ionized plasma, therefore their performance is affected by gas discharge properties, which, in turn, depend on the pressure and temperature at the actuator placement location. Characterization of actuators is initially performed in a laboratory chamber without external flow. The pressure and temperature at the actuator flight operation conditions need to be simultaneously set in the chamber. A simplified approach is desired. It is assumed that the plasma discharge depends only on the gas density, while other temperature effects are assumed to be negligible. Therefore, tests can be performed at room temperature with chamber pressure set to yield the same density as in operating flight conditions. The needed chamber pressures are shown for altitude flight of an air vehicle and for jet engines at sea-level takeoff and altitude cruise conditions. Atmospheric flight conditions are calculated from standard atmosphere with and without shock waves. The engine data was obtained from four generic engine models; 300-, 150-, and 50-passenger (PAX) aircraft engines, and a military jet-fighter engine. The static and total pressure, temperature, and density distributions along the engine were calculated for sea-level takeoff and for altitude cruise conditions. The corresponding chamber pressures needed to test the actuators were calculated. The results show that, to simulate engine component flows at in-flight conditions, plasma actuator should be tested over a wide range of pressures. For the four model engines the range is from 12.4 to 0.03 atm, depending on the placement of the actuator in the engine. For example, if a DBD plasma actuator is to be placed at the compressor exit of a 300 PAX engine, it has to be tested at 12.4 atm for takeoff, and 6 atm for cruise conditions. If it is to be placed at the low-pressure turbine, it has to be tested at 0.5 and 0.2 atm, respectively. These results have implications for the feasibility and design of DBD plasma actuators for jet engine flow control applications. In addition, the distributions of unit Reynolds number, Mach number, and velocity along the engine are provided. The engine models are non-proprietary and this information can be used for evaluation of other types of actuators and for other purposes.

Ashpis, David E.↗
DOE OSTI2025

Electron temperature gradient instability and transport analysis in NSTX and NSTX-U plasmas

Extensive linear and nonlinear simulations to study electron temperature gradient (ETG) stability and thermal transport in National Spherical Torus Experiment (NSTX) and NSTX-U plasmas were performed using the fully electromagnetic gyrokinetic code CGYRO. Linear simulations were performed to determine ETG thresholds in different discharges, showing that ETG modes in spherical tokamaks can present different scalings compared to conventional aspect-ratio tokamaks. Nonlinear gyrokinetic simulations were conducted for selected cases to calculate electron thermal transport and compare to experimental values. Results are also compared with those of ETG modes in the multi-mode model and the Trapped-Gyro-Landau-Fluid reduced model codes, to better understand their applicability in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗