Fault Coverage-Aware Metrics for Evaluating the Reliability Factor of Solar Tracking Systems

: This paper presents a mathematical approach for determining the reliability of solar tracking systems based on three fault coverage-aware metrics which use system error data from hardware, software as well as in-circuit testing (ICT) techniques, to calculate a solar test factor (STF). Using Euler’s named constant, the solar reliability factor (SRF) is computed to deﬁne the robustness and availability of modern, high-performance solar tracking systems. The experimental cases which were run in the Mathcad software suite and the Python programming environment show that the fault coverage-aware metrics greatly change the test and reliability factor curve of solar tracking systems, achieving signiﬁcantly reduced calculation steps and computation time.


Introduction
The most abundant form of energy on our planet is incontestably the Sun. Despite dropping technology costs and free sunlight, conventional companies have been reluctant to use photovoltaic (PV) panels such as those installed on residences' rooftops for large-scale projects. Concerning this aspect, the main explanation is given by the fact that researchers have developed metrics for producing electricity from coal, gas, oil, and nuclear energy, but not from PV panels, a rapidly emerging technology whose cost depends on the Sun's availability. Fortunately, the solar industry has started to follow a standardized metric that consumers can understand, generally referred to as the levelized cost of energy (LCOE) [1]. In simple terms, LCOE is measured in cents per Kilowatt-hour (kWh), takes into account all the costs of a solar plant over its lifespan, compares them to the total electrical output, and offers an accurate, long-term financial overview. This metric is commonly recognized by the power industry, service companies, and government agencies worldwide. The solar industry has traditionally used an alternative metric, called dollars per watt [2] which compares the initial costs of a solar installation against its potential peak power output. This metric, if taken on its own, lacks the actual energy generated by the solar system and many other variables which are the strong points of the LCOE metrics.
The multitude of LCOE variables [3] can be grouped into five main categories: (a) location (for Sun exposure, seasonality, and local benefits); (b) cost of capital (spanning interest, taxes, and leverage); (c) energy production (solar panel and automated equipment performance evaluation); (d) system costs (similar to dollars per Watt, but covers the total installation costs as well as warranties); (e) operations and maintenance costs (cleaning, monitoring and repairing). Since it is out of scope for this paper to cover all possible categories offered by the LCOE metrics suite, our field of research is restricted to the energy production and maintenance costs domains. Conceptually, regarding the energy production domain, more daylight, more powerful PV panels and more reliable Sun availability translates into more electricity. However, a few aspects are counter-productive, as cold days are better than hot days for solar production, since PV cells overheating can lower energy output [4]. To compensate for the environmental weather conditions, solar tracking systems are used nowadays to steer PV panels according to the Sun's direction, thereby increasing energy production which will significantly boost the return on investment. Regarding operations and maintenance costs, one of the key elements for long-term performance from a solar tracking device is the requirement of high reliability, remote monitoring, proactive maintenance, and service [5]. Thus, proper selection of electrical equipment is crucial for constructing a robust and long-lasting solar tracking system that will maximize solar energy collection. Since solar trackers can be assembled from individual circuits such as microcontroller units (MCUs, Arduino, SMT32 boards), motor drivers (L298N, ULN2003A), DC-DC and stepper motors, it is vital to emphasize that each of these electrical components can become faulty during their operation. Hence, several testing methods were developed to detect software [6], hardware [7] as well as in-circuit [8] errors in dual-axis solar tracking equipment, and to ultimately counteract the negative effects of a malfunctioning or, in more critical scenarios, system failure.
Since LCOE calculation formulas are overly complex for computing the reliability of solar tracking systems, this paper proposes a set of three fault-coverage aware metrics that calculate an STF based on the number of detected errors and determines the SRF that categorizes and certifies the quality of modern, and high-performance dual-axis solar tracking systems.
The paper is organized as follows: Section 2 presents several works related to recent metrics developed to assess the efficiency and reliability of solar panels. Section 3 describes the most common metrics for evaluating the meantime to failure as well as the reliability of automated solar systems according to the well-known Weibull distribution model. Section 4 introduces the proposed metrics for calculating the STF as well as the SRF of solar tracking systems. Section 5 presents the experimental setup and results. Finally, Section 6 concludes this paper.

