modelPEMElectrolyzer_L2

PEMElectrolyzer_L2 Proton exchange membrane electrolyzer with selectable physics submodels

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

Information

1. Purpose of model

This model electrolyzer uses modular physics classes and a Specification record to describe a real-life electrolyzer system. The default model uses physical relationships taken from (Espinosa-López et al, 2018).

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

The user can select or create a model with system specific parameters. The desired input (electric power, current, current density or hydrogen mass flow) can be varied. Physics submodels can be replaced as desired, with essential definitions annotated. The water consumption, hydrogen output, and net- and gross-calorific energy conversion efficiency can be calculated. The user has the option of controlling pressure through gasPortOut and/or temperature through T_input as well.

3. Limits of validity

Original model developed and validated in the range of 20-60 °C with operating pressure of 15-35 bar.

4. Interfaces

epp: electric power port, type can be chosen

gasPortOut: hydrogen outlet

i_dens_set: input for electric current density

i_el_stack_set: input for electric current

P_el_set: input for electric power

m_flow_H2_set: input for hydrogen mass flow

T_input: input for temperature

excessHeatFlowOut: Heat flow rate out port, equal to the cooling power used to regulate temperature at max temp.

temperatureOut: temperatureOut interface equal to operating temperature

5. Nomenclature

(no elements)

6. Governing Equations

Selectable physics equations allow for different governing equations to be used, and consist of equations from (Espinosa-López et al, 2018) by default.

7. Remarks for Usage

(no remarks)

8. Validation

Results have been validated against (Espinosa-López et al, 2018) published figures.

9. References

[1] Manuel Espinosa-López, Philippe Baucour, Serge Besse, Christophe Darras, Raynal Glises, Philippe Poggi, André Rakotondrainibe, and Pierre Serre-Combe. Modelling and experimental validation of a 46 kW PEM high pressure water electrolyser. Renewable Energy, 119, pp. 160-173, 2018. doi: 10.1016/J.RENENE.2017.11.081.

[2] efficiency curve of the inverter taken from the data sheet of SMA "Technische Wirkungsgrade und Derating" URL: https://files.sma.de/dl/1348/WirkungDerat-TI-de-46.pdf page 71, 26.11.2019

[3] J. Webster and C. Bode, “Implementation of a Non-Discretized Multiphysics PEM Electrolyzer Model in Modelica,” in Proceedings of the 13th International Modelica Conference, Regensburg, Germany, March 4–6, 2019, no. 157, pp. 833–840, DOI: 10.3384/ecp19157833.

10. Version History

Model created by John Webster (jcwebste@edu.uwaterloo.ca) in October 2018

Model adjusted for TransiEnt 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
SI.TemperatureT_std298.15STD temperature
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
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
SI.TemperatureT_amb23 + 273.15K, ambient temperature
IntegerwhichInput1Use current, current density, electric power, or mass_flow_H2 as input
BooleanuserSetTempfalseUse T_input as input
Realeta_inv_n0.956nominal efficiency of the inverter
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
Initialization
SI.TemperatureT_op_startT_stdinitial operating temperature of PEM Electrolyzer

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)
Modelica.Blocks.Interfaces.RealInputi_el_stack_set
Modelica.Blocks.Interfaces.RealInputi_dens_set
Modelica.Blocks.Interfaces.RealInputP_el_set
Modelica.Blocks.Interfaces.RealInputm_flow_H2_set
Modelica.Blocks.Interfaces.RealInputT_input
TransiEnt.Basics.Interfaces.General.TemperatureOuttemperatureOut
TransiEnt.Basics.Interfaces.Thermal.HeatFlowRateOutexcessHeatFlowOut

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.TemperatureT_opOperating stack temperature
SI.Efficiencyeta_condVariable for modeling the effiiciency loss of the electrolyseur due to the needed power for the dryer and the water conditioning
SI.EnthalpyFlowRateH_flow_H2_GCVH2 enthalpy flow rate out of electrolyzer, gross calorific value
SI.EnthalpyFlowRateH_flow_H2_NCVH2 enthalpy flow rate out of electrolyzer, net calorific value
SI.Efficiencyeta_NCVEfficiency of the electrolyzer based on NCV
SI.Efficiencyeta_GCVEfficiency of the electrolyzer based on GCV
SI.Currenti_el_stackCurrent across the electrolyzer stack
SI.CurrentDensityi_dens_aOperating current density at anode
SI.VoltageV_el_stackPEM stack voltage
SI.VoltageV_cellPEM cell voltage considering all included physical phenomena
SI.VoltageV_tnor U_tn, Thermoneutral voltage (voltage at which reaction can occur without releasing any heat)
SI.PressuregasPortPressurepressure of hydrogen connected to gasPortOut
SI.Pressurepp_H2OPa, vapour pressure of water vapour, must always be converted to atm
SI.Pressurepp_H2Pa, partial pressure of H2, must always be converted to atm
SI.Pressurepp_O2Pa, partial pressure of O2, must always be converted to atm
SI.Pressurep_catPressure at cathode (H2)
SI.Pressurep_anPressure at anode (O2)
SI.MassFlowRatem_flow_H2H2 mass flow rate out of electrolyzer
SI.MolarFlowRaten_flow_H2OMolar consummation rate of water
SI.MolarFlowRaten_flow_H2Molar production rate of Hydrogen gas
SI.MolarFlowRaten_flow_O2Molar production rate of Oxygen gas
TransiEnt.Producer.Gas.Electrolyzer.Base.Specifications.AREVAGiner46kWSpecification
electrolyzerVoltagevoltage
electrolyzerTemperaturetemperature
electrolyzerPressurespressure
electrolyzerMassFlowmassFlow
Summarysummary
Modelica.Blocks.Sources.RealExpressionT_op_out
Modelica.Blocks.Sources.RealExpressionQ_Flow_Coolingexcess waste heat generated by electrolyzer system, actively cooled by default
Modelica.Blocks.Tables.CombiTable1DsEfficiencyCurve_Inverter

Contents

NameDescription
Outline
Summary
electrolyzerVoltage
electrolyzerTemperature
electrolyzerPressures
electrolyzerMassFlow