A computational method for two-phase flow simulation in an electrolytic cell

Abstract

The present invention discloses a calculation method for fitting two-phase flow simulation of an entire electrolytic cell, and relates to the technical field of electrolytic cell flow simulation. The method performs gas-liquid two-phase flow simulation on the gas and electrolyte mixed discharge section and a single electrolysis chamber in the anode and cathode electrolysis regions of the electrolytic cell, obtaining a pressure variation curve of the gas and electrolyte mixed discharge section and a two-phase pressure distribution cloud map of the single electrolysis chamber in the anode and cathode electrolysis regions; then, liquid phase flow simulation is performed on the gas and electrolyte mixed discharge section and the single electrolysis chamber in the anode and cathode electrolysis regions, fitting the pressure variation curve and the single-phase pressure distribution cloud map to determine the value of the flow resistance of each section; based on the value of the flow resistance of each section, liquid phase flow simulation is performed on the entire electrolytic cell to determine the fitting gas-liquid two-phase flow simulation result of the electrolytic cell. This method can reduce the time required for performing gas-liquid two-phase flow simulation calculations for the entire electrolytic cell.

Description

Technical Field

[0001] The present invention relates to the technical field of electrolytic cell flow simulation, and in particular to a calculation method for electrolytic cell whole cell fitting two-phase flow simulation. Background Art

[0002] However, as the demand for new energy consumption continues to increase, the requirements for electrolyzer size and output are also increasing.

[0003] The market demand for large-sized electrolyzers is constantly increasing. However, due to sealing requirements and size limitations, physical parameters such as mass flow distribution, temperature distribution, and pressure distribution in large electrolyzers are difficult to measure experimentally. Therefore, in order to obtain data such as the mass flow distribution and temperature distribution of the electrolyte inside the electrolyzer, computational fluid dynamics simulation is usually used to perform a three-dimensional simulation of the interior of the electrolyzer. During the water electrolysis process, the electrolyte completes the electrolysis reaction at the anode and cathode electrodes, producing hydrogen and oxygen, which mix with the remaining electrolyte and are discharged from the electrolyzer. The process is a gas-liquid two-phase flow.

[0004] However, due to the large size of the large square electrolytic cell, the simulation of gas-liquid two-phase flow in the entire electrolytic cell has the problem of long simulation calculation time. Summary of the Invention

[0005] Based on this, it is necessary to provide a calculation method for fitting two-phase flow simulation of the entire electrolytic cell to address the above technical problems. This method can reduce the time for gas-liquid two-phase flow simulation calculation of the entire electrolytic cell.

[0006] The present invention adopts the following technical solutions:

[0007] The present invention provides a calculation method for two-phase flow simulation of an electrolytic cell, comprising:

[0008] The electrolytic cell is divided into a grid to obtain the electrolyte inlet section, the cathode and anode electrolysis areas, and the gas and electrolyte mixed discharge section of the electrolytic cell;

[0009] Gas-liquid two-phase flow simulations were performed for the gas and electrolyte mixed discharge section and a single electrolysis chamber in the anode and cathode electrolysis area, respectively. The pressure variation curve of the gas and electrolyte mixed discharge section with position and the two-phase pressure distribution cloud diagram of a single electrolysis chamber in the anode and cathode electrolysis area were obtained.

[0010] During the liquid phase flow simulation of the gas and electrolyte mixed discharge section, the flow resistance value of the gas and electrolyte mixed discharge section is determined by fitting the pressure-position curve by applying flow resistance in the liquid phase flow area.

[0011] In the process of liquid phase flow simulation of a single electrolysis chamber in the cathode electrolysis area, flow resistance is applied in the cathode electrolysis area and the anode electrolysis area respectively until the single-phase pressure distribution cloud and the two-phase pressure distribution cloud are consistent. Figure 1 To determine the flow resistance values ​​of the cathode electrolysis area and the anode electrolysis area;

[0012] According to the values ​​of each flow resistance, the liquid phase flow simulation of the entire electrolytic cell is carried out to determine the fitted gas-liquid two-phase flow simulation results of the electrolytic cell.

[0013] Preferably, the gas and electrolyte mixed discharge section includes a hydrogen and electrolyte mixed discharge section and an oxygen and electrolyte mixed discharge section; and the electrolytic cell is grid-divided, including:

[0014] According to the different flow states and functions of each part of the electrolytic cell, the electrolytic cell is modeled into the electrolyte inlet section, cathode electrolysis area, anode electrolysis area, hydrogen and electrolyte mixed discharge section, and oxygen and electrolyte mixed discharge section;

[0015] The modeled electrolyte inlet section, cathode electrolysis area, anode electrolysis area, hydrogen and electrolyte mixed discharge section, and oxygen and electrolyte mixed discharge section are respectively divided into fluid calculation grids.

[0016] Preferably, a gas-liquid two-phase flow simulation is performed on the gas and electrolyte mixed discharge section to obtain a pressure variation curve of the gas and electrolyte mixed discharge section with position, including:

[0017] Perform gas-liquid two-phase flow simulation on the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section, and extract the pressure values ​​at different points in the flow direction of the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section respectively;

[0018] According to the position of each point and the corresponding pressure value, the pressure variation curve of the hydrogen and electrolyte mixed discharge section and the pressure variation curve of the oxygen and electrolyte mixed discharge section are drawn respectively.

