SEARCH · Engineering Papers
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.
Tunable far-infrared spectroscopy was used to measure pressure broadening of the 118.455/cm rotational lines of OH by H2, He, N2, and O2. Broadening coefficients of H2 and He were measured for the temperature range of 296 to 80 K, and N2 and O2 broadening coefficients were measured for the 296-210 K temperature range.
Using first-principles approaches we find that the Auger recombination in PbSe is anomalous in three distinct ways. First, the direct Auger coefficient is 4 orders of magnitude lower than that of other semiconductors with similar band gaps, a result that can be attributed to the lack of involvement of a heavy-hole band. Second, phonon-assisted indirect Auger recombination prevails, contrary to the common belief that direct Auger is dominant in narrow-gap semiconductors. Third, an unexpectedly weak temperature dependence of the Auger coefficient is observed, which we can now attribute to the indirect nature of the Auger process. The widely accepted explanation of this behavior in terms of an unusual temperature dependence of the band gap is only a secondary effect. Furthermore, our results elucidate the mechanisms underlying the anomalous Auger recombination in IV-VI semiconductors in general, which is critical for understanding and engineering carrier transport.
Densities derived from the drag of five artificial satellites and covering the time interval from 1958 to 1971 have been used to study the variations in the thermosphere arising from changes in solar activity. Exospheric temperatures derived from the densities by means of atmospheric models have been related to the 10.7-cm solar flux, in which the active-area component had been separated from the slow-varying clear-disk component by means of a relation between the flux and sunspot areas. When the day-to-day (27-day) variation is suppressed, the remaining slow temperature variation appears to be linearly related to the clear-disk flux component. If we also relate linearly the day-to-day temperature variations with the day-to-day variations of the flux, we find that the regression coefficient changes from 1 K at sunspot maximum to 2 K at low solar activity; a nearly constant coefficient is obtained by relating the temperature variations to the ratio of the active-area component to the clear-disk component of the flux.
The Seawinds scatterometer on the advanced Earth observing satellite-2 (ADEOS-2) will determine surface wind vectors by measuring the radar cross section. Multiple measurements will be made at different points in a wind-vector cell. When dense clouds and rain are present, the signal will be attenuated, thereby giving erroneous results for the wind. This report describes algorithms to use with the advanced mechanically scanned radiometer (AMSR) scanning radiometer on ADEOS-2 to correct for the attenuation. One can determine attenuation from a radiometer measurement based on the excess brightness temperature measured. This is the difference between the total measured brightness temperature and the contribution from surface emission. A major problem that the algorithm must address is determining the surface contribution. Two basic approaches were developed for this, one using the scattering coefficient measured along with the brightness temperature, and the other using the brightness temperature alone. For both methods, best results will occur if the wind from the preceding wind-vector cell can be used as an input to the algorithm. In the method based on the scattering coefficient, we need the wind direction from the preceding cell. In the method using brightness temperature alone, we need the wind speed from the preceding cell. If neither is available, the algorithm can work, but the corrections will be less accurate. Both correction methods require iterative solutions. Simulations show that the algorithms make significant improvements in the measured scattering coefficient and thus is the retrieved wind vector. For stratiform rains, the errors without correction can be quite large, so the correction makes a major improvement. For systems of separated convective cells, the initial error is smaller and the correction, although about the same percentage, has a smaller effect.
A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). A critical phenomenon of interest with U-Mo monolithic fuel is the large amount of swelling that takes place during operation, particularly at high fission densities. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic and kinetic properties into mesoscale and continuum level fuel performance modeling codes. One such property where there exists incomplete data is the diffusion of relevant species under irradiation. Fuel performance swelling predictions rely on an accurate representation of diffusion in order to determine the rate of fission gas swelling and the local microstructural evolution. In this work, we present molecular dynamics simulations of the radiation driven diffusion of U, Mo and Xe in U-Mo nuclear fuels. Diffusion coefficients for each species are determined over a range of temperatures and compositions. In this work, updated diffusion coefficients are presented that are applicable under irradiation that incorporate both intrinsic and radiation driven diffusion.
A theoretical investigation was conducted to determine optimum experimental conditions for using a lumped-parameter calorimeter to measure heat-transfer coefficients and heating rates. A mathematical model of the transient temperature response of the calorimeter was used with the measured temperature response to predict the heat-transfer coefficient and the rate of heating. A sensitivity analysis was used to determine the optimum transient experiment for simultaneously measuring the heat addition during heating and the convective heat-transfer coefficient during heating and cooling of a lumped-parameter calorimeter. Optimum experiments were also designed for measuring the convective heat-transfer coefficient during both heating and cooling and cooling only.
The Mars Odyssey spacecraft made use of multipass aerobraking to gradually reduce its orbit period from a highly elliptical insertion orbit to its final science orbit. Aerobraking operations provided an opportunity to apply advanced thermal analysis techniques to predict the temperature of the spacecraft's solar array for each drag pass. Odyssey telemetry data was used to correlate the thermal model. The thermal analysis was tightly coupled to the flight mechanics, aerodynamics, and atmospheric modeling efforts being performed during operations. Specifically, the thermal analysis predictions required a calculation of the spacecraft's velocity relative to the atmosphere, a prediction of the atmospheric density, and a prediction of the heat transfer coefficients due to aerodynamic heating. Temperature correlations were performed by comparing predicted temperatures of the thermocouples to the actual thermocouple readings from the spacecraft. Time histories of the spacecraft relative velocity, atmospheric density, and heat transfer coefficients, calculated using flight accelerometer and quaternion data, were used to calculate the aerodynamic heating. During aerobraking operations, the correlations were used to continually update the thermal model, thus increasing confidence in the predictions. This paper describes the thermal analysis that was performed and presents the correlations to the flight data.
The reaction coefficient for hydrogen/deuterium (H/D) exchange and the diffusion of hydrated excess protons within amorphous solid water (ASW) are characterized as a function of temperature. For these experiments, water films are deposited on a Pt(111) substrate at 108 K, and reactions with pre-adsorbed hydrogen atoms produce hydrated protons. Upon heating, protons diffuse within the water, and H/D exchange occurs when they encounter D2O probe molecules deposited in the films. The time-dependent concentration of D2O is monitored with infrared spectroscopy, and it indicates the protons diffusion from the substrate and establish an equilibrium distribution prior to significant H/D exchange for temperatures 114 K ≤T≤ 134 K. By controlling the distance between the D2O molecules and the substrate, we probe the distribution of protons within the film. It decays as x−2 for the examined range of x (12–52 nm) due to the electric field that develops between the diffusing protons and their image charges in the metal substrate. This agrees with the theoretical distance scaling for the equilibrated proton concentration in a dielectric near a metal boundary. From the proton concentration and the measured D2O decay rate, a lower bound for the proton diffusion coefficient ranging from 10−20 m2/s at 114 K to 10−18 m2/s at 134 K is estimated. The diffusion coefficient has an activation energy of 0.40 eV, which is comparable to energies reported for molecular translations and rotations of H2O, suggesting they may play a critical role in the proton diffusion mechanism within ASW.
The intensity of the received signal at Earth in the radio occultations of Titan is attenuated both by refractive defocusing and pressure-induced absorption from N2-N2 and CH4-N2 pairs. Because the absorption strength is different for the two sets of pairs, matching the retrieved absorptivity profile can in principle yield the vertical variation in gaseous methane in the troposphere. There are two factors that make this difficult. The first is the propagation of noise in the phase and amplitude of the received signal in the absorption retrieval. The phase data is first inverted to retrieve vertical profiles of refractivity, from which the refractive defocusing is calculated. This is then subtracted from the observed. intensity attenuation of the received signal to generate a profile of atmospheric absorption. The second problem is the uncertainty in the pressure-induced absorption coefficients. Laboratory data at radio wavelengths is only available near room temperature (see, e.g., [1] for N2-N2), and the extrapolation to the low temperatures in Titan's troposphere is not well established. Ab initio calculations by Borysow et al. [2, 3] provide absorption coefficients at low temperatures and long wavelengths, but their accuracy has come into question. We present examples from Cassini radio occultations of Titan to illustrate the difficulties. For methane mole fractions in the lower troposphere comparable to that inferred from the Huygens probe (approximately 0.05), it will be difficult to separate the contributions of N2-N2 collisions from those of N2-CH4, collisions to the retrieved absorption. However, higher concentrations of CH4 and/or a higher signal-to-noise ratio from a future uplink experiment could result in a successful separation of the two components. However, key to this are highly accurate estimates of the absorption from a combination of laboratory measurements at love temperatures and long wavelengths, and possibly improved theoretical calculations.
Resistance drift in amorphous Ge 2 Sb 2 Te 5 is experimentally characterized in melt-quenched line cells in the range of 300 K to 125 K and is observed to follow the previously reported power-law behavior with drift coefficients in the range of 0.07 to 0.11 in the dark, linearly decreasing with 1/kT. While these drift coefficients measured in the dark are similar to commonly observed drift coefficients (~0.1) at and above room temperature, measurements under light show a significantly lower drift coefficient (0.05 under illumination vs 0.09 in the dark at 150 K). Periodic on/off switching of light shows a sudden decrease/increase in resistance, attributed to photo-excited carriers, followed by a very slow response (~30 min at 150 K) attributed to contribution of electron traps and slow trap-to-trap charge exchanges. A device-level electronic model is used to relate these experimental findings to gradual charging of electron traps in amorphous Ge 2 Sb 2 Te 5 , which gives rise to growth of a potential barrier for holes in time and, hence, resistance drift.
Large quantities of high-temperature air are needed for work with hypersonic flight problems. At temperatures above 2500 degrees Reamur, where conventional heat exchangers have exceeded their material limits, regenerative pebble-bed exchangers may be used with high-temperature refractories. The design of such a heat exchanger requires the use of reliable heat-transfer coefficients for a packed bed. Considerable data are available on the subject, but they spread over two orders of magnitude at any one Reynolds number value. The facility from which the present data were obtained is used at the Lewis Research Center (NASA) for testing air-breathing engine components. The purpose of this work was to obtain heat-transfer data during the initial operation of the bed as a guide to the design of similar equipment. The facility was designed with a conservative estimate of the heat-transfer coefficient, and is shown schematically. Temperatures throughout the packing were measured continuously so that point values of the coefficient might be obtained.
Thin-film spectra of solid NH3 at a resolution of 1 per cm were used to determine its absorption coefficient over the range 50-7000 per cm. The thin films were formed inside a liquid N2 cooled dewar using a variety of substrates and dewar windows. The spectra were recorded with two Fourier spectrometers, one covering the range from 1 to 4 microns and the other from 2.6 to 200 microns. The thickness of the films was measured with a laser interference technique. The absorption coefficients were determined by application of Lambert's law and by a fitting procedure to the observed spectra using thin-film theory. Good agreement was found with the absorption coefficients recently determined by other investigators over a more restricted wavelength range. A metastable phase was observed near a temperature of 90 K and its absorption coefficient is reported. No other major spectral changes with temperature were noted for the range 88-120 K.
An approach for synthesizing buckling results for, and behavior of, thin balanced and unbalanced symmetric laminates that are subjected to uniform heating or cooling and elastically restrained against thermal expansion or contraction is presented. This approach uses a nondimensional analysis for infinitely long, flexurally anisotropic plates that are subjected to combined mechanical loads and is based on useful nondimensional parameters. In addition, stiffness-weighted laminate thermal-expansion parameters and compliance coefficients are derived that are used to determine critical temperatures in terms of physically intuitive mechanical-buckling coefficients. The effects of membrane orthotropy and membrane anisotropy are included in the general formulation. Many results are presented for some common laminates that are intended to facilitate a structural designer's transition to the use of generic buckling design curves. Several curves that illustrate the fundamental parameters used in the analysis are presented, for nine contemporary material systems, that provide physical insight into the buckling response in addition to providing useful design data. Examples are presented that demonstrate the use of generic design curves. The analysis approach and generic results indicate the effects and characteristics of elastically restrained laminate thermal expansion or contraction, membrane orthotropy and anisotropy, and flexural orthotropy and anisotropy in a very general and unifying manner.
An approach for synthesizing buckling results for thin balanced and unbalanced symmetric laminates that are subjected to uniform heating or cooling and elastically restrained against thermal expansion or contraction is presented. This approach uses a nondimensional analysis for infinitely long, flexural anisotropic plates that are subjected to combined mechanical loads. In addition, stiffness-weighted laminate thermal-expansion parameters and compliance coefficients are derived that are used to determine critical temperatures in terms of physically intuitive mechanical-buckling coefficients. Many results are presented for some common laminates that are intended to facilitate a structural designer s transition to the use of the generic buckling design curves. Several curves that illustrate the fundamental parameters used in the analysis are presented, for nine contemporary material systems, that provide physical insight into the buckling response in addition to providing useful design data. Examples are presented that demonstrate the use of the generic design curves.
Ab initio calculation of vibrational excitation and deexcitation cross sections of N2 by low energy resonant electron impact have been carried out. The calculation includes initial target state in v = 0-12 and final target state in v = 0-17. The calculated cross sections are found to agree with available experimental data to within reported experimental error. Effect of target rotation has been investigated. Comparison of vibrational excitation cross sections at J = 0, 50, and 150 shows that the major effect of target rotation is to lower the resonance energy. Vibrational excitation and deexcitation rate coefficients are tabulated for electron temperatures of 1100-55,000 K. The present tabulation supercedes an earlier compilation by the Naval Research Laboratory that was derived from experimental data with a factor-of-two error. The tabulated rate coefficients are suited for flowfield modeling of aeroassisted orbital transfer vehicles.
An interactive model which couples a semi-spectral dynamical model, a radiative transfer code and a two-dimensional chemistry transport model (2-D CTM), is used to assess the atmospheric effects of the High-Speed Civil Transport (HSCT) engine emissions. The residual mean meridional circulation, the zonal-mean temperature and the eddy diffusion coefficients are calculated using zonal means and three longest zonal waves of dynamical variables integrated in the semi-spectral dynamical model. They are used in the 2-D CTM to simulate the distribution of trace gases in the atmosphere. The simulated ozone is sent to the radiative transfer code to calculate the heating rates, which drive the dynamics. This radiative coupling connects the dynamical and photochemical processes and creates feedback when the atmosphere is perturbed. It is found that in most areas the ozone depletion caused by HSCT emissions calculated using the 3-wave model has the features similar to, but with significantly larger magnitude than that calculated by the AER 2-D CTM with prescribed transport parameters and temperature. The difference is mostly due to the differences in the circulation in the two models. The radiative feedback effects are investigated by comparing the ozone depletion calculated with the baseline dynamics and with the dynamics perturbed by the HSCT emissions. The feedback through changes in the residual mean meridional circulation and the eddy diffusion coefficients has moderate effects on the simulated ozone depletion. It reduced the ozone depletion by 20-30% in northern mid and high-latitudes. However, the feedback through changes in the zonal- mean temperature is negligible.
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.
Rate coefficients for the reaction C2H + C2H2 - C4H2 + H are measured over the temperature range 170-350 K. The reactions are carried out in a temperature variable flow cell. C2H radicals are produced by pulsed laser photolysis of C2H2 and a tunable infrared color-center laser is used to probe the transient removal of C2H in absorption to derive the rate coefficients. The rate coefficient is obtained and shown to be independent of temperature over the range 170-350 K. The reaction studied is of central importance for models of the photochemistry of the atmospheres of the outer planets, in particular for the satellite Titan, and the implications of the present results for these models are discussed.