UPFC Location and Performance Analysis in Deregulated Power Systems

We deal with the effect of Unified Power Flow Controller UPFC installation on the objective function of an electricity market. Also this paper proposes a Novel UPFC modelling in OPF which facilities the consideration of the impact of four factors on power market. These include the series transformer impedance addition, the shunt reactive power injection, the in-phase component of the series voltage and the quadrature component of the series voltage. The impact of each factor on the electricity market objective function is measured and then compared with the results from a sensitivity approach. The proposed sensitivity approach is fast so it does not need to repeat OPF solutions. The total impacts of the factors are used to offer UPFC insertion candidate points. It is shown that there is a clear match between the candidate points of the sensitivity method and those proposed by the introduced UPFC modelling in our test case. Furthermore, based on the proposed method, the relation between settings of UPFC series part and active and reactive power spot prices is presented.


Introduction
Limitations in transmission and generation system expansion, such as right-of-way and environmental problems, have made it inevitable to use the current network capacity as much as possible 1 . The competition in a restructured power system leads to it's optimization and new ways for cost reduction. Flexible AC Transmission Systems FACTS , which are developed as a result of recent progress in power electronic technology and communication systems, have opened alternative ways of increasing loadability, better network control and cost reduction. FACTS devices can be used for congestion management 2 , energy loss minimization 3 , power flow control 4 , security enhancement 1 , social welfare maximization 5 and network stability improvement 6 .
To manage power pricing in a PoolCo power market, an ISO implements Optimal Power Flow OPF in which the main objective is to maximize social welfare subject to some network constraints 7, 8 . FACTS settings in steady state applications are determined together with optimal power flow variables in a single unified framework. In some electricity markets, ISO may own all FACTS devices. In this case, it is responsible for both their operation and planning. On the other hand, in some electricity markets, FACTS devices may be owned by different entities that are paid by ISO based on "Ancillary Services" they provide for ISO. In this case, also, ISO controls FACTS devices; but studies related to FACTS planning, which we deal with in this paper, is a subject of interest for FACTS investors.
Among FACTS devices, the Unified Power Flow Controller UPFC is able to simultaneously compensate reactive power and control active and reactive power flows of a transmission line 9 . By employing UPFCs, electricity generation cost and active power losses can be reduced 1, 10 . Using a UPFC in a power market in order to minimize the market cost may lead to the reduction in spot prices of load buses 5 . Both real and reactive power spot prices may subsequently change drastically. An impact on transmission cost allocation in a power market, as a result of UPFC operation, has been reported in 11 . In spite of the above mentioned steady state effects of UPFCs, to the best of our knowledge, no discussion has so far been presented about the desired UPFC settings from an OPF solution and the effect of each of the UPFC functions, including shunt reactive power compensation and active and reactive power flow control.
In this paper, a proposed detailed UPFC modeling, including internal UPFC state and control variables and serial and shunt impedances, is incorporated in the OPF formulation. Through that, the factors influencing the objective function of an electricity market resulting from UPFC installation, namely, the series transformer impedance insertion, the shunt reactive injection, the in-phase component of the series voltage and the quadrature component of the series voltage have taken into account. Also, to measure the impact of these components on improving the objective function of an electricity market, two approaches, namely a differencing method and a sensitivity analysis, are presented. The above four impacts are added in each case to identify the potential points for UPFC installation. The sensitivity approach is fast as it needs to run OPF only once in the base case system without UPFC, to derive the sensitivity coefficients. Therefore, the computational burden in more accurate UPFC allocation techniques such as 12-14 could be significantly decreased if this approach is used to limit the search space. The relation between UPFC series part settings and Locational Marginal Prices LMP is another subject presented in this paper. This paper is organized as follows: In Section 2, optimal power flow and its implementation are presented. Incorporation of the UPFC modelling in OPF is described in Section 3. Then in Section 4, to validate the proposed approaches, a UPFC is placed on all possible points of a test system and the impacts of all pre-mentioned components on improving the objective function of an electricity market are computed by the two approaches. UPFC allocation is also discussed in this section. Finally, concluding remarks are presented in Section 5.

