Next Article in Journal
Numerical Investigation of the Three-Dimensional HCIR Partial Differential Equation Utilizing a New Localized RBF-FD Method
Next Article in Special Issue
Fractal Complexity of a New Biparametric Family of Fourth Optimal Order Based on the Ermakov–Kalitkin Scheme
Previous Article in Journal
Driver Training Based Optimized Fractional Order PI-PDF Controller for Frequency Stabilization of Diverse Hybrid Power System
Previous Article in Special Issue
Convergence Analysis of a New Implicit Iterative Scheme and Its Application to Delay Caputo Fractional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Seventh Order Family of Jarratt Type Iterative Method for Electrical Power Systems

by
Saima Yaseen
1,†,
Fiza Zafar
1,† and
Francisco I. Chicharro
2,*,†
1
Centre for Advanced Studies in Pure and Applied Mathematics, Bahauddin Zakariya University, Multan 60800, Pakistan
2
Instituto de Matemática Multidisciplinar, Universitat Politècnica de València, 46022 València, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2023, 7(4), 317; https://doi.org/10.3390/fractalfract7040317
Submission received: 20 February 2023 / Revised: 31 March 2023 / Accepted: 4 April 2023 / Published: 6 April 2023
(This article belongs to the Special Issue Applications of Iterative Methods in Solving Nonlinear Equations)

Abstract

:
A load flow study referred to as a power flow study is a numerical analysis of the electricity that flows through any electrical power system. For instance, if a transmission line needs to be taken out of service for maintenance, load flow studies allow us to assess whether the remaining line can carry the load without exceeding its rated capacity. So, we need to understand about the voltage level and voltage phase angle on each bus under steady-state conditions to keep the bus voltage within a specific range. In this paper, our goal is to present a higher order efficient iterative method to carry out a power flow study to determine the voltages (magnitude and angle) for a specific load, generation and network conditions. We introduce a new seventh-order three-step iterative scheme for obtaining approximate solution of nonlinear systems of equations. We attain the seventh-order convergence by using four function evaluations which makes it worthy of interest. Moreover, we show its applicability to the electrical power system for calculating voltages and phase angles. By calculating the bus angle and voltage level, we conclude that the performance of the power system is assessed in a more efficient manner using the new scheme. In addition, dynamical planes of the methods applied on nonlinear systems of equations show global convergence.

1. Introduction

Research in the field of computational mathematics is constantly expanding with the development of new numerical schemes or the modification of old ones. A growing percentage of modern mathematical research is focused on understanding nonlinear phenomena and systems. Numerous areas of science like business and engineering involve these nonlinear phenomena as there is nonlinearity in all physical events. Nonlinear equations naturally control a wide range of phenomena, including fluid and plasma mechanics, gas dynamics, elasticity, relativity, chemical processes, combustion, ecology, biomechanics, economic modeling issues, transport theory, and many others that can be seen in the available literature. Moré [1] presented several nonlinear models and the majority of them were stated in the form of G ( r ) = 0 , for a multivariate vector-valued function G : W R n R n . Also, Grosan and Abraham [2] analyzed the relevancy of the nonlinear systems in neurophysiology, kinematics syntheses, chemical equilibrium, combustion, and economic modeling problems. Furthermore, Lin et al. [3] also expressed the applications of the nonlinear system in transportation theory. As another important application of nonlinear systems is electrical power systems. Different studies are implemented on steady-state conditions of the electrical system and the study of power flow (or load flow) is one of them. Power flow analysis is generally used in the design and analysis of power systems. We also study load flow to verify the productivity, stability, and reliability of electrical power transfer from generators to consumers through the grid system. It is the determination of steady-state conditions of a power system, particularly for power generation and load demand. Computational procedures (numerical algorithms) are required to determine the steady-state operating characteristics of a power system network from the given line data and bus data. Voltages for each bus are determined so, that the line flows and losses can be computed. Buses are divided into three categories: slack, generation, and load buses. Table 1 collects the quantities for each bus type.
Since both voltage magnitude and angles are specified, the slack bus is frequently known as the reference bus. In load flow studies, the active power (P) and reactive power (Q) in power systems is balanced by slack bus. The slack bus acts to partly ensure an equal number of variables and constraints, otherwise the system would be over-determined with more equations than unknowns without a specified slack bus. Load buses are also known as PQ buses where both net real and reactive power loads are specified. Buses are referred to as regulated or PV buses because the voltage magnitude is regulated and the net real power is specified. While just the voltage angle is unknown in PV buses whereas both voltage magnitudes and angles are unknown in PQ buses. Finding the unknown variables in the power systems is the first step in solving problems in power flows. In 1956, the first automated digital solution to the power flow problem was given by Ward and Hale [4]. Ever since power flow analysis was initiated in 1956 [4,5,6], it used lots of numerical methods such as the Gauss method [4], Newton’s method [7] and fast decoupled method [8] for calculating the voltages (magnitude and angle) for a specific load, generation, and network condition. The power flow problems involve the computation of voltage magnitude and phase angle at each bus in a power system under the following conditions:
(i)
The system is in a sinusoidal steady-state with balanced three phase steady-state conditions.
(ii)
Constant, linear and lumped-parameter branches are used to make up the transmission network.
(iii)
Demands at each (load) bus are the specified real (P) and reactive power (Q).
(iv)
The specified real power generation at each (generator) bus excluding one generator bus.
A schematic of the power flow problem is represented in Figure 1.
The non-linear equations are dependent on the voltages, phase angle, and how users are interconnected to demand the substations and the electrical power. Newton-Raphson method is widely used to compute the roots of nonlinear equations.
r ( m + 1 ) = r ( m ) G ( r ( m ) ) G ( r ( m ) ) .
During the 1960s, the Newton-Raphson method became the method that was frequently used to study load flow analysis. Different researchers [7,8] computed the voltages by a employing matrix formulation of (1) and establishing a system of linearized equations. The function G ( r ( m ) ) of (1) for systems is described in matrix form. An electrical power system containing the real power and the reactive power of the network is formed by the partial derivatives of non-linear equations. The matrix is formulated to generate the Jacobian matrix but in the case of large systems, strong computation is required. Taking this into consideration, higher order efficient numerical methods are worthy of interest. Different authors have offered high-order, multi-step methods to find solutions of non-linear equations to increase efficiency and reduce the number of iterations [9,10,11]. However, mostly sixth-order iterative methods have been presented by researchers employing weight functions and parameter approaches [12,13,14,15]. Only a few seventh-order convergent iterative schemes are available, e.g., [16,17,18]. Furthermore, in some cases, the extension of iterative schemes from univariate to multivariate is a challenging task as some univariate iterative schemes are easily extendable to multivariate cases while others can not be extended to the multivariate cases or require special algebraic manipulations.
This paper attempts to provide a new three-step seventh order numerical scheme to solve the nonlinear system of equations associated with an electrical power system while taking into account these challenging characteristics.

