Thermal stratification phenomenon simulation method and apparatus, device, medium, and program product

Abstract

The present application relates to the technical field of nuclear reactor thermal hydraulic simulation analysis, and provides a thermal stratification phenomenon simulation method and apparatus, a device, a medium, and a program product. The method comprises: on the basis of an interactive model computation result in a pre-configured program interface, controlling a one-dimensional system program and a three-dimensional computational fluid mechanics program to perform coupling computation, wherein the program interface is configured between the one-dimensional system program and the three-dimensional computational fluid mechanics program, and the model computation result comprises results obtained by performing simulation computation on a liquid metal reactor by the one-dimensional system program and the three-dimensional computational fluid mechanics program; and obtaining a real-time computation result of coupling computation of the one-dimensional system program and the three-dimensional computational fluid mechanics program, wherein the real-time computation result comprises a simulation result of a thermal stratification phenomenon of the liquid metal reactor in a transient operation process. The present application can reduce the computing resource demand and the time cost and improve the efficiency while ensuring that the thermal stratification phenomenon of the liquid metal reactor is accurately simulated.

Description

Simulation method, device, equipment, medium and program product of thermal stratification phenomenon TECHNICAL FIELD

[0001] The present application relates to the technical field of nuclear reactor thermal-hydraulic simulation analysis, in particular to a simulation method, device, equipment, medium and program product of thermal stratification phenomenon. BACKGROUND

[0002] At present, the research on thermal stratification phenomenon of liquid metal reactors is mainly carried out by experimental and numerical simulation methods. Among them, the experimental method is limited by economy and operability, and usually cannot simulate the whole reactor system. At the same time, it is also difficult to measure the data of liquid metal pool temperature field and flow field for the realization of experiment. In addition, the numerical simulation method mainly uses three-dimensional computational fluid dynamics (CFD) program to simulate the complex thermal-hydraulic phenomenon of liquid metal pool. The accuracy of three-dimensional CFD simulation depends on the grid division accuracy. When simulating the whole reactor primary circuit, fine grid division requires huge computing resources.

[0003] SUMMARY

[0004] The main purpose of the embodiments of the present application is to provide a simulation method, device, equipment, medium and program product of thermal stratification phenomenon, which aims to reduce the demand for computing resources and improve the computing efficiency under the premise of accurately simulating the thermal stratification phenomenon of liquid metal reactors.

[0005] To achieve the above purpose, the first aspect of the embodiments of the present application provides a simulation method of thermal stratification phenomenon, which comprises:

[0006] Based on the interaction of model calculation results in the preconfigured program interface, the one-dimensional system program and the three-dimensional computational fluid dynamics program are coupled for calculation; the program interface is configured between the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the model calculation results include the results obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively for simulating and calculating the liquid metal reactor;

[0007] Obtain the real-time calculation results of the one-dimensional system program and the three-dimensional computational fluid dynamics program for coupled calculation, and the real-time calculation results include the simulation results of the thermal stratification phenomenon of the liquid metal reactor in the transient operation process.

[0008] In some embodiments, the model calculation results include first model calculation results and second model calculation results; the method further comprises:

[0009] The first model calculation result is obtained by calculating the first one-dimensional model of the first part of the liquid metal reactor through the one-dimensional system program.

[0010] The second model calculation result is obtained by calculating a three-dimensional model of a second part of the liquid metal reactor through the three-dimensional computational fluid dynamics program.

[0011] The second part is connected with the first part at an inlet of the second part and an outlet of the first part, and the second part is connected with the first part at an outlet of the second part and an inlet of the first part.

[0012] In some embodiments, the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program based on the interaction of the model calculation results in the preconfigured program interface comprises:

[0013] The first model calculation result is transmitted to the three-dimensional computational fluid dynamics program through the program interface when the first model calculation result is obtained by the one-dimensional system program;

[0014] The first model calculation result is taken as an input of the three-dimensional model to calculate the three-dimensional model at a Kth time step through the three-dimensional computational fluid dynamics program; K is a positive integer greater than or equal to 1;

[0015] The second model calculation result is transmitted to the one-dimensional system program through the program interface when the second model calculation result is obtained by the three-dimensional computational fluid dynamics program;

[0016] The second model calculation result is taken as an input of the first one-dimensional model to calculate the first one-dimensional model at a K+1th time step through the one-dimensional system program.

[0017] In some embodiments, after the calculation of the three-dimensional model through the three-dimensional computational fluid dynamics program with the first model calculation result as an input of the three-dimensional model, the method further comprises:

[0018] If the second model calculation result does not converge when the second model calculation result is obtained by the three-dimensional computational fluid dynamics program, the input of the first one-dimensional model is adjusted, and the first model calculation result is obtained by recalculating the first one-dimensional model at a Kth time step through the one-dimensional system program.

[0019] In some embodiments, the transmission of the second model calculation result to the one-dimensional system program through the program interface comprises:

[0020] In a case where the second model calculation result converges, the second model calculation result is transmitted to the one-dimensional system program through the program interface.

[0021] In some embodiments, the second model calculation result includes: a threshold inlet-outlet mass flow rate and a threshold monitoring point location temperature at a connection between an outlet of the second portion and an inlet of the first portion;

[0022] The method further includes:

[0023] In a case where both the threshold inlet-outlet mass flow rate and the threshold monitoring point location temperature are balanced, it is determined that the second model calculation result converges;

[0024] In a case where any one of the threshold inlet-outlet mass flow rate and the threshold monitoring point location temperature is not balanced, it is determined that the second model calculation result does not converge.

[0025] In some embodiments, before the first model calculation result is obtained by calculating, through the one-dimensional system program, the first one-dimensional model of the liquid metal reactor first portion, the method further includes:

[0026] A third model calculation result is obtained by calculating, through the one-dimensional system program, a second one-dimensional model of the liquid metal reactor as a whole.

[0027] In some embodiments, the first model calculation result is obtained by calculating, through the one-dimensional system program, the first one-dimensional model of the liquid metal reactor first portion, including:

[0028] The first one-dimensional model is calculated by the one-dimensional system program with the third model calculation result as an input of the first one-dimensional model, so as to obtain the first model calculation result.

[0029] To achieve the above object, a second aspect of the embodiments of the present application proposes a simulation device for thermal stratification phenomenon, the device comprising:

[0030] A coupling calculation module is configured to control a one-dimensional system program and a three-dimensional computational fluid dynamics program to perform coupling calculation based on interaction of model calculation results in a pre-configured program interface; the program interface is configured between the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the model calculation results include results obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively from simulation calculation of a liquid metal reactor;

[0031] A result obtaining module is configured to obtain real-time calculation results of the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the real-time calculation results include simulation results of thermal stratification of the liquid metal reactor during transient operation.

[0032] To achieve the above object, a third aspect of embodiments of the present application provides a computer device, comprising a memory and a processor, the memory stores a computer program, and the processor implements each step of the method provided in the first aspect when executing the computer program.

[0033] To achieve the above object, a fourth aspect of embodiments of the present application provides a computer readable storage medium, which stores a computer program, and the computer program implements each step of the method provided in the first aspect when executed by a processor.

[0034] To achieve the above object, a fifth aspect of embodiments of the present application provides a computer program product, which comprises a computer program, and the computer program implements each step of the method provided in the first aspect when executed by a processor.

[0035] The simulation method, device, equipment, medium and program product of the thermal stratification phenomenon described above, through the program interface configured between the one-dimensional system program and the three-dimensional computational fluid dynamics program, interact the model calculation results of the liquid metal reactor by the one-dimensional system program and the three-dimensional computational fluid dynamics program in real time, so as to control the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program; in this way, the process of the thermal stratification phenomenon of the liquid metal reactor system under the transient response condition is simulated. In this way, by obtaining the real-time calculation results of the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program, the simulation results of the thermal stratification phenomenon of the liquid metal reactor during the transient operation process can be obtained.

[0036] That is, the simulation method, device, equipment, medium and program product of the thermal stratification phenomenon provided by the present application can realize real-time data exchange of one-dimensional model and three-dimensional model simulation data through the coupling calculation process between the one-dimensional system program and the three-dimensional computational fluid dynamics program through the program interface between the one-dimensional system program and the three-dimensional computational fluid dynamics program, so that the simulation of the thermal stratification phenomenon under the transient response condition of the pool type liquid metal reactor can be realized. Moreover, compared with the analysis mode of separately simulating the liquid metal reactor by using the one-dimensional system program or the three-dimensional computational fluid dynamics program, the present application comprehensively combines the advantages of the three-dimensional computational fluid dynamics program in accurately simulating the thermal stratification phenomenon and the advantages of the system program in quickly simulating the overall transient response of the reactor system, accurately simulates the thermal stratification phenomenon in the transient operation process of the reactor, can save computing resources, improve simulation accuracy and computing efficiency, and further support the performance matching of the in-vessel components of the liquid metal reactor and the development of the thermal-hydraulic performance of the whole reactor. BRIEF DESCRIPTION OF DRAWINGS

[0037] Fig. 1 is a step flow diagram of the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application in some embodiments;

[0038] Fig. 2 is a coupling architecture diagram related to the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application;

[0039] Fig. 3 is a timing diagram related to the model calculation of the liquid metal reactor in the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application;

[0040] Fig. 4 is a structure region division diagram of the liquid metal reactor related to the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application;

[0041] Fig. 5 is another timing diagram related to the model calculation of the liquid metal reactor in the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application;

[0042] Fig. 6 is a step flow diagram of the coupling calculation related to the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application;

[0043] Fig. 7 is a coupling calculation flow diagram related to the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application;

[0044] Fig. 8 is a coupling calculation loop example diagram related to the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application;

[0045] Fig. 9 is a temperature cloud diagram of the heating horizontal tube heating section in the example shown in Fig. 8;

[0046] Fig. 10 is a flow velocity transient curve of the heating horizontal tube heating section in the example shown in Fig. 8;

[0047] Fig. 11 is a schematic structural diagram of a simulation device for thermal delamination in some embodiments of the present application;

[0048] Fig. 12 is a schematic diagram of the internal structure of a computer device according to an embodiment of the present application. DETAILED DESCRIPTION

[0049] For the purpose of the present application, the technical solutions and advantages are more clearly and explicitly described below in combination with the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and not intended to limit the present application.

[0050] It should be noted that although the functional modules are divided in the device schematic diagram, and the logical order is shown in the flowchart, in some cases, the steps shown or described can be performed in a manner different from the module division in the device or the order in the flowchart. The terms "first", "second", etc. in the specification and claims and the above-described drawings are used to distinguish similar objects, and do not necessarily describe a specific order or sequence.

[0051] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the present application belongs. The terms used herein are only for the purpose of describing the embodiments of the present application and are not intended to limit the present application.

[0052] Before introducing the embodiments of the present application, the related technical terms involved in the embodiments of the present application are briefly introduced.

[0053] One-dimensional system program and three-dimensional computational fluid dynamics program.

[0054] Both the one-dimensional system program and the three-dimensional computational fluid dynamics program can be used for simulation analysis of nuclear reactors. The one-dimensional system program is also called a system-level analysis program, which can be mainly used for fast simulation of the overall reactor system. The one-dimensional system program can be written by the user based on the computer programming language Python or directly called. Currently, mature one-dimensional system programs available for direct calling on the market include LOCUST, TRACE, ATHLET, etc. In addition, the three-dimensional computational fluid dynamics program is mainly used for simulation analysis of complex thermal-hydraulic phenomena of the reactor. When the three-dimensional computational fluid dynamics program simulates the reactor, the more fine the grid division, the higher the accuracy of the simulation results, and the higher the cost of computing resources and time of the device. Currently, common three-dimensional computational fluid dynamics programs include STAR-CCM+, OpenFOAM, FLUENT, and CFX, etc.

[0055] Liquid metal reactor.

[0056] Liquid metal reactor, also known as liquid metal fuel reactor, is a reactor using liquid metal or liquid alloy containing fissile material as fuel. The liquid metal fuel has good heat conduction performance, no irradiation swelling problem, simple core structure, and can be operated at high temperature and low pressure, but the metal fuel must be heated to be in a molten state. The liquid metal reactor is composed of a reactor, a primary loop, an intermediate loop, a secondary loop and a propulsion shaft system.

[0057] Thermal stratification phenomenon of liquid metal reactor

[0058] Due to the strong heat conduction of the liquid metal coolant, after adding the liquid metal coolant, the liquid metal reactor is very easy to form thermal stratification phenomenon in the upper liquid metal pool of the reactor. Once the thermal stratification phenomenon occurs in the liquid metal reactor, it will cause the internal components of the reactor to form obvious thermal stress, and affect the coolant flow distribution and system transient response in the reactor, which is very unfavorable for the structure and safe operation of the reactor. Therefore, the research on the thermal stratification phenomenon of the liquid metal reactor is of great importance to the safe operation of the liquid metal reactor, the performance matching of the internal components and the development of the overall thermal hydraulic performance of the reactor.

[0059] Next, the overall concept of the simulation method of the thermal stratification phenomenon provided by the embodiment of the present application is further described.

[0060] At present, the research on the thermal stratification phenomenon of the liquid metal reactor is mainly carried out by experimental and numerical simulation methods. Among them, the experimental method is limited by economy and operability, and usually cannot simulate the whole reactor system, and it is also difficult to measure the data of the temperature field and flow field of the liquid metal pool for the realization of the experiment. In addition, the numerical simulation method mainly uses three-dimensional computational fluid dynamics (CFD) program to simulate the complex thermal hydraulic phenomenon of the liquid metal pool in detail, and the accuracy of the three-dimensional CFD simulation depends on the grid division accuracy. When simulating the whole primary loop of the reactor, fine grid division requires huge computing resources and time cost.

[0061] To this end, the embodiment of the present application provides a simulation method of thermal stratification. The coupling calculation process between the one-dimensional system program and the three-dimensional computational fluid dynamics program is performed through the program interface between the one-dimensional system program and the three-dimensional computational fluid dynamics program, the real-time data exchange of the one-dimensional model and the three-dimensional model simulation data can be realized, and thus the simulation of the thermal stratification under the transient response condition of the pool type liquid metal reactor can be realized. Compared with the analysis mode of simulating the liquid metal reactor by using the one-dimensional system program alone, the embodiment of the present application can establish a complete model of the liquid metal reactor, and thus the simulation accuracy can be improved. Compared with the analysis mode of simulating the liquid metal reactor by using the three-dimensional computational fluid dynamics program alone, the embodiment of the present application does not need to establish a three-dimensional model for each part of the liquid metal reactor, and thus the computing resources can be saved and the computing efficiency can be improved. That is, the embodiment of the present application combines the advantages of the three-dimensional computational fluid dynamics program in accurately simulating the thermal stratification and the advantages of the system program in quickly simulating the overall transient response of the reactor system, accurately simulates the thermal stratification in the transient operation process of the reactor, can save computing resources, improve simulation accuracy and computing efficiency, and further support the performance matching of the in-vessel components of the liquid metal reactor and the development of the thermal-hydraulic performance of the whole reactor.

[0062] It should be understood that the simulation method of thermal stratification provided by the embodiment of the present application can be applied to a terminal, can be applied to a server side, and can be software running in a terminal or a server side. In some embodiments, the terminal can be a smart phone, a tablet computer, a notebook computer, a desktop computer, etc.; the server side can be configured as a separate physical server, can be configured as a server cluster or a distributed system composed of multiple physical servers, can be configured as a cloud server providing basic cloud computing services such as cloud service, cloud database, cloud computing, cloud function, cloud storage, network service, cloud communication, middleware service, domain name service, security service, CDN, and big data and artificial intelligence platform, etc.; and the software can be an application that implements the simulation method of thermal stratification, but is not limited to the above forms.

[0063] Alternatively, the present application can also be used in a plurality of general or special computer system environments or configurations. For example: personal computers, server computers, handheld devices or portable devices, tablet devices, multi-processor systems, microprocessor-based systems, set-top boxes, programmable consumer computer devices, network PCs, minicomputers, mainframe computers, distributed computing environments including any of the above systems or devices, etc. The present application can be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types.

[0064] For the sake of understanding and description, the terminal device applies the simulation method of thermal stratification phenomenon provided by the embodiments of the present application as an example to describe the present application in detail. The implementation of the simulation method of thermal stratification phenomenon provided by the embodiments of the present application in any form of subject can refer to the process of the simulation method of thermal stratification phenomenon provided by the embodiments of the present application.

[0065] Please refer to FIG. 1, which is a step flowchart of the simulation method of thermal stratification phenomenon provided by the embodiments of the present application in some embodiments. It should be understood that although the execution order of some method steps is shown in FIG. 1, the simulation method of thermal stratification phenomenon provided by the embodiments of the present application can of course adopt an execution order of method steps different from that shown in the figure based on different design needs of actual application. That is, the order of the method steps shown in FIG. 1 does not constitute a limitation on the execution logic order of the simulation method of thermal stratification phenomenon provided by the embodiments of the present application, and any reasonable change based on the order of the method steps shown in FIG. 1 should be included in the protection scope of the simulation method of thermal stratification phenomenon provided by the embodiments of the present application.

[0066] As shown in FIG. 1, in some embodiments, the simulation method of thermal stratification phenomenon provided by the embodiments of the present application can include steps S101 to S102.

