Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solution thermodynamics”

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 145 records · Page 8

Chondritic Models of 4 Vesta: Comparison of Data from the Dawn Mission with Predicted Internal Structure and Surface Composition/Mineralogy

While the HEDs provide an extremely useful basis for interpreting data from the Dawn mission, there is no guarantee that they provide a complete vision of all possible crustal (and possibly mantle) lithologies that are exposed at the surface of Vesta. With this in mind, an alternative approach is to identify plausible bulk compositions and use mass-balance and geochemical modelling to predict possible internal structures and crust/mantle compositions and mineralogies. While such models must be consistent with known HED samples, this approach has the potential to extend predictions to thermodynamically plausible rock types that are not necessarily present in the HED collection. Nine chondritic bulk compositions are considered (CI, CV, CO, CM, H, L, LL, EH, EL). For each, relative proportions and densities of the core, mantle, and crust are quantified. This calculation is complicated by the fact that iron may occur in metallic form (in the core) and/or in oxidized form (in the mantle and crust). However, considering that the basaltic crust has the composition of Juvinas and assuming that this crust is in thermodynamic equilibrium with the residual mantle, it is possible to calculate a single solution to this problem for a given bulk composition. Of the nine bulk compositions tested, solutions corresponding to CI and LL groups predicted a negative metal fraction and were not considered further. Solutions for enstatite chondrites imply significant oxidation relative to the starting materials and these solutions too are considered unlikely. For the remaining bulk compositions, the relative proportion of crust to bulk silicate is typically in the range 15 to 20% corresponding to crustal thicknesses of 15 to 20 km for a porosity-free Vesta-sized body. The mantle is predicted to be largely dominated by olivine (greater than 85%) for carbonaceous chondrites, but to be a roughly equal mixture of olivine and pyroxene for ordinary chondrite precursors. All bulk compositions have a significant core, but the relative proportions of metal and sulphide can be widely different. Using these data, total core size (metal+ sulphide) and average core densities can be calculated, providing a useful reference frame within which to consider geophysical/gravity data of the Dawn mission. Further to these mass-balance calculations, the MELTS thermodynamic calculator has been used to assess to what extent chondritic bulk compositions can produce Juvinas-like liquids at relevant degrees of partial melting/crystallization. This work will refine acceptable bulk compositions and predict the mineralogy and composition of the associated solid and liquid products over wide ranges of partial melting and crystallization, providing a useful and self-consistent reference frame for interpretation of the data from the VIR and GRaND instruments onboard the Dawn spacecraft.

Toplis, M. J.↗

COMOC: a finite-element algorithm for the Navier-Stokes equations

Finite-element algorithm devised to facilitate solution of two-dimensional Navier-Stokes equations governing kinematics and thermodynamics of variable-viscosity, compressible, multiple-species fluid. Algorithm has been implemented into existing computer program system.

Baker, A. J.↗

Ocean surface pressure fields from satellite-sensed winds

The University of Washington's planetary boundary layer model is inverted to use remotely sensed satellite scatterometer-derived surface winds as input to calculate maritime surface pressure fields. The analysis of three different synoptic storm situations is performed using the model and is then compared to conventional National Weather Service analyses. Agreement is good. Isolation of the PBL secondary flow, stratification and thermal wind effects in the model revealed that each may be significant under certain conditions. However, the model shows sensitivity to the thermodynamic features only in a general sense and even a neutral stratification solution gives a good approximation. The high density of the scatterometer data produces mesoscale (hundreds of kilometers) dynamic details, which cannot be confirmed by conventional data at this stage.

Brown, Robert A.↗

Extensions of the CSCM methodology for nonequilibrium reacting gas flows

The paper describes the extension of the CSCM flux difference splitting to produce conservative unconditionally stable implicit upwind methods for the extended equation set. The CSCM-S single data level relaxation algorithm is argued to be particularly advantageous for efficiently solving the set of equations. Attention is given to issues of alternative state variable representations and to complementary thermodynamic relations that parameterize and need to be effectively coupled with the solution procedure for the partial differential equations. Within the framework of iteration, consideration is given to more efficient approximate solution procedures than the full block coupled LU decomposition. Results of early numerical experiments are given to illustrate the developments.

