Methane catalytic cracking simulation method, system, terminal and medium

Abstract

The application provides a methane catalytic cracking simulation method, a system, a terminal and a medium, wherein the methane catalytic cracking simulation method is a simulation process executed in cycles, and each simulation process comprises the following steps: based on a reference flow field condition, a reference reaction condition is combined, a flow field calculation is performed, and an updated flow field condition is obtained; based on the updated flow field condition and the reference reaction condition, a material transport reaction rule is combined, and an updated reaction condition is obtained; based on the reference flow field condition and the updated flow field condition, and the reference reaction condition and the updated reaction condition, convergence verification is performed, if the convergence verification is passed, the updated flow field condition and the updated reaction condition are output, and the iterative simulation is stopped; otherwise, the simulation process is re-executed. The application can accurately predict the nucleation, growth, migration and accumulation behavior of carbon particles, so that better simulation effect is achieved.

Description

Technical Field

[0001] This application belongs to the field of computational fluid dynamics simulation technology, and relates to a simulation technology for simulating the cracking of methane in a molten medium to produce hydrogen, and particularly relates to a simulation method, system, terminal and medium for catalytic cracking of methane. Background Technology

[0002] Methane catalytic cracking technology is one of the important methods for industrial hydrogen production. Simulation analysis of methane catalytic cracking is helpful in predicting the temperature distribution and product composition within the reactor, thereby assisting in the design of efficient and low-energy-consumption methane catalytic cracking hydrogen production schemes.

[0003] Currently, CFD (Computational Fluid Dynamics) models are commonly used for simulation. However, carbon particles are generated during the catalytic cracking of methane, and existing methods often simplify carbon particles as gaseous components or ignore them altogether. This makes it impossible to accurately predict the nucleation, growth, migration, and accumulation behavior of carbon particles, resulting in low accuracy in the simulation of methane catalytic cracking. Summary of the Invention

[0004] The purpose of this application is to provide a method, system, terminal and medium for simulating methane catalytic cracking, in order to solve the problem that the accuracy of the simulation of carbon particle behavior in the prior art is not high, which in turn affects the accuracy of the simulation of methane catalytic cracking.

[0005] In a first aspect, this application provides a methane catalytic cracking simulation method for simulating methane catalytic cracking reactions in target areas and iteratively simulating the methane catalytic cracking reaction process in each target area;

[0006] The iterative simulation is a cyclically executed simulation process, with a preset time step between each simulation process; each simulation process includes:

[0007] Based on the reference flow field and the reference reaction, flow field calculations are performed to obtain updated flow field conditions. The reference flow field conditions include reference pressure distribution data, reference velocity distribution data, reference temperature distribution data, reference volume fraction data, and reference gas-liquid interface data. The updated flow field conditions include updated pressure distribution data, updated velocity distribution data, updated temperature distribution data, updated volume fraction data, and updated gas-liquid interface data.

[0008] Based on the updated flow field and the reference reaction conditions, and combined with the laws governing mass transport and reaction, updated reaction conditions are obtained; wherein, the reference reaction conditions include reference component concentration distribution data, reference reaction rate data, reference carbon mass data, and reference carbon motion data; the updated reaction conditions include updated component concentration distribution data, updated reaction rate data, updated carbon mass data, and updated carbon motion data;

[0009] Based on the reference flow field and the updated flow field, as well as the reference response and the updated response, convergence verification is performed. If the convergence verification passes, the updated flow field and the updated response are output, and the iterative simulation is stopped. Otherwise, the updated flow field is used as the new reference flow field and the updated response is used as the new reference response, and the simulation process is re-executed.

[0010] In one embodiment of this application, the step of performing flow field calculations based on a reference flow field and a reference reaction to obtain updated flow field conditions includes:

[0011] Extract the reference velocity distribution data and the reference pressure distribution data, and calculate the intermediate velocity distribution data using the law of conservation of momentum; perform pressure correction based on the intermediate velocity distribution data to obtain the updated pressure distribution data; and perform velocity correction based on the updated pressure distribution data to obtain the updated velocity distribution data.

[0012] Extract the reference reaction rate data and calculate the updated temperature distribution data using the law of energy conservation;

[0013] Based on the updated velocity distribution data, the updated volume fraction data is calculated using the gas-liquid transport law; based on the updated volume fraction data, the interface is reconstructed to obtain the updated gas-liquid interface;

[0014] The updated pressure distribution data, the updated velocity distribution data, the updated temperature distribution data, the updated volume fraction data, and the updated gas-liquid interface data are output as the updated flow field conditions.

[0015] In one embodiment of this application, obtaining the updated reaction information based on the updated flow field and the reference reaction information, combined with the laws governing mass transport and reaction, includes:

[0016] Extract the reference component concentration distribution data and the updated temperature distribution data, and calculate the updated reaction rate data using the laws governing methane catalytic cracking reaction;

[0017] Based on the updated reaction rate data, combined with the reference component concentration distribution data and the updated rate distribution data, the updated component concentration distribution data is calculated using the component transport laws.

[0018] Based on the updated reaction rate data, the reference carbon mass data is updated, and the updated carbon mass data is calculated and obtained.

[0019] Based on the updated carbon mass data and the updated velocity distribution data, the stress on carbon particles is analyzed to update the reference carbon motion data and obtain the updated carbon motion data.

[0020] In one embodiment of this application, updating the reference carbon mass data based on the updated reaction rate data and calculating the updated carbon mass data includes:

[0021] Obtain the reference carbon particle diameter data, and combine it with the updated reaction rate data to obtain the updated carbon particle diameter data;

[0022] Based on the updated carbon particle diameter data, the updated carbon mass data is calculated and obtained.

[0023] In one embodiment of this application, the reference carbon motion data includes reference carbon position data and reference carbon velocity data; the updated carbon motion data includes updated carbon position data and updated carbon velocity data; the step of analyzing the force on carbon particles based on the updated carbon mass data and the updated velocity distribution data to update the reference carbon motion data and obtain the updated carbon motion data includes:

[0024] Based on the updated carbon mass data and the updated carbon velocity data, acceleration data is obtained;

[0025] The reference carbon velocity data is updated based on the acceleration data to obtain the updated carbon velocity data;

[0026] The reference carbon position data is updated based on the updated carbon velocity data to obtain the updated carbon position data;

[0027] Based on the updated carbon velocity data and the updated carbon position data, the updated carbon motion data is obtained.

[0028] In one embodiment of this application, the convergence verification includes:

[0029] Based on the reference volume fraction data and the updated volume fraction data, obtain volume fraction residual data; based on the reference component concentration distribution data and the updated component concentration distribution data, obtain component concentration residual data; based on the reference carbon motion data and the updated carbon motion data, obtain carbon motion residual data.

