Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion rate”

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 415 records · Page 23

Spectrum of the gamma-ray diffuse component observed from HEAO-1

The spectrum of the diffuse X and gamma ray background was measured between 15 keV and 4 MeV with the scintillation detectors aboard the HEAO 1 satellite. The apertures of the detectors were modulated on time scales of hours and the difference in counting rates measured the diffuse component flux. The observed spectrum is presented and compared with other measurements. At least two components are indicated, one below -100 keV and the other above. Possible origins are discussed.

Gruber, D. E.↗

Cooling rate estimations based on kinetic modelling of Fe-Mg diffusion in olivine

A finite one-dimensional kinetic model was developed to estimate the cooling rates of lunar rocks. The model takes into consideration the compositional zonation of olivine and applies Buening and Buseck (1973) data on ion diffusion in olivine. Since the 'as-solidified' profile of a given olivine is not known, a step-function, with infinite gradient, is assumed; the position of this step is based on mass balance considerations of the measured compositional profile. A minimum cooling rate would be associated with the preservation of a given gradient. The linear cooling rates of lunar rocks 12002 and 15555 were estimated by use of the olivine cooling-rate indicator to be 10 C/day and 5 C/day, respectively. These values are lower than those obtained by dynamic crystallization studies (10-20 C/day).

Taylor, L. A.↗

Mechanistic nuclear fuel performance modeling of uranium nitride

Uranium mononitride (UN) is a nuclear fuel candidate for advanced reactor designs and an alternative being considered for light water reactors due to its higher thermal conductivity and uranium density than UO 2 . As with any nuclear fuel, swelling and fission gas release are important factors for safety, while also being some of the hardest phenomena to predict with a high degree of confidence. Getting a grasp on the gas swelling behavior and release is crucial to lower the barrier for UN utilization. An accelerated swelling rate at high temperatures observed experimentally, sometimes referred to as “breakaway swelling,” further complicates the prediction of fuel performance of UN. A mechanistic model has been developed using a multiscale approach to describe the intragranular and intergranular fission gas behavior. Lower-length-scale calculations have been employed to inform models of the gas and self-diffusion behavior, resolution rate, and bubble shape. Leveraging previous work on high burnup UO 2 , two populations of intragranular bubbles are considered; small bulk bubbles and larger bubbles located along dislocations. The dislocation bubbles were found to be crucial to the overall swelling behavior, and the breakaway swelling transition was associated with the transition in the gas atom diffusion mechanism from an irradiation-induced athermal diffusion regime at lower temperatures to an intrinsic thermal equilibrium regime at higher temperatures, accelerating the growth of the dislocation bubbles. Similarly, the threshold for fission gas release was associated with the grain boundary vacancy diffusivity surpassing the gas atom diffusivity at sufficiently high temperatures, allowing the over-pressurized grain boundary bubble to grow in size and interconnect. Using thermo-mechanical models with the fission gas model, two integral fuel pin assessment cases were simulated. Finally, this work demonstrates the ability of a multiscale approach to accelerate the understanding of advanced fuel forms when experimental data is limited.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Creep deformation and cavitation in an additively manufactured Al-8.6Cu-0.4Mn-0.9Zr (wt%) alloy

