Abstract

Fuel cells could be a highly effective and eco-friendly technology to transform chemical energy stored in fuel to useful electricity and thus are presently appraised as a standout among the most encouraging advancements for future energy demand. Solid oxide fuel cells (SOFCs) have several advantages over other types of fuel cells, such as the flexibility of fuel used, high energy conversion, and relatively inexpensive catalysts due to high-temperature operation. The single chambers, wherein the anode and cathode are exposed to the same mixture of fuel, are promising for the portable power application due to the simplified, compact, sealing-free cell structure. The empirical regression models, such as artificial neural networks (ANNs), can be used as a black-box tool to simulate systems without solving the complicated physical equations merely by utilizing available experimental data. In this study, the performance of the newly proposed BSCF/GDC-based cathode SOFC was modeled using ANNs. The cell voltage was estimated with cathode preparation temperature, cell operating temperature, and cell current as input parameters by the one-layer feed-forward neural network. In order to acquire the appropriate model, several network structures were tested, and the network was trained by backpropagation algorithms. The data used during the training, validation, and test are the actual experimental results from our previous study. The optimum conditions to achieve maximum power of the cell were then determined by the genetic algorithm and the developed ANN.

1. Introduction

Fuel cells could be a highly effective and eco-friendly technology to transform chemical energy stored in fuel to useful electricity and thus are presently appraised as a standout among the most encouraging advancements for future energy demand [1]. Among different classes of fuel cells, the solid oxide fuel cell (SOFC) has displayed an extraordinary integration of benefits, including high energy efficiency, broad fuel versatility, significant contamination allowance, and low waste discharge [2, 3]. Currently, yttria-stabilized zirconia (YSZ), a zirconia-based electrolyte, is the most widely utilized electrolyte as it has sufficient oxide-ion conductivity and, furthermore, demonstrates suitable phase durability in both oxidizing and reducing environments [4]. However, in order to lessen the ohmic polarization loss, the zirconia-based electrolyte needs a very high operation temperature ∼1000°C, which provides adequate oxygen conductivity. Therefore, high temperature places rigorous constraints on its accompanying materials. Reducing the operation temperature to 500–800°C, classified as the intermediate temperature range, is supposed to elongate the lifespan of SOFC, minimize its degradation rate, and diversify the range of material selections [5].

On the contrary, reducing operation temperature requires a new pair of electrolyte-cathode, which is more ion conductive than the conventional pair of Y-stabilized ZrO2 and (La1−xSrx)MnO3. Therefore, the nickel cermet has been the unanimous selection for the anode composition. LAMOX is an oxide conductor family based initially on the parent La2Mo2O9 crystal. The lanthanide-doped LAMOX has the fundamental characteristics of an electrolyte at intermediate temperatures such as high ion conductivity [6, 7] and good phase stability [8]. However, the LAMOX composition is only occasionally investigated as the SOFC electrolyte.

On the cathode side, Ba0.5Sr0.5Co1−xFexO3−δ (BSCF) is usually suggested as a candidate in the intermediate temperature range, in perspective of its exceptional catalytic activity toward oxygen reduction combined with high ion conductivity [9], low electrochemical surface-exchange resistance [10], and reasonable fabrication cost [10]. However, the perovskite of BSCF does have its own drawback. The composition has the phase instability problem, similar to many other perovskites with cobalt on the B site [11]. It should be noted that the cathodic loss contributes the significant part of SOFC’s polarization loss. The large cathode polarization loss normally occurs from the fact that the mixed conductor of the cathode perovskite frequently shows the ionic conductivity significantly lower than its electronic conductivity [12, 13]. Hence, a composite cathode of ionic conductor and mixed conductor is adopted to minimize the cathodic loss.

