Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution species”

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 253 records · Page 14

PDF approach for compressible turbulent reacting flows

The objective of the present work is to develop a probability density function (pdf) turbulence model for compressible reacting flows for use with a CFD flow solver. The probability density function of the species mass fraction and enthalpy are obtained by solving a pdf evolution equation using a Monte Carlo scheme. The pdf solution procedure is coupled with a compressible CFD flow solver which provides the velocity and pressure fields. A modeled pdf equation for compressible flows, capable of capturing shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed, and an averaging procedure is developed to provide smooth Monte-Carlo solutions to ensure convergence. Two supersonic diffusion flames are studied using the proposed pdf model and the results are compared with experimental data; marked improvements over CFD solutions without pdf are observed. Preliminary applications of pdf to 3D flows are also reported.

Hsu, A. T.↗

Catalytic methods using molecular oxygen for treatment of PMMS and ECLSS waste streams, volume 2

Catalytic oxidation has proven to be an effective addition to the baseline sorption, ion exchange water reclamation technology which will be used on Space Station Freedom (SSF). Low molecular weight, polar organics such as alcohols, aldehydes, ketones, amides, and thiocarbamides which are poorly removed by the baseline multifiltration (MF) technology can be oxidized to carbon dioxide at low temperature (121 C). The catalytic oxidation process by itself can reduce the Total Organic Carbon (TOC) to below 500 ppb for solutions designed to model these waste waters. Individual challenges by selected contaminants have shown only moderate selectivity towards particular organic species. The combined technology is applicable to the more complex waste water generated in the Process Materials Management System (PMMS) and Environmental Control and Life Support System (ECLSS) aboard SSF. During the phase 3 Core Module Integrated Facility (CMIF) water recovery tests at NASA MSFC, real hygiene waste water and humidity condensate were processed to meet potable specifications by the combined technology. A kinetic study of catalytic oxidation demonstrates that the Langmuir-Hinshelwood rate equation for heterogeneous catalysts accurately represent the kinetic behavior. From this relationship, activation energy and rate constants for acetone were determined.

Akse, James R.↗

Analysis of gas absorption to a thin liquid film in the presence of a zero-order chemical reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

Rajagopalan, S.↗

Probability density function approach for compressible turbulent reacting flows

The objective of the present work is to extend the probability density function (PDF) tubulence model to compressible reacting flows. The proability density function of the species mass fractions and enthalpy are obtained by solving a PDF evolution equation using a Monte Carlo scheme. The PDF solution procedure is coupled with a compression finite-volume flow solver which provides the velocity and pressure fields. A modeled PDF equation for compressible flows, capable of treating flows with shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed. Two super sonic diffusion flames are studied using the proposed PDF model and the results are compared with experimental data; marked improvements over solutions without PDF are observed.

Hsu, A. T.↗

Interaction of transient radiation in nongray gaseous systems

A general formulation is presented to investigate the transient radiative interaction in nongray absorbing-emitting species between two parallel plates. Depending on the desired sophistication and accuracy, any nongray absorption model from line-by-line models to the wide band model correlations can be employed in the formulation to investigate the radiative interaction. Special attention is directed to investigate the radiative interaction in a system initially at a uniform reference temperature and suddenly the temperature of the bottom plate is reduced to a lower but constant temperature. The interaction is considered for the case of radiative equilibrium as well as for combined radiation and conduction. General as well as limiting forms of the governing equations are presented and solutions are obtained numerically by employing the method of variation of parameters. Specific results are obtained for CO, CO2, H2O, and OH. The information on species H2O and OH is of special interest for the proposed scramjet engine application. The results demonstrate the relative ability of different species for radiative interactions.

Tiwari, S. N.↗

Surface-slip equations for multicomponent, nonequilibrium air flow