Optimal Power Flow Implementation
The main objective of an electricity market is to maximize the social welfare which consists of bid prices of generation units and loads 8 . For the sake of simplicity, customers' loads are assumed to be constant. However, consideration of more accurate load models and bid prices of customers are also possible. The mathematical formulation of an optimal power  flow problem can be expressed as where the cost function, f, is the total bid offers of the generators. Note that since we have assumed that the price elasticity of demand is zero, minimizing f is equivalent to maximizing the welfare 8 ; when no UPFC is installed the control variables, w, are active power generations, P G , and reactive power generations, Q G . The state variables, t, include load bus voltages, V L , and load bus angles, θ L . The equality constraints, g, in the optimization are nonlinear AC load flow equations. The inequality constraints, h, are as following.

Gi
Upper and lower active powers of generator-i, Upper and lower reactive powers of generator-i, Upper and lower voltage magnitudes of generator-i, Upper and lower voltage magnitudes of load bus-i.
In this paper the optimal power flow solution is based on separating the control variables, w, from the state variables, t 15 . The algorithm of optimal power flow is shown in Figure 1.

Mathematical Problems in Engineering
Former side End side Figure 2: UPFC equivalent circuit.

Novel UPFC Modelling in OPF
The Acha's UPFC modelling 10, 16 consists of two voltage sources and two impedances representing the series and the shunt converters and transformers in a UPFC, as shown in Figure 2 where i Z s and Z p represent series and shunt transformers leakage impedances, respectively ii I S and I p denote series and shunt converter currents, respectively iii Q P is the net shunt reactive power injected to the former side bus The UPFC control parameters in 10, 16 modelling are the amplitude and the angle of the series converter voltage phasor U s , ϕ s and the amplitude and the angle of the shunt converter voltage phasor U p , ϕ p . However, none of these control parameters are directly effective in active or reactive power flow from UPFC converters. Thus, it makes this model inappropriate to use for a performance analysis. In this paper, that modelling is enhanced to resolve this issue. The UPFC control parameters of proposed model include the in-phase and quadrature components of the series converter voltage U sx , U sy as shown in Figure 3 a , and the in-phase and quadrature components of the shunt converter voltage U px , U py as shown in Figure 3 b . U sx and U px are at the same angle as V 1 while U sy and U py are perpendicular to V 1 . These parameters can mathematically be expressed as U s U sx jU sy × e jδ 1 , U p U px jU py × e jδ 1 .

3.1
Under normal operating conditions of a power system, δ 1 − δ 2 and V 1 − V 2 are small and the resistances of Z s and Z p are small, as well. Thus it can be supposed that, U sx and U sy influence only reactive power flow and active power flow from bus 1 to bus 2, respectively. In other words, the in-phase and the quadrature components of the UPFC series voltage are comparable in operation to a tap changer and a phase-shifter, respectively Figure 3 a . Figure 3: Phasor diagrams of a series converter and b shunt converter voltages.

Mathematical Problems in Engineering
On the other hands U px and U py are responsible for flowing the reactive and active powers, respectively, in the shunt part of the UPFC equivalent circuit in Figure 2.
In order to incorporate the proposed UPFC model into the OPF algorithm, three UPFC parameters, namely, U sx , U sy and U px , should be added to the set of the optimization control variables, w, and at the same time, the only remaining parameter, U py , should be added to the set of the state variables, t. According to U py function, this parameter is incorporated into the Jacobian matrix and mismatch equations of the load flow to satisfy the active power balance equation in the UPFC. Also, the UPFC operational limits given below should be added to the optimization inequalities, h.

UPFC Performance Analysis in a Power Market
The UPFC model is composed of two voltage sources and two impedances, representing physical converters and transformers. To determine how much the installation of a UPFC may affect a power system, we can include the components of the UPFC model, one by one, and accordingly, study the effect of each. First, we import the series and shunt impedances. Note that the resistances of the transformers can be neglected as they are much smaller than their reactances. Then, the series voltage components are enabled. U sx and U sy are independent variables and enabling them has an impact on the system. On the other hand, among the shunt voltage components, U px has a similar behaviour and could be treated similarly. However, U py is a dependent variable and is modified according to the other three control variables, to keep the active power balance in the UPFC. So, once for instance, U sx is enabled, U py would change accordingly and therefore, its effect would be taken into account and as such needs not to be calculated separately. Besides, some part of the shunt reactive power produced due to enabling U px is lost in the shunt impedance, Z p , which can be considered as the main influence of Z p . Therefore, the effects of U px and the shunt impedance can be combined if, instead, we consider the net reactive power injected to the former side bus, Q P .
In brief, the influence of UPFC installation on a power market can be considered as the total impacts of four functions. i The insertion of the UPFC series transformer impedance on the line ii Reactive power injection, Q P , at the former side bus due to U px iii Reactive power flow in the series part due to U sx iv Active power flow in the series part due to U sy The series transformer impedance addition regularly increases the OPF objective function since it increases the line impedance. The next three components are the variables of the optimization and should decrease the objective function.

The Relation between UPFC Series Voltage Components and LMPs
In an electricity power market, when the power price active LMP at the sending end of a transmission line is cheaper than the one at the receiving end, flowing the active power from the sending end to the receiving end is desirable. In this case, given a UPFC installed at the sending end of the line, U sy should be set at a positive value to produce this flow and vice versa. Likewise, this rule is also true about the reactive power. That is to say, if the reactive LMP at the sending end of the transmission line is less than the one at the receiving end, U sx ought to be set at a positive value to cause this flow and vice versa. If the maximum thermal current of the line is reached in the base case, UPFC series parameters are, however, set in a different manner. U sx and U sy settings should be so selected in this case to decrease the line current.

Case Studies and Results Analysis
Validation tests are performed on the six bus 11 lines test system shown in Figure 4  The OPF results of the test system are summarized in Table 1. The OPF cost in this electricity market, f, is 8815.09 $/hr. The reactive power generation of G3, Q G3 , and the current at the receiving ends of line 2-4 and line 1-5, I l4-2 and I l5-1 , are set to their maximum values. In order to allocate UPFCs in a system, all possible locations should be evaluated and the number, location and size of the UPFCs should be determined. The possibility of installing the UPFCs at both ends of all lines in the six bus system is considered which constitutes 22 cases. Then, the optimal power flow problem is solved for these installation cases and finally, the OPF costs are compared.
With reference to the size, the maximum UPFC series voltage can be up to 0.5 pu of the line voltage 17 or even more. This determines the series converter MVA. Also, the series and shunt converters may 17 or may not 12 have the same size. Determination of converter MVAs is a matter of UPFC allocation. In this paper, however, we would like to find some areas as candidate points for UPFC installation by performing a sensitivity analysis. Once the promising candidate points are determined by this approach, more precise UPFC allocation algorithms such as 12, 13 would be necessary to select the final points. So, with regard to the purpose of this paper, the UPFC series and shunt converters are sized into the relatively small fixed value of 4 MVA, equal to 0.04 pu with S base 100 MVA. Since the current ratings of the test system lines are 0.4 pu on average, the UPFC maximum series voltage would be typically 0.1 pu. Also, the assumption of constant sized converters removes the need to calculate the UPFC investment cost.
Apart from the size of converters, other UPFC ratings may vary at different points. The maximum voltage magnitude of the shunt converter, U max p , is always a bit more than the line nominal voltage. Here, it is chosen 1.2 pu in all cases. Since U p is normally about 1.0 pu, the maximum current of the shunt converter, I max p , will be the same as the converter apparent power rating, 0.04 pu, in all cases. The maximum current of the series part, I max s , is practically selected to be equal to the line current thermal limit 17 . Given the nominal power and the maximum current of the series converter, the maximum voltage magnitude of the series converter, U max s , is

4.1
The resistances and the reactances of the coupling transformers are chosen from typical figures based on their voltage level and nominal power.
A differencing method which includes the following steps, applied to all 22 UPFC placement cases.
i By letting three UPFC control parameters free, run optimal power flow and obtain UPFC settings U * px , U * sx and U * sy and OPF cost function, f 4 opf .
ii Use UPFC in zero compensation mode Q p 0, U sx 0 and U sy 0 and obtain OPF cost function, f 1 opf .  Figure 5: UPFC equivalent circuit and power injection model 5 .
iii Use UPFC in the operating mode Q p free, U sx 0 and U sy 0 and obtain OPF cost function, f 2 opf .
iv Use UPFC in the operating mode Q p free, U sx U * sx and U sy 0 and obtain OPF cost function, f 3 opf .
In each step, one of the four UPFC elements effective in changing the OPF cost function is added. Then, the objective function of the step, f k opf , is obtained by the OPF solution. So the OPF cost function alteration caused by adding an element y k , Δf k difr , can be computed as where the OPF cost for the base case system with no UPFC, f 0 opf , is 8815.09 $/hr as given in Table 1. The change in the OPF cost function due to enabling an element y k can also be calculated by a sensitivity analysis, Δf k sens , as shown in 4.3 .
where y * 1 is the series transformer leakage impedance; y * 2 denotes the net reactive power injected by the shunt converter, Q p ; y * 3 and y * 4 are the in-phase and the quadrature components of the series voltage, respectively, obtained in step I; ∂f/∂y k is the OPF cost function sensitivity with respect to the element y k . The sensitivity factors are calculated using the OPF results of the main system with no UPFC.
The sensitivity of OPF objective function with respect to shunt converter reactive power injection, ∂f/∂Q p , is equal to the reactive LMP at the bus to which the UPFC is connected. Active and reactive LMPs are the Lagrangian multipliers of power flow equations in optimal power flow, which are obtained after solving OPF. ∂f/∂U sx and ∂f/∂U sy coefficients can be calculated using Figure 5. Suppose that a UPFC is installed at the sending end of line 1-2. Figure 5 a shows the equivalent circuit of the UPFC 5 in which I t and I q are the in-phase and the quadrature components of the shunt converter current with respect to V 1 .

