Application of Time-Voltage Characteristics in Overcurrent Scheme to Reduce Arc-Flash Incident Energy for Safety and Reliability of Microgrid Protection

The interconnection between diverse Distribution Generations (DGs) that utilize various technologies and complex structure of networks are the most characteristic of modern Distribution Networks (DN). The wide adoption of DGs considerably affects the power flow dynamics in the DN and consequently the fault characteristics. The excessive level of fault currents can pose risks of heat (high temperature) and pressure in accordance to Arc Flash (AF) incident energy in microgrids. This research studies the relationship between AF severity and the solving of coordination problem of Overcurrent Relays (OCRs) in DN, and introduces a novel equation that considers the AF qualities in solving the coordination problem for OCRs. In this study, a novel optimization problem, the AF severity with the optimal coordination of OCRs in DN is presented and the Water Cycle Optimization Method (WCOM) is employed to find the best combination of the OCR’s settings in the DN while considering the AF induced energy. The proposed optimization approach and the novel equation are evaluated with an IEC microgrid and compared with the conventional protection method and Particle Swarm Optimization (PSO) used in optimizing the coordination of OCR in the DN. The optimal settings of the OCR scheme are achieved and examined on the modified IEC microgrid benchmark system. In order to verify the result, an industrial simulation package (ETAP) and OCR (GE Multiin, model-750/760) was used in this work.


Background
Currently, due to their extensive inclination to integrate more sustainable and clean energy sources, together with the development of transmission and distribution technologies, DNs have turned into groups partially independent of the networks [1]. The most straightforward possible description of microgrids is that of a group of DGs connected directly to the load areas. It may be fully or partially fed from the main utility through an intelligent management system [2]. A microgrid's unique nature poses challenges that do not appear in conventional distribution networks [3]. Nevertheless, the need for more reliable, efficient, flexible, and cost-effective networks reinforces the trend of transforming DNs into active networks [4][5][6][7]. One of the challenges that microgrids bring, which are one the most important fields of research nowadays, is the coordination of protection schemes to be more responsive and reliable in case of faults [8][9][10][11]. The renewable energy resources, 2 of 19 or DGs, have converted the conventional DN from a radial system into a multi-looped system. As a result of the fact that workers perform their tasks on charged conductors, workers should have Personal Protective Equipment (PPE) to handle the high degree of AF incident energy [12]. The importance of such measures is demonstrated by the ten-year study conducted by Electricite de France on 120,000 workers, which showed that electrical arc wounds caused more than 77% of the injuries [13]. The excessive energy generated by integrating more energy sources into the network than the equipment's operational limits, may cause equipment damage and longer repair time in addition to costly repair expenses. Such faults shorten the life of the equipment and also limit the capabilities of protection schemes. The ability to operate the DN with DGs with the aim of limiting the impact of AF in the DN will increase the selectivity and protection coordination complexity. Therefore, selecting the optimal protection scheme is becoming one of the most challenging and important tasks to safely and securely operate the DN with DG systems. In this article, the fault current, AF incident energy characteristics and the OCRs coordination challenges in DN equipped with DGs can be described as a novel optimization problem. The solving of the novel optimization problem will help to find and adjust the OCR's settings by considering the AF qualities in order to minimize the total operation time for OCRs and increase the reliability of the DN.

Literature Review
Evaluation of arc risks is one of the most critical factors that needs to be estimated and assessed in case of faults to reduce the amount of energy produced by AF. According to the reports and statistics provided by national associations (NFPA [14] and OSHA [15]), there are 30,000 cases of injury as a result of AF incidents annually in the US. There were also 7000 burn injuries caused by AFs, of which 2000 were admitted to hospitals annually. In addition, 400 deaths annually as a result of injuries from incidents were also counted. According to NFPA, the evaluation and analysis of arc risks and the labeling of AF is one of the necessary prerequisites for selecting the appropriate protection devices when operators work on a charged cable within the boundaries of AF safety zones.
Various solutions have been proposed to reduce AF risk. Among these solutions, there are the proposal to introduce arc-resistant switchgear, subsidiary breaker management assembly, arc reduction switch, and a relay with optical sensor for lights during arc incidence. Such methods are criticized for the expensive cost of upgrading the equipment in sub-stations. Further, several authors have suggested changing the setting of the relays to minimize the energy generated by AF incidents through minimizing the overall operation time of relays [16,17]. The authors in [18] suggested merging OCR coordination considering the AF. In this study [18], the minimum operational time is obtained while preserving an adequate coordination time interval. In another study, the FPA-based method was used to obtain the best possible solution in terms of total operation time for the overcurrent relays [19].
A number of studies have discussed and presented AF problem for DNs. For example, the authors in [20] studied the relationship between cables length and AF. While in [21], the effect of electrode geometry on the limits of the protection of AF were studied, and the author recommended that this effect should be increased according to IEEE 1584-2018 specifications. The effect of AF on the fault calculations via IEC and ANSI were studied in [22], whereas the authors in [23] compared both methodologies integrated and momentary method to improve AF calculation. A simplified hazards analysis of arc-flash using boundary curves were conducted in [24]. In [25,26], current limiting fuses were leveraged to minimize Arc Flash Incident Energy (AFIE), while in [27,28] researchers developed an approach that can detect and mitigate the AF at the speed of light and sound. An MV high speed grounding switch was used for AF mitigation in [29,30]. In [31], the authors improved the protection level in relays to detect AF. In [32,33], researchers used a graphical approach to analyze accident energy and knock, to discover the maximum and minimum value of AFIE. The study in [34] introduces a simplified arc-flash analysis methodology to assess the settings of over-current protective devices and arc rated PPE. In [15], the author suggested taking into account the anti-islanding schemes actions and procedures to eliminate AF as well as focusing on the interconnect transformer, as it increases the percentage of zero sequence sources and temporarily de-sensitizes the feeder protective relay. Most approaches for the OCR coordination problem in relevant research are formulated as a constrained optimization problem [20][21][22][23][24][25][26][27]. However, no study to date has considered AF in the optimization problem with OCRs. In Saad et al. [35], the primary investigation, which used OCRs to minimize the AFIE, was presented by employing the water cycle algorithm. However, the proposed approach in [35] did not consider the benefits of employing definite regions to minimize the total operation time for all OCRs and to minimize the AFIE. The paper [35], did not present and investigate the performance of the proposed approach with a different fault level or fault location. In addition, the authors in [35] did not investigate and compare the performance of the WCOM to common and powerful algorithm such as PSO. This article aims to develop and employ powerful optimization algorithms (WCOM and PSO) due to the limited number of studies that have considered the significance of using these new metaheuristic algorithms. In this paper, the investigation of the benefits of using the proposed optimization methods to achieve the minimum operating time for all relays (primary and backup) and the minimum level of the AFIE in the network, while maintaining an appropriate coordination time (CTI), are extended and presented. The proposed new optimization problem aims to guarantee the mixed coordination between OCRs and AF in a way will not affect the traditional function of OCRs by using WCOM and by being compared to PSO. This paper is one of the first extensive investigations looking into the best combination of the OCR's settings in the DN while considering the AF induced energy for different operations and fault scenarios by using the WCOM and PSO algorithms.

