Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic functions”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

PREWATE: An interactive preprocessing computer code to the Weight Analysis of Turbine Engines (WATE) computer code

The Weight Analysis of Turbine Engines (WATE) computer code was developed by Boeing under contract to NASA Lewis. It was designed to function as an adjunct to the Navy/NASA Engine Program (NNEP). NNEP calculates the design and off-design thrust and sfc performance of User defined engine cycles. The thermodynamic parameters throughout the engine as generated by NNEP are then combined with input parameters defining the component characteristics in WATE to calculate the bare engine weight of this User defined engine. Preprocessor programs for NNEP were previously developed to simplify the task of creating input datasets. This report describes a similar preprocessor for the WATE code.

Fishbach, L. H.↗
NASA NTRS1991

No-vent fill pressurization tests using a cryogen simulant

The results are described of an experimental program which studied the performance of various no-vent fill techniques for tank-to-tank liquid transfer. The tests were performed using a cryogen simulant (Freon-114) and a test bed consisting of a multiple tank/plumbing network that enabled studies of a variety of different inlet flow and active mixing regimes. Several results and conclusions were drawn from the 26 transfer experiments comprising the program. Most notable was the significant improvement in fill performance (i.e., minimized fill time and maximized fill fraction) with increased agitation of the liquid surface. Another was the close correlation between measured condensation rates and those predicted by recent theories which express condensation as a function of turbulent eddy effects on the liquid surface. In most cases, test data exhibited strong agreement with an analytical model which accounts for tank heat transfer and thermodynamics in a 1 g environment.

Schmidt, G. R.↗
NASA NTRS1992

No-vent fill pressurization tests using a cryogen simulant

The results are described of an experimental program which studied the performance of various no-vent fill techniques for tank-to-tank liquid transfer. The tests were performed using a cryogen simulant (Freon-114) and a test bed consisting of a multiple tank/plumbing network that enabled studies of a variety of different inlet flow and active mixing regimes. Several results and conclusions were drawn from the 26 transfer experiments comprising the program. Most notable was the significant improvement in fill performance (i.e., minimized fill time and maximized fill fraction) with increased agitation of the liquid surface. Another was the close correlation between measured condensation rates and those predicted by recent theories which express condensation as a function of turbulent eddy effects on the liquid surface. In most cases, test data exhibited strong agreement with an analytical model which accounts for tank heat transfer and thermodynamics in a 1 g environment.

Schmidt, G.↗
NASA NTRS1994

NLTE model stellar atmospheres with line blanketing near the series limits

In this paper we study the influence of line-merging regions at the intermediate long-wavelength side of a continuum threshold on the computed model atmosphere structure and predicted spectrum. In order to model these regions sufficiently accurately, we have developed two concepts. First, we have extended the occupation probability formalism of Hummer and Mihalas to non-local thermodynamic equilibrium (non-LTE) plasmas. Second, in order to treat the very complicated opacity in the line merging region, we have generalized the concept of opacity distribution functions to treat non-LTE situations. All Rydberg states are consistently included within this framework, so that no arbitrary cutoff of high (LTE) levels is made. We have calculated several pure hydrogen models atmospheres for two effective temperatures, T(sub eff) = 20000 and 35000 K, and discussed the differences between models calculated with various treatments of the line merging. In particular, we have shown that the error in the predicted profiles of Balmer lines resulting from the neglect of line merging is typically of the order of 3-4%, while the errors in the far-UV portion of the Balmer continuum reaches 15-35%. The errors generally decrease with increasing effective temperature. At the same time, the internal accuracy of the models is shown to be about or below 0.5% for all predicted spectral features. We conclude that for interpreting current high-accuracy spectrophotometric observations models including the line merging are necessary, and that the formalism developed in this paper is capable of providing a sufficiently accurate and robust modeling technique.

Hubeny, I.↗
NASA NTRS1997

Integrins, tensegrity, and mechanotransduction