[0067] Step S101, based on the model calculation result in the pre-configured program interface, control the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program.

[0068] It should be noted that the program interface is configured between the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the model calculation result includes the results obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively for simulating and calculating the liquid metal reactor.

[0069] Step S102, obtain the real-time calculation result of the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the real-time calculation result includes the simulation result of the thermal stratification phenomenon of the liquid metal reactor in the transient operation process.

[0070] The terminal device can combine the one-dimensional system program and the three-dimensional computational fluid dynamics program to simulate the thermal stratification phenomenon of the liquid metal reactor in the transient operation process when the terminal device executes the simulation method for the thermal stratification phenomenon provided in the embodiments of the present application. For example, the terminal device can pre-establish a program interface between the one-dimensional system program and the three-dimensional computational fluid dynamics program. During the process of calculating the corresponding model of the liquid metal reactor by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively, the terminal device can interact the model calculation results obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively to calculate the model of the liquid metal reactor through the program interface, so as to perform coupled calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program. In this way, the terminal device can obtain the simulation result of simulating the thermal stratification phenomenon of the liquid metal reactor in the transient operation process by obtaining the real-time calculation result obtained by the coupled calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program.

[0071] In some embodiments, when the terminal device establishes the program interface between the one-dimensional system program and the three-dimensional computational fluid dynamics program, the terminal device can establish an extension interface of the one-dimensional system program as the program interface between the one-dimensional system program and the three-dimensional computational fluid dynamics program, or the terminal device can also establish an extension interface of the three-dimensional computational fluid dynamics program as the program interface between the one-dimensional system program and the three-dimensional computational fluid dynamics program.

[0072] For example, assuming that the three-dimensional computational fluid dynamics program is STAR-CCM+ software, the terminal device can establish a JAVA extension interface of the STAR-CCM+ software as the program interface between the STAR-CCM+ software and the one-dimensional system program. Thus, the model calculation result obtained by the terminal device through the one-dimensional system program for model calculation of the liquid metal reactor can be transmitted to the STAR-CCM+ software through the JAVA extension interface of the STAR-CCM+ software.

[0073] In some embodiments, the terminal device can also set the input file path and the output file path of the one-dimensional system program when establishing the program interface between the one-dimensional system program and the three-dimensional computational fluid dynamics program. The input file path is a file path used by the terminal device to obtain the model parameters (such as the fluid temperature and flow of the liquid metal reactor) of the liquid metal reactor through the one-dimensional system program for model calculation of the liquid metal reactor, and the output file path is a file path used to store the model calculation result obtained by the one-dimensional system program for model calculation. After the terminal device sets the input file path and the output file path of the one-dimensional system program, the terminal device can start the one-dimensional system program to calculate the one-dimensional model corresponding to the liquid metal reactor through the one-dimensional system program.

[0074] Exemplarily, refer to FIG. 2, which is a schematic diagram of a coupling architecture involved in the simulation method of the thermal stratification phenomenon provided in the embodiments of the present application. As shown in FIG. 2, when starting a one-dimensional system program to perform model calculation on the one-dimensional model corresponding to the liquid metal reactor, the terminal device can call a local method in the established JAVA extension interface of the three-dimensional computational fluid dynamics program STAR-CCM+ software, to interact with the dynamic link library NUSOLSYS of the one-dimensional system program in real time through the local method, to pass the model calculation result obtained by the one-dimensional system program performing model calculation to the STAR-CCM+ software, and to realize coupling calculation of the STAR-CCM+ software and the one-dimensional system program.

[0075] It should be noted that the local method called by the terminal device in the JAVA extension interface of the STAR-CCM+ software can be a program method corresponding to the coupling interface code as shown below.

[0076] In the embodiments of the present application, through the process that the terminal device respectively performs model calculation on the model corresponding to the liquid metal reactor by using the one-dimensional system program and the three-dimensional computational fluid dynamics program, the model calculation result obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively performing model calculation on the model of the liquid metal reactor is interacted through the pre-established program interface between the one-dimensional system program and the three-dimensional computational fluid dynamics program, so that the one-dimensional system program and the three-dimensional computational fluid dynamics program perform coupling calculation. In this way, the terminal device can obtain the real-time calculation result obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program performing coupling calculation, and obtain the simulation result of simulating the thermal stratification phenomenon of the liquid metal reactor in the transient operation process.

[0077] That is to say, the simulation method of the thermal stratification phenomenon provided in the embodiments of the present application can combine the one-dimensional system program and the three-dimensional computational fluid dynamics program, and comprehensively utilize the advantages of the three-dimensional computational fluid dynamics program in accurately simulating the thermal stratification phenomenon and the advantages of the system program in quickly simulating the overall transient response of the reactor system, to realize accurate simulation of the thermal stratification phenomenon in the transient operation process of the reactor, that is, to improve the simulation accuracy, save the computing resources, and improve the computing efficiency, so as to further support the performance matching of the liquid metal reactor in-vessel component and the development of the overall thermal-hydraulic performance of the reactor.

[0078] In some embodiments, the model calculation result obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively performing model calculation on the model corresponding to the liquid metal reactor includes a first model calculation result and a second model calculation result.

[0079] Please refer to FIG. 3, which is a timing diagram of model calculation of a liquid metal reactor involved in the simulation method of thermal stratification phenomenon provided in the embodiments of the present application.

[0080] As shown in FIG. 3, in some embodiments, the simulation method of thermal stratification phenomenon provided in the embodiments of the present application can further include steps S301-S302.

[0081] In step S301, the first one-dimensional model of the first part of the liquid metal reactor is calculated by the one-dimensional system program to obtain the first model calculation result.

[0082] In step S302, a three-dimensional model of the second part of the liquid metal reactor is calculated by the three-dimensional computational fluid dynamics program to obtain the second model calculation result.

[0083] The terminal device can divide the liquid metal reactor into a first part P1 and a second part P2 according to the region where thermal stratification occurs in the liquid metal reactor. The inlet of the second part P2 is connected to the outlet of the first part P1, and the outlet of the second part P2 is connected to the inlet of the first part P1. Then, the terminal device establishes a first one-dimensional model M2 for the first part P1 through a one-dimensional system program, and calculates the first one-dimensional model M2 through the one-dimensional system program to obtain a first model calculation result. In addition, the terminal device also establishes a three-dimensional model M3 for the second part P2 through a three-dimensional computational fluid dynamics program, and calculates the three-dimensional model M3 through the three-dimensional computational fluid dynamics program to obtain a second model calculation result.

[0084] Exemplarily, please refer to FIG. 4, which is a structural region division diagram of a liquid metal reactor involved in the simulation method of thermal stratification phenomenon provided in the embodiments of the present application.

[0085] As shown in FIG. 4, the terminal device divides the part of the primary circuit from the steam generator inlet to the main pump, the lower plenum of the reactor core, the outlet of the reactor core, and the part of the secondary circuit, into the first part P1 of the liquid metal reactor according to the region where thermal stratification occurs in the liquid metal reactor, and divides the liquid metal pool in the upper plenum of the reactor pressure vessel from the outlet of the reactor core to the steam generator inlet, into the second part P2 of the liquid metal reactor.

[0086] It should be noted that the first model calculation result includes the flow rate, temperature and other data of the outlet boundary of the first part P1, and the second model calculation result includes the threshold inlet and outlet mass flow rate and threshold monitoring point location temperature and other data of the connection between the outlet of the second part P2 and the inlet of the first part P1. It should be understood that the first model calculation result and the second model calculation result can of course each further include other data not listed here, and the simulation method of the thermal stratification phenomenon provided in the embodiments of the present application does not specifically limit the types of data included in the first model calculation result and the second model calculation result.

[0087] In the embodiment, the terminal device divides the liquid metal reactor into the first part P1 and the second part P2 according to the region where the thermal stratification occurs in the liquid metal reactor, and the inlet of the second part P2 is connected to the outlet of the first part P1, and the outlet of the second part P2 is connected to the inlet of the first part P1. Thus, the terminal device establishes a first one-dimensional model M2 for the first part P1 through the one-dimensional system program, and obtains the first model calculation result by calculating the first one-dimensional model M2 through the one-dimensional system program, and establishes a three-dimensional model M3 for the second part P2 through the three-dimensional computational fluid dynamics program, and obtains the second model calculation result by calculating the three-dimensional model M3 through the three-dimensional computational fluid dynamics program.

[0088] In this way, the first model calculation result and the second model calculation result are interacted between the one-dimensional system program and the three-dimensional computational fluid dynamics program through the program interface, so that the one-dimensional system program and the three-dimensional computational fluid dynamics program are coupled to calculate, and the thermal stratification phenomenon occurring at the interface between the first part P1 and the second part P2 of the liquid metal reactor can be accurately simulated in combination with the one-dimensional system program and the three-dimensional computational fluid dynamics program.

