Improved Methodology for Power Transformer Loss Evaluation: Algorithm Reﬁnement and Resonance Risk Analysis

: The decline in power quality within electrical networks is adversely impacting the energy efﬁciency and safety of transmission elements. The growing prevalence of power electronics has elevated harmonic levels in the grid to an extent where their signiﬁcance cannot be overlooked. Ad-ditionally, the increasing integration of renewable energy sources introduces heightened ﬂuctuations, rendering the prediction and simulation of working modes more challenging. This paper presents an improved algorithm for calculating power transformer losses attributed to harmonics, with a comprehensive validation against simulation results obtained from the Power Factory application and real-world measurements. The advantages of the algorithm are that all evaluations are performed in real-time based on single-point measurements, and the algorithm was easy to implement in a Programmable Logic Controller (PLC). This allows us to receive the exchange of information to energy monitoring systems (EMSs) or with Power factor Correction Units (PFCUs) and control it. To facilitate a more intuitive understanding and visualization of potential hazardous scenarios related to resonance, an extra Dijkstra algorithm was implemented. This augmentation enables the identiﬁcation of conditions, wherein certain branches exhibit lower resistance than the grid connection point, indicating a heightened risk of resonance and the presence of highly distorted currents. Recognizing that monitoring alone does not inherently contribute to increased energy efﬁciency, the algorithm was further expanded to assess transformer losses across a spectrum of Power Factory Correction Units power levels. Additionally, a command from a PLC to a PFCU can now be initiated to change the capacitance level and near-resonance working mode. These advancements collectively contribute to a more robust and versatile methodology for evaluating power transformer losses, offering enhanced accuracy and the ability to visualize potentially critical resonance scenarios.


Introduction
The global transition to renewable energy sources has witnessed a remarkable rise in recent decades, driven by mounting environmental concerns and the pursuit of sustainable energy solutions [1].Solar photovoltaic systems, wind turbines, and other power electronic technologies have emerged as cornerstones of the green energy revolution.In 2022, wind and solar power generated 22% of the EU's electricity, overtaking gas for the first time [2].Additionally, in modern industrial power supply systems, there has been a significant increase in the installed capacity of nonlinear electric sources [3], and total capacity in Europe, only in photovoltaics, has increased by 25% from 167.5 GW in 2021 to 208.9 GW 2022 [4].However, as these sources are integrated into existing electrical grids, complex challenges for power quality appear [5].
The rise of renewable energy sources, often characterized by variable [6] and nonlinear power generation profiles, has brought about a significant transformation in power Energies 2023, 16, 7837 2 of 16 grids [7].Traditional electrical networks, designed for steady-state and sinusoidal operation, are now contending with a new set of operational dynamics [8].These dynamics encompass fluctuations in energy production, voltage irregularities, and a growing prevalence of harmonic distortions [9].Recently, the widespread use of power electronic systems increased in magnitude [10], therefore it has become a key issue in installations [11].Harmonic distortion not only poses threats to the operational lifespan of equipment such as power transformers but also leads to energy losses that impact economic viability [12].
This article delves into the intricate interplay between harmonic distortion and energy losses, with a particular focus on the consequences of harmonic resonance and decreased efficiency of power transformers in power distribution systems.Energy losses are a fundamental issue in the energy sector [13].It is imperative to increase energy efficiency constantly [14].This pursuit is not only beneficial for both industry and the climate [15], but it is also indispensable, as adept energy management can help to reduce costs and increase profit [16].Real-time monitoring and recording of energy consumption data can aid in identifying important energy losses in a facility [17] and allow the implementation of various energy-saving practices to reach up to 8% in energy efficiency [18].Nevertheless, energy monitoring itself does not save energy [19] because data must be not only collected, but also conclusions should be made on how to improve, and this process is significantly more challenging in the presence of harmonics [20].
In the domain of transformer losses attributed to harmonic distortion, the majority of analyses are in alignment with the recommendations outlined in IEEE Std C57.110-2018 recommendations [21][22][23].Building upon similar technological foundations, article [24] uses a comparable methodology to evaluate the maximum loading capacity of transformers in the presence of harmonic influences.Simulations hold a significant role in scientific literature and are indispensable.In the context of transformer losses, it is essential to consider both simulations and real measurements for a comprehensive understanding.Comparisons of transformer losses have been conducted through real measurements in relation to the EMTP-RV model [25,26] and with the Matlab/Simulink model [27].Additionally, simulations using ETAP-Etrax Software have been employed for a big railway system in the article [28].Even though the finite element method (FEM) is an effective technique [29,30] for handling the distribution of losses in various metallic structures [31], it takes a lot of time and processing capacity [32].
Currently, studies on power transformer losses due to harmonics are mainly with simulated data [33] or historical data [34].The disadvantage of those methods are calculations and evaluations, which are performed in one operating point [35].The drawback of these methods lies in the fact that calculations and evaluations are conducted at a single operating point, whereas operating points and the loading profile are subject to constant change.The increasing prevalence of renewables further amplifies their substantial influence in this regard [36].The static load model is suitable for the steady state, as it does not involve load dynamics [37], and this gives errors due to discrete time, or the wide usage is limited.
The most precise assessment of power transformer losses can be achieved through measurements [38], particularly in controlled laboratory conditions [39].However, translating this methodology into the real industrial environment proves to be more costly due to the need for duplicate measurement devices and poses complexities, primarily attributed to the utilization of medium voltage measurement equipment.The goals of lowering energy waste and increasing awareness of energy efficiency can only be achieved by measuring and controlling energy consumption [40].Internet of Things technologies are currently supporting these efforts [41,42] by making a vast amount of data readily available for analysis [43].In article [44], energy efficiency was improved only by gathering data with a monitoring system and then identifying energy losses.To add in different articles [45], it was presented that giving enterprises a tool to assess and comprehend their actual energy management capabilities is vital to find areas where improvements can be made.
Monitoring alone is insufficient, and it is necessary to detect abnormal energy consumption data in real-time [46].In article [47], real-time monitoring was presented ac-Energies 2023, 16, 7837 3 of 16 cording to IEEE Std C57.110-2018 [48], which evaluates the increase in the temperature of the transformer due to harmonic distortion.Publication [49] presented monitoring of insulation ageing, and neural networks were used for evaluation in publication [50].To approach assessment of losses in sufficient manner, it is important to develop universal and efficient algorithmic methods [51] in terms of implementation convenience as well as acceptable computational complexity [52].The significance of this study lies in its potential to improve energy monitoring [53] systems with a real-time algorithm to mitigate risks promptly.Real-time experiments validate that the proposed method is effective and has good application prospects [46].For this reason, article [54] presented a simplified algorithm implementation in the programmable logic controller, which resulted in sufficient calculation precision compared to the measurement.During the experiment, which was presented in article [54], it was noticed that the power factor correction unit (PFCU) decreased efficiency due to ongoing series resonance compared to the state where the PFCU was disconnected.This condition is serious not only because of the derating of the transformer [55], but it also can lead to explosions of the capacitors [56].The algorithm must not only calculate the transformer's losses due to harmonics but also evaluate harmful conditions such as resonance.
This paper is organized as follows: Section 2 describes of the procedure of calculating resonance frequencies, additional transformer losses due to harmonic presence, and the equation explanations with the schemes for evaluating resistances versus frequency.The results of the calculations and comparison to simulations are provided in Section 3. The paper ends with the conclusion section, where the results and future work are presented.

