Method, device and storage medium for calculating fouling risk of reactor

Abstract

The present application provides a method, device, and storage medium for calculating the fouling risk of a reactor, which relates to the technical field of nuclear reactors. The method for calculating the fouling risk of a reactor provided by the present application includes: step S110, obtaining input parameters; step S120, calculating the mass evaporation rate based on the input parameters; step S130, calculating the fouling thickness of the reactor core based on the mass evaporation rate; step S140, calculating the boron deposition mass based on the fouling thickness of the reactor core; step S150, determining whether the convergence condition is reached. If the convergence condition is reached, step S160 is entered; otherwise, the process returns to step S110; step S160, outputting the calculation result. The present application calculates the fouling thickness of the reactor core and the boron deposition mass based on the input parameters through the risk formation mechanism and important phenomena caused by reactor fouling, so as to evaluate the risks brought by fouling, including CILC risk and CIPS risk.

Description

Technical Field

[0001] This application relates to the technical field of nuclear reactors, and in particular, to a method, device, and storage medium for calculating the fouling risk of a reactor. Background Art

[0002] Axial power shift (CIPS) and local cladding corrosion (CILC) caused by fouling in the reactor core will affect the safety and economy of reactor operation, and are also one of the key challenges in improving the performance of pressurized water reactors, which has attracted more and more attention in the nuclear power industry. The phenomena caused by fouling in the reactor core are important influencing factors for the safety and economy of reactor operation, and the evaluation of the fouling risk in the reactor core may become an important issue that the industry will focus on in the future. Summary of the Invention

[0003] In view of this, this application provides a method, device, and storage medium for calculating the fouling risk of a reactor to realize the fouling risk analysis of the reactor core.

[0004] In a first aspect, this application provides a method for calculating the fouling risk of a reactor, including:

[0005] Step S110: Obtain input parameters;

[0006] Step S120: Calculate the mass evaporation rate based on the input parameters;

[0007] Step S130: Calculate the fouling thickness of the reactor core based on the mass evaporation rate;

[0008] Step S140: Calculate the boron deposition mass based on the fouling thickness of the reactor core;

[0009] Step S150: Determine whether the convergence condition is reached. If the convergence condition is reached, go to step S160; otherwise, return to step S110;

[0010] Step S160: Output the calculation result.

[0011] In a second aspect, this application provides a computing device, including:

[0012] At least one processor; and

[0013] At least one memory storing instructions, which when executed by the at least one processor alone or jointly, cause the computing device to execute the method as described in the first aspect.

[0014] In a third aspect, this application provides a computer storage medium storing instructions, which when executed by at least one processor of a computing device alone or jointly, cause the computing device to execute the method as described in the first aspect.

[0015] Compared with the prior art, the present application has the following advantages:

[0016] The method for calculating the fouling risk of the reactor provided by the present application includes: step S110, obtaining input parameters; step S120, calculating the mass evaporation rate based on the input parameters; step S130, calculating the fouling thickness of the reactor core based on the mass evaporation rate; step S140, calculating the boron deposition mass based on the fouling thickness of the reactor core; step S150, determining whether the convergence condition is reached. If the convergence condition is reached, go to step S160; otherwise, return to step S110; step S160, outputting the calculation result. The present application calculates the fouling thickness of the reactor core and the boron deposition mass based on the input parameters through the risk formation mechanism and important phenomena caused by reactor fouling, so as to evaluate the risks brought by fouling, including CILC risk and CIPS risk. Description of the Drawings

[0017] The accompanying drawings are provided to provide a further understanding of the present application. They are incorporated and constitute a part of the present application. The accompanying drawings illustrate embodiments of the present application and, together with this specification, serve to explain the principles of the present application. In the accompanying drawings:

[0018] Figure 1 is a schematic flowchart of a method for calculating the fouling risk of a reactor provided by an embodiment of the present application;

[0019] Figure 2 is a schematic flowchart of a method for calculating the mass evaporation rate provided by an embodiment of the present application;

[0020] Figure 3 is a schematic flowchart of a method for calculating the fouling thickness of a reactor core provided by an embodiment of the present application;

[0021] Figure 4 is a schematic flowchart of a method for calculating the boron deposition mass provided by an embodiment of the present application;

[0022] Figure 5 is a schematic diagram of fouling grid division provided by an embodiment of the present application;

[0023] Figure 6 is a schematic structural diagram of a computing device provided by an embodiment of the present application. Detailed Embodiments

[0024] To more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the accompanying drawings required for the description of the embodiments. Obviously, the accompanying drawings in the following description are only some examples or embodiments of the present application. For those of ordinary skill in the art, without creative efforts, the present application can also be applied to other similar scenarios based on these drawings. Unless obvious from the language context or otherwise stated, the same reference numerals in the figures represent the same structure or operation.

[0025] As shown in the present application, unless the context clearly indicates an exceptional situation, words such as "a", "an", "one", and / or "the" are not specifically singular and may also include plural. Generally speaking, the terms "comprising" and "including" only indicate the inclusion of the clearly identified steps and elements, and these steps and elements do not constitute an exclusive list. The method or device may also include other steps or elements.

[0026] At the same time, the present application uses specific terms to describe the embodiments of the present application. For example, "one embodiment", "an embodiment", and / or "some embodiments" mean a certain feature, structure, or characteristic related to at least one embodiment of the present application. Therefore, it should be emphasized and noted that the "one embodiment" or "an embodiment" or "an alternative embodiment" mentioned twice or more at different positions in this specification is not necessarily the same embodiment. In addition, certain features, structures, or characteristics in one or more embodiments of the present application can be appropriately combined.

[0027] Unless otherwise specifically stated, the relative arrangements of components and steps, numerical expressions, and values described in these embodiments do not limit the scope of the present application. At the same time, it should be understood that for the convenience of description, the dimensions of each part shown in the drawings are not drawn according to the actual proportional relationship. Technologies, methods, and devices known to those of ordinary skill in the relevant fields may not be discussed in detail, but in appropriate cases, the said technologies, methods, and devices should be regarded as part of the specification. In all the examples shown and discussed here, any specific value should be interpreted as merely exemplary and not as a limitation. Therefore, other examples of the exemplary embodiments may have different values. It should be noted that similar reference numerals and letters represent similar items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further discussed in subsequent drawings.

