Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic assessment”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Estimation of Thermodynamic and Dynamic Contributions to Sea Ice Growth in the Central Arctic Using ICESat-2 and MOSAiC SIMBA Buoy Data

The fine spatial resolution of the ICESat-2 (IS2) satellite altimeter allows monitoring the evolution of sea ice thickness with detailed dynamic information (e.g. ridges and leads). In this study, we first assess the ability of IS2 to estimate thermodynamic ice growth and dynamic thickening during the ice-growing season in the central Arctic Ocean. As an indicator of the thermodynamic ice growth, we use 10 thermistor string-based sea ice mass balance array (SIMBA) buoys deployed at a scale of ~ 50 km from the Icebreaker Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. We collect IS2 data within 20 km buffer distance from the individual buoys, and calculate the mode, median, and mean of the IS2-derived ice thickness. The IS2 modal thickness shows the least bias (−0.169 m) with the buoy ice thickness, representing level ice thickness. In addition, the increasing rate of the IS2 modal thickness is close to the thermodynamic ice growth with a small bias of −0.054 cm/day. However, the increasing rates of the IS2 median and mean thickness are greater than the thermodynamic ice growth by about 0.114 cm/day and 0.198 cm/day, respectively, because they also include ice growth caused by thickness redistribution during dynamic deformation. The dynamic contributions may account for 26.1 ± 10.3% and 34.4 ± 10.1% of the total increase of the IS2 median and mean thickness, respectively. Within a ~ 50 km radius area from the MOSAiC Central Observatory, IS2 measurements exhibit that the ridge fraction increased from <2% in November to ~ 4% in March (~0.029%/day of average increasing rate) and ridge height increased about 0.047 cm/day during the same period. However, lead formation does not show significant contributions to the dynamic ice thickening because leads are temporary features lasting only 2–3 days. Although there are considerable uncertainties in IS2 ice thickness estimation and IS2-buoy thickness comparison, this study emphasizes the importance of combining measurements by IS2 and SIMBA buoys to explain the regional sea ice mass balance with separating the thermodynamic and dynamic contributions.

satellite altimeter↗
DOE OSTI2020

Three-dimensional modeling and simulations of single-crystal and bi-crystal titanium for high-strain-rate loading conditions

High purity single crystal titanium (Ti) under shock wave loading is modeled under both one-dimensional and three-dimensional cylindrical conditions. Cylinder sizes of 10 μm and 20 μm radius are both considered in order to assess influence of boundary conditions. A thermodynamically consistent single crystal model for application to shock conditions is presented. The model accounts for the coupled non-linear elastic, dislocation slip, deformation twinning, and structural phase transformation response of the titanium material. Plate impact experiment results using a copper flyer are used to compare against the simulations for crystals oriented in [0001] and [10$\overline1$1] crystallographic directions. The one-dimensional and threedimensional simulations of the two differently oriented single crystals indicate differences between the one-dimensional and three-dimensional representation, especially for the [10$\overline1$1] oriented single crystal. This orientation breaks the relative orientation symmetry between the crystal and cylinder which otherwise exists for the [0001] oriented single crystal. A significant amount of heterogeneity in the field response of the[10$\overline1$1] oriented simulation was demonstrated due to the highly coupled nature of the deformation. A bi-crystal model composed of both the [0001] and [10$\overline1$1] orientations with the boundary between the two along the axis of the cylinder is also considered for a cylinder model size of 10 μm. The results indicate a strong interaction between the two grains that affects the ω phase volume fraction achieved relative to the single crystal calculations.

36 MATERIALS SCIENCE↗
NASA NTRS1969

Thermophysical properties of sodium

Assessment is given of physical and thermodynamic properties of sodium. FORTRAN subroutine computes enthalphy and entropy of sodium in given state, and composition, molecular weight, volume, and compressibility factor of corresponding vapor. Tabular results for saturated liquid and vapor are presented for a 500-2500 degree F range.

Golden, G. H.↗
NASA NTRS1983

The coal-fired gas turbine locomotive - A new look

Advances in turbomachine technology and novel methods of coal combustion may have made possible the development of a competitive coal fired gas turbine locomotive engine. Of the combustor, thermodynamic cycle, and turbine combinations presently assessed, an external combustion closed cycle regenerative gas turbine with a fluidized bed coal combustor is judged to be the best suited for locomotive requirements. Some merit is also discerned in external combustion open cycle regenerative systems and internal combustion open cycle regenerative gas turbine systems employing a coal gasifier. The choice of an open or closed cycle depends on the selection of a working fluid and the relative advantages of loop pressurization, with air being the most attractive closed cycle working fluid on the basis of cost.

