SEARCH · Engineering Papers
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.
Irradiation induced non-equilibrium point defect populations influence mass transport in oxides, which in turn affects their stability and performance in hostile environments. In this study a strong dose rate dependence is observed.
Ion exchange in Rb-doped KTiOPO 4 has facilitated significant advancements in ferroelectric domain engineering, yet understanding the underlying mechanisms remains in its infancy. We perform time-of-flight secondary ion mass spectrometry analysis on multiple periodically ion-exchanged and periodically poled Rb-doped KTiOPO 4 samples under different temperatures and annealing durations. The results are compared between annealing in air, which involved ex situ annealing before periodic poling, and vacuum annealing conducted in situ after periodic poling. The Rb + diffusion profile after periodic ion exchange forms a tooth-shaped pattern. We show that in situ annealing causes a surface pinning effect on the nonpolar face, limiting Rb + migration along the polar axis at the surface. Once the pinned layer is removed through milling, the underlying Rb + diffusion is distinctively different from the surface. Additionally, the rate of Rb + diffusion during in situ annealing is linear, while the periodic domain structures remain stable after annealing. These results contribute to understanding the ionic diffusion process in a ferroelectric ionic conductor and using ion exchange to tailor the linear and nonlinear properties of KTiOPO 4 .
Triple-conducting materials have been proved to improve the performance of popular protonic ceramic electrolysis cells. However, partially because of the complexity of the water splitting reaction involving three charge carriers, that is, oxygen (O 2– ), proton (H + ), and electron (e – ), the triple-conducting reaction mechanism was not clear, and the reaction conducting pathways have seldom been addressed. In this study, the triple conducting Ruddlesden–Popper phase Pr 1.75 Ba 0.25 NiO 4+δ as an anode on the BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3–δ electrolyte was fabricated and its electroresponses were characterized by electrochemical impedance spectroscopy with various atmospheres and temperatures. The impedance spectra are deconvoluted by means of the distribution of the relaxation time method. The surface exchange rate and chemical diffusivity of H + and O 2– are characterized by electrical conductivity relaxation. The physical locations of electrochemical processes are also identified by atomic layer deposition with a surface inhibitor. A microkinetics model is proposed toward conductivities, triple-conducting pathways, reactant dependency, surface exchange and bulk diffusion capabilities, and other relevant properties. Lastly, the rate-limiting steps and suggestions for further improvement of electrode performance are presented.
Reactant transport is an important consideration in the design of ideal electrode structures. For the oxygen reduction reaction catalyzed by Pt/C in proton exchange membrane fuel cell cathodes, O 2 and H + delivery to Pt surfaces and H 2 O transport away are required. Some Pt nanoparticles may only be accessible via micropores that are too small for ionomer molecules to enter, possibly requiring flooding for H + transport. Here, to test if these “buried” Pt particles can play a role in activity through this proposed transport mechanism, we have performed atomic-scale simulations based on reactive force field molecular dynamics with an emphasis on confinement below 20 Å. Diffusion coefficients as a function of the molar concentration and local environment have been evaluated in water domains confined in two-dimensional graphene nanochannels of various channel heights representing a morphological model for micropores in proton exchange fuel cell cathodes. Our study shows that local atomic-scale structures can strongly modify H + , O 2 , and H 2 O transport rates in flooded micropores less than 20 Å in size. We find that there is a critical crossover in diffusion behavior around the 20 Å spacing with larger pores having bulk-like diffusion properties and confinement below 20 Å monotonically decreases diffusion rates. As pore size decreases, we observe locally dispersed water regions that ultimately strand reactants from long-distance transport. These findings suggest that flooded micropores may in fact be viable transport pathways for relevant reactants and products if the pore walls on opposite sides remain separated by ≥10 Å separation. Furthermore, the confinement effect is so strong that N-doping and C-vacancy defects in the C pore wall have only a minimal impact on diffusion rates and their effects are counterintuitively more apparent at larger spacings. These findings provide valuable insight regarding cathode performance and the role “stranded” catalyst particles may play in fuel cell cathodes.
Alloys and composites of aluminum (Al) have shown promise in regulating and enhancing particle combustion for energy generation applications. Recent work on aluminum-silicon (Al-Si) spherical alloy particles has shown improved combustion at low heating rates through enhanced diffusion accompanying a lower melting temperature. This study extends reactivity analysis to higher heating rates, comparing oxidation of Al-Si with Al. Flame speeds of Al-Si powder mixed with molybdenum trioxide (MoO 3 ) powder (Al-Si+MoO 3 ) exhibited a faster transition to steady propagation relative to Al+MoO 3 . Bomb calorimetry experiments revealed up to 5.8% greater early temperature rise for Al-Si powder. Rapid steady propagation for Al-Si mixture and faster temperature rise for Al-Si particles were attributed to accelerated kinetics evidenced in thermal equilibrium analysis of the mixtures using a differential scanning calorimeter (DSC). Larger DSC exotherms in the early stages of oxidation (i.e., 480–720 °C) correlated with early heat release and promoted steady flame propagation for Al-Si+MoO 3 compared with Al+MoO 3 . Furthermore, in multiple heating rate DSC studies, Al-Si+MoO 3 ignited at a lower heating rate (i.e., 15 ⁰C/min) than Al+MoO 3 (i.e., 20 ⁰C/min). Both mixtures ignited at temperatures less than 625 °C which is above the melting temperature of Al-Si (574 °C) and below Al (659 °C). Thus, alloy particle fuels with enhanced diffusion-controlled kinetics promote steady flame propagation and show promise for energy generation applications. This is especially promising for technologies driven by enhanced condensed phase combustion such as some propellant additives or applications dependent on reliable burn rates such as pyrotechnic time delay formulations and primers.
In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model
Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (T g ) and chain length and thus have the potential for phase separation. Here, one crosslinker has a high T g corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low T g . A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.
Sandström and He reanalyzed our 304L creep data (Acta Materialia 121137, 2025) at 10 MPa/750 °C and proposed that: (i) primary dislocation creep dominates based on the ϕ -model; (ii) measured creep rates deviate from published 304H data and Coble creep predictions using literature grain boundary (GB) diffusivity; and (iii) creep-induced elemental redistribution (CIER) could be explained by dislocation emission/absorption (DEAA) and GB creep (GBC), as opposed to diffusional creep. Here, in this response, we propose that: (i) the creep behavior at 15 MPa/750 °C—for which we demonstrated CIER—is only partially captured by the ϕ -model; (ii) direct comparison between 304L and 304H creep data is inappropriate due to carbides in 304H, and Coble creep predictions estimated from CALPHAD-based bulk diffusivity span experimental rates; and (iii) CIER exhibits distinct fingerprint for diffusional creep characterized by anisotropy and stress orientation dependence whereas neither DEAA nor GBC have been mechanistically substantiated.
Electrical transport in semiconducting and metallic particle suspensions is an enabling feature of emerging grid-scale battery technologies. Although the physics of the transport process plays a key role in these technologies, no universal framework has yet emerged. Here, we examine the important contribution of shear flow to the electrical transport of non-Brownian suspensions. We find that these suspensions exhibit a strong dependence of the transport rate on the particle volume fraction and applied shear rate, which enables the conductivity to be dynamically changed by over 10 7 decades based on the applied shear rate. We combine experiments and simulations to conclude that the transport process relies on a combination of charge and particle diffusion with a rate that can be predicted using a quantitative physical model that incorporates the self-diffusion of the particles.
Binary interactions, especially mass transfer and mergers, can strongly influence the evolution of massive stars and change their final properties and the occurrence of supernovae. Here, we investigate how binary interactions affect predictions of the diffuse flux of neutrinos. By performing stellar population syntheses including prescriptions for binary interactions, we show that the resulting detection rates of the diffuse supernova neutrino background is enhanced by 15%–20% compared to estimates without binary considerations. A source of significant uncertainty arises due to the presently sparse knowledge of the evolution of rapidly rotating carbon-oxygen cores, especially those created as a result of mergers near the white dwarf to core collapse boundary. As a result, the enhancement effect may be as small as a few percent if the effects of rotation in postmerger systems are neglected, or as large as 75% if trends are extrapolated. Our estimates serve to highlight that binary effects can be important.
Nickel-based alloys are known for their excellent corrosion resistance and high-temperature strength. However, at high temperatures, polycrystals experience grain boundary-stress-driven vacancy diffusion leading to diffusion-driven creep. Understanding the fundamentals of defect migration is essential to predict diffusion creep behavior. In this talk, we explore the Nabarro-Herring diffusion creep and the associated solute segregation that is likely to result from the different diffusivities of the solute and solvent elements. Density functional theory calculations are performed to determine the formation, migration and binding energies of the solute and vacancy defects in stress-free and stressed states. The dilute-limit transport coefficients are then calculated using the self-consistent mean field theory as implemented in the Kinetic Cluster Expansion (KineCluE) code. Predictions of diffusion creep strain rates and solute segregation are then presented for various solutes, stress and temperatures. This work opens research avenues for exploration of diffusion creep via characterization of elemental segregation to grain boundaries.
Abstract The drying of apple slices in an electrohydrodynamic (EHD) system with a needle‐plate electrode was characterized and mathematically modeled at voltages of 30, 40, and 50 kV and with distances between needles of 5, 6, and 7 cm. The results showed that the drying rate, effective water diffusion coefficient, and rehydration ratio of the apple slices, and specific energy consumption of drying, all increased with increasing voltage. Initially, there was a period when the drying rate was approximately constant in the low‐voltage EHD, but not in the high voltage conditions. At the same voltage, the drying rate was highest when the distance between needles was 6 cm. The voltage and distance between needles had a significant ( p < .05) interactive effect on the effective water diffusion coefficient and specific energy consumption, but not on the rehydration ratio. General mathematical models for drying can be used to represent the drying behavior of apple slices. The Midilli and Kucuk model and the Logarithmic model had the best fits, with correlation coefficients of 0.998803 and 0.998267, respectively. Predictive functions for the moisture ratio of apple slices based on voltage and distance between needles were established and verified. The results indicated that the Logarithmic model could more accurately predict the moisture ratio of apple slices dried using the EHD method. Practical Applications As a new drying technology, electrohydrodynamic (EHD) drying has unequaled advantages. According to the results, using the appropriate voltage and distance between needles improves not only drying efficiency, but also the rehydration ratio of apple slices, which reflects the quality of dried products. The Logarithmic model can accurately predict the EHD drying behavior of apple slices, which is of great significance in optimizing actual production.
Manipulating CdSeTe bandgrading to enhance photocurrent and carrier lifetime is an essential part of high-performance CdTe photovoltaics (PVs). Here, we examine Se diffusion kinetics in single-crystal and polycrystalline CdTe during deposition, thermal annealing, and CdCl 2 treatments. Se distributions are determined by dynamic secondary-ion-mass spectroscopy and Auger electron spectroscopy depth profiling and coupled with electron backscatter diffraction images of the crystalline structure. Effective bulk and grain boundary diffusion coefficients are determined by analytical models and discussed in the context of processing and film morphology. Se is found to diffuse in CdTe at much higher rates during CdCl 2 treatments than with thermal processing alone. GB diffusion also occurs at a significantly faster rate than bulk diffusion. As a result of these two effects, the near interface bulk and GB Se diffusion during CdCl 2 treatments dominates the bandgrading profiles in CdTe PVs.
Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.
A recently discovered, enhanced Ge diffusion mechanism along the oxidizing interface of Si/SiGe nanostructures has enabled the formation of single-crystal Si nanowires and quantum dots embedded in a defect-free, single-crystal SiGe matrix. Here, we report oxidation studies of Si/SiGe nanofins aimed at gaining a better understanding of this novel diffusion mechanism. Here, a superlattice of alternating Si/Si 0.7 Ge 0.3 layers was grown and patterned into fins. After oxidation of the fins, the rate of Ge diffusion down the Si/SiO 2 interface was measured through the analysis of HAADF-STEM images. The activation energy for the diffusion of Ge down the sidewall was found to be 1.1 eV, which is less than one-quarter of the activation energy previously reported for Ge diffusion in bulk Si. Through a combination of experiments and DFT calculations, we propose that the redistribution of Ge occurs by diffusion along the Si/SiO 2 interface followed by a reintroduction into substitutional positions in the crystalline Si.
Polymeric membranes are ubiquitous in transport applications including gas separation, water desalination, solid-state batteries, and fuel cells. The transport mechanism in certain classes of membranes is well understood. It is generally believed that polymeric membranes used in gas separation follow the solution-diffusion model. However, the transport mechanism in other membranes such as the ones used in water desalination is less clear as both the solution-diffusion and hydrodynamic models have been proposed. In this contribution, we study the structure as well as the water and polymer dynamics of several polymeric membranes using small-angle and quasi-elastic neutron scattering. We demonstrate that on the scale of a few nanometers, water appears to diffuse through a desalination membrane at a rate comparable to bulk water, while water diffuses at a rate that is a factor of 2 slower in an anion-exchange membrane. These results illustrate the importance of studying both the structure and dynamics of membrane materials and have important implications in understanding the transport mechanisms that differentiate the different types of polymeric membranes.
Hydrogen diffusion in metals and alloys plays an important role in the discovery of new materials for fuel cell and energy storage technology. While analytic models use hand-selected features that have clear physical ties to hydrogen diffusion, they often lack accuracy when making quantitative predictions. Machine learning models are capable of making accurate predictions, but their inner workings are obscured, rendering it unclear which physical features are truly important. To develop interpretable machine learning models to predict the activation energies of hydrogen diffusion in metals and random binary alloys, we create a database for physical and chemical properties of the species and use it to fit six machine learning models. Our models achieve root-mean-squared errors between 98–119 meV on the testing data and accurately predict that elemental Ru has a large activation energy, while elemental Cr and Fe have small activation energies. By analyzing the feature importances of these fitted models, we identify relevant physical properties for predicting hydrogen diffusivity. While metrics for measuring the individual feature importances for machine learning models exist, correlations between the features lead to disagreement between models and limit the conclusions that can be drawn. Instead grouped feature importance, formed by combining the features via their correlations, agree across the six models and reveal that the two groups containing the packing factor and electronic specific heat are particularly significant for predicting hydrogen diffusion in metals and random binary alloys. In conclusion, this framework allows us to interpret machine learning models and enables rapid screening of new materials with the desired rates of hydrogen diffusion.
Abstract We report a simulation study on the effects of helium (He) bubble size on the morphological evolution and pattern formation on the surface of tungsten used as a plasma-facing component (PFC) in nuclear fusion devices. We have carried out a systematic investigation based on self-consistent dynamical simulations of surface morphological evolution according to an atomistically-informed, 3D continuum-scale model that captures well the relevant length and time scales of surface nanostructure formation in PFC tungsten. The model accounts for PFC surface diffusion, driven by the biaxial compressive stress originating from the over-pressurized He bubbles in the near-surface region of PFC tungsten as a result of He plasma exposure, combined with the formation of self-interstitial atoms in tungsten that diffuse toward the PFC surface and the flux of surface adatoms generated as a result of surface vacancy-adatom pair formation upon He implantation; this transport of surface adatoms contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. Our detailed analysis reveals that varying the average He bubble size plays an important role in the PFC surface growth kinetics as well as the resulting surface topography. Specifically, we find that an increase in the He bubble size leads to a deceleration in the growth rate of the tungsten nanotendrils that emanate from the PFC surface. We also find that the separation distance between the resulting surface features increases with increasing He bubble size, as well as over time. This coarsening effect is a thermally activated process resulting in an accurate description of the temperature dependence of the average surface feature separation by an Arrhenius relation.