Strategic Sizing and Placement of Distributed Generation in Radial Distributed Networks Using Multiobjective PSO

Distributed generators (DGs) o ﬀ er signi ﬁ cant advantages to electric power systems, including improved system losses, stability, and reduced losses. However, realizing these bene ﬁ ts necessitates optimal DG site selection and sizing. This study proposes a traditional multiobjective particle swarm optimization (PSO) approach to determine the optimal location and size of renewable energy-based DGs (wind and solar) on the Namibian distribution system. The aim is to enhance voltage pro ﬁ les and minimize power losses and total DG cost. Probabilistic models are employed to account for the random nature of wind speeds and solar irradiances. This is used in an algorithm which eventually optimizes the siting and sizing of DGs using the nearest main substation as reference. The proposed method is tested on the Vhungu-Vhungu 11 kV distribution network in Namibia. Four cases were considered: base case with no DG, solar power, wind power, and a hybrid of both wind and solar. Optimal values for each case are determined and analyzed: 0.69.93 kW at 26 km for solar PV-based DG and 100 kW at 42 km for wind-based DG. These ﬁ ndings will serve as a valuable blueprint for future DG connections on the Namibian distribution network, providing guidance for optimizing system performance.


Introduction
Namibia faces a significant electricity demand, with 611 MW needed to meet its requirements.To fulfill this demand, the country heavily relies on imports, obtaining 60% of its electricity from South Africa and Zambia [1].However, Namibia possesses a unique advantage due to its geographic location, as it boasts one of the highest solar irradiation levels in the world, reaching 3000 kWh/m 2 and an estimated maximum solar potential of 344 GW.Additionally, the country's 1600 km coastline experiences a yearly average wind speed of 10 m/s, offering a substantial wind energy potential of 27.201 GW.Furthermore, Namibia's terrain is covered with invader bush, which, if effectively utilized, has the potential to generate up to approximately 100 MW of power [2].As Namibia looks to the future, there is an expectation of integrating numerous distributed generators (DGs) into its distribution networks.The integration of these renewable energy sources into medium voltage networks brings forth numerous benefits.However, it is crucial to emphasize that these benefits can only be fully realized if the size and location of these DGs are carefully determined and optimized.The optimal sizing and siting of DGs serve as the foundation for effective renewable energy planning in the country.The placement of DGs holds particular significance when considering the impact on power losses within the system.Addressing this issue requires careful consideration of the DGs' size and location, taking into account system losses, voltage stability, and overall network performance.Namibia's abundant solar irradiation, coastal wind resources, and potential utilization of invader bush present promising opportunities for renewable energy integration.However, optimizing the size and placement of DGs is paramount to ensure the realization of these benefits while mitigating potential challenges associated with power losses and voltage instability.By implementing effective planning strategies, Namibia can harness its renewable energy potential and pave the way for a sustainable and resilient energy future.
In residential distributed generation (DG) systems, the presence of voltage instability poses a potential challenge.This is primarily due to the limitations of low voltage lines, typically operating at around 11 kV, which are unable to effectively accommodate the penetration of over 500 kW of generated electricity [3].The voltage level at each bus, where both the DG units and loads are connected, further influences the extent of power losses experienced within the system [4].The adverse impacts of DG integration on power losses have been extensively studied and documented [5].To address the challenges associated with DG integration, researchers have solved the optimization problem of determining the optimal size and location of DG units within an existing distribution network [6].By utilizing advanced optimization techniques such as multiobjective particle swarm optimization (PSO), researchers strive to identify the optimal configuration of distributed generation (DG) units that simultaneously minimize power losses and consider additional objectives such as voltage stability and cost.Extensive research and empirical studies have consistently demonstrated that a considerable portion of DGs, typically ranging from 10% to 50%, can be effectively integrated into distribution networks without compromising their overall performance [7].However, achieving optimal integration and minimizing the impact on power losses require careful planning, taking into account factors such as network topology, load demand patterns, and available capacity.
Sizing and siting of distributed generators (DGs) in a radial distribution network pose a complex optimization problem.Single or a multiobjective approach is used.Single objective functions are aimed at minimizing system power losses, improving system voltage profiles, minimizing costs, enhancing system reliability, and more.On the other hand, multiobjective functions combine two or more single objective functions to capture multiple criteria simultaneously.Many articles in literature have worked on optimization techniques for DG integration; their main contribution and research gaps are individually analyzed in Table 1.Some literature focus on technical aspects [8][9][10][11][12][13][14], economic impacts [15][16][17], and environmental considerations [18].Some literature considers a combination of both [16,19,20].Various methods exist in the literature for solving the DG siting and sizing problem, each with its own advantages and disadvantages.Among them, particle swarm optimization (PSO) is widely utilized due to its simplicity, ease of implementation, and ability to explore a broad search space without being constrained by local optima.The weighted sum approach is commonly used in practical applications due to its simplicity and versatility in incorporating various constraints beyond the feasible region.This approach is indicated in [10,11,14,20].This approach involves assigning relative importance to different objectives based on their significance in the overall system performance.For instance, minimizing power losses or improving voltage profiles may be assigned higher priority compared to other objectives.Furthermore, it is crucial to consider the uncertainty associated with solar and wind resources [9,10,12,13,21,22].As solar and wind power generations exhibit intermittent characteristics, their output is subject to variability and uncertainty.Probabilistic or stochastic methods can be employed to model and analyze this uncertainty, enabling robust decisions in DG siting and sizing.
From the existing literature, it is evident that the majority of the researchers have focused on singular components of the objective functions, such as technical, economic, or environmental factors.However, this study takes a comprehensive approach by considering both technoeconomic aspects when addressing the optimal allocation of distributed generator (OADG) problem.One common limitation found in previous research is the assumption of dispatchable DGs, neglecting the stochastic nature of renewable energybased DGs, specifically solar and wind.This assumption can lead to undesirable outcomes, as it fails to account for the intermittent and unpredictable nature of these energy sources.To overcome this limitation, this study acknowledges and incorporates the stochastic behavior of solar and wind resources into the analysis.Moreover, while some studies consider the random nature of renewable energybased DGs, they often model it on a monthly or yearly basis [10,13,15].This approach can potentially impact the accuracy of the results, as it overlooks the time-varying and hourly probabilistic characteristics of wind and solar resources.In contrast, this paper introduces an hourly probabilistic model that effectively captures the random nature of these energy sources, resulting in a more realistic representation of their behavior.The main contributions of this study can be summarized as follows: (i) Development of an hourly probabilistic model: by incorporating an hourly probabilistic model, the study accurately captures the random characteristics of wind and solar resources, enabling a more precise analysis of their impact on the DG allocation problem (ii) Consideration of technoeconomic parameters: in addition to system losses, this study takes into account important technoeconomic parameters such as DG investment cost and voltage profiles.By considering these factors, the proposed approach ensures that the DG sizing and siting decisions are not only technically sound but also economically viable (iii) Introduction of a distance function: to determine the optimal distance of the DG from the slack, a distance function is introduced.This function helps in identifying the most favorable location for the DGs, considering factors such as power flow and voltage stability (iv) Simultaneous siting of wind and solar DGs: unlike previous studies that often focus on either wind or solar DGs, this study considers the simultaneous siting of both types of DGs.This comprehensive approach allows for a more holistic and integrated analysis of renewable energy integration in the distribution system (v) This study employed a traditional weighted multiobjective PSO approach.It utilizes relative weights to optimize multiple objectives effectively

