In this paper, we study a solid oxide fuel cell (SOFC) controlled by a multi-input multi-output (MIMO) compensator, which uses the blower/fan power and cathode inlet temperature as actuators. The usable power of the fuel cell (FC) is maximized by limiting the air flow rate deliberately when an increase in power is demanded. Possible rate bounds on the cathode inlet temperature are also modeled. These bounds could represent the physical limitations (due to slow dynamics of heat exchangers) and/or a control concept for accommodating the power saving objective. Applying proper limits to the amplitude and rate of the actuator signals, and incorporating antiwindup (AW) techniques, can raise the net power of the FC by 16% with negligible effects on the spatial temperature profile.

## Introduction

Conventional electricity generation using fossil fuels is neither efficient nor suitable in terms of the high pollution it produces. Increasingly, fuel cells (FCs) are used in a variety of applications for electricity generation [1,2]. The ultralow emission, zero noise pollution of fuel cells, and their high efficiency help fuel cells as strong candidates for the new generation of power plants. Fuel cells are classified based on their efficiency, operating temperature, and type of electrolyte. A popular class of fuel cells with efficient long-term operation and fuel flexibility is solid oxide fuel cells (SOFC), which is characterized by its solid electrolyte. SOFCs, operating at high temperature, are being used in a wide range of mobile and stationary applications taking advantage of the high quality heat byproduct, useful for reformation [3,4]. One of the main challenges in commercialization of SOFCs is to be able to safely control the system in the presence of disturbances such as power demand variation, and/or nonlinearities such as actuator saturation [5,6]. One of the objectives of this technical note is to ensure power following while controlling the resulting changes in the temperature profile, in transient operation of solid oxide fuel cells. This would require a high performance control technique for minimizing the temperature gradients and, thus, thermal fatigue and material damage.

For high efficiency and low degradation of the fuel cell due to thermal cycling, the fuel cell temperature should remain fairly constant during operation. This has led to increased attention to SOFC thermal control [7,8]. To meet the demands of developing control strategies, in Ref. [9], a control-oriented multi-input multi-output (MIMO) nonlinear thermal model of the SOFC is developed and a temperature controller is proposed. In Ref. [10], a systematic approach to the multivariable robust control of a hybrid fuel cell gas turbine plant is presented, and the transient operation of the system is studied using an empirical model. In a recent work [11], a model of the SOFC in frequency domain is developed for system control design and stability. Such models are useful to understand the fuel cell interface with remote power systems for applications such as those found in aviation, terrestrial vehicles, and naval transport industry. In Ref. [12], a dynamic states estimator is designed to track and predict the behaviors of unmeasurable states inside SOFC using stochastic filtering algorithms. In Ref. [13], dynamic programming strategies are proposed for a hybrid system to achieve an optimal schedule while minimizing fuel consumption. Chaudhari et al. [14] investigate the transient behavior of a 25 kW SOFC and an internal combustion engine hybrid power system and apply model predictive control techniques in order to determine control parameters and set-points. Ansari et al. [15] propose a neural network model for prediction of SOFC performance based on the Levenberg-Marquardt back propagation algorithm.

The use of compressor/turbine generators with fuel cells in a hybrid format is common [5,16]. Here, for a 5 kW SOFC, we use a variable speed blower for cooling purposes and assume that there exists an external reformer. The complex dynamics of the SOFC system is modeled and advanced control techniques are applied in order to achieve stable and safe operation while maintaining system performance. The objective is to maximize usable power of a 5 kW SOFC while minimizing the spatial temperature variations. We start by assuming that a controller is designed for SOFC to primarily meet this thermal control objective. By relying on a high-performance MIMO controller, proposed in Ref. [17], the nonminimum phase like behavior caused by the fan can be addressed; however, a temporary surge in blower power due to the inertia of the blade will result (overshoot in blower power demand). Large overshoots in the blower power lowers the net power of the FC available for external use. That would necessitate conservative operation of the fuel cell. In order to address this issue, we propose to limit the power sent to the blower with a fictitious saturation bound. While all actuators have a limited range of operation due to physical characteristics and constraints, the approach here creates a limit that is lower than the physical limitations, and attempts to leverage the thermal mass of the fuel cell. By preventing the blower from using too much power, more of the power of the fuel cell can be made available for external demand, while keeping a high-performance controller, for tracking smaller changes or when the power demand is reduced.

