modelGasturbine_with_Gasport

Extends from TransiEnt.Producer.Electrical.Conventional.Components.NonlinearThreeStatePlant (Slew Rate limited (=nonlinear), Minimum power limited (shuts down below minimum power), with primary and secondary balancing controller where secondary balancing power is lumped inside), TransiEnt.Producer.Electrical.Base.PartialNaturalGasUnit (Adds a gas interface with a mass flow defining boundary to child components), TransiEnt.Producer.Electrical.Base.PartialCDEUnit (Adds a gas interface with a mass flow defining boundary to child components), TransiEnt.Basics.Icons.GasTurbine.

Information

1. Purpose of model

Combined Cycle Gas turbine model with startup time and setpoint limits.

Preconfigured combined cycle gas plant model with gas port.

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

(no remarks)

3. Limits of validity

4. Interfaces

P_el_set: input for electric power in [W] (setpoint for electric power)

P_SB_set: input for electric power in [W] (secondary balancing setpoint)

epp: active power port (choice of power port)

gasPortIn: inlet for real gas

gasPortIn_2: second gas port - only active if useSecondGasPort==true

gasPortOut_CDE: outlet port for CDE. only active if 'useCDEPort=true'

5. Nomenclature

(no remarks)

6. Governing Equations

The mass flow of gasPortIn is claculated via needed electrical power output and the efficiency of the power plant. If the second gasPort 'gasPortIn_2' is used, the mass flow of the first gas port is reduced depending on the enthalpy flow of the second gas port. The mass flow of the second gas port needs to be defined from outside.

7. Remarks for Usage

(no remarks)

8. Validation

(no remarks)

9. References

(no remarks)

10. Version History

Model created by Oliver Schülting (oliver.schuelting@tuhh.de) on Nov 2018

Parameters