Materials and Methods
One of the vital issues in industrial and distribution power systems is increasing the use of nonlinear loads and increasing power of inverter-based distributed generation, which leads to higher harmonic disturbances in a grid [57].Harmonic currents in the inverter are generated because of usage pulse-width modulation [58].In power systems, harmonics are generally managed by adding filters, circuit detuning, and usage of phase-shifting transformers [59] or inverters with harmonic current compensation compatibility [60].
The power transformer, the main equipment in the utility network [61], plays a very important role because the failure of the transformer can cause serious power outages [20].Under normal operating conditions, transformers are typically not a significant source of harmonics in power systems.Despite historically being the first cause of harmonics, their relationship between primary voltage and currents follows a non-linear magnetization curve, experiencing distortion if located within saturation regions [62].This distortion affects the magnetizing current, but the impact on harmonics is minimal even with a large number of transformers operating [63].In addition to the fact that the transformer is a harmonic source itself, the transformer's efficiency can drop drastically while supplying non-linear loads.
Normally in industrial plants in parallel to transformers, PFCUs are installed, which consist of capacitor banks.The capacitors are not a source of harmonic current but can have the effect of magnifying such currents [64].Capacitor banks in industrial power systems are mainly used for power factor correction purposes.A PFCU with installed capacitor banks is predominantly employed in reactive power compensation technology owing to its highly cost-effective pricing.In various applications of a PFCU within a power system network, the interaction between the inductance of a transformer and the capacitance of a PFCU introduces a noteworthy consideration.Specifically, due to the combinations of these elements, a frequency is invariably present at which the capacitance aligns in either parallel or series resonance with the transformer or the grid.When this condition appears on, or close to, cases of series resonance, the increased harmonic currents will circulate through the power transformer between the utility network and the power factor equipment, causing additional losses to the system [65].Harmonic resonance is a key concern in the power quality analysis of industrial power systems.This event can threaten the reliability of the system and can cause significant issues in industrial plants [66].Under series resonance conditions, there is an augmentation of current through the capacitor bank, leading to supplementary losses.Consequently, this can give rise to adverse effects, including circuit breaker tripping or fuses blowing.Additionally, physicochemical processes occurring in the dielectric may result in accelerated ageing of the insulation, contributing to heightened service costs, permanent damage, and related complications [67].
If the system is considered from the point of view of the distribution network, then the power transformer is in series resonance with the power factor correction unit [68].The resonant frequency f res can be found with Equation (1).
where Q C is the reactive power of the capacitor, S r is the apparent rated power of the transformer, and U k is the impedance of the voltage transformer in percent.
In article [54], its algorithm to evaluate the effect of harmonics on a power transformer was implemented in PLC.This implementation allowed to measure losses and efficiency of the power transformer in real-time and was tested in real conditions, and the results were compared to the measurements of the MV side of the transformer.Observations during the experiment revealed that ongoing resonance has a more pronounced impact on the efficiency of the transformer compared to power factor correction.In certain scenarios, it was determined that it might be more advantageous not to compensate reactive power than to engage in compensation while neglecting the associated risk of resonance.To create an algorithm, firstly, key parameters of the transformer must be calculated.The copper losses of the power transformer P I 2 R in kW are calculated according to Equation (2).
where K is the constant coefficient (for a single-phase transformer, K equals 1; for a threephase transformer, K equals 1.5), R 1 and R 2 are the resistances of primary and secondary windings, respectively, and I 1f and I 2f are the rated current of primary and secondary windings, respectively.When the copper losses of the power transformer are known, the nominal stray losses P TSL-N are found according to Equation (3).
where P LL is the rated load losses of the transformer.According to article [48], the ratio of total stray losses for the oil-filled power transformer, equal to 1000 kVA of power, is 50 percent for the nominal Eddy current losses P EC-N in kW and 50 percent for other stray losses P OTSL in kW; these can be calculated according to Equations ( 4) and (5).
For distorted currents, RMS of current I pu in proportional units is calculated according to Equation (6).
Energies 2023, 16, 7837 5 of 16 where h is the harmonic order, h max is the highest order of measured harmonic, and I h is the per unit rms current at harmonic h.The harmonic loss factor of Eddy current losses F HL-STR is calculated according to Equation (7).
The harmonic loss factor of other stray losses is F HL-STR calculated according to Equation (8).
The loading of the transformer P LLpu per unit is calculated according to Equation (9).
where S L is the apparent power of the load in kVA, and S N is the apparent power of the power transformer in kVA.The losses of the power transformer P TL in W with given harmonic content can be calculated according to Equation (10).(10) where P 0 is the idle losses of the transformer.The active power load P MV in kW on the MV side is calculated according to Equation (11).
P MV = P TL + P L (11) where P L represents the active load of the transformer.
To implement this algorithm, the PLC is needed to do frequency sweep calculations in real-time.To calculate the parameters for the resonance, we need to calculate the reactive resistances of the power factor correction unit and the power transformer.Considering what is presented in Figure 1, a simple system of an industrial distribution network, which has a transformer and a PFCU, the MV busbars are connected to the grid via an HV transformer.
where h is the harmonic order, hmax is the highest order of measured harmonic, and Ih is the per unit rms current at harmonic h.The harmonic loss factor of Eddy current losses FHL-STR is calculated according to Equation (7).
The harmonic loss factor of other stray losses is FHL-STR calculated according to Equation (8).
The loading of the transformer PLLpu per unit is calculated according to Equation (9).
where SL is the apparent power of the load in kVA, and SN is the apparent power of the power transformer in kVA.The losses of the power transformer PTL in W with given harmonic content can be calculated according to Equation (10).
where P0 is the idle losses of the transformer.The active power load PMV in kW on the MV side is calculated according to Equation (11).
where PL represents the active load of the transformer.
To implement this algorithm, the PLC is needed to do frequency sweep calculations in real-time.To calculate the parameters for the resonance, we need to calculate the reactive resistances of the power factor correction unit and the power transformer.Considering what is presented in Figure 1, a simple system of an industrial distribution network, which has a transformer and a PFCU, the MV busbars are connected to the grid via an HV transformer.If active resistance is neglected, the equivalent impedance Z z can be calculated according to Equation (12).
where L s is the inductance of the system, L ST is the inductance of the HV transformer, L TR is the inductance of the MV transformer, f is the frequency, and C PFCU is the capacitance of the power factor correction unit.
To calculate different scenarios of the PLC, there is a need to know how current injections in industrial branches are changed due to different capacitances of the PFCU.First Kirchhoff's law states that the current flowing into a node (or a junction) must be equal to the current flowing out of it.In Figure 2, it is presented that the current injected into the distribution busbar is divided into the grid and to another industrial branch.
Energies 2023, 16, x FOR PEER REVIEW 6 of 16 If active resistance is neglected, the equivalent impedance Zz can be calculated according to Equation (12).
where Ls is the inductance of the system, LST is the inductance of the HV transformer, LTR is the inductance of the MV transformer, f is the frequency, and CPFCU is the capacitance of the power factor correction unit.
To calculate different scenarios of the PLC, there is a need to know how current injections in industrial branches are changed due to different capacitances of the PFCU.First Kirchhoff's law states that the current flowing into a node (or a junction) must be equal to the current flowing out of it.In Figure 2, it is presented that the current injected into the distribution busbar is divided into the grid and to another industrial branch.When the resistance values of the system are known, the ratio of current absorbed by the industrial branch I1% and grid can be calculated according to Equation (13) if the total injected current IT is known.
where RS is the resistance of the network, RST resistance of the HV transformer, and RTR is the resistance of the industrial branch transformer.One of the benefits of algorithm was listed as the possibility to use it only with the measurements on the low voltage side, which is cheaper and safer.For this reason, Equation ( 14) can be used to calculate the harmonic current flown in MV busbars IT in case the measurements are only available on the 0.4 kV side (IB).
Equations ( 13) and ( 14) can be represented in an equivalent scheme of the network, which is provided in Figure 3.When the resistance values of the system are known, the ratio of current absorbed by the industrial branch I 1% and grid can be calculated according to Equation (13) if the total injected current I T is known.
) 2 I T (13) where R S is the resistance of the network, R ST resistance of the HV transformer, and R TR is the resistance of the industrial branch transformer.
One of the benefits of algorithm was listed as the possibility to use it only with the measurements on the low voltage side, which is cheaper and safer.For this reason, Equation ( 14) can be used to calculate the harmonic current flown in MV busbars I T in case the measurements are only available on the 0.4 kV side (I B ).
Equations ( 13) and ( 14) can be represented in an equivalent scheme of the network, which is provided in Figure 3. Equations ( 13) and ( 14) allow us to recalculate the current drawn from the MV g to the industrial plant branch with different powers of the PFCU and to find the corr tion in which PFCU mode is best for the power transformer loss, considering not o displacement power factor but harmonic distortion too.
Harmonic currents will tend to flow into the system capacitance [69], or, in the of resonance, relatively higher currents can flow to the branch that has the smallest sistance.For quick evaluation of the possible current harmonic path, the Dijkstra a rithm can be used.The algorithm can be used for finding the shortest paths from one tex to another in a weighted graph with non-negative weights.In this scenario, the we of the path is equal to the absolute impedance of the branch.The Dijkstra algorithm added to the algorithm for a more user-friendly understanding of impedance relat ships and harmonic current flow.In the test device, the algorithm is used four time calculate harmonic current flow in the following frequencies: 250 Hz, 350 Hz, 550 Hz, 650 Hz. Figure 4 shows a simplified equivalent scheme for the usage of the Dijkstra a rithm.For harmonic current to flow from point A to point B, there are two possible w with different resistances.This algorithm calculates resistance for desired frequenc and, in the case of resonance, the path through ZTR1 will be shown; this means that transformer is exposed to higher current harmonics.Equations ( 13) and ( 14) allow us to recalculate the current drawn from the MV grid to the industrial plant branch with different powers of the PFCU and to find the correlation in which PFCU mode is best for the power transformer loss, considering not only displacement power factor but harmonic distortion too.
Harmonic currents will tend to flow into the system capacitance [69], or, in the case of resonance, relatively higher currents can flow to the branch that has the smallest resistance.For quick evaluation of the possible current harmonic path, the Dijkstra algorithm can be used.The algorithm can be used for finding the shortest paths from one vertex to another in a weighted graph with non-negative weights.In this scenario, the weight of the path is equal to the absolute impedance of the branch.The Dijkstra algorithm was added to the algorithm for a more user-friendly understanding of impedance relationships and harmonic current flow.In the test device, the algorithm is used four times to calculate harmonic current flow in the following frequencies: 250 Hz, 350 Hz, 550 Hz, and 650 Hz. Figure 4 shows a simplified equivalent scheme for the usage of the Dijkstra algorithm.For harmonic current to flow from point A to point B, there are two possible ways, with different resistances.This algorithm calculates resistance for desired frequencies, and, in the case of resonance, the path through Z TR1 will be shown; this means that the transformer is exposed to higher current harmonics.
The Dijkstra algorithm's implementation allows us to monitor frequency sweep changes in bigger systems and detect dangerous resonance scenarios in real time.The Dijkstra algorithm's implementation allows us to monitor frequency sweep changes in bigger systems and detect dangerous resonance scenarios in real time.

