Research Article | Open Access

Volume 2019 |Article ID 5131696 | https://doi.org/10.1155/2019/5131696

Vasilija Sarac, Tatjana Atanasova-Pacemska, "Multiparameter Analysis for Efficiency Improvement of Single-Phase Capacitor Motor", Mathematical Problems in Engineering, vol. 2019, Article ID 5131696, 13 pages, 2019. https://doi.org/10.1155/2019/5131696

# Multiparameter Analysis for Efficiency Improvement of Single-Phase Capacitor Motor

Revised03 Apr 2019
Accepted16 Apr 2019
Published22 May 2019

#### Abstract

Single-phase motors are known for their small power ratings and their usage in various household appliances. Although they are not large electricity consumers, their widespread application contributes to the overall electricity consumption. In addition, standard IEC 60034-30-1:214 defines the efficiency levels for single- and three-phase motors and stipulates the increased electrical efficiency for the electrical motors. Therefore, this paper sets the parametric analysis of permanently split capacitor motor with five different design parameters that have impact on the efficiency of the motor. As an output from the parametric analysis, two different optimized motor models are obtained with increased efficiency. The impact of each parameter on motor efficiency—as well as on the other operating characteristics, like starting torque, overloading capacity, rated current, starting current, total losses, and power factor—is analyzed and adequate conclusions are derived. The obtained motor models are verified with Finite Element Method (FEM) for magnetic flux density distribution.

#### 1. Introduction

Capacitor motors belong to the single-phase motors that have a capacitor in the auxiliary winding usually aimed for motor starting. In the case of the permanently split capacitor motor, there is only one capacitor, placed in the auxiliary winding, which is permanently in operation, during starting as well as during motor running. In general, the single-phase motors have low power ratings but they are widespread in the household appliances and everywhere where single-phase supply is available. This makes the electricity consumption from the single-phase motors a considerable portion of the overall electricity consumption in the households and industry. In March 2014, the standard IEC 60034-30-1 has been published. It defines four efficiency classes for induction motors (single- and three-phase), IE1 standard efficiency, IE2 high efficiency, IE3 premium efficiency, and IE4 super premium efficiency (Figure 1). From January 2017 the legally specified minimum efficiency IE3 must be maintained for power ratings of electrical motors from 0.75 kW up to 375 kW. Improving the motor efficiency has become an imperative for electrical motor producers. Nevertheless, the production costs should be maintained on competitive level through analyzing the various optimization aspects of the motors.

#### 2. Mathematical Theory for Single-Phase Motor Modeling

Despite its simple construction, the mathematical model of the single-phase motors, which describes the various motor operating modes, is quite complex due to the existence of the rotating elliptical electromagnetic field in the motor air gap. This air gap electromagnetic field is highly unsymmetrical; the application of the well-known theory and the mathematical models for the three-phase symmetrical induction motors is improper and inaccurate. Therefore, the unsymmetrical magnetomotive force (mmf) currents and voltages corresponding to the two windings of single-phase induction machines may be decomposed in two symmetrical systems (Figure 2) which are the forward and backward components of the two-phase system [15].

Here following relations are valid:where and are general variables associated with the auxiliary winding and main winding, respectively.

The superposition principle yields the following.

The first step in the motor analysis is to determine all motor parameters: , main stator winding resistance; , main winding leakage reactance; , auxiliary stator winding resistance; , auxiliary stator winding leakage reactance; , magnetizing reactance; , rotor winding resistance; and , rotor winding leakage reactance [15]. Calculation of the motor parameters is based on the motor geometry obtained from the producer.

Based on the method of symmetrical components, impedance parameters , , and are calculated, necessary for obtaining the symmetrical components of the stator currents and electromagnetic torque [15]:where s is the motor slipwhere C is the capacitance in the auxiliary winding and a is the ratio of the winding turns of the main and auxiliary stator winding. The symmetrical components of the currents of the main stator winding are calculated fromwhere Vs is the motor power supply voltage.

The supply line current is found from the following.The currents in the main winding and the auxiliary stator winding are calculated, respectively.Components of the electromagnetic torque-direct and inverse are found:where p is the pair of the motor poles.

