Control of the Air Supply Subsystem in a PEMFC with Balance of Plant Simulation

This paper deals with the design of a control scheme for improving the air supply subsystem of a Proton Exchange Membrane Fuel Cell (PEMFC) with maximum power of 65 kW. The control scheme is evaluated in a plant simulator which incorporates the balance of plant (BOP) components and is built in the aspenONE R © platform. The aspenONE R © libraries and tools allows introducing the compressor map and sizing the heat exchangers used to conduct the reactants temperature to the operating value. The PEMFC model and an adaptive controller were programmed to create customized libraries used in the simulator. The structure of the plant control is as follows: the stoichiometric oxygen excess ratio is regulated by manipulating the compressor power, the equilibrium of the anode-cathode pressures is achieved by tracking the anode pressure with hydrogen flow manipulation; the oxygen and hydrogen temperatures are regulated in the heat exchangers, and the gas humidity control is obtained with a simplified model of the humidifier. The control scheme performance is evaluated for load changes, perturbations and parametric variations, introducing a growing current profile covering a large span of power, and a current profile derived from a standard driving speed cycle. The impact of the control scheme is advantageous, since the control objectives are accomplished and the PEMFC tolerates reasonably membrane damage that can produce active surface reduction. The simulation analysis aids to identify the safe Voltage-Current region, where the compressor works with mechanical stability.


Introduction
The automotive industry has strongly promoted the development of electric vehicles.Recent trends in this industry allow the commerce of vehicles powered by Proton Exchange Membrane Fuel Cells (PEMFC).The first commercial models, Hyundai ix35 FCEV or Tucson Fuel Cell (100 kW) started production at Ulsan Korea in 2013 and then distributed in Europe and other countries in 2015 [1].The Toyota Mirai model (114 kW) was launched in Japan in December 2014 [2].Honda has been developing fuel cell cars for decades, the Honda FCX Clarity (100 kW) model was produced between 2008 and 2014 but a new series is offered for lease from spring 2016 [3].Mercedes-Benz B-Class F-Cell (100 kW) is on the testing process, thus few units have been sold to reduced groups since 2010 [4].In the United States, the sales of vehicles powered by fuel cells are restricted to areas with sufficient hydrogen fueling stations [5].A public-private partnership in California is making one of the first efforts to promote hydrogen powered vehicles [6].In this state, the established hydrogen stations guarantee the autonomy of the fuel cell vehicles market.
Even when the PEMFC technology is highly attractive for portable and vehicular markets, it is not entirely mature in terms of cost, size and marketing level.Nowadays, PEMFC powered vehicle prototypes are mostly in evaluation and testing stages; at this point, a demanding task is to improve the design and the performance of the fuel cell-based power systems.The various challenges of PEMFC researches usually engage thorough analysis coming from comprehensive PEMFC system models and simulators developments.From the viewpoint of dealing with a power system, a PEMFC must operate under appropriate conditions of temperature, pressure, flow and humidity.The coupling of auxiliary equipment forms a fuel cell system and aids to attain the design operating conditions.The coordinated operation of all the components of the fuel cell system is known as the balance of plant (BOP).The dynamics of the auxiliary components have a decisive effect on the fuel cell performance, especially in systems delivering power higher than 10 kW.In particular, electric vehicles need power from 50 to 100 kW, depending on their size and use.
Summarizing, the auxiliary equipment deals with reactants supply and with temperature and humidity management in fuel cell systems.One of the former works on modeling, dynamic analysis, control and optimization of PEMFC systems was reported by Pukrushpan, 2004 [7], it was considered the basis for developing the simulator of the present study.The control of the reactants supply is of relevant importance for conceiving PEMFC power systems.Fang, 2015 [8], for example, modeled a hydrogen injection subsystem, consisting of gaseous fuel injectors that offer greater flexibility for controlling the anode pressure, and consequently for compensating abrupt pressure reduction when purge occurs.But a major concern in the system management is the air supply subsystem; herein the problems of oxygen excess ratio regulation, equalization of anode and cathode pressures and water management control are essentially addressed in PEMFC control studies using a variety of techniques, e.g., nonlinear model-based predictive control of gas pressures was applied to deal with different events of errors in PEMFC systems, along with several types of dynamic load changes [9,10].It has been acknowledged that oxygen excess ratio control is one of the most important actions to enhance protection and performance of fuel cells; thus, well-suited control strategies with this target have been proved, three related works are the following: a sliding mode control achieved by adjusting the compressor supply voltage was assisted with a nonlinear observer predicting the oxygen excess ratio and improved with a novel tuning procedure [11].In a second work, a feedforward fuzzy-PID control was proposed; the model used for designing the oxygen excess ratio controller included cathode and anode mass flow transients, membrane hydration dynamics, as well as the fuel cell BOP simulation [12], the controller was developed to adapt the PID parameters to achieve the regulation of the air flow rate using on-line fuzzy logic optimization loop.A third approach used adaptive control under exigent scenarios [13,14].Similar control objectives are addressed by designing classical PI and PIDs controllers [15] or enhanced PI tuning based on genetic algorithms [16].More recently, fault-tolerant control strategies [17] are applied.Typically, the control problem must deal with non-linear behaviors of the plant; however, control oriented linear models such as state-space models [18] and identification models [19] are also developed.
The BOP control structure also includes the regulation or monitoring of the power delivered by the fuel cell; the management of temperature to treat the problems of feeding reactants temperature regulation, as well as the fuel cell stack cooling to mitigate the heat produced by the exothermic reactions using a heat removing circuit [20][21][22][23][24].A current topic is the development of hybrid systems [25][26][27][28][29][30][31][32]; therefore, the addition of peripheral elements, such as auxiliary power units or backup sources intensifies the challenge of modeling highly coupled multi-physical systems.There are also complete studies concerning the design, improvement or modeling of converters for fuel cell applications (e.g., [33][34][35]).Kim, 2015 [36] proposed a process for implementing new renewable energy systems, and offered an analysis of the key factors affecting system performances of fuel cell systems.Some simulation platforms provide tools and libraries that facilitate the integration of subsystems in specialized simulation environments.Studies about PEMFC systems have presented advanced process simulators, for example, Ziogou, 2011 [37] reported the use of the gPROMS platform to simulate a PEMFC BOP.A well-known software for process engineering applications is the Advanced System for Process Engineering (aspenONE R v8, Aspen Technology, Inc., Burlington, MA, USA).This simulation platform has been widely applied for process design and development on various renewable energy sources areas, some examples are the works on biomass or solid waste utilization for gasification [38][39][40][41].More related works are those of Rabbani, 2013 [42], who took advantage of the Aspen Plus R simulator to model and control a Ballard power module of 21.2 kW, and presented a dynamic analysis for vehicular applications.Also within the Aspen Plus R platform, two simulators of fuel cell systems were built in which a steam reforming subsystem was considered for hydrogen production: Salemme, 2010 and 2013 [43,44] calculated the operative conditions that maximize the energy efficiency of a PEMFC system; and Suther, 2010 [45] developed a simulator with a macro-level model for a Solid Oxide Fuel Cell (SOFC) stack, this model yields the output composition of the exhaust and the produced work, and calculates the heat available for the fuel reformer and the electrochemical properties of the SOFC.
In this paper, a control scheme for a PEMFC with BOP components simulation is presented.The dynamic simulation of the PEMFC BOP is executed in the aspenONE R platform for maximal power solicitations of 65 kW.The built simulator includes a compressor for oxygen supply, two heat exchangers to regulate the temperature of the reactants at the electrodes inputs, and a humidifier to saturate with water the air entering into the fuel cell.The overall system diagram is depicted in Figure 1 and the control structure is entirely described in Section 3. The emphasis is on introducing a practical framework for assessment of the control performance in the reactants supply subsystem, given that: (1) The controllers performance evaluation is made in terms of the interactions and constrains imposed by the auxiliary equipment; thereby preventing unsafe or degraded operating conditions.(2) The operating condition definition and sizing of the auxiliary equipment is made in accordance with the PEMFC capabilities.(3) The static and dynamic characteristic responses of the auxiliary components are captured in the simulator.(4) The simulator facilitates control tests under changes in the system inputs, disturbances and parametric variations.
The aspenONE R was selected as the process simulation platform because this software includes unit operations modeling blocks, dynamic simulation tools and own libraries programming facilities, it is used for academic purposes and in industry projects, counting a number of plant operating challenges.Moreover, it can be interfaced to the Matlab/Simulink platform, enabling the use of control tools and toolboxes.