Literature Review
Even though LCOE metrics became today's standard for choosing reliable and efficient solar tracking systems [9,10], the author in [11] demonstrates that LCOE inaccurately disfavors renewable energy sources. Furthermore, based on his research, the LCOE equation which was formulated by the U.S. National Renewable Energy Laboratory (NERL) overprices solar energy by 27%, and wind energy by 18% as compared to natural gas-based power, showing that the LCOE model will be potentially replaced in the future by new metrics called "present value of the cost of energy" (PVCOE).
Recent advancements in the field of renewable energy show that besides the standard LCOE metrics formulas, other solar-energy oriented metrics have been developed for assessing the performance and reliability of static solar PV systems. Along with the solar tracking methods which are used to optimize the PV panel's position towards the Sun, new techniques such as solar power forecasting were introduced to the market to improve the performance and reliability of solar PV panels. Thus, the authors in [12] present a collection of value-based solar forecasting metrics that can be applied to a variety of input variables (various time horizons, geographic locations, and applications) to improve the accuracy of solar forecasting capabilities. Additionally, the authors implement a framework for a systematic approach to evaluate the sensitivity of the proposed metrics to three types of solar forecasting improvements (uniform improvements, ramp forecasting magnitude improvements, and ramp forecasting threshold adjustment). Their experimental results show that the proposed metrics offer an efficient evaluation of solar power forecasts and highlight the economic and reliability effects of improved solar energy forecasts. Researchers from the U.S. NREL associate the term reliability with stability in their work [13], which involves a study of transient stability, frequency response, control stability, and the impact of concentrating solar power (CSP) on grid reliability. Their investigation demonstrates that large amounts of CSP units do not pose any challenge to grid stability, and that frequency response can be enhanced significantly by frequency sensitive controls on CSP and solar PV installations. Furthermore, another team of research scientists from the same institute highlights the role of reliability and durability in PV system economics [14]. By improving the conventional LCOE metrics formula with the system architecture cost, a capacity factor (system location, orientation, tracking), and a reliability and durability expression, they were able to determine the total life cycle cost divided by the total lifetime energy production.
On the other hand, the authors in [15] propose an optimal modeling and metrics design for solar energy-powered base stations, proving that 4G cellular communication units can benefit from complete clean energy in the future. The authors establish a model for describing the dynamic behavior of solar energy-based base stations by using a stochastic queue model, and further develop three performance metrics (service outage probability, solar energy utilization efficiency, and mean depth of discharge) which will form the key element for their minimization algorithm that takes into consideration the communication reliability, efficiency, and durability. The United Nations Sustainable Development Goals no. 7 indicates a major global task: "Ensure access to affordable, reliable, sustainable, and modern energy for all" [16]. The authors in [17] adhere to this goal and address the problem of reliability cost in decentralized solar power systems in sub-Saharan regions. By employing a multi-step optimization algorithm based on the LCOE metrics, the authors achieve a 90% reduction of energy cost for household solar systems with battery storage in Africa.
Recent efforts in the field of green energy have been also made to introduce the term "availability" to the standards, contract language, and performance models of PV systems. Therefore, the authors in [18] define availability as the operational status of an electrical energy generating technology, which gives a hint about the component health and system reliability (uptime, downtime, and condition stages). Moreover, availability is described as the fraction of time in which an industrial unit is capable of providing service to generate energy, so that reliability in this scenario is purely given by a probabilistic model, which will be extended later in this paper. Following, the authors in [19] present a comprehensive study of a large PV system involving field data, failure times, and repair times for five years. Their reliability program, which strengthens the interaction between reliability model, failure model, and accelerated tests, shows that periodic maintenance and repair operations are vital for system availability since the granularity of the proposed model is given by the level of collected failure data. These scientific findings direct us to a new approach in the reliability domain where the predicted number of component failures becomes extremely important, thus opening the path to developing new metrics based on fault coverage. This paper distinguishes itself from the above-mentioned works by proposing three metrics that take into consideration the fault coverage obtained from the testing methods that were applied to modern dual-axis solar tracking systems. More precisely, an STF is computed that considers software, hardware, and in-circuit errors, thus determining the SRF that defines the quality, reliability, robustness, and durability of mobile solar trackers.

The Weibull Distribution Model for Calculating the Reliability of Dual-Axis Solar Tracking Systems
During the construction of various automation equipment, the requirement to ensure their safe operation has led to the use of certain safety factors in their design. The notions of reliability and safety in operation have emerged and highlighted the concern for finding mathematical models and calculation methods that allow making the most accurate predictions regarding the behavior, over a certain period, of installations and equipment in conditions of known exploitation. Reliability, from the maintenance point of view, is a probability representing the ability of a product to function without defects, involving a slight chance of failure, and can be expressed mathematically as in Equation (1) [20]: where R is the probability that the product will continue to meet the specifications, and F denotes the probability of failure, both variables being time functions, depending on the parameter t. The advantage of relation (1) is that it can be applied to systems that are fully functional or fully defective. In our context, if a chain of solar tracking systems contains irreparable products (most unfavorable scenario), the corresponding reliability curve will look like the one presented in Figure 1. As can be seen in Figure 1, the computation of the reliability factor denoted with R relies on N products (solar tracking systems) and the final moment, denoted with T, in which all the components of the products are completely damaged.

Failure Frequency (Hazard Rate)
The failure frequency, also known as hazard rate, describes how many products are out of service (malfunctioning) in a given period. The mathematical expression associated with the failure frequency is given by Equation (2) [20]: where λ denotes the failure frequency, N the total number of products, T i the corresponding time failures for each of the N products (industrial components), and the Mean Time to Failure (MTF) represents the average time to failure. The reliability equation will be written as in (3) [20]: where N-i designates the number of surviving products and N the total number of products. The reliability curve diagram can be drawn accordingly as presented in Figure 2. The derived linear reliability function can be expressed as in (4) [20]: where a denotes a random coordinate of the X-axis, and b denotes a random value from the Y-axis. First, according to Figure 2, the reliability linear function will be computed at moment t = 0, to determine the coordinate b, as presented in Equation (4a): Secondly, since we know that at t = T (final moment) the reliability will be 0, coordinate a can be computed as presented in relation (4b): According to relations (4a) and (4b), the reliability linear function can be rewritten as in Equation (5): Consider, for a more concrete case, that a test was performed on a chain of 15 PV sensors (in our case, PV cells) that are used to automate solar tracking equipment. The test time is measured in units of days. The measurement times were collected and ordered in Table 1. First, for all the collected failure times, the associated MTF and hazard rate can be determined for the solar tracker's sensor modules. According to Equation (2), the MTF formula can be rewritten as in (6) [20]: In other words, the MTF is given by the total test time divided by the total number of failures during testing phases. The MTF expression is calculated according to the catalog data from Table 1, as in expression (6a): The corresponding failure frequency will be determined from relation (2) as seen in expression (6b): To establish the relative error (E) of our previous calculus, the MTF will be calculated a second time by using the following linear interpolation function of the reliability that depends on the time parameter as presented in relation (7) [20]: Thus, the new MTF formula will be written as in Formula (8): By subtracting the MTF 1 from the initial MTF relation and by dividing the result by the average time to failure, the expression (9) will be obtained as follows: The approximation method helps determine the relative error for determining the MTF, which in our case is 0%, showing that the failure frequency calculus is accurate in establishing the time to failure model for irreparable products. For reparable components, however, failure times can be replaced by maintenance operational times which allow the fixing of damaged products, as can be seen in Figure 3. The mathematical relation for mean downtime (MDT) will be further established, as in Equation (10) [20]: where T Di denotes failures or maintenance times, N F the total number of defects, and T, the total functioning time of N products. Accordantly, the operating time of the service products can be expressed as in the mathematical rule (11): where N represents the number of components, T is the time interval, while N F and MDT having the same meaning as in Equation (10). Consequently, at this point the Mean Time Between Failures (MTBF) formula can be defined, as in Equation (12) [21]: where the MTBF is the total operating time of all components divided by the number of occurring defects. The associated failure frequency (hazard rate) can be thus calculated with the relation (13):

Availability (Efficiency)
Availability is defined as the ability of a system/equipment to be in working order under given conditions, at a given time, or for a given period, assuming that the provision of the necessary external means (maintenance) is ensured [13]. Availability is considered a complex form of assessing the quality of a system/product because it includes both reliability and maintainability. Mathematically, availability can be expressed as in Equation (14) [20]: Based on Equation (12) the formula for the product N·T can be rewritten as in relation (15): According to Equation (15) we can now write the final form of A from (14): where, for convenience, the MTF and MTBF parameters can be used to characterize the automated system. In conclusion, a product or component is available as long as it can provide service, and unavailable when it is undergoing maintenance operations. As seen in Equation (1), the probability model of availability and unavailability can be expressed as in the mathematical rule (17) [20]: More generally, we will write unavailability as in the Equation (18): where unavailability depends only on the MDT divided by the total time (operating time + malfunctioning time). Consider, for a real-life scenario, that a company acquires 30 identical motor drivers that are used in solar tracking systems and perform a test on them for 3 months, recording the failure times (maintenance times), measured in days, in Table 2. According to Table 2, the active number of failures is N F = 40. The total time is calculated for N T = 30 ·90 ·24 = 64800 h = 2700 days. Accordingly, we can compute the MDT parameter with the relation (19): The total failure time will then be calculated with the Formula (20): Furthermore, the MTBF will be determined with relation (21): The failure frequency will be computed with Equation (22): By relying on the previous standard metrics, the general conclusion is that the MTBF is 19.9 days and the MTF is 3.225 days, which means that once every approximately 19 days one motor driver malfunctions and requires 3.225 days for maintenance operations, resulting in an 86% availability rating and 14% unavailability rating.

Typical Variation Form of the Hazard Rate (Bathtub Curve)
The variation in time of the reliability function (R) in regards to the hazard rate (λ) can be expressed mathematically as in Equation (24) [21]: It was already stated that at moment t = 0 all the products function according to their operating specifications so that the relationship between reliability and the instantaneous failure rate is established with the help of the equation system (25) [21]: Generally, the MTF relation will become as presented in (26): The MTF Formula (27) is given consequently by the subgraph area: The graphical representation of Equation (26) can be seen in Figure 4. This graphical representation applies to any solar tracking system (automation equipment) which is described by the run-in period (usually around 6 months) followed by the normal operating period (warranty period, which depends on the quality of the acquired product) and, lastly, the catastrophic defect area where generally, a solar tracking system will require several maintenance human interventions.
Finally, there is a second and last approach for analyzing the typical variation form of the hazard rate, given by the so-called Weibull distribution model, found in relation (28) [21]: According to Equation (28), the relation for the reliability variation can be written in the Equation (29): where β is a positional parameter, η is a scalar parameter, t 0 is the initial moment, and exp represents the exponential form of Euler's constant named value.

Applications
The abovementioned metrics can be applied to a variety of real-life scenarios regarding the calculus of reliability, MTF, MDT, and hazard rate of solar tracking systems. More specifically, to establish accurate predictions about the reliability of the entire solar tracking system, it is of major importance to formulate problems concerning the failure rates of a solar tracker's components. For instance, consider that an MCU (e.g., Arduino UNO) which is used to control solar tracking equipment, has a constant failure rate of 0.1 / year, and we want to determine the reliability, respectively the probability of failure for 1, 2, and 3 years. This problem can be solved by just relying on the equation set (30) as follows: The aforementioned solution shows that the reliability of the Arduino UNO board decreases in years by an average value of 8%, while the chance of failure increases every year by approximately 9%. Another way to analyze the behavior of solar tracking systems is from the defects points of view, according to a Weibull distribution given by the parameters t 0 = 0, η = 625 days and β = 1.5. It is considered that for the given input parameters, a test was performed on 1000 random components that were mass-produced to construct solar trackers. The goal here is to determine the survival (operating) probability of these components after t = 10 4 h, by applying the Weibull equation, as presented in (31): Additionally, for another queue of 1000 random components that were mass-produced in the industrial environment, it is of interest to determine the period after which 100 components became faulty. It is assumed that the component testing was started at the moment t = 0. Therefore, if t 0 = 1000 components, moment t 1 will correspond to the subtraction of 100 components from the initial stack, in other words, t 1 = 1000 -100 = 900 remaining components. The associated reliability R(t 0 ) = 1 (as a general rule) and R(t 1 ) = 0.9. The expression of R(t 1 ) is expanded according to Formula (32): Following, the last obtained expression for t 0 = 0 will be computed as in (33): Finally, a more convenient manner to represent the number of faulty components is through a graph chart which shows the evolution in time of 10 pairs containing 100 components each, until the final moment where all products are completely damaged, as can be seen in Figure 5.

Proposed Fault Coverage-Aware Model for Calculating the Solar Reliability Factor
Since the standard Weibull method relies initially on the computation of time to failure models such as MTF, MDT, and MTBF which are usually applied to individual components, this paper proposes a set of three fault coverage-aware metrics that make use of experimental data gathered from hardware, software and in-circuit oriented testing techniques.
The aim of testing electronic solar tracking equipment is to identify errors that may occur in digital circuits. Reliability analysis is crucial for businesses as well as regular consumers for choosing robust and long-lasting components. To satisfy these needs, the reliability of solar trackers is calculated in two major ways: (a) by using hardware testing equipment such as built-in self-test (BIST) modules and (b) by focusing on the calculus of reliability based on white-box software testing (WBST) algorithms.

Fault-Coverage Aware Metric for Determining the Reliability of Electronic Components with the Help of BIST Devices
BIST architectures are today's preferred hardware testing routines due to their facile deployment and integration in various digital circuits, thus eliminating the need of employing standalone automated test equipment (ATE). A static deployed BIST architecture is generally composed of three components: (a) test pattern generator (TPG) which injects test vectors into the circuit chain, (b) device under test (DUT), sometimes also referred to as circuit under test (CUT), and (c) output response analyzer (ORA) which collects the outputted signals of the CUT and generates corresponding signatures for each of the patterns, as can be seen in Figure 6. Usually, a BIST architecture makes use of a linear feedback shift register (LFSR) as a TPG unit and a multiple input signature register (MISR) as an ORA unit, both hardware devices relying on the deployment of D-type flip flops (FFs). More important, however, is the fact that the total number of deployed FFs in the MISR unit increases the chance of detecting hardware errors, implicitly expanding the fault coverage domain. Therefore, during the BIST routine's testing phase, collected data (signatures) from the experimental cases is used to outline the fault coverage. Further, an STF parameter is defined, which will be expressed mathematically as in Formula (34): where T V represents the number of executed test vectors, N E the number of errors per test case, T P the total number of test patterns, and D the number of identical devices used for detecting errors (which in our case are given by the FF units). Since the generalized STF formula can be adapted to various test scenarios, several STF properties will be listed as follows: (1) The number of errors N E is always equal, lower, or greater than T V in any test scenario regardless of the noise factors that are considered, and can be expressed mathematically as in (35): (2) The number of test vectors T V is always equal to or lower than the total number of test patterns T P given by the expression (36): (3) The total number of test patterns is calculated with the Equation (37): (4) If no errors are detected during the testing routines, the STF parameter will be 0 resulting in 100% reliability and availability of the tested equipment, expressed in Equation (38): To demonstrate the validity of these properties, a few concrete cases are presented further. Regarding the STF, the following considerations are listed: (a) for a total number of T V = 7 (test cases) the test system has successfully identified N E = 10 bit-flip errors. It should be noted, according to property (1) that the number of errors detected may be greater than the number of test cases because multiple errors (burst errors) may occur within a single test vector. Additionally, for (b) a number D = 4 FFs are used in the structure of the MISR in this paper. According to property (3), the total number of test cases T P is calculated based on the Formula (39): Since T V = 7 and T P = 10 property (2) is also satisfied. Finally, by having this hypothetical data the variables from Equation (34) can be adequately substituted as in (40): Supposing that during the detection stage, all the errors were corrected by using the Hamming code algorithms, parameter N E will be reduced to the value 0. The corresponding equation will be expressed as in (41): thus satisfying property (4) of the first set of metrics. In the last stage the SRF will be calculated with the Formula (42): Since the SRF expression is the exponential of the STF equation, the entire relationship can be rewritten as in expression (43): At this point, the SRF of the automated sun-tracking equipment will be computed as in Equation (44): Similarly, if the STF parameter is reduced to 0, the mathematical expression will become as in expression (45): The final mathematical relation connects the SRF parameter with the initially computed STF variable according to property (4), thus showing that hardware components that are not affected by permanent errors can achieve 100% reliability and availability. Any other types of defects caused by transient errors can be considered negligible since their footprint on the electronic components is only temporary and does not cause any serious harm to the overall functionality of the tested devices [6].

