Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Probing the electronic structure and spectroscopy of pyrrolyl and imidazolyl radicals using high-resolution photoelectron imaging of cryogenically cooled anions

High-resolution photoelectron imaging and photodetachment spectroscopy of cryogenically cooled pyrrolide and imidazolide anions are used to probe the electronic structure and spectroscopy of pyrrolyl and imidazolyl radicals. The high-resolution data allow the ground state vibronic structures of the two radicals to be completely resolved, yielding accurate electron affinities of 2.1433 ± 0.0008 eV and 2.6046 ± 0.0006 eV for pyrrolyl and imidazolyl radicals, respectively. Fundamental frequencies for eight vibrational modes of pyrrolyl and ten vibrational modes of imidazolyl are measured, including several nonsymmetric Franck–Condon-forbidden modes. Two electronic excited states are also observed for the two radicals, displaying diffuse spectral features in both systems. The observations of nonsymmetric vibrational modes in the ground states and the diffuse excited state features provide strong evidence for vibronic couplings between the ground state and the two close-by excited states. The 2-pyrrolide isomer is also observed as a minor species from the electrospray ionization source and the electron affinity of 2-pyrrolyl is measured to be 1.6690 ± 0.0030 eV along with five vibrational frequencies. Even though the HOMOs of both pyrrolide and imidazolide anions are p orbitals, photodetachment spectroscopy reveals completely different threshold behaviors for the two anions: a d-wave-dominated spectrum for pyrrolide and an s-wave-dominated spectrum for imidazolide. Furthermore, the current study provides a wealth of electronic and spectroscopic information, which is ideal to compare with more accurate vibronic coupling calculations for these two important radicals, as well as interesting information about the photodetachment dynamics of the two anions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Time-dependent modeling of coupled plasma-wall dynamics

Effects of dynamic wall outgassing on divertor plasma during an ELM-like heat pulse are modelled using newly coupled edge plasma code UEDGE and all reaction-diffusion code FACE. Different divertor regimes are simulated depending on the degree of the plasma detachment. It is shown that hydrogen outgassing from the divertor target plate induced by ELM pulse in the semi-detached regime can trigger deep plasma detachment. The induced detachment does not spontaneously reverse after the pulse due to ineffective implantation of low energy plasma hydrogen in the target material. Furthermore, we demonstrate that plasma neutrals play essential role in heat transport in a detached divertor. The modelled hydrogen dynamics in the wall material demonstrates that the heat pulse leads to enhanced hydrogen transport toward the wall bulk.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗
DOE OSTI2024

Optical pulse-induced ultrafast antiferrodistortive transition in SrTiO 3

The ultrafast dynamics of the antiferrodistortive phase transition in perovskite SrTiO 3 is monitored via time-domain Brillouin scattering. Using femtosecond optical pulses, we initiate a thermally driven tetragonal-to-cubic structural transformation and detect the crystal phase through changes in the frequency of Brillouin oscillations (BO) induced by propagating acoustic phonons. Coupling the measured BO frequency with a spatiotemporal heat diffusion model, we demonstrate that, for a sample kept in the tetragonal phase, deposition of sufficient thermal energy induces a rapid transformation of the heat-affected region to the cubic phase. The initial phase change is followed by a slower reverse cubic-to-tetragonal phase transformation occurring on a timescale of hundreds of picoseconds. Here, we attribute this ultrafast phase transformation in the perovskite to a structural resemblance between atomic displacements of the R-point soft optic mode of the cubic phase and the tetragonal phase, both characterized by anti-phase rotation of oxygen octahedra. The structural relaxation time exhibits a strong temperature dependence consistent with the prediction of the equation of motion describing collective oxygen octahedra rotation based on the energy landscape of the phenomenological Landau theory of phase transitions. Evidence of such a fast structural transition in perovskites can open up new avenues in information processing and energy storage sectors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗
DOE OSTI2026

Multi-physics Preconditioning for Thermally Activated Batteries

Thermal batteries, also known as molten-salt batteries, are single-use reserve power systems activated by pyrotechnic heat generation, which transitions the solid electrolyte into a molten state. The simulation of these batteries relies on multiphysics modeling to evaluate performance and behavior under various conditions. This paper presents advancements in scalable preconditioning strategies for the Thermally Activated Battery Simulator (TABS) tool, enabling efficient solutions to the coupled electrochemical systems that dominate computational costs in thermal battery simulations. We propose a hierarchical block Gauss-Seidel preconditioner implemented through the Teko package in Trilinos, which effectively addresses the challenges posed by tightly coupled physics, including charge transport, porous flow, and species diffusion. The preconditioner leverages scalable subblock solvers, including smoothed aggregation algebraic multigrid (SA-AMG) methods and domain-decomposition techniques, to achieve robust convergence and parallel scalability. Strong and weak scaling studies demonstrate the solver’s ability to handle problem sizes up to 51.3 million degrees of freedom on 2048 processors, achieving near sub-second setup and solve times for the end-to-end electrochemical solve. These advancements significantly improve the computational efficiency and turnaround time of thermal battery simulations, paving the way for higher-resolution models and enabling the transition from 2D axisymmetric to full 3D simulations.

