Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Physically constrained 3D diffusion for inverse design of fiber-reinforced polymer composite materials

Designing fiber-reinforced polymer composites (FRPCs) with a tailored nonlinear stress-strain response is crucial for applications such as energy absorption in crash structures, flexible robotics, and impact-resistant protective gear. However, the inherent complexities of composite materials and the multitude of parameters involved, render traditional design and optimization methods inadequate for achieving effective inverse design of composites. In this paper, we present an AI-based inverse design framework that effectively and efficiently generates FRPCs with targeted nonlinear stress-strain responses. We introduce a physically constrained diffusion model (PC3D_Diffusion) capable of managing the complexities of composite materials and producing detailed, high-quality designs. We propose a loss-guided, learning-free approach to generate physically feasible microstructure designs by explicitly enforcing physical constraints during the generation process. For training purposes, 1.35 million FRPC samples were created, and their corresponding stress-strain curves were computed using established physics-based computational models. The results show that PC3D_Diffusion consistently generates high-quality designs with tailored mechanical behaviors, while guaranteeing compliance with the physical constraints. PC3D_Diffusion advances FRPC inverse design and may facilitate the inverse design of other 3D materials, offering potential applications in industries reliant on materials with custom mechanical properties.

Xu, Pei [Clemson Univ., SC (United States)]↗
DOE OSTI2019

Multigrid treatment of implicit continuum diffusion

Implicit treatment of diffusive terms of various differential orders common in continuum mechanics modeling, such as computational fluid dynamics, is investigated with spectral and multigrid algorithms in non-periodic 2D domains. In doubly periodic time dependent problems these terms are handled efficiently by spectral methods, but in non-periodic systems solved with distributed memory parallel computing and 2D domain decomposition, this efficiency is lost for a large number of processors. Here, we built and present here a multigrid algorithm for these types of problems that outperforms a spectral solution employing the highly optimized FFTW library. This solver is suitable for high performance computing and may be able to efficiently treat implicit diffusion of arbitrary order by introducing auxiliary equations of lower order. We test these solvers for fourth and sixth order diffusion with harmonic test functions as well as turbulent 2D magnetohydrodynamic simulations. It is also shown that an anisotropic operator without mixed-derivative terms improves model accuracy and speed, and we examine the impact that the various diffusion operators have on the energy, the enstrophy, and the qualitative aspect of a simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Enabling fast charging of lithium-ion batteries through secondary- /dual- pore network: Part I - Analytical diffusion model

Battery performance is strongly correlated with electrode microstructural properties. Enabling fast charging of lithium-ion batteries requires improved through-plane ionic diffusion that can be achieved through, among other strategies, structured electrodes with a secondary- or dual-pore network (SPN). In this work, an analytical model investigates the impact of such an SPN on ionic diffusion with a composite electrode, considering various pore-channel geometries and comparing to standard electrodes with identical gravimetric- and volumetric-specific theoretical capacities. Relevant SPN design parameters and tortuosity coefficients are identified according to three optimization objectives that aim to balance the improved overall through-plane diffusion, thanks to the coarse aligned channels, and degraded in-plane diffusion because of the porous matrix densification required to maintain gravimetric- and volumetric-specific theoretical capacities. The model indicates that a relatively low amount of SPN is required and that electrodes with high through-plane tortuosity and low in-plane tortuosity benefit most from such architecture.

25 ENERGY STORAGE↗
DOE OSTI2022

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Diffusive mass transport in agglomerated glassy fallout from a near-surface nuclear test

Aerodynamically-shaped glassy fallout is formed when vapor phase constituents from the nuclear device are incorporated into molten carriers (i.e. fallout precursor materials derived from soil or other near-field environmental debris). The effects of speciation and diffusive transport of condensing constituents are not well defined in models of fallout formation. Previously we reported observations of diffuse micrometer scale layers enriched in Na, Fe, Ca, and 235 U, and depleted in Al and Ti, at the interfaces of agglomerated fallout objects. Here in this paper, we derive the timescales of uranium mass transport in such fallout as it cools from 2500 K to 1500 K by applying a 1-dimensional planar diffusion model to the observed 235 U/ 30 Si variation at the interfaces. By modeling the thermal transport between the fireball and the carrier materials, the time of mass transport is calculated to be <0.6 s, <1 s, <2 s, and <3.5 s for fireball yields of 0.1 kt, 1 kt, 10 kt, and 100 kt respectively. Based on the calculated times of mass transport, a maximum temperature of deposition of uranium onto the carrier material of ~2200 K is inferred (1σ uncertainty of ~200 K). We also determine that the occurrence of micrometer scale layers of material enriched in relatively volatile Na-species as well as more refractory Ca-species provides evidence for an oxygen-rich fireball based on the vapor pressure of the two species under oxidizing conditions. These results represent the first application of diffusion-based modeling to derive material transport, thermal environments, and oxidation-speciation in near-surface nuclear detonation environments.

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

