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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Experimental quantification of the impact of heterogeneous mix on thermonuclear burn

In inertial confinement fusion, deuterium–tritium (DT) fuel is brought to densities and temperatures where fusion ignition occurs. However, mixing of the ablator material into the fuel may prevent ignition by diluting and cooling the fuel. MARBLE experiments at the National Ignition Facility provide new insight into how mixing affects thermonuclear burn. These experiments use laser-driven capsules containing deuterated plastic foam and tritium gas. Embedded within the foam are voids of known sizes and locations, which control the degree of heterogeneity of the fuel. Initially, the reactants are separated, with tritium concentrated in the voids and deuterium in the foam. During the implosion, mixing occurs between the foam and gas materials, leading to DT fusion reactions in the mixed region. Here, it is shown that by measuring the ratios of DT and deuterium–deuterium neutron yields for different macropore sizes and gas compositions, the effects of mix heterogeneity on thermonuclear burn may be quantified, supporting an improved understanding of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Integrated diagnostic analysis of inertial confinement fusion capsule performance

A conceptual model is developed for typical inertial confinement fusion implosion conditions that integrates available diagnostic information to determine the stagnation properties of the interior fill and surrounding shell. Assuming pressure equilibrium at peak compression and invoking radiative and equation-of-state relations, the pressure, density, and electron temperature are obtained by optimized fitting of the experimental output to smooth, global functional forms. Typical observational data that may be used includes x-ray self-emission, directional neutron time-of-flight signals, neutron yield, high-resolution x-ray spectra, and radiographic images. This approach has been validated by comparison with radiation-hydrodynamic simulations, producing semi-quantitative agreement. Model results implicate poor kinetic energy coupling to the hot core as the primary cause of the observed low thermonuclear burn yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Ping-pong modes and higher-periodicity multipactor

Multipactor is a vacuum discharge driven by secondary electron emission. Multiple period multipactors have long been known to exist but have been studied less extensively. In a period-n multipactor, electrons undergo multiple impacts in one rf period, with the synchronous phase alternating periodically between multiple values. A novel resonant form is proposed that combines one- and two-surface impacts within a single period, provided the total transit time is an odd number of rf half-periods and the product of secondary yields exceeds unity. For low fD products, the simplest such mode is shown to significantly increase the upper electric field boundary of the multipacting region and lead to overlap of higher-order bands. The results agree nicely with 3-D particle-in-cell code simulations. An alternative, map-based method is introduced for analyzing higher-periodicity multipactor. Practical implications of the findings are discussed, including consequences for multipactor suppression strategies using a dc magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Alpha channeling with high-field launch of lower hybrid waves

Although lower hybrid waves are effective at driving currents in present-day tokamaks, they are expected to interact strongly with high-energy particles in extrapolating to reactors. In the presence of a radial alpha particle birth gradient, this interaction can take the form of wave amplification rather than damping. While it is known that this amplification more easily occurs when launching from the tokamak high-field side, the extent of this amplification has not been made quantitative. Here, by tracing rays launched from the high-field-side of a tokamak, the required radial gradients to achieve amplification are calculated for a temperature and density regime consistent with a hot-ion-mode fusion reactor. These simulations, while valid only in the linear regime of wave amplification, nonetheless illustrate the possibilities for wave amplification using high-field launch of the lower hybrid wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Kinetic simulations of laser parametric amplification in magnetized plasmas

Laser pulse compression using magnetized resonance near the upper-hybrid frequency is promising for achieving higher output intensity in regimes previously thought impossible using unmagnetized plasmas. Using one dimensional particle-in-cell simulations, we verify that, by partially replacing plasma with an external transverse magnetic field of megagauss scale, the output pulse can be intensified by a factor of a few, due to the increased allowable amplification time despite a decreased growth rate. Further improvement is impeded by the generation of an electromagnetic wakefield, to which the amplified pulse loses more energy than it does in the unmagnetized case. This limitation can however be circumvented by the use of a stronger pump. In contrast to unmagnetized compression, the magnetized amplification remains efficient when the pump intensity is well above the wavebreaking threshold, until a higher phase-mixing threshold is exceeded. Here, this surprising resilience to wavebreaking in magnetized plasma is of great benefit for magnetized compression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Evidence of electron microbunching in laser-driven modulated downramp injection and prospects for beam-driven implementation

