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 469 records · Page 26

Normal stability of slow manifolds in nearly periodic Hamiltonian systems

Kruskal [J. Math. Phys. 3, 806 (1962)] showed that each nearly periodic dynamical system admits a formal U(1) symmetry, generated by the so-called roto-rate. We prove that such systems also admit nearly invariant manifolds of each order, near which rapid oscillations are suppressed. We study the nonlinear normal stability of these slow manifolds for nearly periodic Hamiltonian systems on barely symplectic manifolds—manifolds equipped with closed, non-degenerate 2-forms that may be degenerate to leading order. In particular, we establish a sufficient condition for long-term normal stability based on second derivatives of the well-known adiabatic invariant. We use these results to investigate the problem of embedding guiding center dynamics of a magnetized charged particle as a slow manifold in a nearly periodic system. Here, we prove that one previous embedding and two new embeddings enjoy long-term normal stability and thereby strengthen the theoretical justification for these models.

97 MATHEMATICS AND COMPUTING↗

Efficient quantum algorithm for dissipative nonlinear differential equations

Significance Nonlinear differential equations appear in many domains and are notoriously difficult to solve. Whereas previous quantum algorithms for general nonlinear differential equations have complexity exponential in the evolution time, we give the first quantum algorithm for dissipative nonlinear differential equations that is efficient provided the dissipation is sufficiently strong relative to nonlinear and forcing terms and the solution does not decay too rapidly. We also establish a lower bound showing that differential equations with sufficiently weak dissipation have worst-case complexity exponential in time, giving an almost tight classification of the quantum complexity of simulating nonlinear dynamics. Furthermore, numerical results for the Burgers equation suggest that our algorithm may potentially address complex nonlinear phenomena even in regimes with weaker dissipation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation of complexity dynamics in a DC glow discharge magnetized plasma using recurrence quantification analysis

Recurrence is an ubiquitous feature which provides deep insights into the dynamics of real dynamical systems. A suitable tool for investigating recurrences is recurrence quantification analysis (RQA). It allows, e.g., the detection of regime transitions with respect to varying control parameters. We investigate the complexity of different coexisting nonlinear dynamical regimes of the plasma floating potential fluctuations at different magnetic fields and discharge voltages by using recurrence quantification variables, in particular, DET, L{sub max}, and Entropy. The recurrence analysis reveals that the predictability of the system strongly depends on discharge voltage. Furthermore, the persistent behaviour of the plasma time series is characterized by the Detrended fluctuation analysis technique to explore the complexity in terms of long range correlation. The enhancement of the discharge voltage at constant magnetic field increases the nonlinear correlations; hence, the complexity of the system decreases, which corroborates the RQA analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Picosecond Thermodynamics in Underdense Plasmas Measured with Thomson Scattering

The rapid evolution of electron density and temperature in a laser-produced plasma were measured using collective Thomson scattering. Unprecedented picosecond time resolution, enabled by a pulse-front-tilt compensated spectrometer, revealed a transition in the plasma-wave dynamics from an initially cold, collisional state to a quasi-stationary, collisionless state. The Thomson-scattering spectra were compared with theoretical calculations of the uctuation spectrum using either a conventional Bhatnagar Gross Krook (BGK) collision operator or the rigorous Landau collision terms: the BGK model overestimates the electron temperature by 50% in the most-collisional conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the study of hydrodynamic instabilities in the presence of background magnetic fields in high-energy-density plasmas

Blast-wave-driven hydrodynamic instabilities are studied in the presence of a background B-field through experiments and simulations in the high-energy-density (HED) physics regime. In experiments conducted at the Laboratoire pour l’utilisation des lasers intenses (LULI), a laser-driven shock-tube platform was used to generate a hydrodynamically unstable interface with a prescribed sinusoidal surface perturbation, and short-pulse x-ray radiography was used to characterize the instability growth with and without a 10-T B-field. The LULI experiments were modeled in FLASH using resistive and ideal magnetohydrodynamics (MHD), and comparing the experiments and simulations suggests that the Spitzer model implemented in FLASH is necessary and sufficient for modeling these planar systems. These results suggest insufficient amplification of the seed B-field, due to resistive diffusion, to alter the hydrodynamic behavior. Although the ideal-MHD simulations did not represent the experiments accurately, they suggest that similar HED systems with dynamic plasma-β (=2μ0ρv2/B2) values of less than ∼100 can reduce the growth of blast-wave-driven Rayleigh–Taylor instabilities. These findings validate the resistive-MHD FLASH modeling that is being used to design future experiments for studying B-field effects in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Laser fluence dependence on emission dynamics of ultrafast laser induced copper plasma

