Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two-temperature”

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.

53 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Closure theory for high-collisionality multi-ion plasmas

A general formalism is developed to construct and solve a system of linearized moment equations for parallel and perpendicular closures in high-collisionality plasmas. It is applicable for multiple ion species with arbitrary masses, temperatures, charges, and densities. The convergence of closure coefficients is evaluated by increasing the number of moments from 2 to 32 for scalar, vector, and rank-2 tensor moments. As an example, the complete set of closure coefficients for a deuterium-carbon plasma over the entire Hall parameter range is presented. Furthermore, the closure coefficients at various temperature ratios show that the one-temperature closure coefficients can differ significantly from the two-temperature coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Spectroscopic diagnostics of high-temperature plasma in stellar coronae using Fe XXIV-XXIV K-shell lines with XRISM

The RS CVn type binary star GT Muscae was observed during its quiescence using the Resolve X-ray microcalorimeter spectrometer onboard XRISM. The main and satellite lines of the Fe XXIV-XXIV K-shell transitions were resolved for the first time from stellar sources. We conducted line ratio analysis to investigate any deviations from collisional ionization equilibrium and Maxwell electron energy distribution with a single temperature. By using five combinations of direct excitation lines and dielectronic recombination satellite lines in three line complexes (Fe He$\alpha$, Ly$\alpha$, and He$\beta$), we found that the plasma is well characterized by two-temperature thermal plasmas with temperatures of 1.7 and 4.3 keV, which is consistent with thermal broadening of Fe xxv, and the broad-band fitting results in the 1.7–10 keV band. Other forms of deviation from a single-temperature plasma, such as different ionization and electron temperatures or the $\kappa$ distribution for the electron energy distributions, are not favored, which is reasonable for stellar coronae at quiescence. This study demonstrates the utility of the Fe K-shell line ratio diagnostics to probe plasma conditions using X-ray microcalorimeters.

X-rays: stars↗
DOE OSTI2021

Role of element-specific damping in ultrafast, helicity-independent, all-optical switching dynamics in amorphous (Gd,Tb)Co thin films

Ultrafast control of the magnetization in ps timescales by fs laser pulses offers an attractive avenue for applications such as fast magnetic devices for logic and memory. However, ultrafast helicity-independent all-optical switching (HI-AOS) of the magnetization has thus far only been observed in Gd-based, ferrimagnetic amorphous (a-) rare earth-transition metal (a-RE-TM) systems, and a comprehensive understanding of the reversal mechanism remains elusive. Here, we report HI-AOS in ferrimagnetic a-Gd 22-x Tb x Co 78 thin films, from x=0 to 18, and elucidate the role of Gd in HI-AOS in a-RE-TM alloys and multilayers. Increasing Tb content results in increasing perpendicular magnetic anisotropy and coercivity, without modifying magnetization density, and slower remagnetization rates and higher critical fluences for switching but still shows picosecond HI-AOS. Simulations of the atomistic spin dynamics based on the two-temperature model reproduce these results qualitatively and predict that the lower damping on the RE sublattice arising from the small spin-orbit coupling of Gd (with L=0) is instrumental for the faster dynamics and lower critical fluences of the Gd-rich alloys. Annealing a-Gd 10 Tb 12 Co 78 leads to slower dynamics which we argue is due to an increase in damping. These simulations strongly indicate that accounting for element-specific damping is crucial in understanding HI-AOS phenomena. The results suggest that engineering the element-specific damping of materials can open up new classes of materials that exhibit low-energy, ultrafast HI-AOS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Calculation of electron-ion temperature equilibration rates and friction coefficients in plasmas and liquid metals using quantum molecular dynamics