[0028] In addition, although the terms used in the present application are selected from well-known and commonly used terms, some of the terms mentioned in the specification of the present application may be selected by the applicant according to his or her judgment, and their detailed meanings are described in the relevant parts of this description. In addition, it is required to understand the present application not only through the actual terms used, but also through the meaning implied by each term.

[0029] In this application, flowcharts are used to illustrate the operations performed by devices or equipment according to the embodiments of this application. It should be understood that the operations before or below do not necessarily need to be executed precisely in sequence. Instead, they can be executed in reverse order or simultaneously, and steps can be processed in parallel. At the same time, other operations can be added to these processes, or one or more steps can be removed from these processes.

[0030] During the operation of the reactor, corrosion products (such as nickel, iron, etc.) on the coolant system pipes will be released into the coolant and deposited on the surface of the fuel assembly cladding to form dirt. Boron and lithium hydroxide accumulate in the pore-like gaps of the dirt. When the boron accumulates to a certain extent, it will precipitate onto the dirt, resulting in a decrease in the power of the upper part of the fuel assembly, and thus the power peak will shift downward, that is, the CIPS phenomenon occurs. If the dirt deposition is thick at a local position, the coolant will not be able to flow through the dirt to cool the surface of the fuel cladding, resulting in a relatively high temperature of the fuel cladding. If the fuel cladding is at a very high temperature for a long time, it will accelerate the corrosion and failure of the fuel cladding, thus leading to the CILC phenomenon. When the reactor power is increased, the average temperature of the reactor coolant is increased, or the core power peak factor is increased, the CIPS / CILC phenomenon is very likely to occur, and it usually occurs in the high-power fuel assemblies in the core.

[0031] In the embodiments of this application, during the process of calculating the dirt risk of the reactor, the calculation object (such as the core, steam generator, etc.) has a large space. For the accuracy of the calculation results, the calculation object can be divided into multiple nodes, and calculations are performed for each node separately. Then, in the embodiments of this application, performing calculations on the calculation object means performing calculations on multiple nodes of the divided calculation object.

[0032] Figure 1 is a schematic flow diagram of a method for calculating the dirt risk of a reactor provided by an embodiment of this application. As Figure 1 shown, the method for calculating the dirt risk of the reactor includes the following steps:

[0033] Step S110: Obtain input parameters.

[0034] The input parameters include some parameters required for subsequent calculations. In actual operations, the input parameters can be provided by the user through an input text file, the calculation results of other modules (such as computational fluid dynamics software, etc.), restart file information, etc.

[0035] In some embodiments, the values of some parameters to be solved can be iteratively updated during the iterative calculation process, such as the dirt thickness, etc. The initial value needs to be set at each time step. For example, the initial value of the input parameters at the current time step can be set according to the calculation results of the previous time step.

[0036] Step S120: Calculate the mass evaporation rate based on the input parameters.

[0037] Embodiments of the present application calculate the mass evaporation rate based on input parameters, which can improve the calculation accuracy and thus improve the accuracy of fouling risk analysis. In some embodiments, please refer to Figure 2 , calculating the mass evaporation rate based on the input parameters includes the following steps:

[0038] Step S201, construct the momentum equation.

[0039] Construct the momentum equation based on the heat flux density transferred by heat conduction, the mass evaporation rate, the vapor-phase enthalpy value of the fouling surface, and the liquid-phase enthalpy value of the fouling surface. An example of the relationship formula of the momentum equation is as follows:

[0040] (1)

[0041] In the formula: Q is the heat flux density transferred by heat conduction, with the unit of W / m 2 ;

[0042] x is the coordinate axis perpendicular to the fouling surface;

[0043] G evap is the mass evaporation rate, with the unit of kg / m 2 ·s;

[0044] H l is the liquid-phase enthalpy value of the fouling surface, with the unit of J / kg;

[0045] H v is the vapor-phase enthalpy value of the fouling surface, with the unit of J / kg.

[0046] Among them, Q, H l , H v are input parameters.

[0047] Step S202, construct the energy equation.

[0048] Construct the energy equation based on the vapor pressure inside the fouling, the liquid pressure inside the fouling, the liquid surface tension inside the fouling, the contact angle of the fouling, and the critical pore radius of the fouling. An example of the relationship formula of the energy equation is as follows:

[0049] (2)

[0050] In the formula: P v is the vapor pressure inside the fouling, with the unit of Pa;

[0051] P l is the liquid pressure inside the fouling, with the unit of Pa;

[0052] is the liquid surface tension inside the fouling, with the unit of N / m;

[0053] is the fouling contact angle, in °;

[0054] is the critical pore radius of fouling, in m.

[0055] Among them, and are input parameters.

[0056] Step S203: Construct a heat transfer equation.

[0057] Based on the internal temperature of fouling, the saturation temperature of the liquid inside the fouling, the mass evaporation rate, and the latent heat of vaporization of the liquid inside the fouling, construct a heat transfer equation. An example of the relationship of the heat transfer equation is as follows:

[0058] (3)

[0059] In the formula: h v is the volume heat transfer coefficient, in J / m 3 ·s·°C;

[0060] T crud is the internal temperature of fouling, in °C;

[0061] T sat is the saturation temperature of the liquid inside the fouling, in °C;

[0062] G evap is the mass evaporation rate, in kg / m 2 ·s;

[0063] h fg is the latent heat of vaporization of the liquid inside the fouling, in J / kg.

[0064] Among them, h v , T crud , T sat , h fg are input parameters.

[0065] Step S204: Calculate the mass evaporation rate based on the momentum equation, the energy equation, and the heat transfer equation.

[0066] In actual operation, step S204 may include the following steps:

[0067] 1) Set the boundary value on the fouling side;

[0068] 2) Set the upper limit value of the fouling layer, and set the upper limit value of the fouling layer not exceeding 500;

[0069] 3) Solve the momentum equation. When the solution of the relationship (1) converges, proceed to the next calculation;

[0070] 4) Solve the energy equation;

[0071] 5) Solve the heat transfer equation. When the solution of relation (3) converges, otherwise return to step 3);

[0072] 6) When the calculation reaches the target condition, such as when the calculation has reached the dry-out position, end the process; otherwise, return to step 2).

[0073] Step S130: Calculate the fouling thickness of the reactor core based on the mass evaporation rate.

