SEARCH · Engineering Papers

Results for “isotherm model parameters”

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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Scalar dissipation and mixing in turbulent reacting flows

The mixing term that appears in probability density function-based formulations of the statistical behavior of turbulent reacting flows was studied with the use of direct numerical simulations of an irreversible, second-order, isothermal chemical reaction between initially segregated reactants in decaying, homogeneous turbulence. Three-dimensional, time-dependent pseudospectral calculations with 64-cubed Fourier modes were made as well as were Monte Carlo simulations of the joint concentration pdf equations for coalescence-dispersion (C/D) models. It was found that the development of the scalar microscales depends largely on the initial concentration distributions (both isotropic and nonisotropic distributions were used) and is relatively independent of the reaction rate constant. Joint concentration probability density functions were compared to predictions of D/D models that make use of parameters evaluated from the nonreacting or inert system.

Leonard, Andy D.↗

NASA NTRS2010

Herschel-ATLAS: Dust Temperature and Redshift Distribution of SPIRE and PACS Detected Sources Using Submillimetre Colours

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500/microns using both PACS and SPIRE. We fit isothermal modified-blackbody spectral energy distribution (SED) models in order to extract the dust temperature of sources with counterparts in GAMA or SDSS with either a spectroscopic or a photometric redshift. For a subsample of 331 sources detected in at least three FIR bands with significance greater than 30 sigma, we find an average dust temperature of (28 plus or minus 8)K. For sources with no known redshifts, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters and compare to colours of observed sources to establish the redshift distribution of those samples. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3sigma, we find an average redshift of 2.2 plus or minus 0.6.

Amblard, A.↗

NASA NTRS1994

The effects of physical aging at elevated temperatures on the viscoelastic creep on IM7/K3B

Physical aging at elevated temperature of the advanced composite IM7/K3B was investigated through the use of creep compliance tests. Testing consisted of short term isothermal, creep/recovery with the creep segments performed at constant load. The matrix dominated transverse tensile and in-plane shear behavior were measured at temperatures ranging from 200 to 230 C. Through the use of time based shifting procedures, the aging shift factors, shift rates and momentary master curve parameters were found at each temperature. These material parameters were used as input to a predictive methodology, which was based upon effective time theory and linear viscoelasticity combined with classical lamination theory. Long term creep compliance test data was compared to predictions to verify the method. The model was then used to predict the long term creep behavior for several general laminates.

Gates, Thomas S.↗

DOE OSTI2022

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

NASA NTRS1987

Calculation of thermomechanical fatigue life based on isothermal behavior

The isothermal and thermomechanical fatigue (TMF) crack initiation response of a hypothetical material was analyzed. Expected thermomechanical behavior was evaluated numerically based on simple, isothermal, cyclic stress-strain - time characteristics and on strainrange versus cyclic life relations that have been assigned to the material. The attempt was made to establish basic minimum requirements for the development of a physically accurate TMF life-prediction model. A worthy method must be able to deal with the simplest of conditions: that is, those for which thermal cycling, per se, introduces no damage mechanisms other than those found in isothermal behavior. Under these assumed conditions, the TMF life should be obtained uniquely from known isothermal behavior. The ramifications of making more complex assumptions will be dealt with in future studies. Although analyses are only in their early stages, considerable insight has been gained in understanding the characteristics of several existing high-temperature life-prediction methods. The present work indicates that the most viable damage parameter is based on the inelastic strainrange.

Halford, Gary R.↗

NASA NTRS1987

Calculation of thermomechanical fatigue life based on isothermal behavior

The isothermal and thermomechanical fatigue (TMF) crack initiation response of a hypothetical material was analyzed. Expected thermomechanical behavior was evaluated numerically based on simple, isothermal, cyclic stress-strain-time characteristics and on strainrange versus cyclic life relations that have been assigned to the material. The attempt was made to establish basic minimum requirements for the development of a physically accurate TMF life-prediction model. A worthy method must be able to deal with the simplest of conditions: that is, those for which thermal cycling, per se, introduces no damage mechanisms other than those found in isothermal behavior. Under these assumed conditions, the TMF life should be obtained uniquely from known isothermal behavior. The ramifications of making more complex assumptions will be dealt with in future studies. Although analyses are only in their early stages, considerable insight has been gained in understanding the characteristics of several existing high-temperature life-prediction methods. The present work indicates that the most viable damage parameter is based on the inelastic strainrange.

Halford, G. R.↗