Mathematical Problems in Engineering 9
Injecting powers P 1S , Q 1S , P 2S and Q 2S in Figure 5 b is equal to UPFC insertion on the sending end of line 1-2 in Figure 5 a . These powers can be represented in terms of U sx , U sy and I q as The consequence of these power injections in changing OPF cost function can be estimated by LMPs. In order to compute ∂f/∂U sx , the variables U sy and I q in 4.4 -4.7 are set to zero and the chain rule is used as where ALMP 1 , RLMP 1 , ALMP 2 and RLMP 2 are the active and the reactive LMPs at line 1-2 both ends. Also, ∂I l1-2 /∂U sx and ∂I l2-1 /∂U sx are the derivatives of the current through line 1-2 with respect to U sx ; likewise, λ I l1-2 and λ I l2-1 are the Lagrangian multipliers of the maximum current constraints at the sending and the receiving ends of line 1-2. However, the last two terms in 4.8 may seem to be irrelevant. The reason these terms are added can be explained as follows: UPFC series voltage causes the line current to change. This change, when the maximum line current is binding, produces a second change in OPF cost which can be estimated using the maximum current Lagrangian multiplier. In our test case, nonetheless, the maximum current multiplier only at the receiving end of line 2-4 and at the receiving end of line 1-5 is nonzero. ∂I l1-2 /∂U sx in 4.8 can be simply derived based on the definition of I l1-2 .
where S l1-2 , P 12old and Q 12old are the apparent, active and reactive powers of line 1-2 while U sx , U sy and I q are set to zero. A similar procedure can be employed to calculate ∂f/∂U sy i.e., let U sx and I q be zero and use an equation similar to 4.8 .  ∂f/∂X s can be calculated by substituting power injections P 1S , Q 1S , P 2S and Q 2S for the series transformer impedance, X s , as shown in Figure 6. In a similar way to 4.8 , we obtain for X s

4.10
where ∂I l1-2 /∂X s is calculated by an equation similar to 4.9 . After calculating the injection powers in Figure 6, the derivatives in 4.10 are obtained as:

4.11
The values of ∂f/∂y k , y * k and Δf k sens for k 1, . . . , 4 are shown in Tables 2 through 5

Line Impedance Increase
From Table 2, it is evident that inserting the UPFC series transformer at either the sending end or the receiving end of a line produces roughly similar change in the OPF cost function. Furthermore, in most cases 13 cases out of 22 , the OPF cost function increases when the UPFC series transformer impedance is inserted.

Shunt Reactive Power Injection
By reviewing Δf 2 difr in Table 3 and comparing the results of UPFC insertion on all lines connected to a particular bus, it may be concluded that connecting the UPFC to a certain bus, no matter on which line, would approximately lead to the same amount of shunt reactive compensation. For example, in installing the UPFC at the receiving end of line 2-3 and the sending ends of lines 3-5 and 3-6 in which the UPFC is connected to bus 3, Q p takes very close values of 2.49, 3.85 and 3.63 $/hr, respectively. Consequently, the results of Table 3 are grouped according to the buses not the lines. Also, it can be seen that whenever a UPFC is connected to one of the load buses, the OPF sets the shunt converter current, I p , to its maximum value, that is 0.04 pu. These cases are marked by * in Table 3. This is due to the fact that by producing reactive power through a UPFC, active power loss as a result of reactive power flow on transmission lines would decrease.

