Electrodynamics of Reactive Power in the Space of Inter-Substation Zones of AC Electrified Railway Line

In railway traction, the definition of “electromagnetic field” is functionally connected to the concept of the reactive power consumed by the electric rolling stock, and characterized by the running and standing electromagnetic waves in the space of the inter-substation zones from the site of the AC traction system. Such a definition is established and theoretically justified by the theory of electromagnetic fields. This article uses the methodology of this theory, in particular, a method for power balance estimation in electromagnetic fields based on Maxwell’s equations, as well as methods for the analysis of running and standing electromagnetic waves based on the theory of reflection, propagation and transmission of plane harmonic waves. The research considers the regularities of standing electromagnetic waves in the space of inter-substation zones of electric traction systems, which occur due to the incomplete reflection of incident waves from the contact wire and metal parts of the roof surface and the frontal part of the body of the electric rolling stock. The flow of electricity to the roof surface and the frontal part of the body of an electric locomotive is considered. The possibility of using existing methods to reduce wave reflections and thereby to effectively compensate for reactive power in the space of inter-substation zones is discussed.


Introduction
The fundamental role of reactive power is widely recognized in the electric power industry, and the ambiguity of the concept and formulae for determining reactive power is well known. In five years from now, we will celebrate the 100th anniversary of the "reactive power" concept [1,2], which has brought about a long-lasting discussion that has continued to the present day; let us review some of them.
The works of Constantin J. Budeanu [2] and Stanisław S. Fryze [3,4] were the initial studies at the first stage of the development of the reactive power definition, which lasted until the 1980s. Comprehensive review, classification and analysis of the methods developed during this period, and the names of their authors, are given in monographs by Leszek S. Czarnecki [5] and Volodymyr E. Tonkal [6].
The beginning of the second stage in the development of the power theory was laid in [7][8][9], proved by Hirofumi Akagi et al. in 1982. The main direction in this work is to obtain mathematical relationships and develop algorithms for the control systems of active filters when solving problems of compensating the reactive component of the first harmonic of the current and suppressing its higher harmonics. As a result, the author proposed the so-called "generalized theory of the instantaneous reactive power in three-phase circuits", which is frequently known as the "p-q theory". Unfortunately, in this theory, the authors do not explain reactive power's physical meaning, considering it as a calculated value when analyzing three-phase three-wire systems without a neutral wire. The improved p-q theory was proposed in [10], in which the limitations of its original version were overcome.
As opposed to the p-q theory, based on a plane, the power theory, based on the use of a spatial Cartesian coordinate system, was proposed by Fang Z. Peng et al. in [11] for the three-phase four-wire systems. This theory is applicable for all the following conditions of operation: sinusoidal or non-sinusoidal, as well as balanced or not.
Further development of the previous theory is given in the studies [12,13]. Considering the essence and degree of generalization of the obtained version, the authors proposed calling it the "p-q-r theory". The supply voltage and consumer currents are converted into rotating rectangular spatial pqr coordinates, where the reactive powers q q and q r are calculated as scalar and vector products of voltage and current vectors. A comparative analysis of the considered three power theories is given in the review article [14].
In recent years, the actual publications, including [15][16][17], state the precise physical foundations of the fundamental theory of instantaneous active and reactive power.
Compared with previous power theories, the mathematical relationships of which are intended to determine the instantaneous active and reactive power in three-phase circuits with valve converters, the expressions of instantaneous reactive power in [18][19][20] are proposed together with a physical interpretation. It allowed determining the reactive power in single-phase circuits and individual phases of three-phase systems, including circuits equipped with converters.
From the history of reactive power development, it is seen that the main subject of long-lasting discussions and scientific research is the ambiguity of the concept and formulae for its definition in electrical circuits with non-sinusoidal voltages and currents. From the time of the first work of C. J. Budeanu and up to the 1980s, such ambiguity has already led to the emergence of more than ten methods for determining and calculating reactive power, proposed by many authors. Moreover, in each of these methods, various researchers proposed their approaches to assessing reactive power. For a simple comparison, reactive power was calculated for the DE1 electric locomotive (Ukraine) according to crucial theories. The calculation results are shown in [20] and Table 1. The below theories and expressions were used: -by S. S. Fryze Q F = S 2 − P 2 , -by differential -by integral where U-pantograph voltage, I-locomotive current, k-harmonic number, ϕ-phase shift, S-total (apparent) power, P-active power, S-total (apparent) power, Q-reactive (inactive) power, Q d -the sum of the reduced reactive powers of all of the circuit elements whose voltage and current have the same harmonics, Q i -the sum of the reduced reactive powers of all of the circuit elements whose voltage and current have different harmonics (according to O. A. Maevski).
The different values of reactive powers in Table 1, as well as the different signs, show clear evidence of the ambiguity of the reactive power defined by these theories. In our opinion, it has not been excluded because of numerous reasons and, perhaps, the most important of them is a state of the problem of "reactive power", of which the solution is based on the "circuit" approach, that is, on the theory of electric and magnetic circuits. However, the theory of circuits is known to historically originate from the nucleation and development of a more complex fundamental theory of electromagnetic fields. Many physical concepts, such as the theory of electric and magnetic circuits, like most laws, are based precisely on the formulas and postulates of the field theory. The "circuit" principles for the analysis of processes in electrical devices are widespread, but also provide numerous obligatory assumptions and simplifications of electromagnetic processes, which determine their nature and mechanisms only approximately, and in some cases give erroneous results. Moreover, electric power processes in all electrotechnical devices and systems are the processes of creation, transformation and propagation of electromagnetic fields and the interaction of fields with electric charges. Therefore, to solve electric power problems in electric transport systems, a "field" approach is needed. This approach is based on the fact that electromagnetic energy is transmitted from traction substations to the electric rolling stock (electric locomotives, electric trains, trams, metro trains) of electric transport systems through the airspace of inter-substation zones, by means of electromagnetic waves (i.e., by the movement of the electromagnetic field) rather than through the traction network (i.e., catenary system [21] and rails [22]). In this case, the amount and direction of the transmitted energy are determined by the Poynting vector → S , known from the theory of electrical engineering and radio engineering [23], as follows: H-the values of the vectors of the strengths of electric and magnetic fields at points within the space of inter-substation zones.
In this case, the incident (direct) running electromagnetic waves fall at the surface of the contact wire, rails and body of the electric rolling stock, and form reverse (reflected) waves, which play a double role in the consumption of active and reactive power. Firstly, the flow of active energy decreases, since it is equal to the difference between the energy fields of the incident and reflected waves. Secondly, when combined with the incident wave under certain conditions, the reflected waves cause the appearance of standing waves, which are a characteristic feature of reactive power. Classical approaches [23][24][25] describing the processes of the complete reflection of incident waves, and thus the formation of standing waves, which currently prevail in theoretical electrical engineering, are limited to the "circuit" approach (according to the theory of electrical circuits) and mainly to the analysis of the operating modes of a long line in the limiting cases of its load, as follows: idle, short circuit, reactive load, etc. To solve the problem set in this article, a "field" approach must be applied, in which the feeder or substation zones of the electrified section are considered simple conducting systems with special properties for the propagation of electromagnetic waves from traction substations to the electric rolling stock, which is associated with currents flowing in the catenary system and rails and, finally, the charges arising on them.
A review of the existing investigations has shown that there are currently no scientific publications on this issue. Some exceptions are articles [26][27][28][29], which, using elements of the electromagnetic field theory, make a physical explanation for the occurrence of reactive power Q and distortion power D in electrical circuits with nonsinusoidal voltages and currents. They mostly use the well-known expression of the Umov-Poynting theorem in its complex form to express the balance of powers over volume V of the periodic electromagnetic field, and conclude that the real part of this expression describes active power P and the imaginary part describes reactive power Q. Finally, the authors obtain expressions for the E, D, H and B quantities of the electromagnetic field using EMF, voltage and current in an external nonlinear electric circuit. Ultimately, the authors adhere to the concept of C. J. Budeanu that the total power is determined by the following ratio: thus, the processes of wave propagation of energy in the electromagnetic field are not considered in these works.
The need for a physical interpretation of the inactive components of the total power using the Poynting vector is also emphasized by the results of the research presented in [24]. The author correctly describes the mechanism of the flow of electricity that propagates in the dielectric space around the wires of a three-phase power line. According to this research, the electrical energy of running electromagnetic waves, dissipated in the conductors of the supply line, is distributed over the phases of the load or accumulated and returned by inductions or capacitances. The author divides these waves into two groups. The first group consist of waves with a unidirectional flow of energy from source to load. These waves, making oscillations (called "active"), have a nonzero average value of the transferred energy. The waves of the second group, making oscillations (called "inactive"), do not transfer energy. They spread in the dielectric space around the supply line conductors.
As well as in the previous research, research [30] states that a plane electromagnetic wave is characterized by two flows of energy and two, active and reactive, densities. As these energy flows are inseparable at the phase shift angle ϕ = 0. If ϕ = 0, the reactive components of the energy flow will disappear. The author mathematically and physically substantiated the algorithm for separating the active and reactive components of the energies of electromagnetic waves, and expressed them in the following way: where Π a and Π r active and reactive components of Poynting vector, A-power, αattenuation coefficient, β-phase coefficient, ω-angular frequency, ρ-wave impedance, z-coordinate in the direction of which the Pointing vector is moved. The importance of such division of the energy of electromagnetic waves is shown by the author in an example of the propagation speed of wave energy in a medium with losses. Quite an original consideration of the components of the total power and their balance is shown in [31]. This article proves that in order to balance complex power, the following applies: that is, the balance of this expression for both sides, time t − js must be considered. This is due to the fact that there cannot exist a balance of reactive energy in the "active" (real) time. This balance must be considered at "reactive" time t − js. The authors try to show the stated idea on the example of time-dependent, i.e., parametric, RLC circuits. The fulfilment of the power balance in the electromagnetic field is also studied in the processes of dissipation and absorption of the energy of electromagnetic waves. In particular, a finite-size coated sphere scatterer is considered in [32], which is illuminated by incident electromagnetic waves. The percentage of absorbed and scattered energies determined by the Poynting vector is analyzed. The study found that the absorbed energy, called active energy, is about 25% of the energy of the incident electromagnetic field. It is stated that the imaginary part of the complex power balance expression is the dissipated power, and its occurrence is due to the phase shift between the incident electromagnetic field and polarization currents. Finally, it concludes that the diffuser must be optimally designed to eliminate the reactive power.
From the above analysis follows that the authors of the publications, recognizing the presence of active and inactive components of the energy of electromagnetic waves, do not consider the processes of wave reflection, and do not associate the physical interpretation of reactive power with the emergence and suppression of standing electromagnetic waves in the dielectric space of the electromagnetic field.

