Theoretical and Experimental Study to Determine Voltage Violation, Reverse Electric Current and Losses in Prosumers Connected to Low-Voltage Power Grid

: The impact of PV generation distributed in a low voltage transmission line depends on many factors: The distribution lines and PV generators characteristics, its location, operational control, local meteorological conditions, electricity consumption proﬁle, and the electricity cost variation. An atypical and challenging behavior of photovoltaic distributed generation (DG) insertion in consumer units (CUs), implies in some circumstances, as the reverse directionality of the power ﬂow between the load equipped with a photovoltaic system generator and the electrical grid, when a CU contains a distributed generation and low power consumption, the power ﬂow will be directed to the power electric grid. In this work, the modeling of a low-voltage real feeder was performed, setting the variables of the system under real operating conditions. As result, voltage levels variability throughout the feeder, the electrical losses, and the asymmetry between the phases were observed. Through simulation scenarios, the occurrence of voltage increase under di ﬀ erent penetration scenarios of distributed generation was veriﬁed and there was a 10% increase in reference voltage as well as the occurrence of higher electrical losses by reverse current, reaching 1200% more with a DG penetration, in the massive presence of the photovoltaic generator. The mitigatory action used in this work was able to attenuate the negative impacts to the feeder circuit, ensuring the integrity grid and the consumer unit.


Introduction
Nowadays, the technological advance increasingly demands electricity to supply its needs, whether at a residential, industrial, or commercial level. Even with the environmental awareness and efforts to rationalize the use of electricity, there is a strong growth of electronic equipment, such as cellphones, televisions, air conditioning, and others in urban environments. Energy complementarity is a way of aggregating power from different sources to the electrical matrix, supplying the demand on moments of scarcity. Photovoltaic solar energy emerges as one of the options to mitigate this problem [1]. The development of this energetic modality in Brazil occurred with the creation of the Resolution 414/2012, added and renumbered to 687/2015, along with the creation of module 8 of distribution procedure. These two normative documents are responsible for regulating energy insertion, generation of electricity connected to the network, and monitoring the quality of the electricity [2]. Due to this standardization process and to the enormous potential for solar resource utilization in Brazil, in the year 2017, from the 10,561 connections performed until the month of May of that year [3], 10,453 came The data mentioned in the table above also refer to the type of service extension of the consumer unit (CU), that in turn depends on the installed demand in the residence; therefore, Table 2 associates the relation between installed power and service voltage. Table 2. Service voltage in function of local demand.

Connection Installed Power Voltage
One-phase P ≤ 15 kW 220 V Three-phase 15 ≤ P ≤ 75 kW 380 V The values shown in Table 2 vary according to the regions and the electrical utilities. The data mentioned in the respective table came from the local electrical utility regulations [14]. The penetration of distributed generation (DG) was conditioned to the limits of power established by local regulation of the electricity company [15], where each DG may vary the power up to the one-phase maximum limit P ≤ 15 kW (100%).

Influence of Photovoltaic Systems in the Voltage of a CU
The electric power distribution system was initially conceived to deliver power to a consumer radially. With the creation of the normative resolution 687/2015 from the National Electric Energy Agency (ANEEL), the point of view regarding electricity consumers is under new definition, classifying a CU as an electricity prosumer, which ensures that it goes on to produce locally its own electricity [2]. Through the illustrations in Figure 1a, it is possible to visualize the unidirectional flow dynamics (conventional) and in Figure 1b, the reverse flow in a generic section of a feeder. The images simplify the connection of multiple CUs to a feeder knot (V ) through the extension impedance Z . The section impedance, Z , interconnects the bars to the distribution substation. Depending on the type of feeder, the line impedance (Z ) and the service extension impedance are represented as active and reactive components (R + jX), characterized by the cable type and the total distance of the network section. Similarly, the demand (P ) and photovoltaic generators (P ) are  The images simplify the connection of multiple CUs to a feeder knot (V 2 ) through the extension impedance Z 2n . The section impedance, Z L1 , interconnects the bars to the distribution substation. Depending on the type of feeder, the line impedance (Z L1 ) and the service extension impedance are represented as active and reactive components (R + jX), characterized by the cable type and the total distance of the network section. Similarly, the demand (P UC ) and photovoltaic generators (P PV ) are represented by the complex power given by S = P + jQ. The direction of the current flow (I 2n ) is variable, changing its direction when there is a photovoltaic generator in the CUs. When P PV = 0, there is a voltage drop in the service extension (V sn ) given by Equation (1). V sn = i 2n ·Z sn = i 2n · r sn + jx sn = V rn + jV xn (1) When P PV P UC , a power in reverse direction will emerge by the impedance Z sn , the current flowing through the section is adopted as −I 2n ; from that, Equation (2) is obtained considering the voltage module on the side of the CU. V sn = −i 2n ·Z sn = −i 2n · r sn + jx sn = − P PV V 2n ·r sn + Q UC V 2n ·x sn (2) The negative sign in the first portion merely indicates that the increase of voltage will occur on the side of the consumer. Inserting the real component of Equation (2) in the respective term of Equation (3a), it is possible to calculate the voltage value in the side of the consumer when the power supplied by the photovoltaic system exceeds the local demand.
Analyzing Equation (3a), considering that the demand profile is constant throughout the day, too much power increase from the photovoltaic generator will inevitably cause an increase in the voltage on the load bar. The phenomenon mentioned can be seen in Figure 2.

Real Feeder Identification
A low-voltage distribution system can operate with levels of three-phase 380 V and one-phase 220 V, and commonly presents an aerial topology consisting in bars and conductors, carrying power to four wires in the three-phase star model, when there is double voltage. Circuit energization is conditioned to a substation with primary feeding of 13,800 V in the delta module and maximum power of 75 kVA. Figure 3 below illustrates a real feeder.

