Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion rate”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Windows(Registered Trademark)-Based Software Models Cyclic Oxidation Behavior

Oxidation of high-temperature aerospace materials is a universal issue for combustion-path components in turbine or rocket engines. In addition to the question of the consumption of material due to growth of protective scale at use temperatures, there is also the question of cyclic effects and spallation of scale on cooldown. The spallation results in the removal of part of the protective oxide in a discontinuous step and thereby opens the way for more rapid oxidation upon reheating. In experiments, cyclic oxidation behavior is most commonly characterized by measuring changes in weight during extended time intervals that include hundreds or thousands of heating and cooling cycles. Weight gains occurring during isothermal scale-growth processes have been well characterized as being parabolic or nearly parabolic functions of time because diffusion controls reaction rates. In contrast, the net weight change in cyclic oxidation is the sum of the effects of the growth and spallation of scale. Typically, the net weight gain in cyclic oxidation is determined only empirically (that is, by measurement), with no unique or straightforward mathematical connection to either the rate of growth or the amount of metal consumed. Thus, there is a need for mathematical modeling to infer spallation mechanisms. COSP is a computer program that models the growth and spallation processes of cyclic oxidation on the basis of a few elementary assumptions that were discussed in COSP: A Computer Model of Cyclic Oxidation, Oxidation of Metals, vol. 36, numbers 1 and 2, 1991, pages 81-112. Inputs to the model include the selection of an oxidation-growth law and a spalling geometry, plus oxide-phase, growth-rate, cycle-duration, and spall-constant parameters. (The spalling fraction is often shown to be a constant factor times the existing amount of scale.) The output of COSP includes the net change in weight, the amounts of retained and spalled oxide, the total amounts of oxygen and metal consumed, and the terminal rates of weight loss and metal consumption.

Smialek, J. L.↗
DOE OSTI2024

Morphology selection in dealloying: A phase field study of the coupling among kinetic mechanisms

A phase field model is used to investigate how the corrosion rate and morphology of a dealloying metal is controlled by the kinetic rates of bulk solid diffusion, interfacial diffusivity, and interface dissolution. A specific focus is investigating the manner by which variations in these kinetic parameters synergistically influence the onset of a morphological instability at the solid-liquid interface, which is known to precede rapid dealloying corrosion towards highly topologically complex structures such as bicontinuous porous solids. This work demonstrates that this instability can be suppressed by either an increase of interdiffusion in the solid alloy, or by a decreasing rate of dissolution at the solid-liquid interface. When the instability is suppressed, the result is a relatively planar, passivated dealloying front. Importantly, however, coupled changes in kinetic parameters can also promote the morphological instability and sustain fine-scale ligament and pore formation. Finally, the relevance of the solid-liquid interfacial diffusion for determining the rate of dealloying and ligament formation is highlighted. By elucidating the mechanistic influence of bulk and interfacial kinetic factors on the dealloyed morphology, this work helps to rationalize the contrasting morphologies observed in various dealloying systems.

Diffusion↗
NASA NTRS1985

Simulations of magnetic-flux transport in solar active regions

The evolution of several observed solar active regions is simulated by solving a transport equation for magnetic flux at the photosphere. The rates of rotation, meridional flow, and diffusion of the flux are determined self-consistently in the calculations. The findings are in good quantitative agreement with previous measures of the rotation rate and diffusion constant associated with photospheric magnetic fields. Although the meridional velocities are consistent in direction and magnitude with recently reported poleward flows, relatively large uncertainties in the velocity determinations make this result inconclusive.

Devore, C. R.↗
NASA NTRS1995

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗
NASA NTRS1987

Combustion rate limits of hydrogen plus hydrocarbon fuel: Air diffusion flames from an opposed jet burner technique

