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,099 records · Page 61

SUPPRESSION OF ELECTRON THERMAL CONDUCTION IN THE HIGH β INTRACLUSTER MEDIUM OF GALAXY CLUSTERS

Understanding the thermodynamic state of the hot intracluster medium (ICM) in a galaxy cluster requires knowledge of the plasma transport processes, especially thermal conduction. The basic physics of thermal conduction in plasmas with ICM-like conditions has yet to be elucidated, however. We use particle-in-cell simulations and analytic models to explore the dynamics of an ICM-like plasma (with small gyroradius, large mean free path, and strongly sub-dominant magnetic pressure) driven by the diffusive heat flux associated with thermal conduction. Linear theory reveals that whistler waves are driven unstable by electron heat flux, even when the heat flux is weak. The resonant interaction of electrons with these waves then plays a critical role in scattering electrons and suppressing the heat flux. In a 1D model where only whistler modes that are parallel to the magnetic field are captured, the only resonant electrons are moving in the opposite direction to the heat flux, and the electron heat flux suppression is small. In 2D or more, oblique whistler modes also resonate with electrons moving in the direction of the heat flux. The overlap of resonances leads to effective symmetrization of the electron distribution function and a strong suppression of heat flux. The results suggest that thermal conduction in the ICM might be strongly suppressed, possibly to negligible levels.

79 ASTRONOMY AND ASTROPHYSICS↗

Anomalous hydrodynamics effective actions from holography

We derive an effective action for charged plasmas with an anomalous (abelian) global current charge current using holography. The holographic description is captured by the dynamics of an Einstein-Maxwell-Chern-Simons theory in an asymptotically AdS spacetime. The ’t Hooft anomaly contribution, which is encoded in the Chern-Simons term, contributes at the Gaussian order in the effective action only in the momentum diffusion sector, where it leads to chiral shear waves. However, as we demonstrate, beyond the Gaussian order, there is non-trivial imprint of the parity-violating anomaly term in sound and charge diffusion dynamics as well.

79 ASTRONOMY AND ASTROPHYSICS↗

First-principles equation-of-state table of silicon and its effects on high-energy-density plasma simulations

Using density-functional theory–based molecular-dynamics simulations, we have investigated the equation of state for silicon in a wide range of plasma density and temperature conditions of ρ=0.001–500g/cm 3 and T=2000–10 8 K. With these calculations, we have established a first-principles equation-of-state (FPEOS) table of silicon for high-energy-density (HED) plasma simulations. When compared with the widely used SESAME-EOS model (Table 3810), we find that the FPEOS-predicted Hugoniot is ~20% softer; for off-Hugoniot plasma conditions, the pressure and internal energy in FPEOS are lower than those of SESAME EOS for temperatures above T ≈ 1–10 eV (depending on density), while the former becomes higher in the low- T regime. The pressure difference between FPEOS and SESAME 3810 can reach to ~50%, especially in the warm-dense-matter regime. Implementing the FPEOS table of silicon into our hydrocodes, we have studied its effects on Si-target implosions. When compared with the one-dimensional radiation-hydrodynamics simulation using the SESAME 3810 EOS model, the FPEOS simulation showed that (1) the shock speed in silicon is ~10% slower; (2) the peak density of an in-flight Si shell during implosion is ~20% higher than the SESAME 3810 simulation; (3) the maximum density reached in the FPEOS simulation is ~40% higher at the peak compression; and (4) the final areal density and neutron yield are, respectively, ~30% and ~70% higher predicted by FPEOS versus the traditional simulation using SESAME 3810. All of these features can be attributed to the larger compressibility of silicon predicted by FPEOS. Furthermore, these results indicate that an accurate EOS table, like the FPEOS presented here, could be essential for the precise design of targets for HED experiments.

36 MATERIALS SCIENCE↗

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma Imaging, LOcal Measurement, and Tomographic Experiment (PILOT): A Mission Concept for Transformational Multi-Scale Observations of Mass and Energy Flow Dynamics in Earth’s Magnetosphere