The characteristic emission features of a laser-produced plasma strongly depend strongly on the laser fluence. We investigated the spatial and temporal dynamics of neutrals and ions in femtosecond laser (800 nm, ~40 fs, Ti:Sapphire) induced copper plasma in vacuum using both optical emission spectroscopy (OES) and spectrally resolved two-dimensional (2D) imaging methods over a wide fluence range of 0.5-77.5 J/cm 2 . 2D fast gated monochromatic images showed distinct plume splitting between the neutral and ions especially at moderate to higher fluence ranges. OES studies at low to moderate laser fluence regime confirm intense neutral line emission over the ion emission whereas this trend changes at higher laser fluence with dominance of the latter. This evidences a clear change in the physical processes involved in femtosecond laser matter interaction at high input laser intensity. Furthermore, the obtained ion dynamics resulting from the OES, and spectrally resolved 2D imaging are compared with charged particle measurement employing Faraday cup and Langmuir probe and results showed good correlation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

THE CONTRIBUTION OF CORONAL JETS TO THE SOLAR WIND

Transient collimated plasma eruptions in the solar corona, commonly known as coronal (or X-ray) jets, are among the most interesting manifestations of solar activity. It has been suggested that these events contribute to the mass and energy content of the corona and solar wind, but the extent of these contributions remains uncertain. We have recently modeled the formation and evolution of coronal jets using a three-dimensional (3D) magnetohydrodynamic (MHD) code with thermodynamics in a large spherical domain that includes the solar wind. Our model is coupled to 3D MHD flux-emergence simulations, i.e., we use boundary conditions provided by such simulations to drive a time-dependent coronal evolution. The model includes parametric coronal heating, radiative losses, and thermal conduction, which enables us to simulate the dynamics and plasma properties of coronal jets in a more realistic manner than done so far. Here, we employ these simulations to calculate the amount of mass and energy transported by coronal jets into the outer corona and inner heliosphere. Based on observed jet-occurrence rates, we then estimate the total contribution of coronal jets to the mass and energy content of the solar wind to (0.4–3.0)% and (0.3–1.0)%, respectively. Our results are largely consistent with the few previous rough estimates obtained from observations, supporting the conjecture that coronal jets provide only a small amount of mass and energy to the solar wind. We emphasize, however, that more advanced observations and simulations (including parametric studies) are needed to substantiate this conjecture.

79 ASTRONOMY AND ASTROPHYSICS↗

The ion polytropic coefficient in a collisionless sheath containing hot ions

The fluid approach has been widely used to study plasma sheath dynamics. For a sheath containing hot ions whose temperature is greater than the electron's, how to truncate the fluid hierarchy chain equations while retaining to the fullest extent of the kinetic effects is always a difficult problem. In this paper, a one-dimensional, collisionless sheath containing hot ions is studied via particle-in-cell simulations. By analyzing the ion energy equation and taking the kinetic effects into account, we have shown that the ion polytropic coefficient in the vicinity of the sheath edge is approximately constant so that the state equation with the modified polytropic coefficient can be used to close the hierarchy chain of the ion fluid equations. The value of the polytropic coefficient strongly depends on the hot ion temperature and its concentration in the plasma. The semi-analytical model is given to interpret the simulation results. As an application, the kinetic effects on the ion saturation current density in the probe theory are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Multistep and Degenerate Variational Integrators for Applications in Plasma Physics

Geometric integrators yield high-fidelity numerical results by retaining conservation laws in the time advance. A particularly powerful class of geometric integrators is symplectic integrators, which are widely used in orbital mechanics and accelerator physics. An important application presently lacking symplectic integrators is the guiding center motion of magnetized particles represented by non-canonical coordinates. Because guiding center trajectories are foundational to many simulations of magnetically confined plasmas, geometric guiding center algorithms have high potential for impact. The motivation is compounded by the need to simulate long-pulse fusion devices, including ITER, and opportunities in high performance computing, including the use of petascale resources and beyond. This dissertation uses a systematic procedure for constructing geometric integrators - known as variational integration - to deliver new algorithms for guiding center trajectories and other plasma-relevant dynamical systems. These variational integrators are non-trivial because the Lagrangians of interest are degenerate - the Euler-Lagrange equations are first-order differential equations and the Legendre transform is not invertible. The first contribution of this dissertation is that variational integrators for degenerate Lagrangian systems are typically multistep methods. Multistep methods admit parasitic mode instabilities that can ruin the numerical results. These instabilities motivate the second major contribution: degenerate variational integrators. By replicating the degeneracy of the continuous system, degenerate variational integrators avoid parasitic mode instabilities. The new methods are therefore robust geometric integrators for degenerate Lagrangian systems. These developments in variational integration theory culminate in one-step degenerate variational integrators for non-canonical magnetic field line flow and guiding center dynamics. The guiding center integrator assumes coordinates such that one component of the magnetic field is zero; it is shown how to construct such coordinates for nested magnetic surface configurations. Additionally, collisional drag effects are incorporated in the variational guiding center algorithm for the first time, allowing simulation of energetic particle thermalization. Advantages relative to existing canonical-symplectic and non-geometric algorithms are numerically demonstrated. All algorithms have been implemented as part of a modern, parallel, ODE-solving library, suitable for use in high-performance simulations.

