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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Storage and Diffusion of Carbon and Nitrogen in Crustal Materials

Nitrogen is a primary constituent of Earth’s atmosphere and an essential component of many biological processes. Recent developments in understanding N incorporation into silicate minerals suggest that the solid earth and lithosphere could provide a substantial reservoir for the planetary N budget. Furthermore, determining the concentration and distribution of low abundances of N in silicates presents a significant challenge in microanalysis. We tested capabilities for high spatial resolution quantification of low levels of N using EPMA with the goal of developing standards to be used in coupled N and noble gas analyses of silicate minerals. We investigated the concentration and distribution of N in a suite of ammonium-bearing silicate minerals that reflect the dominant N-bearing phases of the lithosphere: tobelite, buddingtonite, ammonioleucite, tsaregorodtsevite, and hyalophane. These minerals form during NH 4 -rich alteration of K-bearing silicate minerals, and the bonding environment of N in each mineral is likely to differ based on crystal structure. Using an optimized analytical routine for N quantification enables sensitivity approaching 0.02 weight% N and high-resolution quantitative mapping of sub-micron variations in N contents.

58 GEOSCIENCES↗
DOE OSTI2022

Uptake of N 2 O 5 by aqueous aerosol unveiled using chemically accurate many-body potentials

The reactive uptake of N 2 O 5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N 2 O 5 in water. The free energy profile highlights that N 2 O 5 is selectively adsorbed to the liquid–vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction–diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction–diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N 2 O 5 in more complex solutions.

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

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

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

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Band structures and Z 2 invariants of two-dimensional transition metal dichalcogenide monolayers from fully relativistic Dirac-Kohn-Sham theory using Gaussian-type orbitals

