Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic Modeling”

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 289 records · Page 16

Thermal Design for Extra-Terrestrial Regenerative Fuel Cell System

The Advanced Exploration Systems (AES) Advanced Modular Power Systems (AMPS) Project is investigating different power systems for various lunar and Martian mission concepts. The AMPS Fuel Cell (FC) team has created two system-level models to evaluate the performance of regenerative fuel cell (RFC) systems employing different fuel cell chemistries. Proton Exchange Membrane fuel cells PEMFCs contain a polymer electrolyte membrane that separates the hydrogen and oxygen cavities and conducts hydrogen cations (protons) across the cell. Solid Oxide fuel cells (SOFCs) operate at high temperatures, using a zirconia-based solid ceramic electrolyte to conduct oxygen anions across the cell. The purpose of the modeling effort is to down select one fuel cell chemistry for a more detailed design effort. Figures of merit include the system mass, volume, round trip efficiency, and electrolyzer charge power required. PEMFCs operate at around 60 C versus SOFCs which operate at temperatures greater than 700 C. Due to the drastically different operating temperatures of the two chemistries the thermal control systems (TCS) differ. The PEM TCS is less complex and is characterized by a single pump cooling loop that uses deionized water coolant and rejects heat generated by the system to the environment via a radiator. The solid oxide TCS has its own unique challenges including the requirement to reject high quality heat and to condense the steam produced in the reaction. This paper discusses the modeling of thermal control systems for an extraterrestrial RFC that utilizes either a PEM or solid oxide fuel cell.

thermal modeling↗

High-resolution SOFIA/EXES Spectroscopy of SO2 Gas in the Massive Young Stellar Object MonR2 IRS3: Implications for the Sulfur Budget

Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the ν3 rovibrational band of SO2 obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO2 gas in the massive YSO MonR2 IRS3. This yields a SO2/H abundance lower limit of 5.6 ± 0.5 × 10(exp -7), or >4% of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of b < 3.20 km s(exp -1). The small line widths and high temperature (Tex = 234 ± 15 K) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov–Gerasimenko), which is followed by SO2 formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO2 from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.

Infrared ISM↗

Computational Tools to Design Molten Salt Electrolytes for High-Temperature Environments

Current Li-ion batteries are designed for a small operating window of 5 °C to 55 °C. Modifications to the electrolyte for operation in Mars atmosphere extended this range on the lower bound to ~ -40 °C. These operating-temperature ranges are far from those presented in high-temperature environments, such as the Venus surface, where temperatures are around 450 °C. Protecting the state-of-the-art Li-ion batteries require insulations that decreases the volumetric capacity and limiting the operational time-window. These unique challenges require a paradigm shift in materials used for designing high-temperature batteries. Molten salt electrolyte-based batteries offer a plausible route to designing high-temperature batteries. ZEBRA batteries are known to be one of the safest energy storage devices operating at 270-350 °C. Nitrate based eutectics, operating at 150 °C have also been used in Li-O2 batteries. In this study, we will examine molten-salt electrolyte transport and electrochemical properties using first-principles computations and benchmark against experiments. Further, we will present thermodynamics-based models for designing and predicting melting point of molten salt mixtures. Based on these computational tools, new molten salt mixtures designed with desired operating temperatures and electrochemical windows will be presented. Implications of these new eutectics in the context of high-temperature environment exploration will also be discussed.

Balachandran Radhakrishnan↗

Molten Salt Electrolytes: Computational Analysis of Transport, Electrochemical Properties and Designing New Mixtures

Electric aircraft propulsion has gained a traction over the last decade due to possible high-energy density electrochemistries with reasonable cycle life and identification of non-flammable chemistries that can guarantee much better safety. Commercial Li-ion batteries cannot reach such high capacities because of weight limitations that arise from the widely used intercalation electrodes. Also, use of organic electrolytes make them highly susceptible to fire on exposure to air and humidity. Currently, use of Lithium metal anode along with conversion electrochemistries such as Li-O2 and Li-S are being pursued to achieve such high energy densities. Safer inorganic solid-state electrolytes, recently discovered Water-in-Salt electrolytes, molten salt electrolytes are some of the alternatives being pursued for a safer/non-flammable battery. Of these safer alternatives, molten-salt electrolytes offer some attractive properties: liquid at operating temperatures resulting in better electrode-wetting, stable interface with Li-metal anodes and good ionic conductivity; their only drawback being high operating temperatures. In this talk, we will discuss some of the computational studies, driven via atomistic simulations, of the properties of molten-salt electrolytes including transport mechanism and interface stability. We will also discuss predicting electrochemical properties of these high-temperature electrolytes. Further, we will discuss a thermodynamic approach to predicting new molten-salt eutectics with lower melting points.

