UDC 621.9:681.54
doi: 10.20998/2074-272X.2016.3.05
A.V. Chaban, V.R. Levoniuk, I.M. Drobot, A.F. Herman
MATHEMATICAL MODEL OF ELECTROMAGNETIC PROCESSES IN LEHERA LINE AT OPEN-CIRCUIT OPERATION
Purpose. The work proposed for the modeling of transients in Lehera line uses a modified Hamilton-Ostrogradskiy principle. The above approach makes it possible to avoid the decomposition of a single dynamic system that allows you to take into account some subtle hidden movements. This is true for systems with distributed parameters, which in the current work we are considering. Methodology. Based on our developed new interdisciplinary method of mathematical modeling of dynamic systems, based on the principle of modified Hamilton-Ostrogradskiy and expansion of the latter on the non-conservative dissipative systems, build mathematical model Lehera line. The model allows to analyze transient electromagnetic processes in power lines. Results. In this work the model used for the study of transients in the non-working condition Lehera line. Analyzing the results shows that our proposed approach and developed based on a mathematical model is appropriate, certifying physical principles regarding electrodynamics of wave processes in long power lines. Presented in 3D format, time-space distribution function of current and voltage that gives the most information about wave processes in Lehera line at non-working condition go. Originality. The originality of the paper is that the method offinding the boundary conditions of the third kind (Poincare conditions) taking into account all differential equations of electric power system, i.e. to find the boundary conditions at the end of the line involves all object equation. This approach enables the analysis of any electric systems. Practical value. Practical application is that the wave processes in lines affect the various kinds of electrical devices, proper investigation of wave processes is the theme of the present work. References 12, figures 12.
Key words: mathematical modeling, Hamilton-Ostrogradskiy principle, Euler-Lagrange equation, electric power system, power line with distributed parameters.
В работе, на основе обобщенного междисциплинарного (интердисциплинарного) метода математического моделирования, основанного на модификации интегрального вариационного принципа Гамильтона-Остроградского, предложена математическая модель двухпроводной длинной линии электропередач, которая работает на холостом ходу. Представлены результаты компьютерной симуляции переходных процессов в виде рисунков, которые анализируются.
Библ. 12. рис. 12.
Ключевые слова: математическое моделирование, принцип Гамильтона-Остроградского, уравнение Эйлера-Лагранжа, электроэнергетическая система, линия электропередач с распределенными параметрами.
Introduction. Mathematical modeling of complex electrical systems today is an important technical problem. With mathematical simulation device can analyze electromagnetic and electromechanical transient processes in electrical facilities and systems, not using for the latest expensive full-scale experiments. No exception here and electricity.
In the current work as an example Lehera system analysis [1] uses a long line with distributed parameters that runs on direct current. We know that these lines have found their proper place in power systems around the world. Transmission of electricity in this way: it reduces losses in the lines (due absence of skin effect phenomenon) makes possible association between local power systems that operate with varying frequency and reduces the cost of construction for large distances due to fewer wires and auxiliary fittings, etc.
Unfortunately, in our country for some reason, however, and economic, in 2014 were brought down the only DC line Volgograd-Donbass, which has been designed for a voltage of 800 kV. But in highly developed foreign countries mentioned lines are not only effective, but commissioned new due to certain advantages mentioned above. Here we can mention the following lines: Line Pacific DC power 1400 MW, ± 400 kV voltage, length of 1362 km for the transmission of electricity from hydropower plants in Oregon grid in Los Angeles; power transmission line HPP «Xiangjiaba» - China's Shanghai ± 800 kV voltage guarantees transmitting 6400 MW over a distance of 2000 km; Canada three transmission line
length of about 900 km, built by HPP Nelson River, located in the Arctic Circle, to the city Winnipeg in the South of country. Epps was the third power of 2000 MW at a voltage of ± 500 kV; Brazil put into operation two chains of the line of HPP Itaipu of throughput of 3,150 MW at a voltage of ± 600 kV. The length of each circuit of about 800 km, and others [2].
Analysis of last investigations. Among a number of scientific papers devoted to the analysis of transient processes in power systems look at some of them, next to the theme of this work.
In [3] developed a mathematical model of two- and three wire power line AC to study transient processes and phenomena of overvoltage in a 500 kV line. Based on the software code ATP-EMTP transients were calculated and investigated the phenomenon over during the emergency state line.
A practical approach in the study of transient electromagnetic processes in [4] is represented. After describing the many cases of simulation modeling for selected items grid requirements are presented. Also, a comparative analysis of studies of transient electromagnetic processes in the correct and incorrect model grid is made.
In [5] the mathematical model of electromagnetic transients in electrical systems that is based on discrete nodal equations in phase coordinates and implicit numerical integration methods, which enables modeling transients with symmetric and asymmetric switching and injuries in electrical networks of any configuration.
© A.V. Chaban, V.R. Levoniuk, I.M. Drobot, A.F. Herman
The work [6] covers a wide range of analysis and derivatives re-established processes in electricity under the original angle. Materials book is based on the classical approach to the modeling of electric-energy-systems. Unfortunately, apparently because of the limitation of the volume of the book, the latter did not present the results of computer simulation of wave processes in lines.
The goal improvement on the base of utilization of variation approaches, methods of mathematical modeling of transients in Lehera line which operates on a non-working pace, and due to this more correctly simulation of wave processes.
Variation model of the Lehera line. To build mathematical models of objects under consideration with a high level of adequacy to properly use basic fundamental laws applied physics, applied in the relevant fields of science [9]. In our case this is electrodynamics [1, 6, 7].
Mathematical modeling is usually using two approaches. The first - a classic approach based on the law of conservation of energy and the second - variation based on minimizing functional of the system [9]. Each of these approaches has its advantages and drawbacks, but when used properly leads to reliable results [8]. In other words, the roads leading to the final model is differ but obtained the result - the same. Usually, choosing the right approach to modeling is a proper of researcher.
We offer analysis of transients in line Lehera use modified Hamilton-Ostrogradskiy principle (variation technique) [9]. The above approach allows you to avoid decomposition of single dynamic system, and to obtain initial state equation energetic exclusively on a single approach, enhanced by constructing Lagrange function [9]. In other words, the proposed way allows you to build dynamic systems models based on interdisciplinary approaches. This is especially true for systems with distributed parameters, and in that long transmission lines, as in the equations of the facility it is necessary to consider: electrostatic effects (arc phenomenon), thermodynamic effects (conductors heating, especially during melting ice) mechanical impacts on wires, in particular, various oscillations (especially resonant and close to resonance (beat fluctuations) processes) and others. In the current work we do not consider the above mentioned effects, but these effects, we plan to further consider our investigations actually for that we offer this approach.
A key element of the principle of modified Hamil-ton-Ostrogradskiy is extended non-conservative Lagran-gian. We present its analytical form [8, 10]:
L* = T* - P* + ®* - d*, (1)
where L* is the modified Lagrange function, T * is the kinetic co-energy, P* is the potential energy, O* is the dissipation energy, D* is the energy of external on-potential forces.
We have already mentioned that the line Lehera generally seen as a system with distributed parameters [10, 11]. Then the elements of the modified Lagrange function will not feature power, and their respective densities [1]. So, functional of action by Hamilton-Ostrogradskiy will have a form [9]:
i2 f 1 S = jl L* + JLldl dt, I = JLldl, here L = 0, (2)
h V i J i
where S is the action by Hamilton-Ostrogradskiy, Ll is the
linear density of the modified Lagrange function, I is the
energetic functional.
We write components of the expanded Lagrange
function (mean linear density) [9]:
cT_ 3x
_ T =
lQ CP * _
8x
_ P, =■
J_
2C0
Q2,Q _8Q = i ; (3)
8t
CO' 8x
_O, =O, 3-OlB =
R
0 n2 go
2C
2 ^x
dz , (4)
y\t=z
where i(x, t) is the current in the line, R0, go, C0, L0 are the line parameters, is the external energy dissipation, is the internal energy dissipation, Q(x,t) is the charge of the line.
A ! ;,
i(0 j\-/ iy-> .Aif1" L«dx
A'».,
is =0
m i
+ ' /jh'ïj
a fi -t- %Mrl y \ Y q
&,<** I
dt
X
/
Fig. 1. Electric circuit of the Lehera line at open-circuit operation
It is important to note that in equation (4) minus sign appears! This is due to the fact that the function external dissipation depends on the leakage current that flow between the line wires. Obviously, the electric transmission line during the transfer of energy from the source to the consumer consumes the energy dissipated in in space. In other words, energy is transferred exclusively via the electromagnetic field lines and wires only indicate the direction of electromagnetic wave propagation [1].
Taking into account the equations (3), (4) the energy functional will look like [9]:
I =
lL0 Qt2 -
i
2C0
-jlRo Qt2 - 202 Qx2
2C2
dZdl. (5)
At=z
We write the variation of the energy functional (5) and equate it to zero
if t \
l iv l
LoQt + RojQt\t=Tdr o
1 t
- Qx + jO j Qxit=Tdz
SQt -
A
(6)
oo
Qx \dl = 0
Next, for each element of integrand expressions we use the rule of integration by parts, also known as Gauss-Ostrogradskiy theorem. Then, for the first bracket will be [9]:
2
o
-Jf \L0Qt + R j Q,
dr
t\t=r
0 J
dQdl + Qt
(7)
and for the second one:
(
fA ± Q +12 J q
Jar cQx + c2JQx
rdr dQdl + Qx , (8)
i V " ~u 0 J
where Qt, Qx are the boundary conditions for the functional (5).
From here we can write
d=fU-^+12 fa 2Q
J C n2 J
l
C0 dx2 C02 0 dx2 |t=r
dr - L
0
a 2q at2
(9)
-R0 — J— dr \dQdl + ^ = 0, ^ = +^x.
at 0 at t=r
It is easy to see that variation of energy functional can be zero only when equality to zero of integrand or variations function of the charge of the line. As 8Q never can be equal to zero [9, 12] the energy functional (9) obtains a stationary value only in the case when integrand equals zero, i.e. at presence of the Euler-Poisson equation [9, 12]
_L+10 fa2Q c ax2 c,2 J --2
ax2 dr-L-R0? = 0. (10)
1 a2q r a2q d --y = ^^tt+R
aQ
C ax2 u at2 u at
Taking into account the expression [1]
-TT. (11)
j
J at2 It=r
dr=aQ-, fat
at at
dr = Q
(12)
we obtain finally the commonly known telegraph equation [1, 10]
^ = L0C0 ^ + (R0C0 + g0L0 ^ + g0RQ . (13) ar2 at2 at
The equation of the Lehera line (telegraph equation) is written for the function of the charge of the line. However, it is easily transformed to the common telegraph equation
a2i
a2i
^=LC^+((C + gnL) + gAA*=(Qui). (14)
ax2 at2 at
Experience shows that for more optimal description of physical processes in the line it is useful to use as a general function a function of voltage, i.e. 1 = u(x,t) [1, 10].
We rewrite (14) in such a way:
( A
i Li
a 2u ax2
-(( L0 + C0R0 V - g0Ru
1=(15)
The most important problem solving equations (15) is to determine the initial (v(x,t)|i=to) and boundary (u(x,t)|x=0 and u(x,t)|x=l) conditions. As to the first, then the problem is solved in the accustomed way (they calculate their from previous research or take zero). The main prob-
lem is to find the boundary conditions. In general, the voltage at the beginning u(x,t)|x=0 and at the end u(x,t)|x=l. Of the line are unknown. In the particular case (on the current work) voltage is known at the beginning of the line, while at the end of the line - no. Actually finding this voltage we loan.
We write equations (11) in such a way (taking into account Qx(x,t) = C0u(x,t)):
-M^ = R0i(x,t)+ L ^. (16)
ax at
Further, for the power transmission line we write equations (15), (16) in the discrete space (we discretize them by the line method)
dVj dt
(C0L0 )-1
j-1 - 2uj + uj +1
10L0 +
(Ax)2 V"0 0 (17)
+ C0R0)vj -g0R0uj )i u1 = u(x,t)x=0, uN = u(x,t)x=l
di,
uj+1 uj-1 2Ax
= Rij + L0-L; (18)
dU :
-j- = Vj, j = 2,..., N -1. (19)
dt
We rewrite equations (17), (18) for the N-th node of discretization in the correspondence with Fig. 1 in such a form:
dv
N
1
We write for the equation (10) an expression of steady-state connections [1, 10]
dt C0L0
1
(Ax)2
(uN-1 - 2uN + uN+1
- (0L0 + C0R0 )vN - g0R0uN ]
dN+1
- u
N-1
2Ax
= 0,
(20)
(21)
where uN+1 is the discretization node voltage function at the fictitious layer [10] which will be found from the equation (21). Then,
uN+1 = uN-1 .
(22)
Taking into account (20) and (22) we write the final equation of the long line for the N-th node
dvN dt
2
C0 L0 (Ax)
2U N-1
2
V C0 L0 (Ax)2
10 R0 Cn L0
g0 L0 + C0 R0 C0 L0
N ■■
du n dt
■ = vN
(23)
(24)
An important functional dependence that is interesting for the potential users is current value in elements of the Lehera line. Its calculation is possible in such a way. Discretizing equation (16) by the line method (right derivative) we have:
u,+1 - u, Ax
= R^j +
A dt
(25)
From here we obtain finally
di
it=frj- u'+' >-70 'j,j=1.....N - L (26)
Compatible integration is subject to this system of differential equations: (17), (19), (23), (24), (26).
t=r
u
N
Computer simulation results. Computer simulation is carried out for the Lehera line at DC at open-circuit operation. The line has the following parameters: Ro = o.86-io-1 Q/km, Lo = o.134-io-12 H/km, Co = o.85-io-8 F/km, go = o.375-io-7 S/km, length of the line l = 6oo km. The line is supplied by the DC voltage w(x,t)|x=o = 4oo kV.
In Fig. 2-4 spatial distribution of the electromagnetic wave as functional dependences of currents (i) and voltages (2) is presented. From these Figures we can see physical basics of electromagnetic processes in the long Lehera line. Let us analyze these processes.
Fig. 2. Distributions of current (1) and voltage (2) in the line at t = 0.002 s
0 0.04 0.08 0.12
Fig. 5. Representation of the voltage transient function in the central point of the line
Recall that although the line is at open-circuit operation, leakage currents and currents in cell lines will be present. Actually the reason - is capacitive currents between the wires line. Obviously, at the end of the transmission line current will equal zero because the line is unloaded.
Fig. 3. Distributions of current (1) and voltage (2) in the line at t = 0.004 s
Fig. 2 shows the spatial distribution of functions of current and at time 0.002 s. Analyzing the mentioned Figure it is easy to see that the function of voltage begins to decrease, and the central line rises sharply upward. A stream function in the same place on the contrary - falls.
Fig. 6. Representation of the current transient function in the central point of the line
Fig. 3 shows the same, but in time 0.004 s. If at time 0.002 s (see Fig. 2) function of voltage increased to 460 kV but at time 0.004 s this growth was 750 kV. Voltage has almost doubled. As for the current, they fell by almost four times.
Fig. 4. Distributions of current (1) and voltage (2) in the line at t = 0.1 s
Fig. 7. Representation of the voltage transient function at the end of the line
0 0.04 0.08 0.12
Fig. 8. Representation of the current transient function at the end of the line
Fig. 4 shows again the same as in Fig. 2 and 3 at the time when the transition process is almost completed. Of the Figure shows that deviation of functions of voltage and current almost took a minimum value. In other words, the amplitude of the electromagnetic wave due dissipative process significantly decreased. Oscillation process practically attenuates.
600 x, km
Fig. 11. Temporal-spatial distribution of the current function at t e [0; 0.02] s
Fig. 9 and 10 represent the line voltage as a function of time and spatial coordinates. These figures are presented in 3D format. Notably relatively high information content of these figures, which is that the spatial and temporal distribution of creating three-dimensional space. It is advisable to analyze the Figures mentioned in comparison with Fig. 2-5 and 7.
600^0 x, km
Fig. 9. Temporal-spatial distribution of the voltage function at t £ [o; o.o2] s
Fig. 5-8 show transient functional dependences of voltage and current (temporal distribution). The first two Figures concern central node of the line for voltage and the central segment of the line for line. The second two Figures - the penultimate node of the line and penultimate discrete circuit of the line.
Through a comparative analysis of the above Figures it is easy to see that the function of voltage (see Fig. 5 and 7) changes little. See quite a different picture (see Fig. 6 and 8) regarding currents. The current changes almost 8 times. This is because the line power transmission line is (open-circuit operation).
x, km
Fig. io. Temporal-spatial distribution of the voltage function at t £ [o.o2; o.Q4] s
Fig. 12. Temporal-spatial distribution of the current function at t e [0.02; 0.04] s
Fig. 11 and 12 show the same as the previous two, but the function of current. As we see the function of voltage and current are in opposite phase. Because the nature of the stress associated with the electric field and the magnetic current we can make a conclusion of space perpendicular fields E and B which maintains the classic electrodynamics [1]. The presented figures it is advisable to analyze in comparison with Fig. 2-4, 6, 8.
Conclusions.
1. Variation approaches to modeling processes in long power lines make it possible to avoid the decomposition of a unified system, while the final form of the equation of state exclusively from single energy approach by building Lagrange expanded function.
2. An important point in the solution of differential equations of state of long line is a search of boundary conditions that often is veiled, incorrectly set, and the use of boundary conditions Neumann and Poincare boundary conditions. Finding these conditions entails full engagement of system of differential equations of studied object including transformers, reactors, compensation devices, etc., which greatly complicates the calculation of transients in a long line.
3. Experience shows that during the analysis of local power systems as the best option long line the telegraph equation it is advisable to write in a function of voltage. In the case of modeling of local energy systems where they use electromagnetic model elements of these systems
(¥ - type and A - type) have difficulty with utilization of known method of nodal voltages which makes impossible to determine voltage at the beginning and end of the line and therefore it is impossible to correctly solve the equation. All this calls into question the degree of adequacy o eventual results that are obtained by known engineering program Mathematica, MatLab, etc., a especially the use of these programs becomes impossible when considering the circuit-field model elements. In this case, each actual task we must use appropriate apparatus of mathematical modeling.
4. Based on the results of computer simulation we can make a number of conclusions:
• voltage function has the greatest amplitude of oscillations at the end of the line when the current function - at the beginning of the line;
• spatial distribution of functional dependence of the line of sending (Fig. 2, 3 and 4) confirms the physical principles regarding electrodynamics of wave processes in long lines of power supply;
• presented in 3D format temporal-spatial distribution of functions of current and voltage provides most information on wave processes in Lehera line at open-circuit operation.
Materials of this work will be used in further studies that will cover long three-phase power lines with various kinds and types of loads.
REFERENCES
1. Shimoni K. Teoreticheskaya elektrotekhnika [Theoretical Electrical Engineering]. Moscow, Mir Publ., 1956. 773 p. (Rus).
2. Available at: https://www.energetika.in.ua (Accessed 25 July 2014).
3. Nayir A. Simulation of transient processes on overvoltage in electric transmission lines using ATP-EMTP. Turkish Journal of Electrical Engineering & Computer Sciences, 2013, no.21, pp. 1553-1556. doi: 10.3906/elk-1108-8.
4. Sowa P., Kumala R., Luszcz K. Modeling of power system components during electromagnetic transients. International
Journal of Innovative Science, Engineering & Technology, 2014, vol.1, iss.10, pp. 715-719.
5. Vepryk Yu.N., Minchenko A.A. Switching overvoltage in 750 kV power line. Elektrotekhnika i elektromekhanika - Electrical engineering & electromechanics, 2009, no.4, pp. 17-20. (Rus). doi: 10.20998/2074-272X.2009.4.04.
6. Kyrylenko O.V. Matematychne modelyuvannya v elektroen-erhetytsi [Mathematical modeling in the power]. Lviv, Lviv Polytechnic National University Publ., 2010. 608 p. (Ukr).
7. Neyman L.R., Demirchyan K.S. Teoreticheskie osnovy elek-trotekhniki. V 2-kh t. T. 1 [Theoretical bases of electrical engineering. In 2 vols. Vol. 1]. Leningrad, Energoizdat Publ., 1981, p. 536. (Rus).
8. Uayd D., Vudson G. Elektromekhanicheskoye preobrazova-niye energii [Electromechanical energy conversion]. Leningrad, Energiia Publ., 1964. 539 p. (Rus).
9. Chaban A.V. Pryntsyp Hamiltona-Ostrohradskoho v elek-tromekhanichnykh systemakh [The principle of Hamilton-Ostrogradskii in electromechanical systems]. Lviv, Taras Soroka Publ., 2015. 488 p. (Ukr).
10. Chaban A.V Matematychne modeliuvannia kolyvnykh prot-sesiv v elektromekhanichnykh systemakh [Mathematical modeling of oscillating processes in electromechanical systems]. Lviv, Taras Soroka Publ., 2008. 328 p. (Ukr).
11. Zeveke G.V., Ionkin P.A., Netushil A.V., Strakhov S.V. Osnovy teorii tsepey [Fundamentals of circuit theory]. Moscow, Energiia Publ., 1975. 752 p. (Rus).
12. Vasydzu K. Variatsionnyye metody v teorii uprugosti i plas-tichnosti [Variational methods in the theory of elasticity and plasticity]. Moscow, Mir Publ., 1987. 542 p. (Rus).
Received 22.01.2016
A.V. Chaban', Doctor of Technical Science, Professor,
V.R. Levoniuk1, Assistant Lecturer,
I.M. Drobot1, Senior Instructor,
A.F. Herman1, Senior Instructor,
1 Lviv National Agrarian University,
1, V. Velykoho Str., Dubliany, Lviv Region, 30831, Ukraine,
e-mail: atchaban@gmail.com, bacha1991@ukr.net
How to cite this article:
Chaban A.V., Levoniuk V.R., Drobot I.M., Herman A.F. Mathematical model of electromagnetic processes in Lehera line at open-circuit operation. Electrical engineering & electromechanics, 2oi6, no.3, pp. 3o-35. doi: io.2o998/2o74-272X.2oi6.3.o5.