Bipolar plate runner structure design method

Abstract

The invention provides a bipolar plate flow channel structure design method which comprises the following steps: S1, based on a bipolar plate snakelike flow channel design formula, a plurality of dimensionless parameter sets are obtained by adjusting parameters in the formula, and each dimensionless parameter set represents one snakelike flow channel structure of a bipolar plate; s2, according to each dimensionless parameter set, establishing a bipolar plate model of a corresponding runner structure; performing analogue simulation on each bipolar plate model to obtain a flow velocity distribution cloud picture of each bipolar plate model; s3, determining a target bipolar plate model according to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud picture of each bipolar plate model; and according to the target bipolar plate model, obtaining design data of the to-be-designed bipolar plate runner structure. The optimization design parameters of the snake-shaped flow channel can be quickly and accurately determined, and the design efficiency is high.

Description

Technical Field

[0001] The invention relates to the technical field of fuel cells and water electrolysis hydrogen production, and in particular to a bipolar plate flow channel structure design method. Background Art

[0002] The flow channel design of the bipolar plate is particularly important for fuel cell stacks and PEM water electrolysis hydrogen production stacks. When the reaction is in progress, both of the above devices carry out mass and heat transfer of gas-liquid two-phase flow in the flow channel and diffusion layer. In order to ensure the normal and efficient electrochemical reaction in the stack, the gas-liquid two-phase flow needs to be evenly distributed on the membrane electrode side, that is, the gas-liquid two-phase flow needs to be evenly distributed in the diffusion layer and the bipolar plate flow channel. The most effective way to make the water distribution uniform is to design an efficient flow channel.

[0003] At present, the more popular bipolar plate flow channel design shapes include: parallel flow channels, nipple or square convex structure flow channels, serpentine flow channels, etc. These flow channels are widely used due to their simple structure and easy processing and manufacturing. However, when these flow channels are used in fuel cell stacks and PEM water electrolysis hydrogen production stacks, there will be problems of uneven distribution of gas-liquid two-phase flow in the flow channels to varying degrees.

[0004] For the bipolar plate serpentine flow channel, due to the limited width of the inlet main pipe, the large number of flow channel branches and the wide bipolar plate, the flow field width will be much larger than the width of the inlet main pipe, resulting in the gas-liquid two-phase flow being difficult to be evenly distributed on the activation surface in various tributaries of the flow channel, which will lead to different electrochemical reaction rates. Insufficient water vapor will cause the reaction rate to decrease or even stop, and even cause local overheating and dry burning, which will irreversibly reduce the life of the battery stack. In addition, in the PEM water electrolysis hydrogen production device, the nature of the oxygen at the anode makes it difficult for oxygen to escape as easily as hydrogen. Uneven water distribution in the bipolar plate flow channel will make it easier for oxygen to adsorb on the anode surface, making it difficult for water to contact the anode to continuously undergo oxygen evolution reaction, thereby greatly reducing the reaction rate.

[0005] Therefore, it is necessary to optimize the design of the bipolar plate serpentine flow channel to improve the uniformity of water distribution in the flow channel, thereby improving the working efficiency of the fuel cell stack and the PEM water electrolysis stack device. However, the existing design method of the bipolar plate serpentine flow channel structure is difficult to quickly and accurately determine the optimized design parameters of the serpentine flow channel in order to achieve the desired uniformity of water distribution in the flow channel. The design process is redundant and cumbersome, and the workload is large.

[0006] Therefore, an improved bipolar plate flow channel structure design method is urgently needed. Summary of the invention

[0007] In view of the above-mentioned shortcomings and deficiencies of the prior art, the present invention provides a bipolar plate flow channel structure design method, which can quickly and accurately determine the optimized design parameters of the serpentine flow channel and has high design efficiency.

[0008] In order to achieve the above object, the main technical solutions adopted by the present invention include:

[0009] An embodiment of the present invention provides a method for designing a bipolar plate flow channel structure, comprising the following steps:

[0010] S1. Based on the bipolar plate serpentine flow channel design formula, multiple dimensionless parameter sets are obtained by adjusting the parameters in the formula, and each dimensionless parameter set represents a serpentine flow channel structure of the bipolar plate;

[0011] The design formula of the bipolar plate serpentine flow channel is:

[0012] L=[n·d+(n-1)·m]·(2t+1)+2·t·w

[0013] Wherein, L is the total width of the serpentine flow channel; n is the number of parallel channels of the serpentine flow channel; d is the width of a parallel channel; m is the interval width between adjacent parallel channels in the serpentine flow channel; t is the number of bends of the serpentine flow channel, and the serpentine flow channel bends upward and then bends downward as one bend; w is the interval width between the serpentine flow channels;

[0014] The dimensionless parameter set is denoted as (n, d, m, t, w);

[0015] S2. According to each dimensionless parameter set, a bipolar plate model corresponding to the flow channel structure is established; each bipolar plate model is simulated to obtain a flow velocity distribution cloud map of each bipolar plate model;

[0016] S3. Determine a target bipolar plate model according to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud diagram of each bipolar plate model; and obtain design data of the bipolar plate flow channel structure to be designed according to the target bipolar plate model.

[0017] Optionally, in S1, based on the bipolar plate serpentine flow channel design formula, the parameters in the formula are adjusted according to the bipolar plate size and a pre-set parameter value range to obtain multiple dimensionless parameter sets.

[0018] Optionally, in S2, each bipolar plate model is simulated using Ansys, COMSOL or SolidWorks software.

[0019] Optionally, in S2, the K-epsilon standard turbulence model is selected, the Coupled coupling solution algorithm is adopted, the convergence standard of the iterative residual is set to 1*e-3, and the simulation is run for each bipolar plate model to obtain the flow velocity distribution cloud map of each bipolar plate model.

[0020] Optionally, the input parameters for running the simulation of each bipolar plate model are the same; the input parameters include the bipolar plate model liquid inlet, the bipolar plate model liquid outlet, fluid viscosity, temperature, initial flow rate and pressure.

[0021] Optionally, a plane located in the middle of the plane where the top of the flow channel is located and the plane where the bottom of the flow channel is located is used as the generation section of the flow velocity distribution cloud map.

[0022] Optionally, in S2, each bipolar plate model is simulated to obtain a flow velocity distribution cloud map, a pressure distribution map and a streamline map of each bipolar plate model;

[0023] Correspondingly, in S3, the target bipolar plate model is determined based on the target water distribution uniformity, target pressure distribution uniformity of the bipolar plate flow channel to be designed, and the flow velocity distribution cloud map, pressure distribution map, and streamline map of each bipolar plate model, and the design data of the bipolar plate flow channel structure to be designed is obtained based on the target bipolar plate model.

[0024] Optionally, the target bipolar plate model is determined according to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud diagram of each bipolar plate model, including:

[0025] According to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud map of each bipolar plate model, a bipolar plate model that meets the target water distribution uniformity is obtained. If there is only one bipolar plate model that meets the target water distribution uniformity, then the bipolar plate model is used as the target bipolar plate model. If there are more than two bipolar plate models that meet the target water distribution uniformity, the flow velocity distribution cloud maps of the two or more bipolar plate models are compared, and the bipolar plate model whose average flow velocity meets the requirements is selected as the target bipolar plate model.

[0026] Optionally, the flow channel corner in the target bipolar plate model is a right angle, and S3 also includes: changing the flow channel corner in the target bipolar plate model to a rounded corner, obtaining a rounded corner target bipolar plate model, simulating the rounded corner target bipolar plate model, obtaining a flow velocity distribution cloud map of the rounded corner target bipolar plate model, and determining the design data of the bipolar plate flow channel structure to be designed according to the flow velocity distribution cloud map of the target bipolar plate model and the flow velocity distribution cloud map of the rounded corner target bipolar plate model.

