CN117668423A - Method, system, equipment and storage medium for constructing air-cooled fuel cell model - Google Patents
Method, system, equipment and storage medium for constructing air-cooled fuel cell model Download PDFInfo
- Publication number
- CN117668423A CN117668423A CN202410137521.5A CN202410137521A CN117668423A CN 117668423 A CN117668423 A CN 117668423A CN 202410137521 A CN202410137521 A CN 202410137521A CN 117668423 A CN117668423 A CN 117668423A
- Authority
- CN
- China
- Prior art keywords
- air
- equation
- fuel cell
- cooled fuel
- 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 85
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 230000008878 coupling Effects 0.000 claims abstract description 17
- 238000010168 coupling process Methods 0.000 claims abstract description 17
- 238000005859 coupling reaction Methods 0.000 claims abstract description 17
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 238000007664 blowing Methods 0.000 claims abstract description 8
- 239000007789 gas Substances 0.000 claims description 71
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 45
- 239000007788 liquid Substances 0.000 claims description 30
- 230000003197 catalytic effect Effects 0.000 claims description 17
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 15
- 239000001301 oxygen Substances 0.000 claims description 15
- 229910052760 oxygen Inorganic materials 0.000 claims description 15
- 238000003487 electrochemical reaction Methods 0.000 claims description 13
- 239000001257 hydrogen Substances 0.000 claims description 13
- 229910052739 hydrogen Inorganic materials 0.000 claims description 13
- 239000012528 membrane Substances 0.000 claims description 12
- 238000009792 diffusion process Methods 0.000 claims description 11
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 7
- 238000004134 energy conservation Methods 0.000 claims description 6
- 229920000554 ionomer Polymers 0.000 claims description 6
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 6
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 claims description 6
- 239000007787 solid Substances 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 239000003792 electrolyte Substances 0.000 claims description 3
- 150000002431 hydrogen Chemical class 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 229910052697 platinum Inorganic materials 0.000 claims description 3
- 239000011521 glass Substances 0.000 claims 1
- 229920000642 polymer Polymers 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 6
- 238000010276 construction Methods 0.000 description 4
- 238000012546 transfer Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04305—Modeling, demonstration models of fuel cells, e.g. for training purposes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04007—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids related to heat exchange
- H01M8/04014—Heat exchange using gaseous fluids; Heat exchange by combustion of reactants
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M8/00—Fuel cells; Manufacture thereof
- H01M8/04—Auxiliary arrangements, e.g. for control of pressure or for circulation of fluids
- H01M8/04298—Processes for controlling fuel cells or fuel cell systems
- H01M8/04992—Processes for controlling fuel cells or fuel cell systems characterised by the implementation of mathematical or computational algorithms, e.g. feedback control loops, fuzzy logic, neural networks or artificial intelligence
-
- 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/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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Mathematical Optimization (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Electrochemistry (AREA)
- Manufacturing & Machinery (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Chemical & Material Sciences (AREA)
- Sustainable Energy (AREA)
- Sustainable Development (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Algebra (AREA)
- Combustion & Propulsion (AREA)
- Bioinformatics & Computational Biology (AREA)
- Analytical Chemistry (AREA)
- Operations Research (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Databases & Information Systems (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Automation & Control Theory (AREA)
- Fuzzy Systems (AREA)
- Medical Informatics (AREA)
- Fuel Cell (AREA)
Abstract
The invention discloses a method, a system, equipment and a storage medium for constructing an air-cooled fuel cell model, and relates to the technical field of fuel cells. Comprising the following steps: establishing a calculation domain according to the real three-dimensional geometric structure of the fan and the air-cooled fuel cell stack, and dividing a fan blade area in the calculation domain into a rotation area; setting two different air supply modes of air suction and air blowing simulated by different rotation directions of a rotation area, and simulating turbulent flow of the rotation area under the different air supply modes by using a RANS model to obtain a control equation set of turbulent flow; coupling a control equation set of turbulent flow with a gas component mass fraction equation of the rotating area to obtain a control equation set of turbulent flow with multiple gas components; and the control equation set is coupled with a multi-physical field transmission control equation of the air-cooled fuel cell stack to obtain a model equation set of the air-cooled fuel cell. The method can improve the accuracy of the calculation result and provide accurate theoretical guidance for structural design optimization of the air-cooled fuel cell.
Description
Technical Field
The present invention relates to the field of fuel cell technologies, and in particular, to a method, a system, an apparatus, and a storage medium for constructing an air-cooled fuel cell model.
Background
The air-cooled fuel cell has the advantages of simple system, small volume and the like, and has great prospect in the fields of small unmanned aerial vehicles and the like. An air-cooled fuel cell is a highly complex and versatile system involving multiple physical field transport and electrochemical reaction coupling. With the continuous rising of the cost of the air-cooled fuel cell simulated by the real materials, the air-cooled fuel cell is simulated by a properly simplified model, so that the cost can be reduced while certain accuracy is maintained, and the deep understanding and technical progress of the air-cooled fuel cell are promoted.
In the prior art, the simulation means for the air-cooled fuel cell mostly only focuses on the fuel cell, firstly, describes the electrochemical reaction in the fuel cell through an electrode dynamics model, describes the heat and mass transfer process in the fuel cell through thermodynamic and fluid dynamics equations, and finally combines the electrochemical reaction and the heat and mass transfer process to establish a complete control equation set. And the control equation set is established for the electrochemical reaction process and the multi-physical-field thermal mass transmission process in the fuel cell to carry out numerical solution, so as to obtain a model equation of the air-cooled fuel cell.
The defects of the prior art are as follows: in the process of obtaining the air-cooled fuel cell model equation, the influence of a flow field generated by rotation of a fan on the hydrothermal state in the fuel cell is not considered, so that a calculation result has obvious deviation from an actual situation.
Disclosure of Invention
In view of the foregoing, it is desirable to provide a method, a system, an apparatus, and a storage medium for constructing an air-cooled fuel cell model.
The embodiment of the invention provides a method for constructing an air-cooled fuel cell model, which comprises the following steps:
establishing a calculation domain according to the real three-dimensional geometric structure of the fan and the air-cooled fuel cell stack, and dividing a fan blade area in the calculation domain into a rotation area;
setting two different air supply modes of air suction and air blowing simulated by different rotation directions of a rotation area, and simulating turbulent flow of the rotation area under the different air supply modes by using a RANS model to obtain a control equation set of the turbulent flow; coupling a control equation set of turbulent flow with a gas component mass fraction equation of the rotating area to obtain a control equation set of turbulent flow with multiple gas components;
and coupling a multi-physical field transmission control equation of the air-cooled fuel cell stack through a control equation set of turbulent flow of multiple gas components to obtain a model equation set of the air-cooled fuel cell.
Additionally, the computational domain further includes a fixed region comprising: an air inlet region, a fan grill region, an air cavity region, an air cooled fuel cell stack region, and an air outlet region; the air inlet area absorbs external air, sequentially passes through the rotating area, the fan grid area, the air cavity area and the air-cooled fuel cell stack area, and is discharged through the air outlet area.
In addition, the mass fraction equation of the gas component is as follows:
wherein,for time (I)>For the density of the mixed gas>For liquid water saturation, +.>For the speed of the mixed gas>For porosity->Is hydrogen, oxygen, water vapor mass fraction, +.>For effective gas diffusion coefficient, +.>Is a source term of the gas composition equation.
Additionally, the multi-gas component turbulent flow control equation set includes:
mass conservation equation of the mixed gas:
mixed gas momentum conservation equation:
turbulence energy equation:
turbulent energy dissipation ratio equation:
wherein,for time (I)>For the density of the mixed gas>For liquid water saturation, +.>For the speed of the mixed gas>For porosity->For the source term of the conservation equation of mass of the mixed gas, < ->For gas pressure +.>Is aerodynamic viscosity>For the source term of the conservation of momentum equation, +.>For turbulent motion, add>For turbulent viscosity>Is a turbulent energy item->And->Is two in numberEmpirical coefficient (F)>Being the plancut number of the turbulent energy dissipation ratio.
In addition, the model equation set of the air-cooled fuel cell includes:
the hydraulic equation:
wherein,for time (I)>Is of liquid water density->For porosity->For liquid water saturation, +.>Is of intrinsic permeability->For relative permeability, ++>Is hydraulic and is->Is aerodynamic viscosity>Is a hydraulic equation source term; the relationship between capillary pressure and liquid water saturation is constructed using the Leverett-J equation, where capillary pressure is the difference between air pressure and hydraulic pressure:
wherein,for capillary pressure>And->Is pneumatic and hydraulic respectively,)>Is the surface tension coefficient +.>For contact angle, ++>Is of intrinsic permeability->Is liquid water saturation;
membrane state water equation:
wherein,is dry film density, < >>For time (I)>For the water content in the form of a film,/->For electrolyte volume fraction,/a->Is electroosmotic drag coefficient->Is ion current density, +.>For Faraday constant, EW is film equivalent mass, ">Is the effective membrane state water diffusion coefficient +.>Is a source term of a membrane state water equation;
electron potential equation:
ion potential equation:
in the method, in the process of the invention,and->Respectively electron potential and ion potential, +.>And->Effective electron conductivity and effective ion conductivity, respectively,>and->Respectively an electron potential equation and an ion potential equation source term;
energy conservation equation:
wherein,for time (I)>、/>And->Mixed gas, liquid water and solid density, respectively, < >>For liquid water saturation, +.>And->The speeds of the mixed gas and the liquid water are respectively +.>For porosity->Is specific heat capacity->For temperature, < >>For effective heat conductivity, +.>The subscripts g, l, and s are mixed gas, liquid water, and solid, respectively, for the energy conservation equation source term.
In addition, the electron potential equation source termThe method comprises the following steps:
the ion potential equation source termThe method comprises the following steps:
wherein,for the anodic catalytic layer electrochemical reaction rate, +.>Electrochemical reaction rate of the cathode catalytic layer;
electrochemical reaction rate of anode catalytic layerAnd the electrochemical reaction rate of the cathode catalytic layer->The calculation formula of (2) is as follows:
wherein a and c are an anode and a cathode respectively,in order to exchange the current density of the current,/>is a temperature correction coefficient>For effective reaction area per unit mass of platinum, < >>And->Hydrogen and oxygen partial pressures, respectively,/->And->Hydrogen and oxygen permeation henry coefficients, respectively, +.>And->Reference concentrations of hydrogen and oxygen, respectively, +.>Is a universal gas constant>For temperature, < >>For exchange coefficients +.>To activate the overpotential +.>For caking effective coefficient, +.>For the porosity of the catalytic layer->For the caking radius>To agglomerate the outer ionomer film thickness +.>For the oxygen diffusion coefficient in ionomers, < >>The specific surface area of the agglomerate.
In addition, an air-cooled fuel cell model building system includes:
the calculation domain establishing module is used for establishing a calculation domain according to the real three-dimensional geometric structure of the fan and the air-cooled fuel cell stack and dividing a fan blade area in the calculation domain into a rotation area;
the simulation module is used for setting two different air supply modes of air suction and air blowing of different rotation directions of the rotation area, simulating turbulent flow of the rotation area under the different air supply modes through the RANS model, and obtaining a control equation set of the turbulent flow; coupling a control equation set of turbulent flow with a gas component mass fraction equation of the rotating area to obtain a control equation set of turbulent flow with multiple gas components;
and the coupling module is used for coupling a control equation set of turbulent flow of multiple gas components with a multiple physical field transmission control equation of the air-cooled fuel cell stack to obtain a model equation set of the air-cooled fuel cell.
In addition, the computer equipment comprises a memory and a processor, wherein the memory stores a computer program, and the processor realizes the steps of the air-cooled fuel cell model building method when executing the computer program.
In addition, a storage medium has stored thereon a computer program which, when executed by a processor, implements the steps of the air-cooled fuel cell model building method described above.
Compared with the prior art, the method, the system, the equipment and the storage medium for constructing the air-cooled fuel cell model have the following beneficial effects:
according to the invention, a calculation domain is established according to the real three-dimensional geometric structure of the fan and the air-cooled fuel cell stack, and a fan blade area in the calculation domain is divided into a rotation area; setting two different air supply modes of air suction and air blowing simulated by different rotation directions of a rotation area, and simulating turbulent flow of the rotation area under the different air supply modes by using a RANS model to obtain a control equation set of the turbulent flow; coupling a control equation set of turbulent flow with a gas component mass fraction equation of the rotating area to obtain a control equation set of turbulent flow with multiple gas components; and the control equation set is coupled with a multi-physical field transmission control equation of the air-cooled fuel cell stack to obtain a model equation set of the air-cooled fuel cell.
Compared with the prior art, the technical scheme combines the control equation set obtained by simulating turbulent flow generated by rotation of the rotation region with the multi-physical-field transmission control equation set obtained by the gas-water-heat-electricity transmission process in the internal working process of the fuel cell to construct the model equation set of the air-cooled fuel cell, so that the accuracy of a calculation result can be improved, and accurate theoretical guidance is provided for structural design optimization of the air-cooled fuel cell.
Drawings
FIG. 1 is a flow chart of a method of modeling an air-cooled fuel cell in accordance with one embodiment;
FIG. 2 is a graph of fan static pressure, stack and air cavity pressure drop for different fan-to-stack distances for an air-cooled fuel cell model construction method provided in one embodiment;
FIG. 3 is a graph of air flow rate into each single cell of a stack for different fan-to-stack distances for an air-cooled fuel cell model construction method provided in one embodiment;
FIG. 4 is a graph showing the average temperature of each single-cell catalytic layer in a stack at different fan-stack distances for an air-cooled fuel cell model building method according to one embodiment;
fig. 5 is a graph of operating voltages of each single cell in a stack for different fan-stack distances in an air-cooled fuel cell model construction method according to an embodiment.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
In one embodiment, a method for constructing an air-cooled fuel cell model is provided, as shown in fig. 1, and the method includes:
the calculation domain comprises an air inlet area, a fan blade rotation area, a fan grid area, an air cavity area, an air-cooled fuel cell pile area and an air outlet area, wherein the pile is formed by connecting a plurality of single cells in series, and each cell comprises a polar plate/flow field, a gas diffusion layer, a microporous layer, a catalytic layer, a proton exchange membrane and other components. The calculation domain is divided into a rotation region and a fixed region by adopting a multiple reference system model, wherein the region where the fan blades are located is the rotation region, the central axis of the rotation region is set as a rotation axis, the actual rotation speed of the fan is the rotation speed of the rotation region, and the simulation of two different air supply modes of air suction and air blowing of the fan can be realized by setting different rotation directions, namely clockwise rotation and anticlockwise rotation. The fixed area includes: an air inlet region, a fan grill region, an air cavity region, an air cooled fuel cell stack region, and an air outlet region; the air inlet area absorbs external air, sequentially passes through the rotating area, the fan grid area, the air cavity area and the air-cooled fuel cell stack area, and is discharged through the air outlet area.
The turbulence flow of the rotating area is simulated by adopting the RANS model, and the unsteady flow field generated by the rotation of the fan blade can be converted into the unsteady flow field by adopting the simulation calculation method, so that the method has good calculation precision and calculation efficiency. Meanwhile, the RANS model equation is coupled with a gas component mass fraction equation to obtain a control equation set of turbulent flow of multiple gas components, wherein the control equation set is as follows:
mass conservation equation of the mixed gas:
mixed gas momentum conservation equation:
turbulence energy equation:
turbulent energy dissipation ratio equation:
the mass fraction equation of the gas component is as follows:
wherein,for time (I)>Representing the density of the mixed gas>Indicating liquid water saturation, ++>Indicating the speed of the mixed gas, +.>Indicating porosity->Representing the source term of the conservation equation of mass of the mixed gas, < +.>Indicating gas pressure, +.>Represents aerodynamic viscosity, +.>Representing the source term of the conservation of momentum equation, +.>Represents the mass fraction of hydrogen, oxygen and water vapor, < >>Indicating the effective gas diffusion coefficient,/-, for>Representing the source term of the gas composition equation, +.>Representing turbulence energy->Represents turbulent viscosity->Represents a turbulent energy item, < ->And->Representing two empirical coefficients>Prandtl number, which represents the turbulent energy dissipation rate.
Based on a control equation of turbulent flow of multiple gas components, coupling a multi-physical field transmission control equation of a fuel cell stack to obtain a fuel cell model, wherein the specific equation is as follows:
the hydraulic equation:
wherein,for time (I)>Is of liquid water density->For porosity->For liquid water saturation, +.>Is of intrinsic permeability->For relative permeability, ++>Is hydraulic and is->Is aerodynamic viscosity>Is a hydraulic equation source term. And constructing a relation between capillary pressure and liquid water saturation by using a Leverett-J equation, wherein the capillary pressure is a difference between air pressure and liquid pressure and is as follows:
wherein,for capillary pressure>And->Is pneumatic and hydraulic respectively,)>Is the surface tension coefficient +.>For contact angle, ++>Is of intrinsic permeability->Is liquid water saturation.
Membrane state water equation:
wherein,is dry film density, < >>For time (I)>For the water content in the form of a film,/->For electrolyte volume fraction,/a->Is electroosmotic drag coefficient->Is ion current density, +.>For Faraday constant, EW is film equivalent mass, ">Is to haveEffective membrane state water diffusion coefficient->Is a source term of a membrane state water equation.
Electron potential equation:
ion potential equation:
wherein,and->Respectively representing electron potential and ion potential, +.>And->Representing the effective electron conductivity and the effective ion conductivity, respectively,/-respectively>And->The electron potential equation and the ion potential equation source terms are represented, respectively.
Wherein,and->The electrochemical reaction rates of the anode and cathode catalytic layers are calculated by a Butler-Volmer equation.
Wherein a and c are an anode and a cathode respectively,representing the exchange current density, ">Representing the temperature correction coefficient, ">Represents the effective reaction area of platinum per unit mass, < >>And->Respectively represent hydrogen and oxygen partial pressures,/->And->Respectively represent the hydrogen and oxygen permeation Henry coefficients,>and->Respectively representing hydrogen and oxygen reference concentrations, +.>Representing the general gas constant, +.>Indicate temperature,/->Representing the exchange coefficient>Represents an activation overpotential, +.>Indicating caking effective coefficient, +.>Representing the porosity of the catalytic layer>Indicates the caking radius>Represents the thickness of the ionomer film outside the agglomerate,/->Represents the oxygen diffusion coefficient in the ionomer,>indicating the specific surface area of the agglomerate.
Energy conservation equation:
wherein,for time (I)>、/>And->Mixed gas, liquid water and solid density, respectively, < >>For liquid water saturation, +.>And->The speeds of the mixed gas and the liquid water are respectively +.>For porosity->Is specific heat capacity->For temperature, < >>For effective heat conductivity, +.>The subscripts g, l, and s are mixed gas, liquid water, and solid, respectively, for the energy conservation equation source term.
All the equations are three-dimensional full-size models of the air-cooled fuel cells of the coupled fans, and the distribution and performance of the hydrothermal state in the air-cooled fuel cell stacks can be obtained through numerical solution calculation. The boundary conditions corresponding to the control equations are set as follows: the boundary conditions of the air inlet and the air outlet in the calculation domain are set to be constant pressure, and the pressure, the temperature and the relative humidity of the air inlet and the air outlet are the same as the external environment of the air-cooled galvanic pile. The outer boundary of the calculation domain is set as an air forced convection heat transfer boundary condition. In the fuel cell stack, the boundary condition at the end face of the anode plate of each single cell is set as the operating current density of the stack, and the boundary condition at the end face of the cathode plate is set as the reference potential. The anode side gas flow rate was calculated using the following formula:
the mass fraction of each gas component at the inlets of the cathode and the anode is calculated by the following formula:
wherein,represents anode intake mass flow,/>Indicating the fuel cell operating current density,/->Represents the anode intake stoichiometry,/->Representing membrane electrode area, +.>Indicating the number of single cells in the stack, +.>Represents the anode inlet hydrogen concentration,/-, and>represents anode inlet pressure, +.>Indicating the relative humidity of the anode intake air,/-)>Represents the water saturation vapor pressure, < >>Is Faraday constant, +.>Representing the general gas constant, +.>Representing the molar mass of the gas component.
In one embodiment, an air-cooled fuel cell model building system is provided, the system comprising:
the calculation domain establishing module is used for establishing a calculation domain according to the real three-dimensional geometric structure of the fan and the air-cooled fuel cell stack and dividing a fan blade area in the calculation domain into a rotation area;
the simulation module is used for setting two different air supply modes of air suction and air blowing of different rotation directions of the rotation area, simulating turbulent flow of the rotation area under the different air supply modes through the RANS model, and obtaining a control equation set of the turbulent flow; coupling a control equation set of turbulent flow with a gas component mass fraction equation of the rotating area to obtain a control equation set of turbulent flow with multiple gas components;
and the coupling module is used for coupling a control equation set of turbulent flow of multiple gas components with a multiple physical field transmission control equation of the air-cooled fuel cell stack to obtain a model equation set of the air-cooled fuel cell.
A computer device comprises a memory and a processor, wherein the memory stores a computer program, and the processor realizes the steps of the air-cooled fuel cell model construction method when executing the computer program.
A storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the above-described air-cooled fuel cell model building method.
Example 1
The fan is a direct current fan of the model 9GA0624P1J03 of San Ace, the weight of the fan is 160 g, the diameter of the fan is 56 mm, the diameter of the hub is 35 mm, the fan comprises 5 blades, and the thickness of the fan is 38 mm. The rated voltage of the fan is 24V, pulse width modulation is adopted for speed regulation, and when the duty ratio of the fan is 100%, the rotating speed of the fan is 17500 r min -1 The fan power consumption was 18W. The fuel cell stack consisted of 30 single cells stacked in series, each cell having an area of 50 mm ×200 mm, a cathode and anode plate thickness excluding the flow field area of 2 mm, an anode flow field depth of 0.4 mm, and a cathode flow field depth of 1.5 mm. The gas diffusion layer, microporous layer, anode catalytic layer, cathode catalytic layer, and proton exchange membrane thicknesses were 0.1 mm, 0.03 mm, 0.005 mm, 0.01 mm, and 0.01 mm, respectively. The fan is 25 or 50 cm from the stack.
The fuel cell operating current density was set to 0.4A cm -2 The ambient pressure is 1 atm, the ambient temperature is 25 ℃, the ambient humidity is 40%, the anode inlet pressure is 1 atm, the anode inlet humidity is 40%, the anode inlet temperature is 25 ℃, and the anode inlet stoichiometric ratio is 5.0.
Based on the parameters, solving the three-dimensional full-size model of the air-cooled fuel cell stack of the coupled fan. Fig. 2 shows the fan inlet static pressure, the air-cooled stack pressure drop and the air cavity pressure drop calculated by the simulation method at the fan and stack distances of 25 and 50 cm, fig. 3 shows the wind speed distribution entering each single cell at the fan and stack distances of 25 and 50 cm, fig. 4 shows the average temperature of the cathode catalytic layer of each cell in the stack at the fan and stack distances of 25 and 50 cm, and fig. 5 shows the output voltage of each cell in the stack at the fan and stack distances of 25 and 50 cm.
The above examples illustrate only a few embodiments of the invention, which are described in detail and are not to be construed as limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention. Accordingly, the scope of protection of the present invention is to be determined by the appended claims.
Claims (9)
1. The method for constructing the air-cooled fuel cell model is characterized by comprising the following steps of:
establishing a calculation domain according to the real three-dimensional geometric structure of the fan and the air-cooled fuel cell stack, and dividing a fan blade area in the calculation domain into a rotation area;
setting two different air supply modes of air suction and air blowing simulated by different rotation directions of a rotation area, and simulating turbulent flow of the rotation area under the different air supply modes by using a RANS model to obtain a control equation set of the turbulent flow; coupling a control equation set of turbulent flow with a gas component mass fraction equation of the rotating area to obtain a control equation set of turbulent flow with multiple gas components;
and coupling a multi-physical field transmission control equation of the air-cooled fuel cell stack through a control equation set of turbulent flow of multiple gas components to obtain a model equation set of the air-cooled fuel cell.
2. The method of claim 1, wherein the calculation domain further comprises a fixed area comprising: an air inlet region, a fan grill region, an air cavity region, an air cooled fuel cell stack region, and an air outlet region; the air inlet area absorbs external air, sequentially passes through the rotating area, the fan grid area, the air cavity area and the air-cooled fuel cell stack area, and is discharged through the air outlet area.
3. The method for constructing an air-cooled fuel cell model according to claim 1, wherein the mass fraction equation of the gas component is:
wherein,for time (I)>For the density of the mixed gas>For liquid water saturation, +.>For the speed of the mixed gas>For porosity->Is hydrogen, oxygen, water vapor mass fraction, +.>For effective gas diffusion coefficient, +.>Is a source term of the gas composition equation.
4. The method of constructing an air-cooled fuel cell model of claim 1, wherein the system of multi-gas component turbulent flow control equations comprises:
mass conservation equation of the mixed gas:
mixed gas momentum conservation equation:
turbulence energy equation:
turbulent energy dissipation ratio equation:
wherein,for time (I)>For the density of the mixed gas>For liquid water saturation, +.>For the speed of the mixed gas>For porosity->For the source term of the conservation equation of mass of the mixed gas, < ->For gas pressure +.>Is aerodynamic viscosity>For the source term of the conservation of momentum equation, +.>For turbulent motion, add>For turbulent viscosity>Is a turbulent energy item->And->For two empirical coefficients +.>Being the plancut number of the turbulent energy dissipation ratio.
5. The method for constructing an air-cooled fuel cell model according to claim 1, wherein the model equation set of the air-cooled fuel cell comprises:
the hydraulic equation:
wherein,for time (I)>Is liquid water density，/>For porosity->For liquid water saturation, +.>For the purpose of the inherent permeability, the polymer,for relative permeability, ++>Is hydraulic and is->Is aerodynamic viscosity>Is a hydraulic equation source term; the relationship between capillary pressure and liquid water saturation is constructed using the Leverett-J equation, where capillary pressure is the difference between air pressure and hydraulic pressure:
wherein,for capillary pressure>And->Is pneumatic and hydraulic respectively,)>Is a surface tension systemCount (n)/(l)>As a contact angle of the glass,is of intrinsic permeability->Is liquid water saturation;
membrane state water equation:
wherein,is dry film density, < >>For time (I)>For the water content in the form of a film,/->For electrolyte volume fraction,/a->Is electroosmotic drag coefficient->Is ion current density, +.>For Faraday constant, EW is film equivalent mass, ">Is the effective membrane state water diffusion coefficient +.>Is a source term of a membrane state water equation;
electron potential equation:
ion potential equation:
wherein,and->Respectively electron potential and ion potential, +.>And->Effective electron conductivity and effective ion conductivity, respectively,>and->Respectively an electron potential equation and an ion potential equation source term;
energy conservation equation:
wherein,for time (I)>、/>And->Mixed gas, liquid water and solid density, respectively, < >>In order to achieve a liquid water saturation,and->The speeds of the mixed gas and the liquid water are respectively +.>For porosity->Is specific heat capacity->For temperature, < >>For effective heat conductivity, +.>The subscripts g, l, and s are mixed gas, liquid water, and solid, respectively, for the energy conservation equation source term.
6. The method for constructing an air-cooled fuel cell model according to claim 5, wherein,
the electron potential equation sourceItemsThe method comprises the following steps:
the ion potential equation source termThe method comprises the following steps:
wherein,for the anodic catalytic layer electrochemical reaction rate, +.>Electrochemical reaction rate of the cathode catalytic layer;
electrochemical reaction rate of anode catalytic layerAnd the electrochemical reaction rate of the cathode catalytic layer->The calculation formula of (2) is as follows:
wherein a and c are an anode and a cathode respectively,to exchange current density, +.>Is a temperature correction coefficient>For effective reaction area per unit mass of platinum, < >>And->Hydrogen and oxygen partial pressures, respectively,/->And->Hydrogen and oxygen permeation henry coefficients, respectively, +.>And->Reference concentrations of hydrogen and oxygen, respectively, +.>Is a universal gas constant>For temperature, < >>For exchange coefficients +.>To activate the overpotential +.>For caking effective coefficient, +.>For the porosity of the catalytic layer->For the caking radius>To agglomerate the outer ionomer film thickness +.>For the oxygen diffusion coefficient in ionomers, < >>The specific surface area of the agglomerate.
7. An air-cooled fuel cell model building system, comprising:
the calculation domain establishing module is used for establishing a calculation domain according to the real three-dimensional geometric structure of the fan and the air-cooled fuel cell stack and dividing a fan blade area in the calculation domain into a rotation area;
the simulation module is used for setting two different air supply modes of air suction and air blowing of different rotation directions of the rotation area, simulating turbulent flow of the rotation area under the different air supply modes through the RANS model, and obtaining a control equation set of the turbulent flow; coupling a control equation set of turbulent flow with a gas component mass fraction equation of the rotating area to obtain a control equation set of turbulent flow with multiple gas components;
and the coupling module is used for coupling a control equation set of turbulent flow of multiple gas components with a multiple physical field transmission control equation of the air-cooled fuel cell stack to obtain a model equation set of the air-cooled fuel cell.
8. A computer device comprising a memory and a processor, the memory having stored therein a computer program, characterized in that the processor, when executing the computer program, implements the steps of the air-cooled fuel cell model building method of any one of claims 1-6.
9. A storage medium having stored thereon a computer program, wherein the computer program, when executed by a processor, implements the steps of the air-cooled fuel cell model building method of any one of claims 1-6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410137521.5A CN117668423B (en) | 2024-02-01 | 2024-02-01 | Method, system, equipment and storage medium for constructing air-cooled fuel cell model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410137521.5A CN117668423B (en) | 2024-02-01 | 2024-02-01 | Method, system, equipment and storage medium for constructing air-cooled fuel cell model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117668423A true CN117668423A (en) | 2024-03-08 |
| CN117668423B CN117668423B (en) | 2024-05-03 |
Family
ID=90081111
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410137521.5A Active CN117668423B (en) | 2024-02-01 | 2024-02-01 | Method, system, equipment and storage medium for constructing air-cooled fuel cell model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117668423B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118263467A (en) * | 2024-04-09 | 2024-06-28 | 扬州大学 | Fuel cell with periodic structure and band gap regulating and controlling method thereof |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100330439A1 (en) * | 2009-06-26 | 2010-12-30 | Fischel Halbert | Electrochemical cells utilizing Taylor Vortex Flows |
| US20140255812A1 (en) * | 2013-03-08 | 2014-09-11 | Global Energy Science, Llc (A California Limited Liability Company) | Taylor vortex flow electrochemical cells utilizing particulate electrolyte suspensions |
| CN108763737A (en) * | 2018-05-25 | 2018-11-06 | 天津大学 | Proton Exchange Membrane Fuel Cells full battery D beam element method for establishing model |
| CN113297756A (en) * | 2021-07-27 | 2021-08-24 | 中国汽车技术研究中心有限公司 | Fuel cell transient modeling method adopting explicit format updating algorithm |
| CN114937787A (en) * | 2022-04-28 | 2022-08-23 | 上海交通大学内蒙古研究院 | Air-cooled fuel cell test system |
| CN115310376A (en) * | 2022-06-29 | 2022-11-08 | 山东交通学院 | Fuel cell air compressor parameter setting method based on load excitation amplitude analysis |
-
2024
- 2024-02-01 CN CN202410137521.5A patent/CN117668423B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100330439A1 (en) * | 2009-06-26 | 2010-12-30 | Fischel Halbert | Electrochemical cells utilizing Taylor Vortex Flows |
| US20140255812A1 (en) * | 2013-03-08 | 2014-09-11 | Global Energy Science, Llc (A California Limited Liability Company) | Taylor vortex flow electrochemical cells utilizing particulate electrolyte suspensions |
| CN108763737A (en) * | 2018-05-25 | 2018-11-06 | 天津大学 | Proton Exchange Membrane Fuel Cells full battery D beam element method for establishing model |
| CN113297756A (en) * | 2021-07-27 | 2021-08-24 | 中国汽车技术研究中心有限公司 | Fuel cell transient modeling method adopting explicit format updating algorithm |
| CN114937787A (en) * | 2022-04-28 | 2022-08-23 | 上海交通大学内蒙古研究院 | Air-cooled fuel cell test system |
| CN115310376A (en) * | 2022-06-29 | 2022-11-08 | 山东交通学院 | Fuel cell air compressor parameter setting method based on load excitation amplitude analysis |
Non-Patent Citations (2)
| Title |
|---|
| SWARN S. KALSI等: ""Fuel Cell Plant Concept with Dynamic Characteristics of a Rotating Generator"", 《2006 IEEE PES POWER SYSTEMS CONFERENCE AND EXPOSITION》, 5 February 2007 (2007-02-05), pages 1 - 4 * |
| 伊亭等: ""联合循环中高温燃料电池建模综述"", 《华东电力》, 30 September 2001 (2001-09-30), pages 17 - 20 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118263467A (en) * | 2024-04-09 | 2024-06-28 | 扬州大学 | Fuel cell with periodic structure and band gap regulating and controlling method thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117668423B (en) | 2024-05-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kurnia et al. | Progress on open cathode proton exchange membrane fuel cell: Performance, designs, challenges and future directions | |
| CN117668423B (en) | Method, system, equipment and storage medium for constructing air-cooled fuel cell model | |
| CN112909303B (en) | Method for establishing transient real-time model of proton exchange membrane fuel cell | |
| CN106960970B (en) | Proton Exchange Membrane Fuel Cells internal water distribution forecasting method | |
| CN107908209B (en) | Fuel cell test platform and temperature and humidity control system thereof | |
| CN101098021B (en) | Membrane humidifier for a fuel cell | |
| CN113297756B (en) | Fuel cell transient modeling method adopting explicit format updating algorithm | |
| CN115050985B (en) | Leaf vein runner bionic air-cooled fuel cell with lung air inlet mode and method thereof | |
| CN109815644A (en) | Fuel battery engines gas handling system control method based on map and Linear-Quadratic Problem | |
| Bhatt et al. | Polymer exchange membrane (PEM) fuel cell: a review | |
| CN101210750A (en) | Method for driving air-conditioner by utilizing fuel battery waste heat | |
| CN115312815B (en) | Electrochemical performance calculation method for air-cooled proton exchange membrane fuel cell stack | |
| Tong et al. | Research on influencing factors of heat transfer enhancement fins in fuel cell cooling channel | |
| CN114744254A (en) | Modeling method of hydrogen circulating pump in fuel cell system | |
| CN110203406B (en) | Unmanned aerial vehicle and portable power equipment | |
| CN102956908B (en) | Humidifier for proton exchange membrane fuel cell | |
| CN108682885B (en) | Design method for geometric dimension of porous electrode of microfluid fuel cell | |
| CN204011605U (en) | A kind of fuel cell cold start-up heater | |
| CN115394368A (en) | Battery performance prediction method based on catalytic layer carbon corrosion | |
| Park et al. | Analysis of a shell-and-tube type gas-to-gas membrane humidifier for an automotive polymer electrolyte membrane fuel cell power system | |
| Fujinuma et al. | Basic Study for Model Construction of The Water Recovery System in Polymer Electrolyte Fuel Cells | |
| Prajapati et al. | Novel custom-made design and analysis of hydrogen PEM fuel cell for high endurance small UAV applications | |
| CN115513489B (en) | Fuel cell heat dissipation subsystem simulation method | |
| CN116207304B (en) | Modeling method of membrane humidifier for fuel cell, membrane humidifier module and model | |
| CN208028156U (en) | Heat-insulated and heat-exchange device for unmanned plane fuel cell under low temperature environment |
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 |