^{1}

^{*}

^{1}

^{1}

^{1}

^{1}

^{2}P-EnergyLab, EA 4079, University of La Reunion, Saint Denis, France

This paper proposes a parameterized nonlinear model-based predictive control (NMPC) strategy to tackle the oxygen excess ratio regulation challenge of a proton exchange membrane fuel cell. In practice, the most challenging part regarding NMPC strategies remains the on-line implementation. In fact, NMPC strategies, at least in their basic form, involve heavy computation to solve the optimization problem. In this work, a specific parameterization of control actions has been designed to address this limitation and achieve on-line implementation. To assess the effectiveness and relevance of the proposed strategy, the controller has been implemented on-line, experimentally validated on a real fuel cell and compared to the built-in controller. Performance of the parameterized NMPC controller in terms of setpoint tracking accuracy, disturbances rejection and computational cost, have tested under several control scenarios. Experimental results have shown the excellent tracking capability, disturbances rejection ability and low computational cost of the NMPC controller, regardless of the operating conditions. Moreover, compared to the built-in controller the proposed strategy has demonstrated better disturbances rejection capability. Overall, the proposed parameterized NMPC controller appears as an excellent candidate to address the oxygen excess ratio regulation issue.

For the last decade, to reduce greenhouse gas emissions and fossil fuel dependence, numerous renewable energy technologies have been studied. In regard to renewable energies, one of the primary drawbacks is the variability of the supply flows, which raises the key issue of energy storage to counteract the intermittent nature of their conversion. In this context, the hydrogen vector represents a promising alternative, as long as the hydrogen is produced from a renewable energy technology (e.g. electrolytic hydrogen obtained using photovoltaic energy), and stored to optimize potential gaps and surplus of intermittent production. Downstream the hydrogen chain, electricity is generated from fuel cells, genuine zero-emission power generators. Due to their high power density and low operating temperature, proton exchange membrane fuel cells (PEMFC) have proved to be the most suitable fuel cell technology for both transportation and stationary applications [

However, several bolts still remain to be removed to improve their reliability and energy conversion, reduce their cost or extend their lifetime. Among them, one of the most important is related to their control. Indeed, global efficiency improvement, optimal hydrogen and air consumption, and reliable and accurate power response remain challenging control goals.

Numerous control strategies have been reported in the literature for PEMFC systems, ranging from PID controllers [

Several works dealing with benefits of predictive control strategies regarding PEMFC control have been reported. Wu et al. [

Regarding nonlinear model-based predictive control strategies, the modeling step remains essential and mandatory for the control design. In the literature, a wide variety of models, either mechanistic or black-box, with different levels of complexity, can be consulted. Regarding mechanistic approaches, numerous one-dimensional and multi-dimensional models have been proposed [

To ensure optimal performance of PEMFC, one parameter, namely oxygen excess ratio, requires special attention. The oxygen excess ratio or stoichiometric ratio represents the ratio of inlet oxygen flow to reacted oxygen flow and is widely used to guarantee safety and to reach a high performance. Arce et al. [

In this paper, a real-time implementable nonlinear model-based predictive control (NMPC) strategy is developed to tackle the oxygen excess ratio regulation issue. In practice, the most challenging part regarding NMPC strategies remains the on-line implementation. In fact, NMPC strategies, at least in their basic form, involve heavy computation to solve the optimization problem. In this work, a specific parameterization of control actions is designed to address this limitation and achieve on-line implementation. This parameterization allows to significantly decrease the optimization problem dimension and so the computational time required to solve the optimization problem. Regarding real-time control goal, and according to its short computational time compared to classical mechanistic model and its low sensitivity to noise, an artificial neural network (ANN) model is used as internal model to predict the PEMFC behavior.

The rest of this paper is organized as follows: The experimental setup of the fuel cell unit is presented in Section 2. Section 3 is dedicated to the control design—the ANN model predictor is designed and experimentally validated—and the real-time implementable NMPC scheme dedicated to PEMFC oxygen excess ratio control is detailed. Experimental validation of the proposed control strategy is presented in Section 4. Eventually, conclusions and prospects are drawn in Section 5.

In the present study, a 50 cm^{2} single home fuel cell assembling, with a commercial membrane electrode assembly (MEA) from Paxitech, is used. The fuel cell is assembled with Teflon glass fiber as seal joints and clamped at 11 Nm bolt torque. MEA used for this work is composed of symmetric electrodes (0.5 mgPt·cm^{−2}, Pt 40 wt%/C) with a Nafion® N115 membrane. Air and pure hydrogen are used as oxidant and fuel, respectively. All experiments are performed on a Fuel Cell Test Station (FCTS) manufactured by Fuel Cell Technologies, Inc.

The FCTS provides control of temperature, flow rate, pressure and humidity for both feed gases. The bench integrates an AC-Impedance measurement system and a programmable DC electronic load. Inlet gas flow rates are controlled using two mass flow controllers, located before the humidifier. Humidification and fuel cell temperatures are controlled using temperature controllers. Reactant gases are humidified through external water tanks, and desired gas humidification value is reached through the regulation of the water temperature in the humidifiers. Output pressures are regulated using backpressure controllers designed to operate between 0 and 400 kPa. An electronic load measures cell output voltage and current (

The Supervisory Control and Data Acquisition (SCDA) system relies on a proprietary graphical user interface composed of a set of modular virtual instruments (VI) developed with LabView^{TM} graphical software. This makes the experimental unit a flexible environment, allowing to perform a wide range of tests based on various configuration settings. Since VI’s can easily be modified and/or replaced, alternative monitoring and control strategies can be integrated to the SCDA system. As regards to innovative control strategies, a set of VIs makes it possible to replace the original built-in controller by model-based approaches.

Regarding NMPC control strategies, the first step is to design a model of the system. This model is expected to predict the system behavior several steps ahead. Among the various PEMFC models that can be consulted in the literature, a very few are dedicated to control purposes. In real-time control context, most of complex and heavy computations mechanistic models cannot be considered. Recently, several works emphasize the interest of ANN to model PEMFC systems. Here, due to its short computational time, its low sensitivity to noise and its reliability an ANN model is designed and used as predictive model.

The proposed control strategy involves an ANN model of the PEMFC used to predict the system output several steps ahead.

The oxygen excess ratio, namely λ O 2 , denotes the ratio of oxygen supply to oxygen consumed in the cathode. λ O 2 is a function of the inlet air flow rate at the cathode m ˙ a i r , the load current I, the pressure P c a , i n , the temperature T c a , i n , the relative humidity ϕ c a , i n , the number n of cells of the fuel cell stack, the Faraday constant F, the oxygen mole fraction y O 2 , c a , i n and the molar masses M O 2 , M N 2 and M v of oxygen, nitrogen and vapor, respectively.

λ O 2 = m ˙ a i r I f ( P c a , i n , T c a , i n , ϕ c a , i n , n , F , y O 2 , c a , i n , M O 2 , M N 2 , M v ) (1)

However, del Real et al. [

λ O 2 = f ( m ˙ a i r , I ) (2)

In this paper, an ANN model is designed to predict λ O 2 several steps ahead, once the load current I and the inlet air flow rate m ˙ a i r are available. In practice, the load current is a measured disturbance, whereas the inlet air flow rate is the manipulated variable.

Regarding real-time implementation goal, the computational speed of the model is a key component. The model is expected to predict the system output with sufficient accuracy, while avoiding heavy computations. To address this problem, several ANN architectures have been investigated to determine which one provides the best trade-off between prediction accuracy and computational speed. Finally, a fully connected three layers network is designed. The hidden layer has five neurons with a tangent sigmoid activation function φ 1 , whereas the output layer has one neuron using a linear activation function φ 2 .

Identification (training phase) and validation of the ANN model are performed on two different sets of data collected on the PEMFC system presented in Section 2 and sampled with a 3 seconds period. Regarding the poor extrapolation performance of ANN models, the data set used for the training phase has to be chosen wisely. In this aim, the training data set, used to identify the ANN model, covers the whole operating conditions and ranges from minimum to maximum inputs/outputs values. Experiments are performed for different values of load current (ranging from 5A to 20A) and inlet air flow rate (ranging from 500 sccm to 2800 sccm), which leads to oxygen excess ratio values ranging from 5 to 11.

The adjustable parameters of the ANN model (so-called weights) are specified by the vector θ and are determined from the training data set Z N :

Z N = { [ u ( t ) , y ( t ) ] | t = 1 , ⋯ , N } (3)

where u ( t ) and y ( t ) are the set of inputs and the corresponding desired outputs, respectively.

The training phase aims to determine a mapping from the set of training data to the set of possible weights, so that the ANN model produces predictions y ^ ( t ) :

Z N → θ ^ (4)

The prediction error approach, which is the strategy applied here, is based on the introduction of a measure of closeness in terms of a mean square error criterion:

V N ( θ , Z N ) = 1 2 N ∑ t = 1 N [ y ( t ) − y ^ ( t | θ ) ] T [ y ( t ) − y ^ ( t | θ ) ] (5)

The weights are then found as:

θ ^ = arg min θ V N ( θ , Z N ) (6)

The training phase is performed using Levenberg-Marquart algorithm with an error goal of 0.001. The validation results are presented in

Two criteria, namely root mean square error ( RMSE = ∑ i = 1 n ( y i − y ^ i ) 2 / n ) and absolute maximal error ( AME = max i ( | y i − y ^ i | ) ), are considered to assess the one step ahead prediction performance of the ANN model. n represents the number of data, y and y ^ denote experimental and simulated data, respectively.

With an RMSE of 0.003 and an AME under 0.1 the predictive performance accuracy of the ANN model is more than sufficient, especially since it is dedicated to control purposes (

In this study, a parameterized NMPC controller is used to control the PEMFC oxygen excess ratio ( λ O 2 ), once the manipulated variable is the air inlet flow rate ( m ˙ a i r ) and the load current (I) is a measured disturbance.

NMPC is an optimization-based multivariable constrained control method that uses a nonlinear model to predict the future behavior of the process. Classical

Training | Validation | |
---|---|---|

RMSE | 4.726 × 10^{−4} | 0.003 |

AME | 0.079 | 0.98 |

NMPC is defined by solving on-line an optimization problem that leads to an optimal sequence of future control actions that minimizes a given cost function over a certain prediction horizon. NMPC offers a great flexibility since it can handle input/state constraints, process nonlinearities, and even economical or operating constraints [

In this work, to address this problem and reach real-time control goal, a parameterized NMPC strategy is designed. When dealing with classical NMPC strategies, the predictive horizon sets the optimization problem dimension, which could lead to huge problem dimension. The main feature of parameterized NMPC strategy consists in decreasing the optimization problem dimension using a particular parameterization of the control sequence u:

u : = ( u ( p ) ( k + 1 ) , ⋯ , u ( p ) ( k + N y ) ) (7)

where u ( p ) ( k + 1 ) is the control action at time k + 1, N_{y} the predictive horizon and p a set of parameters.

p ^ = arg min p ( J ( u ( p ) ) ) (8)

where the optimal solution p ^ is the set of parameters that minimizes the cost function J while meeting the problem constraints.

In this context, as the problem dimension depends on the dimension of p, itis straightforwardly linked to the parameterization choice and decoupled from the predictive horizon N_{y}. In this manner, a well-structured parameterization scheme may significantly decrease the optimization problem dimension, which is a quite important feature to fulfilreal-time requirements [

Here, since the system is open-loop stable, an approximation of the steady control action is used in the parameterization of the control profile to decrease the degree of freedom of the control profile.

In the first step, several experiments with the fuel cell in open loop are carried out. Even if λ O 2 depends on several variables, Gruber et al. [

In a second step, these data are gathered to generate a stationary map correlating the load current, the inlet air flow rate and the oxygen excess ratio.

In a last step, based on this stationary map a simplified model is identified to approximate the steady state control action as a function of the oxygen excess ration and the load current:

u * = a 1 λ O 2 I (9)

a 1 is a constant parameter that has to be determined off-line.

In this work, a specific parameterization of the control actions based on the approximation of the steady state control action u * and a single scalar parameter p is proposed:

u ( p ) ( k + i ) = u ∗ ( k + i ) + p E ∗ ( k + i ) for i ∈ { 1 , ⋯ , N y } (10)

p ∈ [ − 1 , 1 ] is a scalar parameter, u ∗ ( k + i ) = a 1 y s p ( k + i ) I ( k + i ) is the approximation of the steady control action obtained in stationary steady regime for λ O 2 ( k + i ) and I ( k + i ) , and E ∗ ( k + i ) is the trajectory tracking error obtained using u ∗ ( k + i ) . Specifically, E ∗ ( k + i ) = y s p ( k + i ) − y ^ A N N ∗ ( k + i ) and y ^ A N N ∗ ( k + i ) is the ANN model output using u ∗ ( k + i ) and I ( k + i ) as inputs.

In this manner, the control action u ( p ) ( k + i ) can be seen as a dynamic correction of the approximation of the control action in steady state. In the present case, this is a simple but efficient way to significantly reduce the degree of freedom of the control action profile and so to reduce the computation effort required to solve the optimization problem.

Here the cost function is written to minimize the setpoint tracking error and to limit the variations of the control variable:

J ( u ( p ) ) = ∑ i = 1 N y γ ( y s p ( k + i ) − ( y ^ A N N ( k + i ) + e m ( k ) ) ) 2 + ∑ i = 1 N u β ( u ( p ) ( k + i ) − u ( p ) ( k + i − 1 ) ) 2 (11)

y s p ( k + i ) and y ^ A N N ( k + i ) denote the desired value of oxygen excess ratio and the oxygen excess ratio estimated by ANN model, respectively. u ( p ) ( k + i ) = m ˙ a i r ( p ) ( k + i ) is the control action at time k + i, and γ and β are two weighting parameters. N_{y} and N_{u} are the control horizon and the predictive horizon, respectively. The prediction horizon corresponds to the future time interval used to compute the process output predictions with the process model. The control horizon corresponds to the time interval when present and future control actions are computed. The control horizon is shorter or equal to the prediction horizon. For time intervals between N u and N y the control actions are constant and set equal to u ( p ) ( k + N u ) .

Note that to improve the robustness of the approach against plaint/model mismatch, the error modeling e m ( k ) = y ^ A N N ( k ) − y ( k ) at time k is used over the predictive horizon to correct the ANN model outputs y ^ A N N ( k + i ) . y ( k ) represents the oxygen excess ratio measured at time k. Indeed, the cell performance changes along its life, which necessary involves plant/model mismatch. Therefore, to assure an optimal control of the PEMFC along its life, the controller robustness against plant/model mismatch is crucial.

In this context, the optimal value p ^ is obtained by minimizing the cost function J, which is subject to constraints on manipulated variable:

p ^ = arg min p ( J ( u ( p ) ) ) (12)

Subject to

u ( p ) ( k + i ) ∈ [ u min , u max ] for i = 1 , ⋯ , N y (13)

Here, the optimization problem that lead to the sequence of future control actions relies entirely on a single scalar parameter, which dramatically reduces the computational effort compared to a classical NMPC strategy [

To assess the performance of the proposed control strategy in terms of tracking capability, disturbances rejection, and robustness against plant-model mismatch, series of experiments are performed on the PEMFC system presented in Section 2. In this aim, several VI’s of the SCDA system are modified to replace the built-in controller by the parameterized NMPC controller. With this setup, the optimal control action, computed on-line in Matlab® environment, is sent through a specific VI to be applied to the experimental unit.

Here, to exemplify the controller performance three control scenarios are considered: The first scenario illustrates the tracking capability of the controller for a constant value of load current. The second scenario is designed to study the controller ability to cope with disturbances. Eventually, the third scenario evaluates the controller performance in terms of setpoint tracking accuracy in presence of disturbances. A comparison between the built-in controller and the proposed control strategy is made for the second scenario. Note that this comparison cannot be performed for the first and third scenario due to technical constraints. In fact, the built-in controller of the SCDA system, integrated by Fuel Cell Technologies, is designed to operate at constant oxygen excess ratio. In other word, when using the original setup, the oxygen excess ratio value has to be set before the experiments and cannot be modified during the experiments. The built-in controller is a Proportional Integral (PI) controller. However, the SCDA system provided by Fuel Cell Technologies, Inc does not allow to access the controller parameters or to modify them.

In the sequel, the inlet air flow rate is taken as manipulated variable, whereas the load current is a measured disturbance. It is assumed that all other variables, required to ensure safety operation of the fuel cell, are properly controlled. The cell temperature and the sampling period are set to 75˚C and 3 s respectively. The optimization problem is solved on-line using Levenberg-Marquardt algorithm. All experiments are performed using control parameters listed in

Parameter | Value |
---|---|

a_{1} | 16.57 |

N y | 30 |

N u | 30 |

γ | 1 |

β | 1 |

The tracking capability of the controller is evaluated using an oxygen excess ratio setpoint that covers the whole operating condition. In this control scenario, the load current is set to 15A and kept constant during the entire experiments.

To assess the controller performance in terms of disturbances rejection, a set of step changes on the load current is performed. The dynamics of the PEMFC is strongly correlated to the power level. Therefore, to verify that the controller performs accurately whatever the power level, the current steps is chosen to cover the whole operating conditions. Several experiments, with different sets of step changes in load current and different oxygen excess ration values, are carried out to compare the performance of the original built-in controller and the proposed NMPC strategy.

In each and every case, the NMPC controller performed better than the built-in controller.

The NMPC controller demonstrates that it can maintain the oxygen excess ratio to the desired level, even in presence of disturbance. Moreover, compared to the built-in controller, the proposed controller exhibits significantly better disturbances rejection capability. Indeed, for five steps in current (t = 48 s, t = 138 s, t = 351 s, t = 438 s and t = 528 s) an important overshot is observed with the built-in controller, whereas the proposed controller efficiently rejects the disturbances.

This control scenario is designed to assess the controller performance when the oxygen excess ratio setpoint is modified on-line according to the load current. This control scenario, which is a common working scenario in automotive applications, is significantly important. Indeed, several works demonstrated that updating the oxygen excess ratio reference according to the load current could significantly increase the PEMFC efficiency while avoiding oxygen starvation [

The proposed controller demonstrates excellent tracking capability even in presence of disturbances (

In order to highlight the benefits of the proposed parameterized NMPC scheme compared to a classical NMPC strategy, both strategies are implemented in Matlab environment and their performance are evaluated in simulation. In this aim, a mechanistic model, which has been presented and validated in a previous work [

(step changes on load current—variable oxygen excess ratio setpoint) for different values of the control horizon N u and the predictive horizon N y : test 1: N y = N u = 10 , test 2: N y = N u = 20 and test 3: N y = N u = 30 .

The computational effort is evaluated using the mean computational time (MCT) and the maximal computational time ( C T ∞ ). This last criterion represents the worst case in terms of computational time.

M C T = ∑ i = 1 n C T ( i ) n (14)

C T ∞ = max i ( C T ( i ) ) (15)

where n is the number of data and C T ( i ) the time in second required to solve the control problem at time t = i . Note that all simulations are run in Matlab environment on a computer (2.7GHz, RAM@32Go).

Note that due to its too high computation cost, the classical NMPC controller has not been experimentally tested on the real fuel cell.

Regarding the sampling period of 3 seconds, the computational cost of the parameterized NMPC controller appears to be totally acceptable.

It is important to emphasize that in addition to the presented examples, numerous experiments have been carried out in the whole range of operating conditions. In every case, the proposed NMPC controller demonstrated highly satisfactory results in terms of tracking capability and disturbances rejection. Moreover, the computational time of the proposed parameterized NMPC controller is significantly lower than the one of a classical NMCP controller.

Classical NMPC | Parameterized NMPC | |
---|---|---|

test 1: N y = N u = 10 | 0.804 | 0.102 |

MCT (s) | 1.394 | 0.122 |

CT_{∞} (s) | ||

test 2: N y = N u = 20 | ||

MCT (s) | 3.723 | 0.201 |

CT_{∞} (s) | 4.918 | 0.215 |

test 3: N y = N u = 30 | ||

MCT (s) | 9.158 | 0.301 |

CT_{∞} (s) | 11.446 | 0.322 |

Classical NMPC | Parameterized NMPC | |
---|---|---|

test 1: N y = N u = 10 | 0.804 | 0.102 |

MCT (s) | 1.394 | 0.122 |

CT_{∞} (s) | ||

test 2: N y = N u = 20 | ||

MCT (s) | 3.723 | 0.201 |

CT_{∞} (s) | 4.918 | 0.215 |

test 3: N y = N u = 30 | ||

MCT (s) | 9.158 | 0.301 |

CT_{∞} (s) | 11.446 | 0.322 |

The oxygen excess ratio is considered as a performance variable of the system and its regulation is an important issue since this parameter determines the safety of the fuel cell. In this paper, to address oxygen excess ratio control challenge, a parameterized NMPC strategy has been developed. In a first stage, due to its short computational time and its low sensitivity to noise an artificial neural network (ANN) model has been designed. The oxygen excess ratio is a function of inlet air flow rate, load current, relative humidity of air at the cathode inlet, stack temperature and inlet pressure at the cathode. However, regarding real-time control goal a simplified model has been proposed. This ANN model, used as predictor in the control strategy, is expected to predict the oxygen excess ratio of the PEMFC several steps ahead, once the load current and the inlet air flow rate are available. The validation procedure has been performed on experimental data and the model has shown good performance in terms of prediction accuracy. In a second stage, a parameterized NMPC approach has been designed. This approach, based on a particular parameterization of the control sequence, has led to a low-dimensional optimization problem. Indeed, the optimization problem that leads to the sequence of future control actions relies entirely on a single scalar parameter, which dramatically reduces the computational effort compared to a classical NMPC strategy.

Eventually, the controller has been implemented on-line and experimentally validated on a real fuel cell. Numerous control scenarios have been experimentally conducted to evaluate the controller performance, especially in terms of setpoint tracking accuracy, disturbances rejection and computational cost. These control scenarios gather all the possible scenarios in which the system would have to work. In each and every case, the controller demonstrated highly satisfactory results since it tracked efficiently the desired oxygen excess ratio value while compensating disturbances, regardless of the operating conditions. Besides, experimental comparison demonstrated that the proposed controller had much better disturbances rejection capability than the built-in controller. Eventually, to emphasize the benefit of the proposed controller in terms of computational time efficiency, a comparison has been performed between the proposed NMPC controller and a classical NMPC controller. In all case scenarios, the computational cost of the proposed parametrized NMPC controller was significantly lower than the one of the classical NMPC controller.

Overall, the proposed parameterized NMPC controller appears as an excellent candidate to address the oxygen excess ratio regulation issue.

The authors declare no conflicts of interest regarding the publication of this paper.

Damour, C., Grondin, D., Benne, M., Grondin-Perez, B. and Chabriat, J.-P. (2020) Experimental Evaluation of Parameterized Nonlinear MPC Applied to PEM Fuel Cell. Engineering, 12, 99-116. https://doi.org/10.4236/eng.2020.122010