Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA DYNAMICS”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Theory of gradient drift instabilities in low-temperature, partially magnetised plasmas

A fluid dispersion theory in partially magnetised plasmas is analysed to examine the conditions under which large-wavelength modes develop in Penning-type configurations, that is, where an electric field is imposed perpendicular to a homogeneous magnetic field. The fluid dispersion relation assuming a slab geometry shows that two types of low-frequency, gradient drift instabilities occur in the direction of the E×B and diamagnetic drifts. One type of instability, observed when the equilibrium electric field and plasma density gradient are in the same direction, is similar to the classic modified Simon–Hoh instability. A second instability is found for conditions in which (i) the diamagnetic drift is in the direction opposite to the E×B drift and (ii) the magnitude of the diamagnetic drift is sufficiently larger than the electron thermal speed. The present fluid dispersion theory suggests that the rotating spokes driven by such fluid instabilities propagate in the same direction as the diamagnetic drift, which can be in the same direction as or opposite to the E×B drift, depending on the plasma conditions. This finding may account for the observation, in some plasma devices, of the rotation of large-scale structures in both the E×B and -E×B directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Validation of magnetized gas-jet experiments to investigate the effects of an external magnetic field on laser-plasma instabilities

Laser-plasma instabilities (LPI) play a detrimental role in energy coupling to the target in inertial confinement fusion (ICF). The recent development of applied strong magnetic fields for use in ICF and laboratory astrophysics experiments has opened opportunities to investigate the role of external magnetic fields on LPIs. Recent numerical studies have shown that stimulated Raman scattering (SRS) can be mitigated by external magnetic fields in the kinetic regime of the instability and warrant systematic experimental studies to validate modelling. To this end, we design experiments at the OMEGA-EP laser facility to investigate the effect of an external perpendicular B-field of 5–30 T on the backscattered light spectrum from a gas-jet target. We present measurements from a proof-of-principle experiment, where the backward-SRS (BSRS) is in the kinetic regime, for which the magnetization is expected to produce the greatest effects on instability growth. New simulations motivated by this experiment are used to inform the proposal of an upgraded experimental design. Our simulation predictions show that the new design is suited to experimentally demonstrating BSRS mitigation by an external magnetic field, despite the magnetization effects on the hydrodynamics, as well as the inherent temporal and spatial variations in plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quasilinear theory of Brillouin resonances in rotating magnetized plasmas

Both spin and orbital angular momentum can be exchanged between a rotating wave and a rotating magnetized plasma. Through resonances the spin and orbital angular momentum of the wave can be coupled to both the cyclotron rotation and the drift rotation of the particles. It is, however, shown that the Landau and cyclotron resonance conditions which classically describe resonant energy–momentum exchange between waves and particles are no longer valid in a rotating magnetized plasma column. In this case a new resonance condition which involves a resonant matching between the wave frequency, the cyclotron frequency modified by inertial effects and the harmonics of the guiding centre rotation is identified. A new quasilinear equation describing orbital and spin angular momentum exchanges through these new Brillouin resonances is then derived, and used to expose the wave-driven radial current responsible for angular momentum absorption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A kinetic model for pure electron plasma compression via the rotating wall technique

An approximate model for pure electron plasma compression is developed for the case where the rotating wall (RW) electric field couples to the ${E}\times{B}$rotation and axial bounce motion of the electrons. The key assumption in the model is that, throughout the compression, the plasma remains in a slowly evolving thermal equilibrium defined by the plasma temperature and angular momentum. Linearized drift kinetic theory is employed to derive an expression for torque exerted by the RW field on the plasma through coupling to the resonant plasma particles, and averaging is used to find the torque that both compresses and heats the plasma. The evolution equations for the angular velocity and temperature of the plasma include the compression and heating from the torque and cooling from cyclotron radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma potential shaping using end-electrodes in the Large Plasma Device

We perform experiments in the Large Plasma Device at the University of California, Los Angeles, studying how different end-electrode biasing schemes modify the radial potential profile in the machine. We impose biasing profiles of different polarities and gradient signs on a set of five concentric electrodes placed 12 m downstream from the plasma source. We find that imposing concave-down profiles (negative potential radial gradient) on the electrodes creates radial potential profiles halfway up the plasma column that are comparable to those imposed on the electrodes and a few electron temperature in height, regardless of the biasing polarity. On the other hand, imposing concave-up profiles (positive potential radial gradient) leads to non-monotonic radial potential profiles. This observation can be explained by the current drawn through the electrodes and the parallel plasma resistivity, highlighting their important role in controlling the rotation of plasma. Concave-down plasma potential profiles, obtained by drawing electrons on the axis, are predicted to drive azimuthal drift velocities that can approach significant fractions of the ion sound speed in the central region of the plasma column.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
DOE OSTI2020