[0019] Preferably, a gas-liquid two-phase flow simulation is performed on a single electrolysis chamber in the anode and cathode electrolysis regions to obtain a two-phase pressure distribution cloud diagram of the single electrolysis chamber in the anode and cathode electrolysis regions, including:

[0020] Gas-liquid two-phase flow simulation is performed on a single electrolysis chamber in the cathode electrolysis area and a single electrolysis chamber in the anode electrolysis area to obtain a two-phase pressure distribution cloud map of a single electrolysis chamber in the cathode electrolysis area and a two-phase pressure distribution cloud map of a single electrolysis chamber in the anode electrolysis area.

[0021] Preferably, the two-phase flow control equations in the gas-liquid two-phase flow simulation include the gas-liquid two-phase flow mass conservation equation and the gas-liquid two-phase flow momentum conservation equation;

[0022] The mass conservation equation for gas-liquid two-phase flow is:

[0023]

[0024] The momentum conservation equation for gas-liquid two-phase flow is:

[0025]

[0026] Where, l is the liquid phase, g is the gas phase; t is the time, s is the volume fraction, ρ is the fluid density, u is the fluid velocity, and P is the pressure. is the stress-strain tensor, S m is the fluid mass source term, S u is the fluid momentum source term.

[0027] Preferably, the liquid phase flow control equations in the liquid phase flow simulation include the liquid phase mass conservation equation and the liquid phase momentum conservation equation;

[0028] Liquid mass conservation equation:

[0029]

[0030] Liquid phase momentum conservation equation:

[0031]

[0032] Among them, the flow resistance applied in the liquid phase flow region represents the regulation of the fluid momentum source term.

[0033] Preferably, the liquid phase flow region is a porous medium region; the liquid phase flow control equations in the liquid phase flow simulation include the porous medium liquid phase flow mass conservation equation and the porous medium liquid phase flow momentum conservation equation;

[0034] Mass conservation equation for liquid flow in porous media:

[0035]

[0036] Momentum conservation equation for liquid flow in porous media:

[0037]

[0038] S ul =-(Aμu i +Bρ|u|u i )

[0039] Where ε is the porosity, μ is the fluid dynamic viscosity, A is the viscous drag coefficient, and B is the inertial drag coefficient. The fluid momentum source term in the liquid phase is adjusted by the viscous drag coefficient and the inertial drag coefficient.

[0040] Preferably, the fitting gas-liquid two-phase flow simulation result includes the flow and heat exchange parameter values ​​of each electrolytic chamber in the electrolytic cell; according to the values ​​of each flow resistance, the liquid phase flow simulation of the entire electrolytic cell is performed to determine the fitting gas-liquid two-phase flow simulation result of the electrolytic cell, including:

[0041] The electrolyte inlet section, the cathode and anode electrolysis areas, and the gas and electrolyte mixed discharge section are geometrically and grid-combined to form the entire electrolytic cell;

[0042] The flow resistance values ​​of each part are set to the corresponding parts of the entire electrolytic cell, and the liquid phase flow simulation of the entire electrolytic cell is performed to obtain the flow and heat transfer parameter values ​​in each electrolytic chamber.

[0043] The present invention provides a computing device for simulating two-phase flow in an electrolytic cell, comprising:

[0044] A partitioning module is used to divide the electrolytic cell into grids to obtain the electrolyte inlet section, the cathode and anode electrolysis areas, and the gas and electrolyte mixed discharge section of the electrolytic cell;

[0045] The first simulation module is used to simulate the gas-liquid two-phase flow in the gas and electrolyte mixed discharge section and the single electrolysis chamber in the anode and cathode electrolysis area, respectively, to obtain the pressure variation curve of the gas and electrolyte mixed discharge section with position, and the two-phase pressure distribution cloud diagram of the single electrolysis chamber in the anode and cathode electrolysis area;

[0046] The second simulation module is used to determine the value of the flow resistance of the gas and electrolyte mixed discharge section by fitting the pressure-position variation curve by applying flow resistance in the liquid phase flow region during the liquid phase flow simulation of the gas and electrolyte mixed discharge section;

[0047] The third simulation module is used to simulate the liquid phase flow in a single electrolysis chamber in the cathode electrolysis area by applying flow resistance in the cathode electrolysis area and the anode electrolysis area respectively until the single-phase pressure distribution cloud and the two-phase pressure distribution cloud are consistent. Figure 1 To determine the flow resistance values ​​of the cathode electrolysis area and the anode electrolysis area;

[0048] The fourth simulation module is used to perform liquid phase flow simulation on the entire electrolytic cell according to the values ​​of each flow resistance, and determine the fitted gas-liquid two-phase flow simulation results of the electrolytic cell.

[0049] The present invention provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the calculation method for the above-mentioned electrolytic cell whole-tank fitting two-phase flow simulation is implemented.

[0050] The present invention provides a computer 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 calculation method for the two-phase flow simulation of the entire electrolytic cell is implemented.

[0051] At least one of the above technical solutions adopted by the present invention can achieve the following beneficial effects:

[0052] In the present invention, the electrolytic cell is divided into multiple parts, and liquid-phase flow simulation is performed on each part respectively to fit the simulation results of the gas-liquid two-phase flow simulation, and the flow resistance of each part is determined so that the liquid-phase flow simulation corresponding to the flow resistance is consistent with the effect of the gas-liquid two-phase flow simulation. Therefore, the flow resistance corresponding to each part can be set to the corresponding part, and the liquid-phase flow simulation of the entire electrolytic cell is performed through the flow resistance of the corresponding part to simulate the gas-liquid two-phase flow simulation of the electrolytic cell, thereby improving the efficiency of the gas-liquid two-phase flow simulation and reducing the time for gas-liquid two-phase flow simulation calculation of the entire electrolytic cell. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] The drawings described herein are used to provide a further understanding of the present invention and constitute a part of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0054] Figure 1 A schematic flow chart of a calculation method for two-phase flow simulation of an electrolytic cell provided by the present invention;

