modelGasturbine

Model of a gasturbine with three states (halt / startup / running), pyhsical constraints (Pmin,Pmax,Pgradmax) and first order dynamics
Diagram of Gasturbine

Extends from Base.PartialTurbine (Abstract model of a turbine defining nothing but the interfaces), TransiEnt.Basics.Icons.GasTurbine.

Information

1. Purpose of model

Turbine model modeled by

- minimum / maximum power ouput

- maximum power gradient

- on / off status dependent on schedule value without time delay. Turbine shuts down if scheduled value is below minimum power

- sends on/off status to momentum of inertia statistics

2. Level of detail, physical effects considered, and physical insight

(no remarks)

3. Limits of validity

Note, that no statistics are involved!

4. Interfaces

mpp: mechanical power port

P_target: input for electric power in W

P_spinning_set: input for electric power in W (setpoint for spinning reserve power)

isGeneratorRunning: BooleanOutput

5. Nomenclature

(no elements)

6. Governing Equations

(no equations)

7. Remarks for Usage

The input P_spinning_set is supposed to be gradient limited by limtations of the primary balancing offer mechanism which has normally a higher gradient limit than the rest of the plant.

The input P_target is the sum of secondary balancing setpoint and scheduled set point

8. Validation

(no remarks)

9. References

(no remarks)

10. Version History

Model created by Pascal Dubucq (dubucq@tuhh.de) on 01.10.2014

Model modified by Robert Flesch (flesch@xrg-simulation.de) in Feb 2021: added feedback of actual power to control

Parameters

TypeNameDefaultDescription
SI.PowerP_n1e6Nominal power of plant
RealP_min_star0.3Dimensionless minimum power (=20% of nominal power)
RealP_max_star1Dimensionless maximum power (=Nominal power)
RealP_grad_max_star0.1/60Dimensionless maximum power gradient per second (12% of P_nom per minute)
RealP_turb_init0Initial or guess value of turbine power (in p.u.)
SI.TimeT_plant10Turbine first order dynamic
SI.Timet_startup360Startup time (no output during startup)
BooleanuseSlewRateLimitertruechoose if slewRateLimiter is activated
Numerical
RealTd1e-3The higher Nd, the closer y follows u
BooleanuseThreshfalseUse threshould for suppression of numerical noise
Realthres1e-7If abs(u-y)< thres, y becomes a simple pass through of u. Increasing thres can improve simulation speed. However to large values can make the simulation unstable. A good starting point is the choice thres = tolerance/1000.
Realthres_hyst1e-10Threshold for hysteresis for switch from halt to startup (chattering might occur, hysteresis might help avoiding this)
SI.Timet_eps10Threshold time for transitions

Connectors

TypeNameDefaultDescription
TransiEnt.Basics.Interfaces.General.MechanicalPowerPortmpp (from PartialTurbine)
TransiEnt.Basics.Interfaces.Electrical.ElectricPowerInP_target (from PartialTurbine)
Modelica.Blocks.Interfaces.BooleanOutputisGeneratorRunning
TransiEnt.Basics.Interfaces.Electrical.ElectricPowerInP_spinning_setSetpoint for spinning reserve power

Components

TypeNameDefaultDescription
Modelica.Units.SI.ActivePowerP_T (from PartialTurbine)mpp.tau*der(mpp.phi)
SimCentersimCenter
Modelica.Blocks.Math.Gainnormalize
Modelica.Blocks.Math.GaindeNormalize
Boundaries.Mechanical.PowerMechanicalBoundary
RealP_set_star-normalize.yNormalized, always positive
RealP_is_star-deNormalize.uNormalized, always positive
Modelica.Blocks.Continuous.FirstOrderplantDynamic
TransiEnt.Components.Turbogroups.OperatingStates.ThreeStateDynamicoperationStatus
Modelica.Blocks.Sources.BooleanExpressionbooleanExpression