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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evidence for the 3D Radiative Effects of Boundary-Layer Clouds From Observations of Direct and Diffuse Surface Solar Fluxes

Numerical experiments have revealed the importance of horizontal transport of light in the presence of clouds (“3D effects”), with consequences for climate, weather, and solar resource availability predictions. Yet, analysis of 3D effects from observations remain sparse because of the difficulty to isolate the effect of horizontal transport in radiation measurements. In this study, we provide observational evidence for 3D effects based on the direct-diffuse partition of surface solar fluxes. It is compared to outputs from the ecRad radiative transfer scheme run on retrieved cloud profiles. Furthermore, the direct-beam calculation takes careful account of the field-of-view of the pyrheliometer to ensure consistency between observed and modeled direct fluxes. Only the solver that accounts for 3D effects is able to reproduce the observed mean direct-diffuse partition as a function of solar zenith angle and cloud cover, in particular at large solar zenith angles where cloud sides intercept most of the direct beam.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Larger than uncorrelated vacancy diffusion contributions in chemically disordered crystalline materials

Arising from a variational approach to compute diffusion coefficients, we introduce “superkinetic kinosons”, which are contributions to the diffusion flux of mobile defects and atomic species by single jump mechanisms that can exceed their uncorrelated values. In vacancy-mediated diffusion in crystalline random and high entropy materials, these contributions primarily arise from intermittent jumps of slow moving atoms that remove the vacancy out of correlation traps and enable a quantification of such effects on the diffusion coefficients. They can be significant contributors to diffusion, even when faster moving atoms have infinite percolation networks. Furthermore, their existence and importance provide fundamental insights about underlying aspects of vacancy diffusion such as formation of localized correlation traps and deviation of vacancy diffusion from percolation behavior.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗
NASA NTRS1992

Helium diffusion in the sun

We calculate improved standard solar models using the new Livermore (OPAL) opacity tables, an accurate (exportable) nuclear energy generation routine which takes account of recent measurements and analyses, and the recent Anders-Grevesse determination of heavy element abundances. We also evaluate directly the effect of the diffusion of helium with respect to hydrogen on the calculated neutrino fluxes, on the primordial solar helium abundance, and on the depth of the convective zone. Helium diffusion increases the predicted event rates by about 0.8 SNU, or 11 percent of the total rate, in the chlorine solar neutrino experiment, by about 3.5 SNU, or 3 percent, in the gallium solar neutrino experiments, and by about 12 percent in the Kamiokande and SNO solar neutrino experiments. The best standard solar model including helium diffusion and the most accurate nuclear parameters, element abundances, and radiative opacity predicts a value of 8.0 SNU +/- 3.0 SNU for the C1-37 experiment and 132 +21/-17 SNU for the Ga - 71 experiment, where the uncertainties include 3 sigma errors for all measured input parameters.

Bahcall, J. N.↗
DOE OSTI2024

HTO and selenate diffusion through compacted Na-, Na–Ca-, and Ca-montmorillonite

Radionuclide transport in smectite clay barrier systems used for nuclear waste disposal is controlled by diffusion, with adsorption significantly retarding transport rates. While a relatively minor component of spent nuclear fuel, 79 Se is a major driver of the safety case for spent fuel disposal due to its long half-life (3.3×10 5 yr) and its low adsorption to clay (K D < 10 L/kg), thus a thorough understanding of Se diffusion through clay is critical for understanding the long-term safety of spent fuel disposal systems. Through-diffusion experiments with tritiated water (HTO, conservative tracer) and Se(VI) were conducted with a well-characterized, purified montmorillonite source clay (SWy-2) under a constant ionic strength (0.1 M) and three different electrolyte compositions: Na + , Ca 2+ , and a Na + -Ca 2+ mixture at pH 6.5 in order to probe the effects of electrolyte composition and interlayer cation composition on clay microstructure, Se(VI) aqueous speciation, and ultimately diffusion. Further, the results were modeled using a reactive transport modeling approach to determine values of porosity (ε), D e (effective diffusion coefficient), and K D (distribution coefficient for adsorption). HTO diffusive flux was higher in Ca-montmorillonite (D e =1.68×10 -10 m 2 s -1 ) compared to Na-montmorillonite (De=7.83×10 -11 m 2 s -1 ). This increase in flux is likely due to a greater degree of clay layer stacking in the presence of Ca 2+ compared to Na + , which leads to larger inter-particle pores. Overall, the Se(VI) flux was much lower than the HTO flux due to anion exclusion, with Se(VI) flux following the order Ca (D e = 1.03×10 -11 m 2 s -1 ) > Na–Ca (D e = 2.12×10 -12 m 2 s -1 ) > Na (D e = 1.28×10 -12 m 2 s -1 ). These differences in Se(VI) flux are due to a combination of factors, including (1) larger accessible porosity in Ca-montmorillonite due to clay layer stacking and smaller electrostatic effects compared to Na-montmorillonite, (2) larger accessible porosity for neutral-charge CaSeO4 species which makes up 32% of aqueous Se(VI) in the pure Ca system, and (3) possibly higher Se(VI) adsorption for Ca-montmorillonite. Through a combination of experimental and modeling work, this study highlights the compounding effects that electrolyte and counterion compositions can have on radionuclide transport through clay. Diffusion models that neglect these effects are not transferable from laboratory experimental conditions to in situ repository conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

