The Effect of Unbalanced Impedance Loads on the Short-Circuit Current

: Conventional short-circuit studies often neglect the load current because the short-circuit current (SCC) ﬂowing from generators is much greater than the SCC that is affected by various loading conditions. As distributed or clustered loads that are unbalanced in phases are connected to the grid, they can also change the magnitude and phase angle of the SCC, despite their small capacities. Thus, the objective of this study is to present algorithms that are able to analyze such an impedance unbalanced load. For this purpose, this study initially derives an SCC model of the unbalanced impedance load in phases. Since the proposed SCC model requires the pre-fault voltage, it uses a power-ﬂow analysis algorithm that iteratively calculates the current that is to be injected and the pre-fault voltage, using the bus impedance matrix. Then, the proposed SCC calculation algorithm transforms the unbalanced loads into equivalent impedances, using the pre-fault voltage, and adds them to sequence networks as input data, using the proposed SCC model. The proposed algorithms are veriﬁed in various case studies. As a result, the proposed SCC calculation algorithms are more accurate, because they do not neglect unbalanced loads.


Introduction
The load current in short-circuit studies is often neglected because the magnitude of the short-circuit current (SCC) flowing from rotating generators is much greater than the magnitude of the SCC that is affected by loads. As distributed or clustered loads that are unbalanced in phases are connected to the grid, they can significantly change the magnitude and phase angle of the SCC, despite their small capacities. Thus, the load current should not be neglected in short-circuit studies. Therefore, many studies have presented various methods for calculating the SCC. For example, in [1], the equivalent circuits of a faulted network, which includes renewable energy sources based on voltage source converters, are proposed. The SCC problem is also solved by the iterative current compensation method with the bus impedance matrix [2]. Since the topology of the distribution system is usually either radial or weakly meshed, the Thevenin equivalent circuit is also solved by the backward and forward sweep (BFSW) method [3]. To synthesize the stochastic variation in the SCC, in [4], the probabilistic simulations (e.g., Monte Carlo simulations) are performed. In [5,6], the decomposition method of the bus impedance matrices to bus-current-injection to branch-current (BIBC), and branch-current to bus-voltage (BCBV) matrices, which were originally proposed for power-flow calculation, are also applied to the SCC calculation. Moreover, to apply the SCC calculation algorithm to the real-time system, a hybrid compensation method that uses the Thevenin equivalent impedance matrix, based on phase representation, is presented in [7].
While the previous studies are based on balanced three-phase systems, the SCC of the multi-phase distribution systems is calculated by representing the multi-phase distribution system as the equivalent three-phase system in [8]. Similarly, when an internal fault occurs in the twelve-phase transmission A distribution system is usually neither balanced nor fully transposed. To the best of our knowledge, none of the previous studies have derived an SCC model of the unbalanced wye-or delta-connected loads, so this study proposes the following equivalent short-circuit models of such an unbalanced load. The proposed SCC model is based on an unbalanced network. Thus, it is applicable to either an unbalanced or untransposed system.

Delta-Connected Unbalanced Impedance Loads
A delta-connected unbalanced impedance load is connected to between phases A and B, in Figure 1. Each line-to-line voltage is as follows: If applying V ABC = AV 012 to (1), the sequence voltage is: If multiplying A −1 to each side, the sequence voltage is as follows: The positive-and negative-sequence components of (3) are as follows: The sequence component of the line current is as follows: Combining Equations (4) and (5) yields the following: Equation (6) is represented by the positive-, negative-, and zero-sequence networks in Figure 2.  (2) If multiplying A −1 to each side, the sequence voltage is as follows: The positive-and negative-sequence components of (3) are as follows: (1 1 120 ) 1 3 The sequence component of the line current is as follows: Combining Equations (4) and (5) yields the following: Equation (6) is represented by the positive-, negative-, and zero-sequence networks in Figure 2.

Wye-Connected Unbalanced Impedance Loads
A wye-connected single-phase impedance load is connected to between phases A and N, in Figure 3. Each phase voltage is as follows: If multiplying A −1 to each side, the sequence voltage is as follows: The positive-and negative-sequence components of (3) are as follows: (1 1 120 ) 1 3 The sequence component of the line current is as follows: Combining Equations (4) and (5) yields the following: Equation (6) is represented by the positive-, negative-, and zero-sequence networks in Figure 2.

