Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusivity”

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 253 records · Page 14

Diffusion mechanisms in Ir-coated Re for high-temperature, radiation-cooled rocket thrusters

Materials used for radiation-cooled rocket thrusters must be capable of surviving under extreme conditions of high temperatures and oxidizing environments. Thruster chambers were developed using chemical-vapor-deposited (CVD) Re coated with CVD Ir on the inside surface which is exposed to hot combustion gases. Ir serves as an oxidation barrier protecting the Re which maintains structural integrity at high temperatures. In order to predict and extend the performance limits of these Ir-coated Re thrusters, the diffusion kinetics of CVD materials at temperature are studied. Thruster end ring sections were examined using electron microprobe analysis both before and after exposure to high temperature vacuum environments. The resulting elemental maps for Re, Ir, and Mo in the near-surface region allow identification of diffusion mechanisms operating at these temperatures. Line scans for Ir and Re were fit using a diffusion model to extract relevant diffusion constants. The fastest diffusion process is seen to be grain boundary diffusion with Re diffusing down grain boundaries in the Ir overlayer. The measured dependence of the diffusion rate on temperature will allow prediction of operating lifetimes for these thrusters.

Hamilton, J. C.↗

An analog study of a shock-position diffuser control on a supersonic turbojet engine

Shock position was assumed to be controlled by air bleed through a diffuser-exit bypass door operated by a second-order servomotor. Pressure disturbance resulting from changes in bypass-door position had whereas a disturbance in engine speed had a large effect on diffuser pressure and could impose severe requirements on diffuser control response. The controlled diffuser response improved with faster servomotor response and smaller diffuser dead time and lag. The effect of diffuser dead time on diffuser response was greater than the effect of diffuser lag.

Novik, David↗

Radiant extinction of gaseous diffusion flames

The absence of buoyancy-induced flows in microgravity significantly alters the fundamentals of many combustion processes. Substantial differences between normal-gravity and microgravity flames have been reported during droplet combustion, flame spread over solids, candle flames, and others. These differences are more basic than just in the visible flame shape. Longer residence time and higher concentration of combustion products create a thermochemical environment which changes the flame chemistry. Processes such as flame radiation, that are often ignored under normal gravity, become very important and sometimes even controlling. This is particularly true for conditions at extinction of a microgravity diffusion flame. Under normal-gravity, the buoyant flow, which may be characterized by the strain rate, assists the diffusion process to transport the fuel and oxidizer to the combustion zone and remove the hot combustion products from it. These are essential functions for the survival of the flame which needs fuel and oxidizer. Thus, as the strain rate is increased, the diffusion flame which is 'weak' (reduced burning rate per unit flame area) at low strain rates is initially 'strengthened' and eventually it may be 'blown-out'. Most of the previous research on diffusion flame extinction has been conducted at the high strain rate 'blow-off' limit. The literature substantially lacks information on low strain rate, radiation-induced, extinction of diffusion flames. At the low strain rates encountered in microgravity, flame radiation is enhanced due to: (1) build-up of combustion products in the flame zone which increases the gas radiation, and (2) low strain rates provide sufficient residence time for substantial amounts of soot to form which further increases the flame radiation. It is expected that this radiative heat loss will extinguish the already 'weak' diffusion flame under certain conditions. Identifying these conditions (ambient atmosphere, fuel flow rate, fuel type, etc.) is important for spacecraft fire safety. Thus, the objective is to experimentally and theoretically investigate the radiation-induced extinction of diffusion flames in microgravity and determine the effect of flame radiation on the 'weak' microgravity diffusion flame.

Atreya, Arvind↗

Comments on the Diffusive Behavior of Two Upwind Schemes

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and locally one-dimensional finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2.5 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a speedup of 29 over finite volume.

Wood, William A.↗

Laser Velocimeter Measurements of the Flow Within a Vane-Island Diffuser

