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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Numerical simulation to study the transient self focusing of laser beam in plasma

In this paper, we present the numerical simulation for the coupled system of equations governing the dynamics of laser and Ion Acoustic Wave (IAW) in a collisionless plasma, when the coupling between the waves is through ponderomotive non-linearity. The nonlinear evolution of the laser beam is studied when the pump laser is perturbed by a periodic perturbation. By changing the perturbation wave number, we have studied its effect on the nonlinear evolution pattern of laser beam. In order to have a physical insight into the nonlinear dynamics of laser beam evolution in time and space, we have studied the laser and IAW spectra containing spatial harmonics. The magnitude of these harmonics changes with time and leads to time dependent localization of laser beam in spatial domain. The nonlinear dynamics of this localization is investigated in detail by using simulation and a semi-analytical model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗
DOE OSTI2023

Photon Doppler Velocimetry to Spatially Resolve Plasma Density in a Power Flow Gap

The understanding of power flow plasmas is important as we look towards next generation pulsed power (NGPP) as current losses could prohibit the goals of that facility. Therefore, it is important to have accurate diagnostics of the plasma parameters on the current machines, which can be used to help inform and improve simulations. Having these plasma parameters will help validate models and simulations to provide confidence when they are expanded to conditions relevant to NGPP. One important plasma parameter that can be measured is the electron density, which can be measured by photonic Doppler velocimetry (PDV). A PDV system has several key advantages over other interferometers by measuring relatively low densities (> 1 × 10 15 cm -2 ) with both spatial and temporal resolution. Experiments were performed on the Mykonos pulsed power machine, which is a 1 MA sub scale machine in which recent platforms have been developed to explore current densities relevant to the inner magnetically insulated transmission line (MITL) on the Z machine. Experiments were performed on two different platforms, the thin foil platform and the Mykonos parallel plate platform (MP3). In addition, a combination of both single-point and multi-point measurements were used. The single-point measurements proved to be very promising, providing a clear increase in density at about 70 ns into the current rise on thin foil experiments up to about 5 × 10 17 cm -3 before the probe stopped providing signal. While we did also see returns from multi-point measurements on both platforms, the signals were not as easy to interpret due to strong background effects. However, they do show initial promise for this diagnostic to measure density at several points across a 1 mm gap. These measurements provide insights in how to improve the diagnostic so that it can provide useful information on power flow relevant experiments.

47 OTHER INSTRUMENTATION↗
NASA NTRS1988

Simulation of electrostatic turbulence in the plasma sheet boundary layer with electron currents and bean-shaped ion beams

Plasma data from ISEE-1 show the presence of electron currents as well as energetic ion beams in the plasma sheet boundary layer. Broadband electrostatic noise and low-frequency electromagnetic bursts are detected in the plasma sheet boundary layer, especially in the presence of strong ion flows, currents, and steep spacial gradients in the fluxes of few-keV electrons and ions. Particle simulations have been performed to investigate electrostatic turbulence driven by a cold electron beam and/or ion beams with a bean-shaped velocity distribution. The simulation results show that the counterstreaming ion beams as well as the counterstreaming of the cold electron beam and the ion beam excite ion acoustic waves with a given Doppler-shifted real frequency. However, the effect of the bean-shaped ion velocity distributions reduces the growth rates of ion acoustic instability. The simulation results also show that the slowing down of the ion bean is larger at the larger perpendicular velocity. The wave spectra of the electric fields at some points of the simulations show turbulence generated by growing waves.

Nishikawa, K.-I.↗
DOE OSTI2016

