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 1,783 records · Page 99

Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP

Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The stability of expanding reactive shocks in a van der Waals fluid

Despite the extensive literature accumulated since the pioneering works of D'yakov and Kontorovich in the 1950s, the stability of steady shocks is still an open question when realistic boundary conditions are accounted. The consideration of a supporting mechanism, which is indeed a necessary condition for shock steadiness, modifies the perturbation shock dynamics in the unstable range. The Noh problem is a suitable example to form steady expanding shocks. This configuration is of great interest to the high-energy-density-physics community because of its direct application to inertial confinement fusion and astrophysics, for which the stagnation of a supersonically converging material via an accretion shock front is ubiquitous. Here in this work, we extend the generalized Noh problem, both base-flow solution and linear stability analysis, to conditions where endothermic or exothermic transformations undergo across the shock. Within the spontaneous acoustic emission conditions found for a van der Waals gas, we find that cylindrical and spherical expanding shocks become literally unstable for sufficiently high mode numbers. Counterintuitively, the effect of exothermicity or endothermicity across the shock is found to be stabilizing or destabilizing, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics and symmetries of quarks as probes to phases of QCD. Final Technical Report

In the period of August 1, 2017-July 31, 2019, supported by the DOE Grant No. DE-SC0018209, we have performed the comprehensive theoretical research on broad aspects of QCD plasma that is created in relativistic heavy-ion collisions, with particular emphasis on the dynamics and symmetries of quarks in the hot and dense QCD matter. Our research has also extended to the study of dynamical aspects of possible QCD critical point, that has been searched for experimentally in the Beam Energy Scan (BES) program currently going on in BNL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced Dual Confocal Measurement System

A measurement instrument utilizing dual, chromatic, confocal, distance sensors has been jointly developed by General Atomics (GA) and Sandia National Laboratories (SNL) for thickness and flatness measurement of target components used in dynamic materials properties (DMP) experiments on the SNL Z-Machine (Z). Compared to previous methods used in production of these types of targets, the tool saves time and yields a 4x reduction in thickness uncertainty which is one of the largest sources of error in equation of state (EOS) measurements critical to supporting the NNSA’s Stockpile Stewardship program and computer modeling of high energy density experiments. It has numerous differences from earlier instruments operating on the dual confocal sensor principle to accommodate DMP components including larger lateral travel, longer working distance, ability to measure flatness in addition to thickness, built-in thickness calibration standards for quickly checking calibration before and after each measurement, and streamlined operation. Thickness and flatness of 0.2mm-3.3mm thick sections of diamond machined copper and aluminum can be measured to “sub-micron” accuracy. Sections up to 6mm thick can be measured with as-yet undermined accuracy. Furthermore, samples must have one surface which is flat to within 300µm, lateral dimensions of no more than 50mm x 50mm, and height less than 40mm.

36 MATERIALS SCIENCE↗

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗

Perspectives on relativistic electron–positron pair plasma experiments of astrophysical relevance using high-power lasers

The study of relativistic electron–positron pair plasmas is both of fundamental physics interest and important to understand the processes that shape the magnetic field dynamics, particle acceleration, and radiation emission in high-energy astrophysical environments. Although it is highly desirable to study relativistic pair plasmas in the laboratory, their generation and control constitutes a critical challenge. Significant experimental and theoretical progress has been made over recent years to explore the use of intense lasers to produce dense relativistic pair plasma in the laboratory and study the basic collective plasma processes associated with these systems. Important challenges remain in terms of improving the number of pairs, system size, and control over the charge neutrality required to establish laboratory platforms that can expand our understanding of relativistic pair plasma and help validate underlying models in conditions relevant to high-energy astrophysical phenomena. Furthermore, we highlight recent progress in this field, discuss the main challenges, and the exciting prospects for studying relativistic pair plasmas and astrophysics relevant instabilities in the laboratory in the near future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Latitude variation of recurrent MeV-energy proton flux enhancements in the heliocentric radial range 11 to 20 AU and possible correlation with solar coronal hole dynamics

Recurrent low energy (not less than 0.5 MeV) proton flux enhancements, reliable indicators of corotating plasma interaction regions in interplanetary space, have been observed on the Voyager 1 and 2 and Pioneer 11 spacecraft in the heliographic latitude range 2 deg S to 23 deg N and the heliocentric radial range 11 to 20 AU. After a period of rather high correlation between fluxes at different latitudes in early 1983, distinct differences develop. The evolution of the fluxes appears to be related to the temporal and latitudinal dynamics of solar coronal holes, suggesting that information about the latitudinal structure of solar wind stream sources propagates to these distances.

Christon, S. P.↗

Plasma waves and continuum radiation in planetary magnetospheres

Voyager data on whistler mode waves and electron cyclotron harmonic emissions are analyzed to understand the interaction of the waves with the dynamics of the electrons. The occurrence and characteristics of Jovian whistler mode chorus and the interactions with the plasma in and near the Io torus are emphasized. Bernstein waves, especially those near the upper hybrid or plasma frequency are discussed to provide insight into an important plasma diagnostic tool and to compare the relevant portions of the electron distribution function at Earth, Jupiter, and Saturn. The nonthermal continuum radiation common to the magnetospheres of Earth, Jupiter, and Saturn is considered. Because of the very low frequency of these radio waves and their close association with upper hybrid resonance emissions, continuum radiation is often associated more closely with the plasma wave spectrum of a planetary magnetosphere than with the planet's radio spectrum.

Kurth, W. S.↗

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Contribution of low-energy ionospheric protons to the plasma sheet