Lombard, C. K.↗

Operator split methods in the numerical solution of the finite deformation elastoplastic dynamic problem

The spatial formulation of the elastoplastic dynamic problem for finite deformations is considered. A thermodynamic argument leads to an additive decomposition of the spatial rate of deformation tensor and allows an operator split of the evolutionary equations of the problem into elastic and plastic parts. This operator split is taken as the basis for the definition of a global product algorithm. In the context of finite element discretization the product algorithm entails, for every time step, the solution of a nonlinear elastodynamic problem followed by the application of plastic algorithms that operate on the stresses and internal variables at the integration points and bring in the plastic constitutive equations. Suitable plastic algorithms are discussed for the cases of perfect and hardening plasticity and viscoplasticity. The proposed formalism does not depend on any notion of smoothness of the yield surface and is applicable to arbitrary convex elastic regions, with or without corners. The stabiity properties of the global product algorithm are shown to be identical to those of the algorithm used for the integration of the nonlinear elastodynamic problem. Numerical examples illustrate the accuracy of the method.

Pinsky, P. M.↗

The computation of non-equilibrium chemically-reacting flows

Numerical techniques for the analysis of flows with nonequilibrium thermodynamics and chemistry are developed and demonstrated. Particular attention is given to the formulation of the thermodynamic and chemical gas models, the derivation of the governing equations, and the solution procedures. Results for typical problems involving shock tubes and supersonic diffusers are presented in extensive graphs and briefly characterized.

Grossman, B.↗

Integration of satellite and surface data using a radiative-convective oceanic boundary-layer model

A mixing-line boundary-layer model is used to retrieve cloud-top height from satellite-derived cloud-top temperatures, using 700-hPa National Meteorological Center (NMC) analyses and the Comprehensive Ocean and Atmosphere Data Set (COADS) surface data as supporting datasets. Results are compared with the fixed-lapse-rate method of retrieving boundary-layer depth from sea surface temperatures (SST) and cloud-top temperatures. A radiative-convective equilibrium boundary-layer model is used to retrieve boundary-layer structure given SST and surface wind, satellite cloud-top temperatures and cloud fraction, and the 700-hPa NMC thermodynamic analyses. Good agreement is found between the COADS data and the model solutions for low-level temperature and moisture. This suggests that equilibrium boundary-layer models may be of use over remote oceans in the retrieval of boundary-layer structure.

Betts, Alan K.↗

Theoretical spectra of circumstellar dust shells around carbon-rich asymptotic giant branch stars

Realistic modeling of circumstellar dust shells around evolved stars has to include a physical treatment of the interaction among hydrodynamics, thermodynamics, radiative transfer, chemistry and dust formation and -growth. A self-consistent solution of this problem is presented in the case of stationary, spherical symmetric dust-driven winds. The resulting shell structure and the mass-loss rate are completely determined by the three fundamental stellar parameters stellar mass M(stellar), stellar luminosity L(stellar) and effective temperature T(sub eff) and by the element abundances epsilon(sub i). A detailed calculation of the transport coefficients of the dust component by means of the particle size distribution function and the solution of the non-grey radiative transfer problem provide realistic synthetic spectra of the dust shell models. We discuss the dependence of the resulting spectra on the stellar parameters in terms of infrared two color diagrams. Application of these model calculations to the prominent infrared object IRC +10216 yields best agreement with the observed spectrum and the visibility data at maximum light for the stellar parameters M(stellar) = 0.7 solar mass, L(stellar) = 2.4 x 10(exp 4) solar luminosity, T(stellar) = 2010K and a carbon to oxygen ratio of epsilon(sub c)/epsilon(sub o) = 1.40, which corresponds to a mass-loss rate of M-dot = 8 x 10(exp -5) solar mass/yr. In this model only amorphous carbon grains are considered as the main opacity source. From this model a distance to IRC +10216 of d = 170pc is deduced. The total mass contained in the circumstellar dust shell implies and initial main sequence mass of M(sub ZAMS) greater than or = 1.3 solar mass for IRC +10216.

Winters, J. M.↗

Real-Gas Effects on Binary Mixing Layers

