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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Effects of negative triangularity shaping on energetic particle driven Alfvén eigenmodes in DIII-D

Shape variations from positive to negative triangularity may lead to improved performance regimes in a tokamak with a reduction in microturbulence as well as improved thermal confinement. The motivation of this investigation is to model and compare the neutral beam driven Alfvén eigenmode (AE) activity in two DIII-D discharges with positive and negative triangularity shaping of plasma. The simulations are performed using the linear version of the FAR3d code, which solves the reduced MHD equations for thermal plasma with addition of moment equations for the energetic ion density and parallel velocity with appropriate Landau closure relations. Overall, our numerical results indicate that for similar physical parameters, the unstable AEs observed in the negative triangularity case have lower growth rates as compared to the positive triangularity regime. Our findings may be useful to analyze the influence of the reverse-D like geometry on the AE instabilities in DIII-D and may lead to better configurations for minimizing fast ion losses in a tokamak device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2001

Etching in Chlorine Discharges Using an Integrated Feature Evolution-Plasma Model

Etching of semiconductor materials is reliant on plasma properties. Quantities such as ion and neutral fluxes, both in magnitude and in direction, are often determined by reactor geometry (height, radius, position of the coils, etc.) In order to obtain accurate etching profiles, one must also model the plasma as a whole to obtain local fluxes and distributions. We have developed a set of three models that simulates C12 plasmas for etching of silicon, ion and neutral trajectories in the plasma, and feature profile evolution. We have found that the location of the peak in the ion densities in the reactor plays a major role in determining etching uniformity across the wafer. For a stove top coil inductively coupled plasma (ICP), the ion density is peaked at the top of the reactor. This leads to nearly uniform neutral and ion fluxes across the wafer. A side coil configuration causes the ion density to peak near the sidewalls. Ion fluxes are thus greater toward the wall's and decrease toward the center. In addition, the ions bombard the wafer at a slight angle. This angle is sufficient to cause slanted profiles, which is highly undesirable.

Hwang, Helen H.↗

DOE OSTI2017

Skin Effect Simulation for Area 11 Dense Plasma Focus Hot Plate

Two arc flashover events occurred at the DPF Area 11 facility. These flashover events happened in the same location on the bank current delivery plates. The damage from one of these events can be seen on the left-hand side of Figure 1. Since the flashovers occurred in the same area of the bank, and the reliability of the bank is important for future DPF experiments, a failure analysis effort was initiated. Part of this failure analysis effort was an effort to understand the physical reasons behind why the flashover happened, and why it happened in the same place twice. This paper summarizes an effort to simulate the current flow in the bank in order to understand the reasons for the flashover.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

A two-fluid study of oblique tearing modes in a force-free current sheet

Kinetic simulations have demonstrated that three-dimensional reconnection in collisionless regimes proceeds through the formation and interaction of magnetic flux ropes, which are generated due to the growth of tearing instabilities at multiple resonance surfaces. Since kinetic simulations are intrinsically expensive, it is desirable to explore the feasibility of reduced two-fluid models to capture this complex evolution, particularly, in the strong guide field regime, where two-fluid models are better justified. With this goal in mind, this paper compares the evolution of the collisionless tearing instability in a force-free current sheet with a two-fluid model and fully kinetic simulations. Our results indicate that the most unstable modes are oblique for guide fields larger than the reconnecting field, in agreement with the kinetic results. The standard two-fluid tearing theory is extended to address the tearing instability at oblique angles. The resulting theory yields a flat oblique spectrum and underestimates the growth of oblique modes in a similar manner to kinetic theory relative to kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

Current carrying properties of double layers and low frequency auroral fluctuations

A summary is given of a systematic parametric study in which the ion-transit time effects on double-layer fluctuations are investigated in a series of numerical simulations involving several lengths of finite extent plasmas. Typical values of ion drift velocities both smaller and larger than the ion-thermal speed on the high potential end of the simulation plasma are given. It is shown that the link between the electron-current fluctuations and ion-transit time is the Langmuir criterion for the existence of double layers. The time period of the recurring phenomena depends on the ion dynamics; for ions entering the simulation plasma in such a way that the ion drift velocity is less than approximately the ion thermal speed, the time period is governed by the ion transit time from the location of ion injection to the location of the double layer.

Singh, N.↗

DOE OSTI2020

Optimized continuum x-ray emission from laser-generated plasma

We study continuum x-ray emission from hot plasma at the National Ignition Facility (NIF). We find that the x-ray yield in the multi-keV photon energy range is larger in Ti than in Ag or Au. This apparent paradox is due to Ti K-shell vacancies generated by the extraordinary energy density achieved by the NIF lasers. This is supported by direct observations of large continuum enhancement above the Ti K-series limit due to both free–bound (recombination) emission and strong Ly α (H-like) emission. Detailed calculations agree well with our measurements and support our conclusions.  X-ray emission from matter pervades the universe and plays an important role in the high-energy density (HED) regime (~10 5 J/cm 3 ). For example, the emission and absorption of x-rays controls the thermal state and evolution of stellar interiors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

Plasma contactors for electrodynamic tether

The role plasma contactors play in effective electrodynamic tether operation is discussed. Hollow cathodes and hollow cathode-based plasma sources have been identified as leading candidates for the electrodynamic tether plasma contactor. Present experimental efforts to evaluate the suitability of these devices as plasma contactors, conducted concurrently at NASA Lewis Research Center and Colorado State University, are reviewed. These research programs include the definition of preliminary plasma contactor designs, and the characterization of their operation both as electron emitters and electron collectors to and from a simulated space plasma. Results indicate that ampere-level electron currents, sufficient for electrodynamic tether operation, can be exchanged between hollow cathode-based plasma contactors and a dilute plasma.

Patterson, Michael J.↗

DOE OSTI2025

Addressing key physics problems in high-energy-density plasmas with a novel kinetic simulation capability

Many important physical processes in inertial confinement fusion (ICF) and dense Z-pinch (DZP) experiments require a kinetic (velocity-space-dependent) description. Conventional particle-in-cell (PIC) methods are poorly suited for high-energy-density (HED) plasmas, due to restrictive time-step constraints and the inability to conserve energy. In a previous LDRD (21-FS-048), we demonstrated that a fully implicit PIC formulation overcomes these limitations: it conserves energy even when coupled with Coulomb collision models and can be solved efficiently with large grid cells and large time steps. Thus, it is feasible to use this method to study kinetic effects in ICF and DZP plasmas on hydro-like time and spatial scales. In this follow-on LDRD, we advanced this methodology into a high-fidelity tool for production-scale simulations and used it to answer key questions relevant to ICF and DZP experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Observation of laser ablation of silicon as a function of pulse length at constant fluence via time-resolved x-ray spectroscopy

We investigate the ablation of silicon as a function of laser pulse length at a constant fluence using time-resolved x-ray spectroscopy data obtained from OMEGA EP experiments at the University of Rochester's Laboratory for Laser Energetics. Our targets consisted of three-layer planar structures composed of Si (50 μm), Cu (25 μm), and SiO 2 (500 μm) layers. The Si layer was irradiated by a 351-nm laser with varying pulse widths of 250 ps, 500 ps, 1 ns, and 10 ns while maintaining a constant fluence of ~27.9 kJ/cm 2 . Electron temperatures and densities of the ablated plasma were determined by analyzing the time-resolved x-ray spectroscopy data through a comparison of experimental measurements with synthetic results obtained from Si atomic calculations in a steady state and non-local thermodynamic equilibrium. These calculations were computed using PrismSPECT. Additionally, radiation-hydrodynamics simulations with FLASH are used to generate simulated plasma-density and plasma-temperature profiles, which are then compared with the experimental measurements. Our analyses reveal that increasing the laser pulse length at a constant fluence results in a decrease in electron temperatures and densities. Furthermore, the longer pulses with lower intensities lead to deeper ablation regions before reaching the peak ablation but lower ionization balances in the silicon layer. Furthermore, these findings emphasize the critical role of laser pulse length in plasma ablation and shock generation for laser-impulse studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

Yang, Y. S.↗

DOE OSTI2018

Wavelet methods for studying the onset of strong plasma turbulence

Recent simulations have demonstrated that coherent current sheets dominate the kinetic-scale energy dissipation in strong turbulence of magnetized plasma. Wavelet basis functions are a natural tool for analyzing turbulent flows containing localized coherent structures of different spatial scales. Here, wavelets are used to study the onset and subsequent transition to fully developed turbulence from a laminar state. Originally applied to neutral fluid turbulence, an iterative wavelet technique decomposes the field into coherent and incoherent contributions. In contrast to Fourier power spectra, finite time Lyapunov exponents, and simple measures of intermittency such as non-Gaussian statistics of field increments, the wavelet technique is found to provide a quantitative measure for the onset of turbulence and to track the transition to fully developed turbulence. The wavelet method makes no assumptions about the structure of the coherent current sheets or the underlying plasma model. Temporal evolution of the coherent and incoherent wavelet fluctuations is found to be highly correlated (a Pearson correlation coefficient of >0.9) with the magnetic field energy and plasma thermal energy, respectively. The onset of turbulence is identified with the rapid growth of a background of incoherent fluctuations spreading across a range of scales and a corresponding drop in the coherent components. This is suggestive of the interpretation of the coherent and incoherent wavelet fluctuations as measures of coherent structures (e.g., current sheets) and dissipation, respectively. The ratio of the incoherent to coherent fluctuations Ric is found to be fairly uniform in the turbulent state across different plasma models and provides an empirical threshold of ~0.1 for turbulence onset. The utility of this technique is illustrated through examples. First, it is applied to the Kelvin–Helmholtz instability from different simulation models including fully kinetic, hybrid (kinetic ion/fluid electron), and Hall MHD simulations. Second, the wavelet diagnostic is applied to the development of turbulence downstream of the bowshock in a global magnetosphere simulation. Finally, the wavelet technique is also shown to be useful as a de-noising method for particle simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Nonlinear wavenumber shift of large amplitude Langmuir waves

Nonlinear particle-in-cell simulation is carried out to investigate the nonlinear behavior of the Langmuir wave launched with a fixed frequency in a uniform plasma. It is found that in the strong driving case, the launched wave propagates in a phase velocity larger than that predicted by the linear theory; there appears a nonlinear down-shift of wavenumber. The phase velocity of the nonlinear wave and the down-shift of the wavenumber are demonstrated to be determined by the velocity of nonlinearly accelerated resonant electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗

DOE OSTI2017

Nonlinear fishbone dynamics in spherical tokamaks

Linear and nonlinear kinetic-MHD hybrid simulations have been carried out to investigate linear stability and nonlinear dynamics of beam-driven fishbone instability in spherical tokamak plasmas. Realistic NSTX parameters with finite toroidal rotation were used. The results show that the fishbone is driven by both trapped and passing particles. The instability drive of passing particles is comparable to that of trapped particles in the linear regime. The effects of rotation are destabilizing and a new region of instability appears at higher q min (>1.5) values, q min being the minimum of safety factor profile. In the nonlinear regime, the mode saturates due to flattening of beam ion distribution, and this persists after initial saturation while mode frequency chirps down in such a way that the resonant trapped particles move out radially and keep in resonance with the mode. Correspondingly, the flattening region of beam ion distribution expands radially outward. A substantial fraction of initially non-resonant trapped particles become resonant around the time of mode saturation and keep in resonance with the mode as frequency chirps down. On the other hand, the fraction of resonant passing particles is significantly smaller than that of trapped particles. Our analysis shows that trapped particles provide the main drive to the mode in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Prediction of Breakdown in Air and Solid Dielectrics: A Complete Plasma Model from Discharge Initiation to Flashover (Final Technical Report)

A major challenge for most of the technologies of today, and the technologies of the future, is the compact and reliable design of electrical insulation systems for their electrical components. Although controlled low-temperature plasmas (LTPs) or non-thermal plasmas in air at atmospheric pressure are of interest for a wide range of technologies, the appearance of such discharges is considered highly undesirable in insulation systems because they are precursors of a breakdown or flashover, where complete insulation failure takes place. For non-thermal plasmas, often only one mechanism, or the transition from one discharge mechanism to another, has typically been studied to date. On the other hand, in the case of thermal plasma applications, the conditions following the occurrence of an electric arc or flashover have always been considered. Here, there is a big technical gap: there is no paper describing a complete plasma model simulating all discharge mechanisms and their transitions from discharge initiation to electrical arc and flashover. To address this gap, a complete plasma model for both thermal and non-thermal plasmas should be developed using a hydrodynamic approach. The project aimed to build the model. The work has the potential to address industrial applications and is of great importance to higher voltage systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Reconstruction and interpretation of ionization asymmetry in magnetic confinement via synthetic diagnostics

Abstract Strong poloidal refueling asymmetry in the DIII-D tokamak is inferred from line radiation measurements. Synthetic diagnostics in neutral transport modeling coupled to gyrokinetic simulations illuminate implications for the plasma flow profile in the scrape-off layer of single-null beam-driven discharges. Recycling occurs primarily either on the inner or outer divertor legs, depending on the toroidal magnetic field direction. By reversing the toroidal magnetic field, the observed line radiation asymmetry is nearly eliminated or reversed. It is determined that, while relatively simple physics can describe the observed ionization asymmetry, predicting the overall brightness of the hydrogenic Lyman- α signal requires detailed simulation of the plasma and resulting turbulence. To this end, kinetic plasma simulations fully coupled to comprehensive neutral transport calculations—a novel capability—provide first-principles reproduction of Lyman- α observations on DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Gyrokinetic studies of trapped electron mode turbulence in the Helically Symmetric eXperiment stellarator

Gyrokinetic simulations of plasma microturbulence in the Helically Symmetric eXperiment are presented. Using plasma profiles relevant to experimental operation, four dominant drift wave regimes are observed in the ion wavenumber range, which are identified as different flavors of density-gradient-driven trapped electron modes. For the most part, the heat transport exhibits properties associated with turbulence driven by these types of modes. Additionally, long-wavelength, radially localized, nonlinearly excited coherent structures near the resonant central flux surface, not predicted by linear simulations, can further enhance flux levels. Integrated heat fluxes are compatible with experimental observations in the corresponding density gradient range. Despite low shearing rates, zonal flows are observed to regulate turbulence but can be overwhelmed at higher density gradients by the long-wavelength coherent structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Monte-Carlo fluid approaches to detached plasmas in non-axisymmetric divertor configurations

Fluid transport modeling in three-dimensional boundaries of toroidal confinement devices is reviewed with the emphasis on a Monte-Carlo approach to simulate detached plasmas. The loss of axisymmetry in such configurations presents a major challenge for numerical implementation of the standard fluid model widely applied to fusion experimental devices. A large-scale effort has been made to address this problem under complementary aspects including different magnetic topologies and numerical techniques. Here in this paper, we give a brief review of the different strategies pioneered and the challenges involved. A more detailed description is provided for the Monte-Carlo code—EMC3-Eirene, where the physics model and the basic idea behind the applied Monte-Carlo method are presented. The focus is put on its applications to detachment studies for stellarators and tokamaks. Here, major achievements and difficulties encountered are described. Model limitations and further development plans are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