Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion coefficients”

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,297 records · Page 72

Micromechanical thermal analysis of interphase region in a titanium aluminide MMC

The high reactivity between the fiber and matrix in silicon carbide/titanium aluminide MMCs leads to the formation of brittle reaction products at the fiber/matrix (F/M) interface. Also, the high thermal expansion coefficient mismatch between the fiber and matrix leads to high tensile residual stresses at the F/M interface, and this can lead to premature cracking during cooldown. One solution to these problems is the use of a metallic fiber coating like Ta which acts as an F/M diffusion barrier and reacts with the matrix to form a beta stabilized compliant layer. A finite element micromechanics analysis was performed to study the effects of Ta and beta interphase layers on the thermal residual stresses during consolidation. A 5-micron-thick beta layer reduced cool-down stresses by 8 percent compared to a 2 percent reduction computed for a Ta layer of the same thickness. Plastic yielding in the Ta was not effective in reducing cool-down stresses. Compliant alpha-2 particles next to the stiffer gamma particles reduced stresses in the gamma particles by less than 2 percent. A simple closed form analysis was developed to calculate thermal residual stresses in a fiber/interphase/matrix system. A 2-micron-thick Ag interphase layer was found to reduce residual stresses by about 11 percent.

Naik, R. A.↗

Effect of Film-Hole Shape on Turbine Blade Film Cooling Performance

The detailed heat transfer coefficient and film cooling effectiveness distributions as well as tile detailed coolant jet temperature profiles on the suction side of a gas turbine blade A,ere measured using a transient liquid crystal image method and a traversing cold wire and a traversing thermocouple probe, respectively. The blade has only one row of film holes near the gill hole portion on the suction side of the blade. The hole geometries studied include standard cylindrical holes and holes with diffuser shaped exit portion (i.e. fanshaped holes and laidback fanshaped holes). Tests were performed on a five-blade linear cascade in a low-speed wind tunnel. The mainstream Reynolds number based on cascade exit velocity was 5.3 x 10(exp 5). Upstream unsteady wakes were simulated using a spoke-wheel type wake generator. The wake Strouhal number was kept at 0 or 0.1. Coolant blowing ratio was varied from 0.4 to 1.2. Results show that both expanded holes have significantly improved thermal protection over the surface downstream of the ejection location, particularly at high blowing ratios. However, the expanded hole injections induce earlier boundary layer transition to turbulence and enhance heat transfer coefficients at the latter part of the blade suction surface. In general, the unsteady wake tends to reduce film cooling effectiveness.

Han, J. C.↗

Preliminary Investigation of a New Type of Supersonic Inlet

A supersonic inlet with supersonic deceleration of the flow entirely outside of the inlet is considered. A particular arrangement with fixed geometry having a central body with a circular annular intake is analyzed, and it is shown theoretically that this arrangement gives high pressure recovery for a large range of Mach number and mass flow and therefore is practical for use on supersonic airplanes and missiles. For some Mach numbers the drag coefficient for this type of inlet is larger than the drag coefficient for the type of inlet with supersonic compression entirely inside, but the pressure recovery is larger for all flight conditions. The differences in drag can be eliminated for the design Mach number. Experimental results confirm the results of the theoretical analysis and show that pressure recoveries of 95 percent for Mach numbers of 1.33 and 1.52, 92 percent for a Mach number of 1.72, and 86 percent for a Mach number of 2.10 are possible, with the configurations considered. If the mass flow decreases, the total drag coefficient increases gradually and the pressure recovery does not change appreciably. The results of this work were first presented in a classified document issued in 1946.

DIFFUSERS, SUPERSONIC↗

Diffuse reflectance of the ocean - The theory of its augmentation by chlorophyll a fluorescence at 685 nm

The radiative transfer equation is modified to include the effect of fluorescent substances and solved in the quasi-single scattering approximation for a homogeneous ocean containing fluorescent particles with wavelength independent quantum efficiency and a Gaussian shaped emission line. The results are applied to the in vivo fluorescence of chlorophyll a (in phytoplankton) in the ocean to determine if the observed quantum efficiencies are large enough to explain the enhancement of the ocean's diffuse reflectance near 685 nm in chlorophyll rich waters without resorting to anomalous dispersion. The computations indicate that the required efficiencies are sufficiently low to account completely for the enhanced reflectance. The validity of the theory is further demonstrated by deriving values for the upwelling irradiance attenuation coefficient at 685 nm which are in close agreement with the observations.

