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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

The influence of ions and the induced secondary emission on the nanosecond high-gradient microwave breakdown at metal surface

The mechanism of ultrafast breakdown at metal/vacuum interface in the high-power microwave waveguides is studied. In order to realize the nanosecond discharge, the required ambient gas pressure above the metal surface is approximately calculated as high as several Torr due to the low ionization-rate for high-energy electrons and short pulse. The local high pressure may come from the evaporated microscopic protrusions due to Joule heating and gas desorption. Besides, ions accelerated by the ambient space charge field could obtain sufficient high energy to collide and sputter the metal atoms to increase the ambient pressure. The positive feedbacks during the rapid discharge are studied by particle-in-cell simulation. The relatively high-energy ions could generate secondary electrons. It is shown that, as the positive feedback, the secondary electrons induce the gas desorption and stronger ionization, resulting in ion and electron density increasing as well as sheath field further increasing. As a result, more higher-energy ions bombard metal surface, leading to higher secondary electron yield and higher density plasma generated to cut off the microwave transmission finally. These nonlinear courses realize the ultrafast discharge in waveguides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Nonlinear shearing modes approach to the diocotron instability of a planar electron strip

The nonlinear evolution of the diocotron instability of a planar electron strip is investigated analytically by means of the nonlinear shearing mode for the solution of the initial and boundary value problems. The method is based on the sheared spatial coordinates which account for the motion of electron flow in the electrostatic field of the unstable diocotron modes in addition to the unperturbed sheared motion of the electron flow on the transformed shear coordinates. The time evolutions are studied by the solution of the initial and boundary value problems. The obtained solutions for the perturbed electrostatic potential include two nonlinear effects—the effect of the distortion of the boundaries of the planar electron strip and the effect of the coupling of the sheared nonmodal diocotron modes. It was proved by a two-dimensional particle-in-cell simulation that the developed theory is valid as long as the distortion of the boundaries of the basic shear flow does not change the frequency and growth rate of the linear diocotron instability in the transformed coordinates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Conservative discontinuous Galerkin interpolation: Sheared boundary conditions

Local studies of accretion disks and laboratory magnetized plasmas employ analytical coordinate mappings that introduce sheared boundary conditions (BCs). Here, we present a discontinuous Galerkin (DG) algorithm to apply such BCs based on projections and quadrature-free integration. The procedure is high-order accurate, preserves moments exactly and works in multiple dimensions. Tests of increasing complexity are provided, beginning with translations of one and two dimensional fields, followed by 3D and 5D simulations with sheared (twist-shift) BCs. Results show that the algorithm is (p + 1)-order accurate in the DG representation and (p + 2) -order accurate in the cell averages, with p being the order of the polynomial basis. Quantification of the algorithm's hyperdiffusion and discussion of aliasing errors are given. This technique enables conservative local simulations of plasma turbulence with DG, not possible until now.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Transport and losses of fusion-born alpha particles in the presence of tearing modes using the new Toroidal Accelerated Particle Simulator (TAPaS)

Here, the transport and losses of fusion-born alpha particles are studied in the presence of a single-helicity tearing mode, characterized by (m = 2, n = 1). The analysis is performed by means of the recently developed Toroidal Accelerated Particle Simulator (TAPaS). Although such modes have usually been believed to result only in a local flattening of the radial profiles, it is shown that the density profile can exhibit a global modification leading to significant losses of alpha particles. This is due to the fact that, although the magnetic field does not exhibit any chaotic behaviour, the trajectories of alpha particles do, as revealed by their Poincaré maps. Such results are in qualitative agreement with past observations and simulations of energetic particles generated by neutral beam injection in TFTR, DIII-D and AUG tokamaks. In-depth analysis is carried out to characterize the impact of the tearing mode on the transport and losses of fusion-born alpha-particles with a realistic density profile. The impact of the amplitude is evidenced. Moreover, the effect of the island rotation frequency is assessed based on a detailed analysis of the linear resonances in phase-space, in agreement with the simulation results. Finally, the probability density function of the exit time has been computed and the transport of alpha particles has been found to be anomalous.

anomalous transport↗
NASA NTRS1980

Some features of inverted-V events as seen from simulated double layers

Results from a numerical simulation of a double layer show some features similar to those of inverted-V events. The strong heating of thermal and precipitating electrons is observed along with extremely low frequency fluctuations found during inverted-V events. It is suggested that after the acceleration of auroral electrons by the double layer, the precipitating free electrons are heated by the nonlinear effects of the electron beam plasma instability. Fluctuations and pulsations of auroral electron fluxes during auroral events are caused by a relaxation type of oscillation. The finite extent of one dimensional plasma is simulated by solving the Vlasov and Poisson equations as an initial and boundary value problem.

Singh, N.↗
NASA NTRS2000

Magnetic Flux Compression Concept for Nuclear Pulse Propulsion and Power

The desire for fast, efficient interplanetary transport requires propulsion systems having short acceleration times and very high specific impulse attributes. Unfortunately, most highly efficient propulsion systems which are within the capabilities of present day technologies are either very heavy or yield very low impulse such that the acceleration time to final velocity is too long to be of lasting interest, One exception, the nuclear thermal thruster, could achieve the desired acceleration but it would require inordinately large mass ratios to reach the range of desired final velocities. An alternative approach, among several competing concepts that are beyond our modern technical capabilities, is a pulsed thermonuclear device utilizing microfusion detonations. In this paper, we examine the feasibility of an innovative magnetic flux compression concept for utilizing microfusion detonations, assuming that such low yield nuclear bursts can be realized in practice. In this concept, a magnetic field is compressed between an expanding detonation driven diamagnetic plasma and a stationary structure formed from a high temperature superconductor (HTSC). In general, we are interested in accomplishing two important functions: (1) collimation of a hot diamagnetic plasma for direct thrust production; and (2) pulse power generation for dense plasma ignition. For the purposes of this research, it is assumed that rnicrofusion detonation technology may become available within a few decades, and that this approach could capitalize on recent advances in inertial confinement fusion ICF) technologies including magnetized target concepts and antimatter initiated nuclear detonations. The charged particle expansion velocity in these detonations can be on the order of 10 (exp 6)- 10 (exp 7) meters per second, and, if effectively collimated by a magnetic nozzle, can yield the Isp and the acceleration levels needed for practical interplanetary spaceflight. The ability to ignite pure fusion micro-bursts with reasonable levels of input energy is an equally challenging scientific problem. It remains to be seen, however, whether an effective ignition driver can be developed which meets the requirements for practical spaceflight application (namely high power density, compactness, low weight, and low cost). In this paper, system level performance and design issues are examined including generator performance, magnetic flux compression processes, magnetic diffusion processes, high temperature superconductor (HTSC) material properties, plasmadynamic processes, detonation plasma expansion processes, magnetohydrodynamic instabilities, magnetic nozzle performance, and thrust production performance. Representative generator performance calculations based on a simplified skin layer formulation are presented as well as the results of exploratory small-scale laboratory experiments on magnetic flux diffusion in HTSC materials. In addition, planned follow-on scientific feasibility experiments are described which utilize high explosive detonations and high energy gas discharges to simulate the plasma conditions associated with thermonuclear micro-detonations.

Litchford, Ronald J.↗
DOE OSTI2018

Alternative model of space-charge-limited thermionic current flow through a plasma

It is widely assumed that thermionic current flow through a plasma is limited by a “space-charge-limited” (SCL) cathode sheath that consumes the hot cathode's negative bias and accelerates upstream ions into the cathode. In this paper, we formulate a fundamentally different current-limited mode. In the “inverse” mode, the potentials of both electrodes are above the plasma potential, so that the plasma ions are confined. The bias is consumed by the anode sheath. There is no potential gradient in the neutral plasma region from resistivity or presheath. The inverse cathode sheath pulls some thermoelectrons back to the cathode, thereby limiting the circuit current. Thermoelectrons entering the zero-field plasma region that undergo collisions may also be sent back to the cathode, further attenuating the circuit current. In planar geometry, the plasma density is shown to vary linearly across the electrode gap. A continuum kinetic planar plasma diode simulation model is set up to compare the properties of current modes with classical, conventional SCL, and inverse cathode sheaths. SCL modes can exist only if charge-exchange collisions are turned off in the potential well of the virtual cathode to prevent ion trapping. With the collisions, the current-limited equilibrium must be inverse. Inverse operating modes should therefore be present or possible in many plasma devices that rely on hot cathodes. Evidence from past experiments is discussed. Finally, the inverse mode may offer opportunities to minimize sputtering and power consumption that were not previously explored due to the common assumption of SCL sheaths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Gyrokinetic particle simulation of a field reversed configuration

