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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Viscoplastic Model Development to Account for Strength Differential: Application to Aged Inconel 718 at Elevated Temperature

The magnitude of yield and flow stresses in aged Inconel 718 are observed to be different in tension and compression. This phenomenon, called the Strength differential (SD), contradicts the metal plasticity axiom that the second deviatoric stress invariant alone is sufficient for representing yield and flow. Apparently, at least one of the other two stress invariants is also significant. A unified viscoplastic model was developed that is able to account for the SD effect in aged Inconel 718. Building this model involved both theory and experiments. First, a general threshold function was proposed that depends on all three stress invariants and then the flow and evolution laws were developed using a potential-based thermodynamic framework. Judiciously chosen shear and axial tests were conducted to characterize the material. Shear tests involved monotonic loading, relaxation, and creep tests with different loading rates and load levels. The axial tests were tension and compression tests that resulted in sufficiently large inelastic strains. All tests were performed at 650 C. The viscoplastic material parameters were determined by optimizing the fit to the shear tests, during which the first and the third stress invariants remained zero. The threshold surface parameters were then fit to the tension and compression test data. An experimental procedure was established to quantify the effect of each stress invariant on inelastic deformation. This requires conducting tests with nonproportional three-dimensional load paths. Validation of the model was done using biaxial tests on tubular specimens of aged Inconel 718 using proportional and nonproportional axial-torsion loading. These biaxial tests also helped to determine the most appropriate form of the threshold function; that is, how to combine the stress invariants. Of the set of trial threshold functions, the ones that incorporated the third stress invariant give the best predictions. However, inclusion of the first stress invariant does not significantly improve the model predictions. The model shows excellent predictive capability for nonproportional load paths. Additionally, it reduces to the well-known models of Mises Drucker and Drucker-Prager. The requisite experiments involve reasonably simple load paths in the axial-shear stress plane and hence can be performed on a variety of different materials: be they metallic, geological. polymeric, ceramic or granular. The general form of the threshold function allows representation of inelastic deformation in a range of materials.

Iyer, Saiganesh↗

NASA NTRS1974

Nonequilibrium boundary layer at a stagnation point for a hydrogen-helium stream over ablating graphite

The nonequilibrium axisymmetric stagnation point boundary layer over an ablating graphite surface is considered. The external stream is a high temperature mixture of hydrogen and helium. Variable thermodynamic and transport properties are assumed. Lennard-Jones potential model is used to calculate the transport coefficients of each species. Although the mixture rules for viscosity of the gas mixture are used, the weighting functions are more sophisticated than those commonly employed. For the conductivity of the mixture, generalized Wassiljewa coefficients are used. Seven species with 28 dissociation/recombination reactions are considered. Hansen's model for the dissociation rate constants is employed. The recombination rate constants are obtained by invoking detailed balance principles assisted by the JANAF thermodynamic data and the Hansen-Pearson thermodynamic data for C3.

Liu, T.-M.↗

NASA NTRS2004

Application of Simulated Annealing and Related Algorithms to TWTA Design