Wye-Connected Unbalanced Impedance Loads
A wye-connected single-phase impedance load is connected to between phases A and N, in Figure 3. Each phase voltage is as follows:  A wye-connected single-phase impedance load is connected to between phases A and N, in Figure 3. Each phase voltage is as follows: where The sequence voltage is as follows: The sequence current is as follows: Combining Equations (8) and (9) yields the following: Equation (10) is represented by the positive-, negative-, and zero-sequence networks in Figure 4.
where A = [1 1 1;1 1∠240° 1∠120°;1 1∠120° 1∠240°]. The sequence voltage is as follows: The sequence current is as follows: Combining Equations (8) and (9) yields the following: Equation (10) is represented by the positive-, negative-, and zero-sequence networks in Figure 4.  If an unbalanced wye-or delta-connected impedance load is connected to other phases (e.g., B-N, C-N, B-C, or C-A), the equivalent model in Figures 2 and 4 can be replaced by the corresponding impedance that is connected to other phases. The conventional short-circuit studies have often where A = [1 1 1;1 1∠240° 1∠120°;1 1∠120° 1∠240°]. The sequence voltage is as follows: The sequence current is as follows: Combining Equations (8) and (9) yields the following: Equation (10) is represented by the positive-, negative-, and zero-sequence networks in Figure 4.  If an unbalanced wye-or delta-connected impedance load is connected to other phases (e.g., B-N, C-N, B-C, or C-A), the equivalent model in Figures 2 and 4 can be replaced by the corresponding impedance that is connected to other phases. The conventional short-circuit studies have often  If an unbalanced wye-or delta-connected impedance load is connected to other phases (e.g., B-N, C-N, B-C, or C-A), the equivalent model in Figures 2 and 4 can be replaced by the corresponding impedance that is connected to other phases. The conventional short-circuit studies have often ignored the load currents, because the magnitude of the SCC is higher than the magnitude of the load current. However, the proposed method transforms the unbalanced impedance loads to the equivalent SCC models and adds them to each sequence network, so that it can calculate the SCC more accurately. That is, it can take the effect of unbalanced loads on the SCC into account. However, the proposed unbalanced load model is only valid for the impedance loads. To take other load types (e.g., constant current and power load types) into account, the following approximation method is also used.

Transformation of Constant Current and Power Loads to Equivalent Impedance
The current that is flows to the constant power and current loads at the nominal voltage can be close to the current that is flows to constant impedance loads. Furthermore, the rated capacity of the constant power and current loads is often specified by active and reactive power in kVA or MVA. Thus, the constant power and current loads can be approximated by the constant impedance load in short-circuit analyses [21,22,27].
where V = the pre-fault voltage in per unit (p.u.), and S = the complex power scheduled at the load in p.u. However, since the voltage of the unbalanced loads during a fault is usually different from the pre-fault voltage, the proposed approximation has inherent errors to the actual SCC. Thus, this study compares the SCC that is approximated by the proposed method to the SCC that is calculated by the following superposition rule.

Superposition
The voltage of the faulted bus (V f ) in the SCC equations, presented in Table A1 of the Appendix A, is often assumed as 1∠0 • p.u. If ignoring the load in Figure 5, the SCC (I f ,no−load ) is as follows: The load current before the fault is as follows: As a result of the assumption of V f (e.g., 1∠0 • p.u.), the bus voltage at the slack generator should be increased by the following: Thus, the actual SCC that takes the load current into account can be determined by the following superposition rule: Note that the superposition rule based on Equation (15) can be applied for line-to-line (LL), line-to-line ground (LLG), and three-phase faults.

Pre-Fault Voltage Calculation
The proposed equivalent SCC model required the pre-fault voltage in order to more accurately calculate the SCC, so a power-flow calculation algorithm was presented in [31,32]. The method represented the power system network as the bus impedance matrix, or Zbus, and then calculated the power flow without any further processing of the impedance matrix. Indeed, Zbus represented the relationship between the injected currents and the induced voltages by the following:  bus V Z I . (16) To build the bus impedance matrix, in Figure 6, this study introduced the following four rules: (1) The slack node needed to be connected to the reference (e.g., the ground). If a slack node was connected to the reference, the initial Zbus is built by the following: (2) The new nodes needed to be connected to the slack node. If the new node j, with an impedance of z p.u., was connected to the previous node i, the following rule is used: For example, if the new nodes 2 and 3, in Figure 6, were connected to the slack bus (e.g., node 1), Zbus was modified by the following: (3) The new node needed to be connected to the existing two nodes. If the new node, with an impedance of z p.u., was connected to the existing two nodes (e.g., nodes i and j), the following rule modified the existing Zbus: For example, if a new node, with an impedance of 0.4j p.u., was connected to the existing nodes (e.g., nodes 1 and 3, thus i = 1 and j = 3), Zbus was modified by the following:

Pre-Fault Voltage Calculation
The proposed equivalent SCC model required the pre-fault voltage in order to more accurately calculate the SCC, so a power-flow calculation algorithm was presented in [31,32]. The method represented the power system network as the bus impedance matrix, or Z bus , and then calculated the power flow without any further processing of the impedance matrix. Indeed, Z bus represented the relationship between the injected currents and the induced voltages by the following: To build the bus impedance matrix, in Figure 6, this study introduced the following four rules: (1) The slack node needed to be connected to the reference (e.g., the ground). If a slack node was connected to the reference, the initial Z bus is built by the following: (2) The new nodes needed to be connected to the slack node. If the new node j, with an impedance of z p.u., was connected to the previous node i, the following rule is used: where Z and Z (v−1) ii = the element of the ith row and ith column, at the previous step. For example, if the new nodes 2 and 3, in Figure 6, were connected to the slack bus (e.g., node 1), Z bus was modified by the following: (3) The new node needed to be connected to the existing two nodes. If the new node, with an impedance of z p.u., was connected to the existing two nodes (e.g., nodes i and j), the following rule modified the existing Z bus : For example, if a new node, with an impedance of 0.4j p.u., was connected to the existing nodes (e.g., nodes 1 and 3, thus i = 1 and j = 3), Z bus was modified by the following: (4) The Kron reduction of the existing node. If the existing node was connected to the reference through an impedance of z p.u., the following Kron reduction rule was used: where n = the node connected to the reference, i and j = 1, . . . , n.
If the existing node (e.g., n = 3) was connected to the reference, through an impedance of 0.1j p.u., final Z bus is as follows: These four rules built the bus impedance matrices that calculated the power flow and SCC of the positive-, negative-, and zero-sequence networks of the test feeders.
(4) The Kron reduction of the existing node. If the existing node was connected to the reference through an impedance of z p.u., the following Kron reduction rule was used: where n = the node connected to the reference, i and j = 1, ... , n.
If the existing node (e.g., n = 3) was connected to the reference, through an impedance of 0.1j p.u., final Zbus is as follows:  The method presented in [31,32] calculated the currents that were to be injected and the voltages that were induced by the currents. If the actual pre-fault voltages were known at the first iteration, the method could determine the actual currents. But, the pre-fault voltages that were estimated at the first step were based on the initial nominal voltages. Thus, the steps that determined the injected currents and the induced voltages, or Equation (16), were iterated until they achieved the convergence in the following mismatch: However, the method presented in [31,32] did not model the tap-changing transformers and P-V buses. Therefore, this study decomposed the system into the series and shunt components, in order to model the tap-changing transformers. The detailed modeling methods of P-V buses and tapchanging transformers were done beyond this study, so they will be presented in Part II of this paper.

Seven-Bus Radial System
To verify the proposed iterative bus impedance method that was able to calculate not only the power flow in the normal steady state, but also the SCC that was affected by loads, a distribution system with unbalanced loads and a distributed generation (DG) source is modeled in Figure 7. The system included a substation transformer that was connected in a delta-grounded wye at 115/12.47 line-to-line kV and had a tap-changing position of eight, lines, a distribution transformer that was  The method presented in [31,32] calculated the currents that were to be injected and the voltages that were induced by the currents. If the actual pre-fault voltages were known at the first iteration, the method could determine the actual currents. But, the pre-fault voltages that were estimated at the first step were based on the initial nominal voltages. Thus, the steps that determined the injected currents and the induced voltages, or Equation (16), were iterated until they achieved the convergence in the following mismatch: ∆V However, the method presented in [31,32] did not model the tap-changing transformers and P-V buses. Therefore, this study decomposed the system into the series and shunt components, in order to model the tap-changing transformers. The detailed modeling methods of P-V buses and tap-changing transformers were done beyond this study, so they will be presented in Part II of this paper.