25 ENERGY STORAGE↗
DOE OSTI2021

Effective field theory of stochastic diffusion from gravity

Planar black holes in AdS have long-lived quasinormal modes which capture the physics of charge and momentum diffusion in the dual field theory. How should we characterize the effective dynamics of a probe system coupled to the conserved currents of the dual field theory? Specifically, how would such a probe record the long-lived memory of the black hole and its Hawking fluctuations? We address this question by exhibiting a universal gauge invariant framework which captures the physics of stochastic diffusion in holography: a designer scalar with a gravitational coupling governed by a single parameter, the Markovianity index. We argue that the physics of gauge and gravitational perturbations of a planar Schwarzschild-AdS black hole can be efficiently captured by such designer scalars. We demonstrate that this framework allows one to decouple, at the quadratic order, the long-lived quasinormal and Hawking modes from the short-lived ones. It furthermore provides a template for analyzing fluctuating open quantum field theories with memory. In particular, we use this set-up to analyze the diffusive Hawking photons and gravitons about a planar Schwarzschild-AdS black hole and derive the quadratic effective action that governs fluctuating hydrodynamics of the dual CFT. Along the way we also derive results relevant for probes of hyperscaling violating backgrounds at finite temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Efficient Reformulation of Linear and Nonlinear Solid-Phase Diffusion in Lithium-ion Battery Models using Symmetric Polynomials: Mass Conservation and Computational Efficiency

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms, along with suitable discretization or approximations for the solid-phase diffusion equation. The solid-phase diffusion equation represents the main computational burden for typical pseudo-2-dimensional (p2D) models since these equations in the pseudo r -dimension must be solved at each point in the computational grid. This substantially increases the complexity of the model as well as the computational time. Traditional approaches towards simplifying solid-phase diffusion possess certain significant limitations, especially in modeling emerging electrode materials which involve phase changes and variable diffusivities. A computationally efficient representation for solid-phase diffusion is discussed in this paper based on symmetric polynomials using Orthogonal Collocation and Galerkin formulation (weak form). A systematic approach is provided to increase the accuracy of the approximation (p form in finite element methods) to enable efficient simulation with a minimal number of semi-discretized equations, ensuring mass conservation even for non-linear diffusion problems involving variable diffusivities. These methods are then demonstrated by incorporation into the full p2D model, illustrating their advantages in simulating high C-rates and short-time dynamic operation of Lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Poisson_CCD: A dedicated simulator for modeling CCDs

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Representing methane emissions from wet tropical forest soils using microbial functional groups constrained by soil diffusivity

Tropical ecosystems contribute significantly to global emissions of methane (CH 4 ), and landscape topography influences the rate of CH 4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH 4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH 4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH 4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O 2 ) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O 2 concentrations increased in the valley. We simulated the dynamics of CH 4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH 4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O 2 , acetate, and associated changes in soil pH with topography regulated simulated CH 4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH 4 production (acetate, CO 2 , and H 2 ) and oxidation (O 2 and CH 4 ) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH 4 . A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH 4 emissions. This study enhanced the predictive capability for CH 4 emissions associated with complex topography and drought in wet tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

CEFR simulation using diffusion code system RAST-F

In this study, the RAST-F code system, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver, was used for the neutronic simulation of the CEFR start-up experiments. The numerical solution of the RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. The RAST-F solution of the selected experimental simulations was compared against the measurement data. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of k{sub eff} and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental data within two-sigma of experimental uncertainty. The good agreement of these results indicates that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Tuning the Magnetic Properties of CrI 3 Using Ni Thin Film Deposition for Applications in Spintronic Devices