[0027] The beneficial effects of the present invention are:

[0028] The bipolar plate flow channel structure design method proposed in the present invention characterizes the bipolar plate serpentine flow channel structure through a bipolar plate serpentine flow channel design formula, can accurately describe the geometric shape, size and layout of the serpentine flow channel, and uses the formula to perform parameterized design and optimization of the flow channel structure. It can perform numerical simulation and simulation analysis on a computer to improve the flow properties of the fluid, greatly reduce the number and cost of experiments, shorten the design cycle, and improve design efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] Figure 1 A schematic diagram of a flow chart of a bipolar plate flow channel structure design method according to a specific embodiment of the present invention;

[0030] Figure 2 It is a schematic diagram of characterizing the bipolar plate serpentine flow channel according to the bipolar plate serpentine flow channel design formula according to a specific embodiment of the present invention;

[0031] Figure 3 This is a schematic diagram of the structure of the flow channel of the n2d1m1t6w1 type bipolar plate according to Example 1 of the present invention;

[0032] Figure 4 It is a schematic diagram of the structure of the flow channel of the N7D0.5M0.5T3W0.5 bipolar plate according to Example 1 of the present invention;

[0033] Figure 5 This is a schematic diagram of the structure of the flow channel of the n8d1m1t1w1 type bipolar plate according to Example 1 of the present invention;

[0034] Figure 6 A flow velocity distribution cloud diagram of the flow channel of the n2d1m1t6w1 type bipolar plate according to Example 1 of the present invention;

[0035] Figure 7 The flow velocity distribution cloud diagram of the flow channel of the N3D0.5M0.5T7W0.5 bipolar plate according to Example 1 of the present invention;

[0036] Figure 8 Flow velocity distribution cloud diagram of the flow channel of the n8d1m1t1w1 type bipolar plate according to Example 1 of the present invention;

[0037] Fig. 9 Schematic diagram of the structure of the rounded n8d1m1t1w1 bipolar plate flow channel according to Example 1 of the present invention;

[0038] Fig.10 Flow velocity distribution cloud diagram of the rounded corner n8d1m1t1w1 bipolar plate flow channel according to Example 1 of the present invention;

[0039] Fig.11 is a cell voltage-current density curve according to Example 1 of the present invention, which includes n2d1m1t6w1 type bipolar plate flow channels, n7d0.5m0.5t3w0.5 type bipolar plate flow channels, n8d1m1t1w1 type bipolar plate flow channels, rounded n8d1m1t1w1 type bipolar plate flow channels and cube d1m1 type flow channels;

[0040] Fig.12 It is a schematic diagram of the structure of a cube d1m1 type flow channel according to comparative example 1 of the present invention;

[0041] Fig.13 is a flow velocity distribution cloud diagram of a cube d1m1 type flow channel according to comparative example 1 of the present invention;

[0042] Fig.14 It is a schematic structural diagram of a horizontal d1m1 type flow channel according to Comparative Example 2 of the present invention;

[0043] Fig.15 It is a schematic structural diagram of a vertical d1m1 type flow channel according to Comparative Example 2 of the present invention;

[0044] Fig.16 is a flow velocity distribution cloud diagram of the horizontal d1m1 type flow channel according to Comparative Example 2 of the present invention;

[0045] Fig.17 Graph showing the flow velocity distribution of the vertical d1m1 type flow channel according to Comparative Example 2 of the present invention. DETAILED DESCRIPTION

[0046] In order to better explain the present invention and facilitate understanding, the present invention is described in detail below through specific implementation modes in conjunction with the accompanying drawings.

[0047] like Figure 1 As shown, the present invention provides a method for designing a bipolar plate flow channel structure, comprising the following steps:

[0048] S1. Based on the bipolar plate serpentine flow channel design formula, multiple dimensionless parameter sets are obtained by adjusting the parameters in the formula, and each dimensionless parameter set represents a serpentine flow channel structure of the bipolar plate.