A series of detailed non-intrusive laser velocimeter flow field surveys were performed to validate computational fluid dynamic (CFD) codes for rotating pump components. Three component laser-2-focus (L2F) laser velocimetry was used to acquire data within a flow passage of a diffuser at design flow conditions. The thirteen vane-island type diffuser was designed to allow diffusion in the axial direction. Multiple laser velocimeter surveys were performed throughout several diffuser passages and at multiple diffuser depths (hub to shroud). The axial (3D) component of the velocity at all locations was small. The radial and meridional velocity components nearest the diffuser hub increased in magnitude relative to the flow nearest the shroud as the flow progressed through the diffuser. A continuity check across the diffuser throat and based on the meridional velocity component yielded a match with a facility flow meter of 98.0%.

Ferguson, Thomas V.↗

Suppression of Soot Formation and Shapes of Laminar Jet Diffusion Flames

Laminar nonpremixed (diffusion) flames are of interest because they provide model flame systems that are far more tractable for analysis and experiments than practical turbulent flames. In addition, many properties of laminar diffusion flames are directly relevant to turbulent diffusion flames using laminar flamelet concepts. Finally, laminar diffusion flame shapes have been of interest since the classical study of Burke and Schumann because they involve a simple nonintrusive measurement that is convenient for evaluating flame shape predictions. Motivated by these observations, the shapes of round hydrocarbon-fueled laminar jet diffusion flames were considered, emphasizing conditions where effects of buoyancy are small because most practical flames are not buoyant. Earlier studies of shapes of hydrocarbon-fueled nonbuoyant laminar jet diffusion flames considered combustion in still air and have shown that flames at the laminar smoke point are roughly twice as long as corresponding soot-free (blue) flames and have developed simple ways to estimate their shapes. Corresponding studies of hydrocarbon-fueled weakly-buoyant laminar jet diffusion flames in coflowing air have also been reported. These studies were limited to soot-containing flames at laminar smoke point conditions and also developed simple ways to estimate their shapes but the behavior of corresponding soot-free flames has not been addressed. This is unfortunate because ways of selecting flame flow properties to reduce soot concentrations are of great interest; in addition, soot-free flames are fundamentally important because they are much more computationally tractable than corresponding soot-containing flames. Thus, the objectives of the present investigation were to observe the shapes of weakly-buoyant laminar jet diffusion flames at both soot-free and smoke point conditions and to use the results to evaluate simplified flame shape models. The present discussion is brief.

Xu, F.↗

Diffusion Characteristics of Upwind Schemes on Unstructured Triangulations

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and dimensionally-split finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the dimensionally-split finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2-3 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a 20-25 speedup over finite volume.

Wood, William A.↗

Microgravity Turbulent Gas-Jet Diffusion Flames

A gas-jet diffusion flame is similar to the flame on a Bunsen burner, where a gaseous fuel (e.g., propane) flows from a nozzle into an oxygen-containing atmosphere (e.g., air). The difference is that a Bunsen burner allows for (partial) premixing of the fuel and the air, whereas a diffusion flame is not premixed and gets its oxygen (principally) by diffusion from the atmosphere around the flame. Simple gas-jet diffusion flames are often used for combustion studies because they embody the mechanisms operating in accidental fires and in practical combustion systems. However, most practical combustion is turbulent (i.e., with random flow vortices), which enhances the fuel/air mixing. These turbulent flames are not well understood because their random and transient nature complicates analysis. Normal gravity studies of turbulence in gas-jet diffusion flames can be impeded by buoyancy-induced instabilities. These gravitycaused instabilities, which are evident in the flickering of a candle flame in normal gravity, interfere with the study of turbulent gas-jet diffusion flames. By conducting experiments in microgravity, where buoyant instabilities are avoided, we at the NASA Lewis Research Center hope to improve our understanding of turbulent combustion. Ultimately, this could lead to improvements in combustor design, yielding higher efficiency and lower pollutant emissions. Gas-jet diffusion flames are often researched as model flames, because they embody mechanisms operating in both accidental fires and practical combustion systems (see the first figure). In normal gravity laboratory research, buoyant air flows, which are often negligible in practical situations, dominate the heat and mass transfer processes. Microgravity research studies, however, are not constrained by buoyant air flows, and new, unique information on the behavior of gas-jet diffusion flames has been obtained.

Source record↗

Effects of g-Jitter on Diffusion in Binary Liquids