Real Feeder Identification
A low-voltage distribution system can operate with levels of three-phase 380 V and one-phase 220 V, and commonly presents an aerial topology consisting in bars and conductors, carrying power to four wires in the three-phase star model, when there is double voltage. Circuit energization is conditioned to a substation with primary feeding of 13,800 V in the delta module and maximum power of 75 kVA. Figure 3 below illustrates a real feeder.  Figure 2 illustrates voltage behavior in a low voltage service extension, the occurrence of overvoltage due to the excess of reverse power flow, and structural physical conditions of the electrical circuit. The voltage rise may spread throughout the feeder in case more DGs are connected to the main knot of the feeder, as well as if the weather conditions are favorable to overcome nominal values of photovoltaic generators.

Real Feeder Identification
A low-voltage distribution system can operate with levels of three-phase 380 V and one-phase 220 V, and commonly presents an aerial topology consisting in bars and conductors, carrying power to four wires in the three-phase star model, when there is double voltage. Circuit energization is conditioned to a substation with primary feeding of 13,800 V in the delta module and maximum power of 75 kVA. Figure 3 below illustrates a real feeder. The distribution lines were modeled taking into consideration the layout and real characteristics of the cables and spacing between the bars (all raised in the field), and so the nodal admittance matrices were defined from the information of their sequential impedance, resistance, and reactance. The standardization of bars and section of the network was declared in one of the subroutines from the main script of the code. Geometry modeling of the sustaining structures and isolation of the conductors follows references found in the regulation of the local energy distributor [16]. The feeder supplies energy to 202 residential consumers, distributed in the 20 bars of the system, where 189 are connected to the one-phase circuit and the others connected to the three-phase circuit.

Simulation Data Acquisition
Monitoring PV plants is an urgent activity for practical implementation due to the information that allows owners to maintain, operate, and control these systems. The work done by Beránek [17] focuses on designing and testing a unique monitoring system with third-generation software for photovoltaic arrays. The distribution lines were modeled taking into consideration the layout and real characteristics of the cables and spacing between the bars (all raised in the field), and so the nodal admittance matrices were defined from the information of their sequential impedance, resistance, and reactance. The standardization of bars and section of the network was declared in one of the subroutines from the main script of the code. Geometry modeling of the sustaining structures and isolation of the conductors follows references found in the regulation of the local energy distributor [16]. The feeder supplies energy to 202 residential consumers, distributed in the 20 bars of the system, where 189 are connected to the one-phase circuit and the others connected to the three-phase circuit.

Simulation Data Acquisition
Monitoring PV plants is an urgent activity for practical implementation due to the information that allows owners to maintain, operate, and control these systems. The work done by Beránek [17] focuses on designing and testing a unique monitoring system with third-generation software for photovoltaic arrays.
Simulation standardization was performed in the OPENDSS environment in temporal mode. In this work, there was the concern to visualize the voltage behavior at the service extension of a consumer unit (CU) facing the insertion of photovoltaic generators, in order to confront tolerance limits pre-established by regulation. For that, a residence located in the feeder at a final section of the network (section 08) was considered, at the last bar (B9), properly instrumented and with data acquisition system (SCADA). Through the SCADA records, it was possible to determine the input profiles for the simulation. Figure 4 illustrates the composition of the experimental environment, as well as the layout of the instrumentation used. Simulation standardization was performed in the OPENDSS environment in temporal mode. In this work, there was the concern to visualize the voltage behavior at the service extension of a consumer unit (CU) facing the insertion of photovoltaic generators, in order to confront tolerance limits pre-established by regulation. For that, a residence located in the feeder at a final section of the network (section 08) was considered, at the last bar (B9), properly instrumented and with data acquisition system (SCADA). Through the SCADA records, it was possible to determine the input profiles for the simulation. Figure 4 illustrates the composition of the experimental environment, as well as the layout of the instrumentation used.  The power systems elements used in OPENDSS, such as PVSystem and Load, presented models that used the monitored variable illustrated in Figure 4b as input values. Data acquisition was adjusted to register data every one minute, enabling the assembling of the solar irradiation loadshapes, demand, and operating temperature of the photovoltaic module.  The power systems elements used in OPENDSS, such as PVSystem and Load, presented models that used the monitored variable illustrated in Figure 4b as input values. Data acquisition was adjusted to register data every one minute, enabling the assembling of the solar irradiation loadshapes, demand, and operating temperature of the photovoltaic module.

Solar Irradiance Data and Operating Temperature of PV Panel
The energy quality from the photovoltaic systems is closely linked to the quantity of incident solar energy on the PV panels surface, the operating temperature of the cell, and panel efficiency. The PV System element of the OPENDSS requires irradiance and panel temperature input data as input quantities in its mathematical model to estimate daily electrical power delivered to the power system to which it is connected. The solar irradiation profile was collected through a pyranometer and a temperature sensor installed together to the photovoltaic modules, both connected to the Scientific Campbell CR1000 datalogger.
Irradiation and temperature loadshapes operate together with the basis quantities of irradiation and operating temperature of the cell, being adopted in such quantities as 1 kW/m 2 and 25 • C, respectively.