Gyrokinetic particle simulations of interactions between energetic particles and magnetic islands induced by neoclassical tearing modes

Interactions between energetic particles (EPs) and neoclassical tearing mode (NTM) islands in the DIII-D tokamak are studied using the global gyrokinetic toroidal code (GTC). GTC simulations find that the EP radial profile is partially flattened within the magnetic island regions and that there are stochastic regions in the particle phase space. Radial particle flux is induced mainly around the magnetic island regions and decreases with time to almost zero when the initial EP distribution achieves a new steady-state in the absence of EP sources. Stochastic regions of magnetic field lines induced by the superposition of multiple islands have weak effects on the particle flux when the width of stochastic regions is smaller than the EP drift orbit width. In conclusion, the perturbed parallel EP current induced by the magnetic islands has weak stabilizing effects on the linear growth rate of the NTM instability in this DIII-D experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hot electrons between cold walls

We consider electron cooling in a collisionless plasma slab between two cold and freely emitting walls. Numerical calculations suggest a counterintuitive behavior of this system: the cooling rate slows down and eventually stops, leaving the system with a significant fraction of its initial thermal energy. Analytical treatment within the Vlasov–Poisson model reveals a set of steady-states with a two-component distribution of electrons: the primary electrons trapped within the potential wells and the secondary electrons forming the counterstreaming beams. We show that such steady-states are linearly stable with respect to one-dimensional perturbations. Establishment of a particular steady-state depends on initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Neural network representability of fully ionized plasma fluid model closures

The closure problem in fluid modeling is a well-known challenge to modelers aiming to accurately describe their systems of interest. Over many years, analytic formulations in a wide range of regimes have been presented but a practical, generalized fluid closure for magnetized plasmas remains an elusive goal. Here, as a first step toward constructing a novel data-based approach to this problem, we apply ever-maturing machine learning methods to assess the capability of neural network architectures to reproduce crucial physics inherent in known magnetized plasma closures. We find encouraging results, indicating the applicability of neural networks to closure physics, but also arrive at recommendations on how one should choose appropriate network architectures for the given locality properties dictated by the underlying physics of the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Laser-heated capillary discharge waveguides as tunable structures for laser-plasma acceleration

Laser-heated capillary discharge waveguides are novel, low plasma density guiding structures able to guide intense laser pulses over many diffraction lengths and have recently enabled the acceleration of electrons to 7.8 GeV by using a laser-plasma accelerator (LPA). These devices represent an improvement over conventional capillary discharge waveguides, as the channel matched spot size and plasma density can be tuned independently of the capillary radius. This has allowed the guiding of petawatt-scale pulses focused to small spot sizes within large diameter capillaries, preventing laser damage of the capillary structure. High performance channel-guided LPAs require control of matched spot size and density, which experiments and simulations reported here show can be tuned over a wide range via initial discharge and laser parameters. In this paper, measurements of the matched spot size and plasma density in laser-heated capillary discharges are presented, which are found to be in excellent agreement with simulations performed using the MHD code MARPLE. Strategies for optimizing accelerator performance are identified based on these results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Laser transport and backscatter in low-density SiO 2 and Ta2O 5 foams

Experiments using a single 527 nm wavelength beam interacting with sub- and supercritical density SiO 2 and Ta 2 O 5 foams examined laser propagation and backscatter from laser–plasma instabilities such as Stimulated Brillouin Scattering (SBS). Two densities of each material were examined, and multiple diagnostics were used to characterize the propagation and backscatter. For 5 mg/cc SiO 2 (n e /n c = 0.375), the laser propagation distance was well approximated by treating the foam as a gas. However, for the 2 mg/cc SiO 2 foam (n e /n c = 0.15), the same model over-predicts the propagation distance by ~40%. Furthermore, existing analytical theories on propagation through subcritical foams were able to account for this difference. The laser heat wave propagated ~1/2 as far in Ta 2 O 5 than SiO 2 foams with similar electron density. We showed that this difference is due to the increased radiation losses in the higher Z foam. The fraction of backscattered light scales linearly with incident laser intensity for the range of intensities examined. Ta 2 O 5 foams had significantly lower levels of backscatter (1–3%) than the SiO 2 (4–8%), which is consistent with estimates of large Landau damping due to the presence of the oxygen atoms. The measured fraction of SBS backscattered laser energy for a 2 mg/cc SiO 2 foam shot was ~4 times lower than predicted by simulations assuming a gas-like foam. We found that we needed to assume increased ion heating such that T i /T e ~1.2–1.5 in the plasma to agree with the measured SBS reflectivity. Analytical models of laser-heated foams predict preferential heating of the ions as has been observed in previous experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Microtearing instabilities and electron thermal transport in low and high collisionality NSTX discharges