Results
To check the accuracy of implemented the algorithm, the same conditions as those presented in article [54] were modeled in a Power Factory environment.The model consisted of two MV voltage connection points, which were the first branch representing the consumer with the transformer and PFCU and the second branch representing the current distortion source.The model of the Power Factory environment is shown in Figure 5.

Results
To check the accuracy of implemented the algorithm, the same conditions as those presented in article [54] were modeled in a Power Factory environment.The model consisted of two MV voltage connection points, which were the first branch representing the consumer with the transformer and PFCU and the second branch representing the current distortion source.The model of the Power Factory environment is shown in Figure 5.The Dijkstra algorithm's implementation allows us to monitor frequency sweep changes in bigger systems and detect dangerous resonance scenarios in real time.

Results
To check the accuracy of implemented the algorithm, the same conditions as those presented in article [54] were modeled in a Power Factory environment.The model consisted of two MV voltage connection points, which were the first branch representing the consumer with the transformer and PFCU and the second branch representing the current distortion source.The model of the Power Factory environment is shown in Figure 5.The input data for the simulation are listed in Table 1.It was noticed that additional iterations were needed to improve the accuracy of the algorithm because the transformer loading increased due to the harmonics losses.The results after the improved algorithm are listed in Table 3.When the model was tested, it allowed us to do the simulations with more various scenarios, such as different powers in the PFCUs.Figure 6 shows us that, due to resonance even with a relatively small THD of current, the most efficient scenario does not fully compensate for the reactive power.For this reason, it is useful to make calculations in real-time, and PLC allows for the possibility to receive real data from the PFCU; in the case of a resonance scenario, the PLC can send a command to the PFC to change its power.
Energies 2023, 16, x FOR PEER REVIEW 9 of 16 The input data for the simulation are listed in It was noticed that additional iterations were needed to improve the accuracy of the algorithm because the transformer loading increased due to the harmonics losses.The results after the improved algorithm are listed in Table 3.When the model was tested, it allowed us to do the simulations with more various scenarios, such as different powers in the PFCUs.Figure 6 shows us that, due to resonance even with a relatively small THD of current, the most efficient scenario does not fully compensate for the reactive power.For this reason, it is useful to make calculations in realtime, and PLC allows for the possibility to receive real data from the PFCU; in the case of a resonance scenario, the PLC can send a command to the PFC to change its power.To evaluate the resistances and frequency sweep of the system's simplified model, Power Factory was used.A simplified model with a neglected load is provided in Figure 7.To evaluate the resistances and frequency sweep of the system's simplified model, Power Factory was used.A simplified model with a neglected load is provided in Figure 7.In Figure 8, network reactance results according to the frequency sweep in the Power Factory are shown.To evaluate the resistances and frequency sweep of the system's simplified model, Power Factory was used.A simplified model with a neglected load is provided in Figure 7.According to Equation (13), the current ratios between branches are provided in Figure 9.It can be seen the harmonic current to branch is exceeding the harmonic current to grid only in the resonance scenario.The results of Equation ( 14), recalculating the total distorted current in MV busbars, are represented in Table 4.The results of Equation ( 14), recalculating the total distorted current in MV busbars, are represented in Table 4.When the currents circulating in MV busbars are known, it allows us to recalculate the current drawn by the branch in different power of PFCU IB'.The results are provided in Table 5.The table shows that the harmonic distortion rises due to the increasing power of the PFCU; this is because the resonant frequency is shifting towards lower frequencies, according to Equation (1), and resonant frequency can be calculated.The results of the resonant frequency calculations due to different powers of the PFCU are shown in Figure 10.Recalculated current IB' allows us to find the most efficient scenario for the power transformer using Equation (10).The results are shown in Figure 11.Recalculated current IB' allows us to find the most efficient scenario for the power transformer using Equation (10).The results are shown in Figure 11.Recalculated current IB' allows us to find the most efficient scenario for the power transformer using Equation (10).The results are shown in Figure 11.The working regime of the PFCU with 100 kVAr of power leads to decreasing load losses of 48 percent compared to a fully compensated working mode of 140 kVAr.This leads to lower usage of the PFCU due to lower load and lower current harmonic quantities, which are shown in Figure 12.The working regime of the PFCU with 100 kVAr of power leads to decreasing load losses of 48 percent compared to a fully compensated working mode of 140 kVAr.This leads to lower usage of the PFCU due to lower load and lower current harmonic quantities, which are shown in Figure 12.Recalculated current IB' allows us to find the most efficient scenario for the power transformer using Equation (10).The results are shown in Figure 11.The working regime of the PFCU with 100 kVAr of power leads to decreasing load losses of 48 percent compared to a fully compensated working mode of 140 kVAr.This leads to lower usage of the PFCU due to lower load and lower current harmonic quantities, which are shown in Figure 12.

