modelStack

Diagram of StackOpen the diagram
▶ Check in your browser

Parameters

TypeNameDefaultDescription
IntegerN10Number of volumes of each channel
IntegerNseries1Number of modules connected in series to build a stack of SOFCs
IntegerNparallelFluid1Number of modules in fluid dynamic parallel
IntegernX10Number of species in the fluid model
BooleanconstantDiffusionfalseTrue if diffusion coefficient for concnetration losses should be considered constant
BooleanleakfalseTrue if a linear loss proportional to the current should be included to introduce a fuel leakeage effect
BooleannonLinConctrueTrue if concentration losses should be considered as non linear
Boolean[N]initialEquationfill(true, N)If True
Types.SpecificHeatCapacitycmPEN400Specific heat capacity of the PEN
Types.SpecificHeatCapacitycmPlate500Specific heat capacity of the interconnecting plate
Types.MolarEnergyEanodeAnodic energy of activation for activation loss
Types.MolarEnergyEcathodeCathodic energy of activation for activation loss
Types.SpecificConductivityAreakAnodePre-exponential factor for anodic activation loss
Types.SpecificConductivityAreakCathodePre-exponential factor for cathodic activation loss
Types.PerUnitporosityPorosity coefficient for anode and cathode electrodes
Types.DensityrhoPEN6600PEN density
Types.DensityrhoPlate8000Interconnecting plate density
Types.SpecificConductivityLengthpreExp3.34e4Pre-exponential factor for the evaluation of membrane conductivity (original = 3.34e4)
Types.ElCurrent[N]iPENmaxzeros(N)Max current for each volume for leakage losses
Types.ElPotentialvLeak00Leakage loss at zero current
Types.LengthH1e-3Channel height
Types.LengthW0.05Channel width
Types.LengthtotLength0.05Channel total length
Types.LengthtauAE400e-6Anode electrode thickness
Types.LengthtauCE16e-6Cathode electrode thickness
Types.LengthtauSE4e-6Solid electrolyte thickness
Types.LengthtauI3e-3Interconnecting plate thickness
Modelica.Units.SI.FaradayConstantFModelica.Constants.FFaraday constant C/mol
Types.TemperatureT_start1000General start temperature
Types.Pressurep_start27.5e5
Types.Densityrho_start0.22
Types.CurrentDensityj_start1.5
Types.CurrentDensityj0Anode_start15000
Types.CurrentDensityj0Cathode_start6000
Types.MassFlowRatew_start1e-6
Types.Temperature[N + 1]T_start_Anode{T_start for i in 1:N + 1}
Types.Temperature[N + 1]T_start_Cathode{T_start for i in 1:N + 1}
Types.Temperature[N]T_start_PEN{T_start for i in 1:N}
Types.Temperature[N]T_start_Plate_Anode{T_start for i in 1:N}
Types.Temperature[N]T_start_Plate_Cathode{T_start for i in 1:N}
Types.PerUnit[N + 1,nX]X_start_Anodeones(N + 1, nX)/nX
Types.PerUnit[N + 1,nX]X_start_Cathodeones(N + 1, nX)/nX
Types.Pressure[N]p_start_Anodeones(N)*p_start
Types.Pressure[N]p_start_Cathodeones(N)*p_start
Types.Density[N + 1]rho_start_Anodeones(N + 1)*rho_start
Types.Density[N + 1]rho_start_Cathodeones(N + 1)*rho_start
Types.MassFlowRate[N + 1]w_start_Anodeones(N + 1)*w_start
Types.MassFlowRate[N + 1]w_start_Cathodeones(N + 1)*w_start
Types.CurrentDensity[N]j_startPENones(N)*j_start
Types.CurrentDensity[N]j0Anode_startPENones(N)*j0Anode_start
Types.CurrentDensity[N]j0Cathode_startPENones(N)*j0Cathode_start
Types.Pressurep_start_homotopy_cathode_inp_start_Cathode[1]
Types.Pressurep_start_homotopy_anode_inp_start_Anode[1]
Types.Pressurep_start_homotopy_cathode_outp_start_Cathode[end]
Types.Pressurep_start_homotopy_anode_outp_start_Anode[end]
Types.TemperatureT_start_homotopy_cathode_inT_start_Cathode[1]
Types.TemperatureT_start_homotopy_anode_inT_start_Anode[1]
Types.TemperatureT_start_homotopy_cathode_outT_start_Cathode[end]
Types.TemperatureT_start_homotopy_anode_outT_start_Anode[end]
Types.MassFraction[nX]X_start_homotopy_cathode_inX_start_Cathode[1, :]
Types.MassFraction[nX]X_start_homotopy_anode_inX_start_Anode[1, :]
Types.MassFraction[nX]X_start_homotopy_cathode_outX_start_Cathode[end, :]
Types.MassFraction[nX]X_start_homotopy_anode_outX_start_Anode[end, :]
Types.PressuredeltapNomAnode
Types.PressuredeltapNomCathode
Types.EfficiencyetaDCAC0.98
Types.MolarMass[10]MMDiffusionLosses.Utilities.ConstantsFluidAndSolid.MM
RealQlossFrac0

Connectors

TypeNameDefaultDescription
SOFCPoliMi.Interfaces.FlangeAanodeIn
SOFCPoliMi.Interfaces.FlangeBanodeOut
SOFCPoliMi.Interfaces.FlangeAcathodeIn
SOFCPoliMi.Interfaces.FlangeBcathodeOut
SOFCPoliMi.Interfaces.FuelCellInterfaces.ElPincathodePin
Interfaces.FuelCellInterfaces.ElPinanodePin
Modelica.Blocks.Interfaces.RealOutputTinAnode
Modelica.Blocks.Interfaces.RealOutputTinCathode
Modelica.Blocks.Interfaces.RealInputQloss

Components

TypeNameDefaultDescription
SOFCPoliMi.Components.FuelCell.Module[N]module
SOFCPoliMi.Components.FuelCell.FlowMultiplieranodeInletMultiplier
SOFCPoliMi.Components.FuelCell.FlowMultiplieranodeOutletMultiplier
SOFCPoliMi.Components.FuelCell.FlowMultipliercathodeInletMultiplier
SOFCPoliMi.Components.FuelCell.FlowMultipliercathodeOutletMultiplier
SOFCPoliMi.Components.FuelCell.VoltageMultipliervoltageMultiplier
Types.ElPotentialvStackanodePin.v - cathodePin.v
Types.ElPotentialvChannelmodule[1].anodeElPin.v - module[1].cathodeElPin.v
Types.PowerpElStackvStack*anodePin.i*etaDCAC
Types.PowerenBalance-Qloss*QlossFrac + NparallelFluid*(module[1].anodeChannel.entFlowIn + module[1].cathodeChannel.entFlowIn - (module[N].anodeChannel.entFlowOut + module[N].cathodeChannel.entFlowOut)) - pElStack
Types.CurrentDensitycurrentDensityactualCurrent/(totLength*W)
Types.ElCurrentactualCurrentabs(voltageMultiplier.single.i)Single channel current
Types.MoleFraction[module[1].anodeChannel.nX]XinAnodemodule[1].anodeChannel.Xin
Types.ElCurrentwholeConsCurrent2*F*module[1].anodeChannel.wIn*(XinAnode[2]/MM[2] + XinAnode[4]/MM[4] + 4*XinAnode[5]/MM[5] + 7*XinAnode[6]/MM[6] + 10*XinAnode[7]/MM[7])
Types.PerUnitUfactualCurrent/wholeConsCurrent
Types.MassFlowRatetotalWO2sum(module[i].anodeChannel.wO2 for i in 1:N)*NparallelFluid
SOFCPoliMi.Components.Initializers.HomotopyInitializerXhomotopyInitializerCathodeIn
SOFCPoliMi.Components.Initializers.HomotopyInitializerXhomotopyInitializerAnodeIn

Contents

NameDescription
Fluid
ReactionRates
PENOhmRes