Fault Detection and Isolation for a Cooling System of Fuel Cell via Model-based Analysis

: The development of fuel cell electric vehicles in recent years has led to increased interest in the use of fuel cells as sources of renewable energy. To achieve successful commercialization of fuel cell vehicles, it will be necessary to guarantee the safety, reliability, and lifetime of fuel cell systems by predictive fault detection and isolation (FDI). In this study, the parity equation, an observer, and a Kalman ﬁlter are employed together to compare the characteristics of FDI, focusing on the sensors of the thermal management system. Residuals corresponding to the di ﬀ erence between temperature outputs of linear models under driving cycles and nonlinear temperature outputs are used to isolate faults. Then, assessment of three model-based sensor FDI schemes is used to isolate sensor faults using the Cumulative Sum Control Chart (CUSUM) method. Generated residuals are evaluated by CUSUM to detect the presence of a sensor fault. As a result, isolated sensor faults are assessed. residuals parity Kalman ﬁlter, a lumped fuel cell system.


Introduction
The development of fuel cell electric vehicles (FCEVs) in recent years has led to increased interest in the use of fuel cells as sources of renewable energy for fuel cell vehicles, owing to their high efficiency, low pollution, and good power density [1][2][3][4][5][6]. Proton exchange membrane fuel cell (PEMFC) systems are widely applied to fuel cell vehicles because of their quick fuelling and extended running time [7,8].
To achieve commercialization, however, it will be necessary to solve technical issues related to the safety, reliability, and lifetime of fuel cell systems [9,10]. For commercialization, automakers around the world are targeting 5000 h of operation for transportation applications by 2020. However, according to a research survey, 78% of fuel cell vehicles did not exceed 2000 h of durability. Especially, since the reliability and durability affect human safety during vehicle operation, many researchers have performed durability tests and faults diagnosis. In addition, fault detection and fault diagnosis have become more important at the system level, in particular because fuel cell systems such as air supply, fuel supply, thermal management, power converter, etc., are dependently operated, and thus a fault in any subsystem can damage the entire fuel cell vehicle system. Therefore, it is necessary to detect initial faults early, before a minor fault becomes a serious fault.
A few researchers have attempted to address the FDI of fuel cell systems using model-based methods. Pei et al. investigated water fault diagnosis-based pressure drop using different indicators. They reviewed recent studies on the use of pressure drop, as found in calculation results, to diagnose water fault [11]. Shao et al. presented an ANN (artificial neural network) ensemble method to diagnose The design of the nonlinear dynamic fuel cell models assumes: 1. The fuel on the anode side contains hydrogen and water vapour and the air on the cathode side contains oxygen, nitrogen, and water vapour; 2. In the humidifier model, the injected water is vaporized without heat losses; 3. Both the fuel and air were assumed to be ideal gases; 4. The temperature in the fuel cell stack is uniform; and 5. The electrochemical reaction is quasi-steady.

Air Supply Model
Air is supplied to the stack with a blower. A semiempirical model is applied to build a model of a turbo blower. The experimental data were obtained from Franklin et al. [27], and the performance curve and efficiency of the blower are shown in Figure 2a,b. The experimental data for the blower were provided in the form of a look-up table. The blower is operated by a DC motor. The model of the dynamic electric motor was derived from the specification data [24]. In addition, the PI controller regulates the motor speed to maintain an oxygen excess ratio of 2.2.  The design of the nonlinear dynamic fuel cell models assumes: 1.
The fuel on the anode side contains hydrogen and water vapour and the air on the cathode side contains oxygen, nitrogen, and water vapour; 2.
In the humidifier model, the injected water is vaporized without heat losses; 3.
Both the fuel and air were assumed to be ideal gases; 4.
The temperature in the fuel cell stack is uniform; and 5.
The electrochemical reaction is quasi-steady.

Air Supply Model
Air is supplied to the stack with a blower. A semiempirical model is applied to build a model of a turbo blower. The experimental data were obtained from Franklin et al. [27], and the performance curve and efficiency of the blower are shown in Figure 2a,b. The experimental data for the blower were provided in the form of a look-up table. The blower is operated by a DC motor. The model of the dynamic electric motor was derived from the specification data [24]. In addition, the PI controller regulates the motor speed to maintain an oxygen excess ratio of 2.2.
Air is supplied to the stack with a blower. A semiempirical model is applied to build a model of a turbo blower. The experimental data were obtained from Franklin et al. [27], and the performance curve and efficiency of the blower are shown in Figure 2a,b. The experimental data for the blower were provided in the form of a look-up table. The blower is operated by a DC motor. The model of the dynamic electric motor was derived from the specification data [24]. In addition, the PI controller regulates the motor speed to maintain an oxygen excess ratio of 2.2.

Hydrogen Supply Model
The anode fuel flow includes the hydrogen flow rate and the humid vapour. The inlet anode mass flow rate into the stack is calculated as where J is the current density, sto H2 is the hydrogen stoichiometry ratio, A is the active area through the cell, and n is the number of cells. The relative humidity φ in,an at the anode inlet is calculated as where y H 2 O,in,an is the mass fraction of vapor at the anode side, M H 2 O is the molecular weight of water, and p sat is the saturation pressure.

Dynamic Model
A dynamic model of the fuel cell stack is designed to implement the sensor fault scenario. The stack model assumes a lumped mass for the dynamic response of thermal energy input. The lumped stack obtains heat from the electrochemical reaction ( . H rea ) and the heat is rejected via convection heat transfer to the surroundings ( . Q surr ), involving gases exhausted from the anode channel and the cathode channel ( . Q gas ), and convection in the water coolant ( . Q cool ). The stack model is as follows [28]: where m sta , and c p , sta are the mass and the specific heat of fuel cell stack, and .
T sta is the stack temperature. A dynamic cooling reservoir model is also developed. The water in the fuel cell stack and cooling circuit is considered as a lumped mass. The reservoir is assumed to be placed at the inlet of the fuel cell stack. The governing equation of the coolant reservoir, then, is m c,Rv c p,Rv where m c,Rv , c p , RV is the mass of reservoir, and specific heat of reservoir, T RV is the reservoir temperature. T mix is the sum of coolant outlet temperature through the stack and the radiator. The total mass of water present in the reservoir, radiator, and cooling channel of the stack is 7 kg. A bypass valve and radiator fan are used to control the coolant inlet temperature. The dynamic thermal management system with fault diagnosis is shown in Figure 3.

Electrochemical Reaction Model
The fuel cell stack produces electric power and heat by chemical reaction. When the fuel cell stack is requested to supply power, electric power is determined by the current and the electric potential. The electric potential of the fuel cell stack is calculated by the following equation: where Vsta is stack voltage, R is electric resistance,  is water content, and Vact is activation loss; the Nernst voltage, Vnern, is described as follows [ In Equation (5), the activation loss Vact can be expressed as follows [26,29].
where  is a tuning factor that can be adjusted by experimental data (,  and is the oxygen concentration at the cathode. The Ohmic loss is expressed by the ionic resistance in the membrane. This loss term is applied in Equation (5) in terms of the water content (). The electric resistance in the membrane can be expressed as follows [29]:  (8) where z is membrane thickness.

State Variables for Design of FDI
For sensor FDI, the state variables in our study are the temperature of the fuel cell stack and the coolant inlet temperature. The selected state variables can be represented in nonlinear state space form.

Electrochemical Reaction Model
The fuel cell stack produces electric power and heat by chemical reaction. When the fuel cell stack is requested to supply power, electric power is determined by the current and the electric potential. The electric potential of the fuel cell stack is calculated by the following equation: where V sta is stack voltage, R is electric resistance, λ is water content, and V act is activation loss; the Nernst voltage, V nern , is described as follows [25]: In Equation (5), the activation loss V act can be expressed as follows [26,29].
where ε is a tuning factor that can be adjusted by experimental data (ε 1 = −0.5708, ε 2 = 0.0017, ε 3 = 0.000017903, ε 4 = −0.0001) and c O 2 is the oxygen concentration at the cathode. The Ohmic loss is expressed by the ionic resistance in the membrane. This loss term is applied in Equation (5) in terms of the water content (λ). The electric resistance in the membrane can be expressed as follows [29]: where z is membrane thickness.

State Variables for Design of FDI
For sensor FDI, the state variables in our study are the temperature of the fuel cell stack and the coolant inlet temperature. The selected state variables can be represented in nonlinear state space form.
The coolant flow rate and the opening area of the bypass valve are employed as the manipulated variables: The output variables are the same as the state variables: Disturbances are the result of fluctuation of current and ambient temperature: Thus, the state and measurement equations in the nonlinear system are rewritten using Equation (9): where f ( x(t), u(t)) is the nonlinear system function and g( x(t), u(t)) is the nonlinear measurement function; w(t) and v(t) are the Gaussian white noise and the covariance of the measurement noise, respectively.

Scheme of Model-Based FDI
The fault detection and isolation algorithms of the temperature sensors applied in our study are a parity space, a state observer, and a Kalman filter, shown in Figure 4. Sensor faults in a fuel cell system are diagnosed by detection algorithms, in parallel with the nonlinear fuel cell system model, based on inputs such as the current and the bypass valve fraction. The residual is a clue allowing us to judge sensor faults; residuals are generated by comparing the outputs of the nonlinear fuel cell system to the outputs of the parity space, state observer, and Kalman filter. When a fault is detected, the residual signal shows trends different from the regular pattern. Since the trends vary according to the algorithms, it is necessary to evaluate sensor faults by established criteria. In Figure 4, while the parity space algorithm produces residuals independently from the inputs, the state observer and Karman Filter deliver their output to the healthy model. Those outputs are used to evaluate signal faults.

Fault Definition Model
The stack temperature and coolant inlet temperature are key control factors for the thermal management system (TMS) of a fuel cell. Thus, this study considers sensor faults regarding the stack temperature and coolant inlet temperature. In Equation (14), fs is the sensor fault in the stack temperature sensor, and the coolant inlet sensor of the fuel cell.
In this study, two typical sensor faults are presented. Fault types are stuck and scaling, due to sensor core degradation and external shock or vibrations [30]. A stuck fault has a constant value caused by flaws resulting from sensor degradation during normal operation. Scaling is assumed to be in a form that does not change. Thus, the stuck and scaling faults can be considered additive faults, and can be expressed by the following equation:

Fault Definition Model
The stack temperature and coolant inlet temperature are key control factors for the thermal management system (TMS) of a fuel cell. Thus, this study considers sensor faults regarding the stack temperature and coolant inlet temperature. In Equation (14), f s is the sensor fault in the stack temperature sensor, and the coolant inlet sensor of the fuel cell.
In this study, two typical sensor faults are presented. Fault types are stuck and scaling, due to sensor core degradation and external shock or vibrations [30]. A stuck fault has a constant value caused by flaws resulting from sensor degradation during normal operation. Scaling is assumed to be in a form that does not change. Thus, the stuck and scaling faults can be considered additive faults, and can be expressed by the following equation: where x is a state variable that represents the stack temperature or the coolant inlet temperature, y is the corresponding sensor measurement, and f s is the sensor fault variable.

Design of Parity Space Algorithm
The parity space approach is an open-loop method, one that utilizes analytical redundancy. It is based on the discrete state space model expressed by Equation (14). This equation is established to link the sensor outputs and inputs, as follows [31]: where T, H u , H w , H v , and H f are the related matrices, and, Y(k) is the temperature sensor output.
In this study, a first-order discrete linear dynamic model is used, consisting of two inputs, two disturbances, and two outputs. The state variables are directly output. The detailed matrices T, H u , U, Y, and F are as follows:

Design of State Observer Algorithm
A linearized model was derived to implement the sensor FDI for thermal management of the fuel cell. The linear fuel cell system was developed in our previous study [31]. The nonlinear model was linearized by a Gaussian Jacobian method along the reference operating point. The vector .
x(t), which depends on the states x 1 and x 2 , is approximated by After the nonlinear state space form is redesigned into a continuous state space form, the state space representation of dynamic temperature that includes the sensor fault term is rewritten as follows: .
where F is assumed to be a known matrix, and f s is the magnitude of the sensor faults. The system matrix form is defined as follows: The state observer approach uses the output error between the estimated output and the model output. Hence, the linear fuel cell model was used to design the state observer. Among other approaches, the Luenberger state observer was designed for FDI. The linear fuel cell system, which can be described by the state space model in Equation (18), was used. Then, the state observer was designed to reconstruct the estimation state variable based on the measured inputs and outputs, as in the following equation:ˆ. x Using the observer gain, Equation (20) yields the implementation form of the state observer. It is assumed that the linear fuel cell system is observable.

Design of Kalman Filter Algorithm
A Kalman filter algorithm for the fuel cell system requires the discretization of the continuous state. Hence, the input u is assumed to be piecewise constant over the sampling period.
Thus, the input value is changed at discrete-time intervals kT, where upon Equation (13) becomes With the system matrix A d = e AT , B d = ( T o e Aτ dτ)B, C d = C, and D d = D. While the Kalman filter algorithm is used to generate the residual, the noise covariance is adjusted by the residual-based temperature. The Kalman filter algorithm requires system dynamic modes, initialization, and a measurement model based on Equation (14). It is assumed that an initial estimate of the state variable is known at time t k . Then, the corresponding error covariance matrices can be calculated as: The gain of the steady-state Kalman filter is as follows: The updated state estimate measurement is then

Residual Generation between Nonlinear System and Fault Free Model
The residual was generated using the difference between the output of the nonlinear and linear fuel cell models and the parity space algorithm, state observer algorithm, and Kalman filter algorithm.
where r i is the residual of the temperature i, and y noi andŷ loi represent the output signals measured in the nonlinear model and the estimated output temperature, which are defined by Equations (18), (21), and (26).

Fault of Sensor Evaluation with Residual and Isolation
Evaluation of the residual of the fault diagnosis scheme is important, and helps to identify the presence of a fault. A nonzero residual can be generated because of model error or sensor noise even when the system is in a normal state. This may induce missed or false fault detections. The number of false alarm signals was minimized by applying the Cumulative sum (CUSUM). Therefore, the log-likelihood ratio of the residual is defined as: where r(k) is the residual, L(.) is the log-likelihood ratio, and u and σ are the average value and standard deviation, respectively.
Then, CUSUM can be expressed as: Eventually, fault detection is defined as: where H is the residual threshold calibrated through the mean value. A result of f = 1 indicates a fault.

Validation of Stack Model by Polarization Curve
The I-V characteristic curve of the fuel cell is generally used to verify the fuel cell system model. In this study, the test bench for the unit cell was constructed to conduct model verification to confirm the accuracy of the model in terms of the static characteristics. Figure 5a,b shows the test bench for the unit cell and the validation results. The test facilities were composed of a unit cell, an MFC to measure the flow rates, a humidifier to humidify the gases, an air fan to reject the heat in the unit cell, and other sensors. The active area of the fuel cell stack is 25 cm 2 , the stack temperature is 70 • C, and the pressure of anode and cathode is 3 atm. To verify the designed fuel cell model, a turning coefficient was selected. Details of the tuning coefficients can be found in our previous research [25]. As shown in Figure 5b, I-V performance curves between the experimental data and the model over time are in good agreement. measure the flow rates, a humidifier to humidify the gases, an air fan to reject the heat in the unit cell, and other sensors. The active area of the fuel cell stack is 25 cm 2 , the stack temperature is 70 °C, and the pressure of anode and cathode is 3 atm. To verify the designed fuel cell model, a turning coefficient was selected. Details of the tuning coefficients can be found in our previous research [25]. As shown in Figure 5b, I-V performance curves between the experimental data and the model over time are in good agreement.

Fault Scenarios
The ability to detect fault of the parity, the observer, and the Kalman filter models proposed in this study is illustrated as the fault detection. Fault scenarios such as stuck and scaling are considered and implemented using the nonlinear fuel cell simulation. Table 1 summarizes the set of faults implemented in this study as a case study.

Fault Scenarios
The ability to detect fault of the parity, the observer, and the Kalman filter models proposed in this study is illustrated as the fault detection. Fault scenarios such as stuck and scaling are considered and implemented using the nonlinear fuel cell simulation. Table 1 summarizes the set of faults implemented in this study as a case study.

Fault Detection and Isolation System
The first purpose of this section is to evaluate the detection performance of the proposed FDI scheme. The FDI scheme was tested under a realistic driving cycle; the Highway Fuel Economy Test (HWFET) was selected. The HWFET cycle, which has a high velocity profile, was converted and scaled based on a value of current density of 0.7 A/cm 2 . Before loading the HWFET profile, the temperature of the stack and the coolant inlet were set to 25 • C to ensure each sensor had the same initial state. The HWFET profile of the current, based on the current density, is plotted in Figure 6a. The corresponding temperature sensors of the stack and the coolant inlet in the fuel cell system are shown in Figure 6b In the simulation test, before implementing the fault scenarios, the alarm residual threshold under normal conditions needs to be determined. The thresholds were selected by calculating the standard deviation using the step load with three cases. The first case changed the load from 0.6 to 0.8 to 1.0, and 1.2 A/cm 2 . The second case changed the load from 0.8 to 1.0 to 0.6, and 0.9 A/cm 2 . The final cases changed the load from 1.0 to 0.8 to 0.6, and 0.4 A/cm 2 . The value of the standard deviation was derived by obtaining the mean value of the three residuals. Then, the values of the residual tolerance of the stack temperature sensor and the coolant inlet temperature sensor were obtained by multiplying the standard deviation by the allowable factor. The threshold was selected to be slightly greater than the maximum value of the allowed deviation in the no-fault tests. As a result, the residual threshold of the stack temperature was set at 5, −5; that of the coolant inlet temperature was also set at 5, −5, as shown in Table 2.  Figure 7, with the actual stack and coolant inlet temperature sensor outputs, residuals, and corresponding detection and isolation results that were recorded. First, the real stack temperature and real coolant inlet temperature were determined; it can be observed from Figure 7a,b that the real stack temperature value increased after 450 s, while the real coolant inlet temperature value was maintained at 60 • C despite the occurrence of downshoot at 450 s. When the stack sensor stuck fault occurs, the water pump does not operate because the stack temperature sensor output is lower than the desired temperature control value (70 • C). This explains why the real stack temperature increased.
It can be seen in Figure 7c that the responses of the parity and Kalman filter residuals of the stack temperature to this fault are highly sensitive, whereas the estimator residual responds only slightly to this fault. Moreover, it can be seen from Figure 7d that the responses of the parity and Kalman Filter residuals of the coolant inlet temperature to these faults are slight. On the other hand, the estimator residual does not respond to this fault. According to the residual characteristics listed in Table 2, fault f 1 will occur when the stack temperature sensor responds to this fault. The CUSUM test result L k of the residual confirms the presence of this fault, as shown in Figure 7e,f. This is because the stack CUSUM value, L k (r(s)), exceeds the allowed threshold value at the time 450 s. On the other hand, the coolant threshold H does not cross the allowed threshold value. Therefore, there is no sensor fault in the coolant inlet temperature sensor. Finally, the isolated fault f 1 is shown in Figure 7g.

Fault Detection and Isolation with Coolant Inlet Temperature Sensor Stuck
Similarly, coolant inlet temperature sensor detection and isolation in the case of a coolant inlet sensor stuck fault was evaluated under HWFET cycles based on the current density. Taking a fault resulting from the coolant inlet temperature sensor being stuck at 10 • C upon injection into the coolant inlet at time 450 s, the detection and isolation results are plotted in Figure 8. Figure 8a,b shows the value measured by the temperature sensor and the real temperature value of the stack and coolant inlet, respectively. Figure 8a,b indicates that the real stack temperature increased after 450 s, because the radiator fan stopped. Because the real coolant inlet temperature increases, the stack temperature increases continuously, even though the water pump speed increases. Thus, a faulty coolant inlet temperature sensor can cause critical damage to the fuel cell stack.
In addition, Figure 8c,d shows the generated residual values r obtained with the three respective methods. The figure shows that the parity and the Kalman filter residuals exhibit a shift drop after 450 s, after which the value is maintained until the end of the cycle. On the other hand, the estimator residual undergoes a distinct drop at approximately 450 s; then, the residual immediately becomes similar to what it was before fault injection. Therefore, the results of the CUSUM test, shown in Figure 8e,f, show that the threshold was exceeded at 450 s, indicating a fault in the coolant inlet temperature sensor. Consequently, the CUSUM result L k of the residual r i confirms this result as the threshold H is crossed at 450 s, as shown in Figure 8e. That is why even if the residuals in the stack are zero, the actual and sensor values are the same. This means that the stack sensor has not failed.
The isolated fault signal of the temperature sensor is plotted in Figure 8g.

Fault Detection and Isolation with Coolant Inlet Temperature Sensor Stuck
Similarly, coolant inlet temperature sensor detection and isolation in the case of a coolant inlet sensor stuck fault was evaluated under HWFET cycles based on the current density. Taking a fault resulting from the coolant inlet temperature sensor being stuck at 10 °C upon injection into the coolant inlet at time 450 s, the detection and isolation results are plotted in Figure 8. Figure 8a,b shows the value measured by the temperature sensor and the real temperature value of the stack and coolant inlet, respectively. Figure 8a,b indicates that the real stack temperature increased after 450 s, because the radiator fan stopped. Because the real coolant inlet temperature increases, the stack temperature increases continuously, even though the water pump speed increases. Thus, a faulty coolant inlet temperature sensor can cause critical damage to the fuel cell stack.
In addition, Figure 8c,d shows the generated residual values r obtained with the three respective methods. The figure shows that the parity and the Kalman filter residuals exhibit a shift drop after 450 s, after which the value is maintained until the end of the cycle. On the other hand, the estimator residual undergoes a distinct drop at approximately 450 s; then, the residual immediately becomes similar to what it was before fault injection. Therefore, the results of the CUSUM test, shown in Figure  8e,f, show that the threshold was exceeded at 450 s, indicating a fault in the coolant inlet temperature sensor. Consequently, the CUSUM result Lk of the residual ri confirms this result as the threshold H is crossed at 450 s, as shown in Figure 8e. That is why even if the residuals in the stack are zero, the actual and sensor values are the same. This means that the stack sensor has not failed.
The isolated fault signal of the temperature sensor is plotted in Figure 8g.
(a) This section presents an evaluation of the coolant inlet temperature sensor. The sensor experiences 20% off-scaling in terms of degrees of temperature at 450 s. The detection and isolation results are plotted in Figure 9. Figure 9a,b shows the temperature sensor measurements and the real temperature values of the stack and coolant inlet, respectively. From the temperature screening results, it can be seen that the value of the stack temperature sensor is the same as the real temperature. On the other hand, the real value of the coolant inlet temperature increased when a scaling fault occurred at the coolant inlet temperature sensor. However, the coolant inlet temperature was maintained by the radiator fan controller, even during sensor faults. Figure 9c,d shows the speed of the water pump and the radiator fan. Under normal conditions, the motor speed of the radiator fan decreased from its operating speed after 450 s. This is why the motor controller detected the scaled value of the coolant inlet temperature from the sensor. Hence, the controller reduced the motor speed because of the low heat capacity in the cooling water and, consequently, the real temperature of the coolant inlet increased. This, in turn, caused the water pump motor speed to increase in response to the high temperature at the coolant inlet.
It can be seen from Figure 9e,f that the parity and Kalman filter residuals of the stack responded to this fault in the form of a peak, in contrast to the estimator residual, which responded only slightly. Furthermore, Figure 9e shows that the residuals of the stack temperature hardly responded to this fault. Especially, the maximum peer value cannot be exceeded in the case of the three residuals. That is why the temperature measured by the stack temperature sensor recovered immediately because of the water pump. On the other hand, the estimator residual of the coolant inlet temperature did not respond to this fault, whereas the parity and the Kalman filter residuals exceeded the threshold value. According to the residual characteristics in Table 2, fault f 2 will be detected when the coolant inlet temperature sensor responds to this fault. Unlike the stuck case, the residuals of the stack and the coolant inlet exist in the form of a peak. This is why the controllers of the water pump and the radiator fan continue operating in the case of a sensor scaling fault. The CUSUM test results indicate that a fault occurred in the coolant inlet temperature sensor, as shown Figure 9g,h. The coolant inlet temperature L k (r(c)) crossed the allowed threshold value at 450 s, and L k (r(c)) was maintained. On the other hand, the coolant threshold did not cross the specified threshold value. Thus, a sensor fault in the coolant inlet temperature sensor was indicated. Finally, the isolated fault signal f2 is shown in Figure 9i.

Conclusions
To evaluate three methods of fault detection and isolation under a real driving load profile, a parity equation, a state observer, and a Kalman filter method have been designed. Also, this paper presents a comparison of three residuals derived from the parity equation, state observer, and Kalman filter, as well as a method for detection and isolation of faults in a cooling system temperature sensor in a lumped fuel cell system.
1. The parity equation, state observer, and Kalman filter methods for fault detection were developed in this study. The temperature outputs of a nonlinear fuel cell system were monitored; then, estimated linear model outputs were compared with those of the nonlinear fuel cell system to generate a residual. The generated residuals were evaluated by the CUSUM method to identify the presence of a sensor fault.

Conclusions
To evaluate three methods of fault detection and isolation under a real driving load profile, a parity equation, a state observer, and a Kalman filter method have been designed. Also, this paper presents a comparison of three residuals derived from the parity equation, state observer, and Kalman filter, as well as a method for detection and isolation of faults in a cooling system temperature sensor in a lumped fuel cell system.

1.
The parity equation, state observer, and Kalman filter methods for fault detection were developed in this study. The temperature outputs of a nonlinear fuel cell system were monitored; then, estimated linear model outputs were compared with those of the nonlinear fuel cell system to generate a residual. The generated residuals were evaluated by the CUSUM method to identify the presence of a sensor fault.