Contribution of the Paper
The current of the AF in DN systems varies considerably depending on the fault current type and the DN configuration. The variation being mostly due to the impedance of the arc itself. The AF analysis can be used in two main purposes to help identify specifications with respect to relay settings to minimize AFIE levels within the equipment capability at different network topologies. Although overcurrent relays need to be coordinated to have sufficient Coordination Time Interval (CTI) between main and backup pairs, AF analysis necessitates that the total clearing time must be as minimal as possible to minimize the intensity of the AFIE. Given this, most OCR coordination procedures with inverse Time-Current Characteristics (TCC)s are in compliance with the AF analysis to simplify acquiring the optimal setting considering arc-flash analysis in DN systems The contributions of this paper are listed below: • A novel optimization problem to present the AF severity with the optimal coordination of OCRs in DN. In this problem, additional constraints are suggested to be placed on the optimization problem in order to reduce AF incident energy. Moreover, this paper suggests a new constraint to the optimization problem to comply with the operational limitation of industrial relays. In this work, a defined region was the optimization problem with the AFIE categories. This aims to find the optimum setting of OCRs that considers the AFIE categories; definite regions and DGs would provide a significantly more secure environment for workers and equipment. The proposed new optimization problem aims to guarantee the mixed coordination between OCRs and AF in away will not affect the traditional function of OCRs by using WCOM, as described in following subsection.

•
A new optimal coordination scheme of OCRs in DN that incorporates the DGs and the AF qualities for different scenarios is presented and solved by a new optimization method (WCOM) and compared to the powerful algorithm PSO.

•
The tripping time of OCR installed on the DG bus is discussed by presenting a new coordination strategy (non-standard) for the OCR model. The proposed relay A comparative investigation is performed between the standard and non-standard techniques for OCR coordination, which does not include the AF, and the approach suggested considers the AF within the OCR coordination optimization model under different fault scenarios. This will provide the DN operators and engineers an indicator about the impact of AF on the protection and operation of power systems.

•
The proposed OCR scheme approach is tested and verified by the use of ETAP industrial software, to provide sufficient analysis and prove the robustness of the new approach.

Outline of Paper
This paper is outlined as follows: Section 2 introduces the problem statement and the mathematical model for the AF protection problem. Sections 3 and 4 present the novel optimal protection coordination scheme for the OCRs. The case study and the analysis of results are discussed in Section 5. The summary and conclusions of this paper are presented and discussed in Section 6.

Problem Description: Arc-Flash Protection
Arcing short-circuits in LV and MV switchgear are one of the hazard places in the event of a fault, which causes several operational hazards and damage to equipment in addition to individual casualties. A somewhat considerable tripping time of the overcurrent scheme commonly used in DNs may cause the perpetuation of the arcing energy. Therefore, reducing the tripping time, even to a few hundreds of milliseconds, will maintain the safety requirements in sub-stations [36]. The inverse TCCs and AF categories shown in Figure 1 are an indicator to identify AF hazard categories. Figure 1 shows the relationship between AFIE and OCRs for a DN without DG. If the AF categories are not considered when coordinating overcurrent relays, then category 0 is basically located under the fault defined region in the TCCs. Therefore, it is possible that when a fault occurs, the AFIE will be higher than the incident energy assigned to category 0 and thus will increase the risk of exposing equipment and operators to risk.