Microtearing mode (MTM) real frequency, growth rate, magnetic fluctuation amplitude and resulting electron thermal transport are studied in systematic NSTX scans of relevant plasma parameters. Here, the dependency of the MTM real frequency and growth rate on plasma parameters, suitable for low and high collision NSTX discharges, is obtained by using the reduced MTM transport model. The plasma parameter dependencies are compared and found to be consistent with the results obtained from MTM using the Gyrokinetic GYRO code. The scaling trend of collision frequency and plasma beta is found to be consistent with the global energy confinement trend observed in the NSTX experiment. The strength of the magnetic fluctuation is found to be consistent with the gyrokinetic estimate. In earlier studies, it was found that the version of the Multi-Mode (MM) anomalous transport model, which did not contain the effect of MTMs, provided an appropriate description of the electron temperature profiles in standard tokamak discharges and not in spherical tokamaks. When the MM model, which involves transport associated with MTMs, is incorporated in the TRANSP code and is used in the study of electron thermal transport in NSTX discharges, it is observed that the agreement with the experimental electron temperature profile is substantially improved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

ITER cold VDEs in the limit of perfectly conducting walls

Recently, it has been shown that a vertical displacement event (VDE) can occur in ITER even when the walls are perfect conductors, as a consequence of the current quench (CQ) [A. H. Boozer, Phys. Plasmas 26, 114501 (2019)]. We used the extended-MHD code M3D-C1 with an ITER-like equilibrium and induced a CQ to explore cold VDEs in the limit of perfectly conducting walls, using different wall geometries. In the case of a rectangular first wall with the side walls far away from the plasma, we obtained very good agreement with the analytical model developed by Boozer that considers a top/bottom flat-plates wall. We show that the solution in which the plasma remains at the initial equilibrium position is improved when bringing the side walls closer to the plasma. When approximating the ITER first wall as a perfect conductor, the plasma remains stable at the initial equilibrium position far beyond the value predicted by the flat-plates wall limit. When considering an opposite limit in which only the inner shell of the ITER vacuum vessel acts as a perfect conductor, the plasma is displaced during the CQ, but the edge safety factor remains above 2 longer in the current decay compared to the flat-plates wall limit. In all the simulated cases, the VDE is found to be strongly dependent on the plasma current, in agreement with a similar finding in the flat-plates wall limit, showing an important difference with hot VDEs in which the CQ is not a necessary condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Example of exponentially enhanced magnetic reconnection driven by a spatially bounded and laminar ideal flow

In laboratory and natural plasmas of practical interest, the spatial scale Δ d at which magnetic field lines lose distinguishability differs enormously from the scale a of magnetic reconnection across the field lines. In the solar corona, plasma resistivity gives a/Δ d ~ 10 12 , which is the magnetic Reynold number R m . The traditional resolution of the paradox of disparate scales is for the current density j associated with the reconnecting field Brec to be concentrated by a factor of Rm by the ideal evolution, so j ~ B rec /µ 0 Δ d . A second resolution is for the ideal evolution to increase the ratio of the maximum to minimum separation between pairs of arbitrarily chosen magnetic field lines, Δ max /Δ min , when calculated at various points in time. Reconnection becomes inevitable where Δ max /Δ min ~ R m . A simple model of the solar corona will be used for a numerical illustration that the natural rate of increase in time is linear for the current density but exponential for Δ max /Δ min . Reconnection occurs on a time scale and with a current density enhanced by only ln(a/Δ d ) from the ideal evolution time and from the current density B rec /µ 0 a. In both resolutions, once a sufficiently wide region, Δ r , has undergone reconnection, the magnetic field loses static force balance and evolves on an Alfv´enic time scale. The Alfv´enic evolution is intrinsically ideal but expands the region in which Δ max /Δ min is large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Observations of a plectonemic configuration in a stable magnetized plasma jet

Astrophysical jets are collimated high-speed outflows emerging from spinning and accreting matter around celestial objects and may spontaneously result from self-organized processes. Magnetic self-organization is commonly observed in laboratory plasma physics experiments; however, they require close-fitting flux conservers to constrain and stabilize the toroidal or cylindrical structures. Here we report the first observations of a long, stable, free-boundary plasma jet far from chamber walls, embedding a double-helix magnetic structure resembling a force-free plectonemic Taylor state. The jets arise from an experimental setup that mimics an accretion disk and has no close-fitting solid flux conserver. The results support the hypothesis that self-organization could be a universal, intrinsic explanation for jet formation, collimation, and stability and may help explain double-helix features in celestial observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