The microgravity environment offers the potential to measure the binary diffusion coefficients in liquids without the masking effects introduced by buoyancy-induced flows due to Earth s gravity. However, the background g-jitter (vibrations from the shuttle, onboard machinery, and crew) normally encountered in many shuttle experiments may alter the benefits of the microgravity environment and introduce vibrations that could offset its intrinsic advantages. An experiment during STS-85 (August 1997) used the Microgravity Vibration Isolation Mount (MIM) to isolate and introduce controlled vibrations to two miscible liquids inside a cavity to study the effects of g-jitter on liquid diffusion. Diffusion in a nonhomogeneous liquid system is caused by a nonequilibrium condition that results in the transport of mass (dispersion of the different kinds of liquid molecules) to approach equilibrium. The dynamic state of the system tends toward equilibrium such that the system becomes homogeneous. An everyday example is the mixing of cream and coffee (a nonhomogeneous system) via stirring. The cream diffuses into the coffee, thus forming a homogeneous system. At equilibrium the system is said to be mixed. However, during stirring, simple observations show complex flow field dynamics-stretching and folding of material interfaces, thinning of striation thickness, self-similar patterns, and so on. This example illustrates that, even though mixing occurs via mass diffusion, stirring to enhance transport plays a major role. Stirring can be induced either by mechanical means (spoon or plastic stirrer) or via buoyancy-induced forces caused by Earth s gravity. Accurate measurements of binary diffusion coefficients are often inhibited by buoyancy-induced flows. The microgravity environment minimizes the effect of buoyancy-induced flows and allows the true diffusion limit to be achieved. One goal of this experiment was to show that the microgravity environment suppresses buoyancy-induced convection, thereby mass diffusion becomes the dominant mechanism for transport. Since g-jitter transmitted by the shuttle to the experiment can potentially excite buoyancy-induced flows, we also studied the effects of controlled vibrations on the system.

Duval, Walter M. B.↗

Generating a Simulated Fluid Flow Over an Aircraft Surface Using Anisotropic Diffusion

A fluid-flow simulation over a computer-generated aircraft surface is generated using a diffusion technique. The surface is comprised of a surface mesh of polygons. A boundary-layer fluid property is obtained for a subset of the polygons of the surface mesh. A pressure-gradient vector is determined for a selected polygon, the selected polygon belonging to the surface mesh but not one of the subset of polygons. A maximum and minimum diffusion rate is determined along directions determined using a pressure gradient vector corresponding to the selected polygon. A diffusion-path vector is defined between a point in the selected polygon and a neighboring point in a neighboring polygon. An updated fluid property is determined for the selected polygon using a variable diffusion rate, the variable diffusion rate based on the minimum diffusion rate, maximum diffusion rate, and angular difference between the diffusion-path vector and the pressure-gradient vector.

Rodriguez, David L.↗

Atomistic and mesoscale simulations to determine effective diffusion coefficient of fission products in SiC

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) particles serves as the barrier to prevent escape of fission products produced in the fuel kernel. Knowing the diffusion coefficient of fission products through SiC is critical to determining whether fission gas can escape from the particle. It has been observed in experiments that Ag accumulated in grain boundaries and triple junctions in SiC. It is hypothesized that grain boundary diffusion is the primary pathway by which fission products penetrate the SiC layer. In this report, the effective diffusion coefficient of the fission product Ag through the grain boundary network is calculated using a combination of atomistic and phase-field methods. The grain boundary diffusion coefficient is calculated using molecular dynamics simulations. The bulk diffusion coefficient is determined using a combination of density functional theory and nudged elastic band methods. An effective diffusion coefficient is calculated, accounting for the grain structure using a phase-field method. The effective diffusion coefficient will be incorporated into Bison and fission product release calculations are compared to available experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Sluggish Diffusion of Cations in CeO 2 Probed through Molecular Dynamics and Metadynamics Simulations