Equations are presented for the surface slip (or jump) values of species concentration, pressure, velocity, and temperature in the low-Reynolds-number, high-altitude flight regime of a space vehicle. These are obtained from closed-form solutions of the mass, momentum, and energy flux equations using the Chapman-Enskog velocity distribution function. This function represents a solution of the Boltzmann equation in the Navier-Stokes approximation. The analysis, obtained for nonequilibrium multicomponent air flow, includes the finite-rate surface catalytic recombination and changes in the internal energy during reflection from the surface. Expressions for the various slip quantities have been obtained in a form which can readily be employed in flow-field computations. A consistent set of equations is provided for multicomponent, binary, and single species mixtures. Expression is also provided for the finite-rate species-concentration boundary condition for a multicomponent mixture in absence of slip.

Gupta, Roop N.↗

Distribution and metabolism of quaternary amines in salt marshes

Quaternary amines such as glycine betaine (GBT) are common osmotically active solutes in much of the marine biota. GBT is accumulated by various bacteria, algae, higher plants, invertebrates, and vertebrates in response to salinity or water stresses; in some species, GBT occurs at tens to hundreds of millimolar concentrations and can account for a significant fraction of total nitrogen. Initial studies suggest that GBT is readily converted to two potential methane precursors, trimethylamine (TMA) and acetate, in anoxic sediments. TMA is apparently the most important methane precursor in surface sediments containing sulfate reducing bacteria. In salt marshes, the bulk of the methane formed may be due to the metabolism of TMA rather than other substrates. Current research is focussed on testing this hypothesis and on determining the role of quaternary amino osmoregulatory solutes in methane fluxes from marine environments. Preliminary studies have dealt with several problems: (1) determination of GBT concentrations in the dominant flora and fauna of salt marshes; (2) synthesis of radiolabelled GBT for metabolic studies; and (3) determination of fates of BGT in marine sediments using radiotracers. Both GC and HPLC techniques have been used to assay GBT concentrations in plant and animal tissues. S. alterniflora is probably the only significant source of GBT (and indirectly of methane) since the biomass and distribution of most other species is limited. Current estimates suggest that S. alterniflora GBT could account for most of the methane efflux from salt marshes.

King, Gary M.↗

Lumping and splitting: Toward a classification of mineral natural kinds

How does one best subdivide nature into kinds? All classification systems require rules for lumping similar objects into the same category, while splitting differing objects into separate categories. Mineralogical classification systems are no exception. Our work in placing mineral species within their evolutionary contexts necessitates this lumping and splitting because we classify “mineral natural kinds” based on unique combinations of formational environments and continuous temperature-pressure-composition phase space. Consequently, we lump two minerals into a single natural kind only if they: (1) are part of a continuous solid solution; (2) are isostructural or members of a homologous series; and (3) form by the same process. A systematic survey based on these criteria suggests that 2310 (~41%) of 5659 IMA-approved mineral species can be lumped with one or more other mineral species, corresponding to 667 “root mineral kinds,” of which 353 lump pairs of mineral species, while 129 lump three species. Eight mineral groups, including cancrinite, eudialyte, hornblende, jahnsite, labuntsovite, satorite, tetradymite, and tourmaline, are represented by 20 or more lumped IMA-approved mineral species. A list of 5659 IMA-approved mineral species corresponds to 4016 root mineral kinds according to these lumping criteria. The evolutionary system of mineral classification assigns an IMA-approved mineral species to two or more mineral natural kinds under either of two splitting criteria: (1) if it forms in two or more distinct paragenetic environments, or (2) if cluster analysis of the attributes of numerous specimens reveals more than one discrete combination of chemical and physical attributes. A total of 2310 IMA-approved species are known to form by two or more paragenetic processes and thus correspond to multiple mineral natural kinds; however, adequate data resources are not yet in hand to perform cluster analysis on more than a handful of mineral species. We find that 1623 IMA-approved species (~29%) correspond exactly to mineral natural kinds; i.e., they are known from only one paragenetic environment and are not lumped with another species in our evolutionary classification. Greater complexity is associated with 587 IMA-approved species that are both lumped with one or more other species and occur in two or more paragenetic environments. In these instances, identification of mineral natural kinds may involve both lumping and splitting of the corresponding IMA-approved species on the basis of multiple criteria. Based on the numbers of root mineral kinds, their known varied modes of formation, and predictions of minerals that occur on Earth but are as yet undiscovered and described, we estimate that Earth holds more than 10000 mineral natural kinds.

