A Fuzzy Adaptive PID Controller Design for Fuel Cell Power Plant

Yuxiao Qin 1, Li Sun 1,* ID , Qingsong Hua 2,* and Ping Liu 3 1 Key Lab of Energy Thermal Conversion and Control of Ministry of Education, Chien-Shiung Wu College, Southeast University, Nanjing 210096, China; 213151403@seu.edu.cn 2 School of Mechanical and Electrical Engineering, Qingdao University, Ningxia Road 308, Qingdao 266071, China 3 School of Electrical Engineering and Intelligentization, Dongguan University of Technology, Dongguan 523808, China; pingliu81@163.com * Correspondence: sunli12@seu.edu.cn (L.S.); qihu@qdu.edu.cn (Q.H.); Tel.: +86-180-0120-6210 (L.S.)


Introduction
A fuel cell is an efficient electrochemical device which translates chemical energy into electrical energy [1].It is widely used as a power supply for vehicles and large power plants [2,3].Because of its high efficiency and low emission, it is likely to replace conventional fossil energy [4].Solid oxide fuel cells (SOFCs) are the most popular type of fuel cells because they have high energy efficiency and long-term stability [5,6].Additionally, SOFCs are cost efficient because they do not require precious metals [6].
Although SOFCs have many notable advantages, many challenges still exist.To get a high-quality power supply from SOFCs, the output voltage must be accurately controlled.However, due to the nonlinearity, slow dynamics, and limited constraint of SOFCs [7], controlling the output voltage can be extremely hard.Furthermore, actuator saturation may occur due to limited fuel flow, and improper control may cause the fuel utilization to fall out of the operating threshold [8].As a result, an efficient control scheme for the output voltage of SOFCs has become a primary focus for fuel cell research.
Conventional proportional-integral-derivative (PID) controllers dominate in power plant control systems because of their reliability and simplicity [9,10].A new survey of over 100 boiler-turbine units in Guangdong Province, China shows that over 90% use PID controllers [10].However, studies have shown that conventional PID controllers are not suitable for the control of SOFCs [11], and a comparable H∞ control strategy [12] has also shown to be unsatisfying.To deal with the control difficulties of SOFCs, more advanced control strategies should be applied.With the characteristics of high robustness and satisfying control performance in nonlinear systems, model predictive control (MPC) strategies are drawing more and more attention from researchers.Various MPC strategies based on different models [13][14][15] have been developed for the control of SOFCs, and they all have satisfactory simulation results.However, to apply MPC methods, a huge amount of calculation is required, which will lead to the demand of high-performance computers.This adds to the difficulty of realizing this type of control scheme.
Considering the drawbacks of the MPC strategy, we choose to focus on fuzzy logic control.Fuzzy logic control (FLC) is a kind of control scheme based on the fuzzy logic rules created by experts in the field.Fuzzy logic control is robust, easy to tune, and can control multiple input and output sources [16].In addition, with a series of well-designed fuzzy logic rules, FLC can reach a good control performance with a small amount of calculation.For these reasons, fuzzy logic control schemes are commonly used for nonlinear systems such as fuel cells [16,17] and wind turbine systems [18].The proposed controller combines the FLC and PID controls, and the role of the fuzzy controller is to modify the proportion gain K p and the integral gain K i , at the same time.Due to the constraint of the fuel flow rate, an anti-windup structure based on the back-calculation method [19] is applied to attenuate the phenomenon of actuator saturation.In addition, the proposed controller applies a differential forward algorithm to reduce the "snap back" phenomenon of the derivative action and to improve the robustness of the controller at the same time [20].
To summarize, this paper (1) identifies the SOFC model and uses the identified model as a control object to optimize conventional PID controllers; (2) combines a fuzzy logic control scheme and a PID control scheme to design our proposed fuzzy adaptive PID controller; and (3) develops an anti-windup structure based on a back-calculation method to reduce saturation time and overshoot.This paper is organized as follows: In Section 2, the model of the SOFC system is demonstrated.Section 3 identifies the original model, and tunes the parameter of the PID controller.The designing of the fuzzy controller is described in Section 4. Section 5 compares the results of the dynamic response of normal PID controller and the proposed controller, and conclusions are drawn in Section 6.

