Wavy, Tapered Channels and an AI Ensemble Push Fuel Cells to Record Power

Hydrogen fuel cells promise clean electricity with nothing but water as exhaust, yet their real-world performance has long been throttled by an unglamorous component: the network of channels etched into the plates that feed hydrogen and oxygen to the reaction sites. A new study published in Results in Engineering by Mehrdad Ghasabehi and Mehrzad Shams reports a strikingly effective redesign of that plumbing, combining tapered main channels with wavy sub-channels in a proton exchange membrane fuel cell (PEMFC). Using a physics-informed surrogate model and an ensemble of five optimization algorithms, the researchers arrived at a configuration delivering a peak power density of 0.913 watts per square centimeter, a 45.6 percent improvement over a conventional parallel flow field and a 17.9 percent gain over a tapered design without waves, all while keeping the energy spent on pumping reactants to a mere 1.86 percent of the cell’s output.

The core problem the team attacked is flow maldistribution. In a classic parallel flow field, each channel experiences a different hydraulic resistance, so gas preferentially rushes through the outer channels while the central ones starve. The consequences cascade through the entire cell: oxygen-deprived dead zones form in the middle of the active area, local current density collapses, and water droplets accumulate in the sluggish regions, partially flooding the porous gas diffusion layer. In the conventional design simulated by the researchers, dead zones showed oxygen concentrations more than 30 percent below the average, and localized hotspots reached volumetric heat generation rates of 3.6 times 10 to the ninth watts per cubic meter, roughly 2.6 times the catalyst layer average.

The modified design solves this with two cooperating mechanisms. First, the main channels taper: they are wide at both the inlet and outlet manifolds and progressively vary in width between them, defined by a tapering ratio that compares the width difference to the end width. Counterintuitively, a larger tapering ratio produces a lower pressure profile along the main channels, consistent with Bernoulli’s principle, because the cross-section does not simply shrink. Tapering also rearranges the sub-channel connections, redistributing hydraulic losses and suppressing the preferential edge-channel flow. In the simulations, mid-channel flow speeds rose from 0.13 to 0.77 meters per second, and velocities across different channels became far more uniform, reducing central flow starvation by more than a factor of four.

Second, the sub-channel walls are shaped into smooth, connected circular arcs that repeat every wavelength, creating a converging-diverging wave pattern with no sharp corners. As gas flows through these undulations, centrifugal forces drive it outward and generate pairs of stationary, counter-rotating vortices, akin to the Dean vortices familiar from curved pipe flow. The enhanced design exhibited vortex pairs with a vorticity magnitude of 3000 per second inside every sub-channel, whereas the simple parallel design showed almost none. These vortices perform double duty: they mix fresh oxygen from the channel core toward the gas diffusion layer, and they sweep along the layer surface, repeatedly stripping away the stagnant boundary layer that would otherwise throttle diffusion.

The payoff is visible in the species and water maps. Oxygen concentration at the gas diffusion layer-catalyst layer interface rose 38 percent, from 0.0045 to 0.0062 kilomoles per cubic meter, and the spatial variation dropped from more than 30 percent below average in the dead zones to under 5 percent. Under limiting-current conditions, where concentration losses dominate, the limiting current density climbed 40 percent, from 1.5 to 2.1 amperes per square centimeter. Water saturation fell from 0.13 to 0.11 while membrane water content stayed constant, because higher local gas speeds and vortex shear detach droplets and push them downstream, and the diverging wave sections provide sudden pressure drops that aid water removal without demanding much extra compressor power.

Simulating every conceivable design with a full three-dimensional, multiphase computational fluid dynamics model would be prohibitively expensive, so the researchers built a surrogate. They ran a Central Composite Design of 81 simulation cases, spanning five operating variables: anode and cathode stoichiometries, pressure, wave amplitude, and tapering ratio, each evaluated at three voltages of 0.3, 0.5, and 0.7 volts, plus 20 random unseen cases for testing. The key innovation is the Physics-Informed Response Surface Methodology, or PIRSM. Conventional response surfaces fit a single low-order polynomial to the whole polarization curve, but a fuel cell’s voltage-current relationship spans three physically distinct regimes, activation, ohmic, and concentration losses, and a lone polynomial cannot capture them all. In one demonstration, a conventional fit deviated by 17.7 percent at open-circuit voltage.

