Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffuse fraction”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Chain Dynamics in a Dilute Magnetorheological Fluid

The structure, formation, and dynamics of dilute, mono-dispersive ferrofluid emulsions in an external magnetic field have been investigated using dynamic light scattering techniques. In the absence of the magnetic field, the emulsion particles are randomly distributed and behave like hard spheres in Brownian motion. An applied magnetic field induces a magnetic dipole moment in each particle. Dipolar interactions between particles align them into chains where correlation functions show two decay processes. The short-time decay shows the motion of straight chains as a whole where the apparent chain length increases with the applied magnetic field and the particle volume fraction. Good scaling results are obtained showing that the apparent chain length grows with time following a power law with exponent of 0.6 and depends on the applied field, particle volume fraction, and diffusion constant of the particles. The long-time decay in the correlation function shows oscillation when the chains reach a certain length with time and stiffness with threshold field This result shows that chains not only fluctuate, but move in a periodic motion with a frequency of 364 Hz at lambda = 15. It may suggest the existence of phonons. This work is the first step in the understanding of the structure formation, especially chain coarsening mechanism, of magnetorheological (MR) fluids at higher volume fractions.

Liu, Jing↗
NASA NTRS2001

Estimating Net Primary Productivity Using Satellite and Ancillary Data

The net primary productivity (C) or annual rate of carbon accumulation per unit ground area by terrestrial plant communities is the difference of the rate of gross photosynthesis (A(sub g)) and autotrophic respiration (R) per unit ground area. Although available observations show that R is a large and variable fraction of A(sub g), viz., 0.3 to 0.7, it is generally recognized that much uncertainties exist in this fraction due to difficulties associated with the needed measurements. Additional uncertainties arise when these measurements are extrapolated to regional or global land surface using empirical equations, for example, using regression equations relating C to mean annual precipitation and air temperature. Here, a process-based approach has been taken to calculate A(sub g) and R using satellite and ancillary data. A(sub g) has been expressed as a product of radiation use efficiency, magnitude of intercepted photosynthetically active radiation (PAR), and normalized by stresses due to soil water shortage and air temperature away from the optimum range. A biophysical model has been used to determine the radiation use efficiency from the maximum rate of carbon assimilation by a leaf, foliage temperature, and the fraction of diffuse PAR incident on a canopy. All meteorological data (PAR, air temperature, precipitation, etc.) needed for the calculation are derived from satellite observations, while a land use, land cover data (based on satellite and ground measurements) have been used to assess the maximum rate of carbon assimilation by a leaf of varied cover type based on field measurements. R has been calculated as the sum of maintenance and growth components. The maintenance respiration of foliage and live fine roots at a standard temperature of different land cover has been determined from their nitrogen content using field and satellite measurements, while that of living fraction of woody stem (viz., sapwood) from the seasonal maximum leaf area index as determined from satellite observations. These maintenance respiration values were then adjusted to that corresponding to air temperature according to a prescribed non-linear variation of respiration with temperature. The growth respiration has been calculated from the difference of Ag and maintenance respiration, according to the two-compartment model. The results of calculations will be reported for 36 consecutive months (1987-1989) over large contiguous areas (ca. 10(exp 5) sq km) Of agricultural land and tropical humid evergreen forests, and compared with available field data.

Choudhury, B. J.↗
DOE OSTI2021

Environment-sensitivity functions for gross primary productivity in light use efficiency models

The sensitivity of photosynthesis to environmental changes is essential for understanding carbon cycle responses to global climate change and for the development of modeling approaches that explains its spatial and temporal variability. We collected a large variety of published sensitivity functions of gross primary productivity (GPP) to different forcing variables to assess the response of GPP to environmental factors. These include the responses of GPP to temperature; vapor pressure deficit, some of which include the response to atmospheric CO 2 concentrations; soil water availability (W); light intensity; and cloudiness. These functions were combined in a full factorial light use efficiency (LUE) model structure, leading to a collection of 5600 distinct LUE models. Each model was optimized against daily GPP and evapotranspiration fluxes from 196 FLUXNET sites and ranked across sites based on a bootstrap approach. The GPP sensitivity to each environmental factor, including CO 2 fertilization, was shown to be significant, and that none of the previously published model structures performed as well as the best model selected. From daily and weekly to monthly scales, the best model's median Nash-Sutcliffe model efficiency across sites was 0.73, 0.79 and 0.82, respectively, but poorer at annual scales (0.23), emphasizing the common limitation of current models in describing the interannual variability of GPP. Although the best global model did not match the local best model at each site, the selection was robust across ecosystem types. The contribution of light saturation and cloudiness to GPP was observed across all biomes (from 23% to 43%). Temperature and W dominates GPP and LUE but responses of GPP to temperature and W are lagged in cold and arid ecosystems, respectively. The findings of this study provide a foundation towards more robust LUE-based estimates of global GPP and may provide a benchmark for other empirical GPP products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding and control of Zener pinning via phase field and ensemble learning