Combustion of H2/hydrocarbon (HC) fuel mixtures may be considered in certain volume-limited supersonic airbreathing propulsion applications. Effects of HC addition to H2 were evaluated, using a recent argon-bathed, coaxial, tubular opposed jet burner (OJB) technique to measure the extinction limits of counterflow diffusion flames. The OJB flames were formed by a laminar jet of (N2 and/or HC)-diluted H2 mixture opposed by a similar jet of air at ambient conditions. The OJB data, derived from respective binary mixtures of H2 and methane, ethylene, or propane HCs, were used to characterize BLOWOFF and RESTORE. BLOWOFF is a sudden breaking of the dish-shaped OJB flame to a stable torus or ring shape, and RESTORE marks sudden restoration of the central flame by radial inward flame propagation. BLOWOFF is a measure of kinetically-limited flame reactivity/speed under highly stretched, but relatively ideal impingement flow conditions. RESTORE measures inward radial flame propagation rate, which is sensitive to ignition processes in the cool central core. It is concluded that relatively small molar amounts of added HC greatly reduce the reactivity characteristics of counterflow hydrogen-air diffusion flames, for ambient initial conditions.

Pellett, Gerald L.↗
DOE OSTI2021

Molecular dynamics simulation of the effect of cooling rate on the structure and properties of lithium disilicate glass

The effects of cooling rate on the structure and properties of lithium disilicate (LS2) glass are investigated using molecular dynamics (MD) computer simulations. The evolution of structural features such as pair distribution function, bond angle distribution, and Li coordination distribution are determined, and correlated with dynamic and static properties to elucidate the effects of cooling rate. The density, elastic moduli, and diffusion coefficient are found to be highly sensitive to cooling rate, whereas the Si pair distribution function and bond angle distribution are weakly affected by the cooling rate. Additionally, the changes of Si-O-Si bond angle and Li coordination number suggest the formation of Li cluster at lower cooling rates. Furthermore, by comparing results from other simulations and reported experiments, we confirm that the increase of cooling rate leads to an increase of conductivity and a decrease of density. Finally, at very high cooling rates, we find that all the atoms in LS2-glass do not relax simultaneously, but do so in two distinct configurations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Parabolic oxidation kinetics of a plutonium alloy at room temperature

Auger electron spectroscopy (AES) was used to characterize the oxide scale formed on an ~0.5 wt. % Ga plutonium metal alloy. Analysis of changes to the characteristic Auger transitions for Pu and O during sputter depth profiles were used in combination with previously measured sputter yields to determine the thickness of the oxide scale and provide insight into its chemistry. From this data, the parabolic kinetics for the Pu-Ga alloy were determined as a function of relative humidity, RH = 14–92%, and oxygen partial pressure, p(O 2 ) = 31–157 Torr, at 25°C. For each relative humidity investigated, the parabolic rate constant was found to be independent of oxygen partial pressure. Furthermore, power-law fits to the relative humidity effect on the rate constant was found to have a power of n = 1.265, indicating the diffusion species which rate-limits the reaction isn’t molecular water, which subsequent dissociates at the oxide-metal interface, or hydroxyl from dissociation at the gas-surface interface.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Effects of a temperature-dependent rheology on large scale continental extension

The effects of a temperature-dependent rheology on large-scale continental extension are investigated using a thin viscous sheet model. A vertically-averaged rheology is used that is consistent with laboratory experiments on power-law creep of olivine and that depends exponentially on temperature. Results of the calculations depend principally on two parameters: the Peclet number, which describes the relative rates of advection and diffusion of heat, and a dimensionless activation energy, which controls the temperature dependence of the rheology. At short times following the beginning of extension, deformation occurs with negligible change in temperature, so that only small changes in lithospheric strength occur due to attenuation of the lithosphere. However, after a certain critical time interval, thermal diffusion lowers temperatures in the lithosphere, strongly increasing lithospheric strength and slowing the rate of extension. This critical time depends principally on the Peclet number and is short compared with the thermal time constant of the lithosphere. The strength changes cause the locus of high extensional strain rates to shift with time from regions of high strain to regions of low strain. Results of the calculations are compared with observations from the Aegean, where maximum extensional strains are found in the south, near Crete, but maximum present-day strain rates are largest about 300 km further north.

Sonder, Leslie J.↗
NASA NTRS1989

Effects of a temperature-dependent rheology on large-scale continental extension

