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 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} DEGOV1 block diagram. The speed deviation is summed with the reference REF, which is V60 times the droop R, and drives a governor lead-lag with time constants T1, T2 and T3, then an actuator lead-lag of gain K with T4 and T6, then a lag T5. The result passes through an integrator limited between TMIN and TMAX, a Pade approximation of the engine dead time TD, and a lag TE, giving the mechanical torque Tm, from which Pm is Tm times omega. DEGOV1 block diagram. The speed deviation is summed with the reference REF, which is V60 times the droop R, and drives a governor lead-lag with time constants T1, T2 and T3, then an actuator lead-lag of gain K with T4 and T6, then a lag T5. The result passes through an integrator limited between TMIN and TMAX, a Pade approximation of the engine dead time TD, and a lag TE, giving the mechanical torque Tm, from which Pm is Tm times omega.
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)

Internal parameters: V60 (selected reference input, mechanical torque or electrical power according to SWM) and REF (V60 · R), both set at initialisation

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 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 ENTSOE_simp block diagram. The speed deviation omega minus one is subtracted from the load reference C, which is Tm times R, divided by the droop R, and drives a first-order lag of time constant T1 limited between VMIN and VMAX, then a lead-lag with zero T2 and pole T3, giving the mechanical power Pm and hence the torque Pm divided by omega. ENTSOE_simp block diagram. The speed deviation omega minus one is subtracted from the load reference C, which is Tm times R, divided by the droop R, and drives a first-order lag of time constant T1 limited between VMIN and VMAX, then a lead-lag with zero T2 and pole T3, giving the mechanical power Pm and hence the torque Pm divided by omega.
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 GAST model represents a gas turbine unit, following the IEEE Committee Report (1973) and the convention common to industry tools. 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 GAST block diagram. The speed deviation is subtracted from the load reference LR, divided by the droop R, and meets a minimum gate that also takes the temperature limit. The gated signal passes a limiter between VMIN and VMAX, a lag T1 for the valve and a lag T2 for the compressor discharge to give p3. The damping term Dt times the speed deviation is added to give the mechanical power. A parallel path takes p3 through a lag T3 and the expression AT plus KT times AT minus p4 to form the temperature limit fed back to the minimum gate. GAST block diagram. The speed deviation is subtracted from the load reference LR, divided by the droop R, and meets a minimum gate that also takes the temperature limit. The gated signal passes a limiter between VMIN and VMAX, a lag T1 for the valve and a lag T2 for the compressor discharge to give p3. The damping term Dt times the speed deviation is added to give the mechanical power. A parallel path takes p3 through a lag T3 and the expression AT plus KT times AT minus p4 to form the temperature limit fed back to the minimum gate.
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 TGOV1 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} TGOV1 block diagram. The speed deviation is subtracted from the load reference C, which is Tm times R, divided by the droop R, and drives a first-order lag T1 limited between VMIN and VMAX, then a lead-lag with zero T2 and pole T3 representing the steam chest and reheater. The damping term Dt times the speed deviation is added to the turbine output to give the mechanical power Pm, and the torque is Pm divided by omega. TGOV1 block diagram. The speed deviation is subtracted from the load reference C, which is Tm times R, divided by the droop R, and drives a first-order lag T1 limited between VMIN and VMAX, then a lead-lag with zero T2 and pole T3 representing the steam chest and reheater. The damping term Dt times the speed deviation is added to the turbine output to give the mechanical power Pm, and the torque is Pm divided by 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 TGOV1 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 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 HYGOV block diagram. One minus omega, less the permanent droop R times the gate demand offset, passes a lag Tf to give the error e. That error is divided by r times Tr, rate-limited to plus or minus VELM and integrated between the gate limits Gmin and Gmax. Adding e over r and clipping gives the gate demand c, which drives a servo of time constant Tg to give the actual gate g. The water column relates the flow Q and head h through the water-starting time TW, and the turbine power is At times h times Q minus the no-load flow. HYGOV block diagram. One minus omega, less the permanent droop R times the gate demand offset, passes a lag Tf to give the error e. That error is divided by r times Tr, rate-limited to plus or minus VELM and integrated between the gate limits Gmin and Gmax. Adding e over r and clipping gives the gate demand c, which drives a servo of time constant Tg to give the actual gate g. The water column relates the flow Q and head h through the water-starting time TW, and the turbine power is At times h times Q minus the no-load flow.
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)