Guiding Center↗

Effect of Electronic Inertia on the Gravito-Electrostatic Sheath Structure Formation

The gravito-electrostatic sheath (GES) model, previously proposed to address the fundamental issues on the surface emission mechanism of outflowing solar plasma on the basis of plasma−wall interaction processes with inertialess electrons on both bounded and unbounded scales, is reformulated in the light of active electron inertial response amid geometrical curvature effects. We accordingly derive the electron population distribution law considering both weak electron inertia and geometrical curvature effects in a new analytic construct coupled with the GES structure equations in a closed form. The analysis shows that the GES characteristics and hence plasma outflow dynamics are noticeably affected because of electron inertia. As a consequence of the electron inertia inclusion in contrast with the previous GES formalism, it is found that the GES width gets reduced (–5%), the sheath boundary gets contracted (–7%), the net current density at the surface gets reduced (–25%), the GES potential enhances (+17%), the transonic horizon decreases (‒35%), self-gravity enhances (+2%), and so forth. The obtained results are in fair accord with the existing model predictions centered around both the earlier GES formalisms and standard fluid-kinetic predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optical spectroscopy and modeling of uranium gas-phase oxidation: Progress and perspectives

We report that studies related to U gas-phase oxidation through plasma- and thermo-chemistry are important for many fields, including environmental monitoring, forensic analysis, debris analysis in a weapon detonation event, and nucleation physics. Recently, significant efforts have been made to understand the chemical pathways involved in the progression from U atoms to diatoms (UO) and polyatomic molecules (UxOy), employing optical spectroscopy tools and computational modeling. In many studies, laser ablation of U or a U-containing flow reactor are used as a highly resource-efficient, repeatable, tunable, and lab-scale testbed for studying gas-phase oxidation in U plasmas. The spectroscopic analysis of high-temperature gas-phase oxidation of U is challenging due to the congested U spectra, resolution limitations of instrumentation, and the numerous chemical reaction pathways possible. This article focuses on the current understanding and challenges related to studying U plasma chemistry, specifically U gas-phase oxidation and molecular formation, via optical spectroscopy of plasmas and associated computational and spectral modeling. The physical and chemical processes involved in the evolution from U atoms to U oxide molecules to nanoparticles and agglomerates (i.e., debris) are discussed in the context of optical spectroscopic studies. The article concludes by highlighting opportunities for future research efforts based on existing knowledge published in the literature.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Is there an exact magnetic moment for charged particle motion in a time-dependent, homogeneous magnetic field?

The non-perturbative guiding-centre model provides an exact alternative to full-orbit simulations of charged particle dynamics in situations where traditional guiding-centre theory may fail. We demonstrate that the charged particle motion in a homogeneous, time-varying magnetic field is a solvable example of the non-perturbative guiding-centre model. This entails showing that the exact magnetic moment of Qin and Davidson can be constructed to be asymptotic to the adiabatic invariant series of Kruskal. In contrast to the perturbative invariant, the exact invariant contains information about parametric resonances. These resonances destroy the conservation of the usual magnetic moment over very long times. This refutes some previous claims about the all-time invariance of the magnetic moment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Review of the second charged-particle transport coefficient code comparison workshop

We report the results of the second charged-particle transport coefficient code comparison workshop, which was held in Livermore, California on 24–27 July 2023. This workshop gathered theoretical, computational, and experimental scientists to assess the state of computational and experimental techniques for understanding charged-particle transport coefficients relevant to high-energy-density plasma science. Data for electronic and ionic transport coefficients, namely, the direct current electrical conductivity, electron thermal conductivity, ion shear viscosity, and ion thermal conductivity were computed and compared for multiple plasma conditions. Additional comparisons were carried out for electron–ion properties such as the electron–ion equilibration time and alpha particle stopping power. Overall, 39 participants submitted calculated results from 18 independent approaches, spanning methods from parameterized semi-empirical models to time-dependent density functional theory. In the cases studied here, we find significant differences—several orders of magnitude—between approaches, particularly at lower temperatures, and smaller differences—roughly a factor of five—among first-principles models. We investigate the origins of these differences through comparisons of underlying predictions of ionic and electronic structure. The results of this workshop help to identify plasma conditions where computationally inexpensive approaches are accurate, where computationally expensive models are required, and where experimental measurements will have high impact.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characterization of self-magnetic pinch (SMP) radiographic diode performance on RITS-6 at Sandia National Laboratories. I. Diode dynamics, DC heating to extend radiation pulse

