A Multiphysics Coupling Analysis Method and System for Industrial-Grade Water Electrolysis for Hydrogen Production

Abstract

A multiphysics coupling analysis method for industrial-grade water electrolysis hydrogen production is presented. This method first identifies the key physical fields affecting the electro-hydrogen energy conversion efficiency in the industrial-grade water electrolysis hydrogen production device. Then, based on the mechanical structure and electrochemical reaction mechanism of the industrial-grade water electrolysis hydrogen production device, it establishes the governing equations for each key physical field. Next, it combines the governing equations of each key physical field with the physical processes to construct a multiphysics coupling model. Finally, it solves the multiphysics coupling model to obtain the field distribution of each key physical field. This invention, based on a multiphysics model of a compact structure of an industrial-grade water electrolysis hydrogen production device, explores the coupling mechanism between various physical fields, thereby guiding the design and operation of industrial-grade water electrolysis hydrogen production devices.

Description

Technical Field

[0001] This invention belongs to the field of new energy technology, specifically relating to an industrial-grade multi-physics field coupling analysis method and system for hydrogen production by water electrolysis. Background Technology

[0002] Hydrogen energy boasts numerous advantages, including high energy density, excellent combustion performance, diverse utilization pathways, zero carbon emissions, and minimal environmental pollution, making it an ideal secondary energy source. Renewable energy-based water electrolysis for hydrogen production is an effective pathway for converting clean electricity into green hydrogen energy, playing a central role in the green hydrogen energy industry.

[0003] With the vigorous development of new energy power construction, efficient hydrogen production from water electrolysis devices is particularly important. Current research on the characteristics of electro-hydrogen production systems rarely considers the aforementioned compact structure, making it difficult to apply to industrial-scale devices; furthermore, the mechanisms by which multiple physical fields affect reaction performance are not yet clear, making it impossible to quantitatively guide the design of operating parameters. Summary of the Invention

[0004] The purpose of this invention is to address the aforementioned problems in the existing technology by providing an industrial-grade multiphysics coupling analysis method and system for hydrogen production through water electrolysis.

[0005] To achieve the above objectives, the technical solution of the present invention is as follows:

[0006] A multiphysics coupling analysis method for industrial-grade water electrolysis hydrogen production includes the following steps:

[0007] Step A: Identify the key physical fields affecting the energy conversion efficiency of electricity to hydrogen in an industrial-grade water electrolysis hydrogen production unit;

[0008] Step B: Establish the control equations for each key physical field based on the mechanical structure and electrochemical reaction mechanism of the industrial-grade water electrolysis hydrogen production device;

[0009] Step C: Construct a multiphysics coupling model by combining the governing equations of each key physical field with the physical processes;

[0010] Step D: Solve the multiphysics coupling model to obtain the field distribution of each key physics field.

[0011] The key physical fields include the current density field, flow field, temperature field, and substance concentration field.

[0012] In step B, the governing equation for the current density field is:

[0013]

[0014] In the above formula, This is an electrode polarization overvoltage. For gases, the universal constant is... For temperature, It is Faraday's constant. It is a hyperbolic sine function. , These are current density and exchange current density, respectively.

[0015] The governing equations for the flow field are:

[0016]

[0017] In the above formula, Let be the fluid density, and t be the time. For gradient operators, It is a velocity vector. For pressure, For fluid dynamic viscosity, It is the identity matrix. It is the gravitational acceleration vector. This is the surface tension term;

[0018] The governing equation for the temperature field is:

[0019]

[0020] In the above formula, For constant pressure heat capacity, Thermal conductivity, It serves as the heat source for the temperature field;

[0021] The governing equation for the concentration field of the substance is:

[0022]

[0023] In the above formula, Let be the flux of substance i. Let be the diffusion coefficient of substance i. , Here, the concentration and charge number of substance i are respectively. For electromobility, The potential of the electrolyte. For diffusion and electromigration rates, This refers to the electrolyte flow rate.

[0024] In step C, the multiphysics coupling model includes:

[0025]

[0026]

[0027]

[0028]

[0029]

[0030]

[0031]

[0032] In the above formula, The rate at which the substance is deposited from the electrode surface. The current density at the electrode surface is... It is Faraday's constant. The number of charges on a substance. , These are the electrode equilibrium potential and the standard electrode equilibrium potential, respectively. , These represent the concentrations of oxides and products, respectively. For standard concentration, For constant pressure heat capacity, Thermal conductivity, As the heat source of the temperature field, For temperature, Let be the flux of substance i. Let be the diffusion coefficient of substance i. , Here, the concentration and charge number of substance i are respectively. For electromobility, The potential of the electrolyte. For diffusion and electromigration rates, Electrolyte flow rate , These are Joule heat sources and electrochemical heat sources, respectively. The electrode active surface area is denoted by m, which represents the cathode or anode. , These are the electrode current density vector and the electrolyte current density vector, respectively. , These are the electrode potential and the electrolyte potential, respectively. The electrode potential is... This represents the current density of the local double layer of the electrode.

[0033] Step D uses COMSOL Multiphysics software to solve the multiphysics coupling model.

[0034] An industrial-grade multiphysics coupling analysis system for hydrogen production by water electrolysis includes a key physics field determination module, a governing equation establishment module, a multiphysics coupling model construction module, and a coupling model solution module.

[0035] The key physical field determination module is used to determine the key physical fields that affect the energy conversion efficiency of electricity-hydrogen in industrial-grade water electrolysis hydrogen production devices.

[0036] The control equation establishment module is used to establish control equations for each key physical field based on the mechanical structure and electrochemical reaction mechanism of the industrial-grade water electrolysis hydrogen production device.

[0037] The multiphysics coupling model construction module is used to construct a multiphysics coupling model by combining the control equations of each key physical field and the physical process.

[0038] The coupling model solving module is used to solve the multiphysics coupling model and obtain the field distribution of each key physical field.

[0039] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0040] This invention discloses a multi-physics coupling analysis method for industrial-grade water electrolysis hydrogen production. First, it identifies the key physical fields affecting the electro-hydrogen energy conversion efficiency in the industrial-grade water electrolysis hydrogen production device. Then, based on the mechanical structure and electrochemical reaction mechanism of the industrial-grade water electrolysis hydrogen production device, it establishes the governing equations for each key physical field. Next, it combines the governing equations of each key physical field with the physical processes to construct a multi-physics coupling model. Finally, it solves the multi-physics coupling model to obtain the field distribution of each key physical field. This method, based on a multi-physics model under the compact structure of an industrial-grade water electrolysis hydrogen production device, explores the coupling mechanism between various physical fields, thereby guiding the design and operation of industrial-grade water electrolysis hydrogen production devices. Attached Figure Description

[0041] Figure 1 This is a flowchart of the present invention.

[0042] Figure 2 This diagram illustrates the coupling relationships between key physical fields.

[0043] Figure 3 This is a geometric diagram of the electrolytic cell described in Example 1.

[0044] Figure 4 This is a diagram showing the electrolyte flow rate distribution.

[0045] Figure 5 This is a diagram showing the hydrogen ion concentration distribution in the cathode electrolyte under steady-state conditions.

[0046] Figure 6 This is a graph showing the hydrogen ion concentration distribution in the anolyte under steady-state conditions.

[0047] Figure 7 This is a diagram showing the hydrogen concentration distribution at the cathode under steady-state conditions.

[0048] Figure 8This is a vector diagram showing the distribution of electrolyte potential and current density in the electrolyte.

[0049] Figure 9 This is a diagram showing the temperature and isothermal surface distribution during stable operation of the electrolytic cell.

[0050] Figure 10 This is a system structure diagram as described in Example 2. Detailed Implementation

[0051] The present invention will now be described in further detail with reference to specific embodiments and accompanying drawings.

[0052] Example 1:

[0053] This embodiment uses Figure 3 The industrial-grade water electrolysis electrolyzer for hydrogen production shown is the subject of analysis. This electrolyzer includes:

[0054] Connecting conductors: Composed of highly conductive metallic conductors, they serve to connect and conduct electricity when multiple electrolytic cells are connected in series or parallel.

[0055] Electrolyte channel: The electrolyte flows into the electrolytic cell through this channel, and the reaction products also flow out through the outlet through this channel.

[0056] Porous current collector: A current collector composed of a porous conductive medium. Its porous structure allows the electrolyte to penetrate into the porous catalyst layer to participate in the electrode reaction, while also limiting the bubble radius of the gas-phase products.

[0057] Porous catalyst layer: Electrochemical reactions occur in the catalyst layer, reactants are consumed, and gaseous products are generated at the same time.

[0058] Membrane electrode assembly layer: The assembly structure of the membrane and electrode. The membrane is used to prevent the products from the anode and cathode from mixing.

