CN112800643A - Multi-physical-field coupling calculation simplification method for corrugated-runner fuel cell - Google Patents
Multi-physical-field coupling calculation simplification method for corrugated-runner fuel cell Download PDFInfo
- Publication number
- CN112800643A CN112800643A CN202011615085.6A CN202011615085A CN112800643A CN 112800643 A CN112800643 A CN 112800643A CN 202011615085 A CN202011615085 A CN 202011615085A CN 112800643 A CN112800643 A CN 112800643A
- Authority
- CN
- China
- Prior art keywords
- fuel cell
- cutting section
- flow channel
- bipolar plate
- area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000000446 fuel Substances 0.000 title claims abstract description 50
- 238000004364 calculation method Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000008878 coupling Effects 0.000 title claims abstract description 17
- 238000010168 coupling process Methods 0.000 title claims abstract description 17
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 17
- 239000012530 fluid Substances 0.000 claims abstract description 65
- 239000012528 membrane Substances 0.000 claims abstract description 41
- 230000002457 bidirectional effect Effects 0.000 claims abstract description 27
- 238000009792 diffusion process Methods 0.000 claims abstract description 27
- 238000006243 chemical reaction Methods 0.000 claims abstract description 20
- 238000004458 analytical method Methods 0.000 claims abstract description 11
- 238000012937 correction Methods 0.000 claims abstract description 10
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 24
- 239000001257 hydrogen Substances 0.000 claims description 23
- 229910052739 hydrogen Inorganic materials 0.000 claims description 23
- 239000000110 cooling liquid Substances 0.000 claims description 21
- 239000000463 material Substances 0.000 claims description 9
- 230000000737 periodic effect Effects 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000006555 catalytic reaction Methods 0.000 claims description 6
- 239000012495 reaction gas Substances 0.000 claims description 6
- 230000017525 heat dissipation Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000010205 computational analysis Methods 0.000 claims description 2
- 239000002826 coolant Substances 0.000 description 7
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 239000003054 catalyst Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003487 electrochemical reaction Methods 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/30—Hydrogen technology
- Y02E60/50—Fuel cells
Abstract
The invention discloses a corrugated flow channel fuel cell multi-physical field coupling calculation simplification method, which comprises the following steps: establishing a geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell; establishing a grid with completely consistent upper cutting section and lower cutting section of the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell, performing Boolean operation on the bipolar plate of the full-area single cell, and respectively obtaining the fluid flow of each channel in a reaction region after three-cavity fluid is distributed by a distribution region; carrying out structural mechanics analysis on the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell to obtain new membrane electrode diffusion layers in different deformation states at different positions after the bidirectional corrugated flow channel bipolar plate is compressed; and combining the new membrane electrode diffusion layer with the new air fluid domain to replace the membrane electrode structure and the air fluid domain of the established geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell, and obtaining a correction model of the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell.
Description
Technical Field
The invention relates to the technical field of fuel cells, in particular to a multi-physical-field coupling calculation simplification method for a corrugated flow channel fuel cell.
Background
Proton exchange membrane fuel cells have recently taken an important position in the energy application field of the global future as a clean, high-efficiency electrochemical energy conversion device using renewable energy, hydrogen energy. The proton exchange membrane fuel cell has complex processes such as multi-scale, multi-phase, multi-mass field, multi-component heat and mass transfer and the like, and besides the diagnosis and analysis of the complex physical process in the proton exchange membrane fuel cell by using experimental means, a reasonable calculation model is established by means of mathematical means so as to deeply understand various transmission mechanisms of the complex physical process in the cell, so that the proton exchange membrane fuel cell is widely applied at present.
According to different dimensions established by the calculation model, the calculation model of the proton exchange membrane fuel cell can be divided into a one-dimensional model, a two-dimensional model and a three-dimensional model. One-dimensional models are often used for simple qualitative analysis and can output some simple rules in the transmission process inside the battery. In order to take account of the calculation efficiency and the accuracy, the one-dimensional model is superposed in another dimension to obtain a two-dimensional model, or the other two dimensions are superposed simultaneously to obtain a three-dimensional model. The three-dimensional model can establish a real single-cell geometric structure calculation domain, can solve a real fluid flow and component transmission process, can mutually couple and solve the real fluid flow and component transmission process with other physical processes, and is widely applied to the design and development of proton exchange membrane fuel cells.
Because the three-dimensional model takes the real flow field structure into consideration and involves a multi-physical field coupling solving process, the actual full-area single cell structure is established to analyze the influence of the flow field structure design on the battery performance, the model calculation amount can be greatly increased, and the calculation stability can be correspondingly reduced. Therefore, due to the limitation of computing resources, reasonable simplification work cannot be avoided in the process of establishing the three-dimensional computing model.
Especially, application number (CN107145658A) in the prior art discloses a numerical simulation method for designing flow field parameters of a bipolar plate of a proton exchange membrane fuel cell, which includes firstly performing simulation analysis on stamping formability of the bipolar plate through Dynaform to obtain a proper stamping forming range, then performing flow field design of the bipolar plate, and finally simulating output performance of the cell through CFD to obtain a flow field size of the bipolar plate with the best performance. The influence of punch forming of the bipolar plate is considered in the process, but the structural force is not considered to analyze the compression condition of the MEA, the compression characteristic of the MEA can directly influence the flow field condition of the bipolar plate, and further the performance of the battery is influenced.
In addition, in the prior art, part of calculation models are low-dimensional models, and the real geometric structure of the battery is often abstracted into physical parameters such as length, thickness, height and the like, so that the real flowing process in the battery and the spatial distribution rule of each physical quantity cannot be reflected. The other part of the model is a three-dimensional model, a full-area single cell is often simplified into a single straight flow channel model for qualitative calculation analysis, the influence brought by a real flow field structure is not considered, and the influence of structural force on the compression characteristic of the MEA and the flow field structure is not considered in the established multi-physical field coupling model.
Disclosure of Invention
According to the problems in the prior art, the invention discloses a method for simplifying the multi-physical-field coupling calculation of a corrugated flow channel fuel cell, which specifically comprises the following steps:
establishing a geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell;
carrying out structural mechanics analysis on the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell to obtain new membrane electrode diffusion layers in different deformation states at different positions after the bidirectional corrugated flow channel bipolar plate is compressed;
combining the new membrane electrode diffusion layer with the bipolar plate, and obtaining new air fluid domains at different positions of the bipolar plate due to the embedding influence of the membrane electrode diffusion layer through Boolean operation; combining the new membrane electrode diffusion layer with the new air fluid domain to replace the membrane electrode structure and the air fluid domain of the established geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell, and obtaining a correction model of the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell;
establishing corresponding non-common node interfaces of the correction model on a cathode and anode diffusion layer and cathode and anode reaction gas interface, a cathode and anode reaction gas and cathode and anode cooling liquid interface, a cathode and anode cooling liquid central interface and a cathode and anode cooling liquid and cathode and anode collector plate interface respectively; setting periodic boundary conditions between corresponding layers of materials on the upper cutting section and the lower cutting section respectively, and establishing a grid with the completely consistent upper cutting section and lower cutting section of the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell;
performing Boolean operation on the bipolar plate of the full-area single cell to respectively obtain three-cavity fluid domains of an air cavity, a hydrogen cavity and a cooling liquid cavity of the full-area single cell, and respectively performing three-cavity fluid distribution calculation to respectively obtain the fluid flow of each channel in a reaction region after the three-cavity fluid is distributed by a distribution region;
and respectively adding three-cavity fluid inlet input boundary conditions, three-cavity pressure outlet boundary conditions, current collecting plate upper and lower boundary potential boundary conditions and heat dissipation boundary conditions to the correction model, respectively establishing periodic boundary conditions between the materials of the upper cutting section and the lower cutting section, and performing flow, heat transfer, component transmission and electrochemical multi-physical-field coupling computational analysis on the correction model.
Further, when a simplified geometric model of the full-area bidirectional corrugated flow channel fuel cell is established: firstly, placing two bipolar plate geometric models according to actual assembly positions, removing a distribution area and reserving a reaction area; respectively establishing an upper cutting section and a lower cutting section in a flow passage of a reaction area, keeping the integrity of three air flow passages of the bipolar plate, ensuring the periodicity of the upper cutting section and the lower cutting section in a full-area single cell, and respectively cutting the anode bipolar plate and the cathode bipolar plate by using the upper cutting section and the lower cutting section; performing Boolean operation on the cut bipolar plate, establishing corresponding air fluid domain, hydrogen fluid domain and cooling liquid fluid domain, establishing other layers of materials of the proton exchange membrane fuel cell between the upper cutting section and the lower cutting section according to the actual structure size of the single cell,
wherein, other layers comprise a proton exchange membrane, an anode catalysis layer, a cathode catalysis layer, an anode microporous layer, a cathode microporous layer, an anode diffusion layer, a cathode diffusion layer, an anode current collecting plate and a cathode current collecting plate.
Further, when the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell is subjected to structural mechanics analysis: the air fluid domain, the hydrogen fluid domain and the cooling liquid fluid domain are removed, the air side bipolar plate and the hydrogen side bipolar plate are arranged as rigid bodies, periodic boundary conditions are arranged on the upper cutting section and the lower cutting section, the hydrogen side bipolar plate is used as a supporting surface, and a certain assembling force is exerted on the air side bipolar plate, so that a new membrane electrode diffusion layer with deformation conditions at different positions after the two-way corrugated flow channel bipolar plate is compressed is obtained.
By adopting the technical scheme, the method for simplifying the coupling calculation of the multi-physical field of the corrugated flow channel fuel cell provided by the invention can obtain the distribution condition of each parameter of the complex full-area single cell under the real condition by establishing a flow related to the coupling calculation of the multi-physical field of structural force, fluid, temperature, performance and the like and carrying out calculation analysis on the simplified model, and can carry out physical quantity analysis on the complex full-area single cell by applying the established model, thereby achieving better balance between the calculation efficiency and the accuracy, being used for rapidly analyzing the characteristic analysis of each parameter of the complex proton exchange membrane fuel cell and providing a reliability basis for the flow field structure design of the bipolar plate.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a schematic view of a bipolar plate of a PEM fuel cell according to the present invention;
FIG. 2 is a top view of a portion of a channel in a reaction area of a bi-directional corrugated channel bipolar plate according to the present invention;
FIG. 3 is a simplified schematic view of the direction of flow of the mold of the present invention;
FIG. 4 is a simplified schematic diagram of a model cell according to the present invention in the thickness direction.
In the figure: 1. a hydrogen bridging zone; 2. a coolant bridge zone; 3. an air bridge zone; 4. a three-chamber fluid distribution region; 5. a reaction zone; 6. an air flow passage; 7. a hydrogen gas flow channel; 8. cutting the section upwards; 9. cutting a section downwards; 10. an air inlet minimum flow boundary; 11. an air inlet mean flow boundary; 12. an air inlet maximum flow boundary; 13. a hydrogen gas inlet boundary; 14. a proton exchange membrane; 15. an anode catalyst layer; 16. a cathode catalyst layer; 17. an anodic microporous layer; 18. a cathode microporous layer; 19. an anode diffusion layer; 20. a cathode diffusion layer; 21. a hydrogen gas fluid domain; 22. an air-fluid domain; 23. an anode bipolar plate; 24. cathode bipolar plate.25. anode coolant fluid field; 26. a cathode coolant fluid field; 27. an anode current collector plate; 28. and a cathode collector plate.
Detailed Description
In order to make the technical solutions and advantages of the present invention clearer, the following describes the technical solutions in the embodiments of the present invention clearly and completely with reference to the drawings in the embodiments of the present invention:
as shown in fig. 1, a bipolar plate of a pem fuel cell is schematically illustrated, wherein one side of the bipolar plate is a hydrogen fluid region, the other side of the bipolar plate is an air fluid region, and the middle forming region is a coolant fluid region. The flowing condition of the three-cavity fluid of the air cavity, the hydrogen cavity and the cooling liquid cavity of the bipolar plate is that the fluid passes through the bridge plate area, namely the hydrogen bridge plate area 1, the cooling liquid bridge plate area 2 and the air bridge plate area 3 from the inlet, is redistributed by the fluid distribution area 4, then passes through each flow channel of the reaction area 5, and flows out of the single pool by the outlet bridge plate area after converging by the outlet distribution area. Fig. 2 is a top view of a part of the channels in the reaction zone of the bi-directional corrugated channel bipolar plate, wherein the solid lines represent the air channels 6, and the dotted lines represent the hydrogen channels 7, wherein the corresponding positions of air and hydrogen are present, so that the electrochemical reaction is stronger.
Under the condition that a complete single-flow-channel model on two sides cannot be established, the invention discloses a detailed establishing process for simplifying a geometric model of a full-area bidirectional corrugated flow-channel fuel cell. After the flow of each flow channel in the reaction area is acted by the distribution area, the flow of each flow channel is inconsistent, so that the model provided by the application needs to ensure the integrity of 3 air flow channels, and the flow of the 3 air flow channels is respectively specified as the minimum flow, the average flow and the maximum flow of each flow channel in the reaction area after distribution by the distribution area so as to calculate and analyze the influence condition of the distribution area on each parameter of the battery.
Firstly, placing two bipolar plate geometric models according to actual assembly positions, removing a distribution area and reserving a reaction area, as shown in figure 3, respectively establishing an upper cutting section 8 and a lower cutting section 9 in a flow channel of the reaction area, reserving the integrity of 3 air flow channels in the model, and ensuring the periodicity of the upper cutting section and the lower cutting section in a full-area single cell.
The anode bipolar plate 23 and the cathode bipolar plate 24 are cut by the upper cut section and the lower cut section. On the basis of the cut bipolar plate and the upper and lower cut sections, an air fluid domain 22, a hydrogen fluid domain 21, an anode coolant fluid domain 25 and a cathode coolant fluid domain 26 are respectively established through Boolean operation. Then, establishing other layers of materials of the proton exchange membrane fuel cell between the upper cutting section and the lower cutting section according to the actual structure size of the single cell, wherein the other layers are as follows: a proton exchange membrane 14, an anode catalysis layer 15, a cathode catalysis layer 16, an anode microporous layer 17, a cathode microporous layer 18, an anode diffusion layer 19, a cathode diffusion layer 20, an anode current collecting plate 27, and a cathode current collecting plate 28, and a schematic longitudinal section is shown in fig. 4.
The simplified full-area bidirectional corrugated flow channel fuel cell geometric model comprises a cathode and anode cooling liquid fluid domain, wherein the cathode and anode cooling liquid fluid domain is the central region of a cathode and anode bipolar plate, the established simplified full-area bidirectional corrugated flow channel fuel cell geometric model is subjected to structural mechanics calculation, and new membrane electrode diffusion layers in deformation states at different positions after the bidirectional corrugated flow channel bipolar plate is compressed are obtained. The specific calculation flow is as follows: the air fluid domain, the hydrogen fluid domain and the cooling liquid fluid domain are removed, the air side bipolar plate and the hydrogen side bipolar plate are arranged as rigid bodies, periodic boundary conditions are arranged on the upper cutting section and the lower cutting section, the hydrogen side bipolar plate is used as a supporting surface, and a certain assembling force is exerted on the air side bipolar plate, so that a new membrane electrode diffusion layer with deformation conditions at different positions after the two-way corrugated flow channel bipolar plate is compressed is obtained. The deformation result of the membrane electrode structure, namely the new membrane electrode diffusion layer is combined with the bipolar plate, and new air fluid domains at different positions of the bipolar plate, which are influenced by the membrane electrode structure embedding, can be obtained through Boolean operation. And replacing the membrane electrode structure and the air fluid domain in the built geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell with a new membrane electrode diffusion layer and a new air fluid domain, and correcting the original model to consider the influence of structural force on each parameter of the cell, namely obtaining a corrected model.
The established correction models are respectively used for establishing corresponding non-common node interfaces at the cathode and anode diffusion layer and cathode and anode reaction gas interface, the cathode and anode reaction gas and cathode and anode cooling liquid interface, the cathode and anode cooling liquid central interface and the cathode and anode cooling liquid and cathode and anode collector plate interface, so that the model grids can be processed into structural grids, and the grid quality and the calculation precision are improved. Periodic boundary conditions between corresponding layers of materials are respectively set on the upper cutting section and the lower cutting section, then a grid with completely consistent model upper cutting section and model lower cutting section is established, and the data between the upper cutting section and the lower cutting section can be conveniently transmitted without difference.
Performing Boolean operation on the bipolar plate of the full-area single cell to respectively obtain three-cavity fluid domains of the air cavity, the hydrogen cavity and the cooling liquid cavity of the full-area single cell, then respectively performing three-cavity fluid distribution calculation to respectively obtain the fluid flow of each channel in the reaction region after the three-cavity fluid is distributed by the distribution region. The model of the present application is a counter-current reaction situation, as shown in fig. 3, the average flow value of each channel in the reaction region is respectively taken as three boundary conditions of hydrogen inlet 13 and coolant inlet 12, and the minimum flow value, the average flow value and the maximum flow value of each channel in the reaction region are respectively taken as three boundary conditions of air inlet, namely, the minimum flow boundary 10 of air inlet, the average flow boundary 11 of air inlet and the maximum flow boundary 12 of air inlet.
And respectively adding three-cavity fluid inlet input boundary conditions, three-cavity pressure outlet boundary conditions, current collecting plate upper and lower boundary potential boundary conditions and heat dissipation boundary conditions to the corrected model, respectively establishing periodic boundary conditions among materials of upper and lower cutting sections, and then loading commercial software such as a PEMFC (proton exchange membrane fuel cell) module in FLUENT (flash fuel cell) for calculation. In order to increase the calculation convergence and stability, a certain water content is given during initialization, a flow equation is started during calculation, and other equations such as a component equation, an energy equation and electrical property are started in sequence after the fluid is calculated and stabilized.
The invention relates to a simplifying method for multi-physical field coupling calculation of a corrugated flow channel fuel cell, which can obtain the distribution condition of each parameter of a complicated full-area single cell under the real condition by establishing a flow related to the multi-physical field coupling calculation of structural force, fluid, temperature, performance and the like, wherein the simplifying model can take the influence of a single cell distribution area and a real flow field structure into consideration, and can be closer to various transmission mechanisms of a proton exchange membrane fuel cell under the real condition by taking the influence of the structural force on the compression characteristic and the flow field of a membrane electrode MEA (membrane electrode assembly) into consideration in the coupling calculation of the multi-physical field.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art should be considered to be within the technical scope of the present invention, and the technical solutions and the inventive concepts thereof according to the present invention should be equivalent or changed within the scope of the present invention.
Claims (3)
1. A corrugated flow channel fuel cell multi-physical field coupling calculation simplification method is characterized by comprising the following steps:
establishing a simplified full-area bidirectional corrugated flow channel fuel cell geometric model with an upper cutting section and a lower cutting section;
carrying out structural mechanics analysis on the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell to obtain new membrane electrode diffusion layers in different deformation states at different positions after the bidirectional corrugated flow channel bipolar plate is compressed;
combining the new membrane electrode diffusion layer with the bipolar plate, and obtaining new air fluid domains at different positions of the bipolar plate due to the embedding influence of the membrane electrode diffusion layer through Boolean operation; combining the new membrane electrode diffusion layer with the new air fluid domain to replace the membrane electrode structure and the air fluid domain of the established geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell, and obtaining a correction model of the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell;
establishing corresponding non-common node interfaces of the correction model on a cathode and anode diffusion layer and cathode and anode reaction gas interface, a cathode and anode reaction gas and cathode and anode cooling liquid interface, a cathode and anode cooling liquid central interface and a cathode and anode cooling liquid and cathode and anode collector plate interface respectively; setting periodic boundary conditions between corresponding layers of materials on the upper cutting section and the lower cutting section respectively, and establishing a grid with the completely consistent upper cutting section and lower cutting section of the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell;
performing Boolean operation on the bipolar plate of the full-area single cell to respectively obtain three-cavity fluid domains of an air cavity, a hydrogen cavity and a cooling liquid cavity of the full-area single cell, and respectively performing three-cavity fluid distribution calculation to respectively obtain the fluid flow of each channel in a reaction region after the three-cavity fluid is distributed by a distribution region;
and respectively adding three-cavity fluid inlet input boundary conditions, three-cavity pressure outlet boundary conditions, current collecting plate upper and lower boundary potential boundary conditions and heat dissipation boundary conditions to the correction model, respectively establishing periodic boundary conditions between the materials of the upper cutting section and the lower cutting section, and performing flow, heat transfer, component transmission and electrochemical multi-physical-field coupling computational analysis on the correction model.
2. The corrugated flow channel fuel cell multiphysics coupling computation simplification method of claim 1, further characterized by: when a geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell is established: firstly, placing two bipolar plate geometric models according to actual assembly positions, removing a distribution area and reserving a reaction area; respectively establishing an upper cutting section and a lower cutting section in a flow passage of a reaction area, keeping the integrity of three air flow passages of the bipolar plate, ensuring the periodicity of the upper cutting section and the lower cutting section in a full-area single cell, and respectively cutting the anode bipolar plate and the cathode bipolar plate by using the upper cutting section and the lower cutting section; performing Boolean operation on the cut bipolar plate, establishing corresponding air fluid domain, hydrogen fluid domain and cooling liquid fluid domain, establishing other layers of materials of the proton exchange membrane fuel cell between the upper cutting section and the lower cutting section according to the actual structure size of the single cell,
wherein, other layers comprise a proton exchange membrane, an anode catalysis layer, a cathode catalysis layer, an anode microporous layer, a cathode microporous layer, an anode diffusion layer, a cathode diffusion layer, an anode current collecting plate and a cathode current collecting plate.
3. The corrugated flow channel fuel cell multiphysics coupling computation simplification method of claim 1, further characterized by: when the geometric model of the simplified full-area bidirectional corrugated flow channel fuel cell is subjected to structural mechanics analysis: the air fluid domain, the hydrogen fluid domain and the cooling liquid fluid domain are removed, the air side bipolar plate and the hydrogen side bipolar plate are arranged as rigid bodies, periodic boundary conditions are arranged on the upper cutting section and the lower cutting section, the hydrogen side bipolar plate is used as a supporting surface, and a certain assembling force is exerted on the air side bipolar plate, so that a new membrane electrode diffusion layer with deformation conditions at different positions after the two-way corrugated flow channel bipolar plate is compressed is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011615085.6A CN112800643B (en) | 2020-12-30 | 2020-12-30 | Multi-physical field coupling calculation simplifying method for corrugated flow passage fuel cell |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011615085.6A CN112800643B (en) | 2020-12-30 | 2020-12-30 | Multi-physical field coupling calculation simplifying method for corrugated flow passage fuel cell |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112800643A true CN112800643A (en) | 2021-05-14 |
| CN112800643B CN112800643B (en) | 2024-03-08 |
Family
ID=75804687
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011615085.6A Active CN112800643B (en) | 2020-12-30 | 2020-12-30 | Multi-physical field coupling calculation simplifying method for corrugated flow passage fuel cell |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112800643B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113782106A (en) * | 2021-07-29 | 2021-12-10 | 空气动力学国家重点实验室 | Solving method for high enthalpy dissociation gas surface catalytic reaction rate for CFD |
| CN115048818A (en) * | 2022-08-15 | 2022-09-13 | 潍柴动力股份有限公司 | Bipolar plate strength simulation model building method and device and simulation model building equipment |
| CN117448858A (en) * | 2023-10-18 | 2024-01-26 | 三一氢能有限公司 | Flow field structure and electrolytic tank |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060089803A1 (en) * | 2002-12-27 | 2006-04-27 | Riken | Method and device for numberical analysis of flow field of non-compressive viscous fluid, directly using v-cad data |
| CN101000962A (en) * | 2006-01-11 | 2007-07-18 | 董辉 | Preparation method of wedge dipole plate for fuel cell |
| JP2012237652A (en) * | 2011-05-12 | 2012-12-06 | Mitsubishi Motors Corp | Method of estimating temperature of solid |
| CN111079337A (en) * | 2019-12-23 | 2020-04-28 | 畔星科技(浙江)有限公司 | Multi-physical-field coupling simulation method for proton exchange membrane fuel cell |
| CN111090964A (en) * | 2019-12-30 | 2020-05-01 | 新源动力股份有限公司 | Full-area single-pool model building method combining fluid and performance simulation |
-
2020
- 2020-12-30 CN CN202011615085.6A patent/CN112800643B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060089803A1 (en) * | 2002-12-27 | 2006-04-27 | Riken | Method and device for numberical analysis of flow field of non-compressive viscous fluid, directly using v-cad data |
| CN101000962A (en) * | 2006-01-11 | 2007-07-18 | 董辉 | Preparation method of wedge dipole plate for fuel cell |
| JP2012237652A (en) * | 2011-05-12 | 2012-12-06 | Mitsubishi Motors Corp | Method of estimating temperature of solid |
| CN111079337A (en) * | 2019-12-23 | 2020-04-28 | 畔星科技(浙江)有限公司 | Multi-physical-field coupling simulation method for proton exchange membrane fuel cell |
| CN111090964A (en) * | 2019-12-30 | 2020-05-01 | 新源动力股份有限公司 | Full-area single-pool model building method combining fluid and performance simulation |
Non-Patent Citations (1)
| Title |
|---|
| 罗鑫: "PEM燃料电池流场传质模拟分析", 《中国硕士论文全文库，工程科技I辑》, 15 April 2019 (2019-04-15) * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113782106A (en) * | 2021-07-29 | 2021-12-10 | 空气动力学国家重点实验室 | Solving method for high enthalpy dissociation gas surface catalytic reaction rate for CFD |
| CN113782106B (en) * | 2021-07-29 | 2024-04-16 | 空气动力学国家重点实验室 | Solving method for high enthalpy dissociation gas surface catalysis reaction rate for CFD |
| CN115048818A (en) * | 2022-08-15 | 2022-09-13 | 潍柴动力股份有限公司 | Bipolar plate strength simulation model building method and device and simulation model building equipment |
| CN115048818B (en) * | 2022-08-15 | 2023-01-06 | 潍柴动力股份有限公司 | Method and device for building bipolar plate strength simulation model and simulation model building equipment |
| CN117448858A (en) * | 2023-10-18 | 2024-01-26 | 三一氢能有限公司 | Flow field structure and electrolytic tank |
| CN117448858B (en) * | 2023-10-18 | 2024-04-19 | 三一氢能有限公司 | Flow field structure and electrolytic tank |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112800643B (en) | 2024-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112800643B (en) | Multi-physical field coupling calculation simplifying method for corrugated flow passage fuel cell | |
| Aiyejina et al. | PEMFC flow channel geometry optimization: a review | |
| CN107145658B (en) | Numerical simulation method for flow field parameter design of bipolar plate of proton exchange membrane fuel cell | |
| CN111079337B (en) | Multi-physical field coupling simulation method for proton exchange membrane fuel cell | |
| Chen et al. | Physics field distributions within fuel cell stacks with manifolds penetrating through the plane zone and open outlet features | |
| US8728682B2 (en) | Fuel cell having a fastening structure including a specified configuration of elastic members | |
| CN115101787A (en) | Flow simulation method and device of fuel cell, electronic equipment and medium | |
| Peng et al. | Design, optimization, and fabrication of slotted-interdigitated thin metallic bipolar plates for PEM fuel cells | |
| CN114741877A (en) | Fuel cell multi-physical field coupling simulation analysis method under heterogeneous compression state | |
| Qiu et al. | Channel dimensional error effect of stamped bipolar plates on the characteristics of gas diffusion layer contact pressure for proton exchange membrane fuel cell stacks | |
| Huang et al. | Performance improvement in a proton exchange membrane fuel cell with an innovative flow field design | |
| Boucetta et al. | A review on the performance and modelling of proton exchange membrane fuel cells | |
| CN116845286A (en) | Optimized cooling method for three-dimensional non-isothermal proton exchange membrane fuel cell based on heavy truck | |
| EP4239511A1 (en) | Design assistance system, design assistance method, and design assistance program | |
| Barzegari et al. | Developed endplate geometry for uniform contact pressure distribution over PEMFC active area | |
| Miftah et al. | Study effect of stress in the electrical contact resistance of bipolar plate and membrane electrode assembly in proton exchange membrane fuel cell: A review | |
| CN112666462A (en) | Method and device for evaluating PEMFC performance | |
| CN116207304B (en) | Modeling method of membrane humidifier for fuel cell, membrane humidifier module and model | |
| Ansari | Reduced-order modeling of PEM fuel cell based on POD and PODI: an efficient approach toward combining highest accuracy with real-time performance | |
| Kalaiyarasan et al. | Design and Analysis of Bipolar Plate for 500w Pem Fuel Cell Stack | |
| Chen et al. | Error of single-phase proton exchange membrane fuel cell model based on Brinkman-Darcy's law in different flow fields | |
| CN114117855A (en) | Fuel cell bipolar plate compression resistance analysis method based on finite element technology | |
| Robalinho et al. | COMSOL Modelling and Simulation of PEM Fuel Cell’s Flow Channels | |
| Gass et al. | An advanced 1D physics-based model for PEM hydrogen fuel cells with enhanced overvoltage prediction | |
| Azzam et al. | Design and Comparison Analysis of Various Flow Configurations in Bipolar Plates via Numerical Simulation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |