Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic consistency”

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 127 records · Page 7

Surface film effects on drop tube undercooling studies

The effects of various gaseous atmospheric constituents on drop-tube solidified samples of elemental metals were examined from a microstructural standpoint. All specimens were prepared from the purest available elements, so effects of impurities should not account for the observed effects. The drop-tube gas has a definite effect on the sample microstructure. Most dramatically, the sample cooling rate is effected. Some samples receive sufficient cooling to solidify in free fall while others do not, splating at the end of the drop tube in the sample catcher. Gases are selectively absorbed into the sample. Upon solidification gas can become less soluble and as a result forms voids within the sample. The general oxidation/reduction characteristics of the gas also affect sample microstructures. In general, under the more favorable experimental conditions including reducing atmospheric conditions and superheatings, examination of sample microstructures indicates that nucleation has been suppressed. This is indicated by underlying uniform dendrite spacings throughout the sample and with a single dendrite orientation through most of the sample. The samples were annealed yielding a few large grains and single or bi-crystal samples were commonly formed. This was especially true of samples that were inadvertently greatly superheated. This is in contrast with results from a previous study in which surface oxides were stable and contained numerous sites of nucleation. The number of nucleation events depends upon the surface state of the specimen as determined by the atmosphere and is consistent with theoretical expectations based upon the thermodynamic stability of surface oxide films. Oxide-free specimens are characterized by shiny surfaces, with no observable features under the scanning electron microscope at 5000X.

Ethridge, E. C.↗

Enthalpies of a binary alloy during solidification

The purpose of the paper is to present a method of calculating the enthalpy of a dendritic alloy during solidification. The enthalpies of the dendritic solid and interdendritic liquid of alloys of the Pb-Sn system are evaluated, but the method could be applied to other binaries, as well. The enthalpies are consistent with a recent evaluation of the thermodynamics of Pb-Sn alloys and with the redistribution of solute in the same during dendritic solidification. Because of the heat of mixing in Pb-Sn alloys, the interdendritic liquid of hypoeutectic alloys (Pb-rich) of less than 50 wt pct Sn has enthalpies that increase as temperature decreases during solidification.

Poirier, D. R.↗

Vapor-liquid equilibrium thermodynamics of N2 + CH4 - Model and Titan applications

A thermodynamic model is presented for vapor-liquid equilibrium in the N2 + CH4 system, which is implicated in calculations of the Titan tropospheric clouds' vapor-liquid equilibrium thermodynamics. This model imposes constraints on the consistency of experimental equilibrium data, and embodies temperature effects by encompassing enthalpy data; it readily calculates the saturation criteria, condensate composition, and latent heat for a given pressure-temperature profile of the Titan atmosphere. The N2 content of condensate is about half of that computed from Raoult's law, and about 30 percent greater than that computed from Henry's law.

Thompson, W. R.↗

Polarimetric signatures of sea ice. 1: Theoretical model

Physical, structral, and electromagnetic properties and interrelating processes in sea ice are used to develop a composite model for polarimetric backscattering signatures of sea ice. Physical properties of sea ice constituents such as ice, brine, air, and salt are presented in terms of their effects on electromagnetic wave interactions. Sea ice structure and geometry of scatterers are related to wave propagation, attenuation, and scattering. Temperature and salinity, which are determining factors for the thermodynamic phase distribution in sea ice, are consistently used to derive both effective permittivities and polarimetric scattering coefficients. Polarmetric signatures of sea ice depend on crystal sizes and brine volumes, which are affected by ice growth rates. Desalination by brine expulsion, drainage, or other mechanisms modifies wave penetration and scattering. Sea ice signatures are further complicated by surface conditions such as rough interfaces, hummocks, snow cover, brine skim, or slush layer. Based on the same set of geophysical parameters characterizing sea ice, a composite model is developed to calculate effective permittivities and backscattering covariance matrices at microwave frequencies to interpretation of sea ice polarimetric signatures.

Nghiem, S. V.↗

COLD-SAT Dynamic Model Computer Code

COLD-SAT Dynamic Model (CSDM) computer code implements six-degree-of-freedom, rigid-body mathematical model for simulation of spacecraft in orbit around Earth. Investigates flow dynamics and thermodynamics of subcritical cryogenic fluids in microgravity. Consists of three parts: translation model, rotation model, and slosh model. Written in FORTRAN 77.