Liddle, S. G.↗
NASA NTRS2012

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.↗
NASA NTRS1986

Shock-tube pyrolysis of acetylene - Sensitivity analysis of the reaction mechanism for soot formation

The impact of thermodynamic parameters on the sensitivity of model predictions of soot formation by shock-tube pyrolysis of acetylene were assessed analytically. The pyrolysis process was treated as having three components: initiation, the initial pyrolysis stages; cyclization, formation of larger molecules and radicals and small aromatic molecules; and polymerization, further growth of aromatic rings. Rate equations are reviewed for each component. Thermodynamic effects were assessed by varying the C2H-H and C2H3-H bond energies and the Ct-(Ct) group additivity value. Any change in the C2H-H bond energy had a significant impact on the temperature and the maximum amount of the soot yield. The findings underscore the necessity of using accurate thermodynamic data for modeling high-temperature chemical kinetics.

Frenklach, M.↗
DOE OSTI2024

Thermal Property Modeling and Assessment of the Physical Properties of FLiNaK

Here, the thermodynamic properties of the LiF–NaF–KF system and its subsystems were re-evaluated using the CALPHAD method, focusing on the deviations from ideality in the heat capacity [C p (T)] of mixtures, temperature-dependence of mixing enthalpies (Δ mix H), vapor pressures of mixtures, and enthalpies of fusion (ΔH fus ) using reported values and measurements. The results of this work confirmed a pseudoternary eutectic at 732.9 K and 46.5LiF–11.5NaF–42KF mol % (FLiNaK) with determined Cp(T) = 69.822 + 0.000679T + 2,137,152T –2 J mol –1 K –1 (732.9 < T < 1200 K) and ΔH fus = 18,051 J mol –1 in good agreement to the experimental values. The resulting accurate thermodynamic representation can be used to find relevant excess thermochemical properties in the whole pseudoternary space and provide insights into their composition and temperature dependence. A detailed evaluation of select thermophysical properties was also conducted to determine recommended density, thermal diffusivity, thermal conductivity, vapor pressures, and boiling point with calculated uncertainties over the temperature range of interest for FLiNaK–salt energy applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2003

Internal Flow Thermal/Fluid Modeling of STS-107 Port Wing in Support of the Columbia Accident Investigation Board

As part of the aero-thermodynamics team supporting the Columbia Accident Investigation Board (CAB), the Marshall Space Flight Center was asked to perform engineering analyses of internal flows in the port wing. The aero-thermodynamics team was split into internal flow and external flow teams with the support being divided between shorter timeframe engineering methods and more complex computational fluid dynamics. In order to gain a rough order of magnitude type of knowledge of the internal flow in the port wing for various breach locations and sizes (as theorized by the CAB to have caused the Columbia re-entry failure), a bulk venting model was required to input boundary flow rates and pressures to the computational fluid dynamics (CFD) analyses. This paper summarizes the modeling that was done by MSFC in Thermal Desktop. A venting model of the entire Orbiter was constructed in FloCAD based on Rockwell International s flight substantiation analyses and the STS-107 reentry trajectory. Chemical equilibrium air thermodynamic properties were generated for SINDA/FLUINT s fluid property routines from a code provided by Langley Research Center. In parallel, a simplified thermal mathematical model of the port wing, including the Thermal Protection System (TPS), was based on more detailed Shuttle re-entry modeling previously done by the Dryden Flight Research Center. Once the venting model was coupled with the thermal model of the wing structure with chemical equilibrium air properties, various breach scenarios were assessed in support of the aero-thermodynamics team. The construction of the coupled model and results are presented herein.

Sharp, John R.↗
DOE OSTI2026

Report priority gaps in high temperature thermodynamic data (Interim Progress Report)

This interim progress report (Level 4 Milestone Number M4SF-26LL010203023) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite Host Rock Properties & Processes SF-26LL01020302. Our focus is to assess gaps in data availability and understanding for radionuclide thermodynamics within the context of a “hot repository” concept and expand SUPCRT-NE database development to address higher temperatures needed for a DPC DGR disposal concept. The database is intended to inform the argillite GDSA baseline model. The leading European thermochemical database (Thermochimie) is only applicable to temperatures below 80°C. Thus, a US effort to integrate and expand upon other international thermodynamics database efforts is needed, particular if a “hot repository” concept moves forward.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

The solubility of La hydroxide and stability of La 3+ and La hydroxyl complexes at acidic to mildly acidic pH from 25 to 250 °C

The mobility of rare earth elements (REE) in natural hydrothermal systems can be assessed using geochemical modeling, which requires reliable thermodynamic data of relevant aqueous species. In this study, we evaluate the controls of pH and temperature on La speciation and the role of hydroxyl complexes in REE transport at hydrothermal conditions. Batch-type hydrothermal solubility experiments were conducted using synthetic La hydroxide powders equilibrated in perchloric acid-based aqueous solutions at temperatures between 150 and 250 °C and starting pH of 2 to 5. The La hydroxide solubility is retrograde with temperature and displays a strong pH dependence with a decrease in La concentrations from acidic to mildly acidic pH spanning between 3 and 5 orders of magnitude (e.g. log La molality of –2.5 to –7.2 at 250 °C). Thermodynamic optimizations using GEMSFITS allow to retrieve the standard partial molal Gibbs energies for the La 3+ aqua ion and the formation constants for the La hydroxyl species (i.e., LaOH 2+ , La(OH) 2 + , La(OH) 3 0 ) between 25 and 250 °C. A comparison between the experimentally derived thermodynamic properties with the calculated values from the Helgeson-Kirkham-Flowers equation of state parameters indicates an increased divergence with temperature. Discrepancies in standard partial molal Gibbs energies range between ~ 1 – 12 kJ/mol and result in a predicted La hydroxide solubility differing by up to 3 orders of magnitude at 250 °C. Speciation calculations indicate a higher stability of La 3+ and LaOH 2+ over the other La hydroxyl species in the studied pH range of 3.4 to 6. Here, the optimized thermodynamic properties for La aqueous species have important implications for modeling the solubility of REE minerals such as monazite and the mobility of REE in hydrothermal systems.

58 GEOSCIENCES↗
DOE OSTI2022

Evaluating Tropical Precipitation Relations in CMIP6 Models with ARM Data

Abstract A set of diagnostics based on simple, statistical relationships between precipitation and the thermodynamic environment in observations is implemented to assess phase 6 of the Coupled Model Intercomparison Project (CMIP6) model behavior with respect to precipitation. Observational data from the Atmospheric Radiation Measurement (ARM) permanent field observational sites are augmented with satellite observations of precipitation and temperature as an observational baseline. A robust relationship across observational datasets between column water vapor (CWV) and precipitation, in which conditionally averaged precipitation exhibits a sharp pickup at some critical CWV value, provides a useful convective onset diagnostic for climate model comparison. While a few models reproduce an appropriate precipitation pickup, most models begin their pickup at too low CWV and the increase in precipitation with increasing CWV is too weak. Convective transition statistics compiled in column relative humidity (CRH) partially compensate for model temperature biases—although imperfectly since the temperature dependence is more complex than that of column saturation. Significant errors remain in individual models and weak pickups are generally not improved. The conditional-average precipitation as a function of CRH can be decomposed into the product of the probability of raining and mean precipitation during raining times (conditional intensity). The pickup behavior is primarily dependent on the probability of raining near the transition and on the conditional intensity at higher CRH. Most models roughly capture the CRH dependence of these two factors. However, compensating biases often occur: model conditional intensity that is too low at a given CRH is compensated in part by excessive probability of precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗
NASA NTRS1982

Turbofan forced mixer-nozzle internal flowfield. Volume 1: A benchmark experimental study

An experimental investigation of the flow field within a model turbofan forced mixer nozzle is described. Velocity and thermodynamic state variable data for use in assessing the accuracy and assisting the further development of computational procedures for predicting the flow field within mixer nozzles are provided. Velocity and temperature data suggested that the nozzle mixing process was dominated by circulations (secondary flows) of a length scale on the order the lobe dimensions which were associated with strong radial velocities observed near the lobe exit plane. The 'benchmark' model mixer experiment conducted for code assessment purposes is discussed.

Paterson, R. W.↗
NASA NTRS1991

Thermal diffusion in metal-sulfide liquids - Early results

Experiments were carried out to evaluate the Soret effect in liquid Fe-Ni-S-P alloys in order to gain a better understanding of the physical and thermodynamic properties of metallic liquids and to assess the possibility of systematic errors in Czochralski growth techniques. The metal-sulfide liquids were found to show a substantial Soret effect and, contrary to previous expectation, can be as large as that seen in the silicate system. The segregation is largely produced by S-Fe interactions. The P redistribution occurs to reduce activity gradients in P consequent upon S gradients, and P segregation can be approximately predicted in Fe-Ni-S liquids using the activity model of Jones and Malvin (1990). It is inferred that the sense of the Soret segregation, with P going to the cold end and S going to the hot end, is in accordance with the prediction of Malvin et al. (1986), who speculated that the crystal-pulling experiments of Sellamuthu and Goldstein (1985) was influenced by the Soret effect.