Diffusion in random networks

The ensemble averaging technique is applied to model mass transport by diffusion in random networks. The system consists of an ensemble of random networks, where each network is made of pockets connected by tortuous channels. Inside a channel, fluid transport is assumed to be governed by the one-dimensional diffusion equation. Mass balance leads to an integro-differential equation for the pocket mass density. The so-called dual-porosity model is found to be equivalent to the leading order approximation of the integration kernel when the diffusion time scale inside the channels is small compared to the macroscopic time scale. As a test problem, we consider the one-dimensional mass diffusion in a semi-infinite domain. Because of the required time to establish the linear concentration profile inside a channel, for early times the similarity variable is xt $-$1/4 rather than xt $-$1/2 as in the traditional theory. We found this early time similarity can be explained by random walk theory through the network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A Method for Including Diffusive Effects in Texture Evolution

This work describes a model for representing grain misorientation in continuum models for texture evolution. The model relies on a connection, derived here, between the fluctuation of the reorientation velocity in discrete models of texture evolution and an isotropic diffusion operator that can be used to represent misorientation effects in continuum models for texture evolution. Previous work has developed a method for computing the fluctuation of the reorientation velocity in discrete models. These values, along with the average reorientation velocity, are used to set up inhomogeneous continuum advection and diffusion operators. Here, the complete homogenized continuum model uses the discrete harmonics to represent and evolve the orientation distribution function, now including a diffusion operator representing the effects of intragranular misorientation. Results from the new model are compared to fully-resolved crystal plasticity finite element simulations and experimental data to demonstrate the accuracy of the proposed diffusion operator and the overall method for simulating texture evolution.

(B) crystal plasticity↗
DOE OSTI2017

Molecular dynamics simulations of NMR relaxation and diffusion of bulk hydrocarbons and water

Molecular dynamics (MD) simulations are used to investigate 1 H nuclear magnetic resonance (NMR) relaxation and diffusion of bulk n-C 5 H 12 to n-C 17 H 36 hydrocarbons and bulk water. The MD simulations of the 1H NMR relaxation times T 1,2 in the fast motion regime where T 1 = T 2 agree with measured (de-oxygenated) T2 data at ambient conditions, without any adjustable parameters in the interpretation of the simulation data. Likewise, the translational diffusion DT coefficients calculated using simulation con gurations are well-correlated with measured diffusion data at ambient conditions. The agreement between the predicted and experimentally measured NMR relaxation times and diffusion coefficient also validate the force elds used in the simulation. The molecular simulations naturally separate intramolecular from intermolecular dipole-dipole interactions helping bring new insight into the two NMR relaxation mechanisms as a function of molecular chain-length (i.e. carbon number). Comparison of the MD simulation results of the two relaxation mechanisms with traditional hard-sphere models used in interpreting NMR data reveals important limitations in the latter. With increasing chain length, there is substantial deviation in the molecular size inferred on the basis of the radius of gyration from simulation and the tted hard-sphere radii required to rationalize the relaxation times. This deviation is characteristic of the local nature of the NMR measurement, one that is well-captured by molecular simulations.

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

Revisiting the diffusion mechanism of helium in UO 2 : A DFT+ U study

The understanding of migration properties of helium atoms after their generation through α-decay of actinides in spent nuclear fuels is important for the safety of nuclear fuel storage and disposal. The diffusion of helium in UO 2 is revisited by using the DFT+U simulation methodology employing the “U-ramping” method to address the issue of metastable energy states. A novel diffusion mechanism by helium interstitials, the “asymmetric hop” mechanism, is reported and compared to other diffusion mechanisms including an oxygen vacancy mediated mechanism and available experimental diffusion data. We show that the new mechanism is the dominant one over a wide temperature range.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Diffusion in undoped and Cr-doped amorphous UO 2