NASA NTRS1990

Application of the sine-Poisson equation in solar magnetostatics

Solutions of the sine-Poisson equations are used to construct a class of isothermal magnetostatic atmospheres, with one ignorable coordinate corresponding to a uniform gravitational field in a plane geometry. The distributed current in the model (j) is directed along the x-axis, where x is the horizontal ignorable coordinate; (j) varies as the sine of the magnetostatic potential and falls off exponentially with distance vertical to the base with an e-folding distance equal to the gravitational scale height. Solutions for the magnetostatic potential A corresponding to the one-soliton, two-soliton, and breather solutions of the sine-Gordon equation are studied. Depending on the values of the free parameters in the soliton solutions, horizontally periodic magnetostatic structures are obtained possessing either a single X-type neutral point, multiple neural X-points, or solutions without X-points.

Webb, G. M.↗

NASA NTRS2017

Early Hydrodynamic Escape Limits Rocky Planets to Less Than or Equal to 1.6 Earth Radii

In the past decade thousands of exoplanet candidates and hundreds of confirmed exoplanets have been found. For sub-Neptune-sized planets, those less than approx. 10 Earth masses, we can separate planets into two broad categories: predominantly rocky planets, and gaseous planets with thick volatile sheaths. Observations and subsequent analysis of these planets show that rocky planets are only found with radii less than approx. 1.6 Earth radii. No rocky planet has yet been found that violates this limit. We propose that hydrodynamic escape of hydrogen rich protoatmospheres, accreted by forming planets, explains the limit in rocky planet size. Following the hydrodynamic escape model employed by Luger et al. (2015), we modelled the XUV driven escape from young planets (less than approx.100 Myr in age) around a Sun-like star. With a simple, first-order model we found that the rocky planet radii limit occurs consistently at approx. 1.6 Earth radii across a wide range of plausible parameter spaces. Our model shows that hydrodynamic escape can explain the observed cutoff between rocky and gaseous planets. Fig. 1 shows the results of our model for rocky planets between 0.5 and 10 Earth masses that accrete 3 wt. % H2/He during formation. The simulation was run for 100 Myr, after that time the XUV flux drops off exponentially and hydrodynamic escape drops with it. A cutoff between rocky planets and gaseous ones is clearly seen at approx. 1.5-1.6 Earth radii. We are only interested in the upper size limit for rocky planets. As such, we assumed pure hydrogen atmospheres and the highest possible isothermal atmospheric temperatures, which will produce an upper limit on the hydrodynamic loss rate. Previous work shows that a reasonable approximation for an upper temperature limit in a hydrogen rich protoatmosphere is 2000-3000 K, consistent with our assumptions. From these results, we propose that the observed dichotomy between mini-Neptunes and rocky worlds is simply explained by an early episode of thermally-driven hydrodynamic escape when host stars have saturated XUV fluxes.

Lehmer, O. R.↗

DOE OSTI2021

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗

DOE OSTI2020

The lattice parameter – composition relationship of the body centered cubic uranium-niobium alloys

In this paper, in-situ time-of-flight neutron diffraction experiments were performed on uranium-niobium alloy with 6 wt% Nb to study the lattice parameter-composition relationship for the body centered cubic (bcc) alloys. Based on lattice-parameter measurements of Nb-supersaturated γs phase over the range of 450 – 790 °C, an improved Vegard's-type relationship was established for determination of Nb concentrations at elevated temperatures. Neutron diffraction data were also collected as a function of aging time over the 450 - 600 °C range when γs phase decomposed isothermally into the mixture of orthorhombic α-U and Nb-rich bcc phase, as well as at room-temperature on the alloys aged ex-situ at 500 °C up to five years. From these measurements, the average Nb concentrations in the bcc phase were determined based on the Rietveld refinements of weight fraction and mass conservation relations (lever rule). Over the 15at% – 78at% range of Nb concentrations that correspond to different experimental aging times, the lattice parameters at constant temperatures exhibit a nonlinear S-shaped variation with Nb concentration, and the associated excess volumes of mixing can be described by a subregular solution model of the Redlich–Kister type of polynomial. Over the full range of composition the S-shaped deviation from Vegard's law can be modeled using a combination of an elastic continuum model and a perturbation to the radii of the solute atoms in the solvent, suggesting that electronic interactions between solute and solvent atoms could play an important role in the compositional dependence of lattice parameter for the γ-phase U-Nb alloys. While Vegard's law is a straightforward and reasonably good approximation for the bcc solid solutions in the U-Nb system, the Nb concentrations determined from the weight fraction refinements of diffraction data provide internally consistent, mass-conserving estimates of solute redistribution for the monotectoid reaction.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Corrosion Testing of Refractory inContact with Molten Glasses Designed for Waste Vitrification - VSL Touchpoint Matrix Glasses