The BSCF, gadolinium-doped ceria (GDC), and La1.8Dy0.2Mo2O9 (LDM) powders, therefore, were synthesized and characterized for fuel cell evaluation. The effect of different preparation conditions on the fuel cell behavior was found in this study. On the contrary, modeling and simulation of SOFCs are useful approaches to investigate effects of various design parameters and operating conditions on SOFC performances and also to help in SOFC improvements. The obtained results can be used in optimization, control, improvement, and design of the fuel cell system [1416]. Currently, there are plenty of models that have been published to contribute more knowledge of the fuel cell phenomena. The approaches for fuel cell modeling can be classified into two categories: theoretical and empirical [15, 17, 18]. In the theoretical mathematical approach, the spatial dimensions of the models vary from simplified zero (0-D) [19, 20] and one (1-D) [18, 2123] to more complex two (2-D) [2427] and three (3-D) [2830], with different characteristics and targets to different research purposes. The mathematical models, which are generally founded based on momentum, mass, and energy conservation laws, necessitate the information about many parameters and properties of the fuel cell system, the complicated equations, and the time-consuming numerical techniques. Therefore, using theoretical models makes the problem too complicated and impractical for some engineering applications.

Another approach, namely, empirical or data-driven modeling, might be more practical for the fuel cell user from the viewpoint the SOFC behavior can be inferred without deeply understanding the internal components [15, 17]. The empirical approaches can be applied to simulate the object’s behavior without an algorithmic solution, merely by utilizing available experimental data. These approaches include least-squares support vector machine (LS-SVM) [31] and the Hammerstein model [32, 33]. Among them, artificial neural network (ANN) presents some advantages such as high nonlinearity, matches experimental data with a low degree of error, and does not require any specific analytical equations and system descriptions. ANN is a black-box model, whereas numerical models allow for the examination of parameter distributions inside the cell. Numerical models may also be more accurate than ANN models if, for example, the operating conditions change significantly during dynamic operation. However, utilization of ANNs for modeling SOFCs appears a very promising way due to its rapid computation. Genetic algorithms (GAs) are a class of metaheuristic searching optimization techniques inspired by the mechanics of natural evolution and genetics [34]. Compared to conventional optimization methods, GAs are remarkably simple, robust, and able to handle the nondifferentiable, discontinuous, or multimodal problems. GAs have been successfully applied in a wide range of applications in various engineering, manufacturing, and management areas [35].

In this study, the performance of the BSCF/GDC-based cathode SOFC was modeled using ANNs. The cell voltage was estimated with cathode preparation temperature, cell operating temperature, and cell current as input parameters by the one-layer feed-forward neural network. In order to acquire the appropriate model, several network structures were tested, and the network was trained by backpropagation algorithms. The data used during the training, validation, and test have been obtained from the actual experimental results from our previous study [36]. The optimum condition to achieve maximum power was then determined by genetic algorithms and the developed ANN.

2. Experiments and Methodology

2.1. Experiments

The description of experimental configuration could be found in our study [36, 37]. The brief explanation of the cell is summarized as follows. The porous structure of the anode-supported cell was fabricated through the following steps. Commercial NiO (>99.5%), 1000°C calcined GDC, and 1000°C calcined LDM powders were mixed in a weight ratio of 4.8 : 3.2 : 0.8. Then, the mixture was mixed with 15 wt.% of graphite flash, uniaxial pressed, and then sintered for 2 h at 1175°C to form an anode with a size of ∼0.4 mm in thickness and ∼13.2 mm in diameter.

La1.8Dy0.2Mo2O9 (LDM) electrolyte was applied to the polished porous anode pellet by the spin-coating method with a desired thickness of ∼60 μm. As mentioned above, BFSC cannot be applied directly onto LDM due to chemical reaction. Therefore, it needs an additional barrier layer to prevent undesirable reaction.

A barrier layer of 0.5 wt.% Fe2O3-doped Gd0.1Ce0.9O1.95 (i-GDC) was applied on the LDM layer for one time to form NiO + GDC/LDM/i-GDC button. The button was then co-sintered at 1225°C for 5 hours with a heating and cooling rate of 2°C min−1.

To complete the final cell of NiO + GDC/LDM/i-GDC/BSCF + GDC, the surface of the i-GDC layer was screen-printed on by the blend of composite cathode and then sintered in air at 1000°C (cell a), 1015°C (cell b), 1025°C (cell c), or 1050°C (cell d) for 2 h. In this study, three single-cell sets were fabricated, all of which use a cathode gradient. The composite layer of 50 wt.% BSCF and 50 wt.% GDC was applied onto the i-GDC; meanwhile, the composite layer of 90 wt.% BSCF and 10 wt.% GDC layer was exposed to the air atmosphere.

