Calculation method for analyzing three-dimensional temperature field in pressurized water reactor guide tube absorber

Abstract

The invention provides a calculation method for analyzing a three-dimensional temperature field in a pressurized water reactor guide pipe absorber, and relates to the field of nuclear power plant thermotechnical hydraulic engineering design, and the calculation method comprises the following steps: firstly, calculating the coolant temperature in a guide pipe according to a preset guide pipe coolant calculation formula and guide pipe information; secondly, obtaining a composition structure of an absorber material, and establishing a one-dimensional cylinder radial heat conduction model and a model boundary condition about the composition structure based on the composition structure of the absorber material; then parameter information of the composition structure of the absorber material is obtained, and temperature field information of the absorber is calculated based on a one-dimensional cylinder radial heat conduction model according to the parameter information and the coolant temperature in the guide pipe; and finally, judging whether the temperature field information of the absorber meets a preset temperature limit value requirement or not, and obtaining the cooling effect of the guide pipe on the absorber based on a judgment result. According to the embodiment, the temperature distribution of the internal structure of the guide pipe absorber can be accurately obtained, so that data support is provided for thermal hydraulic design of the absorber in the guide pipe.

Description

Technical Field

[0001] The present disclosure relates to the field of thermal hydraulic design of nuclear power plants, and in particular to a calculation method for analyzing a three-dimensional temperature field in a guide tube absorber of a pressurized water reactor. Background Art

[0002] Thermal hydraulic design plays a pivotal role in the design and operation of nuclear reactors. It is not only directly related to the cooling efficiency and operational safety of the reactor, but also has a profound impact on the economy and reliability of the reactor. Thermal hydraulic design covers two core aspects: hydraulics and critical heat flux density characteristics. It aims to ensure that the fuel rods are fully cooled by optimizing the fluid flow and heat transfer process, while accurately controlling the flow rate by the guide tubes, maintaining a reasonable pressure drop and lifting force to meet complex operating requirements. However, in actual operation, the close coupling between thermal hydraulic design and the mechanical design of fuel assemblies makes the optimization between the two a complex and time-consuming iterative process. In particular, the management of the flow by the guide tubes must ensure sufficient cooling flow to meet the core cooling requirements and prevent the coolant from boiling in the guide tubes.

[0003] Currently, although a variety of methods and technologies have been applied in practice for the management of bypass flow in guide tubes of pressurized water reactors, they are still insufficient in providing accurate three-dimensional temperature distribution information of absorber materials in guide tubes. This information is crucial for accurately evaluating the cooling performance of guide tubes and predicting potential thermal-hydraulic risks. Due to the complex flow and heat transfer processes in guide tubes, traditional calculation and analysis methods often find it difficult to fully capture the details of temperature distribution, resulting in certain deviations in design results. In addition, the lack of accurate three-dimensional temperature field data may also make it difficult to achieve the best optimization between thermal-hydraulic design and fuel assembly mechanical design, thereby affecting the overall performance of the reactor. Summary of the invention

[0004] The disclosed embodiments at least provide a calculation method for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor, which can accurately obtain the temperature distribution of the internal structure of the guide tube absorber, thereby providing data support for the thermal-hydraulic design of the absorber in the guide tube.

[0005] The embodiment of the present disclosure provides a calculation method for analyzing the three-dimensional temperature field in a guide tube absorber of a pressurized water reactor, comprising:

[0006] Calculate the coolant temperature in the guide tube according to a preset guide tube coolant calculation formula and guide tube information;

[0007] Acquire the composition structure of the absorber material, and establish a one-dimensional cylindrical radial heat conduction model and model boundary conditions about the composition structure based on the composition structure of the absorber material;

[0008] Acquiring parameter information of the composition structure of the absorber material, and calculating the temperature field information of the absorber based on the one-dimensional cylindrical radial heat conduction model according to the parameter information and the temperature of the coolant in the guide tube;

[0009] It is determined whether the temperature field information of the absorber meets the preset temperature limit requirement, and based on the determination result, the cooling effect of the guide tube on the absorber is obtained.

[0010] The embodiment of the present disclosure provides a computing device for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor, comprising:

[0011] A coolant solution module, used to calculate the coolant temperature in the guide tube according to a preset guide tube coolant calculation formula and guide tube information;

[0012] A model building module, used for obtaining the composition structure of the absorber material, and building a one-dimensional cylindrical radial heat conduction model and model boundary conditions about the composition structure based on the composition structure of the absorber material;

[0013] A temperature field solving module, used for obtaining parameter information of the composition structure of the absorber material, and calculating the temperature field information of the absorber based on the one-dimensional cylindrical radial heat conduction model according to the parameter information and the temperature of the coolant in the guide tube;

[0014] The cooling effect determination module is used to determine whether the temperature field information of the absorber meets the preset temperature limit requirements, and obtain the cooling effect of the guide tube on the absorber based on the judgment result.

[0015] An embodiment of the present disclosure provides a computer device, including: a processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the computer device is running, the processor and the memory communicate via the bus, and when the machine-readable instructions are executed by the processor, a calculation method for three-dimensional temperature field analysis in an absorber of a guide tube of a pressurized water reactor as described in any possible implementation manner described above is performed.

[0016] An embodiment of the present disclosure provides a computer-readable storage medium having a computer program stored thereon. The computer program is executed by a processor to implement a calculation method for analyzing a three-dimensional temperature field in a pressurized water reactor guide tube absorber as described in any of the above possible implementations.

[0017] The calculation method for three-dimensional temperature field analysis in the guide tube absorber of a pressurized water reactor provided in the embodiment of the present disclosure first calculates the temperature of the coolant in the guide tube according to the preset guide tube coolant calculation formula and guide tube information; then, obtains the composition structure of the absorber material, and establishes a one-dimensional cylindrical radial heat conduction model and its boundary conditions based on it; uses the model, combines the parameter information of the absorber material and the coolant temperature, calculates the temperature field information inside the absorber; finally, by judging whether the temperature field information of the absorber meets the preset temperature limit requirements, the cooling effect of the guide tube on the absorber is evaluated. In this way, the present disclosure can obtain the temperature distribution inside the guide tube absorber more accurately and comprehensively, thereby providing data support for the thermal hydraulic design of the absorber in the guide tube.

[0018] In order to make the above-mentioned objectives, features and advantages of the present disclosure more obvious and easy to understand, preferred embodiments are specifically cited below and described in detail with reference to the attached drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the technical solutions of the embodiments of the present disclosure, the following briefly introduces the drawings required to be cited in the embodiments. The drawings herein are incorporated into the specification and constitute a part of the specification. These drawings illustrate embodiments consistent with the present disclosure and are used together with the specification to illustrate the technical solutions of the present disclosure. It should be understood that the following drawings only illustrate certain embodiments of the present disclosure and should not be regarded as limiting the scope. For ordinary technicians in this field, other relevant drawings can be obtained based on these drawings without creative work.

[0020] Figure 1 A flow chart of a calculation method for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor provided by an embodiment of the present disclosure is shown;

[0021] Figure 2 A schematic diagram of an infinite length tube model in a calculation method for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor provided by an embodiment of the present disclosure is shown;

[0022] Figure 3 A schematic diagram of an absorber structure in a calculation method for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor provided by an embodiment of the present disclosure is shown;

[0023] Figure 4 A schematic diagram of the structure of a computing device for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor provided by an embodiment of the present disclosure is shown;

[0024] Figure 5 A schematic diagram of the structure of a computer device provided by an embodiment of the present disclosure is shown. DETAILED DESCRIPTION

[0025] In order to make the purpose, technical scheme and advantages of the embodiments of the present disclosure clearer, the technical scheme in the embodiments of the present disclosure will be clearly and completely described below in conjunction with the drawings in the embodiments of the present disclosure. Obviously, the described embodiments are only part of the embodiments of the present disclosure, rather than all of the embodiments. The components of the embodiments of the present disclosure generally described and shown in the drawings here can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present disclosure provided in the drawings is not intended to limit the scope of the present disclosure for protection, but merely represents the selected embodiments of the present disclosure. Based on the embodiments of the present disclosure, all other embodiments obtained by those skilled in the art without making creative work belong to the scope of protection of the present disclosure.

[0026] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, further definition and explanation thereof is not required in subsequent drawings.

[0027] The term "and / or" herein only describes an association relationship, indicating that three relationships may exist. For example, A and / or B may represent the following three situations: A exists alone, A and B exist at the same time, and B exists alone. In addition, the term "at least one" herein represents any combination of at least two of any one or more of a plurality of. For example, including at least one of A, B, and C may represent including any one or more elements selected from the set consisting of A, B, and C.

[0028] In the design of the reactor, thermal hydraulic design is an extremely important content. Its design content includes two aspects: hydraulics and critical heat flux density characteristics. It is mainly to meet the requirements of sufficient cooling of fuel rods, flow restrictions by guide tubes, related pressure drop and lifting force requirements, etc. In order to meet the requirements of thermal hydraulics and mechanical properties, there is an iterative process between thermal hydraulic design and fuel assembly mechanical design to achieve the best compromise between thermal hydraulic design and mechanical design. For the bypass flow of the guide tube, this part of the flow is used to cool the control rods, burnable poison rods, resistance plug rods and other related component rods of the reactor core, and leaks through the end plug of the guide tube. Generally, it is required that the total bypass flow of the guide tube does not exceed the design limit on the one hand, so as to ensure that the core is fully cooled; on the other hand, it is required that the coolant in the guide tube does not boil.

[0029] Research has found that although a variety of methods and technologies have been applied in practice for the management of bypass flow in pressurized water reactor guide tubes, they are still insufficient in providing accurate three-dimensional temperature distribution information of absorber materials in guide tubes. This information is crucial for accurately evaluating the cooling performance of guide tubes and predicting potential thermal-hydraulic risks. Due to the complex flow and heat transfer processes in guide tubes, traditional calculation and analysis methods often find it difficult to fully capture the details of temperature distribution, resulting in certain deviations in design results. In addition, the lack of accurate three-dimensional temperature field data may also make it difficult to achieve the optimal optimization between thermal-hydraulic design and fuel assembly mechanical design, thereby affecting the overall performance of the reactor.

[0030] Based on the above research, a calculation method for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor is provided in an embodiment of the present disclosure. First, the temperature of the coolant in the guide tube is calculated according to a preset guide tube coolant calculation formula and guide tube information; then, the composition structure of the absorber material is obtained, and based on this, a one-dimensional cylindrical radial heat conduction model and its boundary conditions are established; using this model, combined with the parameter information of the absorber material and the coolant temperature, the temperature field information inside the absorber is calculated; finally, by judging whether the temperature field information of the absorber meets the preset temperature limit requirements, the cooling effect of the guide tube on the absorber is evaluated.

[0031] In the disclosed embodiment, the absorber temperature field information calculated by combining the one-dimensional cylindrical radial heat conduction model with the composition and structural parameter information of the absorber material can more accurately and comprehensively obtain the temperature distribution inside the guide tube absorber, thereby providing data support for the thermal-hydraulic design of the absorber in the guide tube.

[0032] To facilitate understanding of this embodiment, the execution subject of the calculation method for three-dimensional temperature field analysis in the absorber of the guide tube of a pressurized water reactor provided by the embodiment of the present disclosure is first introduced in detail. The execution subject of the calculation method for three-dimensional temperature field analysis in the absorber of the guide tube of a pressurized water reactor provided by the embodiment of the present disclosure is a computer device. The computer device can be a server. Among them, the server can be an independent physical server, or a server cluster or distributed system composed of multiple physical servers, or a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud storage, big data and artificial intelligence platforms.

[0033] The following is a detailed description of the calculation method for analyzing the three-dimensional temperature field in the absorber of the guide tube of a pressurized water reactor provided in the embodiment of the present application in conjunction with the accompanying drawings. Figure 1 As shown, it is a flow chart of a calculation method for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor provided by an embodiment of the present disclosure, and the method includes the following S101 to S104:

[0034] S101, calculating the coolant temperature in the guide tube according to a preset guide tube coolant calculation formula and guide tube information.

[0035] It is understandable that in a pressurized water reactor (PWR), the guide tube is an important part of the nuclear reactor control system. The guide tube is used to provide precise guidance and positioning for the control rod, ensuring that the control rod can be accurately and stably inserted or withdrawn from the reactor core. Among them, the coolant temperature in the guide tube is an important parameter, which directly affects the thermal hydraulic performance and safety of the reactor. The temperature of the coolant in the guide tube is determined by the power of the guide tube, and the power of the guide tube includes: guide tube wall power, control rod power and coolant power.

[0036] Specifically, the preset guide tube coolant calculation formula includes a coolant power calculation formula and a coolant enthalpy calculation formula, and the coolant power calculation formula includes a guide tube wall power calculation formula, a control rod power calculation formula, and a coolant power calculation formula; the guide tube information may include guide tube wall information, control rod information, and coolant information. Therefore, when calculating the coolant temperature in the guide tube according to the preset guide tube coolant calculation formula and the guide tube information, the following (a) to (b) may be included:

[0037] (a) calculating the guide tube wall power based on the guide tube wall power calculation formula and the guide tube wall information; calculating the control rod power based on the control rod calculation formula and the control rod information; and calculating the coolant power based on the coolant power calculation formula and the coolant information;

[0038] (b) According to the guide tube wall power, control rod power and coolant power, the enthalpy value of the guide tube segment is calculated based on the coolant enthalpy value calculation formula, and the coolant temperature is determined according to the enthalpy value.

[0039] Here, the guide tube wall information may include the volume heat release rate of the guide tube wall, the radial length of the inner wall of the guide tube, the radial length of the outer wall of the guide tube, and the thickness of the guide tube wall. The guide tube wall power is obtained from the surface heat flux density of the inner wall of the guide tube. The long tube model of the guide tube is Figure 2 Taking the infinite length tube model shown in the figure as an example, the model uses a cylindrical coordinate system, does not consider the axial heat conduction of the tube wall, and assumes that the radial heat conduction of the tube is uniform, so the heat conduction of the tube is only a one-dimensional equation of its tube diameter. Therefore, the heat conduction differential equation containing the heat source guided to the tube wall can be expressed as:

[0040]

[0041] Its boundary condition is: when r = r i Sometimes, there are When r = r oSometimes, there is Where Φ represents the volume heat release rate of the guide tube wall, expressed in kW / m 3 ; r i The radial length of the inner wall of the guide tube, in m; r o The radial length of the outer wall of the guide tube, in m; h i The surface heat transfer coefficient of the inner wall of the guide tube is expressed in kW / (m 2 .K);h o The surface heat transfer coefficient of the outer wall of the guide tube is expressed in kW / (m 2 .K);t fi The temperature of the fluid in the guide tube, in K; t fo The temperature of the fluid outside the guide tube is in K. λ is the thermal conductivity of the guide tube wall in kW / (mK).

[0042] Here, the heat conduction differential equation with the heat source guided to the tube wall is integrated twice to obtain:

[0043]

[0044] C1 and C2 can be obtained according to the above boundary conditions:

[0045]

[0046] From this solution we can get:

[0047]

[0048] Where:

[0049]

[0050] The temperature distribution of the guide tube wall is:

[0051]

[0052] The heat flux density of the guide tube wall to the guide tube is:

[0053]

[0054] For example, after obtaining the heat flux density of the guide tube wall to the guide tube, the guide tube wall power can be obtained by solving the guide tube wall power calculation formula, and the guide tube wall power calculation formula is as follows:

[0055] Q 1,i =q×Ph×z;

[0056] Among them, Q 1,i It is represented by the guide tube wall power of the i-th block; q is represented by the heat flux density of the guide tube wall, in W / m2 ; Ph represents the heating circumference of the inner wall of the guide tube, in m; z represents the calculated length of the node, in m.

[0057] For example, the control rod information may include the volume heat release rate of the control rod absorber material, the directness of the control rod absorber material, and the outer diameter of the control rod, and the control rod calculation formula includes:

[0058]

[0059] Among them, Q 2,i Expressed as the control rod power of the i-th node; Φ abs Indicates the volume heat release rate of the control rod absorption material, in kW / m 3 ;D abs Indicates the diameter of the control rod absorber material, in m; Φ clad Indicates the volume heat release rate of the control rod cladding metal, in kW / m 3 ;D R Represents the outer diameter of the control rod, in meters.

[0060] For example, the coolant information may include the inner diameter of the guide tube, the outer diameter of the absorber, and the coolant volume heat release rate, etc. The coolant power calculation formula is expressed as:

[0061]

[0062] Among them, Q 3,i It is expressed as the coolant power of the i-th block; D i It is expressed as the inner diameter of the guide tube of the i-th block, in m; D R Expressed as the outer diameter of the absorber, in m; Φ wat Expressed as the volume heat release rate of the coolant, in kW / m 3 ; NABS=1 means an absorber is inserted; NABS=0 means no absorber is inserted.

[0063] It can be understood that after obtaining the guide tube wall power, control rod power and coolant power, the enthalpy of the guide tube segment can be calculated according to the coolant enthalpy calculation formula. The enthalpy is used to determine the coolant temperature, which usually involves looking up the coolant physical property table according to the enthalpy. If the enthalpy is higher than the liquid saturation enthalpy of water, the temperature can be taken as the saturation temperature. Here, the coolant enthalpy calculation formula includes:

[0064]

[0065] Q i =Q 1,i +Q 2,i +Q 3,i ;

[0066] Among them, H i It is expressed as the outlet enthalpy of the ith node, in kJ / kg, H i-1 It is expressed as the outlet enthalpy value of the previous node of the i-th node, in kJ / kg, with the initial value H0 being 0; Q i represents the total power absorbed by the coolant of the i-th node; W i It is expressed as the flow rate of the i-th block in the guide pipe, in kg / s. Here,

[0067] It can be understood that after determining the coolant temperature according to the enthalpy value, it is also necessary to calculate the absorber outer wall temperature according to the coolant temperature based on the absorber outer wall temperature solution formula, and the absorber outer wall temperature solution formula is expressed as:

[0068]

[0069] Among them, T w Expressed as the outer wall temperature of the absorber, in °C; T f It is expressed as the coolant temperature in °C; q is expressed as the heat flux density of the guide tube wall in W / m 2 ; h represents the absorber surface heat transfer coefficient, unit is m 2 ℃ / W. Here, the calculation formula for h is:

[0070]

[0071] Where Nu is the Nusselt number, which is a dimensionless number; D eIB Expressed as equivalent diameter, unit is m.

[0072] S102, obtaining a composition structure of an absorber material, and establishing a one-dimensional cylindrical radial heat conduction model and model boundary conditions related to the composition structure based on the composition structure of the absorber material.

[0073] Here, see Figure 3 As shown, Figure 3 Taking the absorber structure as an example, the composition structure of the absorber material proposed in the present disclosure includes a first material, a second material, a third material, a core block and a cladding; according to the composition structure of the absorber material, a corresponding one-dimensional cylindrical radial heat conduction model can be constructed, which is expressed as:

[0074]

[0075] Where r is the radial coordinate in m; k(T) is the thermal conductivity in W / M / ℃; Φ is the volume heat release rate in W / m 3 .

[0076] Exemplarily, the model boundary conditions constituting the structure include a first material boundary condition, a second material boundary condition, a third material boundary condition, a core block boundary condition, and a cladding boundary condition.

[0077] Here, the first material boundary condition includes:

[0078]

[0079] Among them, r material_1 Indicates the radius of the first material, in m; T Omaterial_1 Expressed as the temperature of the first material in °C.

[0080] Here, the second material boundary condition includes:

[0081]

[0082] Among them, q material_1 Expressed as the heat flux density of the first material in W / m 2 ; r material_2 Expressed as the radius of the second material, in m; T Omaterial_2 Expressed as the temperature of the second material in °C.

[0083] Here, the third material boundary condition includes:

[0084]

[0085] Among them, r xinkuai Expressed as the radius of the core block, in m; q xinkuai Expressed as the heat flux density of the core in W / m 2 ; r material_3 Expressed as the radius of the third material, in m; T Omaterial_3 Expressed as the temperature of the third material in °C.

[0086] Here, the core boundary conditions include:

[0087]

[0088] Among them, q material_2 Expressed as the heat flux density of the second material in W / m 2 ; T Oxinkuai It is expressed as the temperature of the core block in °C.

[0089] Here, the shell boundary conditions include:

[0090]

[0091] Among them, q qixi Expressed as the heat flux density of the air gap in W / m2 ; r baoqiao It is expressed as the radius of the cladding, in m; T w Expressed as the outer wall temperature of the absorber in °C.

[0092] Exemplarily, the model boundary conditions of the constituent structure also include air gap boundary conditions, which can be expressed as:

[0093]

[0094] Among them, q material_3 Expressed as the heat flux density of the third material in W / m 2 ; r qixi Expressed as the radius of the air gap, in m; T Oqixi Expressed as air gap temperature in °C.

[0095] In some other embodiments, the composition structure of the absorbent material may only include material 1, a shell, a core block or other composition forms, which are not specifically limited here.

[0096] S103, acquiring parameter information of the composition structure of the absorber material, and calculating the temperature field information of the absorber based on the one-dimensional cylindrical radial heat conduction model according to the parameter information and the temperature of the coolant in the guide tube.

[0097] It can be understood that after obtaining the one-dimensional cylindrical radial heat conduction model and the model boundary conditions about the component structure, the parameter information of the component structure of the absorber material is obtained, which may include the radius, thermal conductivity, volume heat release rate, etc. corresponding to each component structure. Then, the temperature field information of the absorber is calculated based on the one-dimensional cylindrical radial heat conduction model according to the parameter information and the absorber outer wall temperature calculated according to the guide tube coolant temperature.

[0098] Exemplarily, before calculating the temperature field information of the absorber based on the one-dimensional cylindrical radial heat conduction model, it is also necessary to determine the temperature field information calculation formula based on the one-dimensional cylindrical radial heat conduction model and the model boundary conditions about the component structure. The temperature field information calculation formula includes the first material temperature information calculation formula, the second material temperature information calculation formula, the third material temperature information calculation formula, the core block temperature information calculation formula, the cladding temperature information calculation formula and the air gap temperature information calculation formula.

[0099] Specifically, the first material temperature information calculation formula includes:

[0100]

[0101] Among them, T m1 Indicates the temperature distribution of the first material, in °C; Φ material_1Indicates the volume heat release rate of the first material, in kW / m 3 ;k material_1 It represents the thermal conductivity of the first material, in W / M / °C.

[0102] And, the second material temperature information calculation formula includes:

[0103]

[0104] Among them, T m2 Indicates the temperature distribution of the second material, in °C; Φ material_2 Indicates the volume heat release rate of the second material, in kW / m 3 ;k material_2 It represents the thermal conductivity of the second material, in W / M / °C.

[0105] And, the third material temperature information calculation formula includes:

[0106]

[0107] Among them, T m3 Indicates the temperature distribution of the third material, in °C; Φ material_3 Indicates the volume heat release rate of the third material, in kW / m 3 ;k material_3 Indicates the thermal conductivity of the third material, in W / M / °C.

[0108] And, the core block temperature information calculation formula includes:

[0109]

[0110] Among them, T xin Indicates the core block temperature distribution, in °C; Φ xinkuai Indicates the volume heat release rate of the core block, in kW / m 3 ;k xinkuai Indicates the thermal conductivity of the third material, in W / M / °C.

[0111] And, the calculation formula of cladding temperature information includes:

[0112]

