modelPEN
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| Boolean | explicit | true | |
| Boolean | constantDiff | false | |
| Boolean | nonLinConc | true | |
| Boolean | leak | false | |
| Integer | nX | 10 | Number of species in the mixture 1: H2O, 2: H2, 3: CO2, 4: CO, 5: CH4, 6: O2 |
| Modelica.Units.SI.FaradayConstant | F | Modelica.Constants.F | Faraday constant C/mol |
| Types.SpecificHeatCapacityMol | R | Modelica.Constants.R | Universal gas constant per unit mol |
| Types.Area | S | W*L | Channel contact surface with electrodes |
| Types.SpecificHeatCapacity | cm | 400 | Specific heat capacity of the PEN |
| Types.HeatCapacity | Cm | M*cm | Heat capacity of the PEN |
| Types.Length | L | 5.1e-4 | Channel length |
| Types.Mass | M | rho*V | PEN mass |
| Types.Density | rho | 6600 | PEN density |
| Types.Area | A | tauAE*W | Electrode cross section |
| Types.Length | tauAE | 50e-6 | Anode electrode thickness |
| Types.Length | tauCE | 50e-6 | Cathode electrode thickness |
| Types.Length | tauSE | 150e-6 | Solid electrolyte thickness |
| Types.Volume | V | S*(tauAE + tauCE + tauSE) | PEN volume |
| Types.Length | W | 5.56e-3 | Channel width |
| Real | beta | 0.5 | Transfer coefficient |
| Types.DiffusionCoefficient | dEffAE | 3.66e-5 | Diffusivity coefficient for the anode electrode |
| Types.DiffusionCoefficient | dEffCE | 1.37e-5 | Diffusivity coefficient for the cathode electrode |
| Types.MolarEnergy | Eanode | 140e3 | |
| Types.MolarEnergy | Ecathode | 137e3 | |
| Types.SpecificConductivityArea | kAnode | 6.54e11 | |
| Types.SpecificConductivityArea | kCathode | 2.35e11 | |
| Types.Temperature | T_start | 920 | |
| Types.MassFraction[nX] | X_start_anode | ones(nX)/nX | |
| Integer | dominantSpecAnode | integer(Modelica.Math.Vectors.find(max(X_start_anode), X_start_anode)) | |
| Types.MassFraction[nX] | X_start_cathode | ones(nX)/nX | |
| Integer | dominantSpecCathode | integer(Modelica.Math.Vectors.find(max(X_start_cathode), X_start_cathode)) | |
| Types.Pressure | p_start | 101325 | |
| Types.MolarMass[nX] | MM | DiffusionLosses.Utilities.ConstantsFluidAndSolid.MM | |
| Types.Pressure[nX] | pX_anode_start | p_start/MM[dominantSpecAnode]*X_start_anode.*MM | |
| Types.Pressure[nX] | pX_cathode_start | p_start/MM[dominantSpecCathode]*X_start_cathode.*MM | |
| Types.CurrentDensity | j_start | 1.5 | |
| Types.CurrentDensity | j0Anode_start | 15000 | |
| Types.CurrentDensity | j0Cathode_start | 6000 | |
| Types.Temperature | Tnom | 950 | Nominal value for PEN temperature |
| Modelica.Units.SI.DimensionlessRatio | logVal | 1 | Homotopy parameter to simplify vOCP |
| Types.PerUnit | porosity | 0.35 | Porosity coefficient for anode and cathode electrodes |
| Types.ElPotential | vRef1 | 1.15 | |
| Types.ElPotential | vRef2 | 0.8 | |
| Types.ElCurrent | iRef1 | 0 | |
| Types.ElCurrent | iRef2 | 3 | |
| Types.ElCurrent | iPENmax | ||
| Types.ElPotential | vLeak0 | ||
| Types.Pressure | pH2Ofixed | 14e5 | |
| Types.Pressure | pH2fixed | 2.5e5 | |
| Types.Pressure | pO2fixed | 3.5e5 |
Connectors
| Type | Name | Default | Description |
|---|---|---|---|
| Interfaces.HeatPort | cathodeHeatPort | ||
| Interfaces.HeatPort | anodeHeatPort | ||
| Interfaces.FuelCellInterfaces.ElPin | anodePin | ||
| Interfaces.FuelCellInterfaces.ElPin | cathodePin | ||
| Interfaces.FuelCellInterfaces.PressureInputSignal[nX] | partialPressuresAnode | ||
| Interfaces.FuelCellInterfaces.PressureInputSignal | absolutePressureAnode | ||
| Interfaces.FuelCellInterfaces.PressureInputSignal | absolutePressureCathode | ||
| Interfaces.FuelCellInterfaces.PressureInputSignal[nX] | partialPressuresCathode | ||
| Interfaces.FuelCellInterfaces.RHOROutputSignal | rHOR | ||
| Interfaces.HeatPort | radCathode | ||
| Interfaces.HeatPort | radAnode |
Components
| Type | Name | Default | Description |
|---|---|---|---|
| Types.SpecificMolarGibbsFreeEnergy | dg0HOR | H2 oxidation reaction rate | |
| Types.ElCurrent | i | PEN current | |
| Types.CurrentDensity | j | PEN current density relative to channel-PEN contact area | |
| Types.CurrentDensity | j0Anode | ||
| Types.CurrentDensity | j0Cathode | ||
| Types.Power | elP | Power output | |
| Types.ElPotential | vOCP | PEN Open Circuit voltage (Nernst law) | |
| Types.ElPotential | vPEN | PEN voltage | |
| Types.ElPotential | vActAnode | Voltage loss due to activation phenomena in the anode | |
| Types.ElPotential | vActCathode | Voltage loss due to activation phenomena in the cathode | |
| Types.ElPotential | vConc | Voltage loss due to concentration effects | |
| Types.ElPotential | vOhm | Voltage loss due to ohmic effects | |
| Types.HeatFlowRate | Qanode | ||
| Types.HeatFlowRate | Qcathode | ||
| Types.Temperature | T | ||
| Types.Pressure | pAnode | ||
| Types.Pressure | pCathode | ||
| Types.Pressure | pH2 | Partial pressure H2 in the channel | |
| Types.Pressure | pH2O | Partial pressure H2O in the channel | |
| Types.Pressure | pO2 | Partial pressure O2 in the channel | |
| Types.Pressure | pH2TPB | Triple phase boundary H2 partial pressure | |
| Types.Pressure | pH2OTPB | Triple phase boundary H2O partial pressure | |
| Types.Pressure | pO2TPB | Triple phase boundary O2 partial pressure | |
| Types.MolarEnergy | RT | Variable expressing the product R*T: in homotopy lambda0 step it is R*Tnom | |
| Resistance | rOhm | ||
| HOR | hor | ||
| DiffusionLosses.ConcResitanceDGM | concLosses | ||
| Components.FuelCell.DiffusionLosses.DiffusivityCoefficients | D | ||
| Types.DiffusionCoefficient | dEffH2O | if constantDiff then dEffAE else D.D_H2O_p | |
| Types.DiffusionCoefficient | dEffH2 | if constantDiff then dEffAE else D.D_H2_p | |
| Types.DiffusionCoefficient | dEffO2 | if constantDiff then dEffCE else D.D_O2_p | |
| Components.HeatTransferModels.RadiativeHeatTransfer | heatRadAnode | ||
| Components.HeatTransferModels.RadiativeHeatTransfer | heatRadCathode | ||
| Types.ElPotential | vLeak | if leak then vLeak0*(1 - i/iPENmax) else 0 | |
| Real | logValCheck | logBound(pH2O/pH2/sqrtBound(pO2/101325)) | |
| Real | logValCheck2 | -RT/2/F*logBound(pH2O/pH2/sqrtBound(pO2/101325)) | |
| Real | vHomotopyLinCheck | (i - iRef1)/(iRef2 - iRef1)*(vRef2 - vRef1) + vRef1 | |
| Real | logBoundCHeck | logBound(pH2OTPB/pH2O*pH2/pH2TPB) + 0.5*logBound(pO2/pO2TPB) | |
| Real | logBoundWorse | log(1e-9) + 0.5*log(1e-9) | |
| Types.Pressure | pH2Oh | ||
| Types.Pressure | pH2h | ||
| Types.Pressure | pO2h |
Contents
| Name | Description |
|---|---|
| Resistance | |
| HOR |