A data processing method and related apparatus

Abstract

The application discloses a data processing method and a related device, and relates to the technical field of electrochemistry. The method comprises the following steps: constructing a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell; performing grid division on the three-dimensional physical model based on a free triangle grid and a mapping grid to obtain a grid three-dimensional physical model; and performing multiple steady-state solving and multiple transient-state solving on the grid three-dimensional physical model to obtain target data. In this way, the change characteristics of the material transmission, electrochemical reaction and temperature field in the methane-assisted solid oxide electrolysis cell can be researched, and the target data such as the current distribution data, the voltage distribution data, the gas concentration distribution data of the cathode side and the anode side, the pressure distribution data, the velocity distribution data, and the temperature distribution data of the solid and the fluid of the methane-assisted solid oxide electrolysis cell can be obtained; and the method has the advantages of simple calculation process, fast convergence speed and accurate results.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of electrochemistry, and in particular to a data processing method and related device. BACKGROUND

[0002] Methane-assisted solid oxide electrolysis cell (CA-SOEC) refers to an electrolysis cell obtained by introducing methane (CH4) on the anode side of a solid oxide electrolysis cell (SOEC). The methane-assisted solid oxide electrolysis cell can effectively reduce the electrolysis energy consumption while preparing high-grade synthesis gas, and realizes the optimal utilization of energy and resources. How to study the change characteristics of the material transport, electrochemical reaction and temperature field in the methane-assisted solid oxide electrolysis cell has become one of the technical problems to be solved in the technical field of electrochemistry. SUMMARY

[0003] Based on the above problems, the present application provides a data processing method for studying the change characteristics of the material transport, electrochemical reaction and temperature field in the methane-assisted solid oxide electrolysis cell to obtain target data, which has the advantages of simple calculation process, fast convergence speed and accurate calculation results.

[0004] The embodiments of the present application disclose the following technical solutions:

[0005] The first aspect of the present application provides a data processing method, comprising:

[0006] constructing a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell;

[0007] dividing the three-dimensional physical model into a meshed three-dimensional physical model based on a free triangular mesh and a mapping mesh;

[0008] performing multiple steady-state solving and multiple transient-state solving on the meshed three-dimensional physical model to obtain target data; the target data includes current distribution data, voltage distribution data, gas concentration distribution data on the cathode side, gas pressure distribution data on the cathode side, gas velocity distribution data on the cathode side, gas concentration distribution data on the anode side, gas pressure distribution data on the anode side, gas velocity distribution data on the anode side, and temperature distribution data of the solid and fluid of the methane-assisted solid oxide electrolysis cell.

[0009] In an optional implementation manner, the constructing a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell comprises:

[0010] determining a plurality of geometric dimensions and a plurality of physical parameters of the methane-assisted solid oxide electrolysis cell;

[0011] constructing an initial three-dimensional physical model based on the plurality of geometric dimensions;

[0012] determining a plurality of physical fields corresponding to the initial three-dimensional physical model; the physical fields are used to describe the physical environment inside or around the three-dimensional physical model;

[0013] determining a target physical field corresponding to each of the physical fields based on a target working condition and the plurality of physical parameters;

[0014] constructing the three-dimensional physical model based on the plurality of target physical fields and the initial three-dimensional physical model.

[0015] In an optional implementation, the plurality of physical fields include a chemical physical field, a secondary current distribution physical field, a physical field of concentration mass transfer in the porous medium on the cathode side, a physical field of concentration mass transfer in the porous medium on the anode side, a physical field of Brinkman equation on the cathode side, a physical field of Brinkman equation on the anode side, and a physical field of solid and fluid heat transfer.

[0016] In an optional implementation, the plurality of geometric dimensions include an anode support layer thickness, an anode catalyst layer thickness, a cathode support layer thickness, a cathode catalyst layer thickness, an electrolysis cell thickness, an anode bipolar plate thickness, a cathode bipolar plate thickness, a size of an anode gas channel, and a size of a cathode gas channel.

[0017] In an optional implementation, the target working condition is determined by a working pressure, a working temperature, an anode inlet gas, a cathode inlet gas, a cathode inlet flow rate, and an anode inlet flow rate of the methane-assisted solid oxide electrolysis cell.

[0018] In an optional implementation, the plurality of physical parameters include constant-pressure heat capacity, thermal conductivity, electrical conductivity, and density of solid materials in the methane-assisted solid oxide electrolysis cell, porosity, permeability of porous electrode materials, and diffusion coefficient, density, viscosity, thermal conductivity, and constant-pressure heat capacity of cathode gas, and diffusion coefficient, density, viscosity, thermal conductivity, and constant-pressure heat capacity of anode gas.

[0019] In an optional implementation, the plurality of times of steady-state solving and the plurality of times of transient-state solving on the meshed three-dimensional physical model to obtain target data include:

[0020] performing a first transient-state solving on the meshed three-dimensional physical model to obtain an initial current distribution;

[0021] based on the initial current distribution, performing multiple steady-state solving on the meshed three-dimensional physical model to obtain a first intermediate calculation result;

[0022] based on the first intermediate calculation result, performing second-time transient-state solving on the meshed three-dimensional physical model to obtain a second intermediate calculation result;

[0023] based on the second intermediate calculation result, performing last-time steady-state solving on the meshed three-dimensional physical model to obtain the target data.

[0024] The second aspect of the present application provides a data processing device, the device comprising:

[0025] a first model construction module, configured to construct a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell;

[0026] a second model construction module, configured to perform meshing on the three-dimensional physical model based on a free triangle mesh and a mapped mesh to obtain a meshed three-dimensional physical model;

[0027] a target settlement result acquisition module, configured to perform multiple steady-state solving and multiple transient-state solving on the meshed three-dimensional physical model to obtain target data; the target data comprising secondary current distribution data, concentration substance transfer data in a porous medium on a cathode side, concentration substance transfer data in a porous medium on an anode side, Brinkman equation data on the cathode side, Brinkman equation data on the anode side, and solid and fluid heat transfer data.

[0028] The third aspect of the present application provides a computer readable storage medium, having a computer program stored thereon, the program being executed by a processor to implement the steps of the method introduced in any implementation manner of the first aspect.

[0029] The fourth aspect of the present application provides an electronic device, comprising:

[0030] a memory having a computer program stored thereon;

[0031] a processor configured to execute the computer program in the memory to implement the steps of the method introduced in any implementation manner of the first aspect.

[0032] Compared with the prior art, the present application has the following beneficial effects:

[0033] The application discloses a data processing method, comprising: constructing a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell; performing meshing on the three-dimensional physical model based on a free triangle mesh and a mapping mesh to obtain a meshed three-dimensional physical model; and performing multiple steady-state solving and multiple transient-state solving on the meshed three-dimensional physical model to obtain target data. The target data comprises current distribution data, voltage distribution data, gas concentration distribution data on the cathode side, gas pressure distribution data on the cathode side, gas velocity distribution data on the cathode side, gas concentration distribution data on the anode side, gas pressure distribution data on the anode side, gas velocity distribution data on the anode side, and temperature distribution data of the solid and the fluid. Thus, the method can be used for researching the change characteristics of the material transmission, electrochemical reaction and temperature field in the methane-assisted solid oxide electrolysis cell, and obtaining the target data. The method has the advantages of simple calculation process, fast convergence speed and accurate calculation result. BRIEF DESCRIPTION OF DRAWINGS

[0034] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0035] Figure 1 A flow chart of a data processing method provided by the embodiments of the present application;

[0036] Figure 2 An external schematic diagram of a methane-assisted solid oxide electrolysis cell provided by the embodiments of the present application;

[0037] Figure 3 A schematic diagram of an electrochemical reaction in a methane-assisted solid oxide electrolysis cell provided by the embodiments of the present application;

[0038] Figure 4 A schematic diagram of meshing on a three-dimensional physical model provided by the embodiments of the present application;

[0039] Figure 5 A schematic diagram of a polarization curve provided by the embodiments of the present application;

[0040] Figure 6 A structural schematic diagram of a data processing device provided by the embodiments of the present application. DETAILED DESCRIPTION

[0041] With the rapid development of new energy, renewable energy such as wind power and photovoltaic has gradually increased its proportion in the process of building new power systems. However, the intermittency and volatility of renewable energy can easily cause imbalance between power generation and power consumption, and thus the access of energy storage technology is urgently needed to maintain the smooth operation of the power system. Solid oxide electrolysis cell (SOEC) has attracted widespread attention due to its high efficiency, high reaction kinetics, flexible reactant selectivity, no noble metal catalyst and low emission.

[0042] Methane-assisted solid oxide electrolysis cell (CA-SOEC) refers to an electrolysis cell obtained by introducing methane (CH4) on the anode side of a solid oxide electrolysis cell (SOEC). CA-SOEC can effectively reduce the energy consumption of electrolysis while producing high-grade syngas, achieving optimal utilization of energy and resources. CA-SOEC has great potential in reducing the operating cost of water electrolysis for hydrogen production, improving the efficiency of water electrolysis for hydrogen production, and promoting the commercialization of high-temperature electrolysis technology. It is an indispensable part of future new energy networks.

[0043] In the traditional scheme, two methods, theoretical calculation and experimental research, can be used to study the variation characteristics of material transport, electrochemical reaction and temperature field in CA-SOEC. Due to the long research period, high cost and difficulty in measuring some physical quantities, theoretical calculation has become an important means for CA-SOEC research. In the traditional scheme, the research on CA-SOEC can only construct a two-dimensional physical model, and cannot simulate the coupling system of multiple physical fields in CA-SOEC, that is, the traditional scheme cannot conduct high-precision research on CA-SOEC.