[0113] Among them, T bao Indicates the cladding temperature distribution in °C; Φ baoqiao Indicates the volume heat release rate of the cladding, in kW / m 3 ;k baoqiao Indicates the thermal conductivity of the cladding, in W / M / ℃.

[0114] And, the air gap temperature information calculation formula includes:

[0115]

[0116] Among them, T qi Indicates the air gap temperature distribution, in °C; Φ qixi Indicates the volume heat release rate of the air gap, in kW / m 3 ;k qixi Indicates the thermal conductivity of the air gap, in W / M / ℃.

[0117] Exemplarily, the temperature field information of the absorber, namely the first material temperature information, the second material temperature information, the third material temperature information, the core block temperature information, the cladding temperature information and the air gap temperature information can be calculated according to the above temperature field information calculation formula.

[0118] S104, judging whether the temperature field information of the absorber meets the preset temperature limit requirement, and obtaining the cooling effect of the guide tube on the absorber based on the judgment result.

[0119] It is understandable that when judging the cooling effect of the cooling tube on the absorber, it is necessary to judge whether the first material temperature information, the second material temperature information, the third material temperature information, the core block temperature information, the cladding temperature information and the air gap temperature information meet the corresponding preset temperature limit requirements. Among them, the preset temperature limit requirements corresponding to each type of component structure are for each material or area, and a reasonable preset temperature limit is set according to its material properties, working environment and safety standards. If all temperature information are lower than the preset temperature limit requirements, it means that the cooling effect of the guide tube on the absorber is good; otherwise, it means that the cooling effect of the guide tube on the cooling body does not meet the cooling requirements, which may be due to insufficient cooling medium flow, excessive temperature, cooling system design defects or abnormal heat generation inside the absorber. At this time, an alarm should be triggered immediately, and emergency measures (such as increasing the cooling medium flow, reducing system power, etc.) may need to be taken to prevent potential safety risks.

[0120] The calculation method for three-dimensional temperature field analysis in the guide tube absorber of a pressurized water reactor provided in the embodiments of the present disclosure calculates the absorber temperature field information by combining a one-dimensional cylindrical radial heat conduction model with the composition and structural parameter information of the absorber material, and can more accurately and comprehensively obtain the temperature distribution inside the guide tube absorber, thereby providing data support for the thermal-hydraulic design of the absorber in the guide tube.

[0121] Those skilled in the art will appreciate that, in the above method of specific implementation, the order in which the steps are written does not imply a strict execution order and does not constitute any limitation on the implementation process. The specific execution order of the steps should be determined by their functions and possible internal logic.

[0122] Based on the same inventive concept, the presently disclosed embodiment also provides a calculation device for three-dimensional temperature field analysis inside the pressurized water reactor guide tube absorber corresponding to the calculation method for three-dimensional temperature field analysis inside the pressurized water reactor guide tube absorber. Since the principle of solving the problem by the device in the presently disclosed embodiment is similar to the calculation method for three-dimensional temperature field analysis inside the pressurized water reactor guide tube absorber mentioned above in the presently disclosed embodiment, the implementation of the device can refer to the implementation of the method, and the repeated parts will not be repeated.

[0123] Reference Figure 4 FIG. 4 is a schematic diagram of a computing device 400 for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor provided by an embodiment of the present disclosure, wherein the device comprises:

[0124] A coolant solution module 401 is used to calculate the coolant temperature in the guide tube according to a preset guide tube coolant calculation formula and guide tube information;

[0125] A model building module 402 is used to obtain the composition structure of the absorber material, and to build a one-dimensional cylindrical radial heat conduction model and model boundary conditions about the composition structure based on the composition structure of the absorber material;

[0126] The temperature field solving module 403 is used to obtain parameter information of the composition structure of the absorber material, and calculate the temperature field information of the absorber based on the one-dimensional cylindrical radial heat conduction model according to the parameter information and the temperature of the coolant in the guide tube;

[0127] The cooling effect determination module 404 is used to determine whether the temperature field information of the absorber meets the preset temperature limit requirement, and obtain the cooling effect of the guide tube on the absorber based on the determination result.

[0128] In some possible embodiments, the preset guide tube coolant calculation formula includes a coolant power calculation formula and a coolant enthalpy calculation formula, the coolant power calculation formula includes a guide tube wall power calculation formula, a control rod power calculation formula and a coolant power calculation formula; the guide tube information includes guide tube wall information, control rod information and coolant information; the coolant solution module 401 is specifically used to:

[0129] Calculate the guide tube wall power based on the guide tube wall power calculation formula and the guide tube wall information; calculate the control rod power based on the control rod calculation formula and the control rod information; and calculate the coolant power based on the coolant power calculation formula and the coolant information;

[0130] According to the guide tube wall power, control rod power and coolant power, the enthalpy value of the guide tube segment is calculated based on the coolant enthalpy value calculation formula, and the coolant temperature is determined according to the enthalpy value;

[0131] The guide tube wall power calculation formula includes:

[0132] Q 1,i =q×Ph×z;

[0133] Among them, Q 1,i It is represented by the guide tube wall power of the i-th block; q is represented by the heat flux density of the guide tube wall, in W / m 2 ; Ph represents the heating circumference of the inner wall of the guide tube, in m; z represents the calculated length of the node, in m;

[0134] The control rod calculation formula includes:

[0135]

[0136] Among them, Q 2,i Expressed as the control rod power of the i-th node; Φ abs Indicates the volume heat release rate of the control rod absorption material, in kW / m 3 ;D abs Indicates the diameter of the control rod absorber material, in m; Φ clad Indicates the volume heat release rate of the control rod cladding metal, in kW / m 3 ;D R represents the outer diameter of the control rod, in m;

[0137] The coolant power calculation formula includes:

[0138]

[0139] Among them, Q 3,i It is expressed as the coolant power of the i-th block; D i It is expressed as the inner diameter of the guide tube of the i-th block, in m; D R Expressed as the outer diameter of the absorber, in m; Φ wat Expressed as the volume heat release rate of the coolant, in kW / m 3 ; NABS = 1 means that an absorber is inserted; NABS = 0 means that no absorber is inserted;

