Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic consistency”

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

Metastable sound speed in gas-liquid mixtures

A new method of calculating speed of sound for two-phase flow is presented. The new equation assumes no phase change during the propagation of an acoustic disturbance and assumes that only the total entropy of the mixture remains constant during the process. The new equation predicts single-phase values for the speed of sound in the limit of all gas or all liquid and agrees with available two-phase, air-water sound speed data. Other expressions used in the two-phase flow literature for calculating two-phase, metastable sound speed are reviewed and discussed. Comparisons are made between the new expression and several of the previous expressions -- most notably a triply isentropic equation as used, a triply isentropic equation as used, among others, by Karplus and by Wallis. Appropriate differences are pointed out and a thermodynamic criterion is derived which must be satisfied in order for the triply isentropic expression to be thermodynamically consistent. This criterion is not satisfied for the cases examined, which included two-phase nitrogen, air-water, two-phase parahydrogen, and steam-water. Consequently, the new equation derived is found to be superior to the other equations reviewed.

Bursik, J. W.↗

The major volume /density/ of solid oxygen in equilibrium with vapor

Data from the literature on the molar volume of solid oxygen have been compiled and critically analyzed. A correlated and thermodynamically consistent set of molar volumes, including the volume changes at the various solid phase transitions, is presented. Evidence for the existence of a delta-solid phase is reviewed. Uncertainties in the data and in the recommended set of values are discussed.

Roder, H. M.↗

Morphological instability in rapid directional solidification

Mullins and Sekerka (1964) showed for fixed temperature gradient that the planar interface is linearly stable for all pulling speeds V above some critical value, the absolute stability limit. Near this limit, where solidification rates are rapid, the assumption of local equilibrium at the interface may be violated. Here, nonequilibrium effects are incorporated into a linear stability analysis of the planar front by allowing the segregation coefficient and interface temperature to depend on V in a thermodynamically consistent way. The absolute stability limit of the cellular mode is modified. A new oscillatory state is formed which, in the absence of latent heat, has a critical wavenumber of zero; by itself this instability would lead to the formation of solute bands in the solid. This mode has its own absolute-stability limit determined by solute trapping and kinetics. Under certain conditions, there exists a window of stability above the steady absolute-stability boundary and below the oscillatory-stability boundary; here the planar segregation-free state is restabilized.

Merchant, G. J.↗

Phase equilibrium of methane and nitrogen at low temperatures - Application to Titan

Since the vapor phase composition of Titan's methane-nitrogen lower atmosphere is uniquely determined as a function of the Gibbs phase rule, these data are presently computed via integration of the Gibbs-Duhem equation. The thermodynamic consistency of published measurements and calculations of the vapor phase composition is then examined, and the saturated mole fraction of gaseous methane is computed as a function of altitude up to the 700-mbar level. The mole fraction is found to lie approximately halfway between that computed from Raoult's law, for a gas in equilibrium with an ideal solution of liquid nitrogen and methane, and that for a gas in equilibrium with pure liquid methane.

Kouvaris, Louis C.↗

A Unified Theory of Non-Ideal Gas Lattice Boltzmann Models

A non-ideal gas lattice Boltzmann model is directly derived, in an a priori fashion, from the Enskog equation for dense gases. The model is rigorously obtained by a systematic procedure to discretize the Enskog equation (in the presence of an external force) in both phase space and time. The lattice Boltzmann model derived here is thermodynamically consistent and is free of the defects which exist in previous lattice Boltzmann models for non-ideal gases. The existing lattice Boltzmann models for non-ideal gases are analyzed and compared with the model derived here.

Luo, Li-Shi↗

Thin Interface Asymptotics for an Energy/Entropy Approach to Phase-Field Models with Unequal Conductivities

Karma and Rapped recently developed a new sharp interface asymptotic analysis of the phase-field equations that is especially appropriate for modeling dendritic growth at low undercoolings. Their approach relieves a stringent restriction on the interface thickness that applies in the conventional asymptotic analysis, and has the added advantage that interfacial kinetic effects can also be eliminated. However, their analysis focussed on the case of equal thermal conductivities in the solid and liquid phases; when applied to a standard phase-field model with unequal conductivities, anomalous terms arise in the limiting forms of the boundary conditions for the interfacial temperature that are not present in conventional sharp-interface solidification models, as discussed further by Almgren. In this paper we apply their asymptotic methodology to a generalized phase-field model which is derived using a thermodynamically consistent approach that is based on independent entropy and internal energy gradient functionals that include double wells in both the entropy and internal energy densities. The additional degrees of freedom associated with the generalized phased-field equations can be chosen to eliminate the anomalous terms that arise for unequal conductivities.

McFadden, G. B.↗

Simulation of Delamination Under High Cycle Fatigue in Composite Materials Using Cohesive Models