We currently do not understand the fundamental physical processes that govern mass and energy flow through the Earth’s magnetosphere. Knowledge of these processes is critical to understanding the mass loss rate of Earth’s atmosphere, as well as for determining the role that a planetary magnetic field plays in atmospheric retention, and therefore habitability, for Earth-like planets beyond the solar system. Mass and energy flow processes are challenging to determine at Earth in part because Earth’s planetary magnetic field creates a complex “system of systems” composed of interdependent plasma populations and overlapping spatial regions that perpetually exchange mass and energy across a broad range of temporal and spatial scales. Further, the primary mass carrier in the magnetosphere is cold plasma (as cold as ~0.1 eV), which is invisible to many space-borne instruments that operate in the inner magnetosphere. The Plasma Imaging LOcal and Tomographic experiment (PILOT) mission concept, described here, provides the transformational multi-scale observations required to answer fundamental open questions about mass and energy flow dynamics in the Earth’s magnetosphere. PILOT uses a constellation of spacecraft to make radio tomographic, remote sensing, and in-situ measurements simultaneously, fully capturing cold plasma mass dynamics and its impact on magnetospheric systems over an unprecedented range of spatial and temporal scales. This article details the scientific motivation for the PILOT mission concept as well as a potential mission implementation.

79 ASTRONOMY AND ASTROPHYSICS↗

Ionic transport in high-energy-density matter

Ionic transport coefficients for dense plasmas have been numerically computed using an effective Boltzmann approach. Here, we developed a simplified effective potential approach that yields accurate fits for all of the relevant cross sections and collision integrals. These results have been validated with molecular-dynamics simulations for self-diffusion, interdiffusion, viscosity, and thermal conductivity. Molecular dynamics has also been used to examine the underlying assumptions of the Boltzmann approach through a categorization of behaviors of the velocity autocorrelation function in the Yukawa phase diagram. By using a velocity-dependent screening model, we examine the role of dynamical screening in transport. Implications of these results for Coulomb logarithm approaches are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improved modeling of the solid-to-plasma transition of polystyrene ablator for laser direct-drive inertial confinement fusion hydrocodes

We report that the target performance of laser direct-drive inertial confinement fusion (ICF) can be limited by the development of hydrodynamic instabilities resulting from the nonhomegeneous laser absorption at the target surface, i.e., the laser imprint on the ablator. To understand and describe the formation of these instabilities, the early ablator evolution during the laser irradiation should be considered. In this work, an improved modeling of the solid-to-plasma transition of a polystyrene ablator for laser direct-drive ICF is proposed. This model is devoted to be implemented in hydrocodes dedicated to ICF which generally assume an initial plasma state. The present approach consists of the two-temperature model coupled to the electron, ion and neutral dynamics including the chemical fragmentation of polystyrene. The solid-to-plasma transition is shown to significantly influence the temporal evolution of both free electron density and temperatures, which can lead to different shock formation and propagation compared with an initial plasma state. The influence of the solid-to-plasma transition on the shock dynamics is evidenced by considering the scaling law of the pressure with respect to the laser intensity. The ablator transition is shown to modify the scaling law exponent compared with an initial plasma state.

36 MATERIALS SCIENCE↗

Analysis of double stub tuner control stability in a phased array antenna with strong cross-coupling

Active stub tuning with a fast ferrite tuner (FFT) has greatly increased the ability of fusion ion cyclotron range of frequency (ICRF) systems (50–100 MHz) by allowing for the antenna system to respond dynamically to changes in the plasma load impedance such as during the L–H transition or edge localized modes (ELMs). A high power waveguide double-stub tuner is under development for use with the Alcator C-Mod lower hybrid current drive (LHCD) system at 4.6 GHz. The amplitude and relative phase shift between adjacent columns of an LHCD antenna are critical for control of the launched n || spectrum. Adding a double-stub tuning network will perturb the phase and amplitude of the forward wave particularly if the unmatched reflection coefficient is high. This effect can be compensated by adjusting the phase of the low power microwave drive for each klystron amplifier. Cross-coupling of the reflected power between columns of the launcher must also be considered. The problem is simulated by cascading a scattering matrix for the plasma provided by a linear coupling model with the measured launcher scattering matrix and that of the FFTs. The solution is advanced in an iterative manner similar to the time-dependent behavior of the real system. System performance is presented under a range of edge density conditions from under-dense to over-dense and a range of launched n || . Simulations suggest power reflection coefficients (Γ 2 ) of less than 1% with no contamination of the n || spectrum. Instability of the FFT tuning network can be difficult for certain plasma conditions and relative phasings, but reducing the control gain of the FFT network stabilizes the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Heavy and light flavor jet quenching in heavy-ion collisions in a perturbative QCD approach