2. Derivation of the Scheme

Let us consider the systems of nonlinear equations G ( r ) = 0 for a multivariate vector-valued function G : W R n R n . Then, we can define the divided difference [19,20,21] as:
[ y , x ; G ] j i = ( G j [ y 1 , y 2 , , y i 1 , y i , x i + 1 , , x n ] G j [ y 1 , y 2 , , y i 1 , x i , x i + 1 , , x n ] ) / ( y i x i ) , 1 i , j n ,
where the index j denotes the jth function and the index i represents the nodes. In the development of our scheme, we employed the weight function technique involving divided differences. Our scheme comprises of three steps which are listed below:
v m = r ( m ) G r ( m ) 1 G r ( m ) , z m = v m P u ( m ) G r ( m ) 1 G v m , r ( m + 1 ) = z m S u ( m ) + Q t ( m ) G r ( m ) 1 G z m ,
where
u ( m ) = I G r ( m ) 1 [ r ( m ) , v m ; G ] ,
and
t ( m ) = I G r ( m ) 1 [ v m , z m ; G ] P u ( m ) ,
where P , S , Q : A n × n ( R ) Γ ( R n ) with A n × n be the set of n × n real matrices and Γ ( R n ) be the set of linear operators from R n t o R n . The order of convergence of scheme (2) turns out to be seventh under the conditions on the weight functions stated below.
Theorem 1. 
Let us suppose that G : W R n R n is a Frechet differentiable function in W containing simple root Υ . We consider that G r is continuous and nonsingular at Υ . Then, the convergence is guaranteed and the numerical scheme (2) has seventh-order convergence if the initial guess r ( 0 ) is sufficiently close to the root and the following conditions are satisfied:
P ( 0 ) = I , P ( 0 ) = 2 I , P ( 0 ) = 0 , P ( 0 ) < , S ( 0 ) = I , S ( 0 ) = 2 I , S ( 0 ) = 2 I , S ( 0 ) < , Q ( 0 ) = 0 , Q ( 0 ) = I , Q ( 0 ) < .
Proof. 
Let us consider that e ( m ) = r ( m ) Υ is the error in the m t h iteration. The Taylor’s series expansion of the function G ( r ( m ) ) and its first order derivative G ( r ( m ) ) with the assumption G ( Υ ) 0 leads us to
G ( r ( m ) ) = G ( Υ ) ( e ( m ) + C 2 ( e ( m ) ) 2 + C 3 ( e ( m ) ) 3 + C 4 ( e ( m ) ) 4 + C 5 ( e ( m ) ) 5 + C 6 ( e ( m ) ) 6
+ C 7 ( e ( m ) ) 7 + O ( ( e ( m ) ) 8 ) ) ,
where
C i = 1 i ! G ( Υ ) 1 G i ( Υ ) , i = 2 , 3 ,
and
G ( r ( m ) ) = G ( Υ ) ( I + 2 C 2 e ( m ) + 3 C 3 ( e ( m ) ) 2 + 4 C 4 ( e ( m ) ) 3 + 5 C 5 ( e ( m ) ) 4 + 6 C 6 ( e ( m ) ) 5
+ 7 C 7 ( e ( m ) ) 6 + O ( ( e ( m ) ) 7 ) ) .
G ( r ( m ) ) = G ( Υ ) ( 2 C 2 + 6 C 3 e ( m ) + 12 C 4 ( e ( m ) ) 2 + 20 C 5 ( e ( m ) ) 3 + 30 C 6 ( e ( m ) ) 4
+ 42 C 7 ( e ( m ) ) 5 + O ( ( e ( m ) ) 6 ) ) .
G ( r ( m ) ) = G ( Υ ) ( 6 C 3 + 24 C 4 e ( m ) + 60 C 5 ( e ( m ) ) 2 + 120 C 6 ( e ( m ) ) 3 + 210 C 7 ( e ( m ) ) 4
+ O ( ( e ( m ) ) 5 ) ) .
G ( 4 ) ( r ( m ) ) = G ( Υ ) ( 24 C 4 + 120 C 5 e ( m ) + 360 C 6 ( e ( m ) ) 2 + 840 C 7 ( e ( m ) ) 3
+ O ( ( e ( m ) ) 4 ) ) .
Subsequently, for
v ( m ) = r ( m ) ( G ( r ( m ) ) ) 1 G r ( m ) ,
applying Taylor’s series to (8) we get
v ( m ) = C 2 ( e ( m ) ) 2 + ( 2 C 3 2 C 2 2 ) ( e ( m ) ) 3 + ( 3 C 4 7 C 2 C 3 + 4 C 2 3 ) ( e ( m ) ) 4 + i = 5 7 A i ( e ( m ) ) i
+ O ( ( e ( m ) ) 8 ) ,
where
A i = A i ( C 2 , C 3 , C 6 ) , 5 i 7 .
Next, we apply Taylor’s series expansion to the function
G ( v ( m ) ) = G ( Υ ) ( C 2 ( e ( m ) ) 2 + ( 2 C 3 2 C 2 2 ) ( e ( m ) ) 3 + ( 3 C 4 7 C 2 C 3 + 4 C 2 3 ) ( e ( m ) ) 4
+ i = 5 7 B i ( e ( m ) ) i + O ( ( e ( m ) ) 8 ) ) ,
where
B i = B i ( C 2 , C 3 , C 7 ) , 5 i 7 .
By using (4)–(7), we obtain
[ r ( m ) , v m ; G ] = G ( v m ) G ( r m ) v m r m
= G ( r m ) + G ( r m ) 2 ! ( v m r m ) + G ( r m ) 3 ! ( v m r m ) 2 + O ( e ( m ) ) 3
= G ( Υ ) ( I + C 2 e ( m ) + ( C 3 + C 2 2 ) ( e ( m ) ) 2 + ( 2 C 2 3 + C 4 + 3 C 3 C 2 ) ( e ( m ) ) 3 + O ( e ( m ) ) 4 ) .
Also, we expand u ( m ) = I G r ( m ) 1 [ r ( m ) , v m ; G ] using Taylor’s series expansion
u ( m ) = C 2 e ( m ) + ( 2 C 3 3 C 2 2 ) ( e ( m ) ) 2 + i = 3 6 D i ( e ( m ) ) i + O ( ( e ( m ) ) 7 ) ,
where
D i = D i ( C 2 , C 3 , C 6 , C 7 ) , 3 i 6 .
Also, Taylor’s expansion of function P about the zero matrix
P ( u ( m ) ) = P ( 0 ) + P ( 0 ) C 2 e ( m ) + ( 2 P ( 0 ) C 3 3 P ( 0 ) C 2 2 + 1 2 P ( 0 ) C 2 2 ) ( e ( m ) ) 2 + i = 3 6 E i ( e ( m ) ) i
+ O ( ( e ( m ) ) 7 ) ,
where
E i = E i ( C 2 , C 3 , C 6 , C 7 ) , 3 i 6 .
The next substep
z ( m ) = v ( m ) P ( u ( m ) ) G ( r ( m ) ) 1 G ( v ( m ) ) ,
becomes,
z ( m ) = ( C 2 P ( 0 ) C 2 ) ( e ( m ) ) 2 + ( 2 C 3 2 C 2 2 2 P ( 0 ) C 3 + 4 P ( 0 ) C 2 2 P ( 0 ) C 2 2 ) ( e ( m ) ) 3
+ ( 14 P ( 0 ) C 2 C 3 4 P ( 0 ) C 2 C 3 + 7 P ( 0 ) C 2 3 3 P ( 0 ) C 4 13 P ( 0 ) C 2 3 1 2 P ( 0 ) C 2 3
+ 3 C 4 7 C 2 C 3 + 4 C 2 3 ) ( e ( m ) ) 4 + i = 5 7 F i ( e ( m ) ) i + O ( ( e ( m ) ) 8 ) ,
where
F i = F i ( C 2 , C 3 , C 6 , C 7 , P ( 0 ) , P ( 0 ) , P ( 0 ) , P ( 0 ) , P i v ( 0 ) ) , 5 i 7 .
Using the conditions
P ( 0 ) = I , P ( 0 ) = 2 I ,
(2), becomes
z ( m ) = ( C 2 C 3 + 5 C 2 3 1 2 P ( 0 ) C 2 3 ) ( e ( m ) ) 4 + i = 4 7 H i ( e ( m ) ) i + O ( e ( m ) ) 8 .
H i = H i ( C 2 , C 3 , C 6 , C 7 , P ( 0 ) , P ( 0 ) , P i v ( 0 ) ) , 4 i 7 .
Similarly, we obtain the following expression by expanding Taylor’s series of the function G ( z ( m ) )
G ( z ( m ) ) = G ( Υ ) ( C 2 C 3 + 5 C 2 3 1 2 P ( 0 ) C 2 3 ) ( e ( m ) ) 4 + i = 5 7 H i ( e ( m ) ) i + O ( ( e ( m ) ) 8 ) ) ,
From (9) the weight function S about the zero matrix is given by:
S ( u ( m ) ) = S ( 0 ) + S ( 0 ) C 2 e ( m ) + ( 2 C 3 S ( 0 ) 3 C 2 2 S ( 0 ) + 1 2 C 2 2 S ( 0 ) ) ( e ( m ) ) 2 + i = 3 7 K i ( e ( m ) ) i
+ O ( ( e ( m ) ) 8 ) ) ,
where
K i = K i ( C 2 , C 3 , C 6 , C 7 , S ( 0 ) , S ( 0 ) , S ( 0 ) , S i v ( 0 ) ) , 3 i 7 .
We have the following expression of the operator [ v m , z m ; G ]
[ v m , z m ; G ] = G ( Υ ) ( I + C 2 2 ( e ( m ) ) 2 + ( 2 C 2 C 3 2 C 2 3 ) ( e ( m ) ) 3 + i = 4 6 L i ( e ( m ) ) i + O ( ( e ( m ) ) 7 ) ) ,
where
L i = L i ( P ( 0 ) P ( 0 ) , P ( 4 ) ( 0 ) , C 2 , C 3 , C 6 , C 7 ) , 3 i 6 .
Applying Taylor’s series to t ( m ) = I G r ( m ) 1 [ v m , z m ; G ] P u ( m ) , we get
t ( m ) = ( C 3 + 5 C 2 2 1 2 P ( 0 ) C 2 2 ) ( e ( m ) ) 2 + i = 3 6 L i ( e ( m ) ) i + O ( ( e ( m ) ) 7 ) .
Thus,
Q ( t ( m ) ) = Q ( 0 ) + Q ( 0 ) ( C 3 + 5 C 2 2 ) ( e ( m ) ) 2 + i = 3 6 M i ( e ( m ) ) i + O ( ( e ( m ) ) 7 ) ,
where
M i = M i ( C 2 , C 3 , C 6 , C 7 , Q ( 0 ) , Q ( 0 ) ) , 3 i 6 .
Consequently, the last step
r ( m + 1 ) = z m ( S ( u ( m ) ) + Q ( t ( m ) ) ) ( G ( r m ) ) 1 G z m ,
is given by:
r ( m + 1 ) = ( C 2 C 3 5 C 2 3 + P ( 0 ) C 2 2 ) ( I + S ( 0 ) + Q ( 0 ) ) ( e ( m ) ) 4 + i = 5 7 N i ( e ( m ) ) i + O ( ( e ( m ) ) 8 ) ,
where
N i = N i ( C 2 , C 3 , C 6 , C 7 , P ( 0 ) , S ( 0 ) , S ( 0 ) , S ( 0 ) , Q ( 0 ) , Q ( 0 ) , Q ( 0 ) ) , 5 i 7 .
Applying
P ( 0 ) = 0 , P ( 0 ) < , S ( 0 ) = I , S ( 0 ) = 2 I , S ( 0 ) = 2 I , S ( 0 ) < , Q ( 0 ) = 2 I , Q ( 0 ) = I , Q ( ( 0 ) ) < ,
to (10), we finally have
e ( m + 1 ) = 1 6 ( C 3 C 2 2 + 5 C 2 5 ) ( 96 C 2 2 + S ( 0 ) C 2 2 P ( 0 ) C 2 2 + 24 C 3 ) ( e ( m ) ) 7 + O ( ( e ( m ) ) 8 ) .
This error analysis reveals that the proposed scheme (2) reaches the seventh order of convergence. It completes the proof. □

2.1. Special Cases

Next, we take some special cases of our proposed scheme (2), which are as follows:
Case 1: If we take the weight functions P ( u ( m ) ) , S ( u ( m ) ) and Q ( t ( m ) ) of the following form:
P ( u ( m ) ) = ( a 0 I + a 1 u ( m ) + a 2 ( u ( m ) ) 2 ) 1 ,
S ( u ( m ) ) = ( b 0 I + b 1 u ( m ) + b 2 ( u ( m ) ) 2 ) 1
Q ( t ( m ) ) = c 0 I + c 1 t ( m ) + c 2 ( t ( m ) ) 2 ,
with
a 0 = I , a 1 = 2 I , a 2 = 4 I ,
b 0 = I , b 1 = 2 I , b 2 = 5 I ,
c 0 = 2 I , c 1 = I , c 2 = c 2 ,
then c 2 = 0.5 I , we have a seventh-order scheme, named as F S 1 which is given below
v ( m ) = r ( m ) G r ( m ) 1 G r m , z ( m ) = v ( m ) I 2 u ( m ) + 4 ( u ( m ) ) 2 1 . G r ( m ) 1 G v ( m ) , r ( m + 1 ) = z ( m ) I 2 u ( m ) 5 ( u ( m ) ) 2 1 + 2 I + t ( m ) + 0.5 ( t ( m ) ) 2 . G r ( m ) 1 G z m .
Case 2: If the weight functions P ( u ( m ) ) , S ( u ( m ) ) and Q ( t ( m ) ) are of the following forms:
P ( u ( m ) ) = ( a 0 + a 1 u ( m ) + a 2 ( u ( m ) ) 2 ) 1 ,
S ( u ( m ) ) = b 0 I + b 1 u ( m ) + b 2 ( u ( m ) ) 2
Q ( t ( m ) ) = c 0 I + c 1 t ( m ) + c 2 ( t ( m ) ) 2 ,
along with
a 0 = I , a 1 = 2 I , a 2 = 4 I ,
b 0 = I , b 1 = 2 I , b 2 = I ,
c 0 = 2 I , c 1 = I , c 2 = c 2 ,
then c 2 = 0.5 I , we obtain the following seventh-order scheme namely F S 2
v ( m ) = r ( m ) G r ( m ) 1 G r ( m ) , z ( m ) = v ( m ) I 2 u ( m ) + 4 ( u ( m ) ) 2 1 G r ( m ) 1 G v ( m ) , r ( m + 1 ) = z ( m ) I + 2 u ( m ) + ( u ( m ) ) 2 + 2 I + t ( m ) + 0.5 ( t ( m ) ) 2 . G r ( m ) 1 G z m .
Case 3: If the weight functions P ( u ( m ) ) , S ( u ( m ) ) and Q ( t ( m ) ) are of the following forms:
P ( u ( m ) ) = a 0 + a 1 u ( m ) + a 2 ( u ( m ) ) 2 ,
S ( u ( m ) ) = ( b 0 + b 1 u ( m ) + b 2 ( u ( m ) ) 2 ) 1
and
Q ( t ( m ) ) = c 0 I + c 1 t ( m ) + c 2 ( t ( m ) ) 2 ,
with
a 0 = I , a 1 = 2 I , a 2 = 0 ,
b 0 = I , b 1 = 2 I , b 2 = 5 I ,
c 0 = 2 I , c 1 = I , c 2 = c 2 ,
then c 2 = 0.1 I , we have the following seventh-order scheme named F S 3
v ( m ) = r ( m ) G r ( m ) 1 G r ( m ) , z ( m ) = v ( m ) I + 2 u ( m ) G r ( m ) 1 G v ( m ) , r ( m + 1 ) = z ( m ) I 2 u ( m ) 5 ( u ( m ) ) 2 1 + 2 I + t ( m ) + 0.1 ( t ( m ) ) 2 . G r ( m ) 1 G z m .

2.2. Computational Cost

Two important factors must be considered to evaluate an iterative method’s effectiveness: The number of evaluations of functions ( F E ) and the number of operations (products–quotients) performed at each iteration. Therefore, our aim is to show the effectiveness of the suggested and existing scheme. To achieve this goal, we employ the computational efficiency index C I = p 1 d + o p as well as the efficiency index I= p 1 d defined in [22], where p denotes the order of convergence, d is the number of evaluations performed at each iteration, and op is the number of product-quotient calculations performed per iteration. The number of evaluations of functions for each G , G and first-order divided difference [ . , . ; G ] , at each iteration is n , n 2 and n ( n 1 ) respectively, and these values are used to construct the efficiency index I. Similarly, to determine the computational efficiency index C I , the same computations are required. To calculate a scheme’s efficiency in this manner, the computing effort of each iteration should be considered. When computing the computational efficiency index C I [22], we take into account the fact that using LU decomposition to solve l linear systems require the number of 1 3 n 3 + l n 2 1 3 n products-quotients. In addition, n 2 products are computed in the case when a matrix is multiplied with a vector, and n 2 quotients are needed for the computation of the first-order divided differences operator. In Table 2, the computational efficiency indices of our methods are compared with the seventh-order method presented by Abad et al. [17] named as AC is given below:
v ( m ) = r ( m ) ( G ( r ( m ) ) ) 1 G ( r ( m ) ) , z ( m ) = v ( m ) T ( u ( m ) ) [ r ( m ) , v ( m ) ; G ] 1 G ( v ( m ) ) , r ( m + 1 ) = z ( m ) W ( u ( m ) ) [ v ( m ) , z ( m ) ; G ] 1 G ( z ( m ) ) ,
where
T ( u ( m ) ) = I + u ( m ) , W ( u ( m ) ) = I + ( u ( m ) ) 2 .
Figure 2 and Figure 3 display the computational efficiency index for the existing method and systems of size ranging from 2 to 20 with weight functions F S 1 , F S 2 and F S 3 respectively. Let us notice that both figures demonstrate that the computational efficiency index behavior for the weight functions F S 1 , F S 2 and F S 3 are better for systems of size greater than 10. The term 1 3 n 3 in computational cost is because of the existence of only one type of linear system to be solved at each iteration with the same matrix of coefficients G ( r ) .

3. Numerical Results

3.1. Equations for Load Flow Analysis

The traditional power flow equation is defined in terms of real power ( P ) , reactive power ( Q ) and voltage ( V ) as:
R i = P i + j Q i ,
where P i is the real power injection (positive) or withdrawal (negative) at bus n , Q i is the reactive power injection or withdrawal at bus n and R i is the net complex power injection.
R i = V i I * = V i ( k = 1 n Y i k V k ) * ,
so the power injections can be written as,
R i * = P i j Q i = V i * I = V i * k = 1 n Y i k V k .
The transformation from rectangular to polar coordinates for the complex voltage is given as:
V i = | V i | δ i , V k = | V k | δ k , Y i k = | Y i k | θ i k ,
Hence,
R i * = P i j Q i = ( | V i | δ i ) k = 1 n ( | Y i k | θ i k ) ( | V k | δ k ) ,
R i * = k = 1 n | Y i k V k V i | ( θ i k + δ k δ i ) ,
where θ be the voltage angle and δ be the current angle. Now separating real and imaginary components, we obtain the equations for real and reactive power as follows:
P i = | V i | k = 1 n | Y i k V k | cos ( θ i k + δ k δ i ) ,
and
Q i = | V i | k = 1 n | Y i k V k | sin ( θ i k + δ k δ i ) .
Now, we want to verify the numerical results of our iterative method so, we consider the power flow problem and compare the numerical results of our scheme namely F S 1 , F S 2 and F S 3 with respect to the number of iterations m, the absolute residual error of the corresponding function G r ( m ) , absolute error in two consecutive iterations r ( m ) r ( m 1 ) and computational order of convergence [23] expressed as δ ln G r ( m + 1 ) G r ( m ) ln G r ( m ) G r ( m 1 ) . All numerical results are obtained by using Maple 13 in a PC with specifications: Intel(R) Core(TM) i5 CPU 6300U @ 2.40 GHz, 2.50 GHz (64-bit Operating System) Microsoft Windows 10 Professional and 8 GB RAM with a precision of 300 digits. We compare our method with the seventh-order method presented by Abad et al. [17] named as AC.
Example 1. 
Let us consider 3 bus power system with the generator at bus 1 and bus 3 where the unit for measurement is per-unit (p.u.). The magnitude of the voltage at bus1 is adjusted to 1.05 pu. Voltage magnitude at bus 3 is fixed at 1.04 pu with a real power generation of 2.0 pu. A load consisting of 4.0 pu and 2.5 pu is taken from bus 2. Line impedances can be shown in Figure 4.
As the numbers of buses are three so, Y b u s can be written as:
Y b u s = Y 11 Y 12 Y 13 Y 21 Y 22 Y 23 Y 31 Y 32 Y 33 = 53.851 68 . 198 0 22.360 116 . 56 0 31.6228 108 . 434 0 22.360 116 . 56 0 58.137 63 . 434 0 35.777 116 . 565 0 31.6228 108 . 434 0 35.777 116 . 565 0 67.2309 67 . 249 0
By expanding (15) and (16) we get:
P 2 = V 2 k = 1 3 Y 2 k V k cos ( θ 2 k + δ k δ 2 ) ,
P 3 = V 3 k = 1 3 Y 3 k V k cos ( θ 3 k + δ k δ 3 ) ,
Q 2 = V 2 k = 1 3 Y 2 k V k sin ( θ 2 k + δ k δ 2 ) .
After substituting values of Y b u s in (17)–(19) we get a system of three nonlinear equations with three unknowns given as:
P 2 = 23.478630 V 2 cos ( 2.0344 + δ 2 ) + 26.00005 V 2 2 + 37.20818 V 2 cos ( 2.0344 δ 3 + δ 2 ) + 4 , P 3 = 34.53209 cos ( 1.89254 + δ 3 ) + 37.208184 V 2 cos ( 2.034444788 + δ 2 δ 3 ) + 26.12160874 , Q 2 = 23.47863 V 2 sin ( 2.03444 + δ 2 ) + 52.0000 V 2 2 + 37.20818 V 2 sin ( 2.034444788 δ 3 + δ 2 ) + 2.5 .
The value r ( 0 ) = ( 0 , 0 , 1 ) t is taken as the initial value and the solution of the nonlinear system for the given problem is Υ = ( 0.047062483774 , 0.00870636568 , 0.971677789143 ) t .
Remark 1. 
We have calculated the voltage magnitude and phase angles of the three-bus power system with greater accuracy. The efficiency index (I) and computational efficiency index (CI) [22] of the above example indicate that the four methods present an efficiency index of I = 7 ( 1 / 30 ) . However, they have marginal differences in the computational efficiency indexes, with values 7 ( 1 / 117 ) , 7 ( 1 / 119 ) and 7 ( 1 / 128 ) for the AC, FS 1 -FS 3 and FS 2 methods, respectively. However, these small differences are reflected in the accuracy of the results as evident from Table 3.

3.2. Stability of the Methods

The iterative methods under study will be employed to solve two academic nonlinear systems with the objective of analyzing the proximity of the initial estimate to the solution. In this regard, their corresponding dynamical planes are represented.
Dynamical planes are commonly used to analyze the stability of a method. To solve nonlinear equations, a set of initial estimates z 0 C ^ is selected, representing the real and imaginary parts of z 0 on the abscissae and ordinates, respectively. Each root is assigned a color, and each initial point is represented by the color to whose root it converges. To solve nonlinear systems of two equations, a set of initial guesses x ( 0 ) = x 1 x 2 ( 0 ) is selected, and the representation follows the same rules as in the case of nonlinear equations.
A set of 400 equally spaced initial estimates in the square { x ( 0 ) R 2 : 2 x 1 , x 2 2 } is taken. The convergence is set when x ( k + 1 ) x ( k ) < 10 6 . If this condition is not achieved after 50 iterations, divergence is interpreted. The information about the number of iterations is also included by scaling the brightness of the points: the darkest scaling, the most number of iterations. The roots of the nonlinear system of equations are represented with white stars. The generation of the dynamical planes follows the guidelines given in [24].
The discussion begins with the uncoupled system of nonlinear equations
x 1 2 1 = 0 , x 2 2 1 = 0 ,
whose solutions are ( 1 , 1 ) , ( 1 , 1 ) , ( 1 , 1 ) , and ( 1 , 1 ) . Purple, green, orange and blue colors represent the initial points whose orbit converges to the roots ( 1 , 1 ) , ( 1 , 1 ) , ( 1 , 1 ) , and ( 1 , 1 ) , respectively. Black color represents divergence, and white stars represent the roots. Figure 5 collects the dynamical planes of the methods AC, FS1, FS2 and FS3.
Global convergence can be observed at Figure 5 applying methods AC, FS1 and FS2, with slower convergence for those initial guesses that are close to the axis. Method FS3 has wide convergence except for initial guesses that are close to the axis, since the denominators are very close to zero at these points.
The second system of nonlinear equations to discuss is
x 1 x 2 + x 1 x 2 1 = 0 , x 1 x 2 x 1 + x 2 1 = 0 ,
whose solutions are ( 1 , 1 ) , and ( 1 , 1 ) . Blue and orange colors represent the initial points whose orbit converges to the roots ( 1 , 1 ) , and ( 1 , 1 ) , respectively. Black color represents divergence, and white stars represent the roots. Figure 6 collects the dynamical planes of the methods AC, FS1, FS2 and FS3.
Global convergence can be again observed at Figure 6 applying methods AC, FS1 and FS2, but the line x 1 + x 2 = 0 . This band is broader in method FS3. The darkness of this line is caused by the proximity of the denominator operator to zero, leading to regions of slow convergence or even divergence. The basins of attraction of Figure 6b,c are less intricate than Figure 6a, improving the stability of methods FS1 and FS2 with respect to AC.
Let us remark that both FS1 and FS2 methods allow using values of the initial estimate more distant from the solution.

4. Conclusions

We developed a new seventh-order scheme for a solving nonlinear system of equations using weight functions. We analyzed a power flow problem and calculated phase angle and voltage magnitude using the new schemes. Numerical results of the power flow-based problem are further supported when we compare our scheme with existing families of the same domain [17]. In addition, methods FS1 and FS2 show global convergence, allowing initial guesses further away from the solution.

Author Contributions

Conceptualization, S.Y. and F.Z.; methodology, S.Y.; software, S.Y.; validation, F.Z.; formal analysis, S.Y. and F.Z.; investigation, S.Y. and F.Z.; writing—original draft preparation, S.Y.; writing—review and editing, F.I.C.; visualization, F.I.C.; supervision, F.Z.; project administration, F.I.C.; funding acquisition, F.I.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank the anonymous reviewers for their suggestions and comments that have improved the final version of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moré, J.J. A collection of nonlinear model problems. In Computational Solution of Nonlinear Systems of Equations, Lectures in Applied Mathematics; Allgower, E.L., Georg, K., Eds.; American Mathematical Society: Providence, RI, USA, 1990; Volume 26, pp. 723–762. [Google Scholar]
  2. Grosan, C.; Abraham, A. A new approach for solving nonlinear equations systems. IEEE Trans. Syst. Man Cybernet Part A Syst. Humans 2008, 38, 698–714. [Google Scholar] [CrossRef]
  3. Lin, Y.; Bao, L.; Jia, X. Convergence analysis of a variant of the newton method for solving nonlinear equations. Comput. Math. Appl. 2010, 59, 2121–2127. [Google Scholar] [CrossRef] [Green Version]
  4. Ward, J.B.; Hale, H.W. Digital computer solution of power flow problems. AIEE Trans. 1956, 75, 398–404. [Google Scholar]
  5. Iwamoto, S.; Tamura, Y. Load flow calculation method for ill-conditioned power systems. IEEE Trans. Power Appar. Syst. 1981, 100, 1736–1743. [Google Scholar] [CrossRef]
  6. Stott, B. Review of load-flow calculation methods. Proc. IEEE 1974, 62, 916–929. [Google Scholar] [CrossRef]
  7. Tinney, W.F.; Hart, C.E. Power flow solution by Newton’s method. IEEE Trans. Power Appar. Syst. 1967, PAS-86, 1449–1460. [Google Scholar] [CrossRef]
  8. Stott, B.; Alsac, O. Fast decoupled load flow. IEEE Trans. Power Appar. Syst. 1974, PAS-93, 859–869. [Google Scholar] [CrossRef]
  9. Changbum, C.; Neta, B. A third-order modification of Newton’s method for multiple roots. Appl. Math. Comput. 2009, 211, 474–479. [Google Scholar]
  10. Dzunic, J.; Petkovic, M.S. A family of three-point methods of Ostrowski’s type for solving nonlinear equations. J. Appl. Math. 2012, 2012, 425867. [Google Scholar] [CrossRef] [Green Version]
  11. Soleymani, F.; Vanani, S.K. Numerical solution of nonlinear equations by an optimal eighth-order class of iterative methods. Ann. Univ. Ferrara 2013, 59, 159–171. [Google Scholar] [CrossRef]
  12. Behl, R.; Argyros, I.K. A new higher order iterative scheme for the solutions of nonlinear systems. Mathematics 2020, 8, 271. [Google Scholar] [CrossRef] [Green Version]
  13. Behl, R.; Sarría, I.; González, R.; Magreñán, A.A. Highly efficient family of iterative methods for solving nonlinear models. J. Comput. Appl. Math. 2019, 346, 110–132. [Google Scholar] [CrossRef]
  14. Kansal, M.; Cordero, A.; Bhalla, S.; Torregrosa, J.R. New fourth and sixth-order classes of iterative methods for solving systems of nonlinear equations and their stability analysis. Numer. Algor. 2021, 87, 1017–1060. [Google Scholar] [CrossRef]
  15. Lee, M.Y.; Kim, Y.I. Development of a family of Jarratt-like sixth-order iterative methods for solving nonlinear systems with their basins of attraction. Algorithms 2020, 13, 303. [Google Scholar] [CrossRef]
  16. Amiri, A.; Cordero, A.; Darvishi, M.T.; Torregrosa, J.R. Stability analysis of a parametric family of seventh-order iterative methods for solving nonlinear systems. Appl. Math. Comput. 2018, 323, 43–57. [Google Scholar] [CrossRef]
  17. Abad, M.; Cordero, A.; Torregrosa, J.R. A family of seventh order schemes for solving nonlinear systems. Bull. Math. Soc. Sci. Math. Roum. Tome 2014, 57, 133–145. [Google Scholar]
  18. Behl, R.; Arora, H. A novel scheme having seventh-order convergence for nonlinear systems. J. Comput. Appl. Math. 2022, 404, 113301. [Google Scholar] [CrossRef]
  19. Grau-Sánchez, M.; Noguera, M.; Amat, S. On the approximation of derivatives using divided difference operators preserving the local convergence order of iterative methods. J. Comput. Appl. Math. 2013, 237, 363–372. [Google Scholar] [CrossRef]
  20. Moysi, A.; Ezquerro, J.A.; Hernández-Verón, M.Á.; Magreñán, Á.A. On the set of initial guesses for the secant method. Math. Meth. Appl. Sci. 2023, 2023, 1–14. [Google Scholar] [CrossRef]
  21. Ortega, J.M.; Rheinboldt, W.G. Iterative Solutions of Nonlinear Equations in Several Variables; Academic Press: New York, NY, USA, 1970. [Google Scholar]
  22. Cordero, A.; Torregrosa, J.R. On interpolation variants of Newton’s method for functions of several variables. J. Comput. Appl. Math. 2010, 234, 34–43. [Google Scholar] [CrossRef] [Green Version]
  23. Soleymani, F.; Sharifi, M.; Shateyi, S.; Haghani, F.K. Iterative methods for nonlinear equations or systems and their applications. J. Appl. Math. 2014, 2014, 705375. [Google Scholar] [CrossRef]
  24. Chicharro, F.I.; Cordero, A.; Torregrosa, J.R. Drawing dynamical and parameters planes of iterative families and methods. Sci. World J. 2013, 780153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic of the power flow problem.
Figure 1. Schematic of the power flow problem.
Fractalfract 07 00317 g001
Figure 2. Comparison of Computational Efficiency for 2 n 10 .
Figure 2. Comparison of Computational Efficiency for 2 n 10 .
Fractalfract 07 00317 g002
Figure 3. Comparison of Computational Efficiency for 10 n 20 .
Figure 3. Comparison of Computational Efficiency for 10 n 20 .
Fractalfract 07 00317 g003
Figure 4. Power System.
Figure 4. Power System.
Fractalfract 07 00317 g004
Figure 5. Dynamical planes of solving system (20) with the methods (a) AC, (b) FS1, (c) FS2, and (d) FS3.
Figure 5. Dynamical planes of solving system (20) with the methods (a) AC, (b) FS1, (c) FS2, and (d) FS3.
Fractalfract 07 00317 g005
Figure 6. Dynamical planes of solving system (21) with the methods (a) AC, (b) FS1, (c) FS2, and (d) FS3.
Figure 6. Dynamical planes of solving system (21) with the methods (a) AC, (b) FS1, (c) FS2, and (d) FS3.
Fractalfract 07 00317 g006
Table 1. Bus Classification Used In Power Flow Problems.
Table 1. Bus Classification Used In Power Flow Problems.
Bus TypeFixed QuantitiesVariable Quantities
Slackvoltage magnitude, voltage anglereal power, reactive power
PQreal power, reactive powervoltage magnitude, voltage angle
PVreal power, voltage magnitudereactive power, voltage angle
Table 2. Efficiency Index(I) and Computational Efficiency Index(CI).
Table 2. Efficiency Index(I) and Computational Efficiency Index(CI).
MethodsFE of GFE in [ . , . ; G ] FE in G Total FEICI
A C 321 3 n 2 + n 7 1 3 n 2 + n 7 1 n 3 + 10 n 2
F S 1 321 3 n 2 + n 7 1 3 n 2 + n 7 1 1 3 n 3 + 12 n 2 + 2 3 n
F S 2 321 3 n 2 + n 7 1 3 n 2 + n 7 1 1 3 n 3 + 13 n 2 + 2 3 n
F S 3 321 3 n 2 + n 7 1 3 n 2 + n 7 1 1 3 n 3 + 12 n 2 + 2 3 n
Table 3. Comparison of methods for the load flow problem.
Table 3. Comparison of methods for the load flow problem.
Cases m r ( m ) r ( m 1 ) G r ( m ) δ
F S 1 1 4.706246 e ( 2 ) 1.464301 e ( 6 )
2 2.639187 e ( 8 ) 1.199179 e ( 38 ) 5.10062
3 1.665677 e ( 40 ) 4.493960 e ( 198 ) 7.33718
F S 2 1 4.706248 e ( 2 ) 7.325844 e ( 7 )
2 2.713283 e ( 8 ) 2.067714 e ( 38 ) 4.78636
3 6.219538 e ( 40 ) 1.170039 e ( 196 ) 7.42478
F S 3 1 4.706248 e ( 2 ) 1.539407 e ( 6 )
2 2.497081 e ( 8 ) 2.147565 e ( 37 ) 4.92188
3 3.587427 e ( 39 ) 5.018327 e ( 192 ) 7.40976
A C 1 4.706248 e ( 2 ) 2.850094 e ( 6 )
2 3.482259 e ( 8 ) 1.527893 e ( 36 ) 5.04385
3 2.406074 e ( 38 ) 1.218126 e ( 187 ) 7.37009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yaseen, S.; Zafar, F.; Chicharro, F.I. A Seventh Order Family of Jarratt Type Iterative Method for Electrical Power Systems. Fractal Fract. 2023, 7, 317. https://doi.org/10.3390/fractalfract7040317

AMA Style

Yaseen S, Zafar F, Chicharro FI. A Seventh Order Family of Jarratt Type Iterative Method for Electrical Power Systems. Fractal and Fractional. 2023; 7(4):317. https://doi.org/10.3390/fractalfract7040317

Chicago/Turabian Style

Yaseen, Saima, Fiza Zafar, and Francisco I. Chicharro. 2023. "A Seventh Order Family of Jarratt Type Iterative Method for Electrical Power Systems" Fractal and Fractional 7, no. 4: 317. https://doi.org/10.3390/fractalfract7040317

APA Style

Yaseen, S., Zafar, F., & Chicharro, F. I. (2023). A Seventh Order Family of Jarratt Type Iterative Method for Electrical Power Systems. Fractal and Fractional, 7(4), 317. https://doi.org/10.3390/fractalfract7040317

Article Metrics

Back to TopTop