Seven-Bus Radial System
To verify the proposed iterative bus impedance method that was able to calculate not only the power flow in the normal steady state, but also the SCC that was affected by loads, a distribution system with unbalanced loads and a distributed generation (DG) source is modeled in Figure 7. The system Energies 2018, 11, 1447 9 of 21 included a substation transformer that was connected in a delta-grounded wye at 115/12.47 line-to-line kV and had a tap-changing position of eight, lines, a distribution transformer that was connected in a grounded wye-delta with a capacity of 10 MVA at 12.47/4.16 kV, and a DG source with a capacity of 1 MVA that injected active power to the grid as a voltage source. It also included three unbalanced loads, with a capacity of (a) 10 MVA at a unity PF in phase a, (b) 10 MVA at a lagging PF of 0.9 in phase b, and (c) 5 MVA at a lagging PF of 0.8 in phases b-c. Note that this case study initially used the constant power loads at a base of 100 MVA and it was a heavily unbalanced system. connected in a grounded wye-delta with a capacity of 10 MVA at 12.47/4.16 kV, and a DG source with a capacity of 1 MVA that injected active power to the grid as a voltage source. It also included three unbalanced loads, with a capacity of (a) 10 MVA at a unity PF in phase a, (b) 10 MVA at a lagging PF of 0.9 in phase b, and (c) 5 MVA at a lagging PF of 0.8 in phases b-c. Note that this case study initially used the constant power loads at a base of 100 MVA and it was a heavily unbalanced system. To model the tap-changing transformers, this study modified the power-flow calculation method that was presented in [31,32], by decomposing the system in Figure 8 into the series and shunt components. That is, the tap changer that was installed on the primary side of the transformer was converted to the equivalent π model and was set at a tap of eight (so, a = 1, b = 1 − 0.00625 × 8 = 0.95, and Yex,pt = 0). Using the presented four rules of building the impedance matrix, this study built the following bus impedance matrix, with a size of 7 by 7, which included only the series impedances between the buses in Figure 8, 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0021 + 0.1429i -0.0025 + 0.1276i -0. bus Z  0072 + 0.1122i -0.0072 + 0.1122i -0.0072 + 0.1122i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.0025 + 0.1276i 0.0181 + 0.1860i 0.0087 + 0.1645i 0.0087 + 0.1645i 0.0087 + 0.1645i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.0072 + 0.1122i 0.0087 + 0.1645i 0.0246 + 0.2169i 0.0246 + 0.2169i 0.0246 + 0.2169i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.0072 + 0.1122i 0.0087 + 0.1645i 0.0246 + 0.2169i 0.0546 + 0.2969i 0.0546 + 0.2969i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.0072 + 0.1122i 0.0087 + 0.1645i 0.0246 + 0.2169i 0.0546 + 0.2969i 0.0846 + 0.3769i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i p.u.

. Short-Circuit Current Calculation
To determine the SCC of the proposed case study, this study applied the pre-fault voltage of a bus (to which the loads were connected) to Equations (6), (10), and (11). In other words, the first two loads that were connected in a wye were transformed to the equivalent impedance by the following:  .3615 The third load that was connected in a delta was transformed to the equivalent impedance by the following: 22   The method iteratively calculated the injected currents and induced voltages, using the matrices (e.g., Equations (25) and (26)) until the convergence was achieved. This study compared the power-flow calculation results to the verified methods (e.g., the Newton-Raphson, Gauss-Seidel, and fast decoupled methods), in Table 1. The proposed method showed the same exact pre-fault voltages as those that were determined by the other methods, and successfully modeled the tap-changing transformer with an off-nominal tap ratio, which could be seen in the voltage of bus 2 (e.g., 1.03446∠−1.528 • p.u.). Thus, the pre-fault voltage could be used in the proposed SCC model, which transformed the load to the equivalent impedance in Equation (11).

Heavily-Meshed Network
To verify the proposed methods in a heavily-meshed network, the IEEE 30-bus test system was modeled, which was a power system network that was presented from American Electric Power in December 1961 [33,34], in Figure 11. It included 30 buses, a slack generator (represented as a slack bus on bus 1), five generators (represented as a P-V bus on buses 2, 5, 8, 11, and 13), two shunt capacitors, and 21 loads. To verify the proposed short-circuit analysis method in such a heavilymeshed network, the following conditions were assumed.
(1) Since the zero-sequence impedance data of the system was not available, they were set to as much as 2.75 times the corresponding positive-sequence impedance data [35]. (2) Since the detailed connection data of the transformers was not available, they were connected in a grounded wye-grounded wye configuration and all of the loads were unbalanced in a delta configuration. (3) All of the initial and target voltages of the generators were set to 1.0∠0° p.u., but it could have been specified differently, if necessary. The other data were set to the same data as was presented in [33], including the tap ratios of the tap-changing transformers and reactive power limits of the P-V buses.