[0030] If the volume fraction residual data does not exceed the preset upper limit of the volume fraction residual data, the component concentration residual data does not exceed the preset upper limit of the component concentration residual data, and the carbon motion residual data does not exceed the preset upper limit of the carbon motion residual data, then the convergence verification is passed; otherwise, the convergence verification is not passed.

[0031] In one embodiment of this application, during the first execution of the simulation process:

[0032] Linear distribution calculation is performed based on a preset top pressure value to obtain the reference pressure distribution data;

[0033] And / or, based on a preset speed value, obtain the reference speed distribution data;

[0034] And / or, based on a preset reaction rate value, obtain the reference reaction rate data;

[0035] And / or, based on preset bubble size, preset bubble shape and preset bubble position, the interface is reconstructed to obtain the reference volume fraction data and the reference gas-liquid interface;

[0036] And / or, based on the reference gas-liquid interface, obtain the gas region and the liquid region; use a preset gas temperature value as the temperature value of the gas region and a preset liquid temperature value as the temperature value of the liquid region to obtain the reference temperature distribution data;

[0037] And / or, based on the preset component concentration values ​​corresponding to each component, obtain the reference component concentration distribution data;

[0038] And / or, based on a preset mass value, obtain the reference carbon mass data;

[0039] And / or, based on a preset velocity value, combined with the gas region, obtain the reference carbon motion data.

[0040] Secondly, this application provides a methane catalytic cracking simulation system for iteratively simulating methane catalytic cracking reactions in a target region, including a flow field acquisition module, a reaction acquisition module, and a convergence judgment module;

[0041] The flow field acquisition module is used to perform flow field calculations based on a reference flow field and a reference reaction to obtain updated flow field conditions. The reference flow field conditions include reference pressure distribution data, reference velocity distribution data, reference temperature distribution data, reference volume fraction data, and a reference gas-liquid interface. The updated flow field conditions include updated pressure distribution data, updated velocity distribution data, updated temperature distribution data, updated volume fraction data, and updated gas-liquid interface.

[0042] The reaction status acquisition module is used to acquire updated reaction status based on the updated flow field status and the reference reaction status, combined with the laws governing mass transport and reaction. The reference reaction status includes reference component concentration distribution data, reference reaction rate data, reference carbon mass data, and reference carbon motion data. The updated reaction status includes updated component concentration distribution data, updated reaction rate data, updated carbon mass data, and updated carbon motion data.

[0043] The convergence judgment module is used to perform convergence verification based on the reference flow field situation and the updated flow field situation, as well as the reference reaction situation and the updated reaction situation. If the convergence verification is passed, the updated flow field situation and the updated reaction situation are output, and the iterative simulation is stopped; otherwise, the updated flow field situation is used as the new reference flow field situation, and the updated reaction situation is used as the new reference reaction situation, and the process returns to the flow field situation acquisition module and the reaction situation acquisition module.

[0044] Thirdly, this application provides a terminal, including: a processor and a memory, wherein the memory and the processor are communicatively connected;

[0045] The memory is used to store computer programs, and the processor is used to execute the computer programs stored in the memory to enable the terminal to perform the methane catalytic cracking simulation method as described above.

[0046] Fourthly, this application provides a computer storage medium storing a computer program that, when executed by a processor, implements the methane catalytic cracking simulation method as described above.

[0047] As described above, this application provides a method, system, terminal, and medium for simulating methane catalytic cracking. By simulating the generation, mass change, and movement of carbon particles, it accurately predicts the nucleation, growth, migration, and accumulation behavior of carbon particles, which is beneficial for more accurate prediction of temperature distribution and product composition within the reactor, thereby achieving better simulation results. Attached Figure Description

[0048] Figure 1The diagram shown is a schematic flow chart of a single simulation process in a methane catalytic cracking simulation as described in an embodiment of this application.

[0049] Figure 2 The diagram shown is a flowchart illustrating a flow field calculation method according to an embodiment of this application.

[0050] Figure 3 The diagram shown is a flowchart illustrating a method for obtaining update response status as described in an embodiment of this application.

[0051] Figure 4 The diagram shown is a schematic flowchart of a reference carbon mass data update process as described in an embodiment of this application.

[0052] Figure 5 The diagram shown is a simulation illustration illustrating the relationship between methane cracking and the molten medium as described in an embodiment of this application.

[0053] Figure 6 The diagram shown is a simulation illustration illustrating the relationship between methane cracking and temperature change as described in an embodiment of this application.

[0054] Figure 7 The diagram shown is a simulation illustration of the relationship between methane cracking and pressure changes as described in an embodiment of this application.

[0055] Figure 8 The diagram shown is a flow chart of a methane catalytic cracking simulation system as described in an embodiment of this application.

[0056] Figure 9 The diagram shown is a structural schematic of a terminal as described in an embodiment of this application.

[0057] Explanation of reference numerals in the attached figures

[0058] 41: Flow field acquisition module; 42: Response acquisition module; 43: Convergence judgment module; 50: Terminal; 51: Processor; 52: Memory; 521: Operating system; 522: Application program; 53: User interface; 54: Network interface; 55: Bus system. Detailed Implementation

[0059] The following specific examples illustrate the implementation of this application. Those skilled in the art can easily understand other advantages and effects of this application from the content disclosed in this specification. This application can also be implemented or applied through other different specific embodiments, and various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of this application. It should be noted that, unless otherwise specified, the following embodiments and features in the embodiments can be combined with each other.

[0060] It should be noted that the illustrations provided in the following embodiments are only schematic representations of the basic concept of this application. Therefore, the drawings only show the components related to this application and are not drawn according to the actual number, shape and size of the components in the actual implementation. In the actual implementation, the form, quantity and proportion of each component can be arbitrarily changed, and the layout of the components may also be more complex.

[0061] Existing methods for simulating methane catalytic cracking typically simplify carbon particles to gaseous components or ignore carbon particles altogether, resulting in low accuracy in simulating methane catalytic cracking.

[0062] To address the technical problems existing in the prior art, the following embodiments of this application provide a methane catalytic cracking simulation method, system, terminal, and medium. By simulating the mass and movement of carbon particles, the catalytic cracking of methane is made more accurate, which helps to better reflect the methane catalytic cracking situation and thus assists in designing and obtaining an efficient and low-energy-consumption methane catalytic cracking hydrogen production scheme.

[0063] The technical solutions in the embodiments of this application will be described in detail below with reference to the accompanying drawings.

