Skip to content

IEEE Governor Models

This page documents the IEEE-standard turbine-governor models available in RAMSES, all using the tor_ prefix. These models implement standardised transfer-function blocks for diesel, gas-turbine, steam-turbine, hydraulic, and simplified ENTSO-E governors.


The tor_DEGOV1 model represents the diesel engine governor as defined in the IEEE Committee Report on governors. It consists of three main stages:

1. Speed governor (electrical control / lead–lag)

The speed error is formed from the measured speed deviation and a reference signal derived from either mechanical power (PP) or mechanical torque (TmT_m), selected by the switch SWMSWM:

V60=(1SWM)Tm+SWMPV_{60} = (1 - SWM) \cdot T_m + SWM \cdot P REF=V60RREF = V_{60} \cdot R

The governor speed error drives a second-order lead–lag controller (transfer function approximating a two-pole, two-zero system) with time constants T1T_1, T2T_2, T3T_3:

Ggov(s)1+sT3(1+sT1)(1+sT2)G_{gov}(s) \approx \frac{1 + s T_3}{(1 + s T_1)(1 + s T_2)}

2. Actuator (fuel injection system)

The governor output passes through a lead–lag block with gain KK and time constants T4T_4, T6T_6:

Gact(s)=K1+sT41+sT6G_{act}(s) = K \cdot \frac{1 + s T_4}{1 + s T_6}

followed by a first-order lag T5T_5, and an integrator with limits [TMIN,TMAX][T_{MIN},\, T_{MAX}].

3. Engine dead time and torque output

The actuator output is subject to engine dead time TDT_D, approximated by a second-order Padé:

esTD10.5sTD+0.0833s2TD21+0.5sTD+0.0833s2TD2e^{-sT_D} \approx \frac{1 - 0.5 s T_D + 0.0833 s^2 T_D^2}{1 + 0.5 s T_D + 0.0833 s^2 T_D^2}

The final mechanical torque is computed as:

Tm=11+sTEVengineT_m = \frac{1}{1 + s T_E} \cdot V_{engine} Pm=TmωP_m = T_m \cdot \omega

The droop characteristic is:

Δω+R(V60REF)=0    error=Δω+RVact\Delta \omega + R \cdot (V_{60} - REF) = 0 \implies \text{error} = \Delta\omega + R \cdot V_{act}
ParameterDescription
SWMInput switch: 0 = mechanical torque (TmT_m), 1 = electrical power (PP)
T1Governor time constant — lead numerator (s)
T2Governor time constant — lag denominator (s)
T3Governor time constant — second-order denominator (s)
KActuator gain
T4Actuator lead time constant (s)
T5First-order actuator lag (s)
T6Actuator lag time constant (s)
TMINMinimum fuel/torque limit (pu)
TMAXMaximum fuel/torque limit (pu)
TDEngine dead time (s)
RSpeed droop (pu/pu)
TEEngine time constant (s)
SYNC_MACH g1 g1 1. 1. 0. 0. 600. 570. 6. 0. 2.05
XT 0.15 2.2 0.3 0.2 2. 0.4 0.2 0.1 6.0257 0. 7.00 0.05 1.5 0.05
EXC GENERIC1 3.0618 -0.1 1. 0. 100. -1. -20.0 10. 120. 5. 12.5 0.1 0. 5.
1 75. 15. 0.22 0.012 0.22 0.012 -0.1 0.1
TOR DEGOV1 0 ! SWM: 0 = torque input, 1 = power input
0.02 ! T1: governor lead time constant (s)
0.02 ! T2: governor lag time constant (s)
0.20 ! T3: governor second-order denominator (s)
1.0 ! K: actuator gain
0.25 ! T4: actuator lead time constant (s)
0.009 ! T5: first-order actuator lag (s)
0.038 ! T6: actuator lag time constant (s)
0.0 ! TMIN: minimum fuel/torque limit (pu)
1.05 ! TMAX: maximum fuel/torque limit (pu)
0.01 ! TD: engine dead time (s)
0.05 ! R: speed droop (pu/pu)
0.05 ; ! TE: engine time constant (s)

The tor_ENTSOE_simp model is a simplified ENTSO-E speed governor suitable for primary frequency response studies. It is derived from the ENTSO-E recommendations for equivalent turbine-governor representation.

The governor equation forms a speed error with permanent droop RR:

Rp˙1C+(ω1)=0,C=TmRR \cdot \dot{p}_1 - C + (\omega - 1) = 0, \quad C = T_m \cdot R

This error drives a first-order lag (representing the valve or steam chest) with limits [VMIN,VMAX][V_{MIN},\, V_{MAX}]:

dp2dt=1T1(p1p2),p2[VMIN,VMAX]\frac{d p_2}{dt} = \frac{1}{T_1}(p_1 - p_2), \quad p_2 \in [V_{MIN},\, V_{MAX}]

The turbine mechanical power is modelled by a lead–lag transfer function with time constants T2T_2 (zero) and T3T_3 (pole):

Pm(s)=1+sT21+sT3p2(s)P_m(s) = \frac{1 + s T_2}{1 + s T_3} \cdot p_2(s)

The final mechanical torque output accounting for shaft speed is:

Tmω=PmT_m \cdot \omega = P_m
ParameterDescription
RPermanent speed droop (pu/pu)
T1Governor/valve time constant (s)
VMINMinimum valve position / lower output limit (pu)
VMAXMaximum valve position / upper output limit (pu)
T2Turbine lead time constant — zero (s)
T3Turbine lag time constant — pole (s)

Internal parameter: C = Tm · R (initialised from the load-flow mechanical torque)

SYNC_MACH g1 g1 1. 1. 0. 0. 400. 360. 6. 0. 2.05
XT 0.15 2.2 0.3 0.2 2. 0.4 0.2 0.1 6.0257 0. 7.00 0.05 1.5 0.05
EXC GENERIC1 3.0618 -0.1 1. 0. 100. -1. -20.0 10. 120. 5. 12.5 0.1 0. 5.
1 75. 15. 0.22 0.012 0.22 0.012 -0.1 0.1
TOR ENTSOE_simp 0.05 ! R: permanent droop (pu/pu)
0.50 ! T1: governor/valve time constant (s)
0.0 ! VMIN: minimum valve position (pu)
1.1 ! VMAX: maximum valve position (pu)
0.0 ! T2: turbine lead time constant (s)
10.0 ; ! T3: turbine lag time constant (s)

The tor_GAST model represents a gas turbine unit, following the IEEE Committee Report (1973) and the NEPLAN/PowerWorld convention. It captures the speed governor, compressor discharge lag, and fuel/combustion dynamics.

Speed governor with droop and load reference:

error=R(VactLR)+(ω1)=0\text{error} = R \cdot (V_{act} - LR) + (\omega - 1) = 0

where LR=TmLR = T_m (load reference, initialised from mechanical torque) and RR is the droop.

The speed error drives a value position through a limiter [VMIN,VMAX][V_{MIN},\, V_{MAX}] and a first-order lag T1T_1:

dp2dt=1T1(p1p2)\frac{d p_2}{dt} = \frac{1}{T_1}(p_1 - p_2)

The compressor discharge is modelled by lag T2T_2:

dp3dt=1T2(p2p3)\frac{d p_3}{dt} = \frac{1}{T_2}(p_2 - p_3)

The mechanical power (including speed-dependent damping DTURBD_{TURB}) is:

Pm=p3+DTURB(ω1)P_m = p_3 + D_{TURB} \cdot (\omega - 1)

The exhaust temperature is tracked by a lag T3T_3:

dp4dt=1T3(p3p4)\frac{d p_4}{dt} = \frac{1}{T_3}(p_3 - p_4)

A temperature limit is applied via:

p5=AT+KT(ATp4)p_5 = A_T + K_T \cdot (A_T - p_4)