TypeNameDefaultDescription
IntegernSubgrids (from PartialElectricPowerPlant)simCenter.iDetailedGridFor calculation of statistics in subgrids (=1 local grid, e.g. hamburg, =2 surrounding grid, e.g. UCTE grid
Modelica.Units.SI.PowerP_min (from NonlinearThreeStatePlant)P_min_star*P_el_n
SI.PowerP_init (from NonlinearThreeStatePlant)Start value for P
Physical Constraints
Modelica.Units.SI.PowerP_el_n (from PartialElectricPowerPlant)300e6Nominal power of plant
SI.PowerP_n (from PartialBalancingPowerProvider)300e6Nominal power of plant
RealP_min_star (from NonlinearThreeStatePlant)0.2Fraction of nominal power (=20% of nominal power)
RealP_max_star (from NonlinearThreeStatePlant)1Fraction of nominal power (=100% of nominal power)
RealP_grad_max_star (from NonlinearThreeStatePlant)0.12/60Fraction of nominal power per second (12% per minute)
Base.PartloadEfficiency.PartloadEfficiencyCharacteristicEfficiencyCharLine (from NonlinearThreeStatePlant)TransiEnt.Producer.Electrical.Base.PartloadEfficiency.ConstantEfficiency()choose characteristic efficiency line
SI.TimeH (from NonlinearThreeStatePlant)6Time constant of plant (For easy definition of J)
SI.InertiaJ (from NonlinearThreeStatePlant)P_el_n*2*H/(2*simCenter.f_n*Modelica.Constants.pi)^2Lumped moment of inertia of whole power plant
SI.Timet_startup (from NonlinearThreeStatePlant)0Startup time (P=0 during startup)
SI.Timet_min_operating (from NonlinearThreeStatePlant)0Minimum operation time
BooleansmoothShutDown (from NonlinearThreeStatePlant)trueshut down process will be smoothed - power gradient will be limit by 'P_grad_operating'
SI.TimeMinimumDownTime (from NonlinearThreeStatePlant)0Minimum time the plant needs to be shut down before starting again
SI.TimeT_plant (from NonlinearThreeStatePlant)0.63/P_grad_max_starTime constant of first order model representing plant dynamics
SI.MassFractionCO2_Deposition_Rate (from NonlinearThreeStatePlant)0Fraction of CO2 that is deposited via CCS
Base.CCS.NoCCSCCS_Characteristics (from NonlinearThreeStatePlant)Base.CCS.NoCCS()Choose characteristic efficiency losses due to CCS
Integerquantity (from NonlinearThreeStatePlant)1amount of power plant blocks into which nominal power is split
BooleanuseSecondGasPort (from PartialNaturalGasUnit)false
Statistics
Modelica.Units.SI.Efficiencyeta_total (from PartialElectricPowerPlant)0.6Total efficiency of plant for emission calculation
EnergyResourcetypeOfResource (from PartialElectricPowerPlant)EnergyResource.ConventionalType of energy resource for global model statistics
PrimaryEnergyCarriertypeOfPrimaryEnergyCarrier (from PartialElectricPowerPlant)PrimaryEnergyCarrier.BlackCoalType of primary energy carrier for co2 emissions global statistics
BooleanintegrateElPower (from PartialElectricPowerPlant)simCenter.integrateElPowertrue if electric powers shall be integrated
BooleanintegrateCDE (from PartialElectricPowerPlant)simCenter.integrateCDEtrue if CDE shall be integrated
BooleancalculateCost (from PartialElectricPowerPlant)simCenter.calculateCosttrue if costs shall be calculated
EnergyResourcetypeOfBalancingPowerResource (from PartialBalancingPowerProvider)EnergyResource.ConventionalType of energy resource for global model statistics
IntegernSubgrid (from NonlinearThreeStatePlant)1Index of subgrid for moment of inertia statistics
Block control › Frequency Control
BooleanisPrimaryControlActive (from PartialBalancingPowerProvider)true
BooleanisSecondaryControlActive (from PartialBalancingPowerProvider)false
BooleanisExternalSecondaryController (from PartialBalancingPowerProvider)trueFalse, provider generates its own control setpoint (only for lumped plants)
SI.Timet_SB_act (from PartialBalancingPowerProvider)simCenter.t_SB_actMaximum reaction time for SB
Modelica.Units.SI.PowerP_el_grad_max_SB (from PartialBalancingPowerProvider)0.02*P_n/60Maximum power gradient for secondary balancing (default: 2%/min)
Initialization
Booleanset_P_init (from NonlinearThreeStatePlant)trueUse of P_init or calculation of P_init from grid (CPP)
RealP_init_set (from NonlinearThreeStatePlant)P_el_nInitial or guess value of turbine power
Expert Settings
TransiEnt.Components.Turbogroups.Base.GradientLimitingChoicespowerGradLimChoice (from NonlinearThreeStatePlant)TransiEnt.Components.Turbogroups.Base.GradientLimitingChoices.GradLimInFirstOrderoptions of gradient limitation
BooleanfixedStartValue_w (from NonlinearThreeStatePlant)falseWhether or not the start value of the angular velocity of the plants mechanical components is fixed
SI.TimeTd_GradientLimiter (from NonlinearThreeStatePlant)simCenter.TdTime constant of integrator
BooleanuseThresh (from NonlinearThreeStatePlant)simCenter.useThreshUse threshould for suppression of numerical noise
Realthres (from NonlinearThreeStatePlant)simCenter.thresIf 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. For too small values, the power plant is never turned off.
Realthres_hyst (from NonlinearThreeStatePlant)1e-10Threshold for hysteresis for switch from halt to startup (chattering might occur, hysteresis might help avoiding this)
BooleanuseHomotopyVarSlewRateLim (from NonlinearThreeStatePlant)simCenter.useHomotopytrue if homotopy shall be used in variableSlewRateLimiter in turbine (true might avoid circular equalities, false might help if 0/0 occurs in power plant)
Fundamental Definitions
TILMedia.VLEFluidTypes.BaseVLEFluidmediumGas (from PartialNaturalGasUnit)simCenter.gasModel1
BooleanuseCDEPort (from PartialCDEUnit)false
TILMedia.VLEFluidTypes.BaseVLEFluidmedium_CDE_deposited (from PartialCDEUnit)simCenter.gasModel1
Configuration
BooleanuseGasPort (from PartialNaturalGasUnit)trueTrue if gas port shall be used
Fuel properties
BooleanuseConstantHoCfalseif ticked, a constant heat of combustion for the fuel can be defined
Modelica.Units.SI.SpecificEnthalpyHoC_gas40e6heat of combustion of natural gas
Numerical Stability
BooleanuseLeakageMassFlowfalseConstant leakage gas mass flow of 'm_flow_small' to avoid zero mass flow
Modelica.Units.SI.MassFlowRatem_flow_smallsimCenter.m_flow_smallleakage mass flow if useLeakageMassFlow=true

Connectors

TypeNameDefaultDescription
TransiEnt.Basics.Interfaces.Electrical.ActivePowerPortepp (from PartialElectricPowerPlant)
TransiEnt.Basics.Interfaces.Electrical.ElectricPowerInP_el_set (from PartialDispatchablePowerPlant)Electric power setpoint
TransiEnt.Basics.Interfaces.Electrical.ElectricPowerInP_SB_set (from PartialBalancingPowerProvider)Secondary balancing setpoint
Modelica.Blocks.Interfaces.BooleanInputuseCCS (from NonlinearThreeStatePlant)true, if CCS is used
TransiEnt.Basics.Interfaces.Gas.RealGasPortIngasPortIn (from PartialNaturalGasUnit)
TransiEnt.Basics.Interfaces.General.MassFlowRateOutm_flow_gas (from PartialNaturalGasUnit)
TransiEnt.Basics.Interfaces.Gas.RealGasPortIngasPortIn_2 (from PartialNaturalGasUnit)
TransiEnt.Basics.Interfaces.Gas.RealGasPortOutgasPortOut_CDE (from PartialCDEUnit)
TransiEnt.Basics.Interfaces.Gas.EnthalpyFlowRateOutH_flow_gasEnthalpy flow rate of fuel

Components

TypeNameDefaultDescription
TransiEnt.SimCentersimCenter (from PartialElectricPowerPlant)
TransiEnt.ModelStatisticsmodelStatistics (from PartialElectricPowerPlant)
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectElectricPowercollectElectricPower (from PartialElectricPowerPlant)
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectGwpEmissionsElectriccollectGwpEmissions (from PartialElectricPowerPlant)
TransiEnt.Components.Statistics.Collectors.LocalCollectors.PowerPlantCostcollectCosts (from PartialElectricPowerPlant)
SI.ActivePowerP_el_is (from PartialElectricPowerPlant)-epp.P
RealP_star (from PartialElectricPowerPlant)P_el_is/P_el_n
SI.Efficiencyeta (from PartialElectricPowerPlant)
SI.Timet_fullload (from PartialElectricPowerPlant)E_total_generation/P_el_nfull load time in hours
SI.EnergyE_total_generation (from PartialElectricPowerPlant)Start value for generated electricity statistics
Booleanis_running (from PartialElectricPowerPlant)For continuous plants always true, for discontinous depending on state
SI.EnthalpyFlowRateQ_flow_fuel_is (from PartialElectricPowerPlant)P_el_is/eta
SI.MassFlowRatem_flow_CDE (from PartialElectricPowerPlant)fuelSpecificCO2Emissions.m_flow_CDE_per_Energy*abs(P_el_is)/eta - m_flow_gas_CDE_deposited
SI.MassFlowRatem_flow_gas_CDE_deposited (from PartialElectricPowerPlant)0*time
Realdelta_f_star (from PartialElectricPowerPlant)(epp.f - simCenter.f_n)/simCenter.f_n
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectElectricPowergeneratedSBP (from PartialBalancingPowerProvider)
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectElectricPower[2]collectSBPOffer (from PartialBalancingPowerProvider)
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectElectricPowergeneratedPBP (from PartialBalancingPowerProvider)
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectElectricPower[2]collectPBPOffer (from PartialBalancingPowerProvider)
TransiEnt.Producer.Electrical.Controllers.PrimaryBalancingControllerprimaryBalancingController (from PartialBalancingPowerProvider)
TransiEnt.Components.Sensors.MechanicalFrequencygridFrequencySensor (from PartialBalancingPowerProvider)
TransiEnt.Producer.Electrical.Base.ControlPower.IdealControlOffercontrolPowerModel (from PartialBalancingPowerProvider)
RealP_SB_set_star (from PartialBalancingPowerProvider)P_SB_set_internal/P_n
TransiEnt.Components.Electrical.Machines.ActivePowerGeneratorGenerator (from NonlinearThreeStatePlant)
TransiEnt.Components.Mechanical.TwoStateInertiaWithIdealClutchMechanicalConnection (from NonlinearThreeStatePlant)
TransiEnt.Components.Electrical.Machines.ExcitationSystemsVoltageController.DummyExcitationSystemExciter (from NonlinearThreeStatePlant)
TransiEnt.Components.Turbogroups.FirstOrderThreeStateTurbineTurbine (from NonlinearThreeStatePlant)
Modelica.Blocks.Logical.Switchswitch1 (from NonlinearThreeStatePlant)
Modelica.Blocks.Sources.RealExpressionCO2_Deposited_set (from NonlinearThreeStatePlant)just for visualisation on diagram layer
Modelica.Blocks.Sources.RealExpressionZero (from NonlinearThreeStatePlant)just for visualisation on diagram layer
Modelica.Blocks.Math.Gaingain (from NonlinearThreeStatePlant)
Modelica.Blocks.Math.SumP_set_total (from NonlinearThreeStatePlant)Schedule plus agc setpoint
Base.PartloadEfficiency.PartloadEfficiencyPartloadEfficiency (from NonlinearThreeStatePlant)
Modelica.Blocks.Sources.RealExpressionp_is (from NonlinearThreeStatePlant)
Modelica.Blocks.Nonlinear.VariableLimitervariableLimiter (from NonlinearThreeStatePlant)
Modelica.Blocks.Sources.RealExpressionP_max_with_CCS (from NonlinearThreeStatePlant)
Modelica.Blocks.Sources.RealExpressionP_min_with_CCS (from NonlinearThreeStatePlant)
SI.MassFlowRatem_flow_gas_CDE_deposited_set (from NonlinearThreeStatePlant)m_flow_CDE/(1 - CO2_Deposition_Rate)*CO2_Deposition_Rate
RealP_set_star (from NonlinearThreeStatePlant)-Turbine.P_target/P_el_n
RealP_set_star_sched (from NonlinearThreeStatePlant)-P_el_set/P_el_n
Realdelta_P_star (from NonlinearThreeStatePlant)(P_el_set + P_el_is)/P_el_n
Modelica.Blocks.Math.Gaingain1 (from NonlinearThreeStatePlant)
TransiEnt.Components.Boundaries.Gas.BoundaryRealGas_hxim_flowmassFlowSink (from PartialNaturalGasUnit)
TransiEnt.Components.Sensors.RealGas.NCVSensorvleNCVSensor (from PartialNaturalGasUnit)
TransiEnt.Components.Boundaries.Gas.BoundaryRealGas_pTxmassFlowSink_2 (from PartialNaturalGasUnit)
TransiEnt.Components.Sensors.RealGas.NCVSensorvleNCVSensor_2 (from PartialNaturalGasUnit)
Modelica.Blocks.Sources.RealExpressionm_flow_gas_2_var (from PartialNaturalGasUnit)
Modelica.Blocks.Math.ProductH_flow_secondGasPort (from PartialNaturalGasUnit)
Modelica.Blocks.Sources.RealExpressionm_flow_gas_2_const (from PartialNaturalGasUnit)
Modelica.Blocks.Routing.RealPassThroughm_flow_gas_2 (from PartialNaturalGasUnit)
Modelica.Blocks.Sources.RealExpressionNCV_2_const (from PartialNaturalGasUnit)
Modelica.Blocks.Routing.RealPassThroughNCV_2 (from PartialNaturalGasUnit)
TransiEnt.Components.Boundaries.Gas.BoundaryRealGas_hxim_flowmassFlowSink_CDE_deposited (from PartialCDEUnit)
Modelica.Blocks.Sources.RealExpressionm_flow_set_CDE_deposited (from PartialCDEUnit)just for visualisation on diagram layer
SI.MassFlowRatem_flow_gas_CDE_deposited_gasPort (from PartialCDEUnit)0*time
Modelica.Units.SI.MassFlowRatem_flow_CDE_calculatedByInput
Modelica.Units.SI.MassFlowRatem_flow_small_intern
Modelica.Units.SI.SpecificEnthalpyHoC_gas_actual