Temperature field analysis method and system for natural air cooling of control rod drive mechanism

Abstract

The present invention provides a temperature field analysis method and system for natural air cooling of a control rod drive mechanism. The analysis method comprises: establishing a control rod drive mechanism model, and setting different external air temperatures to perform simulation to obtain corresponding surface heat fluxes and first surface temperatures of the control rod drive mechanism model under different thermal loads; establishing a natural air cooling model, and performing simulation on the basis of the surface heat fluxes to obtain second surface temperatures of the control rod drive mechanism model; and on the basis of an interpolation method, determining a balance point temperature between the first surface temperatures and the second surface temperatures, wherein the control rod drive mechanism model comprises an internal flow channel, and the natural air cooling model is an integrated head package natural air cooling model. On the basis of the principles of sensitivity analysis, thermosyphon of a single drive mechanism model and natural air cooling of a full-field model are subjected to two-way coupling analysis, and finally, a drive mechanism surface heat flux and temperature field calculation result is obtained.

Description

A natural air cooling temperature field analysis method and system for a control rod drive mechanism

[0001] CROSS-REFERENCE TO RELATED APPLICATIONS

[0002] This application claims priority to Chinese patent application No. 202310969953.8, filed on August 2, 2023, entitled “A method and system for analyzing the temperature field of a control rod drive mechanism with natural air cooling,” and the entire contents of that application are incorporated herein by reference. Technical Field

[0003] The present invention relates to the technical field of temperature field analysis, and in particular to a natural air cooling temperature field analysis method and system for a control rod drive mechanism. Background Art

[0004] The statements in this section merely mention background art related to the present invention and do not necessarily constitute prior art.

[0005] The control rod drive mechanism (CRDM) is a step-by-step lifting mechanism in nuclear power reactors. By energizing the yoke coils in a certain sequence, it can lift, hold, or insert the control rod assembly, thereby changing the depth of the control rod assembly in the reactor core. This is used to implement the reactor's startup, shutdown, power regulation, load tracking, and emergency shutdown protection under accident conditions. It is a critical device for controlling the safe operation of the reactor.

[0006] The control rod drive mechanism (CRDM) needs to dissipate a large amount of heat during normal operation. This heat primarily comes from the heating of the magnetic yoke coils and the thermosiphoning effect of the coolant inside the core. To prevent the CRDM's magnetic yoke coils from overheating and ensure normal operation, the CRDM needs to be cooled. Cooling generally uses natural air cooling or forced air cooling. Forced air cooling uses an induced draft device, using a fan to force air over the surface of the drive mechanism to remove heat. This method has a strong cooling capacity and is widely used. Natural air cooling uses natural convection to remove heat and cool the CRDM. This cooling method does not require a fan, has a simple structure, and reduces energy loss, thereby improving economic efficiency. Because natural convection has a relatively weak heat transfer capacity, current reactors generally use low-power drive mechanisms when using natural air cooling.

[0007] To ensure that components such as the CRDM's yoke coils do not exceed temperature limits under natural ventilation conditions, a temperature field analysis based on fluid dynamics (CFD) is performed during the CRDM design process to determine ideal flow path layout, temperature settings, and other structural parameters. However, the CRDM's internal and external surfaces present complex flow field environments and boundary conditions. Its internal flow channels experience reactor coolant thermosiphoning, while its exterior is cooled by natural ventilation and thermal radiation. Furthermore, the CRDM itself is subject to a series of complex physical phenomena, including heat generation and heat conduction. Furthermore, the CRDM has numerous components and a complex structure. Conventional CFD simulations are affected by these numerous factors, resulting in extremely complex computational models and a massive mesh size, making it difficult to strike a balance between computational resources and accuracy.

[0008] Summary of the Invention

[0009] In view of this, the present invention provides a control rod drive mechanism natural air cooling temperature field analysis method and system. Based on the sensitivity analysis principle, a bidirectional coupling analysis process of thermosiphon of a single model of the control rod drive mechanism and natural air cooling of the full-field model is adopted to finally obtain the surface heat flux density and temperature field calculation results of the control rod drive mechanism.

[0010] In order to solve the above problems, the present invention provides the following technical solutions:

[0011] In a first aspect, the present invention provides a method for analyzing a temperature field of a control rod drive mechanism with natural air cooling, comprising the following steps:

[0012] Step S1: establishing a control rod drive mechanism model, performing simulations by setting different external air temperatures, and obtaining surface heat flux density and first surface temperature corresponding to the control rod drive mechanism model under different thermal loads;

[0013] Step S2: establishing a natural air cooling model, performing simulation based on the surface heat flux density, and obtaining a second surface temperature of the control rod drive mechanism model;

[0014] Step S3: determining a balance point temperature between the first surface temperature and the second surface temperature based on an interpolation method;

[0015] in,

[0016] The control rod drive mechanism model includes an internal flow channel;

[0017] The natural air cooling model is an integrated stack top assembly natural air cooling model.

[0018] Preferably, the control rod drive mechanism model is a two-dimensional or three-dimensional single control rod drive mechanism model.

[0019] Preferably, in step S1, the following physical parameters and boundary conditions are set:

[0020] Setting the heating value of the yoke coil in the control rod drive mechanism model;

[0021] Set the outside air temperature;

[0022] Setting a heat transfer coefficient of the surface of the control rod drive mechanism model;

[0023] Setting the metal thermal conductivity of the control rod drive mechanism model;

[0024] Setting the bottom temperature of the control rod drive mechanism model to the reactor upper head coolant temperature;

[0025] The physical properties of the coolant at different temperatures are set to simulate the thermal siphon caused by the difference in hot and cold densities.

[0026] Preferably, the natural air cooling model is a full-field natural air cooling three-dimensional flow channel model.

[0027] Preferably, in step S2, the heat flux density of the outer surface of the control rod drive mechanism model is set in the natural air cooling model according to the surface heat flux density.

[0028] Preferably, it also includes

[0029] Set the temperature of the insulation layer on the reactor pressure vessel head;

[0030] Set the reactor pressure vessel boundary air temperature to the containment ambient temperature;

[0031] Set the air property function that changes with temperature to simulate the natural convection parameters of air caused by temperature difference;

[0032] The thermal radiation emissivity of the control rod drive mechanism model, the upper head insulation layer of the reactor pressure vessel, and the shroud surface is set.

[0033] Preferably, in step S2,

[0034] According to the simulation results of the natural air cooling model, the flow field and temperature field of the natural air cooling model are determined to obtain the second surface temperature of the control rod drive mechanism model.

[0035] In a second aspect, the present invention provides a control rod drive mechanism natural air cooling temperature field analysis system, comprising:

[0036] a control rod drive mechanism thermal load module configured to: establish a control rod drive mechanism model, perform simulations by setting different ambient air temperatures, and obtain surface heat flux density and first surface temperature corresponding to the control rod drive mechanism model under different thermal loads;

[0037] A natural air cooling module is configured to: establish a natural air cooling model, perform simulation according to the surface heat flux density, and obtain a second surface temperature of the control rod drive mechanism model;

[0038] a coupling module configured to: determine a balance point temperature between the first surface temperature and the second surface temperature based on an interpolation method;

[0039] in,

[0040] The control rod drive mechanism model includes an internal flow channel;

[0041] The natural air cooling model is an integrated stack top assembly natural air cooling model.

[0042] In a third aspect, the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the above-mentioned control rod drive mechanism natural air cooling temperature field analysis method are completed.

[0043] In a fourth aspect, the present invention provides a computer device comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of the above-mentioned control rod drive mechanism natural air cooling temperature field analysis method are implemented.

[0044] Compared with the prior art, the present invention has the following beneficial effects:

[0045] The method and system for analyzing the natural air cooling temperature field of a control rod drive mechanism provided by the present invention simulate and calculate the thermosiphon model representing the internal flow channel of the control rod drive mechanism and the natural air cooling model representing the external flow channel separately, thereby reducing the computing resource usage and avoiding the computing resources occupied by traditional fluid mechanics analysis.

[0046] This method simulates the CRDM model at different external air temperatures to obtain the surface heat flux and first surface temperature of the CRDM model under different thermal loads. Using the surface heat flux as a variable, a natural air cooling simulation is performed to obtain the second surface temperature of the CRDM model. The simulation results of the CRDM model and the natural air cooling model are then coupled, using the surface heat flux and second surface temperature as variables, to determine the equilibrium point temperature between the first and second surface temperatures, representing the desired equilibrium state. This method provides more accurate results than traditional fluid dynamics analysis methods that rely on numerous assumptions.

