Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion models”

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 1,135 records · Page 63

Evaluation of Technetium Getters to Improve the Performance of Cast Stone

Cast Stone has been selected as the preferred waste form for solidification of aqueous secondary liquid effluents from the Hanford Tank Waste Treatment and Immobilization Plant (WTP) process condensates and low-activity waste (LAW) melter off-gas caustic scrubber effluents. Cast Stone is also being evaluated as a supplemental immobilization technology to provide the necessary LAW treatment capacity to complete the Hanford tank waste cleanup mission in a timely and cost effective manner. One of the major radionuclides that Cast Stone has the potential to immobilize is technetium (Tc). The mechanism for immobilization is through the reduction of the highly mobile Tc(VII) species to the less mobile Tc(IV) species by the blast furnace slag (BFS) used in the Cast Stone formulation. Technetium immobilization through this method would be beneficial because Tc is one of the most difficult contaminants to address at the U.S. Department of Energy (DOE) Hanford Site due to its complex chemical behavior in tank waste, limited incorporation in mid- to high-temperature immobilization processes (vitrification, steam reformation, etc.), and high mobility in subsurface environments. In fact, the Tank Closure and Waste Management Environmental Impact Statement for the Hanford Site, Richland, Washington (TC&WM EIS) identifies technetium-99 ( 99 Tc) as one of the radioactive tank waste components contributing the most to the environmental impact associated with the cleanup of the Hanford Site. The TC&WM EIS, along with an earlier supplemental waste-form risk assessment, used a diffusion-limited release model to estimate the release of different contaminants from the WTP process waste forms. In both of these predictive modeling exercises, where effective diffusivities based on grout performance data available at the time, groundwater at the 100-m down-gradient well exceeded the allowable maximum permissible concentrations for 99 Tc. (900 pCi/L). Recent relatively short-term (63 day) leach tests conducted on both LAW and secondary waste Cast Stone monoliths indicated that 99 Tc diffusivities were at or near diffusivities where the groundwater at the 100-m down-gradient well would exceed the allowable maximum permissible 99 Tc concentrations. There is, therefore, a need and an opportunity to improve the retention of Tc in the Cast Stone waste form. One method to improve the performance of the Cast Stone waste form is through the addition of “getters” that selectively sequester Tc inside Cast Stone.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Characterization and Modeling of Grain Boundary Chemistry Evolution in Ferritic Steels under Irradiation

Ferritic/martensitic (FM) steels such as HT-9, T-91 and NF12 with chromium concentrations in the range of 9-12 at.% Cr and high Cr ferritic steels (oxide dispersion strengthened steels with 12-18% Cr) are receiving increasing attention for advanced nuclear applications, e.g. cladding and duct materials for sodium fast reactors, pressure vessels in Generation IV reactors and first wall structures in fusion reactors, thanks to their advantages over austenitic alloys. Predicting the behavior of these alloys under radiation is an essential step towards the use of these alloys. Several radiation-induced phenomena need to be taken into account, including phase separation, solute clustering, and radiation-induced segregation or depletion (RIS) to point defect sinks. RIS at grain boundaries has raised significant interest because of its role in irradiation assisted stress corrosion cracking (IASCC) and corrosion of structural materials. Numerous observations of RIS have been reported on austenitic stainless steels where it is generally found that Cr depletes at grain boundaries, consistently with Cr atoms being oversized in the fcc Fe matrix. While FM and ferritic steels are also subject to RIS at grain boundaries, unlike austenitic steels, the behavior of Cr is less clear with significant scatter and no clear dependency on irradiation condition or alloy type. In addition to the lack of conclusive experimental evidence regarding RIS in F-M alloys, there have been relatively few efforts at modeling RIS behavior in these alloys. The need for predictability of materials behavior and mitigation routes for IASCC requires elucidating the origin of the variable Cr behavior. A systematic detailed high-resolution structural and chemical characterization approach was applied to ion-implanted and neutron-irradiated model Fe-Cr alloys containing from 3 to 18 at.% Cr. Atom probe tomography analyses of the microstructures revealed slight Cr clustering and segregation to dislocations and grain boundaries in the ion-irradiated alloys. More significant segregation was observed in the neutron irradiated alloys. For the more concentrated alloys, irradiation did not affect existing Cr segregation to grain boundaries and segregation to dislocation loops was not observed perhaps due to a change in the dislocation loop structure with increasing Cr concentration. Precipitation of α’ was observed in the neutron irradiated alloys containing over 9 at.% Cr. However ion irradiation appears to suppress the precipitation process. Initial low dose ion irradiation experiments strongly suggest a cascade recoil effect. The systematic analysis of grain boundary orientation on Cr segregation was significantly challenged by carbon contamination during ion irradiation or by existing carbon and therefore carbide formation at grain boundaries (sensitization). The combination of the proposed systematic experimental approach with atomistic modeling of diffusion processes directly addresses the programmatic need for developing and benchmarking predictive models for material degradation taking into account atomistic kinetics parameters