Simulated Annealing (SA) is a stochastic optimization algorithm used to search for global minima in complex design surfaces where exhaustive searches are not computationally feasible. The algorithm is derived by simulating the annealing process, whereby a solid is heated to a liquid state and then cooled slowly to reach thermodynamic equilibrium at each temperature. The idea is that atoms in the solid continually bond and re-bond at various quantum energy levels, and with sufficient cooling time they will rearrange at the minimum energy state to form a perfect crystal. The distribution of energy levels is given by the Boltzmann distribution: as temperature drops, the probability of the presence of high-energy bonds decreases. In searching for an optimal design, local minima and discontinuities are often present in a design surface. SA presents a distinct advantage over other optimization algorithms in its ability to escape from these local minima. Just as high-energy atomic configurations are visited in the actual annealing process in order to eventually reach the minimum energy state, in SA highly non-optimal configurations are visited in order to find otherwise inaccessible global minima. The SA algorithm produces a Markov chain of points in the design space at each temperature, with a monotonically decreasing temperature. A random point is started upon, and the objective function is evaluated at that point. A stochastic perturbation is then made to the parameters of the point to arrive at a proposed new point in the design space, at which the objection function is evaluated as well. If the change in objective function values (Delta)E is negative, the proposed new point is accepted. If (Delta)E is positive, the proposed new point is accepted according to the Metropolis criterion: rho((Delta)f) = exp((-Delta)E/T), where T is the temperature for the current Markov chain. The process then repeats for the remainder of the Markov chain, after which the temperature is decremented and the process repeats. Eventually (and hopefully), a near-globally optimal solution is attained as T approaches zero. Several exciting variants of SA have recently emerged, including Discrete-State Simulated Annealing (DSSA) and Simulated Tempering (ST). The DSSA algorithm takes the thermodynamic analogy one step further by categorizing objective function evaluations into discrete states. In doing so, many of the case-specific problems associated with fine-tuning the SA algorithm can be avoided; for example, theoretical approximations for the initial and final temperature can be derived independently of the case. In this manner, DSSA provides a scheme that is more robust with respect to widely differing design surfaces. ST differs from SA in that the temperature T becomes an additional random variable in the optimization. The system is also kept in equilibrium as the temperature changes, as opposed to the system being driven out of equilibrium as temperature changes in SA. ST is designed to overcome obstacles in design surfaces where numerous local minima are separated by high barriers. These algorithms are incorporated into the optimal design of the traveling-wave tube amplifier (TWTA). The area under scrutiny is the collector, in which it would be ideal to use negative potential to decelerate the spent electron beam to zero kinetic energy just as it reaches the collector surface. In reality this is not plausible due to a number of physical limitations, including repulsion and differing levels of kinetic energy among individual electrons. Instead, the collector is designed with multiple stages depressed below ground potential. The design of this multiple-stage collector is the optimization problem of interest. One remaining problem in SA and DSSA is the difficulty in determining when equilibrium has been reached so that the current Markov chain can be terminated. It has been suggested in recent literature that simulating the thermodynamic properties opecific heat, entropy, and internal energy from the Boltzmann distribution can provide good indicators of having reached equilibrium at a certain temperature. These properties are tested for their efficacy and implemented in SA and DSSA code with respect to TWTA collector optimization.

Radke, Eric M.↗

NASA NTRS1988

Structure of a viscoplastic theory

The general structure of a viscoplastic theory is developed from physical and thermodynamical considerations. The flow equation is of classical form. The dynamic recovery approach is shown to be superior to the hardening function approach for incorporating nonlinear strain hardening into the material response through the evolutionary equation for back stress. A novel approach for introducing isotropic strain hardening into the theory is presented, which results in a useful simplification. In particular, the limiting stress for the kinematic saturation of state (not the drag stress) is the chosen scalar-valued state variable. The resulting simplification is that there is no coupling between dynamic and thermal recovery terms in each evolutionary equation. The derived theory of viscoplasticity has the structure of a two-surface plasticity theory when the response is plasticlike, and the structure of a Bailey-Orowan creep theory when the response is creeplike.

Freed, Alan D.↗

NASA NTRS1996

Importance of the Natural Terrestrial Environment with Regard to Advanced Launch Vehicle Design and Development

The terrestrial environment is an important forcing function in the design and development of the launch vehicle. The scope of the terrestrial environment includes the following phenomena: Winds; Atmospheric Thermodynamic Models and Properties; Thermal Radiation; U.S. and World Surface Environment Extremes; Humidity; Precipitation, Fog, and Icing; Cloud Characteristics and Cloud Cover Models; Atmospheric Electricity; Atmospheric Constituents; Vehicle Engine Exhaust and Toxic Chemical Release; Occurrences of Tornadoes and Hurricanes; Geological Hazards, and Sea States. One must remember that the flight profile of any launch vehicle is in the terrestrial environment. Terrestrial environment definitions are usually limited to information below 90 km. Thus, a launch vehicle's operations will always be influenced to some degree by the terrestrial environment with which it interacts. As a result, the definition of the terrestrial environment and its interpretation is one of the key launch vehicle design and development inputs. This definition is a significant role, for example, in the areas of structures, control systems, trajectory shaping (performance), aerodynamic heating and take off/landing capabilities. The launch vehicle's capabilities which result from the design, in turn, determines the constraints and flight opportunities for tests and operations.

Pearson, S. D.↗

NASA NTRS1993

Numerical Model for Isothermal and Non-Isothermal Crystallization of Liquids and Glasses

A new numerical model of isothermal and non-isothermal first order phase transformations, such as the crystallization of liquids and glasses, is presented. This model computes directly the volume fraction transformed, taking into account time-dependent nucleation rates and cluster-size-dependent growth velocities. The model is applied to the crystallization of lithium disilicate glass, using the appropriate kinetic and thermodynamic parameters. The model is used (1) to determine the validity of common methods for computing the volume fraction transformed as a function of time in isothermal experiments when a time-dependent nucleation rate is expected, (2) to simulate non-isothermal differential scanning calorimetry (DSC) studies of glass devitrification as a function of scan rate, and (3) to compute the effect of preannealing on the DSC peak parameters. A novel behavior of the nucleation rate with scan rate is predicted, arising because the relaxation of the cluster distribution cannot be described by a single relaxation time. Comparisons of the calculations with experimental data on this glass demonstrate the validity of the model.