Operation and Simulation of the PEMFC
The operating principle of a mono-cell is shown in Figure 2. The process begins with the hydrogen and oxygen supply at fixed temperature and pressure.Hydrogen is uniformly distributed through the flow channel for reaching the fuel cell anode, and the air supplied by a compressor enters through the distribution channels of the cathode.At the anode, the molecules of H 2 are diffused through the pores of the electrodes until reaching the catalytic layer where they are linked to the catalyst, then the molecules are dissociated into atoms.That is, molecular hydrogen (H 2 ) is converted into protons (H + ) and electrons (e − ) by the oxidation reaction occurring at the anode.H + ions migrate through the electrolyte by jumps.These ions are temporarily attached to negatively charged radicals, contained in the electrolyte membrane; then, these ions are detached and dragged to the cathode by the water contained in the mono-cell.This mechanism gives the name to this type of fuel cell, since proton H + exchange occurs between the electrodes throughout the membrane.On the other hand, electrons are conducted through an external electrical circuit for powering a load.The electron current is the main product of the fuel cell.In parallel, mechanisms that occur at the cathode are: H + ions migrate to the cathode through the membrane, the electrons arriving there by the external circuit, and oxygen molecules O 2 transported to the surface are linked to the catalyst to form water by a reduction reaction.Water is a by-product, together with the heat caused by the exothermic reactions.A mono-cell generates a voltage not greater than 1.23 V, and it typically is lower than 1.0 V, whereas a set of mono-cells electrically connected in series provides a greater voltage and the power required by the load, this arrangement forms the fuel cell stack.
The performance of a fuel cell stack is determined by the current (I)-voltage (V) relationship which is measured at the terminals of the fuel cell at constant temperature and pressure.This I-V relationship characterizes the steady state performance of the PEMFC and is known as the polarization curve.The dynamic response is characterized by the double layer capacitance effect (electrical dynamics), transients due to continuous thermal, flow, pressure or humidity changes, among others mechanisms occurring on the fuel cell stack and the auxiliary equipment, which conduct to different I-V paths on the polarization curves mapping (at different pressures).In this study, the PEMFC polarization curves of the model proposed by Pukrushpan [7] was used to build the PEMFC library.The mathematical model was programmed in Aspen Custom Modeler R and thereafter it was incorporated as a library into the process flowsheet simulation in Aspen Dynamics R .This model considers the performance of a mono-cell; but the stack consists of n mono-cells connected in series, so it is assumed that the voltage losses due to the stack assemblage are negligible and the fuel cell output voltage v st is calculated using (1): The voltage of a mono-cell v f c depends on the thermodynamic voltage or Nerst voltage E at open circuit, as well as the voltage drops or voltage overpotentials due to activation loss v act , ohmic loss v ohm , and concentration loss v conc , which make the process irreversible.The voltage of each single cell v f c is calculated using (2): ( The PEMFC model is summarized in (1)- (15).With respect to the operating variables, an ideal stack temperature T st is assumed to be regulated at 80 • C, i is the current density, p H 2 an and p O 2 ca are the partial pressures of hydrogen and oxygen, p ca is the pressure at the cathode, p an is the pressure at the anode.To estimate the ohmic overpotential, R ohm is the internal electrical resistance of a mono-cell, σ m is the membrane conductivity, t m is the membrane thickness, λ m is the membrane water content.
The values v 0 and v a depend on the partial pressures of the reactants and the temperature of the stack, p sat stands for the saturation pressure of water vapor and, together with v 0 and v a are calculated from nonlinear regression.Particularly, the membrane conductivity varies with the membrane water content and fuel cell temperature, this dependence was determined empirically for Nafion 117 membrane in terms of the model constants, reported also in [7]: c 1 = 10, c 3 = 2, b 2 = 350, b 11 = 0.005139, b 12 = 0.00326.Conversely, the concentration looses depend on the temperature and the reactant partial pressure, hence, the concentration overpotential was also determined empirically, in this case, in terms of c 2 and c 3 , with c 2 subject to the conditional statements of (12).Here, v con is determined also by i max = 2.2, which is the current density than causes abrupt voltage drop in the concentration region.Finally, the number of mono-cells considered is n = 381, P st is the total power delivered by the cell, P cp is the power consumed by the compressor and P net is the net system power delivery.The general stack voltage model takes into account the effect of the operating temperature and pressure.Figure 3 shows the polarization curves and power curves that characterize the steady performance of the PEMFC.
The PEMFC model is formulated based on the material and energy balances in the supply and discharge lines at the anode and the cathode [7].It is assumed that the fuel cell stack operates at constant temperature and a cooling circuit removes heat from the exothermic reactions.Under load variations, the fuel cell operates in a variable pressure mode between 139.9 and 294.7 kPa; the polarization curves in this pressure range are observed in Figure 3. Variable-pressure operation is typically carried-out with a blower (lower power systems) or a compressor (higher power systems), but the anode and cathode pressures should be in equilibrium.

Design and Simulation of the PEMFC BOP Components
The performance of a fuel cell is determined by the polarization curves, the electrical efficiency and the power curves.The equipment included in the simulator was sized and specified with a basis on these performance features.A special emphasis is on the air supply management, since the fuel cell aging is directly related to the air flow rate condition.In effect, oxygen starvation is attributed among others, to factors such as insufficiency in oxygen supply due to sudden increases in power requirements, non-uniform oxygen distribution or poor control of anode and cathode pressure equilibrium [46].A further difficulty taking place at the cathode if oxygen is depleted is that proton (H + ) reduction occurs instead of oxygen reduction [47].Successive oxygen starvation or depletion processes may produce permanent damage or degradation of the fuel cell, which contributes to the reduction of its lifetime.These problems are closely handled with the air supply management and control tasks and some solutions are treated in the control scheme evaluation section.
The oxygen required by the fuel cell is provided by the air supply of an Allied Signal centrifugal compressor, whose operation is represented by performance curves or characteristic curves gathered in a compressor map [48].The characteristic curves reflect the correlation between the compression ratio (outlet pressure/inlet pressure), and the air mass flow that the compressor delivers (Figure 4).The left boundary is the surge line defining the minimum flow limit, at which, for a given speed, the compressor becomes mechanically unstable; whereas the right boundary is the choke or stonewall line, defining the maximum flow that the compressor can deliver.The areas encountered out of these boundaries marked with blue lines in Figure 4 represent the conditions at which problems of pumping and detachment arise.When the air flow is below the surge limit, flow reversal occurs, then noisy and vibrating operation leads to mechanical damage.If the flow is supersonic (the gas velocity is higher than the sound velocity), then, the air demand becomes greater than the compressor capacity.Therefore, the described boundaries outline the compressor capacity and the operating area, free of mechanical instabilities.
The compressor efficiency (Figure 5) is also determined by the pressure ratio and the mass flow, its maximum value is 80%.The region marked in yellow is the most efficient zone and defines the conditions of mechanical stability.
The simulation of the compressor in Aspen Plus R was done using the COMPR modeling block that includes an isentropic compressor model.The parameter that determines the feasibility of approximating a real process with an idealized process is the isentropic or adiabatic efficiency.The COMPR module calculates the discharge conditions from the performance and efficiency curves; therefore, these must be introduced in the module COMPR.The compressor operates with a mechanical efficiency of 80%.The maximum current demand that the fuel cell can tolerate corresponds to an air flow demand of 0.0815 kg/s; under these conditions, the compressor operates at 100,000 RPM and the air comes out at a temperature of 173 • C.
On the other hand, at small air flows, the compressor operates close to the surge limit, loosens the capacity of increasing the outlet pressure, becomes unstable (Figure 4), and the efficiency is considerably reduced (as it can be seen in Figure 5).Taking into account these drawbacks, the PEMFC system is not able to operate at small flow rate conditions, consequently, the simulation of the PEMFC BOP cannot be conducted at startup conditions or under small flow conditions.Although the PEMFC electric performance could operate a wider range of action, the compressor capacities limit the power delivery range.In addition, extremely low reactant flow rate into the fuel cell leads to water evaporation and flooding shifts to the cathode input and likely to air starvation, which in turn, yields undesired output voltage transients and fuel cell performance degradation [46].These difficulties confirm the necessity of operating in an intermediate range of air flow rate, or equivalently, an intermediate current range in the fuel cell.Heat exchangers modify the temperature of the air delivered by the compressor and the temperature of the hydrogen feeding the fuel cell from a pressurized tank; so that the gases can reach the nominal operating temperature of the fuel cell stack.Two tube and shell heat exchangers were sized with the aid of the module Aspen Exchanger Design and Rating R .
The gas flow to treat, the operating pressure and the temperature gradients in heat exchangers are determined at every time by the voltage-current variation in the PEMFC terminals and by the compressor dynamic response.Accordingly, the specification of the input streams is necessary to correctly size heat exchange equipment and then to model it in the process flowsheet.The heat exchanger at the cathode side was specified for a maximum input air flow of 0.0815 kg/s at 173 • C and 3.5 bar, while the heat exchanger at the anode side was sized for changing the temperature of a maximum hydrogen flow of 0.0484 kg/s, entering at 104 • C and 3.5 bar.In both cases the cooling fluid is water, circulating at 25 • C and 3.5 bar.
The main issues of the geometric design from Aspen Exchanger Design and Rating R are provided in Table 1.The heat transfer areas are 0.2419 m 2 and 0.2742 m 2 for the air and hydrogen heaters, respectively.The simulation of the heat exchangers was done in Aspen Plus R using the HEATX module that offers rigorous modeling for shell and tube heat exchangers, with accurate thermal and hydraulic prediction, and establishes the thermodynamic conditions of phase changes if necessary.The design parameters of Table 1 were introduced in the modeling library to simulate the steady state performance; afterwards, the dynamic simulation was carried out in the Aspen Dynamics R platform, including into the flowsheet the rest of the PEMFC BOP components.
The humidifier simulation was simplified using the MIXER module, the modeling block simulates an adiabatic mixer that combines multiple streams into one.It is a rough representation of the humidifier.To saturate the air with water, an air stream is combined with a water vapor stream, so that the humidity controller furnishes instantaneously the necessary vapor content to the oxygen entering into the fuel cell stack, but does not consider the thermal or diffusive transport dynamics; however, the applied humidity controller action prevents the dehydration problem.On the other hand, the fuel cell flooding affects mass transport and produces greater concentration losses, predominantly at the cathode.Excess of water obstructs the pores of the gas diffusion layer and significantly reduces the transport rate of the reactants to the catalysts active sites.As a final consequence, flooding leads to starvation and greater potential drops.The considered voltage stack model includes to some extent the flooding due to the stack pressure increase, since the ohmic and the concentration overpotentials depend on the water content (14) and on the stack pressure ( 11)- (15), respectively.Nevertheless, the simulation results are not analyzed with respect to voltage variation due to the fuel cell flooding; since a water balance derived from basic principles of gas-phase transport is needed to enable the analysis of the effect of the water flooding in the cathode.Valves were included in the simulator, thus the actuators modeling is considered for different control loops.The VALVE module was used to simulate pressure drops in adiabatic valves.This module also calculates the temperature and phase of the output current.The PEMFC BOP simulation was done in the Aspen Dynamics R environment.Figure 6 shows the flowsheet of the open loop system.

Control Scheme of the PEMFC System
The proposed control scheme involves the following control loops (see Figures 1 and 7): The oxygen excess ratio λ O 2 is regulated with the manipulation of the power supplied to the compressor P cp ; the pressure tracking at the anode p an is achieved by manipulating the hydrogen flow rate W an,in , considering the cathode pressure p ca signal as the reference value; the temperature of the reactants entering at the fuel cell stack T ca,in and T an,in are regulated in the heat exchangers with the manipulation of the cooling water flow rate; the relative humidity φ ca is regulated by varying the flow rate of the water in the humidifier.
Summarizing, the models of the compressor, the heat exchangers, the humidifier as well as the control valves and the PID controllers used in different control loops are available in the simulator libraries, whereas the PEMFC model and the adaptive controller designed to regulate the temperature of the air in the cathode side were programmed in Aspen Custom Modeler R .The resulting libraries are exported to Aspen Dynamics R in order to be integrated in the plant flowsheet.

Pressure Controller Design for the Anode Side
Small changes of hydrogen flow rate W an,in affect immediately the anode pressure p an .Fluctuating loads are associated to important changes on the hydrogen requirement, specially when the fuel cell operates in variable pressure mode; this is the case of the PEMFC system under study.Herein, pressure differences between the electrodes arise and may produce polymeric membrane breaking.To prevent this problem, it is necessary to minimize the pressure difference between the cathode and the anode areas.The controller objective is the anode pressure p an (t) tracking to follow the pressure changes at the cathode p ca (t).
A control loop between p an (t) and p ca (t) is proposed, see Figure 8.The cathode pressure is the reference, the controller C(s) is a PI and G(s) is the transfer function between the anode pressure p an and the flow of hydrogen W an,in .System identification is used to determine G(s).As the input, a pseudo random binary signal (PRBS), representing W an,in variations is applied.Then, the output p an is recorded.With the input-output data a ARX model structure is estimated by the least square error method, in this case, a first-order transfer function ( 16) is enough to obtain good fit.
To analyze the dynamics of the anode pressure in function of the hydrogen flow variation, a simulation was done considering variable load (variable I st ), regulated oxygen excess stoichiometric ratio (regulated λ O 2 ) and non-equilibrated pressure condition between the electrodes areas (p ca = p an ) as the initial state.C(s) is a PI controller, C(s) = K p + K i /s, then the transfer function in closed loop p an /p ca is a second order: Because pressure in the anode must follow the pressure in the cathode, then the closed loop steady state gain must be equal to 1.We choose the pole placement technique to calculate the controller in Figure 8.Because the transfer function in a closed loop is modeled as being a second order, we need two poles p 1 and p 2 .Given the unity gain feedback system, it becomes true that: From ( 17) and ( 18) we solve for the gains k p and k i as follows: k i = p 1 p 2 /3.364 × 10 6 and k p = −(p 1 p 2 + 8.032)/3.364× 10 6 .One can remark that both gains are divided by 3.364 × 10 6 , thus it is expected to obtain gains with small values.On the other hand, small gains compensate the high steady state value derived from the relationship between anode pressure and hydrogen flow rate ( 16), (4.188 × 10 5 ).
In order to select the poles, it should be considered again that the PEMFC works with variable pressure.Therefore, it is advisable to be prudent, and not significantly modify the natural dynamics produced by the pole at 8.032, so p 1 = 10 is selected, just a little bit faster.The second pole p 2 should not intervene significantly, for this reason it is placed as far as possible; on the other hand its distance is proportional to the gains k p and k i .With this in mind, p 2 is chosen far enough to have small values of k p and k i , but manageable yet, p 2 = 340.Thus k p = 1.0166 × 10 −4 and k i = 0.001.

Air Temperature Controller Design for the Cathode Heat Exchanger
The temperature of the air at the cathode side and the temperature of hydrogen at the anode side are regulated to the fuel cell operating temperature (80 • C) in order to prevent damage to the polymeric membrane; this is achieved by varying the cooling water flow rate in the heat exchangers.
Due to the nonlinear dynamics of heat exchangers, and because the air delivered by the compressor at the cathode side exhibits significant temperature fluctuations, it could be advantageous to use an adjustable controller, then a Model Reference Adaptive Control (MRAC) is proposed to regulate the air temperature.The idea is to have a closed loop controller with parameters that can be updated by an adjustment mechanism that compares the output of the plant with a model.
The MRAC is built as follows: The correction signal is the error between the temperature of the output stream in the air heat exchanger simulator and the output of the reference model.Since this is a regulation loop, the second signal is a constant of 80 • C and the error becomes e(t) = T ca,in − 80.The control signal is the cold water flow through the jacket of the heat exchanger W incold1 , it is composed of a feedback part, with a setting parameter θ 2 multiplied by the temperature in the cathode T ca,in , and another feedforward part with a setting parameter θ 1 multiplied by the reference, W incold1 = 80θ 1 + θ 2 T ca,in .The parameters θ are adjusted with the MIT rule, the gradients are calculated with the partial derivatives of the error with respect to the parameters θ.As mentioned before, there is not a model for θ 1 because it is a regulation loop, conversely, for θ 2 a first order model (19) with unitary state stable gain is used.
Afterwards, the cathode temperature control is set up as in Figure 9.The first order model is adjusted by translating the time constant of the response measured in the plant, into a change of T ca,in , as a result of the manipulation of the cold water flow rate W incold1 (a m = 6).In this procedure, the gain of the controller γ is chosen by considering the numerical relationship between the variables in the previous controller (16).In that case, the gain is then small, thus a gain with the same order of magnitude is chosen γ = −3 × 10 −6 .The control tool in Aspen Dynamics R handle only PID controllers.The application of different control techniques must be addressed with external toolbox or developing proprietary libraries.Hence, the adaptive controller was programmed in Aspen Custom Modeler R to create the MRAC library, afterwards, it was integrated into the simulator.

Other Regulatory Controllers
Three additional control loops are necessary, all of them in regulation mode, with no major challenge in their design.To deal with the regulatory control problem, classic PIs were tuned based on the response to a step-type excitation using the methods of Ziegler-Nichols, Cohen-Coon and Chien, Hrones and Reswick (CHR).Subsequently, to confront the various configurations, four performance indices were estimated: the Integral of Squared Error (ISE), the Integral of Absolute Error (IAE), the Integral of Time Multiply Squared Error (ITSE) and the Integral of Time Multiply Absolute Error (ITAE).The preferred configuration was the one leading to the best indices, although slight differences were encountered when comparing the controllers performances.
The first PI controller is configured to regulate the excess oxygen ratio regulation.The oxygen excess ratio λ O 2 is defined in (14), in this expression W O 2 ,in is the oxygen mass flow rate in the fuel cell cathode side and W O 2 ,react is the oxygen reacting to produce water at the anode side.According to the stoichiometric balance, λ O 2 relates the oxygen flow accessing in excess into the cathode with respect to the oxygen required to carry out the reaction.The fact of maintaining λ O 2 at a value of 2 contributes to the total consumption of hydrogen.The regulation of λ O 2 at this reference value is subjected to effectively modify the air flow rate by manipulating the compressor power.The selected PI has the gains k p = 3.1 × 10 4 and k i = 0.03, deduced by the Cohen-Coon method.
The second PI controller is configured to regulate the heat exchanger anode temperature.Unlike the cathode, the anode feed comes from a hydrogen storage tank, and the intermediate heat exchanger input is subjected to unimportant temperature variations.The chosen PI has the gains k p = 5 and k i = 2 deduced by the CHR method.
The controller to regulate the relative humidity is calculated as the relation between the partial pressure and the saturation pressure of the water steam.The signal error for this control loop was calculated by using (20).Where y H 2 O,ca,in is the water steam composition at the input of the cathode p ca,in is the water steam composition at the input of the cathode.The PI has the following gains: k p = 2 and k i = 10.

Results and Discussion
The PEMFC system under study is a Multi-Input Multi-Output (MIMO) dynamical system.The control problem is treated with decoupled controllers; nevertheless, attempts were done to overcome the non-lineal behaviors and interactions between the BOP components and between the various control loops of the system.The control scheme provides a suitable solution to the regulatory and tracking objectives, it should prevent instabilities, damage and degradation or should serve as analysis tool to determine the PEMFC system operating boundaries.Therefore, the control scheme effectiveness is interpreted in that sense from tests carried out with two distinctive load profiles: (1) The first profile is the same load current used in [7], which consists of progressive increments of current that enables the PEMFC operation over the entire possible voltage-current region.This load profile serves to evaluate the control scheme under operational changes, system disturbances and parametric variations.The proposed tests should evidence that the fuel cell can operate in the safe and efficient region, at the same time that the auxiliary subsystems do.(2) The second load is a standard driving cycle, 300 s of the motorway Artemis cycle for maximum speeds of 130 km/h (ARTMW130) were considered.The Driving Cycle Simulink Block [49] was used to produce the speed profile, representative of driving conditions.Then, the power load P Motor was estimated from typical small vehicle characteristics as stated in ( 21) [50].Finally, the current demand is deduced from this profile by considering the performance of the PEMFC.
The resultant load profile serves to evaluate the control scheme under faster dynamics and covering a different power range with respect to the first load.The proposed tests should highlight the efficacy to attain the regulatory and tracking objectives and should be useful to determine the operating boundaries of the system, based on the performances of auxiliary subsystems. where: 225 kg/m 3 ).S = Front surface (2.5 m 2 ).