[0049] The design formula of the bipolar plate serpentine flow channel is:

[0050] L=[n·d+(n-1)·m]·(2t+1)+2·t·w

[0051] In the formula, Figure 2 As shown, L is the total width of the serpentine flow channel; n is the number of parallel channels of the serpentine flow channel; d is the width of a parallel channel; m is the spacing width between adjacent parallel channels in the serpentine flow channel; t is the number of bends of the serpentine flow channel, and a serpentine flow channel that bends upward and then downward is counted as one bend; w is the spacing width between serpentine flow channels.

[0052] The dimensionless parameter set is denoted as (n, d, m, t, w).

[0053] In this way, the bipolar plate serpentine flow channel design formula proposed in the present invention is used to characterize the bipolar plate serpentine flow channel structure, which can accurately describe the geometric shape, size and layout of the serpentine flow channel. With the help of this formula, the flow channel structure can be parameterized and optimized to improve the flow properties of the fluid. With the help of this formula, it is convenient to perform numerical simulation and simulation analysis on a computer, thereby greatly reducing the number and cost of experiments, shortening the design cycle, and improving design efficiency.

[0054] Preferably, based on the bipolar plate serpentine flow channel design formula, the parameters in the formula are adjusted according to the bipolar plate size and the pre-set parameter value range to obtain multiple dimensionless parameter sets. The pre-set parameter value range is generally reasonably set by the designer after consulting the literature, patents and experimental data on the design of the bipolar plate serpentine flow channel. By adjusting the parameters in the formula within the pre-set parameter value range, the parametric design of the flow channel structure is made more reasonable and targeted, thereby improving the accuracy and efficiency of the design.

[0055] S2. According to each dimensionless parameter set, a bipolar plate model of a corresponding flow channel structure is established; each bipolar plate model is simulated to obtain a flow velocity distribution cloud diagram of each bipolar plate model.

[0056] Among them, Ansys software was used, K-epsilon standard turbulence model was selected, Coupled coupling solution algorithm was used, the convergence standard of iterative residual was set to 1*e-3, simulation was run for each bipolar plate model, and the velocity distribution cloud map of each bipolar plate model was obtained. In this way, the flow behavior of the fluid in the bipolar plate flow channel structure can be accurately simulated.

[0057] Preferably, the plane located in the middle of the plane where the top of the flow channel is located and the plane at the bottom of the flow channel is used as the section for generating the velocity distribution cloud map. In this way, the velocity distribution cloud map obtained is the most accurate. The velocity distribution cloud map generated by the section at other positions may cause turbulence deviation due to the plane being close to the plane at the bottom of the flow channel or the diffusion layer.

[0058] The input parameters for each bipolar plate model simulation are the same, including the bipolar plate model liquid inlet, bipolar plate model liquid outlet, fluid viscosity, temperature, initial flow rate and pressure. This makes it easier to compare the flow rate distribution cloud map of each bipolar plate model and determine the target bipolar plate model.

[0059] Optionally, each bipolar plate model is simulated using COMSOL or SolidWorks software.

[0060] S3. Determine a target bipolar plate model according to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud diagram of each bipolar plate model; and obtain design data of the bipolar plate flow channel structure to be designed according to the target bipolar plate model.

[0061] Specifically, the target bipolar plate model is determined according to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud map of each bipolar plate model, including: according to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud map of each bipolar plate model, a bipolar plate model that meets the target water distribution uniformity is obtained; if there is only one bipolar plate model that meets the target water distribution uniformity, the bipolar plate model is used as the target bipolar plate model; if there are more than two bipolar plate models that meet the target water distribution uniformity, the flow velocity distribution cloud maps of the two or more bipolar plate models are compared, and the bipolar plate model whose average flow velocity meets the requirements is selected as the target bipolar plate model.