Gordon, H. R.↗

An objective method for determining the generalized transport tensor for two-dimensional Eulerian models

An objective method for deriving the components of a generalized transport tensor for a two-dimensional model is presented. Representative meridional and vertical velocities and thermodynamic scalars at a uniform grid are used to reduce the problem to the solution of two flux equations for two unknowns. One unknown is the stream-function, coefficient of an antisymmetric tensor, which corrects the Eulerian mean motions for Stokes drift; the other is a time constant, which converts the deviatory velocity tensor to a symmetric transport tensor. The complete asymmetric tensor, called a transport tensor, has a divergence which yields both advection and diffusion by the deviatory velocities. Advantages and disadvantages of Lagrangian and Eulerian averages are discussed, and meridional-vertical velocity correlations are provided.

Danielsen, E. F.↗

The effect of dissipative inhomogeneous medium on the statistics of the wave intensity

One of the main theoretical points in the theory of wave propagation in random medium is the derivation of closed form equations to describe the statistics of the propagating waves. In particular, in one dimensional problems, the closed form representation of the multiple scattering effects is important since it contributes in understanding such problems like wave localization, backscattering enhancement, and intensity fluctuations. In this the propagation of plane waves in a layer of one-dimensional dissipative random medium is considered. The medium is modeled by a complex permittivity whose real part is a constant representing the absorption. The one dimensional problem is mathematically equivalent to the analysis of a transmission line with randomly perturbed distributed parameters and a single mode lossy waveguide and the results can be used to study the propagation of radio waves through atmosphere and the remote sensing of geophysical media. It is assumed the scattering medium consists of an ensemble of one-dimensional point scatterers randomly positioned in a layer of thickness L with diffuse boundaries. A Poisson impulse process with density lambda is used to model the position of scatterers in the medium. By employing the Markov properties of this process an exact closed form equation of Kolmogorov-Feller type was obtained for the probability density of the reflection coefficient. This equation was solved by combining two limiting cases: (1) when the density of scatterers is small; and (2) when the medium is weakly dissipative. A two variable perturbation method for small lambda was used to obtain solutions valid for thick layers. These solutions are then asymptotically evaluated for small dissipation. To show the effect of dissipation, the mean and fluctuations of the reflected power are obtained. The results were compared with a lossy homogeneous medium and with a lossless inhomogeneous medium and the regions where the effect of absorption is not essential were discussed.

Saatchi, Sasan S.↗

Cr, Mn, and Ca distributions for olivine in angritic systems: Constraints on the origins of Cr-rich and Ca-poor core olivine in angrite LEW87051

Angrite meteorites are a type of basaltic achondrites that are noted for their very old cyrstallization ages (4.55 b.y.) and unusual chemical and mineralogical properties. In spite of great interest, only four angrites have been found. LEW87051 is the smallest one which weighs 0.6 g. It is a porphyritic rock with coarse subhedral to euhedral olivines set in a fine-grained groundmass which clearly represents a crystallized melt. The largest uncertainty about the petrogenesis of LEW87051 is the relationship between the large olivine crystals and the groundmass. Prinz et al. suggests that olivines are xenocrysts, while McKay et al. proposed a fractional cyrstallization model based on experimental studies. However, the crystals have Cr-rich and Ca-poor cores which do not match experimental olivines. Although Jurewicz and McKay tried to explaine the zoning of the rim by diffusion, some features are not explained. There also exists a definite composition boundary of Fe(2+) and MnO between the core and the rim. To clarify the origin of these olivines, we have performed experiments using LEW87051 analogs to measure the effects of oxygen fugacity on distribution coefficients of various elements in an angritic system.

Mikouchi, T.↗

Molecular Sieve Bench Testing and Computer Modeling