In this work, we discuss a method to calculate with quantum molecular dynamics simulations the rate of energy exchanges between electrons and ions in two-temperature plasmas, liquid metals, and hot solids. Promising results from this method were recently reported for various materials and physical conditions. Like other ab initio calculations, the approach offers a very useful comparison with the experimental measurements and permits an extension into conditions not covered by the experiments. The energy relaxation rate is related to the friction coefficients felt by individual ions due to their nonadiabatic interactions with electrons. Each coefficient satisfies a Kubo relation given by the time integral of the autocorrelation function of the interaction force between an ion and the electrons. These Kubo relations are evaluated using the output of quantum molecular dynamics calculations in which electrons are treated in the framework of finite-temperature density functional theory. The calculation presents difficulties that are unlike those encountered with the Kubo formulas for the electrical and thermal conductivities. In particular, the widely used Kubo-Greenwood approximation is inapplicable here. Indeed, the friction coefficients and the energy relaxation rate diverge in this approximation since it does not properly account for the electronic screening of electron-ion interactions. The inclusion of screening effects considerably complicates the calculations. We discuss the physically motivated approximations we applied to deal with these complications in order to investigate a widest range of materials and physical conditions. Unlike the standard method used for the electronic conductivities, the Kubo formulas are evaluated directly in the time domain and not in the energy domain, which spares one from needing to introduce an extraneous undetermined numerical parameter to account for the discrete character of the numerical density of states. We highlight interesting properties of the energy relaxation rate not shared by other electronic properties, in particular its self-averaging character. We then present a detailed parametric and convergence study with the numerical parameters, including the system size, the number of bands and k-points, and the physical approximations for the dielectric function and the exchange-correlation energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

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↗
DOE OSTI2023

Collisional-energy-cascade model for nonthermal velocity distributions of neutral atoms in plasmas

Nonthermal velocity distributions with much greater tails than a Maxwellian have been observed for radical atoms in plasmas for a long time. Historically, such velocity distributions have been modeled by a two-temperature Maxwell distribution. In this paper, I propose a model based on collisional energy cascade, which has been studied in the field of granular materials. In the collisional energy cascade, a particle ensemble undergoes energy input at the high-energy region, entropy production by elastic collisions among particles, and energy dissipation. For radical atoms, energy input may be caused by the Franck-Condon energy of molecular dissociation or charge-exchange collision with hot ions, and the input energy is eventually dissipated by collisions with the walls. Here, I show that the steady-state velocity distribution in the collisional energy cascade is approximated by the generalized Mittag-Leffler distribution, which is a one-parameter extension of the Maxwell distribution. This parameter indicates the degree of the nonthermality and is related to the relative importance of energy dissipation over entropy production. This model is compared with a direct molecular dynamics simulation for simplified gaseous systems with energy input, as well as some experimentally observed velocity distributions of light radicals in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermal effects mediating the flow induced by laser-induced optical breakdown

Here, a short focused laser pulse can generate a hot plasma, which expands then contracts and can eject a hot jet, the speed and direction of which is sensitive to the details of the plasma kernel. The coupling of thermal and chemical nonequilibrium plasma mechanisms with this hydrodynamic development is assessed with detailed simulations of a two-temperature, three-species plasma model. Time scales for electron recombination, thermal relaxation, and diffusion are compared to that of the plasma expansion to anticipate conditions in which these mechanisms might affect the vorticity generation that leads to the ultimate flow pattern. The effect of these mechanisms are analyzed through comparison with corresponding inert-gas and equilibrium models. Thermal-nonequilibrium effects are found to be weak due to rapid relaxation of the heavy-particle and electronic temperatures. In contrast, chemical equilibration occurs at a rate comparable to the expansion and thereby enhance both it and subsequent hydrodynamic mechanisms as the energy stored in ion formation during the preceding breakdown is released by electron recombination. Thermal conduction, enhanced by high-temperature free electrons, weakens the ejection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗
DOE OSTI2026

Influence of conical wire array geometry on plasma flow and temperature profiles of radiatively cooled jets

The influence of conical wire array geometry on the formation and dynamics of pulsed-power driven plasma jets is investigated. In the experiments, the jet becomes isolated from the inflows as it passes through an aperture, allowing the study of its intrinsic evolution for different conical angles. Here, our results show that, regardless of the array opening angle, the jets are supersonic, highly collisional, and exhibit an exponential axial density decay with a characteristic scale length of 𝐿 𝑛 ≈ 3 mm, significantly shorter than the overall length of the jet. In contrast, axial velocity systematically increases with larger array opening angles. The near invariance of the density profile is attributed to a compensating mechanism between geometric divergence and axial acceleration, consistent with an asymptotically steady flow regime. Additionally, temperature measurements reveal ion-electron thermal decoupling near the base of the jet, with equilibration downstream. Analysis of radiative cooling and collisional energy exchange timescales indicates that the plasma evolves in a radiatively stable regime and that the two-temperature profile is fully described by the competition between both methods. These findings provide valuable information into the internal structure and evolution of pulsed-power plasma jets in regimes dominated by geometric shaping and radiative cooling, demonstrating experimental control over flow acceleration and collimation relevant to both laboratory and astrophysical contexts.

