Method and apparatus for predicting chemical form and chemical reaction of multiphase system

Abstract

The present application provides a method and an apparatus for predicting the chemical form and chemical reaction of a multiphase system. The method comprises: obtaining key conditions corresponding to a multiphase system, the key conditions comprising: temperature, pressure, multiple element types, and respective contents thereof; by using a Gibbs free energy minimization method, predicting to obtain chemical form information under the key conditions, the chemical form information comprising: various pure substance types and respective contents thereof when the multiphase system reaches thermodynamic equilibrium; and, on the basis of the chemical form information, predicting to obtain a chemical reaction that will occur in the multiphase system under the key conditions. The present application can improve the accuracy and efficiency of predicting the chemical reaction of the multiphase system on the basis of achieving the prediction of the chemical form of the elements in the multiphase system.

Description

Method and device for predicting chemical forms and chemical reactions of multiphase systems Technical Field

[0001] The present application relates to the technical field of multiphase systems in the fields of chemical engineering, energy, environment, etc., and in particular to a method and device for predicting the chemical forms and chemical reactions of a multiphase system. Background Art

[0002] The chemical form of an element affects its physical and chemical properties, including melting point, volatility, adsorption / desorption behavior, diffusion processes, and chemical reactions. In particular, the interactions between elements and materials (adsorption, desorption, diffusion, deposition, corrosion, etc.) depend on the chemical form of the element. Therefore, the study of element chemical form and chemical reactions is crucial in fields as diverse as nuclear energy systems, chemical engineering, geochemistry, and metallurgy, supporting advancements in other areas such as materials design, energy utilization, and environmental protection.

[0003] Currently, research on the chemical speciation of multiphase systems is primarily based on thermodynamic equilibrium calculations to determine the stable chemical speciation within the system. Some studies have used thermodynamic calculations to investigate the complex chemical reactions in water-cooled reactor nuclear fuel and to predict the chemical speciation of fission products during accidents at different oxygen potentials and temperatures. Using the example of a Cs-Sr-Ag-IO five-component system in the primary circuit of a high-temperature gas-cooled reactor (HTR-PM) under equilibrium core conditions, existing research has been able to calculate the chemical speciation of some fission products under both normal operating conditions and accident scenarios.

[0004] However, as shown in Figures 1 and 2, current methods are limited to software calculations, which analyze how the content of each element in a multiphase system changes with temperature. Traditional two-dimensional phase diagrams can describe systems with relatively simple element types, but when the system contains a large number of elements, existing phase diagrams cannot provide important information about all elements. Currently, no research has screened the chemical forms of elements for possible chemical reactions, and the description of chemical processes is unclear. No specific method has been proposed to predict the chemical forms and chemical reactions of elements. Summary of the Invention

[0005] In response to at least one problem in the prior art, the present application proposes a method and device for predicting the chemical forms and chemical reactions of a multiphase system, which can improve the accuracy and efficiency of predicting the chemical reactions of a multiphase system on the basis of predicting the chemical forms of elements in the multiphase system.

[0006] In order to solve the above technical problems, this application provides the following technical solutions:

[0007] In a first aspect, the present application provides a method for predicting chemical forms and chemical reactions of a multiphase system, comprising:

[0008] Obtaining key conditions corresponding to a multiphase system, including temperature, pressure, multiple element types and their respective contents;

[0009] Predicting chemical form information under the critical conditions using a Gibbs free energy minimization method, the chemical form information including: various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium;

[0010] Based on the chemical form information, the chemical reactions that will occur in the multiphase system under the key conditions are predicted.

[0011] In one embodiment, there are multiple key conditions, and the key influencing parameter of the multiphase system is at least one of the impurity element type among the multiple element types, the temperature, and the pressure; and the value of the key influencing parameter corresponding to each key condition is different;

[0012] Correspondingly, after obtaining the key conditions corresponding to the multiphase system, the following is also included:

[0013] The chemical form information under each key condition is predicted by the Gibbs free energy minimization method;

[0014] Based on the chemical form information under each key condition, the chemical reaction that will occur in the multiphase system under the key condition is predicted.

[0015] In one embodiment, after predicting the chemical form information under each key condition by the Gibbs free energy minimization method, the method further includes:

[0016] Selecting the type of pure substance and its content containing the preset selected element from the chemical form information under each key condition;

[0017] generating a two-dimensional phase diagram corresponding to the preset selected element according to the type and content of the pure substance containing the preset selected element;

[0018] The two-dimensional phase diagram is divided into multiple regions, each region corresponds to a set of pure substance groups, and each region corresponds to a different pure substance group. Each pure substance group is obtained based on the pure substance type containing the preset selected element. The value of the phase boundary between each region is the preset phase boundary value. The two-dimensional phase diagram also includes: a change curve of each pure substance type containing the selected element.

[0019] In one embodiment, selecting the pure substance type and content of the preset selected element from the chemical form information under each key condition includes:

[0020] Performing interpolation processing based on each key condition and the chemical form information corresponding thereto to obtain interpolated key conditions and the chemical form information corresponding thereto;

[0021] From the chemical form information under each interpolated key condition, the pure substance type and content of the preset selected element are selected.

[0022] In a second aspect, the present application provides a device for predicting chemical forms and chemical reactions of a multiphase system, comprising:

[0023] An acquisition module is used to obtain key conditions corresponding to the multiphase system, wherein the key conditions include temperature, pressure, multiple element types and their respective contents;

[0024] A first calculation module is configured to predict chemical form information under the key conditions using a Gibbs free energy minimization method, the chemical form information including various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium;

[0025] The first prediction module is used to predict the chemical reaction that will occur in the multiphase system under the key conditions based on the chemical form information.

[0026] In one embodiment, there are multiple key conditions, and the key influencing parameter of the multiphase system is at least one of the impurity element type among the multiple element types, the temperature, and the pressure; and the value of the key influencing parameter corresponding to each key condition is different;

[0027] Correspondingly, the device for predicting the chemical form and chemical reaction of a multiphase system further includes:

[0028] A second calculation module is used to predict the chemical form information under each key condition by using a Gibbs free energy minimization method;

[0029] The second prediction module is used to predict the chemical reaction that will occur in the multiphase system under each key condition based on the chemical form information under the key condition.

[0030] In one embodiment, the device for predicting chemical forms and chemical reactions of a multiphase system further comprises:

[0031] An extraction module, configured to select the type of pure substance containing a preset selected element and its content from the chemical form information under each key condition;

[0032] a phase diagram drawing module, configured to generate a two-dimensional phase diagram corresponding to the preset selected element according to the type and content of the pure substance containing the preset selected element;

[0033] The two-dimensional phase diagram is divided into multiple regions, each region corresponds to a set of pure substance groups, and each region corresponds to a different pure substance group. Each pure substance group is obtained based on the pure substance type containing the preset selected element. The value of the phase boundary between each region is the preset phase boundary value. The two-dimensional phase diagram also includes: a change curve of each pure substance type containing the selected element.

[0034] In one embodiment, the extraction module includes:

[0035] an interpolation processing unit, configured to perform interpolation processing based on each key condition and the chemical form information corresponding thereto, to obtain the interpolated key conditions and the chemical form information corresponding thereto;

[0036] The extraction unit is used to select the pure substance type and content of the preset selected element from the chemical form information under each key condition after interpolation processing.

[0037] In a third aspect, the present application provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method for predicting the chemical form and chemical reaction of the multiphase system is implemented.

[0038] In a fourth aspect, the present application provides a computer-readable storage medium having computer instructions stored thereon, which, when executed by a processor, implement the method for predicting the chemical form and chemical reaction of the multiphase system.

[0039] As can be seen from the above technical solution, the present application provides a method and device for predicting the chemical form and chemical reaction of a multiphase system. The method includes: obtaining key conditions corresponding to the multiphase system, the key conditions including temperature, pressure, multiple element types and their respective contents; predicting chemical form information under the key conditions by using the Gibbs free energy minimization method, the chemical form information including various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium; based on the chemical form information, predicting the chemical reaction that will occur in the multiphase system under the key conditions, which can improve the accuracy and efficiency of predicting the chemical reaction of the multiphase system on the basis of predicting the chemical form of the elements in the multiphase system; specifically, not only can the chemical form composition information be provided in the phase diagram, but also the relative content of each chemical form can be provided along the X-axis, so that the phase diagram information is richer, which is beneficial for subsequent analysis, and can clearly describe the composition of the chemical form under the multiphase system and clearly show the changing trend of the chemical form content of each element. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0041] FIG1 is a schematic diagram of a curve showing the change of CsI in a multiphase system with temperature in an example of the prior art;

[0042] FIG2 is a schematic diagram of a curve showing the change of Cs2I2 with temperature in a multiphase system according to an example of the prior art;

[0043] FIG3 is a logic diagram of a method for predicting chemical forms and chemical reactions of a multiphase system in an embodiment of the present application;

[0044] FIG4 is a schematic diagram of a first flow chart of a method for predicting chemical forms and chemical reactions of a multiphase system in an embodiment of the present application;

[0045] FIG5 is a schematic diagram of a second flow chart of the method for predicting chemical forms and chemical reactions of a multiphase system in an embodiment of the present application;

[0046] FIG6 is a schematic diagram of a third flow chart of the method for predicting chemical forms and chemical reactions of a multiphase system in an embodiment of the present application;

[0047] FIG7 is an OT phase diagram of iodine in an example of the present application;

[0048] FIG8 is a schematic diagram of a fourth flow chart of the method for predicting chemical forms and chemical reactions of a multiphase system in an embodiment of the present application;

[0049] FIG9 is a schematic diagram showing the distribution of Gibbs free energy of chemical reaction of iodine in a multiphase system according to an example of the present application;

[0050] FIG10 is a schematic diagram of the structure of a device for predicting chemical forms and chemical reactions of a multiphase system according to an embodiment of the present application;

[0051] FIG11 is a schematic block diagram of the system structure of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0052] In order to enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0053] In order to solve the problems existing in the above-mentioned prior art, the embodiment of the present application provides a chemical reaction detection method and device for a multiphase system. As shown in FIG3 , the content of pure substances that reach thermodynamic equilibrium can be calculated by minimizing the Gibbs free energy according to the element mass and chemical environment of the multiphase system; the chemical form of the element under certain conditions is obtained by analyzing the data; finally, based on the chemical form obtained by analysis, the chemical reactions that may exist in the process are searched and split into simple chemical reactions; combined with the Gibbs free energy change less than 0 and the thermodynamic equilibrium condition, the most likely chemical reaction in the multiphase system is finally determined. Specifically, a method can be provided for batch processing of a multiphase system based on a thermodynamic calculation method or software and for studying and analyzing the chemical forms and chemical reactions obtained by the processing. The chemical form of the multiphase system can be analyzed and the chemical reactions that may occur in the multiphase system can be summarized and analyzed. The thermodynamic equilibrium calculation can be completed in batches by calling the equilibrium calculation module of FactSage or other thermodynamic software; a two-dimensional phase diagram of the system can be drawn, and the error of the phase diagram can be discussed when other parameters are fixed; complex reactions and simple reactions that may occur in the multiphase system can be determined; the influence of the chemical form and chemical reaction of the element after the key influencing parameters are changed can be discussed. This approach uses thermodynamic calculations to determine the chemical speciation of elements in multiphase systems. Based on this calculation and combined with the Gibbs free energy criterion, it can predict the chemical reactions that are likely to occur in these systems. This research deepens our understanding of chemical processes in multiphase systems and provides an important foundation for studying the reaction mechanisms and chemical behavior of elements in these systems. This approach can be widely applied in industrial settings such as chemical engineering, energy, and the environment. For example, in nuclear energy systems, when studying element-material interactions (adsorption, desorption, deposition, etc.), the adsorption and deposition effects on pipelines, as well as their transport behavior, can vary significantly when the elements exist in different chemical speciations. Determining the chemical speciation of elements through this approach allows for targeted research on key issues, reduces trial-and-error costs, and more accurately describes the physicochemical behavior of elements. Furthermore, in nuclear energy systems, effluent monitoring and control are essential. Determining the chemical speciation present in the system allows for the design of targeted measurement methods, monitoring systems, and control measures, improving the accuracy and efficiency of nuclear power plant monitoring.

