modelPEMElectrolyzer_L1

PEMElectrolyzer_L1 Proton exchange membrane electrolyzer
Diagram of PEMElectrolyzer_L1

Extends from TransiEnt.Producer.Gas.Electrolyzer.Base.PartialElectrolyzer (partial class for electrolyzer).

Information

1. Purpose of model

This is a model for an electrolyzer with a replaceable efficiency curve and replaceable dynamic behavior.

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

The efficiency curve, the dynamic behaviour and the wanted input (electric power or hydrogen mass flow) can be chosen. The water consumption is calculated using a constant factor.

3. Limits of validity

(no remarks)

4. Interfaces

epp: electric power port, type can be chosen

gasPortOut: hydrogen outlet

P_el_set: input for electric power

m_flow_H2_set: input for hydrogen mass flow

5. Nomenclature

(no elements)

6. Governing Equations

The hydrogen mass flow or the electric power is calculated depending on the given input and chosen efficiency curve and dynamic behavior.

Heat outcoupling is calculated via the efficiency curve and a simplified correlation for heat losses depending on the electrical power, the operating temperature and the ambient temperature. The default value for the temperature coefficient is based on [1].

7. Remarks for Usage

Via parameter 'useLeakageMassFlow' a small mass flow of 'm_flow_small' is always flowing out of gas port (to avoid Zero-Mass-Flow problems).

The model works with any fluids which contain hydrogen as the last component.

8. Validation

Tested in the check models "TransiEnt.Producer.Gas.Electrolyzer.Check.TestPEMElectrolyzer_L1_Charline" and "TransiEnt.Producer.Gas.Electrolyzer.Check.TestPEMElectrolyzer_L1_Dynamics"

9. References

[1] Espinosa-López, M. et al., Modelling and experimental validation of a 46 kW PEM high pressure water electrolyzer, 2018

10. Version History

Model created by Carsten Bode (c.bode@tuhh.de) in March 2017

Model generalized for different electrical power ports by Jan-Peter Heckel (jan.heckel@tuhh.de) in July 2018

Model modified by Oliver Schülting (oliver.schuelting@tuhh.de) on Nov 2018: added useLeakageMassFlow

Model modified by Oliver Schülting (oliver.schuelting@tuhh.de) on Nov 2019: added simplified heat model to simulate thermal losses heat outcoupling

Model adjusted for base class by Jan Westphal (j.westphal@tuhh.de) in dec 2019

Parameters

TypeNameDefaultDescription
SI.MassFraction[medium.nc - 1]xi_out (from PartialElectrolyzer)zeros(medium.nc - 1)
SI.SpecificEnergy[:]NCV_H2 (from PartialElectrolyzer)TransiEnt.Basics.Functions.GasProperties.getRealGasNCVVector(medium, medium.nc)Net calorific value of hydrogen at 25 C and 1 bar
SI.SpecificEnergy[:]GCV_H2 (from PartialElectrolyzer)TransiEnt.Basics.Functions.GasProperties.getRealGasGCVVector(medium, medium.nc)Gross calorific value of hydrogen at 25 C and 1 bar
SI.SpecificEnthalpyh0 (from PartialElectrolyzer)TILMedia.Internals.VLEFluidConfigurations.FullyMixtureCompatible.VLEFluidFunctions.specificEnthalpy_pTxi(medium, 1e5, 298.15)Specific enthalpy at 25 C and 1 bar
EnergyResourcetypeOfResource (from PartialElectrolyzer)EnergyResource.ConsumerType of energy resource for global model statistics
Coolant
BooleanuseFluidCoolantPort (from PartialHeatProvision)falsechoose if fluid port for coolant shall be used
BooleanuseHeatPort (from PartialHeatProvision)falsechoose if heat port for coolant shall be used
BooleanexternalMassFlowControl (from PartialHeatProvision)falsechoose if coolant mass flow is defined by input
BooleanuseVariableCoolantOutputTemperature (from PartialHeatProvision)falsechoose if temperature of cooland output shall be defined by input
SI.TemperatureT_out_coolant_target (from PartialHeatProvision)500 + 273.15output temperature of coolant - will be limited by temperature which is technically feasible
SI.LinearTemperatureCoefficientk_therm_relative0.325436e-3temperature conductivity for heat losses relative to electrical power
SI.TemperatureT_amb273.15 + 15ambient temperature for heat loss calculation
Fundamental Definitions
TILMedia.VLEFluidTypes.BaseVLEFluidmedium (from PartialElectrolyzer)simCenter.gasModel3Medium model
SI.ActivePowerP_el_n (from PartialElectrolyzer)Nominal power of the electrolyzer
SI.ActivePowerP_el_max (from PartialElectrolyzer)Maximum power of the electrolyzer
SI.TemperatureT_out (from PartialElectrolyzer)283.15Hydrogen output temperature
BooleanuseHomotopysimCenter.useHomotopytrue if homotopy method is used during initialization
RealspecificWaterConsumption10Mass of water per mass of hydrogen
SI.Efficiencyeta_n0.75Nominal efficiency refering to the GCV (min = 0, max = 1)
SI.Efficiencyeta_scale0Sets a with increasing input power linear degrading efficiency coefficient (min=0,max=1)
IntegerwhichInput1use P_el_set or m_flow_H2_set as input
Replaceable Components
BooleanusePowerPort (from PartialElectrolyzer)trueTrue if power port shall be used
Statistics
BooleanintegrateH2Flow (from PartialElectrolyzer)falsetrue if hydrogen mass flow shall be integrated
BooleanintegrateElPower (from PartialElectrolyzer)simCenter.integrateElPowertrue if electric powers shall be integrated
BooleancalculateCost (from PartialElectrolyzer)simCenter.calculateCosttrue if cost shall be calculated
TransiEnt.Basics.Units.MonetaryUnitPerEnergyCspec_demAndRev_el (from PartialElectrolyzer)simCenter.Cspec_demAndRev_freeSpecific demand-related cost per electric energy
RealCspec_demAndRev_other (from PartialElectrolyzer)simCenter.Cspec_demAndRev_other_freeSpecific demand-related cost per cubic meter water
Numerical Stability
BooleanuseLeakageMassFlowfalseConstant leakage gas mass flow of 'm_flow_small' to avoid zero mass flow
SI.MassFlowRatem_flow_smallsimCenter.m_flow_smallleakage mass flow if useLeakageMassFlow=true

