Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PARTIAL PRESSURE”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Direct Observation of Morphological and Chemical Changes during the Oxidation of Model Inorganic Ligand-Capped Particles

Functionalization and volatilization are competing reactions during the oxidation of carbonaceous materials and are important processes in many different areas of science and technology. Here, we present a combined ambient pressure X-ray photoelectron spectroscopy (APXPS) and grazing incidence X-ray scattering (GIXS) investigation of the oxidation of oleic acid ligands surrounding NaYF 4 nanoparticles (NPs) deposited onto SiO x /Si substrates. While APXPS monitors the evolution of the oxidation products, GIXS provides insight into the morphology of the ligands and particles before and after the oxidation. Our investigation shows that the oxidation of the oleic acid ligands proceeds at O 2 partial pressures of below 1 mbar in the presence of X-rays, with the oxidation eventually reaching a steady state in which mainly CH x and –COOH functional groups are observed. The scattering data reveal that the oxidation and volatilization reaction proceeds preferentially on the side of the particle facing the gas phase, leading to the formation of a chemically and morphologically asymmetric ligand layer. This comprehensive picture of the oxidation process could be obtained only by combining the X-ray scattering and APXPS data. The investigation presented here lays the foundation for further studies of the stability of NP layers in the presence of reactive trace gases and ionizing radiation and for other nanoscale systems where chemical and morphological changes happen simultaneously and cannot be understood in isolation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Consistency of fatigue crack growth behavior of pipeline and low-alloy pressure vessel steels in gaseous hydrogen

Here, this study investigates the fatigue crack growth rate (FCGR) behavior of pipeline and low-alloy pressure vessel steels in high-pressure gaseous hydrogen. Despite a broad range of yield strengths and microstructures ranging from ferrite/pearlite, acicular ferrite, bainite, and martensite, the FCGR in gaseous hydrogen remained consistent (falling within a factor of 2–3). Steels with higher fractions of pearlite, typical of older vintage pipeline steels, exhibited modestly lower crack growth rates in gaseous hydrogen compared to steels with lower fractions of pearlite. Crack growth rates in these materials exhibit a systematic dependence on stress ratio and partial pressure of hydrogen, as captured in the recently published fatigue design curves in ASME B31 code case 220 for pipeline steels and ASME BPVC code case 2938 for pressure-vessel steels.

08 HYDROGEN↗
DOE OSTI2024

Nonthermal hydrogen plasma-enabled ambient, fast lignin hydrogenolysis to valuable chemicals and bio-oils

The reduction of fossil fuel resources and the ongoing surge in global energy demand have captured the interest of researchers worldwide, prompting a focus on developing renewable energy sources. For this reason, biomass conversion has emerged as a crucial pathway for renewable fuel production. Lignin, constituting 10–35% of woody biomass, represents a significant and largely untapped sustainable feedstock. Despite the potential of lignin, a substantial portion of this lignocellulosic residue remains unused, with approximately 60% considered waste. This study addresses the challenge of underutilized lignin by introducing an innovative approach to its hydrogenolysis. Despite their potential, existing hydrogenolysis methods face obstacles such as complexity, high cost, and the need for high temperatures or pressures. Herein we report a noncatalytic nonthermal hydrogen plasma method for lignin hydrogenolysis, conducted under ambient temperature and pressure conditions. Our method proves to be highly effective in breaking lignin bonds, achieving complete conversion, and generating valuable gaseous and bio-oil products including methane and aromatic dimers and monomers obtained from guaiacyl and syringyl units within the lignin structure. Our results showed an increase in gaseous products, especially methane, and aromatic monomer yields, as well as a reduction in total bio-oil and biochar yields and lignin functional groups by increasing reaction time, input power, and H2 partial pressure. This research confirms the considerable promise of utilizing noncatalytic nonthermal hydrogen plasma-assisted hydrogenolysis as an effective technique for producing gaseous and liquid fuels from lignin.

Pishva, Parsa↗
DOE OSTI2024

The atomic-level structure and stability of interfaces of Pt nanoparticles in alumina: An experimental and computational evaluation

