Curve-fitting techniques and applications to thermodynamics.
Mathematical methods for incorporating thermodynamic data into least squares fit
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.
Mathematical methods for incorporating thermodynamic data into least squares fit
Unlike the combustion of homogeneous gas mixtures, there are practically no reliable fundamental data (i.e., laminar burning velocity, flammability limits, quenching distance, minimum ignition energy) for the combustion of heterogeneous dust suspensions. Even the equilibrium thermodynamic data such as the constant pressure volume combustion pressure and the constant pressure adiabatic flame temperature are not accurately known for dust mixtures. This is mainly due to the problem of gravity sedimentation. In normal gravity, turbulence, convective flow, electric and acoustic fields are required to maintain a dust in suspension. These external influences have a dominating effect on the combustion processes. Microgravity offers a unique environment where a quiescent dust cloud can in principle be maintained for a sufficiently long duration for almost all combustion experiments (dust suspensions are inherently unstable due to Brownian motion and particle aggregation). Thus, the microgravity duration provided by drop towers, parabolic flights, and the space shuttle, can all be exploited for different kinds of dust combustion experiments. The present paper describes some recent studies on microgravity combustion of dust suspension carried out on the KC-135 and the Caravelle aircraft. The results reported are obtained from three parabolic flight campaigns.
Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. One of the main reasons for this problem are the possible unrealistic assumptions that models/reanalyses make about the subgrid-scale (SGS) variability of meteorological properties, as well as the relationship between SGS variability and grid-scale (GS) cloud properties. We utilize cloud and thermodynamic data of low level (primarily) liquid clouds collected from two aircraft campaigns conducted over the Beaufort Sea to better understand and characterize this problem. Examining data from the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign and the 1998 First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)–Arctic Cloud Experiment (ACE) reveals that GS cloud water variability is closely related with SGS distribution of total water (i.e. water vapor + cloud water). We investigate two related approaches to prediction of GS cloud properties from SGS variability: the critical saturation ratio method, and the critical relative humidity method. We find significant correlation between GS cloud water and SGS supersaturation when the critical saturation ratio is set at 100%, as well as a notable relationship between GS cloud water and the width of the SGS total water distribution. Critical relative humidity also compares well with GS cloud water. However, we also find that the assumptions of a static critical saturation ratio of 100% to be unrealistic, as well as a fixed SGS distribution width. Empirical calculations from the ARISE data show a large sensitivity of these SGS variables to GS relative humidity, and so a SGS parameterization allowing them to vary according to GS thermodynamic properties may result in more realistic GS cloud water values.
Silicates are a common class of materials that are often exposed to high temperatures. The behavior of these materials needs to be understood for applications as high temperature coatings in material science as well as the constituents of lava for geological considerations. The vaporization behavior of these materials is an important aspect of their high temperature behavior and it also provides fundamental thermodynamic data. The application of Knudsen effusion mass spectrometry (KEMS) to silicates is discussed. There are several special considerations for silicates. The first is selection of an appropriate cell material, which is either nearly inert or has well-understood interactions with the silicate. The second consideration is proper measurement of the low vapor pressures. This can be circumvented by using a reducing agent to boost the vapor pressure without changing the solid composition or by working at very high temperatures. The third consideration deals with kinetic barriers to vaporization. The measurement of these barriers, as encompassed in a vaporization coefficient, is discussed. Current measured data of rare earth silicates for high temperature coating applications are discussed. In addition, data on magnesium-iron-silicates (olivine) are presented and discussed.
Thermodynamic properties for a Jet-A fuel were determined by Shell Development Company in 1970 under a contract for NASA Lewis Research Center. The polynomial fit necessary to include Jet-A fuel (liquid and gaseous phases) in the library of thermodynamic properties of the NASA Lewis Chemical Equilibrium Program is calculated. To verify the thermodynamic data, the temperatures of mixtures of liquid Jet-A injected into a hot nitrogen stream were experimentally measured and compared to those calculated by the program. Iso-octane, a fuel for which the thermodynamic properties are well known, was used as a standard to calibrate the apparatus. The measured temperatures for the iso-octane/nitrogen mixtures reproduced the calculated temperatures except for a small loss due to the non-adiabatic behavior of the apparatus. The measurements for Jet-A were corrected for this heat loss and showed excellent agreement with the calculated temperatures. These experiments show that this process can be adequately described by the thermodynamic properties fitted for the Chemical Equilibrium Program.
Rare earth element exploration requires a quantitative understanding of factors governing their mobilization and economic concentration. However, the behavior of rare earth elements in carbonate- bearing hydrothermal fluids associated with carbonatite-hosted deposits is poorly understood, and conflicting mechanisms of rare earth transport by anionic ligands and alkali behavior have been described. Here, we report quantitative data to characterize the role of carbonate-bearing solutions in the hydrothermal mobilization of neodymium. Solubility studies of neodymium phosphate were performed at temperatures ranging from 100 to 600 °C in carbonate-bearing solutions. The thermodynamic data determined for the predominant complex were used to model the separation of neodymium from thorium in a simple flow-through system based on fluid and mineral compositions characteristic of carbonatite deposits. Our data suggest that neodymium transport is controlled by the stability of the carbonate species NdCO 3 OH o , and at temperatures of 500–600 °C, the concentrations of neodymium in solutions can reach ~1000 ppm.
Analysis of phase diagrams and cation distributions within crystalline solutions as means of obtaining thermodynamic data on rock forming crystalline solutions is discussed along with some aspects of partitioning of elements in coexisting phases. Crystalline solutions, components in a silicate mineral, and chemical potentials of these components were defined. Examples were given for calculating thermodynamic mixing functions in the CaW04-SrW04, olivine-chloride solution, and orthopyroxene systems.
During September-October 2001, the East Pacific Investigation of Climate Processes in the Coupled Ocean-Atmosphere System (EPIC-2001) ITCZ field campaign focused on studies of deep convection in the warm-pool region of the East Pacific. In addition to the TAO mooring array, observational platforms deployed during the field phase included the NOAA ship RN Ronald H. Brown, the NSF ship RN Horizon, and the NOAA P-3 and NCAR C-130 aircraft. This study combines C-band Doppler radar, rawinsonde, and surface heat flux data collected aboard the RN Brown to describe ITCZ convective structure and rainfall statistics in the eastern Pacific as a function of 3-5 day easterly wave phase. Three distinct easterly wave passages occurred during EPIC-2001. Wind and thermodynamic data reveal that the wave trough axes exhibited positively correlated U and V winds and a slight westward phase tilt with height. A relatively strong (weak) northeasterly deep tropospheric shear followed the trough (ridge) axis. Temperature and humidity perturbations exhibited mid-to upper level cooling (warming) and drying (moistening) in the northerly (trough and southerly) phase. At low levels warming (cooling) occurred in the northerly (southerly) phase with little change in the relative humidity, though mixed layer mixing ratios were larger during the northerly phase. When composited, radar, sounding, lightning and surface heat flux observations suggest the following systematic behavior as a function of wave phase: approximately zero to one quarter wavelength ahead of (behind) the wave trough in northerly (southerly) flow, larger (smaller) CAPE, lower (higher) CIN, weaker (stronger) tropospheric shear, higher (lower) conditional mean rain rates, higher (lower) lightning flash densities, and more (less) robust convective vertical structure occurred. Latent and sensible heat fluxes reached a minimum in the northerly phase and then increased through the trough, reaching a peak during the ridge phase (leading the peak in CAPE). From a radar echo coverage perspective, larger areas of light rain and slightly larger (10%) area averaged rain rates occurred in the vicinity of, and just behind, the trough axes in southerly flow. Importantly, the transition in convective structure observed across the trough axis when considered with the relatively small change in area mean rain rates suggests the presence of a transition in the vertical structure of diabatic heating across the easterly waves examined. The inferred transition in heating structure is supported by radar diagnosed divergence profiles that exhibit convective (stratiform) characteristics ahead of (behind) the trough.
Thermodynamic data from interferometric density profile studies and light-scattering experiments near the critical isochore of Xe, CO2 and SF6 provide a basis for examining the hypothesized two-scale-factor universality for the correlation function of fluids near the gas-liquid critical point. For the investigation, three-scale-factor universality is assumed, with Ising-like critical exponent values obtained through the renormalization group technique. The two thermodynamic scale factors are found from the density profiles, while the scale factor for the correlation length is obtained from the light-scattering data.
Sixty-nine soundings were taken at San Nicolas Island during the FIRE marine stratocumulus experiment. The archived thermodynamic data at 5mb vertical resolution were reduced to 25mb vertical resolution by averaging. Ceilometer data showed the time variation of cloud-base. A subset of 20 soundings taken during a 2 day period of light winds were analyzed. The soundings are grouped into 6 hr. averages (with 5 soundings in each category) to examine the diurnal variation of the thermodynamic structure. Consistent with the solar heating during the daytime, cloud-base rises during the daytime hours and the cloud layer thins. The subcloud layer also becomes more stable, suggesting some uncoupling of the cloud layer from the surface.
Current studies at NASA Glenn on oxide thermodynamics are discussed. Previous studies on the vaporization of B2O3 in reducing atmospheres led to inconsistent studies when B was used as a reductant. It is shown that liquid B2O3 does not wet B and a clear phase separation was noted in the Knudsen cell. This problem was solved by using FeB and Fe2B to supply a different and constant activity of B. The thermodynamic data thus derived are compared to quantum chemical composite calculations. A major problem in high temperature mass spectrometry is the determination of accurate ionization cross sections, particularly for molecules. The method of Deutsch and Mark shows promise and some sample calculations are discussed. Finally current studies on the thermodynamics of rare earth silicates are discussed. Here the problems are obtaining a measurable signal from SiO2 vaporization and non-equilibrium vaporization. The use of a Ta reducing agent provides a stronger signal, which is related to silica activity. The Whitman-Motzfeld relation adapted to KEMS measurements is applied to obtain equilibrium pressures.
Silica and sapphire are very stable oxide materials. However, at high temperatures they are both susceptible to chemical degradation by three major mechanisms: (1) Deposition of fluxing oxides (2) Non-reactive vaporization (3) Reactive vaporization, particularly by water vapor. In order to model and predict each of these, accurate thermodynamic data and information on the environment are needed. Experimental and theoretical methods used at NASA for identifying the generated gas species and obtaining the relevant thermochemical data are discussed. Modeling and prediction of degradation are discussed for different environments and flow conditions. In addition, surface devitrification or etch pitting is likely under some conditions.
This report documents the library of thermodynamic data used with the NASA Glenn computer program CEA (Chemical Equilibrium with Applications). This library, containing data for over 2000 solid, liquid, and gaseous chemical species for temperatures ranging from 200 to 20,000 K, is available for use with other computer codes as well. The data are expressed as least-squares coefficients to a seven-term functional form for C((sup o)(sub p)) (T) / R with integration constants for H (sup o) (T) / RT and S(sup o) (T) / R. The NASA Glenn computer program PAC (Properties and Coefficients) was used to calculate thermodynamic functions and to generate the least-squares coefficients. PAC input was taken from a variety of sources. A complete listing of the database is given along with a summary of thermodynamic properties at 0 and 298.15 K.
Actinide-bearing molten salts for use as fuels are an essential part of next generation molten salt reactors. Yet, numerous multicomponent salt mixtures are underdeveloped or have not been investigated. Here, this study, based on a combination of experimental and modeling techniques, is dedicated to determining and understanding a variety of properties of the ternary system of NaCl–UCl 3 –PuCl 3 , which represents a scenario for burnup of NaCl–UCl 3 fuel, at two compositions (∼10 and 5 mol % PuCl 3 in eutectic NaCl–UCl 3 pseudobinary) and a range of temperatures. Evaluation of the heat flow and mass loss data showed the 0.61NaCl–0.30UCl 3 –0.09PuCl 3 salt had a melting temperature of 551 ± 5 °C. Two additional thermal effects were observed occurring at approximately 410 and 494 °C. The transition occurring at 410 °C may be due to the presence of oxide in the salt. Extrapolation of thermodynamic data indicates the transition occurring at 494 °C is due to the formation of a liquid phase. Experimental testing determined the density of this system is a linear function of temperature and can be represented by the equation ρ = 4.014–0.0010T(°C), R 2 = 0.992. Additionally, by using atomistic modeling, we found that increasing the PuCl 3 content from 5 to 10 mol % led to the formation of larger Pu 3+ clusters and slower transport of ions.
This paper presents some recent results germane to airborne windshear system design and certification. We first discuss the data analysis procedure and the associated caveats. The relative frequency, severity, and duration of microburst hazards at the various locations is important for determining the tradeoffs between safety and operational impact of false alerts which are encompassed in detection system thresholds. We then consider radar/lidar design issues such as reflective in microbursts and the vertical structure of outflows. Finally, we provide recent surface thermodynamic data associated with microbursts.
Earth-based, earth-orbital, and spacecraft observational data are used in the present evaluation of Venus atmosphere-surface interactions to quantitatively characterize the reactions between C, H, S, Cl, F, and N gases and plausible surface minerals. Calculation results are used to predict stable minerals and mineral assemblages on the Venus surface, in order to ascertain which (if any) of the atmospheric gases are buffeted by mineral assemblages. Chemical equilibrium calculations using extant thermodynamic data on scapolite minerals predict that carbonate-bearing scapolite and sulfate meionite are unstable on the surface of Venus, while chloride-bearing scapolite is stable.
As advances are being made regarding the performance of nuclear fuels, uranium carbides, and composites, such as (U,Zr)C and UO2 + UCx, have recently gained significant interest for deployment in nuclear space propulsion and high temperature gas-cooled reactors, respectively. However, the phase equilibria of several fission products in carbide systems remain unknown and may impact the overall fuel performance, specifically for particle nuclear fuels that are designed for commercial nuclear energy. Furthermore, comprehensive thermodynamic data on Rare Earth (RE) carbides, such as the Nd-C and Ce-C binary systems, remain limited. Presented in this study are the synthesis methods and characterizations of several Nd-C and Ce-C compositions. The findings from this research provide insights on the stability of RE-C binaries that form in irradiated nuclear fuels and address a critical knowledge gap in the current state of thermodynamics for two key RE-C systems.
Chemical equilibrium models used currently to interpret observations of Jupiter are reexamined using new data defining thermal profiles, which are substantially different from those used in the previous models. A model is developed for the chemical reactions controlling the composition of the upper troposphere on Jupiter, specifically the cloud-forming region from 10 bar to 0.1 bar, which includes, for the first time, the effects of aqueous chemistry on the composition and the vertical distribution of many measurable species in the atmosphere, identifying the factors influencing their abundances above the H2O cloud. The thermodynamic data for potential condensates on Jupiter, i.e., NH3(s), NH4SH(s), (NH4)2S(s), and H2S(s), are reexamined, recognizing the lack of data on sulfides for the temperature range of interest on Jupiter. Vertical profiles of mixing ratios for CO2, H2S, NH3, and H2, obtained for several assumed bulk abundances with respect to solar, are presented.