[0055] Figure 2 A schematic flow chart of another calculation method for two-phase flow simulation of an electrolytic cell provided by the present invention;

[0056] Figure 3 A schematic structural diagram of an alkaline electrolytic cell provided by the present invention;

[0057] Figure 4 A schematic diagram of the fitting of the gas-liquid two-phase flow simulation and the liquid phase flow simulation of the oxygen and electrolyte discharge section in an electrolytic cell provided by the present invention;

[0058] Figure 5 A schematic diagram showing the change in pressure of oxygen and electrolyte two-phase flow with position in a single-phase fitting of an electrolyte provided by the present invention;

[0059] Figure 6 A pressure distribution cloud diagram of two-phase flow in a single electrolysis chamber of an anode in an electrolytic cell provided by the present invention and a pressure distribution cloud diagram of the electrolysis chamber after fitting the gas-liquid two-phase flow;

[0060] Figure 7A schematic diagram of a computing device for simulating two-phase flow in an entire electrolytic cell provided by the present invention;

[0061] Figure 8 A schematic diagram of a computer device for realizing a calculation method for fitting two-phase flow simulation of an entire electrolytic cell provided by the present invention.

[0062] Description of reference numerals:

[0063] ①, hydrogen and electrolyte discharge section; ②, oxygen and electrolyte discharge section; ③, electrolyte inlet section; ④, anode electrolysis area; ⑤, cathode electrolysis area. DETAILED DESCRIPTION

[0064] To make the objectives, technical solutions, and advantages of the present invention more clear, the technical solutions of the present invention will be clearly and completely described below in conjunction with specific embodiments of the present invention and corresponding drawings. Obviously, the embodiments described are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0065] In recent years, China has increasingly prioritized renewable energy. With the introduction of the national "dual carbon" strategy, installed capacity for photovoltaic and wind power generation has increased annually. However, these methods are subject to volatility and face challenges in grid connection. Integrating new power systems with energy storage technologies is necessary to increase the grid utilization rate of renewable energy. Therefore, research in new energy and energy storage technologies is urgent. Hydrogen, as an excellent secondary energy source, offers advantages such as high calorific value and pollution-free reaction products. Furthermore, hydrogen is widely available, making it a key component of new energy storage systems. Among various hydrogen production technologies, water electrolysis is a mature and environmentally friendly technology. It is categorized into alkaline water electrolysis (ALK), proton exchange membrane (PEM), solid oxide electrolysis (SOEC), and anion exchange membrane (AEM) electrolysis. The current mature water electrolysis hydrogen production products mainly include proton exchange membrane water electrolysis hydrogen production and alkaline water electrolysis hydrogen production equipment. Solid oxide water electrolysis hydrogen production and anion exchange membrane water electrolysis hydrogen production are still mostly in the laboratory stage.

[0066] However, with the increasing demand for renewable energy consumption, the requirements for electrolyzer size and output are also increasing. However, the gas-liquid two-phase simulation results of large-scale electrolyzers are often difficult to converge, and model debugging is time-consuming.

[0067] Based on this, the present invention provides a calculation method for fitting two-phase flow simulation of the entire electrolytic cell. This method can solve the problems of slow simulation speed and difficult convergence of gas-liquid two-phase flow in the entire electrolytic cell, and provide guidance for the structural design of large electrolytic cells.

[0068] The technical solutions provided by various embodiments of the present invention are described in detail below with reference to the accompanying drawings.

[0069] Figure 1 The figure is a flow chart of a calculation method for fitting two-phase flow simulation of an electrolytic cell in the present invention, which specifically includes the following steps:

[0070] S101, meshing the electrolytic cell to obtain an electrolyte inlet section, a cathode and anode electrolysis areas, and a gas and electrolyte mixed discharge section of the electrolytic cell.

[0071] Among them, the gas and electrolyte mixed discharge section includes the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section; the anode and cathode electrolysis areas include the cathode electrolysis area and the anode electrolysis area; the electrolytic cell is meshed, including: according to the different flow states and functions of each part of the electrolytic cell, the electrolytic cell is modeled into the electrolyte inlet section, the cathode electrolysis area, the anode electrolysis area, the hydrogen and electrolyte mixed discharge section, and the oxygen and electrolyte mixed discharge section; the modeled electrolyte inlet section, cathode electrolysis area, anode electrolysis area, hydrogen and electrolyte mixed discharge section, and oxygen and electrolyte mixed discharge section are meshed for fluid calculation respectively.

[0072] The electrolyte inlet section includes the channels and flow channels leading to the cathode and anode electrolysis regions. The cathode and anode electrolysis regions are each composed of multiple electrolysis chambers, each of which contains electrodes, support structures, catalyst layers, and diaphragms. The hydrogen and electrolyte mixed discharge section includes all exhaust holes extending from the cathode electrolysis region and the total flow channel for hydrogen and electrolyte from all electrolysis chambers. The oxygen and electrolyte mixed discharge section includes all exhaust holes extending from the anode electrolysis region and the total flow channel for oxygen and electrolyte from all electrolysis chambers.

