Assessment of Load Losses Caused by Harmonic Currents in Distribution Transformers Using the Transformer Loss Calculator Software

: Transformer load losses cause various adverse effects, such as derating, a decreased lifetime, and greenhouse gas emissions. In this paper, the load losses caused by non-linear loads on distribution transformers are analyzed. For this study, the load loss expressions provided by the IEEE Standard C57.110 and ANSI/UL 1561-1562 were adapted to the usual case where the transformer currents differ in each phase. The novel load loss expressions adapted from the IEEE Standard C57.110 were applied using the software known as the “Transformer Loss Calculator” (TLC), implemented with LabVIEW. For the application of new load loss expressions, carbon dioxide (CO 2 ) emissions were determined by multiplying the load losses by the emission factors of each country. The experimental results are based on the recordings made by a FLUKE 435 Series II analyzer on the second of two 1000 kVA transformers, feeding real residential distribution networks with very differently distorted loads. An analysis of these transformers shows that the annual energy losses and CO 2 emissions obtained from the adapted load loss expressions could be more than 5% of those determined by the original IEEE and ANSI Standard expressions. Due to these percentage loss and emission differences, it is advisable to use the TLC software in transformer monitoring instruments.


Introduction
The growth of the demand for electrical energy in recent decades means that it is necessary to use increased numbers of transformers in electrical networks [1].These static electrical machines are used to supply electrical power by adapting high-voltage (HV) supply voltages to the values required by customers' Low-Voltage (LV) electrical installations.
However, transformers consume power.They suffer load losses [2] due to the circulation of currents through their windings, as well as core losses due to magnetic hysteresis and Eddy current phenomena.The latter losses are little affected by harmonics [3,4].The losses are seen in the heat energy, resulting in the internal temperature of transformers being raised and leading to two important effects: The first is the deterioration of the transformer and a reduction in its lifetime [1,[5][6][7][8][9].The second effect is the impact of the load losses of the transformer on the environment, measured as carbon dioxide (CO 2 ) emissions into the atmosphere [10][11][12].These emissions are proportional to the energy consumed, that is, to the losses in the case of transformers, with the proportionality coefficient being the emission factor, provided by the grid managers of each country or region [13][14][15].
The presence of non-sinusoidal currents caused by non-linear loads further increases the load losses of the transformers compared to linear loads.Currently, non-linear loads have a large presence in power grids.Discharge lamps in street lighting and fluorescent lighting, as well as LED lamps, motor drivers, and welding devices, among other load electronic devices, are examples of non-linear loads.Likewise, the integration of renewable energies in the electrical network, as is the case with the implementation of solar energy, even in the domestic sphere, has increased the circulation of non-sinusoidal currents both in the lines and in the distribution transformers of LV networks.Based on the above, it can be inferred that the losses caused by non-sinusoidal currents in transformers are very important.
Furthermore, the losses caused by non-sinusoidal currents do not evolve in the same way as those produced by sinusoidal currents.The latter depend on the Root Mean Square (RMS) value of currents, while the losses caused by non-sinusoidal currents also depend on the frequency of the harmonics of these currents [16][17][18][19][20][21][22][23].The IEEE Standard C57.110-2018 [24,25] establishes that, in the presence of non-sinusoidal currents, the load losses (P LL ) of the transformers are the sum of three others: Each of these losses is due to a different physical phenomenon.The first is the wellknown Joule effect, which occurs due to the circulation of currents through the ohmic resistances of the primary and secondary windings.The losses caused by this phenomenon (P JL ) are known as copper losses (or ohmic losses), and their values are proportional to the squares of the RMS values of the currents flowing through the windings.Winding Eddy Current Losses (P ECL ) occur due to the skin and proximity phenomena [3] and depend on the squares of the RMS values of the currents and the frequencies of their harmonics.Other Stray Losses (P OSL ) occur due to the Eddy currents induced by stray magnetic fluxes in the tank and other metallic parts of the transformer [3].These losses are also dependent on the frequencies of the harmonics, but to a lesser extent than Winding Eddy Current Losses (P ECL ).In dry-type transformers, as listed in the ANSI/UL Standard 1561-1562 [26,27], Other Stray Losses (P OSL ) are considered to be negligible [19,26,27].This consideration has also been extended to oil-immersed transformers in papers [16,17,22,28] and transformer manufacturers [29].
The expressions of the load losses were originally established by the IEEE Standard C57.110 and ANSI/UL 1561 and 1562 based on the assumption that the above phenomena affect the three phases of the transformer equally.All of the references we have consulted also make this assumption, even in applications where the transformer currents do not have the same RMS values and harmonic content in each phase, as determined by the Total Harmonic Distortion (THDi z %) of each phase (z = A, B, C): where the order of the harmonic h z = f hz / f 1 is the ratio of the frequency of that harmonic ( f hz ) and the fundamental frequency ( f 1 = 50 − 60 Hz), h z,max is the order of the highest frequency harmonic used in the analysis, and I hz is the RMS value of the harmonic current of order h z and phase z.However, the last two phenomena (Eddy currents and Other Stray Losses) are very dependent on the frequencies of the currents; therefore, the expressions of the load losses included in the IEEE and ANSI standards do not provide correct values in the usual case where the transformer currents are different in each phase.
To avoid these limitations, in this article, the expressions of load losses for three-phase transformers were developed, which are generally applicable, even when the currents are different in each phase of the transformer.This is the main novelty of this paper.These equations are established in Section 2 of this article (Materials and Methods) via the modification/adaptation of the original expressions included in the IEEE Standard C57.110-2018 [24] and ANSI/UL Standards 1561 and 1562 [26,27].
The Section 2 of this article also describes the operational sequence and source code of the "Transformer Loss Calculator" (TLC) software, which enables the values of transformer load losses to be determined according to the IEEE Standard C57.110, with the separation of the values corresponding to the fundamental (active, reactive, and unbalance) and nonfundamental (distortion) frequency currents, which were obtained by applying the Fourier Series and Fortescue's theorem [30].
In the Section 3 of this article, we present the RMS values and the THDi z % of the currents recorded using a FLUKE 435 Series II analyzer in the second of two real 1000 kVA distribution transformers, which feed two Low-Voltage (LV) residential areas, which are very differently distorted loads, located in the Valencian Community (Spain).These results are applied in the Section 4 of this article to obtain the amount of daily and annual losses of each transformer according to the new approach established in this paper and the original method, which is included in the Standard.The separation between the losses at fundamental and non-fundamental frequencies, provided by the TLC software, has allowed for conclusions to be drawn in this section regarding (1) the effects of harmonics, both in terms of the value of losses and CO 2 emissions (measured in kg), and (2) the threshold in the value of THDi z %, above which the use of harmonic filtering is feasible from the point of view of environmental improvement.In addition, the limitations of our novel approach are analyzed in this section.
The Section 5 of this article presents our conclusions.

Materials and Methods
This section describes the procedure used to determine the environmental impacts (measured in kg of CO 2 emitted into the atmosphere) of the circulation of non-sinusoidal currents in the windings of three-phase distribution transformers.
This study is based on the fact that carbon dioxide emissions are proportional to the internal energy consumption (i.e., losses) of transformers.The proportionality factor is the emission factor (EF) of each region or country.The values of the losses caused by harmonic currents were obtained using the TLC software, which separately determines the values of the losses corresponding to the fundamental and non-fundamental frequency currents according to the IEEE Standard C57.110 [24].A novelty of the TLC software is that the load loss expressions defined by the IEEE Standard C57.110 have been extended to the case in which the currents have different values and harmonic contents (THDi z %) in each phase of the transformer windings.This is usually the case in distribution networks, so the TLC software is more useful than the direct application of the loss expressions included in the Standard.
The expressions of transformer load losses according to the ANSI standards are not implemented in the TLC software version 1.For this reason, in this section of the manuscript, we also establish the expressions for transformer load losses according to [26,27], adapted to the case in which the winding currents do not have the same RMS values and harmonic contents (THDi z %) in each phase.

Transformer Load Losses Adapted from IEEE C57.110 and ANSI/UL 1561-1562 Standards
The expressions of the load losses of three-phase transformers are developed in this section by modifying the IEEE Standard C57.110-2018 [24] and ANSI/UL Standards 1561 and 1562 [26,27], and they are adapted to the case in which the non-sinusoidal currents flowing through their windings have different RMS values and harmonic contents (THDi z %) in each phase.
The rated load losses (P LL−R ) of the transformers operating at the fundamental frequency (50-60 Hz) satisfy the following expression [24,25]: The composition of this equation is as follows: As the novelty of this paper, the equations for transformer load losses based on the IEEE Standard C57.110-2018 [24] were modified and adapted to the usual case in which the secondary currents have different RMS values and THDi z % in each phase (z = A, B, C), as follows: where P z JL represents the ohmic losses, P z ECL represents the Winding Eddy Current Losses, and P z OSL represents the Other Stray Losses of each phase (z = A, B, C).
Due to the Joule effect, the total ohmic losses (P JL ) in the transformer windings are obtained, under these conditions, using the following expression: The composition of this equation is as follows: • R DC is the short-circuit ohmic resistance of the transformer, calculated when direct currents are flowing through the windings of the transformer; • I hz is the RMS value of the harmonic current of order h z = f hz / f 1 and phase z = A, B, C.
Under these conditions, the Total Winding Eddy Current Losses (P ECL ) can be calculated as follows: where P ECL−R is the Rated Winding Eddy Current Losses, and the factor is the ratio of the Winding Eddy Current Losses caused by the harmonics in one of the phases (z = A, B, C) and the nominal losses, without Eddy and Other Stray effects.I hz is the RMS value of the harmonic current of order h z = f hz / f 1 circulating through that phase, and h z,max is the order of the highest harmonic frequency considered in the calculation.The Total Other Stray Losses (P OSL ) with non-sinusoidal currents of different RMS values and THDi z % in each phase of the transformer windings can be obtained as follows: where P OSL−R is the Rated Other Stray Losses, and the factor Is related to the Other Stray Losses caused by harmonics in the phases z = A, B, C to the nominal losses without Eddy and Other Stray effects.Likewise, I hz , h z = f hz / f 1 , and h z,max have the same meanings, as indicated above.
In the case where the manufacturer of the transformer provides the value of the Rated Total Stray Losses (P TSL−R ) but the values of P ECL−R and P OSL−R are not known, the IEEE Standard C57.12.90 TM [31] recommends considering that P ECL−R = 0.33•P TSL−R for oil-immersed transformers and P ECL−R = 0.66•P TSL−R for dry-type transformers.The ANSI/UL Standards 1561 and 1562 [26,27] consider Other Stray Losses (P OSL ) to be negligible in dry-type transformers, meaning that they do not directly affect the heating and lifetime of these transformers [29].For this reason, proceeding in the same way as in the previous section, the expressions of the total load losses (P LL ) proposed by the ANSI standards have been adapted to the habitual case in which non-sinusoidal currents have a different RMS value and THDi z % in each phase and are calculated as the sum of the ohmic losses (P z JL ) and the Winding Eddy Current Losses (P z JL ) of each phase (z = A, B, C), as follows: where, due to the Joule effect, the total ohmic losses (P JL ) can be obtained using Equation ( 4), and the Winding Eddy Current Losses (P ECL ) can be calculated as follows: where P ECL−R is the Rated Winding Eddy Current Losses, and the factor is defined for each phase (z = A, B, C) and coincides in value with F z ECL .

Transformer Loss Calculator Software
This section describes the operational sequence and source code of the Transformer Loss Calculator (TLC) software version 1, which was implemented with the LabVIEW platform (2015 version).LabVIEW software is used in numerous applications, such as in medical, chemical, and electrical disciplines, as well as others [32][33][34][35][36][37][38][39].

Operational Sequence of the TLC Software
The operational sequence of the TLC software version 1 depicted in Figure 1 comprises the following steps: Stage (1).Data input, either of the signals (voltages and currents) recorded on each phase of the secondary transformer by means of a physical measurement system or, alternatively, by the manual input of voltage and current values.In the first case, the software is used for the analysis of real transformers.In the second case, the software is used for simulation analysis.The rated values of the transformer are also entered in this step.
Stage (2).Fourier analysis of the signals introduced in stage (1).Stages (3) and (4).Obtaining the Complex Root Mean Square (CRMS) values of the fundamental frequency (50 Hz) and non-fundamental frequency (up to 10,150 Hz) voltages and currents.
Stage (5).Application of Fortescue's Theorem to fundamental frequency voltages and currents.
Stage (6).Obtaining the CRMS values of the components of the fundamental frequency positive-sequence voltages and currents.Stage (7).Determination of the CRMS values of the unbalanced voltage and current components.
Stage (8).Calculation of the CRMS values of the active currents, positive sequence, and fundamental frequency.
Stage (9).Obtaining the CRMS values of the reactive currents, those being the positivesequence and fundamental frequency currents.
Stage (14).Display on a screen the absolute and relative values (in %) of the total load losses of the transformer and individual values corresponding to each of the current components.Stage (5).Application of Fortescue's Theorem to fundamental frequency voltages a currents.
Stage (6).Obtaining the CRMS values of the components of the fundamental f quency positive-sequence voltages and currents.
Stage (7).Determination of the CRMS values of the unbalanced voltage and curre components.
Stage (8).Calculation of the CRMS values of the active currents, positive sequen and fundamental frequency.
Stage (9).Obtaining the CRMS values of the reactive currents, those being the po tive-sequence and fundamental frequency currents.

Programing Modules of TLC Software
The TLC software was developed in a programming module (Figure 2).We used the programming language indicated in Figure 3, which comprises the following stages:

Programing Modules of TLC Software
The TLC software was developed in a programming module (Figure 2).We used the programming language indicated in Figure 3, which comprises the following stages:  Stage (15): Data Input.At this stage, the samples of voltages and currents recorded in each secondary phase of the transformer using a physical measurement system are incorporated, or the CRMS values of these quantities are manually entered, in order to carry out analysis via simulation.In addition, the values of the transformer rated power ( ), load losses ( and  ), and no-load secondary line-to-line voltages ( ) are entered manually.
With the values entered at this stage, the nominal secondary current is obtained as follows: The TLC software was developed in a programming module (Figure 2).We used the programming language indicated in Figure 3, which comprises the following stages:  Stage (15): Data Input.At this stage, the samples of voltages and currents recorded in each secondary phase of the transformer using a physical measurement system are incorporated, or the CRMS values of these quantities are manually entered, in order to carry out analysis via simulation.In addition, the values of the transformer rated power ( ), load losses ( and  ), and no-load secondary line-to-line voltages ( ) are entered manually.
With the values entered at this stage, the nominal secondary current is obtained as follows: Stage ( 15): Data Input.At this stage, the samples of voltages and currents recorded in each secondary phase of the transformer using a physical measurement system are incorporated, or the CRMS values of these quantities are manually entered, in order to carry out analysis via simulation.In addition, the values of the transformer rated power (S R ), load losses (P JL−R and P TSL−R ), and no-load secondary line-to-line voltages (V s0 ) are entered manually.
With the values entered at this stage, the nominal secondary current is obtained as follows: Stage ( 16): Fourier Analysis.At this stage, the CRMS values of the transformer line-to-neutral voltages and line currents of fundamental (V 1z , I 1z ) and non-fundamental (V hz , I hz ) frequencies are obtained for each secondary phase (z = A, B, C).These values can be displayed on the software's screen (Figure 4) or transferred to Excel.
Stage ( 16): Fourier Analysis.At this stage, the CRMS values of the transformer lineto-neutral voltages and line currents of fundamental (  ,  ̅ ) and non-fundamental ( ,  ̅ ) frequencies are obtained for each secondary phase (,  = , , ).These values can be displayed on the software's screen (Figure 4) or transferred to Excel.
with  = 1∠120° as well as the CRMS values of the unbalanced line-to-neutral voltages ( ,  ,  ) and currents ( ̅ ,  ̅ ,  ̅ ) of each phase of the transformer secondary, Stage ( 18): RMS Currents.At this stage, the following RMS values of the transformer secondary currents and its components are calculated as follows:

•
Secondary current (): • Fundamental frequency current ( ) and its components active ( ), reactive ( ), and unbalanced ( ): • Non-fundamental frequency current ( ): Stage ( 19): Power Losses.The total load losses, as well as those caused by each of the secondary current components, are obtained at this stage.They are obtained in accordance Stage ( 17): Symmetrical Components.Fortescue's theorem is applied at this stage to the fundamental frequency voltages and currents obtained in the previous stage to determine the CRMS values of the fundamental frequency positive-sequence secondary line-to-neutral voltages (V + ) and line currents (I + ) as follows: with a = 1∠120 • as well as the CRMS values of the unbalanced line-to-neutral voltages (V Au , V Bu , V Cu ) and currents (I Au , I Bu , I Cu ) of each phase of the transformer secondary, Stage ( 18): RMS Currents.At this stage, the following RMS values of the transformer secondary currents and its components are calculated as follows:

•
Secondary current (I): • Fundamental frequency current (I 1 ) and its components active (I a ), reactive (I r ), and unbalanced (I u ): • Non-fundamental frequency current (I H ): Stage ( 19): Power Losses.The total load losses, as well as those caused by each of the secondary current components, are obtained at this stage.They are obtained in accordance with the expressions adapted from the IEEE Standard C57.110-2018, as seen in Section 2.1, as follows:

•
Total load losses: The terms corresponding to the fundamental frequency harmonics (h z = 1) in Equation ( 20) represent the transformer fundamental frequency load losses, which have three components: • Losses due to active currents (P LLa ), • Losses due to reactive currents (P LLr ), • Losses due to unbalanced currents (P LLu ), The load losses due to harmonics (P LLH ) are the fourth summand of (20): The relative values of the load losses are also determined at this stage, in %, caused by active currents (ε a %), reactive currents (ε r %), unbalanced currents (ε u %), and harmonic currents (ε H %), such as Step (20): Visualization.The RMS values of the currents, as well as the absolute and relative values of the load losses that these currents produce in the transformers, are displayed on a screen of the TLC software for visualization.

Load Losses Adapted from Standard ANSI/UL Standard 1561/1562
The TLC software operates with the load loss expression adapted from the IEEE Standard C57.110-2018.However, we also planned to incorporate the ANSI/UL Standards 1561 and 1562.In this section, based on what is indicated in Section 2.1.2,the equations that provide the load losses of dry-type and oil-immersed three-phase transformers are established when feeding non-linear loads of different values and non-linearity in each phase.
The load losses caused by the circulation of the currents of different RMS values and harmonic contents through the windings of three-phase transformers can be expressed, according to (10), as follows: The losses originated at the fundamental frequency (h z = 1) are where I 1 is the RMS value of the secondary currents of the fundamental frequency, determined by Equation (18), and I R is the transformer nominal secondary current, defined by Equation (13).These losses can be decomposed into active, reactive, and unbalanced losses, in a similar way to what was achieved in Equation (21).
The load losses caused by harmonic currents are expressed, according to the ANSI/UL Standards 1561 and 1562, as follows: Although the ANSI/UL Standards 1561 and 1562 refer to dry-type transformers, where Other Stray Losses (P OSL ) have practically no effect on winding losses, in our opinion, the last three equations could also be applied to the calculation of load losses in oil-immersed transformers.In Section 4, we will outline that the values obtained with Equations ( 27)-( 29) are close to those obtained with the IEEE Standard C57.110-2018 equations.

Results
The expressions of the load losses of the feeding loads of three-phase transformers of different values and non-linearity in each phase, developed in Sections 2.2.2 and 2.2.3, are applied in this section to determine the daily and annual load losses and CO 2 emissions of two oil-immersed Dyn11 distribution transformers when they supply the Low-Voltage (LV) installations of two very differently distorted residential areas.They were obtained from the manufacturer, ORMAZABAL, and their nominal characteristics are outlined in Table 1.

•
The first zone, supplied by transformer 1, is eminently residential, made up of dwellings, and has values of THDi z %, which are not too high in each phase (z = A, B, C);

•
The second zone, supplied by transformer 2, is also residential but has shops, with cold stores, LED, and discharge lighting, as well as small industrial companies, with high values of THDi z % in each phase.
The RMS values and harmonic contents of the secondary line currents of both transformers were recorded using a FLUKE 435 Series II network analyzer between the months of October and November 2022.

Losses and Daily CO 2 Emissions in Transformer 1
Table 2 shows the RMS values of the fundamental frequency line currents, as well as the THDi z % values of these currents registered by the FLUKE analyzer in each phase (z = A, B, C), of the secondary transformer 1 at one-hour intervals throughout the day of 13 November 2022 (a public holiday).Table 3 summarizes the RMS values of the fundamental frequency currents and the  % values of the line currents recorded using the FLUKE analyzer in each phase ( = , , ) of secondary transformer 1 at one-hour intervals throughout the day of 14 November 2022 (working day).Table 3 summarizes the RMS values of the fundamental frequency currents and the THDi z % values of the line currents recorded using the FLUKE analyzer in each phase (z = A, B, C) of secondary transformer 1 at one-hour intervals throughout the day of 14 November 2022 (working day).The waveforms and harmonic diagram of the currents at 16:02 h are represented in Figure 7a,b, respectively.Here, it can be seen that the currents supplied by transformer 1 to these installations on working days are also slightly distorted.The waveforms and harmonic diagram of the currents at 16:02 h are represented in Figure 7a,b, respectively.Here, it can be seen that the currents supplied by transformer 1 to these installations on working days are also slightly distorted.
The values of the load losses in transformer 1, calculated throughout the day of 14 November 2022 for the fundamental frequency currents and the first 24 harmonics of nonfundamental frequency (harmonic currents), adapting the IEEE Standard C57.110 (IEEE-1, IEEE-H

Losses and Daily CO2 Emissions from Transformer 2
The RMS values of the fundamental frequency currents recorded on 9 October 2022 (a public holiday) in each phase ( = , , ) of secondary transformer 2, as well as the  % values, are summarized in Table 4.

Losses and Daily CO 2 Emissions from Transformer 2
The RMS values of the fundamental frequency currents recorded on 9 October 2022 (a public holiday) in each phase (z = A, B, C) of secondary transformer 2, as well as the THDi z % values, are summarized in Table 4.        Table 5 summarizes the RMS values of the fundamental frequency and the  % values of the currents registered in each phase of secondary transformer 2 throughout 10 October 2022 (working day).Table 5 summarizes the RMS values of the fundamental frequency and the THDi z % values of the currents registered in each phase of secondary transformer 2 throughout 10 October 2022 (working day).The waveform and harmonic diagram of the current of phase A of transformer 2, recorded at 11:00 on 10 October 2022, are represented in Figure 11a,b.It can be seen in these figures that the currents continue to be highly distorted on weekdays.

Discussion
In the previous section, the values of the load losses of two distribution transformers supplying residential facilities in a town near the city of Valencia (Spain) were presented.The losses were obtained on working days and holidays by entering the measurements recorded using the FLUKE 435 Series II analyzer in the TLC software, whose algorithms use expressions for load losses according to the IEEE Standard C57.110, but this was adapted to the usual case where the currents have different RMS values and  % in each phase.Likewise, in this paper, the load losses were also calculated according to the ANSI/UL Standards 1561 and 1562, adapting the expressions of this standard, as indicated in Section 2.2.3.
Although both installations are residential, there are important differences in terms of the linearity of the loads in each of them.Practically all of the installations powered by transformer 1 are residential, the loads are evenly distributed among the three phases, and the supplied currents have a very low harmonic content, with a  % value of around 10% during central hours on weekdays and holidays (Tables 2 and 3).The installations powered by transformer 2 are mainly residential, but there is also a significant number of shops and small industries whose loads introduce high harmonic contents, meaning there are currents with a  % greater than 70% at many times on weekdays and public holidays (Tables 4 and 5).The loads connected to transformer 2 are not evenly distributed, with phase C being the most heavily loaded and with the lowest harmonic content.It can also be seen that currents with a higher  % occur the smaller the supplies of transformer 2 are, which is a symptom of the presence of important non-linear loads that are

Discussion
In the previous section, the values of the load losses of two distribution transformers supplying residential facilities in a town near the city of Valencia (Spain) were presented.The losses were obtained on working days and holidays by entering the measurements recorded using the FLUKE 435 Series II analyzer in the TLC software, whose algorithms use expressions for load losses according to the IEEE Standard C57.110, but this was adapted to the usual case where the currents have different RMS values and THDi z % in each phase.Likewise, in this paper, the load losses were also calculated according to the ANSI/UL Standards 1561 and 1562, adapting the expressions of this standard, as indicated in Section 2.2.3.
Although both installations are residential, there are important differences in terms of the linearity of the loads in each of them.Practically all of the installations powered by transformer 1 are residential, the loads are evenly distributed among the three phases, and the supplied currents have a very low harmonic content, with a THDi z % value of around 10% during central hours on weekdays and holidays (Tables 2 and 3).The installations powered by transformer 2 are mainly residential, but there is also a significant number of shops and small industries whose loads introduce high harmonic contents, meaning there are currents with a THDi z % greater than 70% at many times on weekdays and public holidays (Tables 4 and 5).The loads connected to transformer 2 are not evenly distributed, with phase C being the most heavily loaded and with the lowest harmonic content.It can also be seen that currents with a higher THDi z % occur the smaller the supplies of transformer 2 are, which is a symptom of the presence of important non-linear loads that are permanently connected in phases A and B to that transformer.
Table 6 summarizes the values of the energy losses obtained in transformer 1 on 13 November (public holiday) and 14 (working days), and in transformer 2 on 9 October (public holiday) and 10 (working day) of 2022, using the adapted and original load loss expressions for the standards.The values for the original load loss expressions were obtained using the average RMS values of the currents in the three phases.It can be seen in Table 6 there is not a proportional relationship between the distortion rate of the currents (THDi z %) and the losses they produce in the transformers.Thus, the losses in transformer 1, due to the harmonics, are barely 4.3% of the total, even though the average distortion rate of the currents in that transformer (THDi z %) is around 20%.In contrast, the losses due to harmonics in transformer 2, with average distortion rates of 40-50%, are comparable to the losses caused by fundamental frequency currents (Table 6).The differences between the daily energy losses obtained according to the adapted IEEE C57.110 and ANSI/UL 1561-1562 standards constitute the Other Stray Losses.From the values indicated in Table 6, it can be seen that the Other Stray Losses in transformer 1 occur mainly at the fundamental frequency.
In Table 7 are summarized the lost energies in transformers 1 and 2 during 2022, considering that in that year there were 248 working days and 117 public holidays (with Saturdays and Sundays included as public holidays).The values of CO 2 emissions, in kg, produced by these transformers in 2022 (Table 8), were obtained by multiplying the energy losses indicated in Table 7 by the emission factor, which is 0.154 kg CO 2 /kWh consumed in the Valencian Community [13].
As seen in Tables 7 and 8, using adapted expressions for the load losses developed in this paper, it is verified that the elimination of harmonic currents, by using suitable harmonic filters, could reduce the losses by 183.68 kWh and reduce the CO 2 emissions by 28.29 kg (3.2% of the total) in transformer 1 in 2022, of which only 13.56 kWh and 2.09 kg of CO 2 emissions correspond to the Other Stray Losses (0.24% of the total).In transformer 2, with high current distortion rates, the reduction in losses and CO 2 emissions could be significant, at 10,053.18 kWh and 1548.19 kg (45.11% of the total), respectively, although only 726.91 kWh and 111.95 kg would be due to Other Stray Losses (3.26% of the total).
From Tables 7 and 8, it can also be noted that, in the case of transformer 2, the total energy losses and CO 2 emissions obtained annually from the adapted load loss expressions is 5.92% greater than that obtained according to the original expressions for the Standard IEEE C57.110 (5.55% for the ANSI/UL Standards 1561-1562).However, our approach, like that used by the standards, has limitations.There are small differences between the loss values obtained with our approach and the actual transformer losses derived from the measurement procedure.To obtain the real value of the losses in the transformers, it is necessary to use two wattmeters: one arranged in the primary and the other in the secondary.The actual load losses would be obtained by subtracting the no-load losses from the difference between the active powers measured by the wattmeters in the primary and secondary transformers.This measurement procedure is accurate but more laborious and expensive to implement, especially when installing the wattmeter on the high-voltage side.
Like that used by the standards, our approach is based on measuring only one of the transformer windings.With this, we are making an approximation, since the effects of the circulation of vacuum currents through the primary coils have been neglected.As the vacuum currents are mainly of a fundamental frequency and their effective values do not usually change during the operation of the transformers, the maximum error that can be made in the measurement of the losses with the use of the TLC software can be estimated based on the square of the quotient between the RMS value of the no-load current and the RMS value of the nominal current of the transformer (at the nominal frequency) multiplied by 100.The manufacturers provide the value of this quotient.The load loss maximum errors determined in the 1000 kVA distribution transformers analyzed in the practical part of this article, when using our approach, would be less than 0.011%; therefore, the limitations regarding the use of the new expressions of the losses established in the article can be dismissed.

Conclusions
Transformers are important vectors in energy waste and CO 2 emissions in electrical networks.These effects are aggravated when these machines feed non-linear loads.For this reason, new expressions have been established in this article to calculate the load losses of three-phase transformers obtained by the modification/adaptation of the original expressions included in the EEE Standard C57.110-2018 and ANSI/UL Standards 1561 and 1562.The limitations of these new expressions, derived from the fact that the calculation of the losses is carried out based on only one winding current of the transformer and not on both at the same time, as that would be more exact, can be neglected due to the minor errors introduced in the value of the losses (less than 0.011% in 1000 kVA transformers).
The novel expressions of the load losses obtained by adapting the original IEEE Standard C57.110 have been implemented using the TLC calculation version 1 software.Comparisons between the annual values of energy losses obtained with this software and those resulting from the direct application of the IEEE standard have determined differences of up to 5.92% in favor of the former.This shows that it is advisable to incorporate the TLC software in monitoring instruments to operate three-phase transformers.
The study of load losses and CO 2 emissions carried out on real transformers supplying residential areas with Low-Voltage distribution networks has revealed the following:

•
The values of losses and CO 2 emissions obtained according to the IEEE C57.110 standard are always higher than those calculated using the ANSI/UL Standards 1651 and 1562 because the latter standards do not include the effects of Other Stray Losses.However, the differences between the results of those standards are not too large; they are less than 7%, even with a high THDi z %.

•
The current distortion rates (THDi z %), between 10 and 20%, cause only very slight load losses and CO 2 emissions, which are less than 5% of the total.

•
The current distortion rates (THDi z %) of 40-60% produce high losses and CO 2 emissions, which can be higher than 40% of the total.
Thus, THDi z % can be used as a helpful parameter to inform the necessity of using harmonic filters to improve sustainability, reducing losses and CO 2 emissions from transformers.
Finally, in the future, we intend to improve the TLC software and provide information about other parameters that affect the aging and derating of transformers.

2. 1 . 2 .
Load Losses Adapted from the ANSI/UL Std.1561 and 1562 . Fourier analysis of the signals introduced in stage (1).Stages (3) and (4).Obtaining the Complex Root Mean Square (CRMS) values of t fundamental frequency (50 Hz) and non-fundamental frequency (up to 10,150 Hz) vo ages and currents.

Figure 1 .
Figure 1.Operational sequence of the TLC software version 1.

Figure 1 .
Figure 1.Operational sequence of the TLC software version 1.

Figure 3 .
Figure 3. Source code of software TLC.

Figure 3 .
Figure 3. Source code of software TLC.

Figure 3 .
Figure 3. Source code of software TLC.

Figure 4 .
Figure 4.A partial view of the harmonic voltage and current screen reported by the TLC software.Stage (17): Symmetrical Components.Fortescue's theorem is applied at this stage to the fundamental frequency voltages and currents obtained in the previous stage to determine the CRMS values of the fundamental frequency positive-sequence secondary lineto-neutral voltages ( ) and line currents ( ̅ ) as follows:

Figure 4 .
Figure 4.A partial view of the harmonic voltage and current screen reported by the TLC software.

Figure
Figure 5a,b show the waveforms and the harmonic diagram of the currents at 14:24 h.These currents are slightly distorted.

Figure
Figure 5a,b show the waveforms and the harmonic diagram of the currents at 14:24 h.These currents are slightly distorted.Figure 6a,b show the values of transformer 1′s load losses, which were obtained every hour on 13 November 2022 for the first 25 harmonic currents, adapting the IEEE Standard C57.110 (IEEE-1, IEEE-H) and the ANSI/UL Standards 1561-1562 (ANSI-1, ANSI-H).The differences between the values of the losses, according to both standards (IEEE and ANSI), shown in Figure 6a,b, are the Other Stray Losses ( ).

FigureFigure 6 .
Figure 6a,b show the values of transformer 1's load losses, which were obtained every hour on 13 November 2022 for the first 25 harmonic currents, adapting the IEEE Standard C57.110 (IEEE-1, IEEE-H) and the ANSI/UL Standards 1561-1562 (ANSI-1, ANSI-H).The differences between the values of the losses, according to both standards (IEEE and ANSI), shown in Figure 6a,b, are the Other Stray Losses (P OSL ).Environments 2023, 10, 177 12 of 23

Figure
Figure 9a,b show, respectively, the waveform and the harmonic diagram of the Aphase secondary current of transformer 2 at 2:00 on 9 October 2022.It can be observed from these figures that these currents are strongly distorted.

Figure
Figure 9a,b show, respectively, the waveform and the harmonic diagram of the Aphase secondary current of transformer 2 at 2:00 on 9 October 2022.It can be observed from these figures that these currents are strongly distorted.Figure 10a,b show the values of losses in transformer 2 throughout 9 October 2022 (a public holiday) for the fundamental frequency currents and the first 24 harmonics of nonfundamental frequency, adapting the IEEE Standard C57.110 and ANSI/UL Standard 1651-1562.

Figure
Figure 10a,b show the values of losses in transformer 2 throughout 9 October 2022 (a public holiday) for the fundamental frequency currents and the first 24 harmonics of non-fundamental frequency, adapting the IEEE Standard C57.110 and ANSI/UL Standard 1651-1562.
The waveform and harmonic diagram of the current of phase A of transformer 2, recorded at 11:00 on 10 October 2022, are represented in Figure11a,b.It can be seen in these figures that the currents continue to be highly distorted on weekdays.The values of the load losses calculated in transformer 2 throughout 10 October 2022 (working day) for the fundamental frequency currents and the first 24 harmonics of nonfundamental frequency are represented in Figure12a,b, adapting the IEEE Standard C57.110 and ANSI/UL Standard 1561-1562.

Figure 11 .
Figure 11.A-phase line current of transformer 2 at 11:00 (10 October 2022): (a) waveform; (b) harmonic diagram.The values of the load losses calculated in transformer 2 throughout 10 October 2022 (working day) for the fundamental frequency currents and the first 24 harmonics of nonfundamental frequency are represented in Figure 12a,b, adapting the IEEE Standard C57.110 and ANSI/UL Standard 1561-1562.
2 R are the nominal cooper (Joule) losses, as a function of the ohmic short-circuit resistance of the transformer (R DC ) and the RMS value of the secondary rated current (I R ); • P ECL−R are the Rated Winding Eddy Current Losses; • P OSL−R are the Rated Other Stray Losses; • P TSL−R = P ECL−R + P OSL−R are the Rated Total Stray Losses.

Table 1 .
Rated values of the ORMAZABAL transformers 1 and 2.

Table 2 .
Fundamental (50 Hz) RMS and THDi z % values of the secondary line currents of transformer 1 measured each hour starting from 0:02 on 13 November 2022 (public holiday).

Table 2 .
Fundamental (50 Hz) RMS and  % values of the secondary line currents of transformer 1 measured each hour starting from 0:02 on 13 November 2022 (public holiday).

Table 3 .
Fundamental (50 Hz) RMS and  % values of the secondary line currents of transformer 1 measured each hour starting from 0:02 on 14 November 2022 (working day).

Table 3 .
Fundamental (50 Hz) RMS and THDi z % values of the secondary line currents of transformer 1 measured each hour starting from 0:02 on 14 November 2022 (working day).

Table 4 .
Fundamental (50 Hz) RMS and  % values of the secondary line currents of transformer 2 measured each hour starting from 0:00 on 9 October 2022 (public holiday).

Table 4 .
Fundamental (50 Hz) RMS and THDi z % values of the secondary line currents of transformer 2 measured each hour starting from 0:00 on 9 October 2022 (public holiday).

Table 6 .
Daily energy losses, in kWh, caused by fundamental and harmonic currents in transformers 1 and 2, according to the adapted and original (within parentheses) expressions for the IEEE C57.110 and ANSI/UL 1561-1562 standards.

Table 7 .
Total energy losses, in kWh, caused in the year 2022 by the fundamental and harmonic currents in transformers 1 and 2 according to the adapted and original (within parentheses) expressions for the IEEE C57.110 and ANSI/UL 1561-1562 standards.

Table 8 .
Total CO 2 emissions, in kg, caused in the year 2022 by the fundamental and harmonic load losses in transformers 1 and 2 according to the adapted and original (within parentheses) expressions for the IEEE C57.110 and ANSI/UL 1561-1562 standards.