Next Article in Journal
Static Voltage Stability Zoning Analysis Based on a Sensitivity Index Reflecting the Influence Degree of Photovoltaic Power Output on Voltage Stability
Next Article in Special Issue
Optimization and Control of New Power Systems under the Dual Carbon Goals: Key Issues, Advanced Techniques, and Perspectives
Previous Article in Journal
Advanced Optimisation and Forecasting Methods in Power Engineering—Introduction to the Special Issue
Previous Article in Special Issue
Cooperative Game-Based Collaborative Optimal Regulation-Assisted Digital Twins for Wide-Area Distributed Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intelligent Digital Twin Modelling for Hybrid PV-SOFC Power Generation System

1
State Grid Henan Electric Power Research Institute, Zhengzhou 450000, China
2
Anhui Jiyuan Software Co., Ltd., Hefei 230000, China
3
Faculty of Electric Power Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Submission received: 3 February 2023 / Revised: 14 March 2023 / Accepted: 15 March 2023 / Published: 17 March 2023
(This article belongs to the Special Issue Optimization and Control of New Energy Systems)

Abstract

:
Hydrogen (H2) energy is an ideal non-polluting renewable energy and can achieve long-term energy storage, which can effectively regulate the intermittence and seasonal fluctuation of solar energy. Solid oxide fuel cells (SOFC) can generate electricity from H2 with only outputs of water, waste heat, and almost no pollution. To solve the power generation instability and discontinuity of solar photovoltaic (PV) systems, a hybrid PV-SOFC power generation system has become one feasible solution. The “digital twin”, which integrates physical systems and information technology, offers a new view to deal with the current problems encountered during smart energy development. In particular, an accurate and reliable system model is the basis for achieving this vision. As core components, the reliable modelling of the PV cells and fuel cells (FCs) is crucial to the whole hybrid PV-SOFC power generation system’s optimal and reliable operation, which is based on the reliable identification of unknown model parameters. Hence, in this study, an artificial rabbits optimization (ARO)-based parameter identification strategy was proposed for the accurate modelling of PV cells and SOFCs, which was then validated on the PV double diode model (DDM) and SOFC electrochemical model under various operation scenarios. The simulation results demonstrated that ARO shows a more desirable performance in optimization accuracy and stability compared to other algorithms. For instance, the root mean square error (RMSE) obtained by ARO are 1.81% and 13.11% smaller than that obtained by ABC and WOA algorithms under the DDM of a PV cell. Meanwhile, for SOFC electrochemical model parameter identification under the 5 kW cell stack dataset, the RMSE obtained by ARO was only 2.72% and 4.88% to that of PSO for the (1 atm, 1173 K) and (3 atm, 1273 K) conditions, respectively. By establishing a digital twin model for PV cells and SOFCs, intelligent operation and management of both can be further achieved.

1. Introduction

With the rapid progress of society and the development of different industries, the increasing depletion of traditional energy resources and serious environmental deterioration are bound to promote the exploitation of new energy candidates. To effectively alleviate the pressure on the energy supply chain, changes in the energy structure are imperative, thus renewable and clean energy technologies attract worldwide attention. To further promote a low-carbon lifestyle and energy production structure in China and even the world, China has officially proposed a “dual carbon” target in 2020 [1]. In recent years, the application of renewable energy to cope with climate change and promote energy reform has become the general consensus and concerted action around the world [2]. Among them, solar energy has gained large-scale application thanks to its merits, e.g., abundant storage, clean, zero-emission, low cost, and easy installation. Thus, it is considered as one of the most mature renewable energies to be developed [3]. For PV plants, policy-based subsidy schemes or tax exemptions are crucial for their large-scale promotion, thus many countries have accordingly introduced relevant subsidy policies. For instance, in China, for distributed PV plants, a full electricity subsidy policy has been implemented. In detail, for distributed PV plants commissioned before 1 January 2018, the national subsidy is implemented in accordance with RMB 0.42/kWh, and for those commissioned after 1 January 2018, the national subsidy standard is implemented in accordance with RMB 0.37/kWh. In addition, distributed PV power generation systems are not required to pay the various funds and surcharges levied with the tariff, as well as system standby capacity charges and other related grid connection service fees for the self-consumption of electricity. For Germany, according to the Annual Tax Act 2022 approved by the German government, Germany will exempt rooftop PV from income tax and the corresponding value added tax from 2023 onwards if it meets the conditions. All the aforementioned initiatives are aimed at promoting the widespread use of PV power generation. However, considering that solar energy also has the disadvantages of power generation instability and discontinuity that are affected by weather, season, and other factors (e.g., night) [4], it is imperative to employ a complement system to photovoltaic (PV) system to enhance the reliability and continuity of electricity generation. As a green power technology of the 21st century, fuel cell (FC) can realize power generation without combustion [5]. FC can generate electricity from hydrogen (H2) with only the outputs of water, waste heat, and almost no pollution to the environment. At the same time, H2 energy is an ideal non-polluting renewable energy and can achieve long-term energy storage, which can regulate the intermittence and seasonal fluctuation of solar energy. Therefore, the H2 energy produced by the electrolysis of the water electrolyzer is stored in a hydrogen storage device [6], which provides a new energy storage method for PV systems.
Meanwhile, a tricky problem of the FC system power generation is to satisfy the demand for H2 supply. H2 is a clean and efficient energy carrier used for power generation, however, it is not a free primary energy in nature, thus it is crucial to seek a highly efficient and sustainable solution for its supply. Generally, H2 is produced through the electrolysis in the water electrolyzer, while its main problem is the quite high energy demand to support this process [6]. Thus, if the FC system can be combined with the PV system, the electrolyzer can be powered by the excess electricity produced by the PV system and H2 can be produced by water electrolysis. Therefore, a hybrid PV-FC power generation system consisting of the PV system, the electrolyzer, FCs, and hydrogen storage units has become a viable solution [7], which can effectively overcome the intermittent and stochastic nature of power generation from PV systems, as well as the instability of the output power from FCs [8]. In particular, solid oxide FC (SOFC) has been widely used in many areas including portable power supply, distributed power/thermal power system, high-performance power supply, and large-scale power station due to its high energy conversion efficiency and highly flexible fuel requirements [9]. Therefore, this study chooses SOFC as the research object of FC.
This hybrid PV-SOFC power generation system works on the principle that when there is sufficient light, the excess power produced by the PV system is supplied to water electrolysis to produce H2 in the electrolysis tank and stored in H2 storage units. When the power supply of the PV system is insufficient, the SOFCs use the stored H2 to generate electricity [10], enabling efficient use of energy and reducing energy wastage without any pollution. This hybrid system can be used not only as a stand-alone power generation system, but also as a hybrid DC microgrid system in modern power systems to optimize the power supply and energy flow patterns [11]. In general, such hybrid PV-SOFC power generation systems have a wide range of applications that can effectively alleviate the energy crisis and global environmental problems. However, there are many obstacles to the practical application of this hybrid PV-SOFC system, such as unstable energy conversion efficiency. In particular, as the core components of this hybrid PV-SOFC system, the precise modelling of SOFCs and PV cells is of the utmost significance since it is the basis to analyze and predict the output characteristics of the hybrid system and is the most fundamental step for subsequent performance analysis and state monitoring [12].
In recent years, the rapid advancement of the Internet of Things, 5G communications, artificial intelligence, and automation has provided technical support for the transformation of the energy industry and promoted the information and digital transformation of the energy system [13]. The “digital twin”, which integrates physical systems and information technology, provides new ideas and tools to deal with the current problems faced by the development of smart energy [14]. By establishing a digital twin model for PV cells and SOFCs, intelligent operation and management of both can be achieved. The “digital twin” was firstly used in the aviation sector, and then quickly spread to industrial manufacturing, smart cities, and infrastructure engineering. The core idea is to build a virtual model that can quickly, dynamically, and accurately simulate the whole life cycle of a physical entity based on its physical information and operational data, combined with the theoretical models accumulated in the research field [15]. The “digital twin” models of PV cell and SOFC are the twin of the physical objects of both, receiving data from the physical object in real-time and constantly updating the model to keep it in line with the physical object’s operational state. To realize the digital twin of the PV cell and SOFC models, the models that can reflect the physical reality need to be built, which is also the focus and difficulty of the current research since high reliability and accuracy are required regarding the models. Therefore, this work focused on the accurate and fast identification of unknown model parameters, with the aim of providing a feasible methodology for building a more accurate and reliable “digital twin” model for both.
For PV cells, until now, a variety of PV cell models have been proposed to describe their highly non-linear characteristics to accurately reflect their output current-voltage (I-V) and power-voltage (P-V) characteristics, in which double diode model (DDM) has received the most widespread application thanks to its most proper balance between model complexity and accuracy [16]. For SOFCs, the most representative SOFC model is the electrochemical model [17], which is also investigated in this study. Modelling accuracy depends on the model parameters, which means the reliable model parameter identification is the basis for realizing the “digital twin” of model. Regarding PV and SOFC modelling, the main difficulty is the identification of the unknown parameters of their models.

