Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion coefficient”

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 523 records · Page 29

Controlled release of entrapped nanoparticles from thermoresponsive hydrogels with tunable network characteristics

Thermoresponsive hydrogels have been studied intensively for creating smart drug carriers and controlled drug delivery. Understanding the drug release kinetics and corresponding transport mechanisms of nanoparticles (NPs) in a thermoresponsive hydrogel network is the key to the successful design of smart drug delivery systems. Here, we construct a mesoscopic model of rigid NPs entrapped in a hydrogel network in an aqueous solution, where the hydrogel network is formed by cross-linked semiflexible polymers of thermoresponsive poly(N-isopropylacrylamide) (PNIPAM). By varying the temperature crossing the lower critical solution temperature of PNIPAM, we can significantly change the hydrogel network characteristics. We systematically investigate how the matrix porosity and the nanoparticle size affect the transport kinetics of NPs at different temperatures. Quantitative results on the mean-squared displacement and the van Hove displacement distributions of NPs show that all NPs entrapped in the smart hydrogels undergo subdiffusion at both low and high temperatures. For a coil state, the transport of NPs in the hydrogels can be enhanced by decreasing the matrix porosity of the polymer network and NPs' size. However, when the solution temperature is increased above the critical temperature, the hydrogel network collapses following the coil-to-globule transition, with the NPs tightly trapped in some local regions inside the hydrogels. Consequently, the NP diffusion coefficient can be reduced by two orders of magnitude, or the diffusion processes can even be completely stopped. These findings provide new insights for designing controlled drug release from stimuli-responsive hydrogels, including autonomously switch on/off drug release in response to physical and chemical stimuli.

42 ENGINEERING↗

Mg and native defects in cubic silicon carbide from first principles

The diffusion of Mg defects in 3C-SiC is studied using the density functional theory. Mg has the highest burn-in rate as a transmutant in 3C-SiC when it is placed in high-energy neutron irradiation environment of a fusion reactor. The presence and evolution of transmutant defects impact thermal and mechanical properties of this important structural material. This study is focused on understanding the structure, stability, and evolution of Mg defects and the interaction of Mg with native defects in 3C-SiC. Our calculations of diffusion coefficients for different Mg defects suggest that Mg is likely to diffuse faster in pristine 3C-SiC than in the damaged one, in agreement with earlier experimental observations.

36 MATERIALS SCIENCE↗

High-frequency electrostatic waves in the magnetosphere.

High-frequency electrostatic microinstabilities in magnetospheric plasmas are considered in detail. Rather special plasma parameters are found to be required to match the theoretical wave spectrum with satellite observations in the magnetosphere. In particular, it is necessary to have a cold and a warm species of electrons such that (1) the warm component has an anomalous velocity distribution function that is nonmonotonic in the perpendicular component of velocity and is the source of free energy driving the instabilities, (2) the density ratio of the cold component to the hot component is greater than about 0.01, and (3) the temperature ratio of the two components for cases of high particle density is no less than 0.1. These requirements and the corresponding instability criteria are satisfied only in the trapping region; this is also the region in which the waves are most frequently observed. The range of unstable wavelengths and an estimate of the diffusion coefficient are also obtained. The wave are found to induce strong diffusion in velocity space for low-energy electrons during periods of moderate wave amplitude.

Young, T. S. T.↗

Multi-scale modeling of wastage layer formation in metallic fuel cladding

Fuel-cladding chemical interaction (FCCI) is a major concern for U-Zr metallic fuels' performance, primarily due to the formation of a brittle layer (wastage) in the cladding. This brittle layer, resulting from intermetallic compounds between cladding constituents Fe, Cr, and lanthanide fission products, significantly impacts the cladding's mechanical integrity. Recent efforts focus on developing a mechanistic modeling framework to understand lanthanide production, transport to the fuel-cladding interface, and phase transformation to intermetallic phases. A multi-scale computational approach has been used to calculate lanthanide transport rates, with atomistic calculations determining Nd diffusivities through the solid fuel matrix and along pore surfaces. These diffusivities inform a mesoscale model to determine an effective diffusion coefficient, accounting for porosity and infiltration with bond sodium. This effective diffusivity is used in engineering-scale simulations via the BISON fuel performance code, which has been validated against EBR-II and FFTF reactor experiments.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Transport coefficients of partially ionized argon.

Properties of partially ionized Ar computed via Chapman-Enskog-Burnett expressions noting ambipolar diffusion coefficient, electron-atom momentum, electrical and thermal conductivity, etc