Zener pinning refers to the dispersion of fine particles which influences grain size distribution via movement of grain boundaries in a polycrystalline material. Grain size distribution in polycrystals has a significant impact on their properties including physical, chemical, mechanical, and optical to name a few. We explore the use of Phase-field modeling and machine-learning techniques to understand and improve the control of grain size distribution via Zener pinning in polycrystalline materials. We develop a machine learning model that determines the relative importance of various parameters to exercise microstructure control via Zener pinning. Our workflow combines high-throughput phase-field simulations and machine learning to address the computational bottlenecks associated with large-scale simulations as well as identify features necessary for microstructure control in polycrystals. A random forest (RF) regression model was developed to predict grain sizes based on five Phase-field model parameters, achieving an average prediction error of 0.72 nm for the training data and 1.44 nm for the test data. The importance of the input parameters is analyzed using the SHapley Additive exPlanations (SHAP) approach which reveals that diffusivity, volume fraction, and particle diameter are the most important parameters in determining the final grain size. These findings will allow us to select the best second-phase particles, optimize grain size distributions and thus design microstructures with the desired properties. The developed method is a highly versatile and generalizable approach that can be used to assess the combined effects of individual features in the presence of multiple variables.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of microphase separation on the limiting current density in hybrid organic-inorganic copolymer electrolytes

Hybrid organic-inorganic block copolymer electrolytes are of interest to enable batteries containing lithium metal anodes. The conductive block is a standard polymer electrolyte of poly(ethylene oxide) and the mechanically rigid block is an inorganic poly(acryloisobutyl polyhedral oligomeric silsesquioxane) polymer. Here, in this paper, we compare a poly(acryloisobutyl polyhedral oligomeric silsesquioxane)-b-poly(ethylene oxide)-b-poly(acryloisobutyl polyhedral oligomeric silsesquioxane) (POSS-PEO-POSS) triblock copolymer and a poly(ethylene oxide)-b- poly(acryloisobutyl polyhedral oligomeric silsesquioxane) (PEO-POSS) diblock copolymer mixed with lithium bis(trifluoromethanesulfonyl)imide salt. We have experimentally measured the limiting current density in lithium symmetric cells containing hybrid organic-inorganic electrolytes at 90 °C. The cells were polarized at a large range of applied current density. The diblock copolymer electrolyte exhibited a clear plateau in cell potential at all current densities below the limiting current density. At low applied current density, the triblock copolymer electrolyte also exhibited a clear plateau in cell potential. At currents approaching the limiting current density, the triblock copolymer electrolyte exhibited an underdamped potential profile. The cell potential did not reach a plateau at current densities above the limiting current in both systems. The diblock and triblock copolymer electrolytes were fully characterized using electrochemical methods to determine the ionic conductivity, cation current fraction, salt diffusion coefficient, and open circuit voltage as a function of salt concentration. Cell potential and salt concentration as functions of position in the cell at various current densities were calculated using Newman's concentrated solution theory. The theoretical limiting current density was calculated to be the current density at which salt is depleted at the cathode. We see quantitative agreement between experimental measurements and theoretical predictions for the limiting current density in the diblock copolymer electrolyte which has an ordered structure at all salt concentrations, while the experimental limiting current density is lower than the theoretical prediction for the triblock copolymer electrolyte, which exhibits a disordered morphology at high salt concentrations.

25 ENERGY STORAGE↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2024

Radiation induced non-linear oscillations in ITER baseline scenario plasmas in DIII-D