Operation Principle of SOFC
Like other fuel cells, SOFCs convert chemical energy into electric energy.On the anode, there is a continuous fuel flow such as H 2 , CH 4 , and as for the cathode, there is a continuous flow of O 2 or air.The O 2 gets the electron and becomes O 2 , and creates a concentration gradient.The O 2− is delivered from the cathode to the anode, and then reacts with the fuel to transport the electron [21].In this way, a circuit loop is formed.The formulas of the electrochemical reactions are given as follows: Reaction on the anode: 2H

Model Description
The dynamic model of SOFC we explored in this paper is derived from Reference [22], which is also applied by other scientists in the study of SOFCs [23,24].This model takes both the electrochemical and the thermal aspects into account.Its structure is shown in Figure 1. Figure 2 shows the detailed structure of the thermal block.

Model Description
The dynamic model of SOFC we explored in this paper is derived from Reference [22], which is also applied by other scientists in the study of SOFCs [23,24].This model takes both the electrochemical and the thermal aspects into account.Its structure is shown in Figure 1. Figure 2 shows the detailed structure of the thermal block.The controlled variable of the system is the output voltage of the SOFC Vfc (V), the manipulate variable of the system is the fuel flow rate qf (mol/s), while the disturbance variable is the current of the external load Ifc (A).out T is the output temperature,  1. also applied by other scientists in the study of SOFCs [23,24].This model takes both the electrochemical and the thermal aspects into account.Its structure is shown in Figure 1. Figure 2 shows the detailed structure of the thermal block.The controlled variable of the system is the output voltage of the SOFC Vfc (V), the manipulate variable of the system is the fuel flow rate qf (mol/s), while the disturbance variable is the current of the external load Ifc (A).out T is the output temperature,  The controlled variable of the system is the output voltage of the SOFC V fc (V), the manipulate variable of the system is the fuel flow rate q f (mol/s), while the disturbance variable is the current of the external load I fc (A).T out is the output temperature, q in O 2 (mol/s) is the flow rate of the oxygen, q H 2 (mol/s) is the flow rate of the hydrogen, and p H 2 (Pa), p H 2 O (Pa), p O 2 (Pa) represent for the partial pressure of the hydrogen, steam, and oxygen, respectively.The other parameters of the model are listed in Table 1.The partial pressures can be estimated by the electrochemical relationships as follows [25]: V 0 (V) is the voltage of the electricity generated by the transport of the electron in the cell, it is calculated with the Nerst's equation: Using the activation loss, concentration loss, and the ohmic loss, the output voltage can be calculated as: where V fc (V) stands for the output voltage of SOFC, V act (V), V ohm (V), V conc (V) is the activation voltage loss, the ohmic voltage loss and the concentration voltage loss, respectively.They can be represented as a function of temperature by the following equations: where T ( • C) denotes the temperature of the fuel cell.The Fourier number F o and the source term number S 0 are the two major parameters to demonstrate the thermal equations in a dimensionless form [23].They are expressed as: where t (s) denotes the relaxation time which represents the time for the system to attain 90% of the new steady-state value.∆T( • C) denotes the temperature change during the period of relaxation time.
The relationship between F o and S 0 can be defined as [23]: Based on Equations ( 9)-( 11), the expression of T out can be modified as follow: where T ( • C) is the current temperature of the fuel cell, T int ( • C) is the initial temperature of the fuel cell at the time when no load is added, and dt is the step time of Simulink.Fuel utilization factor u f is also one of the important indicators of SOFC.It reveals the proportion of the acted fuel flow to the input fuel flow: q in H 2 in the input of the fuel flow, and q out H 2 is the output of the fuel flow.Usually, u f is constrained to the range of 0.7-0.9 for high efficiency and safety.

Problem Description
The aim of the control design is to maintain the output voltage to keep up with the rated voltage as accurately as possible, and decrease the voltage variation when the current loads change.The top three problems in control design are listed below:

•
Nonlinearity of the System: owning to Nerst's equation in (4), the SOFC system is characterized as nonlinearity [26], and this deteriorates the control performance of the controller when the working conditions drift off from the ideal working conditions.

•
Fuel Flow is Restricted: fuel flow must be restricted between 0 and 2 mol/s, which may cause actuator saturation and the dynamic property may be deteriorated.

•
Hysteresis of the System: the rated voltage and the current load change rapidly, while the effect of the fuel flow on the output voltage is comparatively slow.

System Identification and Controller Tuning
Nowadays, conventional PID controllers are still wildly utilized in the control of industrial process for their high reliability and simplicity [10].The block diagram of a conventional PID controller is shown in Figure 3.The control equation is given as follows: where u(t) stands for the control action, K p , K i , and K d are the proportional integral and derivative gains.e(t) is the tracking error signal.
where ( ) stands for the control action, Kp, Ki, and Kd are the proportional integral and derivative gains.( ) is the tracking error signal.In this section, we aim to linearize the SOFC model and use the identified model for the purpose of controller tuning.We will use a conventional PID controller to control the SOFC system, and find the parameter of the PID controller with the best step response curve.

System Identification
Given that the SOFC model is a nonlinear system, it's hard to tune the controller with the original model.To get the proper parameter of the PID controller, the first step is to linearize the model.We use the System Identification Toolbox of Matlab to identify the system, and we get the identified model as follows: Figure 4 shows the step response of the original model and the identified model.From the response curve, we can tell that the identified model behaves remarkably similar to the original model.Thus, the identification is successful, and we can use the identified model as the linearized model for the controller tuning.In this section, we aim to linearize the SOFC model and use the identified model for the purpose of controller tuning.We will use a conventional PID controller to control the SOFC system, and find the parameter of the PID controller with the best step response curve.

System Identification
Given that the SOFC model is a nonlinear system, it's hard to tune the controller with the original model.To get the proper parameter of the PID controller, the first step is to linearize the model.We use the System Identification Toolbox of Matlab to identify the system, and we get the identified model as follows: G(s) = 2.417s + 0.8276 s 3 + 0.716s 2 + 0.1808s + 0.006914 (15) Figure 4 shows the step response of the original model and the identified model.From the response curve, we can tell that the identified model behaves remarkably similar to the original model.Thus, the identification is successful, and we can use the identified model as the linearized model for the controller tuning.
( ) ( ) ( ) (14) where ( ) stands for the control action, Kp, Ki, and Kd are the proportional integral and derivative gains.( ) is the tracking error signal.In this section, we aim to linearize the SOFC model and use the identified model for the purpose of controller tuning.We will use a conventional PID controller to control the SOFC system, and find the parameter of the PID controller with the best step response curve.

System Identification
Given that the SOFC model is a nonlinear system, it's hard to tune the controller with the original model.To get the proper parameter of the PID controller, the first step is to linearize the model.We use the System Identification Toolbox of Matlab to identify the system, and we get the identified model as follows: Figure 4 shows the step response of the original model and the identified model.From the response curve, we can tell that the identified model behaves remarkably similar to the original model.Thus, the identification is successful, and we can use the identified model as the linearized model for the controller tuning.

Controller Tuning
We use the pidTuner Toolbox of Matlab to tune the PID controller and get several groups of proportional, integral, and derivative gains.Then we apply these parameters to the PID controller and record the results of the step response of each group of parameters.The simulation results are shown in Figure 5.It is easy to tell that the overshoot of the output is hard to eliminate with the control of only a conventional PID controller.If the overshoot is slashed excessively, the settling time will be improperly long.
We use the pidTuner Toolbox of Matlab to tune the PID controller and get several groups of proportional, integral, and derivative gains.Then we apply these parameters to the PID controller and record the results of the step response of each group of parameters.The simulation results are shown in Figure 5.It is easy to tell that the overshoot of the output is hard to eliminate with the control of only a conventional PID controller.If the overshoot is slashed excessively, the settling time will be improperly long.

Fuzzy Control Design
The conventional PID controller is ideal for the control of the linear system.However, due to the nonlinearity of the SOFC system, the control performance will deteriorate with the variation of the working conditions [18].To this end, we can use PID controllers with variable parameters to overcome the change of the operating point.Based on the variable parameters, the proposed fuzzy adaptive PID controller was designed.The basic block diagram of the system is shown in Figure 6, and the detailed structure of the controller is shown in Figure 7.
The equation of the fuzzy adaptive PID is designed as follows: The detailed data of the overshoot and settling time of the step responses are shown in Table 2.After synthesizing all of the results, we finally set the parameter of the conventional PID controller as follows to reduce both the overshoot and the settling time:

Fuzzy Control Design
The conventional PID controller is ideal for the control of the linear system.However, due to the nonlinearity of the SOFC system, the control performance will deteriorate with the variation of the working conditions [18].To this end, we can use PID controllers with variable parameters to overcome the change of the operating point.Based on the variable parameters, the proposed fuzzy adaptive PID controller was designed.The basic block diagram of the system is shown in Figure 6, and the detailed structure of the controller is shown in Figure 7.
The equation of the fuzzy adaptive PID is designed as follows: where u(t) stands for the control action, K p , K i , and K d are the proportional integral and derivative gains.e(t) is the tracking error signal. ( ] ( ) ( ) ( ) where ( )

Fuzzy PI Controller
The input of the controller are error e and the change of error, they are given as follows: ( ) ( ) ( ) ( ) ( ) ( ) where ( ) k r denotes the rated voltage, and ( ) denotes the output voltage.The outputs are the modifying value of the proportion gain Kp1 and the modifying value of the integral gain Ki1.The fuzzy controller includes three stages: fuzzification, fuzzy logic judgement, and defuzzification.

Fuzzy PI Controller
The input of the controller are error e and the change of error, they are given as follows: ( ) ( ) ( ) ( ) ( ) ( ) where ( ) k r denotes the rated voltage, and ( ) denotes the output voltage.The outputs are the modifying value of the proportion gain Kp1 and the modifying value of the integral gain Ki1.The fuzzy controller includes three stages: fuzzification, fuzzy logic judgement, and defuzzification.

Fuzzy PI Controller
The input of the controller are error e and the change of error, they are given as follows: where r(k) denotes the rated voltage, and y(k) denotes the output voltage.The outputs are the modifying value of the proportion gain K p1 and the modifying value of the integral gain K i1 .The fuzzy controller includes three stages: fuzzification, fuzzy logic judgement, and defuzzification.

Fuzzification
Fuzzification is the process of transforming crisp sets into linguistic fuzzy sets by using fuzzy membership functions.To design the fuzzy logic rule, seven fuzzy items are defined: NB, NM, NS, ZO, PS, PM, and PB, which stands for negative big, negative medium, negative small, zero, positive small, positive medium, and positive big, respectively.The fuzzy membership function of NB and PB are Pi type, while the fuzzy membership function for the rest of the fuzzy items are triangular type.The universe of discourse of e, ec, K p1 and K i1 are all [-3, 3], and the gain for e, ec, K p1 and K i1 are 0.5, 0.5, 0.1 and 0.04, respectively.The curves are displayed in Figures 8 and 9.
membership functions.To design the fuzzy logic rule, seven fuzzy items are defined: NB, NM, NS, ZO, PS, PM, and PB, which stands for negative big, negative medium, negative small, zero, positive small, positive medium, and positive big, respectively.The fuzzy membership function of NB and PB are Pi type, while the fuzzy membership function for the rest of the fuzzy items are triangular type.The universe of discourse of e, ec, Kp1 and Ki1 are all [-3, 3], and the gain for e, ec, Kp1 and Ki1 are 0.5, 0.5, 0.1 and 0.04, respectively.The curves are displayed in Figures 8 and 9.