•
A new optimal coordination scheme of OCRs in DN that incorporates the DGs and the AF qualities for different scenarios is presented and solved by a new optimization method (WCOM) and compared to the powerful algorithm PSO.

•
The tripping time of OCR installed on the DG bus is discussed by presenting a new coordination strategy (non-standard) for the OCR model. The proposed relay characteristic is compared with the traditional TCC in terms of the minimum operation time for a notable drop in the bus voltage. • A comparative investigation is performed between the standard and non-standard techniques for OCR coordination, which does not include the AF, and the approach suggested considers the AF within the OCR coordination optimization model under different fault scenarios. This will provide the DN operators and engineers an indicator about the impact of AF on the protection and operation of power systems.

•
The proposed OCR scheme approach is tested and verified by the use of ETAP industrial software, to provide sufficient analysis and prove the robustness of the new approach.

Outline of Paper
This paper is outlined as follows: Section 2 introduces the problem statement and the mathematical model for the AF protection problem. Sections 3 and 4 present the novel optimal protection coordination scheme for the OCRs. The case study and the analysis of results are discussed in Section 5. The summary and conclusions of this paper are presented and discussed in Section 6.

Problem Description: Arc-Flash Protection
Arcing short-circuits in LV and MV switchgear are one of the hazard places in the event of a fault, which causes several operational hazards and damage to equipment in addition to individual casualties. A somewhat considerable tripping time of the overcurrent scheme commonly used in DNs may cause the perpetuation of the arcing energy. Therefore, reducing the tripping time, even to a few hundreds of milliseconds, will maintain the safety requirements in sub-stations [36]. The inverse TCCs and AF categories shown in Figure 1 are an indicator to identify AF hazard categories. Figure 1 shows the relationship between AFIE and OCRs for a DN without DG. If the AF categories are not considered when coordinating overcurrent relays, then category 0 is basically located under the fault defined region in the TCCs. Therefore, it is possible that when a fault occurs, the AFIE will be higher than the incident energy assigned to category 0 and thus will increase the risk of exposing equipment and operators to risk.  The fault characteristics and tripping time are different in the DN with and without DG systems. In addition, the AF Category increases the complexity of choosing the optimal OCR setting. For DNs equipped with DGs, the category 0 will be updated to a higher level than the OCR characteristics because of the increased fault current compared with DNs without DGs. Therefore, the category of AFIE will change from 0 to 1, which mean higher risks compared to DNs without DG. In conventional DNs, which do not connect to DGs in the distribution level, the AF incident energy category is less than the inverse TCCs of the designated OCR. However, with the interconnection of DGs, as in Figure 1, the AF incident energy category will increase from the category of 0 to 1, which requires wearing PPE that matches the energy generated by the electric arc in the event of a fault. Thus, coordinating the TCCs of OCRs and the AFIE category is an urgent necessity to ensure that the AFIE category resulting from the fault remains below the fault point on the TCC of OCRs as in Figure 1. For example, the DN with DG will have a fault current where category 0 is located Above the OCR curve, and therefore it was considered in the calculation as shown in Figure 1. However, when DG is connected to the grid, the fault current may change the AFIE category from 0 to 1, which makes the situation more dangerous compared to the DN without DG (as shown in Figure 1).
As indicated by Figure 1, when DG is connected, the OCR may enter the definite region, consequently increasing the AFIE category from 0 to 1 and 2. This introduces higher risk than the previous operation conditions. Therefore, the DN will require having more secure PPEs and clothes. Therefore, finding the optimum settings of OCRs that consider the AFIE categories, definite regions, and DGs would provide a significantly more secure environment for workers and equipment. Several studies have used a light sensor (energy sensor) or differential protection to detect the occurrence of AF [36][37][38][39]. The occurrence of AF with circuit breakers inside sub-stations makes it possible to use the existing voltage transformers inside those substations to obtain voltage readings, which can be implemented in the inverse TCC of the overcurrent relay to obtain the least possible time in case of fault, as well as the possibility of adjusting the TCC within the numerical OCRs that are commonly installed and utilized in the existing DNs. Thus, the total operational time of OCRs will depend on the voltage and current by introducing the possibility of interpreting voltage measurements in the characteristic's equation of the overcurrent relay without the need for changing the current protection group settings.

Calculation of the AFIE
Many factors, such as the magnitude of arcing current, the total clearing time for fault events, and the dimensions of the working area, significantly affect the AFIE. Therefore, minimizing the tripping time of OCRs in line with the coordination between the relays (primary and backup) may be one of the most influential and feasible solutions to reduce the level of AFIE. The mathematical expression in Equation (1) shows the relationship between the AFIE at the specific working dimensions and the total clearing time as follows [19]: where E n is the normalized AFIE over time period, C f is the estimation variable (1.5 for LV and 1 for HV), t is the total tripping time, W D is the operating distance, and x is an exponent factor [19,39]. Thus, AFIE can be measured for a voltage level up to 15 kV by applying the following equation: where the z 1 and z 2 are configurations based on the grounded systems, I arc is the arcing current, and G p is the gap distance between conductors. Equation (3) describes the arcing current for a voltage level higher than 15 kV.
I f is three-phase short circuit current in kA. As can be seen from Figure 1, it is related to the AFIE category. The incident energy can be higher in active DNs with high DG-share, based on the network topology and nature of the fault. This alone highlights the need to revise the OCR coordination strategy taking into account the arc-flash analysis.

Optimal Coordination Problem of OCRs and AFIE
The proposed model in this article aims to minimize AFIE. In this work, with the help of AFIE considerations achieved from AF analysis, OCR coordination will reduce the overall operational time while maintaining coordination between the relays (primary and successive backup). More operational details for coordinating OCR are found in [40,41]. In the following subsection, the objective function that describes the AFIE problem and the optimization model constrains are presented and discussed.

Formulation of OCR Coordination Problem
The constrained coordination problem (objective function, OF) of OCR in DNs can be represented in the following Equation (4): where the number of OCRs form j = 1 to n, t j,k is the tripping time of relay j at short-circuit current location k. The weight that defines the chance of event of a fault on a selected line is presented by W. Here, the weight value is assumed to be one [19,41], and the fault incident has an equal probability on each line. The total tripping time of the OCRs is limited by the relays and network operation constraints as regards.

Coordination Criteria
The selection of CTI is fundamentally depended on the tripping time for the primary OCR, overshooting time, and the clearing time of breaker. The time coordination constraints can be demonstrated in Equation (5): where t 2,k and t 1,k are the tripping time of backup and primary OCR for a short-circuit current at location k, respectively. The ∆t is the CTI, which is regularly produced in range of 0.2-0.5 s [42] and it set as 0.3 s in this study.

Relay Operating Time Constraints
The tripping time of overcurrent relays is expressed as a constrained variable to obtain the limits for the highest and lowest tripping time, t max and t min , respectively, for OCR j as follows:

Proposed Current Multiplying Setting (CMS) of Industrial Relay Characteristics
This paper also introduces the addition of the operational limitation constraint of industrial relay characteristics to the formulation of inverse OCR coordination. The highest limit for the inverse characteristics of OCR, CMS max , was defined as high as 100, while the lowest possible limit, CMS min , was chosen as 1.1. The following equation describes this relation: 3.

The TMS and PS Constraints
TMS is the primary variable in managing the OCR scheme coordination time and PS is the plug setting for the OCR. Accordingly, the value of the TMS and PS need to be between the lowest and maximum value of TMS and PS, TMS min , PS min , TMS max and PS max , respectively, as following: 3.1.5. The Relay Constraints Considering AFIE In this paper, a new objective function for OCRs considering AFIE is formulated by creating a voltage-current-based time inverse for the OCR model. The reduction in total operational time of this suggested is based on a reduction in the bus voltage when fault occurs near DG's bus-bar when DG connected with the distribution system to develop relay variables including TMS, pickup current, and OCR parameters. OCR scheme is essentially employed to protect low-voltage DNs. The highest and lowest tripping time are determined by calculating the operation time of OCRs, which is complex and non-linear equation and can be expressed for Standard Inverse (SI) characteristic as follows: In order to minimize the tripping time of the OCR scheme, an exponential term was added to the traditional SI equation, and thus less the AFIE. The new OCR equation after modification become as follows: where v j,k is pre-phase fault voltage measured in per unit on OCR bus voltage for relay j and fault location k. q is a fixed variable and is selected as two [19,40].

Addressing Water Cycle Optimization Method (WCOM) in Solving the Complex Overcurrent Relays Coordination Optimization Problem
The WCOM is an innovative heuristic optimization algorithm that is inspired by nature from the water cycle in seas and rivers [43]. In general, The WCOM is a powerful optimization solver as presented in [44] and it can be developed by The Optimization Toolbox in 295 MATLAB/ SIMULINK (MathWorks, Inc., Natica, MA, USA) [45]. This article aims to develop and employ a new powerful optimization algorithm (WCOM) as limited studies have considered the significant of using the new metaheuristic algorithms compared to common algorithms such as PSO. In this paper, the proposed new optimization method (WCOM) aims to achieve the minimum operating time for all relays (primary and backup) and the minimum level of the AFIE in the network while maintaining an appropriate coordination time (CTI). The proposed new optimization problem aims to guarantee the mixed coordination between OCRs and AF in away will not affect the traditional function of OCRs by using WCOM and compared to PSO, as described in following subsection. The new optimization algorithm (WCOM) and common algorithm PSO [10] have been used and employed in this work for solving different operation mode in microgrid. WCOM is mainly carried out through an iteration process. Firstly, WCOM will initially have random parameters and solution for solving the optimization problem in Equation (4). Then, the solution position will be updated based on the best results in previous step. In the next step, the swapping between the solutions locations and importing random new solution aims to achieve the optimum solution and avoid local optimal solutions. Finally, the WCOM algorithm will achieve the best possible solution at the end of the iteration process. The following steps show the main procedures of WCOM: Step 1: set the main parameter of WCOM: the main parameters of WCOM have been determined by testing each of the parameters over a wide range of values, as shown in  Table 1, then the optimal (best) parameter value was chosen to obtain the simulation results in this article [46].
Step 2: generating a random initial population of solutions (streams).
Step 3: determining the number of the number of streams and the decision variables.
Step 4: calculating the objective function value for each possible solution by using the overcurrent optimization problem (Equation (4)).
Step 5: updating the solution position and selecting the new location for the best results in previous step.
Step 6: swapping between the solutions locations to achieve the optimum solution and avoid local optimal solutions.
Step 7: adding a new random solution to avoid local and saturation in choosing the optimal solution.
Step 8: the WCOM is an iteration algorithm, the above process from step 2 is repeated until the maximum number of iterations are achieved.

Illustration of the Proposed Strategy to Reduce Arc-Flash Energy in Over Current Protection Scheme
An efficient arc-flash reduction strategy can significantly reduce the energy level, reduce the possibility of electrical hazard risk accident, decrease equipment damage, and maintain electrical continuity. During an AF event, the generated energy is dependent on voltage, current, and duration of the arc event. Some proposed options recommend reducing system voltages or fault current level through a grounding system. Still, the most effective and least expensive strategy is to minimize the fault-clearing time by using adaptive setting groups in numerical relays and utilizing wireless communications already available in DNs to avoid harming operators and technicians.
This article aims to present a strategy for addressing the coordination problems of protective relays, which include assessing AF in the grid. In the proposed approach, the minimum operating time for all relays (primary and backup) were obtained. As a result, the AFIE levels were reduced in the network while maintaining an appropriate coordination time (CTI). The proposed coordination strategy (voltage OCR with nonstandard scheme) is described in Figure 2 and shows the standard inverse TCCs of overcurrent relays when a fault occurs, and their relation to the AF Categories. In the event of a fault (I) where the AF categories not taken into account, the AF category level will increase from level 0 to level 1 and thus poses a risk to the life of the operators and damage to the equipment. The employment of proposed optimization scheme (POSH) will ensure that the AF category does not increase from level 0 to 1 by providing less tripping time. As a result, the POSH characteristics will be located under the lowest AF category, which is a less risky situation for operators and equipment. In integrating DGs into the distribution network, the possibility of increasing the AF categories will become more significant. When fault (I) occurs and AF analysis is not considered in overcurrent coordination, the AF event will become more dangerous since the fault current enters the relay's definite region. As a result, the AF category of the fault will change from level 1 to level 2. This situation requires wearing safer protective clothes. Consequently, consideration of AF assessment in overcurrent coordination is essential to ensure operators' safety and avoid equipment damage. When POSH characteristics are utilized, these characteristics will ensure a short tripping time so the fault current will not reach the definite region of the relay. Therefore, IE will not increase so high that it may cause losses in lives and equipment.

Work-Flow and Procedures of the Proposed Strategy
In this work, an optimal solution approach for OC relays coordination problem considering AF events was developed and compared to the traditional coordination scheme. The WCOM as described in previous section is used to find the optimal TMS values for all OC relays considering AF events [41]. The proposed approach in this work is tested by using two DN scenarios. The first DN operation scenario (Mode 1). The second scenario (Mode 2) is a microgrid fed by the utility and with DGs. The proposed optimal solution approach aims to adjust the TMS value for all the OC relays and the relays coordination is selected by considering AF events as shown in the work-flow. Figure 3 presents the main steps involved in the proposed novel approach for acquiring the optimal coordination of overcurrent relays. In the beginning, the ETAP program was used to implement the IEC microgrid benchmark by specifying and defining the data of each of the transmission lines, Current Transformer Ratio (CTR), and network topology. Execution for the simulated system's load flow was made to initially identify the proper CTRs and adjust the relay settings. Then, three-phase fault simulations were carried out at different locations in the network in the two simulated modes (Mode 1 and Mode 2). The acquired results through load flow and short-circuit calculations in ETAP were exported to MATLAB/SIMULINK software to solve the optimization problem. The inputs and settings for WCOM were specified to create the population based on the upper and lower limits and the constraints designated in overcurrent coordination problems. Computation for the objective function in each raindrop was conducted to check the convergence of the optimization problem. The extracted optimal setting is inserted into the ETAP program to implement and evaluate the obtained overcurrent relay setting. The calculation of Incident Energy is based on IEEE1584 while AFIE level is determined by NEPA 70E 2009. This work presents a comparative analysis for the proposed inverse time In integrating DGs into the distribution network, the possibility of increasing the AF categories will become more significant. When fault (I) occurs and AF analysis is not considered in overcurrent coordination, the AF event will become more dangerous since the fault current enters the relay's definite region. As a result, the AF category of the fault will change from level 1 to level 2. This situation requires wearing safer protective clothes. Consequently, consideration of AF assessment in overcurrent coordination is essential to ensure operators' safety and avoid equipment damage. When POSH characteristics are utilized, these characteristics will ensure a short tripping time so the fault current will not reach the definite region of the relay. Therefore, IE will not increase so high that it may cause losses in lives and equipment.