PIRSM instead builds three separate second-order response surfaces, each trained only on data from one loss regime, then aggregates them into a third-order polynomial that transitions smoothly between regions and correctly predicts zero output at open-circuit voltage. For the parasitic power ratio, a cube-root transformation motivated by the nonlinear scaling of pressure drop with velocity prevents non-physical predictions such as negative pressures. The resulting models achieved R-squared values of 0.99 against the simulation data for both objectives, with root mean square errors of 0.017 watts per square centimeter for power density and 0.005 for the parasitic ratio, and an R-squared of roughly 0.96 on the held-out test set. The underlying CFD model itself was validated against experimental polarization data with a root mean square error of 0.0222 amperes per square centimeter, and a grid convergence index of 0.58 percent confirmed numerical reliability.

On top of the surrogate sits the second major novelty: an ensemble optimization framework. Rather than relying on a single metaheuristic, the team ran five multi-objective algorithms from distinct families, three evolutionary methods (NSGA-II, SPEA II, and MOEA/D) and two swarm-based methods (multi-objective grey wolf optimization and multi-objective particle swarm optimization). Their Pareto fronts, each reflecting a different search bias, were pooled and filtered by a dominance criterion, in which one solution dominates another only if it is no worse on every objective and strictly better on at least one. Statistical testing over 30 runs per algorithm, using pairwise Wilcoxon signed-rank tests with Bonferroni correction, established a clear hierarchy: NSGA-II was the superior performer, winning 86.7 percent of rank-1 placements with corrected p-values below 0.001, while SPEA II occupied a consistent second tier. The ensemble’s pooled front captured a wider range of trade-offs than any single algorithm, and because the algorithms run in parallel, the wall-clock cost is the maximum of individual run times rather than their sum.

The winning configuration, selected from the ensemble Pareto front using the TOPSIS method, features a tapering ratio of 2.23, a wave amplitude of 0.47 millimeters, a pressure of 4 atmospheres, stoichiometries of 3 on both anode and cathode, and an operating voltage of 0.47 volts. It achieves an oxygen concentration of 0.01396 kilomoles per cubic meter at the catalyst layer interface, roughly 210 percent higher than the simple parallel design, and an oxygen uniformity index of 0.8174, up 26.5 percent. The limiting current density reaches 2.35 amperes per square centimeter, a 56.7 percent increase over the baseline, and the pressure drop falls by a factor of 6.5 relative to the simple parallel layout. The parasitic power ratio of 0.0186 represents an 84.5 percent reduction compared with the conventional design, although the simpler tapered variant without waves retains slightly lower pumping losses.

The researchers caution that their predictions rest on simulations within the tested parameter ranges, and that experimental validation of the waved-tapered geometry remains an essential next step. The smooth arcs of the wave pattern avoid sharp corners and should benefit mechanical integrity, but repeated converging-diverging sections may introduce cyclic stresses, so accelerated stress tests, thermal cycling, and degradation experiments are recommended to confirm durability. Even so, the study demonstrates a compelling template for fuel cell engineering: couple geometric innovations that exploit real flow physics with surrogate models constrained by electrochemistry, and let a diverse committee of optimizers, rather than any single algorithm, chart the trade-off frontier. As hydrogen economies scale up, designs like this one, squeezing more power from every square centimeter while spending almost nothing on parasitic pumping, could help make fuel cells not just clean, but genuinely competitive.

Subject of Research: Multi-objective optimization of a tapered parallel flow field with wavy sub-channels for proton exchange membrane fuel cells using physics-informed response surface methodology

Article Title: Ensemble multi-objective optimization of modified tapered parallel flow field PEM fuel cell using physics-informed response surface methodology

Article References: Ghasabehi, M., & Shams, M. (2026). Ensemble multi-objective optimization of modified tapered parallel flow field PEM fuel cell using physics-informed response surface methodology. Results in Engineering, 32, Article 113366. https://doi.org/10.1016/j.rineng.2026.113366

Image Credits: AI Generated

DOI: Not provided

Keywords: PEM fuel cell, flow field design, tapered channels, wavy sub-channels, response surface methodology, multi-objective optimization, NSGA-II, Pareto front, computational fluid dynamics, power density, parasitic losses, hydrogen energy