A global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates

Abstract

The present invention discloses a global optimization design and comprehensive evaluation method for a fuel cell metal bipolar plate, which relates to the technical field of fuel cell design, and includes: determining the hydrogen molar flow rate, air molar flow rate and heat flow rate of the metal bipolar plate according to the working parameters of the metal bipolar plate; determining the initial geometric parameters of the metal bipolar plate according to the hydrogen molar flow rate, air molar flow rate and heat flow rate; determining the initialization bipolar plate configuration of the metal bipolar plate according to the initial geometric parameters; establishing a three-dimensional model according to the initialization bipolar plate configuration, performing simulation, obtaining simulation results, and determining evaluation indicators based on the simulation results; optimizing the geometric parameters of the metal bipolar plate according to the evaluation indicators to obtain the final geometric parameters of the metal bipolar plate; and determining the final bipolar plate configuration of the metal bipolar plate based on the final geometric parameters, so as to achieve optimal performance.

Description

Technical Field

[0001] The present invention relates to the technical field of fuel cell design, and in particular to a global optimization design and comprehensive evaluation method, device, computer equipment and storage medium for a fuel cell metal bipolar plate. Background Art

[0002] A fuel cell is a power generation device that converts chemical energy in fuel into electrical energy through redox reactions using oxygen or other oxidants. The most common fuel is hydrogen. In the prior art, the design of metal bipolar plates for fuel cells is mainly based on target process requirements and user experience. The above design method is difficult to achieve optimal performance. Summary of the invention

[0003] The embodiment of the present invention provides a global optimization design and comprehensive evaluation method for a fuel cell metal bipolar plate, aiming to solve the problem that the existing design method of the fuel cell metal bipolar plate cannot achieve the optimal performance.

[0004] The embodiment of the present invention provides a global optimization design and comprehensive evaluation method for a fuel cell metal bipolar plate, which includes:

[0005] Obtaining operating parameters of a metal bipolar plate of a fuel cell;

[0006] Determining the hydrogen molar flow rate, air molar flow rate and heat flow rate of the metal bipolar plate according to the operating parameters;

[0007] Determining initial geometric parameters of the metal bipolar plate according to the hydrogen molar flow rate, the air molar flow rate and the heat flow rate;

[0008] Determining an initialization bipolar plate configuration of the metal bipolar plate according to initial geometric parameters of the metal bipolar plate;

[0009] Establishing a three-dimensional model according to the initialized bipolar plate configuration, performing simulation, obtaining simulation results, and determining evaluation indicators based on the simulation results;

[0010] Optimizing the geometric parameters of the metal bipolar plate according to the evaluation index to obtain the final geometric parameters of the metal bipolar plate;

[0011] A final bipolar plate configuration of the metal bipolar plate is determined based on the final geometric parameters.

[0012] A further technical solution is that the working parameters include working current density, hydrogen excess coefficient, reaction active area and number of single cells, and the hydrogen molar flow rate is expressed by the following formula: Calculate, where m1 is the hydrogen molar flow rate, i is the operating current density, N0 is the number of cells, Aact is the reaction active area, S1 is the hydrogen excess coefficient and F is the Faraday constant.

[0013] A further technical solution is that the working parameters also include an oxygen excess coefficient, and the air molar flow rate is expressed by the following formula: Calculate, where S2 is the oxygen excess coefficient.

[0014] A further technical solution is that the working parameters also include the active area of ​​the single cell and the battery voltage, and the heat flow is calculated by the following formula Q=m1ΔH-iAV, wherein ΔH=286KJ / mol, A is the active area of ​​the single cell, and V is the battery voltage.

[0015] A further technical solution is that the initial geometric parameters of the metal bipolar plate include initial design parameters of the hydrogen flow channel and the hydrogen manifold of the metal bipolar plate, wherein the initial design parameters of the hydrogen flow channel and the hydrogen manifold are determined according to the hydrogen molar flow rate.

[0016] A further technical solution is that the initial geometric parameters of the metal bipolar plate also include initial design parameters of the air flow channel and the air manifold of the metal bipolar plate, wherein the initial design parameters of the air flow channel and the air manifold are determined according to the air molar flow rate.

[0017] A further technical solution is that the initial geometric parameters of the metal bipolar plate also include initial design parameters of a distribution area of ​​the metal bipolar plate, wherein the initial design parameters of the distribution area are determined based on the hydrogen molar flow rate and the air molar flow rate.