The atomic-level structure of interfaces between Pt and a transition form of Al2O3 were studied using a combination of electron microscopy and first principles calculations. A model system of Pt nanoprecipitates in Al2O3 were formed in sapphire wafers via high-energy ion implantation of Pt followed by thermal annealing at 1000 °C in air. The Pt nanoparticles took the form of tetrahedra and truncated tetrahedra primarily bound by {111}Pt facets. The high prevalence of these facets motivated the development of density functional theory (DFT) based models of (111)Pt interfaces with six different chemical terminations of ( 2 ¯ 01 ) θ-alumina. The atomic-level structure of the Pt/Al2O3 interfaces was characterized with aberration-corrected scanning transmission electron microscopy (STEM) and the experimental images were compared to STEM image simulations of the DFT models. The model interface with Pt bonded to oxygen-terminated θ-Al2O3, with the Pt located on top of the O and with an underlying layer of octahedral Al, provided the best match to the experimental images. This interfacial termination is also the most stable for the thermal annealing conditions used based on thermodynamic calculations of the interfacial energy as a function of temperature and oxygen partial pressure. This experimentally verified model provides a basis for improving models of Pt/γ-alumina interfaces.

Clauser, A L↗
DOE OSTI2025

Catalytic hydrogenation of HMF to BHMF over copper catalysts

2,5-Bis(hydroxymethyl)furan (BHMF) is a bio-derived building block for polyester production, obtained via the hydrogenation of 5-hydroxymethylfurfural (HMF). First-principles thermodynamic equilibrium calculations indicate that this reaction is not thermodynamically limited under relevant conditions (e.g., 100 °C and high H 2 partial pressure). In this work, crude HMF was employed as the feedstock for BHMF synthesis. Initially, acidic impurities and humins were removed from unrefined HMF through filtration using a packed bed of γ-alumina. A comprehensive study of the filtration process is presented, including filtration kinetics, breakthrough curve analysis, and mathematical modeling. The purified HMF was subsequently hydrogenated over a 10 wt% CuZrO 2 catalyst, using ethanol as the reaction solvent. Batch reactions were first performed for collection of kinetic data to guide the transition to continuous flow operation. Kinetic data was collected in a fixed bed reactor at varying contact time, time on stream, temperature, and HMF concentration. This data was used to develop a kinetic model for HMF hydrogenation. Maximum BHMF production rates were achieved at 130 °C, accompanied by minor formation of byproducts from BHMF ring-opening reactions. The BHMF selectivity was 100 % at 100 °C although with lower reaction rates. Furthermore, catalyst stability tests revealed a loss of up to 50 % in catalytic activity within the first 24 h, likely due to the adsorption of HMF-derived oligomers that are not easily removed by filtration.

Crude HMF filtration↗
DOE OSTI2024

Kinetic assessment of pulp mill-derived lime mud calcination in high CO 2 atmosphere

The chemical pulping of biomass involves the recycling of calcium through the calcination of lime mud, which is mostly comprised of calcium carbonate (CaCO 3 ). Lime mud decomposes under elevated temperatures to generate calcium oxide (CaO) and carbon dioxide (CO 2 ), the kinetics of which are strongly influenced by the CO 2 partial pressure and temperature. Oxy-fuel combustion and electrified lime kilns for lime mud calcination are intriguing methods to decarbonize this highly polluting operation within the biomass pulping industry. However, the high CO 2 concentration in oxy-fuel and electrified calcination processes alters the kinetics and overall reactivity of lime mud. For the first time, a model-fitting method is used to determine the kinetic parameters for lime mud calcination under a wide range of temperatures (550 °C–1250 °C) and under different concentrations of CO2 (0 %, 15 %, 50 %, and 90 %). A kinetic model is developed that accurately predicts the reaction rates as a function of temperature and CO 2 concentration. The apparent activation of energy for lime mud calcination is elevated under a high CO 2 environment. Relative to inert gas (N 2 , Ar), the temperature window for calcination is much smaller under high CO 2 environments. The presence of Na in lime mud does not seem to affect calcination under a high CO 2 environment. Finally, particle size variation does not have a significant effect on calcination under a high CO 2 environment.

