1. 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
® 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® 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® 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
® 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® 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.
2. Dynamic Simulation of the PEMFC System
2.1. 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 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 () is converted into protons () and electrons () by the oxidation reaction occurring at the anode. 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 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: ions migrate to the cathode through the membrane, the electrons arriving there by the external circuit, and oxygen molecules 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® and thereafter it was incorporated as a library into the process flowsheet simulation in
Aspen Dynamics®. 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
is calculated using (
1):
The voltage of a mono-cell
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
, ohmic loss
, and concentration loss
, which make the process irreversible. The voltage of each single cell
is calculated using (
2):
The PEMFC model is summarized in (1)–(15). With respect to the operating variables, an ideal stack temperature
is assumed to be regulated at 80
,
i is the current density,
and
are the partial pressures of hydrogen and oxygen,
is the pressure at the cathode,
is the pressure at the anode. To estimate the ohmic overpotential,
is the internal electrical resistance of a mono-cell,
is the membrane conductivity,
is the membrane thickness,
is the membrane water content.
The values
and
depend on the partial pressures of the reactants and the temperature of the stack,
stands for the saturation pressure of water vapor and, together with
and
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]:
= 10,
= 2,
= 350,
= 0.005139,
= 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
and
, with
subject to the conditional statements of (12). Here,
is determined also by
= 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,
is the total power delivered by the cell,
is the power consumed by the compressor and
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.
2.2. 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 (
) 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® 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 .
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®.
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 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 and 3.5 bar. In both cases the cooling fluid is water, circulating at 25 and 3.5 bar.
The main issues of the geometric design from
Aspen Exchanger Design and Rating® are provided in
Table 1. The heat transfer areas are 0.2419
and 0.2742
for the air and hydrogen heaters, respectively. The simulation of the heat exchangers was done in
Aspen Plus® 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® 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® environment.
Figure 6 shows the flowsheet of the open loop system.
3. Control Scheme of the PEMFC System
The proposed control scheme involves the following control loops (see
Figure 1 and
Figure 7): The oxygen excess ratio
is regulated with the manipulation of the power supplied to the compressor
; the pressure tracking at the anode
is achieved by manipulating the hydrogen flow rate
, considering the cathode pressure
signal as the reference value; the temperature of the reactants entering at the fuel cell stack
and
are regulated in the heat exchangers with the manipulation of the cooling water flow rate; the relative humidity
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®. The resulting libraries are exported to Aspen Dynamics® in order to be integrated in the plant flowsheet.
3.1. Pressure Controller Design for the Anode Side
Small changes of hydrogen flow rate affect immediately the anode pressure . 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 tracking to follow the pressure changes at the cathode .
A control loop between
and
is proposed, see
Figure 8. The cathode pressure is the reference, the controller
is a PI and
is the transfer function between the anode pressure
and the flow of hydrogen
.
System identification is used to determine
. As the input, a pseudo random binary signal (PRBS), representing
variations is applied. Then, the output
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
), regulated oxygen excess stoichiometric ratio (regulated
) and non-equilibrated pressure condition between the electrodes areas (
) as the initial state.
is a PI controller,
, then the transfer function in closed loop
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
and
. Given the unity gain feedback system, it becomes true that:
From (17) and (18) we solve for the gains and as follows: and . One can remark that both gains are divided by , 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), ().
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 , so is selected, just a little bit faster. The second pole 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 and . With this in mind, is chosen far enough to have small values of and , but manageable yet, . Thus and .
3.2. 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 ) 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
and the error becomes
. The control signal is the cold water flow through the jacket of the heat exchanger
, it is composed of a feedback part, with a setting parameter
multiplied by the temperature in the cathode
, and another feedforward part with a setting parameter
multiplied by the reference,
. 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
because it is a regulation loop, conversely, for
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
, as a result of the manipulation of the cold water flow rate
(
= 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
.
The control tool in Aspen Dynamics® 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® to create the MRAC library, afterwards, it was integrated into the simulator.
3.3. 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 is defined in (14), in this expression is the oxygen mass flow rate in the fuel cell cathode side and is the oxygen reacting to produce water at the anode side. According to the stoichiometric balance, 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 at a value of 2 contributes to the total consumption of hydrogen. The regulation of at this reference value is subjected to effectively modify the air flow rate by manipulating the compressor power. The selected PI has the gains and , 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 and 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
is the water steam composition at the input of the cathode
is the water steam composition at the input of the cathode. The PI has the following gains:
and
.
5. 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®-Aspen Exchanger Design and Rating®-Aspen Dynamics®-Aspen Custom Modeler® 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.