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

The study, published in the journal Automotive and Engine Technology, addresses a long-standing trade-off in fuel cell simulation. Highly detailed stack models, which resolve the electrochemistry, heat transfer and fluid dynamics inside individual cells, are computationally expensive, making them impractical for full-vehicle or system-level investigations where thousands of operating seconds must be simulated. On the other end of the spectrum, reduced-order system models run quickly but often strip away the physical fidelity needed to capture how electrochemical processes respond to rapid changes in load, air supply and temperature. The new framework bridges this gap by running a physics-based stack model and a dynamic balance-of-plant model in parallel, exchanging data at every time step so that electrochemical, thermal and auxiliary losses are all accounted for within a single simulation environment.

Technically, the framework is implemented in AVL CRUISE™ M, a system simulation platform widely used in powertrain engineering, and MATLAB®, with the two tools synchronized during the simulation. The researchers evaluated their approach under representative transient operating conditions, the kind of rapidly fluctuating duty cycles that mirror real-world driving, including scenarios aligned with standard vehicle test cycles such as the Worldwide Harmonized Light Vehicles Test Cycle. A central contribution of the work is its systematic assessment of solver synchronization and coupling strategies, examining how the choice of data exchange scheme affects numerical stability, simulation accuracy and computational performance. This matters because co-simulation introduces a fundamental risk: if two solvers exchange information too infrequently or inconsistently, the coupled system can become unstable or produce physically meaningless results, particularly during sharp transients such as sudden acceleration demands.

The headline result is striking: the coupled framework achieved a real-time factor of 1.74, meaning the simulation ran faster than the physical process it represents. For automotive engineers, this figure is significant beyond mere convenience. Real-time capability is the entry ticket to hardware-in-the-loop testing, where a physical controller, such as an actual fuel cell management electronic control unit, is connected to the virtual plant and must respond to simulated signals within strict timing limits. It also opens the door to digital twin implementations, in which a running fuel cell system in a vehicle or test bench is mirrored by a live simulation model that can predict behavior, flag anomalies and support predictive maintenance. A framework that is both physically grounded and fast enough for real-time operation therefore has immediate practical value.

Equally important is what the framework reveals about where the energy actually goes. By integrating the stack model with the balance-of-plant model, the researchers performed an integrated loss analysis that links electrochemical stack losses with system-level inefficiencies. The results show that overall system efficiency declines from 49.6 percent at low load to 38.8 percent at high load, and that the relative contribution of individual loss mechanisms shifts dramatically with operating conditions. At low loads, activation losses associated with the sluggish kinetics of the oxygen reduction reaction at the cathode dominate, while at higher current densities ohmic losses from membrane resistance and concentration losses caused by mass transport limitations grow in importance. Simultaneously, the auxiliary components consume an increasing share of the produced power, since the air compressor must work harder to deliver the oxygen the stack demands at high current.

This load-dependent loss structure has direct implications for how fuel cell vehicles are designed and operated. A fuel cell powertrain that is efficient during gentle city driving may shed several percentage points of efficiency during highway cruising or hill climbing, not because the electrochemistry itself becomes dramatically worse, but because the balance of plant consumes more and the stack crosses into loss-dominated regimes. With a validated co-simulation environment, engineers can quantify these trade-offs before a single physical prototype is built. They can explore, for example, how a more efficient compressor, a revised humidification strategy or an altered operating temperature might shift the efficiency map, and how those choices interact with transient behavior during real driving profiles rather than idealized steady-state points.

The modular architecture of the framework is another distinguishing feature. Unlike existing co-simulation approaches that focus primarily on stack-internal behavior, the new design treats the stack and the balance of plant as interchangeable modules with well-defined interfaces. This means individual components can be swapped, refined or replaced without rebuilding the entire simulation. A researcher could substitute a more detailed compressor model, test a different membrane-electrode-assembly characterization, or plug in an alternative cooling concept while keeping the rest of the environment intact. Such flexibility is essential in industrial practice, where simulation tools and model detail levels evolve continuously across development programs and where different teams may be responsible for different subsystems.