Bollenbacher, G.↗

Polarimetric Signatures of Sea Ice: Theoretical Model - Part 1

Physical, structural, and electromagnetic properties and interrelating processes in sea ice are used to develop a composite model for polarimetric backscattering signatures of sea ice. Physical properties of sea ice constituents such as ice, brine, air, and salt are presented in terms of their effects on electromagnetic wave interactions. Sea ice structure and geometry of scatterers are related to wave propagation, attenuation, and scattering. Temperature and salinity, which are determining factors for the thermodynamic phase distribution in sea ice, are consistently used to derive both effective permittivities and polarimetric scattering coefficients. Polarimetric signatures of sea ice depend on crystal sizes and brine volumes, which are affected by ice growth rates. Desalination by brine expulsion, drainage, or other mechanisms modifies wave penetration and scattering. Sea ice signatures are further complicated by surface conditions such as rough interfaces, hummocks, snow cover, brine skim, or slush layer. Based on the same set of geophysical parameters characterizing sea ice, a composite model is developed to calculate effective permittivities and backscattering covariance matrices at microwave frequencies for interpretation of sea ice polarimetric signatures.

Nghiem, S. V.↗

Plasticity and Damage Modeling of Stress Asymmetry and Dynamic Behavior of AFS Additive Manufactured Aluminum Alloy 2219

The Solid State Additive Manufacturing (AM) process referred as MELD that fabricated the samples in this study, provides a new path for repairing, coating, joining and additive manufacturing metals and metal matrix composites. This research will be the first application of a physics-based microstructure dependent internal state variable (ISV) plasticity and damage material model to capture the mechanical response of an AM Aluminum Alloy (AA) 2219 via the MELD process. In this research, a microstructure-based internal state variable (ISV) plasticity-damage model was used to capture the mechanical behavior of AFS 2219 aluminum alloy. Aeroprobe Corporation, creator and patent holder for the MELD process, fabricated the material by pushing a solid filler rod of AA2219-T861 material through a hollow rotating tool onto an AA2219 T851 plate substrate. As feedstock, solid or powder precursor metals are pushed through a nonconsumable rotating cylindrical tool. Herein, added layers are deposited and metallurgically bonded to substrate material or previously deposited layers by the heat generated from the rotating tool through plastic deformation of the filler material. Once a layer has been added, the tool height increases, and starts the deposition of the next layer. This process results in beneficial properties such as grain refinement, homogenization and reduced porosity (fully dense). This process will experience temperatures similar to those in the weld nugget zone (WNZ) in friction stir welding (FSW), ranging from 0.6-0.9 Tm, with Tm being the melting point of the material. MELD is highly scalable with AA deposition rates reaching over 1000 cm3/hr, which allows for MELD being used for repairs, coatings, and building components. A motivating factor driving the research for physics-based history dependent material modeling of MELD components is the ability to accurately capture the stress-state and strain rate dependence in the material caused by variations in material microstructure from the MELD processing of new or repaired components. The ISV model incorporates microstructural content and is consistent with continuum level kinematics, kinetics, and thermodynamics. These features allow the ISV model to capture large deformations at the structural scale using the kinematic and isotropic hardening, while microscale damage is obtained from the microstructural features. The benefits of the ISV model arise from the inclusion of structure-property relationships identified from microstructural characterization and experimentation. The Bauschinger effect (BE) is an important concept, vital in the accurate prediction of cyclic stress-strain response of ductile materials such as metals. The ISV model has been successfully used to capture the behavior and damage, and the BE of different aluminum alloys and steels. The ISV model uses kinematic and isotropic hardening to help capture deformations of the material at the macro scale. To understand this hardening relationship, calculating the kinematic and isotropic hardening relationship in the material is warranted for a high-fidelity model. Electron Backscattered Diffraction (EBSD) was used to characterize the as-fabricated microstructure, where a fully-dense equiaxed grain morphology with average grain size of 2.5 μm was observed. Microhardness mapping of the as-built structures, monotonic tension and compression experiments at both quasi-static (0.001/s) strain rates, tension-followed-by-compression and compression-followed-by-tension experiments were performed to obtain the set of plasticity and damage constants necessary to capture strain rate and stress state behavior of this additive material. To calibrate the plasticity-damage model, a single set of constants were determined to capture the different stress states the MELD AA2219. One set of the constants was determined from experimental true stress-strain curves for the tension and compression data. Additionally, microstructural information and data from the open literature were used as the other model constants. This research is a first of its kind for AFS AA2219, includes correlating the ISV model to the monotonic experimental results that capture the isotropic and kinematic plasticity mechanical response.