The configuration of the fuel cell system of NiO + GDC/LDM/i-GDC/GDC-BSCF is shown in Figure 1(a). Two Au wires were connected to the single cell for current collecting, using the Keithley 238 source-measure unit. There is a coil arranged before the cell, which extends the heating time for the inflow. Both the coil and the cell were placed inside a quartz tube, together with a thermocouple which monitors the cell temperature. To activate the anode on-site, an extra Ni + GDC disk was attached above the anode. A sketch of the button cell was represented in Figure 1(b). The dimensions of the system are summarized in Figure 1(c).

The cell was operated in a mixture of methane/air with the flow rate of 650 sccm and the molar ratio of 1.5 : 1. The data were collected on the current density and the cell voltage of the fuel cell button for each combination of the sintering temperature of the cathode (1000, 1015, 1025, and 1050°C) and the cell operating temperatures (625, 650, 675, and 700°C).

2.2. Artificial Neural Network Models

Artificial neural network (ANN) is an information processing method inspired by the biological nervous systems. An ANN consists of a number of highly interconnected nodes, called neurons, which are able to receive and transmit the data [38, 39]. Feed-forward multilayer architectures, the simplest form of ANNs, are composed of an input layer, one or more hidden layers, and an output layer. The number of neurons in the input and output layers is equal to the number of inputs and outputs in the system to be modeled. Each neuron is linked to all neurons in the following layer by methods of weighted connections. The input of a neuron is the weighted aggregate of the data from its linked neurons. The response or output of a neuron is the transformation of its input using an appropriate activation function. The neurons will equally transmit their responses to all the neurons in the next layer. Input neurons do not activate their inputs; they just receive and send them to all the neurons in the first hidden layer. The number of neurons in the hidden layers is an important parameter which impacts the precision of the estimation. If too few neurons are used, the ANN might not have enough flexibility to satisfactorily approximate the data. However, if larger than the necessary number of neurons is included, the network then becomes overfitting, and it can memorize the data exactly but fail to generalize to a new situation [40]. In general, there is no rule for the number of neurons in the hidden layer, and it is usually determined by experiments. Figure 2 shows a feed-forward artificial neural network (3-5-1) structure with the input of operating temperature, sintering temperature, and current.

The activation function is used to calculate the output response of a neuron and can take any form, linear or nonlinear. The popular activation function, which is also employed in this study, is the logistic sigmoid [41]. The general form of the function is given as .

In order to enhance the performance of the network, the supplied input-output data should be normalized. In this study, the input data (xi) are scaled to the normalized value (xnorm) as follows:where xmax and xmin are the maximum and minimum of the experimental data, respectively.

The values of mean squared error (MSE) and coefficient of determination (R2) are frequently employed to evaluate the performance of the ANN. They were, respectively, calculated as follows [42]:where n is the total data points, ti is the predicted value from the networks, ai is the experimental response, and and are the average of predicted value and actual value, respectively.

There are several important factors affecting fuel cell performance, including cathode and anode structure, electrolyte material and thickness, cell temperature, and inlet and outlet gas compositions. In this study, two significant factors, namely, cathode sintering temperature and cell operating temperature, were examined. The range for sintering temperature is from 1000°C to 1050°C, and the range for the operating temperature is from 625°C to 700°C. As depicted in our earlier report, the sintering temperature affects the structure of the obtained cathode. The range for the investigated parameters is explained as follows. At the cathode sintering temperature below 1000°C, the low degree crystallization of BSCF minimizes the redox reaction, thus reducing the fuel cell performance. On the contrary, the phase compatibility between BSCF and GDC is ensured only up to 1050°C. Similarly, at the lower operating temperature of 600°C, there are not enough kinetic rate of ionic conducting, leading to poor performance of the fuel cell. However, at operating temperature over 700°C, the BSCF crystal becomes instable under the methane/air atmosphere, which limits the performance of the fuel cell.

Therefore, the network architecture is designed with three inputs, one hidden layer, and one output. The inputs are current density, sintering temperature of the cathode, and cell operating temperature. The network was trained using the backpropagation algorithm [43]. The maximum number of epochs and the minimum performance gradient are assigned to 400 and 10−5, respectively. Training stops when the maximum number of epochs is reached or when either the MSE or performance gradient is less than the predetermined goal. Through various trial experiments, the appropriate network structure for the current simulation is decided. The data sets were randomly divided into three subsets, namely, training, validation, and test, each of which contains 70%, 20%, and 10% number of samples, respectively. The validation and test sets are essential for the evaluation of the validation and modeling power of the networks.