Demand Curves Acquisition
The power distribution system studied supplies 202 CUs, framed as residential consumers, being 83% of the clients connected to the one-phase circuit. The load curves adopted were obtained from a demand monitor, manufacturer ISSO and reference DMI P100. For simulation purposes, two load curves from consumer units linked to the busbar 09, shown previously in Figure 3, were measured. Characteristically, the feeder is besieged by single-family residences, where during most of the day, the load curve operates in region of minimal demand. Figure 5 presents two demand profiles, framed in active power. Taking as reference the curves from the graph above, there are two particular situations that are commonly attributed to residential loads. Reflecting the blue curve to a real situation, it is characterized as a field house, eventually occupied. The brown curve presents an urban behavior, where the demand peaks occur at night. It is important to stress that both demands can belong to a single family group and, fostered by normative resolution 482/2012 ANEEL, local power use can be contemplated by a distributed generation installed remotely in another CU.

Problem Formulation and Methodology
For the study of power flow in distribution networks, the combination of four variables referenced in Table 3 was performed, where two of them are input parameters and the other two incognitos. The bars types definitions are also present, as well as their identification in the studied feeder.  Taking as reference the curves from the graph above, there are two particular situations that are commonly attributed to residential loads. Reflecting the blue curve to a real situation, it is characterized as a field house, eventually occupied. The brown curve presents an urban behavior, where the demand peaks occur at night. It is important to stress that both demands can belong to a single family group and, fostered by normative resolution 482/2012 ANEEL, local power use can be contemplated by a distributed generation installed remotely in another CU.

Problem Formulation and Methodology
For the study of power flow in distribution networks, the combination of four variables referenced in Table 3 was performed, where two of them are input parameters and the other two incognitos. The bars types definitions are also present, as well as their identification in the studied feeder. The load flow calculation involves determining the state of the network and consists on knowing the behavior in permanent regime based on the fundamental frequency, bars voltages, currents, and power injected in the elements. Traditionally, power flow calculation for transmission networks is based on interactive numerical methods such as Newton-Raphson and Gaus-Seidel; however, when reflected in the distribution feeders, these devices lose effectiveness due to mesh operation, high R/X (Resistance/Reactance) ratio, and imbalance [18]. To perform the simulation, it is necessary to name the bars; when there are load bars, it is called bar PQ (active power and reactive power); for a source bar or reference, we assign the bar Vθ (voltage module and voltage angle). The simulation performed in this article will be developed in the OPENDSS-EPRI open source environment, where the network modeling is based on the feeder nodal admittance matrix, and can be seen in Equation (3b).

Network Electrical Parameters Modeling
The electrical representations of the transmission lines involve a series of factors and, in general, the study is performed in two configurations: Curved lines and long lines, being selected from the level of precision required by the study [19]. In this article, a long line model was chosen, generally used in large extension transmission lines, and is illustrated in Figure 6  Circuit interpretation in Figure 6 reveals there is a power transfer from the k bar to the m bar, through the TL (transmission line) with line impedance predominantly reactive. The presence of shunts elements and capacitors, are visualized in the extremities of the circuit, which are the result of the electric field effect between phase and guard conductors and the soil.
The mathematical representation of TL parameters is seen in Figure 6, and shown in Equation (4), separated in function of the resistive and reactive components of the network.
Shunt capacitive impedance is tiny, and therefore is represented by its inverse, the susceptance. Shunt admittance is divided by half in each extremity and given by Equation (5). Circuit interpretation in Figure 6 reveals there is a power transfer from the k bar to the m bar, through the TL (transmission line) with line impedance predominantly reactive. The presence of shunts elements and capacitors, are visualized in the extremities of the circuit, which are the result of the electric field effect between phase and guard conductors and the soil.
The mathematical representation of TL parameters is seen in Figure 6, and shown in Equation (4), separated in function of the resistive and reactive components of the network. Shunt capacitive impedance is tiny, and therefore is represented by its inverse, the susceptance. Shunt admittance is divided by half in each extremity and given by Equation (5).
In the feeder analyzed, as shown in Figure 3, the impedance of each network section is variable. This happens by economic strategy, because if the feeder were scaled with a uniform cross section throughout the circuit, it would imply high financial cost. For the purpose of modeling the electrical parameters of the feeder, the real impedance of each network section was considered, as well as the service extension that interconnects each CU to the circuit busbar.
Tables 4 and 5 present the geometrical configuration of the feeder under study. In Appendix A of this document are the impedance values considered in the network modeling, represented by the impedance matrix.

Photovoltaic System Modeling
The PVSystem element of OPENDSS is an energy converter, responsible for injecting active power to the bar, which in turn depends on solar irradiation, operating temperature of the photovoltaic panel, and inverter efficiency. Radatz [18] presented in his work the set of combined equations that result in the mathematical model used by the OPENDSS software to stimulate the electric power generated by a generator. Equation (6) estimates the power produced by the photovoltaic array, at each instant of time, using the Loadshapes presented in Figure 5 as a basis. Equation (6) calculates instantaneously the power produced by the PV panels, where Pmp base and I base are base quantities referring to the maximum power point of the generator and solar irradiance, respectively. I p.u. is the solar irradiance profile normalized by the reference radiation (1kWp/m 2 ) throughout the day, shown previously in Figure 5. Power correction at the output of the PV's is presented by the factor Pmp (t) factor . The power effectively injected in the electrical network is given by Equation (7), also being expressed in function of time by the quantity Pca (t) .
The calculation of the power injected to the grid uses as basis the power produced by the PV panels, being dependent on the efficiency of the inverter to each load range. In this study, the photovoltaic generator modeling obeyed the conditions of a real system, where the electric companies establish a maximum power value for the connection of micro-generators to the grid, and such criteria are pre-established in [2].