Efficient production of high-energy nonthermal particles during magnetic reconnection in a magnetically dominated ion-electron plasma

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux-dominated flows. We investigate nonthermal particle acceleration during magnetic reconnection in a magnetically dominated ion–electron plasma using fully kinetic simulations. For an ion–electron plasma with a total magnetization of ${\sigma }_{0}={B}^{2}/(4\pi n({m}_{i}+{m}_{e}){c}^{2})$, the magnetization for each species is ${\sigma }_{i}\sim {\sigma }_{0}$ and ${\sigma }_{e}\sim ({m}_{i}/{m}_{e}){\sigma }_{0}$, respectively. We have studied the magnetically dominated regime by varying σe = 103–105 with initial ion and electron temperatures ${T}_{i}={T}_{e}=5-20{m}_{e}{c}^{2}$ and mass ratio ${m}_{i}/{m}_{e}=1-1836$. Our results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several (2–3) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1\lt s\lt 2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be ${\gamma }_{{be}}\sim {\sigma }_{e}$ and ${\gamma }_{{bi}}\sim {\sigma }_{i}$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices sp are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. Finally, we discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds, especially the applications for explaining nonthermal high-energy emissions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Impurity transport, turbulence transitions and intrinsic rotation in Alcator C-Mod plasmas

Linear and nonlinear gyrokinetic simulations are used to probe turbulent impurity transport in intrinsically rotating tokamak plasmas. For this simulation-based study, experimental input parameters are taken from a pair of ICRF heated Alcator C-Mod discharges exhibiting a change in the sign of the normalized toroidal rotation gradient at mid-radius (i.e. a change from hollow to peaked intrinsic rotation profiles). The simulations show that there is no change in the peaking of the calcium impurity between the plasmas with peaked and hollow rotation profiles, suggesting that the impurity transport and the shape of the rotation do not always change together. Furthermore, near mid-radius, r/a = 0.5 (normalized midplane minor radius), the linear and nonlinear gyrokinetic simulations exhibit no evidence of a transition from ion temperature gradient (ITG) to trapped electron mode dominance when the intrinsic rotation profile changes from peaked to hollow. Finally, extensive nonlinear sensitivity analysis is performed, and there is no change in the ITG critical gradient or in the stiffness of ion heat transport with the change in the intrinsic toroidal rotation profile shape, which suggests that the shape of the rotation profile is not dominated by the ITG onset in these cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Nitrogen mass transfer models for plasma-based low-energy ion implantation

The nitrogen mass transfer process in plasma-based low-energy ion implantation (PBLEII) is theoretically and experimentally studied in order to explore the process mechanism of PBLEII and therefore to optimize the apparatus design and the process conditions. An electron cyclotron resonance (ECR) microwave discharge generates the nitrogen plasma with a high density of 10{sup 11}–10{sup 12} ions/cm{sup 3}, which diffuses downstream to the process chamber along the divergent magnetic field. The nitrogen ions in the plasma implant into the surface and transport to the matrix of an austenitic stainless steel under the low negative pulsed bias of −2 kV at a process temperature of 400 °C. A global plasma model is used to simulate the ECR microwave plasma discharge for a range of working pressures and microwave powers. The fluid models are adopted to calculate the plasma downstream diffusion, the sheath expansion and the low-energy ion implantation on the surface. A nonlinear kinetic discrete model is established to describe the nitrogen transport in the austenitic stainless steel and the results are compared with the experimental measurements. Under an average implantation current density of 0.3–0.6 mA/cm{sup 2}, the surface nitrogen concentration in the range from 18.5 to 29 at. % is a critical factor for the nitrogen transport in the AISI 304 austenitic stainless steel by PBLEII, which accelerates the implanted nitrogen diffusion inward up to 6–12 μm during a nitriding time of 4 h.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Numerical simulation of microwave amplification in a plasma channel produced in a gas via multiphoton ionisation by a femtosecond laser pulse

This paper examines the evolution of a nonequilibrium plasma channel produced in xenon by a femtosecond KrF laser pulse. We demonstrate that such a channel can be used to amplify microwave pulses over times of the order of the relaxation time of the photoelectron energy spectrum in xenon. Using the slowly varying amplitude approximation, we analyse the propagation and amplification of an rf pulse in a plasma channel, in particular when the rf field influences the electron energy distribution function in the plasma. (interaction of laser radiation with matter. laser plasma)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

High-power laser interaction with low-density C–Cu foams

