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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Validation of nonlinear gyrokinetic simulations of L- and I-mode plasmas on Alcator C-Mod

New validation of global, nonlinear, ion-scale gyrokinetic simulations (GYRO) is carried out for L- and I-mode plasmas on Alcator C-Mod, utilizing heat fluxes, profile stiffness, and temperature fluctuations. Previous work at C-Mod found that ITG/TEM-scale GYRO simulations can match both electron and ion heat fluxes within error bars in I-mode [White PoP 2015], suggesting that multi-scale (cross-scale coupling) effects [Howard PoP 2016] may be less important in I-mode than in L-mode. New results presented here, however, show that global, nonlinear, ion-scale GYRO simulations are able to match the experimental ion heat flux, but underpredict electron heat flux (at most radii), electron temperature fluctuations, and perturbative thermal diffusivity in both L- and I-mode. Linear addition of electron heat flux from electron scale runs does not resolve this discrepancy. These results indicate that single-scale simulations do not sufficiently describe the I-mode core transport, and that multi-scale (coupled electron- and ion-scale) transport models are needed. A preliminary investigation with multi-scale TGLF, however, was unable to resolve the discrepancy between ion-scale GYRO and experimental electron heat fluxes and perturbative diffusivity, motivating further work with multi-scale GYRO simulations and a more comprehensive study with multi-scale TGLF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

1D EDIPIC

EDIPIC codes were developed for simulation of kinetic effects in low-temperature plasmas. They can be used for simulation of plasmas in electrostatic approximation and include many atomistic and plasma-surface interaction processes needed for simulations of low temperature plasmas. These codes include a detailed model of secondary electron emission induced by electrons and ions. The Monte-Carlo model of electron-neutral collisions and the Langevin model of Coulomb collisions for electrons are also implemented.

Sydorenko, Dmytro↗
DOE OSTI2021

EDIPIC-2D

EDIPIC codes were developed for simulation of kinetic effects in low-temperature plasmas. They can be used for simulation of plasmas in electrostatic approximation and include many atomistic and plasma-surface interaction processes needed for simulations of low temperature plasmas. These codes include a detailed model of secondary electron emission induced by electrons and ions. The Monte-Carlo model of electron-neutral collisions and the Langevin model of Coulomb collisions for electrons are also implemented.

Sydorenko, Dmytro↗
DOE OSTI2024

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Enhanced betatron X-rays from axially modulated plasma wakefields

In the cavitation regime of plasma-based accelerators, a population of high-energy electrons trailing the driver can undergo betatron motion. The motion results in X-ray emission, but the brilliance and photon energy are limited by the electrons' initial transverse coordinate. To overcome this, we exploit parametrically unstable betatron motion in a cavitated, axially modulated plasma. Theory and simulations are presented showing that the unstable oscillations increase both the total X-ray energy and average photon energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Studies on counterstreaming plasma expansion

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

Singh, N.↗
DOE OSTI2017

Excitation of low frequency Alfven eigenmodes in toroidal plasmas

Global gyrokinetic simulations find that realistic density gradients of energetic particles can simultaneously excite low frequency Alfven eigenmodes in toroidal geometry, beta-induced Alfven-acoustic eigenmode (BAAE) and beta-induced Alfven eigenmode (BAE), with similar radial mode widths and comparable linear growth rates even though damping rate of BAAE is much larger than BAE in the absence of energetic particles. Here, this surprising result is attributed to non-perturbative effects of energetic particles that modify ideal BAAE mode polarizations and nonlocal geometry effects that invalidate radially local dispersion relation. Dominant mode changes from BAAE in a larger tokamak to BAE in a smaller tokamak due to the dependence of wave-particle resonance condition on the tokamak size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Study of the Interaction of Plasma with Beryllium That Is a Candidate Material for the First Wall of a Fusion Reactor

The experimental conditions for treatment of a candidate material of the first wall, beryllium, by the plasma have been determined. The plasma parameters obtained using the simulator with a plasma-beam setup during the irradiation of samples as well as the results of material studies obtained by scanning electron microscopy and energy-dispersion analysis of the beryllium surface layer have been considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Benchmark studies of the gyro-Landau-fluid code and gyro-kinetic codes on kinetic ballooning modes

A Gyro-Landau-Fluid (GLF) 3 + 1 model has been recently implemented in BOUT++ framework, which contains full Finite-Larmor-Radius effects, Landau damping, and toroidal resonance [Ma et al., Phys. Plasmas 22, 055903 (2015)]. A linear global beta scan has been conducted using the JET-like circular equilibria (cbm18 series), showing that the unstable modes are kinetic ballooning modes (KBMs). In this work, we use the GYRO code, which is a gyrokinetic continuum code widely used for simulation of the plasma microturbulence, to benchmark with GLF 3 + 1 code on KBMs. To verify our code on the KBM case, we first perform the beta scan based on “Cyclone base case parameter set.” We find that the growth rate is almost the same for two codes, and the KBM mode is further destabilized as beta increases. For JET-like global circular equilibria, as the modes localize in peak pressure gradient region, a linear local beta scan using the same set of equilibria has been performed at this position for comparison. With the drift kinetic electron module in the GYRO code by including small electron-electron collision to damp electron modes, GYRO generated mode structures and parity suggest that they are kinetic ballooning modes, and the growth rate is comparable to the GLF results. However, a radial scan of the pedestal for a particular set of cbm18 equilibria, using GYRO code, shows different trends for the low-n and high-n modes. The low-n modes show that the linear growth rate peaks at peak pressure gradient position as GLF results. However, for high-n modes, the growth rate of the most unstable mode shifts outward to the bottom of pedestal and the real frequency of what was originally the KBMs in ion diamagnetic drift direction steadily approaches and crosses over to the electron diamagnetic drift direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Nonlinear evolution of stimulated Raman scattering near the quarter-critical density