[0044] The present application discloses a data processing method, comprising: constructing a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell; performing grid division on the three-dimensional physical model based on a free triangle grid and a mapping grid to obtain a grid-based three-dimensional physical model; performing multiple steady-state solving and multiple transient-state solving on the grid-based three-dimensional physical model to obtain target data. The target data in the present application includes secondary current distribution data, concentration substance transfer data in the porous medium on the cathode side, concentration substance transfer data in the porous medium on the anode side, Brinkman equation data on the cathode side, Brinkman equation data on the anode side, and solid and fluid heat transfer data. In this way, the method in the present application can study the variation characteristics of material transport, electrochemical reaction and temperature field in the methane-assisted solid oxide electrolysis cell and obtain target data; and has the advantages of simple calculation process, fast convergence speed and accurate calculation results.

[0045] To enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present application, and not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present application.

[0046] Figure 1 This is a flowchart illustrating a data processing method provided in an embodiment of this application. (In conjunction with...) Figure 1 As shown, the data processing method disclosed in this application includes:

[0047] S101, construct a three-dimensional physical model of a methane-assisted solid oxide electrolyzer.

[0048] A solid oxide electrolyzer is a reverse-running solid oxide fuel cell that can use solid oxides as an electrolyte to electrolyze water or other compounds (such as carbon dioxide) under high temperature (typically 650–850°C) and applied voltage conditions, thereby converting electrical and thermal energy into chemical energy.

[0049] Introducing methane into the anode side of a SOEC (Sodium Oxide Electrolyte) couples it with the SOEC process to reduce the potential difference between the anode and cathode and improve the efficiency of hydrogen production from water electrolysis. This methane-assisted solid oxide electrolyzer (CA-SOEC) effectively reduces electrolysis energy consumption while producing high-grade syngas, achieving optimal utilization of energy and resources.

[0050] In one alternative approach, the process of constructing a three-dimensional physical model of a methane-assisted solid oxide electrolyzer includes the following steps:

[0051] The first step is to determine multiple geometric dimensions and multiple physical parameters of the methane-assisted solid oxide electrolyzer.

[0052] Figure 2 This is a schematic diagram of the appearance of a methane-assisted solid oxide electrolyzer provided in an embodiment of this application. (Combined with...) Figure 2 As shown, the methane-assisted solid oxide electrolyzer of this application includes: an anode bipolar plate (1), an anode flow channel (2), an anode support layer (3), a cathode support layer (4), a cathode flow channel (5), an anode catalyst layer (6), an electrolyte (7), a cathode catalyst layer (8), and a cathode bipolar plate (9). The anode catalyst layer (6) is also called the anode active layer; the cathode catalyst layer (8) is also called the cathode active layer.

[0053] The multiple geometric dimensions of the methane-assisted solid oxide electrolysis cell in the present application include: an anode support layer thickness, an anode catalytic layer thickness, a cathode support layer thickness, a cathode catalytic layer thickness, an electrolysis cell thickness, an anode bipolar plate thickness, a cathode bipolar plate thickness, an anode gas channel size, and a cathode gas channel size.

[0054] Exemplarily, Figure 2 The width of the methane-assisted solid oxide electrolysis cell shown in the figure is 2 mm, the thickness of the Ni-YSZ anode support layer (3) is 210 μm, the thickness of the Ni-YSZ anode active layer (6) is 15 μm, the thickness of the YSZ electrolyte layer (7) is 20 μm, the thickness of the Ni-YSZ cathode active layer (8) is 15 μm, and the thickness of the Ni-YSZ cathode support layer (4) is 210 μm.

[0055] The multiple physical parameters in the present application include the constant-pressure heat capacity, thermal conductivity, electrical conductivity, and density of the solid materials in the methane-assisted solid oxide electrolysis cell, the porosity, permeability of the porous electrode materials, and the diffusion coefficient, density, viscosity, thermal conductivity, and constant-pressure heat capacity of the cathode gas, and the diffusion coefficient, density, viscosity, thermal conductivity, and constant-pressure heat capacity of the anode gas.

[0056] Exemplarily, the porosity of the porous electrode material is 0.36, the tortuosity factor is 3, and the permeability is 10 -10 m 2 The cathode inlet gas is composed of 80% H2O and 20% H2, and the anode inlet gas is composed of 29% CH4 and 71% H2O.

[0057] Secondly, based on the multiple geometric dimensions of the methane-assisted solid oxide electrolysis cell, an initial three-dimensional physical model is constructed.

[0058] It can be understood that the geometric dimensions of the constructed initial three-dimensional physical model are exactly the same as the geometric dimensions of the methane-assisted solid oxide electrolysis cell.

[0059] Thirdly, multiple physical fields corresponding to the initial three-dimensional physical model are determined.

[0060] The physical field is used to describe the physical environment inside or around the three-dimensional physical model.

[0061] The multiple physical fields corresponding to the initial three-dimensional model of the methane-assisted solid oxide electrolysis cell in the present application include the chemical physical field, the secondary current distribution physical field, the physical field of concentration substance transfer in the porous medium on the cathode side, the physical field of concentration substance transfer in the porous medium on the anode side, the physical field of Brinkman equation on the cathode side, the physical field of Brinkman equation on the anode side, and the physical field of solid and fluid heat transfer.

[0062] In the fourth step, a target physical field corresponding to each of the physical fields is determined based on the target working condition and the plurality of physical parameters.

