Search NASA⌕ Search

SEARCH · Search NASA

Results for “Temperature evolution”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

An in situ ambient and cryogenic transmission electron microscopy study of the effects of temperature on dislocation behavior in CrCoNi-based high-entropy alloys with low stacking-fault energy

Temperature is known to affect deformation mechanisms in metallic alloys. As temperature decreases, the stacking-fault energy in many face-centered cubic (fcc) alloys decreases, resulting in a change of deformation mode from dislocation slip to deformation twinning. Such an impact of temperature can be more complex in compositionally heterogeneous microstructures that exhibit, for example, local concentration fluctuation such as that in multi-principal element alloys. In this work, we compare the dislocation behavior and mechanical properties of a fcc Cr 20 Mn 10 Fe 30 Co 30 Ni 10 high-entropy alloy at ambient and liquid-nitrogen temperatures. We find that a network of stacking faults is formed by uniformly extended dislocations at ambient temperatures with low stacking-fault energy, whereas at lower temperatures, uneven dissociation of dislocations becomes significant, which results in severe dislocation pile-ups together with their pronounced entanglement. Our findings indicate that as the stacking-fault energy decreases with decreasing temperature, the heterogeneity of the distribution of elements becomes more dominant in tuning the local variation of lattice resistance. As a result, the change in dislocation behavior at low temperatures strongly affects microstructural evolution and consequently leads to significantly more pronounced work hardening.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of residual oxygen on superconducting niobium films

The integration of niobium (Nb) into emerging superconducting circuits can enhance their performance and function. However, growth of high purity Nb can be challenging due to its high reactivity with oxygen. Here, in this study, we examine the role of residual oxygen inside the growth chamber in transforming the structural, chemical, and superconducting properties of Nb films. We demonstrate that an increase in unintentional oxygen impurities lowers the superconducting critical temperature of Nb. This evolution coincides with the reduction of Nb crystal domains, which are separated by highly disordered oxygen-rich regions. Moreover, chemical analysis reveals the formation of niobium monoxide within the film during growth. These findings provide a comprehensive picture of how residual oxygen in the growth chamber can affect the properties of the Nb films. This study contributes to the materials science and engineering knowledge of superconducting Nb growth.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2021

Edge localized mode characteristics and divertor heat flux during stationary and transient phase for CFETR hybrid scenario

The study of edge localized mode (ELM) behavior, stationary heat flux and transient heat flux during the ELM crash phase is performed for a China Fusion Engineering Test Reactor (CFETR) 1 GW hybrid mode operation scenario ($R$ = 7.2 m, ${B_{\text{T}}}$ = 6.5 T, ${I_{\text{p}}}$ = 13.78 MA). Modeling and simulation start with a scenario obtained by multi-code integrated modeling on the one modeling framework for integrated tasks framework. Linear stability and nonlinear simulations of ELM dynamics are carried out using the BOUT++ six-field reduced magnetohydrodynamic module, which show a much smaller ELM energy loss (Δ ELM ~ 0.13%) compared to that of a Type-I ELM. Parametric analysis of the weak linear growth rate and small ELM energy loss characteristics shed light on physics corresponding to a grassy ELM regime for CFETR 1 GW hybrid scenario. The transient heat flux on the divertor target during this small ELM phase is investigated using BOUT++. Here we found that upstream radial transport in the scrape-off-layer (SOL) induced by small ELMs is weak, which keeps it in the drift-dominated region. However, the heat flux width is still broadened to ${\lambda _{\text{q}}}$ = 4.64 mm by the increase of separatrix temperature during ELM nonlinear evolution. The impact of transient peak heat load and ELM energy fluence on tungsten melting and net erosion rate of divertor target is evaluated for the first-time using physics-based transport that connects the pedestal with the SOL. Energy fluence caused by a single ELM pulse is below the tungsten melting limit, while tungsten erosion would exceed the material requirements. We conclude that external mitigation methods, such as divertor detachment and advanced divertor geometry are likely needed for the steady state operation of CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modelling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