A new thermodynamically consistent damage model is proposed for the simulation of high-cycle fatigue crack growth. The basis for the formulation is an interfacial degradation law that links Fracture Mechanics and Damage Mechanics to relate the evolution of the damage variable, d, with the crack growth rate da/dN. The damage state is a function of the loading conditions (R and (Delta)G) as well as the experimentally-determined crack growth rates for the material. The formulation ensures that the experimental results can be reproduced by the analysis without the need of additional adjustment parameters.

Camanho, Pedro P.↗

Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates

Progressive damage and failure analysis (PDFA) tools are needed to predict the nonlinear response of advanced fiber-reinforced composite structures. Predictive tools should incorporate the underlying physics of the damage and failure mechanisms observed in the composite, and should utilize as few input parameters as possible. The purpose of the Enhanced Schapery Theory (EST) was to create a PDFA tool that operates in conjunction with a commercially available finite element (FE) code (Abaqus). The tool captures the physics of the damage and failure mechanisms that result in the nonlinear behavior of the material, and the failure methodology employed yields numerical results that are relatively insensitive to changes in the FE mesh. The EST code is written in Fortran and compiled into a static library that is linked to Abaqus. A Fortran Abaqus UMAT material subroutine is used to facilitate the communication between Abaqus and EST. A clear distinction between damage and failure is imposed. Damage mechanisms result in pre-peak nonlinearity in the stress strain curve. Four internal state variables (ISVs) are utilized to control the damage and failure degradation. All damage is said to result from matrix microdamage, and a single ISV marks the micro-damage evolution as it is used to degrade the transverse and shear moduli of the lamina using a set of experimentally obtainable matrix microdamage functions. Three separate failure ISVs are used to incorporate failure due to fiber breakage, mode I matrix cracking, and mode II matrix cracking. Failure initiation is determined using a failure criterion, and the evolution of these ISVs is controlled by a set of traction-separation laws. The traction separation laws are postulated such that the area under the curves is equal to the fracture toughness of the material associated with the corresponding failure mechanism. A characteristic finite element length is used to transform the traction-separation laws into stress-strain laws. The ISV evolution equations are derived in a thermodynamically consistent manner by invoking the stationary principle on the total work of the system with respect to each ISV. A novel feature is the inclusion of both pre-peak damage and appropriately scaled, post-peak strain softening failure. Also, the characteristic elements used in the failure degradation scheme are calculated using the element nodal coordinates, rather than simply the square root of the area of the element.

Pineda, Evan J.↗

Analysis of the Effect of Water Activity on Ice Formation Using a New Theory of Nucleation

In this work a new theory of nucleation is developed and used to investigate the effect of water activity on the formation of ice within super-cooled droplets. The new theory is based on a novel concept where the interface is assumed to be made of liquid molecules trapped by the solid matrix. Using this concept new expressions are developed for the critical ice germ size and the nucleation work, with explicit dependencies on temperature and water activity. However unlike previous approaches, the new theory does not depend on the interfacial tension between liquid and ice. Comparison against experimental results shows that the new theory is able to reproduce the observed effect of water activity on nucleation rate and freezing temperature. It allows for the first time a theoretical derivation of the constant shift in water activity between melting and nucleation. The new theory offers a consistent thermodynamic view of ice nucleation, simple enough to be applied in atmospheric models of cloud formation.

Water Activity↗

Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation

A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments relies upon the global momentum conservation of the fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. Numerical examples illustrate the method's application to predicting bulk fluid motion including lateral propellant slosh in low-g conditions.

Orr, Jeb S.↗

Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation

A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments is based upon a novel approach that relies on the global momentum conservation of the closed fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. A numerical example illustrates the method's application to prediction of bulk fluid behavior during a spacecraft ullage settling maneuver.

Orr, Jeb S.↗

A Three-Dimensional Thermodynamically Based Function for the Progressive Failure of Unidirectional Composites

A thermodynamically-based work potential theory for modeling progressive damage for laminated, unidirectional composites assuming plane stress (2D Schaperystheory) is extended to three dimensioanl (3D). An internal state variable, S, is defined to account for the dissipated energy due to damage evolution in the form of microstructure changes in the matrix. With the stationary of the total work potential with respect to the internal state variable, a thermodynamically-consistent set of evolution equationsis derived. The internal state variable is related to the transverse and shear modulithrough microdamage functions. In the first part of this work, coupon specimensare prepared to conduct experiments to characterize the relations between the internalstate variable and the transverse modulus as well as shear modulus. The informationis subsequent used for the prediction of three point bending test. In the second partof this work, objectivity is studied. Three separate methods utilizing different definitions of a reduced internal state variable or of the order of the polynomials are used to represent the matrix micro-damage functions are employed. The three methods are implemented in a user defined subroutine within a commercial finite element method software package. Results from numerical simulations of a center-notched composites panel are compared. The agreement in the maximum stress predictions among the three methods indicates that objectivity, with respect to the functional form of themicrodamage functions, is satisfied.

carbon fiber reinforced polymer composites↗

