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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF

Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.

Smolentsev, Sergey↗
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↗
DOE OSTI2017

The accurate particle tracer code

The Accurate Particle Tracer (APT) code is designed for systematic large-scale applications of geometric algorithms for particle dynamical simulations. Based on a large variety of advanced geometric algorithms, APT possesses long-term numerical accuracy and stability, which are critical for solving multi-scale and nonlinear problems. To provide a flexible and convenient I/O interface, the libraries of Lua and Hdf5 are used. Following a three-step procedure, users can efficiently extend the libraries of electromagnetic configurations, external non-electromagnetic forces, particle pushers, and initialization approaches by use of the extendible module. APT has been used in simulations of key physical problems, such as runaway electrons in tokamaks and energetic particles in Van Allen belt. As an important realization, the APT-SW version has been successfully distributed on the world’s fastest computer, the Sunway TaihuLight supercomputer, by supporting master–slave architecture of Sunway many-core processors. Here, based on large-scale simulations of a runaway beam under parameters of the ITER tokamak, it is revealed that the magnetic ripple field can disperse the pitch-angle distribution significantly and improve the confinement of energetic runaway beam on the same time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2019

Closed-loop simulation with Grad-Shafranov equilibrium evolution for plasma control system development

Plasma control design increasingly depends on fast simulations able to connect to operational plasma control system (PCS) software for iterative development. GSevolve is a nonlinear free-boundary simulation that evolves the Grad-Shafranov equilibrium including current and pressure profiles, and can connect to all versions of the DIII-D PCS operational on devices around the world. Its ability to simulate the evolution of current continuously from open to closed field lines enables GSevolve to simulate full-discharge control from breakdown through plasma termination. GSevolve is compliant with the Plasma Control System Simulation Platform (PCSSP), a Matlab/Simulink®-based environment that serves as the primary framework for development of the ITER Plasma Control System (PCS). The code employs a novel response function-updating scheme to achieve fast simulations. The plasma evolves linearly until a significant shape change triggers a nonlinear update of the response while enforcing a prescribed accuracy. Old responses are recalled if the plasma changes back to an old shape. The equilibrium is a function of a state vector containing axisymmetric conductor currents and parameters related to parallel and perpendicular currents, which evolves according to an efficient state space circuit formulation. Control design and validation results are presented for simulations of DIII-D and EAST, connected to their corresponding versions of the DIII-D PCS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulated Night- Shift Schedule Disrupts the Plasma Lipidome and Reveals Early Markers of Cardiovascular Disease Risk

Introduction: The circadian system coordinates daily rhythms in lipid metabolism, storage and utilization. Disruptions of internal circadian rhythms due to altered sleep/wake schedules, such as in night-shift work, have been implicated in increased risk of cardiovascular disease and metabolic disorders. To determine the impact of a night-shift schedule on the human blood plasma lipidome, an in-laboratory simulated shift work study was conducted. Methods: Fourteen healthy young adults were assigned to 3 days of either a simulated day or night-shift schedule, followed by a 24-h constant routine protocol with fixed environmental conditions, hourly isocaloric snacks, and constant wakefulness to investigate endogenous circadian rhythms. Blood plasma samples collected at 3-h intervals were subjected to untargeted lipidomics analysis. Results: More than 400 lipids were identified and quantified across 21 subclasses. Focusing on lipids with low between-subject variation per shift condition, alterations in the circulating plasma lipidome revealed generally increased mean triglyceride levels and decreased mean phospholipid levels after night-shift relative to day-shift. The circadian rhythms of triglycerides containing odd chain fatty acids peaked earlier during constant routine after night-shift. Regardless of shift condition, triglycerides tended to either peak or be depleted at 16:30 h, with chain-specific differences associated with the direction of change. Discussion: The simulated night-shift schedule was associated with altered temporal patterns in the lipidome. This may be premorbid to the elevated cardiovascular risk that has been found epidemiologically in night-shift workers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Transient Relativistic Plasma Grating to Tailor High-Power Laser Fields, Wakefield Plasma Waves, and Electron Injection

In this work, we show the first experiment of a transverse laser interference for electron injection into the laser plasma accelerators. Simulations show such an injection is different from previous methods, as electrons are trapped into later acceleration buckets other than the leading ones. With optimal plasma tapering, the dephasing limit of such unprecedented electron beams could be potentially increased by an order of magnitude. In simulations, the interference drives a relativistic plasma grating, which triggers the splitting of relativistic-intensity laser pulses and wakefield. Consequently, spatially dual electron beams are accelerated, as also confirmed by the experiment.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Gyrokinetic continuum simulation of turbulence in a straight open-field-line plasma

Here, five-dimensional gyrokinetic continuum simulations of electrostatic plasma turbulence in a straight, open-field-line geometry have been performed using a full- discontinuous-Galerkin approach implemented in the Gkeyll code. While various simplifications have been used for now, such as long-wavelength approximations in the gyrokinetic Poisson equation and the Hamiltonian, these simulations include the basic elements of a fusion-device scrape-off layer: localised sources to model plasma outflow from the core, cross-field turbulent transport, parallel flow along magnetic field lines, and parallel losses at the limiter or divertor with sheath-model boundary conditions. The set of sheath-model boundary conditions used in the model allows currents to flow through the walls. In addition to details of the numerical approach, results from numerical simulations of turbulence in the Large Plasma Device, a linear device featuring straight magnetic field lines, are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2015

One-dimensional particle-in-cell simulations of electrostatic Bernstein waves in plasmas with kappa velocity distributions

Electrostatic Bernstein waves that propagate exactly perpendicularly to a static magnetic field in an electron-ion plasma are investigated using one-and-two-halves dimensional particle-in-cell simulations. An ion-to-electron mass ratio of m{sub i}/m{sub e} = 100 is used, allowing sufficient separation of the electron and ion time scales while still accounting for the ion dynamics without resorting to exceptionally long simulation run times. As a consequence of the mass ratio used, both the high frequency electron Bernstein wave and the lower frequency ion Bernstein wave are resolved within a single simulation run. The simulations presented here use isotropic three-dimensional kappa velocity distributions as well as the widely used Maxwellian velocity distribution, and the results from using each of these velocity distributions are analysed and compared. The behaviour of the Bernstein waves is found to be significantly dependent on the spectral index, κ, of the kappa distribution in all frequency domains of the Bernstein waves. In both the Maxwellian and kappa cases, spectral analysis of the electric field (wave) intensities, as a function of ω and k, show very good agreement between the simulation results and the linear dispersion relation for Bernstein waves. This agreement serves to validate the simulation techniques used, as well as the theory of Bernstein waves in plasmas with a kappa velocity distribution. The intensity of the field fluctuations in the simulations containing an abundance of superthermal particles, i.e., where the plasma has a kappa velocity distribution with a low kappa index, is slightly higher compared to the simulations of plasmas with higher kappa values. The plasmas with low kappa values also exhibit a broader region in frequency space of high intensity field fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Implementation and verification of a model linearized multi-species collision operator in the COGENT code

We implement the multi-species linearized model collision operator, based on the approach proposed by Kolesnikov, in the full-f continuum gyrokinetic code COGENT. We describe the operator and use it in several COGENT simulations of highly collisional plasmas. We analyze simulation results with the Braginskii fluid model to illustrate that COGENT recovers both friction and thermal forces. Finally, we simulate the neoclassical transport of heavy impurities with COGENT and illustrate agreement with previously published theoretical results.

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↗
DOE OSTI2017

Tersoff Benchmarking of Be-C-H Interatomic Potential

To demonstrate prolonged fusion power in ITERs (International Thermonuclear Experimental Reactors), the selection of appropriate plasma facing materials (PFMs) is a vital issue that requires the use of interatomic potentials to simulate the plasma interactions with the PFMs. So it is no surprise that the use of computer codes such as LAMMPS to perform studies in order to determine proper PFMs is necessary based on the complexity of such studies. In recent Tokamak advancements, beryllium was chosen as a first wall material due to its low Z (lower energy loss due to plasma contamination), its low tritium inventory and its oxygen gettering quality. With its high melting point and ability to withstand high heat fluxes, beryllium makes a fine choice for a PFM. In this work, Be-C-H (beryllium-carbon-hydrogen) bonds were investigated in a LAMMPS Tersoff simulation, where the interatomic potentials for Be-C-H were converted to LAMMPS Tersoff format, and were validated by calculating the cohesive energy (E{sub coh}), lattice constant (a), and atomic volume (V{sub 0}) for comparison to Density Functional Theory (DFT), Analytical Bond-Order Potentials (ABOP) and experiment potentials. This benchmark serves as the validation of performing Be-C-H bombardment simulations using LAMMPS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