ABSTRACT Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modelling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron–ion or electron–positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann–Vlasov equation. This new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes, not present in previous formulations of two-fluid magnetohydrodynamics. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. We keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The out-of-equilibrium corrections take the form of a collisional 14-moment closure previously described in the context of viscous single fluids. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Calculation of the C 12 + C 12 sub-barrier fusion cross section in an imaginary-time-dependent mean field theory

The 12 C + 12 C sub-barrier fusion cross section is calculated within the framework of a time-dependent Hartree-Fock-based classical model using the Feynman path-integral method. The modified astrophysical S* factor is compared to direct and indirect experimental results. A good agreement with the direct data is found. In the lower-energy region where recent analyses of experimental data obtained with the Trojan horse method (THM) lead to contrasting results, the model predicts a nonresonant S* factor half-way between those results. Low-energy resonances revealed in the THM data are added to the calculation, and the relative reaction rate in the Gamow region is calculated. In particular, including 0 + resonances result in some agreement with the THM data. Here, the role of different resonances is discussed in detail, and their influence on the reaction rate at temperatures relevant to stellar evolution is investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Global Theory of Microtearing Modes in the Tokamak Pedestal

The pedestal of H-mode tokamaks displays strong magnetic fluctuations correlated with the evolution of the electron temperature. The microtearing mode (MTM), a temperature-gradient-driven instability that alters magnetic topology, is compatible with these observations. Here we extend the conventional theory of the MTM to include the global variation of the temperature and density profiles. The offset between the rational surface and the location of the pressure gradient maximum (μ) emerges as a crucial parameter for MTM stability. In conclusion, the extended theory matches observations on the Joint European Torus tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗
DOE OSTI2022

On the combustion of n-butyl acetate droplets

This paper presents an experimental and numerical study of the combustion of isolated n‑butyl acetate droplets in the standard atmosphere. Numerical simulations are reported using a model that incorporates unsteady gas and liquid transport, variable properties, and radiation. Three skeletal mechanisms of n‑butyl acetate, derived from a large detailed mechanism comprised of 819 species and 52,698 reactions, were used in the numerical simulations to evaluate the influence of the kinetic mechanism on burning. The reduced mechanisms comprised 212 species and 5413 reactions, 157 species and 3089 reactions, and 105 species and 1035 reactions. The numerical model did not include soot formation, though qualitatively mild sooting was noted only for droplets larger than 0.7 mm. The numerical predictions were in good agreement with experimental measurements of droplet and flame diameters. Flame extinction was numerically predicted which was attributed to a decrease of the characteristic diffusion time relative to the chemical time as droplet burned. Effects of initial droplet diameter on the evolution of maximum gas temperature (Tmax) and peak mole fractions of CO2 and CO are also examined numerically.

09 BIOMASS FUELS↗
DOE OSTI2024

METHOD FOR EXPOSING IRRADIATED FUEL CLADDING TO RAPID THERMAL TRANSIENTS

Reactor core materials exposed to rapid thermal transients during accidents can experience significant changes in mechanical properties as irradiation hardness is recovered. Understanding material performance during transients using computer modeling requires accurate descriptions of mechanical property evolution with time and temperature. A simple experimental capability was developed to expose irradiated Zircaloy to a simulated thermal transient using immersion in molten tin (at 675°C) to rapidly heat and hold before rapidly cooling with a water quench. Multiple heating-cooling cycles using a nonirradiated tensile specimen showed simplicity and viability. Based on preliminary test results, a robust system was established for routine evaluation of simulated thermal accident scenarios in the Low Activation Materials Design and Analysis (LAMDA) laboratory at Oak Ridge National Laboratory (ORNL). This paper summarizes system development and nonirradiated Zircaloy commissioning test results. Future experiments will conduct tensile tests on irradiated Zircaloy after various rapid temperature profiles.

