RenGovernor#

Renewable turbine governor group.

Common Parameters: u, name, ree, w0, Sn, Pe0

Common Variables: Pm, wr0, wt, wg, s3_y

Available models: WTDTA1, WTDS

WTDTA1#

WTDTA wind turbine drive-train model.

One can set Htfrac to 0 to simulate a single-mass drive train. Htfrac has to be within [0, 1]

User-provided reference speed should be specified in parameter w0. Internally, w0 is set to the algebraic variable wr0.

Note for PSS/E dyr parser:

In PSS/E doc, Freq1 is said to be Hz, but exported data from PSS/E 34 uses per unit. ANDES requires Freq1 in per unit frequency.

Parameters#

Name

Symbol

Description

Default

Unit

Properties

idx

unique device idx

u

\(u\)

connection status

1

bool

name

device name

ree

Renewable exciter idx

mandatory

H

\(H_t\)

Total inertia constant

3

MWs/MVA

non_zero,non_negative,power

DAMP

\(Damp\)

Damp coefficient

0

p.u. (gen base)

power

Htfrac

\(D_{shaft}\)

Turbine inertia fraction (Hturb/H)

0.500

power

Freq1

\(Freq1\)

First shaft torsional resonant frequency, p.u. (Hz)

1

p.u. (Hz)

Dshaft

\(D_{shaft}\)

Shaft damping factor

1

p.u. (gen base)

power

w0

\(\omega_0\)

Default speed if not using a torque model

1

p.u.

reg

0

Sn

\(S_n\)

0

Variables#

Name

Symbol

Type

Description

Unit

Properties

s1_y

\(y_{s1}\)

State

Integrator output

v_str

s2_y

\(y_{s2}\)

State

Integrator output

v_str

s3_y

\(y_{s3}\)

State

Integrator output

v_str

wt

\(\omega_{t}\)

AliasState

Alias of s1_y

wg

\(\omega_{g}\)

AliasState

Alias of s2_y

wr0

\(\omega_{r0}\)

Algeb

speed set point

p.u.

v_str

Pm

\(P_{m}\)

Algeb

Mechanical power

v_str

pd

\(P_{d}\)

Algeb

Output after damping

v_str

wge

\(wge\)

ExtAlgeb

Pe

\(Pe\)

ExtAlgeb

Retrieved Pe of RenGen

Initialization Equations#

Name

Symbol

Type

Initial Value

s1_y

\(y_{s1}\)

State

\(\omega_{r0}\)

s2_y

\(y_{s2}\)

State

\(\omega_{r0}\)

s3_y

\(y_{s3}\)

State

\(\frac{P_{e0}}{\omega_{r0}}\)

wt

\(\omega_{t}\)

AliasState

wg

\(\omega_{g}\)

AliasState

wr0

\(\omega_{r0}\)

Algeb

\(\omega_{0}\)

Pm

\(P_{m}\)

Algeb

\(P_{e0}\)

pd

\(P_{d}\)

Algeb

\(0\)

wge

\(wge\)

ExtAlgeb

Pe

\(Pe\)

ExtAlgeb

Differential Equations#

Name

Symbol

Type

RHS of Equation "T x' = f(x, y)"

T (LHS)

s1_y

\(y_{s1}\)

State

\(- 1.0 P_{d} + \frac{1.0 P_{m}}{y_{s1}} - 1.0 y_{s3}\)

\(2H_t\)

s2_y

\(y_{s2}\)

State

\(- 1.0 Damp \left(- \omega_{0} + y_{s2}\right) + 1.0 P_{d} - \frac{1.0 Pe}{y_{s2}} + 1.0 y_{s3}\)

\(2H_g\)

s3_y

\(y_{s3}\)

State

\(K_{shaft} \left(y_{s1} - y_{s2}\right)\)

\(1.0\)

wt

\(\omega_{t}\)

AliasState

\(0\)

wg

\(\omega_{g}\)

AliasState

\(0\)

Algebraic Equations#

Name

Symbol

Type

RHS of Equation "0 = g(x, y)"

wr0

\(\omega_{r0}\)

Algeb

\(\omega_{0} - \omega_{r0}\)

Pm

\(P_{m}\)

Algeb

\(- P_{m} + P_{e0}\)

pd

\(P_{d}\)

Algeb

\(D_{shaft} \left(y_{s1} - y_{s2}\right) - P_{d}\)

wge

\(wge\)

ExtAlgeb

\(y_{s2} - 1.0\)

Pe

