Virtual Laboratory of Unbalanced Transient Condition in Synchronous Generator

The electrical engineering department at the Sekolah Tinggi Teknologi Nasional (STTNAS), Yogyakarta has recently reconnoitered virtual laboratories for its undergraduate synchronous generator course to complement existing full-scale laboratory equipment. This study explores virtual laboratory development to be treated as an accessorial tool for enhancing instruction. The focus of this synchronous generator course is the dynamictransient behavior of the system after small disturbances as affected by the unbalanced load. The work is mainly carried out through nonlinear simulations under Matlab-Simulink. Results of the first version of the synchronous generator virtual laboratory and details of its development are provided.


Introduction
In Indonesia, significant percentages of the undergraduate students in electrical engineering are working in industry.To meet with thestudents needto expert conception while experimentations support the understanding of the topics, written exercises are required for those students.Real experiments are keys to develop skills dealling with physical processes [1,2].The laboratory education in the case of synchronous generator course is difficult, especially on the specific condition such as transient condition, at the undergraduate level because modernization of machine laboratories need high investment cost.As an alternative, virtual Laboratory can be treated as an accessorial tool of real laboratory to enhance instruction for conventional on-campus students; in which students can enable to improve the skills before going to the actual laboratory, to learn breaking the restriction of ordinary arrangement and enhance the instruction [3,4].
The studies of dynamic synchronous generator date back from the time when the first power systems began to operate, and to start its interconnection process.The excessive advances regarding to such topics have been made since then, with the purpose ofdescribing in detail the synchronous generator dynamics as affected by small disturbances, where the system in ordinarily unbalanced under steady-state conditions.From the transient stability perspective, several studies have already been accomplished during the 1980s in view of the unbalanced system in the dynamic of asynchronous generator following a system disturbance [5][6][7].This transient phenomenon of a power system utilizing synchronous generator as a source of electrical energy is usually analyzed by employed single machine infinite bus takento a balanced power system [8,9].The results are frequency (f) and damping ratio () of the electromechanical model calculated in each studied case of unbalanced scenario.
Relatively little concern has been paid to the dynamic behavior of an unbalanced system considered to small-signal behavior.The model solution is relatively difficult considering analysis given starts with equations containing thephase quantities and transforms them into new equations using α, β components.They are linear differential equations but the coefficients are variable.In [10], the model developments usethe linearized state-space model and are relatedto DG facilities connected to distribution grids.Several simplifications areapplied to this analysis, such as the electric power systemis considered as a balanced three-phase system so it allows being representedby its single-phase equivalent.For all that, it is said that the modeling approach presented here can be applied to any power system, regardless of size and/or voltage level.
Synchronous machines are usually employed as generators connected to a power system.In the current IEEE Standard [11], it is acknowledged that synchronous machines may be accurately modeled by two lumped-parameter equivalent circuits representing the q and the daxis.The number ofrotor damper branches are selected in accordance with the rotor design.In low-order models, these branches correspond to the actual amortisseur windings; higher order models utilize these branches to represent the distributed effects in the rotor iron.To minimize the effort spent on setting up the model, it also takes advantage of simplifications that are offered by the separation in time scales of the different dynamic behaviors and the fact that the severity of a disturbance is usually attenuated as it propagates through the system.For example, the duration of the electrical transient of the network is very short relative to the electromechanical dynamic of the generator; as such, a static representation of the network can be used where longer electromechanical oscillations are primarily of interest.On the other side, the duration of interest may not be long enough to require the inclusion of the dynamic of slower acting component, such as automatic generation control.The model chosen for the generators on the network need not the same; it only has to be compatible with the networks representation in the context of the solution algorithm used.If the type and location of the disturbances are to be known, one can selectively employs a detailed model for generators electrically close to the fault locations, and the choose simpler models for generators further away.For example in [12], the phenomenon of interest is the transients in machines.The electromechanical transient occurs mainly in machines, ignoring the faster transients in lines and load reactances.This simplification is acceptable; this has an advantage of computational efficiency.
With the intention of helping to increase the knowledge about the dynamic behavior of synchronous generator after small disturbances, this study presents a thorough study concerning the dynamic behavior of unbalanced transient state synchronous generator model.The work is mainly carried out through nonlinear simulations under Matlab-Simulink [13][14][15][16][17][18], where the unbalanced power systems are calculated through the EDSA 2000 [19].
This work is organized as follows.A research method is presented on Section 2. Section 3 presents research and analysis, whereas the conclusion followed by the references is presented on Section 4.

Research Method
The system of this study is the single machine infinite bus (SMIB) shown Figure 1.The external line parameters, and , are to be varied to change the electrical strength of the connection between generator and infinite bus.In this case, the infinite bus supposedly represents a large system to which the generator is connected [20,21].In Figure 1, there is only single machine and a simple external network.So, it can easily transform variables of external network to the rotor reference frame of that single machine.If voltage drop across external line is: Then the stator voltage equations with the external impedance included would become where , , , are stator voltages and stator currents of both q-axis and d-axis, respectively.Using transient model, Equation ( 2) and (3) will be changed and the parameters of synchronous generator model will neglect the changes in stator qd0 flux linkages, shown in Table 1.All of objects in Table 1 are the forming components of synchronous generator model.Figure 2 shows unbalanced three-phase synchronous generator that connected to infinite bus using series RL line coupling.The voltage components of infinite bus in the form of An inside the model box of three-phase synchronous generator of Figure 2 is shown in Figure 3.While the inside model box of stator winding is shown in Figure 4.The stator qd0 currents are determined usingthe stator voltage equations that included the external series RL line parameters in series with stator resistance and leakage reactance.And two rotor circuits are represented bythe differential equation in and .Only the physical field winding on teh daxis has an external excitation of ; that for the q-axis, is zero.The dynamic characteristic of SMIB is represented using the voltage that lagging by transient impedance , shown in Figure 5.And the dynamic equation of rotor per unit can be described using the relationship between rotor angle δ and rotor speed , where , , and are mechanic, electric, and damping torques, respectively; while and H are base rotor speed and inertia constant.To change into ( ) and put small disturbance into working point variable and rotor angle into position will yield, where is synchronization torque constant.Applying Laplace transform into Eq.( 5) will result The value of is inversly proportional with total reactance between and .One can access Matlab's GUI facilities to construct a software package of virtual laboratory for studying synchronous generator under unbalanced transient-state condition.As anexample of using Matlab's GUI capabilities, menu and plotting commands are implemented in a script file to provide interactive windows (Figure 6).The main menu, which is displayed after running the file, are shown in Figure 7 and Figure 8.

Results and Analysis
The studied generator is Grati which is one of generating plants of the 500 kV EHV Java-Madura-Bali (Jamali) System shown in Figure 9.The grid consists of 9 generator nodes and 21 load nodes.The Paiton's bus is slack-bus and others are PV buses.The base system capacity is 100,000 MVA [22].Generator ratings and parameters are shown in Table 2.
The simulation of the proposed generator model is carriedout by Matlab.As inputs of this generator model are stator voltages which are derived by analyzing a single-linediagram of 500 kV EHV Jamali System on EDSA 2000 Newton-Raphson (NR) loads flowsoftware.The process of numerical simulation methodcan is presented by the block diagram of Figure 10.
Using EDSA 2000 software program one can get the load flow calculation results from Figure 9. Table 2 and 3    This study is carried out utilizing the created GUI windows.Running the created GUI Mfile, called "AwalGenMod", from Matlab workspace will display the main window shown in Figure 7.The window presented in Figure 8 will appear after clicking on the icon named START.Setting the slider icons of generator's parameters and typical stator voltages and also clicking on the icon named SIMULATE will present the result menu shown in Figure 7.This leads to the following figures which present the real and reactive powers, the current and the voltage, power angle, and electric torque.
In Figure 11, it can see when the synchronous generator connects and delievers energy to the 500 kV EHV Jamali System as much as ( ) p.u, its current in phase could become unbalanced even though all of the grid loads are balanced.Moreover, the generated active and reactive powers will become slightly to oscillate but the rotor speed tends to be constant.And the power angle peak goes up to 1.04 p.u.
As the level of unbalance is increasing up to 7.5%, the currents in phase are still unbalanced.Comparing among figures observably that the increasing of unbalanced percentage of load has an affect on the decreasing and oscillating of the generated active and reactive powers, the power angle and the magnitude of stator current oscillation, and the electric torque.On the contrary, it has no effect on the speed of the machine's rotor; the rotor speed is always constant.

Conclusion
A useful approach for studying unbalanced transient state operation of synchronous generator under has been presented in this paper.Two operation conditions of the synchronous generator, balanced load and 7.5% of unbalanced load, are mathematically modeled then simulated using Matlab.The simulation results state that the increasing level of unbalanced load has significant influence on the parameters ofgenerator dynamic regarding decreasing and oscillating of its magnitudes, except to the rotor speed.The proposed tool is made easy to use by providing anactive link with the simulated models using some of GUI functions.The given

IJEECS
ISSN: 2502-4752  Virtual Laboratory of Unbalanced Transient Condition in Synchronous Generator (Sugiarto K.) 3 qd0 refrence framework are transformed into rotor reference framework using abc to qd0 converter.

Figure 1 .
Figure 1.Single Machine Infinite Bus System

 4 Figure 2 . 14 ]Figure 3 . 5 Figure 4 .Figure 5 .
Figure 2. Unbalanced Transient Generator Model using SMIB with and [14] the normal quadratic square equation, , that has the root values of ( √( ) ) so the equation of undamped frequency will have the values of

Figure 6 .
Figure 6.Architecture of Unbalanced Transient Figure 7.The Main Window of the Developed State Generator Simulator Software Tool

Figure 8 .
Figure 8.The Window of Inserting the Inputs and Displaying the Variable Tests

7
presents inter-phase voltage values of the test generator terminal, before and after loading condition.It is shown that under unbalanced loads condition, the phase angles of terminal generator voltage are deviated from its balanced value.

detik) Tem (p.u) Torsi Listrik Instan
helpfulness of the proposed tool for studying the synchronous generator dynamic under unbalanced transient state condition.