[0018] A further technical solution is that the initial geometric parameters of the metal bipolar plate also include initial design parameters of the cooling channel of the metal bipolar plate, wherein the initial design parameters of the cooling channel are determined according to the heat flow.

[0019] A further technical solution is that the three-dimensional model is established according to the initialization bipolar plate configuration, including:

[0020] Solid meshing is performed according to the initialized bipolar plate configuration to establish a finite element model of the bipolar plate flow channel structure of the hydrogen fuel cell;

[0021] A three-dimensional two-phase flow model of the bipolar plate flow channel structure of a hydrogen fuel cell is established by performing fluid grid division according to the initialized bipolar plate configuration.

[0022] A further technical solution is that the geometric parameters of the metal bipolar plate are optimized according to the evaluation index to obtain the final geometric parameters of the metal bipolar plate, including:

[0023] Taking the evaluation index as the optimization target and the geometric parameters of the metal bipolar plate as the decision variables, a multi-objective optimization algorithm is used to perform multi-objective optimization on the configuration of the metal bipolar plate. When the evaluation index is optimal, an optimal geometric parameter solution set is obtained, and the optimal geometric parameter solution set is used as the final geometric parameters of the metal bipolar plate.

[0024] For example, taking the effective mass transfer coefficient, the coefficient of variation, the synergy angle and the net mass transfer power as optimization targets, and the geometric parameters of the metal bipolar plate as decision variables, the NSGA-III algorithm is used to perform multi-objective optimization on the configuration of the metal bipolar plate. When the effective mass transfer coefficient, the coefficient of variation, the synergy angle and the net mass transfer power are optimized, an optimal geometric parameter solution set is obtained, and the optimal geometric parameter solution set is used as the final geometric parameters of the metal bipolar plate.

[0025] Wherein, the effective mass transfer coefficient is expressed by the following formula: Calculation, where EMTC is the effective mass transfer coefficient, v and are the components of velocity and concentration gradient in the normal direction of the proton exchange membrane, respectively.

[0026] The coefficient of variation is given by the following formula Calculation; where N is the number of flow channels, q k is the inlet mass flow rate of the kth flow channel, is the average mass flow rate of the flow channel.

[0027] The synergy angle is based on the following formula Calculate, where c is the concentration.

[0028] The net mass transfer power is given by the following formula Calculate, where V is the battery voltage, i is the current, ΔP is the loss along the flow channel from the inlet to the outlet, and A is the active area of ​​the single cell.

[0029] An embodiment of the present invention provides a global optimization design and comprehensive evaluation method for a metal bipolar plate of a fuel cell, the method comprising: obtaining the operating parameters of the metal bipolar plate of the fuel cell; determining the hydrogen molar flow rate, air molar flow rate and heat flow rate of the metal bipolar plate according to the operating parameters; determining the initial geometric parameters of the metal bipolar plate according to the hydrogen molar flow rate, the air molar flow rate and the heat flow rate; determining the initialization bipolar plate configuration of the metal bipolar plate according to the initial geometric parameters of the metal bipolar plate; establishing a three-dimensional model according to the initialization bipolar plate configuration, performing simulation, obtaining simulation results, and determining evaluation indicators based on the simulation results; optimizing the geometric parameters of the metal bipolar plate according to the evaluation indicators to obtain the final geometric parameters of the metal bipolar plate; determining the final bipolar plate configuration of the metal bipolar plate based on the final geometric parameters, so as to obtain the optimally designed fuel cell metal bipolar plate, thereby ensuring that the best performance can be achieved. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying any creative work.

[0031] Figure 1 A schematic flow chart of a global optimization design and comprehensive evaluation method for a fuel cell metal bipolar plate provided by an embodiment of the present invention;

[0032] Figure 2 A schematic diagram of the principle of the global optimization design and comprehensive evaluation method of the metal bipolar plate of the fuel cell provided by the embodiment of the present invention;

[0033] Figure 3a A schematic diagram of the optimized structure of the metal bipolar plate for a fuel cell provided in an embodiment of the present invention;

[0034] Figure 3b A schematic diagram of the structure of a fuel cell metal bipolar plate before optimization provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0035] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0036] It should be understood that when used in this specification and the appended claims, the terms "include" and "comprises" indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