[0140] The coolant enthalpy calculation formula includes:

[0141]

[0142] Q i =Q 1,i +Q 2,i +Q 3,i ;

[0143] Among them, H iIt is expressed as the outlet enthalpy of the ith node, in kJ / kg, H i-1 It is expressed as the outlet enthalpy value of the previous node of the i-th node, in kJ / kg, with the initial value H0 being 0; Q i represents the total power absorbed by the coolant of the i-th node; W i It is expressed as the flow rate of the i-th block in the guide pipe, in kg / s.

[0144] In some possible embodiments, the coolant solution module 401 is further used to:

[0145] Calculating the absorber outer wall temperature according to the coolant temperature based on the absorber outer wall temperature solution formula;

[0146] The absorber outer wall temperature solution formula includes:

[0147]

[0148] Among them, T w Expressed as the outer wall temperature of the absorber, in °C; T f It is expressed as the coolant temperature in °C; q is expressed as the heat flux density of the guide tube wall in W / m 2 ; h represents the absorber surface heat transfer coefficient, unit is m 2 ℃ / W.

[0149] In some possible embodiments, the one-dimensional cylindrical radial heat conduction model includes:

[0150]

[0151] Where r is the radial coordinate in m; k(T) is the thermal conductivity in W / M / ℃; Φ is the volume heat release rate in W / m 3 .

[0152] In some possible embodiments, the composition structure of the absorbent material includes a first material, a second material, a third material, a core block and a cladding; the model boundary conditions of the composition structure include a first material boundary condition, a second material boundary condition, a third material boundary condition, a core block boundary condition and a cladding boundary condition;

[0153] The first material boundary condition includes:

[0154]

[0155] Among them, r material_1 Indicates the radius of the first material, in m; T Omaterial_1 It is expressed as the temperature of the first material in °C;

[0156] The second material boundary condition includes:

[0157]

[0158] Among them, q material_1 Expressed as the heat flux density of the first material in W / m 2 ; r material_2 Expressed as the radius of the second material, in m; T Omaterial_2 It is expressed as the temperature of the second material in °C;

[0159] The third material boundary condition includes:

[0160]

[0161] Among them, r xinkuai Expressed as the radius of the core block, in m; q xinkuai Expressed as the heat flux density of the core in W / m 2 ; r material_3 Expressed as the radius of the third material, in m; T Omaterial_3 It is expressed as the temperature of the third material in °C;

[0162] The core block boundary conditions include:

[0163]

[0164] Among them, q material_2 Expressed as the heat flux density of the second material in W / m 2 ; T Oxinkuai It is expressed as the temperature of the core block in °C;

[0165] The cladding boundary conditions include:

[0166]

[0167] Among them, q qixi Expressed as the heat flux density of the air gap in W / m 2 ; r baoqiao It is expressed as the radius of the cladding, in m; T w Expressed as the outer wall temperature of the absorber in °C.

[0168] In some possible embodiments, the model boundary conditions of the constituent structure further include air gap boundary conditions.

[0169] The air gap boundary conditions include:

[0170]

[0171] Among them, q material_3 Expressed as the heat flux density of the third material in W / m2 ; r qixi Expressed as the radius of the air gap, in m; T Oqixi Expressed as air gap temperature in °C.

[0172] In some possible embodiments, the temperature field solving module 403 is further used for:

[0173] Determine a temperature field information calculation formula based on the one-dimensional cylindrical radial heat conduction model and the model boundary conditions about the component structure; the temperature field information calculation formula includes a first material temperature information calculation formula, a second material temperature information calculation formula, a third material temperature information calculation formula, a core block temperature information calculation formula, a cladding temperature information calculation formula and an air gap temperature information calculation formula;

[0174] The first material temperature information calculation formula includes:

[0175]

[0176] Among them, T m1 Indicates the temperature distribution of the first material, in °C; Φ material_1 Indicates the volume heat release rate of the first material, in kW / m 3 ;k material_1 represents the thermal conductivity of the first material, in W / M / °C;

[0177] The second material temperature information calculation formula includes:

[0178]

[0179] Among them, T m2 Indicates the temperature distribution of the second material, in °C; Φ material_2 Indicates the volume heat release rate of the second material, in kW / m 3 ;k material_2 Indicates the thermal conductivity of the second material, in W / M / ℃;

[0180] The third material temperature information calculation formula includes:

[0181]

[0182] Among them, T m3 Indicates the temperature distribution of the third material, in °C; Φ material_3 Indicates the volume heat release rate of the third material, in kW / m 3 ;k material_3 Indicates the thermal conductivity of the third material, in W / M / ℃;

[0183] The core block temperature information calculation formula includes:

[0184]

[0185] Among them, T xin Indicates the core block temperature distribution, in °C; Φ xinkuai Indicates the volume heat release rate of the core block, in kW / m 3 ;k xinkuai Indicates the thermal conductivity of the third material, in W / M / ℃;

[0186] The cladding temperature information calculation formula includes:

[0187]

[0188] Among them, T bao Indicates the cladding temperature distribution in °C; Φ baoqiao Indicates the volume heat release rate of the cladding, in kW / m 3 ;k baoqiao Indicates the thermal conductivity of the cladding, in W / M / ℃;

[0189] The air gap temperature information calculation formula includes:

[0190]

[0191] Among them, T qi Indicates the air gap temperature distribution, in °C; Φ qixi Indicates the volume heat release rate of the air gap, in kW / m 3 ;k qixi Indicates the thermal conductivity of the air gap, in W / M / ℃.

[0192] Based on the same technical concept, the embodiment of the present disclosure also provides a computer device. Figure 5 , which is a schematic diagram of the structure of a computer device 500 provided in an embodiment of the present disclosure, including a processor 501, a memory 502, and a bus 503. The memory 502 is used to store execution instructions, including a memory 5021 and an external memory 5022; the memory 5021 is also called an internal memory, which is used to temporarily store the operation data in the processor 501 and the data exchanged with the external memory 5022 such as a hard disk. The processor 501 exchanges data with the external memory 5022 through the memory 5021.