[0073] Specifically, the electrolytic cell is three-dimensionally modeled and divided into an electrolyte inlet section, a cathode electrolysis area, an anode electrolysis area, a hydrogen and electrolyte mixed discharge section, and an oxygen and electrolyte mixed discharge section according to the functional differences and different flow states of different parts of the electrolytic cell; then, the modeled electrolyte inlet section, cathode electrolysis area, anode electrolysis area, hydrogen and electrolyte mixed discharge section, oxygen and electrolyte mixed discharge section and other parts are respectively divided into fluid calculation grids to obtain the divided electrolyte inlet section, cathode electrolysis area, anode electrolysis area, hydrogen and electrolyte mixed discharge section, and oxygen and electrolyte mixed discharge section.

[0074] S102, performing gas-liquid two-phase flow simulation on the gas and electrolyte mixed discharge section and the single electrolysis chamber in the anode and cathode electrolysis area, respectively, to obtain a pressure variation curve of the gas and electrolyte mixed discharge section with position, and a two-phase pressure distribution cloud diagram of the single electrolysis chamber in the anode and cathode electrolysis area.

[0075] In one embodiment, a gas-liquid two-phase flow simulation is performed on the gas and electrolyte mixed discharge section to obtain a pressure variation curve of the gas and electrolyte mixed discharge section with position, including: performing a gas-liquid two-phase flow simulation on the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section, and extracting the pressure values ​​of the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section at different points in the flow direction; according to the position of each point and the corresponding pressure value, the pressure variation curve of the hydrogen and electrolyte mixed discharge section with position and the pressure variation curve of the oxygen and electrolyte mixed discharge section with position are drawn respectively.

[0076] In one embodiment, a gas-liquid two-phase flow simulation is performed on a single electrolysis chamber in the anode and cathode electrolysis regions to obtain a two-phase pressure distribution cloud map of the single electrolysis chamber in the anode and cathode electrolysis regions, including: a gas-liquid two-phase flow simulation is performed on a single electrolysis chamber in the cathode electrolysis region and a single electrolysis chamber in the anode electrolysis region to obtain a two-phase pressure distribution cloud map of the single electrolysis chamber in the cathode electrolysis region and a two-phase pressure distribution cloud map of the single electrolysis chamber in the anode electrolysis region.

[0077] Among them, the two-phase flow control equations in the gas-liquid two-phase flow simulation include the gas-liquid two-phase flow mass conservation equation and the gas-liquid two-phase flow momentum conservation equation.

[0078] The mass conservation equation for gas-liquid two-phase flow is:

[0079]

[0080]

[0081] The momentum conservation equation for gas-liquid two-phase flow is:

[0082]

[0083] Where, l is the liquid phase, g is the gas phase; t is the time, s is the volume fraction, ρ is the fluid density, u is the fluid velocity, and P is the pressure. is the stress-strain tensor, S m is the fluid mass source term, S u is the fluid momentum source term.

[0084] Therefore, based on formulas (1)-(4), gas-liquid two-phase flow simulation is performed on the hydrogen and electrolyte mixed discharge section, oxygen and electrolyte mixed discharge section, single electrolysis chamber in the cathode electrolysis area, and single electrolysis chamber in the anode electrolysis area. Then, according to different positions in the flow direction of the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section, the pressure value at each position is extracted, and the pressure change curve with position is drawn according to the position information and pressure data. At the same time, the pressure distribution cloud map in the single electrolysis chamber in the cathode electrolysis area and the single electrolysis chamber in the anode electrolysis area is extracted.

[0085] S103, in the process of liquid phase flow simulation of the gas and electrolyte mixed discharge section, by applying flow resistance in the liquid phase flow region, fitting the pressure change curve with position, and determining the value of the flow resistance of the gas and electrolyte mixed discharge section; and in the process of liquid phase flow simulation of a single electrolysis chamber in the anode and cathode electrolysis regions, respectively applying flow resistance in the cathode electrolysis region and the anode electrolysis region until the single-phase pressure distribution cloud diagram and the two-phase pressure distribution cloud diagram are consistent. Figure 1 The flow resistance values ​​of the cathode electrolysis area and the anode electrolysis area are determined.

[0086] Liquid-phase flow simulations were performed for the hydrogen and electrolyte mixing discharge sections and the oxygen and electrolyte mixing discharge sections. By adjusting the flow resistance applied to the liquid-phase flow region, the pressure variation curves for the liquid-phase flow simulations in these sections were fitted to the pressure variation curves for the two-phase flow. The liquid-phase flow region refers to the flow region within the hydrogen and electrolyte mixing discharge sections and the oxygen and electrolyte mixing discharge sections.

[0087] Based on the same principle mentioned above, liquid phase flow simulation is performed on a single electrolysis chamber in the cathode electrolysis area and a single electrolysis chamber in the anode electrolysis area. By applying flow resistance in the cathode electrolysis area and the anode electrolysis area respectively, the same method mentioned above is used to fit the pressure distribution cloud maps in the single electrolysis chamber in the cathode electrolysis area and the single electrolysis chamber in the anode electrolysis area to the two-phase pressure distribution cloud maps during two-phase flow. The flow resistance in the cathode electrolysis area and the flow resistance in the anode electrolysis area corresponding to the fitting of the two-phase pressure distribution cloud maps during two-phase flow are determined as the values ​​of the flow resistance in the cathode electrolysis area and the anode electrolysis area respectively.