36 MATERIALS SCIENCE↗

Transport, charge exchange and loss of energetic heavy ions in the earth's radiation belts - Applicability and limitations of theory

Computer simulations of processes which control the relative abundances of ions in the trapping regions of geospace are compared with observations from discriminating ion detectors. Energy losses due to Coulomb collisions between ions and exospheric neutrals are considered, along with charge exchange losses and internal charge exchanges. The time evolution of energetic ion fluxes of equatorially mirroring ions under radial diffusion is modelled to include geomagnetic and geoelectric fluctutations. Limits to the validity of diffusion transport theory are discussed, and the simulation is noted to contain provisions for six ionic charge states and the source effect on the radiation belt oxygen ion distributions. Comparisons are made with ion flux data gathered on Explorer 45 and ISEE-1 spacecraft and results indicate that internal charge exchanges cause the radiation belt ion charge state to be independent of source charge rate characteristics, and relative charge state distribution is independent of the radially diffusive transport rate below the charge state redistribution zone.

Spjeldvik, W. N.↗

Effect of double layers on magnetosphere-ionosphere coupling

The dynamic aspects of auroral current structures are reviewed with emphasis on consequences for models of microscopic turbulence (MT). A number of models of MT are introduced into a large-scale model of Alfven wave propagation to determine the effect of various models on the overall structure of auroral currents. The effect of a double layer (DL) electric field which scales with the plasma temperature and the Debye length is compared with the effect of anomalous resistivity due to electrostatic ion cyclotron turbulence in which the electric field scales with the magnetic field strength. It is shown that the DL model is less diffusive than the resistive model, indicating the possibility of narrow intense current structures.

Lysak, Robert L.↗

Ocean Turbulence: Two-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 3

In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Canuto, V. M.↗

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.

36 MATERIALS SCIENCE↗

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Geometry of Non-Diffusive Tracer Transport in Gridded Atmospheric Models

A first-order linear and numerically non-diffusive Eulerian transport algorithm, minVAR, was recently developed for preservation of correlations between interrelated tracers during advective transport. The present study extends this work by: (1) providing further investigation of several interesting geometric constructions found in contours of constant minVAR, short for minimum variance, through extension to three coordinate dimensions. These contours capture point-by-point representations of thousands of individual atmospheric aerosol and/or cloud particles as they evolve and are rendered on Eulerian grids at a level of sub-grid resolution limited only by numerical precision; and (2) exploration of geometric similarities between the Arakawa C-grid, used to obtain interpolated values of the wind field at grid scale and minVAR. In particular, we consider interpolation of the u and v horizontal components of wind velocity from grid to sub-grid scales. The last results are motivated by recent applications of the Weather Research and Forecasting (WRF) model applied in the coastal Houston region, where the recent TRacking Aerosol Convection Interactions ExpeRiment (TRACER) field campaign was organized. A unique and fully consistent mapping is obtained between particles moving along meteorological wind trajectories and the non-diffusive, non-dispersive representation of such trajectories on an Eulerian grid.

54 ENVIRONMENTAL SCIENCES↗

The modeling of gas phase permeation through iron and nickel membranes

The gas phase permeation of hydrogen through metal membranes encompasses many kinetic processes. This paper reviews a permeation model which incorporates second order gas-surface reaction kinetics with simple bulk diffusion. The model is used to investigate the effect of this particular surface reaction of steady-state permeation. The dependence of the steady-state permeation flux on temperature, pressure, and thickness of the membrane has been calculated. The model predicts that the bulk controlled steady-state flux will change to a surface limited steady-state flux as either the temperature or thickness of the membrane is reduced. Finally, using independently derived parameters, the model is compared with permeation measurements on iron and nickel membranes.