A minimum gate blocks the valve position: p1=min(pgov,p5)p_1 = \min(p_{gov},\, p_5).

The torque-power conversion:

Tmω=PmT_m \cdot \omega = P_m
ParameterDescription
RPermanent speed droop (pu/pu)
T1Governor/valve time constant (s)
T2Compressor discharge time constant (s)
T3Radiation shield / exhaust temperature lag (s)
ATAmbient temperature load limit (pu)
KTTemperature control loop gain
VMAXMaximum valve position (pu)
VMINMinimum valve position (pu)
DTURBTurbine damping factor (pu torque / pu speed)

Internal parameter: LR = Tm (load reference, set at initialisation)

SYNC_MACH g1 g1 1. 1. 0. 0. 400. 360. 6. 0. 2.05
XT 0.15 2.2 0.3 0.2 2. 0.4 0.2 0.1 6.0257 0. 7.00 0.05 1.5 0.05
EXC GENERIC1 3.0618 -0.1 1. 0. 100. -1. -20.0 10. 120. 5. 12.5 0.1 0. 5.
1 75. 15. 0.22 0.012 0.22 0.012 -0.1 0.1
TOR GAST 0.05 ! R: permanent droop (pu/pu)
0.40 ! T1: governor/valve time constant (s)
0.10 ! T2: compressor discharge time constant (s)
3.00 ! T3: exhaust temperature lag (s)
1.00 ! AT: ambient temperature load limit (pu)
2.50 ! KT: temperature control loop gain
1.05 ! VMAX: maximum valve position (pu)
0.0 ! VMIN: minimum valve position (pu)
0.0 ; ! DTURB: turbine damping factor

The tor_TGOV1D model represents a simple steam turbine governor, based on the IEEE TGOV1 structure with an added damping term. It is one of the most widely used simple governor models in power-system stability studies.

Speed governor with permanent droop:

The speed error (with droop RR) drives a first-order lag T1T_1 with output limits [VMIN,VMAX][V_{MIN},\, V_{MAX}]:

Rp1C+(ω1)=0,C=TmRR \cdot p_1 - C + (\omega - 1) = 0, \quad C = T_m \cdot R dp2dt=1T1(p1p2),p2[VMIN,VMAX]\frac{d p_2}{dt} = \frac{1}{T_1}(p_1 - p_2), \quad p_2 \in [V_{MIN},\, V_{MAX}]

Turbine with reheater and damping:

The valve position passes through a lead–lag block representing the steam chest and reheater:

Pm(s)=1+sT21+sT3p2(s)P_m(s) = \frac{1 + s T_2}{1 + s T_3} \cdot p_2(s)

The total mechanical power includes a speed-dependent damping term DtD_t:

Pm=pturbine+Dt(ω1)P_m = p_{turbine} + D_t \cdot (\omega - 1)

Hence the mechanical torque:

Tm=PmωT_m = \frac{P_m}{\omega}
ParameterDescription
RPermanent speed droop (pu/pu)
T1Steam chest / valve actuator time constant (s)
VMAXMaximum valve position (pu)
VMINMinimum valve position (pu)
T2Lead time constant — reheater zero (s)
T3Lag time constant — reheater pole (s)
DtTurbine damping coefficient (pu torque / pu speed deviation)

Internal parameter: C = Tm · R (initialised from load-flow mechanical torque)

SYNC_MACH g1 g1 1. 1. 0. 0. 400. 360. 6. 0. 2.05
XT 0.15 2.2 0.3 0.2 2. 0.4 0.2 0.1 6.0257 0. 7.00 0.05 1.5 0.05
EXC GENERIC1 3.0618 -0.1 1. 0. 100. -1. -20.0 10. 120. 5. 12.5 0.1 0. 5.
1 75. 15. 0.22 0.012 0.22 0.012 -0.1 0.1
TOR TGOV1D 0.05 ! R: permanent droop (pu/pu)
0.50 ! T1: steam chest/valve actuator time constant (s)
1.05 ! VMAX: maximum valve position (pu)
0.0 ! VMIN: minimum valve position (pu)
2.10 ! T2: lead time constant — reheater zero (s)
7.00 ! T3: lag time constant — reheater pole (s)
0.0 ; ! Dt: turbine damping coefficient