The design of an efficient four-bed molecular sieve (4BMS) CO2 removal system for the International Space Station depends on many mission parameters, such as duration, crew size, cost of power, volume, fluid interface properties, etc. A need for space vehicle CO2 removal system models capable of accurately performing extrapolated hardware predictions is inevitable due to the change of the parameters which influences the CO2 removal system capacity. The purpose is to investigate the mathematical techniques required for a model capable of accurate extrapolated performance predictions and to obtain test data required to estimate mass transfer coefficients and verify the computer model. Models have been developed to demonstrate that the finite difference technique can be successfully applied to sorbents and conditions used in spacecraft CO2 removal systems. The nonisothermal, axially dispersed, plug flow model with linear driving force for 5X sorbent and pore diffusion for silica gel are then applied to test data. A more complex model, a non-darcian model (two dimensional), has also been developed for simulation of the test data. This model takes into account the channeling effect on column breakthrough. Four FORTRAN computer programs are presented: a two-dimensional model of flow adsorption/desorption in a packed bed; a one-dimensional model of flow adsorption/desorption in a packed bed; a model of thermal vacuum desorption; and a model of a tri-sectional packed bed with two different sorbent materials. The programs are capable of simulating up to four gas constituents for each process, which can be increased with a few minor changes.

Mohamadinejad, Habib↗

Transport properties of liquid metal hydrogen under high pressures

A theory is developed for the compressibility and transport properties of liquid metallic hydrogen, near to its melting point and under high pressure. The interionic force law is assumed to be of the screened Coulomb type, because hydrogen has no core electrons. The random phase approximation is used to obtain the structure factor S(k) of the system in terms of the Fourier transform of this force law. The long wavelenth limit of the structure factor S(o) is related to the compressibility, which is much lower than that of alkali metals at their melting points. The diffusion constant at the melting point is obtained in terms of the Debye frequency, using a frequency spectrum analogous with the phonon spectrum of a solid. A similar argument is used to obtain the combined shear and bulk viscosities, but these depend also on S(o). The transport coefficients are found to be about the same size as those of alkali metals at their melting points.

Brown, R. C.↗

Theoretical analysis of infrared radiation shields of spacecraft

For a system of N diffuse, gray body radiation shields which view only adjacent surfaces and space, the net radiation method for enclosures has been used to formulate a system of linear, nonhomogeneous equations in terms of the temperatures to the fourth power of each surface in the coupled system of enclosures. The coefficients of the unknown temperatures in the system of equations are expressed in terms of configuration factors between adjacent surfaces and the emissivities. As an application, a system of four conical radiation shields for a spin stabilized STARPROBE spacecraft has been designed and analyzed with respect to variations of the cone half angles, the intershield spacings, and emissivities.

Shealy, D. L.↗

Post sunset behavior of the 6300 A atomic oxygen airglow emission

A theoretical model of the 6300 A OI airglow emission was developed based on the assumptions that both the charged and neutral portions of the Earth's upper atmosphere are in steady state conditions of diffusive equilibrium. Intensities of 6300 A OI emission line were calculated using electron density true height profiles from a standard C-4 ionosonde and exospheric temperatures derived from Fabry-Perot interferometer measurements of the Doppler broadened 6300 A emission line shape as inputs to the model. Reaction rate coefficient values, production mechanism efficiencies, solar radiation fluxes, absorption cross sections, and models of the neutral atmosphere were varied parametrically to establish a set of acceptable inputs which will consistently predict 6300 A emission intensities that closely agree with intensities observed during the post-sunset twilight period by an airglow observatory consisting of a Fabry-Perot interferometer and a turret photometer. Emission intensities that can only result from the dissociative recombination of molecular oxygen ions were observed during the latter portion of the observational period. Theoretical calculations indicate that contamination of the 6300 A OI emission should be on the order of or less than 3 percent; however, these results are very sensitive to the wavelengths of the individual lines and their intensities relative to the 6300 A OI intensity. This combination of a model atmosphere, production mechanism efficiencies, and quenching coefficient values was used when the dissociative photoexcitation and direct impact excitation processes were contributing to the intensity to establish best estimates of solar radiation fluxes in the Schumann--Runge continuum and associated absorption cross sections. Results show that the Jacchia 1971 model of the upper atmosphere combined with the Ackerman recommended solar radiation fluxes and associated absorption cross sections produces theoretically calculated intensities that more closely agree with the observed intensities than all the other combinations.

Smith, R. E.↗

Analytic model for washout of HCl(g) from dispersing rocket exhaust clouds

The potential is investigated that precipitation scavenging of HCl from large solid rocket exhaust clouds may result in unacceptably acidic rain in the Cape Canaveral, Florida, area before atmospheric dispersion reduces HCl concentrations to safe limits. Several analytic expressions for HCl(g) and HCl(g + aq) washout are derived; a geometric mean washout coefficient is recommended. A previous HCl washout model is refined and applied to a space shuttle case (70 t HCl exhausted up to 4 km) and eight Titan 3 (60 percent less exhaust) dispersion cases. The vertical column density (sigma) decays were deduced by application of a multilayer Gaussian diffusion model to seven standard meteorological regimes for overland advection. The Titan 3 decays of sigma and initial rain pH differed greatly among regimes; e.g., a range of 2 pH units was spanned at x 100 km downwind and t = 2 hr. Environmentally significant pH's .5 for infrequent exposures were shown possible at X = 50 km and t 5 hr for the two least dispersive Titan 3 cases. Representative examples of downwind rainwater pH and G(X) are analyzed. Factors affecting the validity of the results are discussed.

Pellett, G. L.↗

Large-amplitude magnetic variations in quasi-parallel shocks - Correlation lengths measured by ISEE 1 and 2

Wide separations up to more than 1 earth radius between ISEE 1 and 2 during the second half of 1978 have been used to measure the correlation length of magnetic pulsations in quasiparallel shocks. When the two spacecraft were less than a few hundred km apart, magnetic oscillations measured by magnetometers on both spacecraft exhibited virtually identical waveforms, but at distances of several thousand km, the two time series of field variation showed no detailed similarily at all. The correlation coefficients of the pulsations dropped from close to 1.0 for spacecraft separations of less than 100 km to 0.2 for separations of greater than 800 km. A correlation length of several hundred km may be related to the gyroradius of return protons with energy typical of the peaks of diffuse and beam ion distributions.

Greenstadt, E. W.↗

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.↗

Plasma instability in the presence of negative ions

Ion density fluctuations are studied in a diffusion-controlled argon plasma containing, as a dilutant, two electron-attaching species, carbon dioxide and sulfur hexafluoride. It is found that the system becomes increasingly unstable as the concentration of the electron-attaching species is increased. Nonlinear mode-mode couplings have been identified and the coupling coefficients for these interactions have been computed. Turbulent fluctuations are observed to have pronounced three-dimensionality with distinct axial and azimuthal behaviors. The power spectra of these fluctuations are composed of many discrete modes. The importance of changes in the total discharge pressure is dependent on the electron-attachment cross section of the dilutant species. From these results it is shown that the ionization instability is probably responsible for the observed phenomenon.

Johnson, J. A., III↗

Role of pressure diffusion in non-homogeneous shear flows

A non-local model is presented for approximating the pressure diffusion in calculations of turbulent free shear and boundary layer flows. It is based on the solution of an elliptic relaxation equation which enables local diffusion sources to be distributed over lengths of the order of the integral scale. The pressure diffusion model was implemented in a boundary layer code within the framework of turbulence models based on both the kappa-epsilon-(bar)upsilon(exp 2) system of equations and the full Reynolds stress equations. Model computations were performed for mixing layers and boundary layer flows. In each case, the pressure diffusion model enabled the well-known free-stream edge singularity problem to be eliminated. There was little effect on near-wall properties. Computed results agreed very well with experimental and DNS data for the mean flow velocity, the turbulent kinetic energy, and the skin-friction coefficient.

Demuren, A. O.↗

Analysis of Turbulent Flow and Heat Transfer in Noncircular Passages

Previous work on turbulent heat transfer and flow in tubes was generalized and applied to flow in non-circular passages of equilateral triangular and square cross section. Expressions for eddy diffusivity that had been verified for flow and heat transfer in tubes were assumed to apply in general along lines normal to a wall. Velocity distributions, wall shear-stress distributions, and friction factors, as well as wall heat-transfer distributions, wall temperature distributions, and average heat-transfer coefficients were calculated. In addition, results from a previous analysis for axial flow between rods were compared with new experimental data. For calculating wall temperature distributions, uniform heat generation in the passage wall and uniform heat transfer at the outer surface were assumed. The application of the results is restricted to moderately small peripheral wall temperature variations. Calculations were made for Reynolds numbers from 20,000 to 900,000 and Prandtl numbers from 0.73 to 300. Results show that velocities, shear stresses, and heat transfer in the region near the corner were lower than average values and were zero at the corner. Friction factors and average Nusselt numbers were lower than in a tube.

Taylor, Maynard F.↗