Abstract This work shows how the radiation brought about by metals or metal-equivalent radiators such as Kr and Xe produces non-linear dynamics on otherwise stationary β N flattops of DIII-D ITER Baseline Scenario demonstration discharges. The Kr and Xe gases are used to reproduce the radiative loss rates of W in present machines that operate at core temperatures much lower than the expected ITER temperature. Experiments on DIII-D with injection of Kr and Xe, as well as with sources of intrinsic metals reach the range of radiated fraction values expected in the ITER core and experience slow oscillations in temperature and radiated power. In many cases of high radiated fraction, the core temperature decreases enough for the safety factor profile to rise above the 1/1 rational surface, naturally eliminating sawteeth and occasionally producing a persistent helical core. The oscillations can be reproduced by a modified Lotka–Volterra system for temperature and radiated fraction if diffusion and noise are included, which indicates that the interplay between temperature and radiation can be the main cause of the cyclic nature of the system. A new physics based model which includes equations for temperature, density and input power can also reproduce the oscillations observed in the experiments. The present results suggest that the non-linearity of the system can be increased by the inclusion of the inherently non-linear alpha heating term, which is proportional to ∼ n e 2 T i 2 , and obtains oscillations in the model when added to an otherwise more stationary system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Emergence of liquid following laser melting of gold thin films

X-ray structural science is undergoing a revolution driven by the emergence of X-ray Free-electron Laser (XFEL) facilities. The structures of crystalline solids can now be studied on the picosecond time scale relevant to phonons, atomic vibrations which travel at acoustic velocities. In the work presented here, X-ray diffuse scattering is employed to characterize the time dependence of the liquid phase emerging from femtosecond laser-induced melting of polycrystalline gold thin films using an XFEL. In a previous analysis of Bragg peak profiles, we showed the supersonic disappearance of the solid phase and presented a model of pumped hot electrons carrying energy from the gold surface to scatter at internal grain boundaries. This generates melt fronts propagating relatively slowly into the crystal grains. By conversion of diffuse scattering to a partial X-ray pair distribution function, we demonstrate that it has the characteristic shape obtained by Fourier transformation of the measured F ( Q ). The diffuse signal fraction increases with a characteristic rise-time of 13 ps, roughly independent of the incident pump fluence and consequent final liquid fraction. This suggests the role of further melt-front nucleation processes beyond grain boundaries.

36 MATERIALS SCIENCE↗
NASA NTRS1982

A numerical study of the effects of ambipolar diffusion on the collapse of magnetic gas clouds

The gravitational collapse of isothermal, nonrotating magnetic gas clouds have been calculated numerically, including the effects of ambipolar diffusion. The fractional ionization in the clouds is approximated by a power-law function of the gas density, f = K/n to the q-power, where K and q are adjustable parameters. Eleven numerical experiments were run, and the results indicate that the asymptotic character of collapse is determined mainly by the value of q and is largely independent of the other parameters characterizing a cloud (e.g., K, cloud mass). In particular, there is nearly a one-to-one correspondence between q and the slope, x, of the central magnetic field strength-gas density relationship. If q is no more than 0.8, a cloud collapses asymptotically, as though the magnetic field were 'frozen' to the neutral matter. The magnetic field strength at the center of a collapsing cloud is strongly amplified during collapse even for values of q of about 1, despite extremely low values of fractional ionization. A discussion of the theoretical basis for this unexpected behavior is given. Possible implications of our results for the problems of magnetic braking of rotating protostars and star formation in general are also presented.

Black, D. C.↗
NASA NTRS1992

Prediction and measurement of soil bidirectional reflectance

The present model for soil bidirectional reflectance distribution functions in the visible and reflective IR ranges is based on a representation of soil surfaces by a collection of opaque spheres atop a Lambertian horizontal surface. The model is not sensitive to increases in the sphere area index beyond a value of 0.4. Observed soil reflectance factor data for a tilled field are noted to be consistent with those predicted by the model's expression of soil reflectance as a function of the horizontal area obscured by a sphere, the amount of surface in shadow, and relative brightness of a sphere, and the fraction of diffuse illumination.

Irons, James R.↗
NASA NTRS1999

Growing wheat in Biosphere 2 under elevated CO2: observations and modeling

Spring wheat (Triticum aestivum L., cv. Yecora Rojo) was grown in the intensive agricultural biome (IAB) of Biosphere 2 during the l995-l996 winter/spring season. Environmental conditions were characterized by a day/night temperature regime of 27/17 degrees C, relative humidity (RH) levels around 45%, mean atmospheric CO2 concentration of 450 ppmv, and natural light conditions with mean intensities about half of outside levels. Weekly samples of above-ground plant matter were collected throughout the growing season and phenological events recorded. A computer model, CERES-Wheat, previously tested under both field and controlled conditions, was used to simulate the observed crop growth and to help in data analysis. We found that CERES-Wheat simulated the data collected at Biosphere 2 to within 10% of observed, thus suggesting that wheat growth inside the IAB was comparable to that documented in other environments. The model predicts phenological stages and final dry matter (DM) production within l0% of the observed data. Measured DM production rates, normalized for light absorbed by the crop. suggested photosynthetic efficiencies intermediate between those observed under optimal field conditions and those recorded in NASA-Controlled Ecological Life-Support Systems (CELSS). We suggest that such a difference can be explained primarily in terms of low light levels inside the IAB, with additional effects due to elevated CO2 concentrations and diffuse light fractions.