Electromagnetic torque is as follows.Motor input power is as follows.The power factor becomes as follows.The motor output power P2 is obtained from the input power P1 after subtracting all motor losses , which can be divided into several groups: stator and rotor ohmic losses, capacitor losses, iron core losses, and mechanical losses.

Stator ohmic losses are calculated at temperature of 75°C as the sum of the ohmic losses in the main and the auxiliary winding.Capacitor losses can be found fromwhere is the capacitor impedance .

According to [15] rotor ohmic losses at rated load can be approximately calculated aswhere andCore losses are mainly eddy current and hysteresis losses that in analytical form can be expressed as follows.

is the eddy current loss constant.

is the maximum flux density.

is the hysteresis loss constant.

is the material dependant Steinmetz constant (1.6 to 2).

Mechanical losses are mainly covering friction and windage losses. These losses can be determined by driving the motor at the rated speed with no-load or excitation. For the machines that operate at constant or nearly constant speed, these losses are constant. For an accurate calculation of motor losses, it is necessary to have an accurate calculation of motor parameters (resistances and reactances) as well as the currents in all motor windings. In addition, the magnetic flux density in the motor cross-section should be calculated, a parameter that is often more or less accurately predicted in motor analytical calculations. This emphasizes the need for an accurate computer-aided design of the motor capable of handling various design modifications, necessary for improving the motor efficiency factor (η), which can be finally obtained from the ratio of the output power P2 and the input power P1.

#### 3. Computer-Aided Design of the Motor Models

Accurate modeling of the motor in a software program for automated computing of parameters and operating characteristics of the single-phase capacitor motor is an important step in the parametric analysis for efficiency improvement. Exact motor geometry is modeled with respect to the magnetic core dimensions and properties of the core material, i.e., magnetization curve B=f(H), with B being the magnetic flux density and H being the magnetic field strength (Figures 3 and 4).

Properties of core laminations are defined in accordance with lamination type DL 80, RK502.4. Electric properties of the materials are modeled as well. Stator winding, consisting of main and auxiliary winding made of copper, is placed in the main and the auxiliary slots. Rotor winding is of squirrel cage type, made of aluminum. Motor is aimed for constant power operation and it is modeled in the software program Maxwell using available data from the motor’s producer. The accuracy of the developed starting model (BM) is verified by comparing the output results from the model and available data from the producer. Table 1 presents this comparison.

 BM Producer data Rated power (W) 124 124 Rated torque (Nm) 0.416 0.411 Rated speed (rpm) 2846 2880 Rated current (A) 1.35 1.32 Maximum power (W) 210 210 Maximum torque (Nm) 0.83 0.8

Once the starting model has been verified as sufficiently accurate, the parametric analysis is set by varying several motor parameters within the predefined boundaries. The variation ranges of the parameters are presented in Table 2.

 Ranges of variation Step Nr. of conductors in main slot (/) 110÷124 1 Nr. of conductors in auxiliary slot (/) 125÷135 1 Rotor diameter (mm) 65÷65.6 0.1 Motor core length (mm) 32÷37 1 Rotor slot opening (mm) 0÷1.5 0.75

As a first step in the parametric analysis, four varied parameters were defined: the number of the conductors in the main and the auxiliary stator slot, length of magnetic core, and rotor diameter, i.e., the length of the air gap between the stator and the rotor. The parametric analysis resulted in 8085 different combinations that were solved. Out of these solutions, the best result or the best motor model in terms of the maximum efficiency at rated load was selected (model M1). The second parametric analysis was defined by adding one more varied parameter, i.e., the length of the rotor slot opening, and by taking into consideration the previous four varied parameters. This second parametric analysis resulted in 17325 different combinations. Their solutions gave the second motor model (M2) with better efficiency than BM and M1. In the next step, motor models are verified by FEM for magnetic flux density distribution. The flux density B is calculated from the magnetic vector potential A. In order for Maxwell’s equations to be solved, the complete machine cross-section is divided into numerous elements forming the mesh of finite elements (Figure 5). The electromagnetic field problems are solved by solving Maxwell’s equations in a finite region of space with appropriate boundary conditions and user-specific initial conditions in order to obtain solution with guaranteed uniqueness.where E is the electric field, D is the electric displacement equal to εE, ε is the permittivity, B is the magnetic flux density, and H is the magnetic field, equal to the product of material magnetic permeability and magnetic flux density, μB. J is the current density, σE, where σ is the conductivity. ρ is the charge density.

The FEM discretization of the domain of the analyzed object produces a set of matrix differential equations. They are solved with the time decomposition method (TDM). The domain is decomposed along the time axis and all time steps are solved simultaneously instead of solving them time step by time step. The nonlinear matrix equations are linearized for each of the nonlinear iterations. As an output from the FEM model, motor torque is calculated for all motor models. The comparison of the obtained torque from all FEM models with the torque from the parametric analysis verifies the accuracy of the FEM model.

#### 4. Results and Discussions

##### 4.1. Parametric Analysis

The output results from the parametric analysis give an overview of all important motor parameters and data. Motor parameters of all three models (BM, M1, and M2) are presented in Table 3 together with the output results from the parameters variations that resulted in the best motor models in terms of efficiency. All parametric analyses are run for constant power operation; i.e., all motors should have the same output power. The type of the winding is two-layer sinusoidal winding.

 Motor parameter BM M1 M2 Nr. of conductors per layer in main winding - (/) 117 124 124 Nr. of conductors per layer in auxiliary winding - (/) 130 125 125 Core length-L (mm) 35 37 37 Rotor diameter- (mm) 65.4 65.6 65.6 Rotor slot opening –RSO (mm) 1.5 1.5 0 Main winding resistance at 20°C -(Ω) 25.376 27.16 27.16 Auxiliary winding resistance at 20°C-(Ω) 69.14 67.05 67.056 Main winding reactance-(Ω) 16.24 19.18 19.38 Auxiliary winding reactance-(Ω) 24.66 23.8 24 Rotor resistance-(Ω) 22.31 25.66 25.98 Rotor reactance-(Ω) 3.53 4.81 21.525 Magnetizing reactance-(Ω) 387.18 702.3 784.57 Capacitance-C (μF) 6 6 6

For the best motor model in terms of efficiency besides motor parameters, motor operation in three typical operating modes, i.e., the rated load, the locked rotor, and the no-load, is analyzed. Table 4 provides the motor operating characteristics for the rated load. Here, an overview of all motor losses is presented as well. Mechanical losses (frictional and windage losses) were calculated as 9% of the total losses as often they are within range of 8-12 percent of the total losses. The no-load operating characteristics are presented in Table 5 together with the motor locked rotor (starting) characteristics. Detailed presentation of motor operation for different speeds is presented in Figures 6, 7, 8, and 9. Figure 6 presents the input (line) current for different speeds.

 Motor characteristic BM M1 M2 Rated speed (rpm) 2846.48 2831.72 2831.72 Line current (A) 1.35 1.1788 1.159 Main winding current density (A/mm2) 4.422 3.6562 3.608 Auxiliary winding current density (A/mm2) 3.6 3.38 3.284 Main winding current (A) 1.0166 0.84 0.829 Auxiliary winding current (A) 0.3693 0.346 0.336 Capacitor losses (W) 72.378 63.54 59.98 Iron core losses (W) 4.487 3.93 3.13 Stator ohmic losses (Ω) 43.345 33.083 31.9 Rotor ohmic losses (Ω) 15.91 16.603 16.34 Mechanical losses (W) 13.424 13.28 13.28 Total losses (W) 149.545 130.454 124.67 Input power P1 (W) 273.552 254.352 248.7 Output power P2(W) 124 123.89 124.04 Efficiency factor (eta) η(%) 45.332 48.7 49.87 Output torque T (Nm) 0.416 0.4178 0.418 Maximum torque (Nm) 0.83 0.775 0.754 Power factor (/) 0.919 0.98 0.975 Rated slip (/) 0.051 0.056 0.056
 Motor characteristic BM M1 M2 Main winding reactance (Ω) 15.96 18.7 18.7 Auxiliary winding reactance (Ω) 24.31 23.37 23.4 Rotor reactance (Ω) 3.24 4.3 47.59 Magnetizing reactance (Ω) 283.34 524.48 552.5 No-load line current (A) 1.15 0.6799 0.653 No-load current main winding (A) 0.953 0.423 0.392 No -load current auxiliary winding (A) 0.413 0.39 0.386 No-load iron core losses (W) 5.27 4.787 3.82 No-load input power (W) 156.244 118.96 114.97 No-load power factor (W) 0.614 0.795 0.8 No-load speed (rpm) 2985.91 2986.7 2986.7 Starting torque (Nm) 0.116 0.11 0.111 Starting line current (A) 4.077 3.658 3.53 Starting current main winding (A) 3.9325 3.519 3.4 Starting current auxiliary winding (A) 0.2704 0.2712 0.272

From Table 4 it can be observed that motor line current is decreased in both optimized models. The decrease of the line current is mainly the result of the increased resistance of the main winding due to the increase of the number of conductors in the main winding. The decreased line current leads to the decreased input power P1 or electricity consumption (15). The decreased input power affects the motor efficiency factor. The decrease of the current resulted in decreased ohmic losses in the stator windings (17). The parametric analysis is done with the motor output power as a constraint; i.e., output power should be kept constant which implies almost unchanged output torque as the relation between the output power and the torque is directly proportional. Some variations of the rated speed in different motor models that affect the rated output torque are possible, but in these particular models, they are not so emphasized. By taking into consideration that efficiency factor is defined with the ratio of the output and the input power, keeping the output power constant is forcing the decrease of the input power and the motor losses in order to achieve better efficiency factor than in the starting model. The smaller the motor losses, the better the obtained efficiency factor. As the electricity consumption of the motor is decreased, this leads to an energy efficient motor. This is the case with M1 and M2, with both models having smaller total losses than BM.

The influence of the variation of different design parameters on overall motor efficiency is presented in Figure 10. Diagrams are obtained by keeping the starting motor configuration and by varying only one parameter (i.e., Figure 10(a), is variable, =130, L=35 mm, =65.4 mm, RSO=1.5 mm). In Figure 10(b) is variable while all the other parameters remain the same as in the starting model, BM. On the same principle, the other results in Figure 10 are obtained. The efficiency increases with the increase of the conductors in the main winding, the decrease of the main winding current, and the decrease of the copper losses.

Motor output power (P2) is proportional to product D2L where D is the stator bore diameter and L is the core length [15]. Consequently, the increase of motor length gives the increase of the motor output power and the increase of the motor efficiency (η) since output and input power of the motor are related to (22). This type of the motor has increased current magnitude with backward time shift and the increased torque pulsations as the motor approaches the rated speed. To reduce this, the impedance in the auxiliary winding should be reduced and consequently the losses are minimized; i.e., number of conductors in auxiliary winding is decreased (Figure 10(b)).

##### 4.2. FEM Analysis

Parametric analysis serves as a starting point for deriving motor FEM model for verification of proposed models (M1 and M2) in terms of magnetic flux density distribution. Figure 11 presents the distribution of the magnetic flux density for all motor models.

There are some general recommendations for the flux density in the specific parts of the machine. For the stator and the rotor teeth, it should be less than 1.8 T, and for the core (yoke) between 1.3 and 1.5 T [15]. Following these general recommendations, the flux density distribution is within recommended ranges for BM and M1. As for the second optimized model M2, bigger values of the flux density can be observed in the rotor teeth due to the motor construction with closed rotor slots. Yet, this value of the flux density is still below the point of the core saturation with respect to the built-in core material as the saturation point is approximately at 2T (Figure 3). FEM models of the motors are verified in terms of their accuracy by calculating the output torque for one fixed speed. It should be noted that presented diagrams of the torque in Figure 12 do not represent the motor transient characteristics usually plotted during motor acceleration from zero up to the rated speed. The presented diagrams are calculating motor torque at one constant speed, i.e., rated speed, for each of the motor models, within the whole interval. Table 6 presents the comparison between the torque values from the parametric analysis and from the FEM models at rated load operating mode. As the torque in the FEM models has pronounced oscillations, it is calculated as an average value within the time interval from 175 to 200 ms. The presented results in Table 6 and their similarity verify that the FEM models are sufficiently accurate. Therefore, the presented results of the flux density distribution are considered reliable and contribute to the overall estimation of the motor design and the performance.

 Output torque (Nm) Parametric analysis FEM analysis BM 0.4135 0.4271 M1 0.4178 0.4293 M2 0.418 0.4435

#### 5. Conclusion

Single-phase capacitor motors are large electricity consumers if their relative high number in various electrical household appliances is taken into consideration. Therefore, the slight improvement of their efficiency contributes to the overall energy efficiency. The paper presents the impact of several motor parameters such as number of conductors in main and auxiliary winding, machine length, air gap length, and opening of the rotor slots on the efficiency of the permanently split capacitor motor. Parametric analysis was run in Maxwell software which resulted in two improved motor models, M1 and M2, with increased efficiency of 7.4% and 10%, respectively, compared with the starting motor model. The difference between models M1 and M2 is the width of rotor slot opening, with the latter having completely closed rotor slots. The models obtained from the parametric analysis were verified with FEM in terms of the magnetic flux density distribution. The results from the FEM analysis proved that model M1 is well designed with respect to the magnetic flux density distribution as the flux density in machine cross-section is within the recommended intervals for the various parts of the machine. Model M2 has the highest value of flux density in rotor teeth although this value is still below the point of core saturation for the built-in magnetic material. The accuracy of the derived FEM models is verified with the similarity of the obtained results of the output torque from FEM and the parametric analysis. By using the latest advancements in the designing programs for the electrical machines, different motor models can be derived in relatively short computing time, resulting in the best possible motor design in terms of predefined objective function. Alterations in the motor design are easily inputted; the results are quickly obtained giving motor’s designers freedom to explore the different motor variants and to find the best design for a certain application. This paper provides a comprehensible overview of the design parameters of the single-phase permanent split capacitor motor and their impact on the motor efficiency. In order to have comparable results, all motor models are derived with the same capacitor as it affects the motor performance. Authors’ further research will be focused on finding the optimal capacitor value for providing the best starting and operating conditions of the motor.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

The work was supported by the project “Contributions in Mathematical Theory, Mathematical Modeling and Their Application”, financed by “Goce Delcev” University.

#### References

1. “Siemens: Efficiency classes for IEC motor,” https://www.industry.siemens.com/drives/global/en/motor/low-voltage-motor/efficiency-standards/pages/line-motors.aspx. View at: Google Scholar
2. J. Li and J. Shen, “Influence of mechanical parameters on power efficiency of induction motors,” in Proceedings of the 2014 17th International Conference on Electrical Machines and Systems (ICEMS), pp. 2555–2560, Hangzhou, China, October 2014. View at: Publisher Site | Google Scholar
3. D. Žarko, S. Frljić, and S. Stipetić, “Design of premium efficiency (IE3) induction motors using evolutionary optimization and scaling laws,” Przegląd Elektrotechniczny, vol. 92, no. 12, pp. 181–184, 2016. View at: Google Scholar
4. C. Mademlis, I. Kioskeridis, and T. Theodoulidis, “Optimization of single-phase induction motors - Part I: Maximum energy efficiency control,” IEEE Transactions on Energy Conversion, vol. 20, no. 1, pp. 187–195, 2005. View at: Publisher Site | Google Scholar
5. X. Wang, D. Meng, and Y. Xu, “Analysis of single-phase induction motors efficiency optimization control with hall sensor,” Advances Science and Technology Letters, vol. 53, pp. 357–360, 2014. View at: Publisher Site | Google Scholar
6. V. Rashtchi and A. Ghasemian, “Efficiency optimization of induction motor drive using modified particle swam optimization,” in Proceedings of the International Conference on Electrical, Electronics and Instrumentation Engineering, pp. 14–18, Johannesburg, South Africa, 2013. View at: Google Scholar
7. A. M. Sharf and A. A. El-Gammal, “A new technique for maximum efficiency and minimum electrical energy consumption of single phase induction motors based on particle swarm optimization (PSO),” in Proceedings of the 16th National Power Systems Conference, pp. 12–17, Hyderabad, India, 2010. View at: Google Scholar
8. K. Makowski and M. J. Wilk, “Optimization of a single-phase capacitor induction motor by applying a surrogate field-circuit model,” COMPEL - The International Journal for Computation And Mathematics in Electrical And Electronic Engineering, vol. 33, no. 6, pp. 1891–1903, 2014. View at: Publisher Site | Google Scholar
9. S. Szenasi and I. Felde, “Configuring genetic algorithm to solve the inverse heat conduction problem,” Acta Polytechnica Hungarica, vol. 14, no. 6, pp. 133–152, 2017. View at: Google Scholar
10. N. Horiyat, S. M. Shariatmadar, and V. Amir, “Optimization of an experimental unsymmetrical-unbalanced two-phase induction motor using the HGAPSO hybrid technique and the finite element method for increasing efficiency and reducing torque ripple,” Turkish Journal of Electrical Engineering & Computer Sciences, vol. 24, pp. 4154–4173, 2016. View at: Publisher Site | Google Scholar
11. P. K. Choudhary and S. P. Dubey, “Efficiency optimization of induction motor drive in steady-state using artificial neural network,” in Proceedings of the 2016 International Conference on Computation of Power, Energy Information and Commuincation (ICCPEIC), pp. 295–302, Melmaruvathur, Chennai, India, April 2016. View at: Publisher Site | Google Scholar
12. L. Ramesh, S. P. Chowdhury, S. Chowdhury, and Y. H Song, “Efficiency optimization of induction motor using a fuzzy logic based optimum flux search controller,” in Proceedings of the International Conference on Power Electronics, Drives and Energy Systems, pp. 1–6, 2006. View at: Publisher Site | Google Scholar
13. A. G. Yetgin and M. Turan, “Efficiency optimization of slitted–core induction motor,” Journal of Electrical Engineering, vol. 65, no. 1, pp. 60–64, 2014. View at: Publisher Site | Google Scholar
14. A. G. Yetgin and M. Turan, “Efficiency improvement in induction motor by slitted tooth design,” Technical Gazette, vol. 24, no. 5, pp. 1291–1296, 2017. View at: Publisher Site | Google Scholar
15. I. Boldea and S. A. Nasar, The Induction Machines Design Handbook, CRC Press, 2nd edition, 2010.
16. C. E. Feraga and A. Bouldjedi, “Performance of a photovoltaic pumping system driven by a single phase induction motor connected to a photovoltaic generator,” Automatika, vol. 57, no. 1, pp. 163–172, 2016. View at: Publisher Site | Google Scholar
17. M. S. Kurt and A. Fenercioğlu, “Rotor slot distance effects on output parameters in single phase induction motors,” Hittite Journal of Science & Engineering, vol. 5, no. 1, pp. 31–35, 2018. View at: Publisher Site | Google Scholar
18. H. Ahn, K. Kim, J. Choi, E. Kim, and Y. Lim, “Rotor design for an efficient single-phase induction motor for refrigerator compressors,” Energies, vol. 9, no. 3, pp. 1–16, 2016. View at: Publisher Site | Google Scholar
19. G. Y. Zhou and J. X. Shen, “Current harmonics in induction machine with closed-slot rotor,” IEEE Transactions on Industry Application, vol. 53, no. 1, pp. 134–142, 2017. View at: Publisher Site | Google Scholar
20. M. Y. Kaikaa, F. Babaa, A. Khezzar, and M. Boucherma, “Analytical analysis of rotor slot harmonics in the line current of squirrel cage induction motors,” Journal of Electrical Engineering, vol. 57, no. 1, pp. 12–19, 2006. View at: Google Scholar
21. M. Dems, K. Komeza, S. Wiak, and S. F. Coya, “Influence of the closing slots on the additional losses in the induction motor,” COMPEL-The International Journal for Computation and Mathematics in Electrical and Electronic Engineering, vol. 34, no. 2, pp. 531–549, 2015. View at: Publisher Site | Google Scholar
22. T. Tudorache and L. Melcescu, “FEM optimal design of energy efficient induction machines,” Advances in Electrical and Computer Engineering, vol. 9, no. 2, pp. 58–64, 2009. View at: Publisher Site | Google Scholar
23. A. Simion, L. Livadaru, S. Mihai, A. Munteanu, and C.-G. Cantemir, “Induction machine with improved operating performances for electric trucks. A FEM-based analysis,” Advances in Electrical and Computer Engineering, vol. 10, no. 2, pp. 71–76, 2010. View at: Publisher Site | Google Scholar
24. Z. Li, L. Zhang, X. Yu, and Q. Wang, “Modelling and analysis of two-degree of freedom permanent magnet synchronous generator,” EEA - Electrotehnica, Electronica, Automatica, vol. 65, no. 3, pp. 63–69, 2017. View at: Google Scholar

Copyright © 2019 Vasilija Sarac and Tatjana Atanasova-Pacemska. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.