Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective 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 145 records · Page 8

Mechanistic FCCI model for BISON

In this report, multi-scale simulations to develop a mechanistic model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY23 are described. The model assumes that cladding wastage is dominated by the intermetallic phase (Fe,Cr)$_{17}$Nd$_2$, and calculates production and diffusion of Nd in the fuel, diffusion into the cladding, and formation and growth of the intermetallic wastage layer. Atomistic simulations are used to calculate diffusivity of Nd in the (Fe,Cr) system. Mesoscale simulations, informed by previous atomistic calculations, are used to determine the effective diffusivity of Nd through the fuel, accounting for empty and sodium-filled porosity within the fuel. The BISON model incorporates this data and calculates the wastage layer thickness as a function of time, assuming parabolic growth kinetics based on the best available data. The results of the new mechanistic BISON model are compared to existing assessment cases for the X447 experiment in EBR-II, and show good agreement with experimental data and the existing BISON empirical model. To improve predictions of kinetics for future work, a phase-field model that includes fuel, wastage, and cladding phases is developed. Preliminary testing indicates that the kinetics may differ significantly from the assumed parabolic growth law, suggesting direction for future improvement of the model that will allow it to be applied to other fuel designs with greater confidence.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A molecular dynamics study of nuclear quantum effect on the diffusion of hydrogen in condensed phase

In this paper, the quantum effect of hydrogen molecule on its diffusivity is analyzed using Molecular Dynamics (MD) method. The path integral centroid MD (CMD) method is applied for the reproduction method of time evolution of the molecules. The diffusion coefficient of liquid hydrogen is calculated using the Green-Kubo method. The simulation is performed at wide temperature region and the temperature dependence of the quantum effect of hydrogen molecule is addressed. The calculation results are compared with those of classical MD results. As a result, it is confirmed that the diffusivity of hydrogen molecule is changed depending on temperature by the quantum effect. It is clarified that this result can be explained that the dominant factor by quantum effect on the diffusivity of hydrogen changes from the swollening the potential to the shallowing the potential well around 30 K. Moreover, it is found that this tendency is related to the temperature dependency of the ratio of the quantum kinetic energy and classical kinetic energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian Analysis of TRISO Fuel: Quantifying Model Inadequacy, Incorporating Lower-Length-Scale Effects, and Developing Parallel Active Learning Capabilities

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Comparative study of structure formation and mechanical behavior of age-hardened Ti–Nb–Zr and Ti–Nb–Ta shape memory alloys

This work sets out to study the peculiar effects of aging treatment on the structure and mechanical behavior of cold-rolled and annealed biomedical Ti–21.8Nb–6.0Zr (TNZ) and Ti–19.7Nb–5.8Ta (TNT) (at.%) shape memory alloys by means of transmission electron microscopy, X-ray diffractometry, functional fatigue and thermomechanical testing techniques. Dissimilar effects of aging treatment on the mechanical behavior of Zr- and Ta-doped alloys are explained by the differences in the ω-phase formation rate, precipitate size, fraction and distribution, and by their effect on the alloys' critical stresses and transformation temperatures. Even short-time aging of the TNZ alloy leads to its drastic embrittlement caused by “overaging”. On the contrary, during aging of the TNT alloy, formation of finely dispersed ω-phase precipitates is gradual and controllable, which makes it possible to finely adjust the TNT alloy functional properties using precipitation hardening mechanisms. To create in this alloy nanosubgrained dislocation substructure containing highly-dispersed coherent nanosized ω-phase precipitates, the following optimum thermomechanical treatment is recommended: cold rolling (true strain 0.37), followed by post-deformation annealing (600 °C, 15–30 min) and age-hardening (300 °C, 30 min) thermal treatments. It is shown that in TNT alloy, pre-transition diffraction effects (diffuse reflections) can “mask” the β-phase substructure and morphology of secondary phases. - Highlights: • TNZ alloy is characterized by much higher ω-phase precipitation rate than TNT alloy. • Difference in precipitation rates is linked to the difference in Zr and Ta diffusion mobility. • Aging of nanosubgrained TNZ alloy worsens its properties irrespective of the aging time. • Aging time of nanosubgrained TNT alloy can be optimized to improve its properties.

36 MATERIALS SCIENCE↗

Hydrogen effects on thermal diffusivity and electrical resistivity of zircaloy cladding