[0037] It should also be understood that the terms used in this specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the specification of the present invention and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an" and "the" are intended to include plural forms.

[0038] It should be further understood that the term "and / or" used in the present description and the appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0039] As used in this specification and the appended claims, the term "if" may be interpreted as "when" or "upon" or "in response to determining" or "in response to detecting," depending on the context. Similarly, the phrases "if it is determined" or "if [described condition or event] is detected" may be interpreted as meaning "upon determination" or "in response to determining" or "upon detection of [described condition or event]" or "in response to detecting [described condition or event]," depending on the context.

[0040] See also Figure 1 , Figure 1 Schematic diagram of the global optimization design and comprehensive evaluation method of the fuel cell metal bipolar plate provided by the embodiment of the present invention. Figure 1 As shown, the method comprises the following steps:

[0041] S10, obtaining operating parameters of the metal bipolar plate of the fuel cell.

[0042] In a specific implementation, the operating parameters of the metal bipolar plate of the fuel cell can be input by the user.

[0043] Specifically, the operating parameters of the metal bipolar plate include operating current density, hydrogen excess coefficient, reaction active area, number of cells, oxygen excess coefficient, single cell active area and cell voltage. The above parameters can be input by the user.

[0044] S20, determining the hydrogen molar flow rate, air molar flow rate and heat flow rate of the metal bipolar plate according to the operating parameters.

[0045] In a specific implementation, after receiving the operating parameters of the metal bipolar plate of the fuel cell, the hydrogen molar flow rate, air molar flow rate and heat flow rate of the metal bipolar plate are determined according to the operating parameters. The specific calculation method is as follows:

[0046] First, for the hydrogen molar flow rate, the hydrogen molar flow rate is expressed by the following formula: Calculate, where m1 is the hydrogen molar flow rate, i is the operating current density, N0 is the number of cells, Aact is the reaction active area, S1 is the hydrogen excess coefficient and F is the Faraday constant.

[0047] Further, for the air molar flow rate, the air molar flow rate is expressed by the following formula: Calculate, where S2 is the oxygen excess coefficient.

[0048] Further, for the heat flow, the heat flow is calculated by the following formula Q=m1ΔH-iAV, wherein ΔH=286KJ / mol, A is the active area of ​​the single cell, and V is the battery voltage.

[0049] S30, determining initial geometric parameters of the metal bipolar plate according to the hydrogen molar flow rate, the air molar flow rate, and the heat flow rate.

[0050] In a specific implementation, the initial geometric parameters of the metal bipolar plate include the initial design parameters of the hydrogen flow channel, hydrogen manifold, air flow channel, air manifold, distribution area and cooling flow channel of the metal bipolar plate. The initial design parameters may be, for example, dimensional information such as length, width and height.

[0051] Specifically, initial design parameters of the air flow channel and the air manifold are determined according to the air molar flow rate.

[0052] The optional parameter range of the air flow channel is determined according to the air molar flow rate, and then a value is selected from the optional parameter range as the initial design parameter of the air flow channel. The optional parameter range corresponding to different air molar flow rates can be pre-set by those skilled in the art, and the present invention is not specifically limited. For example, the optional parameter range can be preferably 0.4 mm-1 mm.

[0053] The optional parameter range of the air manifold is determined according to the air molar flow rate, and then a value is selected from the optional parameter range as the initial design parameter of the air manifold. The optional parameter range corresponding to different air molar flow rates can be pre-set by those skilled in the art, and the present invention is not specifically limited.

[0054] Furthermore, initial design parameters of the hydrogen flow channel and the hydrogen manifold are determined according to the hydrogen molar flow rate.

[0055] The optional parameter range of the hydrogen flow channel is determined according to the hydrogen molar flow rate, and then a value is selected from the optional parameter range as the initial design parameter of the hydrogen flow channel. The optional parameter range corresponding to different hydrogen molar flow rates can be pre-set by those skilled in the art, and the present invention is not specifically limited. For example, the optional parameter range can be preferably 0.4 mm-1 mm.

[0056] The optional parameter range of the hydrogen manifold is determined according to the hydrogen molar flow rate, and then a value is selected from the optional parameter range as the initial design parameter of the hydrogen manifold. The optional parameter range corresponding to different hydrogen molar flow rates can be pre-set by those skilled in the art, and the present invention is not specifically limited.