As in any other physical system, actuation limitations can cause severe performance degradation (even system failure) [18]. Of course, actuators can be saturated both in terms of the size of the input generated and the rate at which the input can change, i.e., the magnitude and rate actuator saturation, respectively. In this paper, for the first time, the usable power of the fuel cell is maximized by artificially bounding the inlet flow rate, thus, avoiding the overshoots in the blower power. Antiwindup (AW) control techniques are then incorporated in order to guarantee stability and provide a satisfactory performance. The effect of AW on cathode inlet temperature, as the second actuator signal, is studied and possible rate bounds are modeled and compensated using antiwindup magnitude and rate augmentations.

The paper outline is as follows: In Sec. 2, the nonlinear model of a solid oxide fuel cell in coflow configuration is presented followed by the controller, which is designed assuming ideal actuation (Sec. 3). The issue of maximizing usable power of a SOFC by scheduling the actuator's amplitude bounds is discussed in Secs. 4. In Sec. 5, the concept of actuator rate saturation to accommodate the maximized power with acceptable performance is discussed. Section 6 combines both magnitude and rate limitation on actuators and proposes an antiwindup design using the peak-to-peak bound minimization approach to improve the performance of the system constrained by saturation nonlinearities. Simulation results are presented in Sec. 6.2, evaluating the effect of the proposed AW augmentations on system performance. Section 7 concludes the paper.

## Solid Oxide Fuel Cell Model

Solid oxide fuel cells are composed of a solid electrolyte separating the two electrodes, cathode and anode, which are fed by the air and fuel, respectively. In contrast to the typical heat engines running due to a temperature difference, the chemical potential difference is the driving force of the fuel cells. As shown in Fig. 1, the chemical potential difference requires the fuel and air to be apart in anode and cathode, respectively. Under proper operating conditions, negatively charged oxygen ions from the cathode chamber are transferred to the anode chamber, through the separating layer of electrolyte. The release of electrons in cathode then leads to the current and formation of water.

The main product of the fuel cell is the clean desirable electricity directly converted from the chemical energy avoiding the thermal cycle. The voltage available from a single fuel cell is approximately 1 V; thus, a stack of a large number of fuel cells is always of interest.

An integrated nonlinear dynamic model of a SOFC system developed with multiple subsystems is studied here. The model is the same coflow SOFC investigated in several previous efforts including [17,19] and has been evaluated with experimental data. The system is a typical planar coflow SOFC in the 5 kW scale integrated with a variable speed blower. The model has sufficient spatial resolution (i.e., along the flow channels) to capture the effects of power demand variations. A schematic of the SOFC model is shown in Fig. 2 representing four control volumes: cathode and anode gas channels, positive-electrode electrolyte negative-electrode (PEN), and the top and bottom interconnecting plates, which are the same due to periodic boundary condition assumption.

*P*

_{blower}is the motor power supplied to the blower and

*P*

_{impeller}is the loss associated with the impeller, which is given by

where *η* is the blower isentropic efficiency (85%), *γ* is the air-specific heat ratio, *R* is the universal gas constant, *T*_{amb} and *P*_{amb} are the ambient temperature and pressure, respectively. The dynamic nonlinear model of a coflow SOFC studied here is spatially discretized into five nodes along the flow direction (Fig. 2). Each node includes four control volumes: cathode and anode gas channels, PEN, and the interconnecting plate, which adds up to 20 control volumes. A list of the 60 states of the SOFC model is given in Table 2, having N_{2} and O_{2} as anode flow species and $CH4,CO,CO2,H2,H2O,N2$ as cathode flow species.

## Controller Design

One of the main challenges in transient operation of the solid oxide fuel cells is achieving load following with minimum risk of damage. Power variations lead to large changes in temperature profile of the SOFC, which result in undesirable thermal fatigue and consequently serious damage to the structure of the fuel cell [20]. Thermal control design aims at minimizing the spatial temperature variation throughout the SOFC, while maintaining a reasonable performance in load following.