Objective of the Article
The article continues and develops the research presented in [33], dedicated to the "field" approach to substantiate the nature and processes of electricity transmission from traction substations (TS) to the electric rolling stock (ERS).
The primary purpose is to establish and provide theoretical grounds for the "field" functional connection between reactive power consumed by the electric rolling stock, and the occurrence and damping of standing electromagnetic waves in the space of intersubstation zones of electric transport systems, based on the electromagnetic field theory.

1.
For the first time in the Ukraine and Europe, a "field" (i.e., based on the electromagnetic field theory) interpretation of functional connectivity of reactive power, transmitted from a traction substation to the electric rolling stock, with occurrence and dissipation processes of standing electromagnetic waves in inter-substation zones of electric transport systems, is hereby proposed and substantiated.

2.
Dependencies have been found to describe the occurrence of standing waves in the space of inter-substation zones under the incomplete reflection of incident electromagnetic waves from a contact wire, rails, metal and dielectric surfaces of roofs, front and side parts of the body of the ERS.

3.
For the first time, existing methods for the suppression of standing electromagnetic waves are adapted to compensate the reactive power in electric transport systems by applying single-or multi-layered film coatings, with specific electromagnetic properties, on the surfaces of feeder zone devices.

Reactive Power as a Measure of Asymmetry of the Rates of Change of Electric Power in Electrical and Magnetic Fields
Let us consider the starting point of a 560 km long railway line supplied by a 25 kV AC 50 Hz traction power system. Imagine that the space of the specified line assumed to be limited by a traction substation, catenary system, railway vehicle and the rails, has a volume V limited by the enclosing surface A v , as shown in Figure 1. The trains in this section are driven by single-phase AC electric locomotives, types VL80 and VL85 (Ukraine). The locomotives have thyristor power converters supplying the traction motors. The frequency of voltage higher harmonics on the output changes are in the range of 1250 to 1950 Hz, according to [34,35].

Reactive Power as a Measure of Asymmetry of the Rates of Change of Electric Power in Electrical and Magnetic Fields
Let us consider the starting point of a 560 km long railway line supplied by a 25 kV AC 50 Hz traction power system. Imagine that the space of the specified line assumed to be limited by a traction substation, catenary system, railway vehicle and the rails, has a volume limited by the enclosing surface , as shown in Figure 1. The trains in this section are driven by single-phase AC electric locomotives, types VL80 and VL85 (Ukraine). The locomotives have thyristor power converters supplying the traction motors. The frequency of voltage higher harmonics on the output changes are in the range of 1250 to 1950 Hz, according to [34,35]. -enclosing surface, ⃗ and ⃗ -vectors of magnetic field tensions in the space of catenary and rails, ⃗ and ⃗ -normal and tangential vectors of electric field tension ⃗ , ⃗ -vector of electric field tension at the point , and -currents in catenary and rails, ℎ-distance between the contact wire of catenary system and rails, -coordinate axis for height.
Let this space be linear, homogeneous and isotropic with the following electromagnetic properties: permittivity , conductivity and permeability . Then, the alternating electromagnetic field in this space could be described by Maxwell's equations in its complex form according to [23], as follows: where ⃗ and ⃗ are the vectors of electric and magnetic field tensions in the space of inter-substation zones. Since both tensions ⃗ and ⃗ are created by a single energy source, i.e., by a traction substation, electricity is transferred according to Poynting vector , and takes place in one direction from TS to ERS.
Like complex power in the theory of circuits, we form a complex Poynting vector ( Figure 1) as follows: where Let this space be linear, homogeneous and isotropic with the following electromagnetic properties: permittivity ε, conductivity σ and permeability µ. Then, the alternating electromagnetic field in this space could be described by Maxwell's equations in its complex form according to [23], as follows: where → E and → H are the vectors of electric and magnetic field tensions in the space of inter-substation zones.
Since both tensions → E and → H are created by a single energy source, i.e., by a traction substation, electricity is transferred according to Poynting vector S, and takes place in one direction from TS to ERS.
Like complex power S in the theory of circuits, we form a complex Poynting vector ( Figure 1) as follows: where rot Next, let us transform system (1)-(4) as follows: multiply expression (4) by → E, and (2) by → H . Then the following expressions could be written: Energies 2021, 14, 3510 Subtracting (6) from (5) and considering that the difference of the left parts of (5) and (6) is equal to −div → S , we have the following: Then the flow of the Pointing vector through the enclosing surface A V into the volume V of inter-substation zones ( Figure 1) will be equal to the following: By analogy with the theory of circuits and according to the Joule-Lenz law in its differential form [23], the first component on the right side of Expression (8) is the average value of heat losses in volume V of inter-substation zones over the period of the supply voltage, and the second component, i.e., the imaginary part, is the average value of reactive power Q entering into volume V through the enclosing surface A V , as follows: where E and H are the effective values of the harmonic tensions. Given that , and considering Formula (9), the expression of the energy of the "field" (by the electromagnetic field vectors) for the representation of reactive power can be written in the following form: Therefore, reactive power in the space of inter-substation zones by volume V is a measure of the asymmetry of the rates of change in the electric and magnetic components of the energy of the electromagnetic field. Numerically, it is equal to the difference between average energies stored in the magnetic and electric fields in a given volume V, multiplied by 2ω.