Here, we study the propagation of high-power laser beams in micro-structured carbon foams by monitoring the x-ray output from deliberately introduced Cu content. In particular, we characterize this phenomenon measuring absolute time-resolved x-ray yields, time-resolved x-ray imaging, and x-ray spectroscopy. New experimental results for C–Cu foams show a faster heat front velocity than simulation that assumed homogeneous plasma. We suggest the foam micro-structure may explain this trend.

42 ENGINEERING↗
NASA NTRS2008

Plasma Sheet Circulation Pathways

Global simulations of Earth's magnetosphere in the solar wind compute the pathways of plasma circulation through the plasma sheet. We address the pathways that supply and drain the plasma sheet, by coupling single fluid simulations with Global Ion Kinetic simulations of the outer magnetosphere and the Comprehensive Ring Current Model of the inner magnetosphere, including plasmaspheric plasmas. We find that the plasma sheet is supplied with solar wind plasmas via the magnetospheric flanks, and that this supply is most effective for northward IMF. For southward IMF, the innermost plasma sheet and ring current region are directly supplied from the flanks, with an asymmetry of single particle entry favoring the dawn flank. The central plasma sheet (near midnight) is supplied, as expected, from the lobes and polar cusps, but the near-Earth supply consists mainly of slowly moving ionospheric outflows for typical conditions. Work with the recently developed multi-fluid LFM simulation shows transport via plasma "fingers" extending Earthward from the flanks, suggestive of an interchange instability. We investigate this with solar wind ion trajectories, seeking to understand the fingering mechanisms and effects on transport rates.

Moore, Thomas E.↗
DOE OSTI2022

Full radiator-perturber interaction in computer simulations of hydrogenic spectral line broadening by plasmas

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating potential obtained from computer simulations. Such calculations have relied on the multipole expansion of the radiator-plasma interaction often keeping only the dipole term. With increasing density, however, higher multipoles as well as plasma perturbers overlapping the bound electron wave functions are expected to become important. For hydrogenic systems, the atomic matrix elements of the full Coulomb and screened Coulomb interactions are given by analytical formulas. Using these results, a computer simulation approach that accounts for the full radiator-plasma interaction is developed. One benefit is the removal of inherent strong collision divergences in the multipole expansion approximation. Furthermore, it yields the plasma polarization shift produced by perturbers penetrating the wave function of the radiator bound electrons. Here, the model was applied to hydrogenic argon Ly-α , Ly- β , and Ly-$\gamma$ spectral lines in a dense argon plasma at free electron densities of 10 24 or 10 25 cm ⁻3 and temperature of $800$ $eV$ relevant to plasma diagnostic techniques for inertial confinement fusion implosions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

High ionisation fraction plasmas in a low temperature, multidipole cusp plasma

The depletion of neutral helium atoms has been studied in an unmagnetised spherical plasma created by DC discharge in a multidipole confinement field. Knowing the neutral density profile is critical to predicting the equilibrium flow of such plasmas. A model of the emissivity due to electron-impact excitation of neutral atoms in the plasma has been derived and used to fit radiance measurements of several neutral transitions to extract the radial profile of neutral density for plasmas of varying temperature and density. We report a depletion of the core neutral density varying between negligible levels to 80 % of the edge neutral density depending on the input power and fuelling. The corresponding ionisation fraction varies between 30–80 % in the plasma core. A simple neutral diffusion model is sufficient to describe the shape of neutral density profile implied by the radiance measurements. We have used the measurements to include a drag force due to neutral charge-exchange collisions in simulations of driven plasma flow. The simulation predicts a better fit to Mach probe flow measurements when this neutral drag is accounted for. This study shows that accounting for a realistic neutral profile is important to predict the plasma flow geometry and its magnetohydrodynamics (MHD) stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

ELECTRON HEATING IN A RELATIVISTIC, WEIBEL-UNSTABLE PLASMA