Tests under Operational Changes
The first test serves to analyze the performance of the control scheme under changes in the system variables, with special attention to the regulation of the excess oxygen stoichiometric ratio.The load profile consists of progressive increments of current that conduct the PEMFC to operate over the entire possible region (with air flow demands until 0.0815 kg/s).Given a positive load increase, the power consumed by the compressor grows with the oxygen utilization (Figure 10).For a maximum power demand (65 kW), the compressor consumes about 12 kW, i.e., about 18.5% of the power generated by the PEMFC.Under these conditions, the net power (power produced by the fuel cell-compressor consumption) is reduced to about 53 kW.The regulation of the oxygen excess ratio to a value of 2 (Figure 10), successfully conducts the compressor to operate in a flow path-pressure within the zone free of mechanical instability, without exceeding the limits of detachment and pumping (Figure 11).The PEMFC system powers a variable load under isothermal conditions, thereby, the pressure must be modified to guarantee the energy supply.Figure 12 shows the pressure trajectory on the polarization curve of a mono-cell at the nominal fuel cell temperature (80 • C).
It can be seen that the voltage delivered by each mono-cell is around 0.6 V, regardless of the current and power changes.However, the current density covers a range of 0.3 to 1.1 A/cm 2 , which corresponds to a large part of the ohmic region, where the operation of the fuel cell is safe.This region has a good balance between efficiency and power, because at lower currents, the power diminish (Figure 3), whereas at higher currents, the efficiency decreases.With the oxygen consumption, the hydrogen utilization increases to follow the cathode pressure variation.The pressure tracking aids to minimize effectively the pressure difference in the electrodes area, although the cathode pressure is modified in a wide range from 140 to 300 kPa (Figure 13).

