114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199 | def physical_parameters(type_fuel_cell):
"""This function is used to set the physical parameters of the fuel cell system.
Parameters
----------
type_fuel_cell : str
Type of fuel cell system. It can be "EH-31_1.5", "EH-31_2.0", "EH-31_2.25", "EH-31_2.5", "LF",
or "manual_setup".
Returns
-------
Hcl : float
Thickness of the anode or cathode catalyst layer in meters.
epsilon_mc : float
Volume fraction of ionomer in the catalyst layer.
tau : float
Pore structure coefficient in the CL.
Hmem : float
Thickness of the membrane in meters.
Hgdl : float
Thickness of the gas diffusion layer in meters.
epsilon_gdl : float
Anode/cathode GDL porosity.
epsilon_c : float
Compression ratio of the GDL.
Hgc : float
Thickness of the gas channel in meters.
Wgc : float
Width of the gas channel in meters.
Lgc : float
Length of the gas channel in meters.
Aact : float
Active area of the catalyst layer in meters squared.
e : float
Capillary exponent.
Re : float
Electron conduction resistance of the circuit in ohm.m².
i0_c_ref : float
Reference exchange current density at the cathode in A.m-2.
kappa_co : float
Crossover correction coefficient in mol.m-1.s-1.Pa-1.
kappa_c : float
Overpotential correction exponent.
a_slim : float
One of the limit liquid saturation coefficients: the slop of slim function.
b_slim : float
One of the limit liquid saturation coefficients: the intercept of slim function.
a_switch : float
One of the limit liquid saturation coefficients: the slop of s_switch function.
C_dl : float
Volumetric double layer capacitance in F.m-3.
"""
if type_fuel_cell == "manual_setup": # Setup which are not stored in "stored_physical_parameters".
# Fuel cell physical parameters: 𝜔 (which are not controllable by the system)
# Catalyst layer
Aact = 8.5e-3 # m². It is the active area of the catalyst layer.
Hcl = 1e-5 # m. It is the thickness of the anode or cathode catalyst layer.
epsilon_mc = 0.3949198274842546 # It is the volume fraction of ionomer in the CL.
tau = 1.015639135686993 # It is the pore structure coefficient in the CL, without units.
# Membrane
Hmem = 2e-5 # m. It is the thickness of the membrane.
# Gas diffusion layer
Hgdl = 2e-4 # m. It is the thickness of the gas diffusion layer.
epsilon_gdl = 0.7011156494971454 # It is the anode/cathode GDL porosity.
epsilon_c = 0.27052745219052654 # It is the compression ratio of the GDL.
# Gas channel
Hgc = 5e-4 # m. It is the thickness of the gas channel.
Wgc = 4.5e-4 # m. It is the width of the gas channel.
Lgc = 9.67 # m. It is the length of the gas channel.
# Interaction parameters between water and PEMFC structure
e = 5.0 # It is the capillary exponent
# Voltage polarization
Re = 5.694464714060734e-07 # ohm.m². It is the electron conduction resistance of the circuit.
i0_c_ref = 2.787917581303015 # A.m-2. It is the reference exchange current density at the cathode.
kappa_co = 29.793535549174077 # mol.m-1.s-1.Pa-1. It is the crossover correction coefficient.
kappa_c = 1.6136446641573106 # It is the overpotential correction exponent.
a_slim, b_slim, a_switch = 0.0555312850726664, 0.10514269908118055, 0.6365424991141914 # It is the limit
# liquid saturation coefficients.
C_scl = 2e7 # F.m-3. It is the volumetric space-charge layer capacitance.
else: # Stored setup in "stored_physical_parameters".
(Hcl, epsilon_mc, tau, Hmem, Hgdl, epsilon_gdl, epsilon_c, Hgc, Wgc, Lgc, Aact, e, Re, i0_c_ref, kappa_co,
kappa_c, a_slim, b_slim, a_switch, C_scl) = stored_physical_parameters(type_fuel_cell)
return (Hcl, epsilon_mc, tau, Hmem, Hgdl, epsilon_gdl, epsilon_c, Hgc, Wgc, Lgc, Aact, e, Re, i0_c_ref, kappa_co,
kappa_c, a_slim, b_slim, a_switch, C_scl)
|