SEARCH · Search NASA

Results for “diffuse”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Kinetic Monte Carlo Simulation of Oxygen and Cation Diffusion in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia (YSZ) is of interest to the aerospace community, notably for its application as a thermal barrier coating for turbine engine components. In such an application, diffusion of both oxygen ions and cations is of concern. Oxygen diffusion can lead to deterioration of a coated part, and often necessitates an environmental barrier coating. Cation diffusion in YSZ is much slower than oxygen diffusion. However, such diffusion is a mechanism by which creep takes place, potentially affecting the mechanical integrity and phase stability of the coating. In other applications, the high oxygen diffusivity of YSZ is useful, and makes the material of interest for use as a solid-state electrolyte in fuel cells. The kinetic Monte Carlo (kMC) method offers a number of advantages compared with the more widely known molecular dynamics simulation method. In particular, kMC is much more efficient for the study of processes, such as diffusion, that involve infrequent events. We describe the results of kinetic Monte Carlo computer simulations of oxygen and cation diffusion in YSZ. Using diffusive energy barriers from ab initio calculations and from the literature, we present results on the temperature dependence of oxygen and cation diffusivity, and on the dependence of the diffusivities on yttria concentration and oxygen sublattice vacancy concentration. We also present results of the effect on diffusivity of oxygen vacancies in the vicinity of the barrier cations that determine the oxygen diffusion energy barriers.

Good, Brian↗

NASA NTRS1962

Performance Evaluation of Fixed- and Variable-Area Rocket Exhaust Diffusers Using Single and Clustered Nozzles with and Without Gimbaling

An investigation of exhaust diffusers used for altitude simulation in testing rocket engines was conducted with model diffusers and gaseous nitrogen as the working fluid. This investigation was conducted to evaluate the effects on performance of gimbaling clustered nozzles in a single fixed-area exhaust diffuser. A two-nozzle cluster was gimbaled in all attitudes in four different exhaust-diffuser configurations: (1) a straight circular tube, (2) a straight figure-eight tube, (3) a circular tube with a second throat, and (4) a figure-eight tube with a second throat. Gimbaling clustered nozzles had little effect on the operating pressure ratio for either straight-tube or second-throat diffusers, but caused a large increase in the starting-pressure-ratio requirements in some straight-tube exhaust diffusers. The performance of two or four clustered nozzles with no gimbaling was compared with single-nozzle performance on the basis of the ratio of diffuser area to nozzle-throat area. The performance was found to be dependent on the diffuser- to nozzle-throat-area ratio and independent of both the number of primary nozzles and the nozzle-area ratio. A separate investigation was also conducted to evaluate the performance improvement obtainable with a variable-area exhaust diffuser. This type of diffuser achieved a 26-percent reduction in starting pressure ratio, a 40-percent reduction in operating pressure ratio, and a 50-percent reduction in the overall diffuser length over the values attainable with a fixed-area second-throat diffuser.

ROCKET↗

NASA NTRS1989

Tracer versus trace element diffusion - Diffusional decoupling of Sr concentration from Sr isotope composition

Chemical diffusivities for nonalkali major, minor, and selected trace elements have been measured at 10 kbar, 1300 C during interdiffusion of dacitic and rhyolitic melts. Sr tracer diffusivities were measured at the same conditions in melts of dacitic and rhyolitic composition and in a diffusion couple of these two melts undergoing simultaneous chemical diffusion. Chemical diffusivities of all elements studied are functions of the silica content of the melt, but at any one silica content the chemical diffusivities of all elements studied are within a factor of five of the Si and Al diffusivities. Si diffusivity varies from 2.8 x 10 to the -9th sq cm/s at 65 percent SiO2 to 5.1 x 10 to the -10th sq cm/s at 75 percent SiO2. The diffusivities of minor and trace elements are proposed to be controlled by the ability of the local melt structure to accommodate these elements while maintaining LTE. The measured Sr tracer diffusivities are an order of magnitude above Sr chemical diffusivities.

Baker, Don R.↗

NASA NTRS2002