Work-Flow and Procedures of the Proposed Strategy
In this work, an optimal solution approach for OC relays coordination problem considering AF events was developed and compared to the traditional coordination scheme. The WCOM as described in previous section is used to find the optimal TMS values for all OC relays considering AF events [41]. The proposed approach in this work is tested by using two DN scenarios. The first DN operation scenario (Mode 1). The second scenario (Mode 2) is a microgrid fed by the utility and with DGs. The proposed optimal solution approach aims to adjust the TMS value for all the OC relays and the relays coordination is selected by considering AF events as shown in the work-flow. Figure 3 presents the main steps involved in the proposed novel approach for acquiring the optimal coordination of overcurrent relays. In the beginning, the ETAP program was used to implement the IEC microgrid benchmark by specifying and defining the data of each of the transmission lines, Current Transformer Ratio (CTR), and network topology. Execution for the simulated system's load flow was made to initially identify the proper CTRs and adjust the relay settings. Then, three-phase fault simulations were carried out at different locations in the network in the two simulated modes (Mode 1 and Mode 2). The acquired results through load flow and short-circuit calculations in ETAP were exported to MATLAB/SIMULINK software to solve the optimization problem. The inputs and settings for WCOM were specified to create the population based on the upper and lower limits and the constraints designated in overcurrent coordination problems. Computation for the objective function in each raindrop was conducted to check the convergence of the optimization problem. The extracted optimal setting is inserted into the ETAP program to implement and evaluate the obtained overcurrent relay setting. The calculation of Incident Energy is based on IEEE1584 while AFIE level is determined by NEPA 70E 2009. This work presents a comparative analysis for the proposed inverse time OCR approach (POSH) with the conventional Standard Inverse Time Relay (SITR) coordination approach in terms of CTI and overall operational time and overall incident energy.

Simulation Results and Discussion
The presented formulation of the AF severity with the optimal coordination for over current relays in Sections 2 and 3 is evaluated and examined using a modified IEC microgrid with two DGs. This section aims to discussed and show the effectiveness of the proposed new optimal approach (POSH) for OCRs in DNs that incorporate the DGs and the AF qualities for all possible three-phase faults scenarios. In the following subsection, the description of the case study considering the network specifications is presented. Then, the results of the proposed scheme of OCRs and the AF in microgrid are presented and discussed under different IEC network scenarios. Finally, the proposed optimal coordination scheme is tested using the industrial software (ETAP). The case study and the analysis of results are described and discussed in this section.

Case Study
In this paper, the main case study is presented in Figure 4. The proposed case study is a modified IEC microgrid model including two DGs. The base power in the proposed case study has been selected to be 10 MVA and the detailed information of the studied DN is given as by [46]. The case study network contains six OCRs, as presented in Figure 4. Table 2 presents the Plug Setting (PS), Current Transformer Ratio (CTR), and Pickup current (IPP) for each OCR. In order to evaluate the proposed solution, three-phase fault has been simulated at different locations. In addition, this work aims to investigate and present the performance of the proposed solution under different network operation scenarios (with and without DGs). The microgrid system is first considered as a radial network without DG and is protected by the OCR scheme with time-current inverse characteristics for the implementation of the procedure. However, the fault contribution is shared between the utility source and DG when DG is added to the system. The power flow is bidirectional, reducing the capability of the OCR scheme. In order to evaluate the advantages and benefits of the proposed approach in terms of providing minimum Overall Operational Time (OOT) for different fault locations, as shown in Figure 4, with respect to the CTI between primary and secondary OCR pairs, the CMS, fault current, and CMS in different fault locations are presented in Tables 3 and 4.

Simulation Results and Discussion
The presented formulation of the AF severity with the optimal coordination for over current relays in Sections 2 and 3 is evaluated and examined using a modified IEC microgrid with two DGs. This section aims to discussed and show the effectiveness of the proposed new optimal approach (POSH) for OCRs in DNs that incorporate the DGs and the AF qualities for all possible three-phase faults scenarios. In the following subsection, the description of the case study considering the network specifications is presented. Then, the results of the proposed scheme of OCRs and the AF in microgrid are presented and discussed under different IEC network scenarios. Finally, the proposed optimal coordination scheme is tested using the industrial software (ETAP). The case study and the analysis of results are described and discussed in this section.

Case Study
In this paper, the main case study is presented in Figure 4. The proposed case study is a modified IEC microgrid model including two DGs. The base power in the proposed case study has been selected to be 10 MVA and the detailed information of the studied DN is given as by [46]. The case study network contains six OCRs, as presented in Figure  4. Table 2 presents the Plug Setting (PS), Current Transformer Ratio (CTR), and Pickup current (IPP) for each OCR. In order to evaluate the proposed solution, three-phase fault has been simulated at different locations. In addition, this work aims to investigate and present the performance of the proposed solution under different network operation scenarios (with and without DGs). The microgrid system is first considered as a radial network without DG and is protected by the OCR scheme with time-current inverse characteristics for the implementation of the procedure. However, the fault contribution is shared between the utility source and DG when DG is added to the system. The power flow is bidirectional, reducing the capability of the OCR scheme. In order to evaluate the advantages and benefits of the proposed approach in terms of providing minimum Overall Operational Time (OOT) for different fault locations, as shown in Figure 4, with respect to the CTI between primary and secondary OCR pairs, the CMS, fault current, and CMS in different fault locations are presented in Tables 3 and 4.

Results Analysis and Discussion
The results of the proposed POSH and conventional SITR schemes are presented over different network operation scenarios. In this subsection, we will compare the POSH and SITR schemes performance using WCOM optimization algorithm for a DN based on the following configurations: • Mode 1: the utility fed the microgrid (without DG). • Mode 2: the utility and DGs fed the microgrid (with DG).