[0089] In some embodiments, when the terminal device calculates the three-dimensional model M3 through the three-dimensional computational fluid dynamics program, the terminal device takes the first model calculation result obtained by the one-dimensional system program calculating the first one-dimensional model M2 as the input of the three-dimensional model M3, and when the terminal device calculates the first one-dimensional model M2 through the one-dimensional system program for a non-first time step, the terminal device takes the second model calculation result obtained by the three-dimensional computational fluid dynamics program calculating the three-dimensional model M3 as the input of the first one-dimensional model M2. When the terminal device calculates the first one-dimensional model M2 through the one-dimensional system program for a first time step, the input of the first one-dimensional model M2 is the third model calculation result obtained by the terminal device calculating the second one-dimensional model M1 corresponding to the entire liquid metal reactor through the one-dimensional system program.

[0090] Please refer to FIG. 5, which is another time sequence diagram of model calculation of the liquid metal reactor involved in the simulation method of thermal stratification phenomenon provided in the embodiments of the present application.

[0091] As shown in FIG. 5, in some embodiments, before the step S301, the simulation method of thermal stratification phenomenon provided in the embodiments of the present application can further include a step S501.

[0092] In the step S501, the second one-dimensional model of the liquid metal reactor as a whole is calculated by the one-dimensional system program to obtain a third model calculation result.

[0093] Before the terminal device performs the step S301, i.e., before the first one-dimensional model of the first part of the liquid metal reactor is calculated by the one-dimensional system program to obtain the first model calculation result, the terminal device first establishes a second one-dimensional model M1 for the liquid metal reactor as a whole by the one-dimensional system program, and then acquires relevant model parameters (such as fluid temperature and flow of the liquid metal reactor) of the liquid metal reactor by the one-dimensional system program according to a pre-set input file path, and takes the acquired model parameters as input of the second one-dimensional model M1 to calculate the second one-dimensional model M1 to obtain a third model calculation result.

[0094] It should be noted that the third model calculation result includes fluid temperature and flow. It should be understood that, based on different design needs of actual applications, the third model calculation result can of course include other data not listed here, and the simulation method of thermal stratification phenomenon provided in the embodiments of the present application is not specifically limited to the type of data included in the third model calculation result.

[0095] In some embodiments, the step S301 of calculating the first one-dimensional model of the first part of the liquid metal reactor by the one-dimensional system program to obtain the first model calculation result can include:

[0096] The third model calculation result is taken as input of the first one-dimensional model to calculate the first one-dimensional model by the one-dimensional system program to obtain the first model calculation result.

[0097] In the present embodiment, after the terminal device calculates the third model calculation result by the one-dimensional system program for the second one-dimensional model M1, the terminal device takes the third model calculation result as input of the first one-dimensional model M2, and then calculates the first one-dimensional model M2 by the one-dimensional system program to obtain the first model calculation result.

[0098] Exemplarily, the terminal device establishes a second one-dimensional model M1 of the liquid metal reactor as a whole (including the first part P1 and the second part P2) through a one-dimensional system program, and obtains relevant model parameters of the liquid metal reactor from an input file path through the one-dimensional system program, so as to perform simulation through the one-dimensional system program to calculate the second one-dimensional model M1, so that the fluid temperature and flow boundary conditions at the interface position of the first part P1 and the second part P2 in the second one-dimensional model M1 model can be obtained, which are third model calculation results obtained by the one-dimensional system program in the calculation of the second one-dimensional model M1. Then, in the case that the terminal device establishes a first one-dimensional model M2 for the first part P1 through the one-dimensional system program, the terminal device can take the third model calculation results as input of the first one-dimensional model M2, so as to use the third model calculation results as the corresponding import initial boundary conditions of the first one-dimensional model M2, and perform calculation of the first time step of the first one-dimensional model M2 through the one-dimensional system program, so as to obtain system overall parameters (such as outlet boundary flow, temperature and other data) of the first one-dimensional model M2 in the first time step. The system overall parameters are first model calculation results obtained by the one-dimensional system program in the calculation of the first one-dimensional model M2.

[0099] In this embodiment, the terminal device first establishes a second one-dimensional model M1 for the liquid metal reactor as a whole through a one-dimensional system program, and obtains third model calculation results by calculating the second one-dimensional model M1. Then, the terminal device takes the third model calculation results as input of a first one-dimensional model M2, so as to perform calculation of the first time step of the first one-dimensional model M2 through the one-dimensional system program to obtain first model calculation results. In this way, the terminal device can start coupling calculation by combining the one-dimensional system program and the three-dimensional computational fluid dynamics program through the program interface to interact the first model calculation results to the three-dimensional computational fluid dynamics program, so as to realize accurate simulation of the thermal stratification phenomenon occurring at the interface between the first part P1 and the second part P2 of the liquid metal reactor.

[0100] In some embodiments, please refer to FIG. 6, which is a step flow diagram of coupling calculation involved in the simulation method of the thermal stratification phenomenon provided in the embodiments of the present application.

[0101] As shown in FIG. 6, when the terminal device performs the above-mentioned step S101 and controls the one-dimensional system program and the three-dimensional computational fluid dynamics program to perform coupling calculation based on the interaction of model calculation results in the pre-configured program interface, the terminal device can perform steps S601 to S604.

[0102] In step S601, the first model calculation results are transmitted to the three-dimensional computational fluid dynamics program through the program interface when the one-dimensional system program obtains the first model calculation results.

[0103] Step S602, the first model calculation result is calculated by the three-dimensional computational fluid dynamics program as the input of the three-dimensional model to the three-dimensional model K time step calculation; K is a positive integer greater than or equal to 1;

[0104] Step S603, in the case of the second model calculation result obtained by the three-dimensional computational fluid dynamics program, the second model calculation result is transmitted to the one-dimensional system program through the program interface;

[0105] Step S604, the second model calculation result is calculated by the one-dimensional system program as the input of the first one-dimensional model to the first one-dimensional model K+1 time step calculation.

[0106] In the process of the terminal device interacting the first model calculation result and the second model calculation result between the one-dimensional system program and the three-dimensional computational fluid dynamics program through the program interface to control the one-dimensional system program and the three-dimensional computational fluid dynamics program to perform coupled calculation, in the case that the one-dimensional system program calculates the first one-dimensional model M2 to obtain the first model calculation result, the terminal device transmits the first model calculation result to the three-dimensional computational fluid dynamics program through the program interface, so that the terminal device can calculate the three-dimensional model M3 through the three-dimensional computational fluid dynamics program with the first model calculation result as the input of the three-dimensional model M3. Among them, the calculation of the three-dimensional model M3 by the terminal device through the three-dimensional computational fluid dynamics program in the case that the first model calculation result is the input of the three-dimensional model M3 can be the first time step calculation, or it can be the non-first time step calculation. For example, in the case that the terminal device calculates the first one-dimensional model M2 through the one-dimensional system program to obtain the first model calculation result, the terminal device interacts the first model calculation result with the three-dimensional computational fluid dynamics program as the input of the three-dimensional model M3, and the K time step calculation of the three-dimensional model M3 through the three-dimensional computational fluid dynamics program is the first time step calculation, and K is equal to 1 at this time. In addition, in the case that the terminal device calculates the first one-dimensional model M2 through the one-dimensional system program to obtain the first model calculation result, the terminal device interacts the first model calculation result with the three-dimensional computational fluid dynamics program as the input of the three-dimensional model M3, and the K time step calculation of the three-dimensional model M3 through the three-dimensional computational fluid dynamics program is the non-first time step calculation, and K is greater than 1 and equal to N at this time.

[0107] When the terminal device obtains the second model calculation result by performing the Kth time step calculation on the three-dimensional model M3 through the three-dimensional computational fluid dynamics program, the terminal device continues to pass the second model calculation result to the one-dimensional system program through the program interface, so that the terminal device can pass the second model calculation result to the one-dimensional system program as an input of the first one-dimensional model M2 to perform the K+1th time step calculation on the first one-dimensional model M3.

[0108] After the terminal device obtains the first model calculation result by performing the K+1th time step calculation on the first one-dimensional model M2 through the one-dimensional system program, the terminal device continues to pass the first model calculation result to the three-dimensional computational fluid dynamics program through the program interface, and then passes the second model calculation result to the one-dimensional system program through the program interface after obtaining the first model calculation result by performing the K+1th time step calculation on the three-dimensional model M3 through the three-dimensional computational fluid dynamics program to perform the next time step calculation. In this way, after the terminal device passes the first model calculation result and the second model calculation result through the program interface to enable the one-dimensional system program and the three-dimensional computational fluid dynamics program to perform the coupling calculation for a time step length reaching the preset transient calculation requirement time, the terminal device determines that the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program is completed, so that the terminal device can perform the step S103 to obtain the real-time calculation result obtained by the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program, and obtain the simulation result for simulating the thermal stratification phenomenon of the liquid metal reactor in the transient operation process.