Static and Dynamic Effects of Lateral Carrier Diffusion in Semiconductor Lasers

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. It is well known that the value of diffusion coefficients affects the threshold pumping current of a semiconductor laser. At the same time, the strength of carrier diffusion process is expected to affect the modulation bandwidth of an AC-modulated laser. It is important not only to investigate the combined DC and AC effects due to carrier diffusion, but also to separate the AC effects from that of the combined effects in order to provide design insights for high speed modulation. In this presentation, we apply a hydrodynamic model developed by the present authors recently from the semiconductor Bloch equations. The model allows microscopic calculation of the lateral carrier diffusion coefficient, which is a nonlinear function of the carrier density and plasma temperature. We first studied combined AC and DC effects of lateral carrier diffusion by studying the bandwidth dependence on diffusion coefficient at a given DC current under small signal modulation. The results show an increase of modulation bandwidth with decrease in the diffusion coefficient. We simultaneously studied the effects of nonlinearity in the diffusion coefficient. To clearly identify how much of the bandwidth increase is a result of decrease in the threshold pumping current for smaller diffusion coefficient, thus an effective increase of DC pumping, we study the bandwidth dependence on diffusion coefficient at a given relative pumping. A detailed comparison of the two cases will be presented.

Li, Jian-Zhong↗

NASA NTRS2002

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗

NASA NTRS2016

Passive Rocket Diffuser Testing: Reacting Flow Performance of Four Second-Throat Geometries

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure. As one of the nation's largest rocket testing facilities, the performance and design limitations of diffusers are of great interest to NASA's Stennis Space Center. This paper describes a series of tests conducted on four diffuser configurations to better understand the effects of inlet geometry and throat area on starting behavior and boundary layer separation. The diffusers were tested for a duration of five seconds with a 1455-pound thrust, LO2/GH2 thruster to ensure they each reached aerodynamic steady state. The effects of a water spray ring at the diffuser exits and a water-cooled deflector plate were also evaluated. Static pressure and temperature measurements were taken at multiple axial locations along the diffusers, and Computational Fluid Dynamics (CFD) simulations were used as a tool to aid in the interpretation of data. The hot combustion products were confirmed to enable the diffuser start condition with tighter second throats than predicted by historical cold-flow data or the theoretical normal shock method. Both aerodynamic performance and heat transfer were found to increase with smaller diffuser throats. Spray ring and deflector cooling water had negligible impacts on diffuser boundary layer separation. CFD was found to accurately capture diffuser shock structures and full-flowing diffuser wall pressures, and the qualitative behavior of heat transfer. However, the ability to predict boundary layer separated flows was not consistent.

Jones, Daniel R.↗

NASA NTRS2016

Banded Structures in Electron Pitch Angle Diffusion Coefficients from Resonant Wave Particle Interactions

Electron pitch angle (D (alpha)) and momentum (D(pp)) diffusion coefficients have been calculated due to resonant interactions with electrostatic electron cyclotron harmonic (ECH) and whistler mode chorus waves. Calculations have been performed at two spatial locations L = 4.6 and 6.8 for electron energies 10 keV. Landau (n = 0) resonance and cyclotron harmonic resonances n = +/-1, +/-2,...+/-5 have been included in the calculations. It is found that diffusion coefficient versus pitch angle (alpha) profiles show large dips and oscillations or banded structures. The structures are more pronounced for ECH and lower band chorus (LBC) and particularly at location 4.6. Calculations of diffusion coefficients have also been performed for individual resonances. It is noticed that the main contribution of ECH waves in pitch angle diffusion coefficient is due to resonances n = +1 and n = +2. A major contribution to momentum diffusion coefficients appears from n = +2. However, the banded structures in D alpha and Dpp coefficients appear only in the profile of diffusion coefficients for n = +2. The contribution of other resonances to diffusion coefficients is found to be, in general, quite small or even negligible. For LBC and upper band chorus waves, the banded structures appear only in Landau resonance. The Dpp diffusion coefficient for ECH waves is one to two orders smaller than D alpha coefficients. For chorus waves, Dpp coefficients are about an order of magnitude smaller than D alpha coefficients for the case n does not = 0. In case of Landau resonance, the values of Dpp coefficient are generally larger than the values of D alpha coefficients particularly at lower energies. As an aid to the interpretation of results, we have also determined the resonant frequencies. For ECH waves, resonant frequencies have been estimated for wave normal angle 89 deg and harmonic resonances n = +1, +2, and +3, whereas for whistler mode waves, the frequencies have been calculated for angle 10 deg and Landau resonance. Further, in ECH waves, the banded structures appear for electron energies (is) greater than1 keV, and for whistler mode chorus waves, structures appear for energies greater than 2 keV at L = 4.6 and above 200 eV for L = 6.8. The results obtained in the present work will be helpful in the study of diffusion curves and will have important consequences for diffuse aurora and pancake distributions.