\(Pe\)

ExtAlgeb

\(0\)

Services#

Name

Symbol

Equation

Type

Ht2

\(2H_t\)

\(2 D_{shaft} H_{t}\)

ConstService

Hg2

\(2H_g\)

\(2 H_{t} \left(1 - D_{shaft}\right)\)

ConstService

Kshaft

\(K_{shaft}\)

\(\frac{0.5 \cdot 2H_{g} 2H_{t} Freq_{1}^{2}}{H_{t}}\)

ConstService

Blocks#

Name

Symbol

Type

Info

s1

\(s1\)

Integrator

s2

\(s2\)

Integrator

s3

\(s3\)

Integrator

Config Fields in [WTDTA1]

Option

Symbol

Value

Info

Accepted values

allow_adjust

1

allow adjusting upper or lower limits

(0, 1)

adjust_lower

0

adjust lower limit

(0, 1)

adjust_upper

1

adjust upper limit

(0, 1)

WTDS#

Custom wind turbine model with a single swing-equation.

This model is used to simulate the mechanical swing of the combined machine and turbine mass. The speed output is s1_y which will be fed to RenExciter.wg.

PFLAG needs to be set to 1 in exciter to consider speed for Pref.

Parameters#

Name

Symbol

Description

Default

Unit

Properties

idx

unique device idx

u

\(u\)

connection status

1

bool

name

device name

ree

Renewable exciter idx

mandatory

H

\(H_t\)

Total inertia

3

MWs/MVA

non_zero,non_negative,power

D

\(D_{shaft}\)

Damping coefficient

1

p.u.

power

w0

\(\omega_0\)

Default speed if not using a torque model

1

p.u.

reg

0

Sn

\(S_n\)

0

Variables#

Name

Symbol

Type

Description

Unit

Properties

s1_y

\(y_{s1}\)

State

Integrator output

v_str

s3_y

\(y_{s3}\)

State

Unused state variable

wt

\(\omega_{t}\)

AliasState

Alias of s1_y

wg

\(\omega_{g}\)

AliasState

Alias of s1_y

Pm

\(P_{m}\)

Algeb

Mechanical power

v_str

wr0

\(\omega_{r0}\)

Algeb

speed set point

p.u.

v_str

wge

\(wge\)

ExtAlgeb

Pe

\(Pe\)

ExtAlgeb

Retrieved Pe of RenGen

Initialization Equations#

Name

Symbol

Type

Initial Value

s1_y

\(y_{s1}\)

State

\(\omega_{r0}\)

s3_y

\(y_{s3}\)

State

wt

\(\omega_{t}\)

AliasState

wg

\(\omega_{g}\)

AliasState

Pm

\(P_{m}\)

Algeb

\(P_{e0}\)

wr0

\(\omega_{r0}\)

Algeb

\(\omega_{0}\)

wge

\(wge\)

ExtAlgeb

Pe

\(Pe\)

ExtAlgeb

Differential Equations#

Name

Symbol

Type

RHS of Equation "T x' = f(x, y)"

T (LHS)

s1_y

\(y_{s1}\)

State

\(- 1.0 D_{shaft} \left(- \omega_{r0} + y_{s1}\right) + \frac{1.0 \left(P_{m} - Pe\right)}{wge}\)

\(2H\)

s3_y

\(y_{s3}\)

State

\(0\)

wt

\(\omega_{t}\)

AliasState

\(0\)

wg

\(\omega_{g}\)

AliasState

\(0\)

Algebraic Equations#

Name

Symbol

Type

RHS of Equation "0 = g(x, y)"

Pm

\(P_{m}\)

Algeb

\(- P_{m} + P_{e0}\)

wr0

\(\omega_{r0}\)

Algeb

\(\omega_{0} - \omega_{r0}\)

wge

\(wge\)

ExtAlgeb

\(y_{s1} - 1.0\)

Pe

\(Pe\)

ExtAlgeb

\(0\)

Services#

Name

Symbol

Equation

Type

H2

\(2H\)

\(2 H_{t}\)

ConstService

Kshaft

\(K_{shaft}\)

\(1.0\)

ConstService

Blocks#

Name

Symbol

Type

Info

s1

\(s1\)

Integrator

Config Fields in [WTDS]

Option

Symbol

Value

Info

Accepted values

allow_adjust

1

allow adjusting upper or lower limits

(0, 1)

adjust_lower

0

adjust lower limit

(0, 1)

adjust_upper

1

adjust upper limit

(0, 1)