UO 2 fuel pellets are often doped with chromium oxide to obtain favourable properties such as higher density, improved thermal stability, large grain sizes, improved pellet-clad interaction margins, and increased fission gas retention during transients. Chromium has a low solubility limit in UO 2 , with past experimental work reporting solubility limits ranging between 0.004 to 0.06 wt.% Cr. Due to its low solubility, segregation of Cr ions to the grain boundary may occur. Further, the complexity of these boundaries may be high as observed in other ceramics resulting in disordered or amorphous regions along the boundary, affecting a range of material and operational properties of the fuel pellet. To assess these disordered regions, in this work we study amorphous undoped and Cr doped UO 2 systems (containing 10–50 at.% Cr 3+ ) that have been modelled using classical molecular dynamics methods incorporating Cr 3+ into the well-used CRG potential library. Diffusion coefficients, pre-exponential factors, and activation energies for diffusion were computed for oxygen ions, assessing the impact of structure and extrinsic species on migration. Oxygen diffusion was observed to be much faster in the undoped amorphous system compared to its crystalline counterpart. Oxygen diffusion in doped systems decreased with increasing Cr concentration, highlighting the importance of additives to retain fission products and other migratory species.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Simulations of self- and Xe diffusivity in uranium mononitride including chemistry and irradiation effects

A combination of density functional theory and empirical potential atomic scale simulations have been used to determine a model for defect stability and mobility in uranium mononitride (UN), as a function of temperature (T) and N 2 partial pressure (p N 2 ). Using the model, predictions of hypo-stoichiometry under U-rich conditions compare favorably to CALPHAD calculations using the TAF-ID database. Furthermore, our predictions of U and N self-diffusivity are in good agreement with experiments carried out as a function of T at specific partial pressures under thermal equilibrium. The validated atomic scale data have then been implemented within a cluster dynamics method to simulate irradiation-enhanced defect concentrations. All defects and clusters studied have significantly enhanced concentrations, with respect to thermal equilibrium, as T is lowered. The irradiation-enhanced Xe diffusivity is compared to post-irradiation annealing and in-pile experiments. In conclusion, the contributions of various defects and clusters to non-stoichiometry, self-diffusivity, and Xe diffusivity are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Grain boundary self-diffusion and point defect interactions in α -U via molecular dynamics

Though metallic U-Zr fuel has been used in nuclear reactors since the 1960s, many of its fundamental and thermodynamic properties are still unknown. The a-U phase, which has a highly anisotropic crystal structure and physical properties, is present in U-Zr fuel. The character and behavior of a-U grain boundaries will strongly impact fuel thermophysical performance under irradiation. Here, we study the interaction of point defects with grain boundaries, diffusion along grain boundaries, and the predicted diffusional creep behavior of a-U via molecular dynamics. We calculate the segregation energy of vacancies and interstitials to grain boundaries and quantify the biased sink strength of the grain boundaries, and observe that this sink strength is not strongly dependent on the grain boundary orientation. We also find that grain boundary diffusivity is strongly dependent on the grain boundary energy and grain boundary orientation. The presence of point defects within the grain boundary can induce diffusion in grain boundaries with low formation energies and can enhance diffusion in high-energy grain boundaries. We also find that diffusional creep of a-U at prototypical metallic fuel operation conditions is extremely high and could help explain observed metallic fuel swelling behaviors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling Ni redistribution in the hydrogen electrode of solid oxide cells through Ni(OH) 2 diffusion and Ni-YSZ wettability change

Ni redistribution in the hydrogen electrodes of solid oxide cells is an important degradation mechanism. Its driving force is still under debate. This work focuses on the interplay between Ni(OH) 2 diffusion and Ni-YSZ wettability change. With a Ni(OH) 2 diffusion model developed in our previous work, we further employ three models for the Ni-YSZ contact angle to quantify the Ni-YSZ wettability change. The microstructure evolutions in a reconstructed Ni-YSZ electrode are simulated under selected experimental operating conditions. It is shown that the phenomenological model can capture the driving force of Ni spreading/detachment on YSZ surface and the Ni migration reported in experiments can be qualitatively reproduced through the competition between Ni(OH) 2 diffusion and Ni-YSZ wettability change. It is also found that the initial microstructure of the active layer and the microstructure characteristics in the support layer can strongly affect the distribution of steam partial pressure and overpotential, hence the Ni redistribution, in the active layer, which may explain the inconsistencies in experiments. Further, our results show that there are two missing pieces in current theory of Ni redistribution: a mechanism of fast Ni diffusion under humid atmosphere and the physical reason behind the Ni spreading/detachment.

