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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Pink-beam focusing with a one-dimensional compound refractive lens

The performance of a cooled Be compound refractive lens (CRL) has been tested at the Advanced Photon Source (APS) to enable vertical focusing of the pink beam and permit the X-ray beam to spatially overlap with an 80 µm-high low-density plasma that simulates astrophysical environments. Focusing the fundamental harmonics of an insertion device white beam increases the APS power density; here, a power density as high as 500 W mm –2 was calculated. A CRL is chromatic so it does not efficiently focus X-rays whose energies are above the fundamental. Only the fundamental of the undulator focuses at the experiment. A two-chopper system reduces the power density on the imaging system and lens by four orders of magnitude, enabling imaging of the focal plane without any X-ray filter. As a result, a method to measure such high power density as well as the performance of the lens in focusing the pink beam is reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Fast-ion transport in quasisymmetric equilibria in the presence of a resonant Alfvénic perturbation

Significant progress has been made in designing magnetic fields that provide excellent confinement of the guiding-centre trajectories of alpha particles using quasisymmetry (QS). Given the reduction in this transport channel, we assess the impact of resonant Alfvén eigenmodes (AEs) on the guiding-centre motion. The AE amplitudes are chosen to be consistent with experimental measurements and large-scale simulations. We evaluate the drift resonance condition, phase-space island width and island overlap criterion for quasisymmetric configurations. Kinetic Poincaré plots elucidate features of the transport, including stiff transport above a critical perturbation amplitude. Our analysis highlights key departures from the AE-driven transport in tokamaks, such as the avoidance of phase-space island overlap in quasihelical configurations and the enhanced transport due to wide phase-space islands in low magnetic shear configurations. In configurations that are closer to QS, with QS deviations $\delta B/B_0 \lesssim 10^{-3}$ , the transport is primarily driven by the AE, while configurations that are further from QS, $\delta B/B_0 \sim 10^{-2}$ , experience significant transport due to the QS-breaking fields in addition to the AE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A method for evaluating the mean preheat temperature in X-ray driven ablation

A novel method is proposed for evaluating the mean preheat temperature in X-ray driven ablation, based on the equation of state (EOS) of the ablator and the radiation hydrodynamic simulation. The equation of state of plastic (CH) has been discussed in detail. There are two types of planar CH in simulations, with the thick target being 10 μm thicker than the thin target. The difference between the transmission fluxes of the two types of targets can represent the energy absorbed by the last 10 μm of the thick target (or the preheated layer). This energy approximates the internal energy of the preheated layer. The mean preheat temperature of the preheated layer has also been obtained from simulations. The simulation results show that the relationship between the absorbed energy and the mean preheat temperature is similar to the EOS of CH for different conditions (e.g., different values of M-band fraction and radiation temperature) and can be written as ε=2.530×10{sup 11}T{sup ¯1.444} when the mean preheat temperature is below 12 eV. For these cases, the relationship between the surface preheat temperature T{sub S} and the mean preheat temperature T{sup ¯} was T{sub S}=0.63T{sup ¯}. This relation provides the means for demonstrating the proposed method, because the transmission fluxes and the surface preheat temperature T{sub S} can be measured experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Modeling of long range frequency sweeping for energetic particle modes

Long range frequency sweeping events are simulated numerically within a one-dimensional, electrostatic bump-on-tail model with fast particle sources and collisions. The numerical solution accounts for fast particle trapping and detrapping in an evolving wave field with a fixed wavelength, and it includes three distinct collisions operators: Drag (dynamical friction on the background electrons), Krook-type collisions, and velocity space diffusion. The effects of particle trapping and diffusion on the evolution of holes and clumps are investigated, and the occurrence of non-monotonic (hooked) frequency sweeping and asymptotically steady holes is discussed. The presented solution constitutes a step towards predictive modeling of frequency sweeping events in more realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Particle-in-cell simulations of magnetic reconnection in laser-plasma experiments on Shenguang-II facility