Rivera, O. G.↗

Consistent Initial Conditions for the DNS of Compressible Turbulence

Relationships between diverse thermodynamic quantities appropriate to weakly compressible turbulence are derived. It is shown that for turbulence of a finite turbulent Mach number there is a finite element of compressibility. A methodology for generating initial conditions for the fluctuating pressure, density and dilatational velocity is given which is consistent with finite Mach number effects. Use of these initial conditions gives rise to a smooth development of the flow, in contrast to cases in which these fields are specified arbitrarily or set to zero. Comparisons of the effect of different types of initial conditions are made using direct numerical simulation of decaying isotropic turbulence.

Ristorcelli, J. R.↗

Properties of small Ar sub N-1 K/+/ ionic clusters

A self-consistent formalism is developed that, based upon a many-body potential, dynamically determines the thermodynamic properties of ionic clusters without an a priori designation of the equilibrium structures. Aggregates consisting of a single closed shell K(+) ion and N-1 isoelectronic argon atoms were studied. The clusters form crystallites at low temperatures, and melting transitions and spontaneous dissociations are indicated. The results confirm experimental evidence that shows that ionic clusters become less stable with increasing N. The crystallite structures formed by four different clusters are isosceles triangle, skewed form, octahedron with ion in the middle, and icosahedron with the ion in the middle.

Etters, R. D.↗

Equilibrium thermodynamics and stochastic nonlinear acoustic fields

A crystalline solid is considered to consist of a large number of incoherent nonlinear acoustic radiation sources identified with the vibrating particles of the crystalline lattice. Randomization of the field, together with the assumption of a stochastically independent, fluctuating, radiation field at the absolue zero of temperature, leads to an expression of the temperature-dependent radiation field in terms of the zero-point field. The equation is identified with the Planck distribution formula of quantum mechanics in the linear field limit. The thermodynamic state functions are also obtained in terms of the nonlinear acoustic modal energies per unit mass and reduce to the results of the Debye-Einstein stochastic quantum oscillator model in the linear field limit.

Cantrell, J. H.↗

Calculations and curve fits of thermodynamic and transport properties for equilibrium air to 30000 K

A self-consistent set of equilibrium air values were computed for enthalpy, total specific heat at constant pressure, compressibility factor, viscosity, total thermal conductivity, and total Prandtl number from 500 to 30,000 K over a range of 10(exp -4) atm to 10(exp 2) atm. The mixture values are calculated from the transport and thermodynamic properties of the individual species provided in a recent study by the authors. The concentrations of the individual species, required in the mixture relations, are obtained from a free energy minimization calculation procedure. Present calculations are based on an 11-species air model. For pressures less than 10(exp -2) atm and temperatures of about 15,000 K and greater, the concentrations of N(++) and O(++) become important, and consequently, they are included in the calculations determining the various properties. The computed properties are curve fitted as a function of temperature at a constant value of pressure. These curve fits reproduce the computed values within 5 percent for the entire temperature range considered here at specific pressures and provide an efficient means for computing the flowfield properties of equilibrium air, provided the elemental composition remains constant at 0.24 for oxygen and 0.76 for nitrogen by mass.

Gupta, Roop N.↗

Design of closed-cycle MHD generator with nonequilibrium ionization and system

A method is developed to include the nonequilibrium ionization process in the MHD generator duct design equations, and these equations are coupled to the thermodynamic conditions of the closed cycle system. This is used to relate MHD generator size, configuration and gas conditions to the overall thermodynamic efficiency of the system. The system studied consists of an MHD loop (Ar + Cs or He + Cs) topping a steam bottoming plant.

Voshall, R. E.↗

Interior models of Jupiter

