A method and computer program product for optimizing a flow channel cross-section of a hydrogen fuel cell bipolar plate
By establishing a multi-parameter optimization model and performing nested loop traversal search, the cross-section of the bipolar plate flow channel in hydrogen fuel cells is optimized, solving the problem that the flow channel design is difficult to balance area and proportion, improving flow capacity and stability, and making it suitable for mass production of titanium bipolar plates.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- UESTC (SHENZHEN) ADVANCED RES INST
- Filing Date
- 2026-05-08
- Publication Date
- 2026-07-07
Smart Images
Figure CN122133359B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of hydrogen fuel cell technology, and in particular to a method for optimizing the cross-sectional design of bipolar plate flow channels in hydrogen fuel cells and a computer program product. Background Technology
[0002] As a clean energy conversion device, the core component of a hydrogen fuel cell, the bipolar plate, performs multiple functions, including distributing reactant gases, collecting current, discharging generated water, and dissipating heat. The flow channel structure, used to uniformly distribute the reactant fluid and discharge the products, directly affects the gas-liquid two-phase flow characteristics, pressure drop distribution, and electrochemical reaction uniformity within the cell.
[0003] Traditional bipolar plates often employ rectangular or simple trapezoidal flow channels, which suffer from problems such as small flow area, high pressure loss, and high risk of flooding, limiting further improvements in battery performance. In recent years, methods such as topology optimization and biomimetic design have emerged to improve the flow channel configuration, but these often overlook manufacturing process constraints (such as the tensile limit of metal sheets) and actual assembly conditions (such as fixed plate width). Especially for titanium bipolar plates, although they possess excellent corrosion resistance and strength, their molding process is limited by the elongation rate, requiring a reasonable layout of the number and geometry of flow channels within a limited width.
[0004] In the process of realizing this invention, the inventors discovered at least the following problems in the prior art:
[0005] Existing bipolar plate flow channel cross-section designs struggle to balance maximizing the flow channel cross-sectional area with constraints such as aspect ratio and flow channel / rib width ratio, making it difficult to ensure the output power density and operational stability of hydrogen fuel cells. Summary of the Invention
[0006] The purpose of this invention is to provide a method and computer program product for optimizing the flow channel cross-section of a bipolar plate in a hydrogen fuel cell. This addresses the technical problem in existing bipolar plate flow channel cross-section designs where maximizing the flow channel cross-sectional area is difficult to balance with comprehensive constraints such as aspect ratio and flow channel / rib width ratio, thus hindering the assurance of output power density and operational stability of the hydrogen fuel cell. The numerous technical effects of the preferred solutions among the many technical solutions provided by this invention are detailed below.
[0007] To achieve the above objectives, the present invention provides the following technical solution:
[0008] This invention provides a method for optimizing the flow channel cross-section of a hydrogen fuel cell bipolar plate, comprising the following steps: S100: Based on the forming process boundary of the pure titanium bipolar plate, a multi-parameter optimization model of the flow channel area S is established, including the flow channel width a, rib width b, forming angle θ, hypotenuse length d, and the number of flow channels n, and the value range and step size of each parameter are limited; S200: In the multi-parameter optimization model of the flow channel area S, a first constraint is established based on the ratio of the flow channel width a to the rib width b, and a constraint is established based on the aspect ratio h of the flow channel of the pure titanium bipolar plate. The second constraint is to establish the third constraint on the total width W after molding, based on the flow channel width a, rib width b, hypotenuse length d, and the number of flow channels n; S300: Under the first, second, and third constraints, based on the value range and step size of each parameter in the multi-parameter optimization model of flow channel area S, the maximum value of flow channel area S and the corresponding values of flow channel width a, rib width b, molding angle θ, hypotenuse length d, and number of flow channels n are calculated by traversing and searching all parameter combinations in the multi-parameter optimization model of flow channel area S through a five-fold nested loop.
[0009] Preferably, in step S100, the value of the flow channel width a ranges from 0.2 to 2, the value of the rib width b ranges from 0.7 to 1, the value of the forming angle ranges from π / 3 to π / 2, the value of the inclined side length d ranges from 0 to 2, and the value of the number of flow channels n ranges from 20 to 50.
[0010] Preferably, in step S100, the step size of the flow channel width a is 0.05, the step size of the rib width b is 0.05, the step size of the forming angle is π / 180, the step size of the inclined side length d is 0.1, and the step size of the number of flow channels n is 1.
[0011] Preferably, in step S200, the ratio of the width a of the first constraint channel to the width b of the rib is in the range of 0.25≤a / b≤2.5, so that the structural strength and fluid performance of the pure titanium bipolar plate are balanced.
[0012] Preferably, in step S200, the formula for calculating the aspect ratio h of the second constraint channel is h=2. d sinθ / (a+b+4 d The value of cosθ is 0.4 < h < 1.1 to avoid the pure titanium bipolar plate being too flat or too tall.
[0013] Preferably, in step S200, the forming process boundary of the pure titanium bipolar plate is an original width of 50cm, and a maximum width of 67.5cm after stretching based on a stretching ratio of 1.35; the expression for the total width W after forming is W=n (a+b+2d), the third constraint of the total width W is 50≤W≤67.5.
[0014] Preferably, in step S300, the five-fold nested loop traversal search process specifically includes: S310: traversing each channel width a value based on the range and step size of the channel width a; S320: for each channel width a value, traversing each rib width b value based on the range and step size of the rib width b value, and determining whether the ratio of the channel width a value to the rib width b value satisfies the first constraint; S330: for channel width a values and rib width b values that satisfy the first constraint, traversing each forming angle θ value based on the range and step size of the forming angle θ; S340: for each forming angle θ value, traversing each hypotenuse length d value based on the range and step size of the hypotenuse length d; 350: Calculate the runner aspect ratio h and determine whether the runner aspect ratio h satisfies the second constraint; S360: For the runner width a, rib width b, forming angle θ, and hypotenuse length d that satisfy the second constraint, traverse each runner number n value based on the range of runner number n and the step size, and calculate the total width W; S370: If the total width W satisfies the third constraint, calculate the runner area S based on all combinations that satisfy the third constraint. If the new runner area S is greater than the currently recorded maximum value, update the runner area S until all traversals are completed to obtain the optimal solution for the runner area S, and record the corresponding values of runner width a, rib width b, forming angle θ, hypotenuse length d, and runner number n.
[0015] Preferably, in step S300, the formula for calculating the area S of each flow channel is: a d sinθ+d 2 sinθ cosθ, where a d sinθ represents the area of the rectangular portion of each flow channel in the cross-section of the pure titanium bipolar plate, d 2 sinθ cosθ represents the sum of the areas of the two triangular flanks of each flow channel in the cross-section of a pure titanium bipolar plate.
[0016] Preferably, in step S300, the maximum value of the flow channel area S is 1.0132, wherein the flow channel width a is 0.85, the rib width b is 0.70, the forming angle θ is 60 degrees, the inclined side length d is 0.90, and the number of flow channels n is 20.
[0017] A computer program product includes a computer program / instructions that, when executed by a processor, implement the steps of the method for optimizing the flow channel cross-section of a hydrogen fuel cell bipolar plate as described above.
[0018] Implementing one of the above-described technical solutions of the present invention has the following advantages or beneficial effects:
[0019] This invention establishes a multi-parameter optimization model based on the titanium plate forming process boundary, including channel width *a*, rib width *b*, forming angle *θ*, hypotenuse length *d*, and the number of channels *n*. It introduces three constraints: the ratio of channel width *a* to rib width *b*, the channel height-to-width ratio *h*, and the total width *W* after forming. Through traversal search, the feasible solution with the largest cross-sectional area is obtained. While satisfying the titanium forming process constraints and structural stability requirements, the cross-sectional area of a single channel is significantly increased. The larger channel cross-section effectively enhances the flow capacity of reactant gases and cooling media, reduces pressure drop, and improves the internal mass transfer efficiency and hydrothermal management performance of the fuel cell, thereby increasing the output power density and operational stability of the fuel cell. Simultaneously, the optimized channel height-to-width ratio and reasonable rib width support structure balance conductivity, mechanical strength, and processing feasibility, avoiding the flooding or dry film problems easily encountered in traditional designs. This makes it suitable for mass production of metal bipolar plates, demonstrating good engineering application prospects and industrialization value. Attached Figure Description
[0020] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. In the drawings:
[0021] Figure 1 This is a flowchart of an embodiment of the present invention: a method for optimizing the cross-sectional design of a bipolar plate flow channel in a hydrogen fuel cell.
[0022] Figure 2 This is a three-dimensional pure titanium bipolar plate in the method for optimizing the flow channel cross-section of a hydrogen fuel cell bipolar plate according to Embodiment 1 of the present invention. Figure 1 ;
[0023] Figure 3 This is a three-dimensional pure titanium bipolar plate in the method for optimizing the flow channel cross-section of a hydrogen fuel cell bipolar plate according to Embodiment 1 of the present invention. Figure 2 ;
[0024] Figure 4 This is a schematic diagram of the cross-section of a pure titanium bipolar plate in an embodiment of the present invention, namely, a method for optimizing the flow channel cross-section of a hydrogen fuel cell bipolar plate.
[0025] Figure 5 This is a flowchart of the five-fold nested loop traversal search process in the hydrogen fuel cell bipolar plate flow channel cross-section optimization design method of Embodiment 1 of the present invention;
[0026] Figure 6 This is a schematic diagram of all feasible solutions that satisfy the constraints in the optimization design method for the bipolar plate flow channel cross section of a hydrogen fuel cell according to Embodiment 1 of the present invention.
[0027] Figure 7 This is a schematic diagram of the optimal solution set of the Pareto front in a method for optimizing the cross-section of a bipolar plate flow channel in a hydrogen fuel cell, according to Embodiment 1 of the present invention. Detailed Implementation
[0028] To make the objectives, technical solutions, and advantages of the present invention clearer, various exemplary embodiments described below will be referenced to the accompanying drawings, which form part of the exemplary embodiments, illustrating various exemplary embodiments that may be used to implement the present invention. Unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this disclosure. It should be understood that they are merely examples of processes, methods, and apparatuses consistent with some aspects of the present invention disclosed as detailed in the appended claims, and other embodiments may be used, or structural and functional modifications may be made to the embodiments listed herein without departing from the scope and spirit of the present invention.
[0029] In the description of this invention, it should be understood that the terms "center," "longitudinal," "lateral," etc., indicate the orientation or positional relationship based on the accompanying drawings, and are only for the convenience of describing the invention and simplifying the description, and do not indicate or imply that the referred element must have a specific orientation, or be constructed and operated in a specific orientation. The terms "first," "second," etc., are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. The term "multiple" means two or more. The terms "connected" and "linked" should be interpreted broadly, for example, they can be fixed connections, detachable connections, integral connections, mechanical connections, electrical connections, communication connections, direct connections, indirect connections through an intermediate medium, and can be the internal connection of two elements or the interaction relationship between two elements. The term "and / or" includes any and all combinations of one or more of the related listed items. Those skilled in the art can understand the specific meaning of the above terms in this invention according to the specific circumstances.
[0030] To illustrate the technical solution described in this invention, specific embodiments are described below, showing only the parts related to the embodiments of this invention.
[0031] Example 1:
[0032] like Figures 1-4As shown, this invention provides a method for optimizing the flow channel cross-section of a hydrogen fuel cell bipolar plate, including the following steps: S100: Based on the forming process boundary of the pure titanium bipolar plate, a multi-parameter optimization model of the flow channel area S is established, which includes the flow channel width a, rib width b, forming angle θ, hypotenuse length d, and the number of flow channels n. That is, a multi-parameter optimization model with the flow channel area S as the objective function is established, and the value range and step size of each parameter are limited. Thus, the flow channel width a, rib width b, forming angle θ, hypotenuse length d, and the number of flow channels n are discretized values, which can form multiple combinations of the values of each parameter, making it convenient to calculate the flow channel area S by iterative search in the future. S200: In the multi-parameter optimization model of the flow channel area S, the first constraint is set by the ratio of the flow channel width a to the rib width b, the second constraint is set by the flow channel height-to-width ratio h of the pure titanium bipolar plate, and the third constraint is set by the flow channel width a, rib width b, hypotenuse length d and the number of flow channels n. Under the first, second and third constraints, the rationality and manufacturability of the pure titanium bipolar plate structure are ensured. S300: Under the first, second, and third constraints, based on the value range and step size of each parameter in the multi-parameter optimization model of the flow area S, a five-fold nested loop traversal is used to search for all parameter combinations in the multi-parameter optimization model of the flow area S. Nested loop traversal is a powerful tool for handling multi-dimensional data structures and complex traversal scenarios. One loop contains another loop. Each time the outer loop is executed, the inner loop will complete one round. Thus, all combinations of flow width a, rib width b, forming angle θ, hypotenuse length d, and flow number n that satisfy the constraints are calculated, ensuring the reliability of the calculation. The maximum value of the flow area S and the corresponding values of flow width a, rib width b, forming angle θ, hypotenuse length d, and flow number n are obtained. This embodiment establishes a multi-parameter optimization model based on the titanium plate forming process boundary, including channel width *a*, rib width *b*, forming angle *θ*, hypotenuse length *d*, and the number of channels *n*. It introduces three constraints: the ratio of channel width *a* to rib width *b*, the channel height-to-width ratio *h*, and the total width *W* after forming. Through traversal search, the feasible solution with the largest cross-sectional area is obtained. While satisfying the titanium forming process constraints and structural stability requirements, the cross-sectional area of a single channel is significantly increased. The larger channel cross-section effectively enhances the flow capacity of reactant gases and cooling media, reduces pressure drop, and improves the internal mass transfer efficiency and hydrothermal management performance of the fuel cell, thereby improving the output power density and operational stability of the fuel cell. Simultaneously, the optimized channel height-to-width ratio and reasonable rib width support structure balance conductivity, mechanical strength, and processing feasibility, avoiding the flooding or dry film problems easily encountered in traditional designs. This makes it suitable for mass production of metal bipolar plates, demonstrating good engineering application prospects and industrialization value.
[0033] As an optional implementation, in step S100, the value of the flow channel width a ranges from 0.2 to 2, the value of the rib width b ranges from 0.7 to 1, the value of the forming angle θ ranges from π / 3 to π / 2, the value of the hypotenuse length d ranges from 0 to 2, and the value of the number of flow channels n ranges from 20 to 50. The step size of the flow channel width 'a' is 0.05, resulting in values of 0.2, 0.25, 0.3, ..., 2; the step size of the rib width 'b' is 0.05, resulting in values of 0.7, 0.75, 0.8, ..., 1.0; the step size of the forming angle 'θ' is π / 180, or 1°, resulting in values of 60°, 61°, 62°, ..., 90°; the step size of the hypotenuse length 'd' is 0.1, resulting in values of 0, 0.1, 0.2, ..., 2.0; and the step size of the number of flow channels 'n' is 1, resulting in values of 20, 21, 22, ..., 50.
[0034] As an optional implementation, in step S200, the ratio of the first constraint channel width a to the rib width b is in the range of 0.25≤a / b≤2.5, so as to balance the structural strength and fluid performance of the pure titanium bipolar plate and prevent the channel from being too narrow or too wide relative to the rib, so as to ensure the structural rationality and manufacturability of the pure titanium bipolar plate.
[0035] As an optional implementation, in step S200, such as Figure 4 As shown, the formula for calculating the aspect ratio h of the second constrained flow channel is h=2. d sinθ / (a+b+4 d The value of cosθ is 0.4 < h < 1.1. This ratio reflects the proportion of the flow channel “height” to the “periodic width” (one flow channel + one rib + projection of the two oblique sides), avoiding excessive flatness or height which would lead to manufacturing difficulties or excessive flow resistance.
[0036] As an optional implementation, in step S200, the forming process boundary of the pure titanium bipolar plate is an original width of 50 cm and a maximum width of 67.5 cm after stretching based on a stretching ratio of 1.35. The original width and stretching ratio of the pure titanium bipolar plate can be adjusted accordingly according to actual needs; the expression for the total width W after forming is W=n (a+b+2d), the third constraint of the total width W is 50≤W≤67.5, and the goal is to design a trapezoidal flow channel structure with the largest single flow channel cross-sectional area under this width constraint.
[0037] As an optional implementation, in step S300, such as Figure 5As shown, the five-fold nested loop traversal search process specifically includes: S310: Traversing each channel width a value based on the range and step size of the channel width a; S320: For each channel width a value, traversing each rib width b value based on the range and step size of the rib width b value, and determining whether the ratio of the channel width a value to the rib width b value satisfies the first constraint. If not, skipping and directly proceeding to the next rib width b value traversal; S330: For channel width a values and rib width b values that satisfy the first constraint, traversing each forming angle θ value based on the range and step size of the forming angle θ; S340: For each forming angle θ value, traversing each hypotenuse length d value based on the range and step size of the hypotenuse length d; S350: Calculating the channel height-to-width ratio h, and determining whether the channel height-to-width ratio h satisfies the second constraint. If not, ... Skip the current hypotenuse length d value and proceed to the next hypotenuse length d value; S360: For the channel width a, rib width b, forming angle θ, and hypotenuse length d that satisfy the second constraint, traverse each channel number n value based on the range of channel number n and the step size, and calculate the total width W. If the total width W does not satisfy the third constraint, skip it and proceed to the next channel number n value; S370: If the total width W satisfies the third constraint, calculate the channel area S based on all combinations that satisfy the third constraint. If the new channel area S is greater than the currently recorded maximum value, update the channel area S until all traversals are completed to obtain the optimal solution for the channel area S, and record the corresponding values of channel width a, rib width b, forming angle θ, hypotenuse length d, and channel number n, thus realizing the trapezoidal channel structure design for the maximum single channel cross-sectional area of the bipolar plate. The entire nested loop traversal search process can be implemented using tools such as Matlab.
[0038] As an optional implementation, in step S300, the formula for calculating the area S of each flow channel is: a d sinθ+d 2 sinθ cosθ, such as Figure 4 As shown, where a d sinθ represents the area of the rectangular portion of each flow channel in the cross-section of the pure titanium bipolar plate, where the base of the rectangle is 'a' and the height is 'dsinθ', thus the area is 'a'. d sinθ;d2 sinθ cosθ represents the sum of the areas of the two triangles on the side flanks of each flow channel in the cross-section of a pure titanium bipolar plate, such as... Figure 4 As shown, the base of the triangle is dcosθ and the height is dsinθ, therefore the area of each triangle is (d... 2 sinθ cosθ) / 2, and thus finally the formula for calculating the area S of each flow channel is obtained. d sinθ+d 2 sinθ cosθ.
[0039] As an optional implementation, in step S300, based on Figure 6 , Figure 7 Simulation results show that the optimal aspect ratio h of the second constraint channel is 0.465, and the maximum channel area S is 1.013. Based on the traversal process, the specific values of each parameter corresponding to the channel area S being 1.013 can be obtained. Specifically, the maximum channel area S is 1.013, where the channel width a is 0.85, the rib width b is 0.70, the forming angle θ is 60 degrees (i.e., 1.0472 radians), the hypotenuse length d is 0.90, and the number of channels n is 20. At this time, the total width W after forming is 20. (0.85+0.70+2) 0.90) = 67.0cm, which meets the third constraint of the total width W: 50≤W≤67.5.
[0040] The embodiment is merely a specific example and does not indicate that this is the only way to implement the present invention.
[0041] Example 2:
[0042] A computer program product includes a computer program / instructions that, when executed by a processor, implement the steps of the hydrogen fuel cell bipolar plate flow channel cross-section optimization design method described in Embodiment 1. This embodiment establishes a multi-parameter optimization model based on the titanium plate forming process boundary, including flow channel width a, rib width b, forming angle θ, hypotenuse length d, and the number of flow channels n. It introduces three constraints: the ratio of flow channel width a to rib width b, the flow channel height-to-width ratio h, and the total width W after forming. Through traversal search, a feasible solution with the largest cross-sectional area is obtained. Under the premise of satisfying the titanium material forming process constraints and structural stability requirements, the cross-sectional area of a single flow channel is significantly improved.
[0043] The above description is merely a preferred embodiment of the present invention. Those skilled in the art will understand that various changes or equivalent substitutions can be made to these features and embodiments without departing from the spirit and scope of the present invention. Furthermore, under the teachings of the present invention, these features and embodiments can be modified to adapt to specific situations and materials without departing from the spirit and scope of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed herein, and all embodiments falling within the scope of the claims of this application are within the protection scope of the present invention.
Claims
1. A method for optimizing the cross-sectional design of the bipolar plate flow channel in a hydrogen fuel cell, characterized in that, Includes the following steps: S100: Based on the molding process boundary of pure titanium bipolar plates, a multi-parameter optimization model of the flow area S, including the flow channel width a, rib width b, molding angle θ, hypotenuse length d and the number of flow channels n, is established, and the value range and step size of each parameter are limited. S200: In the multi-parameter optimization model of the flow channel area S, the first constraint is set by the ratio of the flow channel width a to the rib width b, the second constraint is set by the flow channel height-to-width ratio h of the pure titanium bipolar plate, and the third constraint is set by the flow channel width a, rib width b, hypotenuse length d and the number of flow channels n. S300: Under the first, second, and third constraints, based on the value range and step size of each parameter in the multi-parameter optimization model of the flow channel area S, the maximum value of the flow channel area S and the corresponding values of the flow channel width a, rib width b, forming angle θ, hypotenuse length d, and number of flow channels n are calculated by traversing and searching all parameter combinations in the multi-parameter optimization model of the flow channel area S through a five-fold nested loop. In step S200, the formula for calculating the aspect ratio h of the second constraint channel is h=2. d sinθ / (a+b+4 d The value of cosθ is 0.4 < h < 1.1 to avoid the pure titanium bipolar plate being too flat or too tall; the forming process boundary of the pure titanium bipolar plate is the original width of 50cm, and the maximum width after stretching is 67.5cm based on the stretching ratio of 1.
35. The expression for the total width W after molding is W=n (a+b+2d), the third constraint of the total width W is 50≤W≤67.5; In step S300, the five-fold nested loop traversal search process specifically includes: S310: Traversing each channel width a value based on the range of channel width a and the step size; S320: For each runner width 'a', iterate through each rib width 'b' value based on the range and step size of the rib width 'b', and determine whether the ratio of runner width 'a' to rib width 'b' satisfies the first constraint; S330: For runner width 'a' and rib width 'b' values that satisfy the first constraint, iterate through each forming angle 'θ' value based on the range and step size of the forming angle 'θ'; S340: For each forming angle 'θ' value, iterate through each hypotenuse length 'd' value based on the range and step size of the hypotenuse length 'd'; S350: Calculate the runner height-to-width ratio 'h', and determine whether the runner height-to-width ratio 'h' satisfies the second constraint. S360: For the runner width a, rib width b, forming angle θ, and hypotenuse length d that satisfy the second constraint, traverse each runner number n value based on the range of runner number n and the step size, and calculate the total width W; S370: If the total width W satisfies the third constraint, calculate the runner area S based on all combinations that satisfy the third constraint. If the new runner area S is greater than the currently recorded maximum value, update the runner area S until all traversals are completed to obtain the optimal solution for the runner area S, and record the corresponding values of runner width a, rib width b, forming angle θ, hypotenuse length d, and runner number n; The formula for calculating the area S of each flow channel is: a d sinθ+d 2 sinθ cosθ, where a d sinθ represents the area of the rectangular portion of each flow channel in the cross-section of the pure titanium bipolar plate, d 2 sinθ cosθ represents the sum of the areas of the two triangular flanks of each flow channel in the cross-section of a pure titanium bipolar plate.
2. The method for optimizing the cross-sectional design of the bipolar plate flow channel in a hydrogen fuel cell according to claim 1, characterized in that, In step S100, the value of the runner width a ranges from 0.2 to 2, the value of the rib width b ranges from 0.7 to 1, the value of the forming angle ranges from π / 3 to π / 2, the value of the inclined side length d ranges from 0 to 2, and the value of the number of runners n ranges from 20 to 50.
3. The method for optimizing the cross-sectional design of the bipolar plate flow channel in a hydrogen fuel cell according to claim 2, characterized in that, In step S100, the step size of the runner width a is 0.05, the step size of the rib width b is 0.05, the step size of the forming angle is π / 180, the step size of the inclined side length d is 0.1, and the step size of the number of runners n is 1.
4. The method for optimizing the cross-sectional design of the bipolar plate flow channel in a hydrogen fuel cell according to claim 1, characterized in that, In step S200, the ratio of the width a of the first constraint channel to the width b of the rib is in the range of 0.25≤a / b≤2.5, so as to balance the structural strength and fluid performance of the pure titanium bipolar plate.
5. The method for optimizing the flow channel cross-section of a hydrogen fuel cell bipolar plate according to claim 1, characterized in that, In step S300, the maximum value of the flow channel area S is 1.0132, where the flow channel width a is 0.85, the rib width b is 0.70, the forming angle θ is 60 degrees, the inclined side length d is 0.90, and the number of flow channels n is 20.
6. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the steps of the method for optimizing the flow channel cross-section of a hydrogen fuel cell bipolar plate as described in any one of claims 1-5.