*u*and

*w*. The linearized reduced order plant is then given by

The *δ* variables represent the variations from the nominal operating conditions. The state vector $\delta xp\u2208\mathbb{R}np$ is the difference between the current state of the system from those at nominal condition, while $\delta up\u2208\mathbb{R}nu$ is the control commands that would be added to the nominal values of the input to obtain the commands to the actuators. Other variables, i.e., the measurement outputs (sensors) $\delta y\u2208\mathbb{R}ny$, and the performance outputs $\delta z\u2208\mathbb{R}nz$, similarly denote variations from the nominal conditions. The exogenous input $\delta w\u2208\mathbb{R}nw$ is the reference signal, representing the change in power demand.

The block diagram of the SOFC control system is shown in Fig. 3. Table 3 summarizes the input and output signals associated with the SOFC model. According to the table and figure, cathode inlet temperature and blower power are inputs to the plant $up1$ and $up2,\u2009up=[up1\u2003up2]\u22a4$, respectively. The anode outlet temperature, plate temperatures at the first, middle, and the last nodes, and the blower shaft speed are the measurement outputs, *y*. The electrolyte temperature at each of the five nodes are the performance outputs, *z*. An outer proportional–integral–derivative loop uses the power tracking signal to obtain the voltage needed and the nominal voltage is then used as the disturbance *w*.

Actuators | Cathode inlet temperature, $Tcath,in$ |

Blower power, P_{blower} | |

Measurement outputs | Anode outlet temperature |

Plate temperature, nodes 1, 3, 5 | |

Blower shaft speed | |

Performance outputs | Electrolyte temperature, nodes 1–5 |

Disturbance | Fuel cell voltage |

Actuators | Cathode inlet temperature, $Tcath,in$ |

Blower power, P_{blower} | |

Measurement outputs | Anode outlet temperature |

Plate temperature, nodes 1, 3, 5 | |

Blower shaft speed | |

Performance outputs | Electrolyte temperature, nodes 1–5 |

Disturbance | Fuel cell voltage |

*C*, with same order as of the plant, is presented by

## Actuator Magnitude Limitations

The fuel cell power is prone to large variations often due to the grid (or microgrid, for example) demand. Therefore, actuator saturation is likely, particularly for a high performance controller. Addressing the problems created by actuation saturation is one of the challenges in FC controller design. As discussed in Ref. [22], fast load following is possible under the assumption of ideal actuation. Under realistic actuation, however, power following degrades and the stability and performance are no longer guaranteed. The main actuation problems are the power needed by the blower to deliver the required cathode air flow, and the proper temperature needed at the cathode inlet. If the power demand change is large, it might lead to saturation problem for the blower power. On other hand, due to the inertia of the blower, the MIMO controller designed for the integrated SOFC and the blower model can cause relatively large overshoots in the blower power and consequently limit the available net power of the FC. Since the power of the blower is provided by the fuel cell, any power not used by the blower can be added to net power supplied by the FC. While we can reduce the overshoot, by using a less aggressive controller, we would like to avoid this solution particularly when power demand is lowered and the FC temperature is prone to a significant drop. This might lead to lower efficiency in ion transport and large temperature variations.

Alternatively here, in order to avoid the overshoot and increase the range of power available, we use an aggressive controller and rely on artificial saturation (enforced by software) that limits the power sent to the blower. This acts on the controller as a standard saturation bound, which can result in performance degradation or instability. The next step is then to develop antiwindup protection schemes, which help maintaining the stability as well as an acceptable performance for the fuel cell under such actuator saturation.

The power demand profile used in simulations in this paper corresponds to ±15% variations in power demand. Figure 4 shows the demand profile composed of a 15% decrease in the power at time *t* = 11,000 s, and then 30% of the nominal power increase at time *t* = 12,000 s. The times chosen for changes in power are far apart from each other and the start-up time, to focus on the dynamics of power set-point change. The simulation results are presented for the nominal conditions, given in Table 4, satisfying the temperature range often used as the most common electrolyte average temperature.

### Blower Power Enforced Saturation.

*u*