[0059] The geometric parameters of the electrolytic cell are shown in Table 1:

[0060] Table 1 Geometric parameters of the electrolytic cell

[0061]

[0062] See Figure 1 A multi-physics coupling analysis method for industrial-grade water electrolysis hydrogen production is performed according to the following steps:

[0063] 1. Identify the key physical fields affecting the electro-hydrogen energy conversion efficiency in an industrial-grade water electrolysis hydrogen production device. The key physical fields include the current density field, flow field, temperature field, and substance concentration field.

[0064] 2. Based on the mechanical structure and electrochemical reaction mechanism of the industrial-grade water electrolysis hydrogen production device, the governing equations for each key physical field are established, among which,

[0065] Current density field:

[0066] The objective of solving the electrochemical reaction process is to determine the current density vector and voltage magnitude at various points in the electrolyte. From the perspective of current formation principles, the directional movement of various charged particles in the electrolytic cell within the electric field generates current; therefore, the governing equations for the current density field include the mass conservation equation and the charge conservation equation. In an electrolytic cell at room temperature, the governing equations for the electrochemical reaction are:

[0067]

[0068]

[0069] In the above formula, It is the current density vector. The conductivity of the diaphragm. It represents the electric potential.

[0070] According to electrode kinetics, when two reversible electrodes are connected by a diaphragm or salt bridge but do not form a closed loop, a stable potential difference will be formed between the two equilibrium electrodes. This potential difference is the reversible voltage of the electrolytic cell. When the electrolytic cell is running, electrode polarization and resistance effects will lead to a decrease in electrolysis efficiency. The electrolysis voltage can be decomposed into:

[0071]

[0072] In the above formula, This is an electrode polarization overvoltage. This refers to the voltage loss caused by the equivalent resistance of the electrolytic cell.

[0073] Given by the Butler-Volmer formula:

[0074] .

[0075] Flow field:

[0076] The momentum transfer processes in an electrolyzer involve the flow of gas and liquid in free channels and porous media. The momentum equations for fluid flow are described by the Navier-Stokes equations:

[0077]

[0078] In the above formula, Let be the fluid density, and t be the time. For gradient operators, It is a velocity vector. For pressure, For fluid dynamic viscosity, It is the identity matrix. It is the gravitational acceleration vector. This is the surface tension term.

[0079] This equation is used to solve for the parameters of free flow of a single-phase fluid, including the velocity vector and pressure. It is used to describe the flow in the region of free flow.

[0080] Temperature field:

[0081] As the electrolysis process proceeds, the temperature of the electrolytic cell gradually increases due to factors such as electrochemical reactions and the equivalent resistance of the cell. This process is accompanied by energy transfer within the cell. There are three modes of heat transfer: conduction, convection, and radiation. Energy transfer due to a temperature gradient is called conduction; energy transfer under forced convection in a flow field is called convection; and energy transfer using photons as carriers is called radiation. This invention only considers the influence of conduction and convection within the electrolytic cell on the temperature distribution. The equation describing this energy transfer process is:

[0082]

[0083] In the above formula, For constant pressure heat capacity, Thermal conductivity, It is the heat source for the temperature field.

[0084] This formula reflects that the heat provided by the heat source to the system is consumed by heat conduction and heat convection, simultaneously causing an increase in the internal temperature of the system. In this study, the heat source term of the heat transfer module comes from the heat of electrochemical reaction and the Joule heat caused by the equivalent resistance of the electrolytic cell; while the dependent variable... This will cause changes in some parameters of other physical fields, such as fluid density. diffusion coefficient Electrode equilibrium potential wait.

[0085] Matter concentration field:

[0086] Mass transfer refers to the diffusion, convection, and electromigration of various substances during electrolysis. In alkaline or proton exchange membrane electrolyzers, the electrolyte often satisfies the dilute solution theory; therefore, the Nernst-Planck equation is used to describe the mass transfer process. The equation used to describe the flux of substance i in the electrolyte is the Nernst-Planck equation:

[0087]

[0088] In the above formula, Let be the flux of substance i. Let be the diffusion coefficient of substance i. , Here, the concentration and charge number of substance i are respectively. For electromobility, The potential of the electrolyte. This refers to the electrolyte flow rate.