Radiographic diodes focus on an intense electron beam to a small spot size to minimize the source area of energetic photons for radiographic interrogation. The self-magnetic pinch (SMP) diode has been developed as such a source and operated as a load for the six-cavity radiographic integrated test stand (RITS-6) inductive voltage adder driver. While experiments support the generally accepted conclusion that a 1:1 aspect diode (cathode diameter equals anode–cathode gap) delivers optimum SMP performance, such experiments also show that reducing the cathode diameter, while reducing spot size, also results in reduced radiation dose, by as much as 50%, and degraded shot reproducibility. Analysis of the effective electron impingement angle on the anode converter with time made possible by a newly developed dose-rate array diagnostic indicates that fast-developing oscillations of the angle are correlated with early termination of the radiation pulse on many of the smaller-diameter SMP shots. This behavior as a function of relative cathode size persists through experiments with output voltages and currents up to 11.5 MV and 225 kA, respectively, and with spot sizes below approximately few millimeters. Since simulations to date have not predicted such oscillatory behavior, considerable discussion of the angle behavior of SMP shots is made to lend credence to the inference. We report there is clear anecdotal evidence that DC heating of the SMP diode region leads to stabilization of this oscillatory behavior. This is the first of two papers on the performance of the SMP diode on the RITS-6 accelerator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

THE KEY ROLE OF SOLAR DYNAMICS IN THE CHROMOSPHERIC HANLE POLARIZATION

The quantum theory of polarized light allows one to model scattering in the solar atmosphere for inferring its properties. This powerful approach has revealed two key long-standing problems in solar physics: the puzzling dilemmas between theory and observations in several anomalously polarized spectral lines and the need for inferring the ubiquitous weak chromospheric magnetic fields, which requires discriminating the Hanle effect in dynamic optically thick plasmas. However, the ever-present dynamics, i.e., the temporal evolution of heatings and macroscopic motions, has been widely disregarded when modeling and interpreting the scattering polarization. This has hindered a consistent theoretical solution to the puzzle while falsifying the Hanle diagnosis. Here, we show that the dynamical evolution is a keystone for solving both problems because its systematic impact allows an explanation of the observations from “anomalous” instantaneous polarization signals. Evolution accounted for, we reproduce amplitudes and (spectral and spatial) shapes of the Ca i 4227 Å polarization at solar disk center, identifying a restrictive arrangement of magnetic fields, kinematics, heatings, and spatio-temporal resolution. We find that the joint action of dynamics, Hanle effect, and low temporal resolutions mimics Zeeman linear polarization profiles, the true weak-field Zeeman signals being negligible. Our results allow reinterpretation of many polarization signals of the solar spectra and support time-dependent scattering polarization as a powerful tool for deciphering the spatio-temporal distribution of chromospheric heatings and fields. This approach may be a key aid in developing the Hanle diagnosis for the solar atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Apparently superluminal superfluids

We consider the superfluid phase of a specific renormalizable relativistic quantum field theory. We prove that, within the regime of validity of perturbation theory and of the superfluid effective theory, there are consistent and regular vortex solutions where the superfluid’s velocity field as traditionally defined smoothly interpolates between zero and arbitrarily large superluminal values. We show that this solution is free of instabilities and of superluminal excitations. We show that, in contrast, a generic vortex solution for an ordinary fluid does develop an instability if the velocity field becomes superluminal. All this questions the characterization of a superfluid velocity field as the actual velocity of “something”.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Zero modes of massive fermions delocalize from axion strings

Massless chiral excitations can arise from the interactions between a fermion and an axion string, propagating along the string and allowing it to superconduct. The properties of these excitations, or zero modes, dictate how the string interacts with light and can thus have important phenomenological consequences. In this paper, we add a nowhere-vanishing Dirac mass for the fermion in the usual model of axion electrodynamics. We find that the zero modes exhibit an interesting phase structure in which they delocalize from the string’s core as the mass increases, up until a critical value past which they disappear. We study this structure from an analytic perspective, with explicit numerical solutions, and via anomaly inflow arguments. Finally, we derive the two-dimensional effective theory of the zero mode and its interactions with the four-dimensional gauge field and show how this effective theory breaks down as the zero modes delocalize.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