[0088] Among them, applying flow resistance in the liquid phase flow region represents adjusting the fluid momentum source term; the fluid momentum source term S can be adjusted u The flow resistance is adjusted in this way, thereby affecting the single-phase flow pressure variation curve with position and the pressure distribution cloud diagram.

[0089] The liquid phase flow control equations in liquid phase flow simulation include the liquid phase mass conservation equation and the liquid phase momentum conservation equation. The liquid phase mass conservation equation is:

[0090]

[0091] The liquid phase momentum conservation equation is:

[0092]

[0093] In an exemplary embodiment, the fluid momentum source term S u The liquid phase flow area can be set to the porous medium area by adjusting the viscous resistance coefficient and the inertial resistance coefficient. ul Parameter size; if the liquid phase flow area is a porous medium area, the liquid phase flow control equations in the liquid phase flow simulation include the porous medium liquid phase flow mass conservation equation and the porous medium liquid phase flow momentum conservation equation.

[0094] The mass conservation equation for liquid flow in porous media is:

[0095]

[0096] Momentum conservation equation for liquid flow in porous media:

[0097]

[0098] S ul =-(Aμu i +Bρ|u|u i ) (9)

[0099] Where ε is the porosity, μ is the fluid dynamic viscosity, A is the viscous drag coefficient, and B is the inertial drag coefficient.

[0100] It should be noted that the gas-liquid two-phase flow simulation and the liquid phase flow simulation can be performed on the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section respectively, or the gas-liquid two-phase flow simulation and the liquid phase flow simulation can be performed on the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section at the same time; similarly, the gas-liquid two-phase flow simulation and the liquid phase flow simulation can be performed on the cathode electrolysis area and the anode electrolysis area respectively, or the gas-liquid two-phase flow simulation and the liquid phase flow simulation can be performed on the cathode electrolysis area and the anode electrolysis area at the same time.

[0101] S104: performing a liquid phase flow simulation on the entire electrolytic cell according to the values ​​of the flow resistances, and determining a fitted gas-liquid two-phase flow simulation result of the electrolytic cell.

[0102] Among them, the fitting gas-liquid two-phase flow simulation results include the flow and heat transfer parameter values ​​of each electrolytic chamber in the electrolytic cell; according to the values ​​of each flow resistance, the liquid phase flow simulation of the entire electrolytic cell is performed to determine the gas-liquid two-phase flow simulation results of the electrolytic cell, including: geometrically and meshing the electrolyte inlet section, the anode and cathode electrolysis areas, and the gas and electrolyte mixed discharge section to form the entire electrolytic cell; setting the flow resistance values ​​of each part to the corresponding parts of the entire electrolytic cell, performing liquid phase flow simulation on the entire electrolytic cell, and obtaining the flow and heat transfer parameter values ​​in each electrolytic chamber.

[0103] Specifically, the geometry and grid of each electrolytic cell part are combined to form the complete electrolytic cell geometry and grid. Then, the value of the flow resistance corresponding to each part is set to the corresponding part. The liquid phase flow simulation of the entire cell is performed to obtain the flow and heat transfer parameters in each electrolytic chamber, completing the two-phase fitting process of the entire cell.

[0104] Combining the partial geometry and grid of each electrolytic cell to form the geometry and grid of the complete electrolytic cell includes: combining the geometry and grid of the hydrogen and electrolyte mixed discharge section, the oxygen and electrolyte mixed discharge section, the cathode electrolysis area, the anode electrolysis area, and the electrolyte inlet section into a complete electrolytic cell.

[0105] The value of the flow resistance includes the resistance coefficient. The recorded resistance coefficient of the hydrogen and electrolyte mixed discharge section and the resistance coefficient of the oxygen and electrolyte mixed discharge section are applied to the corresponding areas of the corresponding complete electrolytic cell, and the resistance coefficient of a single electrolysis chamber extracted from the cathode and anode electrolysis areas are applied to all cathode and anode electrolysis areas. After the corresponding resistance coefficients are applied to the complete electrolytic cell, liquid phase flow simulation is performed to obtain the corresponding flow and heat transfer data, completing the fitting of the two-phase process of the entire cell.

[0106] In an exemplary embodiment, Figure 2 As shown, the present invention also provides a calculation method for fitting two-phase flow simulation of the entire electrolytic cell, where the electrolytic cell is an alkaline electrolytic cell; this embodiment includes the following steps:

[0107] S201, three-dimensional modeling of the electrolytic cell is performed, and the electrolytic cell is divided into an electrolyte inlet section, a cathode electrolysis area, an anode electrolysis area, a hydrogen and electrolyte discharge section, and an oxygen and electrolyte discharge section.

[0108] like Figure 3 As shown, Figure 3 The components of the alkaline electrolyzer are shown; among them, ① is the hydrogen and electrolyte discharge section, ② is the oxygen and electrolyte discharge section, ③ is the electrolyte inlet section, ④ is the anode electrolysis area, and ⑤ is the cathode electrolysis area.

[0109] S202, meshing each part.

[0110] S203 , performing gas-liquid two-phase flow simulation on the hydrogen and electrolyte discharge section and the oxygen and electrolyte discharge section, and extracting the pressure variation curves of the two discharge section structures with respect to position.

[0111] Among them, such as Figure 4 As shown, Figure 4 The schematic diagram of the gas-liquid two-phase flow simulation and liquid phase flow simulation fitting in the oxygen and electrolyte discharge section is shown; Figure 4 The upper middle part shows the two-phase flow of oxygen and electrolyte in the oxygen and electrolyte discharge section, and the lower middle part shows the flow of electrolyte in the discharge section under the application of flow resistance.

[0112] S204 , performing liquid phase flow simulation for the hydrogen and electrolyte discharge section and the oxygen and electrolyte discharge section, applying flow resistance in the flow area, fitting the two-phase flow pressure change curve, and recording the flow resistance parameters of each section.

[0113] Among them, liquid phase flow simulation is performed on the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section. In this embodiment, the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section are set as porous medium models, and the momentum source term brought by the single-phase flow resistance coefficient is adjusted by adjusting the viscous resistance coefficient and the inertial resistance coefficient, so that the pressure distribution curve with position is close to that of the two phases; the flow resistance parameters of each part include the viscous resistance coefficient A1 and the inertial resistance coefficient B1 of the oxygen and electrolyte discharge section, and the viscous resistance coefficient A2 and the inertial resistance coefficient B2 of the hydrogen and electrolyte discharge section. Figure 5 As shown, Figure 5 A schematic diagram showing the change in pressure of the electrolyte single-phase fitting oxygen and electrolyte two-phase flow with position is shown.

[0114] S205, extracting a single electrolysis chamber in the cathode electrolysis area and the anode electrolysis area to perform gas-liquid two-phase flow simulation, and extracting a pressure distribution cloud map in the electrolysis chamber.

[0115] S206, performing liquid phase flow simulation for a single electrolysis chamber in the cathode electrolysis area and the anode electrolysis area, applying flow resistance in the flow area, fitting the pressure cloud map in the two-phase flow electrolysis chamber, and recording the flow resistance parameters of each part.

[0116] Among them, the liquid phase flow simulation is carried out for the electrolysis chamber, and the flow area is set as a porous medium. The momentum source term brought by the single-phase flow resistance coefficient is adjusted by adjusting the viscous resistance coefficient and the inertial resistance coefficient; the flow resistance parameters of each part include the viscous resistance coefficient A3 and the inertial resistance coefficient B3 in the cathode electrolysis chamber, and the viscous resistance coefficient A4 and the inertial resistance coefficient B4 in the anode electrolysis chamber. Figure 6 As shown, Figure 6 The pressure distribution cloud diagram of the two-phase flow in a single electrolysis chamber at the anode and the pressure distribution cloud diagram of the electrolysis chamber after fitting the gas-liquid two-phase flow are shown.

[0117] S207, combining the geometries and meshes of the various parts of the electrolytic cell to form a complete electrolytic cell geometry and mesh.

[0118] S208, setting the recorded flow resistance of each part into the corresponding position, completing the liquid phase flow simulation of the entire tank, and obtaining the flow and heat transfer data in the electrolytic cell.

[0119] Specifically, the various components are combined into a complete electrolyzer, and the hydrogen and electrolyte mixed discharge section, oxygen and electrolyte mixed discharge section, cathode electrolysis area, and anode electrolysis area of ​​the complete electrolyzer are configured as porous media. A1 and B1 are applied to the hydrogen and electrolyte mixed discharge section; A2 and B2 are applied to the oxygen and electrolyte mixed discharge section; A3 and B3 are applied to the cathode electrolysis area; and A4 and B4 are applied to the anode electrolysis area. Electrolyte liquid-phase flow simulation is performed on the entire cell, and flow and heat transfer data within each electrolysis cell are obtained, completing the fitting of the electrolyte single-phase flow simulation to gas and electrolyte gas-liquid two-phase flow simulation.

[0120] In the present invention, when it is necessary to simulate the two-phase flow of the electrolytic cell, the unidirectional flow simulation result can be made consistent with the two-phase flow simulation result by fitting the unidirectional flow simulation, and the flow resistance applied to each part of the electrolytic cell during the unidirectional flow simulation can be extracted. In this way, by setting the resistance of each part to each part of the electrolytic cell, the two-phase flow simulation of the electrolytic cell can be fitted by performing unidirectional flow simulation on the entire electrolytic cell, that is, the result of the unidirectional flow simulation of the electrolytic cell is the fitted two-phase flow simulation result of the electrolytic cell.

[0121] When applying the calculation method for the electrolytic cell two-phase flow simulation provided by the present invention, it is not necessary to Figure 1 The steps are executed in the order shown. The specific execution order of the steps can be determined according to needs, and the present invention does not limit this.

[0122] The above is a calculation method for the whole electrolytic cell fitting two-phase flow simulation provided by one or more embodiments of the present invention. Based on the same idea, the present invention also provides a corresponding calculation device for the whole electrolytic cell fitting two-phase flow simulation, such as Figure 7 shown.

[0123] Figure 7 A schematic diagram of a computing device for simulating two-phase flow in an entire electrolytic cell provided by the present invention is shown. The device 700 includes:

[0124] A division module 701 is used to divide the electrolytic cell into a grid to obtain the electrolyte inlet section, the cathode and anode electrolysis areas, and the gas and electrolyte mixed discharge section of the electrolytic cell;

[0125] The first simulation module 702 is used to perform gas-liquid two-phase flow simulation for the gas and electrolyte mixed discharge section and the single electrolysis chamber in the anode and cathode electrolysis area, respectively, to obtain a pressure variation curve of the gas and electrolyte mixed discharge section with position, and a two-phase pressure distribution cloud diagram of the single electrolysis chamber in the anode and cathode electrolysis area;

[0126] The second simulation module 703 is used to determine the value of the flow resistance of the gas and electrolyte mixed discharge section by fitting the pressure-position variation curve by applying flow resistance in the liquid phase flow region during the liquid phase flow simulation of the gas and electrolyte mixed discharge section;

[0127] The third simulation module 704 is used to simulate the liquid phase flow in a single electrolysis chamber in the cathode electrolysis area by applying flow resistance in the cathode electrolysis area and the anode electrolysis area respectively until the single-phase pressure distribution cloud and the two-phase pressure distribution cloud are consistent. Figure 1 To determine the flow resistance values ​​of the cathode electrolysis area and the anode electrolysis area;

[0128] The fourth simulation module 705 is used to perform liquid phase flow simulation on the entire electrolytic cell according to the values ​​of the flow resistance of each part, and determine the fitted gas-liquid two-phase flow simulation results of the electrolytic cell.

[0129] Regarding the specific limitations of the computing device for fitting the two-phase flow simulation of the entire electrolytic cell, please refer to the limitations of the computing method for fitting the two-phase flow simulation of the entire electrolytic cell above, which will not be repeated here. The various modules in the computing device for fitting the two-phase flow simulation of the entire electrolytic cell can be implemented in whole or in part by software, hardware, and a combination thereof. The above modules can be embedded in or independent of the processor in the computer device in hardware form, or can be stored in the memory of the computer device in software form, so that the processor can call and execute the operations corresponding to the above modules.

[0130] The present invention also provides a computer-readable storage medium, which stores a computer program, which can be used to execute the above Figure 1 The provided calculation method for the two-phase flow simulation of the entire electrolytic cell is used.

[0131] The present invention also provides Figure 8The structural diagram of the computer equipment shown in FIG. Figure 8 As shown in the figure, at the hardware level, the computer device includes a processor, an internal bus, a network interface, a memory, and a non-volatile memory. Of course, it may also include other hardware required for the business. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs it to achieve the above Figure 1 The provided calculation method for the two-phase flow simulation of the entire electrolytic cell is used.

[0132] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided by the present invention can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory or optical memory, etc. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0133] The technical features of the above embodiments can be combined arbitrarily. In order to make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of the present invention.

Claims

1. A calculation method for two-phase flow simulation of an electrolytic cell, characterized in that: include: Dividing the electrolytic cell into grids to obtain an electrolyte inlet section, a cathode and anode electrolysis regions, and a gas and electrolyte mixed discharge section of the electrolytic cell; Performing gas-liquid two-phase flow simulations on the gas and electrolyte mixed discharge section and a single electrolysis chamber in the anode and cathode electrolysis regions, respectively, to obtain a pressure variation curve of the gas and electrolyte mixed discharge section with respect to position, and a two-phase pressure distribution cloud diagram of a single electrolysis chamber in the anode and cathode electrolysis regions; During the liquid phase flow simulation of the gas and electrolyte mixed discharge section, the flow resistance of the gas and electrolyte mixed discharge section is determined by fitting the pressure-position variation curve by applying flow resistance in the liquid phase flow region; During the liquid phase flow simulation of a single electrolysis chamber in the cathode electrolysis region, flow resistances are applied in the cathode electrolysis region and the anode electrolysis region respectively until the single-phase pressure distribution cloud map is consistent with the two-phase pressure distribution cloud map, thereby determining the flow resistance values ​​of the cathode electrolysis region and the anode electrolysis region; According to the values ​​of the flow resistances, liquid phase flow simulation is performed on the entire electrolytic cell to determine the fitted gas-liquid two-phase flow simulation results of the electrolytic cell; The gas-liquid two-phase flow control equations in the gas-liquid two-phase flow simulation include the gas-liquid two-phase flow mass conservation equation and the gas-liquid two-phase flow momentum conservation equation; The mass conservation equation for the gas-liquid two-phase flow is: The momentum conservation equation for the gas-liquid two-phase flow is: in, The subscript is liquid phase, The subscripts are gas phase; t For time, s is the volume fraction, ρ is the fluid density, u is the fluid velocity, is the pressure, is the stress-strain tensor, S m is the fluid mass source term, S u is the fluid momentum source term; The liquid phase flow control equations in the liquid phase flow simulation include the liquid phase mass conservation equation and the liquid phase momentum conservation equation; The liquid phase mass conservation equation: The liquid phase momentum conservation equation: Wherein, applying flow resistance in the liquid phase flow region represents adjusting the fluid momentum source term; The liquid phase flow region is a porous medium region; the liquid phase flow control equations in the liquid phase flow simulation include a porous medium liquid phase flow mass conservation equation and a porous medium liquid phase flow momentum conservation equation; The mass conservation equation for the porous medium liquid phase flow is: Momentum conservation equation for liquid flow in porous media: Wherein, ε is the porosity, μ is the fluid dynamic viscosity, A is the viscous drag coefficient, and B is the inertial drag coefficient; the fluid momentum source term in the liquid phase is adjusted by the viscous drag coefficient and the inertial drag coefficient.