TRANSPORT COEFFICIENT↗

Cluster dynamics simulation of xenon diffusion during irradiation in UO 2

Diffusion of fission gas in UO 2 nuclear fuel impacts several important performance metrics, such as fission gas release, swelling, and thermal conductivity. Current empirical models of fission gas release have significant uncertainty, some of which derives from the bulk diffusion rate and its dependence on, for example, fuel chemistry and irradiation. In this work, we have applied the previously-developed Free Energy Cluster Dynamics (FECD) methodology in the code Centipede to calculate xenon cluster concentrations in UO 2 under intrinsic (high temperature) and irradiation-enhanced (intermediate temperature) conditions in order to develop a model of the xenon diffusion coefficient based on the atomic scale mechanisms responsible for transport. While the diffusion mechanism for xenon in UO 2 is adequately described by the Xe + U 2 O vacancy cluster for intrinsic conditions, a similar process is not capable of capturing measured in-pile fission gas diffusivity at intermediate temperatures. Therefore, a different diffusion mechanism must dominate under this regime. Using calculated atomistic data, we have shown that irradiation-enhanced diffusion at intermediate temperatures occurs via the larger Xe + U 4 O y vacancy clusters, which have lower migration barriers and increase in concentration by several orders of magnitude compared to intrinsic conditions. This mechanism is enabled by the increased uranium vacancy concentration under irradiation due to Frenkel pair production. In addition, the fast migration of uranium interstitials with two attached oxygen interstitials lowers the total uranium interstitial concentration through reactions with sinks. This allows the extended defects, such as Xe + U 4 O y vacancy clusters, to maintain high concentrations by limiting annihilation with attached vacancies. Furthermore, predictions using the Xe + U 4 O y diffusion mechanism are in good agreement with experiment, albeit with some differences in the Arrhenius slope, which we believe may be related to either experimental or model parameter uncertainty. Lastly, an analytical expression suitable for application in fuel performance simulations was derived to capture the predictions of the Centipede simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Grain Boundary Effects on Radiation Damage and Tritium Diffusion in Li–Al–O Ceramics from Molecular Dynamics and Experiments

Understanding irradiation damage and tritium transport in LiAlO2 ceramics is essential for their deployment in tritium-producing burnable absorber rods (TPBARs). Grain boundaries (GBs) play an important role in governing radiation response and tritium transport in LiAlO2 and its secondary phase LiAl5O8. Using molecular dynamics simulations, we investigated defect evolution and tritium diffusion during displacement cascades in single-crystal and bicrystal LiAlO2 and LiAl5O8 at 600 K. The results reveal that GBs suppress the damage in the LiAlO2 bicrystals as compared to single crystals, via interstitial emission. The effect is more pronounced for Li defects, than for Al and O. In LiAl5O8 the damage reduction due to introduction of GB is less pronounced overall. Tritium diffusion coefficients at GBs are higher by factors of 2–10 relative to diffusion inside grain bulk. The effect is particularly pronounced in LiAlO2, where S5, S17, and S25 GBs promote rapid tritium transport, whereas LiAl5O8 exhibits slower diffusion due to reduced free volume at its GBs. These findings evince a trade-off that while GBs mitigate radiation damage by absorbing interstitials, they simultaneously provide fast pathways for tritium migration which is not desirable for TPBAR application. The mechanistic insights gained here establish a foundation for microstructural design strategies to balance radiation tolerance and tritium retention in ceramic breeder materials.

Roy, Ankit↗

Rotational and translational diffusion of liquid n-hexane: EFP-based molecular dynamics analysis

Molecular Dynamics (MD) simulations based on the Effective Fragment Potential (EFP) method are utilized to provide a comprehensive assessment of diffusion in liquid n-hexane. We decompose translational diffusion into components along and orthogonal to the long axis of the molecule. Rotational diffusion is decomposed into tumbling and spinning motions about this axis. Our analysis yields four corresponding diffusion coefficients which are related to diagonal entries in the complete 6 × 6 diffusion tensor accounting for the three rotational and three translational degrees of freedom and for the potential coupling between them. However, coupling between different degrees of freedom is expected to be minimal for a natural choice of the molecular body-fixed axis, so then off-diagonal entries in the tensor are negligible. This expectation is supported by a hydrodynamic analysis of the diffusion tensor which treats the liquid surrounding the molecule being tracked as a viscous continuum. Thus, the EFP MD analysis provides a comprehensive characterization of diffusion and also reveals expected shortcomings of the hydrodynamic treatment, particularly for rotational diffusion, when applied to neat liquids.

