Next Article in Journal
Enhanced Cycling Performance of Spinel LiNi0.5Mn1.5O4 Cathodes through Mg-Mn Hetero-Valent Doping via Microwave Sol-Gel Method
Previous Article in Journal
Energy-Efficient Geopolymer Composites Containing Phase-Change Materials—Comparison of Different Contents and Types
Previous Article in Special Issue
Numerical Investigation of the Performance of a Proton Exchange Membrane Water Electrolyzer under Various Outlet Manifold Structure Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Feature Fusion Method for Life Prediction of Automotive Proton Exchange Membrane Fuel Cell Based on TCN-GRU

1
Hubei Key Laboratory of Advanced Technology for Automotive Components, Wuhan University of Technology, Wuhan 430070, China
2
Hubei Research Center for New Energy & Intelligent Connected Vehicle, Wuhan University of Technology, Wuhan 430070, China
3
National Energy Key Laboratory for New Hydrogen-Ammonia Energy Technologies, Foshan Xianhu Laboratory, Foshan 528200, China
*
Author to whom correspondence should be addressed.
Materials 2024, 17(19), 4713; https://doi.org/10.3390/ma17194713
Submission received: 19 July 2024 / Revised: 13 September 2024 / Accepted: 23 September 2024 / Published: 25 September 2024
(This article belongs to the Special Issue PEMFC Materials: Fabrication, Characterization and Applications)

Abstract

:
The Proton Exchange Membrane Fuel Cell (PEMFC) is a fast-developing battery technology, and the key to its reliability and lifespan improvement lies in the accurate assessment of durability. However, the degradation mechanism of the PEMFC is hard to determine and its internal parameters are highly coupled. Thus, the development of a more accurate life prediction model that meets the actual scenarios is to be investigated urgently. To solve this problem, a multi-feature fusion life prediction method based on the Temporal Convolutional Network-Gated Recurrent Unit (TCN-GRU) is proposed. A TCN algorithm is used as the prediction base model, and two GRU modules are included with the model to strengthen the model’s expression ability and improve its predictive accuracy. Two widely recognized datasets and two operating conditions are utilized for model training and prediction, respectively. Comparisons are made with single-feature parameter models in terms of Root Mean Square Error (RMSE) and the Determination Coefficient (R2). The results show that the prediction accuracy of the TCN-GRU multi-feature fusion model is higher than that of the single-feature models in terms of stability and anti-interference under both operating conditions. The accuracy of the TCN-GRU (three-feature) model is the most optimal in a steady-state condition at 80% of the training set ratio (RMSE = 3.27 × 10−3, R2 = 0.965). Furthermore, with the increase in the input feature parameter, the TCN-GRU model is closer to the real value, which proves once again that the proposed model can meet the accuracy requirements of the life prediction of the PEMFC.

1. Introduction

In recent years, with the aggravation of environmental problems and the energy crisis, more and more countries have realized the importance of clean energy in future energy development, and have started to research and develop infrastructures and energy storage devices around hydrogen energy [1,2,3]. The Proton Exchange Membrane Fuel Cell (PEMFC), a mainstream hydrogen energy power device, is a typical hydrogen fuel cell with high power density, fast startup, long sustainable operation time, and zero pollution and emission [4,5]. Therefore, it has the prospect of wide application in the power field of automobiles, power locomotives, and ships [6].
However, the cost and life of PEMFCs are two important factors limiting their large-scale commercial application [7], and an accurate assessment of durability is key to improving reliability and extending lifespan. Therefore, many scholars have devoted themselves to exploring the life prediction methods for PEMFCs. Nowadays, there are three main methods: model-based method, data-driven method, and hybrid prediction method [8].
The model-based prediction method mainly focuses on modeling from the fuel cell degradation mechanism, which requires less data and has high accuracy. Ou et al. [9] proposed a PEMFC prediction method based on a semi-empirical model to realize the prediction of PEMFC degradation and the estimation of its remaining service life under automotive environmental conditions by introducing electrochemical surface area and equivalent resistance degradation models, respectively. Lechartier et al. [10] presented a combined static and dynamic prediction model and verified the accuracy of the model with experimental data. Ao et al. [11] used a life prediction method based on a frequency-domain Kalman filter (FDKF) and a voltage degradation model. By processing the data in groups, the computation time can be greatly reduced with high accuracy. However, the model-based method requires an in-depth understanding of the aging mechanism of the stack and strong modeling ability. Meanwhile, the internal structure and materials of different stacks are not the same, and the aging mechanism is somewhat different, so some key parameters in the modeling cannot be defined directly, which makes it difficult to establish a complete and accurate mechanism model.
The data-based prediction methods are mainly used to build a system behavior model on PEMFC historical operation data directly for fault diagnosis and life prediction. These methods do not require in-depth analysis of the internal reaction mechanism of the fuel cell, but rather comprehensive analysis of a large amount of historical operation data using a variety of methods, such as statistical modeling, machine learning, deep learning, and hybrid learning [12]. Wu et al. [13] improved on the Relevance Vector Machine (RVM) and introduced the Support Vector Machine (SVM) algorithm for comparison, and the results proved that the improved RVM statistical model can form an adaptation to the prediction process due to the absence of limitations in the kernel function, and the prediction performance is better. Liu et al. [14] compared several life prediction methods based on different structural neural networks, and the results showed that the Adaptive Neuro-Fuzzy Inference System-Fuzzy C-Means (ANFIS-FCM) has the best short-term prediction performance, while the introduction of the Particle Swarm Optimization (PSO) algorithm realizes the automatic adjustment of parameters of ANFIS-FCM. Mezzi et al. [15] designed an echo state network-based variable load prediction method for fuel cells that can be predicted without prior knowledge of the variable load profile. Zuo et al. [16] combined the attention mechanism and GRU for PEMFC prediction, and the results showed that the model has high prediction accuracy on both dynamic and pseudo-stable datasets. However, the data-based prediction method still has the problem of “black box” for the degradation mechanism and state changes of the PEMFC, and there are some limitations in relying on the prediction results to decide the follow-up maintenance measures.
The hybrid prediction methods combine different life prediction methods to improve prediction accuracy by improving individual prediction method weaknesses. There are mainly model–data hybrid-driven methods and data–data hybrid-driven methods [17]. The model–data hybrid-driven method retains the interpretability of the parameters in the model-based method, but also makes full use of the advantages of the two models by blurring them with the help of the data-based method when the mechanism process is not clear. Pan et al. [18] combined model-based AEKF and data-driven NARX with external inputs for predicting PEMFC performance degradation and validated the predictive power of the method using two different datasets. Liu et al. [19] proposed a 2-Stage hybrid prediction method: Stage 1 used an automated machine learning algorithm based on an evolutionary algorithm and an adaptive neuro-fuzzy inference system to achieve voltage degradation prediction in a long time series state. Stage 2 utilized a semi-empirical degradation model based on the predicted data to estimate the remaining lifetime of the battery. The data–data hybrid-driven method can fully utilize the different dimensional information of the algorithm for prediction, thereby improving the prediction accuracy and robustness. Zhu et al. [20] combined Bayesian theory and the self-attention mechanism to propose a B-GRU hybrid model, which can quantify the uncertainty parameters in the prediction process, and experimentally proved that the hybrid algorithm’s computational accuracy is higher than that of the commonly used neural networks when the training data are less than 380 h.
The degradation of the PEMFC’s internal structure is often not accurately characterized by a single feature with poor robustness; the single feature may change due to the performance degradation of the stack, which in turn affects the output characteristics of the stack. And due to the strong coupling of the internal parameters of the PEMFC, it is often difficult to determine the strong correlation parameters using multiple features to predict the degradation trend, which leads to a long prediction time and low prediction accuracy. The Temporal Convolutional Network-Gated Recurrent Unit (TCN-GRU) is an algorithm that combines the ability of the TCN to extract features of time series data with the advantage of GRU for higher computational efficiency in long time series data prediction, which is able to better handle multi-feature data inputs [21], and has been demonstrated in electrical load [22] and lithium-ion battery life prediction [23]. Therefore, this paper proposes a multi-feature fusion method for the life prediction of automotive PEMFCs based on the TCN-GRU. This method is also the first application of the TCN-GRU algorithm to the life prediction of the PEMFC. The correlation of feature parameters in PEMFC experiments is analyzed, and the strength of the correlation between other feature parameters and voltage is projected, so that multiple strongly correlated features are filtered to participate in the fusion life prediction. The TCN-GRU dual-feature and three-feature fusion life prediction models are constructed, and the voltage trend prediction is carried out by utilizing the FCLAB test open dataset and the dynamic cycling operating condition test dataset of the automotive PEMFC system, comparing it with the evaluation indicators to judge the applicability of the models.
The rest of paper is structured as follows. In Section 2, for the characteristics of steady-state operating conditions and dynamic cycling conditions, correlation analysis is performed on the original data, and several strongly correlated feature parameters are filtered and smoothed for noise reduction and reconstruction. Section 3 presents the single-feature parameter and multi-feature parameters fusion prediction models. The life prediction under different operating conditions is carried out in Section 4, and the model evaluation results are analyzed and summarized in Section 5. The conclusions and prospects for future work are summarized in Section 6.