Non-NASA Center↗
NASA NTRS1944

The Surface Structure of Ground Metal Crystals

The changes produced on metallic surfaces as a result of grinding and polishing are not as yet fully understood. Undoubtedly there is some more or less marked change in the crystal structure, at least, in the top layer. Hereby a diffusion of separated crystal particles may be involved, or, on plastic material, the formation of a layer in greatly deformed state, with possible recrystallization in certain conditions. Czochralski verified the existence of such a layer on tin micro-sections by successive observations of the texture after repeated etching; while Thomassen established, roentgenographically by means of the Debye-Scherrer method, the existence of diffused crystal fractions on the surface of ground and polished tin bars, which he had already observed after turning (on the lathe). (Thickness of this layer - 0.07 mm). Whether this layer borders direct on the undamaged base material or whether deformed intermediate layers form the transition, nothing is known. One observation ty Sachs and Shoji simply states that after the turning of an alpha-brass crystal the disturbance starting from the surface, penetrates fairly deep (approx. 1 mm) into the crystal (proof by recrystallization at 750 C).

Boas, W.↗
NASA NTRS2014

Constraints on Transport and Emplacement Mechanisms of Labile Fractions in Lunar Cold Traps

Sustaining the scientific exploration of the Solar System will require a significant proportion of the necessary fuels and propellants, as well as other bulk commodities, to be produced from local raw materials [1]. The viability of mineral production depends on the ability to locate and characterize mineable deposits of the necessary feedstocks. This requires, among other things, a workable understanding of the mechanisms by which such deposits form, which is the subject of Economic Geology. Multiple deposition scenarios are possible for labile materials on the Moon. This paper suggests labile fractions moved diffusely through space; deposits may grow richer with depth until low porosity rock; lateral transport is likely to have occurred with the regolith, at least for short distances; crystalline ice may not exist; the constituent phases could be extremely complex. At present we can constrain the sources only mildly; once on the Moon, the transport mechanisms inherently mix and therefore obscure the origins. However, the importance of expanding our understanding of ore-forming processes on the Moon behooves us to make the attempt. Thus begins a time of new inquiry for Economic Geology.

Rickman, D.↗
DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Heat-Pulse Measurements Reveal Fiber Volume Fractions

Measurements of thermal diffusivities by heat-pulse method constitutes basis of noncontact, nondestructive method of determining fiber volume fractions (FVFs) of samples of composite (matrix/fiber) materials. Modulated radiant heat applied to sample, while resulting modulated infrared emitted is monitored. Phase shift between two modulations indicative of thermal diffusivity and fiber volume fraction of sample. Testing method takes less time, provides data to characterize sample through its thickness, and amenable to scanning for global determination of gradual spatial variations in FVF along sample.

Winfree, William P.↗
DOE OSTI2023

Topology driven and soft phonon mode enabled Na-ion diffusion in quaternary chalcogenides, Na 3 ZnGaX 4 (X = S, and Se)

The compounds Na 3 ZnGaX 4 (X = S, Se) are potential solid electrolyte materials in sodium-based batteries, which have certain advantages over oxide materials and have shown significant ionic conductivity at ambient temperature. In this paper, we bring out atomic-level features of the diffusion process in these new materials using the microscopic techniques of inelastic neutron scattering (INS), quasielastic neutron scattering (QENS), and ab initio molecular dynamics (AIMD) simulations. The insights obtained from these techniques are unique and not available from other macroscopic experiments. Neutron scattering experiments have been performed at temperatures from 100 to 700 K. The simulations have been carried out up to 900 K. We have calculated the phonon spectra and the space–time correlation functions and found good agreement with the results of the neutron scattering experiments. The simulations enable detailed analysis of the atomic-site dependent dynamical information. We observe low-energy phonon modes of ~6 meV involving the vibrations of certain Na atoms in the lattice. This reveals that the Na at 32g Wyckoff sites (Na2) has sufficiently shallow potential among the two available crystallographic sites. This shallow potential facilitates diffusion. Furthermore, the specific structural topology of the network of interconnected zig-zag chains of the Na2 atomic sites provides the low-barrier energy pathways for diffusion. A small fraction of vacancy defects appears essential for diffusion. We further observe that the Na2 atoms undergo jump-like diffusion to the vacant next or the 2nd next neighbour sites at ~4 Å. While the QENS experiments reveal the jump-like diffusion and its time scale, detailed analysis of the AIMD simulations shows that the jumps appear mostly along zig-zag chains of the Na2 sites in the tetragonal ab-plane, as well as between the chains along the c-axis. In conclusion, the net diffusion is essentially 3-dimensional, with little anisotropy despite the anisotropy of the tetragonal crystal structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Equilibrium Fe isotope fractionation between olivine, pyroxene, spinel and MORB glass: Implications for mantle partial melting to generate MORBs