Fault-Coverage Aware Metric for Determining the Reliability of Electronic Components with the Help of WBST Methods
Traditional testing processes aim to generate test vector sets that apply to a specific fault model to optimize coverage while reducing the execution time of the test application. Because the Off-Line test is typically performed after the circuit is assembled as part of a more in-depth manufacturing test, it is often used in maintenance checking for extending the lifespan of the product. Roughly equivalent to off-line testing, on-line testing, also referred to as parallel or concurrent error detection, is a test technique used to permanently verify the integrity of the CUT, as can be seen in Figure 7. One of the preferred online software testing techniques is known under the acronym of WBST routines which periodically checks the CUT for transient or permanent errors. A more particular case of the WBST method is the cloud-based software testing technique which makes use of an Internet of Things (IoT) software TPG for injecting test vectors in MCU devices which are part of solar tracking equipment. By stimulating the inputs of the CUT the outputs are collected in a results gatherer where with the help of an error indicator, the user is capable to establish the fault coverage of the employed WBST method. Moreover, as expected, the previous metric set presented in Section 4.1 for our software test scenarios can be adapted adequately. The STF formula corresponding to the WBST routines will be written as in Equation (46): where N E represents the number of errors per test vector, T V denotes the number of considered test vectors, T P provides the total number of test patterns and B designates the number of breakpoints/software functions that are implemented across the algorithm debugging stage. Besides this minor difference from the proposed hardware metric, a few other properties which are important for demonstrating the validity of the software-based metric will be listed below: (1) The number of errors N E is always equal, lower, or greater than T V in any test scenario regardless of the number of considered breakpoints, and can be expressed mathematically as in statement (47): (2) The number of test vectors T V is always equal to or lower than the total number of test patterns T P given by the expression (48): (3) If no errors are detected during the testing routines the STF parameter will be 0 resulting in 100% reliability and availability of the tested equipment, expressed as in Equation (49): As already anticipated, the second metric removes one of the listed properties from Section 4.1 since there is no clear mathematical relationship between the total number of test patterns T P and the number of breakpoints B. At this point, the validity of the remaining three properties can be demonstrated in a real-life scenario. Regarding the STF, the following considerations are made: (a) for a total number T V = 7 (test cases), the test system has successfully identified N E = 10 calculation errors. Additionally, for (b) a number of T P = 10 test patterns and a number B = 10 breakpoints in the software code, meaning that the test system successfully identified all calculation errors with the deployed software functions. Because properties (1) and (2) are already satisfied from the above statement, let the STF parameter will be calculated as in relation (50): In the last stage, the SRF parameter will be computed with the mathematical rule (51): showing that with almost the same hypothetical set of data from Section 4.1, the reliability of the solar tracking equipment is halved. This is because calculation errors in the software code usually manifest in the behavioral model of the hardware device (incoherent stepper motor rotations leading to improper positioning of the solar panel thus resulting in heavy energy harvesting losses). Along with the satisfying property (3) of the second set of metrics, consider that during the detection stage all the errors were manually corrected by the software engineer, thus minimizing NE to the value 0. The final Equation (52) will become:

Fault-Coverage Aware Metric for Determining the Reliability of Electronic Components with the Help of ICT Methods
ICT is a type of white-box testing procedure where an electrical probe navigates the entire printed circuit board (PCB) area to identify shorts, opens, resistance, capacitance, and other basic electrical components which give a hint of the electronic assembly which was fabricated according to its catalog specifications [8]. The ICT process is conducted with the help of industrial machines such as ATE, bed of nails fixtures, or with a fixtureless ICT setup (one example here being a flying probe device). In work [8], a flying probe inspired in-circuit tester (FPICT) was constructed and deployed to verify the integrity of the entire solar tracking electrical equipment. By implementing a sensorless solution for locating all test points of the PCB, the proposed test device was capable to monitor critical paths and parameters (such as current and voltage drops), thus determining a fault coverage of the malfunctioning CUTs, as can be seen in Figure 8. To demonstrate the flexibility of the metrics set, the mathematical expressions from Sections 4.1 and 4.2 will be adapted to the ICT test scenarios. The STF formula corresponding to the ICT routines will be written as in expression (53): where N E represents the number of errors per test round, T R denotes the number of considered test rounds, N R provides the total number of test routines and P designates the number of probes (nails) that are equipped to the FPICT device. Similarly, as in the previous sections, the properties which are essential for validating the ICT-based metric will be listed below: (1) The number of errors N E is always equal, lower, or greater than T R in any test scenario regardless of the number of probes that are used during testing, and can be expressed mathematically as in expression (54): (2) The number of test rounds T R is always equal to or lower than the total number of test routines N R , given by the expression (55): (3) If no errors are detected during the testing routines the STF parameter will be 0 resulting in 100% reliability and availability of the tested equipment, expressed as in relation (56): Consider a real-life scenario where all possible test points' voltage deviations must be identified. For this purpose, the test system implements a number of T R = 10 rounds for each test point, a total number of N R = 100 test routines, and a number of P = 2 probes, to identify N E = 12 voltage deviations. The STF parameter will be computed in relation (57), based on the previous configuration: With the abovementioned result properties (1) and (2) are successfully satisfied. To demonstrate property (3), a similar but fault-free DUT is considered, leading to the Formula (58): The SRF parameter will be calculated as in Equations (59) and (60): The above described metrics reduce significantly the number of computations executed during the Weibull distribution algorithm, as well as the number of parameters that are involved in the mathematical rules, proving that fault coverage is a promising reliability and availability indicator.

Experimental Setup and Results
This section presents the experimental setup and results as well as a comparison between our work and other related manuscripts from the targeted domain.