[0047] The present invention calculates the inner and outer models of the control rod drive mechanism independently, thus getting rid of the constraints of computing resources and making full-field simulation of the integrated stack top assembly possible.

[0048] It should be understood that the foregoing general description and the following detailed description are merely illustrative and are not restrictive of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following briefly introduces the drawings required for use in the embodiments of the present application.

[0050] FIG1 is a flow chart of a method for analyzing the natural air cooling temperature field of a control rod drive mechanism according to an embodiment of the present invention;

[0051] FIG2 is a schematic diagram of the sensitivity analysis principle according to an embodiment of the present invention;

[0052] FIG3 is a schematic diagram of thermal load at a balance point according to an embodiment of the present invention;

[0053] FIG4 is a schematic diagram of a coupled analysis of a single control rod drive mechanism model and a full-field model according to an embodiment of the present invention;

[0054] FIG5 is a schematic diagram of the surface heat flux density of a control rod drive mechanism model according to an embodiment of the present invention;

[0055] FIG6 is a schematic diagram of a calculation model according to an embodiment of the present invention;

[0056] FIG7 is a schematic diagram of a natural air cooling flow field in an embodiment of the present invention;

[0057] FIG8 is a structural block diagram of a natural air cooling temperature field analysis system for a control rod drive mechanism according to an embodiment of the present invention.

[0058] Reference numerals: 1—CRDM heat load module; 2—natural air cooling module; 3—coupling module; a—CRDM heating analysis curve; b—natural air cooling analysis curve. DETAILED DESCRIPTION

[0059] In order to better understand the technical solution of the present application, the embodiments of the present application are described in detail below with reference to the accompanying drawings.

[0060] It should be noted that the following detailed descriptions are exemplary and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the art to which the present invention belongs.

[0061] The control rod drive mechanism (CRDM) is a step-by-step lifting mechanism in nuclear power reactors. By energizing the yoke coils in a certain sequence, it can lift, hold, or insert the control rod assembly, thereby changing the depth of the control rod assembly in the reactor core. This is used to implement the reactor's startup, shutdown, power regulation, load tracking, and emergency shutdown protection under accident conditions. It is a critical device for controlling the safe operation of the reactor.

[0062] Computational fluid dynamics (CFD) refers to the use of numerical solutions to the differential equations that control fluid flow to derive the discrete distribution of the fluid flow field over a continuous area, thereby approximately simulating the fluid flow conditions and determining the affected design parameters such as flow channels and temperature based on the fluid flow conditions.

[0063] To simplify reactor refueling operations, shorten refueling time, and reduce operator exposure, the integrated pile top assembly (IHP) connects the cables and bridges on the reactor pressure vessel top cover, the core measurement instrument lead-out structure, the control rod drive mechanism cooling fan, etc. with missile shielding and seismic support structures to form an integral structure that does not need to be disassembled during lifting.

[0064] The following example provides a method for analyzing the temperature field of a control rod drive mechanism (CRDM) with natural air cooling. Based on the principle of sensitivity analysis, a bidirectional coupling analysis process of thermosiphoning for a single CRDM model and natural air cooling for a full-field model is employed to ultimately calculate the surface heat flux density and temperature field of the CRDM.

[0065] FIG1 is a flow chart of a method for analyzing the natural air cooling temperature field of a control rod drive mechanism according to an embodiment of the present invention, FIG2 is a schematic diagram of the sensitivity analysis principle according to an embodiment of the present invention, and FIG3 is a schematic diagram of the equilibrium point heat load according to an embodiment of the present invention.

[0066] This embodiment provides a method for analyzing the temperature field of a control rod drive mechanism with natural air cooling, as shown in FIG1 , including the following steps:

[0067] Step S1: establishing a control rod drive mechanism model, performing simulations by setting different external air temperatures, and obtaining the surface heat flux density and first surface temperature corresponding to the control rod drive mechanism model under different thermal loads;

[0068] Step S2: Establishing a natural air cooling model, performing simulation based on the surface heat flux density, and obtaining a second surface temperature of the control rod drive mechanism model;

[0069] Step S3: determining a balance point temperature between the first surface temperature and the second surface temperature based on an interpolation method;

[0070] in,

[0071] The control rod drive mechanism model includes internal flow channels;

[0072] The natural air cooling model is an integrated stack top assembly natural air cooling model.