The effects of a temperature-dependent rheology on large-scale continental extension are investigated using a thin viscous sheet model. A vertically-averaged rheology is used that is consistent with laboratory experiments on power-law creep of olivine and that depends exponentially on temperature. Results of the calculations depend principally on two parameters: the Peclet number, which describes the relative rates of advection and diffusion of heat, and a dimensionless activation energy, which controls the temperature dependence of the rheology. At short times following the beginning of extension, deformation occurs with negligible change in temperature, so that only small changes in lithospheric strength occur due to attenuation of the lithosphere. However, after a certain critical time interval, thermal diffusion lowers temperatures in the lithosphere, strongly increasing lithospheric strength and slowing the rate of extension. This critical time depends principally on the Peclet number and is short compared with the thermal time constant of the lithosphere. The strength changes cause the locus of high extensional strain rates to shift with time from regions of high strain to regions of low strain. Results of the calculations are compared with observations from the Aegean, where maximum extensional strains are found in the south, near Crete, but maximum present-day strain rates are largest about 300 km further north.

Sonder, Leslie J.↗
NASA NTRS1995

Use of the quartz crystal microbalance to determine the monomeric friction coefficient of polyimides

When a thin film of polymer is coated on to a quartz crystal microbalance (QCM), the QCM can be used to detect the rate of increase in weight of the polymer film as the volatile penetrant diffuses into the polymer. From this rate information the diffusion coefficient of the penetrant into the polymer can be computed. Calculations requiring this diffusion coefficient lead to values which approximate the monomeric friction coefficient of the polymer. This project has been concerned with the trial of crystal oscillating circuits suitable for driving polymer coated crystals in an atmosphere of penetrant. For these studies done at room temperature, natural rubber was used as an easily applied polymer that is readily penetrated by toluene vapors, qualities anticipated with polyimides when they are tested at T(g) in the presence of toluene. Three quartz crystal oscillator circuits were tested. The simplest circuit used +/- 5 volt dc and had a transistor to transistor logic (TTL) inverter chip that provides a 180 deg phase shift via a feed back loop. This oscillator circuit was stable but would not drive the crystal when the crystal was coated with polymer and subjected to toluene vapors. Removal of a variable resistor from this circuit increased stability but did not otherwise increase performance. Another driver circuit tested contained a two stage differential input, differential output, wide band video amplifier and also contain a feed back loop. The circuit voltage could not be varied and operated at +/- 5 volts dc; this circuit was also stable but failed to oscillate the polymer coated crystal in an atmosphere saturated with toluene vapors. The third oscillator circuit was of similar construction and relied on the same video amplifier but allowed operation with variable voltage. This circuit would drive the crystal when the crystal was submerged in liquid toluene and when the crystal was coated with polymer and immersed in toluene vapors. The frequency readings obtained when using this oscillating circuit are highly variable. This circuit requires further modification to stabilize frequency readings before its use in studies to determine the diffusion coefficient of penetrant molecules into a polymer film coated on a QCM.

Bechtold, Mary M.↗
NASA NTRS1974

Thermodynamics and kinetics of reactions in protective coating systems

Investigations of fluoride activated packs with Al:Ni ratios greater than 50 a/o prove that the specimen surface is not in equilibrium with the pack at high Al:Ni ratios but that an activity gradient exists between pack and specimen. Therefore, gaseous diffusion and possibly surface reactions play a role in determining the overall rate of Al deposition in such packs. Noticeable differences in coating behavior have been obtained in packs activated with chloride and iodide, and it appears that poorest results are obtained with iodides, better with chlorides, and best with fluorides. A numerical method has been perfected for calculating rates of solid-state diffusion controlled coating formation, allowing for the variation of diffusivity with composition in the NiAl phase. Layer growth rates can now be accurately predicted from a knowledge of the surface and substrate compositions. Furthermore, the correct diffusion profiles are obtained by this method. These differ substantially from the profile obtained when the diffusivity is assumed constant.

Gupta, B.↗
NASA NTRS2002

Effect of Wind Velocity on Flame Spread in Microgravity

A three-dimensional, time-dependent model is developed describing ignition and subsequent transition to flame spread over a thermally thin cellulosic sheet heated by external radiation in a microgravity environment. A low Mach number approximation to the Navier Stokes equations with global reaction rate equations describing combustion in the gas phase and the condensed phase is numerically solved. The effects of a slow external wind (1-20 cm/s) on flame transition are studied in an atmosphere of 35% oxygen concentration. The ignition is initiated at the center part of the sample by generating a line-shape flame along the width of the sample. The calculated results are compared with data obtained in the 10s drop tower. Numerical results exhibit flame quenching at a wind speed of 1.0 cm/s, two localized flames propagating upstream along the sample edges at 1.5 cm/s, a single line-shape flame front at 5.0 cm/s, three flames structure observed at 10.0 cm/s (consisting of a single line-shape flame propagating upstream and two localized flames propagating downstream along sample edges) and followed by two line-shape flames (one propagating upstream and another propagating downstream) at 20.0 cm/s. These observations qualitatively compare with experimental data. Three-dimensional visualization of the observed flame complex, fuel concentration contours, oxygen and reaction rate isosurfaces, convective and diffusive mass flux are used to obtain a detailed understanding of the controlling mechanism, Physical arguments based on lateral diffusive flux of oxygen, fuel depletion, oxygen shadow of the flame and heat release rate are constructed to explain the various observed flame shapes.

Prasad, Kuldeep↗
NASA NTRS1983

The problem of cooling the cold Io torus

Two models are developed for the transport of ions inward from Io's orbit, the first using the assumption that radial diffusion is the dominant transport mechanism while the second uses a combination of diffusive and convective transport. The models include thermal and number density transport, radiation, ionization and pickup of local neutrals, recombination, charge exchange, and Coulomb interactions. It is found that pure diffusive transport can account for the dramatic inward depletion of the torus only by invoking recombination or by postulating a massive increase in the production rate of torus ions sometime prior to the Voyager mission. Recombination is determined to be far too slow to be the cause of the observed density decrease inside of Io. The model combining convection and diffusion is shown to reasonably match the data, but only if a diffusion coefficient 100 times less than that derived from Pioneer observations is used. In addition, it is shown that the Pioneer derived diffusion rate combined with Voyager temperature and density measurements indicate a large nonradiative sink of energy in the inner torus.

Richardson, J. D.↗
NASA NTRS2003

In General, the Total Voltammetric Current from a Mixture of Redox-Active Substances will Not be the Sum of the Currents that Each Substance would Produce Independently at the Same Concentration as in the Mixture

At the potential range where both decamethylferrocene (dMeFc) and ferrocene (Fc) are oxidized with rates controlled by linear diffusion, electrogenerated Fc(+) radicals diffusing outwards from the electrode react quantitatively (K23 C=5.8 x 10(exp 8) with dMeFc diffusing towards the electrode and produce Fc and dMeFc. That reaction replaces dMeFc with Fc, whose diffusion coefficient is higher than that of dMeFc(+), and the total mass-transfer limited current from the mixture is increased by approximately 10%. Analogous observations are made when mass-transfer is controlled by convective-diffusion as in RDE voltammetry. Similar results have been obtained with another, and for all practical purposes randomly selected pair of redox-active substances, [Co(bipy)3](2+) and N - methylphenothiazine (MePTZ); reaction of MePTZ(+) with [Co(bipy)3](2+) replaces the latter with MePTZ, which diffuses faster and the current increases by approximately 20%. The experimental voltammograms have been simulated numerically and the role of (a) the rate constant of the homogeneous reaction; (b) the relative concentrations; and, (c) the diffusion coefficients of all species involved have been studied in detail. Importantly, it was also identified that within any given redox system the dependence of the mass-transfer limited current on the bulk concentrations of the redox-active species is expected to be non-linear. These findings are discussed in terms of their electroanalytical implications.

Leventis, Nicholas↗
NASA NTRS1994

Long-term modulation of galactic cosmic radiation and its model for space exploration

As the human exploration of space has received new attention in the United States, studies find that exposure to space radiation could adversely impact the mission design. Galactic Cosmic Radiation (GCR), with its very wide range of charges and energies, is particularly important for a mission to Mars, because it imposes a stiff mass penalty for spacecraft shielding. Dose equivalent versus shielding thickness calculations, show a rapid initial drop on exposure with thickness, but an asymptotic behavior at a higher shielding thickness. Uncertainties in the radiobiology are largely unknown. For a fixed radiation risk, this leads to large uncertainties in shielding thickness for small uncertainties in estimated dose. In this paper we investigate the application of steady-state, spherically-symmetric diffusion-convection theory of solar modulation to individual measurements of differential energy spectra from 1954 to 1989 in order to estimate the diffusion coefficient, kappa(r,t), as a function of time. We have correlated the diffusion coefficient to the Climax neutron monitor rates and show that, if the diffusion coefficient can be separated into independent functions of space and time: kappa(r,t)=K(t) k(sub 0) beta Pk(sub 1) (r), where beta is the particle velocity and P the rigidity, then (i) The time dependent quantity 1/K(t), which is proportional to the deceleration potential, phi(r,t), is linearly related to the Climax neutron monitor counting rate. (ii) The coefficients obtained from hydrogen or helium intensity measurements are the same. (iii) There are different correlation functions for odd and even solar cycles. (iv) The correlation function for the Climax neutron monitor counting rate for given time, t, can be used to estimate mean deceleration parameter phi(t) to within +/- 15% with 90% confidence. We have shown that k(r,t) determined from hydrogen and/or helium data, can be used to fit the oxygen and iron differential energy spectra with a root mean square error of about +/- 10%, and essentially independent of the particle charge or energy. We have also examined the ion chamber and C-14 measurements which allow the analysis to be extended from the year 1906 to 1990. Using this model we have defined reference GCR spectra at solar minimum and solar maximum. These can be used for space exploration studies and provide a quantitative estimate of the error in dose due to changes in GCR intensities.

Badhwan, G. D.↗
DOE OSTI2026

Role of Wadsley Defects and Cation Disorder to Enhance MoNb 12 O 33 Diffusion

Wadsley-Roth (WR) niobates have emerged as high-rate anode materials that can combine rapid ionic diffusion with good electronic conductivity. WR compounds have been defect-enhanced by limited annealing, however, such materials often contain multiple types of defects. In particular, both Wadsley defects (variable block size) and transition metal disorder have the potential to modify transport rates, however the corresponding effects are not well understood mechanistically. Here, MoNb 12 O 33 (MNO) was calcined at two different temperatures to compare a defect-rich condition (MNO-800) with a proximal order-rich condition (MNO-900) as assessed through XRD, XANES, EXAFS, and STEM characterizations. Galvanostatically cycled lithium half cells of MNO-800 exhibited additional capacity (307 mAhg −1 at 0.1C, 4.66% higher) and improved high-rate capacity of 200 mAhg −1 at 10C. ICI-based overpotential analysis identified solid state diffusion as the dominant rate limiting process where MNO-800 correspondingly exhibited ∼3X faster capacity-weighted diffusivity. A machine-learning interatomic potential was trained to density functional theory and then applied with molecular dynamics (MLIP-MD) to examine the possible roles of Wadsley defects and transition metal disorder. For both defect-types, Li was found to populate and activate fast diffusion paths from window sites at lower extents of lithiation as compared to the order-rich model.

defect↗
NASA NTRS1989

Effects of H2O, CO2, and N2 air contaminants on critical airside strain rates for extinction of hydrogen-air counterflow diffusion flames

Dish-shaped counterflow diffusion flames centered by opposing laminar jets of H2 and clean and contaminant O2/N2 mixtures in an argon bath at 1 atm were used to study the effects of contaminants on critical airside strain. The jet velocities for both flame extinction and restoration are found for a wide range of contaminant and O2 concentrations in the air jet. The tests are also conducted for a variety of input H2 concentrations. The results are compared with those from several other studies.

Pellett, G. L.↗