This paper presents a computational study of real-gas effects on the mean flow and temporal stability of heptane/nitrogen and oxygen/hydrogen mixing layers at supercritical pressures. These layers consist of two counterflowing free streams of different composition, temperature, and density. As in related prior studies reported in NASA Tech Briefs, the governing conservation equations were the Navier-Stokes equations of compressible flow plus equations for the conservation of total energy and of chemical- species masses. In these equations, the expressions for heat fluxes and chemical-species mass fluxes were derived from fluctuation-dissipation theory and incorporate Soret and Dufour effects. Similarity equations for the streamwise velocity, temperature, and mass fractions were derived as approximations to the governing equations. Similarity profiles showed important real-gas, non-ideal-mixture effects, particularly for temperature, in departing from the error-function profile, which is the similarity solution for incompressible flow. The temperature behavior was attributed to real-gas thermodynamics and variations in Schmidt and Prandtl numbers. Temporal linear inviscid stability analyses were performed using the similarity and error-function profiles as the mean flow. For the similarity profiles, the growth rates were found to be larger and the wavelengths of highest instability shorter, relative to those of the errorfunction profiles and to those obtained from incompressible-flow stability analysis. The range of unstable wavelengths was found to be larger for the similarity profiles than for the error-function profiles

Okong'o, Nora↗

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗

An objective method for determining the generalized transport tensor for two-dimensional Eulerian models

An objective method for deriving the components of a generalized transport tensor for a two-dimensional model is presented. Representative meridional and vertical velocities and thermodynamic scalars at a uniform grid are used to reduce the problem to the solution of two flux equations for two unknowns. One unknown is the stream-function, coefficient of an antisymmetric tensor, which corrects the Eulerian mean motions for Stokes drift; the other is a time constant, which converts the deviatory velocity tensor to a symmetric transport tensor. The complete asymmetric tensor, called a transport tensor, has a divergence which yields both advection and diffusion by the deviatory velocities. Advantages and disadvantages of Lagrangian and Eulerian averages are discussed, and meridional-vertical velocity correlations are provided.

Danielsen, E. F.↗

Least squares inversion of Stokes profiles in the presence of velocity gradients

The Auer, Heasley and House Stokes inversion procedure in use at High Altitude Observatory is based on the analytic solution of the equation of transfer for polarized light where the representation of the thermodynamic and magnetic structure of the atmosphere is assumed to have a high degree of invariance, namely, a Milne-Eddington (ME) structure with a constant magnetic field. In the presence of invariance breaking gradients the resultant Stokes profiles are represented only approximately, if at all, by analytic forms. The accuracy of the inversion parameters and their significance as measures of actual structure are explored for the ME and the Landman-Finn sunspot models under the effects of velocity gradients. The resulting field parameters are good to a few percent and prove to be insensitive to the errors committed by the use of a ME-representation, but the resulting ME parameters yield a less precise measure of thermal structure.

Skumanich, A.↗

Estimated heats of fusion of fluoride salt mixtures suitable for thermal energy storage applications

The heats of fusion of several fluoride salt mixtures with melting points greater than 973 K were estimated from a coupled analysis of the available thermodynamic data and phase diagrams. Simple binary eutectic systems with and without terminal solid solutions, binary eutectics with congruent melting intermediate phases, and ternary eutectic systems were considered. Several combinations of salts were identified, most notable the eutectics LiF-22CaF2 and NaF-60MgF2 which melt at 1039 and 1273 K respectively which posses relatively high heats of fusion/gm (greater than 0.7 kJ/g). Such systems would seemingly be ideal candidates for the light weight, high energy storage media required by the thermal energy storage unit in advanced solar dynamic power systems envisioned for the future space missions.

Misra, A. K.↗

Morphological control of inter-penetrating polymer networks

Synthetic organic polymer chemistry has been successful in producing composition of matter with thermal oxidation stability and progressively higher glass transition temperatures. In part, this was done by increasing the steric-hindrance of moieties in the chain of a macromolecule. The resulting polymers are usually quite insoluble and produce molten polymers of very high viscosities. These types of polymers are not easily processed into graphite fiber prepregs by melt or solution impregnation methods. Hence, a technological need exists to produce new knowledge of how to produce polymer-fiber composites from this class of polymers. The concept of freeze drying amic-acid prepolymers with reactive thermoplastic was proposed as a research topic for the ASEE/NASA Summer Faculty Program of 1989 as a means of producing polymer-fiber composites. This process scheme has the thermodynamic attribute that the magnitude of phase separation due to differences in solubility of two organic constituents in solution will be greatly reduced by removing a solvent not by evaporation but by sublimation. Progress to date on evaluating this polymer processing concept is briefly outlined.

