SEARCH · Search NASA
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.
LLEPE (liquid-liquid equilibrium parameter estimator) is a tool for estimating thermodynamic parameters like standard enthalpies, standard molar volumes and Pitzer coefficients, by combining Gibbs minimization with regression to multiphase equilibrium isotherm data. It also includes visualization tools to access data quality and prediction model performance.
H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.
Adsorption is an important phenomenon in surface chemistry, especially for nano-porous shales. In the shale gas-in-place, adsorbed gas could contribute up to 85%. However, adsorption is hard to quantitatively characterize due to different adsorption mechanisms, patterns, surfaces, and pore sizes. Moreover, key thermodynamic parameters, such as the enthalpy of adsorption, are challenging to determine due to uncertainties in adsorbed gas densities used in constructing absolute isotherms. In this paper, we revisit the Brunauer, Emmett, and Teller (BET) model, analytically simplify Ono-Kondo (OK) models for subsurface shales, and compare commonly used mono- and multilayer adsorption models (e.g., Langmuir, supercritical Dubinin-Radushkevich (SDR), supercritical BET (SBET), and simplified OK (OKs) models) with recently proposed pressure-dependent adsorption densities to develop a practical and reliable methodology that can be used in the supercritical state, typical for subsurface black shale conditions.Three independent data sets were used for nitrogen and methane adsorption isotherms at different temperatures. Furthermore, we demonstrate that adsorption predicted by the SDR model is comparable or lower than that predicted by the SBET model, but higher than the amounts predicted by Langmuir and OK models. The nitrogen BET method tends to underestimate the accessible SSA for methane. Measurement of isosteric heat of adsorption is suggested instead of the experimental fitting method, due to the significant difference between the calculated results by two commonly used methods. Experimental fitting and simulation methods are also briefly reviewed to guide future research on shale gas adsorption.
Metal-organic frameworks (MOFs) with coordinatively unsaturated open metal sites (OMSs) are promising sorbent materials in atmospheric water harvesting (AWH) systems at low relative humidity (RH) due to their strong interactions with water molecules. However, accurate computational modelling of water adsorption in those materials are challenging, as standard force fields (FFs) based on the 12-6 Lennard-Jones (L-J) potential cannot properly describe the water-OMS interactions. In biomolecular simulations, the 12-6-4 L-J potential has been successfully used to model metal ion-water interactions in which the extra 1/r 4 -term is for charge-induced dipole electrostatics. In this work, we adopted this strategy and used the 12-6-4 L-J potential to model water-OMS interactions in MOFs for low RH AWH applications. Notably, without any modifications, the parameters used for aqueous metal ions were able to greatly improve the accuracy in predicting water adsorption isotherms in MOF-74, which highlights the simplicity and transferability of this method.
The high density of aluminum nanocrystals (>10 21 m -3 ) that develop during the primary crystallization in Al-based metallic glasses indicates a high nucleation rate (~10 18 m -3 s -1 ). Several studies have been advanced to account for the primary crystallization behavior, but none have been developed to completely describe the reaction kinetics. Recently, structural analysis by fluctuation electron microscopy has demonstrated the presence of the Al-like medium range order (MRO) regions as a spatial heterogeneity in as-spun Al 88 Y 7 Fe 5 metallic glass that is representative for the class of Al-based amorphous alloys that develop Al nanocrystals during primary crystallization. From the structural characterization, an MRO seeded nucleation configuration is established, whereby the Al nanocrystals are catalyzed by the MRO core to decrease the nucleation barrier. The MRO seeded nucleation model and the kinetic data from the delay time (τ) measurement provide a full accounting of the evolution of the Al nanocrystal density (N v ) during the primary crystallization under isothermal annealing treatments. Moreover, the calculated values of the steady state nucleation rates (J ss ) predicted by the nucleation model agree with the experimental results. Moreover, the model satisfies constraints on the structural, thermodynamic, and kinetic parameters, such as the critical nucleus size, the interface energy, and the volume-free energy driving force that are essential for a fully self-consistent nucleation kinetics analysis. In conclusion, the nucleation kinetics model can be applied more broadly to materials that are characterized by the presence of spatial heterogeneities.
Small-volume, high-throughput screening techniques are sought to enable downselection from a large candidate pool of bio-blendstocks to a select few, having physical properties consistent with requirements of downsized, turbo-boosted internal combustion engines. Herein, this work presents a droplet evaporation-based approach to predict heat of vaporization, vapor pressure, diffusion coefficient, and Lennard–Jones parameters for an unknown fuel. Two different schemes, considering the isothermal evaporation of a moving droplet in ambient air, are proposed, which combine droplet velocity and temperature measurements, with some known properties to predict unknown properties. The schemes utilize an inverse solution of a transient model of droplet evaporation solved in an iterative fashion. A baseline scheme, which only requires droplet size change measurements, is evaluated using test data for three liquid fuels, comprising of alkanes and alcohols, as obtained in a temperature-controlled chamber. Results yield temperature-dependent heat of vaporization and vapor pressure predictions within 10 % and 22 %, respectively, of reference values. The advanced scheme, which additionally requires droplet temperature measurement, is numerically evaluated in the current work and will be experimentally validated in future efforts. The advanced scheme is found to significantly improve prediction quality, with deviations less than 2 % and 1 % for heat of vaporization and vapor pressure, while also predicting diffusion coefficient and Lennard–Jones parameters within 5 % and 8 %, respectively. The combined set of approaches, which primarily track droplet evaporation, can be incorporated into a small-volume, high-throughput fuel screening process.
Synchrotron X-ray diffraction has been used to investigate the structure and equation of state (EOS) of hydrazine (N 2 H 4 ) up to 54.3 GPa at 298 K. The diffraction patterns could be fit to a monoclinic unit-cell structure and put strong constraints on previously reported phase transitions documented by vibrational spectroscopy over this pressure range. Pressure–volume ( P–V ) data were fit using a Vinet EOS, yielding parameters: V 0 = 45.2 Å 3 /molecule (fixed), K 0 = 11.8(7) GPa, and K 0 ′ = 6.5(2). Previously measured high-pressure vibrational frequency shifts were used to estimate the vibrational free energy and model P–V–T isotherms from 0 to 1200 K. The results of the P–V–T isotherms are compared to existing shock Hugoniot data on hydrazine and 298 K isotherms for assemblages of possible decomposition products. This comparison suggests dissociation at high density under shock loading. Good correspondence was found between the static lattice EOS as calculated by the model and the previously reported EOS as calculated by density functional theory. Finally, these results resolve existing uncertainties about the EOS and crystal symmetry of hydrazine at high pressure and provide valuable baseline information on this important energetic material.
Visible photon fluxes can influence the rate and selectivity of heterogeneously catalyzed reactions on metal nanoparticle surfaces. Models describing the influence of photon fluxes have typically introduced photon flux dependent apparent thermal kinetic parameters (reaction orders, activation energies, binding energies, etc.). This has relied on empirical fitting of reaction rate data, making mechanistic interpretations of how photon fluxes influence elementary step rates challenging and inconsistent with fundamental descriptions of photochemistry on metal surfaces developed from surface science studies. Using the CO adsorption–desorption quasi-equilibrium reaction on Pt/Al 2 O 3 catalysts as a model system, we measured steady state adsorbed CO (CO*) coverages under isothermal and isobaric (1 mbar CO) conditions as a function of temperature (473–573 K) and of 440 nm photon flux ((0.1–5.2) × 10 3 # hv Pt site –1 s –1 ) using in situ IR spectroscopy. Steady state CO* coverage on Pt was photon flux dependent with increasing photon flux causing decreasing coverage, consistent with photons driving CO* desorption rates faster than thermal CO* desorption rates. However, photon flux dependent CO* coverages were essentially temperature independent, inconsistent with models that describe photon effects using perturbations to apparent thermal kinetic parameters. Instead, 120 steady state CO* coverages as a function of temperature and photon flux are quantitatively described by a kinetic model in which the overall desorption rate is a summation of independent thermal and photon induced CO* desorption rates. Site-resolved analysis reveals distinct kinetic parameters for photon driven desorption of CO* from well-coordinated, under-coordinated, and highly under-coordinated Pt sites, with temperature-dependent apparent quantum efficiencies (AQE) consistent with temperature dependence of vibrational quanta distribution of adsorbed CO. The rigorous kinetic rate laws for independent photon and thermal driven pathways allow for predictive modeling of the influence of photon fluxes on the rates of CO* desorption under catalytic conditions. Further, the analysis provides evidence that steady state continuous wave photon fluxes can drive desorption/adsorption reactions on metal surfaces out of thermal equilibrium, reconciling surface science observations of molecular photodesorption with applied catalysis. The work establishes a general kinetic framework to be considered for photon driven processes on metals, and defines catalyst, reaction, and photon flux characteristic design principles for breaking Sabatier limitations.
Adsorbate molecules present in a reaction mixture may bind to and block catalytic sites. Measurement of the surface coverage of these molecules via adsorption isotherms is critical for modeling and design of catalytic reactions on surfaces. However, it is challenging to measure isotherms in solution in a way that is directly relevant to catalytic activity under reaction conditions, particularly since adsorbates may bind with an enormous range of surface affinity parameters. Here we used the motion of self-propelled catalytic Janus particles, which employ the decomposition of hydrogen peroxide fuel as a propulsion mechanism, to determine the effective surface coverage of thioglycerol, furfural, and ethanol on a platinum surface as a function of concentration in aqueous solution by measuring the decrease in active motion due to the blocking of active sites. For strongly adsorbing thioglycerol, this effective coverage was compared and contrasted to the total adsorbed amount measured using inductively-coupled plasma analysis. Demonstrating the broad applicability of this approach, the surface affinity of the three adsorbates spanned more than four orders of magnitude. For each species, the adsorbate-mediated attenuation of active motion occurred over a wide concentration range and was well-described by a Langmuir isotherm. The strongly interacting thioglycerol had the highest affinity towards the surface (K a = 15.5 ± 4.3 mM –1 ) and fully deactivated the active particle motion at surface saturation. Furfural had an intermediate affinity (K a = 0.42 ± 0.07 mM –1 ) but did not fully block H 2 O 2 access to the surface at apparent saturation, consistent with a maximum fractional surface coverage of θ max = 0.67. Ethanol exhibited even lower affinity (K a = 0.0025 ± 2x10 -4 mM –1 ) and its coverage saturated at only θ max = 0.38. Analysis of isotherms at elevated temperatures enabled direct extraction of the enthalpies of adsorption. The degree of surface coverage at adsorbate saturation appeared to correlate with the relative energies of adsorption for the different adsorbate species and was consistent with adsorbate saturation of one of multiple active site populations towards H 2 O 2 decomposition. Furthermore, computational investigations into solvent effects on furfural adsorption showed good quantitative agreement with the experimental results. This work leverages unique properties of active particles to explore fundamental catalysis questions and demonstrates a novel paradigm for significant and experimentally accessible multidisciplinary research.
The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.
We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of cold dark matter and a strongly dissipative subcomponent. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as atomic dark matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of ∼10 10 M ⊙ systems that vary over the aDM microphysics and the dwarf’s evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts, and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark “clumps,” which are clusters of subresolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction (∼5%) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations.
The leakage and accumulation of plastic in the environment is a significant and growing problem with numerous detrimental impacts and has led to a push toward the design and development of more environmentally benign materials. To this end, we have developed a quantum chemistry-based model for predicting the mobility of polymer materials from molecular structure. Hydrophobicity is used as a surrogate for mobility given that hydrophobic interactions drive much of the partitioning of contaminants in and out of various environmentally relevant compartments. To model polymer hydrophobicity, we adjusted a previously developed Quantum Chemically Calculated Abraham Parameter model to calculate Abraham parameters of small molecules from molecular structure information. The resulting model predicted the octanol-water partition coefficient (K OW ) of polymer repeating units with a root mean square error (RMSE) of 0.48 (log scale). Additionally, the hydrophobicity of high molecular weight polymer materials was captured through solubility parameters and Nile red staining experiments from the literature and predicted with RMSEs of 1.21 (J/cc) 0.5 and 3.42 nm, respectively. Finally, to test the environmental applicability of the model, the relative adsorption capacity of three polymers was predicted and used to unify sorption isotherms across multiple sorbates and polymer sorbents.
Detailed characterization of gas transport and storage in low-permeability and organic-rich shales is associated with an array of challenges due to the complex morphology of the pore space, representing a broad range of pore sizes, combined with the heterogeneous fabric of the shale matrix. These factors, and their interplay during gas transport and sorption, complicate a) the analysis of shale samples at the laboratory scale (~ up to a ft. In length), and b) the prediction of natural gas production and carbon sequestration potential at larger scales. Here, in this work, tandem experiments with inert (helium – He) and adsorbing (krypton - Kr and carbon dioxide - CO 2 ) gases were performed and analyzed to develop an efficient workflow for characterizing and modeling transport and sorption at the laboratory scale. In particular, pressure pulse-decay (PPD) measurements were conducted on an Eagle Ford shale core sample at room temperature using He, Kr, and CO 2 . PPD measurements with He (a non-sorbing gas) were employed to probe the overall porosity, including natural fractures, microcracks, mesopores, and micropores. A triple-porosity model (TPM) was adopted to interpret the gas transport in the shale sample: The pore space is represented by three interacting continua, including macropores (larger fractures), mesopores (including microcracks), and micropores. To facilitate the application of the TPM, a modified analytical approach is introduced to extract effective transport parameters (in terms of the characteristic time for transport at relevant porosity levels) directly from the PPD measurements with He. To validate the modeling approach, model parameters extracted from two He PPD measurements are demonstrated to provide excellent agreement with a 3rd He PPD measurement performed at a higher pressure. The effective transport parameters, extracted from the He experiments, were subsequently converted for application to the Kr and CO 2 PPD experiments, by accounting for relevant transport modes and differences in fluid properties. Excess adsorption isotherms were extracted from the equilibrium pressures of the Kr and CO 2 PPD experiments using the He pore volume as a baseline. These adsorption isotherms were then integrated into the TPM to predict combined gas transport and sorption for Kr and CO 2 , with transport coefficients translated from the He measurements. The predictions for Kr and CO 2 are demonstrated to be in excellent agreement with the experimental observations. This, in turn, demonstrates that the proposed analytical approach provides for an effective characterization of mass transfer rates in shales, that can be applied directly in a TPM representation of mass transfer and sorption.
Abstract There exist many time-temperature parameter (TTP) models for creep rupture prediction of components including the Larson–Miller (LM), Manson–Haferd (MH), Manson–Brown (MB), Orr–Sherby–Dorn (OSD), Manson–Succop (MS), Graham–Walles (GW), Chitty–Duval (CD), Goldhoff–Sherby (GS) models. It remains a challenge to determine which model is “best”, capable of accurate interpolation and physically realistic extrapolation of creep rupture data for a given material. In this study, metamodeling is applied to create a unified TTP metamodel that combines and regresses into twelve TTP models (eight existing and four newly derived). An analysis of the mathematical problems that exist in TTP models is provided. A matlab code is written that can: (1) calibrate the material constants of any of the twelve TTP models (using the metamodel); (2) determine the most suitable stress-parameter function; (3) and report the normalized mean square error (NMSE) of rupture predictions for a given material database. Using the metamodel, and code, a design engineer can make an intelligent selection of the “best” TTP model for creep resistant design. This process is demonstrated using four isotherms of alloy P91 creep rupture data. To assess the influence of material, further validation is performed on alloys Hastelloy X, 304SS, and 316SS. It is determined that the “best” model is dependent on material type and the quality and quantity of available data.
A series of batch experiments were performed to assess the uranium sorption capacity of four mineralogically distinct lithologies from the Negev Desert, Israel, to evaluate the suitability of a potential site for subsurface radioactive waste disposal. The rock specimens consisted of an organic-rich phosphorite, a bituminous marl, a chalk, and a sandstone. The sorption data for each lithology were fitted using a general composite surface complexation model (GC SCM) implemented in PHREEQC. Sorption data were also fitted by a non-mechanistic Langmuir sorption isotherm, which can be used as an alternative to the GC SCM to provide a more computationally efficient method for uranium sorption. This is because all the rocks tested have high pH/alkalinity/calcium buffering capacities that restrict groundwater chemistry variations, so that the use of a GC SCM is not advantageous. The mineralogy of the rocks points to several dominant sorption phases for uranyl (UO 2 2+ ), including apatite, organic carbon, clays, and iron-bearing phases. The surface complexation parameters based on literature values for the minerals identified overestimate the uranium sorption capacities, so that for our application, an empirical approach that makes direct use of the experimental data to estimate mineral-specific sorption parameters appears to be more practical for predicting uranium sorption.
High-entropy alloys (HEAs) represent an interesting alloying strategy that can yield exceptional performance properties needed across a variety of technology applications, including hydrogen storage. Examples include ultrahigh volumetric capacity materials (BCC alloys → FCC dihydrides) with improved thermodynamics relative to conventional high-capacity metal hydrides (like MgH 2 ), but still further destabilization is needed to reduce operating temperature and increase system-level capacity. Here, in this work, we demonstrate efficient hydride destabilization strategies by synthesizing two new Al 0.05 (TiVNb) 0.95–x Mo x (x = 0.05, 0.10) compositions. We specifically evaluate the effect of molybdenum (Mo) addition on the phase structure, microstructure, hydrogen absorption, and desorption properties. Both alloys crystallize in a bcc structure with decreasing lattice parameters as the Mo content increases. The alloys can rapidly absorb hydrogen at 25 °C with capacities of 1.78 H/M (2.79 wt %) and 1.79 H/M (2.75 wt %) with increasing Mo content. Pressure-composition isotherms suggest a two-step reaction for hydrogen absorption to a final fcc dihydride phase. The experiments demonstrate that increasing Mo content results in a significant hydride destabilization, which is consistent with predictions from a gradient boosting tree data-driven model for metal hydride thermodynamics. Furthermore, improved desorption properties with increasing Mo content and reversibility were observed by in situ synchrotron X-ray diffraction, in situ neutron diffraction, and thermal desorption spectroscopy.
Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content
Water absorption plays a key role in the performance of polymeric anion exchange membranes. It influences important properties such as ionic conductivity and mechanical strength and alters their performance as solid electrolytes in hydrogen electrochemical devices for energy conversion. However, computational approaches that address the relationship between the polymer design and the absorption process are scarce. In this work, we introduce a simple thermodynamic model to predict the water absorption isotherms of polyelectrolyte membranes in contact with a water vapor reservoir that incorporates the specific chemical design of the polymers. The model accurately predicts the water content and macrostructural properties of polynorbornene membranes as a function of the water activity and successfully captures the effect of various polymer design parameters. The energy of pairwise attractive interactions predicted by our model provides a means to interpret the absorption process at the molecular level. The model also reveals the most significant favorable and unfavorable contributions to the free energy and indicates that their balance is solely governed by the water volume fraction, regardless of the polymer design. This universal behavior leads to important implications in the search for better ion exchange membranes.