[0057] Furthermore, initial design parameters of the distribution area are determined according to the hydrogen molar flow rate and the air molar flow rate.

[0058] Specifically, firstly, the sum of the hydrogen molar flow rate and the air molar flow rate is calculated to obtain the total flow rate, and the initial design parameters of the distribution area are determined according to the total flow rate.

[0059] The optional parameter range of the distribution area is determined according to the total flow, and then a value is selected from the optional parameter range as the initial design parameter of the distribution area. The optional parameter range corresponding to different total flow rates can be pre-set by those skilled in the art, and the present invention is not specifically limited.

[0060] S40, determining an initialization bipolar plate configuration of the metal bipolar plate according to initial geometric parameters of the metal bipolar plate.

[0061] In a specific implementation, after the initial geometric parameters of the metal bipolar plate are determined, the initialization bipolar plate configuration of the metal bipolar plate can be constructed according to the initial geometric parameters.

[0062] S50, establishing a three-dimensional model according to the initialized bipolar plate configuration, performing simulation, and determining an evaluation index based on the simulation result.

[0063] In a specific implementation, the evaluation index can be set by those skilled in the art, and the present invention does not specifically limit it. For example, the evaluation index can be referred to in Table 1 below.

[0064] Table 1. Evaluation index table

[0065]

[0066] In one embodiment of the present invention, the evaluation indicators are selected as effective mass transfer coefficient, coefficient of variation, synergy angle and net mass transfer power.

[0067] Specifically, a three-dimensional model is established according to the initialized bipolar plate configuration, and simulation is performed to obtain simulation results based on electrochemical characteristics, flow characteristics, and heat transfer characteristics. A secondary analysis of the effective mass transfer coefficient, coefficient of variation, synergy angle, and net mass transfer power is performed on the simulation results to obtain the effective mass transfer coefficient, coefficient of variation, synergy angle, and net mass transfer power.

[0068] Specifically, the modeling process includes two parts: solid meshing is performed according to the initialization bipolar plate configuration to establish a finite element model of the bipolar plate flow channel structure of the hydrogen fuel cell; fluid meshing is performed according to the initialization bipolar plate configuration to establish a three-dimensional two-phase flow model of the bipolar plate flow channel structure of the hydrogen fuel cell. The finite element model and the three-dimensional two-phase flow model can simulate the electrochemical reaction, fluid flow and heat transfer characteristics of the fuel cell to obtain simulation results, which include effective mass transfer coefficient, coefficient of variation, synergy angle, and net mass transfer power. The coefficient of variation is a dimensionless constant that measures the uniformity of the inlet flow rate.

[0069] The effective mass transfer coefficient is given by the following formula Calculation, where EMTC is the effective mass transfer coefficient, v and are the components of velocity and concentration gradient in the normal direction of the proton exchange membrane, respectively.

[0070] The coefficient of variation is given by the following formula Calculation; where N is the number of flow channels, q k is the inlet mass flow rate of the kth flow channel, is the average mass flow rate of the flow channel.

[0071] The synergy angle is based on the following formula Calculate, where c is the concentration.

[0072] The net mass transfer power is given by the following formula Calculate, where V is the battery voltage, i is the current, ΔP is the loss along the flow channel from the inlet to the outlet, and A is the active area of ​​the single cell.

[0073] S60, optimizing and designing the geometric parameters of the metal bipolar plate according to the evaluation index to obtain the final geometric parameters of the metal bipolar plate.

[0074] In a specific implementation, the evaluation index is used as the optimization target, the geometric parameters of the metal bipolar plate are used as the decision variables, and the configuration of the metal bipolar plate is optimized by a multi-objective optimization algorithm. When the evaluation index is optimal, the optimal geometric parameter solution set is obtained, and the optimal geometric parameter solution set is used as the final geometric parameters of the metal bipolar plate. The multi-objective optimization method includes but is not limited to:

[0075] Vector evaluation genetic algorithm, NSGA-II, NAGA-III, IBEA, SMS-EMOA, MOGA, SEEA, SPEA2, MOLS.

[0076] In the embodiment of the present invention, the structural parameters (geometric parameters) of the metal bipolar plate to be optimized can be set by those skilled in the art, and the present invention does not specifically limit them. For example, the structural parameters of the metal bipolar plate to be optimized include but are not limited to the parameters recorded in Table 2.