[0089] The first term on the right-hand side of the equation represents diffusion, the second term represents electromigration of charged matter, and the third term represents convective transfer of matter. Generally, the intensity of convective transfer is much higher than that of diffusion, and the flow velocity in the convection term is the dependent variable in the flow field. Therefore, it is necessary to couple the flow field with the mass transfer process to accurately describe the concentration distribution of matter.

[0090] 3. Construct a multiphysics coupling model by combining the control equations of each key physical field and the physical processes.

[0091] The multiphysics processes within an electrolyzer are not independent but rather exhibit complex coupling relationships. For example, the electrode current density determines the reaction rate, thus affecting the rate of substance formation or consumption and its concentration distribution; conversely, changes in substance concentration also alter the electrode potential, thereby influencing the current density distribution. Furthermore, the fluid flow velocity distribution directly determines the intensity of convective mass transfer and convective heat transfer, while the substance concentration distribution and the electrolyzer's temperature distribution affect fluid properties such as density and dynamic viscosity. The multiphysics coupling relationships within the electrolyzer of this invention are as follows: Figure 2 As shown.

[0092] (1) Current density field - matter concentration field: Faraday's law

[0093]

[0094] In the above formula, The rate at which the substance is deposited from the electrode surface. The current density at the electrode surface is... It is Faraday's constant. The charge number of the substance.

[0095] Faraday's law clarifies the relationship between the flux of substances participating in electrochemical reactions on the electrode surface and the current density, thus establishing a one-way coupling between the current density field and the substance concentration field.

[0096] (2) Matter concentration field - current density field: Nernst equation

[0097]

[0098] In the above formula, , These are the electrode potential and the standard electrode potential, respectively. , These represent the concentrations of oxides and products, respectively. For standard concentration, =1 mol / L.

[0099] The concentration of a substance affects the electrode potential, which in turn affects the electrochemical reaction.

[0100] (3) Current density field-temperature field: electrochemical heat equation

[0101]

[0102] In the above formula, As the heat source of the temperature field, , These are Joule heat sources and electrochemical heat sources, respectively. denoted as the electrode active surface area, where m represents the cathode or anode.

[0103] Joule heat sources include the Joule heat of electrodes and electrolytes, and are defined as follows:

[0104]

[0105] In the above formula, , These are the electrode current density vector and the electrolyte current density vector, respectively. , These are the electrode potential and the electrolyte potential, respectively.

[0106] The electrochemical heat source is generated by electrode polarization overpotential and temperature rise overpotential, which are defined as:

[0107]

[0108] In the above formula, The electrode potential is... The current density of the local electric double layer of the electrode. This is the double-layer potential, which is related to the electrode potential. The difference is the overpotential caused by electrode polarization. The term represents the overpotential caused by temperature rise.

[0109] The coupling relationship between the electrochemical current density field and the temperature field can be established using the above three equations.

[0110] (4) Convective heat and mass transfer

[0111] As shown by the Nernst-Planck equations and the heat transfer equations, both mass transfer and energy transfer processes involve convection terms, meaning that fluid flow influences the concentration and temperature fields. For typical mass transfer processes, the convective mass transfer rate is much greater than that of diffusion; for typical heat transfer processes, the convective heat transfer efficiency is much higher than that of conduction. This indicates that the velocity distribution of the flow field has a significant impact on the concentration distribution of the substance and the temperature distribution of the electrolyzer. Therefore, in modeling, it is necessary to establish a single-phase coupling relationship between the flow field and the mass and heat transfer fields.

[0112] On the other hand, the concentration distribution of substances in the fluid and the distribution of fluid temperature will cause changes in parameters such as fluid density and dynamic viscosity, thereby affecting the flow. However, if the temperature and concentration variations are small within the entire geometric space of the electrolyzer, the changes in fluid dynamic parameters are minimal. Therefore, in this invention, this process is ignored, i.e., only single-phase coupling is established.

[0113] 4. The multiphysics coupling model was solved by simulation using COMSOL Multiphysics software to obtain the field distribution of each key physical field.

[0114] Solution results of the multiphysics coupling model:

[0115] The flow velocity distribution of the electrolyte in the electrolyte channel and porous medium in the electrolytic cell is as follows: Figure 4 As shown. The streamline distribution at the central interface of the electrolytic cell is as follows. Figure 4 As shown in (b), Figure 4(a) shows that in the region where the electrolyte flows freely, the fluid velocity in the center is the fastest, reaching 2.09 mm / s; while the velocity near the wall is 0. This conclusion is consistent with the flow characteristics of viscous fluids. Figure 4 (b) reflects the distribution of streamlines at the center of the electrolyzer. Electrolyte flowing into the electrolyzer from the free channel permeates into the porous medium through osmosis and other processes, reaching the catalyst layer to participate in the electrochemical reaction. Furthermore, due to gravity, a small portion of the anolyte permeates through the porous medium to the cathode, resulting in a denser streamline distribution in the cathode porous medium compared to the anode.

[0116] The hydrogen ion concentration distributions in the cathode and anolyte electrolytes under steady-state conditions are as follows: Figure 5 , Figure 6 As shown. The concentration streamline distributions caused by the convective mass transfer effect are respectively as follows: Figure 5 (b) Figure 6As shown in (d), simulation results show that after hydrogen ions in the cathode electrolyte flow into the electrolyte channel, they are consumed by convection and diffusion as they reach the porous catalyst layer and participate in the electrode reaction. This causes the hydrogen ion concentration in the cathode electrolyte to gradually decrease along the channel. At the anode, hydrogen ions are continuously generated through electrochemical reactions and flow into the anode electrolyte channel with the electrolyte flow and their own diffusion. Therefore, the hydrogen ion concentration in the anode electrolyte gradually increases. Furthermore, a comparison of hydrogen ion concentrations at the anode and cathode reveals that the increase in anode hydrogen ion concentration is numerically equal to the decrease in cathode hydrogen ion concentration, which is a necessary result of charge conservation during electrolysis.

[0117] The concentration distribution of hydrogen in the cathode electrolyte is as follows: Figure 7 As shown, as the reaction proceeds, the hydrogen gas generated in the catalyst layer accumulates in the cathode channel along with the electrolyte flow, causing its concentration to increase along the channel. Through... Figure 7 (e) and Figure 5 (a) Figure 6 (c) The comparison also shows that the amount of hydrogen generated at the cathode is 1 / 2 of the amount of hydrogen ions consumed at the cathode, which also satisfies the law of charge conservation and Faraday's law.

[0118] Electrolyte potential distribution and current density vector distribution in the electrolyte are as follows: Figure 8 As shown in (a) and (b), as the electrolysis reaction proceeds, oxidation occurs at the anode, causing the concentration of oxides in the anolyte to continuously increase and the concentration of reductants to continuously decrease. According to the Nernst equation, the electrode potential at the anode gradually increases. Conversely, reduction occurs at the cathode, causing the concentration of reductants in the cathode electrolyte to continuously increase and the concentration of oxides to gradually decrease. According to the Nernst equation, the electrode potential at the cathode gradually decreases. From the equivalent circuit of the electrolytic cell, it can be seen that the potential difference in the electrolyte gradually decreases as the reaction proceeds, and the current density will also gradually decrease. Therefore, from... Figure 8 It is easy to see that the current density gradually decreases along the direction of the electrolyte channel.

[0119] Temperature distribution during stable operation of the electrolytic cell is as follows Figure 9 As shown in the diagram, the heat sources in the electrolyzer are ohmic losses and exothermic reactions caused by chemical polarization. Due to the high resistance of the electrolyzer membrane, ohmic losses mainly originate from the membrane structure, while the heat source from electrode polarization is the porous catalyst layer. Therefore, the middle section of the electrolyzer is the heat source and has the highest temperature. Furthermore, the changing concentration of substances in the electrolyte channels causes the current along the channels to decrease, which in turn gradually reduces the heat flux of the heat source. Therefore, the temperature is higher near the inlet where the current is higher, and lower near the outlet where the current is lower. The isothermal surface distribution of the electrolyzer is shown in the diagram. Figure 9As shown in (b), the simulation results show that the flow field introduces convective heat transfer, resulting in faster heat transfer in the electrolyte channels and thus causing a depression in the isothermal surface. This is precisely the result of the coupling between the flow field and the temperature field. It should be noted that although the temperature distribution in the electrolyzer is non-uniform, the temperature difference across the entire electrolyzer is very small, stabilizing at around 27 degrees Celsius. This approximates an isothermal region, thus demonstrating that it is reasonable to disregard the changes in hydrodynamic parameters caused by local temperature differences in this model.

[0120] Example 2:

[0121] An industrial-grade multi-physics coupling analysis system for hydrogen production by water electrolysis includes a key physics field determination module 1, a control equation establishment module 2, a multi-physics coupling model construction module 3, and a coupling model solution module 4.