74 ATOMIC AND MOLECULAR PHYSICS↗

High-throughput approach for investigating interdiffusion in medium- and high-entropy alloys

Interdiffusion experiments are usually time-consuming and tedious since diffusion couples must be annealed at several temperatures for a long time. The efforts required to study interdiffusion in multicomponent alloys increase dramatically as multiple diffusion couples are required to cover broad composition ranges and determine the diffusivities of individual elements in different chemical environments. To circumvent this challenge, we present a high-throughput approach applicable to single-phase and compositionally complex alloys, which are assumed to approximate ideal solid solutions. Here, a simple diffusion-multiple experiment combined with a physically based kinetic model is proposed to efficiently determine the diffusion coefficients of the constituent elements in quaternary CrFeCoNi alloys. Compared with tracer diffusivities reported in the literature, the results, thus, obtained do not differ by more than a factor of 2 and were obtained from a single interdiffusion experiment. In contrast, the diffusivities simulated with commercial mobility and thermodynamic databases are strongly overestimated by a factor ranging from 1 to 16. Therefore, our approach enables high-throughput determination of diffusivities and can help in the design of alloys for high-temperature applications where diffusion plays a key role.

36 MATERIALS SCIENCE↗

Accelerated Materials Design for Molten Salt Technologies Using Innovative High-Throughput Methods

The focus of the project is on building an innovative accelerated materials design platform for molten salt technologies using novel high-throughput methods coupled to data analytics. The main objectives is to predict a FeCrMnNi alloy compositional space with better corrosion resistance than stainless steel 316 and identify new molten salt corrosion mechanisms. The project demonstrates the feasibility to use high-throughput methods coupled to data analytics to accelerate alloy design for extreme environments applications. Using a trained and tested machine learning (ML) model, 2000 FeCrMnNi alloy corrosion rate in molten chloride salts were predicted and a compositional field with corrosion rate lower than 316 stainless steel was identified. The ML model interpretability unveiled multiple features of importance in the model prediction. Some features were expected to be of relative significance, such as work function, surface energy and alloy electronegativity, and the ML model interpretability analysis confirmed those. On the other hand, the most important feature is the diffusion coefficient of Ni in the bulk alloy which indicates that a surface diffusion mechanism plays an important role in the overall corrosion mechanism in molten salts.

36 MATERIALS SCIENCE↗

Ring current impoundment of the Io plasma torus

A newly discovered feature in the Io plasma formation that may be described as a ramp separating a high-density plasma ledge on its Jupiterward side from the lower-density radially distended Io plasma disc on its anti-Jupiterward side is observed to coincide with a marked inward decrease in the ring current population. The spatial congruency of the counter-directed maximal gradients in both plasma bodies reveals a profound coupling between them. The existence of the ramp requires a local order-of-magnitude reduction in the diffusion coefficient that governs radial mass transport. It is demonstrated that the diminished diffusive efficiency there is caused by strong pressure gradient inhibition of the interchange instability that underlies mass transport. The Io plasma torus, which is defined as the region of strong ultraviolet emissions, is identified as the plasma ledge. The plasma density in the ledge is high and, incidentally therefore, able to emit strongly because it is impounded against rapid, centrifugal expulsion by the inwardly directed pressure of the ring current at its inner edge.

Siscoe, G. L.↗

Concentrations and behavior of oxygen and oxide ion in melts of composition CaO.MgO.xSiO2

The behavior of oxygen and oxide ion in silicate melts was investigated through their electrochemical reactions at a platinum electrode. Values are given for the diffusion coefficient for molecular oxygen in diopside melt and the activation energy of diffusion. It is shown that molecular oxygen dissociates prior to undergoing reduction and that oxide ion reacts quickly with silicate polymers when it is produced. The concentration of oxide ion is kept low by a buffering effect of the silicate, the exact level being dependent on the silicate composition. Data on the kinetics of reaction of the dissociation of molecular oxygen and on the buffering reactions are provided. It is demonstrated that the data on oxygen in these silicate melts are consistent with those for solid buffers.

Semkow, K. W.↗

The combined effects of stress and moisture on the structural integrity of composites