Contribution
Ref.
(i) The study is aimed at reducing feeder losses and enhancing voltage quality.A Stamford generator, with a rated power of 1350 kW, was integrated into the network (ii) Research gaps: the analysis did not consider time-dependent loads and time-dependent generation.Additionally, renewable energy-based distributed generators (DGs) were not considered [8] (i) The modified voltage index (MVI) method was utilized to determine the optimal placement and size of DG units in order to enhance the voltage stability margin (ii) Research gaps: the study adopted a single DG placement approach, considering one DG at a time and assessing its placement at different buses.Economic variables such as investment and operational costs were not taken into consideration in the analysis [9] (i) The approach is aimed at minimizing power losses in the electrical network while improving voltage stability and network security.To address the stochastic nature of solar irradiance and wind speed, suitable probabilistic models were employed, allowing for realistic representation of their variability in the analysis (ii) Research gaps: the study focused solely on the technical aspects of optimal placement and sizing of solar and wind distributed generators (DGs) in the distribution territory; however, economic factors such as investment costs, operational costs, or financial considerations were not taken into account in the analysis [10] (i) Used weighted multiobjective voltage index to minimize real power loss and enhance the voltage profile within the system (ii) Research gaps: the study considers a general type of distributed generators (DGs) integrated into the network without specifying a particular DG type.The focus of the analysis is solely on technical aspects [11] (i) It employs a multistate modeling approach to account for the uncertain nature of wind and solar resources.The proposed model evaluates deviations in several key parameters, including annual energy losses (AEL), total DG penetration, loss of load expectation (LOLE), and loss of energy expectation (ii) Research gaps: the study primarily focuses on technical aspects and incorporates monthly, seasonal, and yearly models of solar and wind resources.It does not include hourly models, which means that the analysis does not account for the fine-grained variations and fluctuations in solar irradiance and wind speed throughout the day [12] (i) A methodology was proposed to optimize the allocation of different types of renewable distributed generation (DG) units within the distribution system with the objective of minimizing annual energy loss.This methodology involves the utilization of a multistate model for the hourly modeling of renewable energy sources (ii) Research gaps: the proposed methodology only focused solely on the technical aspects of optimally allocating different types of renewable distributed generation (DG) units within the distribution system [13] (i) Crow search algorithm (CSA) was used to determine the optimal size and allocation of distributed generators (DGs).A multiobjective function was formulated to address the objectives of reducing active power losses and improving the voltage profile (ii) Research gaps: the study solely focused on technical aspect only, and the algorithm was applied to IEEE 33-bus without incorporating real data [14] (i) Used time-varying and seasonal optimal placement and sizing of both intermittent renewable energy sources (such as wind energy) and nonintermittent renewable energy sources (such as solar energy).To account for the multistate and hourly probabilistic nature of wind speed and solar irradiance data, the study employed appropriate modelling techniques (ii) Research gaps: does not consider the economic aspects related to cost-benefit analysis or financial feasibility [21] (i) Uses multiobjective management approach that combines network reconfiguration with the allocation and sizing of renewable distributed generations (DGs) with the aim of minimizing active power loss, annual operation costs (including installation, maintenance and active power loss costs), and pollutant gas emissions.The optimization problem is solved by considering the time sequence variance in renewable DGs and load (ii) Research gaps: used simulation models (IEEE 33 bus) without validating the results against real-world data [15] (i) Incorporated both technical and economic aspects, with aim of minimizing power losses and maximizing profit.Additionally, it accounted for the stochastic nature of wind and solar resources (ii) Research gaps: the DG placement study lacked the use of real data and overlooked the significance of considering voltage profiles, which are crucial in DG placement studies [16].
(i) An optimization model was developed to address the allocation of distributed generators (DGs) with the primary objective of minimizing the total planning cost (ii) Research gaps: the optimization model presented in the study focused solely on wind energy; integration challenges with other sources such as wind and solar were not tested and validated [17] 3 Journal of Energy (vi) The study validates the proposed algorithm using real data of wind and solar applied to a practical distribution network as opposed to simulation models mostly used by researchers