[0193] In the embodiment of the present application, the memory 502 is specifically used to store the application code for executing the solution of the present application, and the execution is controlled by the processor 501. That is, when the computer device 500 is running, the processor 501 communicates with the memory 502 through the bus 503, so that the processor 501 executes the application code stored in the memory 502, and then executes the method described in any of the above embodiments.

[0194] Among them, the memory 502 can be, but is not limited to, random access memory (Random Access Memory, RAM), read only memory (Read Only Memory, ROM), programmable read-only memory (Programmable Read-Only Memory, PROM), erasable programmable read-only memory (Erasable Programmable Read-Only Memory, EPROM), electrically erasable read-only memory (Electric Erasable Programmable Read-Only Memory, EEPROM), etc.

[0195] Processor 501 may be an integrated circuit chip with signal processing capabilities. The above-mentioned processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components. The disclosed methods, steps and logic block diagrams in the embodiments of the present invention may be implemented or executed. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0196] It is to be understood that the structure illustrated in the embodiment of the present application does not constitute a specific limitation on the computer device 500. In other embodiments of the present application, the computer device 500 may include more or fewer components than shown in the figure, or combine some components, or split some components, or arrange the components differently. The components shown in the figure may be implemented in hardware, software, or a combination of software and hardware.

[0197] The present disclosure also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the calculation method for three-dimensional temperature field analysis in the absorber of a guide tube of a pressurized water reactor described in the above method embodiment are executed. The storage medium may be a volatile or non-volatile computer-readable storage medium.

[0198] The present disclosure also provides a computer program product that carries a program code. The program code includes instructions that can be used to execute the steps of the calculation method for three-dimensional temperature field analysis in a pressurized water reactor guide tube absorber described in the above method embodiment. For details, please refer to the above method embodiment, which will not be repeated here.

[0199] The computer program product may be implemented in hardware, software or a combination thereof. In one optional embodiment, the computer program product is implemented as a computer storage medium. In another optional embodiment, the computer program product is implemented as a software product, such as a software development kit (SDK).

[0200] Those skilled in the art can clearly understand that, for the convenience and simplicity of description, the specific working process of the system and device described above can refer to the corresponding process in the aforementioned method embodiment, and will not be repeated here. In the several embodiments provided in the present disclosure, it should be understood that the disclosed system and method can be implemented in other ways. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some communication interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0201] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0202] In addition, each functional unit in each embodiment of the present disclosure may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0203] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a non-volatile computer-readable storage medium that is executable by a processor. Based on this understanding, the technical solution of the present disclosure, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present disclosure. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0204] Finally, it should be noted that the above-described embodiments are only specific implementation methods of the present disclosure, which are used to illustrate the technical solutions of the present disclosure, rather than to limit them. The protection scope of the present disclosure is not limited thereto. Although the present disclosure is described in detail with reference to the above-described embodiments, ordinary technicians in the field should understand that any technician familiar with the technical field can still modify the technical solutions recorded in the above-described embodiments within the technical scope disclosed in the present disclosure, or can easily think of changes, or make equivalent replacements for some of the technical features therein; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present disclosure, and should be included in the protection scope of the present disclosure. Therefore, the protection scope of the present disclosure should be based on the protection scope of the claims.

Claims

1. A calculation method for three-dimensional temperature field analysis in a pressurized water reactor guide tube absorber, characterized in that: include: Calculate the coolant temperature in the guide tube according to a preset guide tube coolant calculation formula and guide tube information; Acquire the composition structure of the absorber material, and establish a one-dimensional cylindrical radial heat conduction model and model boundary conditions about the composition structure based on the composition structure of the absorber material; Acquiring parameter information of the composition structure of the absorber material, and calculating the temperature field information of the absorber based on the one-dimensional cylindrical radial heat conduction model according to the parameter information and the temperature of the coolant in the guide tube; It is determined whether the temperature field information of the absorber meets the preset temperature limit requirement, and based on the determination result, the cooling effect of the guide tube on the absorber is obtained.

2. The method according to claim 1, characterized in that The preset guide tube coolant calculation formula includes a coolant power calculation formula and a coolant enthalpy calculation formula, the coolant power calculation formula includes a guide tube wall power calculation formula, a control rod power calculation formula and a coolant power calculation formula; the guide tube information includes guide tube wall information, control rod information and coolant information; The step of calculating the coolant temperature in the guide tube according to the preset guide tube coolant calculation formula and the guide tube information includes: Calculate the guide tube wall power based on the guide tube wall power calculation formula and the guide tube wall information; calculate the control rod power based on the control rod calculation formula and the control rod information; and calculate the coolant power based on the coolant power calculation formula and the coolant information; According to the guide tube wall power, control rod power and coolant power, the enthalpy value of the guide tube segment is calculated based on the coolant enthalpy value calculation formula, and the coolant temperature is determined according to the enthalpy value; The guide tube wall power calculation formula includes: Q 1,i =q×Ph×z Among them, Q 1,i It is represented by the guide tube wall power of the i-th block; q is represented by the heat flux density of the guide tube wall, in W / m 2 ; Ph represents the heating circumference of the inner wall of the guide tube, in m; z represents the calculated length of the node, in m; The control rod calculation formula includes: Among them, Q 2,i Expressed as the control rod power of the i-th node; Φ abs Indicates the volume heat release rate of the control rod absorption material, in kW / m 3 ;D abs Indicates the diameter of the control rod absorber material, in m; Φ clad Indicates the volume heat release rate of the control rod cladding metal, in kW / m 3 ;D R represents the outer diameter of the control rod, in m; The coolant power calculation formula includes: Among them, Q 3,i It is expressed as the coolant power of the i-th block; D i It is expressed as the inner diameter of the guide tube of the i-th block, in m; D R Expressed as the outer diameter of the absorber, in m; Φ wat Expressed as the volume heat release rate of the coolant, in kW / m 3 ; NABS = 1 means that an absorber is inserted; NABS = 0 means that no absorber is inserted; The coolant enthalpy calculation formula includes: Q i =Q 1,i +Q 2,i +Q 3,i Among them, H i It is expressed as the outlet enthalpy of the ith node, in kJ / kg, H i-1 It is expressed as the outlet enthalpy value of the previous node of the i-th node, in kJ / kg, with the initial value H0 being 0; Q i represents the total power absorbed by the coolant of the i-th node; W i It is expressed as the flow rate of the i-th block in the guide pipe, in kg / s.