09 BIOMASS FUELS↗
DOE OSTI2024

Hydrogen-Steam Separation Using Mechanical Vapor Recompression Cycle

Solar thermochemical hydrogen and fuel production is a promising pathways for producing sustainable fuels and chemicals. One of the main challenges in the development of these systems is their low steam conversion extent, dictated by its restrictive thermodynamics requiring extremely high temperatures and low oxygen partial pressure to obtain conversions over 10 \%. While condensing the unreacted steam is technically simple, the latent heat is thus lost, which can be larger than the producd hydrogen higher heating value. We propose to use a modified mechanical vapor recompression cycle, allowing to recover the latent heat by compressing the steam-hydrogen mixture prior to the condensation process, thus creating a temperature different between the hot exhaust and cold inlet streams. We show that this method can recover over 90 % of the latent heat, thus increasing the viability of solar thermochemical hydrogen production cycles even under limited conversion conditions.

hydrogen technology↗
DOE OSTI2024

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗
DOE OSTI2024

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hydrogen embrittlement of Zircaloy-4 fabricated by ultrasonic additive manufacturing

Ultrasonic additive manufacturing (UAM) was successfully applied to the zirconium material system to create a planar geometry. Following fabrication, SS-J3 type tensile specimens of the UAM and wrought Zircaloy-4 with a nominal gage section of 5×1.2×0.75 mm were machined for hydriding studies and mechanical testing. The SS-J3 type tensile specimens were gas-charged with hydrogen using a custom system that precisely controls hydrogen gas flow rate, hydrogen partial pressure, and temperature. To avoid altering the UAM material, the maximum process temperature was limited to 550°C. Using different initial hydrogen gas pressures and flow rates, various hydrogen contents (70–1755 wppm) were achieved for Zircaloy-4 specimens. Tensile testing shows that, regardless of hydrogen content, all UAM specimens measured yield strengths in the range of 557±16 MPa and ultimate tensile strengths of 660±4MPa. However, total elongation clearly decreased as a function of increasing hydrogen content. Further, at the lowest hydrogen content, the total plastic elongation measured 21.5%, and this value decreased to less than 2% when the hydrogen content was increased to 1000 wppm. Cross-sectional optical microscopy images revealed that the hydride distributions are randomly oriented. For the low hydrogen content specimen, these randomly distributed hydrides are isolated from each other. A sandwiched structure, consisting of high-density and low-density layers, was also observed for the specimens with hydrogen content between 200 and 400 wppm. As the hydrogen content increases, the hydrides diffuse into the low-density hydride layers to form a network across the whole specimen. Although the tensile properties of UAM and wrought Zircaloy-4 exhibit the same behavior with increasing hydrogen content, the distinctions in grain orientations led to differences in the hydride orientations at low and intermediate hydrogen concentrations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Type-II superconductivity at 9K in Pb–Bi alloy

In the present work, we report the synthesis of Pb–Bi alloy with enhanced Tc of up to 9K, which is higher than that of Pb. The alloy is synthesized via a solid-state reaction route in the vacuum-encapsulated quartz tube at 700°C in an automated furnace. The synthesized sample is characterized by X-ray Diffraction(XRD) and Energy dispersive X-ray analysis(EDAX) for its phase purity and elemental composition. Rietveld refinement of XRD reveals that the end product is a majority hexagonal Pb 7 Bi 3 , with minor rhombohedral Bi. The electronic transport measurement shows metallic behavior with the Debye temperature of 108K and a superconductivity transition temperature (T c ) below 9K, which is the maximum to date for any reported Pb–Bi alloy, Pb or Bi at ambient pressure. Partial substitution of Bi at the Pb site may modify the free density of electronic states within the BCS model to attain the optimum T c , which is higher by around 2K from the reported T c of Pb. The superconductor phase diagram derived from magneto-transport measurements reveals that the synthesized alloy is a conventional superconductor with an upper critical field (H c2 ) of 3.9 T, which lies well within the Pauli paramagnetic limit. The magnetization measurements carried out following ZFC(Zero Field Cool) protocols infer that the synthesized alloy is a bulk superconductor below 9K. The isothermal M-H(Magnetization vs. Field) measurements performed below T c establish it as a type-II superconductor. Furthermore, the specific heat capacity measurements show that the Pb–Bi alloy is a strongly coupled bulk superconductor below around 9K with possibly two superconducting gaps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