[0074] Calculating the fouling thickness of the reactor core based on the mass evaporation rate includes: obtaining the mainstream nickel solubility and the near-wall nickel solubility of the reactor core; when the near-wall nickel solubility of the reactor core is greater than the mainstream nickel solubility of the reactor core, constructing a mass transfer model of the reactor core based on the mass evaporation rate, the mainstream nickel solubility of the reactor core, the total concentration of corrosion products in the coolant system, and the mass flow rate of fouling deposition transfer in the reactor core; otherwise, constructing a mass transfer model of the reactor core based on the mass evaporation rate, the near-wall nickel solubility of the reactor core, the total concentration of corrosion products in the coolant system, and the mass flow rate of fouling deposition transfer in the reactor core; constructing a mass balance model of the coolant system based on the total concentration of corrosion products in the coolant system, the mass flow rate of fouling deposition transfer in the reactor core, and the total release rate of corrosion products in the coolant system; solving the mass flow rate of fouling deposition transfer in the reactor core based on the mass transfer model of the reactor core and the mass balance model of the coolant system; calculating the fouling thickness of the reactor core based on the mass flow rate of fouling deposition transfer in the reactor core.

[0075] In the embodiments of the present application, by constructing a mass transfer model and a mass balance model, and calculating the fouling thickness based on the mass transfer model and the mass balance model, comprehensively considering the influencing factors of the fouling deposition thickness, the calculation results are more accurate and reasonable.

[0076] In some embodiments, a mass transfer model of the letdown system can also be constructed, specifically including: obtaining the nickel solubility of the letdown system; constructing a mass transfer model of the letdown system based on the nickel solubility of the letdown system and the mass flow rate of fouling deposition transfer in the letdown system. Based on the mass transfer model of the letdown system, the mass flow rate of fouling deposition transfer in the letdown system can be calculated. When constructing the mass balance model, the fouling mass of the letdown system is also considered, that is, a mass balance model of the coolant system is also constructed based on the mass flow rate of fouling deposition transfer in the letdown system, and the mass flow rate of fouling deposition transfer in the reactor core is also solved based on the mass transfer model of the letdown system.

[0077] In some embodiments, a mass transfer model of the steam generator can also be constructed. Specifically, it includes: obtaining the mainstream nickel solubility and the near-wall nickel solubility of the steam generator; when the near-wall nickel solubility of the steam generator is greater than the mainstream nickel solubility of the steam generator, constructing a mass transfer model of the steam generator based on the mass evaporation rate, the mainstream nickel solubility of the steam generator, the total concentration of corrosion products in the coolant system, and the mass flow rate of fouling deposition transfer of the steam generator; otherwise, constructing a mass transfer model of the steam generator based on the mass evaporation rate, the near-wall nickel solubility of the steam generator, the total concentration of corrosion products in the coolant system, and the mass flow rate of fouling deposition transfer of the steam generator. When constructing the mass balance model, the fouling mass of the steam generator is also considered, that is, a mass balance model of the coolant system is also constructed based on the mass flow rate of fouling deposition transfer of the steam generator, and the mass flow rate of fouling deposition transfer of the core is solved based on the mass transfer model of the steam generator.

[0078] In the embodiments of the present application, when constructing the mass balance model, the fouling masses inside the core, inside the steam generator, and inside the letdown system are considered simultaneously, comprehensively considering the locations where fouling is generated and disappears in the primary loop, making the fouling distribution in the entire loop more reasonable and the calculated fouling thickness more accurate. It can be understood that in some other embodiments, the layout in the reactor loop may be different, and the steam generator and / or the letdown system may not be considered, thus omitting the corresponding calculation steps.

[0079] In some embodiments, calculating the fouling thickness of the core based on the mass flow rate of fouling deposition transfer of the core includes: calculating the mass per unit area of fouling of the core based on the mass flow rate of fouling deposition transfer of the core; calculating the fouling thickness of the core based on the mass per unit area of fouling of the core and the fouling density. In the embodiments of the present application, the fouling thickness is directly solved through the given fouling density without separately considering the fouling densities at different positions, reducing the calculation amount.

[0080] In an exemplary embodiment, please refer to Figure 3 , the calculation method of the fouling thickness includes the following steps:

[0081] Step S301, calculate the mainstream nickel solubility of the steam generator.

[0082] The nickel solubility can be obtained by solving a function based on the fluid pH value, and the fluid pH value can be obtained by the table lookup method. In some embodiments, the calculation relationship of the nickel solubility is as follows:

[0083] (4)

[0084] In the formula: Con Ni is the nickel solubility, with the unit of ppb;

[0085] a1, a2, and a3 are coefficients that can be obtained by fitting corrosion experiment data or nuclear power plant operation data, that is, a1, a2, and a3 are input parameters.

[0086] Through the calculation relation (4) of nickel solubility and the mainstream pH value of the steam generator, the mainstream nickel solubility of the steam generator can be calculated. It should be noted that in the above embodiments, the nickel solubility is solved based on the calculation relation related to the fluid pH value. In some other embodiments, it can also be calculated based on other calculation relations of nickel solubility.

[0087] Step S302: Calculate the near-wall nickel solubility of the steam generator.

[0088] Through the calculation relation (4) of nickel solubility and the near-wall fluid pH value of the steam generator, the near-wall nickel solubility Con of the steam generator can be calculated. wall 。

[0089] In some embodiments, the near-wall fluid pH value of the steam generator can be calculated based on the externally input near-wall fluid temperature, boron concentration, and lithium concentration of the steam generator, and this pH value can also be obtained by the look-up table method.

[0090] Step S303: Construct a mass transfer model of the steam generator.

[0091] Based on the nickel solubility of the steam generator, the total concentration of corrosion products in the coolant system, and the mass flow rate of dirt deposition transfer in the steam generator, a mass transfer model of the steam generator is constructed. In an exemplary embodiment, the mass transfer model of the steam generator considers the mass transfer from the mainstream to the near-wall layer caused by turbulence, the mass transfer from the near-wall layer to the wall caused by deposition, etc., and its relation is as follows:

[0092] (5)

[0093] In the formula: Con is the total concentration of corrosion products in the coolant system, and the unit is ppb.

[0094] Con wall is the near-wall nickel solubility of the steam generator, and the unit is ppb. When the near-wall nickel solubility of the steam generator is greater than the mainstream nickel solubility, the near-wall nickel solubility of the steam generator in relation (5) is equal to the mainstream nickel solubility, that is, a mass transfer model of the steam generator is constructed based on the mainstream nickel solubility of the steam generator.

[0095] F is the factor for converting nickel solubility to total concentration.

[0096] is the mass transfer coefficient caused by the mixing of the near-wall layer and the mainstream of the steam generator, and the unit is kg / s·m2 。

[0097] G sed The mass flow rate transferred for fouling deposition in the steam generator, with the unit of kg / s·m 2 ;

[0098] where, F and are input parameters.

[0099] Step S304, calculate the mainstream nickel solubility in the core.

[0100] Based on the fluid temperature, boron concentration, and lithium concentration of the core mainstream input externally, the fluid pH value of the core mainstream can be calculated. Combining with the calculation relationship formula (4) of nickel solubility, the mainstream nickel solubility in the core can be calculated.

[0101] Step S305, calculate the near-wall nickel solubility in the core.

[0102] Based on the fluid temperature, boron concentration, and lithium concentration of the core near-wall input externally, the fluid pH value of the core near-wall can be calculated. Combining with the calculation relationship formula (4) of nickel solubility, the near-wall nickel solubility Con in the core can be calculated. wnode 。

[0103] Step S306, construct the mass transfer model of the core.

[0104] Based on the nickel solubility in the core, the total concentration of corrosion products in the coolant system, and the mass flow rate transferred for fouling deposition in the core, construct the mass transfer model of the core. In an exemplary embodiment, the mass transfer model of the core considers the mass transfer from the mainstream to the near-wall layer caused by boiling at the core, the mass transfer from the mainstream to the near-wall layer caused by turbulence, the mass transfer from the near-wall layer to the wall surface caused by deposition, etc., and its relationship formula is as follows:

[0105] (6)

[0106] In the formula: is the mass transfer efficiency caused by boiling;

[0107] G evap is the mass evaporation rate, with the unit of kg / s·m 2 ;

[0108] Con is the total concentration of corrosion products in the coolant system, with the unit of ppb;

[0109] Con wnodeis the nickel solubility at the near-wall surface of the core, in ppb. When the nickel solubility at the near-wall surface of the core is greater than the mainstream nickel solubility, the nickel solubility at the near-wall surface of the core in relation (6) is equal to the mainstream nickel solubility, that is, the mass transfer model of the core is constructed based on the mainstream nickel solubility of the core.

[0110] F is the factor for converting nickel solubility to total concentration.

[0111] is the mass transfer coefficient caused by the mixing of the near-wall layer of the core and the mainstream, in kg / s·m 2 ;

[0112] G node is the mass flow rate of fouling deposition transfer in the core, in kg / s·m 2 .

[0113] Among them, , G evap , and F are input parameters.

[0114] Step S307: Construct the mass transfer model of the letdown system.

[0115] First, obtain the nickel solubility of the letdown system. Then, construct the mass transfer model of the letdown system based on the nickel solubility of the letdown system and the mass flow rate of fouling deposition transfer in the letdown system. In an exemplary embodiment, based on the fluid temperature, boron concentration, and lithium concentration of the letdown system input externally, the fluid pH value of the letdown system can be calculated. Combining with the calculation relation (4) of nickel solubility, the nickel solubility Con of the letdown system can be calculated. ld . The relation of the mass transfer model of the letdown system is as follows:

[0116] (7)

[0117] In the formula: G ld is the mass flow rate of fouling deposition transfer in the letdown system, in kg / s·m 2 ;

[0118] G sld is the letdown flow rate, in kg / s·m 2 ;

[0119] Con ld is the nickel solubility of the letdown system, in ppb;

[0120] F p is the particulate removal rate;

[0121] F s is the dissolved matter removal rate.

[0122] Among them, Gsld , F p and F s are input parameters.

[0123] Step S308: Construct a mass balance model for the coolant system.

[0124] Based on the total concentration of corrosion products in the coolant system, the mass flow rate of fouling deposition transferred by the steam generator, the mass flow rate of fouling deposition transferred by the reactor core, the mass flow rate of fouling deposition transferred by the letdown system, and the total release rate of corrosion products in the coolant system, construct a mass balance model for the coolant system. In an exemplary embodiment, the mass balance model of the coolant system considers the removal of fouling mass in the letdown system, the total release of corrosion products in the coolant system, the fouling mass in the steam generator, and the fouling mass inside the reactor core. The relationship is as follows:

[0125] (8)

[0126] In the formula: m total is the mass of the coolant in the coolant system, with the unit of kg / m 2 ;

[0127] Con is the total concentration of corrosion products in the coolant system, with the unit of ppb;

[0128] t is the time, with the unit of s;

[0129] G node is the mass flow rate of fouling deposition transferred by the reactor core, with the unit of kg / s·m 2 ;

[0130] G sed is the mass flow rate of fouling deposition transferred by the steam generator, with the unit of kg / s·m 2 ;

[0131] G ld is the mass flow rate of fouling deposition transferred by the letdown system, with the unit of kg / s·m 2 ;

[0132] G dis is the total release rate of corrosion products in the coolant system, with the unit of kg / s·m 2 .

[0133] Among them, m total and G dis are input parameters. It should be noted that when the steam generator is not considered, the parameter G sed in formula (8) can be removed. Similarly, when the letdown system is not considered, the parameter G ld in formula (8) can be removed.

[0134] In some embodiments, the total release rate G of corrosion products in the coolant system may be calculated based on the mainstream pH value in the coolant system. dis An example of the calculation relationship of the total release rate of corrosion products in the coolant system is as follows:

[0135] G dis = f1(PH dis ) (9)

[0136] Where: G dis is the total release rate of corrosion products in the coolant system, with the unit of kg / m 2 ·s;

[0137] f1 is the calculation function of the total release rate of corrosion products in the coolant system;

[0138] PH dis is the mainstream pH value in the coolant system. This parameter is an input parameter.

[0139] It should be noted that in the above embodiments, the total release rate of corrosion products in the coolant system is calculated based on the mainstream pH value in the coolant system. In some other embodiments, it can also be calculated based on other calculation relationships of the total release rate of corrosion products in the coolant system.

[0140] Step S309: Calculate the mass flow rate of fouling deposition transfer.

[0141] Based on the mass transfer model of the steam generator, the mass transfer model of the core, the mass transfer model of the letdown system, and the mass balance model of the coolant system, solve the mass flow rate of fouling deposition transfer of the steam generator and the mass flow rate of fouling deposition transfer of the core. That is, by combining equations (5), (6), (7), and (8), G node can be solved. In some embodiments, G sed and Con can also be solved.