Tests under System Disturbances
The change of temperature in the air delivered by the compressor can represent important perturbations for regulating the temperature of the air entering at the fuel cell cathode (Figure 14).The compressor supplies hot air that reaches until 173 • C at the highest possible air flow.The adaptive controller regulates properly the air temperature in the cathode inlet without exceeding 0.4 • C over the reference value of 80 • C (fuel cell nominal temperature).The controller rejects the disturbances quickly and provides a unitary gain to the regulated variable.In a closed loop, the variation of 0.4 • C cannot be exceeded, since this value is obtained with the maximum disturbance on the temperature at the compressor outlet (Figure 14) that corresponds to the maximum current demanded by the load (Figure 10).However, if exceeded, there is a risk of affecting the humidification of the reactants, causing dehydration of the membrane because at high stack operating temperatures, the reaction of water formation at the cathode may be insufficient to compensate the lack of water; therefore, the membrane becomes drier and the conductivity decreases, leading to higher ionic resistance and temporary larger ohmic losses.A long-term dehydration condition causes irreversible damage to the membrane and develop cracks, along with gas crossover and hot spots (high chemically active areas due to the exothermic reactions), resulting in more gas crossover.Once this process is initiated, a destructive sequence is established, hence, dry membrane condition conducts to shorter the fuel cell life [47].The effect of operating the fuel cell with membrane damage is exemplified in Section 4.3.Moreover, the variation of the air flow rate in the compressor can disturb the relative humidity at the fuel cell cathode input.The humidity controller aids to inject enough steam for controlling the relative humidity.Because the model is quite simplified and the control is perfect, the controller performance is not presented here; instead Figure 15 shows the vapor injected to the air stream allowing perfectly humidified air fed.

Test under Parametric Variations
In order to evaluate the control scheme performance under parametric changes, it is supposed that the polymer membrane is partially damaged, thereby ensuing the system would fail.A reduction of 15% in the effective area is assumed, i.e., decreases from 0.028 to 0.024 m 2 .
The simulation shows that the PEMFC system undergoes natural degradation and delivers reduced voltage and power at its terminals (Figure 16); but, the general degraded behavior does not represent a large deviation of the fuel cell normal operation, since the power delivered by the PEMFC does not decrease more than 6% during the increasing demand of the load, while the regulation of the oxygen excess ratio, effectively allows the compressor to operate with mechanical stability and efficiency despite membrane damage.

Test under Driving Solicitations
The second load profile is representative of a standard driving cycle on the motorway, with maximum speed of 130 km/h.The mechanical parameters used in (21) to compute the power demand of the vehicle motor P Motor correspond to the characteristics of a small vehicle.Figure 17 shows the vehicle solicitations in function of the current and the electric power demanded to the PEMFC, which cover span intervals until 200 A or 40 kW, respectively.Control difficulties may come due to low air flow rates or due to fast and fluctuating input dynamics.Nevertheless, the control scheme must provide useful insight to mark boundaries of the PEMFC system applicability.
The anode pressure tracking control was made with the following results: The anode pressure fits perfectly with the cathode pressure during the whole test even though fast current transients and sudden large changes occurs continuously (Figure 18).Consequently, the risk of polymeric membrane breaking diminish considerably.The control effort is minor, since the hydrogen flow rate varies from 0 to 0.03 kg/s and still follows the fast dynamics.
The air temperature is properly regulated at 80 • C as it was the case of the first load considered, whereby the temperature control efficacy was demonstrated (Figure 14).The hydrogen temperature is also regulated around 80 • C; but the transients are significantly smoother in comparison with the air temperature variation because the hydrogen is delivered by a tank, while the air is delivered by a compressor operating with important transients.For this reason, only the air temperature control was presented.Moreover, the air is without difficulty saturated with water, Figure 19 presents the water vapor injected to the air stream in order to maintain saturation conditions, the required water content varies with the pressure.To conclude the control strategy test under driving standard conditions, the oxygen excess ratio λ O 2 control performance is also presented.The regulatory control of the stoichiometric oxygen excess ratio is successfully achieved, λ O 2 is conserved between 1.95 and 2.08, while the control effort defined by the power consumption of the compressor remains largely, close to 2 kW, with variations achieving until 45%-48% of the maximum compressor power capacity (Figure 20).Summarizing, the control scheme contributes to attain the reactants flow rate, temperature, pressure and humidity, needed by the PEMFC stack to supply the load demand.All the control objectives are attained without overshooting or delays; however, further analysis should consider the auxiliary equipment constrains.Figure 21 shows the fuel cell path relating the air flow rate and the compression ratio over the compressor map.During the driving cycle, the flow conditions exceeding the surge limit are low power demands that conduct to the lowest compressor efficiency region.These difficulties arise due to startups and fast changes in the load.Clearly, a backup source is mandatory to overcome these difficulties, but an energy management strategy is necessary to determine properly the power split between the sources.The PEMFC operation should be restricted by defining a low power limit to overcome low efficient operation of the compressor, in the vicinity of the surge limit.Figure 22 exhibits the PEMFC operation for solicitations over 100 A, for this current demands, the compressor operates 92% of the time between the safety limits.The signals that continue overtaking the surge limits arise due to the fast load transients.Typical energy management strategies consider frequency analysis to determine the dynamic capacity of the fuel cell in order to limit conveniently the the fuel cell intervention.
Considering the tests realized with both current profiles, some concluding issues are that the PEMFC system under study should be limited between 23 to 60 kW, or equivalently from 100 to 250 A. On the other hand, supplementary frequency analysis is mandatory to determine the dynamic capabilities of the fuel cell.The compressor capacity analysis realized herein is often not considered to determine the power split between the sources in hybrid generators.This analysis should improve the power management design, it becomes possible because auxiliary equipment is integrated in the fuel cell system simulator.The contribution may come as it should be expected that dehydration and flooding can be avoided at the same time that air starvation, undesired output voltage transients and fuel cell performance degradation.As it was explained before, special emphasis on the air supply management is important since the fuel cell aging is directly related to the air flow rate conditions, as well as the mechanical instability and damage problems in the compressor.

Conclusions
The developed simulator of a power generation system combines the models of a PEMFC stack, a compressor feeding the fuel cell stack, two heat exchangers modifying the reactants temperature to attain the anode and cathode sides nominal temperature, and a humidifier saturating with water the air entering into the fuel cell stack.The auxiliary equipment was sized based on the PEMFC performance, using design tools from the simulation environment.In addition, two built libraries served for modeling the PEMFC stack behavior and applying the MRAC control into the air temperature control loop.The air supply control scheme showed effective performance in equilibrating the cathode and anode pressure, and for regulating the temperature and humidity of the air and hydrogen input streams.The tests showed that the control scheme, designed for a multi-variable, non-linear and highly coupled system, tolerates well the operating variables changes, disturbances and parametric variation.The consideration of the compressor map enables to restrain the capabilities of the PEMFC system to prevent mechanical instability and damage, and fuel cell aging; although an additional frequency analysis is necessary to determine the power density capacity of the fuel cell.
The control scheme proposed can facilitate the control implementation on real systems, since the controllers can be commercially available.However, it was necessary to incorporate an adaptive controller in the air temperature loop due to the nonlinear dynamics of the heat exchanger.
The set Aspen Plus R -Aspen Exchanger Design and Rating R -Aspen Dynamics R -Aspen Custom Modeler R is advantageous to simulate the BOP, because it makes possible to obtain more realistic behavior of the auxiliary equipment, and to consider thermal, thermodynamic and fluidic phenomena occurring in the system.The built libraries allow to consider the electrochemical phenomena that characterize the fuel cell operation that are not available in the process simulator.

Figure 1 .
Figure 1.Process diagram flowsheet of the PEMFC system with their respective control loops.

Figure 3 .
Figure 3. Fuel cell polarization and power curves.

Figure 4 .
Figure 4. Performance curves of the centrifugal compressor.

Figure 5 .
Figure 5. Compressor efficiency as a function of pressure ratio and mass flow.

Figure 6 .
Figure 6.Flow diagram of PEMFC system in Aspen Dynamics R .

Figure 7 .
Figure 7. PEMFC control scheme, tested in the simulation environment Aspen Dynamics R .

Figure 8 .
Figure 8. Block diagram of the anode pressure controller.
γe a m s+a m T ca,in .

Figure 9 .
Figure 9. Block diagram of the adaptive controller to regulate the air temperature.

Figure 12 .
Figure 12.Current-Voltage path on a mono-cell operating at 80 • C and 139.9 to 294.7 kPa.

Figure 13 .
Figure 13.Minimization of the pressure difference between the cathode and anode, manipulating the hydrogen mass flow.

Figure 14 .
Figure 14.Temperature control air to the cathode inlet, manipulating the flow of cooling water.

Figure 15 .
Figure 15.Water vapor injected to the air stream for maintaining humidity saturation condition.Current increasing steps load case.

Figure 16 .
Figure 16.System response to a parametric variation of 15 % in the effective area of the polymer membrane.

Figure 19 .
Figure 19.Water vapor injected to the air stream for maintaining humidity saturation condition.Standard driving cycle load case.

Figure 20 .
Figure 20.The excess oxygen ratio regulation.

Figure 21 .
Figure 21.Air compression path conforming to the PEMFC solicitations.

Figure 22 .
Figure 22.Air compression path for PEMFC solicitations greater than 100 A.