2. Multi-Feature Correlation Analysis and Data Processing

2.1. Source of Data

The aging process of internal components of the PEMFC such as the proton exchange membrane and bipolar plate is relatively slow. Meanwhile, due to the continuity of the durability experiments and the closure of the power stack, it is not possible to find a reasonable performance index from the internal perspective, so the external index is used as a symbol for characterization [24]. In order to accurately predict the aging trend of the PEMFC in steady-state conditions, the first PEMFC dataset group (FC1) from the 2014 IEEE Prognostics and Health Management (PHM) Data Challenge durability test public dataset is used [25]. The fuel cell test bench used for this data is shown in Figure 1. The specific parameters of this stack are shown in Table 1. The monitoring parameters include air/hydrogen flow and pressure, coolant temperature, output current, output voltage, monomer voltage, and so on.
The dynamic operating conditions data come from the cycling conditions simulation experiment completed on the 60 kW fuel cell bench test platform [27]. The dynamic cycling test conditions are based on the China light-duty vehicle test cycle-passenger car (CLTC-P) driving conditions [28], with a single cycle duration of 1800 s, divided into three speed intervals: low, medium, and high. The three speed ranges are 10–50 km/h (low range), 50–80 km/h (medium range), and 80–114 km/h (high range). The test rig cumulatively tested and sampled data for 96 complete cycles with a data sampling interval of 2 s. The test bench is shown in Figure 2. The schematic of this fuel cell system is shown in Figure 3, which mainly consists of the stack, the air system, the hydrogen system, the thermal management system, and the control system. The specific parameters of this stack are shown in Table 2. In this experiment, the anode of the stack was supplied with hydrogen in a cyclic mode, with the hydrogen pressure maintained at about 15 bar and the pressure controlled by a proportional valve to control the flow area. The unreacted hydrogen is returned to the anode inlet via a hydrogen circulation pump and an ejector. The cathode air is pressurized, cooled, and humidified, and then flows into the cathode flow channel. The cathode inlet pressure and flow rate can be adjusted by the air compressor, backpressure valve, and exhaust throttle. The parameters that can be directly measured on the test platform include: output voltage, output power, cathode/anode inlet/exhaust pressures, air compressor speed, coolant inlet/outlet temperatures, proportional valve and backpressure valve openings, hydrogen circulating pump speed, and so on.

2.2. Correlation Coefficient Analysis Method

Output voltage can be a key feature to characterize the degradation of the PEMFC [29]. Therefore, when selecting other features for joint prediction, it is necessary to filter out the feature parameters that are closely related to the output voltage. There are three main methods for performing correlation analysis: the Pearson correlation coefficient, the Spearman correlation coefficient, and the Kendall correlation coefficient [30].
The Pearson correlation coefficient is applicable to the comparison between two continuous variables and measures the strength and direction of their linear relationship, with a value between −1 and 1. The Pearson correlation coefficient of variables X and Y is denoted as ρ X , Y , as shown in Equation (1) [31]:
ρ X , Y = E ( ( X μ X ) ( Y μ Y ) ) σ X σ Y = E ( X Y ) E ( X ) E ( Y ) E ( X 2 ) E 2 ( X ) E ( Y 2 ) E 2 ( Y )
here, µ is the expected value and σ is the standard deviation.
When ρ X , Y = 1, it indicates that there is a completely positive linear relationship between variables X and Y, i.e., the two variables are completely positively correlated. When ρ X , Y = −1, it means that there is a completely negative linear relationship between variables X and Y, i.e., the two variables are completely negatively correlated. When ρ X , Y = 0, it means that there is no linear relationship between variables X and Y, i.e., the two variables are uncorrelated.
The Spearman correlation coefficient is used to determine the monotonic relationship between two variables, which does not require the data to follow a normal distribution, and it calculates the correlation by ranking the variables, which is applicable in the case of a nonlinear relationship. The rank of a number in a variable is the position of that number in the set of columns from smallest to largest. This method is particularly suitable for rank order data. The Spearman correlation coefficient of variables X and Y is denoted as ρ , which is calculated as shown in Equation (2) [32]:
ρ = 1 6 i = 1 n d i 2 n ( n 2 1 )
here, d i is the rank difference obtained by subtracting Y i from X i and n is the number of data points. This method for the degree of linear relationship is consistent with the Pearson correlation coefficient.
The Kendall correlation coefficient is a statistic used to measure the hierarchical relationship between two variables, which does not depend on the specific values of the data, but rather calculates the correlation based on the hierarchy of the data. The Kendall correlation coefficient is particularly applicable when the data do not satisfy a normal distribution or there are outliers. The Kendall correlation coefficient between two variables is denoted as τ , which is calculated as shown in Equation (3) [33]:
τ = N c N d n × ( n 1 ) / 2
here, N c is the number of two variables that agree in rank, N d is the number of two variables that do not agree in rank, and n is the number of elements.
The feature parameters of the PEMFC are coupled with each other, but do not have a linear relationship with each other, and at the same time, there is no normal distribution law. Therefore, in this paper, the Spearman correlation coefficient and Kendall correlation coefficient are selected as multi-feature filtering methods.

2.3. Multi-Feature Filtering and Data Processing