_{lim}as a known positive constant, and $up2$ to be the blower power. The inputs to the plant

*u*are thus modeled as

_{p}_{;}thus, the actuator signal is

*q*

### Antiwindup Design With Magnitude Bound on Blower Power.

*ν*

_{2}is added only to the second row of the input vector

*u*, which corresponds to the blower power as from Eq. (10)

*w*and

*q*as input signals. Given the linearized SOFC state-space model (7) and the unconstrained controller (8), the closed-loop system with magnitude antiwindup gains can be written as

In order to guarantee the stability of the closed-loop system (16) and establish a performance bound for the AW design, an upper bound for the $L2$ (or Energy) gain *γ*, from the disturbance signal *w* to the performance output *z*, is minimized and the stabilizing AW gain Λ is obtained. The minimization is constrained by a linear matrix inequality ensuring the stability and performance bound *γ* (details omitted due to space limitation). While the algorithm presented is to bound the energy of the performance output (by a fix multiple of the energy of the reference input), other techniques that bound the peak norm (i.e., peak-to-peak or energy-to-peak) are straightforward and can be implemented in a similar fashion with ease [24]. Indeed, we use minimization of the upper bound for peak-to-peak gain in Sec. 6.

While the blower power for the nominal condition is 0.288 kW (from Table 4), increasing the overall power level requires higher air flow rates and thus higher blower power. Furthermore, the overshoot of the blower should also be taken into account. As shown in Fig. 6, the peak blower power (for both overshoot and steady operation at higher power levels) reaches close to 1 kW. Therefore, the maximum net power, which can be achieved form the fuel cell, is limited to 4 kW. Assuming a ±15% variation in the nominal power demand, the nominal value of the fuel cell net power should then be approximately 3.5 kW, in order to respect the 5 kW upper bound of the fuel cell capacity, i.e., $Pnet,max=3.5$ kW. In order to increase the maximum net power available by the fuel cell, the idea of the enforced saturation (bounding the blower power intentionally) is applied, and AW techniques are used to ensure stability and performance. The steady-state value of the blower power (0.66 kW) could be selected as the magnitude bound in order to avoid the excess power in overshoot. Due to the thermal mass of the FC, cutting the power for the short duration of the overshoot would not change temperature profiles significantly, while a %7 increase in the net power could be achieved.

A lower power bound results in more power saving, particularly if the higher power demand is not permanent and a decrease in tracking profile is expected. Therefore, in Fig. 6, the limit of 0.4 kW is selected for the blower power saturation element (recall the nominal value is 0.288 kW).

*u*

_{lim}= 0.4 kW on

*P*

_{blower}. This level implies a 16% increase in the FC available net power. Therefore, by putting an aggressive limit on the blower power and taking advantage of the AW techniques, the same SOFC can be used for 16% higher power demands.

Without antiwindup, as shown in Fig. 8, the 16% additional power made available is at the cost of a large temperature increase while the signal is saturated. By limiting the blower power, not enough air is sent through, causing a large temperature rise along the cell, especially at the outlet. These changes in the temperature profile can have negative impact on the FC and may lead to degradation and thermal fatigue. However, taking advantage of the proposed AW techniques, the error can be reduced to the negligible amount of less than 2 K.

Figure 9 shows the variation of the cathode inlet temperature, the other actuator of the fuel cell system, when it is not subject to any saturation constraint. This figure shows that in order to address the blower power saturation, and reduce the error in the outlet temperature, the augmented controller avoids further heating up of the inlet using the other actuator signal ($Tcath,in$). As shown in Ref. [17], in coflow FC, without saturation constraint, both the air flow rate and cathode inlet temperature are increased in response to higher power demand. The elevated airflow rate is aimed at reducing the average temperature while the higher inlet temperature ameliorates potential thermal gradients along the cell.

In order to accommodate the instantaneous drop in inlet temperature observed in Fig. 9, an actuator with high rate of change is required. However, in the absence of ideal actuation assumption, the rate at which $Tcath,in$ can change is limited by the mechanism used (e.g., heat exchanger, mixing chamber). Under operational conditions, this may lead to actuator rate saturation, causing severe performance degradation. The rate antiwindup design introduced in Sec. 6.1 is used as the remedy, since an artificial rate bound can address this sudden drop in cathode inlet temperature, thus avoiding large temperature gradients. The objective is to extend AW protection schemes, which help maintaining the stability as well as an acceptable performance for the fuel cell under actuator both magnitude and rate saturations.