[0142] Step S310: Calculate the fouling thickness.

[0143] Based on the mass flow rate of fouling deposition transfer calculated in step S309, the fouling thickness can be calculated. In some embodiments, the mass per unit area of fouling can be calculated based on the mass flow rate of fouling deposition transfer, and then the fouling thickness can be calculated based on the mass per unit area of fouling and the fouling density. Exemplarily, the fouling thickness of the core can be calculated through equations (10) and (11). Similarly, the fouling thickness of the steam generator can be calculated using G sed .

[0144] (10)

[0145] (11)

[0146] Where: m node is the mass per unit area of fouling in the core, with the unit of kg / m 2 ;

[0147] G node is the mass flow rate of fouling deposition transfer in the core, with the unit of kg / s·m 2 ;

[0148] t is time, with the unit of s;

[0149] THC node is the fouling thickness in the core, with the unit of m;

[0150] is the fouling density, with the unit of kg / m 3 .

[0151] Among them, are input parameters.

[0152] In the embodiment of the present application when constructing the mass transfer model, the inflow and outflow of coolant and corrosion products in the near-wall coolant layer are considered; the effects of turbulent mixing and boiling cycle cause mass transfer; the dissolution and deposition of fouling on the wall; the increase or decrease of corrosion products in the coolant. These processes make the process of fouling deposition from the mainstream coolant to the wall more detailed and the mass transfer process more accurate.

[0153] Step S140: Calculate the boron deposition mass based on the fouling thickness in the core.

[0154] The embodiment of the present application considers two boron deposition mechanisms, namely boron-lithium compound deposition and boron adsorption, and calculates the total boron deposition mass by combining these two deposition mechanisms to achieve the calculation of the boron deposition mass in the reactor, making the calculation result more accurate.

[0155] Please refer to Figure 4 , the calculation method of boron deposition mass includes the following steps:

[0156] Step S401: Calculate the boron-lithium compound deposition mass inside the fouling.

[0157] The deposition mass of the boron-lithium compound inside the fouling, i.e., the mass of the boron-lithium compound deposited inside the fouling, can be calculated through a certain relational formula. For example, it can be calculated based on the thickness of the boron deposition layer, the boron concentration in the coolant, and the concentration factor. The thickness of the boron deposition layer can be obtained based on the thickness of the fouling. For example, it can be obtained based on the relational formula between the thickness of the boron deposition layer inside the fouling and the thickness of the fouling fitted according to the test data, or based on the concentration of the target substances (such as lithium, boron, etc.) inside the fouling and the thickness of the fouling. That is, the deposition mass of the boron-lithium compound inside the fouling can be calculated based on the concentration of the target substances inside the fouling, the thickness of the fouling in the core, the boron concentration in the coolant, and the concentration factor.

[0158] In some embodiments, obtaining the thickness of the boron deposition layer inside the fouling includes: obtaining the concentration of the target substances inside the fouling; when the concentration of the target substances exceeds the solubility of the deposited boron-lithium compound, obtaining the thickness of the boron deposition layer based on the thickness of the fouling. Among them, obtaining the concentration of the target substances inside the fouling includes: dividing the fouling into multiple layers of grids based on the thickness of the fouling; calculating the concentration of the target substances at the grids. Then, when the concentration of the target substances at the grids exceeds the solubility of the deposited boron-lithium compound, obtaining the thickness of the boron deposition layer based on the thickness of the fouling and the position of the grids. In the embodiments of the present application, the fouling is divided into multiple layers of grids and calculated based on the grid units, making the calculation results more accurate.

[0159] Since the microstructure of the fouling is a capillary porous body, the porous body is considered a group of interconnected pipes. The larger vertical pores are filled with steam, and the smaller pores are filled with liquid. Considering the diffusion of the target substance i (such as lithium, boron, etc.) inside the fouling and the convection caused by the mass flow inside the fouling, mass flow convection equations are constructed for the liquid phase and vapor phase inside the fouling respectively, and then the concentration of the target substances inside the fouling is obtained by solving, making the calculation results more accurate. Specifically, calculating the concentration of the target substances at the grids includes: constructing a mass flow convection equation for the liquid phase inside the fouling based on the diffusion coefficient of the target substance in the liquid phase inside the fouling, the total mass flow rate of the target substance at the grids in the liquid phase inside the fouling, and the mass concentration of the target substance at the grids in the liquid phase inside the fouling; constructing a mass flow convection equation for the vapor phase inside the fouling based on the diffusion coefficient of the target substance in the vapor phase inside the fouling, the total mass flow rate of the target substance at the grids in the vapor phase inside the fouling, and the mass concentration of the target substance at the grids in the vapor phase inside the fouling; solving based on the relational formula between the mass flow convection equation for the liquid phase inside the fouling, the mass flow convection equation for the vapor phase inside the fouling, the total mass flow rate of the target substance at the grids in the liquid phase inside the fouling and the total mass flow rate of the target substance at the grids in the vapor phase inside the fouling, and the relational formula between the mass concentration of the target substance at the grids in the liquid phase inside the fouling and the mass concentration of the target substance at the grids in the vapor phase inside the fouling to obtain the mass concentration of the target substance at the grids in the liquid phase inside the fouling, and taking the obtained mass concentration of the target substance at the grids in the liquid phase inside the fouling as the concentration of the target substances at the grids.

[0160] In an exemplary embodiment, obtaining the thickness of the boron deposition layer inside the fouling includes the following steps:

[0161] Step S1111, divide the grid.

[0162] Please refer to Figure 5 , the fouling deposits on the surface of the fuel assembly cladding. In the figure, q w is the heat flux density on the cladding surface, with the unit of W / m 2 ; q surf is the heat flux density on the fouling surface, with the unit of W / m 2 . Based on the given uniform fouling thickness, divide the grid. The thickness of each layer of the grid is dx. Divide the grid from the fouling surface (close to the fluid side) to the cladding surface in sequence and determine the parameters at each layer of the grid. For example, if the fouling is evenly divided into N - 1 layers of grids, these N - 1 layers of grids share N surfaces, namely surfaces x(1), x(2), …, x(n), x(n + 1), …, x(N). Among them, the thickness of the nth layer of the grid is dx(n). It can be understood that considering the relatively large size of the fuel assembly cladding surface, each layer of the grid can also be divided into multiple grids in the direction perpendicular to the fouling deposition direction, and the calculation is performed for each grid in a loop.

[0163] In the above - mentioned embodiment, the grid is divided from the fouling surface to the cladding surface in sequence. In some other embodiments, the grid can also be divided from the cladding surface to the fouling surface in sequence. The embodiments of the present application do not limit this.

[0164] Step S1112, calculate the concentration of the target substance at the grid. Among them, the target substance can contain lithium, or contain boron, or contain both lithium and boron.

[0165] Considering the diffusion of the target substance i (such as lithium, boron, etc.) inside the fouling and the convection caused by the mass flow inside the fouling, mass - flow convection equations are constructed for the liquid phase and vapor phase inside the fouling respectively. Among them, at the kth layer of the grid, the mass - flow convection equations for the liquid phase and vapor phase inside the fouling are exemplified as follows:

[0166] (12)

[0167] (13)

[0168] (14)

[0169] In the formula: G i,l,k is the total mass flow rate of the target substance i at the kth layer of the grid in the liquid phase inside the fouling, with the unit of kg / s·m 2 ;

[0170] D i,lis the diffusion coefficient of the target substance i in the liquid phase within the fouling, with the unit of m 2 / s;

[0171] is the liquid phase density within the fouling, with the unit of kg / m 3 ;

[0172] A l is the flow area of the part filled with the liquid phase within the fouling, with the unit of m 2 ;

[0173] C i,l,k is the mass concentration of the target substance at the k-th grid in the liquid phase within the fouling, with the unit of ppb;

[0174] G l is the mass flow rate of the liquid phase within the fouling, with the unit of kg / s·m 2 ;

[0175] G i,v,k is the total mass flow rate of the target substance i at the k-th grid in the vapor phase within the fouling, with the unit of kg / s·m 2 ;

[0176] D i,v is the diffusion coefficient of the target substance i in the vapor phase within the fouling, with the unit of m 2 / s;

[0177] is the vapor phase density within the fouling, with the unit of kg / m 3 ;

[0178] A v is the flow area of the part filled with the vapor phase within the fouling, with the unit of m 2 ;

[0179] C i,v,k is the mass concentration of the target substance at the k-th grid in the vapor phase within the fouling, with the unit of ppb;

[0180] G v is the mass flow rate of the vapor phase within the fouling, with the unit of kg / s·m 2 ;

[0181] B is a coefficient.

[0182] Among them, D i,l , , A l , G l , D i,v , , G v , A v and B are input parameters.

[0183] At steady state, the net flux of the target substance i at any position within the fouling is 0, and the following relationship holds at the k-th layer of the grid:

[0184] G i,l,k + G i,v,k = 0 (15)

[0185] By simultaneously solving the relationships (12)-(15), C i,l,k can be obtained. Take the obtained C i,l,k as the concentration of the target substance i at the grid.

[0186] Step S1113: Calculate the thickness of the boron deposition layer.

[0187] Based on the concentration distribution of the target substance i within the fouling, determine whether the target substance i at the grid reaches the critical concentration, that is, whether it exceeds the solubility of the deposited boron-lithium compound, so as to obtain the thickness d B of the boron deposition layer. An example of the calculation formula for the thickness d B of the boron deposition layer is as follows:

[0188] (16)

[0189] In the formula, d B is the thickness of the boron deposition layer, with the unit of m;

[0190] THC node is the fouling thickness of the reactor core, with the unit of m;

[0191] x(j) is the thickness from the fouling surface to the j-th layer of the grid assuming that the concentration of the target substance i at the j-th layer of the grid exceeds the solubility of the deposited boron-lithium compound. If the thickness of each layer of the grid is dx, then the thickness from the fouling surface to the j-th layer of the grid is j·dx, with the unit of m.

[0192] The solubility of the deposited boron-lithium compound can be used as an input parameter or obtained based on the temperature within the fouling. Exemplarily, the solubility of the deposited boron-lithium compound is obtained based on an empirical relationship related to temperature:

[0193] (17)

[0194] In the formula, SI is the solubility of the deposited boron-lithium compound, with the unit of ppb;

[0195] T crud is the temperature within the fouling, with the unit of °C.

[0196] Among them, T crud is an input parameter.

[0197] Calculate the deposition mass of boron-lithium compounds based on the thickness of the boron deposition layer, the boron concentration in the coolant, and the concentration factor. In some embodiments, an example of the calculation formula for the deposition mass of boron-lithium compounds is as follows:

[0198] (18)

[0199] In the formula, M b,1 is the deposition mass of boron-lithium compounds, with the unit of kg;

[0200] k B is the concentration factor;

[0201] C coolant B is the boron concentration in the coolant, with the unit of ppb;

[0202] is the coolant density, with the unit of kg / m 3 ;

[0203] A is the area of the calculated region on the cladding surface, with the unit of m 2 .

[0204] Among them, k B , C coolant B , and A are input parameters.

[0205] Step S402: Calculate the boron adsorption mass on the fouled surface.

[0206] The boron adsorption mass on the fouled surface is the mass of boron adsorbed on the fouled surface, which can be calculated through a certain relationship. For example, the boron adsorption mass can be calculated based on the boron concentration on the cladding surface and the fouling mass in the core. In some embodiments, calculating the boron adsorption mass on the fouled surface includes the following steps:

[0207] Step S1201: Obtain the boron concentration on the cladding surface.

[0208] The boron concentration on the cladding surface can be an input parameter or calculated based on relevant calculation formulas. In an exemplary embodiment, the boron concentration on the cladding surface is calculated based on the enrichment factor and the boron concentration in the coolant. An example of the calculation relationship is as follows:

[0209] (19)

[0210] (20)

[0211] In the formula, C clad B is the boron concentration on the cladding surface, with the unit of ppb;

[0212] C coolant B is the boron concentration in the coolant, with the unit of ppb;

[0213] CF B is the concentration factor;

[0214] G evap is the mass evaporation rate, with the unit of kg / s·m 2 ;

[0215] THC node is the fouling thickness of the core, with the unit of m;

[0216] is the coolant density, with the unit of kg / m 3 ;

[0217] D B is the boron diffusion coefficient in the coolant, with the unit of m 2 / s;

[0218] P is the fouling porosity.

[0219] Among them, C coolant B , , D B , P are input parameters.

[0220] Step S1202: Calculate the boron adsorption mass based on the boron concentration and fouling mass on the cladding surface.