[0109] It should be noted that the transient calculation requirement time T can be defined by the user. It should be understood that T can have different sizes in different feasible implementations based on different design needs of actual applications, and the simulation method for the thermal stratification phenomenon provided in the embodiments of the present application does not limit the specific value of T.

[0110] Exemplarily, please refer to FIG. 7, which is a coupling calculation process schematic diagram involved in the simulation method for the thermal stratification phenomenon provided in the embodiments of the present application.

[0111] As shown in FIG. 7, the terminal device uses the one-dimensional system program to calculate the second one-dimensional model M1 to provide the third model calculation result at the t0 time boundary condition, that is, the terminal device establishes the second one-dimensional model M1 of the liquid metal reactor as a whole through the one-dimensional system program, and performs calculation on the second one-dimensional model M1 through the dynamic link library NUSOLSYS loaded in the one-dimensional system program to perform simulation, so as to obtain the third model calculation result (the fluid temperature and flow boundary conditions of the interface position between the second part P2 and the first part P1 in the second one-dimensional model M1), and then the terminal device uses the third model calculation result as the input of the first one-dimensional model M2, so that in the case that the terminal device establishes the first one-dimensional model M2 of the first part P1 of the liquid metal reactor through the one-dimensional system program, the fluid temperature and flow boundary conditions of the interface position between the second part P2 and the first part P1 in the second one-dimensional model M1 are used as the corresponding import initial boundary conditions of the first one-dimensional model M2. Based on this, the terminal device is coupled for calculation in combination with the one-dimensional system program and the three-dimensional computational fluid mechanics program.

[0112] At the t0 time, the terminal device performs the first time step calculation on the first one-dimensional model M2 through the loading of NUSOLSYS to obtain the first model calculation result (the system overall parameters of the first one-dimensional model M2 at the first time step: outlet boundary flow, temperature data). Then, the terminal device transmits the first model calculation result to STAR-CCM+ through the JAVA extension interface of the established three-dimensional fluid mechanics program STAR-CCM+, so that in the case that the three-dimensional model M3 of the second part P2 of the liquid metal reactor is established using STAR-CCM+, the outlet boundary flow and temperature data obtained by the first one-dimensional model M2 at the first time step are used as the inlet boundary conditions of the three-dimensional model M3 (the first model calculation result is used as the input of the three-dimensional model M3), and the terminal device also performs the first time step calculation on the three-dimensional model M3 through STAR-CCM+ to obtain the second model calculation result (threshold import and export mass flow and threshold monitoring point position temperature).

[0113] At the t1 time, in the case that the terminal device judges that the second model calculation result obtained by the calculation of the three-dimensional model M3 through STAR-CCM+ is converged, the terminal device transmits the second model calculation result to the one-dimensional system program through the JAVA extension interface as the new inlet boundary conditions of the first one-dimensional model M2, and performs the next time step calculation on the first one-dimensional model M2 through the loading of NUSOLSYS. Subsequently, the terminal device transmits the first model result obtained by the calculation of the first one-dimensional model M2 to STAR-CCM+ through the JAVA extension interface as the inlet boundary conditions of the three-dimensional model M3, so that the next time step calculation on the three-dimensional model M3 is continued through STAR-CCM+.

[0114] At the time t2, the terminal device also judges that the second model calculation result obtained by STAR-CCM+ calculating the three-dimensional model M3 converges, and then transmits the second model calculation result to the one-dimensional system program through the JAVA extension interface as the new inlet boundary condition of the first one-dimensional model M2, and calculates the next time step of the first one-dimensional model M2 through NUSOLSYS. Subsequently, the terminal device also transmits the first model calculation result obtained by calculating the first one-dimensional model M2 to STAR-CCM+ through the JAVA extension interface as the inlet boundary condition of the three-dimensional model M3, so as to continue calculating the next time step of the three-dimensional model M3 through STAR-CCM+.

[0115] In the process of continuing the coupled calculation of the one-dimensional system program and STAR-CCM+ by the terminal device performing the same operation as that at the time t1 or t2, the terminal device also judges whether the time step ti reaches the transient calculation required time T. If ti is less than or equal to T, the terminal device continues to perform the same operation as that at the time t1 or t2 to continue the coupled calculation of the one-dimensional system program and STAR-CCM+. If the terminal device judges at the time ti+1 that the time step ti+1 reaches the transient calculation required time (ti+1 is greater than or equal to T), the terminal device confirms that the one-dimensional system program and STAR-CCM+ have completed the entire coupled calculation process, and thus no longer transmits the second model calculation result to the one-dimensional system program to end the calculation.

[0116] In this embodiment, the terminal device uses the one-dimensional system program to calculate the first one-dimensional model to obtain the first model calculation result, and then transmits the first model calculation result to the three-dimensional computational fluid dynamics program through the program interface, so as to calculate the three-dimensional model through the three-dimensional computational fluid dynamics program using the first model calculation result as the input of the three-dimensional model to obtain the second model calculation result. Subsequently, in the case where the time step does not reach the preset transient calculation required time, the terminal device transmits the second model calculation result to the one-dimensional system program through the program interface, so as to use the second model calculation result as the new input of the first one-dimensional model, and thus makes the one-dimensional system program and the three-dimensional computational fluid dynamics program perform the coupled calculation of the next time step. In this way, the one-dimensional system program and the three-dimensional computational fluid dynamics program can be combined to comprehensively use the advantages of the three-dimensional computational fluid dynamics program in accurately simulating the thermal stratification phenomenon and the advantages of the system program in quickly simulating the overall transient response of the reactor system, so as to accurately simulate the thermal stratification phenomenon in the transient operation of the reactor.

[0117] After the terminal device performs the calculation on the three-dimensional model M3 by the three-dimensional computational fluid dynamics program, the terminal device further determines whether the second model calculation result of the three-dimensional model M3 converges. In the case that the second model calculation result does not converge, the terminal device adjusts the model initial value of the first one-dimensional model to perform the model calculation again.

[0118] In some embodiments, after the three-dimensional model of the second part of the liquid metal reactor is calculated by the three-dimensional computational fluid dynamics program to obtain the second model calculation result in step S302, the simulation method of the thermal stratification phenomenon provided by the embodiments of the present application can further include:

[0119] In the case that the second model calculation result is obtained by the three-dimensional computational fluid dynamics program, if the second model calculation result does not converge, the input of the first one-dimensional model is adjusted and the first one-dimensional model is calculated again by the one-dimensional system program to obtain the first model calculation result.

[0120] After the terminal device transmits the first model calculation result to the three-dimensional computational fluid dynamics program through the program interface to take the first model calculation result as the input of the three-dimensional model M3 and performs the calculation on the three-dimensional model M3 by the three-dimensional computational fluid dynamics program, the terminal device determines whether the second model calculation result obtained by the calculation on the three-dimensional model M3 converges. In the case that the second model calculation result does not converge, the terminal device confirms that the input of the three-dimensional model M3 needs to be adjusted and calculated again. Since the input of the three-dimensional model M3 is the first model calculation result, the terminal device adjusts the model initial value of the first one-dimensional model M2 to perform the K-time step calculation on the first one-dimensional model M2 again to obtain a new first model calculation result. Then, the terminal device transmits the new first model calculation result to the three-dimensional computational fluid dynamics program through the program interface as the new input of the three-dimensional model M3, and performs the K-time step calculation on the three-dimensional model M3 by the three-dimensional computational fluid dynamics program to obtain a new second model calculation result. For the new second model calculation result, the terminal device still determines whether the second model calculation result converges. If the second model calculation result still does not converge, the terminal device continues to adjust the model initial value of the first one-dimensional model to perform the model calculation again until the second model calculation result converges.

[0121] In some embodiments, the simulation method of the thermal stratification phenomenon provided by the embodiments of the present application can further include:

[0122] In the case that the threshold inlet and outlet mass flow rate and the threshold monitoring point position temperature are balanced, it is determined that the second model calculation result converges;

[0123] In a case where any one of the threshold inlet-outlet mass flow rate and the threshold monitoring point location temperature is not balanced, it is determined that the second model calculation result is not converged.

[0124] For the second model calculation result, the terminal device monitors whether the threshold inlet-outlet mass flow rate in the second model calculation result is balanced and whether the monitoring threshold monitoring point location temperature in the second model calculation result is balanced. In a case where the terminal device monitors that the threshold inlet-outlet mass flow rate is balanced and also monitors that the threshold monitoring point location temperature is balanced, the terminal device determines that the second model calculation result is converged. If the terminal device monitors that the threshold inlet-outlet mass flow rate is balanced but monitors that the threshold monitoring point location temperature is not balanced, the terminal device determines that the second model calculation result is not converged. Or, if the terminal device monitors that the threshold inlet-outlet mass flow rate is not balanced but monitors that the threshold monitoring point location temperature is balanced, the terminal device also determines that the second model calculation result is not converged. Also or, if the terminal device monitors that the threshold inlet-outlet mass flow rate is not balanced and also monitors that the threshold monitoring point location temperature is not balanced, the terminal device also determines that the second model calculation result is not converged.

[0125] In some embodiments, in the step S603, the second model calculation result is transmitted to the one-dimensional system program through the program interface, which can include:

[0126] In a case where the second model calculation result is converged, the second model calculation result is transmitted to the one-dimensional system program through the program interface.

[0127] In the policy and concept of coupling the one-dimensional system program and the three-dimensional computational fluid dynamics program for calculation, for the second model calculation result, the terminal device only transmits the second model calculation result to the one-dimensional system program through the program interface in a case where the second model calculation result is converged, so that the one-dimensional system program and the three-dimensional computational fluid dynamics program perform coupling calculation of the next time step.

[0128] In this embodiment, whether the second model calculation result of the three-dimensional model M3 is converged is judged by the terminal device. In a case where the second model calculation result is not converged, the terminal device adjusts the model initial value of the first one-dimensional model to re-perform model calculation until the second model calculation result is converged. In this way, the terminal device only transmits the second model calculation result to the one-dimensional system program through the program interface in a case where the second model calculation result is converged, so that the one-dimensional system program and the three-dimensional computational fluid dynamics program perform coupling calculation of the next time step. Therefore, it can be ensured that the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program is performed, and the thermal stratification phenomenon in the transient operation process of the reactor is accurately simulated.

[0129] Next, a preferred embodiment of the terminal device applying the simulation method of the thermal stratification phenomenon provided by the present application is proposed.

[0130] Exemplarily, it is assumed that the three-dimensional computational fluid dynamics program is STAR-CCM+ software. When the terminal device applies the simulation method of the thermal stratification phenomenon provided by the present application, the input file path and the output file path of the one-dimensional system program are set in advance, and the JAVA extension interface of the STAR-CCM+ software is established as the program interface between the STAR-CCM+ software and the one-dimensional system program. In this way, the terminal device can interact with the dynamic link library NUSOLSYS of the one-dimensional system program in real time through the local method in the JAVA extension interface, transfer the model calculation results obtained by the one-dimensional system program to the STAR-CCM+ software, and realize the coupled calculation of the STAR-CCM+ software and the one-dimensional system program.

[0131] Before the terminal device calls the local method in the JAVA extension interface to control the coupled calculation of the STAR-CCM+ software and the one-dimensional system program, the terminal device divides the one-loop part from the steam generator inlet to the main pump, the lower core cavity, the core outlet, and the two-loop part into the first part P1 of the liquid metal reactor according to the region where the thermal stratification of the liquid metal reactor occurs, and divides the liquid metal pool in the upper cavity of the reactor pressure vessel from the core outlet to the steam generator inlet into the second part P2 of the liquid metal reactor.

[0132] After that, in the process of the terminal device calling the local method in the JAVA extension interface to control the coupled calculation of the STAR-CCM+ software and the one-dimensional system program, the terminal device first uses the one-dimensional system program to calculate the second one-dimensional model M1 to obtain the third model calculation result providing the boundary condition at time t0. That is, the terminal device establishes the second one-dimensional model M1 of the liquid metal reactor as a whole (including the first part P1 and the second part P2) through the one-dimensional system program, and performs simulation calculation on the second one-dimensional model M1 by loading the dynamic link library NUSOLSYS in the one-dimensional system program, thereby obtaining the third model calculation result (the fluid temperature and flow boundary conditions at the interface position of the second part P2 and the first part P1 in the second one-dimensional model M1). Then, the terminal device uses the third model calculation result as the input of the first one-dimensional model M2, so that in the case that the terminal device establishes the first one-dimensional model M2 of the first part P1 of the liquid metal reactor through the one-dimensional system program, the fluid temperature and flow boundary conditions at the interface position of the second part P2 and the first part P1 in the second one-dimensional model M1 are used as the corresponding import initial boundary conditions of the first one-dimensional model M2. Based on this, the terminal device performs coupled calculation in combination with the one-dimensional system program and the three-dimensional computational fluid dynamics program.

[0133] At time t0, the terminal device performs the first time step calculation of the first one-dimensional model M2 by loading NUSOLSYS to obtain the first model calculation result (the system overall parameters of the first one-dimensional model M2 at the first time step: outlet boundary flow rate, temperature data). Then, the terminal device transmits the first model calculation result to STAR-CCM+ through the established JAVA extension interface of STAR-CCM+, so that in the case where the three-dimensional model M3 of the second part P2 of the liquid metal reactor is established by using STAR-CCM+, the outlet boundary flow rate and temperature data obtained by the first one-dimensional model M2 at the first time step are taken as the inlet boundary conditions of the three-dimensional model M3 (the first model calculation result is taken as the input of the three-dimensional model M3), and the terminal device also performs the first time step calculation of the three-dimensional model M3 by STAR-CCM+ to obtain the second model calculation result (threshold inlet and outlet mass flow rate and threshold monitoring point position temperature).

[0134] At time t1, in the case where the terminal device judges that the second model calculation result obtained by STAR-CCM+ calculation of the three-dimensional model M3 converges, the terminal device transmits the second model calculation result to the one-dimensional system program through the JAVA extension interface as the new inlet boundary condition of the first one-dimensional model M2, and performs the next time step calculation of the first one-dimensional model M2 by NUSOLSYS. Subsequently, the terminal device transmits the first model result obtained by the calculation of the first one-dimensional model M2 to STAR-CCM+ through the JAVA extension interface as the inlet boundary condition of the three-dimensional model M3, so that STAR-CCM+ continues to perform the next time step calculation of the three-dimensional model M3.

[0135] At time t2, in the case where the terminal device judges that the second model calculation result obtained by STAR-CCM+ calculation of the three-dimensional model M3 converges, the terminal device transmits the second model calculation result to the one-dimensional system program through the JAVA extension interface as the new inlet boundary condition of the first one-dimensional model M2, and performs the next time step calculation of the first one-dimensional model M2 by NUSOLSYS. Subsequently, the terminal device also transmits the first model result obtained by the calculation of the first one-dimensional model M2 to STAR-CCM+ through the JAVA extension interface as the inlet boundary condition of the three-dimensional model M3, so that STAR-CCM+ continues to perform the next time step calculation of the three-dimensional model M3.

[0136] In the process of performing the same operation as that at the time t1 or the time t2 to make the one-dimensional system program and the STAR-CCM+ continue the coupling calculation, the terminal device also judges whether the time step t1 reaches the transient calculation required time T. If t1 is less than or equal to T, the terminal device continues to perform the same operation as that at the time t1 or the time t2 to make the one-dimensional system program and the STAR-CCM+ continue the coupling calculation. If the terminal judges that the time step t1+1 reaches the transient calculation required time (t1+1 is greater than or equal to T) at the time t1+1, the terminal device confirms that the one-dimensional system program and the STAR-CCM+ have completed the entire coupling calculation process, and thus no longer continues to pass the second model calculation result to the one-dimensional system program to end the calculation.

[0137] Finally, the terminal device can obtain the simulation result of simulating the thermal stratification phenomenon of the liquid metal reactor in the transient operation process by obtaining the real-time calculation result obtained by the coupling calculation of the one-dimensional system program and the STAR-CCM+.

[0138] Next, the terminal device applies the simulation method of the thermal stratification phenomenon provided in the embodiments of the present application to a real case of coupling calculation in combination with the one-dimensional system program and the three-dimensional computational fluid dynamics program.

[0139] Please refer to FIG. 8, which is an example schematic diagram of a coupling calculation loop involved in the simulation method of the thermal stratification phenomenon provided in the embodiments of the present application.

[0140] As shown in FIG. 8, the terminal device establishes a second one-dimensional model M1 for the entire liquid metal reactor by the one-dimensional system program, which is a rectangular loop. For the top heating horizontal pipe in the rectangular loop, the terminal device establishes a three-dimensional model M3 by the three-dimensional computational fluid dynamics program, and for the remaining positions in the rectangular loop except the top heating horizontal pipe, such as the pump, the pipeline and the pressure boundary, the terminal device establishes a first one-dimensional model M2 by the one-dimensional system program. In this way, the coupling surface between the one-dimensional system program and the three-dimensional computational fluid dynamics program is located at the inlet and outlet positions of the top heating horizontal pipe. The terminal device takes the loop initial temperature 500 ℃, the heating horizontal pipe length 3 m, the heating horizontal pipe length diameter 10 cm, and the heating section length of the heating horizontal pipe 1 m as the input of the second one-dimensional model M1 to perform calculation on the second one-dimensional model M1. Then, the terminal device further takes the third model calculation result obtained by the calculation on the second one-dimensional model M1 as the input of the first one-dimensional model M2 to perform the first time step calculation on the first one-dimensional model M2, and passes the first model calculation result obtained by the calculation on the first one-dimensional model M2 to the three-dimensional computational fluid dynamics program to perform the above-mentioned coupling calculation process in combination with the one-dimensional system program and the three-dimensional computational fluid dynamics program.

[0141] After the time step of the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program reaches the preset transient calculation requirement time, the terminal device obtains the real-time calculation result obtained by the coupling calculation. The outlet flow rate obtained by the one-dimensional system program performing calculation on the first one-dimensional model M2 is 0.1564 m / s, and the outlet flow rate obtained by the three-dimensional computational fluid dynamics program performing calculation on the three-dimensional model M3 is 0.1579 m / s. The two outlet flow rates are in good agreement, and the outlet flow rate obtained by the three-dimensional computational fluid dynamics program performing calculation on the three-dimensional model M3 is larger, which is probably caused by the volume expansion of the fluid after being heated. In addition, in the simulation result of the terminal device simulating the thermal stratification phenomenon of the liquid metal reactor in the transient operation process, the temperature cloud picture of the heating horizontal tube heating section obtained by the three-dimensional computational fluid dynamics program performing calculation on the three-dimensional model is shown in FIG. 9, and the flow rate transient curve of the heating horizontal tube heating section obtained by the three-dimensional computational fluid dynamics program performing calculation on the three-dimensional model is shown in FIG. 10.

[0142] It should be understood that, although each step in the flowchart involved in each embodiment as described above is displayed in sequence according to the indication of the arrow, these steps are not necessarily executed in sequence according to the indication of the arrow. Unless otherwise explicitly stated herein, the execution of these steps is not strictly limited in sequence, and these steps can be executed in other sequences. Moreover, at least part of the steps in the flowchart involved in each embodiment as described above can include multiple steps or multiple stages, which are not necessarily executed at the same time, but can be executed at different times, and the execution sequence of these steps or stages is not necessarily sequential, but can be executed alternately or alternately with at least part of other steps or steps or stages in other steps.

[0143] Based on the same inventive concept, the embodiment of the present application also provides a reactor core power parameter generation device for implementing the reactor core power parameter generation method. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme described in the above method, so the specific limitations in one or more reactor core power parameter generation device embodiments provided below can refer to the limitations of the reactor core power parameter generation method in the above text, which will not be repeated here.

[0144] In one embodiment, as shown in FIG. 11, a thermal stratification phenomenon simulation device is provided, which includes a coupling calculation module 1201 and a result acquisition module 1201, wherein:

[0145] The coupling calculation module 1201 is configured to control coupling calculation of a one-dimensional system program and a three-dimensional computational fluid dynamics program based on interaction model calculation results in a preconfigured program interface; the program interface is configured between the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the model calculation results include results obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program from simulation calculation of a liquid metal reactor.

[0146] The result acquisition module 1202 is configured to acquire real-time calculation results of the coupling calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the real-time calculation results include simulation results of thermal stratification of the liquid metal reactor in a transient operation process.

[0147] In some embodiments, the model calculation results include first model calculation results and second model calculation results.

[0148] The coupling calculation module 1201 is further configured to obtain the first model calculation results by calculating a first one-dimensional model of a first part of the liquid metal reactor through the one-dimensional system program, and obtain the second model calculation results by calculating a three-dimensional model of a second part of the liquid metal reactor through the three-dimensional computational fluid dynamics program; an inlet of the second part is connected to an outlet of the first part, and an outlet of the second part is connected to an inlet of the first part.

[0149] In some embodiments, the coupling calculation module 1201 is further configured to, in a case where the one-dimensional system program obtains the first model calculation results, transfer the first model calculation results to the three-dimensional computational fluid dynamics program through the program interface; perform Kth time step calculation on the three-dimensional model by taking the first model calculation results as input of the three-dimensional model through the three-dimensional computational fluid dynamics program; K is a positive integer greater than or equal to 1; in a case where the three-dimensional computational fluid dynamics program obtains the second model calculation results, transfer the second model calculation results to the one-dimensional system program through the program interface; and perform K+1th time step calculation on the first one-dimensional model by taking the second model calculation results as input of the first one-dimensional model through the one-dimensional system program.

[0150] In some embodiments, the coupling calculation module 1201 is further configured to, in a case where the three-dimensional computational fluid dynamics program obtains the second model calculation results, if the second model calculation results do not converge, adjust input of the first one-dimensional model and re-perform Kth time step calculation on the first one-dimensional model through the one-dimensional system program to obtain the first model calculation results.

[0151] In some embodiments, the coupling calculation module 1201 is further configured to, in a case where the second model calculation result converges, pass the second model calculation result to the one-dimensional system program through the program interface.

[0152] In some embodiments, the second model calculation result comprises: a threshold inlet-outlet mass flow rate and a threshold monitoring point location temperature at a connection between an outlet of the second part and an inlet of the first part.

[0153] The coupling calculation module 1201 is further configured to, in a case where both the threshold inlet-outlet mass flow rate and the threshold monitoring point location temperature are balanced, determine that the second model calculation result converges; and in a case where any one of the threshold inlet-outlet mass flow rate and the threshold monitoring point location temperature is unbalanced, determine that the second model calculation result does not converge.

[0154] In some embodiments, the coupling calculation module 1201 is further configured to perform a second one-dimensional model calculation of the liquid metal reactor as a whole through the one-dimensional system program to obtain a third model calculation result.

[0155] In some embodiments, the coupling calculation module 1201 is further configured to perform a calculation of the first one-dimensional model through the one-dimensional system program by taking the third model calculation result as an input of the first one-dimensional model to obtain the first model calculation result.

[0156] The above-mentioned various modules in the simulation device for the thermal stratification phenomenon can be realized by software, hardware, and combinations thereof, in whole or in part. The above-mentioned various modules can be embedded in or independent of a processor in a computer device in hardware form, or can be stored in a memory in the computer device in software form, so as to be called and executed by a processor to perform operations corresponding to the above-mentioned various modules.

[0157] In some embodiments, a computer device is provided, which can be a server, and an internal structure diagram of the computer device can be as shown in FIG. 12. The computer device includes a processor, a memory, an input / output interface (I / O), and a communication interface. Among them, the processor, the memory, and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store model parameters related to the liquid metal reactor. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with external terminals through network connection. The computer program is executed by the processor to implement a simulation method of thermal stratification phenomenon.

[0158] Those skilled in the art can understand that the structure shown in FIG. 12 is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. A specific computer device can include more or fewer components than those shown in the figure, or combine certain components, or have a different component arrangement.

[0159] In some embodiments, a computer device is provided, including a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the following steps:

[0160] Based on the model calculation result in the pre-configured program interface, the one-dimensional system program and the three-dimensional computational fluid dynamics program are controlled to perform coupled calculation; the program interface is configured between the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the model calculation result includes results obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively for simulating calculation of the liquid metal reactor;

[0161] The real-time calculation result of the coupled calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program is obtained, and the real-time calculation result includes a simulation result of thermal stratification phenomenon of the liquid metal reactor in a transient operation process.

[0162] In some embodiments, the model calculation result includes a first model calculation result and a second model calculation result; and the processor executing the computer program further implements the following steps:

[0163] The first model calculation result is obtained by calculating the first one-dimensional model of the first part of the liquid metal reactor through the one-dimensional system program;

[0164] The second model calculation result is obtained by calculating a three-dimensional model of a second part of the liquid metal reactor through the three-dimensional computational fluid dynamics program;

[0165] The second part is connected with the first part at an inlet of the second part and an outlet of the first part, and the second part is connected with the first part at an outlet of the second part and an inlet of the first part.

[0166] In some embodiments, the processor, when executing the computer program, further implements the following steps:

[0167] The first model calculation result is transmitted to the three-dimensional computational fluid dynamics program through the program interface when the one-dimensional system program obtains the first model calculation result;

[0168] The first model calculation result is taken as an input of the three-dimensional model to calculate the three-dimensional model at a Kth time step through the three-dimensional computational fluid dynamics program; K is a positive integer greater than or equal to 1;

[0169] The second model calculation result is transmitted to the one-dimensional system program through the program interface when the three-dimensional computational fluid dynamics program obtains the second model calculation result;

[0170] The second model calculation result is taken as an input of the first one-dimensional model to calculate the first one-dimensional model at a (K+1)th time step through the one-dimensional system program.

[0171] In some embodiments, the processor, when executing the computer program, further implements the following steps:

[0172] If the second model calculation result does not converge, the input of the first one-dimensional model is adjusted, and the first model calculation result is obtained by recalculating the first one-dimensional model at a Kth time step through the one-dimensional system program when the three-dimensional computational fluid dynamics program obtains the second model calculation result.

[0173] In some embodiments, the processor, when executing the computer program, further implements the following steps:

[0174] The second model calculation result is transmitted to the one-dimensional system program through the program interface when the second model calculation result converges.

[0175] In some embodiments, the second model calculation result comprises: a threshold inlet-outlet mass flow rate and a threshold monitoring point location temperature at a connection between an outlet of the second portion and an inlet of the first portion; and the processor, when executing the computer program, further implements the following steps:

[0176] In a case where both the threshold inlet-outlet mass flow rate and the threshold monitoring point location temperature are balanced, it is determined that the second model calculation result converges.

[0177] In a case where any one of the threshold inlet-outlet mass flow rate and the threshold monitoring point location temperature is not balanced, it is determined that the second model calculation result does not converge.

[0178] In some embodiments, the processor, when executing the computer program, further implements the following steps:

[0179] A third model calculation result is obtained by performing calculation on a second one-dimensional model of the liquid metal reactor as a whole by the one-dimensional system program.

[0180] In some embodiments, the processor, when executing the computer program, further implements the following steps:

[0181] The first model calculation result is obtained by performing calculation on the first one-dimensional model by the one-dimensional system program with the third model calculation result as an input of the first one-dimensional model.

[0182] In some embodiments, a computer readable storage medium is provided, and the computer readable storage medium stores a computer program. The computer program, when executed by a processor, implements the same steps as the processor in the above computer device when executing the computer program, and details are not repeated here.

[0183] In one embodiment, a computer program product is provided, and the computer program product comprises a computer program. The computer program, when executed by a processor, implements the same steps as the processor in the above computer device when executing the computer program, and details are not repeated here.

[0184] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiments of the methods can be included. Any reference to memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (Read-Only Memory, ROM), magnetic tape, floppy disk, flash memory, optical storage, high-density embedded non-volatile memory, resistive memory (ReRAM), magnetic resistive memory (Magneto resistive Random Access Memory, MRAM), ferroelectric memory (Ferroelectric Random Access Memory, FRAM), phase change memory (Phase Change Memory, PCM), graphene memory, etc. Volatile memory can include random access memory (Random Access Memory, RAM) or external cache memory, etc. As an illustration but not limitation, RAM can be in various forms, such as static random access memory (Static Random Access Memory, SRAM) or dynamic random access memory (Dynamic Random Access Memory, DRAM), etc. The database involved in the embodiments provided in the present application can include at least one of a relational database and a non-relational database. The non-relational database can include a distributed database based on a block chain, etc., without being limited thereto. The processor involved in the embodiments provided in the present application can be a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., without being limited thereto.

[0185] Any combination of the technical features of the above embodiments can be made. In order to make the description simple, all possible combinations of the technical features in the above embodiments are not described, however, as long as the combination of the technical features does not exist contradictory, it should be considered as the scope of the present application.

[0186] The above embodiments only express several implementation manners of the present application, and the description is more specific and detailed, but it should not be understood as a limitation on the scope of the patent of the present application. It should be pointed out that for ordinary skilled in the art, without departing from the concept of the present application, a number of modifications and improvements can be made, which are within the scope of protection of the present application. Therefore, the protection scope of the present application should be subject to the appended claims.

Claims

1. A method of simulating thermal delamination phenomena, characterized by, The method comprises: Based on the interaction model calculation result in the pre-configured program interface, control the one-dimensional system program and the three-dimensional computational fluid dynamics program to carry out coupled calculation; the program interface is configured between the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the model calculation result includes the results obtained by simulating and calculating the liquid metal reactor by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively; Obtain the real-time calculation result of the coupled calculation of the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the real-time calculation result includes the simulation result of the thermal stratification phenomenon of the liquid metal reactor in the transient operation process.

2. The method of claim 1, wherein, The model calculation result includes first model calculation result and second model calculation result; the method further comprises: The first one-dimensional model of the first part of the liquid metal reactor is calculated by the one-dimensional system program to obtain the first model calculation result; The three-dimensional model of the second part of the liquid metal reactor is calculated by the three-dimensional computational fluid dynamics program to obtain the second model calculation result; Wherein, the inlet of the second part is connected with the outlet of the first part, and the outlet of the second part is connected with the inlet of the first part.

3. The method of claim 2, wherein, The method based on the interaction model calculation result in the pre-configured program interface, control the one-dimensional system program and the three-dimensional computational fluid dynamics program to carry out coupled calculation, comprising: In the case that the one-dimensional system program obtains the first model calculation result, the first model calculation result is transmitted to the three-dimensional computational fluid dynamics program through the program interface; The first model calculation result is taken as the input of the three-dimensional model by the three-dimensional computational fluid dynamics program to carry out the calculation of the Kth time step of the three-dimensional model; K is a positive integer greater than or equal to 1; In the case that the three-dimensional computational fluid dynamics program obtains the second model calculation result, the second model calculation result is transmitted to the one-dimensional system program through the program interface; The second model calculation result is taken as the input of the first one-dimensional model by the one-dimensional system program to carry out the calculation of the K+1th time step of the first one-dimensional model.

4. The method of claim 3, wherein, After the first model calculation result is taken as the input of the three-dimensional model by the three-dimensional computational fluid dynamics program to carry out the calculation of the three-dimensional model, the method further comprises: In the case that the three-dimensional computational fluid dynamics program obtains the second model calculation result, if the second model calculation result does not converge, the input of the first one-dimensional model is adjusted, and the first one-dimensional model is recalculated by the one-dimensional system program to obtain the first model calculation result in the Kth time step; The second model calculation result is transmitted to the one-dimensional system program through the program interface, comprising: In the case that the second model calculation result converges, the second model calculation result is transmitted to the one-dimensional system program through the program interface.

5. The method of claim 4, wherein, The second model calculation result comprises: a threshold inlet-outlet mass flow rate and a threshold monitoring point location temperature at a connection between an outlet of the second part and an inlet of the first part; The method further comprises: determining that the second model calculation result converges when both the threshold inlet-outlet mass flow rate and the threshold monitoring point location temperature are balanced; determining that the second model calculation result does not converge when either the threshold inlet-outlet mass flow rate or the threshold monitoring point location temperature is not balanced.

6. The method according to any one of claims 2 to 5, characterized in that, Before the first model calculation result is obtained by calculating the first one-dimensional model of the first part of the liquid metal reactor through the one-dimensional system program, the method further comprises: obtaining a third model calculation result by calculating a second one-dimensional model of the entire liquid metal reactor through the one-dimensional system program; The first model calculation result is obtained by calculating the first one-dimensional model of the first part of the liquid metal reactor through the one-dimensional system program, comprising: calculating the first one-dimensional model through the one-dimensional system program by taking the third model calculation result as an input of the first one-dimensional model, to obtain the first model calculation result.

7. A simulation device of thermal delamination phenomenon, characterized by, The device comprises: a coupling calculation module configured to control a one-dimensional system program and a three-dimensional computational fluid dynamics program to perform coupled calculation based on interaction of model calculation results in a preconfigured program interface; the program interface is configured between the one-dimensional system program and the three-dimensional computational fluid dynamics program, and the model calculation results comprise results obtained by the one-dimensional system program and the three-dimensional computational fluid dynamics program respectively through simulation calculation of a liquid metal reactor; a result acquisition module configured to acquire real-time calculation results of the one-dimensional system program and the three-dimensional computational fluid dynamics program performing coupled calculation, the real-time calculation results comprising simulation results of thermal stratification phenomena of the liquid metal reactor in a transient operation process.

8. A computer device, comprising: The computer device comprises a memory and a processor, the memory stores a computer program, and the processor implements the steps of the method in any one of claims 1 to 6 when executing the computer program.

9. A computer-readable storage medium storing a computer program, the computer program comprising instructions that, when executed by a computer, cause the computer to perform the method of any one of claims 1 to 8. The computer program, when executed by a processor, implements the steps of the method in any one of claims 1 to 6.

10. A computer program product comprising a computer program, characterized in that, The computer program, when executed by a processor, implements the steps of the method in any one of claims 1 to 6. The computer program, when executed by a processor, implements the steps of the method in any one of claims 1 to 6.

Images