Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion behavior”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Effects of Rhenium Addition on the Temporal Evolution of the Nanostructure and Chemistry of a Model Ni-Cr-Al Superalloy: Analysis of the Coarsening Behavior - 2

The temporal evolution of the nanostructure and chemistry of a model Ni-8.5 at.% Cr-10 at.% Al alloy with the addition of 2 at.% Re was studied using transmission electron microscopy and atom-probe tomography in order to measure the number density and mean radius of the y' (LIZ) precipitates and the chemistry of the y'-precipitates and the y (fcc)-matrix. In this article, the coarsening behavior of the y'-precipitates is discussed in detail and compared with the Umantsev-Olson model for multi-component alloys. In addition, the experimental results are evaluated with PrecipiCalc(TradeMark) simulations. The results show that the diffusivities of the solute elements play a major role in the coarsening behavior of the y'-precipitates and that the addition of Re retards the coarsening kinetics and stabilizes the spheroidal morphology of the precipitates by reducing the interfacial energy.

Yoon, Kevin E.↗
DOE OSTI2024

Insights into the structure and dynamics of K + ions at the muscovite–water interface from machine learning potential simulations

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. Here, in this work, we make a first step in this direction by investigating the behavior of the native K + ions at the interface between neat water and the muscovite mica (001) surface with ab-initio -based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K + ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K + diffusion between mica’s ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K + ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Ab-initio molecular dynamics↗
DOE OSTI2023

A Generalization of Threshold-Based and Probability-Based Models of Information Diffusion

Diffusion of information through complex networks is of interest in studies such as propagation prediction and influence maximization, both of which have applications in viral marketing and rumor controlling. There are a variety of information diffusion models, all of which simulate the adoption and spread of information over time. However, there is a lack of understanding of whether, despite their conceptual differences, these models represent the same underlying generative structures. For instance, if two different models utilize different conceptual mechanisms, but generate the same results, does the choice of model matter? A classification of diffusion of information models is developed based on the neighbor knowledge of the model infection requirement and the stochasticity of the model. This classification allows for the identification of models that fall into each respective category. The study involves the analysis of the following agent-based models on directed scale-free networks: (1) a linear absolute threshold model (LATM), (2) a linear fractional threshold model (LTFM), (3) the independent cascade model (ICM), (4) Bass-Rand-Rust model (BRRM) (5) a stochastic linear absolute threshold model (SLATM) (6) a stochastic fractional threshold model (SLFTM), and (7) Dodds–Watts model (DWM). Through the execution of simulations and analysis of the experimental results, the distinctive properties of each model are identified. Our analysis reveals that similarity in conceptual design does not imply similarity in behavior concerning speed, final state of nodes and edges, and sensitivity to parameters. Therefore, we highlight the importance of considering the unique behavioral characteristics of each model when selecting a suitable information diffusion model for a particular application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Investigating intermolecular interactions among CO 2 , water and PEEK-ionene membrane using cryo ToF-SIMS and isotopic labeling

Cryogenic time-of-flight secondary ion mass spectrometry (cryo ToF-SIMS) has emerged as a powerful tool for investigating molecular interactions, speciation, and dynamics in materials for CO 2 capture. In this study, we apply cryo ToF-SIMS to probe interactions between CO 2 , water, and PEEK-ionene membranes—a promising material for direct CO 2 capture due to its selectivity, durability, and efficiency. Despite this potential, the mechanisms governing CO 2 diffusion and the influence of water vapor on CO 2 behavior remain unclear. To address this, we loaded PEEK-ionene membranes with 13 CO 2 and D 2 O and employed cryo ToF-SIMS to visualize the 3D distribution of CO 2 and water within the membrane. While prior studies suggest that 13 CO 2 is absorbed under ambient conditions, our cryo ToF-SIMS analysis revealed no enhancement of the 13 C/ 12 C ratio, suggesting weak CO 2 -membrane interactions. As a result, CO 2 vaporizes even at low temperatures (−140°C) under vacuum conditions. In contrast, D 2 O displayed a relatively homogeneous distribution in the membrane, suggesting stronger water-membrane interactions via hydrogen bonding (18–20 kJ/mol). Interestingly, CO 2 was not detected in D 2 O-loaded membranes, indicating minimal interference from water vapor on CO 2 diffusion. As a comparison, the cryo ToF-SIMS data show that 13 CO 2 can readily react with a basic Na 2 CO 3 aqueous solution to form NaH 13 CO 3 . These findings demonstrate cryo ToF-SIMS as a critical technique for understanding gas-water-membrane interactions, offering insights for membrane functionalization to improve CO 2 capture efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2023

Oxidation of Silicon Particles Suspended in Mullite Matrix With Dual-Mode Vacancy Diffusion of Oxygen