[0077] Table 2. Structural parameters of metal bipolar plates to be optimized

[0078]

[0079] In one embodiment, the evaluation indicators are selected as effective mass transfer coefficient, coefficient of variation, synergy angle and net mass transfer power, and the geometric parameters of the metal bipolar plate are optimized and designed according to the effective mass transfer coefficient, coefficient of variation, synergy angle and net mass transfer power to obtain the final geometric parameters of the metal bipolar plate.

[0080] Specifically, see Figure 2 In one embodiment, the effective mass transfer coefficient, the coefficient of variation, the synergy angle and the net mass transfer power are used as optimization targets, and the geometric parameters of the metal bipolar plate are used as decision variables. The NSGA-III algorithm is used to perform multi-objective optimization on the configuration of the metal bipolar plate. When the effective mass transfer coefficient, the coefficient of variation, the synergy angle and the net mass transfer power are optimized, an optimal geometric parameter solution set is obtained, and the optimal geometric parameter solution set is used as the final geometric parameters of the metal bipolar plate.

[0081] The final geometric parameters of the metal bipolar plate include the final design parameters of the hydrogen flow channel, hydrogen manifold, air flow channel, air manifold, distribution area and cooling flow channel of the metal bipolar plate. The final design parameters may be, for example, the length, width and height of the flow channel, the cross-sectional geometry and other dimensional information, i.e., the geometric parameters.

[0082] Among them, Figure 2 In the embodiment, the flow channel geometric parameters include the geometric parameters of the hydrogen flow channel, the air flow channel and the cooling flow channel. The fuel cell bipolar plate is the metal bipolar plate.

[0083] S70, determining a final bipolar plate configuration of the metal bipolar plate based on the final geometric parameters.

[0084] In a specific implementation, after the final geometric parameters of the metal bipolar plate are determined, the final bipolar plate configuration of the metal bipolar plate can be constructed according to the final geometric parameters.

[0085] See also Figure 3a-Figure 3b , Figure 3a The optimized flow channel design is shown. Figure 3b The flow channel design before optimization is shown. Experiments were conducted based on the above flow channel design, and the experimental results are shown in Table 3 below.

[0086] Table 3. Test results

[0087]

[0088]

[0089] Based on the above results, it can be seen that the EMTC, diffusion coefficient and mass transfer coefficient of the optimized wave flow channel are increased, indicating that the design optimization of the flow channel enhances the mass transfer effect.

[0090] The technical solution of the embodiment of the present invention obtains the working parameters of the metal bipolar plate of the fuel cell; determines the hydrogen molar flow rate, air molar flow rate and heat flow rate of the metal bipolar plate according to the working parameters; determines the initial geometric parameters of the metal bipolar plate according to the hydrogen molar flow rate, the air molar flow rate and the heat flow rate; determines the initialization bipolar plate configuration of the metal bipolar plate according to the initial geometric parameters of the metal bipolar plate; establishes a three-dimensional model according to the initialization bipolar plate configuration, performs simulation, obtains simulation results, and determines evaluation indicators based on the simulation results; optimizes the geometric parameters of the metal bipolar plate according to the evaluation indicators to obtain the final geometric parameters of the metal bipolar plate; determines the final bipolar plate configuration of the metal bipolar plate based on the final geometric parameters, so as to obtain the optimally designed fuel cell metal bipolar plate, ensuring that the best performance can be achieved.

[0091] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0092] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

[0093] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed by the present invention, and these modifications or replacements should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be based on the protection scope of the claims.

Claims

