Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Role of ion temperature on scrape-off layer plasma turbulence

Turbulence in Scrape-off layer (SOL) of tokamak plasma has been studied numerically using interchange modes with the help of electron continuity, quasineutrality, and ion energy equations. Electron temperature is assumed uniform. We have studied dynamics of seeded plasma blob and plasma turbulence to identify the role of ion temperature and its gradient. The ion temperature elongates the blob poloidally and reduces its radial velocity. Initial dipole nature of the plasma blob potential breaks and generates few more dipoles during its propagation in the SOL. Plasma turbulence simulation shows poloidally elongated density and ion temperature structures that are similar to the seeded blob simulation studies. Fluctuations of the density and ion temperature have been presented as function of scale lengths of the density and ion temperature. Reduction of the SOL width and increase of radial electric field have been measured in the presence of the ion temperature. Particle and energy transports have been also presented as the function of the density and ion temperature scale lengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

HIS-TORIC: Extending core ICRF wave simulation to include realistic SOL plasmas

This paper presents a novel approach to incorporating an arbitrarily shaped edge scrape-off-layer (SOL) plasma and an ion cyclotron range of frequency (ICRF) antenna structure into existing core ICRF wave simulation models. We partition the entire computation domain into two sub-domains: a core and an edge region. Simulations in each domain are performed separately with appropriate numerical solvers. For the core, the TORIC ICRF solver [M. Brambilla, Plasma Phys. Contrl. Fusion 41, 1 (1999)] was modified to impose an essential (Dirichlet) boundary condition at its interface with the edge domain. In the edge, a finite element method (FEM) is used to solve a cold collisional plasma model. The domains are then joined together using the continuity boundary condition for the tangential electric and magnetic fields at their interfaces (Hybrid Integration of SOL to TORIC: HIS-TORIC). The model developed here was tested using an ICRH H minority heating scenario on the Alcator C-Mod tokamak [I. Hutchinson et. al., Phys. Plasmas 1, 1511 (1994)]. The simulated pattern of core wave propagation agrees well with a standard TORIC simulation. This approach opens the possibility of using a realistic diverted SOL plasma and a complicated 3D RF antenna together with a rigorous hot core plasma model, while requiring only minimal modification to existing RF codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Self-pinched transport of a high ν/γ electron beam

The self-pinched transport of a 0.5 MeV, 18 kA cylindrical electron beam has been studied experimentally and computationally. The relatively low voltage and high current required for materials surface modification applications leads to complicated beam dynamics as the Alfven limit is approached. Transport and focusing of the high ν/γ beam was done in a sub-Torr, neutral gas-filled, conducting tube in the ion-focused regime. In this regime, beam space charge forces are progressively neutralized to allow focusing of the beam by its self-magnetic field. The beam exhibits stable envelope oscillations as it is efficiently and reproducibly propagated for distances greater than a betatron wavelength. Experimental results follow the trends seen in 2-D particle-in-cell simulations. Results show that the input electron beam can be periodically focused to a peaked profile with the beam half-current radius decreased by a factor of 2.84. This results in an increase of a factor of 8 in beam current density. This focusing is sufficient to produce desired effects in the surface layers of metallic materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. I. Nonlinear and kinetic effects

Results are presented from a series of simulations examining the susceptibility of the cross-beam energy transfer (CBET) instability to nonlinear processes in the context of direct-drive inertial confinement fusion experiments on the OMEGA laser facility. These form the basis for the second paper of this series [A. G. Seaton, L. Yin, R. Follett, B. J. Albright, and A. Le, “Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth,” Phys. Plasmas 29, 042707 (2022)], where we examine the efficacy of increases in laser bandwidth at suppressing CBET. We choose laser and plasma conditions for the simulations that are favorable to CBET and promote nonlinearity. Through a comparison of outputs from the particle-in-cell code vector particle in cell (VPIC) and the linearized fluid code laser-plasma simulation environment (LPSE), a series of nonlinear effects have been identified in the kinetic simulations that include particle trapping, the two-ion wave decay, and ion-acoustic wave self-focusing. These effects produce time-dependent energy transfer, in contrast to the linearized fluid simulations in which a steady state is reached after an initial transient. Ion trapping is shown to allow for increased energy transfer relative to fluid simulations, with the remaining nonlinear processes acting to reduce the energy transfer. Nonlinear dynamics is contrasted for low- and high-intensity beams as well as between speckled and planar beams. For the parameters under consideration, beam profile has a significant effect on nonlinear dynamics, though the greatest sensitivity is to beam intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Simulation of laser interaction with ablative plasma and hydrodynamic behavior of laser supported plasma

When an intense laser beam irradiates on a solid target, ambient air ionizes and becomes plasma, while part of the target rises in temperature, melts, vaporizes, ionizes, and yet becomes plasma. A general Godunov finite difference scheme WENO (Weighted Essentially Non-Oscillatory Scheme) with fifth-order accuracy is used to simulate 2-dimensional axis symmetrical laser-supported plasma flow field in the process of laser ablation. The model of the calculation of ionization degree of plasma and the interaction between laser beam and plasma are considered in the simulation. The numerical simulations obtain the profiles of temperature, density, and velocity at different times which show the evolvement of the ablative plasma. The simulated results show that the laser energy is strongly absorbed by plasma on target surface and that the velocity of laser supported detonation (LSD) wave is half of the ideal LSD value derived from Chapman-Jouguet detonation theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Towards fast and accurate predictions of radio frequency power deposition and current profile via data-driven modelling: applications to lower hybrid current drive

Three machine learning techniques (multilayer perceptron, random forest and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modelling and real-time control applications. The machine learning models use a database of more than 16 000 GENRAY/CQL3D simulations for training, validation and testing. Latin hypercube sampling methods ensure that the database covers the range of nine input parameters ( $n_{e0}$ , $T_{e0}$ , $I_p$ , $B_t$ , $R_0$ , $n_{\|}$ , $Z_{{\rm eff}}$ , $V_{{\rm loop}}$ and $P_{{\rm LHCD}}$ ) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to $\sim$ ms with high accuracy across the input parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Eigenmode characteristics of the double tearing mode in the presence of shear flows

The double tearing mode (DTM) is characterized by two eigen states with antisymmetric or symmetric magnetic island structure, referred to as the even or odd DTM. In this work, we systematically revisit the DTM instabilities in the presence of an antisymmetric shear flow with a focus on eigenmode characteristics as well as the stabilization or destabilization mechanism in a wide parameter region. Both initial value simulation and eigenvalue analysis are performed based on reduced resistive MHD model in slab geometry. A degenerated eigen state is found at a critical flow amplitude v{sub c}. The even (or odd) DTM is stabilized (or destabilized) by weak shear flow below v{sub c} through the distortion of magnetic islands mainly due to the global effect of shear flow rather than the local flow shear. The distortion can be quantified by the phase angles of the perturbed flux, showing a perfect correspondence to the growth rates. As the shear flow increases above v{sub c}, the degenerated eigen state bifurcates into two eigen modes with the same growth rate but opposite propagating direction, resulting in an oscillatory growth of fluctuation energy. It is identified that two eigen modes show the single tearing mode structure due to the Alfven resonance (AR) occurring on one current sheet. Most importantly, the AR can destabilize the DTMs through enhancing the plasma flow exerting on the remaining island. Meanwhile, the local flow shear plays a remarkable stabilizing role in this region. In addition, the eigenmode characteristic of the electromagnetic Kelvin-Helmholtz instability is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

A numerical study of the effect of various reactions, pressure and gas mixture ratio on the density distribution of etchant species (H, Br, Br{sup +}, and HBr{sup +}) in HBr/He plasma

In this study, a fluid model has been used to study the effect of gas mixing ratio and pressure on the density distribution of important etchant species, i.e., hydrogen (H), bromine (Br), Br{sup +}, and HBr{sup +} in HBr/He plasma. Our simulation results show that the densities of active etchant species H, Br, and HBr{sup +} increase with the increase in pressure as well as the HBr fraction in HBr/He mixture. On the contrary, the density of Br{sup +} decreases with the increase in He percentage in HBr/He mixture and with the increase in the pressure. Time averaged reaction rates (of the reactions involved in the production and consumption of these species) have been calculated to study the effect of these reactions on the density distribution of these species. The spatial distribution of these species is explained with the help of the time averaged reaction rates. Important reactions have been identified that contribute considerably to the production and consumption of these active species. The code has been optimized by identifying 26 reactions (out of 40 reactions which contribute in the production and consumption of these species) that have insignificant effect on the densities of H, Br, Br{sup +}, and HBr{sup +}. This shows that out of 40 reactions, only 14 reactions can be used to calculate the density and distribution of the important species in HBr/He plasma discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Nonlinear dynamics of toroidal Alfvén eigenmodes driven by energetic particles

A kinetic simulation code based on a reduced model is developed to study dynamic evolutions of a single toroidicity-induced shear Alfvén eigenmode driven by energetic particles. For zero background damping, it is found that the wave amplitude in nonlinear phase can either saturate for weak energetic particle drives or slowly increase for strong drives. This slow nonlinear growth in strong drive cases is found to be associated with broadening and overlapping of resonances between the wave and trapped particles. For the near-marginal-stability case with a large background damping, the mode nonlinear evolution exhibits strong upward and downward frequency chirping in multiple branches. A hole/clump formation is observed clearly in the corresponding evolution of energetic particle distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Linear simulation of kinetic electromagnetic instabilities in a tokamak plasma with weak magnetic shear

Gyrokinetic simulation and eigenvalue calculation of electromagnetic instabilities are carried out for an experimentally observed low-n mode in weak magnetic shear discharge. With different magnetic shear values, the Gyrokinetic Toroidal Code simulation of the ion temperature gradient mode is consistent with the eigenvalue code (HD7) calculation. Due to the sensitivity of the kinetic ballooning mode (KBM) to global equilibrium, the simulation of the KBM deviates from the eigenvalue results, for the ballooning representation used in HD7 satisfies its spatial scale separations. Under a flat safety factor profile, the KBM is more unstable and its mode structure tends to move with the peak of the ion temperature drive. Further simulation of the KBM in an HL-2A-like equilibrium shows that the β excitation threshold of the mode is lower than 0.2% and the dominant toroidal mode number is n = 4, which is consistent with the measured experimental spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) excited in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a new 3-D relativistic particle-in-cell code, we have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. The simulation has been performed using a long simulation system in order to study the nonlinear stages of the Weibel instability, the particle acceleration mechanism, and the shock structure. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic (HD) like shock structure. In the leading shock, electron density increases by a factor of <_ 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. We discuss the possible implication of our simulation results within the AGN and GRB context. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique. The same technique will be used to calculate radiation from accelerated electrons (positrons) in turbulent magnetic fields generated by Weibel instability.

Nishikawa, K.↗
DOE OSTI2013

Tunability of the terahertz space-charge modulation in a vacuum microdiode

Under certain conditions, space-charge limited emission in vacuum microdiodes manifests as clearly defined bunches of charge with a regular size and interval. The frequency corresponding to this interval is in the terahertz range. In this computational study, it is demonstrated that, for a range of parameters, conducive to generating THz frequency oscillations, the frequency is dependant only on the cold cathode electric field and on the emitter area. For a planar micro-diode of given dimension, the modulation frequency can be easily tuned simply by varying the applied potential. Simulations of the microdiode are done for 84 different combinations of emitter area, applied voltage, and gap spacing, using a molecular dynamics based code with exact Coulomb interaction between all electrons in the vacuum gap, which is of the order 100. It is found, for a fixed emitter area, that the frequency of the pulse train is solely dependant on the vacuum electric field in the diode, described by a simple power law. It is also found that, for a fixed value of the electric field, the frequency increases with diminishing size of the emitting spot on the cathode. Some observations are made on the spectral quality, and how it is affected by the gap spacing in the diode and the initial velocity of the electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Relativistic electron acceleration by oblique whistler waves

Test-particle simulations of electrons interacting with finite-amplitude, obliquely propagating whistler waves are carried out in order to investigate the acceleration of relativistic electrons by these waves. According to the present findings, an efficient acceleration of relativistic electrons requires a narrow range of oblique propagation angles, close to the whistler resonance cone angle, when the wave amplitude is held constant at relatively low value. For a constant wave propagation angle, it is found that a range of oblique whistler wave amplitudes permits the acceleration of relativistic electrons to O(MeV) energies. An initial distribution of test electrons is shown to form a power-law distribution when plotted in energy space. It is also found that the acceleration is largely uniform in electron pitch-angle space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Unconventional ballooning structures for toroidal drift waves