Byun, TS↗
DOE OSTI2021

Using Multigrain Crystallography to Explore the Microstructural Evolution of the α-Olivine to γ-Ringwoodite Transformation and ε-Mg2SiO4 at High Pressure and Temperature

The introduction of multigrain crystallography (MGC) applied in a laser-heated diamond anvil cell (LH-DAC) using synchrotron X-rays has provided a new path to investigate the microstructural evolution of materials at extreme conditions, allowing for simultaneous investigations of phase identification, strain state determination, and orientation relations across phase transitions in a single experiment. Here, we applied this method to a sample of San Carlos olivine beginning at ambient conditions and through the α-olivine → γ-ringwoodite phase transition. At ambient temperatures, by measuring the evolution of individual Bragg reflections, olivine shows profuse angular streaking consistent with the onset of yielding at a measured stress of ~1.5 GPa, considerably lower than previously reported, which may have implications for mantle evolution. Furthermore, γ-ringwoodite phase was found to nucleate as micron to sub-micron grains imbedded with small amounts of a secondary phase at 15 GPa and 1000 °C. Using MGC, we were able to extract and refine individual crystallites of the secondary unknown phase where it was found to have a structure consistent with the ε-phase previously described in chondritic meteorites.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Polymorphism and Negative Linear Compressibility in Pyrazine-d4

Pyrazine (1,4-diazine) is a widely used linker in coordination polymers, molecular magnets, and metal–organic frameworks, yet its solid-state polymorphism is controlled by weak C–H···N interactions. Using variable-temperature and high-pressure neutron diffraction, a previously unrecognized low-temperature transition to phase IV is identified, showing that the form stabilized under pressure can also be accessed below ∼90 K at ambient pressure in the perdeuterated compound. The transition from phase III → IV involves a symmetry-breaking lattice distortion driven by the reorganization of C–H···N hydrogen bonds, accompanied by negative linear compressibility, which is primarily accommodated through the reduction in void-space within the structure. Structural analysis reveals that while the overall monoclinic framework of phase IV is preserved under both temperature and pressure routes, the evolution of hydrogen-bonded network is distinct depending on the pathway taken. These results establish pyrazine as a model molecular solid for studying pressure–temperature polymorphism and demonstrate how weak directional interactions can produce anomalous mechanical responses in organic crystals.

Ridley, Christopher J. [Oak Ridge National Laborat↗
DOE OSTI2022

Simulation of dislocation evolution in microparticle impacts over a wide range of impact velocities

We report the dynamic loading environment created by high velocity microparticle impacts is very difficult to model. The extremely high strain rates, large deformations, and high temperatures affect the flow stress of the material in ways that standard flow stress models cannot capture. Plastic deformation and strain hardening in metals is controlled by the motion of dislocations. Dislocations can be nucleated, stored as forest dislocations, or be annihilated as loading progresses. A comprehensive accounting of dislocation density change is needed to accurately describe rate and temperature dependent dislocation glide and evolution across the wide rage of loading conditions present in microparticle impact problems. Therefore, we implement and apply a newly proposed flow stress model (Hunter and Preston, 2015, 2022) that is then coupled with mobile and immobile dislocation density evolution equations. This model is implemented in Los Alamos National Laboratory’s hydrodynamics code, FLAG, to model copper-on-copper microparticle impacts. This new strength model allows for accurate simulation of particle rebound and flattening across a wide range of impact velocities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Operation above the Greenwald density limit in high performance DIII-D negative triangularity discharges

The density limit in strongly-shaped negative triangularity (NT) discharges is studied experimentally in the DIII-D tokamak. Record-high Greenwald fractions f G are obtained, using gas puff injection only, with values up to near 2, where f G is defined as the ratio of the line-averaged density over n G = I p /(π a 2 ), with I p [MA] the plasma current and a[m] the plasma minor radius. A clear higher operational limit with higher auxiliary power is also demonstrated, with the ohmic density limit about two times lower than with additional neutral beam injection heating. The evolution of the electron density, temperature and pressure profiles are analyzed as well. The core density can be up to twice the Greenwald density and keeps increasing, while the value at the separatrix remains essentially constant and slightly below n G . The edge temperature gradient collapses to near zero and NT plasmas are shown to be resilient to such profiles in terms of disruptivity. We also present the time evolution of the inverse electron pressure scale length with the value at the last closed flux surface (LCFS) decreasing below the value at the normalized radius 0.9 near the density limit, demonstrating the clear drop of confinement starting from the edge. This inverse scale length ‘collapse’ at the LCFS also defines well the characteristic behavior of the kinetic profiles approaching a density limit.

density limit↗
DOE OSTI2021

$A$-type antiferromagnetic order and magnetic phase diagram of the trigonal Eu spin-$\frac{7}{2}$ triangular-lattice compound $\mathrm{EuSn_2As_2}$

The trigonal compound EuSn 2 As 2 was recently discovered to host Dirac surface states within the bulk band gap and orders antiferromagnetically below the N´eel temperature T N = 23.5(2) K from neutron-diffraction measurements. Here the magnetic ground state of single-crystal EuSn 2 As 2 and the evolution of its properties versus temperature T and applied magnetic field H are reported. Included are the zero-field single-crystal neutron diffraction measurements versus T, magnetization M(H, T), magnetic susceptibility χ(H, T) = M(H,T)/H, heat capacity C p (H, T), and electrical resistivity ρ(H, T) measurements. The neutron-diffraction and χ(T) measurements both indicate a collinear A-type antiferromagnetic (AFM) structure below T N , where the Eu 2+ spins S = 7/2 in a triangular ab-plane layer (hexagonal unit cell) are aligned ferromagnetically in the ab plane whereas the spins in adjacent Eu planes along the c axis are aligned antiferromagnetically. The χ(H ab , T) and χ(Hc, T) data together indicate a smooth crossover between the collinear AFM alignment and an unknown magnetic structure at H ≈ 0.12 T. Dynamic spin fluctuations up to 60 K are evident in the χ(T), Cp(T) and ρ(H, T) measurements, a temperature that is more than twice T N . The ρ(H, T) is consistent with a low-carrier-density metal with strong magnetic scattering and does not reflect a contribution of the topological state of the material as reported earlier by ARPES measurements. This observation is consistent with previous ones for other topological insulators where the chemical potential is above the Dirac point so that ARPES readily detects the surface states, whereas resistivity measurements do not. Finally, the magnetic phase diagrams for both H ∥ c and H ∥ ab in the H-T plane are constructed from the T N (H), χ(H, T), Cp(H, T), and ρ(H, T) data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effect of fuel composition and EGR on spark-ignited engine combustion with LPG fueling: Experimental and numerical investigation

This paper presents an experimental and numerical investigation of a spark-ignited (SI) cooperative fuel research (CFR) engine fueled with different liquefied petroleum gas (LPG) fuels and exhaust gas recirculation (EGR). Here, the effects of LPG fuel composition on engine combustion characteristics are initially evaluated at two different compression ratios (CR). Results show normal combustion at CR 7 and heavy knocking combustion at CR 10 for all the tested fuels, with a more substantial impact for the LPG fuel with high proportions of n-butane species. The Livengood-Wu (LW) integral method is then used to analyze the knock occurrence risk of individual fuel based on the reactivity of the tested fuels. The introduction of EGR then demonstrates the potential of knock intensity reduction below the borderline knock limit. A zonal-based kinetic interactions study is also performed to understand the knock mitigation effectiveness of EGR over the pressure–temperature domain relevant to SI engine operation. Finally, a multidimensional, computational fluid dynamics (CFD) simulation model is shown to predict the LPG combustion characteristics and presents the evolution of in-cylinder temperature and chemical species to demonstrate the development of end-gas autoignition events without and with EGR.

02 PETROLEUM↗