[0064] This embodiment provides a method for simulating methane catalytic cracking, used to simulate the methane catalytic cracking reaction within a target region. The target region is a user-pre-defined extracted area for simulation. For example, the target region is a two-dimensional region of a preset size; for instance, the width of the target region is a preset width, such as within the range of 20-50 mm, and the height of the target region is the height of the reaction vessel undergoing the methane catalytic cracking reaction; for example, if the reaction vessel height is 200 mm, then the height of the target region is 200 mm.

[0065] It should be noted that, in order to facilitate the subsequent simulation process, the width of the target area is usually no less than 6 times the initial diameter of the bubble. For example, if the initial diameter of the bubble is set to any value in the range of 2-8mm, then the width of the target area is usually no less than the range of 12-48mm.

[0066] Furthermore, to facilitate calculation, the target region is divided into individual grids based on a preset grid size, allowing for discretization analysis based on each grid, thereby simulating the methane catalytic cracking process. Those skilled in the art can set the specific grid size according to actual needs; this embodiment does not impose specific limitations. For example, the grid size is 0.25mm * 0.25mm.

[0067] The methane catalytic cracking simulation method is essentially an iterative simulation of the methane catalytic cracking reaction process within each target region. The iterative simulation is a cyclical simulation process, with the time interval between each simulation process being a preset time step. For example, the time step is... That is, the time interval between each simulation process is .

[0068] Based on this, at each time interval, the simulation process is re-executed based on the results of the previous simulation, to simulate the methane catalytic cracking reaction in an iterative manner and obtain the reaction results. Specifically, as follows... Figure 1 As shown, a single simulation includes:

[0069] S100 performs flow field calculations based on the reference flow field and the reference reaction to obtain updated flow field conditions.

[0070] Among them, the flow field condition refers to the changes in the motion and thermodynamic state of the gas and liquid during the catalytic cracking reaction of methane.

[0071] Specifically, the reference flow field is the flow field calculated and obtained during the simulation of the previous time step, used to characterize the flow field at the end of the previous time step and at the beginning of the current time step. The reference flow field includes reference pressure distribution data, reference velocity distribution data, reference temperature distribution data, reference volume fraction data, and a reference gas-liquid interface.

[0072] The updated flow field conditions are the flow field conditions calculated during the simulation at the current time step, used to characterize the flow field conditions at the end of the current time step. This includes updating pressure distribution data, velocity distribution data, temperature distribution data, volume fraction data, and the gas-liquid interface.

[0073] It should be noted that during the initial execution of the simulation process, i.e., the simulation process corresponding to the first time step, the reference flow field used for updating is a preset flow field condition. Specifically, during the initial execution of the simulation process:

[0074] The reference pressure distribution data is set as follows: a pressure distribution based on a preset linear distribution of top pressure value, that is, the pressure at the top of the target area is a preset value, for example, 2 bar, and increases linearly with the depth of the target area, wherein the pressure increase gradient is the product of density and gravitational acceleration.

[0075] The reference velocity distribution data is set to 0.

[0076] The reference volume fraction data and reference gas-liquid interface are set as follows: a bubble of a preset size is generated at the bottom center of the target area. The reference volume fraction data and the reference gas-liquid interface are calculated based on this bubble. Those skilled in the art can set the preset size of the bubble according to actual needs. This embodiment does not impose specific limitations. For example, the bubble is circular, and its preset size is a diameter of 2-8 mm.

[0077] The reference temperature distribution data is set as follows: based on the pre-defined reference gas-liquid interface, the target region is divided into a gas region and a liquid region; the gas region is located inside the generated bubble, and the temperature value of all grids within this gas region is set to a preset gas temperature value; the liquid region is located outside the generated bubble, and the temperature value of all grids within this liquid region is set to a preset liquid temperature value. Those skilled in the art can set the gas and liquid temperature values ​​according to actual needs, and this embodiment does not impose specific limitations. It should be noted that since the methane catalytic cracking reaction requires a certain temperature to proceed, the gas and liquid temperature values ​​are set to relatively high values ​​to avoid a large number of time steps during the simulation where heating does not occur, thereby reducing the computational load. For example, the preset gas temperature value is 500 degrees Celsius, and the preset liquid temperature is 1000 degrees Celsius.

[0078] The reaction status refers to the changes in the state and composition of matter caused by chemical reactions during the catalytic cracking reaction of methane.

[0079] Specifically, the reference reaction conditions are the reaction conditions calculated during the simulation of the previous time step, used to characterize the reaction conditions at the end of the previous time step and the beginning of the current time step. These reference reaction conditions include reference component concentration distribution data, reference reaction rate data, reference carbon mass data, and reference carbon motion data. The reference component concentration distribution data includes reference methane concentration distribution data and reference hydrogen concentration distribution data; the reference carbon motion data includes reference carbon position data and reference carbon velocity data.

[0080] It should be noted that, similar to the reference flow field situation, during the initial execution of the simulation process, i.e., the simulation process corresponding to the first time step, the reference response situation used for updating is the preset response situation. Specifically, during the initial execution of the simulation process:

[0081] The reference component concentration distribution data is set as follows: the methane concentration distribution data is a uniform distribution based on a preset methane concentration value. Those skilled in the art can set the preset methane concentration value according to actual needs, and this embodiment does not impose specific limitations here; the reference hydrogen concentration distribution data is 0.

[0082] The reference reaction rate data was set to 0.

[0083] The method for setting the reference carbon mass data and reference carbon motion data is as follows: a carbon particle with a preset diameter is set within any grid, that is, a carbon particle with a preset diameter is randomly selected at any position to generate a carbon particle with a preset diameter. Carbon mass data is then obtained based on the preset diameter, and reference carbon position data is obtained based on this grid. The carbon velocity data is set to 0. Those skilled in the art can set the preset diameter according to actual needs; this embodiment does not impose specific limitations. For example, the preset diameter is 1×10⁻⁻⁻⁴. 9 m.

[0084] Based on the reference flow field, flow field calculations are performed using the laws of fluid mechanics to update the flow field and obtain updated flow field conditions. This allows for the simulation of the motion and thermodynamic states of gas and liquid during the methane catalytic cracking reaction, facilitating the analysis of the methane catalytic cracking reaction.

[0085] In some alternative implementations, such as Figure 2 As shown, the flow field calculation includes:

[0086] S101, extract reference velocity distribution data and reference pressure distribution data, calculate intermediate velocity distribution data using the law of conservation of momentum; perform pressure correction based on intermediate velocity distribution data to obtain updated pressure distribution data; and perform velocity correction based on updated pressure distribution data to obtain updated velocity distribution data.

[0087] The law of conservation of momentum states that the momentum of any given grid is affected by gravity, liquid surface tension and liquid pressure, as well as the momentum transferred from all surrounding grids through convection and diffusion.

[0088] In some alternative implementations, for any grid, the formula for calculating the law of conservation of momentum is expressed as:

[0089]

[0090] in, The calculated grid, i.e., the mesh The intermediate speed, For grid The intermediate velocity of each surrounding grid; For grid The convective diffusion coefficient, For grid The convective diffusion coefficient of each surrounding grid cell can be specifically determined by those skilled in the art based on practical experience. The convection diffusion coefficient of the grid and its surrounding grid is not specifically limited in this embodiment; It is used to reflect the effect of liquid pressure. It should be noted that... This is the reference pressure distribution data obtained at the previous time step; Used to reflect the effects of gravity and liquid surface tension, specifically, , To reflect the effect of gravity, For grid density, , For grid The gas volume fraction in the previous calculation, This represents the density of the gas during the catalytic cracking reaction of methane. This represents the density of the liquid during the catalytic cracking reaction of methane. It is the acceleration due to gravity. For grid Regarding the volume, it should be noted that if the target area is a two-dimensional region, then the mesh... Also a two-dimensional grid, in this case, let the grid... Thickness in units of 1, computational mesh The product of the area and the thickness is the mesh. volume, Used to reflect the effect of liquid surface tension For grid The curvature of the gas-liquid interface can be calculated by those skilled in the art based on a reference gas-liquid interface to obtain the mesh. gas-liquid interface curvature , The gas-liquid interface surface tension coefficient can be set by those skilled in the art based on practical experience, and no specific limitation is made in this embodiment.

[0091] By constructing the above momentum conservation law calculation formula for each grid, the velocity situation of each grid can be calculated, and thus the intermediate velocity distribution data can be obtained.

[0092] It should be noted that the calculated intermediate velocity distribution data actually represents the velocity situation at the end of the current time step, while The reference pressure distribution data obtained in the previous time step corresponds to the end time of the previous time step. Based on this, for the intermediate velocity distribution data obtained through calculation, the pressure distribution needs to be corrected, the updated pressure distribution data corresponding to the end time of the current time step is obtained, and the intermediate velocity distribution data is corrected based on the updated pressure distribution data to obtain the updated velocity distribution data.

[0093] In some alternative implementations, the formula for calculating pressure correction based on intermediate velocity distribution data is as follows:

[0094]

[0095] in, This is the pressure correction amount.

[0096] Based on this, pressure correction values ​​are obtained by performing pressure correction on the intermediate velocity distribution data, and then the pressure distribution data is calculated and updated. Specifically:

[0097]

[0098] in, To update the pressure distribution data.

[0099] Furthermore, velocity correction is performed based on the updated pressure distribution data to obtain the corrected velocity, which is then used as the updated velocity distribution data. Specifically,

[0100]

[0101] in, To update the velocity distribution data.

[0102] S102, extract the reference reaction rate data, and use the law of energy conservation to calculate and update the temperature distribution data.

[0103] The law of energy conservation states that for any given grid, the sum of the energy required for a temperature change within the grid and the energy lost through convection and diffusion is equal to the sum of the energy gained through heat conduction at the grid boundaries and the energy absorbed by the catalytic cracking of methane. It should be noted that the energy absorbed by the catalytic cracking of methane in the current time step is obtained based on the reaction rate, exemplarily the rate of methane consumption.

[0104] In some alternative implementations, for any grid, the energy conservation law is expressed as follows:

[0105]

[0106] in, For grid Specific heat capacity, , This represents the specific heat capacity of the gas during the catalytic cracking reaction of methane. This refers to the specific heat capacity of the liquid during the catalytic cracking reaction of methane. The mesh obtained at the previous time step temperature, The grid updated for the current time step Temperature; For the previous time step, this grid The temperature of each boundary, specifically, is obtained based on the reference temperature distribution data acquired at the previous time step for that grid. The average of the temperatures on both sides of each boundary is taken as the temperature of the corresponding boundary. This grid is the current time step. The velocities at each boundary are specifically obtained based on the updated velocity distribution data for that grid. The velocity on both sides of each boundary is used as the average of the velocities on both sides as the temperature of the corresponding boundary. This grid is the current time step. The direction vectors of each boundary; for this grid The thermal conductivity of each boundary can be set by those skilled in the art based on practical experience; this embodiment does not impose specific limitations here. For this grid The methane consumption rate is specifically obtained based on the reference reaction rate data obtained in the previous time step. It is the enthalpy constant, specifically, .

[0107] Based on this, the updated temperature distribution data is obtained by substituting the reference reaction rate data, reference temperature distribution data, and updated rate distribution data into the calculation formula corresponding to the law of energy conservation.

[0108] S103, based on the updated velocity distribution data, calculate the updated volume fraction data using the gas-liquid transport law; based on the updated volume fraction data, reconstruct the interface to obtain the updated gas-liquid interface.

[0109] Among them, the gas-liquid transport law means that for any grid, the change in gas volume within the grid is equal to the sum of the fluxes of diffusion and convection across the boundaries of that grid.

[0110] In some optional implementations, for any grid, the calculation formula corresponding to the gas-liquid transport law is expressed as:

[0111]

[0112] in, For the current time step grid Volume fraction; For grid The volume fraction of each boundary is specifically obtained by taking the average value of the grid volume fraction on both sides of each interface based on the reference volume fraction data, and then using it as the volume fraction of the corresponding boundary.

[0113] Based on this, by substituting the reference volume fraction data and update velocity fraction data into the calculation formula corresponding to the gas-liquid transport law, the volume fraction corresponding to each grid is calculated to obtain the updated volume fraction data.

[0114] Furthermore, based on the updated volume fractions calculated for each grid, interface reconstruction is performed to obtain the updated gas-liquid interface. Specifically, for each grid containing an interface (i.e., a grid with a volume fraction other than 1 or 0), the volume fraction gradient is calculated based on the volume fractions of that grid and its adjacent grids to determine the interface normal. Then, based on the volume fraction of that grid, the interface is located as a straight line segment passing through that grid, and this straight line segment shares endpoints with the gas-liquid interfaces within its adjacent grids. Finally, several sets of connected line segments are obtained as the updated gas-liquid interface. For example, the Geo-Reconstruct format is used to reconstruct the interface based on the volume fractions of adjacent grids, thereby obtaining the updated gas-liquid interface.

[0115] S104 outputs updated flow field data, including updated pressure distribution data, updated velocity distribution data, updated temperature distribution data, updated volume fraction data, and updated gas-liquid interface data.

[0116] Based on this, the updated flow field situation corresponding to the current time step can be obtained.

[0117] S200 obtains updated reaction information based on updated flow field conditions and reference reaction conditions, combined with the laws of material transport and reaction.

[0118] Among them, the material transport reaction law refers to the reaction changes of components such as methane, hydrogen, carbon particles, and catalyst during the catalytic cracking reaction of methane.