Physics - Plasma physics↗
DOE OSTI2021

Combined Mixed Potential Electrochemical Sensors and Artificial Neural Networks for the Quantificationand Identification of Methane in Natural Gas Emissions Monitoring

Sensors capable of quantifying methane concentration and discriminating between possible sources are needed for natural gas leak detection where multiple spatially overlapping sources including wetlands and agriculture may be present. We report on the fabrication by an additive manufacturing process of a four electrode La 0.87 Sr 0.13 CrO 3 , Indium Tin Oxide (In 2 O 3 90 wt%, SnO 2 10 wt%), Au, Pt mixed potential electrochemical sensor using yttria-stabilized zirconia (YSZ) as a solid electrolyte to natural gas detection. Artificial neural networks (ANNs) are used to automatically decode the possible source and concentration of methane. The ANNs trained on sensor data are capable of correctly discriminating between three sources of methane emissions from simulated mixtures of emissions from cattle, wetlands, or natural gas with >98% accuracy. Quantification error for methane in mixtures of CH 4 in air, CH 4 + NH3 in air, and simulated natural gas is less than 1.5% ppm when a two-temperature dataset is employed.

03 NATURAL GAS↗
DOE OSTI2026

CHEQUP v0.1

CHEQUP (Castro-based Hofi Expansion with QUasineutral Plasma) is a simulation code for modeling the formation of hydrodynamic optical-field-ionized (HOFI) plasma channels, which are used as waveguides in laser-plasma acceleration experiments. This includes experiments performed at LBNL's BELLA facility as well as other laser facilities across the world. CHEQUP extends the open-source Castro hydrodynamics framework with physics modules tailored for modeling HOFI plasma channels -- including multi-species ionization and three-body recombination for mixtures of hydrogen, nitrogen, helium, and argon ; a two-temperature model tracking electron and heavy-species temperatures separately ; and coupling with other codes of the BLAST ecosystem (https://blast.lbl.gov/) such as WarpX, via the openPMD standard. CHEQUP inherits from Castro the ability to run on modern GPU architectures (NVIDIA CUDA, AMD HIP) and supports adaptive mesh refinement (AMR) for efficient multi-scale resolution. Compared to existing tools, CHEQUP would be, to our knowledge, the first open-source code implementing the full HOFI channel formation physics, and the first implementation capable of running on GPUs. This enables significantly faster, large-scale parameter scans critical for the design of next-generation LPA-based accelerators and light sources.

Lehe, Remi [Lawrence Berkeley National Laboratory ↗
DOE OSTI2026

Variation in surface properties as error sources in dynamic temperature measurements

Most temperature measurement methods rely upon the assumption of a single uniform source as part of the analysis in inferring an observed temperature T from the surface of a material. This is true for the radiometric and spectroscopic methods we rely upon in shock physics. However, there are obvious cases where this is a poor assumption that can cause large errors in analysis, such as the measurement of emission from a surface composed of more than one material, such as a rusted surface of Fe and Fe 2 O 3 , or a multi-material surface that could be produced through alloying methods (e.g., mokume gane) or 3D printing. Here we perform an initial examination of the effects of varied surface properties upon the analysis using idealized test data to constrain when we may have better confidence in the resulting calculations. We also examine the structures of the fits to test whether, in principle, comparing the fit with ideal gray bodies would reveal such a structure and enable the researcher to identify when a two-temperature fitting model should be used instead of a single- T model or a wavelength-dependent emissivity model.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Global X-Ray Properties of the Vela and Puppis A Supernova Remnants

The Vela and Puppis A supernova remnants (SNRs) comprise a large emission region of ∼8° diameter in the soft X-ray sky. The HaloSat CubeSat mission provides the first soft X-ray (0.4–7 keV) observation of the entire Vela SNR and Puppis A SNR region with a single pointing and moderate spectral resolution. HaloSat observations of the Vela SNR are best fit with a two-temperature thermal plasma model consisting of a cooler component with kT{sub 1}=0.19{sub −0.01}{sup +0.01} keV in collisional ionization equilibrium and a hotter component with kT{sub 2}=1.06{sub −0.27}{sup +0.45} keV in nonequilibrium ionization. Observations of the Puppis A SNR are best fit with a single-component plane-parallel shocked plasma model with kT=0.86{sub −0.05}{sup +0.06} keV in nonequilibrium ionization. For the first time, we find the total X-ray luminosities of both components of the Vela SNR spectrum in the 0.5–7 keV energy band to be L{sub X}=4.4{sub −1.4}{sup +1.4}×10{sup 34} erg s{sup −1} for the cooler component and L{sub X}=4.1{sub −1.5}{sup +1.8}×10{sup 34} erg s{sup −1} for the hotter component. We find the total X-ray luminosities of the Vela and Puppis A SNRs to be L{sub X}=8.4×10{sup 34} erg s{sup −1} and L{sub X}=6.7{sub −0.9}{sup +1.1}×10{sup 36} erg s{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

NuSTAR Observations of Abell 2163: Constraints on Non-thermal Emission

Since the first non-thermal reports of inverse Compton (IC) emission from the intracluster medium (ICM) of galaxy clusters at hard X-ray energies, we have yet to unambiguously confirm IC emission in observations with newer facilities. RXTE detected IC emission in one of the hottest known clusters, Abell 2163 (A2163), a massive merging cluster with a giant radio halo—the presumed source of relativistic electrons IC scattering CMB photons to X-ray energies. The cluster's redshift (z ~ 0.2) allows its thermal and non-thermal radio emission to fit the field of view of the Nuclear Spectroscopic Telescope Array (NuSTAR), permitting a deep observation capable of confirming or ruling out the RXTE report. The IC flux provides constraints on the average magnetic field strength in a cluster. To determine the global diffuse IC emission in A2163, we fit its global NuSTAR spectrum with four models: single- (1T) and two-temperature (2T), 1T+power-law component (T+IC), and multi-temperature+power law (9T+IC). Each represents different characterizations of the thermal ICM emission, with power-law components added to represent IC emission. We find that the 3–30 keV spectrum can be described by purely thermal emission, with a global average temperature of kT = (11.8 ± 0.2) keV. The IC flux is constrained to <4.0 × 10 –12 erg s –1 cm –2 using the 1T+IC model and <1.6 × 10 –12 erg s –1 cm –2 with the more physical 9T+IC model, both to 90% confidence levels. Here, combining these limits with 1.4 GHz diffuse radio data from the VLA, we find the average magnetic field strength to be >0.22 μG and >0.35 μG, respectively, providing the strongest constraints on these values in A2163 to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

XRISM Reveals Low Nonthermal Pressure in the Core of the Hot, Relaxed Galaxy Cluster A2029

We present XRISM Resolve observations of the core of the hot, relaxed galaxy cluster Abell 2029 (A2029). We find that the line-of-sight bulk velocity of the intracluster medium (ICM) within the central 180 kpc is at rest with respect to the brightest cluster galaxy, with a 3σ upper limit of ∣vbulk∣ < 100 km s−1. We robustly measure the field-integrated ICM velocity dispersion to be σv = 169 ± 10 km s−1, obtaining similar results for both single-temperature and two-temperature plasma models to account for the cluster cool core. This result, if ascribed to isotropic turbulence, implies a subsonic ICM with Mach number M 3D ≈ 0.22 and a nonthermal pressure fraction of 2.6 ± 0.3%. The turbulent velocity is similar to what was measured in the core of the Perseus cluster by Hitomi, but here in a more massive cluster with an ICM temperature of 7 keV, the limit on the nonthermal pressure fraction is even more stringent. Our result is consistent with expectations from simulations of relaxed clusters, but it is on the low end of the predicted distribution, indicating that A2029 is an exceptionally relaxed cluster with no significant impacts from either a recent minor merger or active galactic nucleus activity.

Audard, Marc (ORCID:000000034721034X)↗