Molten Salt Electrolytes↗

Trace Element Partitioning Between Olivine and Melt in Lunar Basalts

Mineral/melt partition coefficients have been widely used to provide insights into magmatic processes. Olivine is one of the most abundant and important minerals in the lunar mantle and mare basalts. Yet, no systematic olivine/melt partitioning data are available for lunar conditions. We report trace element partition data between host mineral olivine and its melt inclusions in lunar basalts. Equilibrium is evaluated using the Fe-Mg exchange coefficient, leading to the choice of melt inclusion-host olivine pairs in lunar basalts 12040, 12009, 15016, 15647, and 74235. Partition coefficients of 21 elements (Li, Mg, Al, Ca, Ti, V, Cr, Mn, Fe, Co, Y, Zr, Nb, Gd, Tb, Dy, Ho, Er, Tm, Yb and Lu) were measured. Except for Li, V, and Cr, these elements show no significant difference in olivine-melt partitioning compared to the data for terrestrial samples. The partition coefficient of Li between olivine and melt in some lunar basalts with low Mg# (Mg# < 0.75 in olivine, or < ~0.5 in melt) is higher than published data for terrestrial * Corresponding author. Email address: youxue@umich.edu 2 samples, which is attributed to the dependence of DLi on Mg# and the lack of literature DLi data with low Mg#. The partition coefficient of V in lunar basalts is measured to be 0.17 to 0.74, significantly higher than that in terrestrial basalts (0.003 to 0.21), which can be explained by the lower oxygen fugacity in lunar basalts. The significantly higher DV can explain why V is less enriched in evolved lunar basalts than terrestrial basalts. The partition coefficient of Cr between olivine and basalt melt in the Moon is 0.11 to 0.62, which is lower than those in terrestrial settings by a factor of approximately 2. This is surprising because previous authors showed that Cr partition coefficient is independent of fO2. A quasi-thermodynamically based model is developed to correlate Cr partition coefficient to olivine and melt composition and fO2. The lower Cr partition coefficient between olivine and basalt in the Moon can lead to more Cr enrichment in the lunar magma ocean, as well as more Cr enrichment in mantle-derived basalts in the Moon. Hence, even though Cr is typically a compatible element in terrestrial basalts, it is moderately incompatible in primitive lunar basalts, with a similar degree of incompatibility as V based on partition coefficients in this work, as also evidenced by the relatively constant V/Cr ratio of 0.039 ± 0.011 in lunar basalts. The confirmation of constant V/Cr ratio is important for constraining concentrations of Cr (slightly volatile and siderophile) and V (slightly siderophile) in the bulk silicate Moon.

partition coefficients↗

A persistently low level of atmospheric oxygen in Earth’s middle age

Resolving how Earth surface redox conditions evolved through the Proterozoic Eon is fundamental to understanding how biogeochemical cycles have changed through time. The redox sensitivity of cerium relative to other rare earth elements and its uptake in carbonate minerals make the Ce anomaly (Ce/Ce*) a particularly useful proxy for capturing redox conditions in the local marine environment. Here, we report Ce/Ce* data in marine carbonate rocks through 3.5 billion years of Earth’s history, focusing in particular on the mid-Proterozoic Eon (i.e., 1.8 – 0.8 Ga). To better understand the role of atmospheric oxygenation, we use Ce/ Ce* data to estimate the partial pressure of atmospheric oxygen (pO2) through this time. Our thermodynamics-based modeling supports a major rise in atmospheric oxygen level in the aftermath of the Great Oxidation Event (~ 2.4 Ga), followed by invariant pO2 of about 1% of present atmospheric level through most of the Proterozoic Eon (2.4 to 0.65 Ga).