Physical forces, such as those due to gravity, play an important role in tissue development and remodeling. Yet, little is known about how individual cells sense mechanical signals or how they transduce them into a chemical response. Rather than listing the numerous signal pathways that have been found to be sensitive to mechanical stimulation, we need to place potential molecular signaling mechanisms within the context of the entire cell. The model presented is based on the concept that cells use tensegrity architecture to organize their cytoskeleton and stabilize their form. Studies with stick and string tensegrity cell models predict that living cells are hard-wired to respond immediately to external mechanical stresses. This hard-wiring exists in the form of discrete cytoskeletal filament networks that mechanically couple specific cell surface receptors, such as integrins, to nuclear matrix scaffolds and to potential transducing molecules that physically associate with the cytoskeleton. If these signaling molecules do function in a "solid-state", then mechanical stresses may be transduced into biochemical responses through force-dependent changes in cytoskeletal geometry or through local alterations in thermodynamic or kinetic parameters. Changes in cytoskeletal tension (prestress) also may play a role in signal amplification and adaptation. Recent experimental results are described which provide direct support for the tensegrity theory.

Review↗
NASA NTRS1976

Background

An analysis was made of the UF6 fueled gas core reactor as a function of cavity reactor criticality and fluid mechanics tests, investigations of uranium optical emission spectra, and radiant heat transfer power plant studies. Data are also given on nuclear and thermodynamic cycle analysis.

Source record↗
NASA NTRS1997

Kinetics and Thermochemistry of Reversible Adduct Formation in the Reaction of Cl((sup 2)P(sub J)) with CS2

Reversible adduct formation in the reaction of Cl((sup 2)P(sub J)) with CS2 has been observed over the temperature range 193-258 K by use of time-resolved resonance fluorescence spectroscopy to follow the decay of pulsed-laser-generated Cl((sup 2)P(sub J)) into equilbrium with CS2Cl. Rate coefficients for CS2Cl formation and decomposition have been determined as a function of temperature and pressure; hence, the equilbrium constant has been determined as a function of temperature. A second-law analysis of the temperature dependence of Kp and heat capacity corrections calculated with use of an assumed CS2Cl structure yields the following thermodynamic parameters for the association reaction: Delta-H(sub 298) = -10.5 +/- 0.5 kcal/mol, Delta-H(sub 0) = -9.5 +/- 0.7 kcal/mol, Delta-S(sub 298) = -26.8 +/- 2.4 cal/mol.deg., and Delta-H(sub f,298)(CS2Cl) = 46.4 +/- 0.6 kcal/mol. The resonance fluorescence detection scheme has been adapted to allow detection of Cl((sup 2)P(sub J)) in the presence of large concentrations of O2, thus allowing the CS2Cl + Cl + O2 reaction to be investigated. We find that the rate coefficient for CS2Cl + O2 reaction via all channels that do not generate Cl((sup 2)P(sub J)) is less than 2.5 x 10(exp-16) cu cm/(molecule.s) at 293 K and 300-Torr total pressure and that the total rate coefficient is less than 2 x 10 (exp -15) cu cm/(molecule.s) at 230 K and 30-Torr total pressure. Evidence for reversible adduct formation in the reaction of Cl((sup 2)P(sub J)) with COS was sought but not observed, even at temperatures as low as 194 K.

Nicovich, J. M.↗
NASA NTRS1999

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

SiC and Si3N4 materials were tested under various turbine engine combustion environments, chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high speed aircraft. Representative CVD, sintered, and composite materials were evaluated in both furnace and high pressure burner rig exposure. While protective SiO2 scales form in all cases, evidence is presented to support paralinear growth kinetics, i.e. parabolic growth moderated simultaneously by linear volatilization. The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were used to map SiO2 volatility (and SiC recession) over a range of temperature, pressure, and velocity. The functional dependency of material recession (volatility) that emerged followed the form: exp(-QIRT) * P(exp x) * v(exp y). These empirical relations were compared to rates predicted from the thermodynamics of volatile SiO and SiO(sub x)H(sub Y) reaction products and a kinetic model of diffusion through a moving, boundary layer. For typical combustion conditions, recession of 0.2 to 2 micron/h is predicted at 1200- 1400C, far in excess of acceptable long term limits.

Smialek, James L.↗
NASA NTRS2012

The Role of Aerosols on Precipitation Processes: Cloud Resolving Model Simulations