It is known that the predictive life of the refractory ceramic liner of nuclear waste glass melters is conservative, as demonstrated by performance of these materials such as in the Defense Waste Processing Facility (DWPF). The motivation for this task is to maximize the useful life of the melters that will be operated at the Waste Treatment and Immobilization Plant (WTP), which will in turn minimize procurement and disposal costs and melter outage times, as well as to identify maximum loadings in the waste glass of those species that corrode melter components. This task was initiated jointly with Pacific Northwest National Laboratory (PNNL) with the objective to develop a methodology and model to enable more accurate prediction of refractory service life under prototypic conditions from laboratory-scale material corrosion tests. Refractory corrosion is generally reported as physical material loss, measured in units of distance (e.g., inch) or as physical material loss rate, measured in units of distance per time (e.g., inch/day). In post-operational melters, the refractory corrosion is measured directly, sometimes reported as corrosion depth. Crucible tests are used in the laboratory to accelerate the refractory corrosion to facilitate a meaningful measurement in a commensurate amount of time. Crucible tests are particularly useful in understanding refractory corrosion across a large glass composition space, where operational testing would be prohibitive. Some of the critical parameters that are known to influence refractory corrosion by molten glass in a crucible test are temperature, system redox, molten salt phases, glass chemistry, and test duration. The majority of data collected for Monofrax® K-3 (hereafter referred to as K-3) corrosion is from crucible tests, but a small amount comes directly from scaled and production melters. Crucible test data has been collected under varying conditions, whereas data collected from operational melters is relatively fewer and represents conditions specific to the melter campaign. The result is that the published data can be grouped and analyzed in multiple ways, not all of which are readily comparable. The Standard Test Method for Isothermal Corrosion Resistance of Refractories to Molten Glass (ASTM C621) outlines the general guidelines used across industry. That method describes a sealed, static test in which the surface area of the refractory coupon and the volume of glass are fixed. A significant portion of the crucible data pertaining to nuclear waste glasses has been collected in a modified configuration; the most notable differences being the surface area of the refractory coupon to volume of the glass and use of a method for bubbling the melt. To our knowledge, the influence of those parameters on the refractory corrosion has not been quantified. In this work, it was determined that static tests and bubbled tests would be performed. Savannah River National Laboratory (SRNL) was tasked with setting up and performing static testing while PNNL was tasked with setting up and performing bubbled testing. Initial activities were performed to establish laboratory methods that reproduce data comparable to existing data sets of K-3 refractory corrosion by low activity waste (LAW) and high-level waste (HLW) glass compositions. Later activities were focused on refining the test parameters to establish a standard test practice to be used between Laboratories and collecting additional data to be used in the enhanced waste glass model development. This document serves primarily to convey the refractory loss measurement results from corrosion testing of K-3 refractory with waste glass compositions developed for use in the WTP melters. The data will be used in the enhanced property/composition models being developed for waste glass vitrification and melter operations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1974

Thermodynamic properties derived from the free volume model of liquids

An equation of state and expressions for the isothermal compressibility, thermal expansion coefficient, heat capacity, and entropy of liquids have been derived from the free volume model partition function suggested by Turnbull. The simple definition of the free volume is used, and it is assumed that the specific volume is directly related to the cube of the intermolecular separation by a proportionality factor which is found to be a function of temperature and pressure as well as specific volume. When values of the proportionality factor are calculated from experimental data for real liquids, it is found to be approximately constant over ranges of temperature and pressure which correspond to the dense liquid phase. This result provides a single-parameter method for calculating dense liquid thermodynamic properties and is consistent with the fact that the free volume model is designed to describe liquids near the solidification point.

Miller, R. I.↗

NASA NTRS1978

Infrared excesses in early-type stars - Gamma Cassiopeiae