Fuzzy Logic Judgment
The fuzzy logic judgment is the core of fuzzy control which is based on fuzzy logic rules.To design the controller, Mamdani Inference Method [27] was applied for the fuzzy judgment.The output membership value is calculated below: The fuzzy rules are established with 49 items, which are listed in Table 3.The rules can be expressed as " If e is Ai and ec is Bi, then Kp1 is Ci and Ki1 is Di", where Ai, Bi, Ci and Di are the fuzzy items of the universe of discourse of e, ec, Kp1 and Ki1, respectively.Fuzzification is the process of transforming crisp sets into linguistic fuzzy sets by using fuzzy membership functions.To design the fuzzy logic rule, seven fuzzy items are defined: NB, NM, NS, ZO, PS, PM, and PB, which stands for negative big, negative medium, negative small, zero, positive small, positive medium, and positive big, respectively.The fuzzy membership function of NB and PB are Pi type, while the fuzzy membership function for the rest of the fuzzy items are triangular type.The universe of discourse of e, ec, Kp1 and Ki1 are all [-3, 3], and the gain for e, ec, Kp1 and Ki1 are 0.5, 0.5, 0.1 and 0.04, respectively.The curves are displayed in Figures 8 and 9.

Fuzzy Logic Judgment
The fuzzy logic judgment is the core of fuzzy control which is based on fuzzy logic rules.To design the controller, Mamdani Inference Method [27] was applied for the fuzzy judgment.The output membership value is calculated below: where cp, ci are the elements in the fuzzy set of Kp1 and Ki1, respectively, and the membership degrees of cp and ci , respectively.The fuzzy rules are established with 49 items, which are listed in Table 3.The rules can be expressed as " If e is Ai and ec is Bi, then Kp1 is Ci and Ki1 is Di", where Ai, Bi, Ci and Di are the fuzzy items of the universe of discourse of e, ec, Kp1 and Ki1, respectively.

Fuzzy Logic Judgment
The fuzzy logic judgment is the core of fuzzy control which is based on fuzzy logic rules.To design the controller, Mamdani Inference Method [27] was applied for the fuzzy judgment.The output membership value is calculated below: where c p , c i are the elements in the fuzzy set of K p1 and K i1 , respectively, and µ Kp1 c p , µ Ki1 (c i ) are the membership degrees of c p and c i , respectively.The fuzzy rules are established with 49 items, which are listed in Table 3.The rules can be expressed as " If e is A i and ec is B i , then K p1 is C i and K i1 is D i ", where A i , B i , C i and D i are the fuzzy items of the universe of discourse of e, ec, K p1 and K i1 , respectively.In the process of defuzzification, the weighted mean was utilized to convert the vector into a single value [28], and then the output of the fuzzy controller was determined.
The calculation of K p1 and K i1 are below: With the output of the Fuzzy controller, the proportion gain K p and the integral gain K i can be calculated as follows: where K p0 is the initial value of the proportion gain, K i0 is the initial value of the integral gain, K 1 and K 2 are both constant which denotes the gain of K p1 and K i1 , respectively.

The Realization of Anti-Windup
In this paper, the fuel flow rate is limited to between 0 and 2 mol/s.Due to this limitation, there might be the problem of actuator saturation.To this end, a method based on back-calculation is applied [19].As shown in Figure 7, a feedback for the integrator is designed.The input of the integrator can be calculated as: where e 1 (t) denotes the input of the integrator, e 0 (t) denotes the tracking error signal, u 0 (t), and u 1 (t) denotes the controlled signal before and after the saturation block, respectively.u 2 (t) is the difference between u 1 (t) and u 0 (t) as (17) shows.From Equations ( 17) and ( 26), we can tell that when the controlled signal is above the upper limit, u 1 (t) will remain the upper limit, and u 1 (t) will be smaller than u 0 (t).This will set u 2 (t) as a negative signal to reduce the integral action.On the contrary, when the controlled signal is under the lower limit, u 1 (t) will remain the lower limit, and u 1 (t) will be bigger than u 0 (t).This will set u 2 (t) as a positive signal to increase the integral action.In this way, the controller is able to get rid of the situation of actuator saturation much faster and reduce the overshoot at the same time.