Creep deformation and cavitation were investigated at 300 ºC in both tension and compression for an additively manufactured Al-8.6Cu-0.5Mn-0.9Zr (wt%) alloy in the as-fabricated state and after various aging treatments (aging at 300 °C/200 h or 350 °C/24 h and overaging at 400 °C/200 h). Creep mechanisms at 300 °C were determined by relating the measured creep response to corresponding microstructural and X-ray computed tomography observations. In compression, alloys in the as-fabricated and two aging conditions exhibited similarly high creep resistance. Overaging (400 °C/200 h) led to substantial coarsening of intragranular θ-Al2Cu precipitates and an expected drop in their Orowan strengthening contribution. In tension, minimum strain rates comparable to those in compression were obtained at any given stress; however, upon accumulation of some plastic strain in the matrix, creep cavities started to form, leading to accelerated tertiary stage creep deformation and rupture. Cavitation occurred exclusively along melt pool boundaries due to locally enhanced diffusion enabled by (i) large grain-boundary area in adjacent fine-grained zones and (ii) localization of creep strain in nearby heat-affected zones. Although cavity growth was initially diffusion-controlled, its rate was determined by matrix creep rate, consistent with constrained cavity growth mechanisms. This study reveals how microstructural complexities induced by the additive manufacturing process affect the creep and cavitation behavior of Al-Cu-Mn-Zr alloys. The underlying creep and cavitation mechanisms uncovered in this study point to pathways that improve the high-temperature properties of additively manufactured alloys.

36 MATERIALS SCIENCE↗

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗

Accelerated molecular dynamics simulations of dislocation climb in nickel

We report the mechanical behavior of materials operating under high temperatures is strongly influenced by creep mechanisms such as dislocation climb, which is controlled by the diffusion of vacancies. However, atomistic simulations of these mechanisms have traditionally been impractical due to the long time scales required. To overcome these time scale challenges, we use Parallel Trajectory Splicing (ParSplice), an accelerated molecular dynamics method, to simulate dislocation climb in nickel. We focus on modeling the activity of a vacancy near a jog on an edge dislocation in order to observe vacancy pipe diffusion and vacancy absorption at the jog. From rigorously constructed trajectories encompassing more than 2000 vacancy absorption events over a simulation time of more than 4μs at 900 K, a comprehensive sampling of available atomistic mechanisms is collated and analyzed further with molecular statics calculations. We estimate average rates for pipe diffusion and vacancy absorption into the jog using data from the dynamic and static calculations, finding very good agreement. Our results strongly suggest that the dominant mechanism for vacancy absorption by jogs is via biased diffusion to the dislocation core followed by fast pipe diffusion to the jog.

36 MATERIALS SCIENCE↗

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The interaction of atomic oxygen with thin copper films

A source of thermal, ground-state atomic oxygen has been used to expose thin copper films at a flux of 1.4 x 10 exp 17 atoms/sq cm s for times up to 50 min for each of five temperatures between 140 and 200 C. Rutherford backscattering spectroscopy was used to characterize the oxide formed during exposure. The observations are consistent with the oxide phase Cu2O. The time dependence and the temperature dependence of the oxide layer thickness can be described using oxide film growth theory based on rate limitation by diffusion. Within the time and temperature ranges of this study, the growth of the oxide layers is well described by the equation L(T,t) = 3.6 x 10 to 8th exp(- 1.1/2k sub B T)t exp 1/2, where L,T, and t are measured in angstroms, degrees Kelvin, and minutes, respectively. The deduced activation energy is 1.10 +/- 0.15 eV, with the attendant oxidation rate being greater than that for the corresponding reaction in molecular oxygen.

Gibson, B. C.↗

Mass transfer coefficients, kL, and air-CO2 ingassing rates in 3.4 m2 and 1-acre raceway ponds.

