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.
Computer code has been developed to provide thermodynamic and transport properties of liquid argon, carbon dioxide, carbon monoxide, fluorine, helium, methane, neon, nitrogen, oxygen, and parahydrogen. Equation of state and transport coefficients are updated and other fluids added as new material becomes available.
The thermodynamic functions, transport properties, and physical properties of both liquid and gaseous hydrogen are presented. The low temperature regime is emphasized. The tabulation of the properties of normal hydrogen in both Si and engineering units is given along with the tabulation of parahydrogen.
The deviation from a Maxwellian velocity distribution caused by excitation of the J = 2 rotational level of parahydrogen (followed by radiative decay) in a molecular hydrogen gas is investigated. It is noted that inelastic collisions deplete the high-velocity tail of the distribution, while elastic collisions tend to refill it; the resulting steady-state distribution has a small depletion in the tail which slightly reduces the cooling rate of the gas. The elastic-collision transport is approximated by a continuous process, a simplified expression is obtained for the inelastic-collision operator, and the Boltzmann equation is then solved analytically. A correction to the inelastic-collision rate is evaluated, and it is found that the relative magnitude of the effect of tail depletion on the inelastic-collision rate is strongly temperature-dependent. The critical molecular density above which the J = 2 level deexcites by superelastic collisions is shown to be a weak function of temperature.
A new method of calculating speed of sound for two-phase flow is presented. The new equation assumes no phase change during the propagation of an acoustic disturbance and assumes that only the total entropy of the mixture remains constant during the process. The new equation predicts single-phase values for the speed of sound in the limit of all gas or all liquid and agrees with available two-phase, air-water sound speed data. Other expressions used in the two-phase flow literature for calculating two-phase, metastable sound speed are reviewed and discussed. Comparisons are made between the new expression and several of the previous expressions -- most notably a triply isentropic equation as used, a triply isentropic equation as used, among others, by Karplus and by Wallis. Appropriate differences are pointed out and a thermodynamic criterion is derived which must be satisfied in order for the triply isentropic expression to be thermodynamically consistent. This criterion is not satisfied for the cases examined, which included two-phase nitrogen, air-water, two-phase parahydrogen, and steam-water. Consequently, the new equation derived is found to be superior to the other equations reviewed.
A theory is developed for isentropic metastable sound propagation in high pressure gas-liquid mixtures. Without simplification, it also correctly predicts the minimum speed for low pressure air-water measurements where other authors are forced to postulate isothermal propagation. This is accomplished by a mixture heat capacity ratio which automatically adjusts from its single phase values to approximately the isothermal value of unity needed for the minimum speed. Computations are made for the pure components parahydrogen and nitrogen, with emphasis on the latter. With simplifying assumptions, the theory reduces to a well known approximate formula limited to low pressure.
Design features of a tested catalytic converter for altering vented cryogenic parahydrogen used as a coolant on spacecraft into a para-ortho equilibrium for channeling to other cooling functions are described. The hydrogen is expected to be stored in either liquid or solid form. A high surface area Ni-on-Si catalyst was selected for tests at an operating pressure of 2 torr at a ratio of 1000 gr catalyst for a gr/sec hydrogen flow. Cylindrical and radial flow geometries were tried and measurements centered on the converter efficiencies at different operating temperatures when the converter was placed in the vent line of the H2 cooler. Efficiencies ranging from 10-100 percent were obtained for varying flow rates. Further testing is necessary to characterize the converter performance under a wider range of operating temperatures and environments.
Energy conversion processes which are potentially important in the outer planets at pressures greater than obut 0.1 bar are reviewed. Generation of buoyancy contrasts by condensation of various constituents is discussed with emphasis on the possible significance of phase changes in substances such as Si and Mg compounds at deep levels. It is demonstrated that, in the absence of nonequilibrium thermodynamic processes, strong kinetic energy generation must accompany the transport of heat out of the high temperature planetary interiors. The possibly dominant role of lagged parahydrogen conversion in the convective transport of heat at levels where T is less than 300 K is discussed. Measurements which may ultimately contribute to a better understanding of energy conversion processes are summarized.
Voyager IRIS spectra of Jupiter are used to derive zonal averages for 270- and 150-mb temperatures, as well as optical depths through the troposphere at two temperatures, ammonia concentrations near the 680-mb level, and the parahydrogen fraction near the 270-mb level. Simple modeling of an axisymmetric circulation incorporating the linear damping of perturbations from a uniform state for both winds and temperature yields results that are consistent with observed thermal wind shears and with the vertical motion field.
It is hypothesized that observations of the Doppler shifting of IR and UV absorption lines may furnish a powerful method for the detection of 4.5-9 min standing acoustic waves trapped in a wave duct beneath the Jupiter tropopause. Similarly, data could be obtained on Jupiter's thermal and density structure, and on the depth to which its zonal winds penetrate. The present model of standing oscillations in the molecular hydrogen envelope gives attention to theoretical eigenfrequencies and to such forcing functions for wave generation as coupling with turbulent and convective motions, thermal instability due to radiative transfer, wave propagation effects in a saturated atmosphere, and ortho-to-parahydrogen conversion. The small contribution that the forcing mechanisms make to velocity amplitudes implies that the Doppler shifting caused by the waves may only be resolvable through the superposition of oscillation records to enhance S/N ratios.
The cumulative effects of an ensemble of moist convective plumes on a conditionally unstable atmosphere are predicted by a model of moist convection on Jupiter in which the heating/cooling and drying/moistening of the environment occur through (1) compensating subsidence, (2) detrainment of updraft air at cloud tops, and (3) the evaporation and melting of falling condensate. Parahydrogen is transported as a passive tracer. Pure moist convective, mixed moist-dry convective, and primarily dry convective regimes are possible, depending on the assumed deep-water abundance, efficiency of condensate evaporation, and initial temperature profile.
The present study describes modifications and additions to the MIPROPS computer program for calculating the thermophysical properties of 17 fluids. These changes include adding new fluids, new properties, and a new interface to the program. The new program allows the user to select the input and output parameters and the units to be displayed for each parameter. Fluids added to the MIPROPS program are carbon dioxide, carbon monoxide, deuterium, helium, normal hydrogen, and xenon. The most recent modifications to the MIPROPS program are the addition of viscosity and thermal conductivity correlations for parahydrogen and the addition of the fluids normal hydrogen and xenon. The recently added interface considerably increases the program's utility.
The Voyager IRIS spectra of the Jovian North Equatorial Belt (NEB) hot spots are reanalyzed using a radiative transfer model which includes the full effects of anisotropic multiple scattering by clouds. The atmospheric model includes the three thermochemically predicted cloud layers, NH3, NH4SH, and H2O. Spectrally dependent cloud extinction is modeled using Mie theory and the refractive indices of NH3 ice, NH4SH ice, water, and H2O ice. The upper tropospheric temperature profile, gas abundances, height-dependent parahydrogen profile, and vertical distribution of NH3 cloud opacity are retrieved from an analysis of the far-infrared (180-1200/cm) IRIS observations. With these properties constrained, the 5-micron (1800-2300/cm) observations are analyzed to determine the atmospheric and cloud structure of the deeper atmosphere (P of greater than 1.5 bars). The results show that the abundance of water is at least 1.5 times solar with 2 times solar (0.00276 mixing ratio relative to H2) providing the best-fit to the Voyager IRIS hot spot observations.
Voyager IRIS observations of the NEB reveal longitudinal variability of 5-micron brightness temperatures of order 100 C. An anisotropic multiple scattering radiative transfer model is used to calculate synthetic spectra for comparison with the IRIS observations. Mie theory is used to model the spectral dependence of cloud extinction from 180 to 2300/cm. Cloud base locations within the model vary with assumed gas abundances according to thermochemical equilibrium. It is found that spatial variations in the abundance profiles of the condensible species, parahydrogen profiles and cloud optical depths can be used as tracers of the local and large-scale dynamics. The variation of cloud opacity is strongly correlated with the variation of relative humidity, which suggests that dynamic depletion of water vapor above the cloud forming level is the most plausible model to explain the spatial variation in the water profile within the NEB.
This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.
The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.
While oxygen will probably be the first fluid liquefied from in-situ developed propellants, hydrogen will also be a sought-after commodity on both the Lunar and Martian surfaces. General hydrogen liquefaction processes are not significantly affected by the atmosphere on Mars and similar processes can be applied on both surfaces. However, hydrogen will require different processes than oxygen (which can be directly liquefied into its final use tank) due to the different spin states the hydrogen nucleus (i.e. parahydrogen and orthohydrogen) can have as well as the significant sensible energy that has to be removed prior to the actual phase change. Due to the low efficiency of 20 K class cryocoolers and high sensible heat removal requirements at higher temperatures, significant benefits can be obtained by precooling the hydrogen gas in a variety of manners. The benefits of three concepts are explored to take advantage of different precooling technics as well as accounting for sensible, latent, and para to ortho conversion heating requirements. Process concepts are shown that can reduce the base input power requirements for heat removal by more than 67% and the radiator area by similar amounts. A majority of the technology development for these processes has already begun, both within NASA and externally, with the few gaps being identified.
While oxygen will probably be the first fluid liquefied from in-situ developed propellants, hydrogen will also be a sought-after commodity on both the Lunar and Martian surfaces. General hydrogen liquefaction processes are not significantly affected by the atmosphere on Mars and similar processes can be applied on both surfaces. However, hydrogen will require different processes than oxygen (which can be directly liquefied into its final use tank) due to the different spin states the hydrogen nucleus (i.e. parahydrogen and orthohydrogen) can have as well as the significant sensible energy that has to be removed prior to the actual phase change. Due to the low efficiency of 20 K class cryocoolers and high sensible heat removal requirements at higher temperatures, significant benefits can be obtained by precooling the hydrogen gas in a variety of manners. The benefits of three concepts are explored to take advantage of different precooling technics as well as accounting for sensible, latent, and para to ortho conversion heating requirements. Process concepts are shown that can reduce the base input power requirements for heat removal by more than 67% and the radiator area by similar amounts. A majority of the technology development for these processes has already begun, both within NASA and externally, with the few gaps being identified.