diffusion↗

NASA NTRS1999

Diffusion Barriers to Increase the Oxidative Life of Overlay Coatings

Currently, most blades and vanes in the hottest section of aero gas turbine engines require some type of coating for oxidation protection. Newly developed single crystal superalloys have the mechanical potential to operate at increasingly higher component temperatures. However, at these elevated temperatures, coating/substrate interdiffusion can shorten the protective life of the coating. Diffusion barriers between overlay coatings and substrates are being examined to extend the protective life of the coating. A previously- developed finite-difference diffusion model has been modified to predict the oxidative life enhancement due to use of a diffusion barrier. The original diffusion model, designated COSIM, simulates Al diffusion in the coating to the growing oxide scale as well as Al diffusion into the substrate. The COSIM model incorporates an oxide growth and spalling model to provide the rate of Al consumption during cyclic oxidation. Coating failure is predicted when the Al concentration at the coating surface drops to a defined critical level. The modified COSIM model predicts the oxidative life of an overlay coating when a diffusion barrier is present eliminating diffusion of Al from the coating into the substrate. Both the original and the modified diffusion models have been used to predict the effectiveness of a diffusion barrier in extending the protective life of a NiCrAl overlay coating undergoing cyclic oxidation at 1100 C.

Nesbitt, James A.↗

NASA NTRS1999

Planar Strain-Rate-Free Diffusion Flames: Initiation, Properties, and Extinction

An effectively strain-rate-free diffusion flame constitutes the most vigorous laminar combustion of initially unmixed reactive gases. Such a diffusion flame is characterized by a relatively long residence time and by a relatively large characteristic length scale. If such a flame were also planar, providing high symmetry, it would be particularly suitable for experimental and theoretical investigations of key combustion phenomena, such as multicomponent diffusion, chemical kinetics, and soot inception, growth, and oxidation. Unfortunately, a planar strain-rate-free diffusion flame is highly disrupted in earth-gravity (e.g., in a counterflow-diffusion-flame apparatus) because of the very rapid onset (approx. 100 ms) of gravity-induced instability. Accordingly, a specially dedicated apparatus was designed, fabricated, and initially checked out for the examination of a planar strain-rate-free diffusion flame in microgravity. Such a diffusion flame may be formed within a hollowed-out squat container (initially configured as 25 cm x 25 cm x 9 cm), with isothermal, noncatalytic, impervious walls. At test initiation, a thin metallic sheet (approx. 1 mm in thickness) that separates the internal volume into two equal portions, each of dimensions 25 cm x 25 cm x 4.5 cm, is withdrawn, by uniform translation (approx. 50 cm/s) in its own plane, through a tightly fitting slit in one side wall. Thereupon, diluted fuel vapor (initially confined to one half-volume of the container) gains access to diluted oxygen (initially with the same pressure, density, and temperature as the fuel, but initially confined to the other half-volume). After a brief delay (approx. 10 ms), to permit limited but sufficient-for-flammability diffusional interpenetration of fuel vapor and oxidizer, burning is initiated by discharge of a line igniter, located along that side wall from which the trailing edge of the separator withdraws. The ignition spawns a triple-flame propagation across the 25 cm x 25 cm centerplane. When a diffusion flame is emplaced in the centerplane, any subsequent travel, and change in temperature, of that planar diffusion flame may be tracked, along with the effectively spatially uniform but temporally evolving pressure within the container. Eventually, nearly complete depletion of the stoichiometrically deficient reactant, along with heat loss to the container surfaces, effects extinction. These data afford an opportunity to check theoretical models of diffusion and chemical kinetics under conditions ranging from intense burning to flame out, or, alternatively, to evolve simple empirical representations of these phenomena. Thus, the project sought to utilize microgravity testing to elucidate commonly encountered phenomenology, arising in the commonly-encountered mode of combustion (whether related to heating, manufacturing, boiling, and propulsion, or to uncontrolled, free-burning fire in structures and wildland vegetation), of those commonly utilized fuels usually categorized as gaseous fuels (such as hydrogen, natural gas, and propane, which are gaseous under atmospheric conditions).

