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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Multiscale Modeling of Deformation, Fracture and Failure of Fusion Materials and Structures Final Report

The main objective of the research performed at UCLA under grant number DE-FG02-03ER54708 was to develop a physical understanding of the underlying mechanisms that control high-temperature deformation and fracture of fusion structural materials. The key aspect is the synergistic interaction between thermo-mechanical damage phenomena (e.g. creep, fatigue, and fracture), and irradiation damage phenomena (helium and hydrogen bubbles, and irradiation-induced defects). Another thrust area of the research was on the development of constitutive equations for the mechanical behavior of irradiated structures in fusion components (mainly the First Wall and Plasma-Facing components), and in utilizing them in descriptions in large-scale Finite Element Modeling of fusion structures. A parallel experimental effort on thermo-mechanical and plasma-induced damage is designed to validate and confront the developed theoretical and computational predictions. The models have been experimentally validated through experimental tests in simulated plasma transients on specially-designed W samples. This grant started at UCLA on November 14, 2003 and has had successive continuations till November 14, 2017. Research performed under the grant was to develop a range of hierarchical models for the post-elastic deformation, fracture and failure of fusion structural materials. The developed multiscale modeling approach has been based on rigorous mathematical, physical, and computational methods at the forefront of computational materials science. At the fundamental and smallest length scale (nm-μm), we developed advanced rate theory and Monte Carlo approaches to model microstructure evolution, non-equilibrium phase transformations, and dislocation-defect interactions. Microscopic and mesoscopic models of radiation hardening, ductile-to-brittle-transition, post-elastic deformation, plastic instabilities, and fracture processes are based on Dislocation Dynamics (DD) and Grain Growth Dynamics (GGD). At the continuum level, we focused on modeling fracture and failure mechanisms for Virtual Integrated Testing (VISTA), and for understanding the limits of both materials-by-design versus structural component design through advancing the new approach of multiscale-multiphysics. The overall research has been strongly coupled with the U.S. national experimental program to ascertain and verify the range of investigated phenomena. In the following, we give a summary of research progress supported by this grant since its inception till November 14, 2017. Detailed research results are available in our publications, listed in the publications section of this report.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Informing ARC divertor design and plasma facing material selection through integrated modeling

This INFUSE 2023 project between UCLA and Commonwealth Fusion Systems used computer modeling to test whether tungsten materials can survive in CFS's ARC fusion reactor divertor. The team simulated plasma conditions and material responses, finding that tungsten-rhenium alloys resist grain growth better than pure tungsten, and that hydrogen buildup depends more on particle flux than temperature. The work helps CFS design durable plasma-facing components for their fusion power plant.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Verification of 5D continuum gyrokinetic code COGENT: Studies of kinetic drift wave instability

COGENT (Continuum Gyrokinetic Edge New Technology) is a kinetic plasma simulation code that is being developed by the edge simulation laboratory (ESL) collaboration. The original version of the code has been developed in a 4D phase space (2D configuration space and 2D velocity space) to address kinetic plasma phenomena in a complex magnetic field geometry including the core, magnetic separatrix, and scrape-off layer regions. This paper is focused on extending the original 4D phase space to a 5D phase space (3D2V) to address full kinetic turbulences in a tokamak edge region. Here, we report on the current status of 5D COGENT, which presently operates in a shear-less, simple slab geometry. As a verification study, we use the problems of collision-less drift wave instability (universal instability) including its modification in the presence of collisional effects. The simulation model includes the gyrokinetic equations for ion and electron species coupled to the long-wavelength limit of the gyro-Poisson equation. Collisional effects are represented by the Krook collision model. Linear analytical results for the drift mode growth rate and real frequency are recovered, and the non-linear stage is modelled and analysed as well. In addition, extensive 5D runs have been performed to address the effects of the drift wave instability on blob/filamentary structures characteristic of a tokamak edge. Finally, a helical-shaped potential perturbation is observed to grow exponentially in time while spinning around the filament axis with electron drift frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Final Technical Report: Magnetic Reconnection in High-Energy Laser-Produced Plasmas

This report describes the final results from the DOE Grant DE-SC0007168, “Fast Magnetic Reconnection in HED Laser-Produced Plasmas.” The recent generation of laboratory high-energy-density physics facilities has opened significant physics opportunities for experimentally modeling astrophysical plasmas. The goal of this proposal is to use these new tools to study fundamental problems in plasma physics and plasma astrophysics. Fundamental topics in this area involve study of the generation, amplification, and fate of magnetic fields, which are observed to pervade the plasma universe and govern its evolution. This project combined experiments at DOE laser facilities with kinetic plasma simulation to study these processes. The primary original goal of the project was to study magnetic reconnection using a new experimental platform, colliding magnetized laser-produced plasmas. However through a series of fortuitous discoveries, the work broadened out to allow significant advancement on multiple topics in laboratory astrophysics, including magnetic reconnection, Weibel instability, and collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗
DOE OSTI2023

Parallel Algebraic Multigrid for Fusion and Higher-Order PDEs