Kelton, K. F.↗

NASA NTRS2003

[Results of Simulation Studies]

Lattice Monte Carlo and off-lattice molecular dynamics simulations of h(sub 1)t(sub 4) and h(sub 4)t(sub l) (head/tail) amphiphile solutions have been performed as a function of surfactant concentration and temperature. The lattice and off-lattice systems exhibit quite different self-assembly behavior at equivalent thermodynamic conditions. We found that in the weakly aggregating regime (no preferred-size micelles), all models yield similar micelle size distributions at the same average aggregation number, albeit at different thermodynamic conditions (temperatures). In the strongly aggregating regime, this mapping between models (through temperature adjustment) fails, and the models exhibit qualitatively different micellization behavior. Incipient micellization in a model self-associating telechelic polymer solution results in a network with a transient elastic response that decays by a two-step relaxation: the first is due to a heterogeneous jump-diffusion process involving entrapment of end-groups within well-defined clusters and this is followed by rapid diffusion to neighboring clusters and a decay (terminal relaxation) due to cluster disintegration. The viscoelastic response of the solution manifests characteristics of a glass transition and entangled polymer network.

Source record↗

NASA NTRS1972

Melting of iron by significant structure theory

From Eyring's method of significant structures, a partition function is derived for liquid iron. The solid at high temperature is described in the Einstein approximation. Magnetic and electronic contributions to the thermodynamic properties of both the liquid and solid phases are considered. The model is compatible with properties (thermal expansion, compressibility, heat capacity, entropy of melting, and volume change on melting) at one atm. The melting temperature at high pressure is found by satisfying the requirement that the Gibbs free energies of the liquid and solid phases are equal at the melting temperature. Under conditions at the earth's core-mantle boundary, the melting temperature of iron is greater than approximately 5000 K, and under inner-outer core conditions the melting temperature is greater than approximately 7000 K. These estimates are consistent with the Lindemann melting law, but not with the Kraut-Kennedy melting law.

Leppaluoto, D. A.↗

NASA NTRS1975

Thermomechanical coupling in fatigue fracture of viscoelastic materials

Subcritical crack growth in a linear viscoelastic material subjected to cyclic loading is investigated starting with the thermodynamic power balance. Physically, it is supposed that the subcritical slow crack growth is due to local weakening of the material in the neighborhood of the crack tip. Under fatigue loading it is assumed that all energy dissipation goes into heat and that this heat build-up is the dominant irreversible process governing crack growth. The cycle averaged temperature distribution around the crack tip is obtained from local application of the first law of thermodynamics (i.e., conservation of energy). The analysis leads to the solution of a nonlinear integro-differential equation for crack length as a function of time which is coupled with the local energy equation. A regular perturbation technique is used to obtain an analytical solution which compares very well with experimental results.

Hufferd, W. L.↗

NASA NTRS1997

A General Reversible Hereditary Constitutive Model: Theoretical Developments - Part 1

Using an internal-variable formalism as a starting point, we describe the viscoelastic extension of a previously-developed viscoplasticity formulation of the complete potential structure type. It is mainly motivated by experimental evidence for the presence of rate/time effects in the so-called quasilinear, reversible, material response range. Several possible generalizations are described, in the general format of hereditary-integral representations for non-equilibrium, stress-type, state variables, both for isotropic as well as anisotropic materials. In particular, thorough discussions are given on the important issues of thermodynamic admissibility requirements for such general descriptions, resulting in a set of explicit mathematical constraints on the associated kernel (relaxation and creep compliance) functions. In addition, a number of explicit, integrated forms are derived, under stress and strain control to facilitate the parametric and qualitative response characteristic studies reported here, as well as to help identify critical factors in the actual experimental characterizations from test data that will be reported in Part II.

Saleeb, A. F.↗

NASA NTRS2018

Particle Size Calibration Testing in the NASA Propulsion Systems Laboratory

The particle size characterization portion of the 2017 Propulsion Systems Laboratory (PSL) Cloud Calibration is described. The work focuses on characterizing the particle size distribution of the icing cloud as a function of simulated atmospheric conditions.These results will aid in upcoming ice crystal and supercooled liquid icing tests in PSL. Measurements acquired with the Phase Doppler Interferometer and High Speed Imaging instruments are presented. Experimental results indicate that the particle size distribution is primarily a function nozzle air and water pressures, and that air speed is not a significant effect for ice crystal clouds in PSL and both thermodynamic conditions and air speed are not significant effects for supercooled liquid water clouds in PSL.

Engine icing↗

