Settings modules

This modul contains some of the required functions for the settings.

EIS_parameters(f_EIS)

This function gives the time parameters for the EIS_current density function.

Parameters:
  • f_EIS (tuple) –

    EIS parameters. It is a tuple containing the power of the initial frequency 'f_power_min_EIS': f_min_EIS = 10**f_power_min_EIS, the power of the final frequency 'f_power_max_EIS', the number of frequencies tested 'nb_f_EIS' and the number of points calculated per specific period 'nb_points_EIS'.
Returns:
  • t_EIS( tuple ) –

    EIS parameters. It is a tuple containing the initial EIS time after stack equilibrium 't0_EIS', a list of time parameters which gives the beginning of each frequency change 't_new_start_EIS', the final time 'tf_EIS', a list of time parameters which gives the estimated time for reaching equilibrium at each frequency 'delta_t_break_EIS', and a list of time parameters which gives the estimated time for measuring the voltage response at each frequency 'delta_t_measurement_EIS'.
Source code in modules/settings_modules.py 
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
def EIS_parameters(f_EIS):
    """This function gives the time parameters for the EIS_current density function.

    Parameters
    ----------
    f_EIS : tuple
        EIS parameters. It is a tuple containing the power of the initial frequency 'f_power_min_EIS':
        f_min_EIS = 10**f_power_min_EIS, the power of the final frequency 'f_power_max_EIS', the number of frequencies
        tested 'nb_f_EIS' and the number of points calculated per specific period 'nb_points_EIS'.

    Returns
    -------
    t_EIS : tuple
        EIS parameters. It is a tuple containing the initial EIS time after stack equilibrium 't0_EIS', a list of time
        parameters which gives the beginning of each frequency change 't_new_start_EIS', the final time 'tf_EIS', a list
        of time parameters which gives the estimated time for reaching equilibrium at each frequency
        'delta_t_break_EIS', and a list of time parameters which gives the estimated time for measuring the voltage
        response at each frequency 'delta_t_measurement_EIS'.
    """

    # Initialisation
    #       Frequencies
    f_power_min_EIS, f_power_max_EIS, nb_f_EIS, nb_points_EIS = f_EIS  # They are the frequency parameters for the EIS
    #                                                                    simulation.
    f = np.logspace(f_power_min_EIS, f_power_max_EIS, num=nb_f_EIS)  # It is the tested frequencies
    nb_period_break_EIS, nb_period_measurement_EIS = 50, 50  # They are the number of temporal periods which are used
    #                                                          for break and for measurement. It is more accurate to use
    #                                                          periods than time as the frequency range is big.
    #       Time parameters
    delta_t_break_EIS = np.array([])  # It is the estimated time for reaching equilibrium at each frequency.
    delta_t_measurement_EIS = np.array([])  # It is the estimated time for measuring the voltage response.

    # Time parameters calculation
    t0_EIS = 120*60  # s. It is the simulation starting time. [0, t0_EIS] is used to let the stack equilibrate to i_EIS.
    t_new_start_EIS = np.array([t0_EIS])  # It is a list of time parameters which gives the beginning of each frequency
    #                                   change.
    for i in range(nb_f_EIS):  # The goal is to measure nb_f_EIS periods of the signal in order to have precise enough values.
        T_i = 1 / (f[i]) # s. It is the period of the signal.
        if i < (nb_f_EIS - 1):
            delta_t_break_EIS = np.concatenate((delta_t_break_EIS, [nb_period_break_EIS * T_i]))
            delta_t_measurement_EIS = np.concatenate((delta_t_measurement_EIS, [nb_period_measurement_EIS * T_i]))
            next_start_EIS = t_new_start_EIS[i] + delta_t_break_EIS[i] + delta_t_measurement_EIS[i]
            t_new_start_EIS = np.concatenate((t_new_start_EIS, [next_start_EIS]))
        else:
            delta_t_break_EIS = np.concatenate((delta_t_break_EIS, [nb_period_break_EIS * T_i]))
            delta_t_measurement_EIS = np.concatenate((delta_t_measurement_EIS, [nb_period_measurement_EIS * T_i]))
            tf_EIS = t_new_start_EIS[-1] + delta_t_break_EIS[-1] + delta_t_measurement_EIS[-1]  # s. It is the
            #                                                                                     simulation ending time

    t_EIS = t0_EIS, t_new_start_EIS, tf_EIS, delta_t_break_EIS, delta_t_measurement_EIS
    return t_EIS