[0054] The details are described in detail through the following embodiments.

[0055] In order to improve the accuracy and efficiency of predicting chemical reactions in a multiphase system based on predicting the chemical forms of elements in the multiphase system, this embodiment provides a method for predicting the chemical forms and chemical reactions in a multiphase system, wherein the execution subject is a device for predicting the chemical forms and chemical reactions in the multiphase system. The device for predicting the chemical forms and chemical reactions in the multiphase system includes, but is not limited to, a server. As shown in FIG4 , the method specifically includes the following contents:

[0056] Step 100: Acquire key conditions corresponding to the multiphase system, wherein the key conditions include temperature, pressure, multiple element types and their respective contents.

[0057] Specifically, the key conditions can be derived from experimental measurements and theoretical calculations. In one example, the temperature is 2.00E+02 (°C) and the pressure is 5.00E-03 (GPa). The various element types and their respective contents corresponding to the multiphase system are as follows:

[0058] (6.18E-08)Kr+(3.82E-07)Xe+(2.58E-10)I+(6.40E-13)Sr+(5.96E-10)Cs+(6.73E-12)Ag+(1.75E-11)Co+(3.88E-13)Fe+(3.73E-15)Cr+(1.81E-15)Mn+(3.33E-13)Ni+(7.93E-12)Rb+(750000)He+(30.75)C+(4.99)H+(6.11)O+(0.20)N. The unit of content is mole.

[0059] Step 200: Predicting chemical form information under the critical conditions using a Gibbs free energy minimization method. The chemical form information includes various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium.

[0060] Specifically, the chemical form of the element in the multiphase system may represent the composition of simple substances and compounds containing the element in the multiphase system.

[0061] For example, the chemical form information under the key conditions includes pure substance types such as H, H₂, He, C, and C₂, with the contents of H, H₂, He, C, and C₂ being (2.88E-20) mol, (4.52E-01) mol, (7.50E+05) mol, (1.17E-67) mol, and (0.00E+00) mol, respectively. Each pure substance can be composed of the aforementioned elements. The pure substance types can include simple substances and compounds.

[0062] Step 300: Based on the chemical form information, predict the chemical reaction that will occur in the multiphase system under the key conditions.

[0063] Specifically, based on the chemical form information, a chemical equation balancing method can be used to predict the chemical reactions that will occur in the multiphase system under the critical conditions. A target pure substance containing a preset selected element can be selected from the various pure substances. Based on the type and content of the target pure substance, a chemical equation balancing method can be used to predict the chemical reactions that will occur in the multiphase system under the critical conditions and that are related to the preset selected element. Furthermore, a simple chemical equation can be screened from the predicted chemical reactions, and the Gibbs free energy values ​​of each pure substance in the simple chemical equation can be obtained and substituted into the simple chemical equation to obtain the Gibbs free energy change of the simple chemical reaction. If the Gibbs free energy change of the simple chemical reaction is less than zero, it is determined that the simple chemical reaction can occur spontaneously. The preset selected element can be at least one of the multiple elements, such as iodine I.

[0064] For example, first, by setting element conditions, all elements related to I are searched and screened. For example, the element conditions are set to I, O, Rb, and H. After searching and screening, pure substances related to I are obtained, including HIO, IRb, etc. The chemical equation balancing algorithm can be used to screen out the chemical reactions including: 1I(g)+1Rb(g)→1IRb(g), 2HIO(g)+1H2(g)→2I(g)+2H2O(g), etc.

[0065] Specifically, by applying the chemical equation balancing algorithm, the chemical equation of the pure substance related to element I after balancing can be obtained; then, according to the chemical equation combination and decomposition algorithm, a simple chemical equation in the system is obtained from the chemical equation of the pure substance related to element I, that is, the sum of the number of reactants and products is less than or equal to four, and most of the time there is only one or two reactants. For the simple chemical equation, the Gibbs free energy value in g0table is substituted into the simple chemical equation (i.e., the screened equation), and the Gibbs free energy change of this simple chemical equation is calculated. If the Gibbs free energy change of the chemical reaction is less than zero, it is determined that the chemical reaction can occur spontaneously. The specific process of applying the chemical equation balancing, combination and decomposition algorithm is as follows:

[0066] The chemical equation can be written as:

[0067] in:

[0068] R j 、P jare components of the reactant and product j respectively; the reactants and products may be pure substances corresponding to the pure substance types in the above chemical form information.

[0069] n rj 、n pj are the stoichiometric coefficients of reactant and product j respectively;

[0070] In a chemical equation, the systematic components of the elements on both sides of the equation must be the same: AN r =BN p

[0071] in:

[0072] In the above formula, m represents the number of system components (elements) of the multiphase system, n r is the number of reactants, n p is the number of generated objects, a ij 、b ij is the ratio of system component j to component i in the reactants and products, N r and N p is a vector expression of the stoichiometric quantities of reactants and products.

[0073] The above formula can be normalized to: CN e =0

[0074] Where: C=(A,-B)N e =(N r ,N p ) T

[0075] n e is the total number of reactants and products.

[0076] r(C) is the rank of the matrix C, when:

[0077] r(C)=min(m,n e ): has a unique zero solution, and the solution is meaningless;

[0078] r(C)=min(m,n e )-1: has infinite non-zero solutions and the only simplest integer ratio solution;

[0079] r(C) <min(m,n e )-1: has infinite non-zero solutions and infinite simplest integer ratio solutions.

[0080] The above is the algorithm for balancing chemical equations.