Gyrokinetic particle simulation of the field-reversed configuration (FRC) has been developed using the gyrokinetic toroidal code (GTC). The magnetohydrodynamic equilibrium is mapped from cylindrical coordinates to Boozer coordinates for the FRC core and scrape-off layer (SOL), respectively. A field-aligned mesh is constructed for solving self-consistent electric fields using a semi-spectral solver in a partial torus FRC geometry. This new simulation capability has been successfully verified and driftwave instability in the FRC has been studied using the gyrokinetic simulation for the first time. Initial GTC simulations find that in the FRC core, the ion-scale driftwave is stabilized by the large ion gyroradius. In the SOL, the driftwave is unstable on both ion and electron scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron Modulational Instability in the Strong Turbulent Regime for an Electron Beam Propagating in a Background Plasma

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. In conclusion, different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A phase-shift-periodic parallel boundary condition for low-magnetic-shear scenarios

Abstract We formulate a generalized periodic boundary condition as a limit of the standard twist-and-shift parallel boundary condition that is suitable for simulations of plasmas with low magnetic shear. This is done by applying a phase shift in the binormal direction when crossing the parallel boundary. While this phase shift can be set to zero without loss of generality in the local flux-tube limit when employing the twist-and-shift boundary condition, we show that this is not the most general case when employing periodic parallel boundaries, and may not even be the most desirable. A non-zero phase shift can be used to avoid the convective cells that plague simulations of the three-dimensional Hasegawa–Wakatani system, and is shown to have measurable effects in periodic low-magnetic-shear gyrokinetic simulations. We propose a numerical program where a sampling of periodic simulations at random pseudo-irrational flux surfaces are used to determine physical observables in a statistical sense. This approach can serve as an alternative to applying the twist-and-shift boundary condition to low-magnetic-shear scenarios, which, while more straightforward, can be computationally demanding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Fluid vs. kinetic magnetic reconnection with strong guide fields

The fast rates of magnetic reconnection found in both nature and experiments are important to understand theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a more detailed derivation of the analytic model presented in a recent letter and present additional simulation results to support the conclusions that the magnetic reconnection rate in this regime is independent of both collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an alternative interpretation for the numerical results presented therein.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Stabilization of high-compression, indirect-drive inertial confinement fusion implosions using a 4-shock adiabat-shaped drive

Hydrodynamic instabilities and poor fuel compression are major factors for capsule performance degradation in ignition experiments on the National Ignition Facility. Using a recently developed laser drive profile with a decaying first shock to tune the ablative Richtmyer-Meshkov (ARM) instability and subsequent in-flight Rayleigh-Taylor growth, we have demonstrated reduced growth compared to the standard ignition pulse whilst maintaining conditions for a low fuel adiabat needed for increased compression. Using in-flight x-ray radiography of pre-machined modulations, the first growth measurements using this new ARM-tuned drive have demonstrated instability growth reduction of ∼4× compared to the original design at a convergence ratio of ∼2. Corresponding simulations give a fuel adiabat of ∼1.6, similar to the original goal and consistent with ignition requirements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical study on Rayleigh-Taylor instabilities in the lightning return stroke

The Rayleigh-Taylor (R-T) instabilities are important hydrodynamics and magnetohydrodynamics (MHD) phenomena that are found in systems in high energy density physics and normal fluids. The formation and evolution of the R-T instability at channel boundary during back-flow of the lightning return stroke are analyzed using the linear perturbation theory and normal mode analysis methods, and the linear growth rate of the R-T instability in typical condition for lightning return stroke channel is obtained. Then, the R-T instability phenomena of lightning return stroke are simulated using a two-dimensional Eulerian finite volumes resistive radiation MHD code. The numerical results show that the evolution characteristics of the R-T instability in the early stage of back-flow are consistent with theoretical predictions obtained by linear analysis. The simulation also yields more evolution characteristics for the R-T instability beyond the linear theory. The results of this work apply to some observed features of the return stroke channel and further advance previous theoretical and experimental work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