Conclusions
Load-flow calculations, while more accurate, are typically executed using historical income data.However, this article establishes that real-time calculations prove to be more effective than load-flow simulations.This effectiveness is attributed to the constantly changing conditions, including, but not limited to, harmonic distortion, voltage levels, powers of the Power Factor Correction Units (PFCUs), and the dynamic behaviors of consumers.To implement the algorithm, PLC was chosen, due to its ability to collect data from power meters, PFCUs, and the possibility to control the PFCUs.To check the accuracy of the algorithm's results, they were compared with the simulation results of Power Factory applications and real measurements.The calculations resulted in a sufficient difference between the calculated values and modelled values, with a relative error of 0.024 percent.The accuracy of the algorithm was improved due to adding additional iterations, and the relative calculation error to simulated data decreased to 0.0108 percent.An additional Dijkstra algorithm was implemented to easier understand or visualize dangerous scenarios of resonance.The Dijkstra algorithm allows for identifying the conditions when some of the branches have lower resistances than the grid connection point, and this means risk of resonance and highly distorted currents.The algorithm operates in a straightforward manner; it checks all the possible directions for the harmonic flow and chooses the path with lower resistance.In normal conditions, the lower resistance path is to the grid, except in resonance conditions.Using the first Kirchhoff's law, the calculation of current harmonic distributions was tested, assuming knowledge of the injected harmonic current.Subsequently, the calculations were reversed to ascertain the extent of harmonics flowing in medium voltage busbars when measurements were solely conducted on the low voltage side.This allowed us to implement the calculation of which power of the PFCU to use; how much harmonic current would be injected into the industrial branch instead of a grid; and to find the most efficient scenario of the PFCU, taking into account not only displacement power factor but harmonic distortion too.The results presented in the results section show that, even with a relatively small THD of current, the transformer losses are affected.The working regime of the PFCU with 100 kVAr of power leads to decreasing load losses of 48 percent compared to a fully compensated working mode of 140 kVAr.After evaluation of transformer losses with different PFCU powers, controlling the PFCU can be achieved only by using industrial communication protocols.The algorithm is easy to implement, and, in different hardware, the only limitations are due to support of the communication between devices.In many cases, only supporting serial communications as a Modbus RTU can be enough, due to its wide usage in industrial equipment.Communication requirements can also be different due to support protocols of EMSs for logging the measurement data.This algorithm should not overlap with the existing control because it is the main target just to control the target power factor, and, in most of the cases, this setting in the controller of the PFCU is fixed.
Future work will focus on calculations for bigger systems and evaluating the effects of voltage harmonic distortion.