Connectors

TypeNameDefaultDescription
Basics.Interfaces.Thermal.FluidPortInfluidPortIn (from PartialHeatProvision)
Basics.Interfaces.Thermal.FluidPortOutfluidPortOut (from PartialHeatProvision)
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_aheat (from PartialHeatProvision)
Basics.Interfaces.General.TemperatureInT_set_coolant_out (from PartialHeatProvision)
TransiEnt.Basics.Interfaces.Electrical.ActivePowerPortepp (from PartialElectrolyzer)
TransiEnt.Basics.Interfaces.Gas.RealGasPortOutgasPortOut (from PartialElectrolyzer)
TransiEnt.Basics.Interfaces.Electrical.ElectricPowerInP_el_setElectric power input (set value)
TransiEnt.Basics.Interfaces.General.MassFlowRateInm_flow_H2_setHydrogen mass flow rate input (set value)

Components

TypeNameDefaultDescription
TransiEnt.SimCentersimCenter (from PartialHeatProvision)
Components.Boundaries.Heat.Heatflow_L1heatFlow_externalMassFlowControl (from PartialHeatProvision)
Modelica.Blocks.Sources.RealExpressionQ_flow_positive (from PartialHeatProvision)
Modelica.Thermal.HeatTransfer.Sources.PrescribedHeatFlowprescribedHeatFlow (from PartialHeatProvision)
Components.Boundaries.Heat.Heatflow_L1_idContrMFlow_tempheatflow_L1_idContrMFlow_temp (from PartialHeatProvision)
Modelica.Blocks.Sources.RealExpressionrealExpression8 (from PartialHeatProvision)
SI.HeatFlowRateQ_flow_heatprovision (from PartialHeatProvision)
SI.TemperatureT_out_coolant (from PartialHeatProvision)min(T_out_coolant_target, T_out_coolant_max)
SI.TemperatureT_out_coolant_max (from PartialHeatProvision)
TransiEnt.ModelStatisticsmodelStatistics (from PartialElectrolyzer)
SI.HeatFlowRateQ_flow (from PartialElectrolyzer)waste heat
SI.PowerP_el (from PartialElectrolyzer)Electric power consumed by the electrolyzer
SI.MassFlowRatem_flow_H2O (from PartialElectrolyzer)water mass flow rate into the electrolyzer
Modelica.Blocks.Sources.RealExpressionrealExpression (from PartialElectrolyzer)
TransiEnt.Components.Boundaries.Electrical.ActivePower.PowerpowerBoundary (from PartialElectrolyzer)
TILMedia.Internals.VLEFluidConfigurations.FullyMixtureCompatible.VLEFluid_pTvleFluidH2 (from PartialElectrolyzer)
SI.Massmass_H2produced H2 mass
SI.HeatFlowRateQ_lossheat losses to environment
SI.Efficiencyeta_NCVEfficiency of the electrolyzer based on NCV
SI.Efficiencyeta_GCVcharline.etaEfficiency of the electrolyzer based on GCV
Summarysummary

Contents

NameDescription
Outlineprotected
Summaryprotected
Dynamics
Charline