Thermal diffusivity measurements on Zirconium-based cladding materials have historically been a challenge due to the difficulty to measure on specimens with curved geometries, including nuclear grade Zircaloy cladding materials. Here, in this work, we first used laser flash analysis method and four-probe configuration method to measure the thermal diffusivity and electrical resistivity of Zircaloy tubes respectively, which show good agreement with Zircaloy plates in this work, as well as previously published data. The consistent results proved the applicability of the laser flash analysis setup and four-probe configuration method for investigating thermal diffusivity and electrical resistivity of Zircaloy tubes. We further investigated the hydrogen effect on thermal diffusivity and electrical resistivity of Zircaloy. Hydrogen plays significant roles in the thermal diffusivity of Zircaloy, which depends on the hydrogen concentration. For higher hydrogen contents (1130 and 1820 wppm in this work), where phase transformation (α-Zr + δ-hydride → α-Zr + β-Zr) occurs at 567 °C, thermal diffusivity decreases as a function of temperature at the α-Zr + δ-hydride phase regime, while increase at higher temperature at the α-Zr + β-Zr phase regime. For low hydrogen concentration, where hydride dissolved into α-Zr matrix phase at higher temperature, the thermal diffusivity is lower than non-hydrided Zircaloy-4 with similar temperature-dependent trends. Such observation demonstrates the hydrogen effects on reducing thermal diffusivity of α-Zr phase. Hydrogen increases the electrical resistivity of Zircaloy. Similar to thermal diffusivity results, the phase transformation causes a reversal of temperature-dependent trends in electrical resistivity results in Zircaloy-4 with higher hydrogen contents.

36 MATERIALS SCIENCE↗

Interaction of Iron with Alumina Refractory Under Flash Ironmaking Conditions

A novel flash ironmaking technology (FIT) based on the direct reduction of iron ore concentrate by a reducing gas such as hydrogen, natural gas, or coal gas has been developed at the University of Utah. A proper choice of refractories is expected to play an important role in the development of FIT and its scale-up. In this work, the interaction of iron with alumina refractory under flash ironmaking conditions has been investigated. A thermodynamic basis for the interaction phenomena has been developed by considering the Fe-Al-O system, and a kinetic model based on solid-state diffusion is proposed to describe the growth of the hercynite (FeAl{sub 2}O{sub 4}) spinel formed as a result of the interaction between Fe, O (from H{sub 2}O and CO{sub 2}), and Al{sub 2}O{sub 3}. Experiments were conducted with Fe powders and pure Al{sub 2}O{sub 3} in the temperature range 1200 to 1500 °C under gas atmospheres relevant to flash ironmaking, and the reacted samples were analyzed using XRD, SEM-EDX, and EPMA techniques. The analyses of reacted samples confirmed the formation of the hercynite (FeAl{sub 2}O{sub 4}) spinel phase and the thickness of the FeAl{sub 2}O{sub 4} spinel as measured using elemental line scans in SEM-EDX and EPMA varied in the range of 0.5 to 1.3 mm. Results showed that the proposed model suitably describes the growth of the hercynite spinel layer, which obeyed the parabolic rate law at all temperatures. The parabolic rate constants were obtained and the diffusion-controlled growth of the FeAl{sub 2}O{sub 4} layer strongly depended on temperature. Furthermore, using the kinetic model, an expression for the effective diffusivity ( (D{sub eff})-bar ) was obtained and its values at the experimental temperatures were determined. The solid-state diffusion was an activated process, as expected, with an activation energy value of 231 kJ/mol.

36 MATERIALS SCIENCE↗

Analysis based on global model of nitrogen plasma produced by pulsed microwave at low pressure

This paper analyzes certain evolution processes in nitrogen plasmas discharged using pulsed microwaves at low pressure. Comparing the results obtained from the global model incorporating diffusion and the microwave transmission method, the temporal variation of the electron density is analyzed. With a discharge pressure of 300 Pa, the results obtained from experiments and the global model calculation show that when the discharge begins the electron density in the plasma rises quickly, to a level above the critical density corresponding to the discharge microwave frequency, but falls slowly when the discharge microwave pulse is turned off. The results from the global model also show that the electron temperature increases rapidly to a peak, then decays after the electron density reaches the critical density, and finally decreases quickly to room temperature when the discharge microwave pulse is turned off. In the global model, the electron density increases because the high electron temperature induces a high ionization rate. The decay of the electron density mainly comes from diffusion effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulating Dopant Diffusion in a Detalied Porous Structure