[0221] In actual operation, based on the relationship between the boron adsorption mass and the boron concentration and fouling mass on the cladding surface fitted according to test data, the mass of boron adsorbed on the fouling surface can be obtained. An example of the relationship is as follows:

[0222] (21)

[0223] Among them, b1, b2, and b3 are constants and are input parameters;

[0224] M crud is the fouling mass, with the unit of kg, and can be calculated based on the fouling mass per unit area m node of the core;

[0225] C clad B is the boron concentration on the cladding surface, with the unit of ppb.

[0226] Step S403: Calculate the boron deposition mass.

[0227] Sum the boron-lithium compound deposition mass and the boron adsorption mass to obtain the boron deposition mass M b , and an example of the relationship is as follows:

[0228] (22)

[0229] Step S150: Determine whether the convergence condition is met. If the convergence condition is met, proceed to step S160; otherwise, return to step S110.

[0230] Determining whether the convergence condition is met may include determining whether the target parameter satisfies the convergence condition, or may include determining whether the time point reaches the end time, or may include both the determination of the target parameter and the determination of the time point. If the target parameter does not satisfy the convergence condition, update the input parameter for iterative calculation; otherwise, proceed to the next step. If the time point reaches the end time, end; otherwise, set the value of the input parameter for subsequent calculation. After each time step calculation is completed, update the current calculated time point. It should be noted that the values of some input parameters are related to time, and time interpolation can be performed on these input parameters, such as the boron and lithium concentrations in the coolant.

[0231] The target parameter may include one or more parameters obtained during the calculation process. In an exemplary embodiment, the target parameters include the calculated mass evaporation rate and the total concentration of corrosion products in the coolant system, making the solution process more accurate. Specifically, determining whether the target parameter satisfies the convergence condition includes: determining whether the error between the calculated mass evaporation rate and the first preset value is less than the first threshold; determining whether the error between the calculated total concentration of corrosion products in the coolant system and the second preset value is less than the second threshold. When the error between the calculated mass evaporation rate and the first preset value is less than the first threshold, the convergence condition is met. When the error between the calculated total concentration of corrosion products in the coolant system and the second preset value is less than the second threshold, the convergence condition is met. When both the mass evaporation rate and the total concentration of corrosion products in the coolant system satisfy the convergence condition, the current time step calculation is completed; otherwise, update the input parameter for iterative calculation of the current time step.

[0232] Step S160: Output the calculation result.

[0233] In some embodiments, the calculation result of the current time step may be output after each time step calculation is completed. In some other embodiments, the calculation result may be output after the time point reaches the end time and all time step calculations are completed. The form of outputting the calculation result can be set as needed.

[0234] Figure 6 It is a schematic structural diagram of a computing device provided by an embodiment of the present application. As Figure 6 shown, the computing device 600 includes one or more processors 610, one or more memories 620 coupled to the processor 610, and one or more communication modules 640 coupled to the processor 610.

[0235] The communication module 640 is used for two-way communication. The communication module 640 has at least one antenna to facilitate communication. The communication interface can represent any interface necessary for communicating with other network elements.

[0236] The processor 610 can be of any type suitable for the local technical network and, by way of non-limiting example, can include one or more of the following: general-purpose computer, dedicated computer, microprocessor, digital signal processor (DSP), and a processor based on a multi-core processor architecture. The computing device 600 can have multiple processors, such as an application-specific integrated circuit chip that is clocked to synchronize with the main processor in a timely manner.

[0237] The memory 620 can include one or more non-volatile memories and one or more volatile memories. Examples of non-volatile memories include, but are not limited to, read-only memory (ROM) 624, electrically programmable read-only memory (EPROM), flash memory, hard disk, compact disc (CD), digital video disc (DVD), and other magnetic and / or optical memories. Examples of volatile memories include, but are not limited to, random access memory (RAM) 622 and other volatile memories that do not persist during a power outage.

[0238] The computer program 630 includes computer-executable instructions to be executed by the associated processor 610. The computer program 630 can be stored in the ROM 624. The processor 610 can perform any appropriate actions and processes by loading the computer program 630 into the RAM 622.

[0239] Embodiments of the present application can be implemented by the computer program 630 such that the computing device 600 can execute any of the disclosed processes discussed with reference to Figures 1-4 Embodiments of the present application can also be implemented by hardware or by a combination of software and hardware.

[0240] In some embodiments, the computer program 630 can be tangibly embodied in a computer-readable medium, which can be included in the computing device 600 (e.g., the memory 620) or other storage devices accessible to the computing device 600. The computing device 600 can load the computer program 630 from the computer-readable medium into the RAM 622 for execution. The computer-readable medium can include any type of tangible non-volatile memory, such as ROM, EPROM, flash memory, hard disk, CD, DVD, etc. The computer program 630 is stored on the computer-readable medium.

[0241] In general, the various embodiments of the present application can be implemented in hardware or dedicated circuits, software, logic, or any combination thereof. Certain aspects can be implemented in hardware, while other aspects can be implemented in firmware or software, which can be executed by a controller, a microprocessor, or other computing devices. Although the various aspects of the embodiments of the present application are shown and described as block diagrams, flowcharts, or using some other graphical representation, it should be understood that, by way of non-limiting example, the blocks, devices, systems, techniques, or methods described herein can be implemented in hardware, software, firmware, dedicated circuits or logic, general hardware or controllers or other computing devices, or some combination thereof.

[0242] The present application also provides at least one computer program product tangibly stored on a non-transitory computer-readable storage medium. The computer program product includes computer-executable instructions, such as instructions included in program modules, which are executed in a device on a target real or virtual processor to perform the methods described above with reference to Figures 1-4 the methods described. In general, program modules include routines, programs, libraries, objects, classes, components, data structures, etc. that perform specific tasks or implement specific abstract data types. In various embodiments, the functions of the program modules can be combined or separated as needed among the program modules. The machine-executable instructions for the program modules can be executed within a local or distributed device. In a distributed device, the program modules can be located in local and remote storage media.

[0243] The program code for performing the methods of the present application can be written in any combination of one or more programming languages. The program code can be provided to a processor or controller of a general-purpose computer, a dedicated computer, or other programmable data processing device, such that when the program code is executed by the processor or controller, the functions / operations specified in the flowchart and / or block diagram are implemented. The program code can be executed entirely on the machine as an independent software package, partially on the machine, partially on the machine, partially on a remote machine, partially on a remote machine, or entirely on a remote machine or server.

