Search NASA⌕ Search
Home/Archive/DOE OSTI

DOE OSTI · 1774873

A comprehensive modeling method for proton exchange membrane electrolyzer development

Ma, Zhiwen·Witteman, Liam·Wrubel, Jacob A.·Bender, Guido
View original source ↗Read source document ↗Export metadata ↓

Abstract

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

Explore related subjects

08 HYDROGENelectrochemical modelinghydrogen productionlow temperature electrolysis water splittingproton exchange membrane electrolysis cell

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Ma, Zhiwen, Witteman, Liam, Wrubel, Jacob A., Bender, Guido. 2021-03-26. A comprehensive modeling method for proton exchange membrane electrolyzer development. https://doi.org/10.1016/j.ijhydene.2021.02.170

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

DOE OSTI2026

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN↗
DOE OSTI2026

SimH 2 : an integrated techno-economic modeling framework for hydrogen pipeline infrastructure and network optimization

Large-scale hydrogen (H 2 ) pipeline transport design and network optimization have seldom been reported due to the lack of a cost model accounting for the relationship between transport cost and hydrogen mass flow rate. Here, this work introduced a system-level cost model for hydrogen pipeline transport at supercritical state and integrated it with an existing CO 2 pipeline network tool, SimCCS, for hydrogen-specific pipeline design and optimization. The Intermountain West (I-West) region of the U.S., historically dependent on fossil fuel-based economies, is chosen to demonstrate the capabilities of our H 2 pipeline cost model and transport network optimization platform called SimH 2 . Two scenarios are examined: one where the pipeline is not allowed to pass through disadvantaged communities and the other where it is permitted. The results highlight that incorporating disadvantaged-community constraints lead to longer pipeline routes and increased transport costs, reflecting the trade-offs involved in equitable infrastructure development. It is demonstrated that the newly developed SimH 2 tool not only enables the efficient design of H 2 transportation pipelines but also optimizes the network by accounting for local terrain and the presence of disadvantaged areas.

08 HYDROGEN↗
DOE OSTI2026

Disordered hydrogen adsorption at the three-fold site of W(110)

Hydrogen often forms disordered phases on metals; however, deciphering the atomic-scale behavior requires experimental and modeling techniques that directly probe the short-range order between neighboring hydrogen adsorbates. Here, to address this challenge, we applied direct recoil spectroscopy (DRS) to investigate hydrogen adsorption on the W(110) surface. We show that the recoiled hydrogen flux measured during DRS is sensitive to the short range order of the adsorbed hydrogen. Ordered hydrogen neighbors are much more efficient at dechanneling incident ions along low-index surface channels, strongly suppressing the recoiled hydrogen flux along these directions. By modeling the DRS measurements with molecular dynamics for ordered and disordered hydrogen phases, we find that at a hydrogen surface coverage of approximately Θ = 0.6, the hydrogen adsorbates are bound to three-fold sites in a disordered phase at room temperature. This finding is consistent with transfer-matrix scaling model calculations that suggest the transition to an ordered hydrogen phase occurs at a greater surface coverage of Θ = 0.83.

08 HYDROGEN↗