Plasma accelerators can generate high-energy, high-brightness electron beams over centimeter-scale distances, offering novel pathways to compact x-ray free-electron lasers. Generating beams pre-bunched at the desired radiation wavelength would significantly enhance longitudinal coherence and reduce saturation length. Plasma density-modulated downramp injection offers an in-situ way to generate such beams with nanometer-scale bunching. Here we report the first experimental evidence of this mechanism in a laser-driven wakefield accelerator, showing that modulated density downramps generate modulated electron energy spectra absent in unmodulated cases. Particle-in-cell simulations reproduce these observations and reveal bunching factors of 0.05 at 0.4 μm, with a compression factor of approximately 7. Building on this demonstration, we propose a beam-driven implementation for FACET-II to generate multi-GeV beams pre-bunched at hundreds of nanometers wavelength, with sub-micrometer emittance, kiloampere peak current, and sub-percent slice energy spread. Two-stage magnetic compression enables tunable bunching from optical to extreme ultraviolet wavelengths while achieving peak currents exceeding 100 kA. Coherent transition radiation calculations confirm diagnostic feasibility. This approach offers a promising path towards compact, high-energy pre-bunched electron sources for advanced photon science applications.

electron microbunching↗
DOE OSTI2015

Numerical simulations of stable explosive-emission center plasma expansion in vacuum

We present the results of our numerical simulations of the stable cathode-flare expansion stage in the vacuum-breakdown switching period. We used the model of two-fluid hydrodynamics with corrections that allowed it to be applied to a consistent description of the electron and ion transport in plasma and the electron transport in a plasma-anode vacuum gap without imposing any joining conditions at the expanding plasma boundary. We have analyzed the dynamics of the density, directed velocity, and temperature profiles for both components and the influence of various factors on the plasma-cloud expansion velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Complexity and Intermittent Turbulence in Space Plasmas

Sporadic and localized interactions of coherent structures arising from plasma resonances can be the origin of "complexity" of the coexistence of non- propagating spatiotemporal fluctuations and propagating modes in space plasmas. Numerical simulation results are presented to demonstrate the intermittent character of the non-propagating fluctuations. The technique of the dynamic renormalization-group is introduced and applied to the study of scale invariance of such type of multiscale fluctuations. We also demonstrate that the particle interactions with the intermittent turbulence can lead to the efficient energization of the plasma populations. An example related to the ion acceleration processes in the auroral zone is provided.

Chang, Tom↗
DOE OSTI2017

Monte Carlo simulation of excitation and ionization collisions with complexity reduction

Kinetic simulation of plasmas with detailed excitation and ionization collisions presents a significant computational challenge due to the multiscale feature of the collisional rates. In the present work, we propose a complexity reduction method based on atomic level grouping for modeling excitation and ionization collisions. High order of accuracy of the reduction method is realized by allowing an internal distribution within each group. We apply the reduction method to the standard Monte Carlo collision algorithm to model an atomic Hydrogen plasma. Numerical results suggest that the stiffness of the collisional kinetics can be significantly reduced with minimal loss in accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Comparison of toroidicity-induced Alfvén eigenmodes and energetic particle modes by gyrokinetic particle simulations

This work reports on linear global gyrokinetic particle simulations of the excitation of toroidicity-induced Alfvén eigenmodes (TAE) and energetic particle modes (EPM), and the comparison between these two modes. The TAE excitation by antenna clarifies the magnetohydrodynamic (MHD) mode structure and the discrete eigenmode exists in the gap between the upper and lower accumulation points. The TAE excitation by fast ions modifies the MHD mode structure because of radial symmetry breaking and the eigenmode frequency moves towards the lower accumulation point. The phase space structure of fast ions shows that both passing and trapped particles contribute to the TAE excitation and that trapped particles dominate the wave-particle resonance in our simulations. The growth rate of TAE is sensitive to the fast ion energy, density, and density gradient, which are also important factors contributing to the transition of the TAE to the EPM. The gyrokinetic particle simulations also confirm the excitation of EPM when the drive is stronger. The frequency of the EPM is determined by the characteristic frequencies of fast ion motion in toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗
DOE OSTI2013

Kinetic aspects of the ion current layer in a reconnection outflow exhaust