Case Study
The main objective of this study was limited to a punctual analysis of section 8 of the feeder, which connects bar 08 to bar 09, and is responsible to supply six one-phase residential consumers and symmetrically distributed between the three phases available by the circuit of the feeder. Table 6 details all the configuration and layout of the CUs. Upon technical requirements to access the power grid [2], a prosumer can produce its local electricity for self-consumption, being able to eventually export its energy surpluses for remote consumption. The combination of some factors as the conditions of solar irradiation, demand curve, and physical characteristics of service extensions can change the regular dynamic of the voltage profile, causing overvoltage at the point of delivery of the electric utility and unbalance in the three-phase network system. Currently, the inverters are equipped with protection systems that ensure system interruption (curtailment) given the occurrence of this phenomeno; in contrast, the unwanted shutdown of the photovoltaic system will imply in efficiency loss and more operational costs.
In order to identify these technical adversities in advance, in Figure 7 is shown a flowchart that performs a study of electrical voltage variability on permanent basis at the point of delivery of the CU, in function of the power installed of the photovoltaic DG. Simply put, the flowchart solves the load flow problem, returning with voltage values at the bars of the power grid and the CUs (Step 3); the values calculated are confronted with the established limits (Step 4); in case the values exceeds, the flowchart will return to Step 2 in order to reset the parameterization, restarting the load flow solution (Step 3). The algorithm is finished in the case of values remain at rated operating conditions (Step 5).
CU, in function of the power installed of the photovoltaic DG. Simply put, the flowchart solves the load flow problem, returning with voltage values at the bars of the power grid and the CUs (Step 3); the values calculated are confronted with the established limits (Step 4); in case the values exceeds, the flowchart will return to Step 2 in order to reset the parameterization, restarting the load flow solution (Step 3). The algorithm is finished in the case of values remain at rated operating conditions (Step 5). The problem is actually solved in the block that equates and analyzes the power flow, that is Step 3, identifying the potentials in the desired bars; with that, in Step 4, the voltage is verified according to the limits already mentioned in [14]. If the voltage limit is exceeded, it is necessary to make decisions that enable the insertion of photovoltaic generator safely to the PES (power electric system).
Simulation was performed through the OPENDSS-EPRI software, open source. The solution of the power flow is obtained through the interactive algorithm for the solution of load flow based on the Newton method. The input parameters for the simulation were obtained from a data acquisition system, which experimentally monitored the operating variables of the photovoltaic generator, the CU and the environmental conditions.
The simulation cases were selected from several combinations between input parameters. Next, the simulation combinations used to obtain the results are presented. The scenario represented by Figure 8 illustrates typical behavior of remote residential consumers, commonly idle during the entire day for a long monthly period, which can be compared to a sporadically inhabited rural property. The profile illustrated also represents single families CUs that are operationally idle during the day, usually having a nightly over-demand. The problem is actually solved in the block that equates and analyzes the power flow, that is Step 3, identifying the potentials in the desired bars; with that, in Step 4, the voltage is verified according to the limits already mentioned in [14]. If the voltage limit is exceeded, it is necessary to make decisions that enable the insertion of photovoltaic generator safely to the PES (power electric system).
Simulation was performed through the OPENDSS-EPRI software, open source. The solution of the power flow is obtained through the interactive algorithm for the solution of load flow based on the Newton method. The input parameters for the simulation were obtained from a data acquisition system, which experimentally monitored the operating variables of the photovoltaic generator, the CU and the environmental conditions.
The simulation cases were selected from several combinations between input parameters. Next, the simulation combinations used to obtain the results are presented. The scenario represented by Figure 8 illustrates typical behavior of remote residential consumers, commonly idle during the entire day for a long monthly period, which can be compared to a sporadically inhabited rural property. The profile illustrated also represents single families CUs that are operationally idle during the day, usually having a nightly over-demand. It is common that residential CUs do not support photovoltaic system, especially because there is no area for its installation (area occupation factor), but thanks to the power sharing mode, a photovoltaic system can integrate the local power system, being generated in another location (remotely). The demand curve above reaches a maximum of 190 W, in power terms, and equals a monthly average consumption of 78 kWh/month.
In addition to the load curve, the profile of solar irradiance for a clear sky day is represented in the same graph, as well as the operational temperature of the photovoltaic panel. The mismatch between the variables is very noticeable, and this is crucial in the results and consequences coming from simulation results. It is common that residential CUs do not support photovoltaic system, especially because there is no area for its installation (area occupation factor), but thanks to the power sharing mode, a photovoltaic system can integrate the local power system, being generated in another location (remotely). The demand curve above reaches a maximum of 190 W, in power terms, and equals a monthly average consumption of 78 kWh/month. In addition to the load curve, the profile of solar irradiance for a clear sky day is represented in the same graph, as well as the operational temperature of the photovoltaic panel. The mismatch between the variables is very noticeable, and this is crucial in the results and consequences coming from simulation results.

Results
In the distribution feeder, it was delimited an interest zone. Figure 9 shows the gradual insertion of PV generators in the CUs. Bar 09 is located downstream from the substation (last bar), configuring a critical spot from an operational point of view. In this bar, six one-phase residential consumers are connected and distributed in the phase sequence already shown in Table 6. It is common that residential CUs do not support photovoltaic system, especially because there is no area for its installation (area occupation factor), but thanks to the power sharing mode, a photovoltaic system can integrate the local power system, being generated in another location (remotely). The demand curve above reaches a maximum of 190 W, in power terms, and equals a monthly average consumption of 78 kWh/month.
In addition to the load curve, the profile of solar irradiance for a clear sky day is represented in the same graph, as well as the operational temperature of the photovoltaic panel. The mismatch between the variables is very noticeable, and this is crucial in the results and consequences coming from simulation results.

