Diesel engine generators are widely used in the world, especially in remote site power systems as distributed generators. A weak distribution feeder with a small diesel engine may suffer from voltage and power fluctuations due to misfiring of the engine cylinder. In this study, new generator model with example engine torque was developed for the electromagnetic transient analysis program for power systems named XTAP. The configuration and verification results of the developed model are presented in the paper. The model is considered to be useful for analyses of small power systems with those diesel engines.
Diesel engine generators are widely used in the world, especially in remote site power systems as distributed generators. A weak distribution feeder with a small diesel engine may suffer from voltage and power fluctuations due to misfiring of the engine cylinder. In order to investigate effects of these phenomena to electric power quality in conditions similar to actual power systems, electromagnetic transient (EMT) analysis is very useful.
An EMT analysis program for power systems named XTAP (eXpandable Transient Analysis Program) was developed by Central Research Institute of Electric Power Industry and is widely used in Japan. The two-stage diagonally implicit Runge-Kutta (2S-DIRK) method is adopted for numerical integration in XTAP
. Convergence performance was also improved
and initializations of whole simulated system are available
, therefore, XTAP can improve numerical stability, convergence performance, and expandability of new models for transient analyses of power systems.
XTAP is concidered to be useful to simulate those small power systems with existing problems. In this study, a new diesel generator model with fluctuating engine torque was developed using XTAP. The configuration and verification results of the model are described in this paper.
2. Developed Diesel Generator Model
- 2.1 Overall structure
The target of this study is small-scale synchronous diesel engine generators to supply local electric power systems.
The newly developed diesel generator model for XTAP is composed in the form of Voltage-Behind-Reactance (VBR) model.
The sub-transient inductances for
-axis are assumed to be non-salient here for simplicity, because the purpose of this study is basic consideration of fluctuating mechanical torque of diesel generator, not of detailed modeling of synchronous generator.
The generator dynamics can be represented by the voltage source behind the sub-transient inductance
when sub-transient inductances have no saliency. The
- axis magnetic flux which make the internal voltage behind
have only rotor magnetic flux component as described below.
In the developed model, the voltage source
is connected to the system through the sub-transient inductance
and the armature resistance
phase domain. Since diesel generators are salient machine for transient time domain, single
-axis and single
-axis damper winding were into consideration. The magnetic saturation characteristics were also included in the model.
- 2.2 Equations for generator dynamics
The equations used in the developed model are given in this section. The variables for the armature circuit are calculated in SI values. For the other variables, per unit representations are used. The base values are the generator rated capacity, the rated terminal voltage, and the rated rotational speed.
Since the target of this study is generators with non-salient sub-transient inductances as mentioned above, we can introduce the sub-transient inductance
(p.u.H) as follows:
-axis sub-transient inductance (p.u.H).
Stator voltage equations in
phase are given by
is the armature winding resistance (Ω);
are the generator terminal voltage in
are the armature current of each phase (A). Symbol
denotes time derivative operator as shown in
in (2) are the generator internal induced voltage behind
(V) and given by
0 transformation matrix;
are the generator internal induced voltage behind
are obtained from
is the rotor angular velocity (p.u.rad/s).
are the magnetic flux linkage behind
(p.u.Wb) and assumed here as
-axis armature current (p.u.A).
Per unit air-gap torque
(p.u.Nm) is expressed by
The swing equations in per unit are given by
is the rotor angular velocity deviation (p.u.),
, the mechanical torque input (p.u.Nm),
, the damping coefficient (p.u.),
, the rotor angle (rad),
, the synchronous rotor angular velocity (rad/s),
, the phase angle (rad), and
, time (s).
- 2.3 Magnetic saturation characteristics
The magnetic saturation characteristics are included in the developed model. The no-load magnetic saturation characteristic curve shown in
was adopted to represent the example characteristics and used with a liner interpolation.
No-load saturation characteristic curve.
No-load saturation characteristic curve.
By introducing the
-axis mutual inductance current
(p.u.A) shown in (14) and (15), the unsaturated values of the
-axis air-gap magnetic fluxes
(p.u.Wb) are written as (16) and (17);
The base values for the per unit system for the air-gap fluxes is the saturated value of the air-gap flux when the machine operates with no-load at the rated rotational speed and the rated generator terminal voltage.
The magnetomotive force
(p.u.A/m) induced by the current
are given by
The internal induced voltage
(p.u.V) considering the magnetic saturation can be obtained by referring the no-load magnetic saturation curve. Then the coefficients involving the saturation can be calculated by (19) and (20).
The saturated value of the
-axis mutual inductance
(p.u.H) and the
-axis air-gap magnetic flux
(p.u.Wb) are derived by (21) and (22), respectively.
In the developed model, it is assumed that the
-axis air-gap magnetic flux decrease in proportional to the saturation coefficient
compared to the unsaturated value. The second terms on the rightest sides of (21) and (22) represent the decrease amount due to the magnetic saturation.
- 2.4 Fluctuating mechanical torque of diesel engine
Waveform of mechanical torque of diesel engine generator and its possibility for voltage flicker was thoroughly investigated in
. It concluded that the misfiring of the diesel engine has caused the flicker in that case study.
Diesel engine is a kind of four-cycle internal combustion engines and the number of cylinders firing in each revolution
of the shaft is given by
is the number of cylinders. As for detail of the deriving the relation, refer to
Among the common examples of the number of the poles, 6 and 8 are mentioned as the most sensitive frequencies for human perception in
. A diesel generator with 6 poles and 12 cylinders is considered in this study.
Some example waveforms of the engine torque are also available on
. The normalized engine torque for two revolutions shown in the Fig. 7 and the Table A.1, as well as the case engine torque with #7 misfiring in the Fig. 9 is used in the next section. The data of the Fig. 9 was derived by assuming the net torque of the misfiring cylinder to be zero; therefore, it is very large disturbance to the generator
3. Simulation Results
The characteristics of the developed model were verified through some numerical simulations by use of XTAP.
- 3.1 Simulation conditions
The simulated system is shown in
. A synchronous diesel generator model was connected to an infinite bus through a step-up transformer and long transmission lines. The step-up transformer was simulated by inductances, the infinite bus was by 3-phase sinusoidal voltage source. The mutual capacitances of transmission lines were neglected. The initial power flow condition and the per unit impedance are described in
. The reference values of line voltage, capacity and frequency are 220 V, 100 kVA, and 60 Hz, respectively.
Schematic diagram of the simulated power system
Generator constants are listed in
. The base value of per unit system is the rated value of the generator: 220 V, 100 kVA, and 60 Hz. Since the one damper winding is considered in each of
-axis circuit as mentioned above, there are no
-axis transient constants as
The generator was equipped with AVR and speed governor (GOV).
exhibits the block diagram of the AVR. It is one of the standard AVR models used for the power system analysis in Japan. The characteristic of the GOV was simply expressed by (24), where
is the initial value of the mechanical torque.
is the main body of the engine torque. The oscillation due to firing of cylinders was simulated as coefficients times to
Block diagram of the AVR used in the simulation
It should be noted that there are many variations in characteristics of actual AVRs and GOVs for diesel engine generators. The individual controller model should be used when this diesel generator model is used for analyses.
- 3.2 Simulation results of fluctuating engine torque
shows the simulated engine torque. The waveform was composed by combination of sinusoidal waves.
Simulation results of engine torque
Fig. 3 (a)
exhibits the composed engine torque accords with the data in
Fig. 3 (b)
shows the engine torque with misfiring of Cylinder #7.
exhibits the simulation results with engine torque oscillation.
Simulation results for the torque variation
Fig. 4 (a)
shows the waveform of the simulated torque. The large sag in the torque is caused by the misfiring of cylinder #7.
Fig. 4 (b)
describes the resulting large disturbance in the generator active power output. The result indicates that the misfiring of cylinders may cause not only voltage flickers but also change in active power output and might leads to stability problem in a remote weak power system.
Fig. 4 (c)
shows the generator terminal voltage. Before the misfiring, small oscillations due to normal mechanical torque fluctuation was successfully simulated. The voltage dropped by the misfiring and largely changed due to the active power imbalance.
Fig. 4 (d)
, the rotor angular velocity variation and the oscillation due to the torque variation was observed. To simulate these oscillations in the shafts is also important to develop generator control systems such as GOVs and excitation control systems with the consideration of shaft noise.
- 3.3 Simulation results of generator variables
exhibits the simulation results of generator variables for 1 to 5 seconds to the engine torque oscillation including cylinder #7’s misfiring.
Simulation results of generator variables.
Fig. 5 (a)
shows the waveform of the applied fluctuating engine torque. The misfiring causes disturbances to the power system as seen in the figures.
The effect of magnetic saturation was checked through the numerical simulations. The red lines in
Fig. 5 (b)
show the results considering the magnetic saturation, and the blue lines, the results without considering saturation.
Among these simulations, the same reference values were given for the generator active power of
=0.9 (p.u.W) and the generator terminal voltage
=1.01 (p.u.V) as seen in
Fig. 5 (b)
There are seen small differences between with and without saturation in the generator reactive power (
Fig. 5 (d)
) and the phase angle (
Fig. 5 (e)
). The difference due to magnetic saturation in the internal electromotive force leaded to these results. In this simulation, the unsaturated case was simulated only by letting the saturation coefficient
=1.0 in order to check the results of applying the coefficient
. Therefore, when the saturation is considered, the phase angle became smaller and the system become slightly stable compared to the unsaturated case. The results will vary if the unsaturated case is simulated by using the saturated values of generator constants.
On the other hand, little difference are seen in the variables involving the balance of the input and output torques such as the active power (
Fig. 5 (b)
) and the rotor angular velocity (
Fig. 5 (f)
) for this case.
Through these numerical simulations, it is confirmed that the developed diesel generator model can simulate fluctuating engine torque and magnetic saturation, and contribute to EMT simulations to asses effects of these phenomena to electric power quality. Investigation of these assessments in small remote sites or islanded power systems with several diesel engine generators using the developed model is a future work.
In this study, the new diesel generator model with example fluctuating engine torque was developed for transient analyses using XTAP. The magnetic saturation characteristics were also included.
It was verified through the numerical simulations and comparison with the previous study that the model can simulate fluctuations in its terminal voltage and active power output due to the variation of engine torque.
The developed model is considered to be useful for analyses of small power systems with diesel generators. Effects of various fluctuating engine torque to electric power quality in remote sites or small islanded power systems should be investigated as a future work.
This work was partly supported by research grant from Japan Power Academy.
Orie Sakamoto She received her B.S., M.S., and Ph.D degree from the University of Tokyo, Tokyo, Japan, in 2001, 2003, and 2006, respectively. Her areas of interest include modeling and control of synchronous generators and power systems including solar and wind power generations. She is a member of IEEJ and IEEE.
“Numerical Integration by the 2-Stage Diagonally Implicit Runge-Kutta Method for Electromagnetic Transient Simulations”
IEEE Trans. Power Delivery
DOI : 10.1109/TPWRD.2008.923397
“A Robust and Efficient Iterative Scheme for the EMT Simulations of Non-linear Circuits”
IEEE Trans. Power Delivery
“A Practical Steady-State Initialization Method for Electromagnetic Transient Simulations”
Proc. International Conference of Power Systems Transients (IPST2011)
Delft, the Netherlands
June 14-17, 2011
Anderson P. M.
“Analysis of a diesel-engine driven generating unit and the possibility for voltage flicker”
IEEE Trans. Energy Conversion
DOI : 10.1109/60.372566