Fault Coverage-Aware Metrics Equation Reduction for Computing the Reliability Test Factor
With regards to our simulation setup, the experimental data is obtained from works [6][7][8], to outline the efficiency and fast computation of the previously described hardware and software fault coverage-aware metrics. Before describing the experimental environment the number of parameters in each equation is narrowed down, to adapt the metrics to the BIST and WBST test scenarios.
First, according to Equation (34) for the BIST routines, it was established that the STF parameter is directly proportional to the number of test vectors multiplied by the number of errors (per test vector) and inversely proportional to the total number of test patterns multiplied by 2 to the power of D (where D denotes the number of FFs devices used during the testing stage). At this point, Equation (34) will be constrained by the properties mentioned in Section 4, a reason for which the product N E · T V will be denoted with E, representing the hardware error factor. According to property (3), Equation (34) can be rewritten as in relation (61): Secondly, by knowing the STF parameter for the WBST routines, the following associations can be realized similarly as in Equation (62): where E designates the error factor for software testing, T P the number of software test patterns, and B the number of breakpoints used during the debugging stage. Finally, by adapting the STF parameter for the ICT routines, the formula can be simplified as in relation (63): where E represents the error factor for ICT, N R the total number of test routines, and P the total number of test probes. Thus, by reducing the equations to the simplest form, we can proceed with the description of the simulation environment and the reliability graph generation of the tested solar tracking system.

Fault Coverage-Aware Metrics for BIST, WBST, and ICT Test Scenarios
Based on our previous efforts to verify the functionality of the solar tracking equipment described in [6][7][8], we can further use the collected experimental data to validate the proposed fault coverage-aware metrics. Since BIST-oriented architectures and WBST algorithms are applied to check various digital devices involved in the automation process of the solar tracking device, the results will be organized similarly to the Weibull distribution model. Even though the experimental results were gathered only from a personalized solar tracking prototype, we firmly sustain the idea that the proposed fault coverage-aware metrics can be generalized for several solar tracking systems available on the market today. Regarding the experimental data obtained from the online built-in self-test (OBIST) architecture presented in [7], the total number of test patterns and detected errors are synthesized in Table 3. According to the data in Table 3, it can be observed that the total number of patterns T P is a constant value (see column 5) and the error factor E is a variable parameter (see column 6). Thus, to generate the graphical representation of the STF parameter, the following equation system (64) must be computed: where D is also a constant value that comprises 16 FFs units involved in the OBIST routines. Therefore, the following results are presented in the equation system (65):  Table 3 (see column 3). In other words, the mathematical relation between the two parameters can be written as seen in Formula (66): The abovementioned equation proves that the defined STF parameter is the fault coverage percentage written in fractional form. By ordering the results in the equation set (64), the graphical representation of the STF factor will be obtained in form of a linear function, as can be seen in Figure 9 (scenario a). As expected, the more errors are detected utilizing hardware testing techniques (OBIST in our case), the more the STF linear function will increase, according to the fault coverage principle.
After computing the equation set (65), it can be immediately noticed that the STF parameter is tightly coupled to the fault coverage (FC) domain in Table 3 (see column 3). In other words, the mathematical relation between the two parameters can be written as seen in formula (66): 100 100 The abovementioned equation proves that the defined STF parameter is the fault coverage percentage written in fractional form. By ordering the results in the equation set (64), the graphical representation of the STF factor will be obtained in form of a linear function, as can be seen in Figure 9 (scenario a). As expected, the more errors are detected utilizing hardware testing techniques (OBIST in our case), the more the STF linear function will increase, according to the fault coverage principle.  Following, the SRF parameter can be computed according to the equation set in (67) where STF 1 , STF 2 , STF 3 , STF 4, and STF 5 were previously calculated in the equation set (64).
Thus, the graphical representation is a linear function that describes the reliability factor of the entire solar tracking system, as shown in Figure 9 (scenario b).
All the presented results were simulated and computed in the Mathcad v15 environment due to its extensive and facile engineering tools.
Furthermore, it was determined that the granularity of each graphical representation can be controlled by adding more test case scenarios to the fault coverage database. Similarly, with regards to the On-line WBST presented in [6], the experimental results in Table 4 were synthesized according to four test batches which contain a huge amount of data packets. For each of the listed data batches depicted in Table 4, we associated a number of 10 breakpoints/software functions that are distributed along with the algorithm to detect as many error types as possible.
Here, it can be observed that the total number of test patterns is T P = 4334 (computed by adding all test vectors together from column 6), and the variable number of errors E is calculated for each data batch. By using the simplified Equation (62), the equation set for calculating the STF parameter in regards to the WBST routines will be constructed as in equation system (68): By solving the equation set (68), it can be immediately noticed that the STF parameter for software test scenarios is completely different from the FC (see column 7). Therefore, the graphical representation of the STF for WBST routines is a curve distribution, as can be seen in Figure 10 (scenario a). Finally, the equation set (69) for the SRF parameter can be computed accordantly: As can be seen in Figure 10 (scenario b), the graphical representation of the STF parameter is a linear distribution and proves that improved fault coverage results in high reliability because software errors can be corrected manually by the software engineer, in comparison to hardware faults which require in certain scenarios the total replacement of the components. Finally, by extending the test scenarios for the FPICT method [8], the experimental dataset presented in Table 5 will be obtained. According to the four multiple point testing batches seen in Table 5, the error factor E can be determined based on the equation set (70): Finally, after solving the equation set (71), the SRF parameter will be obtained by replacing all STF variables in the following equation system (72): As expected, in the case of the ICT method, the graphical representations of the STF parameter ( Figure 11, scenario a) and the SRF parameter ( Figure 11, scenario b) are very similar to the linear distribution of the OBIST technique, both being hardware-oriented testing strategies.