Optimization of a 10–Step Reduced CH₄/O₂ Combustion Mechanism for RDRE Applications

A 10-step, 9-species reduced mechanism for methane/oxygen combustion is evaluated and optimized for use in rotating detonation rocket engine (RDRE) simulations at an operating pressure of 10 bar. The mechanism is assessed against GRI-Mech 3.0 for ignition delay times, adiabatic equilibrium temperature and species, Chapman-Jouguet (CJ) detonation parameters, and laminar flame speed. A sequential Nelder–Mead optimization procedure is applied to the rate constants of reactions R1, R4, and R5, targeting ignition delay agreement over the equivalence ratio range φ = 1.0–1.4 and temperature range T = 1400–2000 K. The optimized mechanism achieves ignition delay agreement within 28% of GRI-Mech 3.0 for the temperature range relevant to RDRE operation (τ ≥ 5 μs). Adiabatic equilibrium temperatures and CJ detonation parameters are reproduced within 0.15% of GRI-Mech 3.0 values, confirming thermodynamic consistency of the mechanism. Laminar flame speed predictions deviate significantly from GRI-Mech 3.0, which is attributed to the absence of HO₂ chemistry in the reduced mechanism; this limitation does not affect RDRE simulation fidelity since laminar flame propagation is of secondary importance in RDREs.

Computational Fluid Dynamics↗

Learning heterogeneous reaction kinetics from X-ray movies pixel-by-pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries 1 and electrocatalysts 2. Experimental characterizations of such materials by operando microscopy produce rich image datasets 3–6, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation 7. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

Chueh, William↗

Organic carbon oxidation state shapes fermentative methanogenic microbiomes and controls greenhouse gas fluxes

Organic compounds with a negative nominal oxidation state of carbon (NOSC) are thermodynamically recalcitrant in anaerobic ecosystems, but few studies have measured the influence of NOSC on carbon degradation rates, gaseous product yields, or microbiome composition. We amended anaerobic rice paddy sediment microcosms with water-soluble monomeric organic carbon compounds varying in NOSC. Consistent with thermodynamic and stoichiometric predictions, negative NOSC compounds are catabolized more slowly but produce more methane per mole of carbon. Negative NOSC microbiomes have higher alpha diversity, more syntrophs and methanogens, and fewer fermentative bacteria. Strikingly, fermentative bacterial taxa display genomically encoded NOSC catabolic preferences both in the lab and field. Negative NOSC-preferring fermenters have longer predicted doubling times, consistent with the thermodynamic recalcitrance of their preferred substrates. We propose that microbial NOSC catabolic preferences may reflect the thermodynamic niche of microorganisms and we anticipate that extending research on microbial catabolic preferences to a greater variety of organic carbon substrates and diverse microbiomes will improve our understanding of microbial carbon cycling and trait evolution.

Hu, Ruiwen↗

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Characteristics of Subcooled Liquid Methane During Passage Through a Spray-Bar Joule-Thompson Thermodynamic Vent System

NASA s Marshall Space Flight Center (MSFC) conducted liquid methane (LCH4) testing in November 2006 using the multipurpose hydrogen test bed (MHTB) outfitted with a spray-bar thermodynamic vent system (TVS). The basic objective was to identify any unusual or unique thermodynamic characteristics associated with subcooled LCH4 that should be considered in the design of space-based TVSs. Thirteen days of testing were performed with total tank heat loads ranging from 720 W to 420 W at a fill level of approximately 90%. During an updated evaluation of the data, it was noted that as the fluid passed through the Joule Thompson expansion, thermodynamic conditions consistent with the pervasive presence of metastability were indicated. This paper describes the observed thermodynamic conditions that correspond with metastability and effects on TVS performance.

Hastings, L. J.↗

Increasing the Thermodynamic Driving Force of the Phosphofructokinase Reaction in Clostridium thermocellum

Glycolysis is an ancient, widespread, and highly conserved metabolic pathway that converts glucose into pyruvate. In the canonical pathway, the phosphofructokinase (PFK) reaction plays an important role in controlling flux through the pathway. Clostridium thermocellum has an atypical glycolysis and uses pyrophosphate (PP i ) instead of ATP as the phosphate donor for the PFK reaction. The reduced thermodynamic driving force of the PP i -PFK reaction shifts the entire pathway closer to thermodynamic equilibrium, which has been predicted to limit product titers. Here, in this study, we replace the PP i -PFK reaction with an ATP-PFK reaction. We demonstrate that the local changes are consistent with thermodynamic predictions: the ratio of fructose 1,6-bisphosphate to fructose-6-phosphate increases, and the reverse flux through the reaction (determined by 13 C labeling) decreases. The final titer and distribution of fermentation products, however, do not change, demonstrating that the thermodynamic constraints of the PP i -PFK reaction are not the sole factor limiting product titer.

59 BASIC BIOLOGICAL SCIENCES↗