A diffusion–deformation model with damage for polymer undergoing rapid decompression failure

Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multiscale, mechanistic modeling of cesium transport in silicon carbide for TRISO fuel performance prediction

Understanding cesium (Cs) transport in TRistructural ISOtropic (TRISO) particle fuel is crucial for predicting fission product release in high-temperature reactors. However, current challenges include significant scatter in diffusivity data and unexplained temperature-dependent diffusion regimes in the silicon carbide layer. This study addresses these challenges by developing a multiscale, mechanistic Cs transport model integrating atomistic simulations and phase field modeling. Our model quantifies temperature and grain size effects on Cs diffusivity, attributing experimentally observed regimes to a transition from bulk-dominated diffusivity at high temperatures to grain boundary-dominated diffusivity at lower temperatures. The model, validated against diffusion measurements and advanced gas reactor (AGR)-1 and AGR-2 post-irradiation fission product release data, enhances the predictive capability of the BISON fuel performance code. This study advances our understanding of Cs release from TRISO particles and its dependence on temperature and silicon carbide grain size, with implications for the safety and efficiency of high-temperature nuclear reactors.

BISON↗
NASA NTRS1972

Instabilities in uranium plasma and the gas-core nuclear rocket engine

The nonlinear evolution of unstable sound waves in a uranium plasma has been calculated using a multiple time-scale asymptotic expansion scheme. The fluid equations used include the fission power density, radiation diffusion, and the effects of the changing degree of ionization of the uranium atoms. The nonlinear growth of unstable waves is shown to be limited by mode coupling to shorter wavelength waves which are damped by radiation diffusion. This mechanism limits the wave pressure fluctuations to values of order delta P/P approximates 0.00001 in the plasma of a typical gas-core nuclear rocket engine. The instability is thus not expected to present a control problem for this engine.

Tidman, D. A.↗
NASA NTRS1971

Instabilities in uranium plasma.

The nonlinear evolution of unstable sound waves in a uranium plasma has been calculated using a multiple time-scale asymptotic expansion scheme. The fluid equations used include the fission power density, radiation diffusion, and the effects of the changing degree of ionization of the uranium atoms. The nonlinear growth of unstable waves is shown to be limited by mode coupling to shorter wavelength waves which are damped by radiation diffusion. This mechanism limits the wave pressure fluctuations to values of order delta P/P equal to about .00001 in the plasma of a typical gas-core nuclear rocket engine. The instability is thus not expected to present a control problem for this engine.

Tidman, D. A.↗
NASA NTRS1987

Impeller shroud to casing leakage flow simulations in the Space Shuttle Main Engine high pressure fuel pump

Quasi-three-dimensional Navier-Stokes calculations were carried out for the Space Shuttle Main Engine high-pressure fuel pump to simulate the impeller shroud to casing leakage flow. This flow geometry was modeled as an axisymmetric cavity flow with a stationary surface representing the casing, and a rotating surface denoting the impeller. A 63 x 81-node mesh provided sufficient resolution in the regions of greatest flow variations and reduced the effects of numerical diffusion. The turbulence field was closed with the high Reynolds number form of the k-epsilon model supplemented with wall functions in the vicinity of the walls. Finally, a parametric study quantified the effects of through mass flow changes on this leakage flow.

Sindir, Munir M.↗
NASA NTRS1996

A Computational-Experimental Development of Vortex Generator Use for a Transitioning S-Diffuser

The development of an effective design strategy for surface-mounted vortex generator arrays in a subsonic diffuser is described in this report. This strategy uses the strengths of both computational and experimental analyses to determine beneficial vortex generator locations and sizes. A parabolized Navier-Stokes solver, RNS3D, was used to establish proper placement of the vortex generators for reduction in circumferential total pressure distortion. Experimental measurements were used to determine proper vortex generator sizing to minimize total pressure recovery losses associated with vortex generator device drag. The best result achieved a 59% reduction in the distortion index DC60, with a 0.3% reduction in total pressure recovery.

Wendt, Bruce J.↗
NASA NTRS2005

Calibration of High Temperature Thermal Conductivity System: New Algorithm to Measure Heat Capacity Using Flash Thermal Diffusivity in Thermoelectric Materials

A viewgraph presentation describing thermoelectric materials, an algorithm for heat capacity measurements and the process of flash thermal diffusivity. The contents include: 1) What are Thermoelectrics?; 2) Thermoelectric Applications; 3) Improving Thermoelectrics; 4) Research Goal; 5) Flash Thermal Diffusivity; 6) Background Effects; 7) Stainless Steel Comparison; 8) Pulse Max Integral; and 9) Graphite Comparison Algorithm.

heat capacity↗
NASA NTRS2013

Kinetic Monte Carlo Simulation of Cation Diffusion in Low-K Ceramics

Low thermal conductivity (low-K) ceramic materials are of interest to the aerospace community for use as the thermal barrier component of coating systems for turbine engine components. In particular, zirconia-based materials exhibit both low thermal conductivity and structural stability at high temperature, making them suitable for such applications. Because creep is one of the potential failure modes, and because diffusion is a mechanism by which creep takes place, we have performed computer simulations of cation diffusion in a variety of zirconia-based low-K materials. The kinetic Monte Carlo simulation method is an alternative to the more widely known molecular dynamics (MD) method. It is designed to study "infrequent-event" processes, such as diffusion, for which MD simulation can be highly inefficient. We describe the results of kinetic Monte Carlo computer simulations of cation diffusion in several zirconia-based materials, specifically, zirconia doped with Y, Gd, Nb and Yb. Diffusion paths are identified, and migration energy barriers are obtained from density functional calculations and from the literature. We present results on the temperature dependence of the diffusivity, and on the effects of the presence of oxygen vacancies in cation diffusion barrier complexes as well.

Good, Brian↗
DOE OSTI2026

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗
NASA NTRS1986

Group-kinetic theory of turbulence

The two phases are governed by two coupled systems of Navier-Stokes equations. The couplings are nonlinear. These equations describe the microdynamical state of turbulence, and are transformed into a master equation. By scaling, a kinetic hierarchy is generated in the form of groups, representing the spectral evolution, the diffusivity and the relaxation. The loss of memory in formulating the relaxation yields the closure. The network of sub-distributions that participates in the relaxation is simulated by a self-consistent porous medium, so that the average effect on the diffusivity is to make it approach equilibrium. The kinetic equation of turbulence is derived. The method of moments reverts it to the continuum. The equation of spectral evolution is obtained and the transport properties are calculated. In inertia turbulence, the Kolmogoroff law for weak coupling and the spectrum for the strong coupling are found. As the fluid analog, the nonlinear Schrodinger equation has a driving force in the form of emission of solitons by velocity fluctuations, and is used to describe the microdynamical state of turbulence. In order for the emission together with the modulation to participate in the transport processes, the non-homogeneous Schrodinger equation is transformed into a homogeneous master equation. By group-scaling, the master equation is decomposed into a system of transport equations, replacing the Bogoliubov system of equations of many-particle distributions. It is in the relaxation that the memory is lost when the ensemble of higher-order distributions is simulated by an effective porous medium. The closure is thus found. The kinetic equation is derived and transformed into the equation of spectral flow.

Tchen, C. M.↗
DOE OSTI2022

Dynamics of Glyceline and Interactions of Constituents: A Multitechnique NMR Study

In this work, the dynamics of the organic components of the deep eutectic solvent (DES) glyceline are analyzed using an array of complementary nuclear magnetic resonance (NMR) methods. Fast-Field Cycling 1 H relaxometry, Pulsed Field Gradient diffusion, Nuclear Overhauser Effect Spectroscopy (NOESY), 13 C NMR relaxation and pressure dependent NMR experiments are deployed to sample a range of frequencies and modes of motion of the glycerol and choline components of the DES. Generally, translational and rotational diffusion of glycerol are more rapid than those of choline while short range rotational motions observed from 13 C relaxation indicate slow local motion of glycerol at low choline chloride (ChCl) content. The rates of glycerol and choline local motions become more similar at higher ChCl. This result taken together with pressure dependent NMR studies show that the addition of ChCl makes it easier to disrupt glycerol packing. Finally, a relatively slow hydroxyl H-exchange process between glycerol and choline protons is deduced from the data. Consistent with this, NOESY results indicate relatively little direct H-bonding between glycerol and choline. These results suggest that the glycerol H-bonding network is disrupted as choline is added, but primarily in regions where there is intimate mixing of the two components. Thus, the local dynamics of most of the glycerol, resembles that of pure glycerol until substantial choline chloride is present.

13C NMR relaxation↗
NASA NTRS1999

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

We develop and test a 1-point closure turbulence model with the following features: 1) we include the salinity field and derive the expression for the vertical turbulent diffusivities of momentum K(sub m) , heat K(sub h) and salt K(sub s) as a function of two stability parameters: the Richardson number R(sub i) (stratification vs. shear) and the Turner number R(sub rho) (salinity gradient vs. temperature gradient). 2) to describe turbulent mixing below the mixed layer (ML), all previous models 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, the three diffusivities have the same functional dependence on R( sub i) and R(sub rho) as 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. The procedure frees the model from adjustable background diffusivities and indeed we employ the same model throughout the entire vertical extent of the ocean. 3) in the local model, the turbulent diffusivities K(sub m,h,s) are given as analytical functions of R(sub i) and R(sub rho). 5) the model is used in an O-GCM and several results are presented to exhibit the effect of double diffusion processes. 6) the code is available upon request.

Canuto, V. M.↗