Optimal Load-Tracking Operation of Grid-Connected Solid Oxide Fuel Cells through Set Point Scheduling and Combined L 1-MPC Control

An optimal load-tracking operation strategy for a grid-connected tubular solid oxide fuel cell (SOFC) is studied based on the steady-state analysis of the system thermodynamics and electrochemistry. Control of the SOFC is achieved by a two-level hierarchical control system. In the upper level, optimal setpoints of output voltage and the current corresponding to unit load demand is obtained through a nonlinear optimization by minimizing the SOFC’s internal power waste. In the lower level, a combined L1-MPC control strategy is designed to achieve fast set point tracking under system nonlinearities, while maintaining a constant fuel utilization factor. To prevent fuel starvation during the transient state resulting from the output power surging, a fuel flow constraint is imposed on the MPC with direct electron balance calculation. The proposed control schemes are testified on the grid-connected SOFC model.


Introduction
The globally exploding capacity of renewable energy during the past decade makes a remarkable contribution to the conservation of fossil-fuel energy resources and an increasing number of distributed generations (DGs).However, the intermittent energy sources, such as solar energy and wind energy connected to the power grid, brings severe shortage of peaking regulation capacity, which has become the major factor that limits the further promotion of these renewable energies.To this end, seeking for new peaking regulation units have become a renewed issue in the field of power and energy engineering in recent years.It is reported that some large scale thermal power plants are incorporated into the peaking units recently through a series of technical reformation [1,2].However, the variable load operation would have a negative impact on the efficiency of thermal power plants and brings extra cost to the maintenance [3].In the recent years, the academia realize that electrolytic hydrogen making combined with hydrogen storage might be a promising solution to address the challenge brought from renewable energies [4,5].As a complementary technology for hydrogen energy storage, solid oxide fuel cells (SOFCs) has drawn a lot of attention these years for its high efficiency of energy conversion and is treated as a candidate for peaking power plant when being connected to a power gird due to its capability for fast power response.
Omitting the preheaters and the fuel processers, Figure 1 gives a brief schematic diagram of a grid-connected SOFC that is discussed in this paper.The fuel (hydrogen) and the oxidant (usually oxygen from compressed air) are fed into the anode and cathode, respectively.The electrochemical reaction takes place at the electrolyte and produces electric potentials.The reaction is When the circuit loop is closed, the reaction proceeds to output electric power to the external load.In practice, SOFCs are usually connected to a micro power grid through voltage source invertors (VSI) for DC-AC transformation.In the grid level, the dispatching system is working to generate load demand signals to the power sources that are connected to the grid.These power sources make response to the load demands, thus the energy are balanced in the power grid and the frequency are maintained [6].Consequently, for grid-connected SOFCs, its output power should track the unit load demand from the dispatching system.With this operation mode, the control problem considered for a grid-connected SOFC lies in: 1.
As peaking power plant, the SOFC system should have fast load variation ability to compensate the power fluctuation from intermittent energies in the gird.In other words, the peaking SOFC systems are expected to have a fast response to the unit load demand.

2.
When the SOFC operates in a steady state, the system efficiency should be optimized in order to exploit for the maximum profit.An explicit for the efficiency is the ratio of the output electrical power to the chemical power that is released by the fuel, which implies that with a specific output power, the maximum of the efficiency is equivalent to the minimum of the heating power that is wasted with the cells' internal resistance.