Mullite bond coats filled with silicon particles have recently been proposed for the next generation of environmental barrier coatings (EBCs). The sacrificial oxidation of the silicon particles eliminates some of the environmental durability issues associated with traditional EBC designs. To investigate this novel design approach, a mathematical formulation and resulting numerical method have been developed for simulating the oxidation of silicon particles embedded in a mullite matrix. The focus is on the oxidation of the silicon in a dry oxygen environment where the oxidant is presumably oxygen. Oxygen transport from the exterior environment through the mullite is assumed to occur via two lattice defect diffusion mechanisms. As a result, the rate of oxygen transport is dictated by the self-diffusivity of the lattice defects. The numerical method is implemented using the COMSOL Multiphysics® program (COMSOL, Inc.). The method is applied to simulate the oxidation of silicon particles in a mullite pellet in recent oxidation experiments. Numerical results compare favorably with the oxidation behavior observed in the oxidation experiments. The numerical method is also used to investigate the effects of various features of the mullite-silicon particle system on the oxidation behavior.

silicon oxidation↗
DOE OSTI2022

A Robust Numerical Treatment of Solid-Phase Diffusion in Pseudo Two-Dimensional Lithium-Ion Battery Models

Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust numerical algorithm for solving solid-phase diffusion equations in physics-based battery models. In this work, we discuss the origins of numerical instabilities that can occur when solving the solid-phase diffusion equations using iterative methods. Then, in order to resolve such issues, we propose a simple numerical treatment to the surface flux term of discretized solid-phase diffusion equations. To demonstrate its numerical robustness, the proposed method is implemented into a pseudo two-dimensional (P2D) physics-based battery model and simulations are conducted at wide ranges of operating conditions. Even with extremely poor initial guesses for the Li+ concentrations of the active materials, computations using the proposed method do not diverge and the their computational speeds are comparable to those with conventional initial guesses. Comprehensive tests of the proposed method are also performed with a dynamic current profile based on US06 driving profile and a multi-stage charging profile with very high initial C-rate (12C).

battery modeling↗
DOE OSTI2023

Wall diffuser velocity effects on American Shad ( Alosa sapidissima ) inside a fishway entrance channel

Attraction water for fishways is typically introduced through a diffuser inside the entrance channel, often through the floor or wall. In the spring of 2019, this laboratory study examined how 151 adult American Shad (Alosa sapidissima) responded to different gross velocities through a wall diffuser inside a full-scale fishway entrance channel. Two velocity conditions were studied, 0.152 m/s and 0.305 m/s, both without turning vanes inside the auxiliary water channel. The fish were tracked using the passive integrated transponder telemetry technique. The results of the experiments showed no difference in American Shad behavior when exhibited to the low and high velocity treatments. Moreover, shad passed the diffuser in roughly 3 out of every 4 attempts, regardless of the treatment. However, the similarity in shad behavior and passage performance is believed to be more of a result of the similarity in flow fields that resulted from the lack of flow guidance devices inside the auxiliary water channel. Furthermore, these findings therefore highlight the importance of properly maintained flow guidance devices, an often-overlooked component of an auxiliary water system.

13 HYDRO ENERGY↗
DOE OSTI2025

Experimental observation of nonlinear relation between pressure and water flux is consistent with the solution-diffusion model

In several recent studies, it has been proposed that the fundamental understanding of penetrant transport in dense polymer membranes occurring via the solution-diffusion model, which has been the generally accepted theoretical framework for describing penetrant transport in such materials for the past several decades, is flawed. An alternate mechanistic framework based on the idea of two-phase flow in a porous medium (i.e., pore-flow) has been broadly advanced instead, with proponents of this approach claiming that the pore-flow theoretical framework provides the necessary mechanistic insight to design novel polymeric membrane materials for emerging applications. In this study, we show experimental results for hydraulic permeation of water that are entirely consistent with the solution-diffusion theory, without modification, for three dense polymeric membranes: crosslinked poly(ethylene glycol diacrylate) (XLPEGDA), Nafion 117 ionomer in the sodium counterion form (Nafion 117-Na), and cellulose acetate (CA). By measuring water flux at transmembrane pressures up to 240 bar, we observe a nonlinear relationship between the transmembrane pressure (TMP) and water flux, J w , for XLPEGDA and Nafion 117-Na, while this relationship is linear for CA. We demonstrate that the behavior of these three materials is described via the solution-diffusion model. According to the solution-diffusion model, flux is, to a good approximation, proportional to the transmembrane concentration difference induced by the pressure difference across the membrane, rather than to TMP itself. Water sorption isotherms are reported for all three materials. They further justify the nonlinear relationship between TMP and J w observed in XLPEGDA and Nafion 117-Na, emphasizing that the nonlinearity in the flux/TMP relationship stems from nonlinearities in the sorption isotherm with pressure. Additionally, the relationship between water flux and TMP can be predicted, a priori, with no adjustable parameters when a predictive model for the diffusion coefficient of water is employed in conjunction with the experimental water sorption isotherms in the solution-diffusion model. Furthermore, our results demonstrate the validity of the solution-diffusion model to describe transport of penetrants in dense polymer membranes, while highlighting the sensitivity of the solution-diffusion model to the many physical and mathematical simplifications commonly applied to the theory in literature.