The tor_hygov model represents a hydraulic turbine and governor following a simplified penstock/water-column formulation, consistent with the IEEE representation for hydro governors (IEEE Std 1110).

Governor with droop and PI controller:

The speed error (with permanent droop RR and transient droop rr) enters a lag TfT_f to produce the error signal ee:

input=1ωR(cP0)\text{input} = 1 - \omega - R \cdot (c - P_0) dedt=1Tf(inpute)\frac{d e}{dt} = \frac{1}{T_f}\left(\text{input} - e\right)

The rate of change of gate demand is computed via the transient droop:

x˙2=erTr\dot{x}_2 = \frac{e}{r \cdot T_r}

This is rate-limited to [VELM,+VELM][-V_{ELM},\, +V_{ELM}], then integrated with gate limits [Gmin,Gmax][G_{min},\, G_{max}]:

dx4dt=x3,x4[Gmin,Gmax]\frac{d x_4}{dt} = x_3, \quad x_4 \in [G_{min},\, G_{max}]

The PI controller output (desired gate) is:

x5=x4+er,c=clip(x5,Gmin,Gmax)x_5 = x_4 + \frac{e}{r}, \quad c = \text{clip}(x_5,\, G_{min},\, G_{max})

The actual gate position follows through the servo:

dgdt=1Tg(cg)\frac{d g}{dt} = \frac{1}{T_g}(c - g)

Penstock / water column:

The head deviation dHdH is determined by the non-linear water flow equation and water-starting time TWT_W:

1(Qg)2=dH1 - \left(\frac{Q}{g}\right)^2 = dH TWdQdt=dHT_W \cdot \frac{dQ}{dt} = dH

Turbine mechanical torque:

Tmω=At(1dH)(QQnl)Dturb(1ω)gT_m \cdot \omega = A_t \cdot (1 - dH)(Q - Q_{nl}) - D_{turb} \cdot (1 - \omega) \cdot g Pm=TmωP_m = T_m \cdot \omega
ParameterDescription
RPermanent droop (pu/pu)
rTransient droop (pu/pu)
TrTransient droop reset time (s)
TfFilter time constant for speed error (s)
TgGate servo time constant (s)
VELMGate velocity limit (pu/s)
GmaxMaximum gate opening (pu)
GminMinimum gate opening (pu)
TwWater starting time (s)
AtTurbine gain (pu power / pu flow)
DturbTurbine damping factor
QnlNo-load water flow (pu)

Internal parameter: Po = (P / At) + Qnl (initial gate setpoint from load-flow)

SYNC_MACH g1 g1 1. 1. 0. 0. 800. 760. 3. 0. 0.95
XT 0.15 1.1 0.25 0.2 0.7 * 0.2 0.1 6.0257 0. 5.00 0.05 * 0.1
EXC GENERIC1 1.8991 -0.1 0. 1. 100. -1. -11 10. 70. 10. 20. 0.1 0. 4.
1 75. 15. 0.2 0.01 0.2 0.01 -0.1 0.1
TOR hygov 0.05 ! R: permanent droop (pu/pu)
0.30 ! r: transient droop (pu/pu)
5.00 ! Tr: transient droop reset time (s)
0.05 ! Tf: filter time constant for speed error (s)
0.20 ! Tg: gate servo time constant (s)
0.20 ! VELM: gate velocity limit (pu/s)
1.00 ! Gmax: maximum gate opening (pu)
0.0 ! Gmin: minimum gate opening (pu)
1.50 ! Tw: water starting time (s)
1.10 ! At: turbine gain (pu power / pu flow)
0.50 ! Dturb: turbine damping factor
0.05 ; ! Qnl: no-load water flow (pu)