A Review on Direct Air Capture of Carbon Dioxide: Sorbent Materials, Process Engineering, Industrial Scale-Up, and Future Perspectives

The relentless accumulation of anthropogenic greenhouse gases has driven atmospheric carbon dioxide concentrations to approximately 426 ppm, necessitating the aggressive deployment of negative-emission technologies to achieve net zero by 2050. Direct air capture (DAC) offers a scalable, location-independent approach to atmospheric carbon removal; however, it is fundamentally constrained by the significant thermodynamic barriers associated with capturing CO 2 from ultradilute ambient conditions, requiring minimum thermodynamic energy inputs substantially higher than those for postcombustion point sources. This comprehensive review critically examines the technological landscape of DAC, focusing on the interdependent triad of sorbent material design, contactor engineering, and regeneration thermodynamics. We evaluate the fundamental boundaries of adsorption, emphasizing that an optimal adsorption enthalpy and isosteric heat of adsorption must balance the high CO 2 uptake capacity with the energetic penalties of sorbent regeneration. A systematic, comparative analysis of state-of-the-art sorbents is presented, encompassing mesoporous silicas, zeolites, carbon-based materials (CBMs), metal−organic frameworks (MOFs), porous organic polymers (POPs), and polymeric membranes. Special attention is devoted to surface functionalization strategies, particularly amine grafting and impregnation, which transition capture mechanisms from physisorption to chemisorption to enhance selectivity under ambient moisture and low partial pressures. Furthermore, we assess the operational merits of various reactor configurations, including gas−solid, gas−liquid, and membrane contactors, alongside regeneration cycles such as temperature, vacuum, pressure, and moisture swing adsorption. Finally, the review bridges fundamental materials science with industrial application by chronicling the scale-up milestones of pioneering entities and providing a strategic roadmap for advancing DAC technology readiness levels toward global deployment.

Adsorption↗
DOE OSTI2025

Methanogenesis and Acetogenesis in Hydrogenotrophy with Carbonate Minerals: Dependence on Mineral Surface Area, Biofilm Growth, and Microbial Community

The production, storage, and use of hydrogen are anticipated to grow substantially to achieve energy and climate goals. Consequently, microbial communities in many terrestrial and subsurface Earth environments could be exposed to elevated hydrogen concentrations. Hydrogen stimulates metabolic processes that reduce aqueous chemical species, such as bicarbonate or sulfate, that can exchange with solid mineral phases, but the controls on microbial hydrogenotrophy with mineral sources of electron acceptors are not fully understood. Herein, we applied laboratory experiments and biogeochemical modeling to study the response of a natural microbial community to an elevated partial pressure of hydrogen in the presence of carbonate minerals of varying composition, solubility, and size. Experimental incubations and simulation results showed that hydrogen consumption by microbial communities was initially dominated by sulfate reduction and, subsequently, transitioned to acetogenesis and methanogenesis. The rates of acetogenesis and methanogenesis were not correlated with the solubility of carbonate minerals. Instead, we observed strong linear correlations between the rates and surface area of carbonate minerals. Methane and acetate production slowed down in all incubations after about 2 weeks of incubation, although biogeochemical modeling predicted that the metabolic processes were not thermodynamically limited. Electron microscopy and infrared spectroscopy showed that biofilms with diverse microorganisms grew on the carbonates during this period. The methane δ13C value significantly increased, consistent with slower growth at elevated pH. This work highlights that microbial communities form biofilm on carbonate mineral surfaces as a response to hydrogen and that biofilm formation could pose a strong kinetic limitation to hydrogenotrophic metabolism utilizing carbonate minerals.

Acetobacterium↗
DOE OSTI2025

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2025

The Kinetics of the Coherent Stage of Hydroxide Corrosion on Li 2 O-Covered LiH