The broader context makes this work timely. As automakers ramp up hydrogen programs for heavy-duty trucks, buses and passenger vehicles, the pressure to compress development cycles while meeting ever-stricter efficiency and durability targets has intensified. Fuel cell durability is highly sensitive to operating conditions: water management, membrane humidity, temperature gradients and load transients all influence degradation mechanisms. A system-level simulation environment that faithfully represents the interplay between the stack and its auxiliaries under realistic transient loads provides a virtual proving ground for control strategies that balance performance against long-term reliability. The authors specifically point to control strategy development, hardware-in-the-loop testing and digital twin implementations as the natural next applications of their framework.

In an era when the hydrogen economy is alternately hyped and dismissed, studies like this one do the unglamorous but essential work of turning chemistry into drivable engineering. By demonstrating that a physics-based stack model and a dynamic balance-of-plant model can be coupled stably, accurately and fast enough for real-time use, the Graz and Darmstadt researchers have given the fuel cell community a practical tool for answering one of its most consequential questions: not just how efficiently a membrane can convert hydrogen into electricity in isolation, but how an entire system behaves when the road demands everything at once. With system efficiency spanning nearly eleven percentage points between low and high load, and with loss mechanisms shifting beneath every change of the accelerator pedal, that system-level view may prove to be exactly what the next generation of hydrogen vehicles needs.

Beyond the headline performance figures, the publication carries several details that will interest specialists in simulation methodology. The framework’s modular design is consistent with the Functional Mock-Up Unit (FMU) standard, an industry convention that allows models developed in different tools to be packaged with standardized interfaces and exchanged between simulation environments. Adherence to such conventions is what makes a co-simulation genuinely cross-tool rather than a bespoke pairing of two specific software packages, and it explains how a physics-based stack model developed in MATLAB can coexist with a commercial powertrain platform like AVL CRUISE M without forcing either side to abandon its native solver.

The abbreviation list of the paper also hints at the physical depth retained inside the stack module. Terms such as area-specific resistance, gas diffusion layer, bipolar plate, and membrane-electrode-assembly indicate that the model resolves the layered structure of an individual cell rather than treating the stack as a black box. Likewise, the explicit distinction between the hydrogen oxidation reaction at the anode and the oxygen reduction reaction at the cathode, alongside side reactions such as hydrogen evolution, suggests that the electrochemical submodel distinguishes loss pathways that simplified models typically lump together.

Another noteworthy aspect is the choice of reference for efficiency. The reported values are given on the basis of the higher heating value, the more conservative convention that accounts for the latent heat of water vaporization. Efficiencies quoted on a higher heating value basis are systematically lower than those quoted on the lower heating value basis often used in combustion engineering, which is worth remembering when comparing the reported 49.6 percent peak against figures from other fuel cell studies.

The institutional collaboration itself is instructive. The work emerged from the Institute of Thermodynamics and Sustainable Propulsion Systems in Graz and the Institute for Internal Combustion Engines and Powertrain Systems in Darmstadt, both groups with deep roots in conventional powertrain simulation now applying that heritage to hydrogen systems. Published open access with funding from Graz University of Technology, and having passed through peer review between May and August 2026, the study is freely available to the development community it aims to serve.

Subject of Research: Integrated co-simulation of PEM fuel cell systems coupling electrochemical stack modeling with balance-of-plant dynamics for automotive applications

Article Title: Integrated co-simulation of PEM fuel cell systems for automotive applications: modeling, coupling, and system behavior

Article References: Strametz, F., & Endres, J. (2026). Integrated co-simulation of PEM fuel cell systems for automotive applications: modeling, coupling, and system behavior. Automotive and Engine Technology. https://doi.org/10.1007/s41104-026-00177-x

Image Credits: AI Generated

DOI: 10.1007/s41104-026-00177-x

Keywords: PEM fuel cell, co-simulation, balance of plant, fuel cell stack modeling, system-level simulation, real-time simulation, loss analysis, AVL CRUISE M, MATLAB, automotive powertrain, digital twin, hardware-in-the-loop