Fault Coverage-Aware Metrics Comparison with Other Related Works
Regarding the novelty of our research, it is worth mentioning that only two state-ofthe-art works [14,19] consider failure models as a relevant impact factor on the reliability determination of PV systems. A summarized comparison between the author's approach and our proposed metrics can be seen in Table 6. First, it is important to outline the fact that the author's studies from [14,19] are focused on failures occurring in static deployed PV systems where the tested equipment is limited to the PV panel and the extra components which are used for solar energy conversion. Our work compensates for this deficiency by taking into account the number of errors that can occur in mobile PV systems, where advanced electrical equipment is used to steer the mechanical parts of the solar tracking system in different directions. Secondly, according to the table data in the Weibull distribution model [14], the actual number of failures of the tested components are determined without knowing the nature of the deterioration, for which the fault-coverage aware metrics offer additional details such as the cause factor given by the existence of software, hardware and ICT errors in the automation equipment. Thirdly, by analyzing the last two columns of Table 6, it is observable that the availability/reliability versus time factor allows us to establish an accurate prediction of the PV system performance for long-term and short-term usage. Since the proposed metrics are applied to testing techniques that are generating the results through simulation, the fault coverage is determined in a reduced time manner. Thus, the SRF parameter describes the performance of the entire solar tracking system for short-term usage, without the ability to make long-term predictions of its real-time behavior. To overcome this issue, as future work we plan to fuse our software, hardware, and ICT testing strategies in a combined testing suite that will be physically deployed to verify the functionality of the automated solar tracking equipment in real-time over several weeks. The final scope is to increase the global reliability factor of the solar tracking system by making use of the proposed fault coverage aware metrics presented in this paper.
Finally, as can be seen in Table 7, this paper presents a comparison between the proposed fault coverage-aware metrics and the Weibull distribution model regarding calculation steps and execution time. This comparison was realized using the Python programming environment. It can be seen that the proposed fault coverage-aware metrics significantly reduce the computation time by 85.91% and the number of calculation steps by 74.81%, showing that our metrics are considerably more efficient when compared to the standard Weibull distribution model.

Conclusions
This paper presents two different mathematical approaches for calculating the reliability of solar tracking systems, which are a particular case of automation systems. For this, first, the general model of the Weibull distribution model is described as a collection of standard metrics such as MTF, MTBF, MDT, and the famous hazard rate, which helps in shaping a probabilistic methodology for calculating the reliability and availability of solar tracking systems.
The theoretical background was completed with concrete examples from real-life scenarios, where individual components were tested to find out their failure frequency. Accordantly, it was observed that today's standard metrics are first, excluding the nature of a component's defects, thus being limited in determining the failure moments of components in automation systems; secondly, to establish an accurate prediction of the reliability and availability, the testing time of standard metrics can vary from days, weeks, months, to even years. To overcome the aforementioned limitations of standard metrics used in the literature, this paper proposes, to the best of our knowledge, for the first time in literature, a set of three novel fault coverage-aware metrics that make use of experimental data gathered from hardware, software as well as ICT-oriented testing techniques. While our metrics are most suitable for BIST architectures, WBST routines, and ICT methods, we firmly believe that they can be potentially extended for several methodologies in the hardware, software, and ICT testing domains.
Finally, when compared to the Weibull distribution model for determining the reliability of solar tracking systems, due to the use of two parameters called STF and SRF, our fault coverage-aware metrics require significantly fewer calculation steps and reduced computation time.