Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma thermodynamics”

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 253 records · Page 14

Near-field models and simulations of pellet ablation in tokamaks

Here, detailed numerical studies of the ablation of a single neon pellet in the plasma disruption mitigation parameter space have been performed. Simulations were carried out using FronTier, a hydrodynamic/low magnetic Reynolds number MHD code with explicit tracking of material interfaces. FronTier's physics models resolve the pellet surface ablation and the formation of a dense, cold cloud of ablated material, the deposition of energy from hot plasma electrons, the expansion of the ablation cloud along magnetic field lines, and the radiation losses. A local thermodynamic equilibrium model based on Saha equations has been used to resolve atomic processes in the cloud and Redlich–Kwong corrections to the ideal gas equation of state for cold and dense gases, which have been used near the pellet surface. The FronTier pellet code is a next generation of the code described. It has been validated against the semi-analytic improved neutral gas shielding model in the 1D spherically symmetric approximation. Main results include quantification of the influence of atomic processes and Redlich–Kwong corrections on the pellet ablation in spherically symmetric approximation and verification of analytic scaling laws in a broad range of pellet and plasma parameters. Using axially symmetric MHD simulations, the properties of ablation channels and the reduction of pellet ablation rates in magnetic fields of increasing strength have been studied. While the main emphasis has been given to neon pellets for the plasma disruption mitigation, selected results on deuterium fueling pellets have also been presented.

97 MATHEMATICS AND COMPUTING↗

Converging towards atomic and nuclear pressures (Final Report)

his report details the development of transformational tools and techniques, allowing the accurate characterization of matter at the most extreme conditions yet studied in the high energy density (HED) domain. Before this work, most experiments quantitatively mapping the nature of matter in the HED realm were performed using planar geometry, which allows for the isolation and measurement of variables in a single thermodynamic state. Such measurements include equation of state variables, optical conductivity, heat transport and more. However, due to a variety of plasma processes, such experiments are limited to below ~10 TPa (100 Mbar) pressures. To achieve higher pressures, convergent experiments are needed. Historically, convergent experiments have not been used for benchmark data because they are integral measurements. That is each part of a convergent target undergoes a time dependent wide range of states, making it difficult to accurately isolate any particular quantity for a given thermodynamic state. This effort developed innovative convergent HED platforms and techniques enabling the exploration of matter from atomic to nuclear scale pressures. This effort also performed pioneering experiments that yielded the first rigorous benchmark data at these extreme conditions. This funding award began with the overarching goal to understand the behavior of matter at extreme (atomic-to-nuclear scale) pressures. The motivation for these goals lie in the fact that every time scientists explore matter beyond the threshold of an atomic unit, there is a fundamental shift in science. The atomic unit for energy, mass, charge, length, and time have each been explored, and each time such a threshold was crossed, a new sequence of discoveries was made resulting in significant awards such as the Nobel Prize. The only unexplored atomic unit is pressure, and this effort set the course for exploring matter at and beyond such pressures. That most of the recently discovered extrasolar planets and stars, as well as matter in the late implosion stages of inertial fusion targets, have deep internal pressures at and beyond atomic pressures amplifies the importance of this effort. Much of the initial work in developing techniques to create these atomic-to-nuclear pressures already existed within the HED community, predominantly at large scale laser facilities such as the National Ignition Facility (NIF) and the Omega60 laser the University of Rochester and although these experiments were routinely performed the ability to extract information about the underlying physical states and processes remained elusive due to the complexity of the experiments, extreme scales in both time and space, and the integrated nature of all the measurements. This work built a rigorous framework so such measurements can routinely be made. This was done in part through the introduction of Bayesian inference techniques into the field of HED science. These techniques allow for the self-consistent extraction of the relevant variables and their explicit and implicit correlations, so as to make use of integrated experimental data to constrain physical models and provide accurate uncertainty bounds for the data. The techniques developed here are fully transparent and proved immediately useful providing quantitative rigorous benchmarks for physical states and processes at some of the most extreme conditions yet explored on Earth. This funding is directly or partly responsible for 7 publications in major peer-reviewed scientific journals, a doctoral thesis, 3 invited talks at major conferences, and the training of a post-doc, 2 graduate students, and 1 undergraduate student. Beyond the effort supported by this award launched the introduction of Bayesian inference into the HED physics community and helped push the usage if modern data-science techniques within the physics community at large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗

Thermal transport modeling of laser-irradiated spheres

Thermal transport of uniformly laser-irradiated spheres of various materials is investigated computationally. One-dimensional simulations of low- to mid-Z materials (Be, Al, and Cu) are performed to evaluate the impact of nonlocal electron transport on experimental observables under laser intensities of relevance to direct-drive inertial confinement fusion. Here we compare thermal transport models of different levels of fidelity: flux-limited Spitzer–Harm diffusion, the Schurtz–Nicolai–Busquet (SNB) reduced-order nonlocal model, and a Fokker–Planck description. Spitzer–Harm diffusion with different flux-limiter factors are compared with different implementations of the SNB model in the HYDRA radiation hydrodynamics code. Under the conditions of interest, the peak heat flux in the thermal front with the SNB model shows good agreement with Fokker–Planck calculations, with the largest errors below 10% at 10 15 W/cm 2 laser intensity. From HYDRA-SNB simulations, two experimentally relevant effects are observed from nonlocal heat transport when compared to flux-limited Spitzer–Harm modeling: coronal temperatures are cooler due to reduced heat fluxes in the expanding plasma and (for mid-Z materials) x-ray emissions are enhanced due to preheating in the dense plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic field effects on laser energy deposition and filamentation in magneto-inertial fusion relevant plasmas

We report on experimental measurements of how an externally imposed magnetic field affects plasma heating by kJ-class, nanosecond laser pulses. The experiments reported here took place in gas cells analogous to magnetized liner inertial fusion targets. We observed significant changes in laser propagation and energy deposition scale lengths when a 12T external magnetic field was imposed in the gas cell. We find evidence that the axial magnetic field reduces radial electron thermal transport, narrows the width of the heated plasma, and increases the axial plasma length. Reduced thermal conductivity increases radial thermal gradients. This enhances radial hydrodynamic expansion and subsequent thermal self-focusing. Our experiments and supporting 3D simulations in helium demonstrate that magnetization leads to higher thermal gradients, higher peak temperatures, more rapid blast wave development, and beam focusing with an applied field of 12T.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling of tungsten filament in gas discharge in H - ion source

Sources of negative ions such as H - are essential elements of proton accelerators and tokamaks. They have limited lifetime. The replacement of an ion source is a costly process causing delays in the operation of the entire machine. The hot tungsten filament is the key element that limits the lifetime of the H - ion source at the Los Alamos Neutron Science Center (LANSCE) facility. An accurate model, which describes the filament physics and provides a reliable estimate on its lifetime, may significantly improve the operation of the facility. In this report we describe such a model, including a comprehensive list of relevant physical processes from first principles. The model can be used to describe different regimes of operation using diagnostics data as input parameters. It has been benchmarked against the data collected during the production cycle at LANSCE and shows good agreement with experimental data.

47 OTHER INSTRUMENTATION↗

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗

An empirical polytrope law for solar wind thermal electrons between 0.45 and 4.76 AU - Voyager 2 and Mariner 10

In this paper empirical evidence is presented that between 0.4 and 5 AU the thermal portion (but not all) of the solar wind electron population obeys a polytrope relation. It is also shown that this functional relationship is a member of a broader class of possible laws required of a steady state, fully ionized plasma whose proper frame electric field is dominated by the polarization electric field. The empirically determined, thermodynamically interesting value of the polytrope index (1.175) is virtually that predicted (1.16) by the theoretical considerations of Scudder and Olbert (1979). Strong, direct, empirical evidence for the nearly isothermal behavior of solar wind electrons as has been indirectly argued in the literature for some time is provided.

Sittler, E. C., Jr.↗

Temperature Inhibition of Plasma-Driven Methane Conversion in DBD Systems

Low-temperature non-thermal plasmas produce highly reactive chemical environments made up of electrons, ions, radicals, and vibrationally excited molecules. These reactive species, when combined with catalysts, can help drive thermodynamically unfavorable chemical reactions at low temperatures and atmospheric pressure. The conversion of methane (CH 4 ) to produce other value-added chemicals is a good model system because of its applicability to a wide range of industries. To effectively create these plasma catalytic systems, a fundamental understanding of the plasma-phase chemistry alone is imperative. While there have been many studies on methane plasmas and how certain operating conditions (i.e., gas composition and power) affect the plasma, there is limited understanding on how changing bulk reaction temperature affects the plasma properties and ensuing plasma chemistry. Here, in this work, we use a dielectric barrier discharge to investigate the effects of temperature on the reaction chemistry and the plasma’s electrical properties in various methane-gas mixtures. Results show that increasing temperature leads to a reduction in methane conversion as well as changes to both the gas and dielectric material pre-breakdown, which manifests itself in temperature-dependent electrical properties of the plasma. Experiments at various temperatures and power show a positive correlation between key electrical plasma properties (average charge and lifetime per filament) and the measured methane conversion as a function of temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectral deconvolution without the deconvolution: Extracting temperature from x-ray Thomson scattering spectra without the source-and-instrument function