[0063] The target working condition is determined by the working pressure, the working temperature, the anode inlet gas, the cathode inlet gas, the cathode inlet flow rate, and the anode inlet flow rate of the methane-assisted solid oxide electrolysis cell.

[0064] For example, the working pressure of the methane-assisted solid oxide electrolysis cell corresponding to the target working condition is 1 atm, the working temperature is 700°C, the cathode inlet gas is composed of 80% H2O and 20% H2, the anode inlet gas is composed of 29% CH4 and 71% H2O, the cathode inlet flow rate is 0.5 m / s, and the anode inlet flow rate is 0.16 m / s.

[0065] In the determination of the target working condition, the physical parameters involved in each physical field are defined, and the target physical parameters corresponding to the physical field can be obtained.

[0066] For example, under the target working condition, the generation method of the target physical field corresponding to the chemical physical field (referred to as a chemical module) is as follows: after adding the thermodynamic system on the cathode and anode sides, selecting the gas system, generating the chemical system, and selecting the mass transfer, the concentration substance transfer is selected. The cathode side substances are hydrogen and water vapor, and the anode side substances are methane, water vapor, carbon monoxide, hydrogen, and carbon dioxide.

[0067] Figure 3 A schematic diagram of an electrochemical reaction in a methane-assisted solid oxide electrolysis cell is provided for the embodiments of the present application. In combination with Figure 3 As shown in the figure, the water vapor at the cathode of the methane-assisted solid oxide electrolysis cell obtains electrons to generate hydrogen and oxygen ions, and the oxygen ions are transmitted to the anode through the electrolyte; the methane and water vapor at the anode react to generate hydrogen and carbon monoxide, which further undergoes an electrochemical reaction with the oxygen ions conducted from the cathode to generate carbon dioxide and water.

[0068] In the anode chemical system, the methane reforming and water vapor shift reaction rates are calculated through reversible reactions. Specifically, as shown in formulas (1) and (2):

[0069]

[0070] In formulas (1) and (2), the letters have the following meanings: R MSR represents the methane water vapor reforming reaction rate, R WGSR represents the reversible water vapor shift reaction rate; k rf represents the reaction rate constant of the reforming reaction rate; k sf represents the reaction rate constant of the reversible water vapor shift reaction rate; K pr represents the reaction equilibrium constant of the reforming reaction rate; Kps K represents a reaction equilibrium constant representing a rate of a reversible water vapor shift reaction; p CO and represent partial pressures of corresponding gases, respectively.

[0071] Exemplarily, under a target working condition, a generation method of a target physical field corresponding to a secondary current distribution physical field (referred to as a secondary current distribution module) is as follows: a charge conservation equation is selected to calculate a current density to construct the target physical field corresponding to the secondary current distribution module.

[0072] An expression of the charge conservation equation is shown in formula (3):

[0073]

[0074] In formula (3), meanings of each letter are as follows: i l represents an ion current density, i s represents an electron current density, represents a gradient operator, σ l represents an ion conductivity, σ s represents an electron conductivity, φ ss represents an ion potential, φ ll represents an electron potential, Q l represents a charge source term of an ion current, Q s represents a charge source term of an electron current.

[0075] The electrode is composed of an electron conductor and an ion conductor, and thus, intrinsic electron conductivity and intrinsic ion conductivity thereof can be corrected by formula (4) and formula (5). Expressions of formula (4) and formula (5) are as follows:

[0076] σ l,eff = θ(1- ∈)τ -1 σ l (4)

[0077] σ s,eff = θ(1- ∈)τ -1 σ s (5)

[0078] In formula (4) and formula (5), θ represents a volume fraction of a conductor, ∈ represents a porosity of the electrode, τ represents a tortuosity factor of the electrode, σ l represents an intrinsic ion conductivity of a material, σ s represents an intrinsic electron conductivity of a material, σ l,eff represents an ion effective conductivity, and σ s,eff represents an electron effective conductivity.

[0079] Wherein, the relationship between the activation polarization of the anode and the cathode and the current density is described by the Butler-Volmer (BV) equation. The expression of BV is shown in formula (6) and formula (7):

[0080]

[0081] In formula (6) and formula (7), j an represents the anode electrode reaction source (A / m 3 ), j ca represents the cathode electrode reaction source (A / m 3 ), Av represents the active specific surface area of the electrode (2.14e 5 m -1 ), R represents the gas constant; T represents the working temperature; F represents the Faraday constant, n is the number of electrons involved in the reaction, that is, the number of transferred charges (generally taken as 2), α an represents the anode charge transfer coefficient (generally taken as 0.5), and α ca represents the cathode charge transfer coefficient (generally taken as 0.5).

[0082] In formula (6) and formula (7), and j ca represent the exchange current densities of the anode and the cathode; and j ca have the same unit, A / m 2 .

[0083] Since the anode simultaneously occurs the electrochemical reactions of hydrogen and carbon monoxide, it is necessary to calculate the exchange current densities of the two electrochemical reactions respectively. Among them, the electrolysis rate of water is 2.2 times of the electrolysis rate of carbon dioxide, so there are formula (8), formula (9) and formula (10).