The capability to model LiOH growth on vacuum-baked LiH is a necessary precursor to making kinetic predictions of undesirable hydrogen outgassing from corroded LiH materials when placed in sealed-system applications. Here, in this work, LiH samples with initial LiOH surface corrosion layers less than 1.1 μm were outgassed by vacuum baking at high temperature to convert LiOH to Li 2 O. Diffuse-reflectance infrared Fourier transform (DRIFT) spectroscopy was used to measure the subsequent LiOH regrowth during exposure to 25–375 ppm of H 2 O vapor at room temperature. Logarithmic kinetics best model the growth of this coherent hydroxide corrosion layer (up to 1.5 μm) on previously vacuum-baked LiH samples. The logarithmic kinetics are attributed to the dissociation of H 2 O on the corrosion layer and the establishment of a H + /OH – electric field across the thickening LiOH corrosion layer. A time-dependent model of LiOH corrosion growth as a function of both initial LiOH thickness before vacuum baking and moisture partial pressures during re-exposure was developed to help researchers better assess the unwanted hydrogen outgassing potential from LiH/LiD materials.

Matt, Sarah M. [Lawrence Livermore National Labora↗
DOE OSTI2026

Pyrite (001) Interface Chemistry is Controlled by a Sulfoxy Termination

Pyrite (FeS 2 ) is the most common sulfide mineral on Earth, forming through inorganic reactions in the crust and oceanic hydrothermal systems and via microbially driven processes in anaerobic sediments. The pyrite–water interface is the site of a wide range of adsorption and reaction processes in Earth systems including oxidation that dramatically affects the geochemistry of surface waters and influences global carbon and oxygen cycles. Mechanistic geochemical models of pyrite interfacial reactivity, however, are limited by the lack of experimentally derived atomistic structures of the reduced and reacting surfaces. Here, in this work, we reveal the atomic-scale structure of the pyrite (001)-water interface that forms at very low oxygen partial pressures, relevant to suboxic environments in Earth. The interface structure and surface speciation were obtained using the crystal truncation rod method supported by ambient-pressure photoelectron spectroscopy and density functional theory calculations. The surface is dominantly composed of disulfide groups bound to a single oxygen atom, forming a sulfoxy group that has no known molecular or bulk mineral analog. This surface is interpreted as the first step in the oxidative dissolution of pyrite. The sulfoxy group is readily protonated through surface acid–base reactions that alter the structure of interfacial water and the free energy of interfacial reactions. Surface iron sites are not oxidized. Surprisingly, this interface can likely develop in equilibrium with bulk pyrite in some reducing and acidic solutions. This termination is therefore likely representative of pyrite surfaces under a vast range of experimental, industrial and Earth conditions.

oxidation↗
DOE OSTI2025

LASSO for CALPHAD Model Selection Enables Data-Efficient Thermodynamic Modeling: An Application in Thermochemical Hydrogen Production Materials

Phenomenological CALPHAD (CALculation of PHAse Diagrams) models, widely used for multicomponent materials, often contain a considerable number of parameters and require fitting using data from a relatively small number of experimental measurements or theoretical calculations. Sometimes these parameters are introduced for the purpose of improving model fits but without clear physical justification, which leads to overparametrized models with poor generalization performance. Automated approaches for optimal model selection based on the available data therefore become critical. Here, in this work, a least absolute shrinkage and selection operator (LASSO)-based approach is developed for model selection by leveraging the linearity of the CALPHAD model with respect to its parameters to convert the model selection and fitting to a LASSO minimization problem. We demonstrate its utility for thermodynamic modeling of thermochemical hydrogen (TCH) production materials using lanthanum strontium manganite (LSM) as an example. Various TCH-relevant properties, including oxygen stoichiometry as a function of oxygen partial pressure, enthalpy of reduction, and entropy of reduction, are successfully predicted with reasonable accuracy using a minimal set of model parameters. Importantly, the model selection and fitting involve minimal human decision; it can therefore be applied to high-throughput DFT defect calculations and yield efficient workflows for TCH material modeling and optimization.

CALPHAD↗