Since the monitoring parameters of the steady-state condition dataset and the dynamic cycling condition dataset are not exactly the same, and the operating environments are different. Therefore, the multi-feature correlation analysis and data processing of the two conditions are performed separately.

2.3.1. Steady-State Condition Multi-Feature Filtering and Data Processing

In the FCLAB FC1 group dataset, the test bench monitors 24 feature parameters. Not all of them affect the degradation of the PEMFC and some of the parameters are repetitive in properties, so the single voltages U 1 U 5 , the current density i , and the current I are excluded. Since hydrogen and air are compressible fluids, the relationship between their pressures and flow rates cannot be directly described by Bernoulli’s equation and they cannot be replaced with each other. Therefore, both pressures and flow rates are filtered into the feature parameters to be analyzed. Finally, a total of 17 feature parameters are retained, as shown in Table 3.
The correlation analysis of 17 feature parameters in the FC1 group dataset is shown in Figure 4. When the absolute value of the correlation coefficient is 0.8–1.0, it is called a strong correlation. When the absolute value of the correlation coefficient is 0.3–0.8, it is called weak correlation. When the absolute value of the correlation coefficient is less than 0.3, there is no correlation between the two values [34]. Comparing with the Figure 4a Spearman results, the ρ values of D o u t A I R and T i n H 2 with U t o t are 0.67 and −0.51, respectively, which are higher than the correlation of other feature parameters, and the same results are found in the Figure 4b Kendall correlation coefficient. Therefore, D o u t A I R , T i n H 2 , and U t o t are selected as the feature parameters for the life prediction of the steady-state operating conditions.
If the original FC1 group dataset is used directly for life prediction, it is computationally intensive and highly susceptible to overfitting, so the dataset needs to be reconstructed. Firstly, interval sampling is performed, and the output voltage values are sampled every 70 s. In order to reduce the effect of noise and improve the prediction accuracy, the interval sampling data are smoothed by using a one-dimensional Discrete Wavelet Transform (DWT) [35,36]. The three selected feature parameters before and after the noise reduction and smoothing are shown in Figure 5. As can be seen in Figure 5, the voltage information after the noise reduction process retains the original data change trend, with fewer data anomalies and being overall smoother, which is conducive to model prediction.

2.3.2. Dynamic Cycling Condition Multi-Feature Filtering and Data Processing

The parameters that can be directly monitored in the 60 kW fuel cell bench test platform include output voltage, output power, and more than ten feature parameters. The feature parameters that are not related to the degradation of the PEMFC are also excluded. The final 16 feature parameters retained are shown in Table 4.
Correlation analysis is performed on 16 feature parameters in the dynamic cycling condition dataset, which are shown in Figure 6. In Figure 6a, it can be found that U m i n , P i n H 2 , P i n A I R , n A C P S p d , and P o u t A I R are strongly correlated with U o u t , but the difference is not obvious. From comparing the correlation coefficient in Figure 6b, it can be found that the absolute values of U m i n and P i n A I R are closer to 1, which indicates a strong correlation between these two variables and U o u t . Therefore, U o u t , U m i n , and P i n A I R are selected as the features parameters for the life prediction of the dynamic cycling conditions.
In the actual operation process, due to the PEMFC being affected by the environment, temperature, operating conditions, its own system, and other factors, the output data collected by the dynamic cycling condition fluctuate, resulting in a large number of outliers [37]. If it is divided by the operating condition trend or data threshold, there may be data overlapping the calculation, which affects the prediction accuracy. Conditional Generative Adversarial Networks (CGANs) is an extension of Generative Adversarial Networks (GANs), which allow the model to introduce conditional information in the process of generating data, so that the generator can generate samples of specific categories or features based on pre-set conditions [38]. Therefore, the original data containing high-frequency noise are reconstructed using the CGAN unsupervised model to make the data preprocessing and feature extraction process adaptive and self-learning. The data before and after processing are shown in Figure 7. It can be found in Figure 7 that all three feature parameters are able to retain the change trend of the original data after data reconstruction, while the noise is significantly suppressed and the curve is smoother.

3. Design of Life Prediction Model

The commonly used deep learning algorithms and their usage scenarios are shown in Table 5.

3.1. Design of Single-Feature Models

The prediction method with output voltage as a single feature is added to facilitate the comparison and analysis of the advantages and disadvantages of the proposed multi-feature life prediction methods. For the steady-state operating condition life prediction, two prediction methods, the traditional RNN and LSTM, are selected for comparison based on the algorithm characteristics. The RNN model is built based on TensorFlow, and in order to prevent the prediction from overfitting, the Dropout technique is introduced. This technique randomly removes some neurons and their single-step connections during the training process and at each iteration, so that the whole training process is “sparsely” sampled, which avoids the same network from being trained repeatedly and ensures the model’s generalization ability [39]. Meanwhile, in order to ensure that the learning rate η can be adjusted adaptively and maintain the adaptability and stability of each parameter, Adam is chosen as the optimizer, with the number of neurons in the hidden layer being 150. The probability of Dropout random deletion is 0.4, the activation function is ReLU, and the number of iterations is 2000. Considering the training speed and generalization ability, the Batch Size is formulated to be 64 at each iteration.
Based on the LSTM algorithm structure, the Dropout technique is also introduced to prevent predictive overfitting. According to the conclusions of Zhang et al. [40], the Learning Rate η of 0.01 and Dropout random deletion probability of 0.4 are selected to train the network. The number of neurons in the hidden layer of the LSTM is 150, the number of iterations is 2000, and the learning rate is decayed for 50 generations, with a learning rate decay rate of 0.2.
For dynamic cycling operating condition life prediction, a CNN-LSTM-based PEMFC voltage prediction model is selected based on the algorithmic characteristics with the same configuration as above [41,42]. The model utilizes the advantages of the data feature extraction capability of the CNN algorithm and the suitability of the LSTM algorithm for dealing with long-term dependencies. The CNN is mainly used to extract the underlying features of the data to achieve dimensionality reduction, and the use of the LSTM network is mainly to grasp the contextual information more deeply and comprehensively, to avoid the loss of contextual information due to the excessive depth of the CNN network; at the same time, it can allow the model to focus more on the temporal attributes of the data.

3.2. Design of Multi-Feature Models

Since the essence of multi-feature fusion prediction is the processing of time series features, TCN is selected as the prediction base model in this paper. In order to improve the prediction accuracy of multi-feature data input, the TCN structure is modified, and the improved TCN multi-input model structure is shown in Figure 8. The TCN activation function is replaced with the Leaky ReLU. The Leaky ReLU function introduces a negative slope on the base of ReLU, as shown in Equation (4). Because the Leaky ReLU function has a non-zero gradient, it allows faster convergence during training and avoids the problems of neurons not being updated and gradient vanishing. In order to enhance the model’s capacity and expressiveness [43], three layers of residual blocks are used in the TCN model to help the model capture the long-term correlation relationships within the time series more accurately, while avoiding too many residual blocks leading to model overfitting.
f ( x ) = x ,   i f   x > 0 0.01 x ,   o t h e r w i s e
Compared to the LSTM model using three gating mechanisms, the GRU structure is simpler and has fewer hyperparameters. Hence, it has higher computational efficiency in long time series data prediction and is easier to adjust the model hyperparameters and structure [44]. Similarly, in order to better capture multiple sequence features and strengthen the model’s expression ability, two GRU modules are added to the TCN-GRU model. The improved TCN-GRU model is shown in Figure 9. After repeated experimental tuning, the appropriate hyperparameters of the TCN-GRU multi-feature input model are determined, as shown in Table 6.

4. Multi-Operating Condition Life Prediction for PEMFCs

In this paper, predictions are made according to the Training Set: Test Set = 4:6 (40%), 6:4 (60%), 8:2 (80%) division. The influence of the experimental environment configuration on the results during the prediction process is minimized as much as possible.

4.1. Life Prediction for Steady-State Operating Condition

According to the multi-feature correlation filtering results in the previous section, U t o t and D o u t A I R are selected as the input values for the dual-feature prediction model. U t o t , D o u t A I R , and T i n H 2 are selected as the input values for the three-feature fusion prediction model.
The datasets of U t o t and D o u t A I R are imported into the input layer of the TCN-GRU model, and the prediction results under the conditions of 40%, 60%, and 80% of the training set ratio are obtained, respectively, as shown in Figure 10. By analyzing the prediction curves in the figures, it can be found that the results of the TCN-GRU dual-feature fusion prediction model are significantly better than those of the LSTM and RNN models, and the dual-feature fusion prediction fitting performance is better as the number of training sets increases.
The datasets of U t o t ,     D o u t A I R , and T i n H 2 are imported into the input layer of the TCN-GRU model, and the prediction results are obtained under the conditions of 40%, 60%, and 80% of the training set ratio, respectively, as shown in Figure 11. Analyzing the prediction curves in the figure, it can be observed that the three-feature fusion prediction results are close to the original data at 40%, 60%, and 80% of the training set ratio, indicating that the TCN-GRU model can characterize the degradation inside the stacks more rapidly and accurately and make a prediction after more stacks’ feature parameters are introduced.

4.2. Life Prediction for Dynamic Cycling Condition

According to the multi-feature correlation filtering results in the previous section, U o u t and U m i n are selected as the input values for the dual-feature prediction model. U o u t , U m i n , and P i n A I R are selected as the input values for the three-feature fusion prediction model.
The datasets of U o u t and U m i n are imported into the input layer of the TCN-GRU model, and the prediction results under the conditions of 40%, 60%, and 80% of the training set ratio are obtained, respectively, as shown in Figure 12. The dynamic cycling condition of the automotive PEMFC has more uncertainties compared with the steady-state operating condition. Theoretically, the multi-feature fusion prediction has better robustness by combining the trends of multi-feature parameters to comprehensively analyze the degradation characteristics of the PEMFC, and weakening the influence of single-feature parameter anomalies and noise. By qualitatively analyzing Figure 12, it can be seen that the dual-feature fusion prediction is more desirable for dynamic life prediction. From the local zoom in the figure, it can be noticed that the dual-feature prediction curve can better reflect the trend of the voltage value compared with the CNN-LSTM.
The datasets of U o u t ,   U m i n , and P i n A I R are imported into the input layer of the TCN-GRU model, and the prediction results are obtained under the conditions of 40%, 60%, and 80% of the training set ratio, respectively, as shown in Figure 13. As can be seen from the figure, the triple-feature fusion prediction is slightly better than the dual-feature fusion prediction, especially when the training set ratio is at 60%. From the local zoom in the figure, the three-feature fusion prediction is closer to the dual-feature fusion prediction, and the fitting effect is better. As the training data volume is larger, the prediction trend for the remaining lifetime is more accurate. At 80% of the training set ratio, the prediction results are almost the same as the original data.

5. Results and Discussion

In this paper, predictions are made according to the Training Set: Test Set = 4:6 (40%), 6:4 (60%), 8:2 (80%) division. The influence of the experimental environment configuration on the results during the prediction process is minimized as much as possible. In order to quantitatively assess model performance to visualize the degree of model fitting and prediction ability, it is necessary to use appropriate model evaluation indicators to select, adjust, interpret, and improve the model [45]. According to the characteristics of operating conditions, this paper adopts Root Mean Square Error (RMSE) and the Determination Coefficient (R2) as the evaluation indexes of model prediction results.

5.1. Analysis of Life Prediction Results for Steady-State Operating Condition

The prediction evaluation indicators RMSE and R2 of the RNN, LSTM, TCN-GRU (dual-feature), and TCN-GRU (three-feature) models are shown in Table 7 and Figure 14. It can be concluded as follows:
  • The TCN-GRU (three-feature) has the highest prediction accuracy, with an RMSE value of 3.27 × 10−3 and an R2 value of 0.965 at 80% of the training set ratio, almost perfectly fitting the original data, which is significantly better than the other two prediction models. It indicates that the multi-feature fusion prediction model’s anti-overfitting and generalization ability is better than the RNN and LSTM models, and it is more suitable for the prediction of long time series and data with multiple parameters.
  • The TCN-GRU (three-feature) prediction results are better than the TCN-GRU (dual-feature), compared with the RMSE value reduced by at least 3.82%, and even reduced by 23.6% in the training set ratio of 80%, which indicates that in the long time series prediction, the more features involved in the fusion of the prediction, the greater the model’s ability to resist the abnormal data points or noise, and the better the stability of the model.
  • Comparing the TCN-GRU multi-feature fusion prediction results horizontally, with the increase in the training set, the RMSE value of the TCN-GRU (three-feature) decreases by 12.3% and 26.02% in turn, the R2 value improves by 1.18% and 2.33%, and the prediction effect becomes better with the increase in the training set, which proves the multi-feature fusion prediction model can meet the requirements of the PEMFC life prediction.

5.2. Analysis of Life Prediction Results for Dynamic Cycling Condition

The prediction evaluation indicators RMSE and R2 of the CNN-LSTM, TCN-GRU (dual-feature), and TCN-GRU (three-feature) models are shown in Table 8 and Figure 15. It can be concluded as follows:
  • The TCN-GRU (three-feature) has the highest prediction accuracy, which is slightly better than the TCN-GRU (dual-feature). Comparing the RMSE values, the TCN-GRU (three-feature) model reduces at least 6.79% compared to the CNN-LSTM model, which indicates that the multi-feature fusion prediction model is more resistant to noise and has better robustness in the long time series prediction of dynamic cycling conditions.
  • Comparing the R2 values, it can be found that the fitting ability of the TCN-GRU multi-feature fusion prediction model has a disconnected improvement compared with the traditional deep learning network. It further shows that a single deep learning network cannot meet the requirements of time series prediction with large data volume, and it is necessary to flexibly combine multiple algorithms, give full play to the advantages of each algorithm, and reasonably build a joint model according to the data volume and data characteristics, in order to achieve better prediction results.

6. Conclusions