Xiao-Ming Liu↗

List of Commercial and Advanced Developmental Niobium-Based Alloys of the Space Age

This report compiles a list of commercial and developmental niobium-based alloys developed during the Space Age (late 1950s through mid-1970s) for extreme-temperature applications, including rocket engine thrust chambers, hypersonic re-entry thermal protection systems, and space fission reactor loops. Niobium (Nb) was widely pursued because it provided the lowest density (~8.6 g/cc) among the primary refractory metals, a high melting temperature (~2470°C), exceptional low-temperature ductility, good formability, and compatibility with liquid alkali metals. An evaluation of physical metallurgy mechanisms, focusing on solid-solution strengthening via heavy refractory solutes (W, Mo, Ta), dual-purpose reactive solutes (Hf, Zr, Ti), and dispersion strengthening using carbides, nitrides, and oxides is presented. Additionally, the report compares Western and Soviet Union metallurgical approaches, explaining how supply chain factors and manufacturing infrastructure influenced element selection, interstitial chemistry, and alloy identification/naming conventions. Cataloging these historical alloy chemical compositions serves as a foundational reference for modern alloy additive manufacturing, thermodynamic CALPHAD modeling, and machine-learning discovery pipelines for next-generation extreme-environment niobium-based alloys.

Physical Metallurgy↗

Modeling ARRM Xenon Tank Pressurization Using 1D Thermodynamic and Heat Transfer Equations

As a first step in understanding what ground support equipment (GSE) is required to provide external cooling during the loading of 5,000 kg of xenon into 4 aluminum lined composite overwrapped pressure vessels (COPVs), a modeling analysis was performed using Microsoft Excel. The goals of the analysis were to predict xenon temperature and pressure throughout loading at the launch facility, estimate the time required to load one tank, and to get an early estimate of what provisions for cooling xenon might be needed while the tanks are being filled. The model uses the governing thermodynamic and heat transfer equations to achieve these goals. Results indicate that a single tank can be loaded in about 15 hours with reasonable external coolant requirements. The model developed in this study was successfully validated against flight and test data. The first data set is from the Dawn mission which also utilizes solar electric propulsion with xenon propellant, and the second is test data from the rapid loading of a hydrogen cylindrical COPV. The main benefit of this type of model is that the governing physical equations using bulk fluid solid temperatures can provide a quick and accurate estimate of the state of the propellant throughout loading which is much cheaper in terms of computational time and licensing costs than a Computation Fluid Dynamics (CFD) analysis while capturing the majority of the thermodynamics and heat transfer.

Xenon storage↗

Effects of precipitation on the thermodynamic structure of the trade wind boundary layer

A model of the thermodynamic structure of the trade wind boundary layer is formulated to include the parameterization of precipitation in relatively shallow clouds. Although the area-averaged simulated precipitation rates are relatively small (less than 1 mm/day), the inclusion of precipitation has an appreciable effect on the predicted thermodynamic structure. The cloud layer structure simulated with precipitation is warmer, drier, and more stable than that simulated without precipitation. The simulated inversion height is lowered by as much as 60 mbar when precipitation is included.

Albrecht, Bruce A.↗

Specific heats and thermodynamic critical fields in Zn-doped YBa2Cu3O(7-x) according to an induced-pairing model

Electronic specific heats and thermodynamic critical fields are calculated in a mean-field version of an induced-pairing model for superconductivity, and compared with results of Loram et al. (1990) on YBa2(Cu(1-y)Zn(y))3O(7-x). This model involves induction of pairing of holes in a wideband by strongly bound electronlike pairs. It is assumed that the planar hole concentration for no Zn addition is close to, but slightly higher than, that for the maximum Tc, and that it increases by 0.015 per planar Cu ion for each increase of y by 0.01. Parameters of the model are taken to be the same as in a previous publication in which energy gaps were discussed, except that an effective hybridization parameter is adjusted for each Zn concentration to give agreement with the observed Tc. Results are presented for y = 0.0, 0.01, and 0.03. The agreement with experiment is good for thermodynamic critical fields, and is fair for specific heats. For specimens with larger y, with relatively low T(c)s, it is argued that the model should be supplemented to include effects of a BCS-type interaction amongst the wideband carriers.