Kuhn, David K.↗

An engineering assessment of gas/surface interactions in free-molecular aerodynamics

Data available on the gas/surface interaction (GSI) phenomenon, particularly in the rarefied to near-free-molecular flow regime are collected, categorized, and analyzed with a purpose of developing a GSI model that could be used as a guide to spacecraft designers. The study shows that there are not enough useful data for building a high-confidence GSI model. However, sufficient results are obtained to suggest that continued reliance on the diffuse GSI model is inappropriate, particularly in the near-free-molecular regime.

Knox, E. C.↗

Communication: Modeling of concentration dependent water diffusivity in ionic solutions: Role of intermolecular charge transfer

The translational diffusivity of water in solutions of alkali halide salts depends on the identity of ions, exhibiting dramatically different behavior even in solutions of similar salts of NaCl and KCl. The water diffusion coefficient decreases as the salt concentration increases in NaCl. Yet, in KCl solution, it slightly increases and remains above bulk value as salt concentration increases. Previous classical molecular dynamics simulations have failed to describe this important behavior even when polarizable models were used. Here, we show that inclusion of dynamical charge transfer among water molecules produces results in a quantitative agreement with experiments. Our results indicate that the concentration-dependent diffusivity reflects the importance of many-body effects among the water molecules in aqueous ionic solutions. Comparison with quantum mechanical calculations shows that a heterogeneous and extended distribution of charges on water molecules around the ions due to ion-water and also water-water charge transfer plays a very important role in controlling water diffusivity. Explicit inclusion of the charge transfer allows us to model accurately the difference in the concentration-dependent water diffusivity between Na{sup +} and K{sup +} ions in simulations, and it is likely to impact modeling of a wide range of systems for medical and technological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusive equilibrium distributions of He(+) in the plasmasphere

Recent satellite observations of thermal ion composition in the near-equatorial plasmasphere have shown that He(+) comprises 5-10 percent typically and occasionally 25 percent or more of the total thermal ion density. A steady state diffusive equillibrium model for the distributions of H(+), He(+) and O(+) along a plasmaspheric flux tube is used to elicit effects that may help explain these observed high He(+) fractional concentrations. The model indicates that both the ionospheric composition and the temperature distribution along the flux tubes are important factors controlling the equatorial He(+) composition, through the phlasma scale height and thermal diffusion effects. Direct comparison of the model results with thermal ion observations by ISEE-1 indicates that the effects incorporated into the model may explain some of the elevated He(+) concentrations. In some instances, however, effects not included in the model may be of importance.

Waite, J. H., Jr.↗

Modeling of Particle Acceleration at Multiple Shocks via Diffusive Shock Acceleration: Preliminary Results

Successful forecasting of energetic particle events in space weather models require algorithms for correctly predicting the spectrum of ions accelerated from a background population of charged particles. We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box. We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (E(sub max)) appropriate for quasi-parallel and quasi-perpendicular shocks and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Parker, L. Neergaard↗

Modulation of Galactic Cosmic Rays in the Inner Heliosphere, Comparing with PAMELA Measurements

We develop a numerical model to study the time-dependent modulation of galactic cosmic rays in the inner heliosphere. In the model, a time-delayed modified Parker heliospheric magnetic field (HMF) and a new diffusion coefficient model, NLGCE-F, from Qin and Zhang, are adopted. In addition, the latitudinal dependence of magnetic turbulence magnitude is assumed to be ∼(1+sin{sup 2}θ)/2 from the observations of Ulysses, and the radial dependence is assumed to be ∼r{sup S}, where we choose an expression of S as a function of the heliospheric current sheet tilt angle. We show that the analytical expression used to describe the spatial variation of HMF turbulence magnitude agrees well with the Ulysses, Voyager 1, and Voyager 2 observations. By numerically calculating the modulation code, we get the proton energy spectra as a function of time during the recent solar minimum, it is shown that the modulation results are consistent with the Payload for Antimatter-Matter Exploration and Light-nuclei Astrophysics measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Atomistic models of In and Ga diffusion in Cu(In,Ga)Se 2

To improve the performance of Cu(In,Ga)Se 2 thin-film photovoltaic devices, a robust understanding of the dominant diffusion pathways of the alloy species In and Ga is needed. Here, the most probable defect complexes and mechanisms for In and Ga diffusion are identified with the aid of density functional theory. The binding energies and migration barriers for these complexes are calculated in bulk CuInSe 2 and CuGaSe 2 . Analytic models and kinetic lattice Monte Carlo simulations are employed to predict the diffusivity of In and Ga under variations in composition and temperature. Here, we find that a model based on coulombic interactions between group III antisites and vacancies on the Cu-sublattice produces results that match well with experiment.

14 SOLAR ENERGY↗

Secondary Waste Form Development and Data Package for Hanford's Integrated Disposal Facility Performance Assessment - 16479

The Hanford Low-Activity Waste (LAW) Vitrification facility is planned to immobilize Hanford tank waste for on-site, near-surface disposal at Hanford's Integrated Disposal Facility (IDF). During this process, liquid and solid secondary wastes will be generated which are also slated for disposal at the IDF. In order to receive authorization to dispose of the waste in the IDF, a Performance Assessment (PA) must be completed in accordance with DOE Order 435.1, Radioactive Waste Management, and its accompanying Manual. The PA must include calculations of potential doses to a member of the public for a 1,000-year period to provide a reasonable expectation that performance objectives will not be exceeded and that the waste will be managed in a manner that is protective of worker and public health and safety, and the environment. Liquid and solid secondary wastes are planned to be immobilized within cementitious waste forms similar to those used elsewhere in the United States and around the world for permanent disposal of low-level radioactive waste (LLW). Data packages and supporting documents are being developed to provide the technical underpinning for the analyses of contaminant releases to the environment from these waste forms for the IDF PA. Traditionally, releases from cementitious waste forms have utilized a diffusion-limited release model using contaminant-specific diffusivities measured in standard leach tests as inputs. More mechanistic modeling techniques are under development to increase the rigor of the PA analyses. In addition, recent testing has shown diffusivities of key contaminants can be reduced by several orders-of-magnitude through optimized formulation of the waste form and use of additives to bind specific contaminants. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Reactive nanolaminate pulsed-laser ignition mechanism: Modeling and experimental evidence of diffusion limited reactions

We irradiated Al/Pt nanolaminates with a bilayer thickness (tb, width of an Al/Pt pair-layer) of 164 nm with single laser pulses with durations of 10 ms and 0.5 ms at 189 W/cm 2 and 1189 W/cm 2 , respectively. The time to ignition was measured for each pulse, and shorter ignition times were observed for the higher power/shorter pulse width. While the shorter pulse shows uniform brightness, videographic images of the irradiated area shortly after ignition show a non-uniform radial brightness for the longer pulse. A diffusion-limited single step reaction mechanism was implemented in a finite element package to model the progress from reactants to products at both pulse widths. Finally, the model captures well both the observed ignition delay and qualitative observations regarding the non-uniform radial temperature.

36 MATERIALS SCIENCE↗

Calculations of hydrogen diffusivity in Zr-based alloys: Influence of alloying elements and effect of stress

This report summarizes the progress on modeling hydrogen diffusivity in Zr-based alloys. The presence of hydrogen (H) can detrimentally affect the mechanical properties of many metals and alloys. To mitigate these detrimental effects requires fundamental understanding of the thermodynamics and kinetics governing H pickup and hydride formation. In this work, we focus on H diffusion in Zr-based alloys by studying the effects of alloying elements and stress, factors that have been shown to strongly affect H pickup and hydride formation in nuclear fuel claddings. A recently developed accelerated kinetic Monte Carlo method is used for the study. It is found that for the alloys considered here, H diffusivity depends weakly on composition, with negligible effect at high temperatures in the range of 600-1200 K. Therefore, the small variation in compositions of these alloys is likely not a major cause of the very different H pickup rates. In contrast, stress strongly affects H diffusivity. This effect needs to be considered for studying hydride formation and delayed hydride cracking.

36 MATERIALS SCIENCE↗