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.
Coefficients of space charge effects for growing oxide films, derived from perturbation treatment, evaluated in terms of boundary concentrations of diffusing species
Simultaneous measurements of ambient O2 and N2 concentrations and temperature were made near winter solstice 1975 by the open source neutral mass spectrometer on Atmosphere Explorer-D. The data, taken at midmorning and all latitudes, are analyzed using scatter plots and correlation coefficients. Results indicate that O2 concentrations at about 200 km is controlled both by temperature variations and wind-induced diffusion. Since the winds are most effective in controlling O2 at lower altitudes, it is reasonable to expect that temperature effects will dominate over wind effects at higher altitudes.
The development of an advanced radiator concept design meeting SP-100 thermoelectric requirements is reported. Carbon-carbon heat pipes are used to produce this lightweight, high-performance radiator. Two feasibility problems had to be solved to enable the design: the production of a carbon-carbon heat pipe tube with integral fins, and the development of a coating that protects the carbon-carbon substrate from 875-K potassium working fluid. Carbon-carbon tubes with integral fins were successfully produced using a T-300 fiber, an angle interlocking weave architecture, and pitch densification. A barrier coating to protect the inside diameter of these tubes which employed the CVD of niobium over a thin rhenium interlayer was developed. The rhenium interlayer proved to be critical to the success of this coating technique by providing a gradation in the coefficient of thermal expansion, carrying a portion of the induced stress load, improving coating adhesion, and providing a partial carbon diffusion barrier.
The purpose of this work is to carry out direct numerical simulations of diffusion controlled combustion with non-unity Lewis numbers for the reactants and products, thus accounting for the differential diffusion effects of the temperature and concentration fields. We use a formulation based on combining the conservation equations in a way to eliminate the reaction terms similar to the method used by Burke and Schumann (1928) for unity Lewis numbers. We present calculations for an axisymmetric fuel jet and for a planar, time evolving mixing layer, leaving out the effects of thermal expansion and variations of the transport coefficients due to the heat release. Our results show that the front of the flame shifts toward the fuel or oxygen sides owing to the effect of the differential diffusion and that the location of maximum temperature may not coincide with the flame. The dependence of the distribution of the reaction products on their Lewis number has been investigated.
The intensity of turbulence and the Lagrangian correlation coefficient in a gaseous hydrogen (GH2)-gaseous oxygen (GOX) rocket combustion chamber are studied. A helium tracer gas was injected at an upstream point on the combustion chamber centerline. Samples collected downstream were analyzed for the spreading of the tracer gas concentration. Using Taylor's turbulent diffusion theory and Gosman's numerical method, turbulence parameters were found. In order to evaluate turbulence field inhomogeneity in the axial direction, Taylor's theory was extended. An exponential function was used to approximate the Lagrangian correlation coefficient. It was found that turbulence intensity varies from 15% near the injector to 4% at the nozzle entrance, and that the Lagrangian coefficient rapidly decreases when the dispersion time difference increases.
A new method of solving the radiative transfer equation is developed in which the scattering and absorption coefficients may have arbitrary variations with depth, and in which both internal (thermal) emission and incident radiation are allowed. Specular and diffuse reflection at both boundaries also is taken into account. The method begins by forming a paired set of coupled first-order differential equations for the symmetric and antisymmetric parts of the radiation field after writing the scattering integral as a numerical quadrature. These differential equations are broken into finite difference form, in which the symmetric and antisymmetric parts of the radiation field are found on alternate grid points. Numerical results for a number of test problems are shown, demonstrating that the method is very fast, that it returns specific intensities and fluxes that are accurate to at least a percent, and that it can be applied to optically thick problems.
The spectrum of sunlight reflected by Jupiter is analyzed by comparing observations of Woodman (1979) with multiple-scattering computations. The analysis yields information on the vertical cloud structure at several latitudes and on the abundance of CH4 and NH3 in the atmosphere of Jupiter. The abundances of CH4 and NH3 suggest that all ices and rocks are overabundant on Jupiter by a factor of 2 or more, providing an important constraint on models for the formation of Jupiter from the primitive solar nebula. The pressure level of the clouds, the gaseous NH3 abundance, the mean temperature profile, and the Clausius-Clapeyron relation suggest that these clouds are predominantly ammonia crystals with the cloud bottom at 600-700 mb. A diffuse distribution of aerosols exists between 150 and 500 mb, and the spectral variation of albedo reflects a changing bulk absorption coefficient of the material composing the aerosols and is diagnostic of the aerosol composition.
The Andre et al. (1976, 1978) third-order closure model, in which the time rate of change terms, the relaxation and rapid effects for pressure-related terms, and the clipping approximation are used along with the quasi-normal closure, is invoked in the study of turbulence in a cloudy layer that is radiatively cooled from above. A spurious oscillation whose greatest amplitude lies near the inversion is shown by analysis to arise from the mean gradient and buoyancy terms of the triple-moment equations. An attempt is made to damp the oscillation through the introduction of diffusion terms into the triple-moment equations. The results obtained are noted to be sensitive to the ad hoc eddy coefficient applied in the third-moment equations.
A theoretical model for computing the displacement damage defect density and the short-circuit current (I sub sc) degradation in proton-irradiated (AlGa)As-GaAs p-n junction solar cells is presented. Assumptions were made with justification that the radiation induced displacement defects form an effective recombination center which controls the electron and hole lifetimes in the junction space charge region and in the n-GaAs active layer of the irradiated GaAs p-n junction cells. The degradation of I sub sc in the (AlGa)As layer was found to be negligible compared to the total degradation. In order to determine the I sub sc degradation, the displacement defect density, path length, range, reduced energy after penetrating a distance x, and the average number of displacements formed by one proton scattering event were first calculated. The I sub sc degradation was calculated by using the electron capture cross section in the p-diffused layer and the hole capture cross section in the n-base layer as well as the wavelength dependent absorption coefficients. Excellent agreement was found between the researchers calculated values and the measured I sub sc in the proton irradiated GaAs solar cells for proton energies of 100 KeV to 10 MeV and fluences from 10 to the 10th power p/square cm to 10 to the 12th power p/square cm.
A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.
A space-based experiment is currently under development to study diffusion-controlled, gas-phase, low temperature oxidation reactions, cool flames and auto-ignition in an unstirred, static reactor. At Earth's gravity (1g), natural convection due to self-heating during the course of slow reaction dominates diffusive transport and produces spatio-temporal variations in the thermal and thus species concentration profiles via the Arrhenius temperature dependence of the reaction rates. Natural convection is important in all terrestrial cool flame and auto-ignition studies, except for select low pressure, highly dilute (small temperature excess) studies in small vessels (i.e., small Rayleigh number). On Earth, natural convection occurs when the Rayleigh number (Ra) exceeds a critical value of approximately 600. Typical values of the Ra, associated with cool flames and auto-ignitions, range from 104-105 (or larger), a regime where both natural convection and conduction heat transport are important. When natural convection occurs, it alters the temperature, hydrodynamic, and species concentration fields, thus generating a multi-dimensional field that is extremely difficult, if not impossible, to be modeled analytically. This point has been emphasized recently by Kagan and co-workers who have shown that explosion limits can shift depending on the characteristic length scale associated with the natural convection. Moreover, natural convection in unstirred reactors is never "sufficiently strong to generate a spatially uniform temperature distribution throughout the reacting gas." Thus, an unstirred, nonisothermal reaction on Earth does not reduce to that generated in a mechanically, well-stirred system. Interestingly, however, thermal ignition theories and thermokinetic models neglect natural convection and assume a heat transfer correlation of the form: q=h(S/V)(T(bar) - Tw) where q is the heat loss per unit volume, h is the heat transfer coefficient, S/V is the surface to volume ratio, and (T(bar) - Tw ) is the spatially averaged temperature excess. This Newtonian form has been validated in spatially-uniform, well-stirred reactors, provided the effective heat transfer coefficient associated with the unsteady process is properly evaluated. Unfortunately, it is not a valid assumption for spatially-nonuniform temperature distributions induced by natural convection in unstirred reactors. "This is why the analysis of such a system is so difficult." Historically, the complexities associated with natural convection were perhaps recognized as early as 1938 when thermal ignition theory was first developed. In the 1955 text "Diffusion and Heat Exchange in Chemical Kinetics", Frank-Kamenetskii recognized that "the purely conductive theory can be applied at sufficiently low pressure and small dimensions of the vessel when the influence of natural convection can be disregarded." This was reiterated by Tyler in 1966 and further emphasized by Barnard and Harwood in 1974. Specifically, they state: "It is generally assumed that heat losses are purely conductive. While this may be valid for certain low pressure slow combustion regimes, it is unlikely to be true for the cool flame and ignition regimes." While this statement is true for terrestrial experiments, the purely conductive heat transport assumption is valid at microgravity (mu-g). Specifically, buoyant complexities are suppressed at mu-g and the reaction-diffusion structure associated with low temperature oxidation reactions, cool flames and auto-ignitions can be studied. Without natural convection, the system is simpler, does not require determination of the effective heat transfer coefficient, and is a testbed for analytic and numerical models that assume pure diffusive transport. In addition, mu-g experiments will provide baseline data that will improve our understanding of the effects of natural convection on Earth.
An advanced engineering model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multispecies, chemically reacting and multidiffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and the Dufour energy transfer terms. In addition to the spectrum of multispecies aspects developed, the model developed in this study is also conservative in axisymmetric flow for both inviscid and viscous flow environments and the boundary conditions employ a viscous, chemically reacting, reference plane characteristics method. 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 and the 25 lbf film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent flow plug and spool nozzle analysis cases performed. Further, the Soret term was shown to represent an important fraction of the diffusion fluxes occurring in a transpiration cooled rocket engine.
Diffusive and radiative effects are incorporated into an analysis for the vaporization time of drops in film boiling. The momentum, energy, and continuity equations are solved with some appropriate simplifications so as to obtain a simple closed form solution for the overall film boiling heat transfer coefficient. Next, a theoretical expression for the droplet vaporization time is developed and compared to the measured vaporization times of water droplets vaporizing into air, argon, nitrogen, and helium. The agreement between experiment and theory is good. Under the helium blanket, the diffusive evaporative component is significant in comparison to the film boiling component.
A Monte Carlo modeling technique is described for mathematically simulating free molecular flows over a concave spherical surface and a concave cylindrical surface of finite length. The half-angle of the surfaces may vary from 0 to 90 degrees, and the incident flow may have an arbitrary speed ratio and an arbitrary angle of attack. Partial diffuse reflection and imperfect energy accommodation for molecules colliding with the surfaces are also considered. Results of heat transfer, drag and lift coefficients are presented for a variety of flow conditions. The present Monte Carlo results are shown to be in very good agreement with certain available theoretical solutions.
This theoretical study addresses the nature of convective motions in a toroidal loop of binary fluid oriented in the vertical plane and heated from below. The boundaries of the loop are impermeable, but gradients of the solute can be set up by Soret diffusion in the direction around the loop. The existence and stability of steady solutions are discussed over the Rayleigh number-Soret coefficient parameter plane. When the Soret coefficient is negative, periodic and chaotic oscillations analogous to those of thermohaline convection are predicted. When the Soret coefficient is positive, relaxation oscillations and low Rayleigh number chaotic motions are found. Both sets of phenomena are predicted to occur for realistic thermosyphon parameters.
A leaf is modeled as a deterministic two-dimensional structure consisting of a network of circular arcs designed to represent the internal morphology of major species. The path of an individual ray through the leaf is computed using geometric optics. At each intersection of the ray with an arc, the specular reflected and transmitted rays are calculated according to the Snell and Fresnel equations. Diffuse scattering is treated according to Lambert's law. Absorption is also permitted but requires a detailed knowledge of the spectral attenuation coefficients. An ensemble of initial rays are chosen for each incident direction with the initial intersection points on the leaf surface selected randomly. The final equilibrium state after all interactions then yields the leaf bidirectional reflectance and transmittance distributions. The model also yields the internal two dimensional light gradient profile of the leaf.
Two turning vane designs were experimentally evaluated for the fan-drive corner (corner 2) coupled to an upstream diffuser and the high-speed corner (corner 1) of the 0.1 scale model of NASA Lewis Research Center's proposed Altitude Wind Tunnel. For corner 2 both a controlled-diffusion vane design (vane A4) and a circular-arc vane design (vane B) were studied. The corner 2 total pressure loss coefficient was about 0.12 with either vane design. This was about 25 percent less loss than when corner 2 was tested alone. Although the vane A4 design has the advantage of 20 percent fewer vanes than the vane B design, its vane shape is more complex. The effects of simulated inlet flow distortion on the overall losses for corner 1 or 2 were small.
Measurements of the velocity characteristics of the flows in two curved diffusers of rectangular cross section with C and S-shaped centerlines are presented and related to measurements of wall heat transfer coefficients along the heated flat walls of the ducts. The velocity results were obtained by laser-Doppler anemometry in a water tunnel and the heat transfer results by liquid crystal thermography in a wind tunnel. The thermographic technique allowed the rapid and inexpensive measurement of wall heat transfer coefficients along flat walls of arbitrary boundary shapes with an accuracy of about 5 percent. The results show that an increase in secondary flow velocities near the heated wall causes an increase in the local wall heat transfer coefficient, and quantify the variation for maximum secondary-flow velocities in a range from 1.5 to 17 percent of the bulk flow velocity.