Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2464-2470 2464 www.etasr.com Singh and Verma: Performance Analysis and Comparison of Symmetrical and Asymmetrical Dual Stator … Performance Analysis and Comparison of Symmetrical and Asymmetrical Dual Stator Induction Generators for Wind Energy Conversion Systems Ramesh Singh Department of Electrical Engineering Delhi Technological University Delhi, India Vishal Verma Department of Electrical Engineering Delhi Technological University Delhi, India Abstract—The emergence of micro/nano level wind generation has opened the research on induction generator (IG) topologies having easier and finer control available in multiphase generators. To establish the suitability of multiphase generators for wind generators the analysis of their performance based on the developed model for the same rating of six phase symmetrical (60˚) and asymmetrical (30˚) IGs is presented. A comparative performance evaluation of grid-excited six phase symmetrical and asymmetrical IGs also known as Dual Stator Induction Generators (DSIGs) is presented through simulation results in MATLAB/SIMULINK environment amidst load perturbations, limited variation of wind speed and perturbations in voltage and frequency of the non stiff micro-grid to which they are connected. Based on the performance indices like flux of direct and quadrature axis, speed variations, terminal voltage drop/rise, range of operational speed variation etc., a comparative analysis with the help of the results is drawn to establish the suitability of asymmetrical multiphase IGs for grid connected wind generators. Keywords-dual stator induction generator; wind energy conversion system I. INTRODUCTION Recent advances and research in the field of power electronics has changed the scenario of power generation from large scale generation to distributed micro, nano and pico generation. Due to their suitability for power generating at variable speed and having rugged brushless structure with maintenance, free operation squirrel cage induction generators (SCIGs) are fast finding place for wind generation. The specific requirements for grid connected generator coupled with wind turbines include better regulation of terminal voltage, better stability, uniform distribution of flux, higher electrical inertia to meet the transients, higher efficiency and cost effectiveness and better low voltage ride through (LVRT) capability. Normally multistage gear train/single stage plenary gear train is employed for coupling SCIG to enable its working in a wide range of speed variation, but larger gear ratio and multistage gear train reduces efficiency to a great extent by incurring losses [1]. Reported literature reveals that full rated back to back voltage source converter (BTB-VSC) may be employed to harvest energy in higher speed range by vector or scalar control of IG side VSC [2]. Such topologies increase the cost of the system as full power passes through the DC link. To overcome such issues DFIGs were introduced where direct grid connection of stator winding was possible for operation amidst large variation of wind speed [2]. Such configuration removes the dependency on gear train thus becomes more efficient and cost effective for large power range [3]. Employing wound rotor induction machine increases the machine cost and requires regular maintenance [3]. Moreover the LVRT capability of DFIG is also very poor [4]. To offset the increased cost and LVRT capability such machines are well suited for coupling with big wind turbines [4]. To make the suitability of DFIG regarding free structure maintenance i.e. getting rid of slip-rings and IG with better LVRT brushless doubly fed induction generator (BDFIG) which has two stator windings with one nested loops rotor winding is reported in literature [5]. BDFIGs have demonstrated better LVRT capability and have shown their suitability to work with variable speed operation [5]. The increased cost of BDFIGs like DFIG leaves the scope for furthering in research with multi stator IGs for wind energy applications. Other reported topology included split phase induction generators (SPIGs) where there are two stator windings with squirrel cage rotor but the poles of both stators are different. Sharing of power between them is made via rotor, making transient response slower and incurring considerable losses [5-7]. Dual stator induction generators (DSIGs) have emerged as suitable alternative to both DFIGs and BDFIGs where the two stator windings are used with squirrel cage rotor, and both stators have the same number of poles and are excited with power frequency, thus establishing transaction of power through both stator windings independently and via rotor. In addition they establish flux coupling between both stators enhancing transient characteristics and efficiency [6, 8, 9]. A narrow speed operation without multi-gear systems remains the only limiting factor of DSIG [8]. The demerit is offset as the power density is high in such machines and rotor is cheaper with normal stator construction. The reliability of power generation is also increased as now two three phase sets deliver power to the connected three phase loads [6, 9, 10]. Normally Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2464-2470 2465 www.etasr.com Singh and Verma: Performance Analysis and Comparison of Symmetrical and Asymmetrical Dual Stator … the DSIG has been proposed with both symmetrical and asymmetrical winding structures for isolated application and with different kinds of load. Therefore, a need to critically evaluate the performance of both kinds of generators under real life conditions emerges. DSIGs in case of low wind speed can act as two winding transformers to transect power between two stator windings connecting two power circuits and in case of adequate wind speed, generate power from wind and supply to the two networks connected across the two windings, thereby reducing the effective loading on either or both feeders. An assembly / cluster of such small rated DSIGs connected in parallel may support additional benefits. In this paper, the performance of DSIGs with symmetrical and asymmetrical stator winding coupled with weak grid is evaluated on the basis of voltage regulation and LVRT capability which in turn is based on establishment of air gap flux. Since the saturation of flux poses an impediment in operation of IGs, a careful study of flux in air gap is also carried out. The connected load is perturbed intentionally to observe and compare the performance of the mentioned above two types of DSIGs with respect to direct and quadrature axis flux and terminal voltage. The fault ride through capability is also evaluated by creating a LG fault. The stability during fault and the voltage recovery performance is also compared. The simulated results under MATLAB environment are presented and evaluated for different operating conditions. II. SYSTEM CONFIGURATION In the proposed scheme a 5 HP, 4 Pole DSIG is connected directly to two different sets of three phase distribution feeders (DFs), DF 1 and DF 2 with appropriate relative phase shifts in respective set of phases. Figure 1 shows the test bench configuration used for comparative evaluation of symmetrical and asymmetrical DSIG. According to the connected IG appropriate relative phase shifts in respective set of phases is provide on the DF1. Reactive power demand of DSIG is supported by the distribution feeders to which the winding sets are connected respectively. Two separate capacitor banks (C= 9 µF) are provided at each stator winding terminal to partially support reactive power requirement and filter the currents for range of current harmonics and restraining them to enter the terminals of DSIG. Different types of passive loads are connected to both set of three phase windings of DSIG to study the effect of load dynamics seen by the proposed configurations. III. MATHEMATICAL MODELING Two balanced three phase windings displaced to each other by 30˚(asymmetrical) and 60˚(symmetrical) are considered for modeling of DSIG under the following assumptions: Air-gap is uniform and the windings are sinusoidally distributed around the air-gap and the magnetic saturation and core-losses are neglected. The voltage equations are written in d-q axis for DSIG in synchronous reference frame as: )()(. 11111 dsqsqssqs pIRV   (1a) )()(. 11111 qsdsdssds pIRV   (1b) )()(. 22222 dsqsqssqs pIRV   (2a) Fig. 1. Block diagram of the considered system. )()(. 22222 qsdsdssds pIRV   (2b) ))(()(. drrqrqrrqr wpIRV   (3a) ))(()(. qrrdrdrrdr wpIRV   (3b) )()(. 212111 qrqsqsmqsqslmqslsqs iiiLiiLiL  (4a) )()(. 212111 drdsdsmdsdslmdslsds iiiLiiLiL  (4b) )()(. 212122 qrqsqsmqsqslmqslsqs iiiLiiLiL  (5a) )()(. 212122 drdsdsmdsdslmdslsds iiiLiiLiL  (5b) )(. 21 qrqsqsmqrlrqr iiiLiL  (6a) )(. 21 drdsdsmdrlrdr iiiLiL  (6b) Since the model is designed for squirrel cage induction generator, therefore, qrV =0 and drV =0. Thus the real and reactive power flowing in the two circuits can be represented as in equations (7a-d). 11111 .. qsqsdsds ivivp  (7a) 11111 .. dsqsqsds ivivq  (7b) 22222 .. qsqsdsds ivivp  (7c) 22222 .. dsqsqsds ivivq  (7d) Where χqs1, χqs2 are stator q-axis flux linkage components, χds1, χds2 are stator d-axis flux linkage components, χ'qr, χ'dr are rotor q-axis and d-axis flux linkage components, iqs1, iqs2 are Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2464-2470 2466 www.etasr.com Singh and Verma: Performance Analysis and Comparison of Symmetrical and Asymmetrical Dual Stator … stator q-axis current components, ids1, ids2 are stator d-axis current components, i'qr, i'dr are rotor q-axis and d-axis current components, Lls is stator leakage inductance, Lm is air gap inductance, Llm is stator mutual leakage inductance, L'lr is rotor leakage inductance, þ is (d/dt), ω and ωr are synchronous flux and rotor speed, Rs and Rr are stator and rotor resistances, vqs1, vds1, vqs2, vds2, v'qr and v'dr are the stator and rotor q-axis and d- axis voltages respectively. The equivalent circuit of DSIG generator in d-q reference frame is shown in Figure 2(a) and 2(b) respectively. The developed modeling equations do not differentiate between the symmetrical and asymmetrical machines since the instantaneous parameters represent the waveforms, therefore more clarity on the subject would be visible through its depiction on phasor diagrams. (a) (b) Fig. 2. (a) Equivalent circuit of DSIG in d-axis (b) Equivalent circuit of DSIG in q-axis IV. COMPARATIVE PHASOR BASED ANALYSIS It is important to analyze both types of DSIG’s from the flux distribution point of view, which can be depicted separately by analysing the effect of loading on both leading and lagging windings and their contribution on direct and quadrature axes. The differences in operating characteristics are visible and strongly influence the results shown with both DSIGs. Figure 3 shows the phasor diagram showing leading and lagging winding currents of asymmetrical and symmetrical DSIG in terms of direct axis and quadrature axis. To study the effect of leading and lagging winding current on one another the projection of one current is taken on the other current and the resultant leading and lagging current phasor is drawn. For simplification and ease of understanding the resultant phasors for both DSIGs are calculated taking leading winding as a reference. The projections of a current phasor are represented by dotted lines with the same color to that of the current phasor. Superscripts ‘*’ and ‘**’ are used to represent the direct axis and quadrature axis projections of lagging winding on leading winding axis and superscripts ‘#’ and ‘##’ are used to represent the direct axis and quadrature axis projections of leading winding on lagging winding axis. Thick solid hashed lines represent the net leading and lagging winding current after incorporating the projections of their respective direct and quadrature current components. Figures 3(a), 3(d) for instance represent the direct and quadrature components for current phasor of leading and lagging winding asymmetrical DSIG. Figure 3(b) represents the projection of lagging winding currents upon leading winding currents of asymmetrical DSIG. Figure 3(c) represents the resultant leading winding current of asymmetrical DSIG incorporating the projections shown in Figure 3(b). Similarly Figures 3(e), 3(f) represent the projections of leading winding currents on the lagging winding currents and the resultant lagging winding currents after incorporating the projections respectively for asymmetrical DSIG. Figures 3(g-l) represent the current phasor for symmetrical DSIG for projection on leading and lagging winding currents. The quadrature and direct axis components of current phasors depicted in Figures 3(c), 3(f), 3(i) and 3(l) are also represented mathematically in equations (8)-(9), (11)- (12), (14)-(15) and (17)-(18) respectively. The net current phasors are given by (10), (13), (16) and (19).  For Asymmetrical DSIG For projection of lagging winding side to leading winding side: 37.20,1 jI ASq  (8) 037.1,1 jI ASd  (9) 37.237.1,1 jI ASNET  (10) For projection of leading winding side to lagging winding side: 19.169.0,2 jI ASq  (11) 18.105.2,2 jI ASd  (12) 37.236.1,2 jI ASNET  (13)  For Symmetrical DSIG For projection of lagging winding side to leading winding side: 37.20,1 jI Sq  (14) 063.0,1 jI Sd  (15) 37.263.0,1 jI SNET  (16) For projection of leading winding side to lagging winding side: 31.055.0,2 jI Sq  (17) 31.233.1,2 jI Sd  (18) 62.278.0,2 jI SNET  (19) Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2464-2470 2467 www.etasr.com Singh and Verma: Performance Analysis and Comparison of Symmetrical and Asymmetrical Dual Stator … Comparing (10) with (16), when projection of lagging winding current component is added with current component of leading winding, it is observed that even for more real power generation the asymmetrical DSIG demands nearly the same reactive power demanded by the symmetrical DSIG, for lesser active power generation. The same is evident from comparison of Figures 3(c) and 3(i) and collaborated by comparing (13) with (19). Asymmetrical DSIG (a) Phasor of Id and Iq for both leading (I1) and lagging (I2)winding currents (b) Projection of Id and Iq of lagging (I2) winding current on leading winding current (I1) (c) Resultant Id and Iq of leading (I1) winding current (d) Phasor of Id and Iq for both leading (I1) and lagging (I2)winding currents (e) Projection of Id and Iq of leading (I1) winding current on lagging winding current (I2) (f) Resultant Id and Iq of lagging (I2) winding current Symmetrical DSIG (g) Phasor of Id and Iq for both leading (I1) and lagging (I2)winding currents (h) Projection of Id and Iq of lagging (I2) winding current on leading winding current (I1) (i) Resultant Id and Iq of leading (I1) winding current (j) Phasor of Id and Iq for both leading (I1) and lagging (I2)winding currents (k) Projection of Id and Iq of leading (I2) winding current on lagging winding current (I1) (l) Resultant Id and Iq of leading (I1) winding current Fig. 3. DSIG phasor diagram: (a)-(f) Asymmetrical DSIG: Projection of Lagging to Leading, (g)-(l) Symmetrical DSIG: Projection of Lagging to Leading Winding Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2464-2470 2468 www.etasr.com Singh and Verma: Performance Analysis and Comparison of Symmetrical and Asymmetrical Dual Stator … V. MATLAB SIMULATION The type of DSIG under consideration is 5 HP, 4 pole, 50 Hz, 415 V. Leading winding side of DSIG is connected to grid with respective phase shift of 30˚ for asymmetrical and 60˚ for symmetrical winding type DSIG. Similarly lagging winding side of DSIG is connected to grid without any phase shift. Both the grid connections are identical except the difference in their phase for connection to leading side. Two distinct equal value three phase capacitor banks are connected to both the winding sets (leading and lagging) of DSIG to partially support the excitation and to filter out any harmonics arriving from load towards generators. Different types of passive loadings (resistive and resistive-inductive loads) are considered connected to both sides with details given in Table I. TABLE I. RATING OF DIFFERENT LOADS Type of Load Value of Load Resistive R= 180 Ω Resistive Inductive Z= 180 Ω (0.8 pf) VI. PERFORMANCE EVALUATION The performance of weak grid coupled DSIGs with symmetrical and asymmetrical stator windings is evaluated comparatively on the basis of voltage regulation and operation under low voltage which in turn is based on establishment of air gap flux. Since the saturation of flux poses as an impediment in operation of IG a careful study of flux in air gap is requisite. The simulated results under MATLAB environment are evaluated for different operating conditions. The connected load is perturbed intentionally to observe and compare the performance of the afore said two types of DSIGs in respect of direct and quadrature axis flux and terminal voltage on the point of common coupling. The utilization of the core and transient characteristic is also been compared by analyzing the direct and quadrature axis flux under steady state and transient conditions during the occurrence of single phase LG fault. Perturbing loads are applied to the terminals of both winding sets of DSIG and their relative effect on their voltage profile and flux are observed and evaluated. To study the effect of cross-coupling, the IGs are loaded in three steps. Initially from t = 0 to t = 0.2s both DSIGs are at no-load. It may be observed from Figure 5 that equal terminal voltages of 340 V are present. At t = 0.2s when both winding sets of the DSIGs are equally loaded with 1.75 kW upf load on each winding set it resulted in the voltage fall from 340 V to 330 V for both symmetrical and asymmetrical DSIGs. Further at t = 0.3s when the loading on leading winding side of both the IGs is increased by 750W each, whereas on lagging side load is kept as it is, a fall in voltage on the leading winding side of the IGs is observed, which varied from 330V to 322V. But as the cross coupling effect which is more prominent for the asymmetrical DSIG the voltage at the lagging winding side increases to 335V, whereas, in the symmetrical DSIG where the effect is not prominent the terminal voltage remain almost same at 330V. Again at t = 0.4s, a load is applied at the lagging winding of both symmetrical as well as asymmetrical DSIG the additional load applied is the preceding perturbation on leading side withdrawn. The terminal voltage of lagging winding side is decreased from 331V to 322V for symmetrical DSIG and from 335V to 328V for asymmetrical DSIG. The effect of cross-coupling in asymmetrical DSIG is once again clearly visible with rise of terminal voltage on the leading winding side i.e. from 322V to 335V. On the other hand, the terminal voltage of leading winding for symmetrical DSIG shows small effect of cross-coupling with a rise of terminal voltage (from 322V to 330V). Fig. 4. MATLAB / Simulink Model of the Proposed System Detailed comparison between symmetrical and asymmetrical DSIG is also carried out by analysis of the flux distribution. It may be observed in Figure 6(a), that the core utilization is very well justified in asymmetrical DSIG where both direct and quadrature axis flux are nearly equal giving a circular shape of d-q flux graph. Whereas, in Figure 6(b) it is evident that the d axis and q axis flux is highly skewed in symmetrical DSIG and the fluxes are not equal, giving it an elliptical shape on the d-q axis indicating its vulnerability to core saturation. Prior to loading the quadrature axis flux values of symmetrical DSIG ranges from -1 to 1 p.u. and -0.75 to 0.75 p.u. for asymmetrical DSIG and, the direct axis flux value ranges from -0.25 to 0.25 p.u. and 0.75 to 0.75 p.u. for symmetrical and asymmetrical DSIG respectively. It may be observed from Figure 6 that after resistive loading on lagging winding side the d-q plot of flux gets slightly compressed vertically and slightly bulges horizontally for both symmetrical and asymmetrical DSIG in the same proportion. The similar effect is observed for both IGs for the loading changes on winding side. To examine the fault ride through capability of DSIG a line to ground fault is created for ten (10) cycles and d-q axis fluxes are plotted in Figure 7 before fault, during fault and during recovery period. Figure 7 shows the plot for direct and quadrature axis flux during pre-fault period for both symmetrical and asymmetrical DSIG. It may be observed from Figure 7 that asymmetrical DSIG shows better stability during fault, which helps during recovery period after fault is removed, compared to symmetrical DSIG. During pre fault conditions the direct and quadrature axis flux values for symmetrical DSIG are observed as (-0.35 to 0.35) and (-0.95 to 0.95) respectively, whereas, for asymmetrical DSIG they vary Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2464-2470 2469 www.etasr.com Singh and Verma: Performance Analysis and Comparison of Symmetrical and Asymmetrical Dual Stator … between (-0.8 to 0.8) and (-0.6 to 0.6) respectively. Quadrature flux of symmetrical DSIG shows more drastic changes when compared with asymmetrical DSIG, whereas, the direct axis flux show relatively lesser decrementing tendency than for the quadrature axis. The symmetrical DSIG witness a decrease from 0.35 to 0.15 p.u. (55% drop) and asymmetrical DSIG witness the decrease from 0.8 to 0.45 p.u. (43% drop) in the direct axis flux, whereas, quadrature axis see a decrease from 0.9 to 0.75 p.u. (15% drop) for symmetrical DSIG and almost no change in quadrature axis flux for asymmetrical DSIG. During recovery period also the tendency of asymmetrical DSIG to get back to normal position is fast as the growth in quadrature axis flux is more rapid than in case of symmetrical DSIG. The symmetrical DSIG shows a slower response compared to asymmetrical DSIG during recovery. The results of the comprehensive analysis are presented in Table II, which depicts the relative merits of asymmetrical DSIG over the symmetrical DSIG for different operating conditions including LVRT and fault conditions. 320 V 320 V 320 V 0.15 0.2 0.25 0.3 0.35 0.4 0.45 320 V Time (s) Loading on both windings Loading on Leading winding Loading on Lagging winding Fig. 5. Terminal Voltages of Symmetrical and Asymmetrical DSIG Fig. 6. Effect of Loading on Air Gap Flux TABLE II. COMPARATIVE CHART OF SYMMETRICAL AND ASYMMETRICAL DSIG UNDER DIFFERENT LOADING AND FAULT CONDITIONS Property Voltage (Cross-Coupling) Current THD Fault Ride Through d-q Flux Graph Symmetrical DSIG Less pronounced Higher Poor Asymmetrical Asymmetrical DSIG More pronounced Lower Good Symmetrical Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2464-2470 2470 www.etasr.com Singh and Verma: Performance Analysis and Comparison of Symmetrical and Asymmetrical Dual Stator … Fig. 7. Effect of LG Fault(10Cycles) on Air Gap Flux VII. CONCLUSION In this paper, we presented a comparative evaluation of operating characteristics of symmetrical and asymmetrical DSIG. The mathematical model of both DSIGs has been developed for comparative analysis. A phasor analysis was also drawn in terms of winding currents for symmetrical and asymmetrical type DSIGs. The comparative performance of the IGs was observed on the basis of demonstrated results for perturbing loads through direct and quadrature axis flux of the air gap and dynamics during recovery of voltage. Demonstrated result also revealed the asymmetrical generator harvests, the advantage of cross coupling which enables the regulation of voltage by magnetizing the effect of reduction in flux. Results demonstrated the superiority of asymmetrical DSIG for wind energy applications. Practically provision of a phase shift of 30˚ is easily attained with the help of star delta transformer thus adds to the merits of asymmetrical generator over the symmetrical one. REFERENCES [1] H. Polinder, F. F. A. Van der Pijl, G.-J. de Vilder, P. J. Tavner, “Comparison of Direct-Drive and Geared Generator Concepts for Wind Turbines”, IEEE Transactions on Energy Conversion, Vol. 21, No. 3, pp. 725-73, 2006 [2] K. Protsenko, D. Xu, “Modeling and Control of Brushless Doubly-Fed Induction Generators in Wind Energy Applications”, IEEE Transactions on Power Electronics, Vol. 23, No. 3, pp. 1191-1197, 2008 [3] F. Bu, Y. Hu, W. Huang, S. Zhuang, K. Shi, “Control Strategy and Dynamic Performance of Dual Stator-Winding Induction Generator Variable Frequency AC Generating System With Inductive and Capacitive Loads”, IEEE Transactions on Power Electronics, Vol. 29, No. 4, pp. 1681- 1692, 2014 [4] T. Long, S. Shao, P. Malliband, E. Abdi, R. A. McMahon, “Crowbarless Fault Ride-Through of the Brushless Doubly Fed Induction Generator in a Wind Turbine Under Symmetrical Voltage Dips”, IEEE Transactions on Industrial Electronics, Vol. 60, No. 7, pp. 2833-2841, 2013 [5] S. Tohidi, H. Oraee, M. R. Zolghadri, S. Shao, P. Tavner, “Analysis and Enhancement of Low-Voltage Ride-Through Capability of Brushless Doubly Fed Induction Generator”, IEEE Transactions on Industrial Electronics, Vol. 60, No. 3, pp. 1146-1155, 2013 [6] R. Boraci, M. Babescu, N. Budisan, A. R. Boraci, “Mathematical Model of The Two Orthogonal Three-Phase Windings Stator Generator”, 5th International Symposium on Applied Computational Intelligence and Informatics, Romania, May 28- 29, 2009 [7] S. Shao, T. Long, E. Abdi, Richard A. McMahon, “Dynamic Control of the Brushless Doubly Fed Induction Generator Under Unbalanced Operation”, IEEE Transactions on Industrial Electronics, Vol. 60, No. 6, pp. 2465-2476, 2013 [8] N. Budisan, O. Prostean, R. Boraci, I. Szeidert, V. Muller, “The Dual Induction Generator for Renewable Energy Conversion Systems. Experimental Results, Problems and Solutions”, IEEE International Joint Conferences on Computational Cybernetics and Technical Informatics, Romania, May 27-29, 2010 [9] R. Singh, V. Verma, “Dual Stator Induction Generator for Rural Electrification”, International Journal of Applied Engineering Research, Vol. 12, No. 16, pp. 5865-587, 2017 [10] F. Bu, L. Zhu, W. Huang, H. Xu, H. Liu, Y. Gao, "Control strategy of AC&DC hybrid generating system based on dual stator-winding induction generator for micro-grid application”, IEEE Transportation Electrification Conference and Expo, pp. 013-017, 2016