Recently, magnetic reconnection has been realized in high-energy-density laser-produced plasmas. Plasma bubbles with self-generated magnetic fields are created by focusing laser beams to small-scale spots on a foil. The bubbles expand into each other, which may then drive magnetic reconnection. The reconnection experiment in laser-produced plasmas has also been conducted at Shenguang-II (SG-II) laser facility, and the existence of a plasmoid was identified in the experiment [Dong et al., Phys. Rev. Lett. 108, 215001 (2012)]. In this paper, by performing two-dimensional (2-D) particle-in-cell simulations, we investigate such a process of magnetic reconnection based on the experiment on SG-II facility, and a possible explanation for the formation of the plasmoid is proposed. The results show that before magnetic reconnection occurs, the bubbles squeeze strongly each other and a very thin current sheet is formed. The current sheet is unstable to the tearing mode instability, and we can then observe the formation of plasmoid(s) in such a multiple X-lines reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Magnetic reconnection process in transient coaxial helicity injection

The physics of magnetic reconnection and fast flux closure in transient coaxial helicity injection experiments in NSTX is examined using resistive MHD simulations. These simulations have been performed using the NIMROD code with fixed boundary flux (including NSTX poloidal coil currents) in the NSTX experimental geometry. Simulations show that an X point is formed in the injector region, followed by formation of closed flux surfaces within 0.5 ms after the driven injector voltage and injector current begin to rapidly decrease. As the injector voltage is turned off, the field lines tend to untwist in the toroidal direction and magnetic field compression exerts a radial J × B force and generates a bi-directional radial E{sub toroidal}×B{sub poloidal} pinch flow to bring oppositely directed field lines closer together to reconnect. At sufficiently low magnetic diffusivity (high Lundquist number), and with a sufficiently narrow injector flux footprint width, the oppositely directed field lines have sufficient time to reconnect (before dissipating), leading to the formation of closed flux surfaces. The reconnection process is shown to have transient Sweet-Parker characteristics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Increase in the energy density of the pinch plasma in 3D implosion of quasi-spherical wire arrays

Results are presented from experimental studies of the characteristics of the soft X-ray (SXR) source formed in the implosion of quasi-spherical arrays made of tungsten wires and metalized kapron fibers. The experiments were carried out at the Angara-5-1 facility at currents of up to 3 MA. Analysis of the spatial distribution of hard X-ray emission with photon energies above 20 keV in the pinch images taken during the implosion of quasi-spherical tungsten wire arrays (QTWAs) showed that a compact quasi-spherical plasma object symmetric with respect to the array axis formed in the central region of the array. Using a diffraction grazing incidence spectrograph, spectra of SXR emission with wavelengths of 20–400 Å from the central, axial, and peripheral regions of the emission source were measured with spatial resolutions along the array radius and height in the implosion of QTWAs. It is shown that the emission spectra of the SXR sources formed under the implosion of quasi-spherical and cylindrical tungsten wire arrays at currents of up to 3 MA have a maximum in the wavelength range of 50–150 Å. It is found that, during the implosion of a QTWA with a profiled linear mass, a redistribution of energy in the emission spectrum takes place, which indicates that, during 3D implosion, the energy of longitudinal motion of the array material additionally contributes to the radiation energy. It is also found that, at close masses of the arrays and close values of the current in the range of 2.4{sup −3} MA, the average energy density in the emission source formed during the implosion of a quasi-spherical wire array is larger by a factor of 7 than in the source formed during the implosion of a cylindrical wire array. The experimental data were compared with results of 3D simulations of plasma dynamics and radiation generation during the implosion of quasi-spherical wire arrays with a profiled mass by using the MARPLE-3D radiative magnetohydrodynamic code, developed at the Keldysh Institute of Applied Mathematics, Russian Academy of Sciences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Kinetic simulations of gas breakdown in the dense plasma focus

We describe the first fully-kinetic, collisional, and electromagnetic simulations of the breakdown phase of a MA-scale dense plasma focus and are shown to agree with measured electrical characteristics, including breakdown time. In the model, avalanche ionization is driven by cathode electron emission and this results in incomplete gas breakdown along the insulator. This reinforces the importance of the conditioning process that creates a metallic layer on the insulator surface. The simulations, nonetheless, help explain the relationship between the gas pressure, the insulator length, and the coaxial gap width. In the past, researchers noted three breakdown patterns related to pressure. Simulation and analytic results show that at low pressures, long ionization path lengths lead to volumetric breakdown, while high pressures lead to breakdown across the relatively small coaxial electrode gap. In an intermediate pressure regime, ionization path lengths are comparable to the insulator length which promotes ideal breakdown along the insulator surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

X-pinch dynamics: Neck formation and implosion