[0084]

[0085] In formula (6) and formula (7), η an is the anode active overpotential, with the unit of V; η ca represents the cathode active overpotential, with the unit of V; γ an represents the anode pre-exponential factor (in the present application, the value of γ an is 7×10 8 A / m 2 ), and γ ca represents the cathode pre-exponential factor (in the present application, the value of γ ca is 1.5×10 9 A / m 2 ), E act,an is the anode activation energy (in the present application, the value of E act,an is 1.3×105 J / mol), E act,ca is the cathode activation energy (E act,ca The value of E 5 J / mol), is the water vapor partial pressure at the three interfaces of the cathode (unit: Pa), is the water vapor partial pressure at the three interfaces of the anode (unit: Pa), is the equilibrium oxygen partial pressure of the anode (unit: Pa), is the equilibrium oxygen partial pressure of the cathode (unit: Pa).

[0086] wherein η an and η ca can be calculated by formula (11), formula (12) and formula (13):

[0087] η an = φ s - φ l - E op (11)

[0088] η ca = φ s - φ l (12)

[0089]

[0090] E op in formula (11), formula (12) and formula (13) represents the start-up voltage of the SOEC electrolysis reaction, which is theoretically only related to temperature, anode and cathode gas composition; and respectively represent the equilibrium oxygen partial pressure on the anode and cathode sides, which can be calculated by using Hydrocarbon System Calculator (HSC).

[0091] Exemplarily, under the target working condition, the generation method of the target physical field corresponding to the physical field of the concentration substance transfer in the porous medium (referred to as the concentration substance transfer module in the porous medium) is: a general diffusion convection equation is selected to obtain the target physical field corresponding to the concentration substance transfer module in the porous medium.

[0092] The expression of the general diffusion convection equation is shown in formula (14):

[0093]

[0094] In formula (14), represents the gradient operator, j idenotes the diffusion flux of gas i, p denotes the density of the gas, u denotes the gas velocity, and i denotes the mass fraction of gas i, R i denotes the reaction source term of gas i.

[0095] For a porous electrode, gas diffusion occurs through free molecule diffusion and Knudsen diffusion, free molecule diffusion dominates in large pores, and Knudsen diffusion becomes significant when the pore size is smaller than the average free path of the molecules, the present model uses the Knudsen-extended Fick model to describe gas transport in a porous electrode, and the gas composition at different positions can be calculated.

[0096] Exemplarily, under the target working condition, the generation method of the target physical field corresponding to the physical field of solid and fluid heat transfer (referred to as a solid and fluid heat transfer module) is: obtaining the target physical field corresponding to the solid and fluid heat transfer module through a general thermodynamic conservation equation.

[0097] The general diffusion convection equation is expressed as formula (15):

[0098]

[0099] In formula (15), C p denotes the constant pressure heat capacity of the material, k denotes the thermal conductivity of the material, Q denotes the heat source term, and the heat source term is composed of an ohmic heat source (S ohm ), an irreversible heat source (S rev ), and a reversible heat source (S irr ).

[0100] Among them, the bipolar plate, the electrolyte layer and the cathode support layer are set as ohmic heat sources, the anode support layer is set as an ohmic heat source and a reversible heat source, and the anode catalyst layer and the cathode catalyst layer are set as ohmic heat sources, irreversible heat sources and reversible heat sources.

[0101] In the fifth step, the three-dimensional physical model is constructed based on the plurality of target physical fields and the initial three-dimensional physical model.

[0102] That is, a plurality of target physical fields are added to the initial three-dimensional physical model to obtain the three-dimensional physical model in the present application.

[0103] In S102, the three-dimensional physical model is meshed based on the free triangular mesh and the mapping mesh to obtain a meshed three-dimensional physical model.

[0104] Figure 4 A schematic diagram of meshing a three-dimensional physical model is provided for the embodiments of the present application. Figure 4As shown, the bipolar plate, gas flow channel and other areas are divided into unstructured grids, and the electrode support layer, electrode catalyst layer, electrolyte and other areas are divided into structured grids. This combination of structured grids and unstructured grids can better adapt to the narrow and long geometric areas of the electrolyte and electrode, reduce the number of grids, and improve the calculation efficiency.

[0105] The specific grid division steps are as follows:

[0106] (1) Select the cathode inlet and anode inlet surface, and divide the cathode, anode flow channel and bipolar plate into free triangular grids;

[0107] (2) Divide the anode catalyst layer, cathode catalyst layer, anode support layer, cathode support layer and electrolyte layer into mapping grids;

[0108] (3) Since the fluid has a large gradient near the wall surface, the free triangular grids at the junction of the flow channel and the bipolar plate are divided into boundary layer grids;

[0109] (4) Take the divided grid surface as the source surface, and scan the remaining part. Since the gas component has a large gradient at the model inlet and outlet, the scanning distribution adopts symmetric distribution, and the cell size ratio is set to 20; then the grid division part is completed.