Jones, J. H.↗
NASA NTRS2006

Changes in Tropical Cyclone Intensity Over the Past 30 Years: A Global and Dynamic Perspective

The hurricane season of 2005 was the busiest on record and Hurricane Katrina (2005) is believed to be the costliest hurricane in U. S. history. There are growing concerns regarding whether this increased tropical cyclone activity is a result of global warming, as suggested by Emanuel(2005) and Webster et al. (2005), or just a natural oscillation (Goldenberg et al. 2001). This study examines the changes in tropical cyclone intensity to see what were really responsible for the changes in tropical cyclone activity over the past 30 years. Since the tropical sea surface temperature (SST) warming also leads to the response of atmospheric circulation, which is not solely determined by the local SST warming, this study suggests that it is better to take the tropical cyclone activities in the North Atlantic (NA), western North Pacific (WNP) and eastern North Pacific (ENP) basins as a whole when searching for the influence of the global-scale SST warming on tropical cyclone intensity. Over the past 30 years, as the tropical SST increased by about 0.5 C, the linear trends indicate 6%, 16% and 15% increases in the overall average intensity and lifetime and the annual frequency. Our analysis shows that the increased annual destructiveness of tropical cyclones reported by Emanuel(2005) resulted mainly from the increases in the average lifetime and annual frequency in the NA basin and from the increases in the average intensity and lifetime in the WNP basin, while the annual destructiveness in the ENP basin generally decreased over the past 30 years. The changes in the proportion of intense tropical cyclones reported by Webster et a1 (2005) were due mainly to the fact that increasing tropical cyclones took the tracks that favor for the development of intense tropical cyclones in the NA and WNP basins over the past 30 years. The dynamic influence associated with the tropical SST warming can lead to the impact of global warming on tropical cyclone intensity that may be very different from our current assessments, which were mainly based on the thermodynamic theory of tropical cyclone intensity.

Wu, Liguang↗
NASA NTRS2011

Inflow/Outflow Boundary Conditions with Application to FUN3D

Several boundary conditions that allow subsonic and supersonic flow into and out of the computational domain are discussed. These boundary conditions are demonstrated in the FUN3D computational fluid dynamics (CFD) code which solves the three-dimensional Navier-Stokes equations on unstructured computational meshes. The boundary conditions are enforced through determination of the flux contribution at the boundary to the solution residual. The boundary conditions are implemented in an implicit form where the Jacobian contribution of the boundary condition is included and is exact. All of the flows are governed by the calorically perfect gas thermodynamic equations. Three problems are used to assess these boundary conditions. Solution residual convergence to machine zero precision occurred for all cases. The converged solution boundary state is compared with the requested boundary state for several levels of mesh densities. The boundary values converged to the requested boundary condition with approximately second-order accuracy for all of the cases.

Carlson, Jan-Renee↗
DOE OSTI2022

Structure and thermodynamics of calcium rare earth silicate oxyapatites, Ca 2 RE 8 (SiO 4 ) 6 O 2 (RE = Pr, Tb, Ho, Tm)

Calcium rare earth silicate oxyapatites, (Ca 2 RE 8 (SiO 4 ) 6 O 2 ), are of interest as components of glass-ceramic nuclear waste forms. To assess their long-term behavior in a geologic repository, it is essential to determine their structure and thermodynamic stability at relevant conditions. Here, in this work, we performed detailed structural and thermodynamic investigations on Ca 2 Pr 8 (SiO 4 ) 6 O 2 , Ca 2 Tb 8 (SiO 4 ) 6 O 2 , Ca 2 Ho 8 (SiO 4 ) 6 O 2 , and Ca 2 Tm 8 (SiO 4 ) 6 O 2 by high energy synchrotron powder X-ray diffraction combined with Rietveld analysis and high temperature oxide melt drop solution calorimetry. Enthalpies of formation from constituent oxides (ΔH f,ox ) were determined to be -765.1 ± 22.8 kJ/mol for Ca 2 Pr 8 (SiO 4 ) 6 O 2 ; -638.9 ± 20.5 kJ/mol for Ca 2 Tb 8 (SiO 4 ) 6 O 2 ; -643.3 ± 10.3 kJ/mol for Ca 2 Ho 8 (SiO 4 ) 6 O 2 ; and -403.2 ± 5.1 kJ/mol for Ca 2 Tm 8 (SiO 4 ) 6 O 2 . These thermodynamic parameters were used in assessing the thermochemical stability of these phases in the presence of water vapor from room temperature to 600 K, as encountered in the subsurface environments of a geological repository.

36 MATERIALS SCIENCE↗