We propose a model that describes the neck formation and implosion in an X-pinch. The process is simulated to go in two stages. The first stage is neck formation. This stage begins with an electrical explosion of the wires forming the X-pinch, and at the end of the stage, a micropinch (neck) is formed in the region where the wires are crossed. The second stage is neck implosion. The implosion is accompanied by outflow of matter from the neck region, resulting in the formation of a “hot spot”. Analytical estimates obtained in the study under consideration indicate that these stages are approximately equal in duration. Having analyzed the neck implosion dynamics, we have verified a scaling which makes it possible to explain the observed dependences of the time of occurrence of an x-ray pulse on the X-pinch current and mass.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Simulation of nonstationary phenomena in atmospheric-pressure glow discharge

Nonstationary processes in atmospheric-pressure glow discharge manifest themselves in spontaneous transitions from the normal glow discharge into a spark. In the experiments, both so-called completed transitions in which a highly conductive constricted channel arises and incomplete transitions accompanied by the formation of a diffuse channel are observed. A model of the positive column of a discharge in air is elaborated that allows one to interpret specific features of the discharge both in the stationary stage and during its transition into a spark and makes it possible to calculate the characteristic oscillatory current waveforms for completed transitions into a spark and aperiodic ones for incomplete transitions. The calculated parameters of the positive column in the glow discharge mode agree well with experiment. Data on the densities of the most abundant species generated in the discharge (such as atomic oxygen, metastable nitrogen molecules, ozone, nitrogen oxides, and negative oxygen ions) are presented.

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

Resistive magnetohydrodynamic reconnection: Resolving long-term, chaotic dynamics

In this paper, we address the long-term evolution of an idealised double current system entering reconnection regimes where chaotic behavior plays a prominent role. Our aim is to quantify the energetics in high magnetic Reynolds number evolutions, enriched by secondary tearing events, multiple magnetic island coalescence, and compressive versus resistive heating scenarios. Our study will pay particular attention to the required numerical resolutions achievable by modern (grid-adaptive) computations, and comment on the challenge associated with resolving chaotic island formation and interaction. We will use shock-capturing, conservative, grid-adaptive simulations for investigating trends dominated by both physical (resistivity) and numerical (resolution) parameters, and confront them with (visco-)resistive magnetohydrodynamic simulations performed with very different, but equally widely used discretization schemes. This will allow us to comment on the obtained evolutions in a manner irrespective of the adopted discretization strategy. Our findings demonstrate that all schemes used (finite volume based shock-capturing, high order finite differences, and particle in cell-like methods) qualitatively agree on the various evolutionary stages, and that resistivity values of order 0.001 already can lead to chaotic island appearance. However, none of the methods exploited demonstrates convergence in the strong sense in these chaotic regimes. At the same time, nonperturbed tests for showing convergence over long time scales in ideal to resistive regimes are provided as well, where all methods are shown to agree. Both the advantages and disadvantages of specific discretizations as applied to this challenging problem are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

The effect of weak collisionality on damped modes and its contribution to linear mode coupling in gyrokinetic simulation

We revisit the characteristics of stable, damped modes originating from the Landau damping by employing a discretized gyrokinetic Vlasov simulation and also eigenvalue analysis in an unsheared slab geometry. By comparing results between gyrokinetic simulation and an eigenvalue analysis, we found that there exists a critical collisionality β{sub c}{sup *} at which the Case-van Kampen (CvK) modes are damped down to the analytically estimated Landau damping rate and an eigenmode consistent with Landau's theory emerges. Consequently, the recurrence phenomenon disappears so that the Landau damping can be properly reproduced. The critical collisionality β{sub c}{sup *} depends on the resolution in velocity space; i.e., a higher (lower) resolution requires a lower (higher) collisionality, while tends to zero (β{sub c}{sup *}→0) as Δv→0. It is found through a reduced model that even in the collisionless case with marginally stable CvK modes, the linear mode coupling between unstable and stable/damped components through a tertiary mode and the resultant energy transfer can be properly calculated such that the stable/damped mode persists as an eigenstate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Simulation of linear and nonlinear Landau damping of lower hybrid waves