## Actuator Rate Saturation

The inlet air temperature is controlled by bypassing the air through a heat exchanger or a mixing chamber with possibly slow dynamics. As the SOFC is a high temperature fuel cell, there might not be a real magnitude bound on the temperature control signal for the cathode inlet. However, the rate at which this temperature can change, in order to satisfy the control commands, can be limited by the heat exchange mechanism or transport delays.

### Rate Model.

In order to study the effects of actuator rate saturation on the fuel cell performance, we first need to have access to the rate signal (not typically available) and apply rate limits. One common approach to model, the rate signal is to insert a first-order filter with gain *K* in the forward loop connecting the controller's output to the plant's input. As shown in Fig. 10, a first-order circuit with a saturation element is added before the actuator. As a result, the actuator signal is guaranteed to be rate bounded.

The gain *K* needs to be selected large enough in order not to affect the dynamics of the original system. Figures 11 and 12 illustrate the effect of rate model with different *K* values on inlet and outlet temperatures, respectively. In fact, gain *K* controls the delay in system response as this is a first-order filter that may influence the FC behavior. Higher *K* corresponds to smaller delay and thus faster response compared to the time constants of the original system. For relatively small values of gain *K*, as shown in Fig. 11, the delay is observed in the inlet of the fuel cell. For $K\u22650.1$ temperature profiles throughout the cell remain unchanged after adding the rate model. Therefore, for the rest of the simulations in this paper, *K* = 0.1 is selected for the gain of the rate model. The signals are available for the antiwindup loop since the loop in Fig. 10 is a part of the compensator.

### Cathode Inlet Temperature Rate Saturation.

Figure 13 shows the block diagram of the SOFC control system with the first-order model added, having the first actuator subject to rate saturation. For the first plant input ($Tcath,in$), we require $|u\u02d9p1|\u2264r$ with $r$ as a known positive constant with unit $K/s$. This model guarantees the cathode inlet temperature signal to be rate bounded.

Figure 14 shows the effect of different rate saturation bounds on the profile of the cathode inlet temperature. When an increase in power demand is requested, the unconstrained controller increases both inlet temperature and the airflow rate in order to minimize the spatial temperature gradients along the cell (dashed curve). However, as shown in the figure, a limited rate of change of $Tcath,in$ causes a lower slope in the inlet temperature rise.

Figure 15 shows the effect of different rate saturation levels on profile of the anode outlet flow temperature. As the rate bound gets tighter, inlet temperature rises slower; thus, a drop in temperature at the outlet is expected. Compared to the cathode inlet temperature (Fig. 14), the anode outlet temperature, Fig. 15, shows less sensitivity to rate saturation since it is located far from the point of restriction and the effects of the bound become less profound toward the end of the cell. The decrease in inlet temperature and the increase in outlet temperature may cause a minor increase in spatial temperature gradients from nominal conditions, particularly at the nodes near the inlet. To address this, we add a rate limited AW with details discussed in Sec. 6.

## Actuators Magnitude and Rate Saturation

where $K=diag(K1,\u2026,Knu1)\u2208\mathbb{R}nu1\xd7nu1$ and *x _{I}* is the integrator's state. In order to minimize the negative effects of rate saturation on fuel cell performance and guarantee the system stability, we rely on antiwindup design for both magnitude- and rate-bounded actuators.

### Antiwindup Design.

*P*

_{blower}) is magnitude bounded, the antiwindup term

*ν*

_{2}is added to the entire input vector

*u*(see Eq. (13) for comparison). Therefore

*x*is defined as a new state variable with dynamics

_{I}*u*

_{1}is the first row of

*u*in Eq. (20). The augmented closed-loop system with state vector $x=[\delta xp\u22a4\u2003xc\u22a4\u2003xI\u22a4]\u22a4$, and

*w*and

*q*as input signals in then presented as