2. Literature Review

For the parameter identification of PV cells and SOFCs in a hybrid PV-SOFC power generation system, many methods have been developed. First, a systematic survey and analysis of previous research on PV cell parameter identification were conducted in literature [18], which covered and discussed a variety of PV cell modelling and parameter identification methods, which can be broadly classified into analytical, deterministic, and meta-heuristic algorithms, while the first two types of methods have their distinct drawbacks, such as low identification accuracy, high computational effort, strong model dependency, and extreme sensitivity to initial conditions and gradient information. Inspiringly, meta-heuristic algorithms can effectively compensate for the defects with the superiorities of strong flexibility, high accuracy, and insensitivity to gradient information, which is considered as the most ideal method to solve the parameter identification problem. Thus far, many meta-heuristic algorithms are utilized for PV cell parameter identification [16], in which a comprehensive review on meta-heuristic algorithms and related variants which have been applied to PV cell parameter identification was undertaken.
In detail, particle swarm optimization (PSO) was applied to extract the solar cell parameters from illuminated current-voltage characteristics for the single and double diode models. An artificial bee colony (ABC) algorithm was devised in reference [19] to accurately identify the solar cells’ parameters, which exhibited a better search capacity to face multi-modal objective functions. An improved whale optimization algorithm (WOA) was developed in reference [20], referred to as IWOA, to accurately extract the parameters of different PV models, in which two prey searching strategies were adopted to effectively balance the local exploitation and global exploration. For accurate modelling of PV module, an efficient parameter estimation technique named grey wolf optimization (GWO) [21] was developed to estimate the unknown parameters of the PV module. GWO showed high efficiency for parameter estimation of PV module under changing weather conditions. In reference [22], a fuzzy adaptive differential evolution algorithm-based method was developed, in which fuzzy selection strategy and adaptive parameter adjustment strategy were introduced to effectively control the crossover probability and mutation factors to avoid the discrimination into the local optimum while improving the convergence of the algorithm. Recently, an optimization algorithm called the musical chairs algorithm was devised to estimate the unknown parameters of PV cells [23]. The idea behind the use of this algorithm is to have a high number of search agents in the beginning to enhance the exploration and continuously reduce this number to enhance exploitation at the end of optimization and reduce the convergence time.
For FC parameter identification and modelling, there are two main categories [17], namely, traditional approaches and meta-heuristic algorithms. For the traditional methods, the fractional order derivative was adopted to analyze the general electrical characteristic of the dynamic model [24]. In the literature [25], the electrochemical impedance spectroscopy analysis was utilized for the electrochemical and thermodynamic information identification of SOFCs. The real-time estimation of SOFC model parameters in continuous time based on a well-studied FC experimental design was also realized in the literature [26].
For meta-heuristic algorithms, many approaches have been proposed and used to cope with the highly non-linear optimization problem of parameter identification for the SOFC. An improved genetic algorithm for a simplified SOFC electrochemical model was proposed in reference [27], which was based on an adapted encoding and control scheme. An improved Jingqiao adaptive differential evolution algorithm was also developed [28] to achieve fast and reliable SOFC parameter extraction using dynamic vector selection and adaptive cross-rate strategies. A bone particle swarm optimization algorithm was combined with a co-evolutionary framework for both high accuracy and robustness during parameter extraction [29]. In the literature [30], a simplified competitive swarm optimizer that used two simplified optimization strategies significantly improved the overall stability and accuracy. A modified version of grey wolf optimization for SOFC model identification was developed in reference [31], in which two scenarios based on temperature and pressure variations were utilized to show the system reliability. In addition, a Levenberg–Marquardt backpropagation algorithm-based parameter identification technique was proposed in literature [32,33] and applied to fast and precise identification of several unknown parameters for SOFC and PEMFC respectively, which showed satisfactory accuracy, speed, and stability. In reference [34], an extreme learning machine-based method was proposed to extract unknown parameters of SOFC model, in which extreme learning machine is applied to overcome two thorny obstacles (e.g., data shortage and noised data) via predicting additional data and updating noised data.
In this study, a parameter identification strategy based on artificial rabbits optimization (ARO) was designed for the accurate modeling of PV cells and SOFCs, whose contributions are outlined as follows:
  • A hybrid PV-SOFC system as an energy source is proposed to enhance energy utilization efficiency and power generation stability, in which the modelling of PV cells and SOFC is investigated since they are the most critical components in the hybrid system. The utilized ARO based parameter identification strategy can effectively improve both the exploration and exploitation ability, and dynamically adjust the proportion between the exploration and exploitation during iterations, which can effectively improve searching efficiency and avoid being trapped in local optimums;
  • For PV battery parameter identification, DDM benchmark PV battery model is used for verification, upon which the effectiveness of ARO for PV model parameter identification is verified. The simulation results indicate that ARO demonstrates higher accuracy and stability compared with other algorithms;
  • For the identification of unknown parameters of SOFC, the representative electrochemical model is used to validate the ARO based parameter identification method under different operating conditions and different datasets. The simulation results show that the ARO can obtain the minimum RMSE and the highest stability under various operating conditions.