Cloud microphysics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e. pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e., 33 bins). Atmospheric aerosols are also described using number density size-distribution functions. The model is tested by studying the evolution of deep cloud systems in the west Pacific warm pool region, the sub-tropics (Florida) and midlatitudes using identical thermodynamic conditions but with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. Results indicate that the low CCN concentration case produces rainfall at the surface sooner than the high CeN case but has less cloud water mass aloft. Because the spectral-bin model explicitly calculates and allows for the examination of both the mass and number concentration of species in each size category, a detailed analysis of the instantaneous size spectrum can be obtained for these cases. It is shown that since the low (CN case produces fewer droplets, larger sizes develop due to greater condensational and collection growth, leading to a broader size spectrum in comparison to the high CCN case. Sensitivity tests were performed to identify the impact of ice processes, radiation and large-scale influence on cloud-aerosol interactive processes, especially regarding surface rainfall amounts and characteristics (i.e., heavy or convective versus light or stratiform types). In addition, an inert tracer was included to follow the vertical redistribution of aerosols by cloud processes. We will also give a brief review from observational evidence on the role of aerosol on precipitation processes.

Tao, Wei-Kuo↗
NASA NTRSDate not supplied

Analysis of Thermoradiative Thermal Energy Conversion

The thermoradiative cell is a new method for converting heat energy to electrical power, first detailed by Strandberg in 2015. The cell is a p-n junction semiconductor device, similar to a photovoltaic cell but thermodynamically operating in the reverse direction, converting the thermal dark current into electrical power while radiating waste heat to space. The power and efficiency can be calculated as a function of bandgap in the Shockley-Queisser detailed-balance limit, in which the thermal emissivity of the cell is due to the recombination of electron-hole pairs, and all other recombination losses are ignored. The current produced is directly proportional to the recombination radiation. The fundamental loss mechanism for the thermoradiative cell is the energy carried by the infrared radiation into space from band-to-band recombination of carriers injected across the junction. In an ideal cell, to maximize the efficiency, the emission energy of these photons would precisely equal the bandgap. This can be achieved, for example, using dielectric filters or meta-material filters to recycle emission at other wavelengths back into the cell. The voltage is proportional to the external bias. These two constraints allow optimization of the optimum bias point for maximum power. Unlike photovoltaic cells, the maximum power operating point is not the same as the maximum efficiency point, and higher efficiency can be achieved at a higher (negative) bias in the ideal case. Incorporating non-ideal losses, however, shifts the maximum efficiency point toward lower bias. Unlike in photovoltaic cells, non-radiative recombination (e.g., Auger losses) will reduce the output current, but will not reduce the conversion efficiency, since the recombination energy is retained in the cell in the form of heat. Since a thermoradiative cell operates by radiating directly to space, the current produced by themoradiative cells will increase as Stefan-Boltzmann radiation; roughly the fourth power of the temperature. Thus, the power produced is highest at high operating temperatures, and, unlike conventional thermal conversion, increasing radiator temperature increases, the efficiency. Thus, the choice of technology will be toward semiconductors resistant to degradation at high temperature.

Thermoradiative↗
NASA NTRS1988

A calorimetric determination of the enthalpy of formation and a description of the defect structure of the ordered beta-phase /Ni, Cu/ /1-x/ Al/x/

In order to describe thermodynamically the defect structure of an ordered B-Hume-Rothery phase, the heat of formation of (Ni,Cu)(1-x)Al(x) was measured at 1100 K as a function of concentration in the range x (sub Al) = 0.4 and 0.55 for three substitution rations x (sub Ni)/x (sub Cu) = infinity; 11; 5. The heat of formation of the NiAl beta-phase is strongly negative. For the stoichiometric composition it is -72.2 kJ/g-atom. On both the nickel-rich side and the aluminum-rich side the magnitude of the enthalpy of formation decreases linearly with concentration. Substitution of nickel for copper decreases the magnitude of the enthalpy of formation over the entire homogeneity range for the phase (Ni,Cu)(1-x)Al(x). The curve for the enthalpy of formation as well as the literature values for the chemical potential of aluminum are described with great accuracy by the disorder model of Wagner-Schottky.

Henig, E. T.↗
NASA NTRS1993

An estimate of the second law thermodynamic efficiency of the various units comprising an Environmental Control and Life Support System (ECLSS)

The objective of this paper is to present an estimate of the second law thermodynamic efficiency of the various units comprising an Environmental Control and Life Support System (ECLSS). The technique adopted here is based on an evaluation of the 'lost work' within each functional unit of the subsystem. Pertinent information for our analysis is obtained from a user interactive integrated model of an ECLSS. The model was developed using ASPEN. A potential benefit of this analysis is the identification of subsystems with high entropy generation as the most likely candidates for engineering improvements. This work has been motivated by the fact that the design objective for a long term mission should be the evaluation of existing ECLSS technologies not only the basis of the quantity of work needed for or obtained from each subsystem but also on the quality of work. In a previous study Brandhorst showed that the power consumption for partially closed and completely closed regenerable life support systems was estimated as 3.5 kw/individual and 10-12 kw/individual respectively. With the increasing cost and scarcity of energy resources, our attention is drawn to evaluate the existing ECLSS technologies on the basis of their energy efficiency. In general the first law efficiency of a system is usually greater than 50 percent. From literature, the second law efficiency is usually about 10 percent. The estimation of second law efficiency of the system indicates the percentage of energy degraded as irreversibilities within the process. This estimate offers more room for improvement in the design of equipment. From another perspective, our objective is to keep the total entropy production of a life support system as low as possible and still ensure a positive entropy gradient between the system and the surroundings. The reason for doing so is as the entropy production of the system increases, the entropy gradient between the system and the surroundings decreases, and the system will gradually approach equilibrium with the surroundings until it reaches the point where the entropy gradient is zero. At this point no work can be extracted from the system. This is called the 'dead state' of the system.

Chatterjee, Sharmista↗
NASA NTRS1976

Ice in space: An experimental and theoretical investigation

The thermodynamics of water ice formation was experimentally investigated under a wide variety of conditions, including those of outer space. This information, and in particular, the lifetime of ice particles as a function of solar distance is an absolute requirement for proper interpretation of photometric profiles of comets. The sublimation of ice particles in a nonequilibrium situation was studied. An oscillating fiber microbalance was used to measure the sublimation rate of water droplets (which were suspended on a long quartz fiber which was oscillating in a vacuum chamber). The influence of particle size, surface temperature, and the index refraction from simulated solar radiation were studied in relation to ice formation. Also examined was the influence of impurities (clathrates) on ice formation. Windows in the vacuum chamber allowed the ice particles to be exposed to a 1 kilowatt xenon arc lamp which was used to simulate solar radiation. Ice is proposed as a possible energy source for comets, as amorphous water ice and ammonia in low temperature and pressure environments demonstrated a clear energy release upon warming. Motion pictures of ice formation were taken and photographs are shown.

Patashnick, H.↗
NASA NTRS1995

Macroscopic and bulk-controlled elastic modes in an interaction of interstitial alcali metal cations within a face-centered cubic crystalline fullerine

The calculations in a majority of previous works for the fulleride (AqC-60) crystals were performed within the framework of the rigid-lattice model, neglecting the distortion relaxation of the host fullerene (C-60) crystal caused by the interstitial alkali-metal (A) cations. However, an each cation is a source of a static distortion field, and the resulting field is a superposition of such fields generated by all cations. This is a reason why the host-crystal distortions depend on the A-cations configurations, i.e. on a type of a spatial bulk distribution of interstitial cations. The given paper seeks to find a functional relation between the amplitudes of the doping-induced structure-distortion waves and of static concentration ones. A semiphenomenological model is constructed here within the scope of statistical-thermodynamic treatment and using the lattice-statistics simulation method(*). In this model the effects due to the presence of q solute A cations over available interstices (per unit cell) on the static inherent reorientation and/or displacements of the solvent molecules from the 'average-lattice' sites' as well as on the lattice parameter a of a elastically-anysotropic 'cubic' C-60 crystal are taken into account.

Tatarenko, Valentine A.↗
NASA NTRS1999

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

Silicon carbide (SiC) and Si3N4 materials were tested in various turbine engine combustion environments chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high-speed aircraft. Representative chemical vapor-deposited (CVD), sintered, and composite materials were evaluated by furnace and high-pressure burner rig exposures. Although protective SiO2 scales formed in all cases, the evidence presented supports a model based on paralinear growth kinetics (i.e., parabolic growth moderated simultaneously by linear volatilization). The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were thus used to map SiO2 volatility (and SiC recession) over a range of temperatures, pressures, and velocities. The functional dependency of material recession (volatility) that emerged followed the form A[exp(-Q / RT)](P(sup x)v(sup y). These empirical relations were compared with rates predicted from the thermodynamics of volatile SiO and SiOxHy reaction products and a kinetic model of diffusion through a moving boundary layer. For typical combustion conditions, recession of 0.2 to 2 micrometers/hr is predicted at 1200 to 1400 C, far in excess of acceptable long-term limits.

Smialek, James L.↗
NASA NTRSDate not supplied

Open-source Numerical Modeling of Solidification Cracking Susceptibility: Application to Refractory Alloy Systems

Introduction. Alloys such as aluminum, nickel-base, and austenitic stainless steels are susceptible to solidification cracking during welding and 3D printing. Compositional optimization is one method used to effectively mitigate solidification cracking of those alloy systems. With the surge in hypersonic and in-space propulsion activities, refractory metals (Nb, Mo, Ta, W, and Re) and their alloy derivatives are increasing in importance due to their extreme high melting point and retention of high-temperature strength; however, their chemistry was most typically optimized to promote ductility during mechanical operations such as drawing and forming. Welding of such alloys has been a challenge due to a number of issues including solidification cracking, atmospheric contamination (O, C, and N), as well as a shift in ductile-to-brittle transition to higher temperature following grain growth induced by welding. Compositional optimization of refractory alloys for solidification cracking resistance in particular is desirable as their usage increases with the advent of advanced manufacturing methods such as 3D printing. This work evaluates the effect of compositional variation in refractory metal systems on the solidification cracking susceptibility with the goals of optimizing existing alloys and joining process techniques, and formulating new alloys with increased solidification cracking resistance. Experimental Procedures. A python code was developed in a Jupyter notebook environment (Michael and Sowards, 2023) to facilitate the calculation of crack susceptibility index proposed by Kou (2015). Composition is entered as a single point, or as a 1-D or 2-D array. The notebook calls pycalphad (Otis and Liu, 2017 and Bocklund et al, 2020) to calculate the evolution of fraction solid as a function of temperature (under either Scheil or equilibrium assumptions) and then evaluates steepness of the fraction solid curve near the terminal stage of solidification to predict solidification cracking resistance. Open source thermodynamic databases available at online repositories are used (van de Walle). The process is setup in an automated fashion to generate plots that show variation in solidification cracking susceptibility according to composition on 1-D line plots or 2-D contour plots. The Jupyter notebook and crack susceptibility algorithm was also integrated with a widely used commercial CALPHAD code for validation and alloy exploration. Results and Discussion. The crack susceptibility model was first validated against a series of refractory alloy compositions evaluated in past work which utilized a specialized Varestraint test built inside a vacuum chamber environment (Lessman and Gold, 1971). The alloys tested in the Varestraint apparatus included T-111 (Ta-8W-2Hf), ASTAR-811C (Ta-8W-1Re-0.7Hf-0.025C), FS-85 (Nb-27Ta-10W-1Zr), T-222 (Ta-9.6W-2.4Hf-0.01C), Ta-10W, B-66 (Nb-5Mo-5V-1Zr), and SCb-291 (Nb-10W-10Ta). The initial test of the model showed a strong correlation with empirical Varestraint data, i.e., a Spearman rank correlation between model predictions and hot cracking measurements was observed to be greater than 0.8. Following the validation, a set of refractory metal binary mixtures was investigated to evaluate sensitivity of Nb, Mo, W, and Ta to C, N, and O content. A series of plots were produced that suggest ppmw ranges of C, N, and O where solidification cracking increases significantly and reaches a maximum. Also comparative ranking of each primary refractory metal to each interstitial was produced. For example C produces greater cracking response in Mo whereas O produces greater cracking response in Ta and Nb. Such compositional values have utility in setting limits on pickup of these interstitial elements during welding and printing rather than using a one-size-fits-all approach. Furthermore, the results have use in determining additive powder recycling requirements, which is especially pertinent for refractory metal powders due to their high cost compared to conventional alloys. Another application created thousands of hypothetical alloys within the nominal specified composition range of two widely used refractory alloys C103 (Nb-10Hf-1Ti) and TZM (Mo-0.5Ti-0.1Zr). The cracking index was calculated for the alloys and results were fed into machine learning regression techniques including Multiple Linear Regression, Ridge Regression, and Lasso Regression to determine relative potency each alloying element had on computed solidification cracking index. A series of linear equations were produced that relate composition of C103 and TZM to solidification cracking index. The crack susceptibility of C103 for example is described by an equation of the form: cracking index ~ O + 0.667*C + 0.635*N + 0.00037*Ta – 0.0008*Hf (in wt.%) From that equation, it is clear that O has strong propensity to induce solidification cracking. Interestingly, Hf is shown to reduce calculated cracking response. Finally, realizing the potential of this method to discover new refractory alloy formulations across the period table that have low solidification cracking sensitivity, the code was applied to new untested alloy systems including W-Zr-C, W-Ta-C, and others. Conclusions. In summary, an open source numerical method has been developed using Python code to calculate Kou’s crack susceptibility index. The method was applied to refractory metals which are inherently difficult to study from a weldability testing standpoint since inert shielding gas is not sufficient and welding is typically done in vacuum, especially in light of findings presented here where oxygen has profound influence on solidification cracking. This work revealed the effect of compositional variations on a series of refractory metals and showed the framework defined here will be useful in 1) the development of new alloys that have improved weldability and 3D printability, 2) placing compositional limits on existing alloys, and 3) ensuring adequate controls of manufacturing processes such as 3D printing where powder reuse is critical. Keywords. pycalphad; Python; refractory metals; solidification cracking. References. B. Bocklund et. al. (2020) http://doi.org/10.5281/zenodo.3630657. S. Kou. (2015) https://doi.org/10.1016/j.actamat.2015.01.034. G.G. Lessmann and R.E. Gold. Welding Journal, issue 1, pp. 1-s – 8-s (1971). F.N. Michael and J.W. Sowards. NASA/TM-20230002218 (2023). R. Otis and Z.-K. Liu. (2017) http://doi.org/10.5334/jors.140. A. Van de Wallle et. al. (2018) https://doi.org/10.1016/j.calphad.2018.04.003.