In this paper, a multi-feature fusion method for life prediction of the PEMFC based on the TCN-GRU is proposed. This model is based on the TCN algorithm that possesses a three-layer residual block, and two GRU modules are added to better capture multiple sequential features and strengthen the model expression capability. To evaluate the model, two commonly used datasets are used for model training under steady-state operating conditions and dynamic cycling conditions, respectively. Strongly correlated feature parameters are imported into the model for prediction after correlation filtering, and RMSE and R2 are chosen as the evaluation indicators of the prediction effect. The main conclusions are summarized as follows:
  • In the steady-state operating condition, the TCN-GRU (three-feature) has the highest prediction accuracy, with RMSE values of 5.04 × 10−3, 4.42 × 10−3, and 3.27 × 10−3, and R2 values of 0.932, 0.943, and 0.965 at 40%, 60%, and 80% of the training set ratio, respectively. Both multi-feature models have better accuracy than the single-feature models. In dynamic cycling conditions, both the TCN-GRU (three-feature) and the TCN-GRU (dual-feature) have better accuracy than the CNN-LSTM. TCN-GRU (three-feature) has the highest prediction accuracy, especially at 60% of training set ratio, with an RMSE value and R2 value of 4.04 × 10−3 and 0.890, respectively. Compared with the single-feature prediction model, the multi-feature model has many advantages such as strong anti-interference, stability, and high prediction accuracy.
  • The results of TCN-GRU multi-feature fusion prediction are compared horizontally. In the steady-state operating condition, the RMSE values of the TCN-GRU (three-feature) decrease by 12.3% and 26.02%, and the R2 values improve by 1.18% and 2.33%, in turn, with the increase in the training set. While in dynamic cycling conditions, the RMSE values of the TCN-GRU (three features) reduce by 22.6% and 29.79%, and the R2 values improve by 0.6% and 3.8%, respectively. It is fully proved that the proposed multi-feature life prediction method has a better prediction effect as the training set increases, and the prediction model can meet the requirements of PEMFC life prediction.
Therefore, our proposed multi-feature fusion prediction method is suitable for PEMFC lifetime prediction in both steady-state operating conditions and dynamic cycling conditions, with better accuracy than single-feature prediction methods. This method is of significance for improving the accuracy of PEMFC life prediction, extending its service period, and enhancing its reliability and stability. While our study shows that the TCN-GRU outperforms the other three algorithms, since the accuracy of the algorithm depends on the input feature parameters, we have neglected the existence of measurement errors and the influence of other factors on the feature parameters in the actual measurements. Therefore, it is our future research direction to improve the accuracy of the TCN-GRU by reducing the influence of other factors on the accuracy of feature parameters. At the same time, we will continue to try to apply this algorithm to different devices and different complex scenarios, to improve the calculation speed and accuracy, and to realize accurate life prediction and faster and earlier fault diagnosis.

Author Contributions

