modelTestStand
Simulate the fuel cell under prescribed conditions
Extends from Modelica.Icons.Example (Icon for runnable examples).
Information
Please see TestConditions regarding how the properties of the environment are used.
To run the cell with pure O2 in the cathode (no N2), set environment.psi_O2_dry to 100%.
Assumptions:
- The outer surface of each end plate has uniform temperature in the yz plane.
- No heat is conducted from the rest of the cell hardware.
- All electronic current passes through the first segment (index
[1, 1]) of each end plate in the yz plane. - There is no shear force on the fluid at either outlet.
- The pressure of each gas species is uniform over each inlet.
- The volumetric flow rates of the gas species are equal at each outlet, where the volumetric flow rate is approximated using the current at the outlet and the density of the gas within the last segment of the cell.
- The outlet pressure is applied to the gas mixture by Dalton's law (additivity of pressure).
- At the outlet, the liquid has the same pressure as the gas (Amagat's law).
- There is no thermal conduction across either outlet.
Components
| Type | Name | Default | Description |
|---|---|---|---|
| Q.Current | zI | Electrical current | |
| Q.Potential | w | load.v*U.V | Potential |
| Q.CurrentAreic | J | zI/cell.caFP.A[Axis.x] | Current density |
| Q.Number | J_Apercm2 | J*U.cm^2/U.A | Current density, in A/cm2 |
| Q.PressureAbsolute | p_an_in | average(average(anSource.gas.H2.boundary.p + anSource.gas.H2O.boundary.p)) | Total pressure at the anode inlet |
| Q.PressureAbsolute | p_ca_in | average(average(caSource.gas.H2O.boundary.p + p_N2_in + caSource.gas.O2.boundary.p)) | Total pressure at the cathode inlet |
| Q.Pressure | Deltap_an | environment.p - p_an_in | Pressure difference down the anode channel |
| Q.Pressure | Deltap_ca | environment.p - p_ca_in | Pressure difference down the cathode channel |
| Q.Potential | Deltaw_O2 | (DataO2.g(caSource[1, 1].gas.O2.boundary.T, caSource[1, 1].gas.O2.boundary.p) - cell.caCL.subregions[1, 1, 1].gas.O2.g)/4 | Voltage loss due to O2 supply |
| Q.Potential | Deltaw_H2O | (DataH2O.g(caSource[1, 1].gas.H2O.boundary.T, caSource[1, 1].gas.H2O.boundary.p) - cell.caCL.subregions[1, 1, 1].gas.H2O.g)/2 | Voltage loss due to H2O removal |
| Q.Potential | Deltaw_H2 | (DataH2.g(anSource[1, 1].gas.H2.boundary.T, anSource[1, 1].gas.H2.boundary.p) - cell.anCL.subregions[1, 1, 1].gas.H2.g)/2 | Voltage loss due to H2 supply |
| Q.Potential | 'Deltaw_e-' | cell.caFP.subregions[1, 1, 1].graphite.'e-'.g_boundaries[1, Side.p] - cell.caCL.subregions[1, 1, 1].graphite.'e-'.g + cell.anCL.subregions[1, 1, 1].graphite.'e-'.g - cell.anFP.subregions[1, 1, 1].graphite.'e-'.g_boundaries[1, Side.n] | Voltage loss due to e- transport |
| Q.Potential | 'Deltaw_H+' | cell.anCL.subregions[1, 1, 1].ionomer.'H+'.g - cell.caCL.subregions[1, 1, 1].ionomer.'H+'.g | Voltage loss due to H+ transport |
| Q.Potential | Deltaw_an | cell.anCL.subregions[1, 1, 1].graphite.'e-Transfer'.Deltag | Anode overpotential |
| Q.Potential | Deltaw_ca | -cell.caCL.subregions[1, 1, 1].graphite.'e-Transfer'.Deltag | Cathode overpotential |
| Cell | cell | ||
| Conditions.ByConnector.BoundaryBus.Single.Sink[cell.n_y,cell.n_z] | anBC | Boundary condition for the anode end plate, except electrical | |
| Conditions.ByConnector.BoundaryBus.Single.Source[cell.anFP.n_x,cell.n_z] | anSource | Source for the anode reactant stream | |
| Conditions.ByConnector.BoundaryBus.Single.Sink[cell.anFP.n_x,cell.n_z] | anSink | Sink for the anode reactant stream | |
| Conditions.ByConnector.BoundaryBus.Single.Source[cell.n_y,cell.n_z] | caBC | Boundary condition for the cathode end plate, except electrical | |
| Conditions.ByConnector.BoundaryBus.Single.Source[cell.caFP.n_x,cell.n_z] | caSource | Source for the cathode reactant stream | |
| Conditions.ByConnector.BoundaryBus.Single.Sink[cell.caFP.n_x,cell.n_z] | caSink | Sink for the cathode reactant stream | |
| Conditions.Environment | environment | Environmental conditions | |
| Modelica.Electrical.Analog.Sources.RampCurrent | load | ||
| Modelica.Electrical.Analog.Basic.Ground | ground | ||
| Conditions.Adapters.MSL.Electronic | anAdapt | Interface with Modelica.Electrical on the anode side | |
| Conditions.Adapters.MSL.Electronic | caAdapt | Interface with Modelica.Electrical on the cathode side | |
| Conditions.Router[cell.anFP.n_x,cell.n_z] | anRouter | Switch to route the anode for reverse flow | |
| Conditions.Router[cell.anFP.n_x,cell.n_z] | caRouter | Switch to route the cathode for reverse flow | |
| TestConditions | testConditions | Test conditions |