Figure 1 .
Figure 1.A simple industrial system consisting of the industrial branch with a transformer and a PFCU.

Figure 1 .
Figure 1.A simple industrial system consisting of the industrial branch with a transformer and a PFCU.

Figure 2 .
Figure 2. Equivalent scheme to represent harmonic currents flow.

Figure 2 .
Figure 2. Equivalent scheme to represent harmonic currents flow.

Figure 3 .
Figure 3. Equivalent scheme to represent harmonic currents flow values.

Figure 3 .
Figure 3. Equivalent scheme to represent harmonic currents flow values.

Figure 4 .
Figure 4. Equivalent scheme to represent logic for Dijkstra algorithm.

Figure 5 .
Figure 5. Model in the Power Factory environment.

Figure 4 .
Figure 4. Equivalent scheme to represent logic for Dijkstra algorithm.

Figure 4 .
Figure 4. Equivalent scheme to represent logic for Dijkstra algorithm.

Figure 5 .
Figure 5. Model in the Power Factory environment.

Figure 5 .
Figure 5. Model in the Power Factory environment.

Figure 6 .
Figure 6.Changes in the transformer losses due to different power of PFCU.

Figure 6 .
Figure 6.Changes in the transformer losses due to different power of PFCU.

Figure 7 .
Figure 7. Model used to evaluate resistances and harmonic flow.In Figure8, network reactance results according to the frequency sweep in the Power Factory are shown.