The rest of this paper is organized as follows: Section 2 illustrates PV cell and SOFC modelling. ARO based parameter identification method is elaborated in Section 3. Case studies results are presented and analyzed in Section 4. At last, conclusions and perspectives are given in Section 5.

3. Hybrid PV-SOFC System Modelling

This section provides a detailed introduction to the mathematical modelling of the two most core components of the hybrid PV-SOFC system, i.e., the PV cell and SOFC, both of which are modelled based on their basic power generation principles. In detail, the DDM of the PV cell and the electrochemical model of SOFC are discussed, meanwhile, the objective functions for performance evaluation are also presented, respectively.

3.1. PV Cell Modelling

The equivalent model of DDM is illustrated in Figure 1. DDM consists of an ideal constant current source, two diodes D1 and D2, a series resistance Rs, and a shunt resistance Rsh. In particular, series resistance represents the total series resistance of material bulk resistance, thin layer resistance, and electrode contact resistance. Meanwhile, shunt resistance is mainly used to characterize the effect of leakage current on edge leakage. Compared with the single diode model (SDM), DDM has a higher precision and can describe and examine I-V and P-V characteristics in a wider range [16].
The output of DDM can be given by [20]:
I = I ph I d 1 I d 2 I sh
where    I sh and    I s indicate the current passing through the shunt resistance  R sh and series resistance  R s I ph denotes the photocurrent. In particular, the current    I d 1 and   I d 2 flowing through  D 1   and    D 2   and    I sh are expressed as:
I d 1 = I 01 [ exp ( q ( V + R s I ) a 1 N s V T ) 1 ]
I d 2 = I 02 [ exp ( q ( V + R s I ) a 2 N s V T ) 1 ]
I sh = V + I R s R sh
where  I 01 and    I 02 indicate the diode saturation current of  D 1   and    D 2 ,   a 1   and    a 2   are ideality factors of  D 1   and    D 2 .
In Equations (2) and (3),    V T is given by [35]:
V T = N s K T q
Hence, the output current of the PV cell based on DDM is presented as [16]:
I = I ph I 01 [ exp ( q ( V + I R s ) a 1 V T ) 1 ] I 02 [ exp ( q ( V + I R s ) a 2 V T ) 1 ] V + I R s R sh
According to Equation (6), in total of seven parameters need to be identified in DDM, namely,  I ph ,   I 01 ,   I 02 ,   R s ,   R sh a 1 , and  a 2 .
In order to effectively and quantitatively assess the performance of parameter identification methods, researchers have proposed a variety of evaluation criteria and objective functions. The computed results of the objective function allow for a comparison of the specific performance of different algorithms, thus effectively improving the shortcomings of one specific algorithm. Thus, the proper selection of the objective function is important during the practical application of parameter identification. The root mean square error (RMSE), which is most commonly applied in this problem [36], is selected as the objective function and the error is calculated using the actual and identified values as follows:
RMSE   ( x ) = 1 N k = 1 N ( f ( V L , I L , x ) ) 2
The error functions  f ( V L , I L , x ) for DDM is described as:
f DDM ( V L , I L , x ) = I ph I 01 [ exp ( q ( V L + I L R s ) a 1 k T pv ) 1 ] I 02 [ exp ( q ( V L + I L R s ) a 2 k T pv ) 1 ] V L + I L R s R sh I L
As can be seen from Equation (8), RMSE (x) should be minimized based on the optimization of the solution vector  x = { I ph ,   I 01 ,   I 02 ,   R s ,   R sh ,   a 1 ,   a 2 } .

3.2. SOFC Modelling

As shown in Figure 2, an SOFC mainly consists of the anode, the cathode, and the electrolyte between them [17]. The typical polarization curve of SOFC has three main distinct stages, namely, activation, concentration, and ohmic polarization, as demonstrated in Figure 3. The total chemical reaction process occurring in SOFC is expressed by the following equation.
At anode:
2 H 2 + 2 O 2 2 H 2 O + 4 e
At cathode:
O 2 + 4 e 2 O 2
In order to achieve simple and effective control of the SOFC, the electrochemical model is widely used, whose output voltage can be expressed as:
V c = N cell ( E o V act V ohm V con )
where  E o indicates the open circuit voltage;  N cell denotes the number of series cells;  V act V ohm , and  V con indicate activation voltage loss, ohmic voltage loss, and concentration voltage loss, respectively.
Based on Butler–Volmer equation, the  V act can be expressed as:
V act = A sin h 1 ( I load 2 I 0 , a ) + A sin h 1 ( I load 2 I 0 , c )
where  A indicates the slope of Tafel line;  I load I 0 , a , and  I 0 , c represent load current density, the exchange current density of anode and cathode, respectively.
Moreover,  V ohm can be described by:
V ohm = I load R ohm
where  R ohm denotes the ionic resistance.
Furthermore,  V con is described by:
V con = B ln ( 1 I load I L )
where  B indicates a constant;  I L represents the limiting current density.
From Equations (11)–(14), the output of SOFC is described by:
V c = N cell ( E o V act V ohm V con ) = N cell ( E o A sin h 1 ( I load 2 I 0 , a ) A sin h 1 ( I load 2 I 0 , c ) + B ln ( 1 I load I L ) I load R ohm )
From Equation (15), electrochemical model has seven unknown parameters (i.e.,  E o ,   A I 0 , a ,   I 0 , c B I L , and  R ohm ), where  I 0 , a > I 0 , c and  I L > I load .
Similar to PV cell, RMSE is defined as an objective function, calculated as:
RMSE ( x ) = 1 N k = 1 N ( error k ) 2 = 1 N k = 1 N ( V m , k V c , k ) 2
where x  = { E o ,   A I 0 , a ,   I 0 , c B I L , R ohm } represents the solution vector,  V m , k and  V c , k represent the voltage measured and the calculated voltage for the kth time, respectively.

4. Artificial Rabbits Optimization

In this section, the proposed background of the ARO and relevant mathematical models for optimization will be introduced.

4.1. Inspiration

The rabbit is a mammalian, herbivorous vertebrate. In the wild, the rabbit is highly stealthy. When it is stationary, its fur color blends in with the surrounding weeds, which also helps them to avoid enemies. Rabbits stay in low, dry bushes or burrows during the day and come out at night to search for food nearby.
ARO has its origins in the natural rabbit’s survival strategies, which mainly contains two aspects: food searching strategy and survival strategy [37]. First, rabbits do not eat the grass near the holes, as the grass grows tall and dense so that the rabbit is less likely to be discovered by enemies and thus will be well protected, so nesting grass is a natural camouflage for rabbits against enemies. Meanwhile, it is also used for food storage, as rabbits will feed on the grass when they cannot find food in the distance. In addition, rabbits have a very wide field of view which is basically used for overhead scanning, thus they can easily search for food in a large searching area. This foraging strategy will be considered as a global exploration during optimization.
Second, there are many natural predators of wild-bred rabbits, including foxes, wild dogs, eagles, and wolves. These natural predators are relatively large, and rabbits tend to hide burrows as a way of avoiding predator attacks when in danger. Rabbits have excellent running ability and are very agile, thus they can stop suddenly in mid-run, and turn sharply to escape pursuit. This survival strategy can effectively increase their chances of survival. This survival strategy will be regarded as the local exploitation during optimization.