Philosophy of mineralogy↗

Numerical solution for perpendicular sonic hydrogen injection into a ducted supersonic airstream

This note discusses a computer program being developed to study the flow field near opposing perpendicular fuel injectors in scramjets. The MacCormack time-split, finite difference relaxation technique was used to solve the full two-dimensional compressible Navier-Stokes equations along with energy and species equations. By using this technique, a program was developed to consider the turbulent nonreacting flow of hydrogen and air in a rectangular duct. A damping term, proportional to the second derivative of pressure and temperature, was used to produce a stable solution behind the hydrogen jet in the neighborhood of the recompression shock. A case using actual conditions encountered in current scramjet design was analyzed, with results agreeing qualitatively with experimental observations.

Drummond, J. P.↗

Equilibrium properties of chemically reacting gases

The equilibrium energy, enthalpy, entropy, specific heat at constant volume and constant pressure, and the equation of state of the gas are all derived for chemically reacting gas mixtures in terms of the compressibility, the mol fractions, the thermodynamic properties of the pure gas components, and the change in zero point energy due to reaction. Results are illustrated for a simple diatomic dissociation reaction and nitrogen is used as an example. Next, a gas mixture resulting from combined diatomic dissociation and atomic ionization reactions is treated and, again, nitrogen is used as an example. A short discussion is given of the additional complexities involved when precise solutions for high-temperature air are desired, including effects caused by NO produced in shuffle reactions and by other trace species formed from CO2, H2O and Ar found in normal air.

Source record↗

Advanced numerics for multi-dimensional fluid flow calculations

In recent years, there has been a growing interest in the development and use of mathematical models for the simulation of fluid flow, heat transfer and combustion processes in engineering equipment. The equations representing the multi-dimensional transport of mass, momenta and species are numerically solved by finite-difference or finite-element techniques. However despite the multiude of differencing schemes and solution algorithms, and the advancement of computing power, the calculation of multi-dimensional flows, especially three-dimensional flows, remains a mammoth task. The following discussion is concerned with the author's recent work on the construction of accurate discretization schemes for the partial derivatives, and the efficient solution of the set of nonlinear algebraic equations resulting after discretization. The present work has been jointly supported by the Ramjet Engine Division of the Wright Patterson Air Force Base, Ohio, and the NASA Lewis Research Center.

Vanka, S. P.↗

Radiative interactions in laminar duct flows

Analyses and numerical procedures are presented for infrared radiative energy transfer in gases when other modes of energy transfer occur simultaneously. Two types of geometries are considered, a parallel plate duct and a circular duct. Fully developed laminar incompressible flows of absorbing-emitting species in black surfaced ducts are considered under the conditions of uniform wall heat flux. The participating species considered are OH, CO, CO2, and H2O. Nongray as well as gray formulations are developed for both geometries. Appropriate limiting solutions of the governing equations are obtained and conduction-radiation interaction parameters are evaluated. Tien and Lowder's wide band model correlation was used in nongray formulation. Numerical procedures are presented to solve the integro-differential equations for both geometries. The range of physical variables considered are 300 to 2000 K for temperature, 0.1 to 100.0 atm for pressure, and 0.1 to 100 cm spacings between plates/radius of the tube. An extensive parametric study based on nongray formulation is presented. Results obtained for different flow conditions indicate that the radiative interactions can be quite significant in fully developed incompressible flows.

Trivedi, P. A.↗

Best approximation to the weighted slab method