[0062] Preferably, S3 also includes: changing the flow channel corner in the target bipolar plate model to a rounded corner, obtaining a rounded target bipolar plate model, simulating the rounded target bipolar plate model, obtaining a flow velocity distribution cloud map of the rounded target bipolar plate model, and determining the design data of the bipolar plate flow channel structure to be designed according to the flow velocity distribution cloud map of the target bipolar plate model and the flow velocity distribution cloud map of the rounded target bipolar plate model.

[0063] In summary, the bipolar plate flow channel structure design method proposed in the present invention characterizes the bipolar plate serpentine flow channel structure through the bipolar plate serpentine flow channel design formula, can accurately describe the geometric shape, size and layout of the serpentine flow channel, and uses the formula to perform parameterized design and optimization of the flow channel structure. It can perform numerical simulation and simulation analysis on a computer to improve the flow properties of the fluid, greatly reduce the number and cost of experiments, shorten the design cycle, and improve design efficiency.

[0064] In order to further optimize the bipolar plate flow channel structure design, in S2, each bipolar plate model is simulated to obtain the flow velocity distribution cloud map, pressure distribution map and streamline map of each bipolar plate model. Correspondingly, in S3, the target bipolar plate model is determined according to the target water distribution uniformity and target pressure distribution uniformity of the bipolar plate flow channel to be designed, and the flow velocity distribution cloud map, pressure distribution map and streamline map of each bipolar plate model, and the design data of the bipolar plate flow channel structure to be designed is obtained according to the target bipolar plate model.

[0065] Among them, the flow velocity distribution cloud map can intuitively observe and compare the distribution and flow velocity of water at each flow channel position, the pressure cross-section diagram can analyze and compare the distribution of water pressure at each flow channel position, and the streamline diagram can intuitively observe and compare the movement of the liquid simulated flow trajectory in the flow channel.

[0066] Example 1

[0067] The order of installation of bipolar plates in fuel cells or PEM water electrolysis devices is: bipolar plate, diffusion layer, catalyst layer, proton exchange membrane, and the other side is a symmetrical structure. The design process of the bipolar plate flow channel structure provided in this embodiment is:

[0068] S1. Based on the bipolar plate serpentine flow channel design formula L = [n·d+(n-1)·m]·(2t+1)+2·t·w, the parameters in the formula are adjusted according to the bipolar plate size and the pre-set parameter value range to obtain multiple dimensionless parameter sets, each dimensionless parameter set characterizing a serpentine flow channel structure of the bipolar plate.

[0069] Among them, the bipolar plate is square, the flow channel depth is 1, the width of the bipolar plate is L<55, and the pre-set parameter value ranges include: 1≤n≤10, 0.5≤d≤2, 0.5≤m≤2, 1≤t≤8, 0.5≤w≤2.

[0070] S2. According to each dimensionless parameter set, a bipolar plate model of the corresponding flow channel structure is established; using Ansys software, the K-epsilon standard turbulence model is selected, the Coupled coupling solution algorithm is used, the convergence standard of the iterative residual is set to 1*e-3, and simulation is run for each bipolar plate model to obtain the flow velocity distribution cloud map of each bipolar plate model.

[0071] Among them, the plane located in the middle of the plane where the top of the flow channel is located and the plane at the bottom of the flow channel is used as the generation section of the flow velocity distribution cloud map; the input parameters for the simulation of each bipolar plate model are the same, namely: the liquid inlet of the bipolar plate model is the lower right corner port, the liquid outlet of the bipolar plate model is the upper left corner port, the fluid is water, the fluid viscosity is the viscosity of water, the temperature is room temperature, the initial flow rate is 50mL / min, and the pressure is standard atmospheric pressure.

