Search NASA⌕ Search

SEARCH · Search NASA

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗
NASA NTRS2001

Effects of Flame Structure and Hydrodynamics on Soot Particle Inception and Flame Extinction in Diffusion Flames

This paper summarizes recent studies of the effects of stoichiometric mixture fraction (structure) and hydrodynamics on soot particle inception and flame extinction in diffusion flames. Microgravity experiments are uniquely suited for these studies because, unlike normal gravity experiments, they allow structural and hydrodynamic effects to be independently studied. As part of this recent flight definition program, microgravity studies have been performed in the 2.2 second drop tower. Normal gravity counterflow studies also have been employed and analytical and numerical models have been developed. A goal of this program is to develop sufficient understanding of the effects of flame structure that flames can be "designed" to specifications - consequently, the program name Flame Design. In other words, if a soot-free, strong, low temperature flame is required, can one produce such a flame by designing its structure? Certainly, as in any design, there will be constraints imposed by the properties of the available "materials." For hydrocarbon combustion, the base materials are fuel and air. Additives could be considered, but for this work only fuel, oxygen and nitrogen are considered. Also, the structure of these flames is "designed" by varying the stoichiometric mixture fraction. Following this line of reasoning, the studies described are aimed at developing the understanding of flame structure that is needed to allow for optimum design.

Axelbaum, R. L.↗
DOE OSTI2023

Pressure effects on reactivity and extinction of n -dodecane diffusion cool flame

Low-temperature combustion at high pressure has gained increasing interest owing to the growing demand for clean and efficient propulsion and energy systems. In particular, cool flame structure, extinction limits, and fuel reactivities at high pressure are critical in affecting the engine performance at near-limit conditions. Here, this paper investigates the effects of pressure on the cool flame extinction limit, structure, radical index, reactivity, and oxygen concentration dependence by experiments, analysis, and modeling. A large n-alkane, namely n-dodecane (n C 12 H 26 ), is selected to study its diffusion cool flame dynamics and reactivity in a high-pressure counterflow burner up to 10 atm. The experimental results show that higher pressure increases the cool flame extinction strain rates, and that the pressure-weighted extinction strain rate (aP) is proportional to the square of pressure, P 2 . A scaling analysis explains the relationship between the dependence of flame structure, heat release rate, and pressure-weighted cool flame strain rate on pressure. Furthermore, radical indexes at different pressures are measured by isolating the thermal and transport effects from the chemical contribution to diffusion cool flame extinction. The radical index clearly shows that the low-temperature reactivity increases with pressure. In addition, due to the critical role of multiple oxygen addition reactions in low-temperature chemistry, the relationship between the cool flame extinction limit and the oxygen concentration is explored. It is found that the cool flame extinction limits are proportional to the n th power of the oxygen concentration, [O 2 ] n , and increasing pressure leads to stronger extinction limit dependence (larger n) on the oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions, while the scaling analysis explains well the experimental results.

33 ADVANCED PROPULSION SYSTEMS↗
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.↗