Nonlinear evolution of stimulated Raman scattering (SRS) near the quarter-critical density is studied using one-dimensional (1D) and two-dimensional (2D) particle-in-cell simulations in homogeneous plasmas. In 1D configuration, with two-plasmon decay (TPD) instability excluded, the system evolves into two regimes distinguished by whether density cavities have been formed or not. At low temperatures, a cavity regime characterised by high absorption and low reflection, and at high temperatures nonlinear frequency shift regime due to particle trapping, are observed. Furthermore, a competition between SRS and TPD in 2D simulations reveals that the nonlinear SRS does play a significant role near the quarter-critical density, whose influences were mostly neglected before.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Blob Formation and Ejection in Coronal Jets due to the Plasmoid and Kelvin–Helmholtz Instabilities

We perform 2D resistive magnetohydrodynamic simulations of coronal jets driven by flux emergence along the lower boundary. The reconnection layers are susceptible to the formation of blobs that are ejected in the jet. Our simulation with low plasma β (Case I) shows that magnetic islands form easily and propagate upward in the jet. These islands are multithermal and thus are predicted to show up in hot channels (335 Å and 211 Å) and the cool channel (304 Å) in observations by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory . The islands have maximum temperatures of 8 MK, lifetimes of 120 s, diameters of 6 Mm, and velocities of 200 km s{sup −1}. These parameters are similar to the properties of blobs observed in extreme-ultraviolet (EUV) jets by AIA. The Kelvin–Helmholtz instability develops in our simulation with moderately high plasma β (Case II) and leads to the formation of bright vortex-like blobs above the multiple high magnetosonic Mach number regions that appear along the jet. These vortex-like blobs can also be identified in the AIA channels. However, they eventually move downward and disappear after the high magnetosonic Mach number regions disappear. In the lower plasma β case, the lifetime for the jet is shorter, the jet and magnetic islands are formed with higher velocities and temperatures, the current-sheet fragments are more chaotic, and more magnetic islands are generated. Our results show that the plasmoid instability and Kelvin–Helmholtz instability along the jet are both possible causes of the formation of blobs observed at EUV wavelengths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

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

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

energetic particles↗
DOE OSTI2022

Modeling deformation and failure in AlSi-polyester abradable sealcoating material using microstructure-based finite element simulation

A plasma-sprayed aluminum-silicon (AlSi)/polyester coating is applied in modern gas turbine engines as an abradable sealcoating to maintain tight clearances between the rotating blades and the static casing. While running the engine, the rotating blades “rub” with the abradable coating, which results in extreme strain rate (up to 10 6 ) dynamics and a high-temperature environment. Due to the difficulty of collecting direct measurements, predictive computational models are important for analyzing the deformation and failure of the abradable material, to help meet the design target of avoiding damage to the blade tip and maintaining high fuel efficiency. In this research, a microstructure-based finite element (FE) computational model was developed to capture the complex mechanical behavior of the AlSi/polyester microstructure. The model is based on a virtual representative-volume-element (RVE) of a metal-polymer microstructure, reconstructed from x-ray computed tomography. It models the plastic deformation of, and damage to, each AlSi and polyester constituents, as well as the failure at their interface. The model was calibrated and validated with uniaxial tension and compression experiments, conducted at two temperatures (298 K and 533 K) at an applied strain rate of 10 3 - 10 4 s -1 . The material exhibited strongly asymmetric tension-compression behavior and a sensitivity to temperature, which was well captured by the model. The model was further applied to investigate changes in mechanical behavior due to variations in constituents’ volume fractions, which provides guidance to the microstructural design of AlSi/polyester abradable materials. The model is expected to facilitate the development of improved abradable materials by bypassing the conventional trial-and-error approach and extensive testing requirements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Ion beam emittance from an ECRIS

Simulation of ion beam extraction from an Electron Cyclotron Resonance Ion Source (ECRIS) is a fully 3 dimensional problem, even if the extraction geometry has cylindrical symmetry. Because of the strong magnetic flux density, not only the electrons are magnetized but also the Larmor radius of ions is much smaller than the geometrical dimension of the plasma chamber (Ø 64 × 179 mm). If we assume that the influence of collisions is small on the path of particles, we can do particle tracking through the plasma if the initial coordinates of particles are known. We generated starting coordinates of plasma ions by simulation of the plasma electrons, accelerated stochastically by the 14.5 GHz radio frequency power fed to the plasma. With that we were able to investigate the influence of different electron energies on the extracted beam. Using these assumptions, we can reproduce the experimental results obtained 10 years ago, where we monitored the beam profile with the help of viewing targets. Additionally, methods have been developed to investigate arbitrary 2D cuts of the 6D phase space. To this date, we are able to discuss full 4D information. Currently, we extend our analysis tool towards 5D and 6D, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1992

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

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

Delva, Magda↗
NASA NTRS1981

Numerical simulation of the generation of electrostatic turbulence in the magnetotail

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

Swift, D. W.↗