Numerical Results for the Proposed POSH and Conventional SITR Schemes
The proposed optimal OCR (POSH) and the conventional SITR schemes are evaluated and compared in terms of CTI, overall relays operational time and overall incident energy using WCOM optimization algorithm are evaluated in this section. Firstly, Table 5 presents the OOT and TMS value for all relays under the three-phase fault condition. In addition, the results were obtained for both network operation scenarios (Mode 1 and Mode 2). The OOT was significantly decreased by using the proposed POSH approach compared to the conventional SITR approach, as presented in Table 5. For example, the OOT decreased from 5.68 S (SITR) to 2.565 S (POSH) in Mode 1 scenario. The OOT increased in Mode 2 compared to Mode 1 for both POSH and SITR approaches. This indicates that the DGs at the network increases the complexity of the coordination problem and OOT for OCRs. However, the POSH approach outperformed the SITR for both modes, where the POSH decreased the OOT by 57.2% and 54.8% at Mode 2 and Mode1, respectively. Figure 5 presents the percentage improvement of using the PITR and compared to SITR schemes in term of OOT for Modes 1 and 2. OOT decreased from 5.68 S (SITR) to 2.565 S (POSH) in Mode 1 scenario. The OOT increased in Mode 2 compared to Mode 1 for both POSH and SITR approaches. This indicates that the DGs at the network increases the complexity of the coordination problem and OOT for OCRs. However, the POSH approach outperformed the SITR for both modes, where the POSH decreased the OOT by 57.2% and 54.8% at Mode2 and Mode1, respectively. Figure 5 presents the percentage improvement of using the PITR and compared to SITR schemes in term of OOT for Modes 1 and 2.  In this work, to evaluate the performance of the proposed POSH and conventional SITR approaches on the clearing of arc faults, Tables 6 and 7 are presented. The main comparison terms for the dealing with arc faults are the clearing time for arc fault (FCT, s), the incident energy of arc (AFIE, Cal/cm 2 ), the level of energy producing during fault (EL), and the boundary of arc flash (AFB, ft), as shown presented in Table 6 for Mode 1 and Table 7 for Mode 2. Firstly, the EL significantly decreased by using the proposed POSH approach compared to the conventional SITR approach. For example, the EL decreased from level 3 (SITR) to level 1 (POSH) in the Mode 1 scenario at relays R2 and R6. Table 7 showed that the EL decreased from level 4 (SITR) to level 0 (POSH) in the Mode 2 scenario at relays R2 and R6. This practice will help to presents an effective solution to deal with the arc fault energy in DNs with DGs and significant enhancement in safety. Secondly, the clearing time for arc fault FCT was significantly decreased by using the proposed POSH approach, as shown in Table 6 for Mode 1 and Table 7 for Mode 2. For example, the FCT decreased from 1.17 S and 0.163 S (SITR) to 0.18 S and 0.02 S (POSH) in In this work, to evaluate the performance of the proposed POSH and conventional SITR approaches on the clearing of arc faults, Tables 6 and 7 are presented. The main comparison terms for the dealing with arc faults are the clearing time for arc fault (FCT, s), the incident energy of arc (AFIE, Cal/cm 2 ), the level of energy producing during fault (EL), and the boundary of arc flash (AFB, ft), as shown presented in Table 6 for Mode 1 and Table 7 for Mode 2. Firstly, the EL significantly decreased by using the proposed POSH approach compared to the conventional SITR approach. For example, the EL decreased from level 3 (SITR) to level 1 (POSH) in the Mode 1 scenario at relays R2 and R6. Table 7 showed that the EL decreased from level 4 (SITR) to level 0 (POSH) in the Mode 2 scenario at relays R2 and R6. This practice will help to presents an effective solution to deal with the arc fault energy in DNs with DGs and significant enhancement in safety. Secondly, the clearing time for arc fault FCT was significantly decreased by using the proposed POSH approach, as shown in Table 6 for Mode 1 and Table 7 for Mode 2. For example, the FCT decreased from 1.17 S and 0.163 S (SITR) to 0.18 S and 0.02 S (POSH) in the Mode 1 scenario for relays R5 and R6, respectively. Finally, the boundary of AF, AFB, and the incident energy of arc (AFIE) decreased in Mode 1 and 2 when using POSH approach compared to SITR. This indicates the ability of the proposed POSH to deal with AF problem and minimize the impact of it. For example, the AFB decreased from 195.5 ft and 12.74 ft (SITR) to 20.03 ft and 3.42 ft (POSH) in Mode 1 scenario for relays R5 and R6, respectively. Its significant to consider the AF analysis in order to solve the OC relays coordination problem in the microgrid network. The reduction in the incident energy of arc (AFIE) shows how the model is applicable to deal with AF problem. Figure 6 presents and compares the AFIE for the proposed POSH and conventional SITR approaches. The results showed that proposed POSH approach considering AF analysis has successfully reduced the AFIE in different fault scenarios and for all relays compared to SITR approach. respectively. Its significant to consider the AF analysis in order to solve the OC relays coordination problem in the microgrid network. The reduction in the incident energy of arc (AFIE) shows how the model is applicable to deal with AF problem. Figure 6 presents and compares the AFIE for the proposed POSH and conventional SITR approaches. The results showed that proposed POSH approach considering AF analysis has successfully reduced the AFIE in different fault scenarios and for all relays compared to SITR approach.

Discussion and Comparison for Different Optimization Algorithms
In this section, we aim to compare the performance of the new WCOM algorithm and the common PSO [10] algorithm under different operational scenarios. The WCOM and PSO have been tested and envaulted with the proposed POSH approach and the conventional SITR approach over the two microgrid operation scenarios (Mode 1 and Mode 2). Firstly, the OOT for all relays in Mode 1 and 2 is presented in Table 8. The WCOM outperformed PSO and obtained the best results of OOT in Mode 1 and 2 by using SITR and POSH approaches. For example, the WCOM achieved an OOT in Mode 2 equal to 6.1 s and 2.613 s compared to 6.21 s and 2.7 s by PSO algorithm for SITR and POSH approaches, respectively. In addition, the compactional cost for the optimization algorithm is presented in Table 8 through the elapsed time term. As shown in Table 6, the WCOM requires less compactional cost to obtain the optimal results compared to the PSO algorithms over all scenarios. The WCOM decreased the elapsed time by 41.8% and 20.8% at Mode 1 and Mode 1 for SITR approach, respectively, compared to PSO algorithm. In order to present the performance of both optimization algorithms (WCOM and PSO) in term of utilization efficiency in CPU, Figure 7 introduces the convergence curves for Mode 1 and 2 with the proposed POSH approach and the conventional SITR approach. The WCOM algorithm outperformed the PSO and achieved higher utilization efficiency in CPU for all cases, as shown in Figure 7. In general, the WCOM curves were more speedy and smoother compared to the PSO algorithm, and the WCOM achieved optimal values for Mode 1 and 2 under the SITR and POSH approaches with the smaller number of iterations.

Evaluation Using the ETAP Software
This section aims to evaluate the proposed POSH and conventional SITR approaches by using powerful industrial software, ETAP software packages. In order to determine and estimate the magnitude of accident energy of AF, a location-specific C-Line was used in ETAP and the normal TCC was considered. As previously described, the specific C-Line indicates that it can be modified within a specific region to minimize the AF energy. The proposed solution in this paper aims to determine and select the optimal C-line configuration with minimum AF energy and relay time operation that leads to optimal OCRs coordination. Figure 8 presents and shows the tripping characteristics and arc energy level outlined by the IEEE Standard 1584TM-2002 for primary OCR (R1). In Figure 8, the lowest line for AFIE category (0 to 4) is presented as the lowest constant energy with the value 1.2 Cal/cm 2 , and this value was obtained by traditional scheme without considering the AF analysis. The proposed OCR curve (POSH) recorded the lowest value of energy with 0.175 Cal/cm 2 . In addition, the SITR approach recorded a 1.12 Cal/cm 2 , which is lower than that of the traditional AFIE category scheme, as shown in Figure 8.

Evaluation Using the ETAP Software
This section aims to evaluate the proposed POSH and conventional SITR approaches by using powerful industrial software, ETAP software packages. In order to determine and estimate the magnitude of accident energy of AF, a location-specific C-Line was used in ETAP and the normal TCC was considered. As previously described, the specific C-Line indicates that it can be modified within a specific region to minimize the AF energy. The proposed solution in this paper aims to determine and select the optimal C-line configuration with minimum AF energy and relay time operation that leads to optimal OCRs coordination. Figure 8 presents and shows the tripping characteristics and arc energy level outlined by the IEEE Standard 1584TM-2002 for primary OCR (R1). In Figure 8, the lowest line for AFIE category (0 to 4) is presented as the lowest constant energy with the value 1.2 Cal/cm 2 , and this value was obtained by traditional scheme without considering the AF analysis. The proposed OCR curve (POSH) recorded the lowest value of energy with 0.175 Cal/cm 2 . In addition, the SITR approach recorded a 1.12 Cal/cm 2 , which is lower than that of the traditional AFIE category scheme, as shown in Figure 8.

Conclusions
In general, AF events generate dangerous level of heat, energy, radiation, and pressure which can cause a human injuries and equipment damage. High-speed arc-flash detection can significantly reduce the energy level and lower the electrical hazard risk, sav-

Conclusions
In general, AF events generate dangerous level of heat, energy, radiation, and pressure which can cause a human injuries and equipment damage. High-speed arc-flash detection can significantly reduce the energy level and lower the electrical hazard risk, saving lives, minimizing equipment damage, and maintaining process continuity. In this work, a complex OC relays coordination problem for a microgrid network with DGs considering the AF events was tested and investigated. This presented two approaches to achieve the minimum overall operation time for all relays in IEC microgrid. The proposed POSH and conventional SITR schemes using the WCOM optimization algorithm were tested and compared over different network operation scenarios. The overall relays operational time and overall incident energy were significantly decreased by using the proposed POSH approach compared to the conventional SITR approach. In addition, ETAP implementation and results show a reduction in arc-flash incidences that demonstrate the effectiveness of the proposed approach. Furthermore, they provide some significant insights into industry representations of prospective AF impact. In this paper, a typical industrial relay is used as an example of how to implement this proposed scheme. Finally, ETAP results quantify the levels to which the proposed approach reduces arc-flash energy and its impact on safety. Certainly, the addition of arc-flash analysis improves overcurrent coordination when compared to the traditional methodologies

Acknowledgments:
We would like to thank the Renewable Energy, Water, and the Environment Center at The Hashemite University for their support.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviation
The