Primitive mid-ocean ridge basalts (MORBs) exhibit Fe isotopic compositions heavier than the upper mantle by +0.074 ± 0.028 ‰ for δ 56 Fe. The processes responsible for this isotopic difference remain unclear. Modeling of Fe isotope fractionation during mantle partial melting requires reliable equilibrium Fe isotope fractionation factors between minerals and melts, for which consistent data are still lacking. Here, in this study, we used Nuclear Resonant Inelastic X-ray Scattering (NRIXS) technique to measure Fe force constants for a MORB glass (ALV 519-4-1) and natural mantle minerals (olivine, orthopyroxene, clinopyroxene, and spinel) to determine the equilibrium Fe isotope fractionation factors between them. The force constants determined in this study, in increasing order, are 167 ± 26 N/m for spinel, 175 ± 17 N/m for olivine, 176 ± 20 N/m for MORB glass, 205 ± 26 N/m for clinopyroxene, and 219 ± 36 N/m for orthopyroxene. We evaluated the previously proposed mechanisms for the heavy Fe isotopic composition of MORBs, including (i) mantle partial melting, (ii) mantle lithological heterogeneity, with pyroxenite in the source, (iii) mantle metasomatism by low-degree melts, and (iv) fractional crystallization of olivine from melts. For (i), we used the pMELTS program to simulate adiabatic decompression melting of mantle peridotites, and calculated Fe isotope fractionation based on Fe 3+ –Fe 2+ equilibrium-controlled fractionation, where Fe 3+ forms stronger bonds and is more incompatible than Fe 2+ . At 10 wt% peridotite melting, corresponding to MORB generation, only +0.03 ‰ Fe isotope fractionation between the melt and the original bulk composition (Δ 56 Fe = δ 56 Fe melt - δ 56 Fe 0 ) was produced, insufficient to account for the observed MORB-upper mantle difference. For (ii), melting of pyroxenites yields smaller Fe isotope fractionation than melting of peridotites, making it unlikely the cause for the MORB-upper mantle isotopic difference. For (iii), both the Fe 3+ /ΣFe ratio and the δ 56 Fe of melts increase with the degree of partial melting, indicating that low-degree melts are not isotopically heavy enough to significantly alter the isotopic composition of lithospheric mantle through metasomatism. For (iv), equilibrium isotope fractionation between olivine and melt is near zero. These results suggest that equilibrium Fe isotope fractionation alone cannot explain the MORB isotopic signature, highlighting the potential role of kinetic isotope fractionation. Using a diffusion model, we calculated kinetic Fe and Mg isotope fractionations associated with (iv) olivine crystallization from a melt, and found that the predicted Fe and Mg isotope fractionations were inconsistent with observations in MORBs. Qualitatively, two processes could have induced kinetic Fe isotope fractionation during MORB generation: (a) Fe-Mg interdiffusion between melt and solid during melt migration and (b) reactive melt-rock interactions during melt focusing. However, a quantitative understanding of their role in modifying the melt isotopic composition remains limited and requires further investigation.

Fe isotopes↗
NASA NTRS1994

The Prediction of Nozzle Performance and Heat Transfer in Hydrogen/Oxygen Rocket Engines with Transpiration Cooling, Film Cooling, and High Area Ratios

An advanced engineering computational model has been developed to aid in the analysis of chemical rocket engines. The complete multispecies, chemically reacting and diffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and Dufour energy transfer terms. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film-cooled nozzle, and a transpiration-cooled plug-and-spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 nozzle and the film-cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. It is demonstrated that thermal diffusion has a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle and was shown to represent a significant fraction of the diffusion fluxes occurring in the transpiration-cooled rocket engine.

Kacynski, Kenneth J.↗