2. The method according to claim 1, characterized in that The gas and electrolyte mixed discharge section includes a hydrogen and electrolyte mixed discharge section and an oxygen and electrolyte mixed discharge section; the grid division of the electrolytic cell includes: According to the different flow states and functions of various parts in the electrolytic cell, the electrolytic cell is modeled into an electrolyte inlet section, a cathode electrolysis area, an anode electrolysis area, a hydrogen and electrolyte mixed discharge section, and an oxygen and electrolyte mixed discharge section; The modeled electrolyte inlet section, the cathode electrolysis area, the anode electrolysis area, the hydrogen and electrolyte mixed discharge section, and the oxygen and electrolyte mixed discharge section are respectively divided into fluid calculation grids.

3. The method according to claim 2, characterized in that A gas-liquid two-phase flow simulation is performed on the gas and electrolyte mixed discharge section to obtain a pressure variation curve of the gas and electrolyte mixed discharge section with position, including: Performing gas-liquid two-phase flow simulation on the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section, and extracting pressure values ​​at different points in the flow direction of the hydrogen and electrolyte mixed discharge section and the oxygen and electrolyte mixed discharge section respectively; According to the positions of each point and the corresponding pressure value, a pressure variation curve of the hydrogen and electrolyte mixed discharge section and a pressure variation curve of the oxygen and electrolyte mixed discharge section are drawn respectively.

4. The method according to claim 2, characterized in that Performing a gas-liquid two-phase flow simulation on a single electrolysis chamber in the anode and cathode electrolysis region to obtain a two-phase pressure distribution cloud map of the single electrolysis chamber in the anode and cathode electrolysis region includes: Gas-liquid two-phase flow simulation is performed on a single electrolysis chamber in the cathode electrolysis area and a single electrolysis chamber in the anode electrolysis area to obtain a two-phase pressure distribution cloud map of the single electrolysis chamber in the cathode electrolysis area and a two-phase pressure distribution cloud map of the single electrolysis chamber in the anode electrolysis area.

5. The method according to claim 1, wherein The fitted gas-liquid two-phase flow simulation result includes the flow and heat exchange parameter values ​​of each electrolytic chamber in the electrolytic cell; performing liquid phase flow simulation on the entire electrolytic cell according to the values ​​of each flow resistance to determine the fitted gas-liquid two-phase flow simulation result of the electrolytic cell, including: The electrolyte inlet section, the cathode and anode electrolysis areas, and the gas and electrolyte mixed discharge section are geometrically and grid-combined to form an electrolytic cell; The flow resistance values ​​of each part are set to the corresponding parts of the entire electrolytic cell, and the liquid phase flow simulation of the entire electrolytic cell is performed to obtain the flow and heat exchange parameter values ​​in each electrolytic chamber.

6. A computing device for two-phase flow simulation of an electrolytic cell, characterized in that: include: A partitioning module is used to divide the electrolytic cell into grids to obtain the electrolyte inlet section, the cathode and anode electrolysis areas, and the gas and electrolyte mixed discharge section of the electrolytic cell; A first simulation module is used to perform gas-liquid two-phase flow simulation on the gas and electrolyte mixed discharge section and a single electrolysis chamber in the anode and cathode electrolysis area, respectively, to obtain a pressure variation curve of the gas and electrolyte mixed discharge section with position, and a two-phase pressure distribution cloud diagram of a single electrolysis chamber in the anode and cathode electrolysis area; the gas-liquid two-phase flow control equations in the gas-liquid two-phase flow simulation include a gas-liquid two-phase flow mass conservation equation and a gas-liquid two-phase flow momentum conservation equation; The mass conservation equation for the gas-liquid two-phase flow is: The momentum conservation equation for the gas-liquid two-phase flow is: in, The subscript is liquid phase, The subscripts are gas phase; t For time, s is the volume fraction, ρ is the fluid density, u is the fluid velocity, is the pressure, is the stress-strain tensor, S m is the fluid mass source term, S u is the fluid momentum source term; a second simulation module for determining the value of the flow resistance of the gas and electrolyte mixed discharge section by fitting the pressure-versus-position curve during a liquid phase flow simulation of the gas and electrolyte mixed discharge section by applying a flow resistance in the liquid phase flow region; the liquid phase flow control equations in the liquid phase flow simulation include a liquid phase mass conservation equation and a liquid phase momentum conservation equation; The liquid phase mass conservation equation: The liquid phase momentum conservation equation: Wherein, applying flow resistance in the liquid phase flow region represents adjusting the fluid momentum source term; The liquid phase flow region is a porous medium region; the liquid phase flow control equations in the liquid phase flow simulation include a porous medium liquid phase flow mass conservation equation and a porous medium liquid phase flow momentum conservation equation; The mass conservation equation for the porous medium liquid phase flow is: Momentum conservation equation for liquid flow in porous media: Wherein, ε is the porosity, μ is the fluid dynamic viscosity, A is the viscous drag coefficient, and B is the inertial drag coefficient; the fluid momentum source term in the liquid phase is adjusted by the viscous drag coefficient and the inertial drag coefficient; a third simulation module for determining the values ​​of the flow resistances in the cathode electrolysis region and the anode electrolysis region by applying flow resistances in the cathode electrolysis region and the anode electrolysis region respectively during a liquid phase flow simulation of a single electrolysis chamber in the anode electrolysis region until the single-phase pressure distribution cloud map is consistent with the two-phase pressure distribution cloud map; The fourth simulation module is used to perform liquid phase flow simulation on the entire electrolytic cell according to the values ​​of each flow resistance, and determine the fitted gas-liquid two-phase flow simulation result of the electrolytic cell.

7. A computer-readable storage medium, characterized in that The storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 5 is implemented.

Images