Figure 8 .
Figure 8. Results of frequency sweep in the Power Factory application.According to Equation(13), the current ratios between branches are provided in Figure9.It can be seen the harmonic current to branch is exceeding the harmonic current to grid only in the resonance scenario.

Figure 7 .
Figure 7. Model used to evaluate resistances and harmonic flow.

Figure 6 .
Figure 6.Changes in the transformer losses due to different power of PFCU.

Figure 7 .
Figure 7. Model used to evaluate resistances and harmonic flow.In Figure8, network reactance results according to the frequency sweep in the Power Factory are shown.

Figure 8 .Figure 8 .
Figure 8. Results of frequency sweep in the Power Factory application.According to Equation(13), the current ratios between branches are provided in Figure9.It can be seen the harmonic current to branch is exceeding the harmonic current to grid only in the resonance scenario.

Energies 2023 , 16 Figure 9 .
Figure 9. Division of the currents to industrial branch and grid.

Figure 9 .
Figure 9. Division of the currents to industrial branch and grid.

Figure 10 .
Figure 10.Resonant frequency versus different powers of the PFCU.

Figure 10 .
Figure 10.Resonant frequency versus different powers of the PFCU.

Figure 11 .
Figure 11.Calculated transformer losses due to different powers of PFCU.

Figure 12 .Figure 11 .
Figure 12.Changes in the harmonic current to branch due to different powers of PFCU.

Figure 10 .
Figure 10.Resonant frequency versus different powers of the PFCU.

Figure 11 .
Figure 11.Calculated transformer losses due to different powers of PFCU.

Figure 12 .Figure 12 .
Figure 12.Changes in the harmonic current to branch due to different powers of PFCU.

Table 1 .
Input data for the simulation.

Table 2
provides the results of the real measurement, the calculation with the PLC, and the simulation with the Power Factory.

Table 2 .
Comparison between real measurements, calculations, and simulation.

Table 3 .
Comparison after additional iteration added to the algorithm.

Table 1 Table 1 .
Input data for the simulation.

Table 2
provides the results of the real measurement, the calculation with the PLC, and the simulation with the Power Factory.

Table 2 .
Comparison between real measurements, calculations, and simulation.

Table 3 .
Comparison after additional iteration added to the algorithm.

Table 4 .
Calculated total distorted current injected to MV busbars.

Table 4 .
Calculated total distorted current injected to MV busbars.

Table 5 .
Calculated total distorted current injected to industrial branch with different powers of PFCU.