pycalphad↗
NASA NTRS1983

The nonlinear viscoelastic response of resin matrix composites

The current paper describes the utilization of a thermodynamic based analytical nonlinear viscoelastic approach to represent lamina properties. Test data to verify the analysis for both transverse and shear properties of a T300/934 composite are presented. Master curves as a function of stress level and temperature are generated. Favorable comparisons between the traditional graphical and the current analytical approaches are shown.

Hiel, C. C.↗
NASA NTRS1975

Thermodynamics and phase separation of dense fully-ionized hydrogen-helium fluid mixtures

The free energy of a hydrogen-helium fluid mixture is evaluated for the temperatures and densities appropriate to the deep interior of a giant planet such as Jupiter. The electrons are assumed to be fully pressure-ionized and degenerate. In this regime, an appropriate first approximation to the ionic distribution functions can be found by assuming hard sphere interactions. Corrections to this approximation are incorporated by means of the perturbation theory of Anderson and Chandler. Approximations for the three-body interactions and the nonlinear response of the electron gas to the ions are included. It is predicted that a hydrogen-helium mixture, containing 10% by number of helium ions, separates into hydrogen-rich and helium-rich phases below about 8000 K, at the pressures relevant to Jupiter (4-40 Megabars). It is also predicted that the alloy occupies less volume per ion than the separated phases. The equations of state and other thermodynamic derivatives are tabulated. Implications of these results are discussed.

Stevenson, D. J.↗