Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion models”

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 1,423 records · Page 79

Oxidative Attack of Carbon/Carbon Substrates through Coating Pinholes

A critical issue with oxidation protected carbon/carbon composites used for spacecraft thermal protection is the formation of coating pinholes. In laboratory experiments, artificial pinholes were drilled through SiC-coatings on a carbon/carbon material and the material was oxidized at 600, 1000, and 1400 C at reduced pressures of air. The attack of the carbon/carbon was quantified by both weight loss and a novel cross-sectioning technique. A two-zone, one dimensional diffusion control model was adapted to analyze this problem. Agreement of the model with experiment was reasonable at 1000 and 1400 C; however results at lower temperatures show clear deviations from the theory suggesting that surface reaction control plays a role.

Jacobson, Nathan S.↗

Experiment 1: Orbital Processing of High-Quality Zn-Alloyed CdTe Compound Semiconductors

Zinc-alloyed cadmium telluride (Zn:CdTe) crystals were grown in unit gravity and in microgravity for comparative analysis. Two crystals were grown on USML-1/STS-50 in 1992, and two additional crystals were grown on USML-2/STS-73 in 1995. The Crystal Growth Furnace (CGF) in the seeded Bridgman-Stockbarger crystal growth geometry was utilized on both missions. Crystals grown on USML-1/STS-50 were found to have solidified with partial wall contact due to the near-absence of the hydrostatic pressure in microgravity, a residual g-vector that was not axial, and the nonwetting sample/ampoule wetting conditions. Crystals grown on USML-2/STS-73 included: a sample/ampoule identical to the USML-1/STS-73 sample/ampoule, with the addition of a restraining spring to simulate hydrostatic pressure internally, and a tapered ampoule which accomplished 2.2 cm of crystal growth without wall contact. Infrared transmission of all ground and flight samples was found to be 63-66 percent, very close to the theoretical 66 percent, suggesting good stoichiometric control. Infrared microscopy confirmed that the primary precipitates were Te and their size (1-10 microns) and density suggested that the flight- and ground-base samples experienced similar thermal histories. Longitudinal macrosegregation, calculated using scaling analysis, was predicted to be low. Nearly diffusion controlled growth was achieved even in unit gravity and macrosegregation data could be fit with a diffusion controlled model. Radial segregation was monitored and was found to vary with fraction solidified, particularly through the shoulder region, where the sample cross section was varying significantly. It was also disturbed in the flight samples in regions where asymmetric wall contact was noted. In regions where a steady state was established, the radial segregation was invariant within our experimental measurement error. Flight samples were found to be much higher in structural perfection than samples processed in unit gravity under identical growth conditions. In regions where solidification had occurred without wall contact, the free surfaces evidenced virtually no twinning, although twins appeared in the flight samples in regions of wall contact and were pervasive in the ground samples. These results were confirmed using optical microscopy and synchrotron x-ray white beam topography. Full-width half-maximum rocking curve widths (FWHM), recorded in arc-seconds, were significantly reduced from 20 a-s (1 g) to 9 a-s (micro-g) for the best regions of the crystals. The 9 a-s (FWHM) rocking curve value in the unconfined flight samples equals the best value reported terrestrially for this material. The ground samples exhibited a fully developed (111)[110] dislocation mosaic structure, whereas dislocations within the flight samples were discrete and no mosaic structure was evident. The defect density was quantitatively reduced from 75,000 (1 g) to 800 (micro-g) +/-50 percent. Dislocation etch pit density results were confirmed using transmission synchrotron white beam and monochromated beam topography. The low defect density is thought to have resulted from the near absence of hydrostatic pressure, which allowed the molten boule to solidify with little or no wall contact. This minimized the transfer of hoop stresses during solidification and post-solidification processing.

Larson, D. J., Jr.↗

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

SiC and Si3N4 materials were tested under various turbine engine combustion environments, chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high speed aircraft. Representative CVD, sintered, and composite materials were evaluated in both furnace and high pressure burner rig exposure. While protective SiO2 scales form in all cases, evidence is presented to support paralinear growth kinetics, i.e. parabolic growth moderated simultaneously by linear volatilization. The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were used to map SiO2 volatility (and SiC recession) over a range of temperature, pressure, and velocity. The functional dependency of material recession (volatility) that emerged followed the form: exp(-QIRT) * P(exp x) * v(exp y). These empirical relations were compared to rates predicted from the thermodynamics of volatile SiO and SiO(sub x)H(sub Y) reaction products and a kinetic model of diffusion through a moving, boundary layer. For typical combustion conditions, recession of 0.2 to 2 micron/h is predicted at 1200- 1400C, far in excess of acceptable long term limits.