Spectrophotometry of the classical Be star Gamma Cas (1-4 microns, with about 2% spectral resolution) is presented. These data, together with existing broad-band observations, are accurately described by simple isothermal LTE models for the IR excess which differ from most previously published work in three ways: (1) hydrogenic bound-free emission is included; (2) the attenuation of the star by the shell is included; and (3) no assumption is made that the shell contribution is negligible in some bandpass. It is demonstrated that the bulk of the IR excess consists of hydrogenic bound-free and free-free emission from a shell of hot ionized hydrogen gas, although a small thermal component cannot be ruled out. The bound-free emission is strong, and the Balmer, Paschen, and Brackett discontinuities are correctly represented by the shell model with physical parameters as follows: a shell temperature of approximately 18,000 K, an optical depth (at 1 micron) of about 0.5, an electron density of approximately 1 trillion per cu cm, and a size of about 2 trillion cm. Phantom shells (i.e., ones which do not alter the observed spectrum of the underlying star) are discussed.

Scargle, J. D.↗

NASA NTRS1991

In-phase and out-of-phase axial-torsional fatigue behavior of Haynes 188 at 760 C

Isothermal, in-phase and out-of-phase axial-torsional fatigue experiments have been conducted at 760 C on uniform gage section, thin-walled tubular specimens of a wrought cobalt-base superalloy, Haynes 188. Test-control and data acquisition were accomplished with a minicomputer. Fatigue lives of the in- and out-of-phase axial-torsional fatigue tests have been estimated with four different multiaxial fatigue life prediction models that were developed primarly for predicting axial-torsional fatigue lives at room temperature. The models investigated were: (1) the von Mises equivalent strain range; (2) the Modified Multiaxiality Factor Approach; (3) the Modified Smith-Watson-Topper Parameter; and (4) the critical shear plane method of Fatemi, Socie, and Kurath. In general, life predictions by the von Mises equivalent strain range model were within a factor of 2 for a majority of the tests and the predictions by the Modified Multiaxiality Factor Approach were within a factor of 2, while predictions of the Modified Smith-Watson-Topper Parameter and of the critical shear plane method of Fatemi, Socie, and Kurath were unconservative and conservative, respectively, by up to factors of 4. In some of the specimens tested under combined axial-torsional loading conditions, fatigue cracks initiated near extensometer indentations. Two design modifications have been proposed to the thin-walled tubular specimen to overcome this problem.

Kalluri, Sreeramesh↗

DOE OSTI2024

High Performance Heat Pipe Power Transient Testing at SPHERE Facility

Microreactors are being researched, designed, and built at Idaho National Laboratory (INL). Microreactors are small reactors defined at less than 20MW of power. These reactor concepts are also being looked at throughout industry for various applications. An important aspect of these reactor designs is economic feasibility i.e. lower overnight capital cost. The driving factors for implementing microreactors are quick setup and takedown, minimal operators, and the ability to manufacture them readily and to fit in mid-sized containers for transport. A specific area of research to aid in successful integration of these factors within the designs is passive heat removal of the core’s thermal power. Interest in heat pipes to achieve this passive heat removal has been shown across multiple industry partners. Because of this interest, INL has developed a test facility to facilitate experimental tests for sodium filled heat pipes. INL has developed the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility to run experiments on high performance, sodium filled heat pipes. As mentioned above, heat pipes are passive heat transfer devices. Radially, heat pipes are broken up into an outer wall, a small annular gap, a wick structure, and a centerline gap. They function by utilizing latent heat transfer. Heat pipes are traditionally separated into three regions, an evaporator (heat input), an adiabatic region, and finally a condenser region (heat removal). As heat is being applied to the evaporator, the working fluid undergoes a phase change to a vapor. This phase change causes a differential pressure across the axial length of the pipe driving flow down the center gap of the heat pipe. The vapor flows down past the adiabatic region to the condenser where the heat is removed. This heat removal forces the working fluid to phase change back to a liquid. The wick structure is then utilized to drive the flow back towards the evaporator by capillary forces. This backflow is aided by the annular gap. Because this heat transfer mechanism functions with latent heat transfer, the heat pipe is close to isothermal down the axial length. Heat pipes can operate under a wide range of working fluids. Considerations for these working fluids are primarily driven by operating temperatures amongst other important factors based around overall performance. Sodium filled heat pipes operate from 450°C up to 900°C. This temperature range works well for the current microreactor designs. In conjunction with this experimental capability, INL has developed a modeling software to simulate heat pipe physics within reactor cores. This modeling software is called Sockeye and functions under the established INL Multiphysics Object Oriented Simulation Environment (MOOSE). SPHERE also supports Sockeye development by providing the modeling team with experimental data on an array of setups and operating parameters to support validation efforts. A power transient experiment was performed utilizing the SPHERE facility to continue to aid with Sockeye development. The testing followed a proposed test plan to ramp up and down the temperature of the heat pipe. Sockeye models steady state heat pipe operation with high accuracy, the data provided by the power transient testing aims to assist with the validation efforts and further enhance transient modeling capability of the tool [2].

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1973

Laminar film boiling on inclined isothermal flat plates.

Laminar film boiling from an inclined flat plate has been investigated analytically. Using the singular perturbation scheme, the complete set of Navier-Stokes equations is solved. The zeroth-order perturbation coinciding with the boundary-layer equations for vertical flat plates governs the problem. The higher-order perturbations become important near the leading edge and for large values of the inclination angle from the vertical. The assumption of zero interfacial velocity shows that, except for fluids having large (rho x mu) ratios, the results can be predicted using the vertical flat plate results by defining a modified Grashof parameter containing a cos phi term. When the interfacial shear is considered, the solutions indicate that for fluids having large (rho x mu) ratios, the heat transfer rates will be larger (approximately 15% maximum) than those predicted by the simplified model using zero interfacial velocity. In general, the inclination decreases the rate of heat transfer as well as the rate of evaporation.

Nagendra, H. R.↗

DOE OSTI2025

Machine-learning and first-principles investigation of lightweight medium-entropy alloys for hydrogen-storage applications

The transition to a low-carbon economy demands efficient and sustainable energy-storage solutions, with hydrogen emerging as a promising clean-energy carrier and with metal hydrides recognized for their hydrogen-storage capacity. Here, we leverage machine learning (ML) to predict hydrogen-to-metal (H/M) ratios and solution energy by incorporating thermodynamic parameters and local lattice distortion (LLD) as key features. Our best-performing ML model provides improvements to H/M ratios and solution energies over a broad class of medium-entripy alloys (easily extendable to multi-principal-element alloys), such as Ti–Nb-X (X = Mo, Cr, Hf, Ta, V, Zr) and Co–Ni-X (X = Al, Mg, V). Ti–Nb–Mo alloys reveal compositional effects in H-storage behavior, in particular Ti, Nb, and V enhance H-storage capacity, while Mo reduces H/M and hydrogen weight percent by 40–50 %. We attributed results in molybdenum-rich alloys to slow hydrogen kinetics, as validated by our pressure-composition-temperature (PCT) isotherm experiments on pure Ti and Ti 5 Mo 95 alloys. Density functional theory (DFT) and molecular dynamics (MD) simulations also confirm that Ti and Nb promote H diffusion, whereas Mo hinders it, highlighting the interplay between electronic structure, lattice distortions, and hydrogen uptake. Notably, our Gradient Boosting Regression model identifies LLD as a critical factor in H/M predictions. Here, to aid material selection, we present two periodic tables illustrating elemental effects on (a) H 2 wt% and (b) solution energy, derived from ML, and provide a reference for identifying alloying elements that enhance hydrogen solubility and storage.

08 HYDROGEN↗

NASA NTRS2016

Predictive Simulation of Gas Adsorption in Fixed-Beds and Limitations due to the Ill-Posed Danckwerts Boundary Condition

The 1-D axially dispersed plug flow model is a mathematical model widely used for the simulation of adsorption processes. Lumped mass transfer coefficients such as the Glueckauf linear driving force (LDF) term and the axial dispersion coefficient are generally obtained by fitting simulation results to the experimental breakthrough test data. An approach is introduced where these parameters, along with the only free parameter in the energy balance equations, are individually fit to specific test data that isolates the appropriate physics. It is shown that with this approach this model provides excellent simulation results for the C02 on zeolite SA sorbent/sorbate system; however, for the H20 on zeolite SA system, non-physical deviations from constant pattern behavior occur when fitting dispersive experimental results with a large axial dispersion coefficient. A method has also been developed that determines a priori what values of the LDF and axial dispersion terms will result in non-physical simulation results for a specific sorbent/sorbate system when using the one-dimensional axially dispersed plug flow model. A relationship between the steepness of the adsorption equilibrium isotherm as indicated by the distribution factor, the magnitude of the axial dispersion and mass transfer coefficient, and the resulting non-physical behavior is derived. This relationship is intended to provide a guide for avoiding non-physical behavior by limiting the magnitude of the axial dispersion term on the basis of the mass transfer coefficient and distribution factor.

Knox, James Clinton↗