X-ray Thomson scattering (XRTS) probes the dynamic structure factor of the system, but the measured spectrum is broadened by the combined source-and-instrument function (SIF) of the setup. In order to extract properties such as temperature from an XRTS spectrum, the broadening by the SIF needs to be removed. Recent work [Dornheim et al. Nat. Commun. 13 , 7911 (2022)] has suggested that the SIF may be deconvolved using the two-sided Laplace transform. However, the extracted information can depend strongly on the shape of the input SIF, and the SIF is in practice challenging to measure accurately. Here, we propose an alternative approach: we demonstrate that considering ratios of Laplace-transformed XRTS spectra collected at different scattering angles is equivalent to performing the deconvolution, but without the need for explicit knowledge of the SIF. From these ratios, it is possible to directly extract the temperature from the scattering spectra, when the system is in thermal equilibrium. We find the method to be generally robust to spectral noise and physical differences between the spectrometers, and we explore situations in which the method breaks down. Furthermore, the fact that consistent temperatures can be extracted for systems in thermal equilibrium indicates that non-equilibrium effects could be identified by inconsistent temperatures of a few eV between the ratios of three or more scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermodynamic effect of the ion sound instability in the ionosphere

During geomagnetic disturbances when the ring current interacts intensely with the plasmasphere, the plasma of this region undergoes a strong heating due to an ion cyclotron instability. This is followed by the transfer of heat along geomagnetic field lines from the heating region to the ionosphere. One of the results of this process is the formation of non-isothermal region (in which T(sub e) greater than 3.4 T(sub i) at ionospheric heights) caused by a rapid cooling the H(+) ions due to their resonant charge exchange with neutral hydrogen. Heat transfer from the top of the flux tube to the ionosphere is investigated using a hydrodynamic model for the ionosphere-plasmasphere coupling. Field-aligned currents, present in the topside ionosphere, are often accompanied by ion sound turbulence. The turbulence scatters electrons, increasing the total electron collision frequency through wave-particle effects. The influence of wave-particle interactions introduces an anomalous component to the total collision frequency, which modifies substantially the heat conduction coefficient of the plasma. As a result, the plasma is heated more intensely above than below this region of ion sound turbulence.

Khazanov, G. V.↗

Transport coefficients of the quark-gluon plasma at the critical point and across the first-order line

A bottom-up Einstein-Maxwell-dilaton holographic model is used to compute, for the first time, the behavior of several transport coefficients of the hot and baryon-rich strongly coupled quark-gluon plasma at the critical point and also across the first-order phase transition line in the phase diagram. The observables under study are the shear and bulk viscosities, the baryon and thermal conductivities, the baryon diffusion, the jet quenching parameter $\hat{q}$, as well as the heavy-quark drag force and the Langevin diffusion coefficients. These calculations provide a phenomenologically promising estimate for these coefficients, given that our model quantitatively reproduces lattice QCD thermodynamics results, both at zero and finite baryon density, besides naturally incorporating the nearly perfect fluidity of the quark-gluon plasma. Here we find that the diffusion of baryon charge, and also the shear and bulk viscosities, are suppressed with increasing baryon density, indicating that the medium becomes even closer to perfect fluidity at large densities. On the other hand, the jet quenching parameter and the heavy-quark momentum diffusion are enhanced with increasing density. The observables display a discontinuity gap when crossing the first-order phase transition line, while developing an infinite slope at the critical point. The transition temperatures associated with different transport coefficients differ in the crossover region but are found to converge at the critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamical control in a prethermalized molecular ultracold plasma: Local dissipation drives global relaxation

Prethermalization occurs as an important phase in the dynamics of isolated many-body systems when coupled degrees of freedom relax in a subspace separated from a state of complete thermodynamic equilibrium by a gap in energy or other conserved quantity. Slow equilibration from a prethermal phase can localize the dynamics of natural and model systems despite high dimensionality and limited disorder. Here, we report the signature of an enduring prethermal regime of arrested relaxation in the molecular ultracold plasma that forms following the avalanche of a state-selected Rydberg gas of nitric oxide. For a wide range of initial conditions, this system enters a critical phase in which a density of NO + and electrons balances a population of Rydberg molecules. Electron collisions mix orbital angular momentum, scattering Rydberg electrons to states of very high ℓ . The rapid predissociation of molecules in states of low- ℓ purifies this high- ℓ ensemble. The angular momentum barrier that separates NO + ions and electrons creates an extraordinary gap between the plasma states of n ≈ ℓ , with measured n > 200 and penetrating states of ℓ = 0 , 1 , and 2. Evolution to a statistically equilibrated state of N and O atoms cannot occur without Rydberg electron penetration, and this gap blocks relaxation for a millisecond or more. A weak radiofrequency (RF) field drives ℓ -mixing electron collisions that erase this gap, causing wholesale dissipation. Remarkably, a local quantum-state transition promoting an exceedingly small fraction of the molecules in the prethermalized ensemble to a predissociative state also acts with a global effect, driving the entire system to a dissipative equilibrium. Using the Lindblad master equation, we illustrate qualitatively similar dynamics for a toy model of an open quantum system that consists of a localized set of spins on which dissipation acts at a single site as a gateway to equilibrium. Published by the American Physical Society 2025

Wang, Ruoxi↗

A novel design of transitional layer structure between reduced activation ferritic martensitic steels and tungsten for plasma facing materials

Plasma-facing components (PFCs) are among the most critical gaps for fusion energy to establish technical and economic feasibility. Tungsten as a first wall/blanket material in PFCs requires to be integrating with reduced activation ferritic martensitic (RAFM) steels as a structural component. Currently, major drawbacks are the requirement of brazing, the formation of a brittle interface, and a large difference between the coefficients of thermal expansion of tungsten and steel. Here, a novel transitional multilayer structure was designed and investigated to join tungsten and RAFM steels using three interlayers. The composition of each interlayer was selected based on computational thermodynamics and diffusion kinetics to ensure a body-centered cubic (bcc) single-phase structure and prevent the formation of a brittle intermetallic phase region in the temperature range of 600–1150 °C. Although the transitional layer structure was designed for additive manufacturing, spark plasma sintering (SPS) as proof of concept was used to bond the individual layers. Interfaces were investigated using scanning and transmission electron microscopy methods but no layered intermetallic phase was observed. Nanoindentation maps across the interface suggest major hardness differences at the interface between tungsten and the vanadium interlayer, as well as the interface between RAFM steel and the FeCrAl interlayer.

36 MATERIALS SCIENCE↗

Electron-scale temperature gradients in kinetic equilibrium: MMS observations and Vlasov–Maxwell solutions

For an ideal gas, the notion of thermodynamic equilibrium requires that the system’s temperature be spatially uniform. We demonstrate that for a collisionless plasma, kinetic equilibrium can be achieved even in the presence of a temperature gradient. Motivated by recent MMS magnetopause observations of thin current layers exhibiting electron-scale temperature gradients, here we present an exact solution to the Vlasov-Maxwell system which accounts for the observed electron temperature gradient and reproduces much of the plasma’s spatial variation along one dimension. Furthermore, the current layer in the Vlasov-Maxwell model is self-consistently supported by a prominent crescent shaped electron velocity distribution that naturally arises due to the temperature transition, which agrees with the electron distribution structures measured by MMS during two asymmetric current layers encountered at the magnetopause.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonequilibrium radiative heating during outer planet atmospheric entry

The contradictory results obtained by investigators assessing the influence of finite-rate ionization on the radiative heating of probes entering the atmospheres of Jupiter and Saturn are discussed. Tiwari and Szema (1979) found that the radiation heating increased, whereas Leibowitz (1973) and Liebowitz and Kuo (1976) found that it decreased in relation to the results obtained when the ionization rate was assumed to be in equilibrium at the local thermodynamic conditions. The study presented here is limited to stagnation shock layers for nonviscous, hydrogen-helium plasmas with cold, nonblowing conditions at the probe heat shield. It is found that the radiative heating comes mainly from the Balmer region of the spectrum, where the shock layer is optically thin.

Nelson, H. F.↗

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