[0073] This example uses conventional 3D modeling software to create solid models. A single CRDM model, including the main flow channels, air gap, and solid regions of the CRDM body, and a naturally cooled model of the integrated reactor top assembly were constructed. The heat dissipation of the CRDM is divided into two components: heat generated by the yoke coils themselves and thermosiphoning of the reactor coolant within the CRDM's internal flow channels. As shown in Figure 2, the lower the CRDM surface temperature (first surface temperature), the greater the temperature difference between the cold and hot sources, the more pronounced the thermosiphon effect, and the greater the heat removed by air during natural cooling, i.e., the greater the surface heat flux density of the CRDM. Simultaneously, with the surface heat flux density as a variable, an increase in the CRDM thermal load corresponds to an increase in the CRDM surface temperature (second surface temperature). This indicates the CRDM thermal sensitivity: the lower the CRDM surface temperature, the greater the temperature difference between the cold and hot sources, and the greater the CRDM heat generation. The CRDM sensitivity to natural cooling: the greater the thermal load, the higher the CRDM surface temperature.

[0074] Based on the changes in the surface heat load and surface temperature of the CRDM, the CRDM heating analysis curve a and the natural air cooling analysis curve b are drawn as shown in Figure 3. The intersection of the two curves is the equilibrium point, that is, the equilibrium state result, which makes the surface temperature and heat load of the CRDM stable. Finding this equilibrium point can obtain the temperature result required for calculation.

[0075] The following describes each step in detail.

[0076] FIG4 is a schematic diagram of the coupled analysis of a single control rod drive mechanism model and a full-field model in an embodiment of the present invention, and FIG5 is a schematic diagram of the surface heat flux density of the control rod drive mechanism model in an embodiment of the present invention.

[0077] Step S1: establishing a control rod drive mechanism model, performing simulation by setting different external air temperatures, and obtaining the surface heat flux density and first surface temperature corresponding to the control rod drive mechanism model under different thermal loads.

[0078] The control rod drive mechanism model includes an internal flow channel, in which a reactor coolant thermal siphon exists. The control rod drive mechanism model is a two-dimensional or three-dimensional single control rod drive mechanism model.

[0079] In the above steps, as shown in Figure 4, the following physical parameters and boundary conditions are set:

[0080] Set the heating value of the yoke coil in the control rod drive mechanism model to a range of 2kW to 4kW;

[0081] Set the outside air temperature, as shown in Figure 2. The outside air temperature range is 100-300°C;

[0082] Set the heat transfer coefficient of the control rod drive mechanism model surface. Usually, a series of temperature values ​​are selected at intervals within the range of 50℃ to 200℃, and the output surface heat flux density is distributed with the change of height.

[0083] Set the metal thermal conductivity of the control rod drive mechanism model to a range of 16 to 24 (W / m·K).

[0084] Set the bottom temperature of the control rod drive mechanism model to the reactor upper head coolant temperature;

[0085] The physical properties of the coolant at different temperatures can be set using a linear interpolation method to simulate the thermal siphon caused by the difference in density between hot and cold temperatures.

[0086] The simulation was performed with the above settings. Based on the simulation results, the surface heat flux density and first surface temperature of the control rod drive mechanism model under different thermal loads were obtained. Figure 5 shows the surface heat flux density of the drive mechanism obtained when the outside air temperature was set to 100°C, ranging from 195.72W / m 2 ~2496.26W / m 2 .

[0087] FIG6 is a schematic diagram of a calculation model in an embodiment of the present invention, and FIG7 is a schematic diagram of a natural air cooling flow field in an embodiment of the present invention.

[0088] Step S2: Establish a natural air cooling model, perform simulation based on the surface heat flux density, and obtain the second surface temperature of the control rod drive mechanism model.

[0089] Among them, the natural air cooling model is a three-dimensional flow channel model of natural air cooling of the integrated top assembly.

[0090] The natural air cooling model mainly includes the external contours of the air domain inside the enclosure and the large external space, the insulation layer of the reactor pressure vessel upper head, the control rod drive mechanism model, and other structures.

[0091] In the above steps, as shown in Figure 4, the following physical parameters and boundary conditions are set:

[0092] According to the surface heat flux density of the control rod drive mechanism model as shown in FIG5 , the heat flux density of the outer surface of the control rod drive mechanism model is set in the natural air cooling model;

[0093] Set the temperature of the insulation layer on the reactor pressure vessel head to 25℃~100℃;