Probabilistic Modelling of Renewable Energy Sources
2.1.Modeling of Wind Turbine Generator.Wind power exhibits a stochastic nature as it is influenced by the random and variable wind speeds.Modeling of such system at a particular hour t is achieved using the Weibull probability distribution function as given in [16] where f v v t is the Weibull PDF, C t is the scale parameter, and k t is the shape parameter at time period t hour.The Weibull PDF parameters at a particular time t hour can be estimated by using the following equations [16]: where σ v t and μ v t represent the standard deviation and the mean wind speeds at specific time t hour.Г is a gamma function at time t hour.

Contribution
Ref.
(i) Multiobjective algorithm was applied to reduce power loss, maximizing the voltage stability index, minimizing voltage deviation, lowering real power loss costs, increasing real power loss savings, and reducing CO 2 emissions (ii) Research gaps: did not incorporate the potential effects of the intermittent nature of renewable distributed generations (DGs).
The intermittent nature of renewable DGs, such as wind and solar, can introduce uncertainty and variability into the power system [18] (i) Three optimization techniques, PSO, variable constraint PSO (VCPSO), and GA algorithms, are applied to find the optimal size and placement of multiple DGs integrated into electrical power network.VCPSO was offered an improved solution for the optimal placement and size of DGs in terms of the accuracy of the global optimality (ii) Draw back: did not address the application of the hybrid PSO algorithm in real-world distribution networks [19] (i) Cost-based analysis was used on distributed generators (DGs), to determine installation costs, operational costs, and maintenance costs.The objective of the analysis was to minimize losses and maximize the loading capability of the system while ensuring that voltage stability is not compromised (ii) Research gaps: assumed a constant factor of 0.95 as power factor for DG operation which is not the case in real-life situations; the power factor of DGs may vary depending on various factors such as load conditions, system requirements, and control strategies [20] (i) Crisscross optimization algorithm and Monte Carlo simulation method (CSO MCS), used to address the optimal distributed generation allocation (ODGA) problem.This method considers the uncertainties associated with wind, solar, and load consumption (ii) Research gaps: does not consider an hourly resolution.This can lead to inaccurate results, as it fails to capture the time-varying nature of renewable energy generation [22] (i) PSO is used on IEEE 33 radial distribution system with different types of voltage-dependent load models.The results reveal that combination of active-reactive power DG is giving better results for power loss reduction and voltage profile improvement (ii) Research gaps: the study did not investigate the integration of renewable energy sources or consider uncertainties related to renewable generation [23] (i) Provided a review of optimization techniques used for DG sizing and placement in a distributed network (ii) Research gaps: did not consider different scenarios and constraints to identify the most appropriate techniques for specific applications [24] (i) MOPSO algorithm has been used to find the optimal solution of DG sizing and locating problem; this was tested on IEEE 33-bus reliability enhancement of the grid which was confirmed (ii) Research gaps: did not consider implementing renewable DGs with uncertain output power, such as PV panels or wind turbines [25] (i) Used multiobjective bat algorithm on IEEE 69-bus; from the obtained results, it is observed that the best localization and sizing of DG unit give more flexibility to the network (ii) Research gaps: did not cost function and the algorithm did not use real data [26] (i) A backtracking search optimization algorithm (BSOA) is developed to enhance voltage profile and reduce real network losses (ii) Research gaps: the study assumed that the output of renewable energy (RE) sources is dispatchable, hence did not consider intermittent nature of renewable energy sources [27] 4 Journal of Energy 2.2.Modeling of PV Generator.Solar irradiances vary throughout the day.Beta probability function is used to describe the uncertain nature of the solar resources [16].Beta PDF at a time period t hour is defined by where α t and β t are the beta parameters.These parameters can be approximated using the following equations [14]: where μ s and α s are the mean and variance of the solar irradiance.

Modeling of Electrical Loads.
Modelling of loads in an electrical network is classified as follows: (i) Constant power loads: there is no relationship between the voltage and the power drawn at each bus (ii) Constant current loads: the voltage varies proportional with the active and reactive powers (iii) Constant impedance loads: the active and reactive powers are directly proportional to the square of the voltage To investigate the effect of changes in load demand over the year, this paper proposes modelling of loads as constant time-varying loads.This model considers the potential increase in loads that may occur during the planning horizon.The active and reactive power demands in year N y are represented by the following [23]: where P L,i N y and Q L,i N y are the active and reactive power loads at bus i after N y years, respectively.P L,i 0 and Q L,i 0 are the active and reactive power loads at bus i for the base year, respectively.r is the annual growth rate in the load.This value is calculated based on the historical data on the load profiles, and N y is the number of years in the planning period.

Power Generation Models
2.4.1.Power Generation by PV Array.PV module power is a function of the temperature (T) and solar irradiance s and is determined by the following equations [12]: where P PV is the power produced by the PV module (W), FF is the fill factor of the module, V is the module voltage output (V), I is the module current output (A), T c is the cell temperature ( °C), T a is the ambient temperature ( °C), T NOCT is the normal operating temperature of the module, T c,ref is the reference temperature (25 °C), s ref is the reference insolation (1000 W/m 2 ), I sc is the short circuit current of the module, V oc is the open circuit voltage of the module (V), V mpp is the voltage at maximum power point (V), I mpp is the current at maximum power point (A), γ sc is the current temperature coefficient (A/ °C), and γ oc is the voltage temperature coefficient (V/ °C).

Power
Generation by Wind Turbine.The output power of a WT at time t hour is given by [24] where P WT v is the power produced by the WT at speed v and P rated , V rated , V ci , and V co are the rated power, rated speed, cut-in wind speed, and cut-out wind speed, respectively.To determine the power produced at each state, the probability of that specific state that gives its overall contribution to the total power produced in hour t should be known.The probabilities of solar insolation and wind speeds at state N s are given by the following equations [12]: where p s N s is the probability of solar irradiance in state N s , p v N s is the probability of solar irradiance in state N s , and s 1 and s 2 are the lower and upper limit insolations for state N s .
w 1 and w 2 are the lower and upper limit wind speeds for state N s .
The power produced in an hour t by the PV module and the WT is the summation of the power produced at each state in that hour and is given by the following 5 Journal of Energy equations [12]: where P t pv is the power produced by the PV module at hour t and P t wt is the power produced by the WT at hour t.

Problem Formulation
The primary objective of this research is to strategically determine the optimal placement and sizing of distributed generators (DGs) to minimize total system losses and overall DG cost and enhance voltage profiles.The problem formulation incorporates a multiobjective function that considers these objectives while adhering to specified constraints [17].To achieve optimal solutions, the algorithm leverages network performance indices.The formulated multiobjective function, considering the use of network performance indices, is expressed as follows: where α 1 , α 2 , α 3 , and α 4 are the weights between [0,1] and α 1 + α 2 + α 3 + α 4 = 1, PLI is the real power loss index, QLI is the reactive power loss index, VPII is the voltage profile improvement index, and CFI is the cost factor index.The weights assigned to each impact index represent the relative importance of each index within the analysis.These weights are determined based on the specific analysis conducted.It is essential to acknowledge that these weight values may vary depending on the unique concerns and preferences of the engineer involved [24].
In this study, the significance of each performance index is assessed through a detailed technical analysis of the network.This analysis carefully evaluates the influence of each objective on the overall system performance.The weights assigned to each index, as presented in Table 2, reflect the outcomes of this analysis.

Objective Functions
3.1.1.Power Loss.The active and reactive power losses on any network are formulated as given by the following equation [16]: where P L and Q L are the active and reactive power losses, I i and R i represent the current and resistance of branch i, N bus is the total number of buses in the network, N y gives the number of years of the specified period, and t is time in hours.

Voltage Profile.
The voltage profile function is derived from voltage magnitude at bus i given in the following equation [25]: where V i is the voltage magnitude at bus i and V rated is the desired steady-state voltage.
3.1.3.Total DG Cost.The total DG cost is the summation of the investment, maintenance, and operation cost.
(1) Investment Cost.The investment cost includes the startup costs such as the amount of money spent on construction and installation of individual DG units.This can be calculated using the following equation [16]: where P Rn is the rated power of DG unit n, C inv is the investment cost of DG n, and N DG is the total number of DGs.
(2) Operation and Maintenance Costs.Operation and maintenance costs are the costs for generation, repairing, and renewable DG equipment.This cost is modelled using the following equation [16]: Perform loadflow and apply voltage, thermal, and power balance limits.
Determine the total power produced by DG and apply DG power limit.
Randomly generate number of DGs and distance of ith particle.
Randomly generate number of DGs and distance of ith particle.
Determine the total power produced by DG and apply DG power limit.
Perform loadflow and apply voltage, thermal, and power balance limits.where Q L,with DG is the reactive power losses after DG placement and Q L,without DG is the reactive power losses of the network without DG.

Voltage Profile Improvement Index. The voltage profile improvement index (voltage stability index
) is an indication of how stable a distribution network is.In this paper, VPII is modeled using equation (20), which compares voltage profile of a network with DG with that without DG [26].

VPII =
VP with DG VP without DG 20 3.2.4.Cost Factor Index.The index of the total DG cost is expressed as given by the following equation [26]: where Cost DG is the cost of the selected DG and Cost max DG is the maximum cost of DG at maximum DG penetration.

Constraints.
In optimization problems, constraints define a region that is feasible within the search space.The minimization problem is subjected to the following equality and equality constraints.

Power Balance Constrain.
Complex power injected at each bus should be equal to complex power drawn at that particular bus.This constraint is shown in the following equations [9]: where P Gi is the active power output of the generator at bus i, P Li is the active power load at bus i, Q Gi is the reactive power output of the generator at bus i, Q Li is the reactive power load at bus i, Y ij and θ ij are the magnitude and angle of the admittance related to bus i and j.

Generation Capacity.
The size of the DGs should be within the following size limits [19]:  where P min DGi and P max DGi are the minimum and maximum power outputs injected at bus i.
(1) Bus Voltage Limits.The voltage at each bus should be maintained within their upper and lower limits.These limits are specified by the network operator [20].
where V i,min and V i,max are the minimum and maximum voltage magnitudes at bus i, respectively.For this work, these limits are taken to be (2) Line Loading Limit.The current through any feeder (I i ) should be within the maximum thermal capacity (I i,max ) of that particular feeder.This limit is given by [20] I i ≤ I i,max 27

Proposed PSO for DG Placement
The proposed algorithm as given in Figure 1 is designed to optimally size and place DGs in a distribution network in order to implement and to minimize total system losses, improve voltage profiles, and minimize total DG cost.PSO parameters such as maximum number of iterations, k max , swarm size, nPop, inertia weight, acceleration coefficients, c 1 and c 2 , velocity limits, v min and v max , and position limits, x min and x max The number of DGs (N) and the distance from the main bus d are randomly generated.
Power produced by two DGs is given by P DG,PV = N × P t pv , 28 Optimal power and site are then evaluated by first performing load flow analysis within the constraints given in Section 3.3.Finally, the objective functions are evaluated using the following steps: (1) Parameter definition: x min , x max , and w (2) Define control variables (number of WT and number of PV, distance from the main bus) (3) Set g best to infinity    The seasons are further classified in 24 hours for one day and 182 days per season.Solar irradiance data is collected for 2 years giving 362 days of analysis, while wind speed is data collected for one year making 182 days of analysis.Solar irradiance was modelled using beta PDF, and the wind speed was fitted with the Weibull PDF.
The solar irradiances each hour t were divided into 10 states with each state adjusted to 1 kW/m 2 .Wind speeds were divided into 14 states with each state adjusted by 1 m/ s.For each hour, the mean and standard deviations, probabilities, and power produced were evaluated.Total power produced in an hour t was determined as the summation of the power produced at different states.The results obtained and the ratings of the PV module and WT are shown in Figures 2 and 3

Simulation Results
As stated during the algorithm's development, the primary objective behind strategically positioning and sizing the DGs is to minimize losses and enhance voltage profiles.To analyze the networks, simulations were conducted using the DIgSILENT Power Factory software.The Newton-Raphson method was employed to assess the load flows.Moreover, the multiobjective particle swarm optimization (PSO) algorithm was created and executed using MATLAB.This is done by considering varying loads and renewable sources.For placement, DGs are injected at the main bus bar (substation), where the number of DG units and the distance (in a range of 5-100 km) from the main substation are randomly generated.The output power of the DGs and distance are adjusted to give a minimum fitness value.This is to balance the total losses and power consumed by loads to  The maximum allowable integration of distributed generation (DG) into the network is limited to 148.87 kW, representing 30% of the peak load.
Case 2. The optimum fitness value was achieved with 189 PV modules rated at 69.93 kW and placed at 26 km from the main substation.
Case 3. The optimum fitness value was achieved with 1 wind turbine rated at 100 kW and placed at 42 km from the main substation.
Case 4. The optimum fitness value was achieved with 179 PV modules rated at 65.12 kW and placed at 26 km from the main substation and 1 wind turbine rated at 100 kW placed at 26 km from the main substation.
For each case, variations in DG penetration were plotted where the optimum size and location of the DGs were achieved based on two constraints: thermal limit and voltage profile limit.Simulation results for increased DG penetration and optimum sizing and siting based on the proposed algorithm are as shown in Figures 5-8 In Figure 5, it is observed that, for each case that thermal limit and voltage profile are within the limit of maximum thermal capacity of 53 A and voltage limit of 0 95 ≤ v i ≤ 1 05, for Case 2 (69 96/496 22 = 14%) and Case 3 (100/496 22 = 20%), the same applies to Case 4 in Figure 6.For optimum siting   11 Journal of Energy given in Figures 7 and 8, the line current decreases with increase in distance and voltage profile increases with increase in distance.So, the optimum site becomes the intersection between the graphs of distance-line current and distancevoltage profile.
5.1.1.Active Power Loss Reduction.Integration of DG units into the system results in significant reduction in power losses with highest of 71.55% (from 99.8 MW to 28.29 MW) recorded in Case2.For Cases 3 and 4, the total power losses were reduced by 71.25% and 71.46%, respectively.The variation of these power losses with respect to time for both summer and winter is presented in Figures 9 and 10.
The proposed algorithm has successfully determined the optimal sizes and locations for the integration of distributed generation (DG) units.The resulting active power losses for each of the four cases are presented in Table 3.

Reactive Power Loss
Reduction.While the distributed generators (DGs) in the system are designed to supply active power exclusively, they also contribute to a substantial reduction in reactive power losses.The network under analysis is a rural distribution network primarily catering to agricultural and residential loads.In such a network, the presence of inductive devices, such as induction motors, leads to reactive power consumption.Consequently, voltage drops and fluctuations may occur, compromising the stability and reliability of the electrical system.The reduction in reactive power losses is visually represented in Figures 11  and 12.
5.1.3.Voltage Profile Improvement.The integration of distributed generation (DG) units has a positive impact on the voltage profiles of the network, as demonstrated in Figure 13.Comparing Case 1 with Cases 2, 3, and 4, it is evident that Case 1 exhibits a relatively poorer voltage profile.With DG penetration, the voltages at each bus remain within the permissible limits of 0.95 p.u. and 1.05 p.u. Significant voltage changes are only observed at buses 5, 11, 14, and 21, which are considered weaker buses.The analysis of voltage variation at these weaker buses is presented in Table 4.Among them, bus 5 is identified as the weakest, and its voltage variation is further examined in Figure 14.Notably, both Case 1 and Case 2 exhibit a 7.9% improvement in voltage profile at bus 5 (from 0.909 p.u to 0.9999 p.u).Although Case 2 shows a greater voltage improvement, the disparity in voltage enhancement across buses in each case is not significant.The higher voltage change observed at bus 5 can be attributed to its larger number of induction motors, resulting in a higher consumption of reactive power.5.1.4.Total Cost Analysis.By using the proposed algorithm on cost analysis and applying the cost factor index as given   12 Journal of Energy in equation ( 21), technical and economical indices are shown in Table 5.
Table 5 shows the operational cost, the investment cost, and the total energy losses after a period of 5 years for each case; it can be shown that the cost of active power loss decreases with DG penetration.This improves economic efficiency eventual reduction in carbon emission.
The total cost given is optimized cost since the DG size and placement are determined based on the optimized algorithm given.This eventually reduces feeder losses, hence reduction in total active power loss.

Conclusions
A weighted multiobjective PSO was presented in this work to site and size renewable energy-based DGs on a distribution network.PSO is implemented to minimize total system losses, improve voltage profiles, and minimize total DG cost.Stochastic nature of solar irradiances and wind speeds was modelled using appropriate probabilistic models.Seasonal model was utilized.To effectively analyze the hourly data for each season, a multistate model was employed.This model divided the data into different states, enabling a comprehensive assessment of the system's performance.The control variables for DG placement and sizing were determined based on the number of DG units, including both photovoltaic (PV) and wind turbine (WT) units.The sizes of the DG units were randomly generated and injected at the slack bus to assess their impact on the system.
The proposed optimization technique was tested on one of Namibian's distribution network, and results obtained showed improved voltage profile, minimized active power loss, and reduced cost of the total power supplied.
To realize the benefit of voltage profile improved, total cost reduction, and minimized power loss, the optimized placement and sizes of DG integration in the network were found: For Case 2, the optimum size of PV was found to be 69.93 kW and placed at a distance of 26 km from the main substation.For Case 3, the optimum size of WT was found to be 100 kW and placed at a distance of 42 km from the main substation.For Case 4>, the optimum size for both PV and wind-based DGs was found to be 65.12 kW and 100 kW, respectively.The PV farm was placed at a distance of 37 km from the main substation and the wind farm at 20 km from the substation.By achieving these positive outcomes, the study underscores the potential benefits of DG integration and emphasizes the importance of carefully determining the appropriate sizes and locations for optimal performance.

Figure 1 :
Figure 1: Algorithm for optimal sizing and placement of DGs on distribution network.

Figure 2 :Figure 3 :
Figure 2: Estimated power produced by PV module at different states (kW).

4. 2 .
Wind, Solar Data, and Load Data Analysis.The proposed algorithm was tested on a Namibian distribution network, Vhungu-Vhungu 11 kV, as detailed in Section 5.

Table 1 :
Literature review of use of MOPSO technique for optimizing DG sizing and placement in a distributed network.
[25]= C oper,maint + C investment 17 3.2.Performance Index Formulation 3.2.1.Real Power Loss Index.The real power loss index based on network with DG placement and without DG placement is modeled as given in the following equation[25]: Print the optimum size and location 4.1.Wind, Solar Data, and Load Data Analysis.Strategic integration of renewable energy in the network appropriate analysis of wind, solar, and load data is important.Seasonal model where a year is divided in winter (May to October) and summer (November to April) is used in data collection.

Table 5 :
Technical and economical indices for 5 years for Vhungu-Vhungu.