The linear physics of lower hybrid waves (LHWs) and their nonlinear interaction with particles through Landau damping are studied with the gyrokinetic electron and fully kinetic ion (GeFi) particle simulation model in the electrostatic limit. Unlike most other wave modes, the LHWs can resonantly interact with both electrons and ions, with the former being highly magnetized and latter nearly unmagnetized around the lower hybrid frequency. Direct interactions of LHWs with electrons and/or ions are investigated for cases with various k{sub ∥}/k,T{sub i}/T{sub e}, and wave amplitudes. In the linear electron Landau damping (ELD), the dispersion relation and the linear damping rate obtained from our simulation agree well with the analytical linear theory. As the wave amplitude increases, the nonlinear Landau effects are present, and a transition from strong decay at smaller amplitudes to weak decay at larger amplitudes is observed. In the nonlinear stage, the LHWs in the long time evolution finally exhibit a steady Bernstein-Greene-Kruskal mode, in which the wave amplitude is saturated above the noise level. While the resonant electrons are trapped in the wave field in the nonlinear ELD, the resonant ions are untrapped in the LHW time scales. The ion Landau damping is thus predominantly in a linear fashion, leading to a wave saturation level significantly lower than that in the ELD. On the long time scales, however, the ions are still weakly trapped. The results show a coupling between the LHW frequency and the ion cyclotron frequency during the long-time LHW evolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Two-dimensional simulations of thermonuclear burn in ignition-scale inertial confinement fusion targets under compressed axial magnetic fields

We report for the first time on full 2-D radiation-hydrodynamic implosion simulations that explore the impact of highly compressed imposed magnetic fields on the ignition and burn of perturbed spherical implosions of ignition-scale cryogenic capsules. Using perturbations that highly convolute the cold fuel boundary of the hotspot and prevent ignition without applied fields, we impose initial axial seed fields of 20–100 T (potentially attainable using present experimental methods) that compress to greater than 4 × 10{sup 4} T (400 MG) under implosion, thereby relaxing hotspot areal densities and pressures required for ignition and propagating burn by ∼50%. The compressed field is high enough to suppress transverse electron heat conduction, and to allow alphas to couple energy into the hotspot even when highly deformed by large low-mode amplitudes. This might permit the recovery of ignition, or at least significant alpha particle heating, in submarginal capsules that would otherwise fail because of adverse hydrodynamic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Transition characteristics and electron kinetics in microhollow cathode discharges

We demonstrate the transition characteristics and electron kinetics of microdischarges in a microgap with a cathode having microstructures using two-dimensional particle-in-cell/Monte Carlo collision (2d3v) simulations. It is found that the microdischarge is inside the hollow cavity at lower pressures, forming microhollow cathode discharges (MHCDs), while the dominant discharge moves outside the cavity at higher pressures, locating above the cathode rectangular protrusion tip. The spatial distributions of the microdischarge parameters (e.g., electric potential and electron density) are presented by capturing the transition characteristics. The electron kinetics of the microdischarges are examined based on the moment analysis of the electron Boltzmann equation. The collisional and non-collisional components of the electron power absorption are decomposed, which were seldomly investigated for microdischarges previously. The results provide fundamental insights into MHCD formation with a structured electrode at varying pressure conditions, which could be beneficial for the design and fabrication of microplasma devices in practical applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Effects of electron recirculation on a hard x-ray source observed during the interaction of a high intensity laser pulse with thin Au targets

The interaction of a high intensity laser pulse on the preplasma of a high-Z solid target produced by the pulse's pedestal generates high-energy electrons. These electrons subsequently penetrate inside the solid target and produce bremsstrahlung photons, generating an x-ray source which can be used for photonuclear studies or to radiograph high area density objects. The source characteristics are compared for targets with thin (20 μm) and thick (100 μm) Au foils on the Omega EP laser at Laboratory for Laser Energetics. Simulations using the particle-in-cell code CALDER show that for a 20 μm thickness Au target, electrons perform multiple round-trips in the target under the effect of the laser ponderomotive potential and the target electrostatic potential. These relativistic electrons have random transverse displacements, with respect to the target normal, attributed to electrostatic fluctuation fields. As a result, the x-ray spot size is increased by a factor 2 for thin target compared to thick targets, in agreement with experimental results. In addition, the computed doses agree with the measured ones provided that electron recirculation in the thin target is taken into account. A dose increase by a factor 1.7 is then computed by allowing for recirculation. In the 100 μm target case, on the other hand, this effect is found to be negligible.

07 ISOTOPE AND RADIATION SOURCES↗