4.2. Modelling of ARO

The food searching and hiding strategies of rabbits result in a switch between the two strategies that are adopted in ARO, and related mathematical modelling is demonstrated in this section.

4.2.1. Detour Food Searching Strategy

As mentioned above, although a rabbit is a herbivore, it does not eat the grass around its nest. This is because the grass around a rabbit’s nest is used for hiding. They only graze randomly in other areas, not in their own area, and we call this food searching behaviour “detour food searching”. In ARO, each rabbit has its own area with grass and burrows, and they tend to visit each other’s locations randomly for food searching. As a result, a mathematical model of rabbit meandering foraging is proposed as follows [37]:
v i ( t + 1 ) = x j ( t ) + R · ( x i ( t ) x j ( t ) ) + r o u n d ( 0.5 · ( 0.05 + r 1 ) ) · n 1 , i , j = 1 , ,   n   and   j i
R = L · c
L = ( e e ( t 1 T ) 2 ) · sin ( 2 π r 2 )
c ( k ) = { 1       if   k = g ( l ) 0                 else                                                 k = 1 ,   ,   d   and   l = 1 ,   ,   r 3 · d
g = r a n d p e r m ( d )
n 1 ~ N ( 0 , 1 )
where  v i ( t + 1 ) denotes the candidate position of the ith rabbit at the time t + 1,  x i ( t ) represents the position of the ith rabbit at the time t, L is the running length which represents the movement pace when performing the detour foraging, while other parameters can be referred to literature [37]. According to Equation (19), the running length L can generate a longer step during the initial iterations, while a shorter step will be generated during the later iterations. In Equation (18), R denotes the running operator which is employed to simulate the running characteristic of rabbits. Figure 4 shows the trajectory of 300 steps of the running operator R in 2-D and 3-D space [37]. From Figure 4, the searching space is probed by a series of sudden turns with the random direction, and with the increase of steps this probe length gradually gets shorter. Taking Figure 4 as an example, the horizontal coordinate x represents the number of iteration steps, and the vertical coordinate y represents the running step length, we can see that as the number of iteration steps x increases, the running length gradually decreases (i.e., the length of the fold lines in each sudden change in the upper right corner of the Figure 4 is obviously shorter than that of the fold lines in the lower left corner).
Equation (17) suggests that searching individuals perform a random search based on each individual’s position. This behaviour allows the rabbit to move to other rabbits’ areas. This particular food searching behaviour of rabbits, i.e., visiting others’ nests rather than their own, greatly facilitates global exploration and ensures the global searching ability.

4.2.2. Random Hiding Strategy

Rabbit burrowing is a habit evolved by rabbits in order to store food, avoid natural enemies, make a nest and settle down, and reproduce in the long-term harsh living environment in the wild. In ARO, a rabbit digs d burrow and selects one of all the burrows at random as a hiding place to hide in order to avoid being predated. Thus, the ith rabbit’s jth burrow is generated based on
b i , j ( t ) = x i ( t ) + H · g · x i ( t ) ,   i = 1 ,   ,   n   and   j = 1 ,   ,   d
H = T t + 1 T · r 4
n 2 ~ N ( 0 , 1 )
g ( k ) = { 1       if   k = = j 0                         else ,   , k = 1 ,   ,   d
where h indicates the hidden parameter that linearly in a randomly perturbed manner during the iterative process. Hence, during the initial stage, burrows are distributed in a larger neighbourhood of the rabbits. This neighbourhood is reduced with iterations increase. To avoid being chased and attacked by the predators, they are declined to randomly select a burrow from their burrows for sheltering, which can be modelled as:
v i ( t + 1 ) = x i ( t ) + R · ( r 4 · b i , r ( t ) x i ( t ) ) ,   i = 1 , ,   n
g r ( k ) = { 1       if   k = = r 5 · d 0                                               else ,   , k = 1 ,   ,   d
b i , r ( t ) = x i ( t ) + H · g r · x i ( t )
where  b i , r ( t ) denotes a randomly selected burrow, r4 and r5 represent two numbers randomly distributed between 0 and 1.

4.2.3. Energy Shrink Strategy

Rabbits tend to search for food in the initial stages of the iteration and often adopt random hiding during later stages of the iteration, which leads to the rabbit’s energy to diminish over time. Hence, an energy factor is applied to characterize the transition between global exploration and local exploitation, as follows:
A ( t ) = 4 ( 1 t T ) ln 1 r
where r denotes the random number between 0 and 1.
The larger the energy factor, the greater the amount of energy and stamina owned by the rabbit to search for food. The smaller value of the energy factor indicated that the rabbits were less energetic and therefore need to hide randomly. In specific, when A(t) > 1, rabbits tend to randomly explore different searching areas for food searching (global exploration); when A(t) ≤ 1, random hiding occurs (local exploitation), as demonstrated in Figure 5. To sum up, a group of random rabbits (candidate solutions) are generated in the searching space. Then, the rabbit will update its position during iterations, either selected randomly in the group or from its burrows. The energy factor A will decrease increase as iterations increase, upon which the detour food searching behavior will be switched to the hiding behavior.
The ARO algorithm has almost the same amount of performing both detour food searching and random hiding in the iterative process, which contribute significantly to balancing exploration and exploitation. Figure 6 shows the schematic diagram of probability calculation of the detour food searching [37].

5. Case Studies

In this section, ARO was applied in several PV cells and SOFC models for unknown parameter identification in the MATLAB/SIMULINK testing environment for performance validation. In particular, the DDM benchmark PV cell model was used for PV cell validation. For SOFC parameter identification, the electrochemical model was used for verification. All case studies are undertaken by Matlab 2022a through a personal computer with IntelR CoreTMi7 CPU at 2.0 GHz and 32 GB of RAM (Santa Clara, CA, USA).

5.1. PV Cell Model Parameter Identification

5.1.1. Design of ARO for PV Cell Parameter Identification

In order to achieve reliable optimization, the optimization variable  x j   is limited to the upper limit  x j max and lower limit  x j min as follows:
x j min x j x j max ,   j = 1 , 2 , J
where J represents the optimizing variables number.
If an individual rabbit violates Equation (31), it will be randomly reset by:
x j = x j min + r 2 ( x j max x j min )
where  r 2   indicates a number between 0 to 1.
The optimization framework for the identification of PV cell parameters based on ARO is shown in Figure 7. As can be seen in Figure 7, the measured I-V data is considered as input to ARO, and then the optimization procedure will be executed using ARO based on a specific PV cell model, and finally the identified parameters will be the output. For DDM parameter estimation, three meta-heuristic algorithms, namely ABC, WOA, and ARO, are applied to compare and validate the performance of ARO. In this section, the most representative DDM of PV cell is used, and the operating conditions are set to a standard operating environment (G = 1000 W/ m 2 and T = 33 °C). 26 sets of experimental I-V data are taken as input, which are acquired from a 57 mm diameter commercial silicon solar cell (R.T.C. France) operating under G = 1000 W/ m 2 and T = 33 °C [20,38], as detailed in Table 1. The maximum iteration number is designed to be 3000, which is found to be sufficient for convergence after several tests, while all methods were run independently for 30 times and the population size of the algorithms was designed to be 50 to obtain statistical results.

5.1.2. Validation on DDM

As tabulated in Table 2, the optimal parameter identification results and RMSE of three algorithms under DDM are demonstrated, which shows that ARO can obtain the smallest RMSE, followed by ABC and WOA, which indicates that ARO has the highest parameter identification accuracy. For example, the RMSE obtained by ARO are 1.81% and 13.11% smaller than that obtained by ABC and WOA algorithms. Moreover, by returning the extracted parameters to the DDM, the simulated data of ARO are in very good agreement with the actual data almost in all data points as shown in Table 2 and Figure 8.
The output I-V and P-V curves obtained by ARO are illustrated in Figure 8. The actual data are acquired from the experiment conducted in literature [38] published in 1986, where a 57 mm diameter commercial silicon solar cell (R.T.C. France) operating under G = 1000 W/m2 at T = 33 °C, with a CBM 8096 microcomputer acting as the controller. Figure 8 demonstrates that the simulated output curves acquired by ARO are basically the same as original actual curves, which can intuitively reflect its high fitting accuracy.
Convergence speed and stability are both important criteria for measuring the performance of an optimization method. The convergence curves of the RMSE of WOA, ABC, and ARO for DDM are demonstrated in Figure 9. Although WOA and ABC converge quickly in the beginning stage, on the other hand, they also easily stagnate quickly and suffer from prematurity, leading to low-quality local optimums. However, ARO tends to have a higher quality global optimum and a more stable convergence process due to its adaptive adjustment between global exploration and local exploitation throughout the iteration, thus avoiding falling into a local optimum.
Figure 10 shows the boxplots for the various algorithms under DDM, visualizing the distribution of RMSE results over 30 independent runs by different methods. ARO algorithm has the smallest upper and lower RMSE bound intervals of all the algorithms, with no outliers exist. This demonstrates that the ARO algorithm is able to achieve high convergence stability and high optimization accuracy for the problem of DDM parameter identification.

5.2. SOFC Model Parameter Identification

In this section, ARO was adopted for the complex optimization problem of SOFC parameter identification, which is typically highly non-linear with multiple constraints. This section evaluated the performance of ARO, GWO, and PSO based on two different V-I datasets under the electrochemical model, respectively. One set of data was collected from a battery stack consisting of 79 cells, as shown in Figure 11.
At the same time, another dataset was selected from a 5 kW cell stack in WCS-SOFC MATLAB/SIMULINK consisting of 96 cells operating under various operating scenarios [39,40], as shown in Figure 12. Table 3 and Table 4 show the dataset size and the parameters range, respectively.

5.2.1. Validation of ARO

This subsection uses the first aforementioned dataset, and the seven unknown parameters of the electrochemical model are set to [Eo, A, Rohm, B, I0,a, I0,c, IL] = [1.15, 0.02, 0.0004, 0.03, 12, 4, 152].
The identification parameters results are illustrated in Table 5 along with the RMSE, in which the RMSE values can give a direct and clear indication of the optimization accuracy. As shown in Table 5, ARO identifies the unknown parameters in a largest similarity to the real measurement data. ARO also achieves the smallest RMSE compared to other algorithms, e.g., the RMSE value achieved by ARO is only 8.67% and 20.23% to that of GWO and PSO. It verifies ARO achieves parameter estimation of the electrochemical model with the highest accuracy.
Meanwhile, Figure 13 demonstrates the boxplot and convergence graphs of various methods under 100 independent runs. Figure 13a illustrates the RMSE results distribution of different algorithms, it can be verified that ARO owns the smallest distribution range with minimal lower and upper bounds. Figure 13b also proves that ARO has a more stable convergence process and higher convergence precision.

5.2.2. Verification on a 5 kW Cell Stack

The effectiveness of ARO was verified on a 5 kW cell stack here. Furthermore, the estimated results under four various operating scenarios are shown in Table 6, which demonstrates that the identification results of ARO can obtain the smallest RMSE under all operating scenarios compared with GWO and PSO. For example, the RMSE obtained by ARO is only 2.72% and 4.88% to that of PSO for the (1 atm, 1173 K) and (3 atm, 1273 K) conditions, respectively.
In addition, Figure 14 shows the boxplots for the various algorithms under different running conditions, which illustrates the RMSE distribution for each algorithm over 20 runs. It is clear that the ARO achieved the smallest range of RMSE distributions under various conditions, and also yielded the smallest lower and upper bounds on RMSE. This validates that ARO can effectively boost convergence stability and identification accuracy.
In addition, Figure 15 depicts convergence curves of different algorithms under various operating scenarios, in which ARO can converge to the minimum RMSE value under all four testing scenarios. In contrast, although PSO also converges quickly in the beginning stage, it stagnates quickly and suffers from prematurity. GWO can converge continuously throughout the whole evolutionary process, while it struggles to obtain high-quality optimal solutions in all scenarios. The convergence result indicated that ARO is able to break the adsorption of local minima and find a more desirable searching direction. Namely, it shows a more stable convergence tendency and achieves higher quality solutions in all scenarios due to its strong exploration mechanism.

6. Conclusions and Perspectives

This paper took the two core components of a hybrid PV-SOFC system, i.e., the PV cell and the SOFC as the research focus with the aim of reliably identifying the unknown parameters in their models and thus achieving accurate modelling, the main contributions of which are outlined as:
  • To deal with parameters identification of the PV cell and SOFC models, this paper applied an ARO based intelligent parameter identification strategy to achieve accurate modelling of both models. The aim was to provide a feasible method for building a more accurate and reliable “digital twin” model of this hybrid power generation system, and ultimately to achieve a highly automated, intelligent, and low-carbon system operation and management system;
  • The proposed ARO-based parameter identification strategy improves the convergence speed and optimization accuracy by designing a dynamic searching mechanism to regulate the searching behaviour during the iterations, i.e., detour foraging and random hiding, so as to better balance global exploration and local exploitation;
  • ARO was applied for parameters estimation for the DDM model of PV cells and the electrochemical model of SOFC, and its effectiveness was fully verified. The simulation results indicated that ARO show higher accuracy and stability in comparison with other algorithms. For instance, for parameter identification of SOFC electrochemical model under the first dataset, the RMSE value achieved by ARO was only 8.67% and 20.23% to that of GWO and PSO. Meanwhile, the convergence curves of the RMSE obtained by different algorithms also proved that ARO acquires a higher quality global optimum and shows a more stable convergence process.
Overall, this study also has some limitations, for instance, the limited diversity of PV cell models and SOFC models are utilized for verification. In addition, all the case studies were only tested under the ideal laboratory simulation environment rather than practical operation conditions. Meanwhile, only the core components of the hybrid PV-SOFC system were studied instead of all components. Therefore, follow-up research is necessary, and future research will focus on the following three areas:
  • For PV cells and SOFC parameter identification, more types of cell models can be applied for validation. Meanwhile, the ability to perform parameter identification online and in real time can be further developed;
  • Future validation of the SOFC models considering degradation mechanisms can bring higher engineering value in terms of health condition monitoring, and fault detection and diagnosis;
  • The research focus of the digital twin includes the creation of digital twin models, physical information fusion, and service applications. This paper focused on the study of building accurate and reliable digital twin models of core components in the hybrid PV-SOFC system. Therefore, further research is required for the digital twin model of the whole system in its entirety and for the application of the methods proposed in this paper in the areas of optimal design, optimal operation, and fault diagnosis.

Author Contributions

Z.G.: writing the original draft and editing. Z.Y., P.N. and C.C.: conceptualization. X.W., J.Z. and X.H.: visualization and contributed to the discussion of the topic. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the science and technology project from State Grid Corporation of China: Research on Digital Twin-based Management and Interaction Technology for Efficient Collaboration of County Energy Internet, grant number (5700-202224202A-1-1-ZN).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Abbreviation
ABCartificial bee colony
AROartificial rabbits optimization
DDMdouble diode mode
FCfuel cell
GWOgrey wolf optimization
H2hydrogen
I-Vcurrent-voltage
PSOparticle swarm optimization
PVphotovoltaic
P-Vpower-voltage
RMSEroot mean square error
WOAwhale optimization algorithm
SDMsingle diode model
SOFCsolid oxide fuel cell
Variables
  a 1 the ideality factors of  D 1
  a 2 the ideality factors of  D 2
  q = 1.6 × 10 19 C the electron charge
  N s the number of series connected PV cells in the PV panel
Vtthe junction thermal voltage
  A the slope of Tafel line
  E o the open circuit voltage
  N cell the number of series cells
  I 01 the diode saturation current of  D 1
  I 02 the diode saturation current of  D 2
  I 0 , a the exchange current density of anode
  I 0 , c the exchange current density of cathode
  I d 1 the current flowing through  D 1
  I d 2 the current flowing through  D 2
  I load the load current density
  I L the limiting current density
  I ph the photocurrent
  I s the current passing through the series resistance
  I sh the current passing through the shunt resistance
  V act the activation voltage loss
  V c , k the voltage calculated voltage for the  k th time
  V con the concentration voltage loss
  V m , k the voltage measured for the kth time
  V ohm the ohmic voltage loss

References

  1. Huang, R.; Zhang, S.; Wang, P. Key areas and pathways for carbon emissions reduction in Beijing for the “Dual Carbon” targets. Energy Policy 2022, 164, 112873. [Google Scholar] [CrossRef]
  2. Zhang, G.; Shi, Y.; Maleki, A.; Rosen, M.A. Optimal location and size of a grid-independent solar/hydrogen system for rural areas using an efficient heuristic approach. Renew. Energy 2020, 156, 1203–1214. [Google Scholar] [CrossRef]
  3. Verbois, H.; Huva, R.; Rusydi, A.; Walsh, W. Solar irradiance forecasting in the tropics using numerical weather prediction and statistical learning. Sol. Energy 2018, 162, 265–277. [Google Scholar] [CrossRef]
  4. Awasthi, A.; Shukla, A.K.; SR, M.M.; Dondariya, C.; Shukla, K.N.; Porwal, D.; Richhariya, G. Review on sun tracking technology in solar PV system. Energy Rep. 2020, 6, 392–405. [Google Scholar] [CrossRef]
  5. Ahmad, M.Z.; Ahmad, S.H.; Chen, R.S.; Ismail, A.F.; Hazan, R.; Baharuddin, N.A. Review on recent advancement in cathode material for lower and intermediate temperature solid oxide fuel cells application. Int. J. Hydrog. Energy 2022, 47, 1103–1120. [Google Scholar] [CrossRef]
  6. Yang, B.; Wang, J.; Yu, L.; Shu, H.; Yu, T.; Zhang, X.; Sun, L. A critical survey on proton exchange membrane fuel cell parameter estimation using meta-heuristic algorithms. J. Clean. Prod. 2020, 265, 121660. [Google Scholar] [CrossRef]
  7. Fathy, A.; Yousri, D.; Alanazi, T.; Rezk, H. Minimum hydrogen consumption-based control strategy of fuel cell/PV/battery/supercapacitor hybrid system using recent approach-based parasitism-predation algorithm. Energy 2021, 225, 120316. [Google Scholar] [CrossRef]
  8. Ghaffari, A.; Askarzadeh, A. Design optimization of a hybrid system subject to reliability level and renewable energy penetration. Energy 2020, 193, 116754. [Google Scholar] [CrossRef]
  9. Singh, M.; Zappa, D.; Comini, E. Solid oxide fuel cell: Decade of progress, future perspectives and challenges. Int. J. Hydrog. Energy 2021, 46, 27643–27674. [Google Scholar] [CrossRef]
  10. Melzi, B.; Kefif, N.; El Haj Assad, M. Modelling and optimal design of hybrid power system photovoltaic/solid oxide fuel cell for a mediterranean City. Energy Eng. 2021, 118, 1767–1781. [Google Scholar] [CrossRef]
  11. Siti Khodijah, A.; Shiba, M.; Obara, S.Y. Design of compensation battery for tidal power-photovoltaics-SOFC microgrids in Ternate and Pulau-Tidore Islands. Int. J. Energy Res. 2021, 45, 2100–2115. [Google Scholar] [CrossRef]
  12. Wang, J.; Yang, B.; Li, D.; Zeng, C.; Chen, Y.; Guo, Z.; Yu, T. Photovoltaic cell parameter estimation based on improved equilibrium optimizer algorithm. Energy Convers. Manag. 2021, 236, 114051. [Google Scholar] [CrossRef]
  13. He, X.; Ai, Q.; Qiu, R.C.; Huang, W.; Piao, L.; Liu, H. A big data architecture design for smart grids based on random matrix theory. IEEE Trans. Smart Grid 2015, 8, 674–686. [Google Scholar] [CrossRef] [Green Version]
  14. Grieves, M. Digital twin: Manufacturing excellence through virtual factory replication. White Pap. 2017. [Google Scholar] [CrossRef]
  15. Singh, M.; Srivastava, R.; Fuenmayor, E.; Kuts, V.; Qiao, Y.; Murray, N.; Devine, D. Applications of Digital Twin across industries: A review. Appl. Sci. 2022, 12, 5727. [Google Scholar] [CrossRef]
  16. Yang, B.; Wang, J.; Zhang, X.; Yu, T.; Yao, W.; Shu, H.; Sun, L. Comprehensive overview of meta-heuristic algorithm applications on PV cell parameter identification. Energy Convers. Manag. 2020, 208, 112595. [Google Scholar] [CrossRef]
  17. Yang, B.; Wang, J.; Zhang, M.; Shu, H.; Yu, T.; Zhang, X.; Sun, L. A state-of-the-art survey of solid oxide fuel cell parameter identification: Modelling, methodology, and perspectives. Energy Convers. Manag. 2020, 213, 112856. [Google Scholar] [CrossRef]
  18. Abbassi, R.; Abbassi, A.; Jemli, M.; Chebbi, S. Identification of unknown parameters of solar cell models: A comprehensive overview of available approaches. Renew. Sustain. Energy Rev. 2018, 90, 453–474. [Google Scholar] [CrossRef]
  19. Oliva, D.; Cuevas, E.; Pajares, G. Parameter identification of solar cells using artificial bee colony optimization. Energy 2014, 72, 93–102. [Google Scholar] [CrossRef]
  20. Xiong, G.; Zhang, J.; Shi, D.; He, Y. Parameter extraction of solar photovoltaic models using an improved whale optimization algorithm. Energy Convers. Manag. 2018, 174, 388–405. [Google Scholar] [CrossRef]
  21. Nayak, B.; Mohapatra, A.; Mohanty, K.B. Parameter estimation of single diode PV module based on GWO algorithm. Renew. Energy Focus 2019, 30, 1–12. [Google Scholar] [CrossRef]
  22. Dang, J.; Wang, G.; Xia, C.; Jia, R.; Li, P. Research on the parameter identification of PV module based on fuzzy adaptive differential evolution algorithm. Energy Rep. 2022, 8, 12081–12091. [Google Scholar] [CrossRef]
  23. Eltamaly, A.M. Musical chairs algorithm for parameters estimation of PV cells. Sol. Energy 2022, 241, 601–620. [Google Scholar] [CrossRef]
  24. Cao, H.; Deng, Z.; Li, X.; Yang, J.; Qin, Y. Dynamic modeling of electrical characteristics of solid oxide fuel cells using fractional derivatives. Int. J. Hydrog. Energy 2010, 35, 1749–1758. [Google Scholar] [CrossRef]
  25. Caliandro, P.; Nakajo, A.; Diethelm, S. Model-assisted identification of solid oxide cell elementary processes by electrochemical impedance spectroscopy measurements. J. Power Sources 2019, 436, 226838. [Google Scholar] [CrossRef]
  26. Jayasankar, B.R.; Huang, B.; Ben-Zvi, A. Receding horizon experiment design with application in sofc parameter estimation. IFAC Proc. Vol. 2010, 43, 541–546. [Google Scholar] [CrossRef] [Green Version]
  27. Yang, J.; Li, X.; Jiang, J.H.; Jian, L.; Zhao, L.; Jiang, J.G.; Xu, L.H. Parameter optimization for tubular solid oxide fuel cell stack based on the dynamic model and an improved genetic algorithm. Int. J. Hydrog. Energy 2011, 36, 6160–6174. [Google Scholar] [CrossRef]
  28. Gong, W.; Cai, Z.; Yang, J.; Li, X.; Jian, L. Parameter identification of an SOFC model with an efficient, adaptive differential evolution algorithm. Int. J. Hydrog. Energy 2014, 39, 5083–5096. [Google Scholar] [CrossRef]
  29. Jiang, B.; Wang, N.; Wang, L. Parameter identification for solid oxide fuel cells using cooperative barebone particle swarm optimization with hybrid learning. Int. J. Hydrog. Energy 2014, 39, 532–542. [Google Scholar] [CrossRef]
  30. Wei, Y.; Stanford, R.J. Parameter identification of solid oxide fuel cell by Chaotic Binary Shark Smell Optimization method. Energy 2019, 188, 115770. [Google Scholar] [CrossRef]
  31. Yang, B.; Chen, Y.; Guo, Z.; Wang, J.; Zeng, C.; Li, D.; Zhang, X. Levenberg-Marquardt backpropagation algorithm for parameter identification of solid oxide fuel cells. Int. J. Energy Res. 2021, 45, 17903–17923. [Google Scholar] [CrossRef]
  32. Yang, B.; Zeng, C.; Wang, L.; Guo, Y.; Chen, G.; Guo, Z.; Zhu, J. Parameter identification of proton exchange membrane fuel cell via Levenberg-Marquardt backpropagation algorithm. Int. J. Hydrog. Energy 2021, 46, 22998–23012. [Google Scholar] [CrossRef]
  33. Yang, B.; Guo, Z.; Yang, Y.; Chen, Y.; Zhang, R.; Su, K.; Zhang, X. Extreme learning machine based meta-heuristic algorithms for parameter extraction of solid oxide fuel cells. Appl. Energy 2021, 303, 117630. [Google Scholar] [CrossRef]
  34. Wang, J.; Xu, Y.P.; She, C.; Xu, P.; Bagal, H.A. Optimal parameter identification of SOFC model using modified gray wolf optimization algorithm. Energy 2022, 240, 122800. [Google Scholar] [CrossRef]
  35. Jordehi, A.R. Parameter estimation of solar photovoltaic (PV) cells: A review. Renew. Sustain. Energy Rev. 2016, 61, 354–371. [Google Scholar] [CrossRef]
  36. Muangkote, N.; Sunat, K.; Chiewchanwattana, S.; Kaiwinit, S. An advanced onlooker-ranking-based adaptive differential evolution to extract the parameters of solar cell models. Renew. Energy 2019, 134, 1129–1147. [Google Scholar] [CrossRef]
  37. Wang, L.; Cao, Q.; Zhang, Z.; Mirjalili, S.; Zhao, W. Artificial rabbits optimization: A new bio-inspired meta-heuristic algorithm for solving engineering optimization problems. Eng. Appl. Artif. Intell. 2022, 114, 105082. [Google Scholar] [CrossRef]
  38. Easwarakhanthan, T.; Bottin, J.; Bouhouch, I.; Boutrit, C. Nonlinear minimization algorithm for determining the solar cell parameters with microcomputers. Int. J. Sol. Energy 1986, 4, 1–12. [Google Scholar] [CrossRef]
  39. Hagan, M.T.; Demuth, H.B.; Beale, M.H. Neural Network Design; PWS Publishing Co.: Boston, MA, USA, 2002. [Google Scholar]
  40. Yang, B.; Li, D.Y.; Zeng, C.Y.; Chen, Y.J.; Guo, Z.X.; Wang, J.B.; Shu, H.C.; Yu, T.; Zhu, J.W. Parameter extraction of PEMFC via Bayesian regularization neural network based meta-heuristic algorithms. Energy 2021, 228, 120592. [Google Scholar] [CrossRef]
Figure 1. Double diode model.
Figure 1. Double diode model.
Energies 16 02806 g001
Figure 2. Electrochemical operation mechanism of SOFC [17].
Figure 2. Electrochemical operation mechanism of SOFC [17].
Energies 16 02806 g002
Figure 3. Polarization curve of SOFC.
Figure 3. Polarization curve of SOFC.
Energies 16 02806 g003
Figure 4. Trajectory of 300 steps of the running operator [37].
Figure 4. Trajectory of 300 steps of the running operator [37].
Energies 16 02806 g004
Figure 5. Searching behavior adjustment based on energy factor A.
Figure 5. Searching behavior adjustment based on energy factor A.
Energies 16 02806 g005
Figure 6. Schematic diagram of probability calculation of detour food searching.
Figure 6. Schematic diagram of probability calculation of detour food searching.
Energies 16 02806 g006
Figure 7. ARO based parameter identification framework.
Figure 7. ARO based parameter identification framework.
Energies 16 02806 g007
Figure 8. Output I-V and P-V curves obtained by ARO: (a) I-V curve; and (b) P-V curve.
Figure 8. Output I-V and P-V curves obtained by ARO: (a) I-V curve; and (b) P-V curve.
Energies 16 02806 g008
Figure 9. Convergence curves obtained by various algorithms for DDM.
Figure 9. Convergence curves obtained by various algorithms for DDM.
Energies 16 02806 g009
Figure 10. Boxplot graph obtained by various algorithms.
Figure 10. Boxplot graph obtained by various algorithms.
Energies 16 02806 g010
Figure 11. Datasets for validation for electrochemical model.
Figure 11. Datasets for validation for electrochemical model.
Energies 16 02806 g011
Figure 12. Datasets for 5 kW cell stack under various operating scenarios.
Figure 12. Datasets for 5 kW cell stack under various operating scenarios.
Energies 16 02806 g012
Figure 13. Boxplot and convergence graphs of various algorithms: (a) boxplot; and (b) convergence graph.
Figure 13. Boxplot and convergence graphs of various algorithms: (a) boxplot; and (b) convergence graph.
Energies 16 02806 g013
Figure 14. Boxplot of different algorithms at various environmental conditions: (a)1 atm, 1173 K; (b) 1 atm, 1273 K; (c) 3 atm, 1173 K; and (d) 3 atm, 1273 K.
Figure 14. Boxplot of different algorithms at various environmental conditions: (a)1 atm, 1173 K; (b) 1 atm, 1273 K; (c) 3 atm, 1173 K; and (d) 3 atm, 1273 K.
Energies 16 02806 g014
Figure 15. Convergence graphs of different algorithms at various environmental conditions: (a) 1 atm, 1173 K; (b) 1 atm, 1273 K; (c) 3a tm, 1173 K; and (d) 3 atm, 1273 K.
Figure 15. Convergence graphs of different algorithms at various environmental conditions: (a) 1 atm, 1173 K; (b) 1 atm, 1273 K; (c) 3a tm, 1173 K; and (d) 3 atm, 1273 K.
Energies 16 02806 g015
Table 1. Benchmark I-V dataset.
Table 1. Benchmark I-V dataset.
I-V Data12345678910111213
VL
(V)
−0.2057−0.1291−0.05880.00570.06460.11850.16780.21320.25450.29240.32690.35850.3873
IL
(A)
0.76400.76200.76050.76050.76000.75900.75700.75700.75550.75400.75050.74650.7385
I-V Data14151617181920212223242526
VL
(V)
0.41370.43730.45900.47840.49600.51190.52650.53980.55210.56330.57360.58330.5900
IL
(A)
0.72800.70650.67550.63200.57300.49900.41300.31650.21200.1035−0.0100−0.1230−0.2100
Table 2. Parameter identification results obtained by various algorithms under DDM.
Table 2. Parameter identification results obtained by various algorithms under DDM.
AlgorithmIph(A)I01(µA)Rs(Ω)Rsh(Ω)a1I02(µA)a2RMSE
ABC0.76040.54500.037252.09781.81040.15111.41961.0037 × 10−3
WOA0.76030.53330.035871.71161.69210.15021.43601.1342 × 10−3
ARO0.76080.22560.036654.62192.00000.21621.45169.8554 × 10−4
Table 3. The size of dataset.
Table 3. The size of dataset.
Datasets for ValidationDatasets for 5 kW Cell Stack
Electrochemical model(1 atm, 1173 K)(1 atm, 1273 K)(3 atm, 1173 K)(3 atm, 1273 K)
3469696565
Table 4. The range of unknown parameters.
Table 4. The range of unknown parameters.
Electrochemical Model
ParametersEo (V)A (V)Rohm (kΩ·cm2)B (V)I0,a (mA/cm2)I0,c (mA/cm2)IL (mA/cm2)
Lower limit0000000
Upper limit1.21113030200
Table 5. Identification results of various algorithms.
Table 5. Identification results of various algorithms.
AlgorithmsIdentified ParametersRMSE
E0ARBI0,aI0,cIL
Actual value1.15000.02000.00040.030012.00004.0000152.00000
GWO1.14790.02980.00000.045320.25586.2900157.49834.9882 × 10−4
PSO1.14950.03160.00000.0457305.4398158.08962.1379 × 10−4
ARO1.14980.02840.00020.0338304.7907153.21754.3264 × 10−5
Table 6. Identification results of various algorithms for 5 kW cell stack dataset.
Table 6. Identification results of various algorithms for 5 kW cell stack dataset.
ConditionsAlgorithmIdentified ParametersRMSE
EoARohmBI0,aI0,cIL
(1 atm, 1173 K)GWO1.08500.08070.00070.164128.208413.3625167.64134.1991 × 10−3
PSO1.0868965910.056100.3998306.79402008.9389 × 10−3
ARO1.08610.02860.00300.075329.95874.6816160.10672.4342 × 10−4
(1 atm, 1273 K)GWO1.08200.05713.73 × 10−80.248819.56958.6351178.37673.9708 × 10−3
PSO1.08210.044900.0366305.29472006.247 × 10−3
ARO1.08120.02830.00230.084929.42724.21366160.37492.2096 × 10−4
(3 atm, 1173 K)GWO1.11490.060900.035122.19038.4946193.50374.5549 × 10−3
PSO1.11220.054700.038812.448512.44852004.8315 × 10−3
ARO1.11400.02950.00300.077529.98814.8380160.37713.2846 × 10−4
(3 atm, 1273 K)GWO1.10320.114300.133229.508325.8728163.57883.3413 × 10−3
PSO1.11120.036200.37448.19478.19472005.7429 × 10−3
ARO1.11200.02880.00230.085429.99994.1656160.39582.8078 × 10−4
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

Guo, Z.; Ye, Z.; Ni, P.; Cao, C.; Wei, X.; Zhao, J.; He, X. Intelligent Digital Twin Modelling for Hybrid PV-SOFC Power Generation System. Energies 2023, 16, 2806. https://0-doi-org.brum.beds.ac.uk/10.3390/en16062806

AMA Style

Guo Z, Ye Z, Ni P, Cao C, Wei X, Zhao J, He X. Intelligent Digital Twin Modelling for Hybrid PV-SOFC Power Generation System. Energies. 2023; 16(6):2806. https://0-doi-org.brum.beds.ac.uk/10.3390/en16062806

Chicago/Turabian Style

Guo, Zhimin, Zhiyuan Ye, Pengcheng Ni, Can Cao, Xiaozhao Wei, Jian Zhao, and Xing He. 2023. "Intelligent Digital Twin Modelling for Hybrid PV-SOFC Power Generation System" Energies 16, no. 6: 2806. https://0-doi-org.brum.beds.ac.uk/10.3390/en16062806

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