25 ENERGY STORAGE↗
DOE OSTI2020

Quantitative nanoscale tracking of oxygen vacancy diffusion inside single ceria grains by in situ transmission electron microscopy

Oxygen vacancy formation and migration in ceria is critical to its electrochemical and catalytic properties in systems for chemical and energy transformation, but its quantification is rather challenging especially at atomic-scale because of disordered distribution. Here we report a rational approach to track oxygen vacancy diffusion in single grains of pure and Sm-doped ceria at -20 °C to 160 °C using in situ (scanning) transmission electron microscopy ((S)TEM). To create a gradient in oxygen vacancy concentration, a small region (~30 nm in diameter) inside a ceria grain is reduced to the C-type CeO 1.68 phase by the ionization or radiolysis effect of a high-energy electron beam. The evolution in oxygen vacancy concentration is then mapped through lattice expansion measurement using scanning nano-beam diffraction or 4D STEM at a spatial resolution better than 2 nm; this allows direct determination of local oxygen vacancy diffusion coefficients in a very small domain inside pure and Sm-doped ceria at different temperatures. Further, the activation energies for oxygen transport are determined to be 0.59, 0.66, 1.12, and 1.27 eV for pure CeO 2 , Ce 0.94 Sm 0.06 O 1.97 , Ce 0.89 Sm 0.11 O 1.945 , and Ce 0.8 Sm 0.2 O 1.9 , respectively, implying that activation energy increases due to impurity scattering. The results are qualitatively supported by density functional theory (DFT) calculations. In addition, our in situ TEM investigation reveals that dislocations impede oxygen vacancy diffusion by absorbing oxygen vacancies from the surrounding areas and pinning them locally. With more oxygen vacancies absorbed, dislocations show extended strain fields with local tensile zone sandwiched between the compressed ones. Therefore, dislocation density should be reduced in order to minimize the resistance to oxygen vacancy diffusion at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Excited-State Diffusion Monte Carlo Calculations: A Simple and Efficient Two-Determinant Ansatz

We perform excited-state variational Monte Carlo and diffusion Monte Carlo calculations using a simple and efficient wave function ansatz. This ansatz follows the recent variation-after-response formalism, accurately approximating a configuration interaction singles wave function as a sum of only two nonorthogonal Slater determinants and further including important orbital relaxation. The ansatz is used to perform diffusion Monte Carlo calculations with large augmented basis sets, comparing to benchmarks from near-exact quantum chemical methods. Here, the significance of orbital optimization in excited-state diffusion Monte Carlo is demonstrated, and the excited-state optimization procedure is discussed in detail. Diffuse excited states in water and formaldehyde are studied, in addition to the formaldimine and benzonitrile molecules.

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

Diffusivity of Carbon Dioxide in Aqueous Solutions under Geologic Carbon Sequestration Conditions

Accurate assessment of the long-term security of geologic carbon sequestration requires knowledge of the mobility of carbon dioxide in brines under pressure and temperature conditions that prevail in subsurface aquifers. In this paper, we report Raman spectroscopic measurements of the rate of CO 2 diffusion in water and brines as a function of pressure, salinity, and concentration of CO 2 . In pure water at 50 ± 2 °C and 90 ± 2 bar, we find the diffusion coefficient, D, to be (3.08 ± 0.03) × 10 –9 m 2 /s, a value that is consistent with a recent microfluidic study but lower than earlier PVT measurements. Under reservoir conditions, salinity affects the mobility of CO 2 significantly and D decreased by 45% for a 4 M solution of NaCl. We find significant differences of diffusivity of CO 2 in brines (0–4 M NaCl), in both the absolute values and the trend compared to the Stokes–Einstein prediction under our experimental conditions. Finally, we observe that D decreases significantly at the high CO 2 concentrations expected in subsurface aquifers (~15% reduction at 0.55 mol/kg of CO 2 ) and provides an empirical correction to the commonly reported D values that assume a tracer concentration dependence on diffusivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