The National Ignition Facility (NIF) at Lawrence Livermore National Laboratory (LLNL) uses highpowered lasers to explore cutting edge fusion technology. Specialized optics, known as fiber lasers, generate and amplify the NIF’s high-powered lasers that facilitate fusion experiments. Fiber lasers are glass fibers that are doped with heavy elements to stimulate light emission. The fiber laser doping processing uses a predeposited layer of porous silica nanoparticles on the surface of the glass to absorb a dopant rich salt solution. Our project simulated the diffusion of dopants into porous silica with prescribed colloidal stacking geometries using COMSOL Multiphysics. With these simulations we identified the porosity and tortuosity of the porous nanoparticle structures as the primary parameters that will dictate how the dopants will diffusive into the porous surface of the glass. Porosity and Tortuosity are geometric parameters that can easily be input into more detailed simulations to estimate an effective diffusivity within the porous layer. Future work should account for surface reactions to accurately simulate the full-scale diffusion of the solution doping process.

36 MATERIALS SCIENCE↗

Dose rate effects in the radiation damage of the plastic scintillators of the CMS hadron endcap calorimeter

We present measurements of the reduction of light output by plastic scintillators irradiated in the CMS detector during the 8 TeV run of the Large Hadron Collider and show that they indicate a strong dose rate effect. The damage for a given dose is larger for lower dose rate exposures. The results agree with previous measurements of dose rate effects, but are stronger due to the very low dose rates probed. Here, we show that the scaling with dose rate is consistent with that expected from diffusion effects.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Adsorptive Transport of Noble Gas Tracers in Porous Media

The transport of noble gas radionuclides in porous media is relevant to the detection of underground nuclear detonations as well as the sequestration of reprocessing off-gases. However, in field tests releasing radioxenon underground, the quantity of radioxenon observed at the surface has fallen well below expectations1. This study examined the diffusivity of noble gases (Kr and Xe) and the inert molecular gas sulfur hexafluoride (SF6) in porous media to observe any unexpected behavior. To replicate the transport of radiogenic noble gases in underground media, a two-bulb gaseous diffusion apparatus was constructed. The two bulbs are connected with a column of 10-30 Ottawa sand and ordinary atmosphere filled both the bulbs and pore spaces. The tracer gases were diluted in an isolated bulb to approximately 1000 ppm. Once released, the gases were allowed to diffuse through the column. Aliquots were withdrawn at regular time intervals from both bulbs and concentrations were quantified using a Shimadzu QP2010 SE gas chromatograph-mass spectrometer. The effective diffusivity was then calculated using a maximum likelihood estimate on the quasi-steady state model.

Paul, Matthew R.↗

Hydrogen Diffusion Coefficient Measures on Thin Film Uranium Oxide

Thin films of uranium oxides, putatively UO 2 and U 3 O 8 , were deposited on palladium/silver (75/25) foil discs (10 μm thick, 10 mm diameter). Methodology for sealing these foils, applying a hydrogen pressure (~ 1 atm) to the oxide side, and measuring the pressure on the permeate side (opposite side of foil) is reported. The experimental apparatus is held at 100 °C to speed diffusion, aided by the pressure differential across the foil (~1 atm (750 torr) on oxide side, initially low vacuum on permeate side, ~1x 10 -4 torr). Early results indicate an effective diffusion coefficient of roughly 7.66 x 10 -17 +/- 2.22 x 10 -17 cm 2 /sec for UO 2 . These values are in line with expectations and prior measures relative to large lag times for system baseline (substrate only). Evaluation of the technique for U 3 O 8 thin films suggests further development will be needed to extend the technique to more oxidized films (U 3 O 8 , UO 3 ). Characterization of the foils (thickness and speciation) post diffusion experiments will be carried out by SIGMA (Eric Tegtmeier and Andy Richards) in FY24 which will sharpen the uncertainty quantification for UO 2 diffusion coefficients, aiding the nucleation model for DRACO.

36 MATERIALS SCIENCE↗

Processes of diffusion in a limited ensemble of charged particles in a static magnetic field

The dynamics of a limited ensemble of charged particles in a spatially uniform static magnetic field has been studied analytically and numerically. Calculations have been performed for particles with various masses and charges in a wide range of the parameters of the systems under study. It has been shown that the transverse diffusion coefficient in an ensemble of particles in the magnetic field for weakly nonideal systems is described by the Townsend formula. Analytical estimates have been proposed for the analysis of the effective diffusion radius of the ensemble of particles owing to their thermal motion. The proposed relations have been tested by the numerical simulation of the problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗

Density functional theory calculations of self- and Xe diffusion in U 3 Si 2