materials↗
DOE OSTI2021

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Composite Behavior of Nanopore Array Large Memristors

Synthetic nanopores were recently demonstrated with memristive and nonlinear voltage-current behaviors, akin to ion channels in a cell membrane. Such ionic devices are considered a promising candidate for the development of brain-inspired neuromorphic computing techniques. In this work, we show the composite behavior of nanopore-array large memristors, formed with different membrane materials, pore sizes, electrolytes, and device arrangements. Anodic aluminum oxide (AAO) membranes with 5 nm and 20 nm diameter pores and track-etched polycarbonate (PCTE) membranes with 10 nm diameter pores are tested and shown to demonstrate memristive and nonlinear behaviors with approximately 107–1010 pores in parallel when electrolyte concentration across the membranes is asymmetric. Ion diffusion through the large number of channels induces time-dependent electrolyte asymmetry that drives the system through different memristive states. The behaviors of series composite memristors with different configurations are also presented. In addition to helping understand fluidic devices and circuits for neuromorphic computing, the results also shed light on the development of field-assisted ion-selection-membrane filtration techniques as well as the investigations of large neurons and giant synapses. Further work is needed to de-embed parasitic components of the measurement setup to obtain intrinsic large memristor properties.

Chemistry↗
DOE OSTI2024

Melting and solidification dynamics during laser melting of reaction-based metal matrix composites uncovered by in-situ synchrotron X-ray diffraction

Laser additive manufacturing (AM) of reaction-based metal matrix composites (MMCs) involves highly complex and non-equilibrium material transformation behavior, including melting, dissolution, precipitation, and solidification. Yet, the dynamics and interplay of these phase transformation processes remain poorly understood, posing substantial challenges in identifying the microstructure formation mechanism, and predicting and controlling the microstructure in the printed parts. Here, in this study, we performed the in-situ X-ray diffraction experiment to characterize the phase evolution dynamics of the 316L + 10 vol.%TiC system during laser melting, which provides direct and quantitative insights of the complex phase reaction and evolution dynamics under rapid heating and cooling conditions relevant to additive manufacturing of reaction-based MMCs. Further in-depth thermodynamic and kinetic calculations revealed that most of the phase evolution behavior observed in the in-situ X-ray diffraction experiment cannot be solely explained by widely used equilibrium thermodynamic models, and diffusion-controlled nonequilibrium dissolution and precipitation kinetics must be considered to elucidate the complex phase evolution behavior, including incomplete TiC dissolution, and three-step TiC precipitation. The three distinct types of precipitates generate unique hierarchical TiC micro- and nanostructures, which enhances the yield strength from 513 MPa to 877 MPa by 71 %, tensile strength from 628 MPa to 1054 MPa by 68 %, and Young's modulus from 193 GPa to 221 GPa by 14 %. The findings of our research provide the knowledge foundation for the design of unique microstructures and advanced MMC materials through laser AM.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Synchrotron-driven Instabilities in Relativistic Plasmas of Arbitrary Opacity

Recent work has shown that synchrotron emission from relativistic plasmas leads the electron distribution to form an anisotropic ring in momentum space, which can be unstable to both kinetic and hydrodynamic instabilities. Fundamental to these works was the assumption that the plasma was optically thin, allowing all emitted radiation to escape. Here, we examine the behavior of these instabilities as the plasma becomes more optically thick. To do this, we extend a recently developed Fokker–Planck operator for synchrotron emission and absorption in mildly relativistic plasmas to ultrarelativistic plasmas. For a given set of plasma parameters, photons emitted by higher-energy electrons tend to be higher frequency, and thus more easily escape the plasma. As a result, the ratio of the photon emission rate (radiative drag) to absorption rate (radiative diffusion) for a given electron is extremely energy dependent. Given this behavior, we determine the critical parameters that control the opacity, and show how the plasma gradually transitions to become more isotropic and stable at higher opacity.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1977

Kinetics data for diffusion of outgas species from RTV 560

A detailed analytical and experimental study was made of the outgassing behavior of RTV 560 silicone rubber. The four outgas species which predominate in the temperature range of 285 K to 425 K were separately identified. The initial concentration of these species in the parent material and their bulk volatilities were determined. The diffusion coefficients and activation energy for diffusion of the two major species were deduced from outgassing rate data. It is shown that by using these data in a diffusion theory model, the outgassing rates of these major species can be predicted for arbitrary geometry and any temperature within the range studied.