Understanding the interior of Jupiter depends upon our knowledge of the thermodynamics and transport properties of hydrogen-helium mixtures at high pressures and temperatures. The current status of this knowledge is reviewed, and attention is given to the metallic-molecular hydrogen transition and the limited solubility of helium in hydrogen. Models of Jupiter are constructed which are consistent with all the observations to date, but which make various assumptions about the thermodynamics and composition of the interior. These models typically consist of a rocky core surrounded by a nearly-solar fluid mixture. In contrast to the models of Podolak and Cameron, a large enhancement of water or helium is not found to be essential. It is concluded that further progress in constructing interior models requires a better understanding of the thermodynamics of dense molecular hydrogen.

Stevenson, D. J.↗

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

Physical processes within the nocturnal stratus-topped boundary layer

There are many physical processes involved in the stratus-topped boundary layer: longwave radiation cooling, entrainment, latent heating, surface heating, solar heating, and drizzling, for example. The manner in which the processes combine to maintain the turbulence within the stratus-topped boundary layer remains an unsolved problem. The large eddy simulation technique is used to examine the first four physical processes mentioned above. First, the contribution of each physical process to the thermodynamic differences between the updraft and downdraft branches of turbulent circulations is examined through a conditional sampling. Second, these mean thermodynamic differences are shown to express well the vertical distributions of heat and moisture fluxes within stratus-topped boundary layers. These provide a method to validate the process partitioning technique. (This technique assumes that the net flux profile can be partitioned into different component-flux profiles according to physical processes, and that each partitioned component flux is linear in height.) In this paper, the heat and moisture fluxes are process partitioned, and each component-flux is found to contribute to the net flux in a way that is consistent with its corresponding process contribution to the mean thermodynamic differences between updrafts and downdrafts. Also, the net flux obtained by summing all component-fluxes agrees very well with that obtained directly from the large-eddy simulations.

Moeng, Chin-Hoh↗

A study of the applicability of nucleation theory to quasi-thermodynamic transitions of second and higher Ehrenfest-order, supplement 3

The work includes an investigation of the applicability of the nucleation theory to second and higher order thermodynamic transitions in the Ehrenfest sense, and a number of significant conclusions relevant to first order transitions, as well. The underlying theoretical method consisted of expanding the Gibbs' free energy in a Maclarin or Taylor series and then using fundamental thermodynamic determinable quantities, and interpreting the results. Work was performed on the existence and interpretation of an interfacial energy between phases in a second order transition in addition to an investigation of the solid-liquid interfacial energy for various polymers. Extensive considerations were devoted to various aspects of a particular polymer, polyvinylidene fluoride (PVDF or PVF2), including an experimetal investigation of the effects of an applied electric field on the morphology of melt crystallization and on the nucleation and growth of polarized domains.

Barker, R. E., Jr.↗

The coupling of radiative transfer to quasi 1-D flows with thermochemical nonequilibrium

Quasi-one-dimensional nonequilibrium nozzle flows with coupled radiative transfer are considered. The strongly coupled formulation of radiation and flowfield leads to a governing set of integro-differential equations. A fully implicit numerical method using the full matrix inversion or block iteration methods is presented to solve these equations. The nonequilibrium gas model consists of two chemical species, molecular and atomic nitrogen. The thermodynamic state of the gas is described by two temperatures, translational-rotational and vibrational, and the thermal radiation is assumed to be governed by the vibrational temperature. In radiative transfer, gases are assumed to be absorbing and emitting, and a detailed spectral dependency of the absorption coefficient is prescribed for a non-gray gas. The numerical solutions of strongly radiating nonequilibrium flows are presented for both gray and non-gray gases.

Gokcen, Tahir↗

An experimental study of trace element partitioning between olivine, orthopyroxene and melt in chondrules - Equilibrium values and kinetic effects

Mineral/melt partition coefficients were measured using an ion microprobe for 32 elements in orthopyroxene and olivine in equilibrium and dynamic crystallization experiments on compositions corresponding to chondrules. The mineral/melt partition coefficients calculated from the measured concentrations for both olivine and orthopyroxene show very little change between equilibrium experiments and dynamic experiments with cooling rates of up to 100 C/h. The results provide a self-consistent set of partition coefficients that can be used in thermodynamic models of equilibrium and kinetic partitioning between olivine, orthopyroxene, and melt. These data can be used in models of partial melting and crystal fractionation in olivine- and orthopyroxene-rich systems, such as chondrules. The results may also be applicable to mantle peridotites, komatiitic and picritic lavas, and ultramafic intrusions.

Kennedy, A. K.↗