An X-band high-impedance relativistic klystron amplifier with an annular explosive cathode

The feasibility of employing an annular beam instead of a solid one in the X-band high-impedance relativistic klystron amplifier (RKA) is investigated in theory and simulation. Small-signal theory analysis indicates that the optimum bunching distance, fundamental current modulation depth, beam-coupling coefficient, and beam-loaded quality factor of annular beams are all larger than the corresponding parameters of solid beams at the same beam voltage and current. An annular beam RKA and a solid beam RKA with almost the same geometric parameters are compared in particle-in-cell simulation. Output microwave power of 100 MW, gain of 50 dB, and power conversion efficiency of 42% are obtained in an annular beam RKA. The annular beam needs a 15% lower uniform guiding magnetic field than the solid beam. Our investigations demonstrate that we are able to use a simple annular explosive cathode immersed in a lower uniform magnetic field instead of a solid thermionic cathode in a complicated partially shielding magnetic field for designing high-impedance RKA, which avoids high temperature requirement, complicated electron-optical system, large area convergence, high current density, and emission uniformity for the solid beam. An equivalent method for the annular beam and the solid beam on bunching features is proposed and agrees with the simulation. The annular beam has the primary advantages over the solid beam that it can employ the immersing uniform magnetic field avoiding the complicated shielding magnetic field system and needs a lower optimum guiding field due to the smaller space charge effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Estimation of China Fusion Engineering Test Reactor performance and burning fraction in different pellet fueling scenarios by a multi-species radial transport model

Tritium self-sufficiency in future deuterium–tritium fusion reactors is a crucial challenge. As an engineering test reactor, the China Fusion Engineering Test Reactor requires a burning fraction of 3% for the goal to test the accessibility to the future fusion plant. To self-consistently simulate burning plasmas with profile changes in pellet injection scenarios and to estimate the corresponding burning fraction, a one-dimensional multi-species radial transport model is developed in the BOUT++ framework. Here several pellet-fueling scenarios are then tested in the model. The results show that the increased fueling depth improves the burning fraction by particle confinement improvement and fusion power increase. Nevertheless, by increasing the depth, the pellet cooling-down may significantly lower the temperature in the core region. Taking the density perturbation into consideration, the reasonable parameters of the fueling scenario in these simulations are estimated as pellet radius ${r}_{{\rm{p}}}=3\,{\rm{mm}},$ injection rate $=\,4\,\mathrm{Hz},$ and pellet injection velocity $=\,1000\mbox{-}2000\,{\rm{m}}\,{{\rm{s}}}^{-1}$ without drift or $450\,{\rm{m}}\,{{\rm{s}}}^{-1}$ with high-field-side drift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Self-consistent evolution of plasma discharge and electromagnetic fields in a microwave pulse compressor

Nanosecond-scale evolution of plasma and RF electromagnetic fields during the release of energy from a microwave pulse compressor with a plasma interference switch was investigated numerically using the code MAGIC. The plasma was simulated in the scope of the gas conductivity model in MAGIC. The compressor embodied an S-band cavity and H-plane waveguide tee with a shorted side arm filled with pressurized gas. In a simplified approach, the gas discharge was initiated by setting an external ionization rate in a layer crossing the side arm waveguide in the location of the electric field antinode. It was found that with increasing ionization rate, the microwave energy absorbed by the plasma in the first few nanoseconds increases, but the absorption for the whole duration of energy release, on the contrary, decreases. In a hybrid approach modeling laser ignition of the discharge, seed electrons were set around the electric field antinode. In this case, the plasma extends along the field forming a filament and the plasma density increases up to the level at which the electric field within the plasma decreases due to the skin effect. Then, the avalanche rate decreases but the density still rises until the microwave energy release begins and the electric field becomes insufficient to support the avalanche process. The extraction of the microwave pulse limits its own power by terminating the rise of the plasma density and filament length. For efficient extraction, a sufficiently long filament of dense plasma must have sufficient time to be formed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