Enabling U sx and U sy
The U sx and U sy compensation results are presented in Tables 4 and 5, respectively. The first and the second row of each line in both tables represent the results of placing UPFC at the sending and the receiving ends of the line, respectively. According to Tables 4 and 5, Δf 3 difr and Δf 4 difr are constantly negative; so, it may be concluded that enabling series voltage components would always cause the objective function of OPF to decrease. Also, it should be noted that the results of the sensitivity analysis, Δf 3 sens and Δf 4 sens , are usually greater than the differencing results, Δf 3 difr and Δf 4 difr ; the exception cases are shown in bold highlighting. Hence, it seems that by moving away from the initial operating point, the compensation slopes of the in-phase and the quadrature components decrease. An important thing to note is that U * sx at the sending end of a line is often very close to −U * sx at its receiving end. For example, the U * sx values in Table 4 for the sending and the receiving ends of line 2-3 are 0.056 and −0.055, respectively. This is also true about U * sy . Thus, moving a UPFC from one end of a line to the other end appears to have low effect on the U * sx or U * sy absolute value. Figure 7 can be used to examine the proposed approach, explaining the relationship between U * sx and U * sy settings and LMPs in an electricity market. It shows both active and reactive LMPs of each system bus inside a box beside the bus . These LMPs are derived from the OPF on the base case system without UPFC. The illustrated arrows at both ends of each line show the directions of the active and reactive power flows as a result of U sy and U sx activation, respectively. Also, the magnitudes of the settings U * sx and U * sy presented in Tables 4 and 5 are shown above each arrow.
The first part of the proposed approach is now applicable to all 22 cases except the four cases of UPFC insertion on line 1-5 and line 2-4, in which the current is set to the maximum value. It is shown that the approach truly predicts all the cases for the U * sx and U * sy settings. Lines 1-5 and 2-4, drawn by bold lines in Figure 7, are operating at their current thermal limit. Thus, the second part of the proposed approach should be evaluated in these cases. Active and reactive powers flow from bus 2 to bus 4 and the chosen U * sx and U * sy values at both ends of this line cause the line current to reduce, verifying the proposed approach. This is also the case in line 1-5 for U * sy ; however, U * sx values in line 1-5 do not follow the approach and are depicted by double line arrows in Figure 7. These violations are not surprising because the OPF problem shows a high degree of nonlinearity. Altogether, it seems that both parts of the approach efficiently predict the relationship between UPFC series voltage components and LMPs in an electricity market.