Cation diffusion in fluorite‐structured CeO 2 , though far slower than anion diffusion, is an important, high‐temperature process because it governs diverse fabrication and degradation phenomena. Herein, cation diffusion is studied by means of classical molecular dynamics and metadynamics simulations. Three different mechanisms are examined: migration involving an isolated cerium vacancy, migration involving a cerium vacancy in a defect associate with an oxygen vacancy, and migration involving a cation divacancy. For each mechanism, defect diffusion coefficients are calculated as a function of temperature, from which the respective activation enthalpy of defect migration is obtained. Through comparisons with experimental cation diffusion data (specifically, of the absolute magnitude of the cation diffusivity as well as its activation enthalpy), it is concluded that cation diffusion takes place predominantly neither by isolated vacancies nor by cation vacancy–oxygen vacancy associates but by cation divacancies.

36 MATERIALS SCIENCE↗

HTO and selenate diffusion through compacted Na-, Na–Ca-, and Ca-montmorillonite

Radionuclide transport in smectite clay barrier systems used for nuclear waste disposal is controlled by diffusion, with adsorption significantly retarding transport rates. While a relatively minor component of spent nuclear fuel, 79 Se is a major driver of the safety case for spent fuel disposal due to its long half-life (3.3×10 5 yr) and its low adsorption to clay (K D < 10 L/kg), thus a thorough understanding of Se diffusion through clay is critical for understanding the long-term safety of spent fuel disposal systems. Through-diffusion experiments with tritiated water (HTO, conservative tracer) and Se(VI) were conducted with a well-characterized, purified montmorillonite source clay (SWy-2) under a constant ionic strength (0.1 M) and three different electrolyte compositions: Na + , Ca 2+ , and a Na + -Ca 2+ mixture at pH 6.5 in order to probe the effects of electrolyte composition and interlayer cation composition on clay microstructure, Se(VI) aqueous speciation, and ultimately diffusion. Further, the results were modeled using a reactive transport modeling approach to determine values of porosity (ε), D e (effective diffusion coefficient), and K D (distribution coefficient for adsorption). HTO diffusive flux was higher in Ca-montmorillonite (D e =1.68×10 -10 m 2 s -1 ) compared to Na-montmorillonite (De=7.83×10 -11 m 2 s -1 ). This increase in flux is likely due to a greater degree of clay layer stacking in the presence of Ca 2+ compared to Na + , which leads to larger inter-particle pores. Overall, the Se(VI) flux was much lower than the HTO flux due to anion exclusion, with Se(VI) flux following the order Ca (D e = 1.03×10 -11 m 2 s -1 ) > Na–Ca (D e = 2.12×10 -12 m 2 s -1 ) > Na (D e = 1.28×10 -12 m 2 s -1 ). These differences in Se(VI) flux are due to a combination of factors, including (1) larger accessible porosity in Ca-montmorillonite due to clay layer stacking and smaller electrostatic effects compared to Na-montmorillonite, (2) larger accessible porosity for neutral-charge CaSeO4 species which makes up 32% of aqueous Se(VI) in the pure Ca system, and (3) possibly higher Se(VI) adsorption for Ca-montmorillonite. Through a combination of experimental and modeling work, this study highlights the compounding effects that electrolyte and counterion compositions can have on radionuclide transport through clay. Diffusion models that neglect these effects are not transferable from laboratory experimental conditions to in situ repository conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗

Mechanistic calculation of the effective silver diffusion coefficient in polycrystalline silicon carbide: Application to silver release in AGR-1 TRISO particles

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission products produced and not retained in the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110m Ag isotope. However, accurately determining the fission gas release rate requires knowing the diffusion coefficient through the SiC layer. In this study, we leverage atomistic calculations of Ag diffusivity in SiC bulk and grain boundaries (GBs) to develop a mesoscale effective Ag diffusion coefficient (D eff ) in SiC. Since GBs serve as pathways for Ag diffusion, D eff is defined as a function of temperature and microstructure variables. In particular, the size of SiC grains in the direction perpendicular to diffusion is shown to significantly affect Ag diffusion. The prediction of the mechanistic, mesoscale approach falls within one order of magnitude of empirical values. The temperature and microstructure-dependent effective Ag diffusivity in SiC is implemented in the fuel performance code Bison with a correction factor to predict Ag release from AGR-1 TRISO fuel particles. We hereby quantify the impact of SiC grain size on Ag release and improve Bison’s predictions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multiscale study of helium diffusion in Ni-Cr alloys: Short-range trapping versus long-range channeling

Ni-Cr alloys are widely employed as structural materials in fast nuclear reactors but are vulnerable to high-temperature helium (He) embrittlement (HTHE) under fast neutron irradiation. A comprehensive understanding of He diffusion in Ni-Cr alloys, which governs the kinetics of HTHE, is therefore essential for developing resilient materials and preventing failure. In this work, we reveal the underlying mechanisms of He diffusion in pure Ni and Ni-Cr alloys by integrating density functional theory (DFT) with atomic kinetic Monte Carlo (AKMC) simulations. Our findings uncover a non-monotonic dependence of He diffusivity on Cr concentration, contradicting the monotonic trends predicted by DFT-parameterized theories. At low Cr concentrations, He diffusion is dominated by short-range trapping, characterized by multiple trapping sites and a distinct mechanism within the first nearest neighbor of Cr, differing from that in pure Ni. At high Cr concentrations, these local traps become interconnected, forming long-range fast diffusion channels that enhance He mobility. The competition between localized trapping and extended channeling results in a diffusivity that first decreases, then increases with rising Cr content. These atomic-scale insights offer critical guidance for the design of radiation-tolerant Ni-based alloys. Moreover, the combined DFT-AKMC methodology and the concept of random walker diffusion through interconnected energy basins present a broadly applicable framework for studying transport phenomena in disordered systems.

36 - MATERIALS SCIENCE↗

Cation–polymer interactions and local heterogeneity determine the relative order of alkali cation diffusion coefficients in PEGDA hydrogels

Current research efforts are focused on endowing polymer membranes with ion–ion selectivity by incorporating ion–polymer interactions into materials to bias the selective partitioning and or diffusivity of one species over another. However, little is known about the impact of such interactions on the mechanisms of ion transport. In this study, we probe the influence of cation–polymer interactions on cation, anion, and salt diffusivity in a model membrane material, poly(ethylene glycol) diacrylate (PEGDA) by modeling concentrated polyethylene oxide solutions via molecular dynamics simulations. These results are compared to published experimental data for LiCl, NaCl, and KCl diffusion in PEGDA. Experimentally, the order of salt and cation diffusion coefficients for LiCl, NaCl, and KCl deviate from the order in aqueous solutions. Here, simulations identify these deviations to arise from cation–polymer coordination in the membrane. Both the fraction of bound cations and the average binding lifetime increases with decreasing cation hydration free energy (moving down the alkali series), leading to different diffusivity trends in the membrane compared to solution. However, to recover the experimentally observed order of diffusivities cations and salt in our simulations, we needed to incorporate membrane heterogeneity explicitly via a polymer charge scaling procedure. Together, our results indicate that cation–polymer interactions, as well as spatial heterogeneity within the membrane, play a critical role in dictating the observed order of alkali cation and salt diffusion coefficients in membranes.

36 MATERIALS SCIENCE↗

($\mathrm{INVITED}$)Counter-ion effect on the diffusion behavior of $\mathrm{Y}$b, $\mathrm{L}$u, and $\mathrm{N}$d ions in $\mathrm{YAG}$ transparent ceramics

The ability to fabricate additively manufactured laser waveguides with sharp dopant concentration interfaces is limited by diffusion of the dopants at the temperatures required to fully densify the material. Compositional analysis of bilayer samples, where each layer was either undoped YAG or YAG doped with Yb, Lu, or Nd, were fabricated such that all combinations were available for testing. Samples were fabricated at both 1750°C and 1850°C to determine the diffusion behavior of each dopant alone and also in the presence of a second dopant. It was found that the experimental concentration profiles exhibited both intragranular (bulk) and grain boundary contributions, and thus fitting to a complementary error function equation required the use of two diffusion coefficients. Nd always diffused further along grain boundaries than the other dopants, due in part to its small segregation coefficient in YAG. It is shown that the presence of Nd as a counter dopant inhibits intragranular diffusion of other dopants while enhancing their grain boundary diffusion. All the observed trends were attributed to a combination of intragranular lattice strain due to: dopant ions replacing yttrium substitutionally, the relative driving forces for the segregation of dopant ions to grain boundaries, and the ability of one dopant to affect the diffusion of a different dopant in the other direction to maintain charge neutrality.

36 MATERIALS SCIENCE↗