Differential Forward Algorithm
Owning to the fact that the differentiator is very sensitive to the variation of the input [20], the proposed controller applies a differential forward algorithm to overcome the sensibility of the differentiator.In the conventional PID, the controller differentiates to the tracking error signal e(t), which is defined in (14).As we can see, when the rated voltage r(t) changes, the output voltage y(t) still remains the same, which causes a sudden variation of e(t).This will lead to a "snap back" of the differential action.
On the contrary, the differential forward algorithm differentiates to the output voltage y(t) directly.Meanwhile, the differentiator applies the first order actual differential algorithm, which can be expressed as: In this way, the variation of the rated voltage is not able to affect the differentiator, and because the output voltage has a relatively slow rate of change, it will not produce a large control value.With this method, the problem of the "snap back" of the differential action is fixed.

Simulation
In this section, to show the comparison of the robustness of the proposed controller and the conventional PID controller, simulation of step response and disturbance response based on the tuned PID controller in Section 3.2 and the proposed controller in Section 4 is implemented, and the results of the simulation are displayed.

Simulation of Step Response
During the process of step response, the current is kept at 300 A, the rated voltage begins with 305 V, and then it is set to increase by 15 V at 50 s, 150 s, 250 s, and 350 s.Then the rated voltage decreases by 30 V at 450 s and 550 s to return to 305 V.The results of the simulation are shown in Figure 10.
directly.Meanwhile, the differentiator applies the first order actual differential algorithm, which can be expressed as: ( ) In this way, the variation of the rated voltage is not able to affect the differentiator, and because the output voltage has a relatively slow rate of change, it will not produce a large control value.With this method, the problem of the "snap back" of the differential action is fixed.

Simulation
In this section, to show the comparison of the robustness of the proposed controller and the conventional PID controller, simulation of step response and disturbance response based on the tuned PID controller in Section 3.2 and the proposed controller in Section 4 is implemented, and the results of the simulation are displayed.

Simulation of Step Response
During the process of step response, the current is kept at 300 A, the rated voltage begins with 305 V, and then it is set to increase by 15 V at 50 s, 150 s, 250 s, and 350 s.Then the rated voltage decreases by 30 V at 450 s and 550 s to return to 305 V.The results of the simulation are shown in Figure 10.The detailed data of the simulation result are displayed in Table 4.The detailed data of the simulation result are displayed in Table 4.As we can tell from the results, the control performance of the PID controller deteriorates with the variation of the working conditions, which is illustrated by the increase of the overshoot and the settling time.As for the proposed controller, the overshoot keeps at 0 for all the working conditions, and the settling time is also much shorter than the conventional PID controller.Meanwhile, with the design of anti-windup structure, there is also a reduction in the saturation time.Especially in working condition 4, the saturation time of the conventional PID controller is two times bigger than the proposed controller.

Simulation of Disturbance Response
During the process of disturbance response, the rated voltage is kept at 305 V, and the current starts with 300 A, then the current is set to increase by 75 A at 50 s, 100 s, 150 s, and 200 s and then decreases by 150 A at 250 s and 300 s to return to 300 A. The results of the simulation are shown in Figure 11.
The detailed data of the simulation result are displayed in Table 5.Similar to the result of step response, the proposed controller shows a better control performance in the control of disturbance response.To illustrate, with the variation of the working condition, the absolute value of the overshoot of the conventional PID controller varies from 0.33 to 1.11, and the overshoot is kept at 0 for the proposed controller.The settling time of the proposed controller is much shorter than the conventional PID controller.In addition, the average of the saturation time of the proposed controller is also shorter than the conventional PID controller.In working condition 4, especially, the saturation time of the proposed controller is three quarters of the saturation time of conventional PID controller.absolute value of the overshoot of the conventional PID controller varies from 0.33 to 1.11, and the overshoot is kept at 0 for the proposed controller.The settling time of the proposed controller is much shorter than the conventional PID controller.In addition, the average of the saturation time of the proposed controller is also shorter than the conventional PID controller.In working condition 4, especially, the saturation time of the proposed controller is three quarters of the saturation time of conventional PID controller.