Hansen, Marion↗

Thermodynamic constraints on the textural evolution of eucrite, EET 90020

Basaltic eucrites, which are possibly parts of the Vestan crust, formed as either lava flows or intrusions, and can offer insights into early crust formation. However most eucrites have experienced thermal metamorphism, which exacerbates the challenges of understanding these samples. A variety of petrogenetic models such as, partial melting of a primitive source, fractional crystallization of magmas emplaced in a primitive crust, and/or partial melting of a eucritic source coupled with melt mixing and assimilation have been proposed in order to explain observed petrologic and chemical characteristics of eucrites. EET 90020 is an unbrecciated eucrite that has experienced significant thermal metamorphism, with temperatures of metamorphic equilibration calculated to be ~840-1040°C. Its petrologic history remains contentious, in part due to differences in trace element analyses [10,17, this work], which has led to multiple metamorphic interpretations [8,10-11]. Here, we combined petrologic observations, chemical analyses, and thermodynamic modeling, to interpret micro-domain textures identified in EET 90020 and better constrain its petrologic history. Ultimately, the development of such textures is a direct result of the geologic processes operating on a young Vesta or similar asteroid. Sample description: We identified three textural domains in EET 90020 (Figure 1); (I) A coarse grain domain with granoblastic plagioclase and pyroxene, (II) a medium grain domain with curved grain boundaries between plagioclase and pyroxene, and (III) a fine grain domain dominated by tridymite that fills in interstitial space between spherical plagioclase. Major element chemistry in phases is consistent across domains. Pyroxene have distinctive high Ca lamellae (Wo40.92En28.1Fs31.0) forming from the low-Ca pigeonite host (Wo3.5En33.7Fs62.8). There are small amounts of fayalitic olivine along pyroxene/oxide grain boundaries, where pyroxenes have reacted with ilmenite/chromite clasts. Plagioclase is anorthitic (Avg ~An88), with little variation (An86-92.5). Trace element bulk rock analyses were collected via ICP-MS and individual mineral analyses with LA-ICP-MS (Figure 2). The bulk sample has an observable enrichment in light rare earth elements, with a notable depletion in Eu. Plagioclase is enriched in LREEs, depleted in HREEs, and has a positive Eu anomaly while pyroxene is enriched in HREEs with a negative Eu anomaly. Variations in trace element analyses is likely influenced by the presence of phosphates because are highly enriched in REEs and have a dramatic effect on the bulk domain compositions. Thermodynamic modeling: The software package, Perple_X, was used to calculate mineral phase equilibria over a range of conditions using a Gibbs free energy minimization approach [12]. Models used all major/minor elements except P2O5. Thermodynamic properties from [13] were used to constrain endmember phase stabilities. Activity models were used to describe mixing in phases with solid solution. Isochemical P-T phase diagrams were constructed for the bulk thin section, and the coarse, medium and fine grain domain compositions. Domain compositions were determined using microprobe analyses of phases identified in each domain and observed modal mineral abundancies. Thermodynamically stable pyroxene compositions were extracted from model results and compared to pyroxene compositions collected via EMPA, allowing us to calculate temperatures for metamorphic equilibrium. Results: Results that use the coarse domain composition (Fig. 3a) demonstrate overlap at T ~ 1020°C indicating high T equilibration. There was no overlap for the medium and fine grain domains, or the bulk composition (e.g., Fig. 3b) indicating that these are not equilibrium assemblages. These results were replicated using data from [8,20]. Discussion: The disequilibrium results from the bulk rock, fine grain and medium grain models imply that EET 90020 was modified during and/or after peak metamorphism (open system behavior). However, this is in contrast to the coarse domain where it appears as if equilibrium was maintained. We suggest that this conflict is explained by highly localized equilibrium occurring at millimeter length-scales, which has been observed in terrestrial metamorphism where aqueous fluid infiltrate [14] or melt loss occurs [15]. It is unlikely that aqueous alteration was significant in modifying the bulk composition of EET 90020, because it lacks hydrated minerals. However, the presence of melt is supported by the presence of an abundance of curved grain boundaries and spherical grains in the medium and fine grain domains This contrasts the coarse domain, where grain boundaries are angular. Given the evidence for disequilibrium in the medium and fine grain domains, it would be inappropriate to draw additional conclusions about textural evolution using thermodynamic models of those domains. However, modeling results of the coarse grain domain, where equilibrium was maintained, can provide useful insights into the geologic evolution of EET 90020. We approximate maximum metamorphic temperatures, using the Ca component in pyroxene, at ~1020°C (Fig. 4), which is consistent with previous two-pyroxene thermometry [11]. At this temperature, up to 10 vol % melt can be produced in the coarse grain domain, which is at the boundary of minimum melt needed in order to segregate and form a melt network [16]. Tridymite and plagioclase would be the first minerals to melt out during heating starting ~990°C. This is consistent with the observed melt textures formed by tridymite and plagioclase in the fine grain domain. We suggest that a partial melting model best explains the textural development of EET 90020. Upon heating, melt generated was either trapped, resulting in no net change in the localized rock composition, or migrated, resulting in disequilibrium between pyroxene and the surrounding matrix. A positive Eu anomaly in plagioclase indicates that all of the plagioclase was melted previously. This is consistent with previous studies that have suggested that EET 90020 represents a residual eucrite that experienced partial melting and subsequent melt loss [8]. We speculate that partial melting could have occurred when EET 90020 was heated in the asteroid’s lower crust or adjacent to an intrusion [7]. Future work will study other, texturally heterogeneous samples to see if EET 90020 presents a unique circumstance or whether crustal partial melting was ubiquitous on the eucrite asteroid.