stored_operating_inputs(type_fuel_cell, voltage_zone)

This function gives the operating inputs which correspond to the given type_fuel_cell.

Parameters:
  • type_fuel_cell (str) –

    Type of fuel cell configuration.

  • voltage_zone (str) –

    Zone of the polarization curve which is considered.
Returns:
  • T_des( float ) –

    Desired fuel cell temperature in Kelvin.

  • Pa_des( float ) –

    Desired anode pressure in Pascal.

  • Pc_des( float ) –

    Desired cathode pressure in Pascal.

  • Sa( float ) –

    Stoichiometric ratio of hydrogen.

  • Sc( float ) –

    Stoichiometric ratio of oxygen.

  • Phi_a_des( float ) –

    Desired anode relative humidity.

  • Phi_c_des( float ) –

    Desired cathode relative humidity.

  • i_max_pola( float ) –

    Maximum current density for the polarization curve.
Source code in modules/settings_modules.py 
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
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
def stored_operating_inputs(type_fuel_cell, voltage_zone):
    """This function gives the operating inputs which correspond to the given type_fuel_cell.

    Parameters
    ----------
    type_fuel_cell : str
        Type of fuel cell configuration.
    voltage_zone : str
        Zone of the polarization curve which is considered.

    Returns
    -------
    T_des : float
        Desired fuel cell temperature in Kelvin.
    Pa_des : float
        Desired anode pressure in Pascal.
    Pc_des : float
        Desired cathode pressure in Pascal.
    Sa : float
        Stoichiometric ratio of hydrogen.
    Sc : float
        Stoichiometric ratio of oxygen.
    Phi_a_des : float
        Desired anode relative humidity.
    Phi_c_des : float
        Desired cathode relative humidity.
    i_max_pola : float
        Maximum current density for the polarization curve.
    """

    # For the ZSW Generic Stack fuel cell
    if type_fuel_cell == "ZSW-GenStack":
        T_des = 68 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.2e5, 2.0e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.6, 1.6  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.398, 0.50  # It is the desired relative humidity.
        y_H2_in = 0.7  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.500e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.700e4 # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "ZSW-GenStack_Pa_1.61_Pc_1.41":
        T_des = 68 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 1.61e5, 1.41e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.6, 1.6  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.398, 0.50  # It is the desired relative humidity.
        y_H2_in = 0.7 # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.200e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 0.700e4 # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "ZSW-GenStack_Pa_2.01_Pc_1.81":
        T_des = 68 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.01e5, 1.81e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.6, 1.6  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.398, 0.50  # It is the desired relative humidity.
        y_H2_in = 0.7  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.500e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.300e4  # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "ZSW-GenStack_Pa_2.4_Pc_2.2":
        T_des = 68 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.4e5, 2.2e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.6, 1.6  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.398, 0.50  # It is the desired relative humidity.
        y_H2_in = 0.7  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.500e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.900e4  # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "ZSW-GenStack_Pa_2.8_Pc_2.6":
        T_des = 68 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.8e5, 2.6e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.6, 1.6  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.398, 0.50  # It is the desired relative humidity.
        y_H2_in = 0.7  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.500e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.900e4  # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "ZSW-GenStack_T_62":
        T_des = 62 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.2e5, 2.0e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.6, 1.6  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.398, 0.50  # It is the desired relative humidity.
        y_H2_in = 0.7  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.500e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.500e4  # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "ZSW-GenStack_T_76":
        T_des = 76 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.2e5, 2.0e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.6, 1.6  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.398, 0.50  # It is the desired relative humidity.
        y_H2_in = 0.7  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.500e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.100e4  # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "ZSW-GenStack_T_84":
        T_des = 84 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.2e5, 2.0e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.6, 1.6  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.398, 0.50  # It is the desired relative humidity.
        y_H2_in = 0.7  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.000e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 0.700e4  # A.m-2. It is the maximum current density for the polarization curve.

    # For EH-31 fuel cell
    elif type_fuel_cell == "EH-31_1.5":
        T_des = 74 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 1.5e5, 1.5e5  # Pa. It is the desired pressure of the fuel gas (at the anode/cathode).
        Sa, Sc = 1.2, 2.0  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.4, 0.6  # It is the desired relative humidity.
        y_H2_in = 1 # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.300e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.700e4  # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "EH-31_2.0":
        T_des = 74 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.0e5, 2.0e5  # Pa. It is the desired pressure of the fuel gas (at the anode/cathode).
        Sa, Sc = 1.2, 2.0  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.4, 0.6  # It is the desired relative humidity.
        y_H2_in = 1  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.500e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.300e4  # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "EH-31_2.25":
        T_des = 74 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.25e5, 2.25e5  # Pa. It is the desired pressure of the fuel gas (at the anode/cathode).
        Sa, Sc = 1.2, 2.0  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.4, 0.6  # It is the desired relative humidity.
        y_H2_in = 1  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 2.800e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.700e4  # A.m-2. It is the maximum current density for the polarization curve.
    elif type_fuel_cell == "EH-31_2.5":
        T_des = 74 + 273.15  # K. It is the desired fuel cell temperature.
        Pa_des, Pc_des = 2.5e5, 2.5e5  # Pa. It is the desired pressures of the fuel gas.
        Sa, Sc = 1.2, 2.0  # It is the stoichiometric ratio (of hydrogen and oxygen).
        Phi_a_des, Phi_c_des = 0.4, 0.6  # It is the desired relative humidity.
        y_H2_in = 1  # It is the molar fraction of H2 in the dry anode gas mixture (H2/N2) injected at the inlet.
        if voltage_zone == "full":
            i_max_pola = 3.000e4  # A.m-2. It is the maximum current density for the polarization curve.
        elif voltage_zone == "before_voltage_drop":
            i_max_pola = 1.600e4  # A.m-2. It is the maximum current density for the polarization curve.

    # For other fuel cells
    else:
        raise ValueError('the type_fuel_cell given is not valid.')

    return T_des, Pa_des, Pc_des, Sa, Sc, Phi_a_des, Phi_c_des, y_H2_in, i_max_pola