Eagles, D. M.↗

Propulsion dynamics of lunar hoppers

A feasibility study was recently completed on a hopping transporter concept for locomotion on the moon. Termed a Lunar Pogo because it would operate similarly to a conventional pogo stick, this vehicle would accelerate up an inclined leg, pick up the leg, and then enter ballistic flight over the moon's surface. Upon recontacting the lunar surface, the Lunar Pogo would decelerate down the leg. Propulsion would be provided by expansion of gas against a piston. Operation would be partially conservative because much of the energy expended by the gas during takeoff would be recovered by compressing the gas during landing. Two models of the ballistics and propulsion have been set up to estimate performance. The simplified first-order model illuminates the thermodynamics of ideal operation. The second-order model, which includes realistic effects such as sliding, provides a better approximation to actual performance. In a one-man Lunar Pogo, an astronaut would typically make 15-m (50-ft) leaps at an average speed of 5-8 km (3-5 miles) per hour.

Meetin, R. J.↗

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↗

Hadron thermodynamics in relativistic nuclear collisions

Various phenomenological models based on statistical thermodynamical considerations were used to fit the experimental data at high P sub T to a two temperature distribution. Whether this implies that the two temperatures belong to two different reaction mechanisms, or consequences of Lorentz-contraction factor, or related in a fundamental way to the intrinsic thermodynamics of Space-Time can only be revealed by further theoretical and experimental investigations of high P sub T phenomena in extremely energetic hadron-hadron collisions.

Ammiraju, P.↗

A Damage Model for the Simulation of Delamination in Advanced Composites under Variable-Mode Loading

A thermodynamically consistent damage model is proposed for the simulation of progressive delamination in composite materials under variable-mode ratio. The model is formulated in the context of Damage Mechanics. A novel constitutive equation is developed to model the initiation and propagation of delamination. A delamination initiation criterion is proposed to assure that the formulation can account for changes in the loading mode in a thermodynamically consistent way. The formulation accounts for crack closure effects to avoid interfacial penetration of two adjacent layers after complete decohesion. The model is implemented in a finite element formulation, and the numerical predictions are compared with experimental results obtained in both composite test specimens and structural components.

Turon, A.↗

Modeling Xenon Tank Pressurization using One-Dimensional Thermodynamic and Heat Transfer Equations

As a first step in understanding what ground support equipment (GSE) is required to provide external cooling during the loading of 5,000 kg of xenon into 4 aluminum lined composite overwrapped pressure vessels (COPVs), a modeling analysis was performed using Microsoft Excel. The goals of the analysis were to predict xenon temperature and pressure throughout loading at the launch facility, estimate the time required to load one tank, and to get an early estimate of what provisions for cooling xenon might be needed while the tanks are being filled. The model uses the governing thermodynamic and heat transfer equations to achieve these goals. Results indicate that a single tank can be loaded in about 15 hours with reasonable external coolant requirements. The model developed in this study was successfully validated against flight and test data. The first data set is from the Dawn mission which also utilizes solar electric propulsion with xenon propellant, and the second is test data from the rapid loading of a hydrogen cylindrical COPV. The main benefit of this type of model is that the governing physical equations using bulk fluid solid temperatures can provide a quick and accurate estimate of the state of the propellant throughout loading which is much cheaper in terms of computational time and licensing costs than a Computation Fluid Dynamics (CFD) analysis while capturing the majority of the thermodynamics and heat transfer.

Thermodynamic tank modeling↗

An Interface Damage Model for the Simulation of Delamination Under Variable-Mode Ratio in Composite Materials

A thermodynamically consistent damage model for the simulation of progressive delamination under variable mode ratio is presented. The model is formulated in the context of the Damage Mechanics (DM). The constitutive equations that result from the variation of the free energy with damage are used to model the initiation and propagation of delamination. A new delamination initiation criterion is developed to assure that the formulation can account for changes in the loading mode in a thermodynamically consistent way. Interfacial penetration of two adjacent layers after complete decohesion is prevented by the formulation of the free energy. The model is implemented into the commercial finite element code ABAQUS by means of a user-written decohesion element. Finally, the numerical predictions given by the model are compared with experimental results.

Turon, Albert↗