A variational approach is used to model the moisture content and distribution profiles which cause swelling in polymers and composites which are used as aircraft structural components. The moisture intrusion stems from humidity, which makes the problem one of diffusion of mass from the environment into a semiinfinite plate with a given diffusion coefficient and thickness. Lagrangian, ordinary differential equations are defined for the mass diffusion, including driving force, mass potential and volume dissipation expressions. Account is taken of the effects of stress on diffusivity, including the dilatational stress caused by inhomogeneous swelling of the material. It is noted that material near the top and bottom of a part undergoes the most swelling, and the evolution of moisture diffusion in the materials is traced. Finally, experimental data are used to define differential equations for predicting the final moisture distribution within the material.

Hiel, C.↗

Orientation and velocity dependence of the nonequilibrium partition coefficient

Monte Carlo simulations based on a Spin-1 Ising Model for binary alloys have been used to investigate the non-equilibrium partition coefficient (k(sub neq)) as a function of solid-liquid interface velocity and orientation. In simulations of Si with a second component k(sub neq) is greater in the [111] direction than the [100] direction in agreement with experimental results reported by Azlz et al. The simulated partition coefficient scales with the square of the step velocity divided by the diffusion coefficient of the secondary component in the liquid.

Beatty, K. M.↗

Adsorption and Surface Diffusion of Metals on α-Al2O3 for Advanced Manufacturing Applications

Abstract The adsorption and diffusion of Mo and Nb adatoms on the $$\alpha$$ α -Al 2 O 3 (0001) surface were explored using density functional theory-based methods. Adsorption energies of Mo and Nb adatoms at minima sites on the surface were determined. Mo and Nb adatoms prefer to adsorb to the same locations on the surface, and larger adsorption energies calculated for Nb compared to Mo indicate that Nb adatom-surface interactions are stronger than Mo. Using minima adsorption sites as initial and final locations for surface diffusion, energy barriers for diffusion were calculated using the nudged elastic band method. Overall, Mo and Nb follow roughly the same diffusion paths. The diffusion pre-factors for Mo and Nb are similar; however, Mo diffusion has a lower energy barrier and thus a larger diffusion coefficient compared with Nb. These results provide insight into the role of surface diffusion of Mo and Nb adatoms during advanced manufacturing processes.

Biaggne, Austin (ORCID:0000000328074960)↗

Hydrogen Diffusion in Slit Pores: Role of Temperature, Pressure, Confinement, and Roughness

Diffusion of hydrogen (H 2 ) is important to understand the leakage risk and transport behavior for H 2 geologic storage. We applied molecular dynamics simulations to investigate the influencing factors of H2 diffusion in the slit pores of calcite, hematite, and quartz, owing to their abundance. It is revealed that the H2 self-diffusion coefficient increases with the temperature, regardless of the type of pore minerals. The diffusion of H 2 in the 20 nm slit pores falls into the bulk diffusion regime when the pressure is 10 MPa. The self-diffusion of H 2 decreases with pressure in all three types of slit pores, following a power law model with the exponents ranging from -0.825 to -0.964. Furthermore, the impact of confinement on H 2 diffusion is more pronounced for the slit pores with stronger interactions with H 2 -like calcite. The role of surface roughness in H 2 diffusion depends on the slit aperture. The rough surface enhances H 2 diffusion in the larger slit pores due to the enlarged effective pore space, whereas it weakens H 2 diffusion in the small slit pores due to stronger adsorption. These findings will fill the knowledge gap on the coupling effect of different factors influencing H 2 diffusion.

08 HYDROGEN↗

Diffusion in polymer solutions.

Diffusion coefficient of ordinary low molecular weight solutes in polymer solutions, noting variation of diffusivity from value in pure solvent

Gainer, J. L.↗

The Fokker-Planck coefficient for pitch-angle scattering of cosmic rays

For the case of homogeneous, isotropic magnetic field fluctuations, it is shown that most theories which are based on the quasi-linear and adiabatic approximation yield the same integral for the Fokker-Planck coefficient for the pitch angle scattering of cosmic rays. For example, despite apparent differences, the theories due to Jokipii and to Klimas and Sandri yield the same integral. It is also shown, however, that this integral in most cases has been evaluated incorrectly in the past. For large pitch angles these errors become significant, and for pitch angles of 90 deg the actual Fokker-Planck coefficient contains a delta function. The implications for these corrections relating cosmic ray diffusion coefficients to observed properties of the interplanetary magnetic field are discussed.

Fisk, L. A.↗