Conceptualization, J.Z., F.Y., C.D., Y.Z., C.Z., J.W. and B.C.; methodology, J.Z., F.Y., C.D. and Y.Z.; software, J.Z., Y.Z., C.Z. and J.W.; validation, J.Z., Y.Z. and J.W.; formal analysis, F.Y. and C.D.; investigation, Y.Z. and C.Z.; resources, J.Z., C.D., Y.Z. and B.C.; data curation, J.Z., C.D., Y.Z. and C.Z.; writing—original draft preparation, J.Z., Y.Z. and C.Z.; writing—review and editing, J.Z., F.Y., C.D. and J.W.; visualization, J.Z., F.Y., C.D., Y.Z., C.Z., J.W. and B.C.; supervision, F.Y., C.D. and B.C.; project administration, F.Y., C.D. and B.C.; funding acquisition, F.Y. and C.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program of China, grant number 2022YFB4003703; Foshan Xianhu Laboratory of the Advanced Energy Science and Technology Guangdong Laboratory, grant number XHRD2024-11233100-01; Key R&D project of Hubei Province China, grant number 2021AAA006.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, W.; Sun, L.; Li, Z.; Fujii, M.; Geng, Y.; Dong, L.; Fujita, T. Trends and future challenges in hydrogen production and storage research. Environ. Sci. Pollut. Res. 2020, 27, 31092–31104. [Google Scholar] [CrossRef]
  2. Wang, J.; Du, C.; Yan, F.; Zhou, Q.; Xu, H. Hierarchical Rewarding Deep Deterministic Policy Gradient Strategy for Energy Management of Hybrid Electric Vehicles. IEEE Trans. Transp. Electrif. 2024, 10, 1802–1815. [Google Scholar] [CrossRef]
  3. Guo, L.; Zhao, L.; Jing, D.; Lu, Y.; Yang, H.; Bai, B.; Zhang, X.; Ma, L.; Wu, X. Solar hydrogen production and its development in China. Energy 2009, 34, 1073–1090. [Google Scholar] [CrossRef]
  4. Zhang, J.; Yan, F.; Du, C.; Li, W.; Fang, H.; Shen, J. Model-based performance optimization of thermal management system of proton exchange membrane fuel cell. Energies 2023, 16, 3952. [Google Scholar] [CrossRef]
  5. Sayed, K.; Abdel-Khalek, S.; Zakaly, H.M.H.; Aref, M. Energy Management and Control in Multiple Storage Energy Units (Battery–Supercapacitor) of Fuel Cell Electric Vehicles. Materials 2022, 15, 8932. [Google Scholar] [CrossRef]
  6. He, B.; Wei, L.; Hu, F.; Dafalla, A.; Guo, J.; Wang, C.; Jiang, F. Analyzing Key Factors Influencing Water Transport in Open Air-Cooled PEM Fuel Cells. Materials 2024, 17, 3267. [Google Scholar] [CrossRef]
  7. Chen, K.; Laghrouche, S.; Djerdir, A. Aging prognosis model of proton exchange membrane fuel cell in different operating conditions. Int. J. Hydrogen Energy 2020, 45, 11761–11772. [Google Scholar] [CrossRef]
  8. Hua, Z.; Zheng, Z.; Pahon, E.; Péra, M.-C.; Gao, F. A review on lifetime prediction of proton exchange membrane fuel cells system. J. Power Sources 2022, 529, 231256. [Google Scholar] [CrossRef]
  9. Ou, M.; Zhang, R.; Shao, Z.; Li, B.; Yang, D.; Ming, P.; Zhang, C. A novel approach based on semi-empirical model for degradation prediction of fuel cells. J. Power Sources 2021, 488, 229435. [Google Scholar] [CrossRef]
  10. Lechartier, E.; Laffly, E.; Péra, M.-C.; Gouriveau, R.; Hissel, D.; Zerhouni, N. Proton exchange membrane fuel cell behavioral model suitable for prognostics. Int. J. Hydrogen Energy 2015, 40, 8384–8397. [Google Scholar] [CrossRef]
  11. Ao, Y.; Laghrouche, S.; Depernet, D.; Chen, K. Proton Exchange Membrane Fuel Cell Prognosis Based on Frequency-Domain Kalman Filter. IEEE Trans. Transp. Electrif. 2021, 7, 2332–2343. [Google Scholar] [CrossRef]
  12. Deng, S.; Zhang, J.; Zhang, C.; Luo, M.; Ni, M.; Li, Y.; Zeng, T. Prediction and optimization of gas distribution quality for high-temperature PEMFC based on data-driven surrogate model. Appl. Energy 2022, 327, 120000. [Google Scholar] [CrossRef]
  13. Wu, Y.; Breaz, E.; Gao, F.; Miraoui, A. A Modified Relevance Vector Machine for PEM Fuel-Cell Stack Aging Prediction. IEEE Trans. Ind. Appl. 2016, 52, 2573–2581. [Google Scholar] [CrossRef]
  14. Liu, H.; Chen, J.; Hissel, D.; Su, H. Short-Term Prognostics of PEM Fuel Cells: A Comparative and Improvement Study. IEEE Trans. Ind. Electron. 2019, 66, 6077–6086. [Google Scholar] [CrossRef]
  15. Mezzi, R.; Yousfi-Steiner, N.; Péra, M.C.; Hissel, D.; Larger, L. An Echo State Network for fuel cell lifetime prediction under a dynamic micro-cogeneration load profile. Appl. Energy 2021, 283, 116297. [Google Scholar] [CrossRef]
  16. Zuo, J.; Lv, H.; Zhou, D.; Xue, Q.; Jin, L.; Zhou, W.; Yang, D.; Zhang, C. Deep learning based prognostic framework towards proton exchange membrane fuel cell for automotive application. Appl. Energy 2021, 281, 115937. [Google Scholar] [CrossRef]
  17. Wang, Y.; Wu, K.; Zhao, H.; Li, J.; Sheng, X.; Yin, Y.; Du, Q.; Zu, B.; Han, L.; Jiao, K. Degradation prediction of proton exchange membrane fuel cell stack using semi-empirical and data-driven methods. Energy AI 2023, 11, 100205. [Google Scholar] [CrossRef]
  18. Pan, R.; Yang, D.; Wang, Y.; Chen, Z. Performance degradation prediction of proton exchange membrane fuel cell using a hybrid prognostic approach. Int. J. Hydrogen Energy 2020, 45, 30994–31008. [Google Scholar] [CrossRef]
  19. Liu, H.; Chen, J.; Hissel, D.; Su, H. Remaining useful life estimation for proton exchange membrane fuel cells using a hybrid method. Appl. Energy 2019, 237, 910–919. [Google Scholar] [CrossRef]
  20. Zhu, W.; Guo, B.; Li, Y.; Yang, Y.; Xie, C.; Jin, J.; Gooi, H.B. Uncertainty quantification of proton-exchange-membrane fuel cells degradation prediction based on Bayesian-Gated Recurrent Unit. eTransportation 2023, 16, 100230. [Google Scholar] [CrossRef]
  21. Hu, Y.; Wang, N.; Lyu, L.; Zhou, X.; Fang, M. Application of Seq2Seq Model Based on TCN-GRU to Multivariate Water Quality Time Series Prediction. In Proceedings of the 2023 12th International Conference of Information and Communication Technology (ICTech), Wuhan, China, 14–16 April 2023; pp. 201–205. [Google Scholar]
  22. Hu, X.; Li, B.; Shi, J.; Li, H.; Liu, G. A Novel Forecasting Method for Short-term Load based on TCN-GRU Model. In Proceedings of the 2021 IEEE International Conference on Energy Internet (ICEI), Southampton, UK, 27–29 September 2021; pp. 79–83. [Google Scholar]
  23. Li, L.; Li, Y.; Mao, R.; Li, L.; Hua, W.; Zhang, J. Remaining Useful Life Prediction for Lithium-Ion Batteries with a Hybrid Model Based on TCN-GRU-DNN and Dual Attention Mechanism. IEEE Trans. Transp. Electrif. 2023, 9, 4726–4740. [Google Scholar] [CrossRef]
  24. Zhang, C.; Hu, H.; Ji, J.; Liu, K.; Xia, X.; Nazir, M.S.; Peng, T. An evolutionary stacked generalization model based on deep learning and improved grasshopper optimization algorithm for predicting the remaining useful life of PEMFC. Appl. Energy 2023, 330, 120333. [Google Scholar] [CrossRef]
  25. Gouriveau, R.; Hilairet, M.; Hissel, D.; Jemei, S.; Jouin, M.; Lechartier, E.; Morando, S.; Pahon, E.; Pera, M.; Zerhouni, N. IEEE PHM 2014 data challenge: Outline, experiments, scoring of results, winners. In Proceedings of the IEEE Conference on. Prognostics Health Management, Spokane, Washington, DC, USA, 22–25 June 2014; pp. 1–6. [Google Scholar]
  26. Chen, K.; Laghrouche, S.; Djerdir, A. Degradation prediction of proton exchange membrane fuel cell based on grey neural network model and particle swarm optimization. Energy Convers. Manag. 2019, 195, 810–818. [Google Scholar] [CrossRef]
  27. Li, X.; Wei, H.; Du, C.; Shi, C.; Zhang, J. Control strategy for the anode gas supply system in a proton exchange membrane fuel cell system. Energy Rep. 2023, 10, 4342–4358. [Google Scholar] [CrossRef]
  28. GB/T 38146.1-2019; China Automotive Test Cycle-Part 1: Light-Duty Vehicles. China Automotive Technology and Research Center Co., Ltd.: Tianjin, China, 2019; p. 36.
  29. Jouin, M.; Gouriveau, R.; Hissel, D.; Péra, M.; Zerhouni, N. Prognostics of PEM fuel cell in a particle filtering framework. Int. J. Hydrogen Energy 2014, 39, 481–494. [Google Scholar] [CrossRef]
  30. Wang, J.; Du, C.; Yan, F.; Duan, X.; Hua, M.; Xu, H.; Zhou, Q. Energy Management of a Plug-in Hybrid Electric Vehicle Using Bayesian Optimization and Soft Actor-Critic Algorithm. IEEE Trans. Transp. Electrif. 2024, 1-1. [Google Scholar] [CrossRef]
  31. Dufera, A.G.; Liu, T.; Xu, J. Regression models of Pearson correlation coefficient. Stat. Theory Relat. Fields 2023, 7, 97–106. [Google Scholar] [CrossRef]
  32. Yu, H.; Hutson, A.D. A robust Spearman correlation coefficient permutation test. Commun. Stat. -Theory Methods 2022, 53, 2141–2153. [Google Scholar] [CrossRef]
  33. Gao, D.; Zhou, Y.; Wang, T.; Wang, Y. A Method for Predicting the Remaining Useful Life of Lithium-Ion Batteries Based on Particle Filter Using Kendall Rank Correlation Coefficient. Energies 2020, 13, 4183. [Google Scholar] [CrossRef]
  34. Lu, Y.; Tang, X.; Wang, H. An Effective Stock Clustering Method Based on Hybrid Correlation Coefficient. In Proceedings of the 2018 IEEE 22nd International Conference on Computer Supported Cooperative Work in Design ((CSCWD)), Nanjing, China, 9–11 May 2018. [Google Scholar]
  35. Ding, R.; Zhang, S.; Chen, Y.; Rui, Z.; Hua, K.; Wu, Y.; Li, X.; Duan, X.; Wang, X.; Li, J.; et al. Application of Machine Learning in Optimizing Proton Exchange Membrane Fuel Cells: A Review. Energy AI 2022, 9, 100170. [Google Scholar] [CrossRef]
  36. Ibrahim, M.; Steiner, N.; Jemei, S.; Hissel, D. Wavelets-based approach for online Fuel Cells Remaining Useful lifetime Prediction. IEEE Trans. Ind. Electron. 2016, 63, 5057–5068. [Google Scholar] [CrossRef]
  37. Luo, J.; Chen, T.; Xiao, F.; Peng, Y. Remaining useful life prediction of PEMFC based on CNN-Birnn model. Int. J. Green Energy 2023, 20, 1729–1740. [Google Scholar] [CrossRef]
  38. Campbell, J.N.A.; Dais Ferreira, M.; Isenor, A.W. Generation of Vessel Track Characteristics Using a Conditional Generative Adversarial Network (CGAN). Appl. Artif. Intell. 2024, 38, 2360283. [Google Scholar] [CrossRef]
  39. Wang, F.-K.; Cheng, X.-B.; Hsiao, K.-C. Stacked long short-term memory model for proton exchange membrane fuel cell systems degradation. J. Power Sources 2020, 448, 227591. [Google Scholar] [CrossRef]
  40. Zhang, C.; Ji, C.; Hua, L.; Ma, H.; Nazir, M.S.; Peng, T. Evolutionary quantile regression gated recurrent unit network based on variational mode decomposition, improved whale optimization algorithm for probabilistic short-term wind speed prediction. Renew. Energy 2022, 197, 668–682. [Google Scholar] [CrossRef]
  41. Peng, Y.; Chen, T.; Xiao, F.; Zhang, S. Remaining useful lifetime prediction methods of proton exchange membrane fuel cell based on convolutional neural network-long short-term memory and convolutional neural network-bidirectional long short-term memory. Fuel Cells 2022, 23, 75–87. [Google Scholar] [CrossRef]
  42. Liu, C.; Shen, J.; Dong, Z.; He, Q.; Zhao, X. Accuracy improvement of fuel cell prognostics based on voltage prediction. Int. J. Hydrogen Energy 2024, 58, 839–851. [Google Scholar] [CrossRef]
  43. Yang, J.; Wang, L.; Zhang, B.; Zhang, H.; Wu, X.; Xu, X.; Deng, P.; Peng, Y. Remaining useful life prediction of vehicle-oriented PEMFC systems based on IGWO-BP neural network under real-world traffic conditions. Energy 2024, 291, 130334. [Google Scholar] [CrossRef]
  44. Tao, Z.; Zhang, C.; Xiong, J.; Hu, H.; Ji, J.; Peng, T.; Nazir, M.S. Evolutionary gate recurrent unit coupling convolutional neural network and improved manta ray foraging optimization algorithm for performance degradation prediction of PEMFC. Appl. Energy 2023, 336, 120821. [Google Scholar] [CrossRef]
  45. Ma, R.; Xie, R.; Xu, L.; Huangfu, Y.; Li, Y. A Hybrid Prognostic Method for PEMFC with Aging Parameter Prediction. IEEE Trans. Transp. Electrif. 2021, 7, 2318–2331. [Google Scholar] [CrossRef]