[0119] The reaction status refers to the changes in the material state and composition of each component during the catalytic cracking reaction of methane.

[0120] The updated reaction status refers to the reaction information calculated during the simulation at the current time step, used to characterize the reaction status at the end of the current time step. This updated reaction status includes updated component concentration distribution data, updated reaction rate data, updated carbon mass data, and updated carbon motion data. Updated component concentration distribution data includes updated methane concentration distribution data and updated hydrogen concentration distribution data; updated carbon motion data includes updated carbon position data and updated carbon velocity data.

[0121] In some alternative implementations, such as Figure 3 As shown, by understanding the laws governing mass transport and reaction, updated reaction information can be obtained, including:

[0122] S210: Extract the reference component concentration distribution data and the current temperature distribution data, and calculate and update the reaction rate data using the laws governing methane catalytic cracking reaction.

[0123] Among them, the methane catalytic cracking reaction law is used to characterize the methane cracking process under the influence of temperature above the activation energy barrier and concentration and the presence or absence of a catalyst.

[0124] Specifically, for any given grid cell, based on the updated volume fraction data, it is determined whether a catalytic reaction has occurred within that grid cell. If the gas volume within that grid cell is small, a catalytic reaction has occurred; if the gas volume within that grid cell is large, a non-catalytic reaction has occurred. For example, using 0.5 as a boundary, if the volume fraction of any grid cell is less than 0.5, a catalytic reaction has occurred within that grid cell; if the volume fraction of any grid cell is not less than 0.5, a non-catalytic reaction has occurred within that grid cell.

[0125] In some optional implementations, the updated reaction rate data is calculated for any given grid as follows:

[0126] when hour:

[0127]

[0128]

[0129]

[0130] in, This represents the non-catalytic reaction rate. The apparent rate constant for non-catalytic reactions, For non-catalytic reactions, specifically, those skilled in the art can set the pre-factor of non-catalytic reactions based on experience. This embodiment does not impose specific limitations here; This is the activation energy for non-catalytic reactions; The molar concentration of methane within this grid. The molar concentration of hydrogen in this grid is specifically extracted based on the reference component concentration distribution data, which is the molar concentration of methane and hydrogen in this grid. The equilibrium constant is expressed in terms of concentration. Let be the equilibrium constant expressed in terms of partial pressures, i.e., It is the thermodynamic equilibrium constant; The universal gas constant is usually taken as . ; The standard pressure is typically 1 bar.

[0131] when hour:

[0132]

[0133]

[0134]

[0135]

[0136] in, For the catalytic reaction rate, This represents the intrinsic reaction rate; This refers to the gas-liquid contact area factor, specifically... ; Let be the apparent rate constant of the catalytic reaction. As a pre-factor for the catalytic reaction, specifically, those skilled in the art can set the pre-factor for the catalytic reaction based on experience. This embodiment does not impose specific limitations here; This is the activation energy of the catalytic reaction.

[0137] Based on the above calculation method for updating reaction rate data, the reaction rate corresponding to each grid is calculated to obtain the updated reaction rate data.

[0138] S220, based on updated reaction rate data, combined with reference component concentration distribution data and updated rate distribution data, uses component transport laws to calculate and obtain updated component concentration distribution data.

[0139] Among them, the component transport law means that for any grid, the change in the mass of each component within the grid and the sum of the mass of each component within the grid through convection are equal to the sum of the mass of each component generated or consumed through reaction within the grid and the mass of each component leaving the grid through diffusion.

[0140] In some optional implementations, for any grid, the calculation formula corresponding to the component transport law is expressed as:

[0141]

[0142] in, For the current time step grid Internal components The mass concentration, specifically for methane, For the current time step grid The mass concentration of methane in the inner layer, For the current time step grid The mass concentration of hydrogen gas inside. For the previous time step grid Internal components The mass concentration, specifically, For the previous time step grid The mass concentration of methane in the inner layer, For the previous time step grid The mass concentration of hydrogen gas inside; For a time step, Components In the grid The mass concentration at each boundary, specifically, methane in the grid Mass concentration at each boundary Hydrogen gas in the grid Mass concentration at each boundary; Components The diffusion coefficient, specifically, Let be the diffusion coefficient of methane. The diffusion coefficient of hydrogen is given. It should be noted that those skilled in the art can set the diffusion coefficients of each component based on practical experience, and this embodiment does not impose any specific limitations. Components In the grid The mass concentration gradient at each boundary, specifically, methane in the grid Mass concentration gradients at each boundary, Hydrogen gas in the grid Mass concentration gradients at each boundary; Components The generation rate or consumption rate, specifically, The rate of methane consumption. The hydrogen production rate is given by the grid; it should be noted that this grid is obtained based on updated reaction rate data. The reaction rate is determined, and the generation or consumption rate of each component is obtained based on this reaction rate.

[0143] Based on this, the component transport equations for methane and hydrogen are constructed according to the calculation formulas corresponding to the component transport laws, and the mass concentrations of methane and hydrogen in each grid at the current time step are solved respectively, so as to obtain updated methane concentration data and updated hydrogen concentration data, and then obtain updated component concentration distribution data.

[0144] S230 updates the reference carbon mass data based on the updated reaction rate data and calculates the updated carbon mass data.

[0145] Specifically, during the catalytic cracking reaction of methane, carbon particles are generated along with hydrogen gas, thereby changing the mass of the carbon particles.

[0146] Based on this, during the simulation process, the reaction rate of methane catalytic cracking is obtained based on the updated reaction rate data, thereby obtaining the carbon particle generation rate, i.e. the carbon particle mass change rate, and then the reference carbon mass data is updated to obtain updated carbon mass data.

[0147] In some alternative implementations, such as Figure 4 As shown, the process of updating the reference carbon mass data includes:

[0148] S231, Obtain reference carbon particle diameter data, and combine it with the updated reaction rate data to obtain updated carbon particle diameter data.

[0149] Specifically, for any given grid, the volume of the carbon particles is obtained based on reference carbon mass data and the density of the carbon particles. Those skilled in the art should be familiar with the density of carbon particles; therefore, this embodiment will not elaborate on it.

[0150] By treating carbon particles as spherical particles, their diameter can be calculated based on their volume. This allows us to obtain reference carbon particle diameter data based on the diameter of the carbon particles in each grid.

[0151] In some optional implementations, the updated carbon particle diameter data is calculated based on the reference carbon particle diameter data and the updated reaction rate data. The calculation method is as follows:

[0152]

[0153] in, For the current time step grid The diameter of the internal carbon particles, For the previous time step grid The diameter of the internal carbon particles, To determine the carbon particle formation rate, specifically, based on updated reaction rate data, a grid is obtained. The corresponding reaction rate of methane catalytic cracking is determined, and based on the reaction rate of methane catalytic cracking, the carbon particle formation rate is obtained. Specifically, ,in Let P be the volume of the current grid. This represents the number of particles in the current grid.

[0154] Based on this, the diameter change of carbon particles in each grid is obtained through the above-mentioned carbon particle diameter update calculation method, which is used as the updated carbon particle diameter data.

[0155] S232, based on the updated carbon particle diameter data, calculate and obtain updated carbon mass data.

[0156] Specifically, based on the diameter of the carbon particles in each grid, the carbon particles are treated as spherical particles. The volume of the updated carbon particles is calculated, and combined with the density of the carbon particles, the mass of the updated carbon particles is further calculated to obtain the carbon particle mass corresponding to each grid after the update, which is used as the updated carbon mass data.

[0157] S240 analyzes the forces acting on carbon particles based on updated carbon mass data and updated velocity distribution data to update the reference carbon motion data and obtain updated carbon motion data.

[0158] The carbon particle motion data includes carbon particle position data and carbon particle velocity data. In other words, the reference carbon motion data includes reference carbon particle position data and reference carbon particle velocity data, and the updated carbon motion data includes updated carbon particle position data and updated carbon particle velocity data.

[0159] Specifically, carbon particles in the target region are affected by fluid drag and gravity, causing them to move within the target region. Based on this, by acquiring the drag force acting on the carbon particles, their motion can be analyzed, thereby obtaining updated carbon motion data.

[0160] In some optional implementations, the carbon particle motion data is updated in the following way:

[0161]

[0162]

[0163]

[0164]

[0165] in, For the current time step grid The velocity of medium carbon particles, For the previous time step grid The velocity of medium carbon particles, For the current time step grid The mass of medium carbon particles, specifically, is extracted from the grid based on updated carbon mass data. The mass of medium carbon particles, For the current time step grid The position of the medium carbon particles For the previous time step grid The position of the medium carbon particles The drag force brought to the fluid. This is the drag coefficient. The drag force is the net volume force acting on the carbon particles, consisting of the difference between gravity and buoyancy. Specifically, those skilled in the art can set the drag coefficient based on experience; this embodiment does not impose specific limitations. The local fluid density is used as the grid density in actual execution. The density, that is , For the area facing the flow, The acceleration due to gravity is taken as 9.8 m / s² in this embodiment, with the direction along the −y direction; The mass of the fluid displaced by the carbon particles corresponds to the magnitude of the buoyant force experienced by the particles in a static fluid. Specifically, a grid is extracted based on updated carbon mass data. The mass of the carbon particles is used to calculate the diameter of the carbon particles in order to calculate the frontal area. For example, .

[0166] Based on this, the drag force on the carbon particles is calculated, and the velocity of the carbon particles is updated based on the drag force, thereby updating the position of the carbon particles. Finally, the position of the carbon particles is updated, and the updated carbon particle position data and updated carbon particle velocity data are obtained, thereby obtaining updated carbon motion data.

[0167] S300 performs convergence verification based on the reference flow field and the updated flow field, as well as the reference response and the updated response. If the convergence verification is successful, the updated flow field and the updated response are output, and the iterative simulation is stopped. Otherwise, the updated flow field is used as the new reference flow field and the updated response is used as the new reference response, and the simulation process is re-executed.

[0168] Convergence verification refers to whether the updated flow field and reaction conditions, compared with the reference flow field and reaction conditions obtained in the previous time step, do not exceed the maximum degree of change, so that the methane catalytic cracking process is in equilibrium.

[0169] Based on this, when the updated flow field and updated reaction conditions pass the convergence verification, the iterative cycle of the simulation process is stopped, and the result of the last round of simulation is output as the final result of the simulation.

[0170] In some alternative implementations, convergence is determined based on the changes in volume fraction, concentration of each component, and carbon movement in the current simulation and the previous simulation.

[0171] Specifically, volume fraction residual data is obtained based on reference volume fraction data and updated volume fraction data. The volume fraction residual data includes the difference in volume fraction between each grid cell at the current time step and the previous time step.

[0172] Based on the reference and updated component concentration distribution data, residual component concentration data are obtained. The component concentrations include methane and hydrogen concentrations. Based on this, residual methane concentration data is obtained, based on the reference and updated methane concentration distribution data. This residual methane concentration data includes the difference in methane concentration between the current and previous time steps for each grid cell. Similarly, residual hydrogen concentration data is obtained, based on the reference and updated hydrogen concentration distribution data. This residual hydrogen concentration data includes the difference in hydrogen concentration between the current and previous time steps for each grid cell.

[0173] Carbon motion residual data is obtained based on reference carbon motion data and updated carbon motion data. The carbon motion includes carbon particle position and carbon particle velocity. Based on this, carbon particle position residual data is obtained based on reference carbon particle position data and updated carbon particle position data. This residual data includes the distance between the current and previous time steps for the same carbon particle. Similarly, carbon particle velocity residual data is obtained based on reference carbon particle velocity distribution data and updated carbon particle velocity distribution data. This residual data includes the difference in velocity between the current and previous time steps for the same carbon particle.

[0174] If the volume fraction residual data does not exceed the preset upper limit of the volume fraction residual, that is, the maximum value of the volume fraction difference in each grid does not exceed the preset volume fraction difference threshold; if the component concentration residual data does not exceed the preset upper limit of the component concentration residual, that is, the maximum value of the methane concentration difference in each grid does not exceed the preset methane concentration difference threshold, and the maximum value of the hydrogen concentration difference in each grid does not exceed the preset hydrogen concentration difference threshold; if the carbon motion residual data does not exceed the preset upper limit of the carbon motion residual, that is, the maximum value of the carbon particle position distance does not exceed the preset distance threshold, and the maximum value of the carbon particle velocity difference does not exceed the preset velocity difference threshold; only when all of the above requirements are met will the convergence verification pass; otherwise, the convergence verification will fail.

[0175] It should be noted that those skilled in the art can set specific values ​​for the volume fraction difference threshold, methane concentration difference threshold, hydrogen concentration difference threshold, distance threshold, and speed difference threshold according to actual needs, and this embodiment does not impose specific limitations here.

[0176] Based on this, the catalytic cracking process of methane can be simulated, and the flow field and reaction conditions of its stable reaction can be output, allowing researchers to obtain information about the catalytic cracking of methane. For details, please refer to [link to relevant documentation]. Figures 5 to 7 The simulation results obtained based on the method of this application show the catalytic cracking of methane under different melting media, temperatures, and pressures.

[0177] like Figure 8 As shown in the figure, this embodiment provides a methane catalytic cracking simulation system for iteratively simulating the methane catalytic cracking reaction in a target area and cyclically simulating the process. Specifically, the methane catalytic cracking simulation system includes a flow field acquisition module 41, a reaction acquisition module 42, and a convergence judgment module 43.

[0178] The flow field acquisition module 41 is used to perform flow field calculations based on the reference flow field and the reference reaction to obtain updated flow field conditions. The reference flow field conditions include reference pressure distribution data, reference velocity distribution data, reference temperature distribution data, reference volume fraction data, and reference gas-liquid interface data. The updated flow field conditions include updating pressure distribution data, updating velocity distribution data, updating temperature distribution data, updating volume fraction data, and updating gas-liquid interface data.

[0179] The reaction status acquisition module 42 is used to acquire updated reaction status based on the updated flow field status and the reference reaction status, combined with the laws of mass transport and reaction. The reference reaction status includes reference component concentration distribution data, reference reaction rate data, reference carbon mass data, and reference carbon motion data; the updated reaction status includes updated component concentration distribution data, updated reaction rate data, updated carbon mass data, and updated carbon motion data.

[0180] The convergence judgment module 43 is used to perform convergence verification based on the reference flow field and the updated flow field, as well as the reference reaction and the updated reaction. If the convergence verification is successful, the updated flow field and the updated reaction are output, and the iterative simulation is stopped. Otherwise, the updated flow field is used as the new reference flow field and the updated reaction is used as the new reference reaction, and the process returns to the flow field acquisition module 41 and the reaction acquisition module 42 to re-execute the next simulation process.

[0181] Based on the same technical concept, the methane catalytic cracking simulation method provided in this embodiment of the invention can be implemented on the terminal side or the server side.

[0182] like Figure 9 The diagram shown illustrates an optional hardware structure of a terminal according to an embodiment of the present invention. This terminal may be a mobile phone, computer device, tablet device, personal digital processing device, factory back-end processing device, etc. The terminal includes at least one processor 51, a memory 52, at least one network interface 54, and a user interface 53. The various components in the device are coupled together via a bus system 55. It is understood that the bus system 55 is used to implement communication between these components. In addition to a data bus, the bus system 55 also includes a power bus, a control bus, and a status signal bus.

[0183] The user interface 53 may include a monitor, keyboard, mouse, trackball, clicker, button, touchpad, or touch screen.

[0184] It is understood that memory 52 can be volatile memory or non-volatile memory, or both. Non-volatile memory can be read-only memory (ROM) or programmable read-only memory (PROM), used as an external cache. By way of example, but not limitation, many forms of RAM are available, such as static random access memory (SRAM) and synchronous static random access memory (SSRAM). The memory characterized in the embodiments of the present invention is intended to include, but is not limited to, these and any other suitable categories of memory.

[0185] In this embodiment of the invention, the memory 52 is used to store various types of data to support the operation of the terminal. Examples of this data include: any executable program for operation on the terminal*, such as the operating system 521 and application programs 522; the operating system 521 contains various system programs, such as the framework layer, core library layer, driver layer, etc., for implementing various basic services and handling hardware-based tasks. Application programs 522 may contain various applications, such as media players, browsers, etc., for implementing various application services. The methane catalytic cracking simulation method provided in this embodiment of the invention can be included in application program 522.

[0186] The methods disclosed in the above embodiments of the present invention can be applied to processor 51, or implemented by processor 51. Processor 51 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by the integrated logic circuit of the hardware in processor 51 or by instructions in the form of software. The processor mentioned above may be a general-purpose processor, a digital signal processor (DSP), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Processor 51 can implement or execute the methods, steps and logic block diagrams disclosed in the embodiments of the present invention. Processor 51 may be a microprocessor or any conventional processor, etc. The steps of the methane catalytic cracking simulation method provided in the embodiments of the present invention can be directly manifested as being executed by a hardware decoding processor, or being executed by a combination of hardware and software modules in the decoding processor. The software modules may be located in a storage medium, which is located in a memory. The processor reads the information in the memory and combines it with its hardware to complete the steps of the aforementioned method.

[0187] In an exemplary embodiment, the terminal* may be used by one or more application-specific integrated circuits (ASICs), DSPs, programmable logic devices (PLDs), or complex programmable logic devices (CPLDs) to execute the aforementioned method.

[0188] This invention also provides a computer-readable storage medium storing a computer program that, when invoked by a processor, implements the methane catalytic cracking simulation method provided by this invention.

[0189] Computer-readable storage media can be tangible devices capable of holding and storing instructions used by an instruction execution device. Computer-readable storage media can be, for example, (but not limited to) electrical storage devices, magnetic storage devices, optical storage devices, electromagnetic storage devices, semiconductor storage devices, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of computer-readable storage media include: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), static random access memory (SRAM), portable compact disc read-only memory (CD-ROM), digital multifunction disc (DVD), memory sticks, floppy disks, and mechanical encoding devices.

[0190] The computer-readable program represented herein can be downloaded from a computer-readable storage medium to various computing / processing devices, or downloaded via a network, such as the Internet, a local area network, a wide area network, and / or a wireless network, to an external computer or external storage device. A network adapter card or network interface in each computing / processing device receives computer-readable program instructions from the network and forwards these instructions to the computer-readable storage medium in the respective computing / processing device.

[0191] In summary, this application simulates the generation, mass change, and motion of carbon particles during the methane catalytic cracking reaction, thereby effectively improving the accuracy of the simulation and facilitating the design of an efficient and low-energy-consumption methane catalytic cracking hydrogen production scheme.

[0192] The descriptions of the processes or structures corresponding to the above figures each have their own emphasis. For parts of a process or structure that are not described in detail, please refer to the relevant descriptions of other processes or structures.

[0193] The above embodiments are merely illustrative of the principles and effects of this application and are not intended to limit this application. Any person skilled in the art can modify or alter the above embodiments without departing from the spirit and scope of this application. Therefore, all equivalent modifications or alterations made by those skilled in the art without departing from the spirit and technical concept disclosed in this application should still be covered by the claims of this application.

Claims

1. A method for simulating methane catalytic cracking, used to simulate methane catalytic cracking reactions in a target region, characterized in that, The catalytic cracking reaction process of methane in each of the target regions was iteratively simulated; The iterative simulation is a cyclically executed simulation process, with a preset time step between each simulation process; each simulation process includes: Based on the reference flow field and the reference reaction, flow field calculations are performed to obtain updated flow field conditions. The reference flow field conditions include reference pressure distribution data, reference velocity distribution data, reference temperature distribution data, reference volume fraction data, and reference gas-liquid interface data. The updated flow field conditions include updated pressure distribution data, updated velocity distribution data, updated temperature distribution data, updated volume fraction data, and updated gas-liquid interface data. Based on the updated flow field and the reference reaction conditions, and combined with the laws governing mass transport and reaction, updated reaction conditions are obtained; wherein, the reference reaction conditions include reference component concentration distribution data, reference reaction rate data, reference carbon mass data, and reference carbon motion data; the updated reaction conditions include updated component concentration distribution data, updated reaction rate data, updated carbon mass data, and updated carbon motion data; Based on the reference flow field and the updated flow field, as well as the reference response and the updated response, convergence verification is performed. If the convergence verification passes, the updated flow field and the updated response are output, and the iterative simulation is stopped. Otherwise, the updated flow field is used as the new reference flow field and the updated response is used as the new reference response, and the simulation process is re-executed.