The parameters of the model from the neural network were stored and used to optimize the power density using genetic algorithms in the next section.

2.3. Optimization by Genetic Algorithms

Genetic algorithms (GAs) are heuristic search algorithms inspired by the evolutionary ideas of natural selection and genetics. Genetic algorithms start with randomly generated chromosomes from the feasible region to produce a population. The fitness of each chromosome was defined by the objective function after converting the chromosomes to natural values. The algorithms iteratively modify a population of individual solutions analogous to processes happening in the nature—selection, reproduction, and mutation. The best fitness of the population comes toward an optimal solution by the principle of “survival of the fittest” [44]. The basic steps of genetic algorithms can be summarized as follows:(1)Randomly generating a primary population of N individuals(2)Producing the next generation by crossing over and mutation among individuals(3)Selecting the new population of N individuals from the parents and generation of (2)(4)Creating the next population by iterating steps (2) and (3) until the maximum number of generations is reached

The algorithm flow of the GA approach is shown in Figure 3.

The objective function in this study is the power density of the fuel cell, which is defined aswhere P is the power density (mW·cm−2), I is the current density (mA·cm−2) and is the input to the ANN model, and V is the cell voltage (V) which is estimated from the ANN model.

The design variables in the optimization problem are the inputs of the artificial neural network model. Their corresponding intervals were discussed in the previous section and are summarized as follows:(i)Sintering temperature of the cathode, [1000–1050] (°C)(ii)Operating temperature of the cell, [625–700] (°C)(iii)Electric current of the cell, [0–1500] (mA·cm−2)

The binary GA was used as the optimization method. A MATLAB™ program which was adapted from the literature [45] is used for the GA computations. The parameters of the GA such as population size and mutation probability values were set to be 100 and 0.10, respectively. The survival of the chromosomes was determined as the roulette wheel selection with elitism; it means that the best solution is always survive to the next population. The number of generations was assigned to be 500.

3. Results and Discussion

3.1. Microstructure of the Button Cell

In order to better understand the structure of the fabricated single cell, the single cell after I-V measurement was observed using the SEM method. The microstructure of the single cell in a cross-sectional image is shown in Figure 4. The SEM image shows an electrolyte thickness of ∼60 μm, total cathode thickness (including 2 layers) of ∼27 μm, and barrier layer i-GDC thickness of 6–8 μm. Moreover, the cross section shows well adherence between the cathode, i-GDC, electrolyte, and anode layers. Note that a crack of electrolyte at the left-hand side was formed by breaking out for the SEM measurement. Finally, the good adherence between layers and appropriate microstructures reveal that it is a good technique to fabricate the single cell for experimental and simulation study.

3.2. Artificial Neural Network Model Validation

An artificial neural network was used for modeling the cell voltage of the fuel cell. Experimental data obtained under different operating conditions were used to train and validate the neural network model. In ANN modeling, the number of hidden neurons is critical because excessive neurons will cause overfitting; however, too small number of neurons may not adequately capture the complexity of the system [46]. After performing some trials, the sufficient number of neurons in the hidden layer was 7. Consequently, a 3-7-1 feed-forward ANN architecture was employed in this study.

The results of the ANN modeling were described using current versus cell voltage graphs (or I-V plots) and were compared with data from experiments. The correlation coefficients (R2) for three subsets (training, validation, and testing) were 99.92%, 99.88%, and 99.87%, respectively. The MSE values for three subsets (training, validation, and testing) were 4.19 × 10−5, 7.14 × 10−5, and 7.04 × 10−5, respectively. The correlation coefficient and MSE of all data were 99.91% and 5.06 × 10−5, respectively.

Figure 5 shows the cell performances from the ANN modeling and the experimental data at the cathode sintering temperature of 1000°C as the operating temperature changes from 625°C to 700°C. The correlation and MSE for this subset are 99.95% and 2.86 × 10−5, respectively.

Figure 6 presents the cell voltage versus current density curves of the ANN predictions and the experimental data at 1015°C of cathode sintering temperature as the operating temperature is varied from 650°C to 700°C. The correlation coefficient and MSE of the model for this subset are 99.83% and 9.74 × 10−5, respectively.