Figure 1. 2014 IEEE PHM Data Challenge durability test bench [25,26].
Figure 1. 2014 IEEE PHM Data Challenge durability test bench [25,26].
Materials 17 04713 g001
Figure 2. 60 kW fuel cell test bench [27].
Figure 2. 60 kW fuel cell test bench [27].
Materials 17 04713 g002
Figure 3. Schematic of the fuel cell system.
Figure 3. Schematic of the fuel cell system.
Materials 17 04713 g003
Figure 4. Matrix of correlation coefficients of feature parameters for steady-state condition: (a) Comparison of Spearman correlation coefficient; (b) Comparison of Kendall correlation coefficient.
Figure 4. Matrix of correlation coefficients of feature parameters for steady-state condition: (a) Comparison of Spearman correlation coefficient; (b) Comparison of Kendall correlation coefficient.
Materials 17 04713 g004
Figure 5. Comparison of DWT smoothing noise reduction results for feature parameters: (a) The original U t o t ; (b) Results of U t o t after smoothing noise reduction; (c) The original D o u t A I R ; (d) Results of D o u t A I R after smoothing noise reduction; (e) The original T i n H 2 ; (f) Results of T i n H 2 after smoothing noise reduction.
Figure 5. Comparison of DWT smoothing noise reduction results for feature parameters: (a) The original U t o t ; (b) Results of U t o t after smoothing noise reduction; (c) The original D o u t A I R ; (d) Results of D o u t A I R after smoothing noise reduction; (e) The original T i n H 2 ; (f) Results of T i n H 2 after smoothing noise reduction.
Materials 17 04713 g005aMaterials 17 04713 g005b
Figure 6. Matrix of correlation coefficients of feature parameters for dynamic cycling condition: (a) Comparison of Spearman correlation coefficient; (b) Comparison of Kendall correlation coefficient.
Figure 6. Matrix of correlation coefficients of feature parameters for dynamic cycling condition: (a) Comparison of Spearman correlation coefficient; (b) Comparison of Kendall correlation coefficient.
Materials 17 04713 g006
Figure 7. Comparison of reconstruction results of feature parameters based on CGAN: (a) The original U o u t ; (b) Results of U o u t after reconstruction; (c) The original U m i n ; (d) Results of U m i n after reconstruction; (e) The original P i n A I R ; (f) Results of P i n A I R after reconstruction.
Figure 7. Comparison of reconstruction results of feature parameters based on CGAN: (a) The original U o u t ; (b) Results of U o u t after reconstruction; (c) The original U m i n ; (d) Results of U m i n after reconstruction; (e) The original P i n A I R ; (f) Results of P i n A I R after reconstruction.
Materials 17 04713 g007
Figure 8. Schematic of the multi-input structure of the TCN model.
Figure 8. Schematic of the multi-input structure of the TCN model.
Materials 17 04713 g008
Figure 9. Schematic of the structure of multi-feature fusion life prediction model based on TCN-GRU.
Figure 9. Schematic of the structure of multi-feature fusion life prediction model based on TCN-GRU.
Materials 17 04713 g009
Figure 10. Results of dual-feature fusion life prediction for steady-state operating condition: (a) 40% of the training set ratio; (b) 60% of the training set ratio; (c) 80% of the training set ratio.
Figure 10. Results of dual-feature fusion life prediction for steady-state operating condition: (a) 40% of the training set ratio; (b) 60% of the training set ratio; (c) 80% of the training set ratio.
Materials 17 04713 g010
Figure 11. Results of three-feature fusion life prediction for steady-state operating condition: (a) 40% of the training set ratio; (b) 60% of the training set ratio; (c) 80% of the training set ratio.
Figure 11. Results of three-feature fusion life prediction for steady-state operating condition: (a) 40% of the training set ratio; (b) 60% of the training set ratio; (c) 80% of the training set ratio.
Materials 17 04713 g011
Figure 12. Results of dual-feature fusion life prediction for dynamic cycling condition: (a) 40% of the training set ratio; (b) 60% of the training set ratio; (c) 80% of the training set ratio.
Figure 12. Results of dual-feature fusion life prediction for dynamic cycling condition: (a) 40% of the training set ratio; (b) 60% of the training set ratio; (c) 80% of the training set ratio.
Materials 17 04713 g012aMaterials 17 04713 g012b
Figure 13. Results of three-feature fusion life prediction for dynamic cycling condition: (a) 40% of the training set ratio; (b) 60% of the training set ratio; (c) 80% of the training set ratio.
Figure 13. Results of three-feature fusion life prediction for dynamic cycling condition: (a) 40% of the training set ratio; (b) 60% of the training set ratio; (c) 80% of the training set ratio.
Materials 17 04713 g013aMaterials 17 04713 g013b
Figure 14. Evaluation of life prediction results for steady-state operating condition: (a) Comparison of RMSE results; (b) Comparison of R2 results.
Figure 14. Evaluation of life prediction results for steady-state operating condition: (a) Comparison of RMSE results; (b) Comparison of R2 results.
Materials 17 04713 g014
Figure 15. Evaluation of life prediction results for dynamic cycling condition: (a) Comparison of RMSE results; (b) Comparison of R2 results.
Figure 15. Evaluation of life prediction results for dynamic cycling condition: (a) Comparison of RMSE results; (b) Comparison of R2 results.
Materials 17 04713 g015
Table 1. 2014 IEEE PHM Data Challenge durability test stack parameters [25,26].
Table 1. 2014 IEEE PHM Data Challenge durability test stack parameters [25,26].
ParametersValueUnit
Active area100cm2
Operating temperature60°C
Pressure1.5bar
Relative humidity50%
Load current70A
Rated current density0.7A/cm2
Table 2. 60 kW fuel cell test bench stack parameters.
Table 2. 60 kW fuel cell test bench stack parameters.
ParametersValueUnit
Active area304cm2
Operating temperature80°C
Pressure2bar
Relative humidity80%
Table 3. Feature parameters for steady-state condition.
Table 3. Feature parameters for steady-state condition.
Feature ParametersSymbolUnit
Stack output voltage U t o t V
Inlet temperature of hydrogen T i n H 2 °C
Outlet temperature of hydrogen T o u t H 2 °C
Inlet temperature of air T i n A I R °C
Outlet temperature of air T o u t A I R °C
Inlet temperature of coolant T i n W A T °C
Outlet temperature of coolant T o u t W A T °C
Inlet pressure of hydrogen P i n H 2 Mbar
Outlet pressure of hydrogen P o u t H 2 Mbar
Inlet pressure of air P i n A I R Mbar
Outlet pressure of air P o u t A I R Mbar
Inlet flow rate of hydrogen D i n H 2 L/min
Outlet flow rate of hydrogen D o u t H 2 L/min
Inlet flow rate of air D i n A I R L/min
Outlet flow rate of air D o u t A I R L/min
Flow rate of coolant D W A T L/min
Air humidity H r A I R F C %
Table 4. Feature parameters for dynamic cycling condition.
Table 4. Feature parameters for dynamic cycling condition.
Feature ParametersSymbolUnit
Stack output voltage U o u t V
Stack output power P o u t W
Inlet temperature of hydrogen T i n H 2 °C
Inlet temperature of air T i n A I R °C
Outlet temperature of air T o u t A I R °C
Inlet temperature of coolant T i n W A T °C
Outlet temperature of coolant T o u t W A T °C
Temperature difference of coolant T W A T D i f °C
Inlet pressure of hydrogen P i n H 2 MPa
Inlet pressure of air P i n A I R MPa
Outlet pressure of air P o u t A I R MPa
Inlet pressure of coolant P i n W A T MPa
Outlet pressure of coolant P o u t W A T MPa
Pressure difference of cathode and anode P i n D i f MPa
Speed of air compressor n A C P S p d rpm
Minimum single voltage U m i n V
Table 5. Comparison of deep learning algorithms.
Table 5. Comparison of deep learning algorithms.
Name of AlgorithmAdvantageDrawbackUsage Scenario
Recurrent Neural Network (RNN)Wide range of applications with internal memory mechanism.Difficult to parallelize training; prone to gradient vanishing or gradient explosion; computationally inefficient.Short-to-medium time series
Convolutional Neural Network (CNN)It is good at capturing local features and can be computed in parallel with strong robustness and generalization ability.The number of parameters is large, and the convolution operation may lead to information loss.Spatial data
Long Short-Term Memory Network (LSTM)Use gating mechanism for easier convergence, stable training, and flexible structure.The model is large and the computational complexity is high and not easy to interpret.Memorizing long time series
Temporal Convolutional Network (TCN)Parallel computation capable; stable gradient propagation; good at capturing long-term dependencies.Input sequences need to be formulated with a fixed length and are sensitive to sequence ordering.Long time series
Gated Recurrent Unit (GRU)The structure is simple, with fewer gating mechanisms, less computational effort, and higher accuracy when the dataset is small.Weak memory, may not be able to capture long-term dependencies, weak processing of complex sequences.Simple long time series
Table 6. Model hyperparameter settings.
Table 6. Model hyperparameter settings.
ParameterValueParameterValue
Neurons30Bach size128
Filters32Learning rate0.001
Sizes3Dropout0.5
Dilations2Epoch30
Table 7. Evaluation of life prediction results for steady-state operating condition.
Table 7. Evaluation of life prediction results for steady-state operating condition.
Name of ModelRMSE (×10−3)R2
40%60%80%40%60%80%
RNN23.9113.2211.350.2770.4030.455
LSTM12.7710.5110.710.5410.6160.647
TCN-GRU (dual-feature)5.245.674.280.9080.9240.938
TCN-GRU (three-feature)5.044.423.270.9320.9430.965
Table 8. Evaluation of life prediction results for dynamic cycling condition.
Table 8. Evaluation of life prediction results for dynamic cycling condition.
Name of ModelRMSE (×10−3)R2
40%60%80%40%60%80%
CNN-LSTM5.606.953.660.6780.6850.789
TCN-GRU (dual-feature)5.304.232.940.8660.8810.887
TCN-GRU (three-feature)5.224.042.910.8850.8900.924
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, J.; Yan, F.; Du, C.; Zhang, Y.; Zheng, C.; Wang, J.; Chen, B. A Multi-Feature Fusion Method for Life Prediction of Automotive Proton Exchange Membrane Fuel Cell Based on TCN-GRU. Materials 2024, 17, 4713. https://doi.org/10.3390/ma17194713

AMA Style

Zhang J, Yan F, Du C, Zhang Y, Zheng C, Wang J, Chen B. A Multi-Feature Fusion Method for Life Prediction of Automotive Proton Exchange Membrane Fuel Cell Based on TCN-GRU. Materials. 2024; 17(19):4713. https://doi.org/10.3390/ma17194713

Chicago/Turabian Style

Zhang, Jiaming, Fuwu Yan, Changqing Du, Yiming Zhang, Chao Zheng, Jinhai Wang, and Ben Chen. 2024. "A Multi-Feature Fusion Method for Life Prediction of Automotive Proton Exchange Membrane Fuel Cell Based on TCN-GRU" Materials 17, no. 19: 4713. https://doi.org/10.3390/ma17194713

APA Style

Zhang, J., Yan, F., Du, C., Zhang, Y., Zheng, C., Wang, J., & Chen, B. (2024). A Multi-Feature Fusion Method for Life Prediction of Automotive Proton Exchange Membrane Fuel Cell Based on TCN-GRU. Materials, 17(19), 4713. https://doi.org/10.3390/ma17194713

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop