Alkaline fuel cell model establishing method and system considering gas phase and liquid phase

Abstract

The invention provides a method and a system for establishing an alkaline fuel cell model considering gas and liquid phases. The method comprises the following steps: dividing an alkaline fuel cell into nine computational domains, namely a cathode / anode plate, a cathode / anode runner, a cathode / anode diffusion layer, a cathode / anode catalyst layer and a diaphragm; respectively establishing an electrochemical model for each computational domain of the alkaline fuel cell, wherein the electrochemical model comprises an electron potential conservation equation and an ion potential conservation equation; establishing a flow field model which comprises a gas phase mass conservation equation, a gas phase momentum conservation equation, a gas phase component conservation equation and a liquid water conservation equation; establishing a temperature field model including an energy conservation equation; setting the operation boundary conditions of the alkaline fuel cell, wherein the operation boundary conditions comprise anode and cathode inlet gas mole fractions, mass flow, temperature and output voltage. According to the three-dimensional multi-physical field model establishing method of the alkaline fuel cell considering the gas phase and the liquid phase, provided by the invention, based on the working principle of the alkaline fuel cell, a modeling method can be provided for characteristic research of the alkaline fuel cell.

Description

Technical Field

[0001] The invention belongs to the field of fuel cells, and in particular relates to a method and system for establishing an alkaline fuel cell model taking gas-liquid two phases into consideration. Background Art

[0002] With the continuous enrichment of human material life and the sustained development of economic level, the increasingly serious energy and environmental problems have become the key factors restricting the sustainable development of human social life. The consumption of non-renewable energy such as traditional mineral fuels is increasing day by day, the reserves are insufficient, and the large amount of pollutants produced during the combustion process has caused serious damage to the environment. Therefore, it is urgent to find and develop clean and renewable energy technologies.

[0003] As a clean and environmentally friendly energy conversion device, fuel cells can directly convert chemical energy into electrical energy, and have the advantages of high efficiency, no pollution, and low noise. The working principle of fuel cells involves complex physical and chemical processes, such as electrochemical reactions, fluid flow, heat and mass transfer, etc. These processes are coupled and influence each other, and jointly determine the performance of alkaline fuel cells. However, most of the existing research focuses on proton exchange membrane fuel cells, and there are fewer studies on alkaline fuel cells. Most of them use 0-dimensional mathematical models to reflect the efficiency characteristics of alkaline fuel cells. The internal operating mechanism cannot be described, and it is difficult to reflect more accurate equipment parameters during the operation of alkaline fuel cells, such as gas concentration, temperature, current density, etc. Summary of the invention

[0004] To solve the above problems, the present invention proposes a method and system for establishing an alkaline fuel cell model taking into account both gas and liquid phases. The method for establishing an alkaline fuel cell model taking into account both gas and liquid phases based on the working principle of alkaline fuel cells can provide a modeling method for studying the characteristics of alkaline fuel cells.

[0005] In order to achieve the above object, the present invention adopts the following technical scheme:

[0006] A method for establishing an alkaline fuel cell model considering gas-liquid two-phases, comprising:

[0007] Construct a three-dimensional geometric model of an alkaline fuel cell to provide a geometric model basis for the subsequent establishment of conservation equations and boundary conditions in electrochemistry, two-phase flow field, and temperature field;

[0008] Make reasonable assumptions about the alkaline fuel cell model and obtain the assumed conditions for the operation of the alkaline fuel cell model for simplification and calculation of various conservation equations;

[0009] Electrochemical models are established for different regions in the geometric model, including the electron potential conservation equation and the ion potential conservation equation, to describe the charge transfer process inside the alkaline fuel cell.

[0010] Two-phase flow field models are established for different regions in the geometric model, including the gas phase mass conservation equation, the gas phase momentum conservation equation, the gas phase component conservation equation, and the liquid water conservation equation, which are used to describe the flow and distribution of gas and liquid water inside the alkaline fuel cell.

[0011] The temperature field models are established for different regions in the geometric model, including the energy conservation equation, to describe the temperature distribution inside the alkaline fuel cell.

[0012] Set various boundary conditions for the electrochemical model, two-phase field model and temperature field model, including anode inlet gas mole fraction, cathode inlet gas mole fraction, anode inlet gas mass flow rate, cathode inlet gas mass flow rate, anode inlet temperature, cathode inlet temperature, anode plate potential, cathode plate potential and electrolyte concentration in the diaphragm;

[0013] Based on the electrochemical model, two-phase field model, temperature field model and various boundary conditions of the electrochemical model, two-phase field model and temperature field model of the three-dimensional geometric model of the alkaline fuel cell, a three-dimensional multi-physics field model of the alkaline fuel cell is established.

[0014] The present invention also provides an alkaline fuel cell modeling system considering gas-liquid two-phase, comprising the following modules:

[0015] The three-dimensional model building module builds a three-dimensional geometric model of an alkaline fuel cell, providing a geometric model foundation for the subsequent establishment of conservation equations and boundary conditions in electrochemistry, two-phase flow field, and temperature field;

[0016] Assumption module, which makes reasonable assumptions for the alkaline fuel cell model to facilitate the simplification and calculation of various conservation equations;

[0017] Electrochemical model building module, which builds electrochemical models for different regions in the geometric model, including the electron potential conservation equation and the ion potential conservation equation, which are used to describe the charge transfer process inside the alkaline fuel cell;

[0018] Two-phase flow field model construction module, which establishes two-phase flow field models for different regions in the geometric model, including gas phase mass conservation equation, gas phase momentum conservation equation, gas phase component conservation equation, and liquid water conservation equation, which are used to describe the flow and distribution of gas and liquid water inside the alkaline fuel cell;

[0019] The temperature field model building module builds temperature field models for different regions in the geometric model, including energy conservation equations, to describe the temperature distribution inside the alkaline fuel cell;

[0020] The boundary condition setting module sets various boundary conditions for the electrochemical model, two-phase field model and temperature field model. The various boundary conditions include anode inlet gas mole fraction, cathode inlet gas mole fraction, anode inlet gas mass flow rate, cathode inlet gas mass flow rate, anode inlet temperature, cathode inlet temperature, anode plate potential, cathode plate potential and electrolyte concentration in the diaphragm.

[0021] Furthermore, the three-dimensional geometric model of the alkaline fuel cell includes nine regions, namely, an anode plate, an anode flow channel, an anode diffusion layer, an anode catalyst layer, a diaphragm, a cathode catalyst layer, a cathode diffusion layer, a cathode flow channel, and a cathode plate.

[0022] Furthermore, the alkaline fuel cell model is run under the following assumptions:

[0023] The fluid flow in alkaline fuel cells is laminar; the cathode and anode reaction gases of alkaline fuel cells are regarded as ideal gases; there is no gas crossover between the anode and cathode; the influence of gravity is not considered; and the plate material is isotropic and uniform.

[0024] Furthermore, the electrochemical model of the alkaline fuel cell is:

[0025] Under the assumed conditions, the electron potential conservation equation is constructed for the plate, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0026] (1)

[0027] in, is the effective electronic conductivity, is the Laplace operator, represents the electron potential, is the source term of the electron potential conservation equation;

[0028] Under the assumed conditions, the ion potential conservation equation is constructed for the diaphragm and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0029] (2)

[0030] in, is the effective ionic conductivity, represents the ionic potential, represents the source term of the ion potential conservation equation;

[0031] The source terms of the electron potential and ion potential conservation equations are:

[0032] (3)

[0033] (4).

[0034] in, , represent the anode and cathode current densities respectively.

[0035] The Butler-Volmer equation is used to express the current density of the electrochemical reaction rates of the cathode and anode reactions:

[0036] (5)

[0037] (6)

[0038] in, represents the liquid water volume fraction, , denote the anode and cathode reference current densities, respectively. , , , are the molar concentrations of hydrogen, oxygen, water, and hydroxide ions, respectively. , , , Respectively represent the reference molar concentrations of hydrogen, oxygen, water, and hydroxide ions, , denote the anode and cathode transfer coefficients, respectively. is Avogadro's constant, is the ideal gas constant, is the temperature, , Respectively represent the anode and cathode activation overpotentials; exp() represents the exponential function;

[0039] The activation overpotential is calculated as:

[0040] (7)

[0041] (8)

[0042] in, , Represent the equilibrium potential of anode and cathode respectively;

[0043] Furthermore, the two-phase flow field model of the alkaline fuel cell is:

[0044] Under the assumed conditions, the gas phase mass conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0045] (9)

[0046] in, represents the porosity of the diffusion layer and the catalytic layer, and the flow channel is regarded as a straight channel. is 1, is the gas phase volume fraction, is the gas phase density, is the gas phase velocity, is the mass source term;

[0047] The mass source term is:

[0048] (10)

[0049] in, is the phase change source term between liquid water and gaseous water, , denote the molar masses of hydrogen and oxygen respectively;

[0050] The phase change source term of liquid water and gaseous water is:

[0051] (11)

[0052] in, , are the evaporation rate of liquid water and the condensation rate of gaseous water, respectively. represents the water vapor partial pressure, represents the saturated vapor pressure, Represents the molar mass of water. The sum of the volume fraction of liquid water and the volume fraction of gaseous water is always 1, that is, ;

[0053] The saturated vapor pressure is calculated as:

[0054] (12)

[0055] Under the assumed conditions, the gas phase momentum conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0056] (13)

[0057] in, represents the gas phase pressure, is the gas phase dynamic viscosity, is the source term of the momentum conservation equation, T represents the transpose, is the gas phase velocity;

[0058] The source term of the momentum conservation equation is expressed as:

[0059] (14)

[0060] in, represents the intrinsic permeability of the catalytic layer and the diffusion layer, Express the relative permeability:

[0061] (15)

[0062] Under the assumed conditions, the gas phase component conservation equations are constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0063] (16)

[0064] in, represents the mass fraction of gas i, i is hydrogen, oxygen, and water vapor, is the effective diffusion coefficient of gas i, is the gas phase component conservation source term;

[0065] The gas effective diffusion coefficient is corrected by the Brugmann model and is expressed as:

[0066] (17)

[0067] in, is the diffusion coefficient of gas i;

[0068] The gas phase component conservation source term is expressed as:

[0069] (18)

[0070] (19)

[0071] (20)

[0072] in, , , are the conservation source terms of hydrogen, oxygen and water vapor components respectively;

[0073] Under the assumed conditions, the liquid water conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0074] (twenty one)

[0075] in, is the density of liquid water, is the liquid phase velocity, is the liquid water source term;

[0076] Calculated by Darcy's law:

[0077] (twenty two)

[0078] (twenty three)

[0079] in, is the liquid phase pressure, is the liquid water dynamic viscosity, is the relative permeability of liquid water, expressed as:

[0080] (twenty four)

[0081] The transport of liquid water in porous media is determined by capillary pressure. drive:

[0082] (25)

[0083] in, is the gas phase pressure, is the capillary pressure, calculated using the LeFevreut-J formula:

[0084] (26)

[0085] in, is the surface tension coefficient of water, is the contact angle;

[0086] In summary, the conservation equation of liquid water is obtained from equations (21) to (26):

[0087] (27)

[0088] in, is the capillary diffusion coefficient:

[0089] (28)

[0090] The source term of the liquid water conservation equation is expressed as:

[0091] (29)

[0092] Under the assumed conditions, the K in the electrolyte of the catalyst layer and the diaphragm in the three-dimensional geometric model of the alkaline fuel cell is constructed. + OH - Concentration conservation equation:

[0093] (30)

[0094] Among them, j represents K+, OH-, represents the concentration of substance j, is the diffusion coefficient of substance j, The charge carried by the substance, is the mobility of substance j;

[0095] Since K+ does not participate in the electrochemical reaction, ;

[0096] Since no electrochemical reaction occurs in the diaphragm and there is no mass transfer across the phase boundary, the source term of the concentration of each substance in the diaphragm is 0;

[0097] The source term of the OH- conservation equation in the catalyst layer is:

[0098] (31).

[0099] Furthermore, the temperature field model of the alkaline fuel cell is:

[0100] Under the assumed conditions, the energy conservation equation is constructed for the plates and diaphragms in the three-dimensional geometric model of the alkaline fuel cell:

[0101] (32)

[0102] in, is the plate density, is the constant-pressure heat capacity of the plate, is the thermal conductivity of the plate, is the heat source term;

[0103] Under the assumed conditions, the energy conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0104] (33)

[0105] in is the average density, , is the average heat capacity at constant pressure, , is the average thermal conductivity, , , are the thermal conductivities of the gas phase and the liquid phase, respectively;

[0106] Heat source term It is expressed as:

[0107] (34)

[0108] Where h is the latent heat coefficient, , are the entropy changes of cathode and anode reactions, respectively.

[0109] Furthermore, the various boundary conditions of the alkaline fuel cell electrochemical model, two-phase field model and temperature field model are:

[0110] The mole fraction of the anode inlet gas and the cathode inlet gas mole fraction are calculated as follows:

[0111] (35)

[0112] (36)

[0113] in, , , , represent the mole fractions of water vapor and hydrogen at the anode, and water vapor and oxygen at the cathode, respectively. is the reference pressure, , Represent the anode and cathode intake humidity respectively;

[0114] The anode inlet gas mass flow rate and cathode inlet gas mass flow rate are calculated as follows:

[0115] (37)

[0116] (38)

[0117] in, , are the anode and cathode mass flow rates, respectively, , are the stoichiometric coefficients of the anode and cathode, respectively. , are the densities of hydrogen and oxygen, The reference working current used for gas flow calculation is generally 1A. , are the concentrations of hydrogen and oxygen, respectively. , are the anode and cathode inlet pressures, , are the inlet temperatures of the anode and cathode, respectively;

[0118] The anode plate potential and cathode plate potential are defined as follows:

[0119] (39)

[0120] (40)

[0121] in, , are the anode and cathode potentials, is the output voltage of the alkaline fuel cell.

[0122] The electrolyte concentration in the separator is .

[0123] The present invention also provides an electronic device, including a memory, a processor and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the above-mentioned alkaline fuel cell modeling method considering gas-liquid two-phases are implemented.

[0124] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps of the above-mentioned alkaline fuel cell modeling method considering gas-liquid two-phases are implemented.

[0125] Beneficial effects:

[0126] Based on the electrochemical, fluid mechanics and thermodynamics principles of alkaline fuel cells, the present invention comprehensively considers various physical and chemical effects inside alkaline fuel cells, establishes electrochemical models, two-phase fluid mechanics models, and temperature field models involved in alkaline fuel cells, clarifies the coupling relationship between the source terms of each conservation equation and each physical field, and then proposes a method for establishing a three-dimensional multi-physical field model of alkaline fuel cells that considers gas and liquid phases. Compared with the traditional 0-dimensional alkaline fuel cell mathematical model, this model can deeply reflect the internal gas concentration, temperature, liquid water volume fraction and other parameters during the operation of the alkaline fuel cell, and can further analyze the operating characteristics of the alkaline fuel cell under different parameters, thereby achieving a comprehensive prediction and optimization of the performance of the alkaline fuel cell, and providing technical support for the study of the characteristics of the alkaline fuel cell. BRIEF DESCRIPTION OF THE DRAWINGS

[0127] Figure 1 It is a schematic diagram of the three-dimensional geometric model of alkaline fuel cell;

[0128] Figure 2 is the gas concentration distribution diagram of alkaline fuel cell;

[0129] Figure 3 is the temperature distribution diagram of alkaline fuel cell;

[0130] Figure 4 This is the distribution diagram of liquid water volume fraction in alkaline fuel cell. DETAILED DESCRIPTION

[0131] To make the purpose and technical solution of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not limited to the present invention.

[0132] The present invention provides a three-dimensional multi-physical field modeling method for an alkaline fuel cell considering gas-liquid two-phase, comprising the following steps:

[0133] Step 1: construct a three-dimensional geometric model of an alkaline fuel cell to provide a geometric model basis for the subsequent establishment of conservation equations and boundary conditions in electrochemistry, two-phase flow field, and temperature field;

[0134] Step 2: Make reasonable assumptions about the three-dimensional geometric model of the alkaline fuel cell to facilitate the simplification and calculation of each conservation equation;

[0135] Step 3, establishing electrochemical models for different regions in the geometric model, including the electron potential conservation equation and the ion potential conservation equation, which are used to describe the charge transfer process inside the alkaline fuel cell;

[0136] Step 4, establishing two-phase flow field models for different regions in the geometric model, including gas phase mass conservation equation, gas phase momentum conservation equation, gas phase component conservation equation, and liquid water conservation equation, which are used to describe the flow and distribution of gas and liquid water inside the alkaline fuel cell;

[0137] Step 5: Establishing temperature field models for different regions in the geometric model, including energy conservation equations, to describe the temperature distribution inside the alkaline fuel cell;

[0138] Step 6. Set various boundary conditions for the electrochemical model, two-phase field model, and temperature field model. The various boundary conditions include anode inlet gas mole fraction, cathode inlet gas mole fraction, anode inlet gas mass flow rate, cathode inlet gas mass flow rate, anode inlet temperature, cathode inlet temperature, anode plate potential, cathode plate potential, and electrolyte concentration in the diaphragm.

[0139] Specifically, the step 1 includes: the three-dimensional geometric model of the alkaline fuel cell includes nine regions, namely, an anode plate, an anode flow channel, an anode diffusion layer, an anode catalyst layer, a diaphragm, a cathode catalyst layer, a cathode diffusion layer, a cathode flow channel, and a cathode plate.

[0140] Specifically, step 2 includes: the assumptions for the operation of the alkaline fuel cell model are: the fluid flow in the alkaline fuel cell is laminar flow; the cathode and anode reaction gases of the alkaline fuel cell are regarded as ideal gases; there is no gas crossover between the anode and the cathode; the influence of gravity is not considered; and the plate material is isotropic and uniform.

[0141] Specifically, the step 3 includes:

[0142] Under the assumed conditions, the electron potential conservation equation is constructed for the plate, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0143] (1)

[0144] in, is the effective electronic conductivity, is the Laplace operator, represents the electron potential, is the source term of the electron potential conservation equation;

[0145] Under the assumed conditions, the ion potential conservation equation is constructed for the diaphragm and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0146] (2)

[0147] in, is the effective ionic conductivity, represents the ionic potential, represents the source term of the ion potential conservation equation;

[0148] The source terms of the electron and ion potential conservation equations are:

[0149] (3)

[0150] (4)

[0151] in, , represent the anode and cathode current densities respectively.

[0152] The Butler-Volmer equation is used to express the current density of the electrochemical reaction rates of the cathode and anode reactions:

[0153] (5)

[0154] (6)

[0155] in, represents the liquid water volume fraction, , denote the anode and cathode reference current densities, respectively. , , , are the molar concentrations of hydrogen, oxygen, water and hydroxide ions, respectively, , , , denote the reference molar concentrations of hydrogen, oxygen, water and hydroxide ions, respectively. , denote the anode and cathode transfer coefficients, respectively. is Avogadro's constant, is the ideal gas constant, is the temperature, , Respectively represent the anode and cathode activation overpotentials; exp() represents the exponential function;

[0156] The activation overpotentials of anode and cathode are:

[0157] (7)

[0158] (8)

[0159] in, , Represent the equilibrium potential of anode and cathode respectively;

[0160] Specifically, the step 4 includes:

[0161] Under the assumed conditions, the gas phase mass conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0162] (9)

[0163] in, represents the porosity of the diffusion layer and the catalytic layer, and the flow channel is regarded as a straight channel. is 1, is the gas phase volume fraction, is the gas phase density, is the gas phase velocity, is the mass source term;

[0164] The mass source term is:

[0165] (10)

[0166] in, is the phase change source term between liquid water and gaseous water, , denote the molar masses of hydrogen and oxygen respectively;

[0167] The phase change source term of liquid water and gaseous water is:

[0168] (11)

[0169] in, , are the evaporation rate of liquid water and the condensation rate of gaseous water, respectively. represents the water vapor partial pressure, represents the saturated vapor pressure, Represents the molar mass of water. The sum of the volume fraction of liquid water and the volume fraction of gaseous water is always 1, that is, ;

[0170] The saturated vapor pressure is calculated as:

[0171] (12)

[0172] Under the assumed conditions, the gas phase momentum conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0173] (13)

[0174] in, represents the gas phase pressure, is the gas phase dynamic viscosity, is the source term of the momentum conservation equation; the superscript T represents the transpose of the matrix, is the gas phase velocity;

[0175] The source term of the momentum conservation equation is expressed as:

[0176] (14)

[0177] in, represents the intrinsic permeability of the catalytic layer and the diffusion layer, Express the relative permeability:

[0178] (15)

[0179] Under the assumed conditions, the gas phase component conservation equations are constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0180] (16)

[0181] in, represents the mass fraction of gas i, i is hydrogen, oxygen, and water vapor, is the effective diffusion coefficient of gas i, is the gas phase component conservation source term;

[0182] The effective diffusion coefficient of gas is corrected by the Bruggemann model and is expressed as:

[0183] (17)

[0184] in, is the diffusion coefficient of gas i;

[0185] The gas phase component conservation source term is expressed as:

[0186] (18)

[0187] (19)

[0188] (20)

[0189] in, , , are the conservation source terms of hydrogen, oxygen and water vapor components respectively;

[0190] Under the assumed conditions, the liquid water conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0191] (twenty one)

[0192] in, is the density of liquid water, is the liquid phase velocity, is the liquid water source term;

[0193] Calculated by Darcy's law:

[0194] (twenty two)

[0195] (twenty three)

[0196] in, is the liquid phase pressure, is the liquid water dynamic viscosity, is the relative permeability of liquid water, expressed as:

[0197] (twenty four)

[0198] The transport of liquid water in porous media is determined by capillary pressure. drive:

[0199] (25)

[0200] in, is the gas phase pressure, is the capillary pressure, calculated using the Leverett-J model formula:

[0201] (26)

[0202] in, is the surface tension coefficient of water, is the contact angle; It represents the intrinsic permeability of the catalytic layer and the diffusion layer;

[0203] In summary, the conservation equation of liquid water is obtained from equations (21) to (26):

[0204] (27)

[0205] in, is the capillary diffusion coefficient:

[0206] (28)

[0207] The source term of the liquid water conservation equation is expressed as:

[0208] (29)

[0209] Under the assumption, the catalyst layer and the diaphragm in the three-dimensional geometric model of the alkaline fuel cell are constructed to form the K + OH - Concentration conservation equation:

[0210] (30)

[0211] Where j represents K + OH - , represents the concentration of substance j, is the diffusion coefficient of substance j, The charge carried by the substance, is the mobility of substance j.

[0212] Due to K + It does not participate in the electrochemical reaction, so ;

[0213] Since no electrochemical reaction occurs in the diaphragm and there is no mass transfer across the phase boundary, the source term of the concentration of each substance in the diaphragm is 0;

[0214] OH in the catalyst layer - The source term of the conservation equation is:

[0215] (31)

[0216] Specifically, the step 5 includes:

[0217] Under the assumed conditions, the energy conservation equation is constructed for the plates and diaphragms in the three-dimensional geometric model of the alkaline fuel cell:

[0218] (32)

[0219] in, is the plate density, is the constant-pressure heat capacity of the plate, is the thermal conductivity of the plate, is the heat source term;

[0220] Under the assumed conditions, the energy conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell:

[0221] (33)

[0222] in is the average density, , is the average heat capacity at constant pressure, , is the average thermal conductivity, , , are the thermal conductivities of the gas phase and the liquid phase, respectively;

[0223] Heat source term It is expressed as:

[0224] (34)

[0225] Where h is the latent heat coefficient, , are the entropy changes of cathode and anode reactions, respectively.

[0226] Specifically, the step 6 includes:

[0227] The various boundary conditions for the electrochemical model, two-phase field model, and temperature field model are set as follows:

[0228] The mole fraction of the anode inlet gas and the cathode inlet gas mole fraction are calculated as follows:

[0229] (35)

[0230] (36)

[0231] in, , , , represent the mole fractions of water vapor and hydrogen at the anode, and water vapor and oxygen at the cathode, respectively. is the reference pressure, , Represent the anode and cathode intake humidity respectively;

[0232] The anode inlet gas mass flow rate and cathode inlet gas mass flow rate are calculated as follows:

[0233] (37)

[0234] (38)

[0235] in, , are the anode and cathode mass flow rates, respectively, , are the stoichiometric coefficients of the anode and cathode, respectively. , are the densities of hydrogen and oxygen, The reference working current used for gas flow calculation is generally 1A. , are the concentrations of hydrogen and oxygen, respectively. , are the anode and cathode inlet pressures, , are the inlet temperatures of the anode and cathode, respectively.

[0236] The anode plate potential and cathode plate potential are defined as follows:

[0237] (39)

[0238] (40)

[0239] in, , are the anode and cathode potentials, is the output voltage of the alkaline fuel cell.

[0240] The electrolyte concentration in the separator is set to .

[0241] By establishing equations according to the above steps, a three-dimensional multi-physics field model of alkaline fuel cells considering gas-liquid two-phase can be constructed to realize the study of the operating characteristics of alkaline fuel cells under different operating parameters.

[0242] Example:

[0243] Build a 3D geometric model of an alkaline fuel cell Figure 1 As shown in the figure, the geometric model includes nine regions, namely, anode plate, anode flow channel, anode diffusion layer, anode catalyst layer, diaphragm, cathode catalyst layer, cathode diffusion layer, cathode flow channel and cathode plate. The relevant dimensions are shown in Table 1:

[0244] Table 1 Dimensional parameters of alkaline fuel cell

[0245]

[0246] In the electrochemical model, the relevant parameter settings are shown in Table 2:

[0247] Table 2 Electrochemical model parameters

[0248]

[0249] In the two-phase flow field model, the relevant parameter settings are shown in Table 3:

[0250] Table 3 Two-phase flow model parameters

[0251]

[0252] In the temperature field model, the relevant parameter settings are shown in Table 4:

[0253] Table 4 Temperature field model parameters

[0254]

[0255] Set the boundary conditions:

[0256] Anode and cathode air humidity , Set to 100%, anode and cathode inlet temperature , At 70℃, the stoichiometric coefficients of anode and cathode are , 2 and 3 respectively, the cathode and anode inlet pressures are 1atm, and the reference pressure is 1atm.

[0257] The calculated anode hydrogen inlet mole fraction is 0.6929, and the mass flow rate is 3.75×10 -8 kg / s, the cathode oxygen inlet mole fraction is 0.6929, and the mass flow rate is 4.29×10 -6 kg / s.

[0258] Electrolyte concentration in the separator It is 6M.

[0259] Substituting the above parameters into the physical field models, the distribution of parameters such as gas concentration, temperature, liquid water volume fraction, etc. inside the alkaline fuel cell under the corresponding boundary conditions can be calculated as follows: Figure 2 , Figure 3 , Figure 4 shown.

[0260] The present invention also provides a three-dimensional multi-physical field modeling system for an alkaline fuel cell considering gas-liquid two-phases, comprising the following modules:

[0261] The three-dimensional model building module builds a three-dimensional geometric model of an alkaline fuel cell, providing a geometric model foundation for the subsequent establishment of conservation equations and boundary conditions in electrochemistry, two-phase flow field, and temperature field;

[0262] Assumption module, which makes reasonable assumptions for the alkaline fuel cell model to facilitate the simplification and calculation of various conservation equations;

[0263] Electrochemical model building module, which builds electrochemical models for different regions in the geometric model, including the electron potential conservation equation and the ion potential conservation equation, which are used to describe the charge transfer process inside the alkaline fuel cell;

[0264] Two-phase flow field model construction module, which establishes two-phase flow field models for different regions in the geometric model, including gas phase mass conservation equation, gas phase momentum conservation equation, gas phase component conservation equation, and liquid water conservation equation, which are used to describe the flow and distribution of gas and liquid water inside the alkaline fuel cell;

[0265] The temperature field model building module builds temperature field models for different regions in the geometric model, including energy conservation equations, to describe the temperature distribution inside the alkaline fuel cell;

[0266] The boundary condition setting module sets various boundary conditions for the electrochemical model, two-phase field model and temperature field model, including the anode and cathode inlet gas mole fraction, mass flow rate, temperature, output voltage, and electrolyte concentration, and finally establishes a three-dimensional multi-physics field model of the alkaline fuel cell.

[0267] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps of the above-mentioned method for modeling a three-dimensional multi-physical field model of an alkaline fuel cell taking into account gas-liquid two phases are implemented.

[0268] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-mentioned method for modeling a three-dimensional multi-physical field model of an alkaline fuel cell taking into account gas-liquid two-phases.

[0269] It should be understood by those skilled in the art that the embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes. The solutions in the embodiments of the present invention may be implemented in various computer languages, for example, object-oriented programming language Java and interpreted scripting language JavaScript, etc.

[0270] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0271] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0272] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0273] Although the preferred embodiments of the present invention have been described, those skilled in the art may make other changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present invention.

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

Claims

1. A method for establishing an alkaline fuel cell model considering gas-liquid two-phase, characterized in that: include: Construct a three-dimensional geometric model of an alkaline fuel cell; Make reasonable assumptions about the three-dimensional geometric model of the alkaline fuel cell and determine the assumed conditions for the operation of the alkaline fuel cell model; Under the assumed conditions, electrochemical models are established for different regions in the three-dimensional geometric model of alkaline fuel cells. The electrochemical models include the electron potential conservation equation and the ion potential conservation equation. Under the assumed conditions, two-phase flow field models are established for different regions in the geometric model. The two-phase flow field models include the gas phase mass conservation equation, the gas phase momentum conservation equation, the gas phase component conservation equation, and the liquid water conservation equation. Under the assumed conditions, temperature field models are established for different regions in the geometric model, and the temperature field model includes energy conservation equations; Set various boundary conditions for the electrochemical model, two-phase field model and temperature field model, including anode inlet gas mole fraction, cathode inlet gas mole fraction, anode inlet gas mass flow rate, cathode inlet gas mass flow rate, anode inlet temperature, cathode inlet temperature, anode plate potential, cathode plate potential and electrolyte concentration in the diaphragm; Based on the electrochemical model, two-phase field model, temperature field model and various boundary conditions of the electrochemical model, two-phase field model and temperature field model of the three-dimensional geometric model of the alkaline fuel cell, a three-dimensional multi-physics field model of the alkaline fuel cell is established.

2. The method for establishing an alkaline fuel cell model considering gas-liquid two-phase according to claim 1, characterized in that: The three-dimensional geometric model of the alkaline fuel cell includes nine regions, namely, an anode plate, an anode flow channel, an anode diffusion layer, an anode catalyst layer, a diaphragm, a cathode catalyst layer, a cathode diffusion layer, a cathode flow channel and a cathode plate.

3. The method for establishing an alkaline fuel cell model considering gas-liquid two-phase according to claim 1, characterized in that: Reasonable assumptions are made for the three-dimensional geometric model of the alkaline fuel cell. The assumptions for the operation of the alkaline fuel cell model include: The fluid flow in alkaline fuel cells is laminar; the cathode and anode reaction gases of alkaline fuel cells are regarded as ideal gases; there is no gas crossover between the anode and cathode; the influence of gravity is not considered; and the plate material is isotropic and uniform.

4. The method for establishing an alkaline fuel cell model considering gas-liquid two-phase according to claim 1, characterized in that: The electrochemical model of an alkaline fuel cell is: Under the assumed conditions, the electron potential conservation equation is constructed for the plate, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （1） in, is the effective electronic conductivity, is the Laplace operator, represents the electron potential, is the source term of the electron potential conservation equation; Under the assumed conditions, the ion potential conservation equation is constructed for the diaphragm and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （2） in, is the effective ionic conductivity, represents the ionic potential, represents the source term of the ion potential conservation equation; The source terms of the electron potential conservation equation and the ion potential conservation equation are: （3） （4） in, , represent the anode and cathode current densities respectively; The Butler-Volmer equation is used to express the current density of the electrochemical reaction rates of the cathode and anode reactions: （5） （6） in, represents the liquid water volume fraction, , denote the anode and cathode reference current densities, respectively. , , , are the molar concentrations of hydrogen, oxygen, water, and hydroxide ions, respectively. , , , Respectively represent the reference molar concentrations of hydrogen, oxygen, water, and hydroxide ions, , denote the anode and cathode transfer coefficients, respectively. is Avogadro's constant, is the ideal gas constant, is the temperature, , Respectively represent the anode and cathode activation overpotentials; exp() represents the exponential function; The activation overpotential is calculated as: （7） （8） in, , represent the equilibrium potentials of the anode and cathode respectively.

5. The method for establishing an alkaline fuel cell model considering gas-liquid two phases according to claim 1, characterized in that: The two-phase flow field model of alkaline fuel cell is: Under the assumed conditions, the gas phase mass conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （9） in, represents the porosity of the diffusion layer and the catalytic layer, and the flow channel is regarded as a straight channel. is 1, is the gas phase volume fraction, is the gas phase density, is the gas phase velocity, is the mass source term; The mass source term is: （10） in, is the phase change source term between liquid water and gaseous water, , denote the molar masses of hydrogen and oxygen respectively; The phase change source term of liquid water and gaseous water is: （11） in, , are the evaporation rate of liquid water and the condensation rate of gaseous water, respectively. represents the water vapor partial pressure, represents the saturated vapor pressure, Represents the molar mass of water. The sum of the volume fraction of liquid water and the volume fraction of gaseous water is always 1, that is, ; The saturated vapor pressure is calculated as: （12） Under the assumed conditions, the gas phase momentum conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （13） in, represents the gas phase pressure, is the gas phase dynamic viscosity, is the source term of the momentum conservation equation, T represents the transpose, is the gas phase velocity; The source term of the momentum conservation equation is expressed as: （14） in, represents the intrinsic permeability of the catalytic layer and the diffusion layer, Express the relative permeability: （15） Under the assumed conditions, the gas phase component conservation equations are constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （16） in, represents the mass fraction of gas i, i is hydrogen, oxygen, and water vapor, is the effective diffusion coefficient of gas i, is the gas phase component conservation source term; The gas effective diffusion coefficient is corrected by the Brugmann model and is expressed as: （17） in, is the diffusion coefficient of gas i; The gas phase component conservation source term is expressed as: （18） （19） （20） in, , , are the conservation source terms of hydrogen, oxygen and water vapor components respectively; Under the assumed conditions, the liquid water conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （21） in, is the density of liquid water, is the liquid phase velocity, is the liquid water source term; Calculated by Darcy's law: （22） （23） in, is the liquid phase pressure, is the liquid water dynamic viscosity, is the relative permeability of liquid water, expressed as: （24） The transport of liquid water in porous media is determined by capillary pressure. drive: （25） in, is the gas phase pressure, is the capillary pressure, calculated using the LeFevreut-J formula: （26） in, is the surface tension coefficient of water, is the contact angle; In summary, the conservation equation of liquid water is obtained from equations (21) to (26): （27） in, is the capillary diffusion coefficient: （28） The source term of the liquid water conservation equation is expressed as: （29） Under the assumption, the catalyst layer and the diaphragm in the three-dimensional geometric model of the alkaline fuel cell are constructed to form the K + OH - Concentration conservation equation: （30） Among them, j represents K+, OH-, represents the concentration of substance j, is the diffusion coefficient of substance j, The charge carried by the substance, is the mobility of substance j; Since K+ does not participate in the electrochemical reaction, ; Since no electrochemical reaction occurs in the diaphragm and there is no mass transfer across the phase boundary, the source term of the concentration of each substance in the diaphragm is 0; The source term of the OH- conservation equation in the catalyst layer is: （31）。 6. The method for establishing an alkaline fuel cell model considering gas-liquid two-phase according to claim 1, characterized in that: The temperature field model of alkaline fuel cell is: Under the assumed conditions, the energy conservation equation is constructed for the plates and diaphragms in the three-dimensional geometric model of the alkaline fuel cell: （32） in, is the plate density, is the constant-pressure heat capacity of the plate, is the thermal conductivity of the plate, is the heat source term; Under the assumed conditions, the energy conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （33） in is the average density, , is the average heat capacity at constant pressure, , is the average thermal conductivity, , , are the thermal conductivities of the gas phase and the liquid phase, respectively; Heat source term It is expressed as: （34） Where h is the latent heat coefficient, , are the entropy changes of cathode and anode reactions, respectively.

7. The method for establishing an alkaline fuel cell model considering gas-liquid two-phase according to claim 1, characterized in that: The various boundary conditions of the electrochemical model, two-phase field model and temperature field model of the alkaline fuel cell are: The mole fraction of the anode inlet gas and the cathode inlet gas mole fraction are calculated as follows: （35） （36） in, , , , represent the mole fractions of water vapor and hydrogen at the anode, and water vapor and oxygen at the cathode, respectively. is the reference pressure, , Represent the anode and cathode intake humidity respectively; The anode inlet gas mass flow rate and cathode inlet gas mass flow rate are calculated as follows: （37） （38） in, , are the anode and cathode mass flow rates, respectively, , are the stoichiometric coefficients of the anode and cathode, respectively. , are the densities of hydrogen and oxygen, The reference working current used for gas flow calculation is generally 1A. , are the concentrations of hydrogen and oxygen, respectively. , are the anode and cathode inlet pressures, , are the inlet temperatures of the anode and cathode, respectively; The anode plate potential and cathode plate potential are defined as follows: （39） （40） in, , are the anode and cathode potentials, is the output voltage of the alkaline fuel cell; The electrolyte concentration in the separator is .

8. The present invention also provides an alkaline fuel cell modeling system taking into account gas-liquid two-phases, characterized in that: Includes the following modules: A three-dimensional model building module is used to build a three-dimensional geometric model of an alkaline fuel cell; Make reasonable assumptions about the three-dimensional geometric model of the alkaline fuel cell and determine the assumed conditions for the operation of the alkaline fuel cell model; Hypothesis module, under hypothetical conditions, establishes electrochemical models for different regions in the three-dimensional geometric model of alkaline fuel cells, and the electrochemical models include the electron potential conservation equation and the ion potential conservation equation; The electrochemical model building module establishes two-phase flow field models for different regions in the geometric model under assumed conditions. The two-phase flow field model includes the gas phase mass conservation equation, the gas phase momentum conservation equation, the gas phase component conservation equation, and the liquid water conservation equation. The two-phase flow field model building module builds temperature field models for different regions in the geometric model under assumed conditions. The temperature field model includes energy conservation equations. The temperature field model construction module sets various boundary conditions for the electrochemical model, two-phase field model and temperature field model. The various boundary conditions include anode inlet gas mole fraction, cathode inlet gas mole fraction, anode inlet gas mass flow rate, cathode inlet gas mass flow rate, anode inlet temperature, cathode inlet temperature, anode plate potential, cathode plate potential, and electrolyte concentration in the diaphragm. The boundary condition setting module establishes a three-dimensional multi-physical field model of an alkaline fuel cell based on the electrochemical model, two-phase field model, temperature field model and various boundary conditions of the electrochemical model, two-phase field model, temperature field model and the three-dimensional geometric model of the alkaline fuel cell.

9. The alkaline fuel cell model establishment system considering gas-liquid two-phase according to claim 8, characterized in that: The three-dimensional geometric model of the alkaline fuel cell includes nine regions, namely, an anode plate, an anode flow channel, an anode diffusion layer, an anode catalyst layer, a diaphragm, a cathode catalyst layer, a cathode diffusion layer, a cathode flow channel and a cathode plate.

10. The alkaline fuel cell model establishment system considering gas-liquid two-phase according to claim 8, characterized in that: The alkaline fuel cell model operates under the following assumptions: The fluid flow in alkaline fuel cells is laminar; the cathode and anode reaction gases of alkaline fuel cells are regarded as ideal gases; there is no gas crossover between the anode and cathode; the influence of gravity is not considered; and the plate material is isotropic and uniform.

11. The alkaline fuel cell model establishment system considering gas-liquid two-phase according to claim 8, characterized in that: The electrochemical model of an alkaline fuel cell is: Under the assumed conditions, the electron potential conservation equation is constructed for the plate, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （1） in, is the effective electronic conductivity, is the Laplace operator, represents the electron potential, is the source term of the electron potential conservation equation; Under the assumed conditions, the ion potential conservation equation is constructed for the diaphragm and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （2） in, is the effective ionic conductivity, represents the ionic potential, represents the source term of the ion potential conservation equation; The source terms of the electron potential and ion potential conservation equations are: （3） （4） in, , represent the anode and cathode current densities respectively; The Butler-Volmer equation is used to express the current density of the electrochemical reaction rates of the cathode and anode reactions: （5） （6） in, represents the liquid water volume fraction, , denote the anode and cathode reference current densities, respectively. , , , are the molar concentrations of hydrogen, oxygen, water, and hydroxide ions, respectively. , , , Respectively represent the reference molar concentrations of hydrogen, oxygen, water, and hydroxide ions, , denote the anode and cathode transfer coefficients, respectively. is Avogadro's constant, is the ideal gas constant, is the temperature, , Respectively represent the anode and cathode activation overpotentials; exp() represents the exponential function; The activation overpotential is calculated as: （7） （8） in, , represent the equilibrium potentials of the anode and cathode respectively.

12. The alkaline fuel cell model establishment system considering gas-liquid two-phase according to claim 8, characterized in that: The two-phase flow field model of the alkaline fuel cell is: Under the assumed conditions, the gas phase mass conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （9） in, represents the porosity of the diffusion layer and the catalytic layer, and the flow channel is regarded as a straight channel. is 1, is the gas phase volume fraction, is the gas phase density, is the gas phase velocity, is the mass source term; The mass source term is: （10） in, is the phase change source term between liquid water and gaseous water, , denote the molar masses of hydrogen and oxygen respectively; The phase change source term of liquid water and gaseous water is: （11） in, , are the evaporation rate of liquid water and the condensation rate of gaseous water, respectively. represents the water vapor partial pressure, represents the saturated vapor pressure, Represents the molar mass of water. The sum of the volume fraction of liquid water and the volume fraction of gaseous water is always 1, that is, ; The saturated vapor pressure is calculated as: （12） Under the assumed conditions, the gas phase momentum conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （13） in, represents the gas phase pressure, is the gas phase dynamic viscosity, is the source term of the momentum conservation equation, T represents the transpose, is the gas phase velocity; The source term of the momentum conservation equation is expressed as: （14） in, represents the intrinsic permeability of the catalytic layer and the diffusion layer, Express the relative permeability: （15） Under the assumed conditions, the gas phase component conservation equations are constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （16） in, represents the mass fraction of gas i, i is hydrogen, oxygen, and water vapor, is the effective diffusion coefficient of gas i, is the gas phase component conservation source term; The gas effective diffusion coefficient is corrected by the Brugmann model and is expressed as: （17） in, is the diffusion coefficient of gas i; The gas phase component conservation source term is expressed as: （18） （19） （20） in, , , are the conservation source terms of hydrogen, oxygen and water vapor components respectively; Under the assumed conditions, the liquid water conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （21） in, is the density of liquid water, is the liquid phase velocity, is the liquid water source term; Calculated by Darcy's law: （22） （23） in, is the liquid phase pressure, is the liquid water dynamic viscosity, is the relative permeability of liquid water, expressed as: （24） The transport of liquid water in porous media is determined by capillary pressure. drive: （25） in, is the gas phase pressure, is the capillary pressure, calculated using the LeFevreut-J formula: （26） in, is the surface tension coefficient of water, is the contact angle; In summary, the conservation equation of liquid water is obtained from equations (21) to (26): （27） in, is the capillary diffusion coefficient: （28） The source term of the liquid water conservation equation is expressed as: （29） Under the assumed conditions, the K in the electrolyte of the catalyst layer and the diaphragm in the three-dimensional geometric model of the alkaline fuel cell is constructed. + OH - Concentration conservation equation: （30） Among them, j represents K+, OH-, represents the concentration of substance j, is the diffusion coefficient of substance j, The charge carried by the substance, is the mobility of substance j; Since K+ does not participate in the electrochemical reaction, ; Since no electrochemical reaction occurs in the diaphragm and there is no mass transfer across the phase boundary, the source term of the concentration of each substance in the diaphragm is 0; The source term of the OH- conservation equation in the catalyst layer is: （31）。 13. The alkaline fuel cell model establishment system considering gas-liquid two-phase according to claim 8, characterized in that: The temperature field model of alkaline fuel cell is: Under the assumed conditions, the energy conservation equation is constructed for the plates and diaphragms in the three-dimensional geometric model of the alkaline fuel cell: （32） in, is the plate density, is the constant-pressure heat capacity of the plate, is the thermal conductivity of the plate, is the heat source term; Under the assumed conditions, the energy conservation equation is constructed for the flow channel, diffusion layer, and catalyst layer in the three-dimensional geometric model of the alkaline fuel cell: （33） in, is the average density, , is the average heat capacity at constant pressure, , is the average thermal conductivity, , , are the thermal conductivities of the gas phase and the liquid phase, respectively; Heat source term It is expressed as: （34） Where h is the latent heat coefficient, , are the entropy changes of cathode and anode reactions, respectively.

14. The alkaline fuel cell model establishment system considering gas-liquid two-phase according to claim 8, characterized in that: The various boundary conditions of the electrochemical model, two-phase field model and temperature field model of the alkaline fuel cell are: The mole fraction of the anode inlet gas and the cathode inlet gas mole fraction are calculated as follows: The calculations are as follows: （35） （36） in, , , , represent the mole fractions of water vapor and hydrogen at the anode, and water vapor and oxygen at the cathode, respectively. is the reference pressure, , Represent the anode and cathode intake humidity respectively; The anode inlet gas mass flow rate and cathode inlet gas mass flow rate are calculated as follows: （37） （38） in, , are the anode and cathode mass flow rates, respectively, , are the stoichiometric coefficients of the anode and cathode, respectively. , are the densities of hydrogen and oxygen, The reference working current used for gas flow calculation is generally 1A. , are the concentrations of hydrogen and oxygen, respectively. , are the anode and cathode inlet pressures, , are the inlet temperatures of the anode and cathode, respectively; The anode plate potential and cathode plate potential are defined as follows: （39） （40） in, , are the anode and cathode potentials, respectively, is the output voltage of the alkaline fuel cell; The electrolyte concentration in the separator is .

15. Provide an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the steps of a method for modeling an alkaline fuel cell model taking into account gas-liquid two phases as described in claims 1-7 are implemented.

16. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of a method for modeling an alkaline fuel cell model taking into account gas-liquid two phases as described in claims 1-7 are implemented.

Images