[0110] S103, multiple steady-state solving and multiple transient-state solving are performed on the gridded three-dimensional physical model to obtain target data.

[0111] In an optional implementation, the multiple steady-state solving and multiple transient-state solving on the gridded three-dimensional physical model to obtain target data includes the following four steps, including:

[0112] First, first transient-state solving is performed on the gridded three-dimensional physical model to obtain an initial current distribution; this step is referred to as transient-state solving initial current distribution.

[0113] Second, based on the initial current distribution, multiple steady-state solving is performed on the gridded three-dimensional physical model to obtain a first intermediate calculation result.

[0114] Specifically, in this step, 6 steady-state calculations are performed. Specifically, through the first steady-state calculation, the steady-state current distribution, the Brinkman equation on the anode side is obtained; through the second steady-state calculation, the Brinkman equation on the cathode side is solved for the secondary current distribution; through the third steady-state calculation, the concentration mass transfer in the porous medium on the anode side is solved for the secondary current distribution; through the fourth steady-state calculation, the concentration mass transfer in the porous medium on the cathode side is solved for the secondary current distribution; through the fifth steady-state calculation, the solid and fluid heat transfer is solved; through the sixth steady-state calculation, the concentration mass transfer in the porous medium on the cathode and anode sides, the Brinkman equation on the cathode and anode sides for the secondary current distribution is solved, and the calculation results obtained in the second step are collectively referred to as the first intermediate calculation results.

[0115] In the third step, based on the first intermediate calculation results, a second transient solution is performed on the grid-based three-dimensional physical model to obtain second intermediate calculation results.

[0116] The second intermediate calculation results include the secondary current distribution, the concentration mass transfer in the porous medium on the cathode and anode sides, the Brinkman equation on the cathode and anode sides, and the solid and fluid heat transfer.

[0117] In the fourth step, based on the second intermediate calculation results, a final steady-state solution is performed on the grid-based three-dimensional physical model to obtain the target data.

[0118] The target data in the present application includes methane-assisted solid oxide cell current distribution data, voltage distribution data, cathode-side gas concentration distribution data, cathode-side gas pressure distribution data, cathode-side gas velocity distribution data, anode-side gas concentration distribution data, anode-side gas pressure distribution data, anode-side gas velocity distribution data, and solid and fluid temperature distribution data. The cathode-side gas includes but is not limited to methane, water vapor, carbon monoxide, carbon dioxide, etc. The anode-side gas includes but is not limited to water vapor and hydrogen.

[0119] It should be noted that in the process of obtaining the target data in the present application, the third step, i.e., based on the first intermediate calculation results, a second transient solution is performed on the grid-based three-dimensional physical model to obtain second intermediate calculation results, which has the advantage of making the calculation results converge and the calculation results more accurate.

[0120] Figure 5 A polarization curve diagram is provided for the embodiments of the present application. As shown in Figure 5 The horizontal axis represents voltage (in V), and the vertical axis represents current (in (A / cm 2 ); Figure 5The small squares in the diagram represent the polarization curves obtained using the data processing method described in this application. Figure 5 The small triangles in the diagram represent the polarization curves obtained experimentally under the same target operating conditions.

[0121] Among them, we obtained Figure 5 The process of representing the polarization curve with small and medium squares involves adding an auxiliary scan in the final step of acquiring the target data. Using the operating voltage as a parameter, the command `range(0, 0.1, 0.8)` is entered into the parameter value list. This will generate a graph showing the relationship between current and operating voltage after calculation. Figure 5 The polarization curve in.

[0122] Combination Figure 5 As shown, the polarization curve obtained by the method in this application is basically the same as the target curve obtained experimentally.

[0123] It should be noted that, Figure 5 The two polarization curves obtained correspond to the same target operating conditions: the operating temperature is 700℃, the cathode inlet gas consists of 80% H2O and 20% H2, the anode inlet gas consists of 29% CH4 and 71% H2O, the cathode inlet flow rate is 0.5 m / s, and the anode inlet flow rate is 0.16 m / s.

[0124] As can be seen from the foregoing embodiments, the data processing method disclosed in this application has the following advantages:

[0125] First, a complete three-dimensional methane-assisted solid oxide electrolysis cell model can be established. Compared with the existing two-dimensional methane-assisted solid oxide electrolysis cell model, the method in this application considers the influence of the ribs on the reversible water-gas shift reaction and the reversible methanation reaction rate, which is closer to actual production.

[0126] Second, since the methane-assisted solid oxide electrolyzer model has many source terms and is complex, and the methane reforming reaction is endothermic, causing a sudden drop in temperature at the inlet and making convergence difficult, this application adopts a step-by-step coupled calculation method that combines transient and steady-state calculations. The entire study is divided into four major calculation steps (two transient calculations and seven steady-state calculations). The results of the transient calculations are used as the initial values ​​for the steady-state calculations, which solves the problem of non-convergence in the multi-physics coupling process caused by the complexity of the model.

[0127] Third, this application combines structured and unstructured meshes, which reduces the number of meshes required while ensuring mesh quality. The model has a fast calculation speed, converges in 30 iterations under fully coupled conditions, and the calculation results are accurate, with an error of less than 5% compared with experimental results under the same working conditions.

[0128] Fourthly, based on the actual operation parameters and results of the methane assisted solid oxide electrolysis cell, a complete three-dimensional model of the methane assisted solid oxide electrolysis cell is established, which considers the processes of mass transfer, heat transfer, motion transfer, chemical reaction and electrochemical reaction, can simulate the temperature distribution, gas concentration distribution, current density distribution and other parameters under different working conditions, and has important guiding significance for the experimental study of the methane assisted solid oxide electrolysis cell.

[0129] Based on the data processing method disclosed in the foregoing embodiments, the application further discloses a data processing device.

[0130] Figure 6 A structural schematic diagram of a data processing device provided by the embodiments of the application is shown. As shown in the figure, the data processing device 600 disclosed by the application comprises: Figure 6

[0131] A first model construction module 601 is configured to construct a three-dimensional physical model of the methane assisted solid oxide electrolysis cell.

[0132] A second model construction module 602 is configured to perform grid division on the three-dimensional physical model based on a free triangular grid and a mapping grid to obtain a grid-based three-dimensional physical model.

[0133] A target settlement result acquisition module 603 is configured to perform multiple steady-state solving and multiple transient-state solving on the grid-based three-dimensional physical model to obtain target data; the target data comprises secondary current distribution data, concentration substance transfer data in the porous medium on the cathode side, concentration substance transfer data in the porous medium on the anode side, Brinkman equation data on the cathode side, Brinkman equation data on the anode side, and solid and fluid heat transfer data.

[0134] In an optional implementation, the first model construction module 601 comprises:

[0135] A parameter acquisition unit is configured to determine a plurality of geometric dimensions and a plurality of physical parameters of the methane assisted solid oxide electrolysis cell.

[0136] An initial model construction unit is configured to construct an initial three-dimensional physical model based on the plurality of geometric dimensions.

[0137] A physical field determination unit is configured to determine a plurality of physical fields corresponding to the initial three-dimensional physical model; the physical fields are used to describe the physical environment inside or around the three-dimensional physical model.

[0138] A target physical field acquisition unit is configured to determine a target physical field corresponding to each of the physical fields based on a target working condition and the plurality of physical parameters.

[0139] ​The three-dimensional physical model acquisition unit is configured to construct the three-dimensional physical model based on the plurality of target physical fields and the initial three-dimensional physical model.

[0140] In an alternative implementation, the target settlement result acquisition module 603 includes:

[0141] The first data acquisition unit is configured to perform first transient-state solving on the meshed three-dimensional physical model to obtain an initial current distribution.

[0142] The second data acquisition unit is configured to perform multiple steady-state solving on the meshed three-dimensional physical model based on the initial current distribution to obtain a first intermediate calculation result.

[0143] The third data acquisition unit is configured to perform second transient-state solving on the meshed three-dimensional physical model based on the first intermediate calculation result to obtain a second intermediate calculation result.

[0144] The fourth data acquisition unit is configured to perform last steady-state solving on the meshed three-dimensional physical model based on the second intermediate calculation result to obtain the target data.

[0145] Based on the data processing method and device provided in the foregoing embodiments, correspondingly, the present application further provides a computer readable storage medium having a computer program stored thereon, and the program is executed by a processor to implement some or all steps of the data processing method mentioned above.

[0146] Based on the data processing method and device provided in the foregoing embodiments, the present application further provides an electronic device, which includes:

[0147] A memory having a computer program stored thereon;

[0148] A processor configured to execute the computer program in the memory to implement some or all steps of the data processing method provided in the foregoing embodiments.

[0149] It should be noted that each of the embodiments of the present specification is described in a progressive manner, and the same or similar parts between the embodiments can be mutually referred to, and each of the embodiments focuses on the differences from other embodiments. In particular, for the device embodiments, since they are basically similar to the method embodiments, they are described more simply, and the relevant parts can be referred to the part of the description of the method embodiments. The above-described device embodiments are only illustrative, and the units described as separate components can or can not be physically separated, and the components indicated as units can or can not be physical units, that is, they can be located in one place, or distributed on multiple network units. Part or all of the modules can be selected to achieve the purpose of the embodiments according to the actual needs. Those skilled in the art can understand and implement it without creative labor.

[0150] The above describes only one specific embodiment of the present application, but the protection scope of the present application is not limited to this. Any skilled person in the art can easily think of changes or replacements within the technical range disclosed in the present application, which should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A data processing method, characterized by, The method comprises: constructing a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell; meshing the three-dimensional physical model based on a free triangle mesh and a mapped mesh to obtain a meshed three-dimensional physical model; performing first transient solving on the meshed three-dimensional physical model to obtain an initial current distribution; performing six times of steady-state solving on the meshed three-dimensional physical model based on the initial current distribution to obtain a first intermediate calculation result; the first intermediate calculation result comprises a secondary current distribution, a Brinkman equation on the anode side, a Brinkman equation on the cathode side, concentration substance transfer in a porous medium on the anode side, concentration substance transfer in a porous medium on the cathode side, and solid and fluid heat transfer; performing second transient solving on the meshed three-dimensional physical model based on the first intermediate calculation result to obtain a second intermediate calculation result; the second intermediate calculation result comprises a secondary current distribution, concentration substance transfer in a porous medium on the anode side, concentration substance transfer in a porous medium on the cathode side, a Brinkman equation on the cathode side, a Brinkman equation on the anode side, and solid and fluid heat transfer; performing last time of steady-state solving on the meshed three-dimensional physical model based on the second intermediate calculation result to obtain target data; the target data comprises current distribution data, voltage distribution data, gas concentration distribution data on the cathode side, gas pressure distribution data on the cathode side, gas velocity distribution data on the cathode side, gas concentration distribution data on the anode side, gas pressure distribution data on the anode side, gas velocity distribution data on the anode side, and temperature distribution data of solid and fluid of the methane-assisted solid oxide electrolysis cell.

2. The method of claim 1, wherein, The method comprises: constructing a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell; determining a plurality of geometric dimensions and a plurality of physical parameters of the methane-assisted solid oxide electrolysis cell; constructing an initial three-dimensional physical model based on the plurality of geometric dimensions; determining a plurality of physical fields corresponding to the initial three-dimensional physical model; the physical fields are used to describe the physical environment inside or around the three-dimensional physical model; determining a target physical field corresponding to each of the physical fields based on a target working condition and the plurality of physical parameters; 3. The method of claim 2, wherein, constructing the three-dimensional physical model based on the plurality of target physical fields and the initial three-dimensional physical model.

4. The method of claim 2, wherein, The plurality of physical fields comprises a chemical physical field, a secondary current distribution physical field, a concentration substance transfer in a porous medium on the cathode side physical field, a concentration substance transfer in a porous medium on the anode side physical field, a Brinkman equation on the cathode side physical field, a Brinkman equation on the anode side physical field, and a solid and fluid heat transfer physical field. The plurality of geometric dimensions comprises an anode support layer thickness, an anode catalyst layer thickness, a cathode support layer thickness, a cathode catalyst layer thickness, an electrolysis cell thickness, an anode bipolar plate thickness, a cathode bipolar plate thickness, an anode gas channel size, and a cathode gas channel size.

5. The method of claim 2, wherein, The target working condition is determined by the working pressure, working temperature, anode inlet gas, cathode inlet gas, cathode inlet flow rate and anode inlet flow rate of the methane-assisted solid oxide electrolysis cell.

6. The method of claim 2, wherein, The plurality of physical parameters include the constant-pressure heat capacity, thermal conductivity, electrical conductivity and density of the solid material in the methane-assisted solid oxide electrolysis cell, the porosity and permeability of the porous electrode material, and the diffusion coefficient, density, viscosity, thermal conductivity and constant-pressure heat capacity of the cathode gas and the diffusion coefficient, density, viscosity, thermal conductivity and constant-pressure heat capacity of the anode gas.

7. A data processing apparatus, characterized by The device comprises: A first model construction module for constructing a three-dimensional physical model of a methane-assisted solid oxide electrolysis cell; A second model construction module for meshing the three-dimensional physical model based on a free triangle mesh and a mapped mesh to obtain a meshed three-dimensional physical model; A target settlement result acquisition module for performing first transient solving on the meshed three-dimensional physical model to obtain an initial current distribution, performing six times of steady-state solving on the meshed three-dimensional physical model based on the initial current distribution to obtain a first intermediate calculation result, the first intermediate calculation result including a secondary current distribution, an anode-side Brinkman equation, a cathode-side Brinkman equation, a concentration substance transfer in a porous medium on the anode side, a concentration substance transfer in a porous medium on the cathode side and solid and fluid heat transfer, performing second transient solving on the meshed three-dimensional physical model based on the first intermediate calculation result to obtain a second intermediate calculation result, the second intermediate calculation result including a secondary current distribution, a concentration substance transfer in a porous medium on the anode side, a concentration substance transfer in a porous medium on the cathode side, a cathode-side Brinkman equation, an anode-side Brinkman equation and solid and fluid heat transfer, and performing last steady-state solving on the meshed three-dimensional physical model based on the second intermediate calculation result to obtain target data, the target data including current distribution data, voltage distribution data, gas concentration distribution data on the cathode side, gas pressure distribution data on the cathode side, gas velocity distribution data on the cathode side, gas concentration distribution data on the anode side, gas pressure distribution data on the anode side, gas velocity distribution data on the anode side and temperature distribution data of the solid and fluid of the methane-assisted solid oxide electrolysis cell.

8. A computer-readable storage medium having stored thereon a computer program, characterized in that, The program is executed by a processor to implement the steps of the method of any one of claims 1-6.

9. An electronic device, comprising: Comprise: A memory having a computer program stored thereon; A processor configured to execute the computer program in the memory to implement the steps of the method of any one of claims 1-6.

Images