1. A global optimization design and comprehensive evaluation method for a fuel cell metal bipolar plate, characterized in that: include: Obtaining operating parameters of a metal bipolar plate of a fuel cell; Determining the hydrogen molar flow rate, air molar flow rate and heat flow rate of the metal bipolar plate according to the operating parameters; Determining initial geometric parameters of the metal bipolar plate according to the hydrogen molar flow rate, the air molar flow rate and the heat flow rate; Determining an initialization bipolar plate configuration of the metal bipolar plate according to initial geometric parameters of the metal bipolar plate; Establishing a three-dimensional model according to the initialized bipolar plate configuration, performing simulation, obtaining simulation results, and determining evaluation indicators based on the simulation results; Optimizing the geometric parameters of the metal bipolar plate according to the evaluation index to obtain the final geometric parameters of the metal bipolar plate; determining a final bipolar plate configuration of the metal bipolar plate based on the final geometric parameters; The step of optimizing the geometric parameters of the metal bipolar plate according to the evaluation index to obtain the final geometric parameters of the metal bipolar plate includes: Taking the evaluation index as the optimization target and the geometric parameters of the metal bipolar plate as the decision variables, a multi-objective optimization algorithm is used to perform multi-objective optimization on the configuration of the metal bipolar plate, and when the evaluation index is optimal, an optimal geometric parameter solution set is obtained, and the optimal geometric parameter solution set is used as the final geometric parameters of the metal bipolar plate; The multi-objective optimization method includes but is not limited to: Vector Evaluation Genetic Algorithm, NSGA-II, NAGA-III, IBEA, SMS-EMOA, MOGA, SEEA, SPEA2, MOLS; The optimization objectives include but are not limited to effective mass transfer coefficient, coefficient of variation, synergy angle and net mass transfer power; The effective mass transfer coefficient is given by the following formula Calculation, where EMTC is the effective mass transfer coefficient, v and are the components of velocity and concentration gradient in the normal direction of the proton exchange membrane, respectively; The coefficient of variation is given by the following formula Calculation; where N is the number of flow channels, q k is the inlet mass flow rate of the kth flow channel, is the average mass flow rate of the flow channel; The synergy angle is based on the following formula Calculate, where c is the concentration; The net mass transfer power is given by the following formula Calculate, where V is the battery voltage, i is the current, ΔP is the loss along the flow channel from the inlet to the outlet, and A is the active area of ​​the single cell; The structural parameters of the metal bipolar plate to be optimized include flow channel size, ridge size, flow channel cross-sectional geometry, flow channel number, distribution area geometry, distribution structure, distribution area size, inlet / outlet cross-sectional shape and inlet / outlet size.

2. The global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates according to claim 1, characterized in that: The working parameters include working current density, hydrogen excess coefficient, reaction active area and number of single cells. The hydrogen molar flow rate is given by the following formula Calculate, where m1 is the hydrogen molar flow rate, i is the operating current density, N0 is the number of cells, Aact is the reaction active area, S1 is the hydrogen excess coefficient and F is the Faraday constant.

3. The global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates according to claim 2, characterized in that: The working parameters also include the oxygen excess coefficient, and the air molar flow rate is expressed by the following formula Calculate, where S2 is the oxygen excess coefficient.

4. The global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates according to claim 3, characterized in that: The working parameters also include the active area of ​​the single cell and the battery voltage. The heat flow is calculated by the following formula Q=m1ΔH-iAV, where ΔH=286KJ / mol, A is the active area of ​​the single cell, and V is the battery voltage.

5. The global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates according to claim 1, characterized in that: The initial geometric parameters of the metal bipolar plate include initial design parameters of a hydrogen flow channel and a hydrogen manifold of the metal bipolar plate, wherein the initial design parameters of the hydrogen flow channel and the hydrogen manifold are determined according to the hydrogen molar flow rate.

6. The global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates according to claim 1, characterized in that: The initial geometric parameters of the metal bipolar plate also include initial design parameters of the air flow channel and the air manifold of the metal bipolar plate, wherein the initial design parameters of the air flow channel and the air manifold are determined according to the air molar flow rate.

7. The global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates according to claim 1, characterized in that: The initial geometric parameters of the metal bipolar plate also include initial design parameters of a distribution area of ​​the metal bipolar plate, wherein the initial design parameters of the distribution area are determined according to the hydrogen molar flow rate and the air molar flow rate.

8. The global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates according to claim 1, characterized in that: The initial geometric parameters of the metal bipolar plate also include initial design parameters of the cooling channel of the metal bipolar plate, wherein the initial design parameters of the cooling channel are determined according to the heat flow rate.

9. The global optimization design and comprehensive evaluation method for fuel cell metal bipolar plates according to claim 1, characterized in that: The establishing of a three-dimensional model according to the initialized bipolar plate configuration comprises: Solid meshing is performed according to the initialized bipolar plate configuration to establish a finite element model of the bipolar plate flow channel structure of the hydrogen fuel cell; A three-dimensional two-phase flow model of the bipolar plate flow channel structure of a hydrogen fuel cell is established by performing fluid grid division according to the initialized bipolar plate configuration.

Images