Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase mapping”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Rapid Identification of Synthetic Routes to Functional Metastable Phases Using X-ray Probed Laser Anneal Mapping (XPLAM) Time-Temperature Quench Maps

Many material systems have known or predicted functional phases that are metastable at standard temperature and pressure. While substantial advances have been made in the high-throughput and combinatorial synthesis of materials with a range of stoichiometries, investigation of thermal processing remains largely the domain of iterative uniform anneals or static gradients. Here we develop X-ray probed laser anneal mapping (XPLAM), a high throughput technique coupling spatially resolved X-ray diffraction with microsecond to millisecond laser gradient anneals to produce temperature–dwell–transformation (TDT) diagrams of the phase as a function of quench time and temperature. In addition to showing regimes where specific metastable phases form preferentially, TDT diagrams provide insight into the submillisecond kinetics of solid–solid phase transitions. This is a unique tool for mapping reaction pathways for metastable phases. As a first demonstration of XPLAM, we study Bi2O3, which has a rich set of polytypes, including the d-phase with an exceptionally high oxygen ion conductivity. We demonstrate the first annealing-driven synthesis of room temperature d-Bi2O3. We expect XPLAM to prove a powerful technique for rapid identification of synthetic routes to metastable phases and to generate the exhaustive data sets required for machine learning-guided exploration of materials processing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structural Manipulation of Phase Transitions by Self‐Induced Strain in Geometrically Confined Thin Films

Strain engineering is a well-known method often used to tune material properties in thin films. The most studied sources of strain are lattice mismatch and differential thermal contraction between the substrate and film. However, in materials which undergo a structural phase transition (SPT), a third and often overlooked source of strain may play a very significant role. If the substrate confines the area of the film, the SPT may induce stress which changes the evolution of the transition. This is a 2D analog of the isochoric phase transition between water and ice, where the freezing point drops below 0 °C. To illustrate this, the prototypical Mott insulator V 2 O 3 which has an SPT coupled to a metal–insulator transition is used to show how self-induced strain can drastically alter structural and electronic properties. This effect provides an elegant approach for mapping the phase diagram of the SPT and the transitions coupled to it. Moreover, the magnitude of self-straining is tunable by modifying the substrate morphology. Furthermore, this effect may be important for numerous materials which exhibit an SPT and are subjected to geometrical constraints.

Kalcheim, Yoav↗
DOE OSTI2026

Thin-Film Stabilization and Magnetism of η -Carbide-Type Iron Nitrides

Transition-metal nitrides in ..eta..-carbide-type structures exhibit unusual bonding motifs and proximity to magnetic instabilities. Yet they remain unexplored in thin-film form due to the difficulty of stabilizing nitrogen-poor ternaries among competing phases. Here, we report the thin-film synthesis and phase-stability mapping of the ..eta..-nitride systems Fe-W-N and Fe-Mo-N. Amorphous Fe-M-N (M = W, Mo) combinatorial libraries deposited by reactive cosputtering crystallize upon rapid thermal annealing, enabling systematic identification of synthesis windows as a function of composition and annealing temperature. Using laboratory powder X-ray diffraction and synchrotron grazing incidence wide-angle X-ray scattering, we establish that Fe 3 Mo 3 N-based ..eta..-carbide phases form over a substantially broader compositional and thermal range than W-based compositions, where ..eta.. structures are stabilized only when the films are Fe-rich. These trends are rationalized using mixed chemical-potential vs composition phase diagrams that capture the narrow nitrogen chemical-potential stability of ..eta..-nitrides. Magnetic measurements reveal that ferromagnetism is induced in Fe-rich Fe 3.54 Mo 2.46 N with a small exchange-bias-like response that is absent in Fe 3 W 3 N-based compositions, highlighting the sensitivity of magnetic behavior to modest deviations from stoichiometry. This work establishes practical thin-film synthesis routes for ..eta..-nitride materials and demonstrates how composition can be tuned to access emergent magnetic phenomena in these complex nitrides.

36 MATERIALS SCIENCE↗
DOE OSTI2026

An approach for fast and accurate simulation of phase change material based thermal energy storage in buildings

Latent heat thermal energy storage (LHTES) has significant potential for mitigating peak electricity demand and enabling load shifting in buildings. Phase Change Material embedded heat exchangers (PCM-HX) can significantly improve energy demand management due to high storage capacity. However, PCM-HX evaluation typically depends on computationally expensive fully transient simulations, posing significant challenges for scalable system- and building-level energy assessments across different climates and system architectures. This paper presents a generalized, accurate, and computationally efficient methodology for simulating building energy systems integrated with LHTES. The PCM-HX transient performance is represented by performance maps generated using a Generalized Resistance-Capacitance Model (GRCM) that enables accurate predictions of arbitrary PCM-HXs at low computational cost. The feasibility of the proposed approach was verified using a case study considering a dual-mode heat pump-thermal energy storage (HP-TES) system simulated in Modelica with Spawn of EnergyPlus™ for a DOE prototype small office building in two locations: Tampa, FL, and International Falls, MN. The PCM-HX performance maps provided accurate predictions of PCM-HX transient behavior, with mean absolute percentage deviations within 2–4% compared to GRCM while also achieving at least 1800× reduction in computational time. Moreover, the HP-TES system achieved energy savings of up to 17.4% in Tampa, FL, and 62.2% in International Falls, MN, demonstrating the broader applicability of the proposed methodology across different climate zones. This work highlights the importance of robust PCM-HX models in enabling accurate and computationally efficient building-level simulations and enabling future research opportunities for investigating optimized HP-TES designs and advanced control strategies for grid-interactive buildings.

Modelica↗
DOE OSTI2025

Imaging Bragg Edge Analysis TooLs for Engineering Structures (iBeatles)

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) provides pulsed neutrons with energies varying from epithermal to cold. In preparation for VENUS, the neutron imaging beamline to be located at beam port 10, we have performed a series of experiments focused on wavelength-dependent radiography and computed tomography for a broad range of applications, from materials science to biological tissues.One of the time-of-flight (TOF) techniques that is of interest to the scientific community is the 2-dimensional mapping of phases and average crystalline plane orientation in samples both ex-situ and during applied stresses such as tensile loading and heating. This technique is known as Bragg edgeimaging and relies on the identification of changes of transmission values, fitting of the edge to measure its displacement, and thus identify the shift in lattice parameter due to stresses. One of the challenges of TOF imaging measurements is the amount of data and the inability to observe Bragg edge shifts in real time during an experiment. Thus, we have been focusing on creating a Python-based interface that allows fast data processing and instantaneous mapping and fitting of the Bragg edges, and their evolution through time. Python libraries and Jupyter notebooks have been implemented to facilitate decision making during an experiment. The advantage of the notebooks is the possibility to guide an experiment as they can quickly process and display Bragg edge data. These notebooks can be used independently, or can be combined in a Python Graphical User Interface (GUI) tool called iBeatles. This interface permits visualization and fitting of the Bragg edges, and ultimately back-projects the fitting results onto the radiographs to display a strain map. Assuming data collection has sufficient statistics, the strain mapping analysis can be performed on a pixel-by-pixel basis. This development is a step forward toward a better user experience at the future VENUS beamline in terms of live feedback and productivity. Analysis that used to take days of switching between different applications can now be done in minutes within the

Bilheux, JeanChristophe [Oak Ridge National Labora↗
DOE OSTI2025

Transitory Topochemical Tailoring of a van der Waals Superconductor

Topochemical intercalation is widely used to access metastable phases with novel electronic properties, but the reverse reaction (deintercalation) typically restores the original state, limiting practical use. Here, in this study, we present a topochemical approach that employs a sacrificial intercalant that thermally decomposes to irreversibly lock in the new electronic state. Using 2-aminobutane as the sacrificial intercalant, we convert the van der Waals (vdW) material 1T-TiSe 2 into a superconductor and the vdW superconductor 2H-NbSe 2 into a nonsuperconducting metal, while preserving the ability to exfoliate the resulting crystals. We find that this transitory intercalation increases the electron density in both materials and partially suppresses the CDW in TiSe 2 . By tuning the thermolysis temperature, we can systematically vary the carrier density in TiSe 2 , enabling us to map its phase diagram. The superconductivity in TiSe 2 is retained in exfoliated flakes, although with a lower critical temperature. This transitory topochemical strategy enables access to new electronic states with precisely tuned carrier densities that are otherwise inaccessible through direct solid-state synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗
DOE OSTI2022

Ambient Noise Attenuation Tomography Reveals an Asymmetric Damage Zone Across San Jacinto Fault Near Anza, California

Here, we perform seismic attenuation tomography of the shallow structure for the San Jacinto Fault in the Ramona Reservation of southern California. The study uses ambient seismic noise recorded by a linear array of 65 3-C sensors across the fault. We extract amplitude decay information, with uncertainty, from noise interferometry functions. To account for strong heterogeneities in the complex shallow fault zone structure, we apply a frequency-dependent amplitude correction for focusing/defocusing effects using inverted phase velocity maps obtained by solving the transport equation. We then estimate the attenuation structure based on the linear station-triplet method for both Love and Rayleigh waves. The attenuation tomography indicates strong attenuation correlated with known San Jacinto Fault surface traces. The Love wave attenuation tomography reveals an asymmetric damage zone that exists primarily on the fault side with faster seismic velocity, consistent with earthquake ruptures on a fault bimaterial interface with preferred propagation direction to the northwest.

58 GEOSCIENCES↗
DOE OSTI2021

Defect chemistry and doping of BiCuSeO

While p-type BiCuSeO is a well-known mid-temperature oxide thermoelectric (TE) material, computations predict that superior TE performance can be realized through n-type doping. In this study, we use first-principles defect calculations to show that Cu vacancies are responsible for the native p-type self doping; yet, we find that BiCuSeO is n-type dopable under Cu-rich growth conditions, where the formation of Cu vacancies is suppressed. We computationally survey a broad suite of 23 dopants and find that only Cl and Br are effective n-type dopants. Therefore, we recommend that future experimental doping efforts utilize phase boundary mapping to optimize the electron concentration and resolve the anomalous p–n–p transitions observed in halogen-doped BiCuSeO. The prospect of n-type doping, as revealed by our defect calculations, paves the path for rational design of BiCuSeO chemical analogues with similar doping behavior and even better TE performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Alfvénic modes excited by the kink instability in PHASMA

Magnetic flux ropes have been successfully created with plasma guns in the newly commissioned PHAse Space MApping (PHASMA) experiment. The flux ropes exhibit the expected m = 1 kink instability. The observed threshold current for the onset of this kink instability is half of the Kruskal–Shafranov current limit, consistent with predictions for the non-line tied boundary condition of PHASMA. The helicity, paramagnetism, and growth rate of the observed magnetic fluctuations are also consistent with kink instability predictions. The observed fluctuation frequency appears to be a superposition of a real frequency due to a Doppler shift of the kink mode arising from plasma flow (~2kHz) and a contribution from a wave mode (~5 kHz). The dispersion of the wave mode is consistent with an Alfvén wave. Distinct from most previ-ous laboratory studies of flux ropes, the working gas in PHASMA is argon. Thus, the ion cyclotron frequency in PHASMA is quite low and the frequency of the Alfvénic mode plateaus at ~0.5 of the ion gyro frequency with increasing background magnetic field strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron-only reconnection and associated electron heating and acceleration in PHASMA

In this work, using incoherent Thomson scattering, electron heating and acceleration at the electron velocity distribution function (EVDF) level are investigated during electron-only reconnection in the PHAse Space MApping (PHASMA) facility. Reconnection arises during the merger of two kink-free flux ropes. Both push and pull type reconnection occur in a single discharge. Electron heating is localized around the separatrix, and the electron temperature increases continuously along the separatrix with distance from the X-line. The local measured gain in enthalpy flux is up to 70% of the incoming Poynting flux. Notably, non-Maxwellian EVDFs comprised of a warm bulk population and a cold beam are directly measured during the electron-only reconnection. The electron beam velocity is comparable to, and scales with, electron Alfvén speed, revealing the signature of electron acceleration caused by electron-only reconnection. The observation of oppositely directed electron beams on either side of the X-point provides “smoking-gun” evidence of the occurrence of electron-only reconnection in PHASMA. 2D particle-in-cell simulations agree well with the laboratory measurements. The measured conversion of Poynting flux into electron enthalpy is consistent with recent observations of electron-only reconnection in the magnetosheath [Phan et al., Nature 557, 202 (2018)] at similar dimensionless parameters as in the experiments. The laboratory measurements go beyond the magnetosheath observations by directly resolving the electron temperature gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dissociative electron attachment studies of gas-phase acetic acid using a velocity map imaging technique

Abstract Advancing instrumentation to explore dissociative electron attachment (DEA) studies allows previously unattainable information to be acquired. Using a newly constructed velocity map imaging spectrometer, we revisited a study on DEA to gas-phase acetic acid. We discuss possible fragmentation channels and compared the corresponding ion yields with previous high electron-energy resolution results. We focus on the channels occurring at higher energies, particularly near 10 eV, and calculate their thermodynamic thresholds. Moreover, we expand previous studies and perform time-sliced imaging near the 10 eV resonance to obtain the kinetic energy distribution of the fragment ions.

Optics↗
DOE OSTI2024

Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER

The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths λ q , input powers, and particle fueling locations. Under H-mode scenario conditions with an upstream separatrix density of 1 x 10 20 m -3 , the most conservative range of λ q extrapolations ( 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power (P SOL = 29 MW) conditions, and 2/3 field with P SOL = 10 MW. Increasing the cross-field SOL diffusivities by 2–10× reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. The simulations include parallel currents but neglect cross-field drifts. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fueling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced with increasing diffusivity. The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Response of the mode Grüneisen parameters with anisotropic compression: A pressure and temperature dependent Raman study of β-Sn

The lattice dynamic response of body-centered tetragonal beta-Sn (I4 1 /amd) under high-pressure and temperature conditions is determined using experimental optical vibration modes. Here, Raman scattering is used to map the phase stability region of beta-Sn to perform mode Gruneisen analysis, and we demonstrate the necessity of an optical intensity calibration for Raman thermometry. The Gruneisen tensor is evaluated along a set of isotherms to address shortcomings of single-mode Gruneisen parameters with respect to anisotropic deformations of this tetragonal structured soft metal. The changes observed here in the Gruneisen tensor as a function of temperature are related to anharmonicity and denote potential criteria for the onset of premelting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Dynamics and nonmonotonic drag for individually driven skyrmions

Here, we examine the motion of an individual skyrmion driven through an assembly of other skyrmions in the absence of quenched disorder. The skyrmion behavior is determined by the ratio of the damping and Magnus terms, known as the intrinsic skyrmion Hall angle. For fixed drive in the damping dominated regime, the effective viscosity decreases monotonically with increasing skyrmion density. In contrast, in the Magnus dominated regime, the velocity varies nonmonotonically with density, and there is a regime in which the skyrmion moves faster with increasing density, as well as a pronounced speed-up effect in which a skyrmion traveling through a dense medium moves more rapidly than for low densities or the single-particle limit. The velocity-force curves in the Magnus-dominated regime exhibit an overshoot effect and other differences from the damping-dominated regime. We find a finite threshold force for skyrmion motion which increases with density as well as a drive-dependent skyrmion Hall angle. We map dynamic phase diagrams showing the threshold for motion, nonlinear flow, speed-up, and saturation regimes. In some cases, increasing the density can reduce the skyrmion Hall angle while producing a velocity boost, which could be valuable for applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Lithium's low-temperature phase transitions: Insights into quantum lattice dynamics and superconductivity

The large lattice dynamics of lithium, driven by its low atomic mass, results in energetically similar structures and significant isotope effects under pressure, posing challenges to current theoretical models. Above 20 GPa and at low temperatures, lithium's electronic properties deviate from simple metallic behavior, with superconductivity emerging in a complex, pressure-dependent manner, alongside an unusual isotope effect. The structural phases of 7 Li reported under these conditions are inconsistent across studies, and the structures of 6 Li remain unexamined. These gaps limit our understanding of the effects of pressure on lithium's electronic properties and the role of quantum lattice effects on its structural behavior under pressure. Here, we integrate experimental and theoretical approaches to investigate the low-temperature structural phase boundaries in lithium isotopes. We map the structural phase diagram of 7 Li from 5 to 55 GPa and 15–75 K, identifying the sequence fcc → ℎ⁡R⁢1 → cI16. A pronounced isotope effect is observed, with 6 Li shifting the fcc → ℎ⁡R⁢1 phase boundary to lower pressures at 15 K. Density functional theory calculations further clarify how these structural changes affect superconducting properties, particularly emphasizing the role of the fcc → ℎ⁡R⁢1 transition in lithium's superconductivity. Furthermore, our findings offer insights into the unique behavior of lithium isotopes under pressure.

36 MATERIALS SCIENCE↗