In this work, we study the flavor dependence of jet quenching in heavy-ion collisions. A next-to-leading-order perturbative QCD framework is applied to account for both quark and gluon contributions to light and heavy flavor hadron productions. A linear Boltzmann transport model coupled to hydrodynamic simulation is utilized to study the nuclear modification of heavy and light flavor jets in the dynamically evolving quark-gluon plasma. Within our jet quenching framework, we obtain a nice description of the nuclear modification factors of charged hadrons, D mesons, B mesons and B–decayed D mesons over a wide range of transverse momentum (8-300 GeV). Our result shows that perturbative QCD is sufficient to explain the color, mass and energy dependence of parton energy loss and jet quenching in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calculation of collisionless pitch-angle scattering of runaway electrons with synchrotron radiation via high-order guiding-centre equation

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Finally, test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimization of laser energy deposition for single-shot high aspect-ratio microstructuring of thick BK7 glass

We investigate the generation of high aspect ratio microstructures across 0.7 mm thick glass by means of single shot Bessel beam laser direct writing. We study the effect on the photoinscription of the cone angle, as well as of the energy and duration of the ultrashort laser pulse. The aim of the study is to optimize the parameters for the writing of a regular microstructure due to index modification along the whole sample thickness. By using a spectrally resolved single pulse transmission diagnostics at the output surface of the glass, we correlate the single shot material modification with observations of the absorption in different portions of the retrieved spectra, and with the absence or presence of spectral modulation. Numerical simulations of the evolution of the Bessel pulse intensity and of the energy deposition inside the sample help us interpret the experimental results that suggest to use picosecond pulses for an efficient and more regular energy deposition. Picosecond pulses take advantage of nonlinear plasma absorption and avoid temporal dynamics effects which can compromise the stationarity of the Bessel beam propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ponderomotive force driven density modifications parallel to B on the LAPD

The ponderomotive force has previously been identified as a possible driver of observed density modifications close to radio frequency actuators during operation. This nonlinear force redistributes density in regions of gradients in the magnitude of an oscillating electric field and describes the influence of the fast time scale RF wave dynamics on slow time scale plasma transport. Depletion of the saturation current (a proxy for the density) measured at the Large Plasma Device (LAPD) was 30–35% during ion cyclotron range of frequencies operation. A coupled 1D plasma transport and cold plasma frequency domain wave solver was developed to self-consistently describe ponderomotive effects and was used to compare with results obtained from the LAPD experiment. The scaled current density driver for the wave model yielded an RF B field in close agreement with two components of the experimental data. However, the 1D parallel model did not accurately reproduce the amplitude or spatial distribution observed in experimental measurements of By. Within the limitations of the 1D model, initial simulation results showed that the ponderomotive force depleted up to 8% for high power (1 MW) and around 1% for the experimental power of 120 kW. This could suggest that the ponderomotive force is not the main driver of density modification for the LAPD experiments presented in this paper. Higher fidelity tools of at least 2D will be required to give a more realistic description of the RF E fields and the effect of the ponderomotive force on the LAPD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A novel method for thermal noise reduction, enabling measurements of broadband, low-amplitude electron temperature fluctuations using individual radiometer channels

A new analysis method has been developed for measurements of broadband, low-amplitude turbulent electron temperature fluctuations in fusion plasmas using individual radiometer channels of a correlation electron cyclotron emission diagnostic. This method takes advantage of differences in the correlation time of thermal noise compared to the correlation time of plasma fluctuations in fusion reactors. The validation of this single-channel method is demonstrated using comparisons with the standard dual-channel radiometer spectral decorrelation method for measurements of turbulent electron temperature fluctuations in the core and edge of low confinement (L), improved confinement (I), and high confinement (H)-mode plasmas at the ASDEX Upgrade tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Formation of a compound flux rope by the merging of two filament channels, the associated dynamics, and its stability