NASA NTRS2019

Particle Size Calibration Testing in the NASA Propulsion Systems Laboratory

Engine Icing↗

NASA NTRS1967

Thermodynamic analysis of dilute ternary systems. 1. The Ag-Au-Sn system

The heats of solution of silver and gold in dilute Ag-Au-Sn alloys have been determined as a function of alloy composition using liquid metal solution calorimetry. From the values obtained, the various enthalpy interaction parameters in the liquid solutions have been calculated. The results are interpreted in terms of the probably bonding changes taking place between silver, gold and tin atoms.

P J Spencer↗

NASA NTRS1972

Lattice model calculation of elastic and thermodynamic properties at high pressure and temperature

The elastic constants and the entire frequency spectrum were calculated up to high pressure for the alkali halides in the NaCl lattice, based on an assumed functional form of the inter-atomic potential. The quasiharmonic approximation is used to calculate the vibrational contribution to the pressure and the elastic constants at arbitrary temperature. By explicitly accounting for the effect of thermal and zero point motion, the adjustable parameters in the potential are determined to a high degree of accuracy from the elastic constants and their pressure derivatives measured at zero pressure. The calculated Gruneisen parameter, the elastic constants and their pressure derivatives are in good agreement with experimental results up to about 600 K. The model predicts that for some alkali halides the Grunesen parameter may decrease monotonically with pressure, while for others it may increase with pressure, after an initial decrease.

Demarest, H. H., Jr.↗

NASA NTRS1973

A general relaxation theory of simple liquids

A relatively simple relaxation theory to account for the behavior of liquids under dynamic conditions was proposed. The general dynamical equations are similar in form to the phenomenological relaxation equations used in theories of viscoelasticity, however, they differ in that all the coefficients of the present equations are expressed in terms of thermodynamic and molecular quantities. The theory is based on the concept that flow in a liquid distorts both the radial and the velocity distribution functions, and that relaxation equations describing the return of these functions to their isotropic distributions, characterizing a stationary liquid, can be written. The theory was applied to the problems of steady and oscillatory shear flows and to the propagation of longitudinal waves. In all cases classical results are predicted for strain rates, and an expression for the viscosity of a liquid, simular to the Macedo-Litovitz equation, is obtained.

Merilo, M.↗

NASA NTRS1991

Thermochemical analysis of chemical processes relevant to the stability and processing of SiC-reinforced Si3N4 composite

Chemical processes relevant to the stability and processing of SiC-reinforced Si3N4 composites have been examined from thermochemical considerations. The thermodynamic stabilities of various interfaces, such as SiC-Si3N4, SiC-Si3N4-Si2ON2, and SiC-Si3N4-SiO2, have been examined as a function of temperature, and the temperatures above which these interfaces become unstable have been calculated. The degradation of SiC during the processing of the composite has been examined. The processing routes considered in this study include the reaction bonded silicon nitride (RBSN) process and the pressure-assisted sintering processes with suitable sintering additives.

Misra, Ajay K.↗

NASA NTRS1997

Effect of Heat Rejection Conditions on Cryocooler Operational Stability

It is well known that cryocooler thermal efficiency is a strong function of heat rejection temperature, roughly following the dependency described by carnot. An equally important and generally overlooked implication of cryocooler heat-rejection thermodynamics is the effect of the heat rejection temparature control mode on cyrocooler performance and operational stability.

Cryocooler↗

NASA NTRS2021

Thermodynamics of the Lu2O3 – SiO2 system and comparison to other Rare Earth Silicates

Environmental barrier coatings are necessary to protect SiC based ceramics and composites from water vapor degradation in harsh engine environments. Currently, rare earth (RE) silicates are the most promising systems to protect SiC based ceramics and composites. This protection is largely based on reduced silica activity in these rare earth silicates which results in a lowered reactivity with water vapor. To that end, previous Knudsen effusion mass spectrometry (KEMS) studies have explored RE = Yb, Y silicates to measure the reduced silica activity and subsequent water vapor reactivity. Similarly, this work employs the KEMS technique to measure the SiO(g) vapor pressure in Lu containing RE silicates to calculate the activity of silica within the monosilicate [log(a(SiO2)) = -2351.1*1/T-1.6731] and the disilicate (log(a(SiO2)) = -4884.0*1/T + 2.2208) as a function of temperature. The enthalpies of formation for Lu monosilicate from the oxides and the elements were calculated to be -45±3 kJ/mol at 1550 K and -2831.1±12 kJ/mol at 298 K, respectively. The measured enthalpy of formation and those found in literature are compared to modeled values from density functional theory and those estimated using electronegativity.

Knudesen Effusion Mass Spectrometry↗