Review Article

State of the Art of Modelling and Design Approaches for Ejectors in Proton Exchange Membrane Fuel Cell

Table 1

Literature survey on ejector simulation and geometric optimization.

Ref.YearTypeMore detailsKey findings

[38]2019NumericalANSYS Fluent, 2D axisymmetric steady flow, SST, 1-phase 1-speciesIncreasing raises ER, after a certain value performance deteriorates
[39]2020NumericalThermodynamic model, annular mixing layerOptimal depends on working condition
[40]2022ExperimentalEjectors with different nozzle and area ratios tested over wide ranges of operating conditionsUp to 4 “optimal” observed; these values were affected by both working condition and area ratios
[41]2022MixExperimental testing and validation of proposed model, ANSYS Fluent, 2D axisymmetric steady flow, SST, 1-phase 1-species, compound choking criterionCompound choking allowed to determine variation in secondary flow choking position under several suboptimal
[42]2019NumericalANSYS Fluent, 2D steady flow, ; model coupled with anodic pressure drop formula, 1-phase 2-speciesOptimal range for (3.00–3.54) and (1–3)
[43]2021Numerical3D steady flow, RNG , 1-phase 2-species is the first geometrical parameter to be optimized
[51]2014NumericalThermodynamic model of ejector coupled with semiempirical stack modelDefinition of two dimensionless parameters to guide ejector design
[52]2017ExperimentalEjector designed and tested at constant load and in fast transient conditionAnode gas recirculation rate ranging from 40% fuel utilization per pass at 25 A stack current to 64% fuel utilization per pass at 160 A stack current
[53]2022NumericalANSYS Fluent, 2D steady flow, SST, 1-phase 2-speciesDetermined order of influence of geometrical parameter:
(i) Low current (110 A)

(ii) Middle current (275 A)

(iii) High current (412.5 A)

Nozzle throat length
[54]2013NumericalThermodynamic modelInlet primary flow temperature affects ejector entrainment ratio and component efficiency
[55]2019MixANSYS Fluent, 2D axisymmetric steady flow, comparison between RNG and SSTRNG model shows higher accuracy than SST; optimal ° and  mm when stack works at its rated power
[56]2020NumericalANSYS Fluent, 2D axisymmetric steady flow, SST, coupled with a pressure drop through anode modelER can be influenced by anode inlet temperature, relative humidity, and differential pressure
[57]2020NumericalOpenFOAM, 3D transient flow, RNG , 2-phase 3-speciesDynamic responses during power variations results from velocity differences between the primary and the secondary flow; increase of nitrogen mass fraction promotes total ER, while it reduces hydrogen ER
[58]2020NumericalCOMSOL, 3D steady flow, coupled with MATLAB/Simulink hydrogen recovery system modelA lower ejector temperature is disadvantageous in removing the moisture content of the recirculated hydrogen gas, thus in practical applications; the hydrogen inlet temperature/pressure must be carefully controlled
[59]2022NumericalThermodynamic model of a 2-phase CPM ejectorER increase from 0.47 to 1.14 as mixing area ratios range from 1.0 to 1.2 under the given conditions
[60]2015NumericalIntegrated lumped parameter-CFD approach, ANSYS Fluent, 2D axisymmetric steady flow, SSTThe model can be used for studying off-design conditions, where ejector component efficiencies are not constant
[61]2016NumericalANSYS Fluent, 2D axisymmetric steady flow, RNG , 2-phase 1-species, optimization through genetic and evolutionary algorithm is the crucial parameter in ejector performance
[62]2021NumericalANSYS Fluent, 2D steady flow, RNG , 1-phase 2-speciesHumidity and temperature of the secondary flow have a noticeable influence on the performance of the ejector
[44]2017NumericalThermodynamic model, hybrid fish swarm algorithmOptimization efficiency increased with respect to genetic algorithm
[45]2018NumericalMultiobjective evolutionary algorithm coupled with a surrogate model based on CFD simulations and are the most important geometrical variable; entrainment ratio can be increased up to 110% and 35%, for air and CO2, respectively
[46]2023NumericalANSYS Fluent, 2D axisymmetric steady flow, response surface methodologyPriority order in optimizing ejector geometry:
[47]2014Numerical2D axisymmetric steady flow, RNG , 1-phase 2-species, artificial neural network and genetic algorithm to obtain optimal geometryOptimal ; optimal ; optimal °; optimal
[48]2021NumericalAutomated CFD workflow, ANSYS Fluent, RNG , 2-phase 1-species, Gaussian process regression machine learning modelThe algorithm can be used to efficiently explore ejector designs with mean average errors between 0.07 and 0.1
[49]2022NumericalANSYS Fluent, 2D steady flow, realizable , optimization via adjoint methodER increased by around 37%
[50]2022NumericalMATLAB and experimental dataset of a steam-centered ejector are applied to train the ANN model of a steam ejector using three different algorithmsLM model yielded the best agreement; the effect of the outlet area ratio is less important with respect to throat area ratio