Results
In the distribution feeder, it was delimited an interest zone. Figure 9 shows the gradual insertion of PV generators in the CUs. Bar 09 is located downstream from the substation (last bar), configuring a critical spot from an operational point of view. In this bar, six one-phase residential consumers are connected and distributed in the phase sequence already shown in Table 6. In the simulation, the PV generation penetration was customized up to the maximum limit allowed in one-phase configuration, according to Table 2, and the DG was gradually inserted into each CU. The results were divided in two simulation cases, assuming the conditions imposed on [14] limiting the installed power for one-phase systems to a maximum 15 kW. The cases were partitioned into: 1. Case 1: P = 5 kW, 33% of maximum limit allowed. 2. Case 2: P = 10 kW, 66% of maximum limit allowed.

Simulation Results
Cases 1 and 2 were divided into three subcases that define the DG penetration criteria: 1. Subcase 1: Only one CU with PV; 2. Subcase 2: Three CU`s with PV; 3. Subcase 3: All the CU`s with PV. In the simulation, the PV generation penetration was customized up to the maximum limit allowed in one-phase configuration, according to Table 2, and the DG was gradually inserted into each CU. The results were divided in two simulation cases, assuming the conditions imposed on [14] limiting the installed power for one-phase systems to a maximum 15 kW. The cases were partitioned into: 1.