Eucrites↗

Variational energy principle for compressible, baroclinic flow. 2: Free-energy form of Hamilton's principle

The first and second variations are calculated for the irreducible form of Hamilton's Principle that involves the minimum number of dependent variables necessary to describe the kinetmatics and thermodynamics of inviscid, compressible, baroclinic flow in a specified gravitational field. The form of the second variation shows that, in the neighborhood of a stationary point that corresponds to physically stable flow, the action integral is a complex saddle surface in parameter space. There exists a form of Hamilton's Principle for which a direct solution of a flow problem is possible. This second form is related to the first by a Friedrichs transformation of the thermodynamic variables. This introduces an extra dependent variable, but the first and second variations are shown to have direct physical significance, namely they are equal to the free energy of fluctuations about the equilibrium flow that satisfies the equations of motion. If this equilibrium flow is physically stable, and if a very weak second order integral constraint on the correlation between the fluctuations of otherwise independent variables is satisfied, then the second variation of the action integral for this free energy form of Hamilton's Principle is positive-definite, so the action integral is a minimum, and can serve as the basis for a direct trail and error solution. The second order integral constraint states that the unavailable energy must be maximum at equilibrium, i.e. the fluctuations must be so correlated as to produce a second order decrease in the total unavailable energy.

Schmid, L. A.↗

A Generalized Fluid Formulation for Turbomachinery Computations

A generalized formulation of the equations of motion of an arbitrary fluid are developed for the purpose of defining a common iterative algorithm for computational procedures. The method makes use of the equations of motion in conservation form with separate pseudo-time derivatives used for defining the numerical flux for a Riemann solver and the convergence algorithm. The partial differential equations are complemented by an thermodynamic and caloric equations of state of a complexity necessary for describing the fluid. Representative solutions with a new code based on this general equation formulation are provided for three turbomachinery problems. The first uses air as a working fluid while the second uses gaseous oxygen in a regime in which real gas effects are of little importance. These nearly perfect gas computations provide a basis for comparing with existing perfect gas code computations. The third case is for the flow of liquid oxygen through a turbine where real gas effects are significant. Vortex shedding predictions with the LOX formulations reduce the discrepancy between perfect gas computations and experiment by approximately an order of magnitude, thereby verifying the real gas formulation as well as providing an effective case where its capabilities are necessary.

Merkle, Charles L.↗