Heavily-Meshed Network
To verify the proposed methods in a heavily-meshed network, the IEEE 30-bus test system was modeled, which was a power system network that was presented from American Electric Power in December 1961 [33,34], in Figure 11. It included 30 buses, a slack generator (represented as a slack bus on bus 1), five generators (represented as a P-V bus on buses 2, 5, 8, 11, and 13), two shunt capacitors, and 21 loads. To verify the proposed short-circuit analysis method in such a heavily-meshed network, the following conditions were assumed.
(1) Since the zero-sequence impedance data of the system was not available, they were set to as much as 2.75 times the corresponding positive-sequence impedance data [35]. (2) Since the detailed connection data of the transformers was not available, they were connected in a grounded wye-grounded wye configuration and all of the loads were unbalanced in a delta configuration. (3) All of the initial and target voltages of the generators were set to 1.0∠0 • p.u., but it could have been specified differently, if necessary. The other data were set to the same data as was presented in [33], including the tap ratios of the tap-changing transformers and reactive power limits of the P-V buses.
The proposed power-flow method showed the same results as those that were calculated by the Gauss-Seidel, Newton-Raphson, and fast decoupled methods, which are presented in the Appendix. Then, this study generated an SLG fault on bus 5. In order to take the effect of the unbalanced loads on the SCC into account, this study calculated the pre-fault voltage by the proposed power-flow algorithm, transformed the loads to the equivalent impedances, and added them to the sequence networks. For example, the load with a capacity of 21.7 + 12.7j MVA that was connected to bus 2 was transformed by the following:  Figure 11. The heavily-meshed IEEE 30-bus test system [33,34].
The proposed power-flow method showed the same results as those that were calculated by the Gauss-Seidel, Newton-Raphson, and fast decoupled methods, which are presented in the Appendix A. Then, this study generated an SLG fault on bus 5. In order to take the effect of the unbalanced loads on the SCC into account, this study calculated the pre-fault voltage by the proposed power-flow algorithm, transformed the loads to the equivalent impedances, and added them to the sequence networks. For example, the load with a capacity of 21.7 + 12.7j MVA that was connected to bus 2 was transformed by the following: where V B2 = 0.977819∠ − 6.035 • p.u. in Table A3 in the Appendix A. Then, this study calculated the SCC by the following: where , and Z f = 0. Table 5 shows the SCC if an SLG fault occurred on bus 5. The actual SCC was calculated by the superposition rule, based on Equation (15). The unbalanced loads increased the magnitude of the SCC, compared with the method that ignored loads (e.g., the conventional method). It was because the unbalanced load conditions could decrease the impedance that was seen from the faulted bus by adding the SCC models of the unbalanced loads. Thus, as either the distributed or clustered loads that were unbalanced in phases were connected to the grid, the phase angle and magnitude of the SCC should have been accurately calculated. In other words, the loads should have been taken into account. As the next validation step, this study presented the magnitude of the SCC when changing the faulted bus from buses 2 to 30, in Figure 12. For example, if an SLG fault occurred on phase a of bus 2, close to the slack generator (e.g., bus 1), the highest SCC of 12.4659∠−69.71 • p.u. would flow. On the other hand, if the fault occurred on the bus (e.g., bus 26) that was remote from the slack generator, the lowest SCC of 0.7160∠−62.62 • p.u. would flow, which was comparable to [27,28,36]. That is, as an SLG fault occurred closer to the slack generator, the higher SCC would flow. Figure 13a,b show a surface map and a contour map of all of the bus voltages when the fault location was changed from buses 2 to 30. In Figure 13a,b, if an SLG fault occurred on phase a of bus 2, the post-fault voltage magnitudes of the slack bus (e.g., bus 1) and the faulted bus (e.g., bus 2) indicated 1.0 p.u. and 0 p.u., respectively. The voltage magnitudes of the buses behind the faulted bus also showed about 0 p.u. If the fault occurred on bus 30, the voltage magnitudes of the slack bus (e.g., bus 1) and the faulted bus (e.g., bus 30) showed 1.0 p.u. and 0 p.u., respectively. In Figure 13b, the phase voltage of the faulted bus was 0 p.u. after an SLG fault. In the contour map, the phase voltage magnitude of 0 p.u. of each faulted bus changed diagonally, as the faulted bus was changed from 2 to 30.

