SEARCH · Search NASA
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.
Bifurcation approximation for binary and thermal diffusion coefficients, mixture viscosity, and thermal conductivity of laminar boundary layer
Ambipolar diffusion is discussed for the case of a weakly-ionized, multicomponent plasma. Shortcomings of some of the previous formulations of this problem are pointed out. In particular, it is shown that the polarization electrostatic field acts to couple the motion of the various ions with the result that the diffusion velocity of each ion depends on the density gradients of all the ions.
A description is presented of the results of a research program directed at an improved understanding of condensate deposition rate phenomena in combustion systems. The conducted experiments make use of real-time optical laser reflectance-interference-polarization techniques in flame environments. The obtained new data and the results of previous gravimetric experiments are employed as a basis for the development of a comprehensive convective diffusion deposition theory, taking into account the assumption of a multicomponent vapor or multisize class particles 'source-free' boundary layer. The theory makes it possible to provide self-consistent salt/ash/soot deposition rate predictions over a wide variety of environmental conditions.
A flux-difference split explicit finite-difference algorithm is presented for solving the parabolized form of the equations governing three-dimensional nonequilibrium chemically reacting flows. The algorithm is based on an explicit noniterative, upwind space-marching scheme developed by Korte, but differs in that the unsteady Riemann problem, rather than the steady Riemann problem, is solved. The algorithm allows either a second or an approximately third-order accurate upwind treatment of the convection terms by employing the unsteady approximate Riemann solver of Roe. The source terms of the species transport equations are treated in either an explicit or implicit manner, and the species diffusion terms are modeled with either a Fickian or a multicomponent model. A validation of the algorithm is performed by comparing computational results with the 2-D Mach 14, 15 degree compression-corner data of Holden. The three-dimensional capability of the algorithm is demonstrated by computing Mach 2.7 flow over a swept wedge scramjet fuel injector, and three-dimensional reacting flow capability is demonstrated by a computing a shock-jet interaction concept for mixing and combustion enhancement.
Experimental and theoretical studies of the GeSe-GeI4 system revealed the multicomponent, multireaction nature of this system and the existence of a diffusion boundary layer; they also yielded equations for the prediction of diffusion limited mass transport rates. Vapor transport and crystal growth experiments of the GeSe-GeI4 and GeSe-Xe systems on earth and in a microgravity environment demonstrated the effects of convection and of microgravity on mass flux and crystal morphology. The surface and bulk morphology of GeSe crystals obtained for both systems in microgravity is considerably improved relative to ground-based specimens. Unexpected crystal growth phenomena were observed for the GeSe-Xenon system on recent Shuttle flights.
With the acquisition of satellite data on the energy spectrum of galactic gamma-radiation, it is clear that such radiation has a multicomponent nature. A calculation of the pulsar gamma ray emission spectrum is used together with a statistical analysis of recent data on 328 known pulsars to make a new determination of the pulsar contribution to galactic gamma ray emission. The contributions from diffuse interstellar cosmic ray induced production mechanisms to the total emission are then reexamined. It is concluded that pulsars may account for a significant fraction of galactic gamma ray emission.
A model for predicting the distribution of liquid fuel droplets and fuel vapor in premixing-prevaporizing fuel-air mixing passages of the direct injection type is reported. This model consists of three computer programs; a calculation of the two dimensional or axisymmetric air flow field neglecting the effects of fuel; a calculation of the three dimensional fuel droplet trajectories and evaporation rates in a known, moving air flow; a calculation of fuel vapor diffusing into a moving three dimensional air flow with source terms dependent on the droplet evaporation rates. The fuel droplets are treated as individual particle classes each satisfying Newton's law, a heat transfer, and a mass transfer equation. This fuel droplet model treats multicomponent fuels and incorporates the physics required for the treatment of elastic droplet collisions, droplet shattering, droplet coalescence and droplet wall interactions. The vapor diffusion calculation treats three dimensional, gas phase, turbulent diffusion processes. The analysis includes a model for the autoignition of the fuel air mixture based upon the rate of formation of an important intermediate chemical species during the preignition period.
In this paper, a new methodology to determine kinetic parameters for simple chemical models and simple transport properties classically used in DNS of premixed combustion is presented. First, a one-dimensional code is utilized to performed steady unstrained laminar methane-air flame in order to verify intrinsic features of laminar flames such as burning velocity and temperature and concentration profiles. Second, the flame response to steady and unsteady strain in the opposed jet configuration is numerically investigated. It appears that for a well determined set of parameters, one- and two-step mechanisms reproduce the extinction limit of a laminar flame submitted to a steady strain. Computations with the GRI-mech mechanism (177 reactions, 39 species) and multicomponent transport properties are used to validate these simplified models. A sensitivity analysis of the preferential diffusion of heat and reactants when the Lewis number is close to unity indicates that the response of the flame to an oscillating strain is very sensitive to this number. As an application of this methodology, the interaction between a two-dimensional vortex pair and a premixed laminar flame is performed by Direct Numerical Simulation (DNS) using the one- and two-step mechanisms. Comparison with the experimental results of Samaniego et al. (1994) shows a significant improvement in the description of the interaction when the two-step model is used.
A new correlation method for binary gaseous diffusion coefficients from very low temperatures to 10,000 K is proposed based on an extended principle of corresponding states, and having greater range and accuracy than previous correlations. There are two correlation parameters that are related to other physical quantities and that are predictable in the absence of diffusion measurements. Quantum effects and composition dependence are included, but high-pressure effects are not. The results are directly applicable to multicomponent mixtures.
Effective numerical treatment of multicomponent viscous flow problems involving the advection of sharp interfaces between materials of differing physical properties requires correction techniques to prevent spurious diffusion and dispersion. We develop a particular algorithm, based on modern shock-capture techniques, employing a two-step nonlinear method. The first step involves the global application of a high-order upwind scheme to a hyperbolic advection equation used to model the distribution of distinct material components in a flow field. The second step is corrective and involves the application of a global filter designed to remove dispersion errors that result from the advection of discontinuities (e.g., material interfaces) by high-order, minimally dissipative schemes. The filter introduces no additional diffusion error. Nonuniform viscosity across a material interface is allowed for by the implementation of a compositionally weighted-inverse interface viscosity scheme. The combined method approaches the optimal accuracy of modern shock-capture techniques with a minimal increase in computational time and memory. A key advantage of this method is its simplicity to incorporate into preexisting codes be they finite difference, element, or volume of two or three dimensions.
UV extinction curves that represent observed variations i diffuse clouds (Zeta Oph, Xi Per, and Sigma Sco) have been derived by combining IUE and Voyager data and have been modeled by a multicomponent size distribution of silicate and carbonaceous giants. We have used theoretical model atmospheres as our reddening-free standards, and, in the case of Zeta Oph, also compare with the use of a stellar standard. An attempt was made to classify the UV spectra via theoretical spectral indices, but still missing line opacity appears to limit the use of this method at the present. The modeling includes five possible grain ingrediants, largge (0.015-0.25 micrometers) and small (0.005 micrometers) silicate grains, similarly large and small carbonaceous grains, and polycyclic aromatic hydrocarbons (PAHs). The observed extinction curves are fitted by the model via nonlinear Chi-square-minimization, varying the fractional amount of cosmic Si and C that are locked up in each igredient. The observed variation from 'steep' to 'flat' UV extinction curves is found to correspond to a modest model variation in the lower limit of the power-law size distribution of the large grains. Graphite is found to be necessary grain ingredient that cannot be replaced by the laboratory amorphous carbon materials. Carbonaceous mantles on the larger silicate grains provide no good fits to the observed curves. Using absorption properties of neutral laboratory PAHs limits PAHs to a few percent of the cosmic C abundance for acceptable fits that also give values for R(sub V) close to the observed values.
This investigation addresses the problem of hydrodynamic stability of a two-component droplet undergoing spherically-symmetrical gasification. The droplet components are assumed to have characteristic liquid species diffusion times that are large relative to characteristic droplet surface regression times. The problem is formulated as a linear stability analysis, with a goal of predicting when spherically-symmetric droplet gasification can be expected to be hydrodynamically unstable from surface-tension gradients acting along the surface of a droplet which result from perturbations. It is found that for the conditions assumed in this paper (quasisteady gas phase, no initial droplet temperature gradients, diffusion-dominated gasification), surface tension gradients do not play a role in the stability characteristics. In addition, all perturbations are predicted to decay such that droplets were hydrodynamically stable. Conditions are identified, however, that deserve more analysis as they may lead to hydrodynamic instabilities driven by capillary effects.
A computer analysis was developed for calculating steady (or unsteady) three-dimensional aircraft component flow fields. This algorithm, called ENS3D, can compute the flow field for the following configurations: diffuser duct/thrust nozzle, isolated wing, isolated fuselage, wing/fuselage with or without integrated inlet and exhaust, nacelle/inlet, nacelle (fuselage) afterbody/exhaust jet, complete transport engine installation, and multicomponent configurations using zonal grid generation technique. Solutions can be obtained for subsonic, transonic, or hypersonic freestream speeds. The algorithm can solve either the Euler equations for inviscid flow, the thin shear layer Navier-Stokes equations for viscous flow, or the full Navier-Stokes equations for viscous flow. The flow field solution is determined on a body-fitted computational grid. A fully-implicit alternating direction implicit method is employed for the solution of the finite difference equations. For viscous computations, either a two layer eddy-viscosity turbulence model or the k-epsilon two equation transport model can be used to achieve mathematical closure.
The International Space Station (ISS) is currently equipped with a complex, heavy, and power consuming system that recovers approximately 50% of O2 from metabolic CO2. Future long duration human missions to the Moon and Mars will necessitate a sustainable and highly efficient metabolic oxygen recovery system capable of yielding a minimum of 75% O2 recovery. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to significantly increase current metabolic O2 recovery efficiency, expand mission sustainability, and reduce complexity of the system. The novel design combines CO2 conversion to O2 along with C2H4 as byproduct and water electrolysis (currently conducted in two separate units) into a single compact unit that runs at standard conditions and is theoretically capable of generating O2 with a theoretical maximum metabolic CO2 conversion of 73% while consuming less than metabolic water. This paper presents a comprehensive multi-physic 3D model developed at MSFC on CO2 conversion to O2 and C2H4 at standard conditions via MFECR. The 3D spatial domain of the model is a replica of the actual MFECR’s 3D drawing generated for the MFECR fabrication and operated to recover O2 from CO2 yielding C2H4 as byproduct. Electrochemical (EC) physics that includes EC multicomponent reaction mechanisms, mass transport, and electrical current density distributions is coupled in the model with all the other physics phenomena involved in the MFECR’s process, such as two-phase flow, free and porous fluid regimes, multicomponent mass transfer, heat transfer, and DC electrical current generation along with Joule heating effect. The EC reaction sections of the MFECR consists of two porous gas diffusion electrodes (GDE) and an electrolyte serpentine channel sandwiched in the middle. The CO2 feeds the cathode serpentine channel and part of the O2 product is fed back to the anode serpentine chamber. An alkaline solution feeds the electrolyte serpentine chamber wetting the GDEs of both, the anode and cathode allowing the OH- ionic transport between them. The EC reactions in the cathode’s GDE yield C2H4 from CO2 and H2 from water while the EC reaction in the anode’s GDE yields O2. The authors will present in this paper the validation of the model using experimental data and the utilization of the validated model in building a reliable simulator that will not only assist the authors on the MFECR design but also the optimization of its operation in the ISS and future spatial human missions.
The behavior of an isolated, stationary, multicomponent gas bubble in a glassmelt containing several dissolved gases is considered. The relevant mass-transport equations are formulated and calculations are performed for the case of two diffusing gases using a quasi-stationary model and a numerical solution of the exact mass-transfer equations. The results obtained from these two approaches are compared. The factors which govern the dissolution or growth of a bubble are thermodynamic and kinetic in origin. The tendency of a bubble to grow or shrink at long times is controlled by departure from overall equilibrium, whereas the short-time bubble dynamics may be dominated by kinetic effects. As a result of the existence of these dual influences, maxima and/or minima occur in the functional dependence of the bubble radius on time.
A nonlinear perturbation theory is formulated for the solution of the multicomponent equations of energy, mass, and momentum conservation in the atmosphere. The theory is three-dimensional and includes the effects of heat conduction and advection, viscosity, ion drag, and diffusion. The theory is described as a superposition of mathematical modes obtained by expanding the physical quantities into vector and spherical harmonics. The coupling between the various modes, both linear and nonlinear, is included. The theory provides a basis for the treatment of the thermosphere and its interaction with the lower atmosphere, where 'mode coupling' is most important. As an example, a comparison is presented between one-dimensional and three-dimensional calculations of the fundamental mode of the diurnal component in the thermosphere. Coupling between the lowest modes is considered to describe the physical conditions of the lower thermosphere where inertia and Coriolis forces become dominant over the ion-drag and viscous forces. In this region, the latitude structures of the temperature, wind field, and diffusively controlled oxygen are shown to change significantly.
A theoretical model is presented which accounts for the electroepitaxial growth kinetics and composition of multicomponent compounds in terms of mass transport in the liquid and phase diagram relationships. The mass transport in the interface is dominated by electromigration in the absence of convection and by diffusion in the presence of convection. The composition of the solid is controlled by the Peltier effect at the growth interface and by the diffusion and mobility constants of the solute components and the growth velocity (current density). Thus, for a given solution composition, the composition of the solid can be varied by varying the current density. For a given current density the composition remains constant even in the case of relatively thick epitaxial layers. All aspects of the model were found to be in good agreement with the growth and composition characteristics of Ga/x-1/Al/x/As layers.
There is increased interest in, and concern about, deposition and corrosion phenomena in combustion systems containing inorganic condensible vapors and particles (salts, ash). To meet the need for a computationally tractable deposition rate theory general enough to embrace multielement/component situations of current and future gas turbine and magnetogasdynamic interest, a multicomponent chemically 'frozen' boundary layer (CFBL) deposition theory is presented and its applicability to the special case of Na2SO4 deposition from seeded laboratory burner combustion products is demonstrated. The coupled effects of Fick (concentration) diffusion and Soret (thermal) diffusion are included, along with explicit corrections for effects of variable properties and free stream turbulence. The present formulation is sufficiently general to include the transport of particles provided they are small enough to be formally treated as heavy molecules. Quantitative criteria developed to delineate the domain of validity of CFBL-rate theory suggest considerable practical promise for the present framework, which is characterized by relatively modest demands for new input information and computer time.