Volume of Reactive Power in the Space of Inter-Substation Zones of Electrified Section
Using Expressions (9) and (10), let us consider the reactive power Q spreading in the space of an inter-substation zone consisting of seven sections with a total length of 560 km. Each inter-substation zone has a one-sided power supply via the catenary of PBSM-95+MF-100 type and the rails of R65 type, which could be specified as follows: PBSM-95 catenary wire-bimetal steel copper wire with a cross-section of 95 m 2 ; MF-100 contact wire-solid copper wire with a cross-section of 100 m 2 , R65 rails-carbon steel rails with a specific weight of 64.88 kg/m [36]. The VL85 electric locomotive with a train is moving along the inter-substation zone at 160 km/h and taking power P = 5369 kW. According to the experimental research of this study, at this speed, the voltage at the pantograph is U = 26.12 kV, the current in the contact wire and the electric locomotive is I c = 208.9 A, and in the rails of the one-side AC traction system is I r = 83.6 A.
Let us define the tensions of the electric and magnetic fields in the space of the inter-substation zones.
The tension of the electric field in a dielectric (i.e., in the air) has the following two components: tangential E τ and normal E n parts, which are shown in Figure 1. The tangential component E τ in the air of inter-substation zones, according to the "conductor- dielectric" boundary condition [23], is equal to the tangential component in conductor E τ c , which can be determined from the differential Ohm's law, as follows: where J is the current density in the equivalent conductor of the catenary system, and σ is the specific electrical conductivity of the material of the contact wire. Given that the cross-section area of an equivalent contact wire is A = 193.3 × 10 −6 m 2 , with a current of I c = 208.9 A, we find that the current density is J = 1.08 × 10 6 A/m 2 , assuming that the conductivity of copper is σ = 5.8 × 10 7 S/m. Then, the tangential component of tension E τ in the air, according to (11), equals 0.186 V/m. The distance between the contact wire and the head of the rail usually varies from h min = 5.75 m to h max = 6.8 m.
Assuming the maximum distance and the above-mentioned voltage U = 26.12 kV, we find the normal component of the tension, which is equal to E n = 3841 V/m. Consequently, as far as E n E τ , we could find that electric field tension in the airspace of inter-substation zones is E = E n = 3841 V/m.
The tension of the magnetic field at any point of the inter-substation space Each of these components has to be defined using the law of total currents in catenary I c and rails I r (Figure 1), as below: where x is the wire-to-rail coordinate, which varies from an equivalent contact wire to an equivalent rail and is mostly equal to the distance between them h = 6.8 m, according to Figure 1; r c and r r are the radii of equivalent contact wires and rails, which are equal to r c = 7.85 × 10 −3 m and r r = 72.7 × 10 −3 m, respectively. So, we can calculate, according to Formulae (12)- (14), the dependence of the resulting tension of the magnetic field H in the inter-substation space on the x coordinate, i.e., on height h. The results are shown in Figure 2 and Table 2. Let us calculate the possible reactive power Q in a volume V of the space of the inter-substation zone. The volume could be found using the following data: the total length of the space of the inter-substation zone of 560 km, its width is taken to be the width of the railway track of 1.524 m, and its height of h = 6.8 m; then, V = 5, 803, 392 m 3 .
Given that the electric field is uniformly distributed over the inter-substation space and magnetic power is distributed according to Figure 2, the reactive power spreading from TS to ERS is found by specified Formula (10), as follows: Finally, based on Expression (15) and Table 3, we can conclude that an electric locomotive in the specified mode with electrical parameters U = 26.12 kV, I r = 208.9 A, I r = 83.6 A, totally could transfer a reactive power of 2930.5 kvar from the feeder zone.  Given that the electric field is uniformly distributed over the inter-substation space and magnetic power is distributed according to Figure 2, the reactive power spreading from TS to ERS is found by specified Formula (10), as follows: Finally, based on Expression (15) and Table 3, we can conclude that an electric locomotive in the specified mode with electrical parameters = 26.12 kV, = 208.9 A, = 83.6 A, totally could transfer a reactive power of 2930.5 kvar from the feeder zone.

Conditions for the Formation of Standing and Mixed Electromagnetic Waves
The alternating electromagnetic field created in the airspace of inter-substation zones is characterized by the harmonic plane running waves of electric and magnetic field tensions that transmit electricity from TS to ERS [23], as follows: E(y, t) = E m in sin(ωt − βy) + E m r f sin(ωt + βy) (16) H(y, t) = H m in sin(ωt − βy) + H m r f sin(ωt + βy) Let us rewrite Expressions (16) and (17) As it is known, a reflected wave will be formed if the incident wave meets another medium (dielectric or conductor) on its way. For a medium with infinite electrical conductivity σ (a lossless medium), when |Γ| = 1, the amplitude of the reflected wave is equal to the amplitude of the incident wave and, as a result of their adding together, a standing wave is formed. In this regard, given that E m = E m in = E m r f , H m = H m in = H m r f , we obtain the following expressions of standing waves: E(y, t) = 2E m cos(βy) sin(ωt) (20) H(y, t) = −2H m sin(βy) cos(ωt) A graphical representation of the standing waves of electric and magnetic fields reflected from the railway vehicle is shown in Figure 3.  Consequently, it is generally accepted that a necessary condition for the oc of standing electromagnetic waves is a complete reflection of the incident wave analyze this issue, considering [23,25,[37][38][39], as follows.
As it is widely recognized, if a plane sinusoidally changing the electromagne moving in medium 1 (in our case, it is the air), with properties , , , encou its way the interface of another medium 2 (parts and devices of the railway vehic properties , , , and this wave falls (the incident wave having tensions perpendicularly to the boundary of the media, it is partially reflected ⃗ , ⃗ w flection ratio Γ. Consequently, it is generally accepted that a necessary condition for the occurrence of standing electromagnetic waves is a complete reflection of the incident waves. Let us analyze this issue, considering [23,25,[37][38][39], as follows.
As it is widely recognized, if a plane sinusoidally changing the electromagnetic wave moving in medium 1 (in our case, it is the air), with properties ε 1 , σ 1 , µ 1 , encounters on its way the interface of another medium 2 (parts and devices of the railway vehicle), with properties ε 2 , σ 2 , µ 2 , and this wave falls (the incident wave having tensions perpendicularly to the boundary of the media, it is partially reflected → E r f , → H r f with a reflection ratio Γ.
where T = η 1 /η 2 -running wave coefficient; η 1 and η 2 -wave complex resistances of media 1 and 2, defined as follows: -for the conducting medium η = jωµ/σ Let us analyze the flow of electricity from the traction substation to the devices of the inter-substation zone.

Energy Flow to the Contact Wire
At the boundary of the contact wire, the electromagnetic wave falling from the air of the inter-substation zone is partly reflected from the wire surface and partly penetrates into it ( Figure 1). Let us determine the properties of the reflected wave.
For air, medium 1, the wave resistance is η air = η 1 = 376.7 Ω. Medium 2 is a conductive medium represented by the catenary system of the PBSM-95 + MF − 100 type. The cross-section area of the equivalent contact wire is 193.3 × 10 −6 m 2 , hence the specific conductivity is σ 2 = σ eq w = 5 × 10 7 S/m and the permeability is µ r2 = 10.
The higher harmonics of the feeder current in the normal operating mode have frequencies up to 750-1950 Hz, which will be taken into account too. The wave resistance of the material of the equivalent wire for f = 1950 Hz is 1.33 Ω according to (24), at the same time for f = 50 Hz it is equal to 0.212 Ω. Then, the reflection ratio on the boundary "air-to-wire" surface, for example, for the current harmonic with a frequency of f = 1950 Hz is defined as Γ 12 ≈ −1 according to Expression (22).
Then, according to Expressions (23)-(24), the field tensions of the reflected wave are as follows: that is, the almost completely reflected wave changes the sign of vector → E, and a similar result is obtained for the rails.

Incidence of the Poynting Vector on the ERS Roof
As defined in [33], the electromagnetic energy supplying the ERS comes from the traction network of the power supply system in two ways (Figure 4). The first way is contact wire 2-roof bushing insulator 4-appliances of the high-voltage cabinet 6-traction motors. The second way is rails 3-bottom part of the ERS-traction motors and the highvoltage cabinet 6. A significantly minor amount of the energy of the electromagnetic waves comes from the airspace 1 of the inter-substation zones to any part of the EPS, as follows: metal body 8 and to the glass of the driver's cabin 7 in (Figure 4). traction network of the power supply system in two ways (Figure 4). The first way is contact wire 2-roof bushing insulator 4-appliances of the high-voltage cabinet 6-traction motors. The second way is rails 3-bottom part of the ERS-traction motors and the highvoltage cabinet 6. A significantly minor amount of the energy of the electromagnetic waves comes from the airspace 1 of the inter-substation zones to any part of the EPS, as follows: metal body 8 and to the glass of the driver's cabin 7 in (Figure 4).  At the outset, let us analyze the basic stream in which the following components of the Poynting vector carrying this stream can be identified (Figure 4):

1.
Basic vector 0.5 → S basic corresponding to the density of energy flow coming into the elements of the traction power circuit of the ERS (i.e., into the high-voltage cabinet) through the roof bushing porcelain insulator 4. To analyze the electromagnetic waves carrying → S basic , the coefficients of reflection Γ 14 and refraction T 14 (i.e., transmission, penetration) of the incident wave on the "air 1-to-porcelain 4" border have to be defined. For perfect dielectrics like air ε r1 = 1.0 and porcelain ε r4 = 6.0, so for Γ and T coefficients we will get the following: As we can see, some part of the energy is reflected from the insulator, as follows: where r f index denotes the reflected value, in index denotes the incident value, pn index denotes the refracted (penetrated, transmitted) value.
The second stream of → S basic (Figure 4) penetrates into the high-voltage cabinet of the ERS from the space around rails 3 (i.e., from the bottom of ERS). The wave propagation and power transmission in this space need more in-depth study, which is not provided in this article.

Vector
→ S roo f is normally pointed to the roof surface ( Figure 4). Its material is a sheet structural steel of 2 mm thick, with the properties of σ = 7 × 10 6 S/m and µ r = 1000. The electromagnetic harmonic waves of various frequencies, incident from air 1 on the border with the surface of metal roof 5 (Figure 4), are partially reflected and partially penetrate into it. They are gradually damped, resulting in electric power losses in the metal roof.
The wave resistance of the roof metal η 5 , found by Expression (24) in the case of the current harmonic with f = 1950 Hz, is 1.33 Ω, and for air it is η 1 = 376.7 Ω. Substituting for η 1 and η 5 in (22), the coefficient of wave reflection from the roof surface being Γ 15 ≈ −1, the refraction (transmission, penetration) coefficient equals T 15 ≈ 0.87 × 10 −5 . A 1950 Hz wave penetrates to the depth of d 5 = 0.115 mm. In the case of basic current harmonic with a frequency of 50 Hz, the wave resistance of the roof metal η 5 is 0.212 Ω, hence the wave penetrates 0.85 mm deep. The metal sheet of the roof is "impenetrable" for the waves, since they are completely damped out within it without reaching an opposite surface. We obtain similar results for the side surfaces of the ERS body.

Incidence of the Waves on the Frontal Part of the Locomotive Body
Electromagnetic waves spread in the middle part of the air space between the contact wire and the rail, and fall on the frontal part of the locomotive body consisting of glass and metallic parts.
The windshield has a thickness of 15 mm and permittivity ε r = 5.5 . . . 10 (let us assume 9.0). At the "air-to-glass" boundary, the waves are partially reflected with a ratio of Γ = −0.5, according to the following: The metal part of the lower frontal part of the ERS body is made of structural steel sheets of 7 mm thick, with parameters σ = 7 × 10 6 S/m and µ r = 1000. Then, the wave resistance of this metal body, according to (25), is equal to η f ront = η 2 = 33.6 × 10 −5 Ω. Substituting η air = η 1 and η f ront = η 2 in (22), the reflection ratio of the waves from the metal part of the body equals Γ 1 = −1.
The performed calculations of the reflection coefficient of electromagnetic waves from electric traction devices (i.e., contact wire, glass, porcelain bushing, as well as metal roofs and sidewalls of the body) show that incomplete reflection of waves from the surfaces of specified parts and devices occurs in the process of electric locomotive operation. Simultaneously, in the space of the inter-substation zone, both active and reactive powers are transferred from TS to ERS. Consequently, there are running and standing waves that superimpose in this space and form mixed electromagnetic waves. Let us show that the necessary condition for standing waves occurrence, namely, the existence of total reflection, is not so strict. The standing waves in the inter-substation space, which characterize the reactive energy [37], are also possible with incomplete reflection, i.e., when the reflection ratio is 0 ≤ Γ < 1.
Given that E m in /H m in = η air = η 1 , where η air is a wave resistance of the airspace of inter-substation zones, let us rewrite Expressions (18) and (19) in the following way: Then, according to the theory of the well-established symbolic method for the analysis of electric circuits, the instantaneous complex values of electric and magnetic field strengths can be found after formal multiplication of Equations (28) and (29) by operator e jωt , as follows: E(y) = E m in 1 + Γe j2βy e j(ωt−βy) Compared to Expressions (18) and (19), where the fields are represented as the sum of the incident and reflected waves, the form of writing of Expressions (30) and (31) corresponds to only one incident wave with an amplitude varying from the coordinate, as follows: These functions are periodic in the inter-substation space and acquire maximum values with (1 + Γ), which are much larger than the amplitude of the incident wave at points called "wave antinodes". Between them, there are "wave nodes", where the amplitude is minimal and equals (1 − Γ) of the incident wave amplitude. The spatial position of these points does not change over time; therefore, for the expressions described by (30) and (31), the definition of "standing wave" is used [40]. Consequently, in a non-perfect conductor (such as the body of an electric locomotive), the reflected wave has a smaller amplitude than the incident wave. So, a small running wave perpendicular to the conductor surface is superimposed on the standing wave, which compensates for power losses in the conductor.
For more advanced analysis, let us consider the dependence of electromagnetic waves type on the coefficient of running wave T. For this, we rewrite the system of Expressions (28) and (29), considering the well-established Euler's formula and Expression (22), then we obtain the following: After some simple transformations of Expressions (34) and (35), and a transition into instantaneous quantities, we obtain the following system: It follows from Expressions (36) and (37) that the expressions for the electric and magnetic field strengths are the sums of running and standing waves. It follows from them that if the coefficient of wave transmission is T = 1, i.e., if the air space of the intersubstation zones is consistent with the materials of the contact wire and the body of the electric rolling stock (η 1 = η 2 ), then standing waves do not arise in the inter-substation space, since there is no wave reflection. If T = 1, i.e., space and materials are not consistent (η 1 = η 2 ), there is an incomplete reflection and both running and standing waves arise in the inter-substation air space, as a result of which a mixed wave is formed. Incomplete reflection is also indicated by the fact that the metals of both the contact wire and the locomotive body are "real" materials, i.e., they have a finite (and not infinite) conductivity. Therefore, it follows that not all the energy of the incident wave is reflected, and some part of it passes into the metal and dissipates in it as heat.
The above suggests that active energy is transmitted from traction substations to locomotives in the inter-substation airspace by running waves. Thus, the energy exchanged between electric and magnetic fields characterizing the standing waves involves the exchange of reactive power too.

Reactive Power Compensation by "Suppression" of Standing Electromagnetic Waves and Discussion of Results
Reactive power consumed by the ERS of different types is quite significant, as in some modes it reaches up to 40-65% of the consumed power [41]. Therefore, the problem of its compensation is significant for electric traction systems. As it follows from the above analysis, this problem has to be solved by "suppression" of standing waves in the space of inter-substation zones. Since standing waves result from reflected waves, reflections of incident waves need to be reduced to compensate for reactive power. Therefore, it is necessary to investigate issues of the suppression of standing waves and how to damp or eliminate reflected waves. As it follows from Expression (22), it is necessary to ensure that wave resistances along both sides of their separation boundary are as close as possible to each other. There are several ways to reduce the reflection coefficient of monochromatic electromagnetic waves [40,42].
The easiest way to reduce reflections is a thin-film method [40]. For the problem solved in this article, this means that a layer of electrically conductive or dielectric coating of certain thickness d with wave resistance η 2 should be applied to the contact wire, rails and the body of the electric rolling stock. In this case, we will get a circuit model consisting of three environments (1,2,3), shown in Figure 5, where η 1 , η 2 and η 3 are wave resistances of the environments, and K 1 , K 2 and K 3 are their wave numbers defined as Therefore, it follows that not all the energy of the incident wave is reflected, and some part of it passes into the metal and dissipates in it as heat.
The above suggests that active energy is transmitted from traction substations to locomotives in the inter-substation airspace by running waves. Thus, the energy exchanged between electric and magnetic fields characterizing the standing waves involves the exchange of reactive power too.

Reactive Power Compensation by "Suppression" of Standing Electromagnetic Waves and Discussion of Results
Reactive power consumed by the ERS of different types is quite significant, as in some modes it reaches up to 40-65% of the consumed power [41]. Therefore, the problem of its compensation is significant for electric traction systems. As it follows from the above analysis, this problem has to be solved by "suppression" of standing waves in the space of inter-substation zones. Since standing waves result from reflected waves, reflections of incident waves need to be reduced to compensate for reactive power. Therefore, it is necessary to investigate issues of the suppression of standing waves and how to damp or eliminate reflected waves. As it follows from Expression (22), it is necessary to ensure that wave resistances along both sides of their separation boundary are as close as possible to each other. There are several ways to reduce the reflection coefficient of monochromatic electromagnetic waves [40,42].
The easiest way to reduce reflections is a thin-film method [40]. For the problem solved in this article, this means that a layer of electrically conductive or dielectric coating of certain thickness with wave resistance should be applied to the contact wire, rails and the body of the electric rolling stock. In this case, we will get a circuit model consisting of three environments (1,2,3), shown in Figure 5, where , and are wave resistances of the environments, and , and are their wave numbers defined as = 2 ⁄ (where = 1,2,3). Air and the contact wire material are considered to be semi-infinite environments, since they are "impermeable" to electromagnetic waves. The incident waves and both meet partial reflections and , respectively, from the first and the second boundaries; the wave passes into the third environment. Air and the contact wire material are considered to be semi-infinite environments, since they are "impermeable" to electromagnetic waves. The incident waves P in1 and P in2 both meet partial reflections P r f 1 and P r f 2 , respectively, from the first and the second boundaries; the wave P in3 passes into the third environment.
The influence of the coating layer (environment 2 in Figure 5) on the formation of reflected waves is analyzed using analogies of the expressions obtained in the theory of electric circuits with distributed parameters for long lines. These analogies allow us to obtain the reflection coefficient from the first boundary of the environments, i.e., in the air of the inter-substation zone, in the form as below [40]: Let us consider under what conditions the reflection coefficient Γ in Expression (39) disappears, i.e., when there are no reflected and, therefore, standing waves.
According to the first variant, if η 1 = η 2 = η 3 , then the real and imaginary parts of the numerator in Expression (39) are equal to zero, and hence Γ = 0. However, such conditions are practically impossible if environment 3 is a contact wire or metal parts of the locomotive body, since η 1 is the wave resistance of the air, η 1 =376.7 Ω, and η 3 is the resistance of the metal, η 3 = 1.33 Ω; as a result, η 1 > η 3 . At the same time, this condition is feasible for dielectric environment 3, that is, the passage porcelain insulator of the locomotive roof or the glass of its frontal part. However, first, it is necessary to choose the dielectric permittivity of porcelain and glass, equal to the dielectric permittivity of the air. This will provide the following condition: η 1 = η 3 . Second, the thickness of environment 2 (coating) must be, regardless of its η 2 , equal to a multiplication of half the wavelength in this environment (d = n·λ 2 /2; n = 1, 2, 3). In this case, the reflection coefficient of the waves from the first boundary in Expression (39) is equal to tan(K 2 ·d) = 0. So, for the first boundary, the following applies: η 2 + η 1 and for the second boundary, the following: Since η 1 = η 3 , the coefficients Γ 1 and Γ 2 have equal amplitudes, but opposite signs, and therefore opposite phases in the reflection planes. However, the passage of the wave reflected from the second boundary, the distance d "back and forth", causes an additional incursion of its phase by the value of 2K 2 d = 2πn. As a result, the waves in the air environment 1, as equal in amplitude, but opposite in phase, cancel out each other and the resulting wave will be absent.
Another way to reduce reflection is based on using the fact that the wave resistance of a non-conductive material is defined as η = µ/ε. By selecting µ and ε, we can get η equal to the air resistance. If the hysteresis loops µ(H) and ε(E) are the same, so that for any pair of tensions H and E ratio µ/ε is the same, then the layer of this absorbing material will be an empty space for the incident wave with a normal incidence.
It is somewhat more complicated to provide a complete or partial reflection of waves from metal parts of the roof and the locomotive body. The most effective way is to make the body with intermediate media with graded or multi-layered (from three to four layers) coatings made of different materials. They could be based on thin films made from metal, magnetodielectric, semiconductor and superplastic materials [40]. In this case, the required value of the reflection coefficient is obtained by a change in thickness of the layers and properties ε, µ, σ of the material. It is necessary to create such an intermediate coating, either graded or multi-layered, in which wave resistance would change exponentially as η 2 = η 1 e βx [40]. It is always possible to choose such value β as to match the impedances of this layer and the air of inter-substation zones. The advantages of multi-layered coatings over single-layered coatings are clear from researches presented in [42,43]. It follows that double-layered coatings have significantly lower reflection Γ (on 12-16%), transmission T (on 37-42%) and absorption A (on 49-58%) coefficients than single-layered coatings.
It should also be noted that the presented-above calculations show results for basic and higher harmonics created by power converters feeding the traction motors. The higher harmonics can be significantly reduced by appropriate configuration of converters and the application of harmonic filters, which exclude waves created by higher frequencies mentioned above in the article. Such a solution substantially improves power quality and reduces reactive energy transferred by higher harmonics. At the same time, the basic component stills transferer reactive power. So, the proposed multi-layered and composite materials could be applied to compensate this power by the way presented above.
These materials have beneficial reflective properties [44][45][46]. It is also necessary to use the fact that a rough (reflexed, corrugated) metal surface reduces the reflection coefficient [43]. The analogical result could be achieved if the same is done with microlayered coatings. Nowadays, layered materials are the most advanced and comparatively cheap solution; their design technique is shown in [47]. We have to note that a more elaborated manufacture of the body of the electric rolling stock, by applying multi-layered reflective coatings, will pay back quickly and noticeably. This conclusion is evident due to the large amounts of reactive power exchanged in inter-substation zones and the everincreasing cost of electricity spent on trains' electric traction.
In conclusion, it must be added that the presented idea in this article must be continued and added with the more comprehensive analysis of the described-below aspects. First, the article considered the field created by the return current flowing in rails because of the high values of magnetic fields. This current concentrates in rails, but its value equals only 40% of the current taken by locomotive. The other part, known as the stray current, enters to the ground and flows through the metal parts and pipelines. Because of the complexity of the process and difficulties to calculate the values in the nearest underground metal infrastructure [48][49][50][51], this current is not considered in the article. Second, for simplification, the influence of eddy currents in the metal body of the electric locomotive is not considered. In our opinion, this issue is very complex and requires separate special studies. Third, as it was written above, the method is in the conceptual phase and its model must be validated through field tests. To establish the adequacy of the "field" approach for reactive power estimation, one should take the following steps. On the first stage: 1.
Select an extended electrified railway line with a few sections, with a length longer than 160 km; 2.
Choose an AC electric locomotive, which will move along this section; 3.
Perform hourly monitoring of field strengths at the surface of the electric loc motive's roof using broadband electric and magnetic field strength meter. Simultaneously, record the voltage across the current collector and locomotive current by oscilloscope and define its harmonic spectrum; 4.
Calculate the volume of reactive energy Q 1 according to Formula (15).
At the second stage of the experimental tests, it is necessary to cover the surface of the roof of the electric locomotive with a single-layer or multi-layer coating film and perform the exact step-by-step procedure for determining the reactive power Q 2 . Finally, after comparing obtained values of the reactive power, if Q 2 > Q 1 , we could conclude that the "field" approach is adequate.

1.
The currently existing methods for analyzing electric power processes in electric traction systems do not consider the fundamental law of the propagation of electromagnetic energy. It does not account accurately the energy spent on train traction; 2.
Electromagnetic harmonic waves spreading in the space of inter-substation zones from the traction substation to electric rolling stock, incident on the surface of the contact wire, as well as on the frontal, lateral and roof surfaces of the body of electric rolling stock, partially reflected and summing up with the incident waves, form standing waves whose energy is reactive; 3.
Reactive power in the space of inter-substation zones is a measure of the asymmetry of the rates of change in the electric and magnetic components of the energy of the electromagnetic field existing in that zone. This power is numerically equal to the difference between the average values of the energies stored in the magnetic and electric fields in the volume of the feeder zone, multiplied by 2ω; 4.
To compensate for the reactive power existing in the space of inter-substation zones, it is necessary to "suppress" (damp) standing waves by applying single-layered or multi-layered film coatings with certain electromagnetic properties to particular parts of the surface of the body of the electric rolling stock; 5.
The stated and theoretically substantiated "field" (based on the theory of the electromagnetic field) approach to the processes of formation, transmission, and "damping" of standing electromagnetic waves, as characteristics of reactive power in intersubstation zones, is applicable to any type of electric locomotives and AC electric traction systems. However, there are no theoretical limitations. Funding: This article was published with financial support under the Open Science Program which was launched at the Warsaw University of Technology as part of the project "Initiative of Excellence-Research University" (IDUB). It was also supported with internal funds of the Electrical Power Engineering Institute.