Conclusions
The main objective of this study is to present a short-circuit analysis algorithm that is able to analyze unbalanced loads. For this purpose, this study derives an SCC model of the unbalanced wyeor delta-connected impedance load. Since the proposed SCC model requires the pre-fault voltage, it uses a power-flow analysis algorithm that iteratively calculates the current that needs to be injected and the pre-fault voltage, using the bus impedance matrix. The case studies of the unbalanced test feeders indicate that the proposed SCC model, which transforms the unbalanced loads to the equivalent short-circuit impedances, shows the least AEs to the actual SCC, compared with the conventional methods. In other words, the proposed short-circuit analysis method could more accurately calculate the SCC, since it did not neglect the unbalanced loads. Thus, the proposed shortcircuit method is able to take the unbalanced wye-or delta-connected loads into account. The proposed SCC model is based on an unbalanced load, so it is applicable to the unbalanced systems. The proposed SCC models are also applicable for the cases either with or without DG. Since the proposed SCC models are more accurate than the conventional methods that neglect the load, the proposed SCC models of delta-or wye-connected unbalanced loads in phases can be incorporated in the power systems analysis programs (e.g., DIgSILENT). However, this study did not examine the effect of unbalanced loads with various capacities that have uncertainty in their capacity and location on the SCC. In addition to the uncertainty, it did not present a method that transforms constant current and power loads into equivalent short-circuit impedance. Implementing these topics in future should provide a more accurate method for short-circuit studies.
Funding: This work was supported by the Inha University Research Grant #57843-1.

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

AE
absolute error BIBC bus-current-injection to branch-current BCBV branch-current to bus-voltage BFSW backward and forward Sweep DG distributed generation

Conclusions
The main objective of this study is to present a short-circuit analysis algorithm that is able to analyze unbalanced loads. For this purpose, this study derives an SCC model of the unbalanced wyeor delta-connected impedance load. Since the proposed SCC model requires the pre-fault voltage, it uses a power-flow analysis algorithm that iteratively calculates the current that needs to be injected and the pre-fault voltage, using the bus impedance matrix. The case studies of the unbalanced test feeders indicate that the proposed SCC model, which transforms the unbalanced loads to the equivalent short-circuit impedances, shows the least AEs to the actual SCC, compared with the conventional methods. In other words, the proposed short-circuit analysis method could more accurately calculate the SCC, since it did not neglect the unbalanced loads. Thus, the proposed short-circuit method is able to take the unbalanced wye-or delta-connected loads into account. The proposed SCC model is based on an unbalanced load, so it is applicable to the unbalanced systems. The proposed SCC models are also applicable for the cases either with or without DG. Since the proposed SCC models are more accurate than the conventional methods that neglect the load, the proposed SCC models of delta-or wye-connected unbalanced loads in phases can be incorporated in the power systems analysis programs (e.g., DIgSILENT). However, this study did not examine the effect of unbalanced loads with various capacities that have uncertainty in their capacity and location on the SCC. In addition to the uncertainty, it did not present a method that transforms constant current and power loads into equivalent short-circuit impedance. Implementing these topics in future should provide a more accurate method for short-circuit studies.
Funding: This work was supported by the Inha University Research Grant #57843-1.

Conflicts of Interest:
The authors declare no conflict of interest.  Table A1 shows the SCC of each sequence network for SLG, LLG, LL, and 3-phase faults. Table A1. Fault current of each sequence network [37].

Type
Positive-Sequence (I + i ) Negative-Sequence (I − i ) Zero-Sequence (I 0 i ) where V f = the Thevenin voltage source, typically 1.0∠0 • p.u., Z + i , Z − i , and Z 0 i = the ith diagonal element of each bus impedance matrix, Z f = the short-circuit impedance. Table A2 shows the SCCs for the SLG, LLG, LL, and 3-phase faults of the seven-bus radial system. Expectedly, the magnitude of the 3-phase SCC is the highest. The proposed method indicates the least AEs to the actual SCC, is determined by the superposition rule based on (15).  The power-flow solutions of the IEEE 30-bus test system, calculated by the proposed method, show the same results as those that are calculated by the Gauss-Seidel, Newton-Raphson, and fast decoupled methods. Note that the reactive power limits of the P-V buses are set to the same data as presented in [33].