Discussion
From the simulation results, we can tell that the proposed controller has the advantage of no

Discussion
From the simulation results, we can tell that the proposed controller has the advantage of no overshoot, less settling time, and less saturation time.Under the shifting working conditions, the working performance of the PID controller deteriorates significantly, while the proposed controller is less affected and performs much better than the conventional PID controller.
In conclusion, the proposed controller is able to handle the nonlinearity of SOFC, and it is expected to be an efficient solution for the problem of controlling the output voltage of SOFCs.

Figure 1 .
Figure 1.The block diagram of the dynamic model of the solid oxide fuel cell (SOFC).

Figure 2 .
Figure 2. The block diagram of the thermal block of the SOFC.

in O q 2 2 Hq 2 H
(mol/s) is the flow rate of the oxygen, (mol/s) is the flow rate of the hydrogen, and ) represent for the partial pressure of the hydrogen, steam, and oxygen, respectively.The other parameters of the model are listed in Table

Figure 1 .
Figure 1.The block diagram of the dynamic model of the solid oxide fuel cell (SOFC).

Figure 1 .
Figure 1.The block diagram of the dynamic model of the solid oxide fuel cell (SOFC).

Figure 2 .
Figure 2. The block diagram of the thermal block of the SOFC.

2 2 Hq 2 H
(mol/s) is the flow rate of the oxygen, (mol/s) is the flow rate of the hydrogen, and ) represent for the partial pressure of the hydrogen, steam, and oxygen, respectively.The other parameters of the model are listed in Table1.

Figure 2 .
Figure 2. The block diagram of the thermal block of the SOFC.

Figure 4 .
Figure 4. Comparison of the step response of the original model and the identified model.

Figure 4 .Figure 4 .
Figure 4. Comparison of the step response of the original model and the identified model.

Figure 5 .
Figure 5. Step response of the identified model.

Figure 5 .
Figure 5. Step response of the identified model.

Figure 6 .
Figure 6.Block diagram of the fuzzy adaptive PID controller.

Figure 7 .
Figure 7. Detailed structure of the proposed controller.

Figure 6 .
Figure 6.Block diagram of the fuzzy adaptive PID controller.

Figure 6 .
Figure 6.Block diagram of the fuzzy adaptive PID controller.

Figure 7 .
Figure 7. Detailed structure of the proposed controller.

Figure 7 .
Figure 7. Detailed structure of the proposed controller.

Figure 8 .
Figure 8. Membership function of the input of fuzzy PID: (a) input with e; (b) input with ec.

Figure 9 .
Figure 9. Membership function of the output of fuzzy PID: (a) output with Kp1; (b) output with Ki1.
cp, ci are the elements in the fuzzy set of Kp1 and Ki1, respectively, and of cp and ci , respectively.

Figure 8 .
Figure 8. Membership function of the input of fuzzy PID: (a) input with e; (b) input with ec.

Figure 8 .
Figure 8. Membership function of the input of fuzzy PID: (a) input with e; (b) input with ec.

Figure 9 .
Figure 9. Membership function of the output of fuzzy PID: (a) output with Kp1; (b) output with Ki1.

Figure 9 .
Figure 9. Membership function of the output of fuzzy PID: (a) output with K p1 ; (b) output with K i1 .

Figure 10 .
Figure 10.Simulation result of step response.

Figure 11 .
Figure 11.Simulation result of disturbance response.

Figure 11 .
Figure 11.Simulation result of disturbance response.

Table 1 .
Nominal parameters of the SOFC system.

Table 2 .
Result of the parameter tuning simulation.

Table 2 .
Result of the parameter tuning simulation.

Table 3 .
The fuzzy logic rules of the controller.

Table 4 .
Detailed data of the simulation result of step response.

Table 4 .
Detailed data of the simulation result of step response.

Table 5 .
Detailed data of the simulation result of disturbance response.

Table 5 .
Detailed data of the simulation result of disturbance response.