stored_physical_parameters(type_fuel_cell)

This function gives the physical parameters which correspond to the given type_fuel_cell.

Parameters:
  • type_fuel_cell (str) –

    Type of fuel cell configuration.
Returns:
  • Hacl( float ) –

    Thickness of the anode catalyst layer in m.

  • Hccl( float ) –

    Thickness of the cathode catalyst layer in m.

  • epsilon_mc( float ) –

    Volume fraction of ionomer in the CL.

  • Hmem( float ) –

    Thickness of the membrane in m.

  • Hgdl( float ) –

    Thickness of the gas diffusion layer in m.

  • epsilon_gdl( float ) –

    Anode/cathode GDL porosity.

  • epsilon_c( float ) –

    Compression ratio of the GDL.

  • Hagc( float ) –

    Thickness of the anode gas channel in m.

  • Hcgc( float ) –

    Thickness of the cathode gas channel in m.

  • Wagc( float ) –

    Width of the anode gas channel in m.

  • Wcgc( float ) –

    Width of the cathode gas channel in m.

  • Lgc( float ) –

    Length of the gas channel in m.

  • Aact( float ) –

    Active area of the cell in m².

  • e( float ) –

    Capillary exponent.

  • 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.
Source code in modules/settings_modules.py 
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
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
def stored_physical_parameters(type_fuel_cell):
    """This function gives the physical parameters which correspond to the given type_fuel_cell.

    Parameters
    ----------
    type_fuel_cell : str
        Type of fuel cell configuration.

    Returns
    -------
    Hacl : float
        Thickness of the anode catalyst layer in m.
    Hccl : float
        Thickness of the cathode catalyst layer in m.
    epsilon_mc : float
        Volume fraction of ionomer in the CL.
    Hmem : float
            Thickness of the membrane in m.
    Hgdl : float
            Thickness of the gas diffusion layer in m.
    epsilon_gdl : float
        Anode/cathode GDL porosity.
    epsilon_c : float
        Compression ratio of the GDL.
    Hagc : float
        Thickness of the anode gas channel in m.
    Hcgc : float
        Thickness of the cathode gas channel in m.
    Wagc : float
        Width of the anode gas channel in m.
    Wcgc : float
        Width of the cathode gas channel in m.
    Lgc : float
        Length of the gas channel in m.
    Aact : float
            Active area of the cell in m².
    e : float
        Capillary exponent.
    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.
    """

    # For the ZSW Generic Stack fuel cell
    if type_fuel_cell == "ZSW-GenStack" or type_fuel_cell == "ZSW-GenStack_Pa_1.61_Pc_1.41" or \
            type_fuel_cell == "ZSW-GenStack_Pa_2.01_Pc_1.81" or type_fuel_cell == "ZSW-GenStack_Pa_2.4_Pc_2.2" or \
            type_fuel_cell == "ZSW-GenStack_Pa_2.8_Pc_2.6" or type_fuel_cell == "ZSW-GenStack_T_62" or \
            type_fuel_cell == "ZSW-GenStack_T_76" or type_fuel_cell == "ZSW-GenStack_T_84":
        # Global
        Aact = 283.87e-4  # m². It is the MEA active area.
        nb_cell = 26  # . It is the number of cell in the stack.
        # Catalyst layer
        Hacl = 8e-6  # m. It is the thickness of the anode catalyst layer.
        Hccl = 17e-6  # m. It is the thickness of the cathode catalyst layer.
        epsilon_mc = 0.25  # It is the volume fraction of ionomer in the CL.
        # Membrane
        Hmem = 15e-6  # m. It is the thickness of the membrane.
        # Gas diffusion layer
        Hgdl = 127e-6  # m. It is the thickness of the gas diffusion layer.
        epsilon_gdl = 0.788  # It is the anode/cathode GDL porosity.
        epsilon_cl = 0.5  # It is the porosity of the catalyst layer, without units.
        epsilon_c = 0.2  # It is the compression ratio of the GDL.
        #   Microporous layer
        Hmpl = 70e-6  # m. It is the thickness of the microporous layer.
        epsilon_mpl = 0.425  # It is the porosity of the microporous layer.
        # Gas channel
        Hagc = 230e-6  # m. It is the thickness of the anode gas channel.
        Hcgc = 300e-6  # m. It is the thickness of the cathode gas channel.
        Wagc = 430e-6  # m. It is the width of the anode gas channel.
        Wcgc = 532e-6  # m. It is the width of the cathode gas channel.
        Lgc = 246.2e-3  # m. It is the length of one channel in the bipolar plate.
        nb_channel_in_gc = 105 # . It is the number of channels in the bipolar plate.
        Ldist = 7.11e-2 # m. It is the length of the distributor, which is the volume between the gas channel and the manifold.
        #   Auxiliaries
        Lm = 25.8e-3  # m. It is the length of the manifold.
        A_T_a = 9.01e-4  # m². It is the inlet/exhaust anode manifold throttle area
        A_T_c = 22.61e-4  # m². It is the inlet/exhaust cathode manifold throttle area
        Vasm, Vcsm = Lm * A_T_a, Lm * A_T_c  # m3. It is the supply manifold volume.
        Vaem, Vcem = Vasm, Vcsm  # m-3. It is the exhaust manifold volume.
        # Interaction parameters between water and PEMFC structure
        e = 4.0  # It is the capillary exponent
        # Voltage polarization
        Re = 1e-06  # ohm.m². It is the electron conduction resistance of the circuit.
        i0_d_c_ref = 14.43  # A.m-2. It is the dry reference exchange current density at the cathode.
        i0_h_c_ref = 1.0  # A.m-2. It is the fully humidified reference exchange current density at the cathode.
        kappa_co = 5 # mol.m-1.s-1.Pa-1. It is the crossover correction coefficient.
        kappa_c = 1.026  # It is the overpotential correction exponent.
        a_slim, b_slim, a_switch = 0.05553, 0.10514, 0.63654  # It is the limit liquid saturation coefficients.
        C_scl = 2e7  # F.m-3. It is the volumetric space-charge layer capacitance.

    # For EH-31 fuel cell
    elif type_fuel_cell == "EH-31_1.5" or type_fuel_cell == "EH-31_2.0" or type_fuel_cell == "EH-31_2.25" or \
            type_fuel_cell == "EH-31_2.5":
        # Global
        Aact = 85e-4  # m². It is the active area of the catalyst layer.
        nb_cell = 1  # . It is the number of cell in the stack.
        # Catalyst layer
        Hacl = 8.593e-6  # m. It is the thickness of the anode catalyst layer.
        Hccl = Hacl  # m. It is the thickness of the cathode catalyst layer.
        epsilon_mc = 0.3986  # It is the volume fraction of ionomer in the CL.
        # Membrane
        Hmem = 16.06e-6  # m. It is the thickness of the membrane.
        # Gas diffusion layer
        Hgdl = 200e-6  # m. It is the thickness of the gas diffusion layer.
        epsilon_gdl = 0.5002  # It is the anode/cathode GDL porosity.
        epsilon_cl = 0.25  # It is the porosity of the catalyst layer, without units.
        epsilon_c = 0.2  # It is the compression ratio of the GDL.
        #   Microporous layer
        Hmpl = 30e-6  # m. It is the thickness of the microporous layer.
        epsilon_mpl = 0.4  # It is the porosity of the microporous layer.
        # Gas channel
        Hagc = 500e-6  # m. It is the thickness of the anode gas channel.
        Hcgc = Hagc  # m. It is the thickness of the cathode gas channel.
        Wagc = 450e-6  # m. It is the width of the anode gas channel.
        Wcgc = Wagc  # m. It is the width of the cathode gas channel.
        Lgc = 144e-3  # m. It is the length of one channel in the bipolar plate.
        nb_channel_in_gc = 67  # . It is the number of channels in the bipolar plate.
        Ldist = 5e-2  # m. It is the estimated length of the distributor, which is the volume between the gas channel and the manifold.
        #   Auxiliaries
        Lm = 2.03  # m. It is the length of the manifold.
        A_T_a = 11.8e-4  # m². It is the inlet/exhaust anode manifold throttle area
        A_T_c = 34.4e-4  # m². It is the inlet/exhaust cathode manifold throttle area
        Vasm, Vcsm = Lm * A_T_a, Lm * A_T_c  # m3. It is the supply manifold volume.
        Vaem, Vcem = Vasm, Vcsm  # m-3. It is the exhaust manifold volume.
        # Interaction parameters between water and PEMFC structure
        e = 4.0  # It is the capillary exponent
        # Voltage polarization
        Re = 1e-06  # ohm.m². It is the electron conduction resistance of the circuit.
        i0_d_c_ref = 14.43  # A.m-2. It is the dry reference exchange current density at the cathode.
        i0_h_c_ref = 1.0  # A.m-2. It is the fully humidified reference exchange current density at the cathode.
        kappa_co = 30.42  # mol.m-1.s-1.Pa-1. It is the crossover correction coefficient.
        kappa_c = 0.4152  # It is the overpotential correction exponent.
        a_slim, b_slim, a_switch = 0.05553, 0.10514, 0.82  # It is the limit liquid saturation coefficients.
        C_scl = 20e6  # F.m-3. It is the volumetric space-charge layer capacitance.

    # For other fuel cells
    else:
        raise ValueError('the type_input given is not valid.')

    return (Hacl, Hccl, epsilon_mc, Hmem, Hgdl, epsilon_gdl, epsilon_cl, epsilon_c, Hmpl, epsilon_mpl, Hagc, Hcgc, Wagc,
            Wcgc, Lgc, nb_channel_in_gc, Ldist, Lm, A_T_a, A_T_c, Vasm, Vcsm, Vaem, Vcem, Aact, nb_cell, e, Re, i0_d_c_ref,
            i0_h_c_ref, kappa_co, kappa_c, a_slim, b_slim, a_switch, C_scl)