*and Λ*

_{m}*are then obtained by solving a convex optimization problem subjected to performance and stability constraints (inspired by the result presented in Ref. [25]). Further details are provided in Appendix B.*

_{r}### Results.

In this section, the behavior of an SOFC under the magnitude and rate actuator bounds discussed in Sec. 6 is studied through simulations. The objective is to achieve a %16 rise in the FC net power by artificially bounding the inlet flow rate, considering the possible physical bounds on the rate of change of the inlet temperature due to slow dynamics of heat exchanger/mixing chamber. The reasonable bound of $0.1\u2009K/s$ is selected as the rate limit with gain *K* = 0.1 based on the discussions in Sec. 5.1.

Figures 18 and 19 show the performance of the optimized system with maximum net power, compared with the original system. Figure 19 shows the temperature gradients from the nominal conditions at each node along the cell. As discussed in Sec. 2, the cell is discretized into five nodes along the flow direction: node 1 located at the inlet toward node 5 placed at the outlet. The required magnitude and rate bounds are applied to the actuators in order to achieve a %16 increase in the net power. The maximized power is available at the cost of at most $5\u2009K$ increase in temperature gradients, from nominal conditions, along the cell.

As an alternative approach, in Ref. [26], a magnitude lower bound was used for cathode inlet temperature, avoiding it to drop significantly due to AW compensation. Here, the rate limit on the inlet temperature, which could be due to physical constraints and/or control design purposes, addresses this issue and keeps the temperature gradients within an acceptable interval. The overall performance here is superior to the one achieved in Ref. [26]. As an example, the rate bound eliminates excessive temperature gradients along the cell.

## Conclusions

The power of the cooling fan of an SOFC is intentionally limited, in order to increase the fuel cell net power available for external use. Possible rate limits on cathode inlet temperature, as the second actuator signal, are also simulated and their effects on the performance of FC are studied. The existing MIMO controller used for spatial temperature control is then augmented using antiwindup techniques to recover the performance of the unconstrained system. A 16% increase in the FC usable power is achieved by applying suitable actuator limits and incorporating antiwindup techniques, only at the cost of less than $5\u2009K$ variations in temperature gradients.

## Funding Data

NSF Grant No. CMMI-1461583.

## Nomenclature

### Appendix A: Antiwindup Design for Multi-Input Multi-Output Systems With a Magnitude Bounded Actuator

where *B*_{22} and *D*_{122} are the second columns of *B*_{2} and *D*_{12}, and $Dcy2$ is the second row of *D _{cy}*. The system matrices are slightly different from the general magnitude AW design due to the fact that here only a subset of input signals are subject to saturation (only the blower power,

*u*

_{2}, and not the cathode inlet temperature,

*u*

_{1}).

**Algorithm.**[Magnitude AW for MIMO systems with partially bounded actuators: Energy gain approach] [22,26]. Consider the plant and the controller introduced in Eqs. (7), (8), and (10), as well as the magnitude limit $m$. Given any solution to the optimization problem

*γ*from

*w*to

*z*.

Inequality (A3) provides the stability condition for the saturated system with antiwindup gain Λ using the standard Lyapunov stability approach.

### Appendix B: Antiwindup Design for Multi-Input Multi-Output Systems With a Magnitude and a Rate-Bounded Actuator

Here, we assume that a possibly conservative estimate of disturbance signal *w*(*t*) is known. Thus, peak-to-peak approach is applied to the current fuel cell saturation problem.

**Algorithm**. [Magnitude and Rate AW for MIMO systems with partially bounded actuators: Peak-to-peak approach] Assume that this system is only exposed to peak-bounded disturbances with known upper bound $wmax\u2208\mathbb{R}$, i.e., $w\u22a4(t)w(t)\u2264wmax2$. Given the saturation model (19) with the magnitude and rate limits $m$ and $r$, let us assume that for a given $0<\alpha <(|Re(\lambda min(A))|/2)$, there exists a solution for

Parameter *α* can be selected by performing a typical line search. The antiwindup gains satisfying the stability and performance are then given by $\Lambda m=XmMm\u22121,\Lambda r=XrMr\u22121$.