Simulation Results
Cases 1 and 2 were divided into three subcases that define the DG penetration criteria: 1. Subcase 1: Only one CU with PV; 2.
Subcase 3: All the CU's with PV.
Case 1 simulation considered the gradual insertion of DG, keeping constant irradiance and demand conditions (according to Figure 8). The power installed for the photovoltaic generation units assumed in the subcases was fixed initially at 5.2 kWp. According to the three subtopics mentioned above, the levels of penetration for the subcases were 16%, 50%, and 100%, respectively. Subcase 3 corresponded to 41.6% of total power from the substation serving this extension, as well as subcase 1, that contemplated only 7% of total power from the substation. Table 7 below presents synthetically the methodological organization for simulation:  Cases 1 and 2 considered the installed power variations of photovoltaic generation; the subcases of photovoltaic penetration are incremented by the topics described in the first paragraph of this article and the power values of 5 kW and 10 kW correspond to cases 1 and 2, respectively.
The graphs in Figure 10 present the dynamic of voltage in the CUs connected to bar 09, for all subcases. Both graphs are limited by a red line expressing the maximum limit of electric voltage established [15].
The graphs in Figure 10 present the dynamic of voltage in the CUs connected to bar 09, for all subcases. Both graphs are limited by a red line expressing the maximum limit of electric voltage established [15]. On the graphs for subcases 1 and 2, it is noticed that the voltage levels are below the allowable limit. In subcase 1, only one PV generator connected to the electrical system was considered, corresponding to 34% of the allowed on one-phase topology. In this first moment, it is possible to see that CU 78 also suffers a minimal voltage increase, the effect arising from electricity generation from PVs connected to CU 75. The CUs interconnected to phases 2 and 3 tend to be more constant throughout the day, presenting a relative average dispersion for phase 1 of approximately 3%. Amplifying the number of PV generators to the electrical system, keeping constant the potency of the PVs (5 kW), a more symmetrical distribution of the voltage profile is noticed, varying only the amplitude of the voltage module, when comparing the CUs with and without PVs (Figure 10c,d). On the graphs for subcases 1 and 2, it is noticed that the voltage levels are below the allowable limit. In subcase 1, only one PV generator connected to the electrical system was considered, corresponding to 34% of the allowed on one-phase topology. In this first moment, it is possible to see that CU 78 also suffers a minimal voltage increase, the effect arising from electricity generation from PVs connected to CU 75. The CUs interconnected to phases 2 and 3 tend to be more constant throughout the day, presenting a relative average dispersion for phase 1 of approximately 3%. Amplifying the number of PV generators to the electrical system, keeping constant the potency of the PVs (5 kW), a more symmetrical distribution of the voltage profile is noticed, varying only the amplitude of the voltage module, when comparing the CUs with and without PVs (Figure 10c,d).
In subcase 2, electrical voltage remained within the allowed limits and the CU 76 and 77 obtained overvoltage in relation to the other CUs, reaching 5% over CU 75, approaching critical value. This fact is due to the conductor employed in the service extension (Table 6), the variable related to the cross section and the length, which directly interfere with the electrical parameters. The last stage of case 1 analysis (subcase 3) takes into consideration all CUs equipped with PV generator, contemplating 100% penetration at bar 09. Observing Figure 10e, the reach of voltage maximum limit for CU 78 is confirmed, without surpassing it. According to regional utility regulation [15], the service extension should be scaled to meet the maximum expected demand, and such fact is partially counterpointed in subcase 3 (Figure 10e). Still, in evidence for subcase 3, the three-phase voltage in bar 09 remains balanced and tolerated within the demanded values (shown in Figure 10f).
Simulation for case 2 also refers to items presented in Table 7, considering the installed PV potency of 10 kW, corresponding to 66% of maximum limit allowed on the network in one-phase configuration.
With all simulations performed for case 2, Even with the increase of the installed power of the DG, the voltage at the connection point remains acceptable, as seen in Figure 11a, other CUs without DG and connected to the same phase sequence are not influenced by overvoltage, the voltage profile may be seen in Figure 11b, the violation of the critical limit for CUs 76 and 77 can be seen in Figure 11c, although such variation is tiny and still ensures feeder security and operability. Voltage variation between phases of CUs in Figure 11c is on average between 1.7% and 1.9%. Still for subcase 2, CUs without DG remained with the voltage values within the required (Figure 11d) The situation is more critical in subcase 3 (Figure 11e), where practically all consumer units surpass voltage critical limit and the CU 78 reaches maximum voltage value, breaking the protection parameter against overvoltage of the inverter. CU 78 surpasses 5% of the critical voltage value, corresponding to 11% of nominal value. Below, the voltage relative variation percentage for critical and nominal limit, 233 and 220 volts, respectively, are shown.
Referring to subcase 3, where all CU's are equipped with distributed generation of the same rated power installed and under the same demand length, the three-phase voltages seen from the mains side remained within the values required by normative resolution, as illustrated by graph of Figure 11f.
It is possible to verify the level of variability of the supply voltage of the CUs (Figure 12) around the critical limit (232 V) and the nominal value (220V), respectively shown by Figure 12a,b. The consumer unit 78 presented the highest amplitude in voltage increase, such effect was motived by the impedance of the service line and the installed power of DG.
Changes in penetration level of photovoltaic generator present large variations in the electrical system phases in CUs that are connected in the same phases, even with all residencies under the same demand and solar irradiation conditions, it is verified that unbalance level (∆U) occurs between CUs 75 and 78 with ∆U of 4%. The Table 8 below quantifies significantly the voltage variations between the CUs connected to the same phases and to distinct phases. Subcase 3 presents higher rate of variation between consumers connected to phase V A . Fitted into the flowchart (Figure 7), the simulation will return to step 2 through the verifying routine (step 5). As a partially mitigating solution, the PVs were partitioned into three-phase configuration, where the power injected in the network by phase was reduced by 30% of PV nominal power. From these considerations, the simulation of load flow was particularly remade for subcase 3 and can be seen graphically below. Referring to subcase 3, where all CU's are equipped with distributed generation of the same rated power installed and under the same demand length, the three-phase voltages seen from the mains side remained within the values required by normative resolution, as illustrated by graph of Figure 11f.
It is possible to verify the level of variability of the supply voltage of the CUs (Figure 12) around the critical limit (232 V) and the nominal value (220V), respectively shown by Figures 12a and 12b. The consumer unit 78 presented the highest amplitude in voltage increase, such effect was motived by the impedance of the service line and the installed power of DG. The mitigating solution adopted to reduce voltages in the CUs was the redefinition of the PVs connection, distributing all the power from the photovoltaic system simultaneously in three phases. The effect is verified in Figure 13; despite remaining above critical limit allowed, maximum overvoltage located in CU 78 were reduced in 4%. Table 9 below compares the phase asymmetry for both cases, showing the subcase 3 belonging to case studies 1 and 2, as well as the results of the mitigating effect. Cheol Jeong [20] related in his work a voltage control strategy that uses the step voltage regulator (SVR), one of the found results that after using the SVR it was possible to release a capacity of 2.7 MW more, maintaining the voltage levels specified. Shin [21] investigated the delivery voltage of the Changes in penetration level of photovoltaic generator present large variations in the electrical system phases in CUs that are connected in the same phases, even with all residencies under the same demand and solar irradiation conditions, it is verified that unbalance level (∆U) occurs between CUs 75 and 78 with ∆U of 4%. The Table 8 below quantifies significantly the voltage variations between the CUs connected to the same phases and to distinct phases. Subcase 3 presents higher rate of variation between consumers connected to phase V . Fitted into the flowchart (Figure 7), the simulation will return to step 2 through the verifying routine (step 5). As a partially mitigating solution, the PVs were partitioned into three-phase configuration, where the power injected in the network by phase was reduced by 30% of PV nominal power. From these considerations, the simulation of load flow was particularly remade for subcase 3 and can be seen graphically below.
The mitigating solution adopted to reduce voltages in the CUs was the redefinition of the PVs connection, distributing all the power from the photovoltaic system simultaneously in three phases. The effect is verified in Figure 13; despite remaining above critical limit allowed, maximum overvoltage located in CU 78 were reduced in 4%. Table 9 below compares the phase asymmetry for both cases, showing the subcase 3 belonging to case studies 1 and 2, as well as the results of the mitigating effect. Cheol Jeong [20] related in his work a voltage control strategy that uses the step voltage regulator (SVR), one of the found results that after using the SVR it was possible to release a capacity of 2.7 MW more, maintaining the voltage levels specified. Shin [21] investigated the delivery voltage of the on-load tap-changers OLTC such that a large-capacity DG would not create unsafe voltage variations when penetrating the distribution network.   It is noticeable from the verification of the reduction of the voltage variation that when the PVs are redefined, the symmetrical distribution of the phases significantly improves the voltage quality in the three-phase circuit, allowing equilibrium. From a technical losses point of view, the impact in electrical losses reduction was noticed in the conductors belonging to the circuit of the power grid section and service extensions to consumer units. From the feeder's operational profile, it is simple to conclude that electrical losses will be relatively low, although relevantly aggravated when there is too much increase in photovoltaic systems connected to the power grid. Figure 14 presents a comparison of technical losses evolution in the studied feeder, were Figure 14a considers the absence of distributed generation and Figure 14b considers the case of massive insertion of photovoltaic generators connected to a single phase Table 9. Simulation subcases for case 2.

Maximum Voltage Drop
Phases-Volts (V)