Fendell, Francis↗

NASA NTRS2001

Reaction Kernel Structure of a Slot Jet Diffusion Flame in Microgravity

Diffusion flame stabilization in normal earth gravity (1 g) has long been a fundamental research subject in combustion. Local flame-flow phenomena, including heat and species transport and chemical reactions, around the flame base in the vicinity of condensed surfaces control flame stabilization and fire spreading processes. Therefore, gravity plays an important role in the subject topic because buoyancy induces flow in the flame zone, thus increasing the convective (and diffusive) oxygen transport into the flame zone and, in turn, reaction rates. Recent computations show that a peak reactivity (heat-release or oxygen-consumption rate) spot, or reaction kernel, is formed in the flame base by back-diffusion and reactions of radical species in the incoming oxygen-abundant flow at relatively low temperatures (about 1550 K). Quasi-linear correlations were found between the peak heat-release or oxygen-consumption rate and the velocity at the reaction kernel for cases including both jet and flat-plate diffusion flames in airflow. The reaction kernel provides a stationary ignition source to incoming reactants, sustains combustion, and thus stabilizes the trailing diffusion flame. In a quiescent microgravity environment, no buoyancy-induced flow exits and thus purely diffusive transport controls the reaction rates. Flame stabilization mechanisms in such purely diffusion-controlled regime remain largely unstudied. Therefore, it will be a rigorous test for the reaction kernel correlation if it can be extended toward zero velocity conditions in the purely diffusion-controlled regime. The objectives of this study are to reveal the structure of the flame-stabilizing region of a two-dimensional (2D) laminar jet diffusion flame in microgravity and develop a unified diffusion flame stabilization mechanism. This paper reports the recent progress in the computation and experiment performed in microgravity.

Takahashi, F.↗

NASA NTRS1946

Performance of the Modified V-1710-93 Engine-Stage Supercharger with a Constant-Area Vaneless Diffuser

As part of an investigation to increase the power output of the V-1710-93 engine at altitude, the engine-stage supercharger was combined with a constant-area vaneless diffuser designed to improve the performance of the engine-stage supercharger at the rated engine operating point. The performance of the modified supercharger was investigated in a variable-component supercharger test rig and compared with that of the standard supercharger with an 8-vaned diffuser. A separate evaluation of the component efficiencies and a study of the flow characteristics of the modified supercharger was made possible by internal diffuser instrumentation. At the volume flow required by the engine for rated operating conditions, the modified supercharger increased the over-all adiabatic efficiency 0.05 and the over-all pressure coefficient 0.035. Furthermore, the capacity of the engine-stage supercharger was increased by replacing the standard 8-vaned diffuser with the vaneless diffuser. The peak over-all adiabatic efficiency for the modified supercharger, however, was 0.05 to 0.07 lower than that of the standard unit over the range of tip speeds investigated. The improved performance of the modified supercharger at rated engine operating conditions resulted from a shift of the point of peak adiabatic efficiency and pressure coefficient of the standard supercharger to a higher flow. The energy loss through the vaneless diffuser was found to be small. Because of the restricted diffuser diameter, however, diffusion was inadequate, which resulted in a relatively small static-pressure rise through the diffuser, high diffuser-exit velocities, and excessive collector-case losses.

Douglas, John E.↗

NASA NTRS2002

Centrifugal Compressor Surge Margin Improved With Diffuser Hub Surface Air Injection

Aerodynamic stability is an important parameter in the design of compressors for aircraft gas turbine engines. Compression system instabilities can cause compressor surge, which may lead to the loss of an aircraft. As a result, engine designers include a margin of safety between the operating line of the engine and the stability limit line of the compressor. The margin of safety is typically referred to as "surge margin." Achieving the highest possible level of surge margin while meeting design point performance objectives is the goal of the compressor designer. However, performance goals often must be compromised in order to achieve adequate levels of surge margin. Techniques to improve surge margin will permit more aggressive compressor designs. Centrifugal compressor surge margin improvement was demonstrated at the NASA Glenn Research Center by injecting air into the vaned diffuser of a 4:1-pressure-ratio centrifugal compressor. Tests were performed using injector nozzles located on the diffuser hub surface of a vane-island diffuser in the vaneless region between the impeller trailing edge and the diffuser-vane leading edge. The nozzle flow path and discharge shape were designed to produce an air stream that remained tangent to the hub surface as it traveled into the diffuser passage. Injector nozzles were located near the leading edge of 23 of the 24 diffuser vanes. One passage did not contain an injector so that instrumentation located in that passage would be preserved. Several orientations of the injected stream relative to the diffuser vane leading edge were tested over a range of injected flow rates. Only steady flow (nonpulsed) air injection was tested. At 100 percent of the design speed, a 15-percent improvement in the baseline surge margin was achieved with a nozzle orientation that produced a jet that was bisected by the diffuser vane leading edge. Other orientations also improved the baseline surge margin. Tests were conducted at speeds below the design speed, and similar results were obtained. In most cases, the greatest improvement in surge margin occurred at fairly low levels of injected flow rate. Externally supplied injection air was used in these experiments. However, the injected flow rates that provided the greatest benefit could be produced using injection air that is recirculating between the diffuser discharge and nozzles located in the diffuser vaneless region. Future experiments will evaluate the effectiveness of recirculating air injection.

Skoch, Gary J.↗

NASA NTRS2010

Vertical Diffusivities of Active and Passive Tracers

The climate models that include a carbon-cycle need the vertical diffusivity of a passive tracer. Since an expression for the latter is not available, it has been common practice to identify it with that of salt. The identification is questionable since T, S are active, not passive tracers. We present the first derivation of the diffusivity of a passive tracer in terms of Ri (Richardson number) and Rq (density ratio, ratio of salinity over temperature z-gradients). The following results have emerged: (a) The passive tracer diffusivity is an algebraic function of Ri, Rq. (b) In doubly stable regimes (DS, partial derivative of T with respect to z > 0, partial derivative of S with respect to z < 0), the passive scalar diffusivity is nearly the same as that of salt/heat for any values of Rq < 0 and Ri > 0. (c) In DC regimes (diffusive convection, partial derivative of T with respect to z < 0, partial derivative of S with respect to z < 0, Rq > 1), the passive scalar diffusivity is larger than that of salt. At Ri = O(1), it can be more than twice as large. (d) In SF regimes (salt fingers, partial derivative of T with respect to z > 0, partial derivative of S with respect to z > 0, Rq < 1), the passive scalar diffusivity is smaller than that of salt. At Ri = O(1), it can be less than half of it. (e) The passive tracer diffusivity predicted at the location of NATRE (North Atlantic Tracer Release Experiment) is discussed. (f) Perhaps the most relevant conclusion is that the common identification of the tracer diffusivity with that of salt is valid only in DS regimes. In the Southern Ocean, where there is the largest CO2 absorption, the dominant regime is diffusive convection discussed in (c) above.

Canuto, V. M.↗

NASA NTRS2017

Diffusion of Redox-Sensitive Elements in Basalt at Different Oxygen Fugacities

The terrestrial planets and moons of our solar system have differentiated over a range of oxygen fugacity conditions. Basalts formed from magmas on the Earth cover a range of more oxidized states (from approximately IW (iron wustite) plus 2 to approximately FMQ (fayalite-magnetite-quartz) plus 3) than crustal rocks from Mars (IW to approximately IW plus 3), and basalts from the Moon are more reduced than both, ranging from IW to IW minus 2. The small body Vesta differentiated around IW minus 4. Characterization of redox sensitive elements' diffusivities will offer insight into behavior of these elements as a function of f (fugacity of) O2 for these planetary bodies. Here, we report a systematic study of the diffusion of redox-sensitive elements in basaltic melts with varying oxygen fugacities (fO2) for trace elements, V, Nb, W, Mo, La, Ce, Pr, Sm, Eu, Gd, Ta, and W. Since fO2 is an intensive variable that is different for the reservoirs of various planets and moons in our solar system, it is important to characterize how changes in redox states will affect diffusion. We conducted experiments in a piston cylinder device at 1300 degrees Centigrade and 1 gigapascal, at the University of Rochester and NASA Johnson Space Center. We buffered some experiments at Ru-RuO2 (FMQ plus 6.00), and conducted other experiments within either a graphite or Mo capsule, which corresponds to fO2s of either FMQ minus1.2, or FMQ minus 3.00, respectively. Characterizing the diffusivities of redox sensitive elements at different fO2s is important because some elements, like Eu, have varying valence states, such as Eu (sup 2 plus) and Eu (sup 3 plus). Differences in charge and ion radii may lead to differences in diffusivities within silicate melts. This could, lead to formation of a Eu anomaly by diffusion, the magnitude of which may be controlled by the fO2. Characterization of trace element diffusion is also important in understanding trace element fractionation. We found, during the course of our investigation, that not only did the diffusivities of the redox sensitive elements change with fO2, but that the diffusivities of all other analyzed elements also changed. This indicates that not only do changes in valence influence trace elements diffusivities but that the structure of melt may have changed with varying oxygen fugacity, probably due to changes in the speciation of the major element Fe.

Szumila, I.↗

NASA NTRSDate not supplied

Thermodynamic and Diffusion Model Estimates on Metamorphic Temperatures and Timescales for Basaltic Eucrite GRA 98098

Introduction: HED meteorites are thought to rep-resent igneous rocks from Vesta’s basaltic crust and preserve evidence of early crustal metamorphism. Determining the temperatures and timescales of thermal metamorphism is important for reconstructing crustal evolution in the early solar system. Here, we study basaltic eucrite Graves Nunataks (GRA) 98098, which has been identified as a highly metamorphosed eucrite [1]. We present new estimates on metamorphic temperatures determined via thermodynamic modeling as well as the initial results from diffusion models constraining timescales of thermal metamorphism. Sample Description: GRA 98098 is an unbrecciated eucrite with a granoblastic plagioclase and pyroxene mineralogy. Millimeter to cm-long lathes of tridymite cross-cut and poikilitically enclose plagioclase and pyroxene [1,this work]. Pyroxene grains have exsolved into Ca-rich (~Wo38En29Fs33) and Ca-poor (~Wo4.5En36Fs59.5) lamellae. Both unzoned and zoned plagioclase grains are observed. Unzoned plagioclase grains are found solely with tridymite laths. These grains have ~An92 compositions. The cores of the zoned plagioclase grains have the same composition and thin, relatively sodic rims (~An67), (Fig. 1). The bulk sample is unusually enriched in highly in-compatible elements and has one of the most fractionated REE patterns reported [1]. Maximum metamorphic temperatures of 985±78°C have been estimated using two-pyroxene thermometry [2]. Methods: Thermodynamic modeling. Thermodynamic models were constructed using the software Perple_X, which employs a Gibbs free energy minimization in order to determine the most stable phase assemblage for a given bulk rock composition [3]. Bulk composition was calculated using mineral com-positions acquired via EMPA (this study) and the observed abundancies present in the thin section. Two bulk compositions were estimated; 1) includes all phases present in the thin section, (assumes that all phases are present during metamorphism), 2) excludes tridymite from the bulk calculation (assumes that tridymite was not present during metamorphism). In order to determine whether metamorphic equilibria was achieved and estimate temperatures of metamorphism, we compared measured pyroxene compositions with thermodynamically predicted compositions [4, Fig. 2]. Diffusion Modeling. Several time-temperature de-pendent diffusion profiles were calculated in order to determine the best match for XAn chemical profiles observed at the edges of the zoned plagioclase (Fig. 3). We assumed that the start condition was a stepwise gradient at the plagioclase/pyroxene interface. We also assumed an average diffusion coefficient (D) and a constant temperature using the equation in [5]. D was determined for two temperatures (T = 1060ºC; near eucrite solidus [6] and T = 985ºC; metamorphism reported in [4]) and then XAn was calculated as a function of distance from plagioclase core to rim using an error function solution to Fick’s second law. Results: Thermodynamic model results are summarized in Fig. 2. For a bulk composition that includes all phases in the thin section, pyroxene endmember compositions plot in the following temperature ranges: Fs ~660-860ºC, En~1000ºC & 1150ºC, and Wo~760-900ºC (Fig. 2a). For a bulk composition that excludes tridymite from the peak metamorphic assemblage (i.e., the bulk composition minus the contribution from tridymite), a temperature range could not be determined for the Fs component of pyroxene. For Wo, T~760-900ºC and En, T~1000ºC & 1150ºC (Fig. 2b). Fig. 3 summarizes the diffusion model results. For T = 1060°C & 985°C, the most appropriate time interval was estimated based on which diffusion curve most matched (solid lines, Fig. 3) the EMPA data. For T = 1060°C, the best looking match was t = 500 ka. For T = 985°C, the best match was t = 7 Ma. Discussion and future work: Temperature estimates from thermodynamic models are not conclusive because the temperature ranges determined for pyroxene endmember stability do not overlap (colored fields in Fig. 2), thus implying that there is disequilibrium between pyroxene crystals and the bulk composition considered [4]. Thus, additional exploration is needed to define a metamorphically equilibrated do-main that accurately records peak temperature. The utility of defining metamorphically equilibrated do-mains to improve the accuracy and level of detail elucidated regarding the petrogenetic history of metamorphose samples has been demonstrated previously [4,7]. We suggest that in the case of Fig. 2a, the thin section composition is not representative of the length scales over which metamorphic equilibrium was achieve and in the case of Fig. 2b, the assumption that tridymite was not present during metamorphism was incorrect. However, results from thermodynamic models can provide insight into the relative timing of mineral and compositional textures. For example from texture alone, it is unclear whether tridymite was igneous in origin and represents the last bits of melt in a crystallizing magma chamber, or if it formed during (and possibly initiated) open system thermal metamorphism. The latter could be consistent with a partial melt hypothesis [8,9] while the former implies that simple fractional crystallization can yield the textures present in GRA 98098. The lack of coincidence be-tween pyroxene endmember compositions in Fig. 2b suggest that the bulk composition minus tridymite was not the assemblage in equilibrium with the pyroxene, suggesting that tridymite was present during metamorphism and formed during igneous crystallization. We conclude that the development of the Na-rich plagioclase rims likely occurred during or immediately after peak thermal metamorphism, because eucrites of similar metamorphic grade and texture have unzoned plagioclase (~An92) [2,4,8], and Na zoning is only observed in the plagioclase not included in the tridymite. This suggests that the zoning formed after tridymite formation, and therefore after igneous crystallization. Thus, the timescales calculated via diffusion modeling possibly represent the time interval over which thermal metamorphism occurred. Cooling rates approximated for the Vestan crust predict that the crust cooled below 300°C around 35-40 Ma after formation[10]. This is consistent with our modeling results that predict formation of the Na rich plagioclase rims occurring at higher temperatures over a period of 0.5 to 7 Ma years. Future work. Additional thermodynamic modeling work will focus on selecting an equilibrated bulk rock domain in which to elucidate metamorphic conditions. Diffusion models currently provide a minimum time-scale, since diffusion slows down as the system cools. Future work with will focus on integrating cooling into the diffusion models and constraining the depth at which thermal metamorphism occurs because it could be used to determine whether the range of time-scales calculated for thermal metamorphism are consistent with the geologic environment.

J S Gorce↗

NASA NTRS2023

Passive Rocket Diffuser Testing: Reacting Flow Performance of Advanced Configurations

Upper stage rocket engines are optimized for operation at the edge of space. Ground testing often requires simulation of high-altitude environments to prevent damage from off-design operation. Passive supersonic diffusers are an effective means of establishing the partial vacuum needed around the nozzle exit. The low cost of diffusers relative to active vacuum sources (ejectors, pumps, etc.) provides a perpetual incentive to improve their performance and expand the envelope of passively testable engines. Second-throat diffusers have been the de facto standard for rocket testing since the 1950s, and numerous studies have focused on parametric optimization of the topology. Yet, even the best second-throat designs may require more driving pressure than an upper stage engine can provide. To overcome conventional limits, engineers at NASA’s Stennis Space Center (SSC) adapted an obscure cold-flow diffuser topology to the extremes of immersion in supersonic combustion products. Aerodynamic experimentation was conducted at SSC’s E-3 test stand to quantify the performance of prototypical hot-fire spike diffuser hardware and enable direct comparison to prior states of the art. A LO X /GH 2 thruster was contoured to approximate the regeneratively-cooled nozzle of an upper-stage engine and fired into 28 downstream diffuser configurations: 16 second-throat, 2 centerbody, and 10 spike. Temperature and pressure were recorded along the test article walls for a range of chamber pressures spanning 3.5 - 5.0 MPa. The second-throat configurations with the lowest start and unstart pressure ratios were shown to outperform equivalent systems in literature and taken as baselines for evaluation of advanced diffuser performance. Compared to these exemplars, centerbody diffusers decreased start pressure ratio by 11% but increased unstart pressure ratio by 15%. Spike diffusers reduced start and unstart pressure ratios by 24% and 9%, respectively.

Rocket↗

NASA NTRS1972

Diffusion of chromium and aluminum in Ni-20Cr and TDNiCr /Ni-20Cr-2ThO2/.

Diffusion coefficients have been measured for Cr51 in fine- and coarse-grained TDNiCr (Ni-20Cr-2ThO2) and in fine-grained Ni-20Cr in the temperature range from 1038 to 1200 C. Selective diffusivities have also been determined for specimens of these alloys which were aluminized to give an initial surface concentration of 5.8 wt % Al. Finally, diffusion coefficients for interdiffusion of aluminum in TDNiCr and Ni-20Cr have been obtained from electron probe microanalysis of the aluminized specimens. For a given grain size and temperature there is no difference in diffusivities for chromium diffusion in TDNiCr or Ni-20Cr. Diffusion coefficients increase with decreasing grain size for both alloys. Comparison of aluminum diffusion data obtained from electron microprobe profiles with radiotracer chromium diffusivities suggests that aluminum diffuses approximately three times faster than chromium in TDNiCr and Ni-20Cr.

Seltzer, M. S.↗

NASA NTRS1973

Directional diffusion coefficients of solar protons inside and outside the bow shock.

The directional diffusion coefficients of low-energy (greater than or equal to 0.3 MeV) solar protons inside and outside the bow shock are examined during the solar flare event of Jan. 24, 1969. The data are derived from simultaneous observations obtained by Explorer 33 inside the magnetosheath and by Explorer 35 in the interplanetary medium. Although the gross properties of the spin-averaged intensities on a diffusion-type plot appear to be the same in both media, the directional intensities show significant variations. It is shown that directional intensities of low-energy protons can be described reasonably well by anisotropic diffusion with an associated diffusion coefficient. Directional diffusion coefficients are found to differ by a factor of as much as three among different directions in space, and from the spin-averaged diffusion coefficient. This suggests that anisotropic diffusion does indeed take place and that so called 'isotropic' diffusion coefficients derived in the past from spin-averaged intensities may actually be directional diffusion coefficients in cases where substantial anisotropies (greater than 50%) exist.

Verzariu, P.↗