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

: Hydrogen (H 2 ) energy is an ideal non-polluting renewable energy and can achieve long-term energy storage, which can effectively regulate the intermittence and seasonal ﬂuctuation of solar energy. Solid oxide fuel cells (SOFC) can generate electricity from H 2 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 identiﬁcation of unknown model parameters. Hence, in this study, an artiﬁcial rabbits optimization (ARO)-based parameter identiﬁcation 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 identiﬁcation 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.


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 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.

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: 1.
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; 2.
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;

3.
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.

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.

PV Cell Modelling
The equivalent model of DDM is illustrated in Figure 1.DDM consists of an ideal constant current source, two diodes D 1 and D 2 , a series resistance R s , and a shunt resistance R sh .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 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.

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.

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]: where  and  indicate the current passing through the shunt resistance  and series resistance  ,  denotes the photocurrent.In particular, the current  and  flowing through  and  and  are expressed as: The output of DDM can be given by [20]: 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 d1 and I d2 flowing through D 1 and D 2 and I sh are expressed as: In Equations ( 2) and ( 3), V T is given by [35]: Hence, the output current of the PV cell based on DDM is presented as [16]: 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: The error functions f (V L , I L , x) for DDM is described as: 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 .

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.In order to achieve simple and effective control of the SOFC, the electrochemical model is widely used, whose output voltage can be expressed as: where  indicates the open circuit voltage;  denotes the number of series cells;  ,  , and  indicate activation voltage loss, ohmic voltage loss, and concentration voltage loss, respectively.
Based on Butler-Volmer equation, the  can be expressed as: where  indicates the slope of Tafel line;  ,  , , and  , represent load current density, the exchange current density of anode and cathode, respectively.Moreover,  ohm can be described by: where  denotes the ionic resistance.Furthermore,  is described by: At anode: At cathode: In order to achieve simple and effective control of the SOFC, the electrochemical model is widely used, whose output voltage can be expressed as: 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: 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: where R ohm denotes the ionic resistance.Furthermore, V con is described by: where B indicates a constant; I L represents the limiting current density.
From Equations ( 11)-( 14), the output of SOFC is described by: 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: 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.

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

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.

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.

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]: 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). (21) where  ⃑  1 denotes the candidate position of the th rabbit at the time t + 1,  ⃑  represents the position of the th 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.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.

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 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: where b i,r (t) denotes a randomly selected burrow, r 4 and r 5 represent two numbers ran- domly distributed between 0 and 1.

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: 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].

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).

Design of ARO for PV Cell Parameter Identification
In order to achieve reliable optimization, the optimization variable  is limited to the upper limit  and lower limit  as follows: where J represents the optimizing variables number.

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).

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 max j and lower limit x min j as follows: where J represents the optimizing variables number.
If an individual rabbit violates Equation ( 31), it will be randomly reset by: 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.
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 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 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.

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

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.  Convergence speed and stability are both important criteria for measu formance of an optimization method.The convergence curves of the RMSE of and ARO for DDM are demonstrated in Figure 9.Although WOA and AB quickly in the beginning stage, on the other hand, they also easily stagnate suffer from prematurity, leading to low-quality local optimums.However, A have a higher quality global optimum and a more stable convergence proce adaptive adjustment between global exploration and local exploitation throu eration, thus avoiding falling into a local optimum.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/m 2 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.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.

RMSE(log)
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.

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.

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.Tables 3 and 4 show the dataset size and the parameters range, respectively.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.Tables 3 and 4 show the dataset size and the parameters range, respectively.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.Tables 3 and 4 show the dataset size and the parameters range, respectively.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.

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 ( Meanwhile, Figure 13 demonstrates the boxplot and convergence graphs of vari methods under 100 independent runs.Figure 13a illustrates the RMSE results distribu of different algorithms, it can be verified that ARO owns the smallest distribution ra with minimal lower and upper bounds.Figure 13b also proves that ARO has a more sta convergence process and higher convergence precision.

Verification on a 5 kW Cell Stack
The effectiveness of ARO was verified on a 5 kW cell stack here.Furthermore, estimated results under four various operating scenarios are shown in Table 6, wh demonstrates that the identification results of ARO can obtain the smallest RMSE un all operating scenarios compared with GWO and PSO.For example, the RMSE obtai by ARO is only 2.72% and 4.88% to that of PSO for the (1 atm, 1173 K) and (3 atm, 127 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.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.

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:

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: 1.
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; 2.
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; 3.
ARO was applied for parameters estimation for the DDM 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: 1.
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; 2.
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; 3.
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.

Energies 2023, 16 , 2806 6 of 21 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 .

Energies 2023 , 21 Figure 5 .
Figure 5. Searching behavior adjustment based on energy factor A.

Figure 5 .
Figure 5. Searching behavior adjustment based on energy factor A.

Figure 5 .Figure 6 .
Figure 5. Searching behavior adjustment based on energy factor A.

Figure 6 .
Figure 6.Schematic diagram of probability calculation of detour food searching.

Energies 2023 ,Figure 8 .
Figure 8.Output I-V and P-V curves obtained by ARO: (a) I-V curve; and (b) P-V curv

Figure 8 .
Figure 8.Output I-V and P-V curves obtained by ARO: (a) I-V curve; and (b) P-V curve.

Figure 9 .
Figure 9. Convergence curves obtained by various algorithms for DDM.

Figure 9 .
Figure 9. Convergence curves obtained by various algorithms for DDM.

Figure 12 .
Figure 12.Datasets for 5kW cell stack under various operating scenarios.

Figure 11 .
Figure 11.Datasets for validation for electrochemical model.

Figure 11 .
Figure 11.Datasets for validation for electrochemical model.

Figure 12 .
Figure 12.Datasets for 5kW cell stack under various operating scenarios.

Figure 12 .
Figure 12.Datasets for 5 kW cell stack under various operating scenarios.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 [E o , A, R ohm , B, I 0,a , I 0,c , I L ] = [1.15,0.02, 0.0004, 0.03, 12, 4, 152].The identification parameters results are illustrated in Table5along with the RMSE, in which the RMSE values can give a direct and clear indication of the optimization accuracy.As shown in Table5, 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.

Table 2 .
Parameter identification results obtained by various algorithms under DDM.

Table 3 .
The size of dataset.

Table 3 .
The size of dataset.

Table 3 .
The size of dataset.

Table 4 .
The range of unknown parameters.

Table 5 .
Identification results of various algorithms.

Table 6 .
Identification results of various algorithms for 5 kW cell stack dataset.