3.
The dynamic of SOFC system takes on a strong nonlinearity under various operating conditions.Since the performance of linear controllers would deteriorate seriously when the system deviates from the nominal design condition, the control system should have the ability to handle the nonlinearity.4.
The system constraints, referring to the limitations of system variables associated with the physical property and safety requirements should be fulfilled during both the transient process and steady-state operation.Among the constraints of SOFC system, fuel utilization (FU) is the most critical one for the system safety.In the first place, an overused fuel condition (FU > 0.9), namely fuel starvation, should be strictly avoided.The reason is that the fuel provides a reducing atmosphere for the electrodes under the operational condition of high temperature inside the SOFC.Lack of fuel gas would break the reducing atmosphere.In this case, the anode materials would be oxidized, triggering a permanent damage of the SOFC [7].Low FU is accepted in the transient process, but a long-term low FU operation dramatically increases the fuel cost, thus causing a low efficiency of SOFC system.Most of the existing research literatures suggest that the FU should be controlled at 0.8 during the steady-state operation and within the range of 0.7 to 0.9 during the transient state for a comprehensive consideration both on safety and economy [8].
The intention of the paper is to put forward an integrated control system design attempting to meet the above requirements for a grid-connected SOFC system.In industrial processes, a classical solution to optimal operation is to adopt a hierarchical strategy in the control system design.The control system is divided into two levels, an upper level, with set point scheduler, and a lower level, with tracking controller.The set point scheduler performs an optimization on the basis of a rigorous steady-state model of the plant and feeds the lower level controller with setpoints corresponding to the optimization results [9,10].It is of course that a complicated model with as more detail as possible is preferred for better optimization results.On the other hand, an excessively complicated model is disadvantageous for controller design, resulting in increasing cost for the implementation of controller and decreasing system reliability.The advantage of hierarchical control structure is that it separates the optimization problem and control problem.The separation of the two problems regarded in optimal operation enables the optimizer and the controller to be designed individually without being concerned by incorporation of each other.Benefiting from the merits discussed above, this hierarchical control is naturally treated as a candidate employed for achieving optimal operation of SOFC.However, as far as the authors' knowledge, such study conducted in a SOFC operation study is still limited.
Making a general survey of mainstream control algorithm for tracking controller in the lower level, Model Predictive Control (MPC) appears to be the most powerful and convenient one to deal with system constraints [11][12][13][14].Some studies on the application of MPC to a SOFC control problem have been performed previously [8,[15][16][17].These designs used varieties of nonlinear models in the MPC formulation to accommodate the nonlinearities.However, extensive online computation for nonlinear MPC limits its suitability for a time-sensitive system as grid-connected SOFCs.Previous works proposed an Active Disturbance Rejection Control (ADRC), which cooperated with MPC to handle the nonlinearity with acceptable computation [18].Despite a given voltage being followed, the FU was not able to be maintained at the expected constant.It is believed that a direct modification of this control algorithm for maintaining a constant FU is difficult, due to the inability of the ADRC algorithm for multi-input-multi-output (MIMO) systems.To meet the requirement for voltage control and constant FU maintaining simultaneously, L1 Adaptive Control (L1AC) will be employed to work in combination with MPC, which will be referred as L1-MPC in this paper.The L1 Adaptive Control is a robust adaptive control algorithm developed from the model reference adaptive control (MRAC) algorithm.In the philosophy of the L1AC algorithm, it treats the nonlinearity, disturbance, and unknown parameters of the system as a lumped uncertainty term in the model.The lumped uncertainty term is estimated with a fast adaptive law and a compensation term is included in L1AC's control law.With the fast adaptation mechanism, the closed-loop system with L1AC approaches a given reference model.In [19][20][21], theoretical analysis has been given to show the guarantee of the transient error between the closed-loop system and reference model.Although the proposer of L1AC mainly shed light on the state feedback algorithms and the existing output feedback L1AC algorithm is only for single-input-single-output (SISO) system [19], it is convenient to extend the SISO L1AC algorithm to a MIMO one.
In this paper, we present a hierarchical control strategy for grid-connected SOFC.In the upper level, set points of output voltage and current corresponding to unit load demand are obtained through solving a nonlinear optimization problem for the maximum efficiency of SOFC in steady-state operation.Steady-state analysis of a tubular SOFC is carried out in advance to formulate the optimization problem.In the lower level, a L1-MPC tracking controller with a combination of MIMO output feedback L1AC and MPC is put forward to steer the output variables to their optimal set points, handling all of the nonlinearities and constraints.To ensure the avoidance of fuel starvation, additional protection logic is designed afterwards.
The rest of this paper is organized as follows.In Section 2, several dynamic SOFC models are first reviewed to choose the most proper dynamic model for the control study in this work.Then, steady-steady analysis of a tubular SOFC is presented.Section 3 establishes the hierarchical control architecture with set point scheduler and the L1-MPC based tracking controller.Simulation results are given in Section 4. Finally, conclusions and future works are presented in Section 5.

Review of SOFC Structures and Models
SOFCs are classified into two major types by their geometrical shape, planar, and tubular.Tubular SOFCs are usually appraised to be more superior to the planar ones for their advantages in sealing and structural integrity.However, the two types of SOFCs are equivalent in fundamental principles.Without the loss of generality, a tubular SOFC is chosen to be discussed in this paper.Figure 2 shows the structure of a tubular SOFC.The air (oxidant) is supplied through the air supply tube (AST) into the cell, and then flows past the annular channel of cathode reversely to the fuel, which flows over the anode channel.The electrolyte is sandwiched between the two electrodes.A reliable dynamic model with proper complexity is essential for the control design of SOFC.Most of the state-of-the-art SOFC papers are based on a benchmark dynamic model that is presented in [22].This model appears to be very popular in the last decade for its good feasibility for a large scale simulation and its acceptable accuracy on the system dynamics.An important deficiency with the model is that its thermal characteristics are omitted.It has been acknowledged that this simplification is reasonable in studying the dynamics of SOFC since the thermal process is very slow and have little influence on the short-term dynamics.However, the temperature of SOFC is a crucial factor to not only the electrodynamics potential, but also the internal resistance, which are both closely associated with the optimal operation of the SOFC system.Therefore, the benchmark model in [22] is no longer appropriate for the research in this work.
Based on a comprehensive studying of the SOFC's principles, Wang and Nehrir established a model that includes both the electrochemical and the thermal characteristics of a tubular SOFC [23].The most significant feature of this model is that it contains detailed mechanism of internal resistance and output voltage.Moreover, with several reasonable assumptions, the heat transfer analysis inside the tube structure is based on the lumped parameter method.These features make the model adequately accurate with reasonable complexity.

Steady-State Analysis
Motivated by the dynamic model presented in [23], the analysis of SOFC under steady-state is with the mass conservation, thermal balance and electrochemical properties, respectively based on the following assumptions: the air at cathode is with large stoichiometric quantity, and 3.
the pressure gradient of the gases is uniform along the channels.
The notations in this section are listed in Table 1.

Mass Conservation
According to Assumption 1, the following equations are obtained: with Faraday's Law, the material balance equations for H 2 , H 2 O, and O 2 in the steady-state are, In the flow channels, the following mass flow stands The mass flows for the gases at outlet of the channels are The fuel utilization is defined by

Thermal Balance
The chemical reaction in a SOFC is carried out in the electrode-electrolyte tier (EET) with the production of heat.The electrode-electrolyte tier is cooled by the fuel and air flows through the channel.In the meanwhile, some of the heat transfers to the air supply tube by means of radiation.A cross-sectional profile of a tubular SOFC and the heat transfer are shown in Figure 3.Note that only half of the profile is drawn here because of the symmetry of the profile.The heat transfer analysis is based on the lumped parameter method.Five lumps are considered here, which are the fuel in the channel, the electrode-electrolyte tier, the air in the annulus, the air supply tube, and the air inside the tube, respectively, from top to bottom in the figure.In the manner of lumped parameter systems, the temperature difference inside each lump is neglected.With the lumped parameter method, the average value of the temperature at two ends is used to represent the gas temperature in the channels: T air,ann = T air,AST,ann + T air,outlet 2 ( 14) The Cell Tube: For the EET, in steady-state the heat produced is completely dissipated in different ways.The following equation stands in this case: The left hand side of the equation is the heat production in the cell, which is the difference between the total chemical energy that is released by the reactant and the actual electrical power output.The right hand side is the heat transfer with the cell in different ways.The heat transfer with other parts of the cell is listed as follows.
The fuel: Air between cell and AST: AST: Air in AST: The heat transfer terms in ( 16) ~( 20) are calculated through heat transfer equations as: Q cell,ann,conv = h cell A cell,inner (T cell − T air,ann ) Q AST,ann,conv = h AST,outer A AST,outer (T air,cell − T AST ) Q AST,inner,conv = h AST,inner A AST,inner (T AST − T air,AST ) Q air,ann, f low = M c M mw,air C air (T air,outlet − T air,AST,ann ) Q air,AST, f low = M c M mw,air C air (T air,AST,ann − T air,inlet )

The Electrochemical Relations
The potential of fuel cells is subject to the Nernst equation: where, E 0 0,cell is the standard electrode potential, the value of which depends on the species of reactant.For the reaction in SOFC, the reference potential is 1.229 V. E cell is actually the open-circuit voltage of the fuel cell.For a SOFC with electrical loads, the output voltage is less than the E cell due to the varieties of voltage drop including activation drop, ohmic resistance drop, and concentration drop.Therefore, the actual output voltage is calculated as The voltage drops are calculated as [23]: Note that the model presented above is for single cell.Since the electric potential provided by the electrochemical reaction in a single fuel cell is very small, numbers of the fuel cells are usually cascaded together (referred to as fuel cell stacks) to acquire a proper voltage level in the practical application.

Control System Design
The total control strategy in this work is proposed to have a hierarchical architecture, with which the control system is generally divided into a supervisory level in the upper and an underlying level in the lower, as shown in Figure 4.In the upper level is a set point (SP) scheduler.The scheduler receives unit load demand from the power grid dispatching system.With the unit load demand, online optimizations are carried out by the SP scheduler with the steady-state model of the fuel cell to find a set of target for voltage and current under which the corresponding steady states of the SOFCs satisfies the load demand with optimal economic performance.To deal with the mismatch between the steady-state model and the real plant, a modification mechanism will be working with the SP scheduler to make the acquired target as precise as possible.
In the lower level, a SOFC controller is working with a voltage source inverter (VSI).With the set point for voltage being provided from the SP scheduler, the SOFC controller steers the SOFC's output voltage to the set point by manipulating the inlet fuel flow to the anode and the inlet air flow to the cathode.At the same time, the fuel utilization should be maintained within the tolerant range (0.7 to 0.9) in transient states and be driven to the proper value (0.8) in steady states.For this purpose, a combined L1-MPC controller is put forward to generate the control signal for the manipulated fuel and air flows, handling the nonlinear dynamics in the SOFCs and taking care of the system constraints with acceptable online computation.The cooperating VSI is integrated with its own controller.The task for the VSI is to control the current flow to the power grid.The specific control algorithm of VSI is out of the scope of this paper and will not be discussed here, but it is worth pointing out that the time constant of the VSI's dynamics (around 0.1 s) is much smaller than that of the SOFC's to a new steady state [24].Therefore, the transient process of the VSI will not be considered in the SOFC controller design.The algorithm discussed in detail for the SOFC controller will only include the upper level set point scheduler with modifier and the level layer L1-MPC combined controller (as marked with dashed box in Figure 4).This part is redrawn in Figure 5 with the signal flow components listed in Table 2.The specifics of the L1-MPC controller will be discussed in this section.

Set Point Scheduler
In the supervisory level, the set points of current and output voltage is acquired from an online nonlinear optimization module using the steady-state model given by ( 1)~(34), and the unit load demand given by the grid dispatching system.With the unit load demand, the working condition of SOFC satisfies and the proper fuel utilization F u,re f = 0.8 (36) The goal of the optimization considered in this work is to minimize the internal power dissipation of the SOFC under given output power.With the notations that y s = [V out , I] T , , the mathematical formulation of the optimization problem is where Φ(y s , x s , u s ) = 0 is the equality constraint constituted by the steady-state model in Section 3. The acquired V out,re f are given to the combined L1-MPC controller together with the expected fuel utilization F u,re f for the regulation of the SOFC, while I re f is fed to the VSI to adjust the output current of the SOFC.For practical SOFC systems, the optimal operation cannot be achieved by simply solving (37) due to the model errors.The model errors come from the structure mismatch during modeling and parameter perturbation during operation.Since the model errors are not avoidable, the nominal solution of (37) is never likely to be optimal in reality.To deal with this problem, an approach used in earlier time is to re-identify the parameters online using the plant measurements in steady states.
The approach is now referred as "two-step approach" in the literatures [25].However, the two-step approach is usually not effective in practice as the parameter identification is easy to converge to a local optimization.Another approach is to solve the optimization problem with a modifier, which is substantially a lumped estimation of model errors, added to the constraints to correct the model values [26].Owing to the avoidance of complex identification and the ability of acquiring an acceptable optimization result in practice, the modifier approach is now widely used in so-called real-time optimization in industry and it is applied in the SOFCs discussed in this paper.
The amendatory optimization problem with modifier is formulated as: where υ(k) denotes the modifier, which is updated with where K υ is a designed parameter that determines the convergence rate of modifier υ.Notice that the elements in vector x are internal states of SOFC and are difficult to measure.In this case, we replace x(k) with x re f (k − 1) and attribute the corresponding error to the lumped modifier υ(k).Before modifier updating, steady-state detection should be conducted to capture the input/output data used for modifier updating.The steady-state detection used to be a much studied topic in process control in the late 20th century and many approaches have been introduced in that period (see [27][28][29]).For gird-connected SOFCs participating in peaking regulation, the output power should change frequently in order to compensate the power deficit induced by intermittent energy sources.In this case, the system will operate at non-steady states for most of the time.In consideration of this characteristics, we apply a method that is given by Cao and Rhinehart [29].The applied approach provided a way to simply use a sequence of filtered data, preventing the requirement for considerable data storage and user expertise.Denoting X as the indicator variable that is used to detect the steady state, the procedure for steady-state detection is carried out iteratively as follows: Step 1: Calculate a first-order filtered term of X by: Step 2: An estimation of mean square deviation is calculated by: Step 3: Use the following equation to get estimation to the squared differences of successive data: Step 4: Taking the ratio of the two estimates determined by ( 41) and (42): In the equations above λ 1 , λ 2 , λ 3 are parameters that satisfy 0 < λ 1 , λ 2 , λ 3 ≤ 1.The criterion for steady state is |R i − 1| ≤ α, where α is a chosen threshold value for the judgment of steady states.
In this work, the output power is chosen as the indicator for steady-state.The reason is that the output power is determined by both the current and voltage.Further, the voltage is determined simultaneously by a large amount of the SOFC parameters, thus it is sensitive to the unsteadiness of the system dynamic.In this sense, the output power might be the most comprehensive and a reliable variable to indicate the achievement of steady states.
When a steady state is detected, the modifier update mechanism with (39) is conducted using corresponding input and output data of the system.The steady state detection procedure can be deployed in an individual module separate to the set point scheduler with higher sampling frequency, or integrated with the set point scheduler.The whole control scheme with the modifier is presented in Figure 4.

Combined L1-MPC Controller
For SOFC tracking controller design, the outputs and inputs of the control problem are denoted as y = [y 1 ,y 2 ] T = [V out ,F u ] T and u = [u 1 ,u 2 ] T = [M H2 ,M air ] T , respectively.Then, the controlled system is considered as: where G(s) is a 2 × 2 transfer function matrix, δ(s) denotes the lumped uncertainty terms, which contains all of the nonlinearity associated with the system, and (s) denote systems, variables, or signals in the s-domain with Laplace transformation.
The system (44) can be rewritten as, where for i = 1, 2, where G ii (s) denotes the (i,i)-th element of G(s), σ i (s) is generalized uncertainty term that satisfies, for i = 1, 2, and δ i (s) denotes the i-th element of δ(s).
The corresponding time-domain form of σ i (s), σ(t) can be written as, The discussion is based on the following assumptions: 1.The diagonal elements of transfer function matrix G(s) are positive real transfer functions.2. There exist constants L > 0 and L 0 > 0, possibly arbitrarily large, such that hold uniformly in t.

The variation rate of σ(t) is bounded,
.
Then, it follows from the basic procedure of L1 adaptive controller, the following output predictor is constructed: for i = 1, 2, where G ii (s) denotes the (i,i)-th element of G(s), and σi (s) is the Laplace form of the estimation σi (t), which is updated by the adaptation law: . σi (t) = ΓProj( σi , −sign(g ii ) y i ) where Proj represents the projection operator [21], y i = ŷi − y i , g ii is the DC gain of the i-th diagonal element in G(s), Γ is the adaptive gain.Then, the control law is presented as where C i (s) is low pass filters with unit DC gain satisfying such that u * (kT MPC ) is from the solution of the following optimization problem: where A, B, and c, are the system matrices of discrete state-space model corresponding to G(s) with sampling time T MPC , u min , and u max are the input constraints, y min and y max are the expected output constraints, µ re f is a reference input corresponding to the output reference y re f with the linear model is the L1 adaptive term at the sample time k, y(k − 1) denotes the measurement of y at the last sample, L is a matrix that keep all of the eigenvalues of matrix (A + L • c) inside the unit circle.
Equation ( 56) is a standard formulation for output feedback MPC with a state observer and can be easily transformed to a quadratic programming (QP) with a few steps of algebraic operation as follows.
The first two equations of (56) makes the prediction of the system output: . . .
In (57), M is a tunable parameter called controlled horizon in MPC theory, with which, the predicted system input u (k + i|k) keeps invariant when i ≥ M.Then, the optimization problem (56) can be transformed to the following QP form with (57): min where, Note that to avoid the FU from exceeding the upper limit and triggering fuel protection, the lower bound for fuel flow (u 1 ) should be updated at each sampling time with the reference current at the moment: where F upper u is the upper limit for FU and is usually set to be 0.9 or is slightly lower than 0.9 for safety.
The optimization problem (56) is solved online at each sampling time with the receding horizon framework.In this manner, the first term of the solution sequence U(k) is applied on the system, i.e., u * (kT MPC ) = u(k|k) (60) Theorem 1.For plant (44) applied with control law in (52), and defining a reference system as with control input Then, the following bound of system transient performance holds: where θ m,i = 4σ 2 b,i + 4 λ max (P i ) λ min (Q i ) σ b,i d σ,i , P i , Q i are positive definite matrices, σ b,i is estimated maximum value for σ i (t), λ max (•) and λ min (•) respectively denote the maximum and minimum eigenvalue of a matrix, and d σ,i is defined with (49).
Remark 1. Defining the following ideal system with no uncertainties: The difference between the output of ideal system (65) and the output of the reference system (58) satisfies Note that δy i (s) is a result of σ re f ,i (s) filtered by the low-pass system G ii (s) cascading with high-pass filter 1 − C i (s).So, if the cut-off frequency of the high-pass filter is much larger than the bandwidth of the low-pass system G ii (s), the cascaded system will result in a no-pass filter and we get δy i (s) → 0 .Thus, it is reasonable to make the approximation that y ≈ y re f ≈ y.

Fuel Starvation Protection
It is very important to point out that the protection constraint imposed in the MPC still cannot guarantee the avoidance of fuel starvation for the following two reasons: 1. Since the sampling time of MPC is usually much longer than the time-scale of power electronic devices.If the unit load demand changes during the interval of two MPC sampling point, MPC may not be able to respond in time before the fuel starvation occurs.
2. Since the action of the valves usually takes seconds, the supplementary fuel lags the instantaneous consumption caused by the leap of the current.
To overcome the first defect, protection logic, as shown in Figure 6, is added to the control system.When the supervisory level produces a new reference current, the current value will be hold in the memory until the next sampling time of MPC.Additionally, a load governor, which is designed as a first-order filter to smooth the current reference is transmitted to the converter system.

Simulation Results
The simulation is conducted on the model introduced in Section 3 with Matlab/Simulink TM with the configuration of a 5 kw SOFC stack.The SOFC stack is composed of 192 cascading basic fuel cells.Without a loss of generality, the simulation is conducted on a single cell.However, the experimental data is amplified by the cell number in the figures drawing for better exhibition.The power tracking performance of the system without SP scheduler is almost as good as that of the system with the proposed control strategy.The reason is that the PID-only control architecture directly controls the output power through changing the current.The architecture in fact avoids the nonlinearity associated with the SOFC.Nevertheless, when checking the waste power of the SOFCs in Figure 9, the system without SP scheduler turns out to be with the highest waste power.The SP + PID system does not have as good economic performance as the one with proposed control strategy since it actually does not track the set points of output voltage and current so well, despite the fact that both of them receives optimal SPs from the SP scheduler.This explanation can be verified by the curves of output voltage and current in Figure 10. Figure 11 presents the fuel utilization curve.Both of the three control methods can handle the fuel utilization in the safe range.The proposed control strategy has a more flexibility action of the manipulated inputs in order to have better tracking performance.As a consequence, the deviation of FU to the expected value is larger than the other two systems in some transient states.Nevertheless, the safe range of FU is always guaranteed because of the protective mechanism employed in this work.Figure 12 shows the manipulated inputs of the SOFCs.The SOFCs are with different steady-state temperature under different load.Since the thermo-dynamic process is so slow that by the end of the simulation, the temperature has actually not reached a new steady state.This causes a continuous regulating of the inlet air flow in order to compensate for the voltage perturbation brought by the temperature change.Based on the simulation results, the features of the proposed control strategy are concluded in Table 3 along with that of the two comparative simulation.By comparing the Simulation 1 (the proposed control strategy) and Simulation 3 (SP Scheduler + PID), the advantage of the proposed L1-MPC in tracking performance and constraint satisfaction is highlighted.It is worth pointing out that the L1 adaptive controller is the key technique that contributes to the tracking performance.It is well-known that the performance of MPC depends on the their tuned parameters.The most contributing parameters for a linear-model based MPC is the weight matrix Q and R, in (56).Trade-offs can be made between tracking performance and robustness by tuning Q and R. When the plant are with nonlinearities, the MPC have to be tuned in a robust way in order to guarantee the stability of the system when the operation point deviates from the nominal one.In this case, the control action is usually conservative, thus the performance is sacrificed.In our work, we add a L1 adaptive control term in the control signal (52), the nonlinearity of the system dynamic is compensated in a fast way.In this case, the MPC parameters can be tuned in more aggressive way, aka, being able to achieve better tracking performance.Similarly, The PID controller in Simulation 3 has to be tuned sufficiently conservative to guarantee with the existence of system nonlinearity, although it has the potential to have a better performance when working in a small range of working conditions.We highly recognize that there exist plenty of advanced PID that can deal with the system nonlinearities.However, it is out of the scope of this paper.We will not make further discussion here although it can be a very attractive work in the future.
The comparison between Simulation 3 and Simulation 2 (Only PID) further shows the effect of set point scheduler to achieve a better economic performance of the SOFCs.Indeed, in Simulation 3, we have implemented an entirely difficult control structure.However, this is the plainest way to make power command tracking with constant fuel utilization maintained.

Conclusions
This paper proposes a procedure to design a hierarchical control strategy for optimal power tracking operation of grid-connected SOFCs.The control strategy consists of a set point scheduler in the upper level and a SOFC controller working together with VSI devices in the lower level.The setpoints scheduler finds the optimal values of output voltage and current minimizing the internal waste power, based on a comprehensive steady-state model of the SOFC.The acquired optimal values are transmitted to the lower-level SOFC controller and VSI device as their set points.To overcome the insufficiency of MPC in dealing with the nonlinearities and uncertainties in the system dynamics, a combined L1-MPC control algorithm is put forward in the SOFC controller design.Additional protection logic is included in the control strategy to prevent SOFC from fuel starvation under unit load demand variation within the SOFC controller's the sampling interval.The simulation result shows the superiority of the proposed method.
A future direction of effort might be improving the economy in the system's transient states.The set point scheduler optimizes the system economy merely in steady states.However, when the micro-grid has frequent demand for peaking regulation, the unit load demand would change continually.In this case, the SOFC system rarely operates in steady states.Actually, this situation is very likely to happen in those micro-grids operating in islanding mode.As a consequence, the economy in transient process becomes crucial to the general economy of the SOFCs.It will be of great interest in the future to apply those approaches to take care of the plant's economy in the transient process.

Anode: 2H 2 +Figure 1 .
Figure 1.Schematic diagram of a solid oxide fuel cells (SOFC) power plant connected to an AC-grid.

Figure 2 .
Figure 2. The cutaway diagram of a tubular solid oxide fuel cell.

Figure 3 .
Figure 3. Heat Transfer in a Tubular solid oxide fuel cell.

Figure 4 .
Figure 4.The hierarchical control structure for optimal SOFC operation.

Figure 5 .
Figure 5.The control scheme designed for optimal SOFC operation.

Figure 6 .
Figure 6.The logical flow chart for fuel starvation protection.

Table 2 .
The components of signal flow in designed control scheme.

Table 3 .
The conducted simulation and expected performance.