Phase Variation (∆U) Distinct Phases (%)
Phase Variation (∆U) Same Phases (%)    It is noticeable from the verification of the reduction of the voltage variation that when the PVs are redefined, the symmetrical distribution of the phases significantly improves the voltage quality in the three-phase circuit, allowing equilibrium. From a technical losses point of view, the impact in electrical losses reduction was noticed in the conductors belonging to the circuit of the power grid section and service extensions to consumer units. From the feeder`s operational profile, it is simple to conclude that electrical losses will be relatively low, although relevantly aggravated when there is too much increase in photovoltaic systems connected to the power grid. Figure 14  With the results shown in the graphs of Figure 14, the aggressive impact of distributed generation in electrical losses of the line each hour is evident. Please note that, at this moment, maximum penetration of photovoltaic systems in the bar under analysis is considered, all with   It is noticeable from the verification of the reduction of the voltage variation that when the PVs are redefined, the symmetrical distribution of the phases significantly improves the voltage quality in the three-phase circuit, allowing equilibrium. From a technical losses point of view, the impact in electrical losses reduction was noticed in the conductors belonging to the circuit of the power grid section and service extensions to consumer units. From the feeder`s operational profile, it is simple to conclude that electrical losses will be relatively low, although relevantly aggravated when there is too much increase in photovoltaic systems connected to the power grid. Figure 14  With the results shown in the graphs of Figure 14, the aggressive impact of distributed generation in electrical losses of the line each hour is evident. Please note that, at this moment, maximum penetration of photovoltaic systems in the bar under analysis is considered, all with With the results shown in the graphs of Figure 14, the aggressive impact of distributed generation in electrical losses of the line each hour is evident. Please note that, at this moment, maximum penetration of photovoltaic systems in the bar under analysis is considered, all with nominal power of 10kWp. Comparing daily losses in the feeder section (section 08) for both cases, in the absence of DG, the losses in section 08 reached 5.11 Wh/day, and in the face of massive insertion, 6244 kWh/day was obtained, around 1200 times higher than normal condition.
The impact of technical losses is also seen when the load flow study is reconsidered in the simulation based on mitigation strategy adopted in this article, based on symmetrical distribution of nominal power of photovoltaic generator by phase. Figure 15 graphically illustrates the comparison between technical losses in the line for high concentration of photovoltaic systems connected to only one phase, and photovoltaic systems symmetrically distributed in the three-phase system.
Taking as reference the network section of the feeder (section 08), it is evident from the graphs above that the losses were severely reduced, around 55.78%. The mitigating strategy divides the installed power from a photovoltaic generator of 10 kWp located in a CU symmetrically in the three phases. For practical purposes, this is possible with the employment of a three-phase inverter or three one-phase micro-inverters, consequently relieving the conductors loading. In the service extensions to the CUs, there was also a reduction of losses, however the extension physical characteristics from each consumer unit vary according to the type of conductor, distance, and installation method, implying relatively different loss values. Losses reduction in the service extensions decreased 89.16%, obtained through an average of the extensions. The following graph (Figure 16) addresses the relative reduction of technical losses by extension against the massive insertion of photovoltaic generators, but strategically based on the mitigating solution.
nominal power of 10kWp. Comparing daily losses in the feeder section (section 08) for both cases, in the absence of DG, the losses in section 08 reached 5.11 Wh/day, and in the face of massive insertion, 6244 kWh/day was obtained, around 1200 times higher than normal condition.
The impact of technical losses is also seen when the load flow study is reconsidered in the simulation based on mitigation strategy adopted in this article, based on symmetrical distribution of nominal power of photovoltaic generator by phase. Figure 15 graphically illustrates the comparison between technical losses in the line for high concentration of photovoltaic systems connected to only one phase, and photovoltaic systems symmetrically distributed in the three-phase system. Taking as reference the network section of the feeder (section 08), it is evident from the graphs above that the losses were severely reduced, around 55.78%. The mitigating strategy divides the installed power from a photovoltaic generator of 10 kWp located in a CU symmetrically in the three phases. For practical purposes, this is possible with the employment of a three-phase inverter or three one-phase micro-inverters, consequently relieving the conductors loading. In the service extensions to the CUs, there was also a reduction of losses, however the extension physical characteristics from each consumer unit vary according to the type of conductor, distance, and installation method, implying relatively different loss values. Losses reduction in the service extensions decreased 89.16%, obtained through an average of the extensions. The following graph ( Figure 16) addresses the relative reduction of technical losses by extension against the massive insertion of photovoltaic generators, but strategically based on the mitigating solution. Partitioning a photovoltaic generator unit from electrical losses reduction point of view is quite livable. It is possible to guarantee that the service extensions cooperate under normalized loading  Taking as reference the network section of the feeder (section 08), it is evident from the graphs above that the losses were severely reduced, around 55.78%. The mitigating strategy divides the installed power from a photovoltaic generator of 10 kWp located in a CU symmetrically in the three phases. For practical purposes, this is possible with the employment of a three-phase inverter or three one-phase micro-inverters, consequently relieving the conductors loading. In the service extensions to the CUs, there was also a reduction of losses, however the extension physical characteristics from each consumer unit vary according to the type of conductor, distance, and installation method, implying relatively different loss values. Losses reduction in the service extensions decreased 89.16%, obtained through an average of the extensions. The following graph ( Figure 16) addresses the relative reduction of technical losses by extension against the massive insertion of photovoltaic generators, but strategically based on the mitigating solution. Partitioning a photovoltaic generator unit from electrical losses reduction point of view is quite livable. It is possible to guarantee that the service extensions cooperate under normalized loading Partitioning a photovoltaic generator unit from electrical losses reduction point of view is quite livable. It is possible to guarantee that the service extensions cooperate under normalized loading conditions, optimizing the load flow. On the other hand, the use of this configuration may not be as economically attractive, for conditioning the use of more devices. A solution to minimize losses can be seen in Gkaidatzis et al. [22] where the size and allocation of DG units was presented in order to minimize losses.