Two-dimensional (2D) materials exhibit a wide range of remarkable phenomena, many of which owe their existence to the relativistic spin-orbit coupling (SOC) effects. To understand and predict properties of materials containing heavy elements, such as the transition-metal dichalcogenides (TMDs), relativistic effects must be taken into account in first-principles calculations. We present an all-electron method based on the four-component Dirac Hamiltonian and Gaussian-type orbitals (GTOs) that overcomes complications associated with linear dependencies and ill-conditioned matrices that arise when diffuse functions are included in the basis. Until now, there has been no systematic study of the convergence of GTO basis sets for periodic solids either at the nonrelativistic or the relativistic level. Here, in this study, we provide such a study of relativistic band structures of the 2D TMDs in the hexagonal (2H), tetragonal (1T), and distorted tetragonal (1T') structures, along with a discussion of their SOC-driven properties (Rashba splitting and Z 2 topological invariants). We demonstrate the viability of our approach even when large basis sets with multiple basis functions involving various valence orbitals (denoted triple- and quadruple-ζ) are used in the relativistic regime. Our method does not require the use of pseudopotentials and provides access to all electronic states within the same framework. Our study paves the way for direct studies of material properties, such as the parameters in spin Hamiltonians, that depend heavily on the electron density near atomic nuclei where relativistic and SOC effects are the strongest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew–Burke–Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew–Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. Finally, the quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

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

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TREERAY/RADPRESSURE, is an extension to the novel radiative transfer method TREERAY implemented in the grid-based Magneto-Hydrodynamics code FLASH. In TREERAY/RADPRESSURE, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M ⊙ . We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

(ISM:) HII regions↗
DOE OSTI2020

A Novel NASICON‐Type Na 4 MnCr(PO 4 ) 3 Demonstrating the Energy Density Record of Phosphate Cathodes for Sodium‐Ion Batteries

Abstract Sodium‐ion batteries (SIBs) have attracted incremental attention as a promising candidate for grid‐scale energy‐storage applications. To meet practical requirements, searching for new cathode materials with high energy density is of great importance. Herein, a novel Na superionic conductor (NASICON)‐type Na 4 MnCr(PO 4 ) 3 is developed as a high‐energy cathode for SIBs. The Na 4 MnCr(PO 4 ) 3 nanoparticles homogeneously embedded in a carbon matrix can present an extraordinary reversible capacity of 160.5 mA h g −1 with three‐electron reaction at ≈3.53 V during the Na + extraction/insertion process, realizing an unprecedentedly high energy density of 566.5 Wh kg −1 in the phosphate cathodes for SIBs. It is intriguing to reveal the underlying mechanism of the unique Mn 2+ /Mn 3+ , Mn 3+ /Mn 4+ , and Cr 3+ /Cr 4+ redox couples via X‐ray absorption near‐edge structure spectroscopy. The whole electrochemical reaction undergoes highly reversible single‐phase and biphasic transitions with a moderate volume change of 7.7% through in situ X‐ray diffraction and ex situ high‐energy synchrotron X‐ray diffraction. Combining density functional theory (DFT) calculations with the galvanostatic intermittent titration technique, the superior performance is ascribed to the low ionic‐migration energy barrier and desirable Na‐ion diffusion kinetics. The present work can offer a new insight into the design of multielectron‐reaction cathode materials for SIBs.

Zhang, Jian↗
DOE OSTI2026

Substrate-Directed Underlayer Growth of Bilayer MoS 2 Revealed by Mo Isotope Labeling

Direct control over the vertical formation sequence and stacking registry in van der Waals (vdW) bilayers is essential for device performance and moiré engineering yet difficult to resolve unambiguously with conventional probes. Here, we use Mo isotope labeling in a two-step chemical vapor deposition process to synthesize bilayer MoS 2 and trace its vertical formation on common substrates. By combining site-selective laser thinning, Raman spectroscopy, time-of-flight secondary ion mass spectrometry, and atomic-resolution scanning transimission electron microscopy (STEM), we find a clear substrate dependence: on SiO 2 /Si, the second layer nucleates and grows beneath the first (underlayer), whereas on sapphire, it forms on top (overlayer). Density functional theory indicates that a larger equilibrium interfacial separation and weaker MoS 2 –substrate interactions on amorphous SiO 2 permit confined interfacial diffusion and underlayer nucleation, whereas stronger interactions and smaller separations on sapphire favor overlayer growth. On SiO 2 , confined epitaxy templates commensurate 2H, 3R, and mixed bilayers, as confirmed by second harmonic generation spectroscopy and STEM. During underlayer coalescence, embedded mirror-twin grain boundaries stitch atomically sharp 2H|3R junctions via alternating 4|8 ring motifs. Molecular-dynamics simulations reveal that these alternating 4|8 motifs accommodate interlayer vdW coupling and locally modulate the stacking registry. These results provide mechanistic insight into confined epitaxial growth and establish isotope labeling as a powerful probe of two-dimensional materials synthesis.

MoS2↗
DOE OSTI2024

Simulating energetic ions and enhanced fusion rates from ion-cyclotron resonance heating with a full-wave/Fokker–Planck model

Reproducing fast-ion enhanced fusion rates from ion-cyclotron resonance heating (ICRH) in tokamaks requires the self-consistent coupling of a full-wave solver and a Fokker–Planck solver, which evolves multiple simultaneously resonant ion species. We introduce a new self-consistent model that iterates the TORIC full-wave solver with the CQL3D Fokker–Planck solver using the integrated plasma simulator (IPS). This model evolves the bounce-averaged ion distribution functions in both parallel and perpendicular velocity-space with a quasilinear radio frequency (RF) diffusion operator valid in the ion finite Larmor radius (FLR) limit and the RF electric fields with the resultant non-Maxwellian FLR dielectric tensor. This produces non-Maxwellian ICRH simulations that are fully self-consistent, fast, and interoperable with integrated modeling frameworks, such as TRANSP/GACODE/IPS-FASTRAN. We demonstrate our model's capabilities by validating it against experimental data in Alcator C-Mod. We then perform the first RF heating simulations of SPARC using self-consistent non-Maxwellian ion distributions to investigate the potential to enhance fusion rates using ion cyclotron resonance heating generated fast ions.

Physics↗
DOE OSTI2023

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Exponential time differencing for the tracer equations appearing in primitive equation ocean models

The tracer equations are part of the primitive equations used in ocean modeling and describe the transport of tracers, such as temperature, salinity or chemicals, in the ocean. Depending on the number of tracers considered, several equations may be added to and coupled to the dynamics system. In many relevant situations, the time-step requirements of explicit methods imposed by the transport and mixing in the vertical direction are more restrictive than those for the horizontal, and this may cause the need to use very small time steps if a fully explicit method is employed. To overcome this issue, we propose an exponential time differencing (ETD) solver where the vertical terms (transport and diffusion) are treated with a matrix exponential, whereas the horizontal terms are dealt with in an explicit way. In this work, we investigate numerically the computational speed-ups that can be obtained over other semi-implicit methods, and we analyze the advantages of the method in the case of multiple tracers.

42 ENGINEERING↗
DOE OSTI2022

Role of electronic energy loss on defect production and interface stability: Comparison between ceramic materials and high-entropy alloys

High-entropy alloys (HEAs) and some complex alloys exhibit desirable properties and significant structural stability in harsh environments, including possible applications in advanced reactors. Energetic ion irradiation is often used as a surrogate for neutron irradiation; however, the impact of ion electronic energy deposition and dissipation is often neglected. Moreover, differences in recoil energy spectrum and density of cascade events on damage evolution must also be considered. In many chemically complex alloys, the mean free path of electrons is reduced significantly, thus their decreased thermal conductivity and slow dissipation of localized radiation energy can have noticeable effects on displacement cascade evolution that is greatly different from metals with high thermal conductivity. In this work, nanocrystalline HEAs of Ni 20 Fe 20 Co 20 Cr 20 Cu 20 and nonequiatomic (NiFeCoCr) 97 Cu 3 , both having much lower room-temperature thermal conductivity than pure Ni or Fe, are chosen as model HEAs to reveal the role that electronic energy loss during ion irradiation has in complex alloys. The response of nanocrystalline HEAs is investigated under irradiation at room temperature using MeV Ni and Au ions that have different ratios of electronic energy to damage energy, which is the energy dissipated in displacing atoms. Different from previously reported amorphization of nanocrystalline SiC, experimental results on these HEAs show that, similar to the process in nanocrystalline oxide materials, both inelastic thermal spikes via electron–phonon coupling and elastic thermal spikes via collisions among atomic nuclei contribute to the overall grain growth. The growth follows a power law dependence with the total deposited ion energy, and the derived value of the power-exponent suggests that the irradiation-induced instability at and near grain boundaries leads to local rapid atomic rearrangements and consequently grain growth. The high power-exponent value can be attributed to the sluggish diffusion and delayed defect evolution arising from the chemical complexity intrinsic to HEAs. Here, this work calls attention to quantified fundamental understanding of radiation damage processes beyond that of simplified displacement events, especially in simulating neutron environments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Moisture-driven CO 2 pump for direct air capture

Direct Air Capture, the chemical separation of carbon dioxide (CO 2 ) from air, is considered a necessary approach to generate negative carbon emissions and in turn limit global warming, though energetic and monetary costs must continue to decline to meet the needed gigaton scale. The promise of energy-efficient and continuous carbon dioxide membrane separation has only recently been investigated with the development of facilitated chemical transport and novel electrochemically driven methods. This work demonstrates the potential for continuous pumping of CO 2 from air using commercial anion exchange membranes driven solely by a water vapor gradient. A reactive transport model that couples ionic transport with moisture sensitive bicarbonate chemistry in charged polymers, also termed moisture swing sorption, was developed alongside experimentally measured equilibrium and transport properties of CO 2 and water vapor in Fumasep FAA-3 to determine the mass transport limits of such a system. The results show that the commercial membrane is kinetically limited by the moisture sensitive chemistry, enabling a CO 2 current of 1.1 μmol/m 2 • s. In diffusion limited materials, a CO 2 pumping flux can reach 11 μmol/m 2 • s in and peak performance is expected with a dry air relative humidity around 50% pumping into a saturated water vapor permeate. In all configurations, water loss becomes the major cost in such a separation and must be managed through optimal polymer design or water recovery mechanisms. The work demonstrates the applicability of a new driving force, the evaporation of water, to meet the energetic demand of CO 2 separation from air.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A high-performance hydroxide exchange membrane enabled by Cu 2+ -crosslinked chitosan

We report ion exchange membranes are widely used to selectively transport ions in various electrochemical devices. Hydroxide exchange membranes (HEMs) are promising to couple with lower cost platinum-free electrocatalysts used in alkaline conditions, but are not stable enough in strong alkaline solutions. Herein, we present a Cu 2+ -crosslinked chitosan (chitosan-Cu) material as a stable and high-performance HEM. The Cu 2+ ions are coordinated with the amino and hydroxyl groups of chitosan to crosslink the chitosan chains, forming hexagonal nanochannels (similar to 1nm in diameter) that can accommodate water diffusion and facilitate fast ion transport, with a high hydroxide conductivity of 67 ms cm -1 at room temperature. The Cu 2+ coordination also enhances the mechanical strength of the membrane, reduces its permeability and, most importantly, improves its stability in alkaline solution (only 5% conductivity loss at 80 degrees C after 1,000 h). These advantages make chitosan-Cu an outstanding HEM, which we demonstrate in a direct methanol fuel cell that exhibits a high power density of 305 mWcm -2 . The design principle of the chitosan-Cu HEM, in which ion transport channels are generated in the polymer through metal-crosslinking of polar functional groups, could inspire the synthesis of many ion exchange membranes for ion transport, ion sieving, ion filtration and more.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