The weighted-slab method is modified so that, although it is still not exact, it gives the 'best' approximation in a minimization sense when energy loss cannot be neglected. In this approximation the species-dependent energy-change term that operates on the path-length distribution is 'averaged' over the slab-model solution for that particular energy and path length.

Jones, Frank C.↗

Evaporation and combustion of LOX under supercritical and subcritical conditions

The objective is to study the evaporation and combustion of LOX under supercritical and subcritical conditions both experimentally and theoretically. In the evaporation studies, evaporation rate and surface temperature were measured when LOX vaporizing in helium environments at pressures ranging from 5 to 68 atm. A Varian 3700 gas chromatograph was employed to measure the oxygen concentration above the LOX surface. For the combustion tests, high-magnification video photography was used to record direct images of the flame shape of a LOX/H2/He laminar diffusion flame. The gas composition in the post-flame region is also being measured with the gas sampling and chromatography analysis. These data are being used to validate the theoretical model. A comprehensive theoretical model with the consideration of the solubility of ambient gases as well as variable thermophysical properties was formulated and solved numerically to study the gasification and burning of LOX at elevated pressures. The calculated flame shape agreed reasonably well with the edge of the observed luminous flame surface. The effect of gravity on the flame structure of laminar diffusion flames was found to be significant. In addition, the predicted results using the flame-sheet model were compared with those based upon full equilibrium calculations (which considered the formation of intermediate species) at supercritical pressures. Except at the flame front where temperature exceeded 2,800 K, the flame-sheet and equilibrium solutions in terms of temperature distributions were in very close agreement. The temperature deviation in the neighborhood of the flame front is caused by the effect of high-temperature dissociation.

Yang, A. S.↗

Application of the regular associated solution model to the Cd-Te and Hg-Te binary systems

The regular associated solution model is used to treat the phase diagrams of the binary II-VI semiconductor alloy systems Hg-Te and Cd-Te. The equations for the species activity coefficients are used without approximations regarding the magnitudes of the various binary interchange energies or the functional dependence on component mol fraction. The values of the four-adjustable parameters required for description of each system are fixed by fitting liquidus data, and the resulting activity coefficients are used to calculate component partial pressures, which are compared with experimental values as an indpendent check of the validity of the model. The results show that the regular associated solution model provides a usefully accurate, but not complete, description for both the Hg-Te and Cd-Te systems. The relationship of this work to previous investigations is discussed.

Kelley, J. D.↗

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗

A finite element solution algorithm for the Navier-Stokes equations

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the steady-state kinematics and thermodynamics of a variable viscosity, compressible multiple-species fluid. For an incompressible fluid, the motion may be transient as well. The primitive dependent variables are replaced by a vorticity-streamfunction description valid in domains spanned by rectangular, cylindrical and spherical coordinate systems. Use of derived variables provides a uniformly elliptic partial differential equation description for the Navier-Stokes system, and for which the finite element algorithm is established. Explicit non-linearity is accepted by the theory, since no psuedo-variational principles are employed, and there is no requirement for either computational mesh or solution domain closure regularity. Boundary condition constraints on the normal flux and tangential distribution of all computational variables, as well as velocity, are routinely piecewise enforceable on domain closure segments arbitrarily oriented with respect to a global reference frame.

Baker, A. J.↗

Non-equilibrium radiation from viscous chemically reacting two-phase exhaust plumes

A knowledge of the structure of the rocket exhaust plumes is necessary to solve problems involving plume signatures, base heating, plume/surface interactions, etc. An algorithm is presented which treats the viscous flow of multiphase chemically reacting fluids in a two-dimensional or axisymmetric supersonic flow field. The gas-particle flow solution is fully coupled with the chemical kinetics calculated using an implicit scheme to calculate chemical production rates. Viscous effects include chemical species diffusion with the viscosity coefficient calculated using a two-equation turbulent kinetic energy model.

Penny, M. M.↗