[0094] Set the reactor pressure vessel boundary air temperature to the containment ambient temperature and its temperature to the air temperature;

[0095] Set the air property function that changes with temperature. Linear interpolation method can be used to simulate the natural convection parameters of air caused by temperature difference.

[0096] The thermal radiation emissivity of the control rod drive mechanism model, the reactor pressure vessel upper head insulation layer, and the shroud surface is set in the range of 0.22 to 0.35.

[0097] The calculation model shown in Figure 6 is obtained through the above settings. The simulation is carried out. According to the simulation results, the flow field and temperature field of the natural air cooling model shown in Figure 7 are determined. Figure 7 is the surface heat flux density result (195.72W / m 2 ~2496.26W / m 2 ) as the input to obtain the natural air cooling flow field results, the flow rate range is: 0m / s ~ 4.38m / s, the second surface temperature of the control rod drive mechanism model is obtained, the second surface temperature range is 300℃ ~ 500℃.

[0098] Step S3: determining a balance point temperature between the first surface temperature and the second surface temperature based on an interpolation method.

[0099] As shown in Figure 3, based on the corresponding first surface temperatures of the CRAD model under different heat loads obtained after simulation in step S1, a CRAD heating analysis curve a is plotted. As the heat load increases within the range of 0 to 4 kW, the corresponding first surface temperature decreases within the range of 0 to 200°C. Based on the corresponding second surface temperatures under different heat loads obtained after simulation in step S2, a natural air cooling analysis curve b is plotted. As the heat load increases within the range of 0 to 4 kW, the corresponding second surface temperature increases within the range of 0 to 250°C. The two curves are compared to a single point, the equilibrium point, at which the two curves have the same surface temperature under the same heat load, resulting in the desired equilibrium state.

[0100] In the above process, the complete internal model of the control rod drive mechanism and the external flow channel model are calculated independently, thus breaking away from the constraints of computing resources and avoiding the problem of insufficient computing resources caused by traditional fluid mechanics analysis, making full-field simulation of the integrated stack top assembly possible.

[0101] By setting different external air temperatures for a single CRDM model, a series of thermal loads and first surface temperatures were obtained for the CRDM model. A sensitivity analysis using the thermal load as a variable was then performed to determine the second surface temperature of the CRDM. The results from the single CRDM model and the full-field model were coupled to obtain the final equilibrium state. This yields more accurate results than traditional computational fluid dynamics analysis methods that rely on numerous assumptions.

[0102] The above process provides a CFD analysis method for coupling the thermosiphon analysis of a single CRDM model with the natural air cooling of the full-field model to obtain the natural air cooling temperature field of the CRDM. In this process, the complete internal model of the CRDM (including metal structural parts and thermosiphon flow channels) and the external flow channel model (including the outer contour of the CRDM and the external air domain) are calculated separately to reduce computing resource usage.

[0103] FIG8 is a structural block diagram of a natural air cooling temperature field analysis system for a control rod drive mechanism according to an embodiment of the present invention.

[0104] As shown in FIG8 , this embodiment provides a control rod drive mechanism natural air cooling temperature field analysis system, including:

[0105] The control rod drive mechanism heat load module 1 is configured to: establish a control rod drive mechanism model, perform simulations by setting different ambient air temperatures, and obtain surface heat flux density and first surface temperature corresponding to the control rod drive mechanism model under different heat loads;

[0106] The natural air cooling module 2 is configured to: establish a natural air cooling model, perform simulation based on the surface heat flux density, and obtain a second surface temperature of the control rod drive mechanism model;

[0107] The coupling module 3 is configured to: determine a balance point temperature between the first surface temperature and the second surface temperature based on an interpolation method;

[0108] in,

[0109] The control rod drive mechanism model includes internal flow channels;

[0110] The natural air cooling model is an integrated stack top assembly natural air cooling model.

[0111] This embodiment provides a computer-readable storage medium, in which a computer program is stored. When the computer program is executed by a processor, the steps of the control rod drive mechanism natural air cooling temperature field analysis method are completed.

[0112] This embodiment provides a computer device including a memory and a processor. The memory stores a computer program. When the processor executes the computer program, the steps of the control rod drive mechanism natural air cooling temperature field analysis method are implemented.

