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.
On-Orbit Cryogenic Refueling: Potential Mission Benefits, Associated Orbital Mechanics, and Fuel Transfer Thermodynamic Modeling Efforts
The placement of cryogenic fuel/propellant depot stations in Earth orbit has the potential to transform the nature and operations for many types of spaceflight missions. Today, spaceflight missions are almost universally required to carry the entire amount of fuel required for the mission, for the entire duration of the mission, from the point of launch. This is the rough equivalent of making a drive from Ohio to California, requiring the traveler to bring along the total sum of gasoline required for the entire trip, without being able to ‘fill-up’ anywhere along the route. Obviously, this framework of travel greatly encumbers the breadth, scope, and efficiency of potential journeys. Cryogenic fuel/propellant depots have not been implemented because many technical, operational, and engineering challenges still exist. These must be overcome prior to the placement of usable on-orbit propellant depots. This thesis investigates three specific engineering challenges related to on-orbit propellant depots, and presents the current state, technological challenges, and ultimate benefits of on-orbit cryogenic refueling.
VESIcal: An Open-source Thermodynamic Model Engine for Mixed Volatile (H2O-CO2) Solubility in Silicate Melts
Modeling the solubility of volatiles in silicate melts is fundamental to the interpretation of volcanic systems and has implications for magma dynamics, eruption style, and material transport between the mantle, crust, and atmosphere. Recent advancements in computational capabilities and access to computing tools has outpaced the functionality and extensibility of previously available modeling platforms. Here we present VESIcal (Volatile Equilibria and Saturation Index calculator), the first comprehensive modeling tool for H2O, CO2, and mixed (H2O-CO2) solubility in silicate melts that: a) allows users access to seven popular models, with easy inter-comparison between models; b) provides universal functionality for all models (e.g., functions for calculating saturation pressures, degassing paths, etc.); c) can process large datasets (1,000’s of samples) automatically; d) can output computed data into an Excel spreadsheet or CSV file for post-modeling analysis; e) integrates plotting capabilities directly within the tool; and f) provides all of this within the framework of a python library, making the tool extensible by the user and allowing any of the model functions to be incorporated into any other code capable of calling python.Here we will provide a demonstration of VESIcal and its capabilities with applications to various volcanic processes affected by volatiles. VESIcal represents the first tool capable of directly comparing multiple solubility models and equations of state. We find that commonly used models predict surprisingly different volatile solubilities, particularly for pure CO2 or mixed CO2-H2O fluids. Even for melt compositions that are well represented in the calibration datasets of multiple models (e.g., MORBs), calculated solubilities for pure CO2 and pure H2O can deviate from one another by factors of >2 leading to 2x deviations in calculated saturation pressures (e.g., 5 to 10 kbar). The solubility of CO2 predicted by different rhyolitic models also differs substantially, overwhelming other sources of uncertainty such as analytical errors on measurements of volatile contents or uncertainties in crustal density profiles. This highlights the importance of model choice when drawing geological conclusions based on volatiles in magmas.
IRIS 2+ : A Comprehensive Database of Stratified Thermodynamic Models in the Low Solar Atmosphere
Explore the source record for details and available documents.
The role of sea ice dynamics in global climate change
The topics covered include the following: general characteristics of sea ice drift; sea ice rheology; ice thickness distribution; sea ice thermodynamic models; equilibrium thermodynamic models; effect of internal brine pockets and snow cover; model simulations of Arctic Sea ice; and sensitivity of sea ice models to climate change.
Efficient Flowline Simulations of Ice Shelf-Ocean Interactions: Sensitivity Studies with a Fully Coupled Model
Thermodynamic flowline and plume models for the ice shelf-ocean system simplify the ice and ocean dynamics sufficiently to allow extensive exploration of parameters affecting ice-sheet stability while including key physical processes. Comparison between geophysically and laboratory-based treatments of ice-ocean interface thermodynamics shows reasonable agreement between calculated melt rates, except where steep basal slopes and relatively high ocean temperatures are present. Results are especially sensitive to the poorly known drag coefficient, highlighting the need for additional field experiments to constrain its value. These experiments also suggest that if the ice-ocean interface near the grounding line is steeper than some threshold, further steepening of the slope may drive higher entrainment that limits buoyancy, slowing the plume and reducing melting; if confirmed, this will provide a stabilizing feedback on ice sheets under some circumstances.
He I lines in B stars - Comparison of non-local thermodynamic equilibrium models with observations
Profiles of He gamma-gamma 4026, 4387, 4471, 4713, 5876, and 6678 have been obtained in 17 stars of spectral type B0-B5. Parameters of the nonlocal thermodynamic equilibrium models appropriate to each star are determined from the Stromgren index and fits to H-alpha line profiles. These parameters yield generally good fits to the observed He I line profiles, with the best fits being found for the blue He I lines where departures from local thermodynamic equilibrium are relatively small. For the two red lines it is found that, in the early B stars and in stars with log g less than 3.5, both lines are systematically stronger than predicted by the nonlocal thermodynamic equilibrium models.
Thermodynamic and Diffusion Model Estimates on Metamorphic Temperatures and Timescales for Basaltic Eucrite GRA 98098
Introduction: HED meteorites are thought to rep-resent igneous rocks from Vesta’s basaltic crust and preserve evidence of early crustal metamorphism. Determining the temperatures and timescales of thermal metamorphism is important for reconstructing crustal evolution in the early solar system. Here, we study basaltic eucrite Graves Nunataks (GRA) 98098, which has been identified as a highly metamorphosed eucrite [1]. We present new estimates on metamorphic temperatures determined via thermodynamic modeling as well as the initial results from diffusion models constraining timescales of thermal metamorphism. Sample Description: GRA 98098 is an unbrecciated eucrite with a granoblastic plagioclase and pyroxene mineralogy. Millimeter to cm-long lathes of tridymite cross-cut and poikilitically enclose plagioclase and pyroxene [1,this work]. Pyroxene grains have exsolved into Ca-rich (~Wo38En29Fs33) and Ca-poor (~Wo4.5En36Fs59.5) lamellae. Both unzoned and zoned plagioclase grains are observed. Unzoned plagioclase grains are found solely with tridymite laths. These grains have ~An92 compositions. The cores of the zoned plagioclase grains have the same composition and thin, relatively sodic rims (~An67), (Fig. 1). The bulk sample is unusually enriched in highly in-compatible elements and has one of the most fractionated REE patterns reported [1]. Maximum metamorphic temperatures of 985±78°C have been estimated using two-pyroxene thermometry [2]. Methods: Thermodynamic modeling. Thermodynamic models were constructed using the software Perple_X, which employs a Gibbs free energy minimization in order to determine the most stable phase assemblage for a given bulk rock composition [3]. Bulk composition was calculated using mineral com-positions acquired via EMPA (this study) and the observed abundancies present in the thin section. Two bulk compositions were estimated; 1) includes all phases present in the thin section, (assumes that all phases are present during metamorphism), 2) excludes tridymite from the bulk calculation (assumes that tridymite was not present during metamorphism). In order to determine whether metamorphic equilibria was achieved and estimate temperatures of metamorphism, we compared measured pyroxene compositions with thermodynamically predicted compositions [4, Fig. 2]. Diffusion Modeling. Several time-temperature de-pendent diffusion profiles were calculated in order to determine the best match for XAn chemical profiles observed at the edges of the zoned plagioclase (Fig. 3). We assumed that the start condition was a stepwise gradient at the plagioclase/pyroxene interface. We also assumed an average diffusion coefficient (D) and a constant temperature using the equation in [5]. D was determined for two temperatures (T = 1060ºC; near eucrite solidus [6] and T = 985ºC; metamorphism reported in [4]) and then XAn was calculated as a function of distance from plagioclase core to rim using an error function solution to Fick’s second law. Results: Thermodynamic model results are summarized in Fig. 2. For a bulk composition that includes all phases in the thin section, pyroxene endmember compositions plot in the following temperature ranges: Fs ~660-860ºC, En~1000ºC & 1150ºC, and Wo~760-900ºC (Fig. 2a). For a bulk composition that excludes tridymite from the peak metamorphic assemblage (i.e., the bulk composition minus the contribution from tridymite), a temperature range could not be determined for the Fs component of pyroxene. For Wo, T~760-900ºC and En, T~1000ºC & 1150ºC (Fig. 2b). Fig. 3 summarizes the diffusion model results. For T = 1060°C & 985°C, the most appropriate time interval was estimated based on which diffusion curve most matched (solid lines, Fig. 3) the EMPA data. For T = 1060°C, the best looking match was t = 500 ka. For T = 985°C, the best match was t = 7 Ma. Discussion and future work: Temperature estimates from thermodynamic models are not conclusive because the temperature ranges determined for pyroxene endmember stability do not overlap (colored fields in Fig. 2), thus implying that there is disequilibrium between pyroxene crystals and the bulk composition considered [4]. Thus, additional exploration is needed to define a metamorphically equilibrated do-main that accurately records peak temperature. The utility of defining metamorphically equilibrated do-mains to improve the accuracy and level of detail elucidated regarding the petrogenetic history of metamorphose samples has been demonstrated previously [4,7]. We suggest that in the case of Fig. 2a, the thin section composition is not representative of the length scales over which metamorphic equilibrium was achieve and in the case of Fig. 2b, the assumption that tridymite was not present during metamorphism was incorrect. However, results from thermodynamic models can provide insight into the relative timing of mineral and compositional textures. For example from texture alone, it is unclear whether tridymite was igneous in origin and represents the last bits of melt in a crystallizing magma chamber, or if it formed during (and possibly initiated) open system thermal metamorphism. The latter could be consistent with a partial melt hypothesis [8,9] while the former implies that simple fractional crystallization can yield the textures present in GRA 98098. The lack of coincidence be-tween pyroxene endmember compositions in Fig. 2b suggest that the bulk composition minus tridymite was not the assemblage in equilibrium with the pyroxene, suggesting that tridymite was present during metamorphism and formed during igneous crystallization. We conclude that the development of the Na-rich plagioclase rims likely occurred during or immediately after peak thermal metamorphism, because eucrites of similar metamorphic grade and texture have unzoned plagioclase (~An92) [2,4,8], and Na zoning is only observed in the plagioclase not included in the tridymite. This suggests that the zoning formed after tridymite formation, and therefore after igneous crystallization. Thus, the timescales calculated via diffusion modeling possibly represent the time interval over which thermal metamorphism occurred. Cooling rates approximated for the Vestan crust predict that the crust cooled below 300°C around 35-40 Ma after formation[10]. This is consistent with our modeling results that predict formation of the Na rich plagioclase rims occurring at higher temperatures over a period of 0.5 to 7 Ma years. Future work. Additional thermodynamic modeling work will focus on selecting an equilibrated bulk rock domain in which to elucidate metamorphic conditions. Diffusion models currently provide a minimum time-scale, since diffusion slows down as the system cools. Future work with will focus on integrating cooling into the diffusion models and constraining the depth at which thermal metamorphism occurs because it could be used to determine whether the range of time-scales calculated for thermal metamorphism are consistent with the geologic environment.
Towards a generalized computational fluid dynamics technique for all Mach numbers
Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.
Development and Integration of an Advanced Stirling Convertor Linear Alternator Model for a Tool Simulating Convertor Performance and Creating Phasor Diagrams
A simple model of the Advanced Stirling Convertors (ASC) linear alternator and an AC bus controller has been developed and combined with a previously developed thermodynamic model of the convertor for a more complete simulation and analysis of the system performance. The model was developed using Sage, a 1-D thermodynamic modeling program that now includes electro-magnetic components. The convertor, consisting of a free-piston Stirling engine combined with a linear alternator, has sufficiently sinusoidal steady-state behavior to allow for phasor analysis of the forces and voltages acting in the system. A MATLAB graphical user interface (GUI) has been developed to interface with the Sage software for simplified use of the ASC model, calculation of forces, and automated creation of phasor diagrams. The GUI allows the user to vary convertor parameters while fixing different input or output parameters and observe the effect on the phasor diagrams or system performance. The new ASC model and GUI help create a better understanding of the relationship between the electrical component voltages and mechanical forces. This allows better insight into the overall convertor dynamics and performance.
Development and Integration of an Advanced Stirling Convertor Linear Alternator Model for a Tool Simulating Convertor Performance and Creating Phasor Diagrams
A simple model of the Advanced Stirling Convertors (ASC) linear alternator and an AC bus controller has been developed and combined with a previously developed thermodynamic model of the convertor for a more complete simulation and analysis of the system performance. The model was developed using Sage, a 1-D thermodynamic modeling program that now includes electro-magnetic components. The convertor, consisting of a free-piston Stirling engine combined with a linear alternator, has sufficiently sinusoidal steady-state behavior to allow for phasor analysis of the forces and voltages acting in the system. A MATLAB graphical user interface (GUI) has been developed to interface with the Sage software for simplified use of the ASC model, calculation of forces, and automated creation of phasor diagrams. The GUI allows the user to vary convertor parameters while fixing different input or output parameters and observe the effect on the phasor diagrams or system performance. The new ASC model and GUI help create a better understanding of the relationship between the electrical component voltages and mechanical forces. This allows better insight into the overall convertor dynamics and performance.
Partitioning of F and Cl Between Apatite and a Synthetic Shergottite Liquid (QUE 94201) at 4 Gpa from 1300 TO 1500 C
Apatite [Ca5(PO4)3(F,Cl,OH)] is present in a wide range of planetary materials. Due to the presence of volatiles within its crystal structure (Xsite), many recent studies have attempted to use apatite to constrain the volatile contents of planetary magmas and mantle sources. In order to use the volatile contents of apatite to accurately determine the abundances of volatiles in coexisting silicate melt or fluids, thermodynamic models for the apatite solid solution and for the apatite components in multicomponent silicate melts and fluids are required. Although some thermodynamic models for apatite have been developed, they are incomplete. Furthermore, no mixing model is available for all of the apatite components in silicate melts or fluids, especially for F and Cl components. Several experimental studies have investigated the apatite-melt and apatite-fluid partitioning behavior of F, Cl, and OH in terrestrial and planetary systems, which have determined that apatite-melt partitioning of volatiles are best described as exchange equilibria similar to Fe-Mg partitioning between olivine and silicate melt. However, McCubbin et al., recently reported that the exchange coefficients vary in portions of apatite compositional space where F, Cl, and OH do not mix ideally in apatite. In particular, solution calorimetry data of apatite compositions along the F-Cl join exhibit substantial excess enthalpies of mixing, and McCubbin et al. reported substantial deviations in the Cl-F exchange Kd along the F-Cl apatite join that could be explained by the preferential incorporation of F into apatite. In the present study, we assess the effect of apatite crystal chemistry on F-Cl exchange equilibria between apatite and melt at 4 GPa over the temperature range of 1300-1500 C. The goal of these experiments is to assess the variation in the Ap-melt Cl-F exchange Kd over a broad range of F:Cl ratios in apatite. The results of these experiments could be used to understand at what composition apatite shifts from a hexagonal unit cell with space group P63/m to a unit cell with monoclinic symmetry within space group P21/b. We anticipate that this transition occurs at >70% chlorapatite based on solution calorimetry data.
Apatite-Melt Partitioning at 1 Bar: An Assessment of Apatite-Melt Exchange Equilibria Resulting from Non-Ideal Mixing of F and Cl in Apatite
The mineral apatite [Ca5(PO4)3(F,Cl,OH)] is present in a wide range of planetary materials. Due to the presence of volatiles within its crystal structure (X-site), many recent studies have attempted to use apatite to constrain the volatile contents of planetary magmas and mantle sources. In order to use the volatile contents of apatite to precisely determine the abundances of volatiles in coexisting silicate melt or fluids, thermodynamic models for the apatite solid solution and for the apatite components in multi-component silicate melts and fluids are required. Although some thermodynamic models for apatite have been developed, they are incomplete. Furthermore, no mixing model is available for all of the apatite components in silicate melts or fluids, especially for F and Cl components. Several experimental studies have investigated the apatite-melt and apatite-fluid partitioning behavior of F, Cl, and OH in terrestrial and planetary systems, which have determined that apatite-melt partitioning of volatiles are best described as exchange equilibria similar to Fe-Mg partitioning between olivine and silicate melt. However, McCubbin et al. recently reported that the exchange coefficients may vary in portions of apatite compositional space where F, Cl, and OH do not mix ideally in apatite. In particular, solution calorimetry data of apatite compositions along the F-Cl join exhibit substantial excess enthalpies of mixing. In the present study, we conducted apatite-melt partitioning experiments in evacuated, sealed silica-glass tubes at approximately 1 bar and 950-1050 degrees Centigrade on a synthetic Martian basalt composition equivalent to the basaltic shergottite Queen Alexandria Range (QUE) 94201. These experiments were conducted dry, at low pressure, to assess the effects of temperature and apatite composition on the partitioning behavior of F and Cl between apatite and basaltic melt along the F-Cl apatite binary join, where there is non-ideal mixing of F and Cl in apatite.
Designing Molten Salt Eutectics
Designing stable electrolytes with target properties is an important challenge in realizing next generation energy storage devices. Molten salt eutectics-based electrolytes are known for their stability with minimal parasitic reactions when compared to traditional organic electrolytes and are an attractive option for different battery chemistries. The operating temperature of the molten salt batteries depends on the melting temperature of the eutectic and hence there is a necessity to discover novel low melting temperature molten salt eutectic mixtures for energy storage applications. In this work we develop a high throughput computational screening approach for molten salt mixtures using thermodynamic modeling and machine learning (ML). COSMO-SAC model and ML approaches were independently developed based on the existing experimental data and these models were further used to predict the eutectic melting temperature and composition of several new binary, ternary, and quaternary mixtures. We show that combining ML and thermodynamic modeling strategies is effective in exploring the vast design space of molten salt mixtures.
I-Xe dating of Allende inclusions - Antiquity and fine structure
Two inclusions from the Allende meteorite exhibit multiple I-Xe isochrons corresponding to an age difference of 3.7 + or - 1.5 m.y. for their constituent minerals. Consequently, I-Xe dating can establish, independently of thermodynamic models, a partial sequence of formation for minerals that incorporated radioactive I-1290. Additionally, the results give an estimate of the cooling rate for the solar nebula when validity of thermodynamic models is assumed.
Propulsion System Simulation Using the Toolbox for the Modeling and Analysis of Thermodynamic System T-MATS
A simulation toolbox has been developed for the creation of both steady-state and dynamic thermodynamic software models. This paper describes the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS), which combines generic thermodynamic and controls modeling libraries with a numerical iterative solver to create a framework for the development of thermodynamic system simulations, such as gas turbine engines. The objective of this paper is to present an overview of T-MATS, the theory used in the creation of the module sets, and a possible propulsion simulation architecture. A model comparison was conducted by matching steady-state performance results from a T-MATS developed gas turbine simulation to a well-documented steady-state simulation. Transient modeling capabilities are then demonstrated when the steady-state T-MATS model is updated to run dynamically.
Propulsion System Simulation Using the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS)
A simulation toolbox has been developed for the creation of both steady-state and dynamic thermodynamic software models. This paper describes the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS), which combines generic thermodynamic and controls modeling libraries with a numerical iterative solver to create a framework for the development of thermodynamic system simulations, such as gas turbine engines. The objective of this paper is to present an overview of T-MATS, the theory used in the creation of the module sets, and a possible propulsion simulation architecture. A model comparison was conducted by matching steady-state performance results from a T-MATS developed gas turbine simulation to a well-documented steady-state simulation. Transient modeling capabilities are then demonstrated when the steady-state T-MATS model is updated to run dynamically.