Multigrid methods play a key role in large-scale scientific simulation because they are among the fastest and most scalable approaches for solving the underlying sparse linear systems of equations that arise from a wide array of Partial Differential Equation (PDE) discretizations. Algebraic multigrid (AMG) is a special type of multigrid method that depends only on the description of the linear system, giving it better portability and broader applicability than geometric multigrid, as it requires no explicit knowledge of the problem geometry. Even though these methods are widely used today, there are still applications where further development is needed. In this report, we focus on PDEs with higher-order terms (e.g., fourth order), concentrating on a PDE that arises in tokamak edge plasma simulations (a tokamak is a machine that confines a plasma using magnetic fields and is believed to be the leading plasma confinement concept for future fusion power plants). General multigrid relaxes a linear system on coarser grids and reverses this process with interpolation, but standard AMG methods struggle with the aforementioned higher-order PDEs. We investigate cyclic coarsening and interpolation heuristics, as well as new iterative approximation methods of refining the solution at each grid to improve the existing multigrid approach. To this end, we ensure that these techniques are transferable to a parallelized setting with LLNL’s supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Fully Implict Magneto-hydrodynamics Simulations of Coaxial Plasma Accelerators

The resistive Magneto-Hydrodynamic (MHD) model describes the behavior of a strongly ionized plasma in the presence of external electric and magnetic fields. We developed a fully implicit MHD simulation tool to solve the resistive MHD governing equations in the context of a cell-centered finite-volume scheme. The primary objective of this study is to use the fully-implicit algorithm to obtain insights into the plasma acceleration and jet formation processes in Coaxial Plasma accelerators; electromagnetic acceleration devices that utilize self-induced magnetic fields to accelerate thermal plasmas to large velocities. We also carry out plasma-surface simulations in order to study the impact interactions when these high velocity plasma jets impinge on target material surfaces. Scaling studies are carried out to establish some basic functional relationships between the target-stagnation conditions and the current discharged between the coaxial electrodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

LPSE : A 3-D wave-based model of cross-beam energy transfer in laser-irradiated plasmas

The new component of the “laser–plasma simulation environment” (LPSE) described here is a practical numerical model that solves the coupled vector equations for the propagation of nearly monochromatic, coherent, lectromagnetic waves in inhomogeneous unmagnetized plasmas. It operates efficiently in the numerically challenging semiclassical regime where characteristic plasma scale lengths are many times greater than the wavelength of the light and solutions are highly oscillatory. Solutions can be obtained in one, two, or three spatial dimensions and time. Additionally, the model includes the effects of nonlinear coupling of electromagnetic waves to the low-frequency plasma perturbations (i.e., ion-acoustic response) that are responsible for stimulated Brillouin scattering. Induced plasma perturbations are assumed to be imposed on a prescribed large-scale inhomogeneous background that includes spatially varying plasma density and flow. Our code is directly relevant to the problem of cross-beam energy transfer in laser-driven inertial confinement fusion. Lastly, it may also be applicable in other areas where eikonal solutions of multicomponent wave equations (or coupled wave equations) are insufficient, such as optical scattering from ultrasound, electron dynamics in quantum devices or in nanoscale light–matter interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗
NASA NTRS1987

Plasma physics abstracts, 1 January - 31 December, 1986

Topics addressed include: ion-cyclotron waves; plasma waves; solar wind lithium releases; bow shock; Pi2 wave bursts; auroral kilometric radiation; ion energization; magnetic field corrections; electric fields; magnetospheric processes; electron acceleration; inner heliosphere; nightside auroral zone; computerized simulation; plasma wave turbulence; and magnetohydrodynamic waves in plasma sheets.

Gurnett, D. A.↗
DOE OSTI2020

Field-emission plasma enhancement of H₂–O₂ micro-combustion

In this paper, a detailed three-dimensional computational model is developed to assess the potential of micron-scale field emission dielectric barrier discharge (FE-DBD) plasma actuators in non-premixed microburners with channel heights in the range of L z = 0.25–0.75 mm. Results for H 2 –O 2 combustion with Reynolds number Re = 100 are studied for different configurations of plasma actuator arrays. Without plasma actuation, results agree with corresponding experiments, although in this case the burner typically suffers from incomplete combustion due to thermal quenching, radical quenching, and slow diffusive mixing. For the range of conditions simulated, plasma actuation increases combustion completeness from 25%–70% to 40%–80% by enhancing mixing, generating radicals, and Joule heating. Ionic wind resulting from the FE-DBD plasma disrupts the otherwise simple laminar flow in the burner, accelerating the growth of a mixing layer and enhancing heat transfer to channel walls. The larger, lower temperature flames reduce thermal gradients (and thereby thermal stresses) in the walls. Radical generation due to electron-impact dissociation extends the reaction zone in the diffusion flame and increases combustion completeness. In all cases, the additional heat release with plasma actuation (100–190 W) is significantly larger than the power needed to sustain the plasma (≲0.1 W), suggesting that plasma actuation may provide a significant advantage in energy generation using realistic portable micro-combustion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1971