Liu, C. K.↗
NASA NTRS1987

Oxygen-18 tracer study of the passive thermal oxidation of silicon

This work focuses on the thermal oxidation of silicon near 1273 K using the double-tracer oxidation method. The results confirm that oxidation occurs by the transport of electrically neutral nonnetwork oxygen through the interstitial space of the vitreous silica (v-SiO2) scale. Simultaneously, self- (or isotopic-) diffusion occurs in the network, resulting in characteristic isotopic fraction distributions near the gas-scale interface. The self-diffusion coefficients calculated from these profiles agree with those reported for tracer diffusion in v-SiO2, and the diffusion coefficient calculated from the scale growth is consistent with reported O2 permeation data. An important parameter that describes the double-oxidation behavior is a ratio relating the scale thickness grown during the second oxidation, the network self-diffusion coefficient for oxygen, and the time of the second oxidation.

Cawley, J. D.↗
DOE OSTI2024

Effects of neutron irradiation on hydrogen isotope permeation on W-Re

Hydrogen and deuterium permeation behavior for neutron damaged W and W-3 %Re was evaluated using a Plasma Driven Permeation device. For neutron damaged W-3 %Re at 773 K and 1473 K, precipitation of Re and Os was observed by atom probe tomography. Neutron irradiation introduced defect damage throughout the sample and enhanced H and D diffusion toward downstream side. The existence of irradiation damage by neutron irradiation was found to increase hydrogen isotope permeation. Moreover, the existence of Re has suppressed the formation of irradiation defects for hydrogen isotope trapping sites. Further, the recovery of irradiation defects by increasing the irradiation temperature led to the structural changes due to the diffusion of Re and Os, which also affected the hydrogen isotope permeation behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of MOF linker rotation and functionalization on methane uptake and diffusion

The flexible degrees of freedom in metal–organic frameworks (MOFs) can have significant effects on guest molecule behavior. However, in the majority of studies applying molecular simulations to MOFs, the framework is assumed to be rigid in order to minimize computational cost. Here, in this work, we assess the significance of this assumption on a representative example of methane uptake and diffusion in UiO-66. We introduce an open-source code to modify MOFs through functionalization and linker rotation and we perform Grand Canonical Monte Carlo and molecular dynamics simulations of methane in each of the functionalized and linker-rotated derivatives of UiO-66. We find that linker rotation moderately influences methane uptake and significantly influences methane diffusion. Our assessment provides ranges of property values that serve as measures of uncertainty of these two properties associated with linker rotation. We further determine that void volume fraction and minimum pore size are the features that govern methane uptake and diffusion, respectively. These findings illustrate the impact of linker rotation on MOFs and provide design principles to guide future investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

Modeling Candle Flame Behavior In Variable Gravity

The burning of a candle, as typical non-propagating diffusion flame, has been used by a number of researchers to study the effects of electric fields on flame, spontaneous flame oscillation and flickering phenomena, and flame extinction. In normal gravity, the heat released from combustion creates buoyant convection that draws oxygen into the flame. The strength of the buoyant flow depends on the gravitational level and it is expected that the flame shape, size and candle burning rate will vary with gravity. Experimentally, there exist studies of candle burning in enhanced gravity (i.e. higher than normal earth gravity, g(sub e)), and in microgravity in drop towers and space-based facilities. There are, however, no reported experimental data on candle burning in partial gravity (g < g(sub e)). In a previous numerical model of the candle flame, buoyant forces were neglected. The treatment of momentum equation was simplified using a potential flow approximation. Although the predicted flame characteristics agreed well with the experimental results, the model cannot be extended to cases with buoyant flows. In addition, because of the use of potential flow, no-slip boundary condition is not satisfied on the wick surface. So there is some uncertainty on the accuracy of the predicted flow field. In the present modeling effort, the full Navier-Stokes momentum equations with body force term is included. This enables us to study the effect of gravity on candle flames (with zero gravity as the limiting case). In addition, we consider radiation effects in more detail by solving the radiation transfer equation. In the previous study, flame radiation is treated as a simple loss term in the energy equation. Emphasis of the present model is on the gas-phase processes. Therefore, the detailed heat and mass transfer phenomena inside the porous wick are not treated. Instead, it is assumed that a thin layer of liquid fuel coated the entire wick surface during the burning process. This is the limiting case that the mass transfer process in the wick is much faster than the evaporation process at the wick surface.

Alsairafi, A.↗
DOE OSTI2021

The influence of preexisting host rock damage on fault network localization

The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.

58 GEOSCIENCES↗