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

Computing eddy-driven effective diffusivity using Lagrangian particles

A novel method to derive effective diffusivity from Lagrangian particle trajectory data sets is developed and then analyzed relative to particle-derived meridional diffusivity for eddy-driven mixing in an idealized circumpolar current. Quantitative standard dispersion- and transport-based mixing diagnostics are defined, compared and contrasted to motivate the computation and use of effective diffusivity derived from Lagrangian particles. We compute the effective diffusivity by first performing scalar transport on Lagrangian control areas using stored trajectories computed from online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) using the Model for Prediction Across Scales Ocean (MPAS-O). Furthermore, the Lagrangian scalar transport scheme is compared against an Eulerian scalar transport scheme. Spatially-variable effective diffusivities are computed from resulting time-varying cumulative concentrations that vary as a function of cumulative area. The transport-based Eulerian and Lagrangian effective diffusivity diagnostics are found to be qualitatively consistent with the dispersion-based diffusivity. All diffusivity estimates show a region of increased subsurface diffusivity within the core of an idealized circumpolar current and results are within a factor of two of each other. The Eulerian and Lagrangian effective diffusivities are most similar; smaller and more spatially diffused values are obtained with the dispersion-based diffusivity computed with particle clusters.

54 ENVIRONMENTAL SCIENCES↗

Performance Evaluation of Fixed- and Variable-Area Rocket Exhaust Diffusers Using Single and Clustered Nozzles with and Without Gimbaling

An investigation of exhaust diffusers used for altitude simulation in testing rocket engines was conducted with model diffusers and gaseous nitrogen as the working fluid. This investigation was conducted to eval­uate the effects on performance of gimbaling clustered nozzles in a single fixed-area exhaust diffuser. A two-nozzle cluster was gimbaled in all attitudes in four different exhaust-diffuser configurations: (1) a straight circular tube, (2) a straight figure-eight tube, (3) a circular tube with a second throat, and (4) a figure-eight tube with a second throat. Gimbaling clustered nozzles had little effect on the operating pressure ratio for either straight-tube or second-throat diffusers, but caused a large increase in the starting-pressure-ratio requirements in some straight-tube exhaust diffusers. The performance of two or four clustered nozzles with no gimbaling was compared with single-nozzle per­formance on the basis of the ratio of diffuser area to nozzle-throat area. The performance was found to be dependent on the diffuser- to nozzle-throat-area ratio and independent of both the number of primary nozzles and the nozzle-area ratio. A separate investigation was also conducted to evaluate the per­formance improvement obtainable with a variable-area exhaust diffuser. This type of diffuser achieved a 26-percent reduction in starting pressure ratio, a 40-percent reduction in operating pressure ratio, and a 50-percent reduction in the overall diffuser length over the values attainable with a fixed-area second-throat diffuser.

ROCKET↗

Methane Diffusion and Adsorption in Shale Rocks: A Numerical Study Using the Dusty Gas Model in TOUGH2/EOS7C-ECBM

Gas production from shale gas reservoirs plays a significant role in satisfying increasing energy demands. Compared with conventional sandstone and carbonate reservoirs, shale gas reservoirs are characterized by extremely low porosity, ultra-low permeability and high clay content. Slip flow, diffusion, adsorption and desorption are the primary gas transport processes in shale matrix, while Darcy flow is restricted to fractures. Understanding methane diffusion and adsorption, and gas flow and equilibrium in the low-permeability matrix of shale is crucial for shale formation evaluation and for predicting gas production. Modeling of diffusion in low-permeability shale rocks requires use of the Dusty gas model (DGM) rather than Fick’s law. The DGM is incorporated in the TOUGH2 module EOS7C-ECBM, a modified version of EOS7C that simulates multicomponent gas mixture transport in porous media. Also included in EOS7C-ECBM is the extended Langmuir model for adsorption and desorption of gases. In this study, a column shale model was constructed to simulate methane diffusion and adsorption through shale rocks. The process of binary CH 4 - N 2 diffusion and adsorption was analyzed. A sensitivity study was performed to investigate the effects of pressure, temperature and permeability on diffusion and adsorption in shale rocks. The results show that methane gas diffusion and adsorption in shale is a slow process of dynamic equilibrium, which can be illustrated by the slope of a curve in CH4mass variation. The amount of adsorption increases with the pressure increase at the low pressure, and the mass change by gas diffusion will decrease due to the decrease in the compressibility factor of the gas. With the elevated temperature, the gas molecules move faster and then the greater gas diffusion rates make the process duration shorter. The gas diffusion rate decreases with the permeability decrease, and there is a limit of gas diffusion if the permeability is less than 1.0×10 -15 m 2 . In conclusion, the results can provide insights for a better understanding of methane diffusion and adsorption in the shale rocks so as to optimize gas production performance of shale gas reservoirs.

04 OIL SHALES AND TAR SANDS↗

Diffusion and Flow of CO{sub 2} in Carbon Anode for Aluminium Smelting

Using a well-made carbon anode with the right porosity characteristics is essential to the successful operation of a Hall–Héroult cell. Bubble generation at this carbon anode and its contribution to the overall voltage drop in aluminum production holds significant potential for reducing overpotential and improving energy savings. This voltage drop is believed to be greatly influenced by the gas diffusion in the anode carbon, for which there are a limited number of measurements to correlate with the carbon anode properties. In the present study, the CO{sub 2} gas diffusion characteristics were first investigated via measuring anode porosity using mercury intrusion porosimetry (MIP) for different samples with different permeability. Second, diffusion experiments were conducted by flowing CO{sub 2} gas into the anode sample at elevated temperatures (up to 960 °C). The impact of temperature, average pore size, and permeability on the diffusion coefficient was investigated. The diffusion coefficient was consequently calculated by curve fitting using diffusion theories. The value obtained varied from 1.38 × 10{sup −6} to 7.89 × 10{sup −6} m{sup 2}/s. A high diffusion coefficient was measured in the carbon sample with larger average pore size and higher permeability. It was found that anomalous diffusion behavior in the temperature range 600 °C to 960 °C was caused by convective flow effects, different diffusion mechanisms, and the Boudouard reaction for these carbon anode samples. Furthermore, by plotting the measured diffusion coefficient against average pore size and permeability, an increasing trend of diffusivity with an increase in pore size and permeability was observed.

36 MATERIALS SCIENCE↗

Diffusion coefficients of transition metals in fcc cobalt

Using first-principles density functional theory (DFT), we calculate the diffusivities of 32 different solute elements—all transition metals, together with Al and Si—in fcc cobalt within the formalism of the five-frequency model. For self-diffusion in fcc cobalt, we compare the accuracy of various approximations to the exchange-correlation energy functional of DFT in estimating the activation energy, and find that only the Perdew-Burke-Ernzerhof (PBE) approximation agrees well with experimental reports and all other functionals largely overestimate it. Our calculations also show that an accurate estimation of the self-diffusion coefficient requires explicit calculation of the effective jump frequency and vacancy formation entropy via phonons. Using accurate self-diffusion data and scaling all solute-related attempt frequencies with respect to the attempt frequency for self-diffusion using a simple relation involving the atomic mass and melting temperature of the solute yields solute diffusivities in excellent agreement with experiments, where such data is available. Additionally, we find that large solutes spontaneously relax toward the nearest neighbor vacancy to relieve the misfit strain, and the extent of this relaxation correlates negatively with the migration energy. Thus, in general, larger solutes have lower migration energies and diffuse faster than smaller solutes in fcc cobalt. However, extremely large solutes, e.g., group III elements Sc, Y, Lu, tend to be trapped in an energy valley located halfway toward the vacancy, and monovacancy mediated diffusion may no longer be valid in such cases. Finally, for all the solutes considered, we systematically tabulate the diffusion-related quantities calculated—diffusion prefactors, migration and activation energies—constructing an extensive and accurate first-principles database for solute diffusion in fcc cobalt.

36 MATERIALS SCIENCE↗

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↗

On the Diffusivity of Moist Static Energy and Implications for the Polar Amplification Response to Climate Warming

Energy balance models (EBMs) have been widely used in a range of climate problems, but the assumption of constant diffusivity in the parameterization of the moist static energy (MSE) flux can be hardly justified. We demonstrate in this study that the diffusive MSE flux can be derived from the basic energy balance equation with a few tolerable assumptions. The estimated diffusivity is both spatially and seasonally dependent, and its midlatitude average is then tested against several scaling theories for the midlatitude eddy diffusivity. The result supports the diffusivity theory of Held and Larichev (1996) modified for the moist atmosphere, affording a dynamics-based parameterization of MSE diffusivity. The implementation of the parameterization in an EBM leads to an interactive MSE diffusivity that accounts for the midlatitude eddy response to climate forcing perturbations. Under a uniform radiative forcing, the EBM with a diffusivity so parameterized produces a weakening of the midlatitude diffusivity and a modestly polar-amplified surface temperature response as an inevitable outcome under the dual constraints of the nonlinear Clausius-Clapeyron relation and the temperature gradient-dependent diffusivity, even in the absence of any poleward amplifying radiative feedbacks. As the consequence of more isothermal temperature and reduced diffusivity, the variance of the midlatitude surface temperature also decreases with warming.

54 ENVIRONMENTAL SCIENCES↗

Effects of Moisture on Diffusion in Unmodified Wood Cell Walls: A Phenomenological Polymer Science Approach

Despite the importance of cell wall diffusion to nearly all aspects of wood utilization, diffusion mechanisms and the detailed effects of moisture remain poorly understood. In this perspective, we introduce and employ approaches established in polymer science to develop a phenomenological framework for understanding the effects of moisture on diffusion in unmodified wood cell walls. The premise for applying this polymer-science-based approach to wood is that wood polymers (cellulose, hemicelluloses, and lignin) behave like typical solid polymers. Therefore, the movement of chemicals through wood cell walls is a diffusion process through a solid polymer, which is in contrast to previous assertions that transport of some chemicals occurs via aqueous pathways in the cell wall layers. Diffusion in polymers depends on the interrelations between free volume in the polymer matrix, molecular motions of the polymer, diffusant dimensions, and solubility of the diffusant in the polymer matrix. Because diffusion strongly depends on whether a polymer is in a rigid glassy state or soft rubbery state, it is important to understand glass transitions in the amorphous wood polymers. Through a review and analysis of available literature, we conclude that in wood both lignin and the amorphous polysaccharides very likely have glass transitions. After developing and presenting this polymer-science-based perspective of diffusion through unmodified wood cell walls, suggested directions for future research are discussed. A key consideration is that a large difference between diffusion through wood polymers and typical polymers is the high swelling pressures that can develop in unmodified wood cell walls. This pressure likely arises from the hierarchical structure of wood and should be taken into consideration in the development of predictive models for diffusion in unmodified wood cell walls.

59 BASIC BIOLOGICAL SCIENCES↗

Influence of strong Coulomb coupling on diffusion in atmospheric pressure plasmas

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion–ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions (x i ≲ 10 —6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions (10 —6 ≲ x i ≲ 10 —2 ) there is a transition from ion-neutral to ion–ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion–ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion–ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions (x i ≳ 10 —2 ), strongly correlated ion–ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. Importantly, the results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

fast neutral gas heating↗

Measurement of moisture-dependent ion diffusion constants in wood cell wall layers using time-lapse micro X-ray fluorescence microscopy

Our future bioeconomy depends on increased utilization of renewable lignocellulosic biomass. Controlling the diffusion of chemicals, such as inorganic ions, within secondary plant cell walls is central to many biomass applications. However, insufficient understanding of intra-cell-wall diffusion within secondary plant cell walls is hindering the advancement of many lignocellulosic biomass applications. In this work, X-ray fluorescence microscopy was used to measure diffusion constants of K + , Cu 2+ , and Cl - diffusing through loblolly pine (Pinus taeda) cell wall layers under 70%, 75%, or 80% relative humidity (RH). Results revealed that diffusion constants increased with RH, the larger Cu 2+ diffused more slowly than the K + , and the Cl - diffusion constant was the same as that for the counter cation, indicating cations and anions diffused together to maintain charge neutrality. Comparison with electrical conductivity measurements showed that conductivity is being controlled by ion mobility over these RH. The results further support that intra-cell-wall diffusion of inorganic ions is a Fickian diffusion process occurring through rubbery amorphous polysaccharides, which contradicts previous assertions that intra-cell-wall diffusion is an aqueous process occurring through water pathways. Researchers can now utilize polymer science approaches to engineer the molecular architecture of lignocellulosic biomass to optimize properties for specific end uses.

59 BASIC BIOLOGICAL SCIENCES↗

Diffusion and mobility of atomic particles in a liquid

The diffusion coefficient of a test atom or molecule in a liquid is determined for the mechanism where the displacement of the test molecule results from the vibrations and motion of liquid molecules surrounding the test molecule and of the test particle itself. This leads to a random change in the coordinate of the test molecule, which eventually results in the diffusion motion of the test particle in space. Two models parameters of interaction of a particle and a liquid are used to find the activation energy of the diffusion process under consideration: the gas-kinetic cross section for scattering of test molecules in the parent gas and the Wigner–Seitz radius for test molecules. In the context of this approach, we have calculated the diffusion coefficient of atoms and molecules in water, where based on experimental data, we have constructed the dependence of the activation energy for the diffusion of test molecules in water on the interaction parameter and the temperature dependence for diffusion coefficient of atoms or molecules in water within the models considered. The statistically averaged difference of the activation energies for the diffusion coefficients of different test molecules in water that we have calculated based on each of the presented models does not exceed 10% of the diffusion coefficient itself. We have considered the diffusion of clusters in water and present the dependence of the diffusion coefficient on the cluster size. The accuracy of the presented formulas for the diffusion coefficient of atomic particles in water is estimated to be 50%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preparation of Diffusion Couples for Irradiation and High Temperature Testing of Representative TRISO PyC/SiC

Understanding the release of fission products from intact tristructural-isotropic (TRISO) coated particle fuel during normal and off-normal service is necessary to ensuring safe and efficient reactor operation. Planar diffusion couples have been fabricated to determine the diffusion kinetics of select fission product species (Ag, Ag+Pd, Eu, and Sr) in pyrocarbon (PyC) and silicon carbide (SiC) materials which are representative of TRISO fuel. Obtaining accurate diffusion coefficients for these systems in representative materials will allow for higher fidelity inputs for fuel performance models. Three different diffusion couple variants have been produced which include a Baseline, PyC, and SiC variant. The Baseline variant serves as an example of representative PyC and SiC layers as the layer properties mimic those found in TRISO fuel from the second Advanced Gas Reactor Fuel Qualification and Development Program. The PyC and SiC variants include variation in the respective layer properties to understand the impact of PyC density and SiC grain size on fission product diffusion. The diffusion couples were designed to allow for exploration of neutron irradiation on diffusion reflecting in reactor behavior as well as high temperature thermal diffusion reflecting diffusive transport during off-normal conditions. The methodology for selecting the diffusion couple design, refinement of the individual variant layer properties, and the complete process for diffusion couple construction are summarized in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Tritium Diffusion in Fe-Al Aluminide Coating Bulk Phases

The theoretical investigations of interstitial tritium diffusion pathways in various Al-rich (>60% at. Al) iron aluminide coating phases found that in FeNiAl 5 tritium prefers to diffuse in between layers, while in Fe 4 Al 13 , tritium prefers to diffuse across layers. In the case of Fe 2 Al x phases, the simulations show that the variable Al content in the channels affects the diffusion behavior of tritium in the channels. For 100% of vacancy, tritium prefers to diffuse at the center of the channel, while for 0% Al vacancy, tritium prefers to diffuse along the channel’s wall. In the case where 50% Al vacancy are present in the channels, tritium can diffuse in a mixed way, with positions oscillating between the channel’s center and wall. However, the large energy barrier involved in this case suggests that tritium diffusion could be correlated with Al diffusion in the channels. While we found that at 600 K interstitial tritium diffuses faster in Fe 2 Al x phases (D Τ ≤10 -11 m 2 .s -1 ), then in Fe 4 Al 13 (D Τ ≈10 -12 m 2 .s -1 ), and finally in FeNiAl 5 (D Τ ≈10 -13 m 2 .s -1 ), we also find that interstitial tritium generally diffuses faster in Fe-Al coating phases than in γ-LiAlO 2 pellets (D Τ ≈10 -14 m 2 .s -1 ).

36 MATERIALS SCIENCE↗

Tracer versus trace element diffusion - Diffusional decoupling of Sr concentration from Sr isotope composition

Chemical diffusivities for nonalkali major, minor, and selected trace elements have been measured at 10 kbar, 1300 C during interdiffusion of dacitic and rhyolitic melts. Sr tracer diffusivities were measured at the same conditions in melts of dacitic and rhyolitic composition and in a diffusion couple of these two melts undergoing simultaneous chemical diffusion. Chemical diffusivities of all elements studied are functions of the silica content of the melt, but at any one silica content the chemical diffusivities of all elements studied are within a factor of five of the Si and Al diffusivities. Si diffusivity varies from 2.8 x 10 to the -9th sq cm/s at 65 percent SiO2 to 5.1 x 10 to the -10th sq cm/s at 75 percent SiO2. The diffusivities of minor and trace elements are proposed to be controlled by the ability of the local melt structure to accommodate these elements while maintaining LTE. The measured Sr tracer diffusivities are an order of magnitude above Sr chemical diffusivities.

Baker, Don R.↗

Static and Dynamic Effects of Lateral Carrier Diffusion in Semiconductor Lasers

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. It is well known that the value of diffusion coefficients affects the threshold pumping current of a semiconductor laser. At the same time, the strength of carrier diffusion process is expected to affect the modulation bandwidth of an AC-modulated laser. It is important not only to investigate the combined DC and AC effects due to carrier diffusion, but also to separate the AC effects from that of the combined effects in order to provide design insights for high speed modulation. In this presentation, we apply a hydrodynamic model developed by the present authors recently from the semiconductor Bloch equations. The model allows microscopic calculation of the lateral carrier diffusion coefficient, which is a nonlinear function of the carrier density and plasma temperature. We first studied combined AC and DC effects of lateral carrier diffusion by studying the bandwidth dependence on diffusion coefficient at a given DC current under small signal modulation. The results show an increase of modulation bandwidth with decrease in the diffusion coefficient. We simultaneously studied the effects of nonlinearity in the diffusion coefficient. To clearly identify how much of the bandwidth increase is a result of decrease in the threshold pumping current for smaller diffusion coefficient, thus an effective increase of DC pumping, we study the bandwidth dependence on diffusion coefficient at a given relative pumping. A detailed comparison of the two cases will be presented.

Li, Jian-Zhong↗

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗

Passive Rocket Diffuser Testing: Reacting Flow Performance of Four Second-Throat Geometries

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure. As one of the nation's largest rocket testing facilities, the performance and design limitations of diffusers are of great interest to NASA's Stennis Space Center. This paper describes a series of tests conducted on four diffuser configurations to better understand the effects of inlet geometry and throat area on starting behavior and boundary layer separation. The diffusers were tested for a duration of five seconds with a 1455-pound thrust, LO2/GH2 thruster to ensure they each reached aerodynamic steady state. The effects of a water spray ring at the diffuser exits and a water-cooled deflector plate were also evaluated. Static pressure and temperature measurements were taken at multiple axial locations along the diffusers, and Computational Fluid Dynamics (CFD) simulations were used as a tool to aid in the interpretation of data. The hot combustion products were confirmed to enable the diffuser start condition with tighter second throats than predicted by historical cold-flow data or the theoretical normal shock method. Both aerodynamic performance and heat transfer were found to increase with smaller diffuser throats. Spray ring and deflector cooling water had negligible impacts on diffuser boundary layer separation. CFD was found to accurately capture diffuser shock structures and full-flowing diffuser wall pressures, and the qualitative behavior of heat transfer. However, the ability to predict boundary layer separated flows was not consistent.

Jones, Daniel R.↗

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