[0072] S3. According to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud diagram of each bipolar plate model, three bipolar plate models that meet the target water distribution uniformity are obtained. The dimensionless parameter sets corresponding to these three bipolar plate models are (2, 1, 1, 6, 1), (7, 0.5, 0.5, 3, 0.5), and (8, 1, 1, 1, 1), respectively. Therefore, these three bipolar plate models are named as n2d1m1t6w1 type bipolar plate flow channel (such as Figure 3 As shown), N7D0.5M0.5T3W0.5 bipolar plate flow channel (as shown Figure 4 As shown), n8d1m1t1w1 bipolar plate flow channel (as shown Figure 5 As shown), the velocity distribution cloud diagram of the flow channel of the n2d1m1t6w1 type bipolar plate is as follows Figure 6As shown, the average flow velocity is about 1.206 m / s, and the velocity distribution cloud diagram of the flow channel of the n3d0.5m0.5t7w0.5 type bipolar plate is as follows Figure 7 As shown, the average flow velocity is about 0.546 m / s. The velocity distribution cloud diagram of the flow channel of the n8d1m1t1w1 bipolar plate is as follows Figure 8 As shown, the average flow velocity is about 0.303 m / s. By comparing the flow velocity distribution cloud diagrams of the three bipolar plate models, the N8D1M1T1W1 bipolar plate flow channel structure with an average flow velocity that meets the requirements (0.3 m / s < average flow velocity < 0.4 m / s) is selected as the target bipolar plate model.

[0073] Change the channel corners in the target bipolar plate model to rounded corners (such as Fig. 9 As shown), a rounded corner target bipolar plate model is obtained, and the rounded corner target bipolar plate model is simulated to obtain a flow velocity distribution cloud diagram of the rounded corner target bipolar plate model (as shown Fig.10 As shown), the average flow rate is about 0.317 m / s. According to the flow rate distribution cloud map of the target bipolar plate model and the flow rate distribution cloud map of the rounded target bipolar plate model, the design data of the bipolar plate flow channel structure to be designed is determined, that is, the design data of the bipolar plate flow channel structure to be designed is the rounded n8d1m1t1w1 type bipolar plate flow channel structure.

[0074] The process of rounded corners is easier to implement than the completely right-angle processing, whether in etching or machining. Since the average flow velocity of the rounded corner n8d1m1t1w1 bipolar plate flow channel structure also meets the requirements (0.3m / s<average flow velocity<0.4m / s), the design data of the bipolar plate flow channel structure to be designed is the rounded corner n8d1m1t1w1 bipolar plate flow channel structure.

[0075] According to the flow velocity distribution cloud diagrams of the above-mentioned N2D1M1T6W1 bipolar plate flow channels, N7D0.5M0.5T3W0.5 bipolar plate flow channels, and N8D1M1T1W1 bipolar plate flow channels, it can be seen that when the number of flow channels increases and the number of bends decreases, the flow velocity of the liquid in the flow channel will be significantly reduced. This conclusion is consistent with the Bernoulli equation formula.

[0076] The N2D1M1T6W1 bipolar plate flow channel, N7D0.5M0.5T3W0.5 bipolar plate flow channel, and N8D1M1T1W1 bipolar plate flow channel were etched, respectively. The substrate used was a 3mm thick TA1 titanium plate, and the selected diffusion layer was titanium felt. The gasket was EPDM rubber. The anode used was an iridium-based catalyst with a loading of 1.5 mg / cm 2 ; The cathode is a platinum-based catalyst with a loading of 0.4 mg / cm 2The selected proton exchange membrane is Chemours N117. The above three bipolar plate flow channels and other components are assembled into electrolyzers, and then the three electrolyzers are tested with constant current. The test parameters are as follows: the electrolyte is pure water with a resistivity of 18.25MΩ·cm, the flow rate is 50mL / min, and the electrolysis temperature is 60℃; the cell voltage-current density curve is as follows Fig.11 shown.

[0077] According to the cell voltage-current density curve, under the same current density, V n8d1m1t1w1 >V n7d0.5m0.5t3w0.5 >V n2d1m1t6w1 , the cell pressure decreases with the increase of liquid flow rate, that is, it decreases with the decrease of n value and the increase of t value, thereby improving the efficiency of the electrolytic cell, which corresponds to the flow rate information of the liquid in the flow channel.

[0078] Comparative Example 1

[0079] Arrange and design the flow channels according to regular convexities, such as Fig.12 As shown, the spacing and the width of the protrusions are both 1, which is defined as a cube d1m1 type flow channel. The simulation calculation of the above embodiment is also performed on this flow channel, and the input and output parameters are the same, and it can be obtained Fig.13 The velocity distribution cloud diagram shown. The uniformity of the velocity distribution in the overall flow channel is lower than that of the flow channel design in Example 1. The liquid flow branches too much at the inlet, and due to different directions and orientations, the flow velocity in each direction of the inlet branch is different, and turbulent vortices are formed at many points. In this case, the gas will stay in the vortex and is not easy to discharge. In the corners of the flow channel, there is an area where the liquid flow cannot be distributed at all, which reduces the effective area. Under the same current input, the effective current density will increase, which will increase the slot pressure and reduce the service life of the electrode.

[0080] Fig.11 The cell pressure-current density curve in also shows that the cell pressure is higher than that of the above-mentioned embodiment. Compared with the flow channel arrangement design of embodiment 1, this flow channel is not a bipolar plate flow channel design with strong practicality.

[0081] Comparative Example 2

[0082] According to the formula L = [n·d+(n-1)·m]·(2t+1)+2·t·w, if t = 0, a completely parallel or vertical flow channel with n as the maximum value is formed, and the flow channel width and the spacing between the flow channels are both 1, such as Fig.14 and Fig.15 As shown, these two flow channel designs are named horizontal d1m1 type and vertical d1m1 type flow channels.

[0083] The simulation calculation of the above embodiment is also performed on this flow channel, and the input and output parameters are the same as those in the embodiment, and it can be obtained Fig.16 and Fig.17 The flow velocity distribution diagram is shown. It can be seen that the flow channel design in this case will make the distribution of the liquid in the flow channel more uneven. Except for the faster flow velocity of the flow channel at the edge, the liquid flow velocity in the parallel or vertical flow channel in the middle part is almost close to 0. Similar to the defects in comparative example 1, the effective electrolysis area will be reduced. Under the same current input, the effective current density will increase, the cell pressure will increase, and the electrode life will also be reduced.

[0084] In the description of the present invention, it should be understood that the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of the features. In the description of the present invention, the meaning of "plurality" is two or more, unless otherwise clearly and specifically defined.

[0085] In the present invention, unless otherwise clearly specified and limited, the terms "installed", "connected", "connected", "fixed" and the like should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium; it can be the internal connection of two elements or the interaction relationship between two elements. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0086] In the present invention, unless otherwise clearly specified and limited, when a first feature is “on” or “below” a second feature, it may be that the first and second features are in direct contact, or that the first and second features are in indirect contact through an intermediate medium. Moreover, when a first feature is “above”, “above” or “above” a second feature, it may be that the first feature is directly above or obliquely above the second feature, or it may simply mean that the first feature is higher in level than the second feature. When a first feature is “below”, “below” or “below” a second feature, it may be that the first feature is directly below or obliquely below the second feature, or it may simply mean that the first feature is lower in level than the second feature.

[0087] In the description of this specification, the description of the terms "one embodiment", "some embodiments", "embodiment", "example", "specific example" or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are contradictory.

[0088] Although the embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and are not to be construed as limitations of the present invention. A person skilled in the art may alter, modify, replace and modify the above embodiments within the scope of the present invention.

Claims