The US DOE 2016 Billion Ton Update (Langholtz 2016) projected that CO2 in flue-gas from power plants or similar sources would limit U.S. algal biofuel potential to under 5 billion gallons gasoline equivalent (gge) per year, due to lack of sufficient land, water and other requirements near such flue gas sources. 2nd generation carbon capture technologies are proposed to overcome this constraint (Davis 2017), and could expand this resource potential nearly 10-fold, assuming CO2 costs near the future, 2025 DOE NETL, goal of $40/metric ton CO2. An alternative pathway is the direct uptake of air CO2 into the algal ponds through an increase in the air-CO2 transfer rate resulting from the reaction of CO2 with hydroxide ions. However, such a ‘chemical enhancement’ in mass transfer depends on a high pH in the culture, and, by extension, the ability of microalgae to maintain high rates of carbon assimilation under such conditions. Abiotic experiments characterized the air-CO2 mass-transfer rate in 1-acre and 3.4 m2 raceway ponds as a function of pH, with rates approaching 10 g C/m2-day at pH 12 in the brackish water tested, equivalent to a biomass productivity potential of near 20 g AFDW/m2-day (@ 0.47 g C/g AFDW). The mass transfer rate was found to be nearly independent of the mass-transfer coefficient, indicating that at turbulence levels achievable in typical raceway ponds, the system is reaction rate, rather than diffusion limited. In biotic trials with an unknown green microalga, biomass productivity in ponds receiving CO2 from only direct surface air-CO2 exchange averaged 5 +/- 0.5 g AFDW/m2-day in August - September central California conditions, about half-that of experimental controls fertilized with supplemental CO2. pH in the air-only treatment reached a maximum of 10.5, supporting a model predicted air-CO2 ingassing rate between 2-3 g C/m2-day, consistent with the observed biomass carbon assimilation rate. Results suggest that unique alkaliphilic strains are required if reliant on direct in-pond air-CO2 transfer for algal inorganic carbon supply, and that reaching an economically viable biomass productivity will require strains that thrive at pH 11 and above.

09 BIOMASS FUELS↗

Physical electronics and surface physics

Current analytical and experimental investigations of the influence of surface processes on the rate of permeation of gases through solids are summarized. A brief description is presented of an analytical model that provides an approximate prediction of the dependence of the permeation rate on the rates of surface processes relative to the rate of the diffusion process. Results of experiments indicate that the desorption process of hydrogen from various materials is influenced by surface impurities to the degree that substantial changes in the spatial distribution of the desorbed molecules are observed when the impurities are varied. Except possibly for H-Nb, the permeation measurements indicate that the desorption process is not influenced to the degree that detectable changes in the permeation rates are observed when the impurities and temperatures are varied. On the basis of analysis, it is unlikely that the desorption process will be the rate-controlling step under the conditions of the present experiments.

Stickney, R. E.↗

PySIDT: Subgraph Isomorphic Decision Trees for Molecular Property Prediction