The magnetospheric transport of low-energy ionospheric ions is examined by means of three-dimensional particle codes. Emphasis is placed on the behavior of polar wind and cleft originating protons. It is demonstrated that, via nonadiabatic motion inside the neutral sheet, these ions can significantly contribute to the populations of the plasma sheet. The importance of this contribution is found to depend critically upon the dynamics of particles originating from the highest latitudes, as these possibly have access to the distant tail. Hence it is shown that polar wind H(+) expelled into the magnetosphere at very low energies (in the electron volt range) preferentially feed the plasma sheet during quiet times, experiencing accelerations up to several kiloelectron volts upon return into the inner magnetosphere. In contrast, during disturbed times, the intensifying magnetospheric convection confines this population to low L shells where it travels in a nearly adiabatic manner. As for the protons originating from the cleft fountain, the simulations reveal that they can be transported up to the vicinity of the distant neutral line in the nightside sector. Via interaction with the neutral sheet, these ionospheric ions are rapidly raised to the characteristic plasma sheet energy range. The density levels contributed by these populations are quite substantial when compared to those measured in situ. These simulations establish an active role of low-energy ionospheric ions in the overall magnetospheric dynamics.

Delcourt, D. C.↗

Dependence of emittance on the length of an electron bunch during laser-plasma acceleration in guiding structures

We report a theoretical analysis and numerical simulation of the dynamics of transverse emittance of an electron bunch during its acceleration in wake fields generated by a laser pulse in a weakly nonlinear mode. Analytical expressions are obtained for the main factors affecting the emittance growth during acceleration and the case is considered when the characteristic transverse size of the injected bunch exceeds the matched radius determined by the focusing force at the injection point, the initial emittance and the electron bunch energy, and the resulting value of emittance is much larger than the initial one. The dynamics of the emittance growth during acceleration as a function of the length of the electron bunch is described, and the length of the bunch is found at which there occurs a complete phase mixing of betatron oscillations of electrons and the emittance increases to its maximum value determined by the bunch parameters and the focusing force at the injection point. Analytical expressions are in good agreement with the results of numerical simulation. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High heat-flux response of high-conductivity graphitic foam monoblocks

Plasma-facing components based on the so-called monoblock design are planned for use in the divertor region of long-pulse plasma devices such as ITER and JT-60SA due to their capacity to handle high heat fluxes with active water cooling. The plasma-facing materials that are preferred for these monoblocks are tungsten for ITER or carbon-carbon fiber reinforced carbon (CFC) for JT-60SA. The requirements for the plasma-facing components include the ability to handle high plasma fluxes resulting in high temperatures. Therefore, reasonably high thermal conductivity is required. In this study, high thermal conductivity graphite foam is explored as a monoblock material. Four different test articles are created, two with the foam brazed to CuCrZr tubes, and two with the foam press fit to the tubes. These monoblocks are then thermally cycled at 8 MW/m 2 in the high heat-flux test facility GLADIS to examine their robustness for long-pulse divertor applications. Experimental results are compared with computational fluid dynamics simulation. Graphite foam shows promise for use in a plasma-facing component, but further development is necessary to address the significant drop in thermal conductivity at high temperatures.

Youchison, Dennis↗

High-energy Radiation and Ion Acceleration in Three-dimensional Relativistic Magnetic Reconnection with Strong Synchrotron Cooling

Abstract We present the results of 3D particle-in-cell simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of nonradiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, ≫ σ , and form a hard power-law energy distribution, f i ∝ γ −1 . This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to either σ in the strong cooling regime or the radiation burnoff limit, γ syn , when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit, ε c ≈ 16 MeV, is only efficiently produced in the strong cooling regime, γ syn < σ . In this regime, we find that the spectral cutoff scales as ε cut ≈ ε c ( σ / γ syn ) and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectra of Magnetic Turbulence in a Relativistic Plasma

Abstract We present a phenomenological and numerical study of strong Alfvénic turbulence in a magnetically dominated collisionless relativistic plasma with a strong background magnetic field. In contrast with the nonrelativistic case, the energy in such turbulence is contained in magnetic and electric fluctuations. We argue that such turbulence is analogous to turbulence in a strongly magnetized nonrelativistic plasma in the regime of broken quasi-neutrality. Our 2D particle-in-cell numerical simulations of turbulence in a relativistic pair plasma find that the spectrum of the total energy has the scaling k −3/2 , while the difference between the magnetic and electric energies, the so-called residual energy, has the scaling k −2.4 . The electric and magnetic fluctuations at scale ℓ exhibit dynamic alignment with the alignment angle scaling close to cos ϕ ℓ ∝ ℓ 1 / 4 . At scales smaller than the (relativistic) plasma inertial scale, the energy spectrum of relativistic inertial Alfvén turbulence steepens to k −3.5 .

79 ASTRONOMY AND ASTROPHYSICS↗

Heavy flavor dynamics across system size at the LHC

One of the fundamental signatures of the Quark Gluon Plasma has been the suppression of heavy flavor (specifically D mesons), which has been measured via the nuclear modification factor, R AA and azimuthal anisotropies, v n , in large systems. However, multiple competing models can reproduce the same data for R AA to v n . In this talk we break down the competing effects that conspire together to successfully reproduce R AA and v n in experimental data using Trento+v-USPhydro+DAB-MOD. Then using our best fit model we make predictions for R AA and v n across system size for 208 PbPb, 129 XeXe, 40 ArAr, and 16 OO collisions. We find that 0–10% centrality has a non-trivial interplay between the system size and eccentricities such that system size effects are masked in v 2 whereas in 30–50% centrality the eccentricities are approximately constant across system size and, therefore, is a better centrality class to study D meson dynamics across system size.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