1. A method for designing a bipolar plate flow channel structure, characterized in that: The following steps are involved: S1. Based on the bipolar plate serpentine flow channel design formula, multiple dimensionless parameter sets are obtained by adjusting the parameters in the formula, and each dimensionless parameter set represents a serpentine flow channel structure of the bipolar plate; The design formula of the bipolar plate serpentine flow channel is: L=[n·d+(n-1)·m]·(2t+1)+2·t·w Wherein, L is the total width of the serpentine flow channel; n is the number of parallel channels of the serpentine flow channel; d is the width of a parallel channel; m is the interval width between adjacent parallel channels in the serpentine flow channel; t is the number of bends of the serpentine flow channel, and the serpentine flow channel bends upward and then bends downward as one bend; w is the interval width between the serpentine flow channels; The dimensionless parameter set is denoted as (n, d, m, t, w); S2. According to each dimensionless parameter set, a bipolar plate model corresponding to the flow channel structure is established; each bipolar plate model is simulated to obtain a flow velocity distribution cloud map of each bipolar plate model; S3. Determine a target bipolar plate model according to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud diagram of each bipolar plate model; and obtain design data of the bipolar plate flow channel structure to be designed according to the target bipolar plate model.

2. The method for designing a bipolar plate flow channel structure according to claim 1, characterized in that: In S1, based on the bipolar plate serpentine flow channel design formula, the parameters in the formula are adjusted according to the bipolar plate size and the pre-set parameter value range to obtain multiple dimensionless parameter sets.

3. The method for designing a bipolar plate flow channel structure according to claim 1, characterized in that: In S2, each bipolar plate model is simulated using Ansys, COMSOL or SolidWorks software.

4. The method for designing a bipolar plate flow channel structure according to claim 1 or 3, characterized in that: In S2, the K-epsilon standard turbulence model is selected, the Coupled coupling solution algorithm is used, the convergence standard of the iterative residual is set to 1*e-3, and the simulation is run for each bipolar plate model to obtain the flow velocity distribution cloud map of each bipolar plate model.

5. The method for designing a bipolar plate flow channel structure according to claim 4, characterized in that: The input parameters for each bipolar plate model simulation run are the same; The input parameters include the bipolar plate model inlet, the bipolar plate model outlet, fluid viscosity, temperature, initial flow rate and pressure.

6. The method for designing a bipolar plate flow channel structure according to claim 1, characterized in that: The plane located in the middle of the plane where the top of the flow channel is located and the plane of the bottom of the flow channel is used as the generating section of the flow velocity distribution cloud map.

7. The method for designing a bipolar plate flow channel structure according to claim 1, characterized in that: In S2, each bipolar plate model is simulated to obtain a flow velocity distribution cloud map, a pressure distribution map and a streamline map of each bipolar plate model; Correspondingly, in S3, the target bipolar plate model is determined based on the target water distribution uniformity, target pressure distribution uniformity of the bipolar plate flow channel to be designed, and the flow velocity distribution cloud map, pressure distribution map, and streamline map of each bipolar plate model, and the design data of the bipolar plate flow channel structure to be designed is obtained based on the target bipolar plate model.

8. The method for designing a bipolar plate flow channel structure according to claim 1, characterized in that: According to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud diagram of each bipolar plate model, the target bipolar plate model is determined, including: According to the target water distribution uniformity of the bipolar plate flow channel to be designed and the flow velocity distribution cloud map of each bipolar plate model, a bipolar plate model that meets the target water distribution uniformity is obtained. If there is only one bipolar plate model that meets the target water distribution uniformity, then the bipolar plate model is used as the target bipolar plate model. If there are more than two bipolar plate models that meet the target water distribution uniformity, the flow velocity distribution cloud maps of the two or more bipolar plate models are compared, and the bipolar plate model whose average flow velocity meets the requirements is selected as the target bipolar plate model.

9. The method for designing a bipolar plate flow channel structure according to claim 1 or 8, characterized in that: The flow channel corners in the target bipolar plate model are right angles, and S3 also includes: The flow channel corner in the target bipolar plate model is changed to a rounded corner to obtain the rounded target bipolar plate model, the rounded target bipolar plate model is simulated to obtain the flow velocity distribution cloud map of the rounded target bipolar plate model, and the design data of the bipolar plate flow channel structure to be designed is determined according to the flow velocity distribution cloud map of the target bipolar plate model and the flow velocity distribution cloud map of the rounded target bipolar plate model.

Images