[0081] The combination and decomposition of chemical reactions is crucial to describing all chemical reactions of nuclides in complex systems. Consider an arbitrary set of chemical equations in a complex system such as:

[0082] Among them, n i,rj and n i,pj is the stoichiometric number of reactants and products j in equation i. The conservation of mass for the above chemical equations can be described as:

[0083] Among them, n s =n r +n p is the total number of components in the complex system. For a given linear space of all chemical equations and simple chemical equations In the linear space, if we select the first chemical equation in the above chemical equation group, that is, select equation 0 as any chemical equation in the entire multiphase system, and select equations 1 to k as all simple chemical equations. The matrix expression is:

[0084] The solution of the above equations is determined by the coefficient matrix N of the linear equation system s The rank r(N s ), as shown below:

[0085] r(N s )=min(k+1,n s ): has a unique zero solution, and the solution is meaningless;

[0086] r(N s )=min(k+1,n s )-1: has infinite non-zero solutions and the only simplest integer ratio solution;

[0087] r(N s ) <min(k+1,n s )-1: has infinite non-zero solutions and infinite simplest integer ratio solutions.

[0088] Assume that equation 0 is any chemical equation in the entire complex system, and chemical equations 1 to k are selected as all basic chemical equations. Then W=Ω means that chemical equation 0 is composed of equations 1 to k, which means that the following equation has a non-zero solution.

[0089] in, make W is composed of β1,β2,…,β k In a certain linear space; let β i is the stoichiometric coefficient of chemical equation i. The following four conditions are usually met:

[0090] Represents the balance of a chemical equation;

[0091] k≥n s -1, which means the number of chemical equations is greater than the number of components in the system;

[0092] This means that each component is present in at least one chemical equation;

[0093] β i There are at most four non-zero coordinates, which means that chemical equation i has at most three products or three reactants.

[0094] Therefore, d(W)=d(Ω)=n s -1 and That is, W=Ω, if k′≥n s -1, it means that the entire reaction system involves a chemical reaction of 1 to 2 reactants and products.

[0095] Then, based on the combination and decomposition of the reactants, the system is screened for simple chemical equations, meaning those with a total of four or fewer reactants and products, most often with only one or two reactants. Other chemical equations are reconstructed using these simple chemical equations.

[0096] For the selected balanced simple chemical reaction equations:

[0097] Where R and P are reactants and products respectively, ν R and ν P are the corresponding stoichiometric coefficients. The molar Gibbs free energy change Δ of the chemical reaction r G m,T for:

[0098] Δ f G m,T (R) and Δ f G m,T (P) are the molar Gibbs free energies of formation of reactants and products, respectively, Δ r G m,T is the molar Gibbs free energy change of the chemical reaction. Substituting the data in g0table into the above formula, we can get the Gibbs free energy change of the screened chemical reactions and determine Δ r G m,T Is it greater than 0? If it is less than 0, it can occur spontaneously.

[0099] In order to further improve the reliability of obtaining key conditions, in one embodiment of the present application, step 100 may include: inputting a first input file and input conditions; applying the first input file and the input conditions to generate a second input file, and each record in the second input file may be equivalent to one of the above-mentioned key conditions.

[0100] In order to determine the chemical reactions occurring under multiple key influencing parameters, and thus facilitate the subsequent analysis of the influence of the key influencing parameters on the chemical form of the multiphase system, in one embodiment, there are multiple key conditions, and the key influencing parameters of the multiphase system are at least one of the impurity elements in the multiple elements, the temperature, and the pressure; the key influencing parameters corresponding to the respective key conditions have different values; accordingly, as shown in FIG5 , after step 100, the method further includes:

[0101] Step 400: Predict the chemical form information under each key condition using the Gibbs free energy minimization method.

[0102] Step 500: Based on the chemical form information under each key condition, predict the chemical reaction that will occur in the multiphase system under the key condition.

[0103] Specifically, each key condition includes: temperature, pressure, multiple elements, and their respective contents; the key influencing parameters corresponding to each key condition can be the same, and the values ​​of the corresponding key influencing parameters can be the same. For example, the key influencing parameters of the multiphase system are temperature and carbon element C. The key influencing parameters of key condition a are: temperature and carbon element, the temperature value is 2.00E+02°C, and the carbon content is 30.75 mol; the key influencing parameters of another key condition b are: temperature and carbon element, the temperature value is 2.00E+02°C, and the carbon content is 29.41 mol; the carbon content in key conditions a and b is different. The values ​​of the key influencing parameters corresponding to each key condition are different, and the other values ​​are the same.

[0104] In order to improve the reliability of the generated phase diagram and facilitate subsequent analysis, as shown in FIG6 , in one embodiment, after step 400, the following steps are further included:

[0105] Step 600: Select the pure substance type and content of the preset selected element from the chemical form information under each key condition.

[0106] Step 700: Generate a two-dimensional phase diagram corresponding to the preset selected element based on the type and content of the pure substance containing the preset selected element; wherein the two-dimensional phase diagram is divided into multiple regions, each region corresponds to a group of pure substance groups, each region corresponds to a different pure substance group, each pure substance group is obtained based on the type of pure substance containing the preset selected element, the phase boundary value between each region is a preset phase boundary value, and the two-dimensional phase diagram also includes: a change curve for each pure substance type containing the selected element.

[0107] Specifically, the preset phase boundary value can be set according to actual conditions, and this application does not impose any restrictions on this. Preferably, it can be 0.01.

[0108] In one example, a two-dimensional phase diagram of iodine is shown in FIG7 , where w1 to w5 represent the change curves corresponding to the pure substances I(g), HI(g), HIO(g), IRb(g), and CsI(g), respectively, with the abscissa representing temperature in degrees Celsius. The left side of the phase diagram represents the ordinate corresponding to the change curve, with the ordinate representing the relative content of each pure substance. v1 to v8 represent pure substance groups, respectively: I(g), HIO(g)+I(g), HI(g)+I(g), CsI(g)+HI(g)+IRb(g), CsI(g)+HI(g)+I(g), CsI(g)+HI(g)+I(g)+IRb(g), CsI(g)+HI(g)+HIO(g)+I(g), and CsI(g)+HI(g)+HIO(g)+I(g)+IRb(g). The right side of the phase diagram shows the corresponding region on the ordinate, which represents the oxygen content in moles (mol). In this example, the phase boundary between regions v1 and v2 is q. The change curve in the two-dimensional phase diagram can be equivalent to the upper phase diagram, and each region in the two-dimensional phase diagram can be equivalent to the lower phase diagram. In the lower phase diagram, stars can be used to represent the reference conditions of two system variables. The HTR-PM system has different operating conditions. Based on one of these operating conditions, the temperature, pressure, and impurity contents of C, H, O, and N are set. The purpose of setting reference conditions is that when considering the changing relationship between two system variables, such as the temperature and O content in Figure 7, the remaining variables (such as pressure, C, H, and N content) serve as reference conditions. The position represented by the star in Figure 7 represents the set reference condition for the system. The dashed line represents the reference condition for the system variable depicted on the Y-axis of the upper phase diagram. For example, the dashed line represents the temperature change of each pure substance calculated under the condition of an O content of 0.35 mol.

[0109] To ensure the integrity of chemical form information, as shown in FIG8 , in one embodiment, step 600 includes:

[0110] Step 601: performing interpolation processing based on each key condition and the chemical form information corresponding thereto to obtain interpolated key conditions and the chemical form information corresponding thereto.

[0111] Step 602: Select the pure substance type and content of the preset selected element from the chemical form information under each interpolated key condition.

[0112] Specifically, based on the interpolation algorithm, each key condition and the chemical form information corresponding thereto, the interpolated key conditions and the chemical form information corresponding thereto may be obtained.

[0113] The interpolation algorithm can be as follows:

[0114] N-dimensional matrix table interpolation function F(x1,...,x N ) is calculated based on the function table value f(x1,…,x N ), expressed as:

[0115] Function f(x1,...,x N ) is an approximate interpolation function F(x1,...,x N ) is:

[0116] First determine the appropriate table range: in interval Now we scale each interval to (0,1), which converts the problem of finding g(y1,...,y N ) matrix table interpolation function G(y1,...,y N ) problem, the interval Ω G ={(y1,…,y N ):0≤y1≤1,..,0≤y N ≤1}.

[0117] The N-dimensional matrix after conversion is the interpolation function G(y1,...,y N ) is determined by the following formula: [g(j1,…,j N ):j1=0 or 1,…,j N =0 or 1]

[0118] Among them, the piecewise linear interpolation F L Computed as a simplex of complexity N! p(1),…,p(N) =[(y1,…,y N ):0≤yp(1) ≤…≤y p(N) ≤1]

[0119] where (p(1), p(2), …, p(N)) is a permutation of integers in the range 1–N. The simplex is formed by the ordering coefficients y1, …, y N OK. L The form is: F L =l0+l1y1+…+l N y N

[0120] Among them, l i is defined as the simplex with vertices s0 and s N The coefficients of f matching at , and the vertex s i At position p(j), j>i is 1, and other positions are 0.

[0121] For Ω F Each z in the range = (x1,...,x N ) T , the error between the interpolation function and the original function is:

[0122] Among them, D θ f(x) is the directional derivative of f(x) in the θ direction.

[0123] To further illustrate this solution, the present application provides a method for detecting a chemical reaction in a multiphase system, which is described in detail as follows:

[0124] S1. Parameter configuration and data preprocessing.

[0125] First, before the calculation begins, set key parameters, including module control parameters (divide the calculation process into two modules, enabling the phase diagram drawing and reaction analysis modules respectively according to the calculation requirements). Second, generate the user input file (Template.equi file) and calculation parameters (input.txt) in the specific format required for FactSage to run. Finally, establish the corresponding file storage path, including input files and output files. The user input file can perform the same function as the first input file described above, and the calculation parameters can perform the same function as the input conditions described above.

[0126] The most critical input conditions in the user input file (i.e., Template.equi file) are the types and contents of the elements considered in the multiphase system. In one example, the input conditions in the user input file include:

[0127] 6.18E-08Kr+3.82E-07Xe+2.58E-10I+6.40E-13Sr+5.96E-10Cs+6.73E-12Ag+1.75E -11Co+3.88E-13Fe+3.73E-15Cr+1.81E-15Mn+3.33E-13Ni+7.93E-12Rb+750000He+<Variable_C> C+<Variable_H> H+<Variable_O> O+<Variable_N> N

[0128] in,<Variable_X> X represents the impurity element X and its content in the multiphase system. X can be C, H, O, or N, and its content is given in the calculation parameters.

[0129] The calculation parameter Y includes: temperature, pressure and various impurity elements. The key influencing parameter of the multiphase system can be at least one of temperature, pressure and various impurity elements. In one example, the correspondence between the calculation parameter Y, the lower limit a of the variation interval, the upper limit b of the variation interval, and the number c of uniform values ​​within the variation interval is shown in Table 1. In this example, the number of uniform values ​​within the temperature variation interval is 26, which can be represented by 26 temperature values ​​taken at fixed intervals between 200°C and 1000°C; the number of uniform values ​​within the pressure variation interval is 1, which can be represented by the pressure remaining unchanged at 0.005GPa; the number of uniform values ​​within the C impurity element variation interval is 24, which can be represented by 24 contents taken at fixed intervals between 0mol and 30.75mol; the number of uniform values ​​within the variation intervals of the H, O, and N impurity elements is all 1, which can be represented by the content remaining unchanged. The contents of the H, O, and N impurity elements are 4.99mol, 6.11mol, and 0.2mol, respectively. The key influencing parameters of the multiphase system are temperature T and impurity element C.

[0130] Table 1

[0131] S2. Call FactSage software to run.

[0132] The FactSage input file *.equi is generated based on the user input file Template.equi and the calculation parameters input.txt. The functions implemented by the FactSage input file can be equivalent to the functions implemented by the second input file mentioned above.

[0133] In the above example, the key influencing parameters of the multiphase system are temperature T and impurity element C. Temperature T takes 26 values ​​and impurity element C takes 24 values. Therefore, 24 FactSage input files can be generated, and the content of impurity element C corresponding to each FactSage input file is different; the content of impurity element C in the same FactSage input file is the same and contains 26 temperature values, that is, there can be 26 records in the same FactSage input file, each record has a different temperature value, and other data are the same. A FactSage input file contains: temperature: (2.00E+02)℃, (2.33E+02)℃, (2.67E+02)℃..., pressure: 5.00E-03GPa, 6.18E-08Kr+3.82E-07Xe+2.58E-10I+6.40E-13Sr+5.96E-10Cs+6.73E-12Ag+1.75E-11Co+3.88E-13Fe+3.73E-15Cr+1.81E-15Mn+3.33E-13Ni+7.93E-12Rb+750000He+1.3C+4.99H+6.11O+0.2N, content unit is mol.

[0134] Apply the FactSage input file and call the FactSage macro module to perform the Equilib equilibrium calculation to obtain the FactSage calculation result .tab file. The FactSage input file corresponds to the .tab file one-to-one. Taking a .tab file as an example, considering the different types of elements in the multiphase system, the pure substances present in the multiphase system will also be different. After calculation, there are 406 pure substances of the elements in this multiphase system, that is, 406 compounds and single substances in total, which will all be output to the .tab file, as shown in Table 2. Here, only one calculation result, that is, a part of a .tab file, is listed for demonstration. Among them, each row of records in Table 2 can be equivalent to the temperature, pressure and content of each pure substance under the above-mentioned key conditions.

[0135] Table 2

[0136] S3. Extract key data from the calculation results to obtain reactant conditions, data for the chemical forms of all elements in the multiphase system under different conditions (datas), and Gibbs free energy data (g0table). For the datas dataset, separate the variables into reactant content and condition data (Xtable), and the content data for all pure substances in the multiphase system under the corresponding conditions (Ytable). This step is crucial for ensuring the integrity of the data obtained by running FactSage. Due to the large computational complexity, FactSage may miss small amounts of data, requiring manual interpolation to complete the data.

[0137] In one example, the reactants can be extracted according to Table 1 as shown in Table 3, with the content unit being mol:

[0138] Table 3

[0139] In this example, the Gibbs free energy values ​​of 406 pure substances at different temperatures were extracted to obtain the Gibbs free energy data g0table, which is listed in Table 4. In Table 4, G is the abbreviation of Gibbs free energy, J is its unit joule, g is gas, and H, H2, He, C, and C2 are substances:

[0140] Table 4

[0141] Xtable summarizes the content of each calculation parameter under all key conditions. Xtable is a 624×6 data table, as shown in Table 5. Only some of them are listed for illustration:

[0142] Table 5

[0143] Ytable summarizes the types and contents of pure substances in a multiphase system under all key conditions. Ytable is a 624×406 data table that can be composed of various .tab files.

[0144] S4. Phase diagram drawing.

[0145] Because complex multiphase systems often involve a large number of pure substances, a threshold for the relative percentage of pure substances must be provided when drawing phase diagrams. This allows for screening based on the concentration of the pure substance to determine its presence in the complex multiphase system. If the relative percentage of a pure substance exceeds the threshold, the pure substance is considered present.

[0146] Set the drawing parameters condition: [250, 0.005, 30.75, 4.99, 6.11, 0.20], which represent temperature, pressure, C, H, O, and N content, respectively.

[0147] Select the X and Y axes: the X axis is temperature T, and the Y axis is O content. After the X and Y axes are selected, the calculated parameters pressure, C, H, and N content are [0.005, 30.75, 4.99, 0.20] in the above drawing parameters.

[0148] Selected element: iodine element I, that is, the OT phase diagram showing element I in the phase diagram.

[0149] Interpolation density: [300, 300].

[0150] Before interpolation, the drawing data is the Xtable and Ytable mentioned above. Based on the interpolation algorithm mentioned above, since the interpolation density selected at this time is [300, 300], and the X-axis and Y-axis are selected, they become cons_Xtable and cons_Ytable after interpolation. cons_Xtable is a 90000×2 data table, and cons_Ytable is a 90000×406 data table. After interpolation, the amount of data increases, which provides convenience for phase diagram drawing.

[0151] A two-dimensional phase diagram has two layers. The lower layer, similar to a traditional phase diagram, provides the chemical species composition of elements in a multiphase system. Two variables (X and Y axes) are varied while other variables in the phase diagram remain constant. Different regions represent different chemical species. The boundary between two chemical species represents the point at which the relative percentage of one component in the adjacent region equals a set threshold, i.e., 0.01. This differs from the zero phase line in traditional phase diagrams. While the multiphase system variable described by the Y axis in the lower layer remains constant, the upper layer depicts the key chemical species of the elements in the multiphase system as a function of the multiphase system variable described by the X axis, i.e., the relative content of each chemical species varies along the X axis. In the color map, different chemical species are represented by different colored lines in the upper layer, and by different colored blocks in the lower layer. In addition, in the lower layer, stars represent the reference conditions for the two multiphase system variables, while dashed lines represent the reference conditions for the multiphase system variable described by the Y axis in the upper layer.

[0152] S5. Chemical reaction analysis.

[0153] According to the preset selected elements, relevant pure substances and chemical reactions are searched, and the possible basic chemical equations are returned. Combined with the g0table obtained above, the potential basic chemical reactions in the considered chemical system are obtained by calculating G0.

[0154] Specifically, after determining the chemical form, a search will be conducted for simple chemical equations, which are defined as chemical equations with one or two products and reactants. Since there are infinite chemical equations for the key conditions in a multiphase system, simple chemical equations will be listed, and other chemical equations in the multiphase system will be constructed based on linear combinations of simple chemical equations. Using the Gibbs free energy of chemical equations, simple chemical equations that can spontaneously occur are screened out from all basic chemical equations in the complex system, as shown in Figure 9. The ordinate represents the chemical reactions that may occur in element I in the multiphase system obtained through calculation, the abscissa represents temperature, and the color blocks represent changes in the Gibbs free energy of the chemical reaction.

[0155] The criteria for possible chemical reactions in complex multiphase systems are as follows:

[0156] a) Changes in the content of a substance within a certain range; b) Changes in the chemical form of the reactants; c) Chemical reactions occur in the direction of ΔG<0.

[0157] Balancing a chemical equation is a method for calculating chemical equations. A balanced chemical equation must obey the law of conservation of mass and accurately represent the mass ratios between reactants and products, providing precise relationships for chemical calculations.

[0158] As can be seen from the above description, this application example provides a method for batch calculations of multiphase systems using the thermodynamic calculation software FactSage and for analyzing the resulting chemical forms and reactions. This method can analyze the chemical forms of multiphase systems and summarize and analyze possible chemical reactions within the system. This method utilizes the equilibrium calculation module of FactSage to perform batch thermodynamic equilibrium calculations; plot a two-dimensional phase diagram for the system; and analyze the errors in the phase diagram when other parameters are fixed. It also identifies possible complex and simple reactions within the chemical system; and discusses the effects of varying key influencing parameters on the chemical forms and reactions of elements. This method can analyze the chemical forms of complex multiphase systems and predict possible chemical reactions within the system. Compared to traditional two-dimensional phase diagrams that only provide information on the chemical phase composition, this method not only provides information on the chemical form composition in the phase diagram, but also provides information on the relative content of each chemical form as it varies along the x-axis. This enriches the phase diagram information and facilitates subsequent analysis. Furthermore, and importantly, this method can also analyze possible chemical reactions within the system. Traditional approaches to phase diagram construction use the zero phase line as the demarcation line, which results in phase diagrams being unable to depict multi-element systems and providing limited information, including content information. However, this method, by setting a threshold, allows the phase demarcation line to be controlled to the 0.01 line or another user-defined demarcation line. This allows phase diagrams to clearly depict the chemical species composition of multiphase systems, providing significantly more information than traditional phase diagrams. Furthermore, by overlaying top-level phase diagrams, this method clearly demonstrates the changing trends in the chemical species content of each element. For chemical equation search, this method, based on a database of pure substances with defined chemical species, first identifies all pure substances with high concentrations associated with the element of interest, along with their thermodynamic data. Using the aforementioned chemical equilibrium algorithm, it then lists possible chemical reactions within the multiphase system. It then calculates the Gibbs free energy changes of the reactions to identify reactions that are likely to occur spontaneously. This method can be used to determine the chemical species and possible chemical reactions in multi-component and multiphase systems with varying elemental concentrations, reflecting physical principles and describing detailed information about the process.

[0159] From a software perspective, in order to improve the accuracy and efficiency of predicting chemical reactions in a multiphase system based on predicting the chemical forms of each element in the multiphase system, the present application provides an embodiment of a device for predicting the chemical forms and chemical reactions of the multiphase system. Referring to FIG10 , the device for predicting the chemical forms and chemical reactions of the multiphase system specifically includes the following contents:

[0160] Acquisition module 01 is used to obtain key conditions corresponding to the multiphase system, including temperature, pressure, multiple element types and their respective contents;

[0161] The first calculation module 02 is configured to predict chemical form information under the key conditions using a Gibbs free energy minimization method, wherein the chemical form information includes various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium;

[0162] The first prediction module 03 is configured to predict, based on the chemical form information, the chemical reaction that will occur in the multiphase system under the key conditions.

[0163] In one embodiment, there are multiple key conditions, and the key influencing parameter of the multiphase system is at least one of the impurity element type among the multiple element types, the temperature, and the pressure; and the value of the key influencing parameter corresponding to each key condition is different;

[0164] Correspondingly, the device for predicting the chemical form and chemical reaction of a multiphase system further includes:

[0165] A second calculation module is used to predict the chemical form information under each key condition by using a Gibbs free energy minimization method;

[0166] The second prediction module is used to predict the chemical reaction that will occur in the multiphase system under each key condition based on the chemical form information under the key condition.

[0167] In one embodiment, the device for predicting chemical forms and chemical reactions of a multiphase system further comprises:

[0168] An extraction module, configured to select the type of pure substance containing a preset selected element and its content from the chemical form information under each key condition;

[0169] a phase diagram drawing module, configured to generate a two-dimensional phase diagram corresponding to the preset selected element according to the type and content of the pure substance containing the preset selected element;

[0170] The two-dimensional phase diagram is divided into multiple regions, each region corresponds to a set of pure substance groups, and each region corresponds to a different pure substance group. Each pure substance group is obtained based on the pure substance type containing the preset selected element. The value of the phase boundary between each region is the preset phase boundary value. The two-dimensional phase diagram also includes: a change curve of each pure substance type containing the selected element.

[0171] In one embodiment, the extraction module includes:

[0172] an interpolation processing unit, configured to perform interpolation processing based on each key condition and the chemical form information corresponding thereto, to obtain the interpolated key conditions and the chemical form information corresponding thereto;

[0173] The extraction unit is used to select the pure substance type and content of the preset selected element from the chemical form information under each key condition after interpolation processing.

[0174] The embodiment of the device for predicting the chemical form and chemical reaction of a multiphase system provided in this specification can be specifically used to execute the processing flow of the embodiment of the method for predicting the chemical form and chemical reaction of a multiphase system described above. Its functions will not be described in detail here, and reference can be made to the detailed description of the embodiment of the method for predicting the chemical form and chemical reaction of a multiphase system described above.

[0175] FIG11 is a schematic diagram of the physical structure of a computer device provided in an embodiment of the present invention. As shown in FIG11 , the electronic device includes: a memory 501, a processor 502, and a computer program stored in the memory 501 and executable on the processor 502. When the processor 502 executes the computer program, the following method is implemented:

[0176] Step 100: Acquire key conditions corresponding to the multiphase system, wherein the key conditions include temperature, pressure, multiple element types and their respective contents.

[0177] Step 200: Predicting chemical form information under the critical conditions using a Gibbs free energy minimization method. The chemical form information includes various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium.

[0178] Step 300: Based on the chemical form information, predict the chemical reaction that will occur in the multiphase system under the key conditions.

[0179] This embodiment provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program. When the computer program is executed by a processor, the following method is implemented:

[0180] Step 100: Acquire key conditions corresponding to the multiphase system, wherein the key conditions include temperature, pressure, multiple element types and their respective contents.

[0181] Step 200: Predicting chemical form information under the critical conditions using a Gibbs free energy minimization method. The chemical form information includes various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium.

[0182] Step 300: Based on the chemical form information, predict the chemical reaction that will occur in the multiphase system under the key conditions.

[0183] It will be understood by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Thus, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

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

[0185] These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a product including an instruction device that implements the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

[0186] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

[0187] Throughout this specification, reference to terms such as "one embodiment," "a specific embodiment," "some embodiments," "for example," "example," "specific example," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0188] The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for predicting chemical forms and chemical reactions of a multiphase system, characterized in that: include: Obtaining key conditions corresponding to a multiphase system, including temperature, pressure, multiple element types and their respective contents; Predicting chemical form information under the critical conditions using a Gibbs free energy minimization method, the chemical form information including: various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium; Based on the chemical form information, the chemical reactions that will occur in the multiphase system under the key conditions are predicted.

2. The method for predicting chemical forms and chemical reactions of a multiphase system according to claim 1, wherein: There are multiple key conditions, and the key influencing parameter of the multiphase system is at least one of the impurity element type among the multiple element types, the temperature, and the pressure; and the key influencing parameter corresponding to each key condition has a different value; Correspondingly, after obtaining the key conditions corresponding to the multiphase system, the following is also included: The chemical form information under each key condition is predicted by the Gibbs free energy minimization method; Based on the chemical form information under each key condition, the chemical reaction that will occur in the multiphase system under the key condition is predicted.

3. The method for predicting chemical forms and chemical reactions of a multiphase system according to claim 2, wherein: After the chemical form information under each key condition is predicted by the Gibbs free energy minimization method, the method further includes: Selecting the type of pure substance and its content containing the preset selected element from the chemical form information under each key condition; generating a two-dimensional phase diagram corresponding to the preset selected element according to the type and content of the pure substance containing the preset selected element; The two-dimensional phase diagram is divided into multiple regions, each region corresponds to a set of pure substance groups, and each region corresponds to a different pure substance group. Each pure substance group is obtained based on the pure substance type containing the preset selected element. The value of the phase boundary between each region is the preset phase boundary value. The two-dimensional phase diagram also includes: a change curve of each pure substance type containing the selected element.

4. The method for predicting chemical forms and chemical reactions of a multiphase system according to claim 3, wherein: The step of selecting the pure substance type and content of the preset selected element from the chemical form information under each key condition includes: Performing interpolation processing based on each key condition and the chemical form information corresponding thereto to obtain interpolated key conditions and the chemical form information corresponding thereto; From the chemical form information under each interpolated key condition, the pure substance type and content of the preset selected element are selected.

5. A device for predicting chemical forms and chemical reactions of a multiphase system, characterized in that: include: An acquisition module is used to obtain key conditions corresponding to the multiphase system, wherein the key conditions include temperature, pressure, multiple element types and their respective contents; A first calculation module is configured to predict chemical form information under the key conditions using a Gibbs free energy minimization method, the chemical form information including various pure substance types and their respective contents when the multiphase system reaches thermodynamic equilibrium; The first prediction module is used to predict the chemical reaction that will occur in the multiphase system under the key conditions based on the chemical form information.

6. The device for predicting chemical forms and chemical reactions of a multiphase system according to claim 5, characterized in that: There are multiple key conditions, and the key influencing parameter of the multiphase system is at least one of the impurity element type among the multiple element types, the temperature, and the pressure; and the key influencing parameter corresponding to each key condition has a different value; Correspondingly, the device for predicting the chemical form and chemical reaction of a multiphase system further includes: A second calculation module is used to predict the chemical form information under each key condition by using a Gibbs free energy minimization method; The second prediction module is used to predict the chemical reaction that will occur in the multiphase system under each key condition based on the chemical form information under the key condition.

7. The device for predicting chemical forms and chemical reactions of a multiphase system according to claim 6, characterized in that: Also includes: An extraction module, configured to select the type of pure substance containing a preset selected element and its content from the chemical form information under each key condition; a phase diagram drawing module, configured to generate a two-dimensional phase diagram corresponding to the preset selected element according to the type and content of the pure substance containing the preset selected element; The two-dimensional phase diagram is divided into multiple regions, each region corresponds to a set of pure substance groups, and each region corresponds to a different pure substance group. Each pure substance group is obtained based on the pure substance type containing the preset selected element. The value of the phase boundary between each region is the preset phase boundary value. The two-dimensional phase diagram also includes: a change curve of each pure substance type containing the selected element.

8. The device for predicting chemical forms and chemical reactions of a multiphase system according to claim 7, wherein: The extraction module comprises: an interpolation processing unit, configured to perform interpolation processing based on each key condition and the chemical form information corresponding thereto, to obtain the interpolated key conditions and the chemical form information corresponding thereto; The extraction unit is used to select the pure substance type and content of the preset selected element from the chemical form information under each key condition after interpolation processing.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the method for predicting the chemical form and chemical reaction of a multiphase system according to any one of claims 1 to 4 is implemented.

10. A computer-readable storage medium having computer instructions stored thereon, characterized in that: When the instructions are executed by a processor, the method for predicting the chemical form and chemical reaction of a multiphase system according to any one of claims 1 to 4 is implemented.

Images