Experimental investigations of a uranium plasma pertinent to a self-sustaining plasma source

The research is pertinent to the realization of a self-sustained fissioning plasma for applications such as nuclear propulsion, closed cycle MHD power generation using a plasma core reactor, and heat engines such as the nuclear piston engine, as well as the direct conversion of fission energy into optical radiation (nuclear pumped lasers). Diagnostic measurement methods and experimental devices simulating plasma core reactor conditions are discussed. Studies on the following topics are considered: (1) ballistic piston compressor (U-235); (2) high pressure uranium plasma (natural uranium); (3) sliding spark discharge (natural uranium); (4) fission fragment interaction (He-3 and U-235); and (5) nuclear pumped lasers (He-3 and U-235).

Schneider, R. T.↗
DOE OSTI2026

A Structure-Preserving Decorated Particle Method for the Vlasov-Poisson System

We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by macro-particles with position and velocity. In contrast, Scovel-Weinstein decorated particles involve additional shape degrees of freedom, while maintaining a finite-dimensional reduction with Hamiltonian structure inherited from the continuum model. Although the original work established this structure three decades ago, its computational potential has remained largely unexplored. We present a practical implementation of the Scovel-Weinstein model and compare it with a standard PIC algorithm. Numerical experiments demonstrate that macro-particles in standard PIC can be replaced by far fewer decorated particles while retaining comparable accuracy. This decorated particle approach offers a new structure-preserving paradigm for kinetic plasma simulation.

65M75, 70H05, 70G65↗
DOE OSTI2024

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm

For decades, the Vlasov–Darwin model has been recognized to be attractive for particle-in-cell (PIC) kinetic plasma simulations in non-radiative electromagnetic regimes, to avoid radiative noise issues and gain computational efficiency. However, the Darwin model results in an elliptic set of field equations that renders conventional explicit time integration unconditionally unstable. We explore a fully implicit PIC algorithm for the Vlasov–Darwin model in multiple dimensions, which overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. The finite-difference scheme for Darwin field equations and particle equations of motion is space–time-centered, employing particle sub-cycling and orbit-averaging. This algorithm conserves total energy, local charge, canonical-momentum in the ignorable direction, and preserves the Coulomb gauge exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. Finally, we demonstrate the accuracy and efficiency properties of the algorithm with various numerical experiments in 2D–3V.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

SIMULTANEOUS IRIS AND HINODE/EIS OBSERVATIONS AND MODELING OF THE 2014 OCTOBER 27 X2.0 CLASS FLARE

We present a study of the X2-class flare which occurred on 2014 October 27 and was observed with the Interface Region Imaging Spectrograph (IRIS) and the EUV Imaging Spectrometer (EIS) on board the Hinode satellite. Thanks to the high cadence and spatial resolution of the IRIS and EIS instruments, we are able to compare simultaneous observations of the Fe xxi 1354.08 Å and Fe xxiii 263.77 Å high-temperature emission (≳10 MK) in the flare ribbon during the chromospheric evaporation phase. We find that IRIS observes completely blueshifted Fe xxi line profiles, up to 200 km s{sup −1} during the rise phase of the flare, indicating that the site of the plasma upflows is resolved by IRIS. In contrast, the Fe xxiii line is often asymmetric, which we interpret as being due to the lower spatial resolution of EIS. Temperature estimates from SDO/AIA and Hinode/XRT show that hot emission (log(T[K]) > 7.2) is first concentrated at the footpoints before filling the loops. Density-sensitive lines from IRIS and EIS give estimates of electron number density of ≳10{sup 12} cm{sup −3} in the transition region lines and 10{sup 10} cm{sup −3} in the coronal lines during the impulsive phase. In order to compare the observational results against theoretical predictions, we have run a simulation of a flare loop undergoing heating using the HYDRAD 1D hydro code. We find that the simulated plasma parameters are close to the observed values that are obtained with IRIS, Hinode, and AIA. These results support an electron beam heating model rather than a purely thermal conduction model as the driving mechanism for this flare.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

ICRF mode conversion in three-ion species heating experiment and in flow drive experiment on the Alcator C-Mod tokamak

In recent three-ion species (majority D and H plus a trace level of 3He) ICRF heating experiments on Alcator C-Mod, double mode conversion on both sides of the 3He cyclotron resonance has been observed using the phase contrast imaging (PCI) system. The MC locations are used to estimate the species concentrations in the plasma. Simulation using TORIC shows that with the 3He level <1%, most RF power is absorbed by the 3He ions and the process can generate energetic 3He ions. In mode conversion (MC) flow drive experiment in D(3He) plasma at 8 T, MC waves were also monitored by PCI. The MC ion cyclotron wave (ICW) amplitude and wavenumber kR have been found to correlate with the flow drive force. The MC efficiency, wave-number k of the MC ICW and their dependence on plasma parameters like Te0 have been studied. Based on the experimental observation and numerical study of the dispersion solutions, a hypothesis of the flow drive mechanism has been proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