3. The method according to claim 2, characterized in that After determining the coolant temperature according to the enthalpy value, the method further comprises: Calculating the absorber outer wall temperature according to the coolant temperature based on the absorber outer wall temperature solution formula; The absorber outer wall temperature solution formula includes: Among them, T w Expressed as the outer wall temperature of the absorber, in °C; T f It is expressed as the coolant temperature in °C; q is expressed as the heat flux density of the guide tube wall in W / m 2 ; h represents the absorber surface heat transfer coefficient, unit is m 2 ℃ / W.

4. The method according to claim 3, characterized in that The one-dimensional cylindrical radial heat conduction model includes: Where r is the radial coordinate in m; k(T) is the thermal conductivity in W / M / ℃; Φ is the volume heat release rate in W / m 3 .

5. The method according to claim 4, characterized in that The composition structure of the absorber material includes a first material, a second material, a third material, a core block and a cladding; the model boundary conditions of the composition structure include a first material boundary condition, a second material boundary condition, a third material boundary condition, a core block boundary condition and a cladding boundary condition; The first material boundary condition includes: Among them, r material_1 Indicates the radius of the first material, in m; T Omaterial_1 It is expressed as the temperature of the first material in °C; The second material boundary condition includes: Among them, q material_1 Expressed as the heat flux density of the first material in W / m 2 ; r material_2 Expressed as the radius of the second material, in m; T Omaterial_2 It is expressed as the temperature of the second material in °C; The third material boundary condition includes: Among them, r xinkuai Expressed as the radius of the core block, in m; q xinkuai Expressed as the heat flux density of the core in W / m 2 ; r material_3 Expressed as the radius of the third material, in m; T Omaterial_3 It is expressed as the temperature of the third material in °C; The core block boundary conditions include: Among them, q material_2 Expressed as the heat flux density of the second material in W / m 2 ; T Oxinkuai It is expressed as the temperature of the core block in °C; The cladding boundary conditions include: Among them, q qixi Expressed as the heat flux density of the air gap in W / m 2 ; r baoqiao It is expressed as the radius of the cladding, in m; T w Expressed as the outer wall temperature of the absorber in °C.

6. The method according to claim 5, characterized in that The model boundary conditions of the constituent structure also include air gap boundary conditions, The air gap boundary conditions include: Among them, q material_3 Expressed as the heat flux density of the third material in W / m 2 ; r qixi Expressed as the radius of the air gap, in m; T Oqixi Expressed as air gap temperature in °C.

7. The method according to claim 6, characterized in that Before calculating the temperature field information of the absorber based on the one-dimensional cylindrical radial heat conduction model, the method includes: Determine a temperature field information calculation formula based on the one-dimensional cylindrical radial heat conduction model and the model boundary conditions about the component structure; the temperature field information calculation formula includes a first material temperature information calculation formula, a second material temperature information calculation formula, a third material temperature information calculation formula, a core block temperature information calculation formula, a cladding temperature information calculation formula and an air gap temperature information calculation formula; The first material temperature information calculation formula includes: Among them, T m1 Indicates the temperature distribution of the first material, in °C; Φ material_1 Indicates the volume heat release rate of the first material, in kW / m 3 ;k material_1 represents the thermal conductivity of the first material, in W / M / ℃; The second material temperature information calculation formula includes: Among them, T m2 Indicates the temperature distribution of the second material, in °C; Φ material_2 Indicates the volume heat release rate of the second material, in kW / m 3 ;k material_2 Indicates the thermal conductivity of the second material, in W / M / ℃; The third material temperature information calculation formula includes: Among them, T m3 Indicates the temperature distribution of the third material, in °C; Φ material_3 Indicates the volume heat release rate of the third material, in kW / m 3 ;k material_3 Indicates the thermal conductivity of the third material, in W / M / ℃; The core block temperature information calculation formula includes: Among them, T xin Indicates the core block temperature distribution, in °C; Φ xinkuai Indicates the volume heat release rate of the core block, in kW / m 3 ;k xinkuai Indicates the thermal conductivity of the third material, in W / M / ℃; The cladding temperature information calculation formula includes: Among them, T bao Indicates the cladding temperature distribution in °C; Φ baoqiao Indicates the volume heat release rate of the cladding, in kW / m 3 ;k baoqiao Indicates the thermal conductivity of the cladding, in W / M / ℃; The air gap temperature information calculation formula includes: Among them, T qi Indicates the air gap temperature distribution, in °C; Φ qixi Indicates the volume heat release rate of the air gap, in kW / m 3 ;k qixi Indicates the thermal conductivity of the air gap, in W / M / ℃.

8. A computing device for three-dimensional temperature field analysis in a guide tube absorber of a pressurized water reactor, characterized in that: include: A coolant solution module, used to calculate the coolant temperature in the guide tube according to a preset guide tube coolant calculation formula and guide tube information; A model building module, used for obtaining the composition structure of the absorber material, and building a one-dimensional cylindrical radial heat conduction model and model boundary conditions about the composition structure based on the composition structure of the absorber material; A temperature field solving module, used for obtaining parameter information of the composition structure of the absorber material, and calculating the temperature field information of the absorber based on the one-dimensional cylindrical radial heat conduction model according to the parameter information and the temperature of the coolant in the guide tube; The cooling effect determination module is used to determine whether the temperature field information of the absorber meets the preset temperature limit requirements, and obtain the cooling effect of the guide tube on the absorber based on the judgment result.

9. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

10. A computer device comprising a storage medium, a processor, and a computer program stored in the storage medium and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

Images