[0244] In the context of the present application, the computer program code or related data can be carried by any suitable carrier to enable a device, apparatus, or processor to perform the various processes and operations described above. Examples of carriers include signals, computer-readable media, etc.

[0245] A computer-readable medium can be a computer-readable signal medium or a computer-readable storage medium. The computer-readable medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. More specific examples of the computer-readable storage medium include an electrical connection having one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0246] Furthermore, although the operations are described in a particular order, this should not be construed as requiring that the operations be performed in the particular order shown or in sequential order, or that all of the illustrated operations be performed, to obtain the desired result. In some cases, multitasking and parallel processing may be advantageous. Similarly, although several specific implementation details are included in the foregoing discussion, these details should not be construed as limitations on the scope of the present application, but rather as descriptions of features specific to particular embodiments. Certain features that are described in the context of separate embodiments may also be implemented in combination in a single embodiment. Conversely, various features that are described in the context of a single embodiment may also be implemented separately or in any suitable subcombination in multiple embodiments.

[0247] Although the present application has been described in language specific to structural features and / or methodological acts, it is to be understood that the application defined in the appended claims is not necessarily limited to the specific features or acts described above. Rather, the specific features and acts described above are disclosed as example forms of implementing the claims.

Claims

1. A method for calculating the fouling risk of a reactor, characterized in that: include: Step S110, obtaining input parameters; Step S120, calculating the mass evaporation rate based on the input parameters; Step S130, calculating the fouling thickness of the core based on the mass evaporation rate; Step S140, calculating the boron deposition mass based on the fouling thickness of the core; Step S150, determine whether the convergence condition is met, if the convergence condition is met, proceed to step S160, otherwise return to step S110; Step S160, outputting the calculation result; Wherein, calculating the mass evaporation rate based on the input parameters comprises: The momentum equation is constructed based on the heat flux density transferred by heat conduction, the mass evaporation rate, the vapor phase enthalpy of the fouling surface, and the liquid phase enthalpy of the fouling surface; An energy equation is constructed based on the vapor pressure in the dirt, the liquid pressure in the dirt, the surface tension of the liquid inside the dirt, the dirt contact angle and the critical pore radius of the dirt; The heat transfer equation is constructed based on the internal temperature of the dirt, the saturated temperature of the liquid inside the dirt, the mass evaporation rate and the latent heat of vaporization of the liquid inside the dirt; Calculating a mass evaporation rate based on the momentum equation, the energy equation, and the heat transfer equation; And the calculating the fouling thickness of the core based on the mass evaporation rate comprises: Obtain the mainstream nickel solubility in the core and the nickel solubility near the wall; In the case where the nickel solubility near the core wall is greater than the mainstream nickel solubility of the core, a mass transfer model of the core is constructed based on the mass evaporation rate, the mainstream nickel solubility of the core, the total concentration of corrosion products in the coolant system, and the mass flow rate of fouling deposition transfer in the core; otherwise, a mass transfer model of the core is constructed based on the mass evaporation rate, the nickel solubility near the core wall, the total concentration of corrosion products in the coolant system, and the mass flow rate of fouling deposition transfer in the core; Constructing a mass balance model of the coolant system based on the total concentration of corrosion products in the coolant system, the mass flow rate of fouling deposits transferred by the core, and the total release rate of corrosion products in the coolant system; Solving the mass flow rate of fouling deposition transfer in the core based on the mass transfer model of the core and the mass balance model of the coolant system; calculating the fouling thickness of the core based on the mass flow rate of fouling deposition transfer of the core; And the calculation of the boron deposition mass based on the fouling thickness of the core includes: Calculating the mass of boron-lithium compound deposition inside the fouling based on the concentration of the target substance inside the fouling, the fouling thickness of the core, the boron concentration in the coolant, and the concentration factor; The boron adsorption mass on the fouling surface is calculated based on the boron concentration on the cladding surface and the mass per unit area of ​​the fouling in the core; The boron deposition mass is obtained by summing the boron-lithium compound deposition mass and the boron adsorption mass.

2. The method according to claim 1, characterized in that The calculating the fouling thickness of the core based on the mass evaporation rate further comprises: Obtain nickel solubility in the downflow system; A mass transfer model of the downflow system is constructed based on the nickel solubility of the downflow system and the mass flow rate of the fouling deposition transferred by the downflow system; A mass balance model of the coolant system is also constructed based on the mass flow rate of fouling deposits transferred by the downcomer system, and the mass flow rate of fouling deposits transferred by the core is solved based on the mass transfer model of the downcomer system.

3. The method according to claim 1, characterized in that The calculating the fouling thickness of the core based on the mass evaporation rate further comprises: Obtain the mainstream nickel solubility and near-wall nickel solubility of the steam generator; In the case where the nickel solubility near the wall of the steam generator is greater than the mainstream nickel solubility of the steam generator, a mass transfer model of the steam generator is constructed based on the mass evaporation rate, the mainstream nickel solubility of the steam generator, the total concentration of corrosion products in the coolant system, and the mass flow rate of fouling deposition transfer of the steam generator; otherwise, a mass transfer model of the steam generator is constructed based on the mass evaporation rate, the nickel solubility near the wall of the steam generator, the total concentration of corrosion products in the coolant system, and the mass flow rate of fouling deposition transfer of the steam generator; A mass balance model of a coolant system is also constructed based on the mass flow rate transferred by the fouling deposits of the steam generator, and the mass flow rate transferred by the fouling deposits of the core is solved based on the mass transfer model of the steam generator.

4. The method according to claim 1, characterized in that The calculating the fouling thickness of the core based on the mass evaporation rate further comprises: The total release rate of corrosion products in the coolant system is calculated based on the mainstream pH value in the coolant system.

5. The method according to claim 1, characterized in that The calculating the fouling thickness of the core based on the mass flow rate of the fouling deposition transfer of the core comprises: Calculating the fouling mass per unit area of ​​the core based on the mass flow rate of fouling deposition transferred by the core; The fouling thickness of the core is calculated based on the fouling mass per unit area and the fouling density of the core.

6. A computing device, characterized in that: include: at least one processor; as well as At least one memory having instructions stored thereon, which, when executed individually or collectively by the at least one processor, cause the computing device to perform the method according to any one of claims 1 to 5.

7. A computer storage medium, characterized in that: The computer storage medium stores instructions, which, when executed individually or collectively by at least one processor of a computing device, cause the computing device to perform the method according to any one of claims 1 to 5.

Images