Experimental Data Result
An experimental array (CU) was developed in this case study to perform the comparison between the simulated and measured data. In the real context, the CU 78 was entirely instrumented storing and generating data used as input parameters in the simulations.
Several levels of massive insertion of systems in bar 09 were considered in the simulation. Subcase 1 analyzed initially was extracted to compare the dynamic of electrical voltage with real data extracted from CU 78. As shown in Figure 17a, the simulated voltage presents a similar behavior to experimental voltage.
The comparison performed between experimental and simulated voltages presented maximum relative error of ±2% (Figure 17b). The real voltage profile obtained from SCADA proves stability below maximum critical limit defined, corroborating with the simulation, and it is noticed that the interaction between the service extension from a CU and reverse power flow will result directly in the voltage behavior. The simulation performed initially in case 1, subcase 1 (Figure 10a) was repeated to verify the interference of the extension in the presence of reverse load flow. Therefore, a new PV was defined, this time connected to another CU (CU 78) and the simulation remade, keeping constant all premises from case 1, comparing with the previous result. The results can be seen in the graphs below. be seen in Gkaidatzis et al. [22] where the size and allocation of DG units was presented in order to minimize losses.

Experimental Data Result
An experimental array (CU) was developed in this case study to perform the comparison between the simulated and measured data. In the real context, the CU 78 was entirely instrumented storing and generating data used as input parameters in the simulations.
Several levels of massive insertion of systems in bar 09 were considered in the simulation. Subcase 1 analyzed initially was extracted to compare the dynamic of electrical voltage with real data extracted from CU 78. As shown in Figure 17a, the simulated voltage presents a similar behavior to experimental voltage.  The comparison performed between experimental and simulated voltages presented maximum relative error of ±2% (Figure 17b). The real voltage profile obtained from SCADA proves stability below maximum critical limit defined, corroborating with the simulation, and it is noticed that the interaction between the service extension from a CU and reverse power flow will result directly in the voltage behavior. The simulation performed initially in case 1, subcase 1 (Figure 10a) was repeated to verify the interference of the extension in the presence of reverse load flow. Therefore, a new PV was defined, this time connected to another CU (CU 78) and the simulation remade, keeping constant all premises from case 1, comparing with the previous result. The results can be seen in the graphs below. Figure 18 portrays the influence of physical aspects and type of conductor employed in the service extension in the divergence that occurs in voltage values (Figure 18a). The larger amplitudes occur when the circuit is reverse loaded, in the period between 09:00 and 16:00, with the effect caused by the mismatch between the demand curve and electric generation. Figure 18b proves the constant behavior of the reverse power flow for both simulations and CUs, where the overvoltage in the simulation with CU 78 is 2% higher regarding the simulation in subcase 1, based on CU 75.  Figure 18 portrays the influence of physical aspects and type of conductor employed in the service extension in the divergence that occurs in voltage values (Figure 18a). The larger amplitudes occur when the circuit is reverse loaded, in the period between 09:00 and 16:00, with the effect caused by the mismatch between the demand curve and electric generation. Figure 18b proves the constant behavior of the reverse power flow for both simulations and CUs, where the overvoltage in the simulation with CU 78 is 2% higher regarding the simulation in subcase 1, based on CU 75.

Conclusions
This article presents a load flow simulation in a secondary electrical power distribution feeder based on OPENDSS open-source tool. An analysis of the electrical voltage variable is performed, seen from the consumer when it is equipped with a photovoltaic electricity generation system. The electrical modeling of a power feeder was performed, as well as the allocation of the loads and the PV generators, matching the simulation scenario with the actual context. The simulation results proved that, in the occurrence of massive penetration of photovoltaic systems in one-phase connections with maximum power installed of 10 kWp, there was an overvoltage effect in the consumer units, countering the local regulations, which allows network connections of photovoltaic systems with one-phase configuration with up to 15 kWp. The solution for the overvoltage problem was mitigated by changing the connection type of PV systems to three-phase, ensuring voltage regulation within limits pre-established by normative resolution; it would also be possible to repower the service extension and keep the configuration. In the simulation considering partial insertion of

Conclusions
This article presents a load flow simulation in a secondary electrical power distribution feeder based on OPENDSS open-source tool. An analysis of the electrical voltage variable is performed, seen from the consumer when it is equipped with a photovoltaic electricity generation system. The electrical modeling of a power feeder was performed, as well as the allocation of the loads and the PV generators, matching the simulation scenario with the actual context. The simulation results proved that, in the occurrence of massive penetration of photovoltaic systems in one-phase connections with maximum power installed of 10 kWp, there was an overvoltage effect in the consumer units, countering the local regulations, which allows network connections of photovoltaic systems with one-phase configuration with up to 15 kWp. The solution for the overvoltage problem was mitigated by changing the connection type of PV systems to three-phase, ensuring voltage regulation within limits pre-established by normative resolution; it would also be possible to repower the service extension and keep the configuration. In the simulation considering partial insertion of PVs, voltage-raising effect spreads in CUs without DG but are connected to the same bar knot. Overvoltage results obtained in the experiment come from a combination of characteristics of service extension, local demand profile, and photovoltaic generator global characteristics, and it is possible to identify that the voltage tends to reach operating critical limit. Theoretical and experimental comparison is based on subcase 1 (considering only one DG insertion with 5 kWp of nominal power connected to the bar of the distribution system), where it is possible to verify a bias between voltage curves, ensuring the fidelity of the result. In general, this work proposes an analysis methodology based on a free software in the regulation restrictions to allow the access of new users to the low-voltage network in a way to ensure stable and incorrupt operation of the electrical system.