With strong gradients in the pedestal of high confinement mode (H-mode) fusion plasmas, gyrokinetic simulations are carried out for the trapped electron and ion temperature gradient modes. A broad class of unconventional mode structures is found to localize at arbitrary poloidal positions or with multiple peaks. It is found that these unconventional ballooning structures are associated with different eigen states for the most unstable mode. At weak gradient (low confinement mode or L-mode), the most unstable mode is usually in the ground eigen state, which corresponds to a conventional ballooning mode structure peaking in the outboard mid-plane of tokamaks. However, at strong gradient (H-mode), the most unstable mode is usually not the ground eigen state and the ballooning mode structure becomes unconventional. This result implies that the pedestal of H-mode could have better confinement than L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fully kinetic simulations of strong steady-state collisional planar plasma shocks

We report on simulations of strong, steady-state collisional planar plasma shocks with fully kinetic ions and electrons, independently confirmed by two fully kinetic codes (an Eulerian continuum and a Lagrangian particle-in-cell). While kinetic electrons do not fundamentally change the shock structure as compared with fluid electrons, we find an appreciable rearrangement of the preheat layer, associated with nonlocal electron heat transport effects. The electron heat-flux profile qualitatively agrees between kinetic- and fluid-electron models, suggesting a certain level of “stiffness,” though substantial nonlocality is observed in the kinetic heat flux. We also find good agreement with nonlocal electron heat-flux closures proposed in the literature. Finally, in contrast to the classical hydrodynamic picture, we find a significant collapse in the “precursor” electric-field shock at the preheat layer leading edge, which correlates with the electron-temperature gradient relaxation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analysis of wave mode content in fully turbulent, moderately collisional plasma laboratory experiment and kinetic simulation

Our major goal for the project was to make magnetic field and plasma density fluctuation measurements using identical diagnostic and acquisition setups on two different laboratory-based turbulent plasma devices: SSX at Swarthmore College and a new plasma source in development at Bryn Mawr College (BMX). We used single-point correlation and k-filtering metrics using these experimental measurements. These metrics were compared to kinetic simulations of the experimental environments in order to determine wave mode content of the laboratory turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Numerical Investigation of Plasma Detachment in Magnetic Nozzle Experiments

At present there exists no generally accepted theoretical model that provides a consistent physical explanation of plasma detachment from an externally-imposed magnetic nozzle. To make progress towards that end, simulation of plasma flow in the magnetic nozzle of an arcjet experiment is performed using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The aim is to compare the computational results with various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver (see Fig. I), and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kuriki and Okada. This geometry is modeled because there is a substantial amount of experimental data that can be compared to the computational results, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how they scale with different input parameters. Further studies will focus on modeling literature experiments both for the purpose of additional code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗
DOE OSTI2018

Linear analyses of peeling-ballooning modes in high beta pedestal plasmas

We present the linear simulations of edge plasma instabilities using the 3-field peeling-ballooning model and gyro-Landau-fluid model under the BOUT++ framework. A series of realistic equilibria of shifted circular geometry are generated by a global equilibrium solver CORSICA, where the Shafranov shift, elongation effects, and bootstrap current are included. The linear growth rate spectrum of the peeling-ballooning modes is shown in a wide range of pressure gradient and parallel current density in the pedestal region. The results show that the bootstrap current stabilizes high beta ballooning modes. Here, the simulations with different fractions of bootstrap current indicate a trend for the existence of the high beta peeling-ballooning mode stability region. Taking the kinetic effects into account, the linear simulations of kinetic peeling-ballooning mode using the gyro-Landau-fluid model show that this region can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