In Figure 7, it can be seen that the cell voltage versus current density curves with various operating temperatures modeled by the ANN agrees well with the experimental data at the sintering temperature of 1025°C. The operating temperature ranges from 625°C to 700°C. The results show that the ANN model is in good agreement with the experimental data with a correlation coefficient of 99.89% and a MSE of 6.58 × 10−5.

Figure 8 shows I-V curves of the cells with the sintering temperature of 1050°C from the ANN and experimental data as the operating temperature changes from 650°C to 700°C. The correlation coefficient is 99.98%, and the MSE is 1.28 × 10−5 for this subset. The results show a good agreement between the model estimated and the experimental data.

From the results, it can be obtained that the ANN model succeeded in predicting the SOFC performance at an acceptable level of precision. It also confirms the capability of the ANN in investigating the effect of various parameters on its behavior.

3.3. Sensitivity Analysis

In order to identify the critical parameters and their degree of importance on the response outputs, a sensitivity analysis should be investigated. The analysis shows that the network output varies corresponding to the change of inputs and identifies the more sensitive parameters, which should be controlled precisely. The results of this analysis would also contribute deeper understanding about the model parameters for the design process.

In this study, the sensitivity analysis based on the weight matrix was chosen because of its simplicity. From the weight magnitude, several expressions have been proposed which have common properties: calculation of the product of the weights that connects the input neuron i and the hidden neuron j and that connects the hidden neuron j and the output neuron k for each of the hidden neurons of the network, which gives the sum of the calculated products. The following equation, reported by Garson [47], is widely used:in which is the sum of the connection weights between the N input neurons and the hidden neuron j. Qik denotes the effect of the input i on the output k, in relation to the other input variables.

The relative significance of each input variable on the output is presented in Figure 9. From the results, it can be obtained that the most important factor affecting the cell voltage is the sintering temperature of the BSCF cathode.

The importance of sintering temperature of the BSCF cathode was also reported in our previous study [36] or it can be seen in Figure 10. From the figure, even though the sintering temperature differs slightly, the cell power performance is significantly affected. Therefore, the sintering temperature should be precisely controlled and carefully decided in order to obtain the cell’s best performance.

3.4. Optimization Results

From the artificial neural network modeling results, genetic algorithm was utilized to search the optimum condition for the power density. The decision parameters were the sintering temperature of the cathode, operating temperature of the furnace, and current density. The boundary constraints for the design parameters were indicated earlier in Section 2.3. The fitness function is the power density.

Figure 11 depicts the fitness values (maximum and mean) of the population versus generation. As indicated in the figure, after about 20 generations, the value of fitness function attained a maximum value and then remained unchanged. After 100 generations, the maximum fuel cell power density of 451.64 mW/cm2 could be achieved at the sintering temperature of 1005°C, operating temperature of 668°C, and current density of 777 mA/cm2.

4. Conclusions

In this study, an artificial neural network was successfully applied to model the behavior of a BSCF-based single-chamber solid oxide fuel cell. The ANN model can adequately predict the SOFC performance without knowledge of numerous physical, chemical, and electrochemical factors and is validated with actual experimental data. The proposed network has a simple structure with one hidden layer and a small number of neurons. The prediction from the ANN model matched the experimental data with a low degree of error. From the trained ANN, the power density was estimated as fitness function, and the genetic optimization algorithm was employed to predict the optimal cell parameters such as sintering temperature of the cathode, operating temperature of the furnace, and current density of the cell. The maximum power density is 451.64 mW/cm2, which could be achieved at the sintering temperature of 1005°C, operating temperature of 668°C, and current density of 777 mA/cm2. The results suggested that the BSCF should be fabricated at 1005°C, and the button cell should be operated at 668°C.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

Part of this study (with different cathode materials-LSCF) was presented at the Third International Conference on Computational Science and Engineering in Ho Chi Minh City, Vietnam, November 2016.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank the Vietnam National Foundation for Science and Technology Development (NAFOSTED) for the financial support under grant number 104.03-2013.20. The authors would also like to thank Prof. Dah-Shyang Tsai (National Taiwan University of Science and Technology) for the kind support during the project.