Smialek, James L.↗

Vapor Growth and Characterization of ZnSeTe Solid Solutions

Six ZnSe(1-x)Te(x) crystals were grown by the physical vapor transport technique. For each of the source material compositions, x = 0.10, 0.20 and 0.30, two crystals were grown - one under the horizontal and the other under the vertical stabilized configurations. The axial and radial compositional uniformity were measured by precision density measurements, wavelength dispersive X-ray spectroscopy (WDS) and optical transmission mappings. The measured radial ZnTe content was quite uniform for all the grown crystals except the horizontally grown crystal for x = 0.30. The WDS results on this crystal indicated a core with uniform ZnTe content, about 0.38, surrounded by a thin region of high ZnTe content with x = 0.8. This feature was confirmed by the SEM back scattering electron images. For the three source compositions the axial compositional variations for the vertically grown crystals were more uniform than that for the horizontally grown crystals. The measured compositions in the crystals grown from source composition of x = 0.10 suggest that the transport mechanism in the system can not be interpreted by a simple one-dimensional diffusion limited model.

Su, Ching-Hua↗

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

Silicon carbide (SiC) and Si3N4 materials were tested in various turbine engine combustion environments chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high-speed aircraft. Representative chemical vapor-deposited (CVD), sintered, and composite materials were evaluated by furnace and high-pressure burner rig exposures. Although protective SiO2 scales formed in all cases, the evidence presented supports a model based on paralinear growth kinetics (i.e., parabolic growth moderated simultaneously by linear volatilization). The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were thus used to map SiO2 volatility (and SiC recession) over a range of temperatures, pressures, and velocities. The functional dependency of material recession (volatility) that emerged followed the form A[exp(-Q / RT)](P(sup x)v(sup y). These empirical relations were compared with rates predicted from the thermodynamics of volatile SiO and SiOxHy reaction products and a kinetic model of diffusion through a moving boundary layer. For typical combustion conditions, recession of 0.2 to 2 micrometers/hr is predicted at 1200 to 1400 C, far in excess of acceptable long-term limits.

Smialek, James L.↗

Flow field topology of transient mixing driven by buoyancy

Transient mixing driven by buoyancy occurs through the birth of a symmetric Rayleigh-Taylor morphology (RTM) structure for large length scales. Beyond its critical bifurcation the RTM structure exhibits self-similarity and occurs on smaller and smaller length scales. The dynamics of the RTM structure, its nonlinear growth and internal collision, show that its genesis occurs from an explosive bifurcation which leads to the overlap of resonance regions in phase space. This event shows the coexistence of regular and chaotic regions in phase space which is corroborated with the existence of horseshoe maps. A measure of local chaos given by the topological entropy indicates that as the system evolves there is growth of uncertainty. Breakdown of the dissipative RTM structure occurs during the transition from explosive to catastrophic bifurcation; this event gives rise to annihilation of the separatrices which drives overlap of resonance regions. The global bifurcation of explosive and catastrophic events in phase space for the large length scale of the RTM structure serves as a template for which mixing occurs on smaller and smaller length scales. Copyright 2004 American Institute of Physics.

Nonlinear Dynamics↗

Oxidation Through Coating Cracks of SiC-Protected Carbon/Carbon

The oxidation of SiC-protected carbon/carbon through machined slots and naturally occurring craze cracks in the SiC was studied. The slot and crack geometries were characterized, and the subsurface oxidation of the carbon/carbon substrate at temperatures of 1000 to 1300 C in air was assessed using weight change, x-ray computed tomography, and optical microscopy of sections. Rate constants were derived from these measurements and compared with a two-step diffusion control model of carbon oxidation. Oxidation kinetic measurements on both the specimens with machined slots and with naturally occurring craze cracks showed good agreement with the model.

Jacobson, Nathan S.↗

Oxidation Through Coating Cracks of SiC-Protected Carbon/Carbon

The oxidation of SiC-protected carbon/carbon through machined slots and naturally occurring craze cracks in the SiC was studied. The slot and crack geometries were characterized, and the subsurface oxidation of the carbon/carbon substrate at temperatures of 1000 to 1300 C in air was assessed using weight change, x-ray computed tomography, and optical microscopy of sections. Rate constants were derived from these measurements and compared with a two-step diffusion control model of carbon oxidation. Oxidation kinetic measurements on both the specimens with machined slots and with naturally occurring craze cracks showed good agreement with the model.

Jacobson, Nathan S.↗

Oxidation Through Coating Cracks of SiC-Protected Carbon/Carbon

The oxidation of SiC-protected carbon/carbon through machined slots and naturally occurring craze cracks in the SiC was studied. The slot and crack geometries were characterized, and the subsurface oxidation of the carbon/carbon substrate at temperatures of 1000 to 1300 C in air was assessed using weight change, x-ray computed tomography, and optical microscopy of sections. Rate constants were derived from these measurements and compared with a two-step diffusion control model of carbon oxidation. Oxidation kinetic measurements on both the specimens with machined slots and with naturally occurring craze cracks showed good agreement with the model.

Jacobson, Nathan S.↗

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations (LES) were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver was used to compute the flow from a 2-inch round nozzle. Three different flow conditions of varying jet Mach numbers and temperature ratios were examined. The LES results showed that the temperature field behaves similar to the velocity field, but with a more rapidly spreading mixing layer. Predictions of mean, mu-bar(sub i), and fluctuating, mu'(sub i), velocities were compared to particle image velocimetry data. Predictions of mean, T-bar, and fluctuating, T', temperature were compared to data obtained using Rayleigh scattering and Raman spectroscopy. Very good agreement with experimental data was demonstrated for the mean and fluctuating velocities. The LES correctly predicts the behavior of the turbulent temperature field, but over-predicts the levels of the fluctuations. The turbulent heat flux was examined and compared to Reynolds-averaged Navier-Stokes (RANS) results. The LES and RANS simulations produced very similar results for the radial heat flux. However, the axial heat flux obtained from the LES differed significantly from the RANS result in both structure and magnitude, indicating that the gradient diffusion type model in RANS is inadequate. Finally, the LES data was used to compute the turbulent Prandtl number and verify that a constant value of 0.7 used in the RANS models is a reasonable assumption.

hot jets↗

Solar Wind Implantation into Lunar Regolith II: Monte Carlo Simulations of Hydrogen Retention in a Surface with Defects and the Hydrogen (H, H2) Exosphere

Recently, the near-infrared observations of the OH veneer on the lunar surface by the Moon Mineralogy Mapper (M3) have been refined to constrain the OH content to 500-750 parts per million (ppm). The observations indicate diurnal variations in OH up to 200 ppm possibly linked to warmer surface temperatures at low latitude. We examine the M3 observations using a statistical mechanics approach to model the diffusion of implanted H in the lunar regolith. We present results from Monte Carlo simulations of the diffusion of implanted solar wind H atoms and the subsequently derived H and H2 exospheres.

Exosphere↗

Analysis of Fibrous Felts for Flexible Ablators Using Synchrotron Hard X-Ray Micro-Tomography

We analyzed the material properties of low-density felts that are used as substrates for new-generation flexible and conformal carbon/phenolic ablators, and compared them with those of a rigid carbon fiber preform that is used to manufacture rigid carbon/phenolic ablators. Micro-tomography measurements were obtained using synchrotron X-rays, allowing the characterization of the materials microstructure at the scale of the fibers. Using the tomography voxels as computational grids, we computed tortuosity and room temperature conductivity. In addition we performed micro-scale simulations of the oxidation of carbon fibers using a random walk model for oxygen diffusion and a probability law to model surface reactions.

Ablation↗

Modeling Validation and Its Design Relevance for NTP Fuel Elements

The development of Nuclear Thermal Propulsion engines depends heavily on modeling and simulation of a number of disciplines—neutronic / fluid / thermal / structural are only a start. Further, rockets involve extreme conditions, particularly high material temperatures in the nuclear fuel. This paper considers validation of the underling physical models and codes and how their accuracy influences several design challenges. Several issues are identified, where models, material data, could enable design solutions. In particular, thermal and CTE mismatch stresses are a concern, and stress predictions would benefit from improved mechanical property data at high temperatures. Nuclear inter-element heating causes propellant flow maldistribution and performance risks, but a simple model provides insights. High temperature chemical diffusion should be modeled to understand fuel element mass loss.

Nuclear Propulsion↗

Modeling Validation and Its Design Relevance for NTP Fuel Elements

The development of Nuclear Thermal Propulsion engines depends heavily on modeling and simulation of a number of disciplines—neutronic / fluid / thermal / structural are only a start. Further, rockets involve extreme conditions, particularly high material temperatures in the nuclear fuel. This paper considers validation of the underling physical models and codes and how their accuracy influences several design challenges. Several issues are identified, where models, material data, could enable design solutions. In particular, thermal and CTE mismatch stresses are a concern, and stress predictions would benefit from improved mechanical property data at high temperatures. Nuclear inter-element heating causes propellant flow maldistribution and performance risks, but a simple model provides insights. High temperature chemical diffusion should be modeled to understand fuel element mass loss.

Nuclear Propulsion↗

Apollo Sample 64455: Petrologic and Geochemical Characterization of a Glass-Coated Impact Melt Rock

Introduction and Background: Impact melts are a key product of the impact cratering process. Impact melts can turbulently entrain pre-existing rock fragments, and melt part of those entrained materials until thermal equilibrium is achieved [1–3]. Lunar impact melt rocks are valuable for determining the precise ages of basin-forming and cratering events, which anchor the cratering chronology and dynamics of the entire solar system [4, and references therein]. Additionally, entrained clasts can inform on the petrology and age of crustal materials that predate the melt-forming impact. In lithologies where no pristine samples are available, entrained materials may hold the record of endogenous lunar volatiles, provided we can untangle their impact histories as high temperature impact melt can heat the clasts and may alter their native volatiles or isotopic signatures [5–6]. Here, we characterize an Apollo 16 sample, 64455,70B, in an effort to understand the physical, microstructural, and geochemical consequences of impact melt contact metamorphism. We acquired optical light and x-ray maps of the sample, as well as quantitative analyses via electron probe microanalysis (EPMA). These tasks allow us to better understand the interaction between impact melts and entrained clasts, ultimately closing knowledge gaps in high-temperature impact processing on planetary surfaces. Sample Description: Sample 64455 is an oblong crystalline rock encased in glass to form an egg-like shape [7–8]. The impact melt coating has altered the outer (1 to 2 mm) rind of the crystalline rock and has preserved schlieren. The interior crystalline rock is likely a crystalline impact melt rock based on its bulk Ir content of 2.25 parts per billion [9]. The composition of the impact melt glass rind cannot be solely derived from the crystalline interior, the glassy rind is more aluminous and less magnesian [10]. Neither the crystalline interior nor the glassy rind were radiometrically age dated, although several studies investigated its cosmogenic isotopes and exposure ages, which reveal the sample was exposed on the lunar surface for ~1–2 Ma [11–13]. Preliminary Results: The interior crystalline rock consists primarily of plagioclase with interstitial pyroxene and olivine. Plagioclase grains are zoned and compositions range from An97–89, while pyroxene compositions generally range from En68Fs19Wo13 to En82Fs15Wo3 and olivine are Fo77 to Fo85. Minor phases include Fe-Ni metal, Ca-phosphates, a K- and Si-rich phase, sulfides, and Fe-phosphides. The glass rind contains approximately 44 wt.% SiO2, 25 wt.% Al2O3, 6 wt.% FeO, 8 wt.% MgO, and 14 wt.% CaO, with other oxides at <1 wt.% each, although the composition of the glass becomes more variable closest to the crystalline rock. The rind consists of rounded plagioclase, pyroxene, and olivine crystals with interstitial melt. In some locations, the rounded crystals are surrounded by closely-packed, blade-like crystals, <5 µm in width. These results indicate that the crystalline interior has partially melted where in contact with the high-temperature impact melt and while further investigation is needed, pyroxene may have preferentially melted relative to olivine and plagioclase. Future Work: We plan to conduct electron backscatter diffraction (EBSD) mapping in order to determine crystallinity and crystal orientation. Next using transmission electron microscopy (TEM) data in conjunction with EBSD and EPMA information, we will model the diffusion of moderately volatile elements between the melt and remnant crystals.

A C Stadermann↗

Does the way we do science foster discovery?

Freedom to explore the unknown is key to scientific discovery. Maximizing modern individualistic measures of scientific productivity like citations and number of publications may impede the progress of science as a whole.

discovery↗

The Hourly GHI, DHI and DNI: Intercomparison of the CERES-Based Data and the NSRDB Data through Comparison with the Ground-Based BSRN Data

The NASA Clouds and the Earth’s Radiant Energy System (CERES) SYN1deg(Ed4.1) provides hourly global horizontal irradiance (GHI), diffuse horizontal irradiance (DHI) and direct horizontal irradiance (DirHI) at 1-degree latitude by 1-degree longitude resolution, and the time span is from March 2000 to near present. While the GHI agrees well with the ground-based Baseline Surface Radiation Network (BSRN) data, the DHI and DirHI are, respectively, positively and negatively biased against the BSRN data, and so is the direct normal irradiance (DNI) derived simply by dividing the DirHI by the cosine of the hourly solar zenith angle (SZA). We found that the biases tend to be well-defined functions of cos(SZA) and the cloud fraction (CLFR). We thus performed a bias-correction on the hourly DHI and DNI in the latitude-cos(SZA)-CLFR phase space, and therefrom, we derived the hourly global tilted irradiance (GTI) at a number of tilt angles using the isotropic diffuse irradiance model, and the results agree well with the GTI derived from the original BSRN near-instantaneous records. Meanwhile, the National Renewable Energy Laboratory (NREL) has produced the National Solar Radiation Database (NSRDB) with finer spatiotemporal resolution using its Physical Solar Model (PSM) and surface-based measured data, model results and satellite-based data as inputs. Now the NSRDB data covers not only the United States and its neighboring regions, but Europe, Africa, Asia and Oceania. The available temporal resolutions are 5, 10, 15, 30 and 60 minutes and the spatial resolutions are 2, 4 and 10 km, depending on the selected region. In this presentation, we compare both the CERES-based and NSRDB GHI, DHI and DNI with their BSRN counterparts and show how finer spatial resolution may give us an advantage in data accuracy and usability.

Taiping Zhang↗

Methane Hydrate Formation from Enhanced Organic Carbon Burial During Glacial Lowstands: Examples from the Gulf of Mexico

Methane hydrates in fine-grained marine sediments are often found within veins and fractures occupying discrete depth intervals that are surrounded by hydrate-free sediments. As they are not connected with gas sources beneath the base of the methane hydrate stability zone (MHSZ), these isolated hydrate-bearing intervals have been interpreted as formed by in situ microbial methane. We investigate here the hypothesis that these hydrate deposits form in sediments that were deposited during glacial lowstands and contain higher amounts of labile particulate organic carbon (POC), leading to enhanced microbial methanogenesis. During Pleistocene lowstands, river loads are deposited near the steep top of the continental slope and turbidity currents transport organic-rich, fine-grained sediments to deep waters. Faster sedimentation rates during glacial periods result in better preservation of POC because of decreased exposure times to oxic conditions. The net result is that more labile POC enters the methanogenic zone and more methane is generated in these sediments. To test this hypothesis, we apply an advection-diffusion-reaction model with a time-dependent deposition of labile POC at the seafloor controlled by glacioeustatic sea level variations in the last 250 kyr. The model is run for parameters estimated at three sites drilled by the 2009 Gulf of Mexico Joint Industry Project: Walker Ridge in the Terrebonne Basin (WR313-G and WR313-H) and Green Canyon near the canyon embayment into the Sigsbee Escarpment (GC955-H). In the model, gas hydrate forms in sediments with higher labile POC content deposited during the glacial cycle between 230 and 130 kyr (marine isotope stages 6 and 7). The corresponding depth intervals in the three sites contain hydrates, as shown by high bulk electrical resistivities and resistive subvertical fracture fills. This match supports the hypothesis that enhanced POC burial during glacial lowstands can result in hydrate formation from in situ microbial methanogenesis. Our results have implications for carbon cycling during glacial/interglacial cycles and for hydrate accumulation in the MHSZ. In particular, once hydrate-bearing intervals formed during glacial periods are buried beneath the MHSZ and dissociate, gas bubbles can rise and recycle microbial methane into the MHSZ.

Malinverno, Alberto↗