Kinetic aspects of the ion current layer at the center of a reconnection outflow exhaust near the X-type region are investigated by a two-dimensional particle-in-cell (PIC) simulation. The layer consists of magnetized electrons and unmagnetized ions that carry a perpendicular electric current. The ion fluid appears to be nonideal, sub-Alfvénic, and nondissipative. The ion velocity distribution functions contain multiple populations, such as global Speiser ions, local Speiser ions, and trapped ions. The particle motion of the local Speiser ions in an appropriately rotated coordinate system explains the ion fluid properties very well. The trapped ions are the first demonstration of the regular orbits in the chaotic particle dynamics [Chen and Palmadesso, J. Geophys. Res. 91, 1499 (1986)] in self-consistent PIC simulations. They would be observational signatures in the ion current layer near reconnection sites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Simulation study of CO 2 laser-plasma interactions and self-modulated wakefield acceleration

3D numerical simulations of the interaction of a powerful CO 2 laser with hydrogen jets demonstrating the role of ionization on the characteristics of induced wakes are presented. Simulations using SPACE, a parallel relativistic Particle-in-Cell code, are performed in support of the plasma wakefield accelerator experiments being conducted at the Brookhaven National Laboratory (BNL) Accelerator Test Facility (ATF). A novelty of the SPACE code is its set of efficient atomic physics algorithms that compute ionization and recombination rates on the grid and transfer them to particles. The influence of ionization on the spectrum of pump laser has been studied for a range of gas densities.Simulations reproduce both Stokes and anti-Stokes shifts in the spectrum of pump laser, similar to those observed in experiments in the spectrum of probe. A good agreement has been achieved with the experiments on the effect of variation in gas density on Stokes/anti-Stokes intensity. In addition, self-injection and trapping of electrons into the self-modulated wakes has been observed and analyzed. Here, the experimentally validated code SPACE will be used for predictive simulations to guide future experiments at BNL ATF.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multi-camera imaging to characterize jet and liner uniformity on the Plasma Liner Experiment (PLX)

At Plasma Liner Experiment, a set of 36 coaxial plasma guns are deployed quasi-uniformly over a 9 ft diameter spherical chamber and are used to form a high-Z spherically compressive plasma liner. Simulations indicate that for the concept to ultimately achieve optimal target density and temperature, a high degree of timing uniformity is required between all guns. Therefore, to aid in quantifying and correcting gun-to-gun nonuniformities, a key diagnostic will consist of up to six fisheye-view CCD cameras positioned inside the main chamber such that each has all plasma guns within its view. The individual cameras can be triggered at different times to determine each plasma jet’s muzzle velocity and structure for different operating conditions. This camera array is currently under development, and the implementation needs and challenges for this camera array are discussed here. Additionally, we detail the analysis methodology for determining jet-to-jet uniformity deviations and how we can correct them, thereby improving overall liner uniformity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Propagation speed of rotation signals for field lines undergoing magnetic reconnection

Reconnection is associated with two bending of the magnetic field lines. Considering the usual plane of a 2D reconnection simulation, the first bending is in-plane and produces the needed topological changes by bringing oppositely directed filed lines in proximity. The second is typical of fast reconnection and is out of plane, leading to the formation of the Hall magnetic field. This second rotation has recently been observed to proceed at superAlfvénic speeds and to carry substantial energy fluxes (Shay et al., Phys. Rev. Lett. 107, 065001 (2011)). We revisit these rotations with a new diagnostics based on dispersing a multitude of virtual probes into a kinetic simulation, akin the approach of multi spacecraft missions. The results of the new diagnostics are compared with the theory of characteristics applied to the two fluid model. The comparison of virtual probes and the method of characteristics confirm the findings relative to the out of plane rotation and uncover the existence of two families of characteristics. Both are observed in the simulation. The early stage of reconnection develops on the slower compressional branch and the later faster phase develops on the faster torsional branch. The superAlfvénic signal is only relevant in the second phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2013

Electron kinetics in a microdischarge in nitrogen at atmospheric pressure

Electron kinetics during a microdischarge in nitrogen at atmospheric pressure is studied using the one-dimensional Particle-in-Cell/Monte Carlo Collisions model. It is obtained that the electron energy distribution function can be divided into three parts, namely, the non-equilibrium low-energy part, the Maxwellian function at moderate energies, and the high-energy tail. Simulation results showed that the role of the high-energy tail of electron energy distribution increases, when the distance between electrodes increases.

36 MATERIALS SCIENCE↗
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↗