Determination of UPFC Installation Candidate Points Using Total Effects of Components
The impacts of the four elements on the OPF cost function in 22 cases are summarized in the stacked column chart shown in Figure 8. There are two columns for each of the 11 transmission lines in the figure. The left and the right columns are associated with UPFC installation at the sending and the receiving ends of the line, respectively. Each column consists of four stacked columns related to the four elements. The first stacked column represents the impact of the series transformer impedance insertion, represented by a vertical arrow. This element in some cases, such as UPFC insertion on both sides of line 2-4, has a positive effect and in some other cases, such as UPFC installation on both sides of line 1-4, has a negative effect on the cost saving. Other elements, however, have always positive effects.
The total compensation of UPFCs can be identified by comparing the total column heights. It can be seen that after enabling the four components, the OPF cost is reduced in all the 22 cases. Another important thing can be inferred from the values for the lines in which one end is a generation bus and the other end is a load bus, including lines 1-4, 1-5, 2-4, 2-6, 3-5 and 3-6. In these cases, it is observed that UPFC installation at the load bus end of the line is more beneficial at the generation bus end. The reason is that the reactive compensation is much more at the load bus end while the other components produce almost the same results at either end.
Six cases out of 22 in which UPFCs have produced the most improvement are marked by * in Figure 8. These six cases are associated with UPFC installation on both ends of lines 1-2, 1-4 and 2-4. Since simultaneous insertion of a UPFC at both ends of a line is unrealistic, candidate points to install UPFCs in the six bus system appear to be the receiving ends of lines 1-2, 1-4 and 2-4. Figure 9 shows the results of the total UPFC cost reductions by both approaches, normalized based on their respective maximum values. It can be seen that both approaches  show the same pattern of compensation at different points of the system. Thus, it confirms the trustworthiness of the sensitivity approach. Furthermore, six points with the highest figures in the differencing and the sensitivity approaches are distinguished by * and marks, respectively, in Figure 9. It is shown that the two approaches offer the same candidate points. Hence, the proposed sensitivity analysis seems to be, effectively, capable of determining the candidate points. From a computational point of view, while the sensitivity method requires only one OPF run and some post studies, the differencing method needs much more calculations, that is, in our test case, 23 OPF runs, one for base case and 22 ones for UPFC installation on all lines. In order to assess how much saving can be obtained through UPFC installation, the cost of UPFC installation must be calculated. The cost of installation of UPFC is taken from Siemens database and reported in 18 given by 4.12 . where C UPFC is the cost of UPFC in US$/kVA and S is the operating range of UPFC in MVA. Therefore, based on the supposed size of UPFC in our case studies, the cost of UPFC installation will be about 749,000 $. This cost will have to be compared with the revenue or benefit that can be derived from UPFC. The revenue derived from UPFC, shown in Figure 8, has the unit of "$/hr" depending on the utilization and level of congestion. In order to compare the cost of FACTS against the anticipated benefits, they have to be converted to a common unit. In this paper, the comparison is made by converting the cost, as well as the benefit or revenue into annuity "$/year" . To compute the annual capital cost and benefit revenue of FACTS, following assumptions have been made: Project lifetime n : 5 years Discount rate r : 10% Average utilization u : 40% Operational cost of FACTS device is neglected.
Annual capital cost of FACTS in $/year can be found as 19 : Thus, the annual capital cost of UPFC in our test case is 197,000$/year. Annual revenue from use of UPFC in $/year can be determined as 19 : The average utilization u gives the percentage of time the UPFC device is considered 100% effective. Since, the demand and supply patterns change during different time period, leading to different price quantity relationship and consequently different setting for FACTS devices. At low load period, the effectiveness of UPFC devices decreases and hence the revenue benefit from use of UPFC decreases. So, to evaluate the benefit of UPFC a utilization factor is considered. Considering the best case, UPFC installation at the receiving end of line 1-4, the annual revenue generated due to UPFC is 217,000 $/year. Consequently, about US$20,000 can be saved each year.

Conclusions
In this paper, a new explicit model for UPFC was proposed in which the parameters were assigned to the active and reactive power flows in the series and shunt parts of the UPFC. Using the proposed model, UPFC settings and power prices in a restructured power market were simultaneously determined to maximize the social welfare. Also based on the proposed model, impact of UPFC installation on the social welfare was considered to be the result of four elements.
By studying the test system with different UPFC positions, the effect of each element on the power market objective function was observed by means of a differencing method. Then, the total UPFC compensations in different cases were compared and suitable UPFC insertion points were suggested. The comparative results obtained by a sensitivity approach showed that two approaches offer almost the same candidate points in our case. Since the results of the sensitivity approach are calculated without repeating OPF solutions, the method is faster than the differencing method. Eventually, based on the functions of UPFC series voltage components, two rules for predicting the sign of these components in an electricity market were proposed and their effectiveness was practically confirmed by case studies. The OPF cost function sensitivity with respect to an element y k V i pu

Mathematical Symbols
The voltage magnitude at but-i Δδ rad The difference between the voltage angles at buses 1 and 2 r 12 pu The resistance of line 1-2 x 12 pu The reactance of line 1-2 g 12 pu The conductance of r 12 jx 12 b 12 pu The susceptance of r 12 jx 12 B pu The shunt susceptance of line 1-2 ALMP 1 $/hr·MW The active LMP at bus-1 RLMP 1 $/hr·MVAR The reactive LMP at bus-1 ALMP 2 $/hr·MW The active LMP at bus-2 RLMP 2 $/hr·MVAR The reactive LMP at bus-2 λ I l1-2 $/hr·pu The Lagrangian multiplier of the maximum current constraint at the sending end of line1-2 λ I l2-1 $/hr·pu The Lagrangian multiplier of the maximum current constraint at the receiving end of line1-2 I t pu The in-phase component of the shunt converter current with respect to V 1 I q pu The quadrature component of the shunt converter current with respect to V 1 S l1-2 MVA The apparent power through line 1-2 at its sending end when the UPFC is disabled P 12old MW The active power flow at the sending end of line 1-2 when the UPFC is disabled Q 12old MVAR The reactive power flow at the sending end of line 1-2 when the UPFC is disabled P 21old MW The active power flow at the receiving end of line 1-2 when the UPFC is disabled Q 21old MVAR The reactive power flow at the receiving end of line 1-2 when the UPFC is disabled