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 19 records

Analytical homogenization techniques applied to the Fickian diffusion: Effective diffusivity coefficient

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. In this work, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing TRISO particles. An extensive survey on effective thermal conductivity modeling was performed in our previous study, considering a random distribution of mono-sized spherical inclusions in a continuous matrix. Using the analogy between heat conduction and the simplified Fickian diffusion (or fission product species conservation), we can use the same analytical homogenization methods to obtain ETC as for the effective diffusivity coefficient (EDC). We performed several numerical experiments at varying conditions to assess the validity of our hypothesis for EDC calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Graph-theoretical KMC Framework for Calculating Effective Diffusivity in TPBAR Components: Effective Diffusivity of Tritium in α-Zr Grain Boundaries

We report the development of a computer simulation tool based on the Graph-theoretical kinetic Monte Carlo (GT-KMC) approach that can be used to simulate diffusion and derive effective diffusivity of species in complex structures of TPBAR components. The GT-KMC framework was implemented in AKSOME, an on-lattice self-learning kinetic Monte Carlo tool developed at PNNL and was previously used to study solute diffusion in metal alloys. This newly developed tool will be helpful in simulating atom diffusion and extracting their effective diffusivities in grain boundaries and interfaces in various TPBAR components. Once benchmarked, the tool will be used to simulate the diffusion of tritium (hydrogen) along grain boundaries in α-Zr using the activation energy barrier data previously obtained from Density Functional Theory calculations.

36 MATERIALS SCIENCE↗

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↗

Predicting Effective Diffusion Coefficients in Mudrocks Using a Fractal Model and Small-Angle Neutron Scattering Measurements

The determination of effective diffusion coefficients of gases or solutes in the water–saturated pore space of mudrocks is time consuming and technically challenging. Yet, reliable values of effective diffusion coefficients are important to predict migration of hydrocarbon gases in unconventional reservoirs, dissipation of (explosive) gases through clay barriers in radioactive waste repositories, mineral alteration of seals to geological CO 2 storage reservoirs and contaminant migration through aquitards. In this study, small angle and very small angle neutron scattering techniques have been utilized to determine a range of transport properties in mudrocks, including porosity, pore size distributions and surface and volume fractal dimensions of pores and grains, from which diffusive transport parameters can be estimated. Using a fractal model derived from Archie's Law, we calculate effective diffusion coefficients from these parameters and compare them to laboratory–derived effective diffusion coefficients for CO 2 , H 2 , CH 4 and HTO on either the same or related mudrock samples. The samples include Opalinus Shale from the underground laboratory in Mont Terri, Switzerland; Boom Clay from a core drilled in Mol, Belgium and a marine claystone cored in Utah, USA. The predicted values were compared to laboratory diffusion measurements. The measured and modeled diffusion coefficients show good agreement, differing generally by less than factor 5. Furthermore, neutron or X–ray scattering analysis is therefore proposed as a novel method for fast, accurate estimation of effective diffusion coefficients in mudrocks, together with simultaneous measurement of multiple transport parameters including porosity, pore size distributions and surface areas, important for (reactive) transport modelling.

58 GEOSCIENCES↗

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

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

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effective diffusion equation in a random velocity field

The effects are studied of assumed random velocity fields on diffusion in a binary fluid. Random velocity fields can result, for example, from the high-frequency components of residual accelerations onboard spacecraft (often called g-jitter). An effective diffusion equation is derived for an average concentration which includes spatial and temporal correlations induced by the fluctuating velocity fields assumed to be Gaussianly distributed. The resulting equation becomes nonlocal, and if correlations between different components of the velocity field exist, it is also anisotropic. The simple limiting case of short correlation times is discussed and an effective diffusivity is obtained which reflects the enhanced mixing caused by the velocity fields. The results obtained in the limit of short correlation times are valid even if the probability distribution of the velocity field is not Gaussian.

Vinals, Jorge↗

Iodine effective diffusion coefficients through volcanic rock: Influence of iodine speciation and rock geochemistry

Accurate prediction of the subsurface transport of iodine species is important for the assessment of long-term nuclear waste repository performance, as well as monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty, given that radioiodine decays into radioxenon. However, the transport of iodine through intact geologic media is not well understood, compromising our ability to assess risk associated with radioiodine migration. The current study's goal is to quantify the matrix diffusion of iodine species through saturated volcanic rock, with particular attention paid to the redox environment and potential speciation changes. Here, diffusion experiments were run for iodide through lithophysae-rich lava, lithophysae-poor lava, and welded tuff, whereas iodate diffusion was studied through welded tuff. Iodine transport was compared with a conservative tracer, HDO, and effective diffusion coefficients were calculated. Likely due to a combination of size and anion exclusion effects, iodine species diffused more slowly than the conservative tracer through all rock types tested. Furthermore, oxidation of iodide to iodate was observed in the lithophysae-poor lava, affecting transport. Results provide much needed data for subsurface transport models that predict radioiodine migration from underground sources, and indicate the pressing need for geochemical and redox interactions to be incorporated into these models.

58 GEOSCIENCES↗

Effective Diffusion Coefficient of Dissolved Metal Ions during 1D Pit Growth

Transport-limited dissolution of metals in corrosion pits is associated with the product of the metal ion diffusion coefficient and total concentration of metal ions (DM+·Ctot). In one dimensional (1D) pit dissolu-tion, the concentration of ionic species varies with location in the pit; at the dissolving pit bottom, the saturation of metal ion and supporting anion is ex-pected under transport control conditions, while the metal ion is nearly depleted at the pit mouth. To describe metal ion transport through pit solution with constantly-changing concentration, an effective diffusion coefficient (Deff) is proposed in this study with consideration of the combined effects of vis-cosity, temperature and electric field inside a grow-ing 1D pit. The modeled Deff·Ctot values were calcu-lated and compared with experimental DM+·Ctot val-ues for S13Cr 1D pit dissolution in different bulk chloride concentrations at 85°C. The modeled Deff·Ctot values were very close to the experimental values, validating the modeling approach suggested in this study.

Jun, Jiheon↗

Numerical Study of Buoyancy and Differential Diffusion Effects on the Structure and Dynamics of Triple Flames

Triple flames arise in a number of practical configurations where fuel and oxidizer are partially premixed, such as in the base of a lifted jet flame. Past experimental studies, theoretical analyses, and numerical modeling of triple flames suggested the potential role of triple flames in stabilizing turbulent flames and in promoting flame propagation. From recent numerical simulations of laminar triple flames, a strong influence of differential diffusion among species and heat on the triple flame structure has been gradually appreciated. This paper reports preliminary numerical results on the influence of gravity and differential diffusion effects on the structure and dynamics of triple flames with a one-step global irreversible chemistry model.

Chen, J. -Y.↗

Convection of tin in a Bridgman system. I - Flow characterization by effective diffusivity measurements

An electrochemical titration method was used to investigate the dynamic states in a cylindrical layer of convecting tin. The liquid tin was contained in a cell, with curved boundaries made of quartz and flat boundaries made of a solid state electrolyte - yttria-stabilized zirconia (YSZ). The electrolyte acted as a window through which a trace amount of oxygen could be pumped in or out by the application of a constant voltage. The concentration at the YSZ interface was monitored by operating the electrochemical cell in the galvanic mode. Experimentally determined effective diffusivities of oxygen were compared with the molecular diffusivity. Dynamic states in the convective flow were thus inferred. Temperature measurements were simultaneously made in order to identify the onset of oscillations from a steady convective regime. The experiments were conducted for two different aspect ratios for various imposed temperature gradients and two different orientations with respect to gravity. Transcritical states were identified and comparison to two-dimensional numerical models were made.

Sears, B.↗

Thermal /Soret/ diffusion effects on interfacial mass transport rates

It is shown that thermal (Soret) diffusion significantly alters convective mass transport rates and important transition temperatures in highly nonisothermal flow systems involving the transport of 'heavy' species (vapors or particles). Introduction of the Soret transport term is shown to result in mass transfer effects similar to those of 'suction' and a homogeneous chemical 'sink'. It is pointed out that this analogy provides a simple method of correlating and predicting thermal diffusion effects in the abovementioned systems.

Rosner, D. E.↗

Double-diffusive effects during solidification

When a fluid contains two diffusing components with different molecular diffusivities, such as heat and solute concentration, convective motion may be generated when potential energy is released owing to differential diffusion. In the solidification of binary alloys and doped semiconductors, temperature and concentration gradients do exist simultaneously. If these gradients are aligned in a suitable manner, convection may ensue. In the case of the binary alloy, such convection may extend through the mushy zone and may be the cause of macrosegregations. In the case of doped semiconductors, the convective motion may cause the nonuniform distribution of dopant in the product. In this paper, the fundamentals of double-diffusive convection will be reviewed, and its effect on solidification will be discussed.

Chen, C. F.↗

Reduction of Porosity in Cementitious Waste Forms to Reduce Effective Diffusivity

Although current practices meet regulations for grout waste forms, development of sequestration materials and techniques for contaminants of concern could expand the use of grout waste forms and reduce risk to established processes from future compositional changes. There are contaminants in waste streams that are not permanently and completely sequestered in current grout waste forms under certain conditions. These contaminants include methyl mercury, 129 I, 99 Tc, and nitrate. Each of these species have minor fractions that are not permanently and completely isolated from the environment and therefore potentially contribute to long-term environmental impacts. An approach of interest includes the technique of reducing porosity to result in reduced leachability of all contaminants simultaneously. The current Savannah River Site (SRS) saltstone formulation was used to evaluate the relationship between leachability and porosity reduced waste forms by altering the water-to-binder ratio. Results indicated decreasing the water-to-binder directly reduced porosity and resulted in the waste form’s ability to more effectively immobilize nitrate. This work was funded by the Department of Energy Office of Environmental Management (DOE-EM) Technology Operations Office (TOO) Project # HQ221818.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Diffusion effects on the helium abundance of the solar transition region and corona

The diffusion of helium in the solar transition region is studied by solving the mass and momentum conservation equations for a hydrogen-helium plasma given a representative temperature profile. Steady state solutions show that two distinct atmospheres may result. In cases where the thermal force on alpha-particles is balanced by the partial pressure gradient force, helium is the dominant coronal species. On the other hand, if it is the frictional force between protons and alpha-particles which balances the thermal force on alpha-particles then hydrogen is the major coronal component. In order to explore which of these solutions are attainable within reasonable time scales, the time-dependent equations are solved, starting from an initial state with a uniform helium abundance of 10 percent. The atmosphere as a whole is close to hydrostatic equilibrium, but due the thermal forces the individual elements are not. This force inbalance leads to a differential flow between species. It is found that this differential flow leads to a significant enhancement of the coronal helium abundance. Even for the relatively shallow temperature gradient used the helium abundance in the lower corona increases to 30 percent over a 24 hr period.

Hansteen, Viggo H.↗

A Method for Including Diffusive Effects in Texture Evolution

This work describes a model for representing grain misorientation in continuum models for texture evolution. The model relies on a connection, derived here, between the fluctuation of the reorientation velocity in discrete models of texture evolution and an isotropic diffusion operator that can be used to represent misorientation effects in continuum models for texture evolution. Previous work has developed a method for computing the fluctuation of the reorientation velocity in discrete models. These values, along with the average reorientation velocity, are used to set up inhomogeneous continuum advection and diffusion operators. Here, the complete homogenized continuum model uses the discrete harmonics to represent and evolve the orientation distribution function, now including a diffusion operator representing the effects of intragranular misorientation. Results from the new model are compared to fully-resolved crystal plasticity finite element simulations and experimental data to demonstrate the accuracy of the proposed diffusion operator and the overall method for simulating texture evolution.

(B) crystal plasticity↗