[0113] The steps or modules involved in the above embodiments correspond to the control rod drive mechanism natural air cooling temperature field analysis method. For detailed implementation, please refer to the relevant description of the control rod drive mechanism natural air cooling temperature field analysis method. The term "computer-readable storage medium" should be understood to mean a single medium or multiple media containing one or more instruction sets; it should also be understood to include any medium capable of storing, encoding, or carrying an instruction set for execution by a processor and causing the processor to perform any method of the present invention.

[0114] The above description is merely a preferred embodiment of the present application and is not intended to limit the present application. Various modifications and variations are possible for those skilled in the art. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present application shall be included within the scope of protection of the present application.

Claims

1. A method for analyzing the temperature field of a control rod drive mechanism with natural air cooling, characterized in that: The following steps are involved: Step S1: establishing a control rod drive mechanism model, and performing simulation by setting different outside air temperatures to obtain the surface heat flux density and the first surface temperature corresponding to the control rod drive mechanism model under different heat loads; Step S2: establishing a natural air cooling model, performing simulation according to the surface heat flux density, and obtaining a second surface temperature of the control rod drive mechanism model; Step S3: determining a balance point temperature between the first surface temperature and the second surface temperature based on an interpolation method; in, The control rod drive mechanism model includes an internal flow channel; The natural air cooling model is an integrated stack top assembly natural air cooling model.

2. The natural air cooling temperature field analysis method for the control rod drive mechanism according to claim 1, characterized in that: The control rod drive mechanism model is a two-dimensional or three-dimensional single control rod drive mechanism model.

3. The natural air cooling temperature field analysis method for a control rod drive mechanism according to claim 1, characterized in that: In step S1, the following physical parameters and boundary conditions are set: Setting the heating value of the yoke coil in the control rod drive mechanism model; Set the outside air temperature; Setting a heat transfer coefficient of the surface of the control rod drive mechanism model; Setting the metal thermal conductivity of the control rod drive mechanism model; Setting the bottom temperature of the control rod drive mechanism model to the reactor upper head coolant temperature; The physical properties of the coolant at different temperatures are set to simulate the thermal siphon caused by the difference in cold and hot densities.

4. The natural air cooling temperature field analysis method for a control rod drive mechanism according to claim 1, characterized in that: The natural air cooling model is a full-field natural air cooling three-dimensional flow channel model.

5. The natural air cooling temperature field analysis method for the control rod drive mechanism according to claim 4, characterized in that: In step S2, According to the surface heat flux density, the heat flux density of the outer surface of the control rod drive mechanism model is set in the natural air cooling model.

6. The natural air cooling temperature field analysis method for the control rod drive mechanism according to claim 5, characterized in that: Also includes Set the temperature of the insulation layer on the reactor pressure vessel head; Set the reactor pressure vessel boundary air temperature to the containment ambient temperature; Set the air property function that changes with temperature to simulate the natural convection parameters of air caused by temperature difference; The thermal radiation emissivity of the control rod drive mechanism model, the upper head insulation layer of the reactor pressure vessel and the shroud surface is set.

7. The natural air cooling temperature field analysis method for the control rod drive mechanism according to claim 6, characterized in that: In step S2, According to the simulation result of the natural air cooling model, the flow field and the temperature field of the natural air cooling model are determined to obtain the second surface temperature of the control rod drive mechanism model.

8. A natural air cooling temperature field analysis system for a control rod drive mechanism, characterized in that: include: The control rod drive mechanism heat load module is configured to: establish a control rod drive mechanism model, perform simulation by setting different outside air temperatures, and obtain the surface heat flux density and the first surface temperature corresponding to the control rod drive mechanism model under different heat loads; A natural air cooling module is configured to: establish a natural air cooling model, perform simulation according to the surface heat flux density, and obtain a second surface temperature of the control rod drive mechanism model; A coupling module is configured to: determine a balance point temperature between the first surface temperature and the second surface temperature based on an interpolation method; in, The control rod drive mechanism model includes an internal flow channel; The natural air cooling model is an integrated stack top assembly natural air cooling model.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by the processor, the steps of the control rod drive mechanism natural air cooling temperature field analysis method according to any one of claims 1 to 7 are completed.

10. A computer device, characterized in that: The method comprises a memory and a processor, wherein a computer program is stored in the memory, and when the processor executes the computer program, the steps of the natural air cooling temperature field analysis method of the control rod drive mechanism according to any one of claims 1 to 7 are implemented.

Images