We present observations of compound flux rope formation, which occurred on 2014 January 1, via merging of two nearby filament channels, the associated dynamics, and its stability using multiwavelength data. We also discuss the dynamics of cool and hot plasma moving along the newly formed compound flux rope. The merging started after the interaction between the southern leg of the northward filament and the northern leg of the southward filament at ≈01:21 UT and continued until a compound flux rope formed at ≈01:33 UT. During the merging, the cool filament plasma heated up and started to move along both sides of the compound flux rope, i.e., toward the north (≈265 km s{sup –1}) and south (≈118 km s{sup –1}) from the point of merging. After traveling a distance of ≈150 Mm toward the north, the plasma cooled down and started to return back to the south (≈14 km s{sup –1}) after ≈02:00 UT. The observations provide a clear example of compound flux rope formation via merging of two different flux ropes and the occurrence of a flare through tether cutting reconnection. However, the compound flux rope remained stable in the corona and had a confined eruption. The coronal magnetic field decay index measurements revealed that both the filaments and the compound flux rope axis lie within the stability domain (decay index <1.5), which may be the possible cause for their stability. The present study also deals with the relationship between the filament's chirality (sinistral) and the helicity (positive) of the surrounding flux rope.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of ELM stability with extended MHD models in JET, JT-60U and future JT-60SA tokamak plasmas

The stability with respect to a peeling–ballooning mode (PBM) was investigated numerically with extended MHD simulation codes in JET, JT-60U and future JT-60SA plasmas. The MINERVA-DI code was used to analyze the linear stability, including the effects of rotation and ion diamagnetic drift (${\omega }_{* {\rm{i}}}$), in JET-ILW and JT-60SA plasmas, and the JOREK code was used to simulate nonlinear dynamics with rotation, viscosity and resistivity in JT-60U plasmas. It was validated quantitatively that the ELM trigger condition in JET-ILW plasmas can be reasonably explained by taking into account both the rotation and ${\omega }_{* {\rm{i}}}$ effects in the numerical analysis. When deuterium poloidal rotation is evaluated based on neoclassical theory, an increase in the effective charge of plasma destabilizes the PBM because of an acceleration of rotation and a decrease in ${\omega }_{* {\rm{i}}}$. The difference in the amount of ELM energy loss in JT-60U plasmas rotating in opposite directions was reproduced qualitatively with JOREK. By comparing the ELM affected areas with linear eigenfunctions, it was confirmed that the difference in the linear stability property, due not to the rotation direction but to the plasma density profile, is thought to be responsible for changing the ELM energy loss just after the ELM crash. A predictive study to determine the pedestal profiles in JT-60SA was performed by updating the EPED1 model to include the rotation and ${\omega }_{* {\rm{i}}}$ effects in the PBM stability analysis. It was shown that the plasma rotation predicted with the neoclassical toroidal viscosity degrades the pedestal performance by about 10% by destabilizing the PBM, but the pressure pedestal height will be high enough to achieve the target parameters required for the ITER-like shape inductive scenario in JT-60SA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Demonstration of TNSA proton radiography on the National Ignition Facility Advanced Radiographic Capability (NIF-ARC) Laser

We report that proton radiography using short-pulse laser drivers is an important tool in high-energy density (HED) science for dynamically diagnosing key characteristics in plasma interactions. Here we detail the first demonstration of target-normal sheath acceleration (TNSA)-based proton radiography the NIF-ARC laser system aided by the use of compound parabolic concentrators (CPCs). The multi-kJ energies available at the NIF-ARC laser allows for a high-brightness proton source for radiography and thus enabling a wide range of applications in HED science. In this demonstration, proton radiography of a physics package was performed and this work details the spectral properties of the TNSA proton probe as well as description of the resulting radiography quality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Suppression of Electron Thermal Conduction by Whistler Turbulence in a Sustained Thermal Gradient

The dynamics of weakly magnetized collisionless plasmas in the presence of an imposed temperature gradient along an ambient magnetic field is explored with particle-in-cell simulations and modeling. Two thermal reservoirs at different temperatures drive an electron heat flux that destabilizes off-angle whistler-type modes. The whistlers grow to large amplitude, δB/B 0 ≃ 1, and resonantly scatter the electrons, significantly reducing the heat flux. Interestingly, the resulting steady-state heat flux is largely independent of the thermal gradient. The rate of thermal conduction is instead controlled by the finite propagation speed of the whistlers, which act as mobile scattering centers that convect the thermal energy of the hot reservoir. The findings are relevant to thermal transport in high-β astrophysical plasmas such as hot accretion flows and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