Uranium silicide, U 3 Si 2 , has been proposed as an advanced nuclear fuel to be used in light water reactors (LWRs). Development of this alternative to the predominant current fuel, UO 2 , is motivated by enhanced accident tolerance as a result of higher thermal conductivity as well as improved fuel cycle economics through increased uranium density. In order to accurately model the fuel performance of U 3 Si 2 , the diffusion rate of point defects, which is related to self-diffusion, and of fission gas atoms must be determined. DFT calculations are used to predict the U and Si point defect concentrations, the corresponding self-diffusivities, the preferred Xe trap site and the Xe diffusivity. Effects of irradiation are not considered. A low defect formation energy and a high entropy for Si interstitials give rise to Si-rich non-stoichiometry at elevated temperatures. Both U and Si self-diffusion and Xe diffusion are anisotropic as a consequence of the tetragonal crystal structure of U 3 Si 2 . Si diffusion occurs by interstitial mechanisms in both the a-b plane and along the c axis, while the U c axis diffusion rate is controlled by a vacancy mechanism. Interstitial diffusion of U is very fast in the a-b plane of the U 3 Si 2 crystal structure. Xe atoms prefer to occupy U vacancy trap sites. The highest Xe diffusion rate occurs by a vacancy mechanism in both the a-b plane and along the c axis. The diffusion rate is similar in the a-b plane and along the c axis. U and Si self-diffusion and Xe diffusion are all faster in U3Si2 than intrinsic U and Xe diffusion in conventional UO 2 nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mesoscopic analyses of the impact of morphology and operating conditions on the transport resistances in a proton-exchange-membrane fuel-cell catalyst layer

Exploring the origins of local transport resistance and characterizing the oxygen transport resistances in the catalyst layer (RCL) are critical for cost reduction. In this paper, a comprehensive mesoscopic model for simulating coupled transport processes of oxygen and water vapor for different structural parameters under different operating conditions in reconstructed microstructures is proposed. The local transport resistance is calculated after achieving the limiting current density and the effective diffusivity of oxygen. The results demonstrate that RCL increases greatly with decreasing platinum loading (LPt) and the transport resistances in other components of the cell dominate for high-loadings. Both the reduced oxygen permeation coefficient in the ionomer thin-film and the adsorption resistance account for the origins of local transport resistance. The local transport resistance increases with the bare carbon ratio for a constant LPt and Pt/C ratio due to the decreased effective ionomer surface, and increases with the I/C ratio due to the increased ionomer thickness and decreased Knudsen diffusivity. Due to the presence of liquid water, a slight decrease followed by an increase of the local transport resistance versus relative humidity is obtained. The contribution of ionomer thin-films to RCL is more sensitive to liquid saturation compared with that of pores. This journal is

Mu, YT↗

In Vivo Protein Dynamics on the Nanometer Length Scale and Nanosecond Time Scale

Selectively labeled GroEL protein was produced in living deuterated bacterial cells to enhance its neutron scattering signal above that of the intracellular milieu. Quasi-elastic neutron scattering shows that the in-cell diffusion coefficient of GroEL was (4.7 ± 0.3) × 10 –12 m 2 /s, a factor of 4 slower than its diffusion coefficient in buffer solution. Furthermore, for internal protein dynamics we see a relaxation time of (65 ± 6) ps, a factor of 2 slower compared to the protein in solution. Comparison to the literature suggests that the effective diffusivity of proteins depends on the length and time scale being probed. Retardation of in-cell diffusion compared to the buffer becomes more significant with the increasing probe length scale, suggesting that intracellular diffusion of biomolecules is nonuniform over the cellular volume. This approach outlined here enables investigation of protein dynamics within living cells to open up new lines of research using “in-cell neutron scattering” to study the dynamics of complex biomolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusion Blurring of GaAs Quantum Wells Grown at Low Temperature

The processes of the diffusion blurring of a periodic system of GaAs quantum wells separated by AlGaAs barriers are studied by photoluminescence spectroscopy. The system is grown by molecular-beam epitaxy at a low temperature (200°C) and additionally doped with Sb and P isovalent impurities. Postgrowth annealing at the temperature 750°C for 30 min induces an increase in the energy corresponding to the photoluminescence peak of the e1–hh1 exciton state in quantum wells because of blurring of the epitaxial GaAs/AlGaAs interfaces due to enhanced Al–Ga interdiffusion in the cation sublattice. For the Al concentration profile defined by linear diffusion into quantum wells, the Schrödinger equation for electrons and holes is solved. It is found that the experimentally observed energy position of the photoluminescence peak corresponds to the Al–Ga interdiffusion length 3.4 nm and to the effective diffusion coefficient 6.3 × 10{sup –17} cm{sup 2} s{sup –1} at the temperature 750°C. This value is found to be close to the corresponding value for GaAs quantum wells grown at low temperatures without additional doping with Sb and P impurities. From the results obtained in the study, it is possible to conclude that enhanced As–Sb and As–P interdiffusion in the anion sublattice only slightly influences the processes of Al–Ga interdiffusion in the cation sublattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