A Diesel Generator Model with Fluctuating Engine Torque Including Magnetic Saturation for Transient Analysis using XTAP

Journal of Electrical Engineering and Technology.
2015.
May,
10(3):
1298-1303

This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/)which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

- Received : November 17, 2014
- Accepted : December 26, 2014
- Published : May 01, 2015

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

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.
d
- and
q
-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
L''
when sub-transient inductances have no saliency. The
d
- and
q
- axis magnetic flux which make the internal voltage behind
L''
have only rotor magnetic flux component as described below.
In the developed model, the voltage source
e''
is connected to the system through the sub-transient inductance
L''
and the armature resistance
R_{a}
in
abc
phase domain. Since diesel generators are salient machine for transient time domain, single
d
-axis and single
q
-axis damper winding were into consideration. The magnetic saturation characteristics were also included in the model.
L''
(p.u.H) as follows:
where
L''_{d}
and
L''_{q}
are the
d
- and
q
-axis sub-transient inductance (p.u.H).
Stator voltage equations in
abc
phase are given by
where
R_{a}
is the armature winding resistance (Ω);
e_{a}
,
e_{b}
, and
e_{c}
are the generator terminal voltage in
abc
phase (V);
i_{a}
,
i_{b}
, and
i_{c}
are the armature current of each phase (A). Symbol
p
denotes time derivative operator as shown in
Variables
e''_{a}
,
e''_{b}
, and
e''_{c}
in (2) are the generator internal induced voltage behind
L''
(V) and given by
Matrix
C
_{abc−dq0}
is the
abc-dq
0 transformation matrix;
e''_{d}
and
e''_{q}
are the generator internal induced voltage behind
L''_{d}
and
L''_{q}
(p.u.V).
e''_{d}
and
e''_{q}
are obtained from
where
ω_{r}
is the rotor angular velocity (p.u.rad/s).
Flux linkage
ϕ_{fkd}
and
ϕ_{fkq}
are the magnetic flux linkage behind
L''_{d}
and
L''_{q}
(p.u.Wb) and assumed here as
where
i_{d}
and
i_{q}
are the
d
- and
q
-axis armature current (p.u.A).
Per unit air-gap torque
T_{e}
(p.u.Nm) is expressed by
The swing equations in per unit are given by
where
s_{g}
is the rotor angular velocity deviation (p.u.),
T_{m}
, the mechanical torque input (p.u.Nm),
D
, the damping coefficient (p.u.),
δ
, the rotor angle (rad),
ω
_{0}
, the synchronous rotor angular velocity (rad/s),
θ
, the phase angle (rad), and
t
, time (s).
No-load saturation characteristic curve.
By introducing the
d
- and
q
-axis mutual inductance current
i_{md}
(p.u.A) and
i_{mq}
(p.u.A) shown in (14) and (15), the unsaturated values of the
d
- and
q
-axis air-gap magnetic fluxes
ϕ_{adu}
and
ϕ_{aqu}
(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
H_{g}
(p.u.A/m) induced by the current
i_{md}
and
i_{mq}
are given by
The internal induced voltage
e_{gs}
(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
d
-axis mutual inductance
L_{mds}
(p.u.H) and the
d
-axis air-gap magnetic flux
ϕ_{ads}
(p.u.Wb) are derived by (21) and (22), respectively.
In the developed model, it is assumed that the
d
- and
q
-axis air-gap magnetic flux decrease in proportional to the saturation coefficient
K
_{s1}
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.
N_{F}
of the shaft is given by
where
N_{C}
is the number of cylinders. As for detail of the deriving the relation, refer to
[4]
and
[5]
.
Among the common examples of the number of the poles, 6 and 8 are mentioned as the most sensitive frequencies for human perception in
[4]
. 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
[4]
. 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
Schematic diagram of the simulated power system
Generator constants are listed in
Table 2
. 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
d
- and
q
-axis circuit as mentioned above, there are no
q
-axis transient constants as
L_{q}
’ and
T_{q}
’.
Generator Constants
The generator was equipped with AVR and speed governor (GOV).
Fig. 2
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
T
_{m0}
is the initial value of the mechanical torque.
T_{m}
is the main body of the engine torque. The oscillation due to firing of cylinders was simulated as coefficients times to
T_{m}
.
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.
Simulation results of engine torque
Fig. 3 (a)
exhibits the composed engine torque accords with the data in
[4]
.
Fig. 3 (b)
shows the engine torque with misfiring of Cylinder #7.
Fig. 4
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.
In
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.
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)
to
(f)
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
P_{g}
=0.9 (p.u.W) and the generator terminal voltage
V_{t}
=1.01 (p.u.V) as seen in
Fig. 5 (b)
and
5 (c)
, respectively.
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
K
_{s1}
=1.0 in order to check the results of applying the coefficient
K
_{s1}
and
K
_{s2}
. 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.
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.

Diesel engine generator
;
Electromagnetic transient analysis
;
Power systems
;
Synchronous generator
;
XTAP

1. Introduction

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
[1]
. Convergence performance was also improved
[2]
and initializations of whole simulated system are available
[3]
, 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
- 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
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.3 Magnetic saturation characteristics

The magnetic saturation characteristics are included in the developed model. The no-load magnetic saturation characteristic curve shown in
Table 1
was adopted to represent the example characteristics and used with a liner interpolation.
No-load saturation characteristic curve.

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 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
[4]
. 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
PPT Slide

Lager Image

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
Fig. 1
. 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
Fig. 1
. The reference values of line voltage, capacity and frequency are 220 V, 100 kVA, and 60 Hz, respectively.
PPT Slide

Lager Image

Generator Constants

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 3.2 Simulation results of fluctuating engine torque

Fig. 3
shows the simulated engine torque. The waveform was composed by combination of sinusoidal waves.
PPT Slide

Lager Image

PPT Slide

Lager Image

- 3.3 Simulation results of generator variables

Fig. 5
exhibits the simulation results of generator variables for 1 to 5 seconds to the engine torque oscillation including cylinder #7’s misfiring.
PPT Slide

Lager Image

4. Conclusion

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.
Acknowledgements

This work was partly supported by research grant from Japan Power Academy.

BIO

Noda T.
,
Takenaka T.
,
Inoue T.
2009
“Numerical Integration by the 2-Stage Diagonally Implicit Runge-Kutta Method for Electromagnetic Transient Simulations”
IEEE Trans. Power Delivery
24
(1)
390 -
399
** DOI : 10.1109/TPWRD.2008.923397**

Noda T.
,
Kikuma T.
2011
“A Robust and Efficient Iterative Scheme for the EMT Simulations of Non-linear Circuits”
IEEE Trans. Power Delivery
26
(2)
1030 -
1038

Noda T.
,
Takenaka K.
“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
Available:

Anderson P. M.
,
Mirheydar M.
1995
“Analysis of a diesel-engine driven generating unit and the possibility for voltage flicker”
IEEE Trans. Energy Conversion
10
(1)
37 -
47
** DOI : 10.1109/60.372566**

Boldea I.
2005
“Synchronous Generators”
CRC Press
3-12 -
3-15

Citing 'A Diesel Generator Model with Fluctuating Engine Torque Including Magnetic Saturation for Transient Analysis using XTAP
'

@article{ E1EEFQ_2015_v10n3_1298}
,title={A Diesel Generator Model with Fluctuating Engine Torque Including Magnetic Saturation for Transient Analysis using XTAP}
,volume={3}
, url={http://dx.doi.org/10.5370/JEET.2015.10.3.1298}, DOI={10.5370/JEET.2015.10.3.1298}
, number= {3}
, journal={Journal of Electrical Engineering and Technology}
, publisher={The Korean Institute of Electrical Engineers}
, author={Sakamoto, Orie}
, year={2015}
, month={May}