Accurate molecular property prediction is important across all fields of chemistry. Deep neural networks (DNNs) have become increasingly popular due to their ability to train automatically, avoiding the incredibly tedious process of constructing and extending traditional property estimation schemes. However, DNNs require large amounts of training data, are challenging to interpret, require large amounts of memory to load even during inference, and have severe difficulties incorporating qualitative chemical knowledge, which are often desired for molecular property prediction tasks. Here, in this study, we present PySIDT (https://github.com/zadorlab/PySIDT), a software for training and running inference on Subgraph Isomorphic Decision Trees (SIDTs). SIDTs are graph-based decision trees made of nodes associated with molecular substructures. Inference is done by descending target molecular structures down the decision tree to nodes with matching subgraph isomorphic substructures and making predictions based on the final (most specific) nodes matched. SIDTs scale down well to dataset sizes much smaller than is feasible for DNNs. As trees of molecular substructures, SIDTs are inherently readable and easy to visualize, making them easy to analyze. They are also straightforward to extend and retrain, facilitate uncertainty estimation, and enable easy integration of expert knowledge. We demonstrate the SIDT approach discussing its application to a diverse range of molecular prediction tasks: rate coefficient estimation, diffusion coefficient estimation, thermochemistry estimation, transition state bond stretch prediction, p K a prediction, stability of molecular structures, stability of surface structures, and prediction of surface lateral interaction energetics. Additionally, we demonstrate the power of the SIDT algorithms in two direct learning curve vanilla comparisons with the popular DNN-based software Chemprop and the popular gradient boosted trees-based software XGBoost on enthalpy of formation and rate coefficient prediction tasks. In particular, in the enthalpy of formation case, vanilla PySIDT is able to outperform vanilla Chemprop and XGBoost across the full range of training/validation set sizes out to 11,560 data points.

Johnson, Matthew Sean [Sandia National Laboratorie↗

Plasma transport in the Io torus - The importance of microscopic diffusion

This paper considers the question of whether the distribution of mass in the Io plasma torus is consistent with the concept of interchange eddy transport. Specifically, the flux tube content exhibits a gradual decrease with increasing radial distance from the source near Io without any evidence for substantial density irregularity associated with the plasma source or loss. Using a simple one-dimensional numerical model to simulate macroscopic interchange eddy transport, it is demonstrated that this smooth equilibrium distribution of mass can occur but only with the inclusion of a minimal level of small scale microscopic mixing at a rate approaching Bohm diffusion. Otherwise, the system exhibits a chaotic appearance which never approaches an equilibrium distribution. Various physical mechanisms for the microscopic diffusion process which is required to provide a sufficiently rapid mixing of material between the macroscopic eddies are discussed.

Mei, YI↗

The application of far IR lasers to the investigation of the optical properties of materials

The effect of diffusion on the saturation and absorption properties of molecular gases is considered. The analysis starts with a simple derivation of the multi-level rate equations. Incorporating the diffusion term, calculations are carried out both for the homogeneous and inhomogeneous broadening cases. Both graphical and analytical solutions are obtained for the absorption coefficient and the saturation intensity. It is shown that the effect of diffusion is characterized by a diffusion function, which varies inversely with the radius of the beam at low pressures and small beam sizes. The analysis is done in terms of two parameters k1a and k2a, which are simple functions of the intensity of the beam I, the effective relaxation rate w, and the effective diffusion constant D the a is designated as the radius of the beam.

Stafsudd, O. M.↗

Optical experiments on thermophoretically augmented submicron particle deposition from 'dusty' high temperature gas flows

A real-time laser reflectivity method and Pt ribbon targets are used to obtain experimental data on the initial deposition rate of MgO(s) particles of approximately 700 nm diameter from otherwise clean combustion products as a function of target temperature (about 950-1450 K) and mainstream gas temperature (about 1500-1600 K). These preliminary data are used to demonstrate the dominant role of thermophoresis (particle drift down a temperature gradient) and to assess the utility of recently developed theoretical methods for predicting and correlating the temperature dependence of thermophoretically augmented convective-diffusion 'dust' deposition rates from flowing hot gases.

Rosner, Daniel E.↗

Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions

Hydrogen isotopes retention and desorption during and after discharges in fusion devices are still not well understood due to the complex device conditions and limitations of in-situ diagnostics and measurements. Here, we simulated well-diagnosed recent experiments at the DIII-D facility to benchmark our ITMC-DYN integrated package of modeling deuterium diffusion, retention, and desorption during and after D discharge irradiation. Modeling results were compared with detail experimental data of D desorption fluxes for various irradiation conditions. We predicted the temporal evolution of free and trapped D distribution in tungsten (W) plasma-facing material (PFM). Effects of key parameters namely diffusion coefficient, recombination rate, trapping energies against different defect types, were examined in these simulations. Existing experimental data of these parameters in literature varies significantly which makes it harder to identify key mechanisms and physics responsible for hydrogen isotope retention and desorption. The purpose of this work is to accurately simulate recent well-diagnosed reactor experiments given the uncertainties in such parameters and identify mechanisms responsible for the retention and desorption. We implemented the best identified diffusion, recombination, and trapping parameters in ITMC-DYN package that integrate both various collisional and thermal processes. We predicted, for example, that sample cooling between discharges in DIII-D operations can significantly affect the spatial distribution of trapped D in W under reactor irradiation conditions. Correct prediction of desorption spectra from samples irradiated during 10 DIII-D discharges showed that up to 35% of D can be retained in high binding energy defects such as vacancy clusters or voids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