The dynamics of two initially unmagnetized relativistic counter-streaming homogeneous ion–electron plasma beams are simulated in two dimensions (2D) using the particle-in-cell (PIC) method. It is shown that current filaments, which form due to the Weibel instability, develop a large-scale longitudinal electric field in the direction opposite to the current carried by the filaments as predicted by theory. This field, which is partially inductive and partially electrostatic, is identified as the main source of net electron acceleration, greatly exceeding that due to magnetic field decay at later stages. The transverse electric field, although larger than the longitudinal field, is shown to play a smaller role in heating electrons, contrary to previous claims. It is found that in one dimension, the electrons become strongly magnetized and are not accelerated beyond their initial kinetic energy. Rather, the heating of the electrons is enhanced by the bending and break up of the filaments, which releases electrons that would otherwise be trapped within a single filament and slow the development of the Weibel instability (i.e., the magnetic field growth) via induction as per Lenz’s law. In 2D simulations, electrons are heated to about one quarter of the initial kinetic energy of ions. The magnetic energy at maximum is about 4%, decaying to less than 1% by the end of the simulation. The ions are found to gradually decelerate until the end of the simulation, by which time they retain a residual anisotropy of less than 10%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Laser-generated plasmas in length scales relevant for thin film growth and processing: simulation and experiment

In pulsed laser deposition, thin film growth is mediated by a laser-generated plasma, whose properties are critical for controlling the film microstructure. The advent of 2D materials has renewed the interest in how this ablation plasma can be used to manipulate the growth and processing of atomically thin systems. For such purpose, a quantitative understanding of the density, charge state, and kinetic energy of plasma constituents is needed at the location where they contribute to materials processes. Here, we study laser-induced plasmas over expansion distances of several centimeters from the ablation target, which is the relevant length scale for materials growth and modification. The study is enabled by a fast implementation of a laser ablation/plasma expansion model using an adaptive Cartesian mesh solver. Simulation outcomes for KrF excimer laser ablation of Cu are compared with Langmuir probe and optical emission spectroscopy measurements. Furthermore, simulation predictions for the plasma-shielding threshold, the ionization state of species in the plasma, and the kinetic energy of ions, are in good correspondence with experimental data. For laser fluences of 1–4 J cm -2 , the plume is dominated by Cu 0 , with small concentrations of Cu+ and electrons at the expansion front. Higher laser fluences (e.g. 7 J cm -2 ) lead to a Cu + -rich plasma, with a fully ionized leading edge where Cu 2+ is the dominant species. In both regimes, simulations indicate the presence of a low-density, high-temperature plasma expansion front with a high degree of ionization that may play a significant role in doping, annealing, and kinetically-driven phase transformations in 2D materials.

2D materials↗
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
NASA NTRS1974

Thermal relaxation of a two dimensional plasma in a dc magnetic field. Part 2: Numerical simulation

The thermal relaxation process for a spatially uniform two dimensional plasma in a uniform dc magnetic field is simulated numerically. Thermal relaxation times are defined in terms of the time necessary for the numerically computer Boltzman H-function to decrease through a given part of the distance to its minimum value. Dependence of relaxation time on two parameters is studied: number of particles per Debye square and ratio of gyrofrequency to plasma frequency.

Hsu, J. Y.↗
NASA NTRS1983

Electromagnetic radiation from beam-plasma instabilities

A computer simulation is developed for the generation of electromagnetic radiation in an electron beam-plasma interaction. The plasma is treated as a two-dimensional finite system, and effects of a continuous nonrelativistic beam input are accounted for. Three momentum and three field components are included in the simulation, and an external magnetic field is excluded. EM radiation generation is possible through interaction among Langmuir oscillations, ion-acoustic waves, and the electromagnetic wave, producing radiation perpendicular to the beam. The radiation is located near the plasma frequency, and polarized with the E component parallel to the beam. The scattering of Langmuir waves caused by ion-acoustic fluctuations generates the radiation. Comparison with laboratory data for the three-wave interactions shows good agreement in terms of the radiation levels produced, which are small relative to the plasma thermal energy.

Pritchett, P. L.↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