2. The method according to claim 1, characterized in that, The process of calculating the flow field based on the reference flow field and the reference reaction to obtain an updated flow field includes: Extract the reference velocity distribution data and the reference pressure distribution data, and calculate the intermediate velocity distribution data using the law of conservation of momentum; perform pressure correction based on the intermediate velocity distribution data to obtain the updated pressure distribution data; and perform velocity correction based on the updated pressure distribution data to obtain the updated velocity distribution data. Extract the reference reaction rate data and calculate the updated temperature distribution data using the law of energy conservation; Based on the updated velocity distribution data, the updated volume fraction data is calculated using the gas-liquid transport law; based on the updated volume fraction data, the interface is reconstructed to obtain the updated gas-liquid interface; The updated pressure distribution data, the updated velocity distribution data, the updated temperature distribution data, the updated volume fraction data, and the updated gas-liquid interface data are output as the updated flow field conditions.

3. The method according to claim 1, characterized in that, The process of obtaining the updated reaction information based on the updated flow field and the reference reaction information, combined with the laws governing mass transport and reaction, includes: Extract the reference component concentration distribution data and the updated temperature distribution data, and calculate the updated reaction rate data using the laws governing methane catalytic cracking reaction; Based on the updated reaction rate data, combined with the reference component concentration distribution data and the updated rate distribution data, the updated component concentration distribution data is calculated using the component transport laws. Based on the updated reaction rate data, the reference carbon mass data is updated, and the updated carbon mass data is calculated and obtained. Based on the updated carbon mass data and the updated velocity distribution data, the stress on carbon particles is analyzed to update the reference carbon motion data and obtain the updated carbon motion data.

4. The method according to claim 3, characterized in that, The step of updating the reference carbon mass data based on the updated reaction rate data and calculating the updated carbon mass data includes: Obtain the reference carbon particle diameter data, and combine it with the updated reaction rate data to obtain the updated carbon particle diameter data; Based on the updated carbon particle diameter data, the updated carbon mass data is calculated and obtained.

5. The method according to claim 3, characterized in that, The reference carbon motion data includes reference carbon position data and reference carbon velocity data; the updated carbon motion data includes updated carbon position data and updated carbon velocity data; the step of analyzing the force on carbon particles based on the updated carbon mass data and the updated velocity distribution data to update the reference carbon motion data and obtain the updated carbon motion data includes: Based on the updated carbon mass data and the updated carbon velocity data, acceleration data is obtained; The reference carbon velocity data is updated based on the acceleration data to obtain the updated carbon velocity data; The reference carbon position data is updated based on the updated carbon velocity data to obtain the updated carbon position data; Based on the updated carbon velocity data and the updated carbon position data, the updated carbon motion data is obtained.

6. The method according to claim 1, characterized in that, The convergence verification includes: Based on the reference volume fraction data and the updated volume fraction data, obtain volume fraction residual data; based on the reference component concentration distribution data and the updated component concentration distribution data, obtain component concentration residual data; based on the reference carbon motion data and the updated carbon motion data, obtain carbon motion residual data. If the volume fraction residual data does not exceed the preset upper limit of the volume fraction residual data, the component concentration residual data does not exceed the preset upper limit of the component concentration residual data, and the carbon motion residual data does not exceed the preset upper limit of the carbon motion residual data, then the convergence verification is passed; otherwise, the convergence verification is not passed.

7. The method according to claim 1, characterized in that, During the first execution of the simulation process: Linear distribution calculation is performed based on a preset top pressure value to obtain the reference pressure distribution data; And / or, based on a preset speed value, obtain the reference speed distribution data; And / or, based on a preset reaction rate value, obtain the reference reaction rate data; And / or, based on preset bubble size, preset bubble shape and preset bubble position, the interface is reconstructed to obtain the reference volume fraction data and the reference gas-liquid interface; And / or, based on the reference gas-liquid interface, obtain the gas region and the liquid region; use a preset gas temperature value as the temperature value of the gas region and a preset liquid temperature value as the temperature value of the liquid region to obtain the reference temperature distribution data; And / or, based on the preset component concentration values ​​corresponding to each component, obtain the reference component concentration distribution data; And / or, based on a preset mass value, obtain the reference carbon mass data; And / or, based on a preset velocity value, combined with the gas region, obtain the reference carbon motion data.

8. A methane catalytic cracking simulation system for iteratively simulating methane catalytic cracking reactions within a target region, characterized in that, It includes a flow field acquisition module, a reaction acquisition module, and a convergence judgment module; The flow field acquisition module is used to perform flow field calculations based on a reference flow field and a reference reaction to obtain updated flow field conditions. The reference flow field conditions include reference pressure distribution data, reference velocity distribution data, reference temperature distribution data, reference volume fraction data, and a reference gas-liquid interface. The updated flow field conditions include updated pressure distribution data, updated velocity distribution data, updated temperature distribution data, updated volume fraction data, and updated gas-liquid interface. The reaction status acquisition module is used to acquire updated reaction status based on the updated flow field status and the reference reaction status, combined with the laws of mass transport and reaction; wherein, the reference reaction status includes reference component concentration distribution data, reference reaction rate data, reference carbon mass data, and reference carbon motion data; the updated reaction status includes updated component concentration distribution data, updated reaction rate data, updated carbon mass data, and updated carbon motion data; The convergence judgment module is used to perform convergence verification based on the reference flow field situation and the updated flow field situation, as well as the reference reaction situation and the updated reaction situation. If the convergence verification is passed, the updated flow field situation and the updated reaction situation are output, and the iterative simulation is stopped; otherwise, the updated flow field situation is used as the new reference flow field situation, and the updated reaction situation is used as the new reference reaction situation, and the process returns to the flow field situation acquisition module and the reaction situation acquisition module.

9. A terminal, characterized in that, include: A processor and a memory, wherein the memory and the processor are communicatively connected; The memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory to cause the terminal to perform the methane catalytic cracking simulation method as described in any one of claims 1 to 7.

10. A computer storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the methane catalytic cracking simulation method as described in any one of claims 1 to 7.

Images