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.
Report presents data on thermo-dynamic properties of 50 chemical elements, isotope (deuterium) of one of elements, and electron gas, all in standard reference states. Data tabulated as functions of temperature; also given in form of least-squares-fit coefficients of two functional forms for heat capacities at constant pressure in standard states, with constants of integration for enthalpy and entropy.
This paper evaluates a thermodynamic ice crystal icing model that has been previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. The model functions between two distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. Freeze-dominated icing occurs when liquid water (from melted ice crystals) freezes and accretes on a surface along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized as having strong adhesion to the surface. The amount of ice accretion is partially dictated by a freeze fraction, which is the fraction of impinging liquid water that freezes. Melt-dominated icing occurs as unmelted ice on a surface accumulates. This melt-dominated icing is characterized by weakly bonded surface adhesion. The amount of ice accumulation is partially dictated by a melt fraction, which is the fraction of impinging ice crystals that melts. Experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. The fraction of mass loss must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface for a given ice growth rate. This mass loss parameter, however, along with the freeze fraction and melt fraction, are the only experimental parameters that are currently not measured directly. Using icing growth rates from ice crystal icing experiments, a methodology that has been previously proposed is used to determine these unknown parameters. This work takes ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) at the Glenn Research Center in 2018 that examined the fundamental physics of ice crystal icing. This paper continues evaluation of the thermodynamic model from a previous effort, with additions to the model that account for sub-freezing temperatures that have been observed at the leading edge of the airfoil during icing. The predicted temperatures were generally in good agreement with measured temperatures. Other key findings include the total wet-bulb temperature being a good first order indicator of whether icing is freeze-dominated (sub-freezing values) or melt-dominated (above freezing). Maximum sticking efficiency values, the fraction of impinging mass that adheres to a surface, was calculated to be about 0.2, and retained this maximum value for a range of melt ratios (0.3 to 0.65 and possibly higher), which is defined as the ratio of liquid water content to total water content. Higher air velocities reduced the maximum sticking efficiency and shifted the icing regime to higher melt ratio values. Finally, the leading edge ice accretion angle was found to be related to ice growth (lower growth rates for smaller angles) and melt ratio (smaller melt ratios resulted in smaller angles, likely due to erosion effects).
For computational simulation of the convection and mixing of stratified fluid in a rotating tank (such as used in Apollo flights) with time-dependent rotation, the Navier-Stokes convection problem was formulated for a circular tank configuration. The final equations results from a general approximate theory for combined forced and contained natural convection in a time-dependent rotating system. The equations are cast in terms of vorticity and stream function in a form convenient for computation, with a transformed coordinate system, and appropriate boundary conditions are derived. Accurate representations for the cryogenic supercritical oxygen thermodynamic properties are used in the computations, and an efficient numerical finite difference scheme and computational method are employed.
The equilibrium properties of a two-dimensional plasma are examined theoretically within the framework of the random-phase approximation. The guiding-center case is of most interest because of the formal equivalence to two-dimensional discrete-line vortex fluids and because of the formal admissibility of negative temperatures. Using a judicious treatment of the periodic boundary conditions, extensive thermodynamics for nonnegative temperatures and a new value for the energy at which the temperature becomes negative are found. An explicit form for the structure function is derived and compared with the result of a Monte Carlo calculation. Good agreement is found for the negative-temperature-threshold prediction and for the general shape of the structure function. An equilibrium energy spectrum is calculated, which for positive temperatures, is the Debye-Hueckel result; for negative temperatures, it predicts an enhanced excitation of the lowest mode, corresponding to macroscopic charge separation.
A model free energy is presented which accurately represents results from 45 high-precision Monte Carlo calculations of the thermodynamics of hydrogen-helium mixtures at pressures of astrophysical and planetophysical interest. The free energy is calculated using free-electron perturbation theory (dielectric function theory), and is an extension of the expression given in an earlier paper in this series. However, it fits the Monte Carlo results more accurately, and is valid for the full range of compositions from pure hydrogen to pure helium. Using the new free energy, the phase diagram of mixtures of liquid metallic hydrogen and helium is calculated and compared with earlier results. Sample results for mixing volumes are also presented, and the new free energy expression is used to compute a theoretical Jovian adiabat and compare the adiabat with results from three-dimensional Thomas-Fermi-Dirac theory. The present theory gives slightly higher densities at pressures of about 10 megabars.
The behavior of sea ice is a strong function of ice strength which, in the central Arctic, depends on the relative fractions of multiyear and first-year ice coverage. A dynamic-thermodynamic sea ice model has been configured to include both multiyear and first-year ice. The model has been run for a 7-year period, 1979-1985, for which a satellite microwave record of multiyear sea ice coverage is available. The simulated concentrations of multiyear ice are slightly higher than the satellite-derived concentrations, although the total coverage of simulated multiyear ice agrees well with the satellite-derived coverage. Both the simulated and the satellite-derived coverages in the early 1980s show large interannual variations that represent potentially large perturbations in the mass balance of Arctic pack ice. These fluctuations are generally consistent in the model and satellite data, and they are generally consistent with year-to-year differences in the large-scale wind forcing. In at least some cases the wintertime distribution of multiyear ice is influenced by the pattern of ice drift as far back as the previous summer. The late December coverage of multiyear ice in the Alaskan region appears to foreshadow summer ice conditions in the same region. The short-term response to an intense synoptic system is considerably larger in the satellite data than in the model results.
The design of new high-energy density molecules as candidates for jet and rocket fuels is an important goal of modern chemical thermodynamics. The NASA Glenn Research Center is home to a database of thermodynamic data for over 2000 compounds related to this goal, in the form of least-squares fits of heat capacities, enthalpies, and entropies as functions of temperature over the range of 300 - 6000 K. The chemical equilibrium with applications (CEA) program written and maintained by researchers at NASA Glenn over the last fifty years, makes use of this database for modeling the performance of potential rocket propellants. During its long history, the NASA Glenn database has been developed based on experimental results and data published in the scientific literature such as the standard JANAF tables. The recent development of efficient computational techniques based on quantum chemical methods provides an alternative source of information for expansion of such databases. For example, it is now possible to model dissociation or combustion reactions of small molecules to high accuracy using techniques such as coupled cluster theory or density functional theory. Unfortunately, the current applicability of reliable computational models is limited to relatively small molecules containing only around a dozen (non-hydrogen) atoms. We propose to extend the applicability of coupled cluster theory- often referred to as the 'gold standard' of quantum chemical methods- to molecules containing 30-50 non-hydrogen atoms. The centerpiece of this work is the concept of local correlation, in which the description of the electron interactions- known as electron correlation effects- are reduced to only their most important localized components. Such an advance has the potential to greatly expand the current reach of computational thermochemistry and thus to have a significant impact on the theoretical study of jet and rocket propellants.
The paper presents a unified structural and thermodynamic theory of simple classical fluids in which the interactions between the particles can be represented by spherically symmetric pairwise potentials. Both the excess entropy and the gross features of the radial distribution function are determined mainly by excluded volume effects, which are in turn governed by a hard-core property intrinsic to any pair potential. The potential beyond this effective hard core is considered relatively weak and can be treated as a perturbation. It is also considered essential to sum all subclasses of diagrams to infinite order. The basic form of a diagrammatic scheme which allows both summation and the determination of the structure and excess entropy is presented, and a statement of universality is derived. The statement has been confirmed within the accuracy of present-day computer simulations, and a possible procedure for calculating both the structure and thermodynamics of every physically conceivable pair potential is presented.
Closed expression for volume as function of pressure and temperature has been verified for over 250 liquids. Equation can assist chemical engineers, soild-state researchers, and others with interest in thermodynamic behavior of liquids.
For a two-fluid model of dilute solutions of He-3 in liquid He-4, a thermodynamic potential is constructed that provides a Lagrangian for deriving equations of motion by a variational procedure. This Lagrangian is defined for uniform velocity fields as a (negative) Legendre transform of total internal energy, and its primary independent variables, together with their thermodynamic conjugates, are identified. Here, similarities between relations in classical physics and quantum statistical mechanics serve as a guide for developing an alternate expression for this function that reveals its character as the difference between apparent kinetic energy and intrinsic internal energy. When the He-3 concentration in the mixtures tends to zero, this expression reduces to Zilsel's formula for the Lagrangian for pure liquid He-4. An investigation of properties of the intrinsic internal energy leads to the introduction of intrinsic chemical potentials along with other intrinsic variables for the mixtures. Explicit formulas for these variables are derived for a noninteracting elementary excitation model of the fluid. Using these formulas and others also derived from quantum statistical mechanics, another equivalent expression for the Lagrangian is generated.
A multiwavelength scanning radiometer was used to measure the angular distribution of scattered radiation deep within a cloud layer at discrete wavelengths between 0.5 and 2.3 microns. The relative angular distribution of the intensity field at each wavelength is used to determine the similarity parameter, and hence single scattering albedo, of the cloud at that wavelength using the diffusion domain method. In addition to the spectral similarity parameter, the analysis provides a good estimate of the optical thickness of the cloud beneath the aircraft. In addition to the radiation measurements, microphysical and thermodynamic measurements were obtained from which the expected similarity parameter spectrum was calculated using accepted values of the refractive index of liquid water and the transmission function of water vapor. An analysis is presented for the results obtained for a 50 km section of clean marine stratocumulus clouds on 10 July 1987. These observations were obtained off the coast of California from the University of Washington Convair C-131A aircraft as part of the First ISCCP Regional Experiment (FIRE). A comparison of the experimentally-derived similarity parameter spectrum with that expected theoretically from the cloud droplet size distribution measured simultaneously from the aircraft is presented. The measurements and theory are in very close agreement for this case of clean maritime clouds.
We have developed a fast-response Nanometer Aerosol Size Analyzer (nASA) that is capable of scanning 30 size channels between 3 and 100 nm in a total time of 3 seconds. The analyzer includes a bipolar charger (P0210), an extended-length Nanometer Differential Mobility Analyzer (Nano-DMA), and an electrometer (TSI 3068). This combination of components provides particle size spectra at a scan rate of 0.1 second per channel free of uncertainties caused by response-time-induced smearing. The nASA thus offers a fast response for aerosol size distribution measurements in high-concentration conditions and also eliminates the need for applying a de-smearing algorithm to resulting data. In addition, because of its thermodynamically stable means of particle detection, the nASA is useful for applications requiring measurements over a broad range of sample pressures and temperatures. Indeed, experimental transfer functions determined for the extended-length Nano-DMA using the Tandem Differential Mobility Analyzer (TDMA) technique indicate the nASA provides good size resolution at pressures as low as 200 Torr. Also, as was demonstrated in tests to characterize the soot emissions from the J85-GE engine of a T38 aircraft, the broad dynamic concentration range of the nASA makes it particularly suitable for studies of combustion or particle formation processes. Further details of the nASA performance as well as results from calibrations, laboratory tests and field applications are presented.
The ground-based research effort under this program is concerned with systematic studies of the effects of variations: (1) of the relative importance of buoyancy-driven convection, and (2) of diffusion and viscosity conditions on crystal properties. These experimental studies are supported by thermodynamic characterizations of the systems, based on which fluid dynamic parameters can be determined. The specific materials under investigation include: the GeSe-GeI4, Ge-GeI4, HgTe-HgI2, and Hg sub (1-x)Cd sub (x) Te-HgI2 systems. Mass transport rate studies of the GeSe-GeI system as a function of orientation of the density gradient relative to the gravity vector demonstrated the validity of flux anomalies observed in earlier space experiments. The investigation of the effects of inert gases on mass flux yielded the first experimental evidence for the existence of a boundary layer in closed ampoules. Combined with a thorough thermodynamic analysis, a transport model for diffusive flow including chemical vapor transport, sublimation, and Stefan flow was developed.
The complex dielectric constant of sea water is a function of salinity at 21 cm wavelength, and sea water salinity can be determined by a measurement of emissivity at 21 cm along with a measurement of thermodynamic temperature. Three aircraft and one helicopter experiments using two different 21 cm radiometers were conducted under different salinity and temperature conditions. Single or multiple ground truth measurements were used to calibrate the data in each experiment. It is inferred from these experiments that accuracies of 1 to 2%/OO are possible with a single surface calibration point necessary only every two hours if the following conditions are met--water temperatures above 20 C, salinities above 10%/OO, and level plane flight. More frequent calibration, constraint of the aircraft's orientation to the same as it was during calibration, and two point calibration (at a high and low salinity level) rather than single point calibration may give even better accuracies in some instances.
The He I lambda 6678 line of early Be stars generally shows violet (V) and red (R) emission whenever hydrogen alpha emission is present, but its use as a diagnostic has been handicapped by a poor understanding of the processes that drive it into emission. In an attempt to address this problem we obtained three series of eschelle spectra of the first two members of the singlet and triplet 2P-nD series of lambda Eri (B2e) during 1992 November 3-5 at Kitt Peak. During these observations lambda 6678 showed substantial emission variability in both the wings and central profile, providing an opportunity to compare its behavior with that of the lambda 4922, lambda 5876, and lambda 4471 lines. We found that the responses of the lines were different in several respects. Whereas the emissions in the V wings of all four lines scaled together, the R wing of the lambda 4922 line invariably responded with increased absorption whenever the R wing of lambda 6678 line showed increased emission. These same trends occurred within the central photospheric profiles. The R-wing behavior shows that much, but not all of the emission in lambda 6678 is caused by matter projected against the stellar disk. The excitation temperatures of the neighboring 2(sup 1) P transitions, lambda 6678 and lambda 4922 must be greater than and less than the photospheric continuum temperature, respectively. We have investigated departures from local thermodynamic equilibrium (LTE) for the He I spectrum in a variety of ad hoc, perturbed model atmospheres. We have found only one way to cause the source function of lambda 6678 to increase so strongly, namely, by increasing the atmospheric temperature in the line formation region to 30,000 - 40,000 K. This effect was discovered by Auer and Mihalas for O3-O4 atmospheric models, but it has not been applied to active B stars. Our models suggest that lambda 6678 emission in Be stars can be used as a sensitive monitor of localized hot spots on these stars' surfaces. The energies involved in heating the active portions of the atmosphere are too high to be produced by gravitational infall. This leaves magnetically induced flares among the few known processes on the surfaces of stars capable of sustaining this energy level.
The thermodynamic properties, density and temperature, as well as transport property parameters involving viscosity, Prandtl number (including diffusion effects), and gaseous radiation absorption coefficients have been correlated as a function of enthalpy at four pressure levels (10 (sup -1), 10 (sup 0), 10, and 10 (sup 2) atmospheres). The correlation formulas are written in a generalized form for which coefficients for a particular property and pressure level are tabulated. The correlation formulas are useful in digital computer programs for non-adiabatic viscous flow problems.
Sources of numerical noise affecting the convergence properties of the Powell's Principal Axis Method of Optimization in the NAVY/NASA gas turbine engine computer code were investigated. The principal noise source discovered resulted from loose input tolerances used in terminating iterations performed in subroutine CALCFX to satisfy specified control functions. A minor source of noise was found to be introduced by an insufficient number of digits in stored coefficients used by subroutine THERM in polynomial expressions of thermodynamic properties. Tabular results of several computer runs are presented to show the effects on program performance of selective corrective actions taken to reduce noise.
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.