Chromium triiodide (CrI 3 ), a van der Waals magnet, has recently been shown to host Ising ferromagnetism down to the monolayer limit. It is a potentially important material in 2D magnet-based applications, such as magnetic sensors and spintronic devices. Prior studies have revealed the coexistence of two different types of interlayer magnetic coupling, the antiferromagnetic (AFM) coupling near the surface and ferromagnetic (FM) coupling in the deep bulk layers, in pristine CrI 3 crystals below the Curie temperature. In this study, we used Ni thin film deposition to tune the surface magnetic states in bulk CrI 3 . A nanometer thickness (4 nm thick) of Ni was deposited on the surface of CrI 3 using electron-beam evaporation to form a Ni/CrI 3 heterostructure. The deposited Ni thin layer forms nanoclusters that completely cover the CrI 3 surface. Magnetic states of CrI 3 before and after Ni deposition are probed by ultralow-frequency magneto-Raman spectroscopy. Instead of seeing three spin wave branches in Raman scattering below 2 T, as in pristine bulk CrI 3 , in Ni/CrI 3 , we observe only a single spin wave branch, which is softened compared to that in pristine CrI 3 and displays a Zeeman shift under an out-of-plane magnetic field up to 7 T. This observation reveals that the AFM layers on pristine CrI 3 crystal surfaces transit into FM layers, so the entire CrI 3 crystal is in the FM state after Ni deposition. First-principles calculations show that Ni atoms tend to diffuse into the CrI 3 lattice, and the FM interlayer coupling has a much lower energy than AFM coupling in the presence of Ni atoms. As a result, our studies show that the magnetic state of CrI 3 can be modified through deposition of a thin metal layer on the surface, offering a route for controlling the spin degree of freedom in van der Waals magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane

Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Hydrogen charging and desorption from microstructural viewpoint: A method for deconvoluting hydrogen desorption spectra and unveiling the hydrogen-microstructure interaction

Understanding the interaction of hydrogen with microstructural features in metallic materials is crucial for designing hydrogen-resistant alloys. Although thermal desorption spectroscopy (TDS) is widely used for investigating the hydrogen binding behavior of various microstructural features, its application to face-centered cubic (fcc) metals and alloys that exhibit low hydrogen diffusivity is limited due to the lumped TDS desorption signals. This paper shows that, by coupling a Sofronis–McMeeking type hydrogen transport model with a microstructure-informed finite-element model, TDS data can be deconvoluted to reveal the underlying adsorption–diffusion–desorption processes, hydrogen diffusivity, and trap-binding energies. In conclusion, the austenitic steel SS316L in solution-annealed condition is used as a demonstration material, and we focused on investigating the interaction of deuterium (hydrogen isotope) with grain boundaries, which is difficult to investigate from experiments alone but critical for design of alloys for hydrogen infrastructure.

Finite element simulation, Polycrystalline Microst↗
DOE OSTI2021

Localization crossover and subdiffusive transport in a classical facilitated network model of a disordered interacting quantum spin chain

Here, we consider the random-field Heisenberg model, a paradigmatic model for many-body localization (MBL) and add a Markovian dephasing bath coupled to the Anderson orbitals of the model's noninteracting limit. We map this system to a classical facilitated hopping model that is computationally tractable for large system sizes and investigate its dynamics. The classical model exhibits a robust crossover between an ergodic (thermal) phase and a frozen (localized) phase. The frozen phase is destabilized by thermal subregions (bubbles), which thermalize surrounding sites by providing a fluctuating interaction energy and so enable off-resonance particle transport. Investigating steady-state transport, we observe that the interplay between thermal and frozen bubbles leads to a clear transition between diffusive and subdiffusive regimes. This phenomenology both describes the MBL system coupled to a bath and provides a classical analog for the MBL transition in the corresponding quantum model, in that the classical model displays long local memory times. It also highlights the importance of the details of the bath coupling in studies of MBL systems coupled to thermal environments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Probing a New Regime of Neutrino Self-Interactions with Astrophysical Neutrinos and the Relativistic Cosmic Neutrino Background

Neutrino self-interactions beyond the standard model have profound implications in astrophysics and cosmology. In this Letter, we study an uncharted scenario in which one of the three neutrino species has a mass smaller than the temperature of the cosmic neutrino background. This results in a relativistic component that significantly broadens the absorption feature on the astrophysical neutrino spectra, in contrast to the sharply peaked absorption expected in the extensively studied scenarios assuming a fully nonrelativistic cosmic neutrino background. By solving the Boltzmann equations for neutrino absorption and regeneration, we demonstrate that this mechanism provides novel sensitivity to sub-keV mediator masses, well below the traditional ∼1–100 MeV range. Future observations of the diffuse supernova neutrino background with Hyper-Kamiokande could probe coupling strengths down to 𝑔 ∼10 −8 , surpassing existing constraints by orders of magnitude. These findings open new directions for discoveries and offer crucial insights into the interplay between neutrinos and the dark sector.

79 ASTRONOMY AND ASTROPHYSICS↗