[0122] The key physical field determination module 1 is used to determine the key physical fields that affect the energy conversion efficiency of electricity-hydrogen in an industrial-grade water electrolysis hydrogen production device.

[0123] The control equation establishment module 2 is used to establish the control equations for each key physical field based on the mechanical structure and electrochemical reaction mechanism of the industrial-grade water electrolysis hydrogen production device.

[0124] The multiphysics coupling model construction module 3 is used to construct a multiphysics coupling model by combining the control equations of each key physical field and the physical process.

[0125] The coupling model solving module 4 is used to solve the multiphysics coupling model and obtain the field distribution of each key physical field.

Claims

1. A multiphysics coupling analysis method for industrial-grade water electrolysis hydrogen production, characterized in that: The analytical method includes the following steps in sequence: Step A: Identify the key physical fields affecting the electro-hydrogen energy conversion efficiency in an industrial-grade water electrolysis hydrogen production device. The key physical fields include the current density field, flow field, temperature field, and substance concentration field. Step B: Based on the mechanical structure and electrochemical reaction mechanism of the industrial-grade water electrolysis hydrogen production device, establish the governing equations for each key physical field. The governing equation for the current density field is: ； In the above formula, This is an electrode polarization overvoltage. For gases, the universal constant is... For temperature, It is Faraday's constant. It is a hyperbolic sine function. , These are current density and exchange current density, respectively. The governing equations for the flow field are: ； In the above formula, Let be the fluid density, and t be the time. For gradient operators, It is a velocity vector. For pressure, For fluid dynamic viscosity, It is the identity matrix. It is the gravitational acceleration vector. This is the surface tension term; The governing equation for the temperature field is: ； In the above formula, For constant pressure heat capacity, Thermal conductivity, It serves as the heat source for the temperature field; The governing equation for the concentration field of the substance is: ； In the above formula, Let be the flux of substance j. Let be the diffusion coefficient of substance j. , Let J represent the concentration and charge number of substance j, respectively. For electromobility, The potential of the electrolyte. For diffusion and electromigration rates, Electrolyte flow rate; Step C: Construct a multiphysics coupling model by combining the governing equations of each key physical field and the physical processes. The multiphysics coupling model includes: ； ； ； ； ； ； ； In the above formula, The rate at which the substance is deposited from the electrode surface. The current density at the electrode surface is... It is Faraday's constant. The number of charges on a substance. , These are the electrode equilibrium potential and the standard electrode equilibrium potential, respectively. , These represent the concentrations of oxides and products, respectively. For standard concentration, For constant pressure heat capacity, Thermal conductivity, As the heat source of the temperature field, For temperature, Let be the flux of substance j. Let be the diffusion coefficient of substance j. , Let J represent the concentration and charge number of substance j, respectively. For electromobility, The potential of the electrolyte. For diffusion and electromigration rates, Electrolyte flow rate , These are Joule heat sources and electrochemical heat sources, respectively. The electrode active surface area is denoted by m, which represents the cathode or anode. , These are the electrode current density vector and the electrolyte current density vector, respectively. , These are the electrode potential and the electrolyte potential, respectively. The electrode potential is... The current density of the local double layer of the electrode; Step D: Solve the multiphysics coupling model to obtain the field distribution of each key physics field.

2. The multiphysics coupling analysis method for industrial-grade water electrolysis hydrogen production according to claim 1, characterized in that: Step D uses COMSOL Multiphysics software to solve the multiphysics coupling model.

3. An industrial-grade multiphysics coupling analysis system for hydrogen production via water electrolysis, characterized in that: The system is used to execute the method of claim 1, including a key physics field determination module (1), a control equation establishment module (2), a multiphysics field coupling model construction module (3), and a coupling model solving module (4). The key physical field determination module (1) is used to determine the key physical fields that affect the energy conversion efficiency of electricity-hydrogen in an industrial-grade water electrolysis hydrogen production device. The control equation establishment module (2) is used to establish the control equations of each key physical field based on the mechanical structure and electrochemical reaction mechanism of the industrial-grade water electrolysis hydrogen production device. The multiphysics coupling model construction module (3) is used to construct a multiphysics coupling model by combining the control equations of each key physical field and the physical process. The coupling model solving module (4) is used to solve the multiphysics coupling model and obtain the field distribution of each key physical field.

Images