Parameter calculation method, device, medium and equipment for large-space component of nuclear reactor

Abstract

The invention discloses a parameter calculation method and device for a large-space component of a nuclear reactor, a medium and equipment. The method comprises the following steps: determining large-space components and other components in a nuclear reactor system; dividing a pressure solving matrix of a plurality of control bodies in the large-space component into a first structure matrix block and a second structure matrix block; solving the first structure matrix block in parallel to obtain vapor phase pressure of each control body in the large-space component; solving the second structure matrix block in parallel to obtain the liquid phase pressure of each control body in the large-space component; and according to the vapor phase pressure and the liquid phase pressure corresponding to the large-space component, iteratively updating the interface pressure parameter of the large-space component until convergence. Therefore, the violent stirring and mixing and heat conduction effects of large-space components of a nuclear reactor system are effectively simulated, and the calculation accuracy of a system analysis program is improved while the large-scale matrix discrete solving difficulty of the system analysis program is reduced.

Description

Technical Field

[0001] The present application relates to the technical field of nuclear reaction systems, and in particular to a method, device, medium and equipment for calculating parameters of large-space components of a nuclear reactor. Background Art

[0002] In nuclear power plants, there are large tank equipment such as core water tanks, injection tanks, or large sodium pools in fast reactors. When complex accident conditions such as LOCA occur, the corresponding equipment will mix cold and hot fluids, and cold and hot fluid stratification may occur, which has a great impact on the system flow and heat transfer process. Therefore, it is necessary to use the nuclear reactor system analysis program to model, analyze and calculate the parameters of each device in the nuclear power plant.

[0003] In the related technology, the entire system is divided into different computing modules, and the pressure matrices of different computing modules are solved in parallel. However, for large-volume tank equipment such as core water supply tanks and injection tanks, or large-space components such as large sodium pools in fast reactors, it is still impossible to accurately simulate the lateral and radial heat conduction and violent mixing effects of large-space components. Summary of the invention

[0004] In view of this, the present application provides a parameter calculation method, device, medium and equipment for large-space components of a nuclear reactor, which can effectively simulate the violent mixing and heat conduction effects of large-space components of a nuclear reactor system, while reducing the difficulty of discrete solution of large matrices of the system analysis program and improving the accuracy of the calculation of the system analysis program.

[0005] According to one aspect of the present application, a method for calculating parameters of a large space component of a nuclear reactor is provided, comprising:

[0006] Categorize nuclear reactor systems and identify large space components and other components in nuclear reactor systems;

[0007] Divide the pressure solution matrix of multiple control bodies in the large space component into a first structure matrix block and a second structure matrix block, wherein the length of the pressure solution matrix is ​​A times the number of control bodies of the large space component, A is a multiple of 2, and the length of the first structure matrix block or the second structure matrix block is the number of control bodies of the large space component;

[0008] Solve the first structure matrix block in parallel to obtain the vapor phase pressure of each control volume in the large space component;

[0009] Solve the second structural matrix block in parallel to obtain the liquid phase pressure of each control volume in the large space component;

[0010] According to the vapor phase pressure and liquid phase pressure corresponding to the large space component, the interface pressure parameters of the large space component are iteratively updated until convergence.

[0011] Optionally, the parameter calculation method of the large space component of the nuclear reactor also includes:

[0012] Establish the vapor and liquid phase conservation equations for large space components;

[0013] Discretize the vapor and liquid phase conservation equations, determine multiple control volumes in the large space component, and obtain the mass-energy conservation equations and momentum conservation equations corresponding to the large space component;

[0014] The mass-energy conservation equations and momentum conservation equations are transformed to establish the pressure solution matrix for multiple control volumes in large space components.

[0015] Optionally, the parameter calculation method of the large space component of the nuclear reactor also includes:

[0016] According to the initial boundary conditions of the mass-energy solution matrix, as well as the vapor pressure, liquid pressure and fluid velocity corresponding to the large space component, the mass-energy conservation equations are solved to obtain other node control parameters of each control body in the large space component except the pressure, where other node control parameters include temperature or cavitation fraction.

[0017] Optionally, the parameter calculation method of the large space component of the nuclear reactor also includes:

[0018] Update the vapor and liquid phase conservation equations based on the equivalent thermal diffusion coefficient of the mass transfer process of large space components;

[0019] Among them, the vapor and liquid phase conservation equations of large space components are expressed as:

[0020]

[0021] Where ρ is density; v is fluid velocity; t is time; z is spatial coordinate; λ is thermal diffusion coefficient; ε is equivalent thermal diffusion coefficient; S φ is the generalized source term; φ is the node control parameter of the control body.

[0022] Optionally, the parameter calculation method of the large space component of the nuclear reactor also includes:

[0023] According to the node control parameters of the control body in the large space component, the physical property parameters of the large space component are updated through the state equation and the physical property function;

[0024] The node control parameters include at least one of pressure, temperature, fluid velocity, and cavitation fraction; the pressure includes vapor phase pressure and liquid phase pressure; and the physical property parameters include at least one of flow rate, mass, and enthalpy.

[0025] Optionally, solving the first structure matrix block in parallel to obtain the vapor phase pressure of each control volume in the large space component includes:

[0026] Determine the coefficients of the pressure solution matrix according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle;

[0027] Using the LU decomposition method, according to the coefficients and initial parameters of the large space components, the matrix solving algorithm is called on different first structure matrix blocks for parallel calculation to obtain the vapor phase pressure of different control volumes in the current iteration cycle;

[0028] The initial parameters include initial inlet pressure, initial outlet pressure, initial inlet flow rate, initial outlet flow rate, initial inlet temperature, and initial outlet temperature.

[0029] Optionally, solving the second structural matrix block in parallel to obtain the liquid phase pressure of each control volume in the large space component includes:

[0030] Determine the coefficients of the pressure solution matrix according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle;

[0031] Using the LU decomposition method, according to the coefficients and initial parameters of the large space components, the matrix solving algorithm is called on different second structure matrix blocks for parallel calculation to obtain the liquid phase pressure of different control volumes in the current iteration cycle;

[0032] The initial parameters include initial inlet pressure, initial outlet pressure, initial inlet flow rate, initial outlet flow rate, initial inlet temperature, and initial outlet temperature.

[0033] Optionally, the parameter calculation method of the large space component of the nuclear reactor also includes:

[0034] Exchange the interface parameters of the connection between the large space component and the target component obtained in the current iteration cycle to update the initial boundary conditions of the pressure solution matrix of the next iteration cycle and / or the initial boundary conditions of the mass-energy solution matrix obtained by discretizing the vapor and liquid phase conservation equations;

[0035] The target component is another component connected to the large space component, and the interface parameter includes at least one of an interface pressure parameter, an interface temperature parameter, and an interface flow parameter.

[0036] According to another aspect of the present application, a parameter calculation device for a large space component of a nuclear reactor is provided, comprising:

[0037] A component classification module is used to classify the nuclear reactor system and identify large space components and other components in the nuclear reactor system;

[0038] A calculation module is used to divide the pressure solution matrix of multiple control bodies in the large space component into a first structure matrix block and a second structure matrix block, wherein the length of the pressure solution matrix is ​​A times the number of control bodies in the large space component, A is a multiple of 2, and the length of the first structure matrix block or the second structure matrix block is the number of control bodies in the large space component; and, the first structure matrix block is solved in parallel to obtain the vapor phase pressure of each control body in the large space component; and, the second structure matrix block is solved in parallel to obtain the liquid phase pressure of each control body in the large space component;

[0039] The parameter updating module is used to iteratively update the interface pressure parameters of the large space component according to the vapor phase pressure and liquid phase pressure corresponding to the large space component until convergence.

[0040] Optionally, the calculation module is also used to establish the vapor and liquid phase conservation equations for the large space component; discretize the vapor and liquid phase conservation equations, determine multiple control bodies in the large space component, and obtain the mass-energy conservation equations and momentum conservation equations corresponding to the large space component; transform the mass-energy conservation equations and momentum conservation equations, and establish the pressure solution matrix for multiple control bodies in the large space component.

[0041] Optionally, the calculation module is also used to solve the mass-energy conservation equations based on the initial boundary conditions of the mass-energy solution matrix, and the vapor pressure, liquid pressure and fluid velocity corresponding to the large space component, to obtain other node control parameters of each control body in the large space component except the pressure, wherein the other node control parameters include temperature or cavitation fraction.

[0042] Optionally, the parameter calculation device for the large space component of the nuclear reactor further includes:

[0043] Equivalent module, used to update the vapor and liquid phase conservation equations according to the equivalent thermal diffusion coefficient of the mass transfer process of large space components;

[0044] Among them, the vapor and liquid phase conservation equations of large space components are expressed as:

[0045]

[0046] Where ρ is density; v is fluid velocity; t is time; z is spatial coordinate; λ is thermal diffusion coefficient; ε is equivalent thermal diffusion coefficient; S φ is the generalized source term; φ is the node control parameter of the control body.

[0047] Optionally, the calculation module is further used to update the physical property parameters of the large space component through the state equation and the physical property function according to the node control parameters of the control body in the large space component;

[0048] The node control parameters include at least one of pressure, temperature, fluid velocity, and cavitation fraction; the pressure includes vapor phase pressure and liquid phase pressure; and the physical property parameters include at least one of flow rate, mass, and enthalpy.

[0049] Optionally, the calculation module is specifically used to determine the coefficients of the first structure matrix block according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle; using the LU decomposition method, according to the coefficients of the first structure matrix block and the initial parameters of the large space component, calling the matrix solution algorithm on different first structure matrix blocks for parallel calculation to obtain the vapor phase pressure of different control bodies in the current iteration cycle;

[0050] The initial parameters include initial inlet pressure, initial outlet pressure, initial inlet flow rate, initial outlet flow rate, initial inlet temperature, and initial outlet temperature.

[0051] Optionally, the calculation module is specifically used to determine the coefficients of the second structure matrix block according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle;

[0052] The LU decomposition method is used to call the matrix solving algorithm on different second structure matrix blocks for parallel calculation according to the coefficients of the second structure matrix block and the initial parameters of the large space component to obtain the liquid phase pressure of different control volumes in the current iteration cycle.

[0053] The initial parameters include initial inlet pressure, initial outlet pressure, initial inlet flow rate, initial outlet flow rate, initial inlet temperature, and initial outlet temperature.

[0054] Optionally, the parameter calculation device for the large space component of the nuclear reactor further includes:

[0055] A condition update module is used to exchange the interface parameters of the connection between the large space component and the target component obtained in the current iteration cycle to update the initial boundary conditions of the pressure solution matrix of the next iteration cycle and / or the initial boundary conditions of the mass-energy solution matrix obtained by discretizing the vapor and liquid phase conservation equations;

[0056] The target component is another component connected to the large space component, and the interface parameter includes at least one of an interface pressure parameter, an interface temperature parameter, and an interface flow parameter.

[0057] According to another aspect of the present application, a readable storage medium is provided, on which a program or instruction is stored. When the program or instruction is executed by a processor, the steps of the parameter calculation method of the large space component of the nuclear reactor are implemented.

[0058] According to another aspect of the present application, a computer device is provided, including a storage medium, a processor, and a computer program stored on the storage medium and executable on the processor, and when the processor executes the program, the steps of the parameter calculation method for the large space component of the nuclear reactor are implemented.

[0059] By means of the above technical scheme, on the one hand, the large space component itself is divided into multiple control bodies with the inlet and outlet of the system as the boundary, and the pressure solution matrix is ​​calculated in parallel at different control body nodes, effectively simulating the violent mixing and heat conduction effects that may occur in different areas of the large space component in the nuclear reactor system, reducing the difficulty of discrete solution of large matrices in the system analysis program, greatly shortening the calculation time, and improving the calculation efficiency of the system analysis program for parameters. On the other hand, the pressure solution matrix of multiple control bodies is further divided into blocks, and the vapor phase pressure and liquid phase pressure are solved respectively after the first structure matrix block and the second structure matrix block are divided, so as to fully consider that the fluids in different phases have different physical properties, and the vapor (gas) phase and liquid phase are solved and calculated independently by targeting the type, so as to avoid the mutual influence between different phases when the entire pressure solution matrix is ​​solved uniformly, which may cause the propagation and accumulation of errors, and make the system calculation more accurate.

[0060] The above description is only an overview of the technical solution of the present application. In order to more clearly understand the technical means of the present application, it can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of the present application more obvious and easy to understand, the specific implementation methods of the present application are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0061] The drawings described herein are used to provide a further understanding of the present application and constitute a part of the present application. The illustrative embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0062] Figure 1 A schematic flow chart of a method for calculating parameters of a large space component of a nuclear reactor provided in an embodiment of the present application is shown;

[0063] Figure 2 A schematic diagram of iterative calculation processes of different components provided in an embodiment of the present application is shown;

[0064] Figure 3 The structural block diagram of the parameter calculation device for the large space component of a nuclear reactor provided in an embodiment of the present application is shown. DETAILED DESCRIPTION

[0065] The present application will be described in detail below with reference to the accompanying drawings and in combination with embodiments. It should be noted that the embodiments and features in the embodiments of the present application can be combined with each other without conflict.

[0066] The embodiments of the present application are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain the present application, and cannot be interpreted as limiting the present application.

[0067] Those skilled in the art will appreciate that, unless expressly stated, the singular forms "a", "an", "said" and "the" used herein may also include plural forms. It should be further understood that the term "comprising" used in the specification of the present application refers to the presence of the features, integers, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or groups thereof. It should be understood that when we refer to an element as being "connected" or "connected" to another element, it may be directly connected or connected to the other element, or there may be intermediate elements. In addition, the "connection" or "connection" used herein may include wireless connection or wireless fusion. The term "and / or" used herein includes all or any unit and all combinations of one or more associated listed items.

[0068] Now, exemplary embodiments according to the present application will be described in more detail with reference to the accompanying drawings. However, these exemplary embodiments can be implemented in a variety of different forms and should not be interpreted as being limited to the embodiments set forth herein. It should be understood that these embodiments are provided to make the disclosure of the present application thorough and complete, and to fully convey the concepts of these exemplary embodiments to those of ordinary skill in the art.

[0069] In this embodiment, a method for calculating parameters of a large space component of a nuclear reactor is provided. Figure 1 As shown, the method includes:

[0070] Step 101, classify the nuclear reactor system to determine large space components and other components in the nuclear reactor system.

[0071] In this embodiment, the nuclear reactor system is divided into large-space components that need to occupy a large space and other components according to the volume of the components, so as to facilitate targeted simulation, analysis and calculation of the large-space components.

[0072] Specifically, the large space components may be large tank equipment such as a core water supply tank, an injection tank, or components such as a large sodium pool in a fast reactor, which are not listed one by one in this application.

[0073] Step 102: dividing the pressure solution matrix of multiple control volumes in the large space component into a first structure matrix block and a second structure matrix block.

[0074] Among them, the length of the pressure solution matrix is ​​A times the number of control volumes of the large space component, A is a multiple of 2, and the length of the first structure matrix block or the second structure matrix block is the number of control volumes of the large space component.

[0075] For example, a large space component is divided into M control volumes, and A is set to 2. Then the 2M×2M pressure solution matrix is ​​generated as follows:

[0076]

[0077] Among them, a i Represents different coefficient values ​​of the pressure solution matrix, and n represents the control body number.

[0078] The 2M×2M pressure solution matrix is ​​divided into 4 blocks, forming 2 first structure matrix blocks of length M and 2 second structure matrix blocks of length M, as follows:

[0079]

[0080] Among them, p represents pressure, k represents the time value, S gi , S ft Represent the vapor phase and liquid phase source terms respectively, i≤n.

[0081] In this embodiment, the inlet and outlet of the large space component of the system are used as the boundary, and the large space component itself is divided into multiple control bodies, so that the pressure solution matrix can be calculated separately at different control body nodes, which not only reduces the calculation scale of the single matrix of the system, but also effectively simulates the violent mixing and heat conduction effects that may occur in different areas of the large space component in the nuclear reactor system. Thereby reducing the difficulty of discrete solution of large matrices in the system analysis program, greatly shortening the calculation time, and improving the calculation efficiency of the parameters of the system analysis program. In addition, for each control body, the pressure solution matrix is ​​further divided into a first structure matrix block and a second structure matrix block, so that the pressure under different phases can be calculated separately through the first structure matrix block and the second structure matrix block, avoiding the mutual influence between different phases when the entire pressure solution matrix is ​​solved uniformly, which may cause the propagation and accumulation of errors, making the system calculation more accurate.

[0082] In one embodiment, before step 102, the parameter calculation method of the large space component of the nuclear reactor also includes: establishing vapor and liquid phase conservation equations for the large space component; discretizing the vapor and liquid phase conservation equations, determining multiple control bodies in the large space component, and obtaining the mass-energy conservation equations and momentum conservation equations corresponding to the large space component; transforming the mass-energy conservation equations and momentum conservation equations, and establishing a pressure solution matrix for multiple control bodies in the large space component.

[0083] Among them, the vapor and liquid phase conservation equations can be expanded into mass conservation equations, momentum conservation equations and / or energy conservation equations according to different output parameters; mass conservation equations: describe how the mass of each fluid changes over time. Momentum conservation equations: describe how the momentum of the fluid changes with flow and interaction. Energy conservation equations: describe the energy transfer and conversion process of the fluid.

[0084] The mass conservation equations are expressed as:

[0085] The momentum conservation equations are expressed as:

[0086] The energy conservation equations are expressed as:

[0087] Where t is time; z is the spatial coordinate; v is the fluid velocity of the two fluids; p is the pressure of the two fluids (vapor phase pressure and liquid phase pressure); g is the gravity of the two fluids plus the fluid velocity; θ is the angle between the flow channel and the horizontal position; F w is the flow resistance of the two fluids; U is the internal energy; q is the heat source energy; ρ m is the density of the mixture of the two fluids.

[0088] In this embodiment, the vapor and liquid phase conservation equations of the large space component are discretized and converted into algebraic equations (mass-energy conservation equations and momentum conservation equations). By transforming the mass-energy conservation equations and momentum conservation equations in the form of two fluids and two pressures, each control body is associated with a pressure equation group, and the pressure equation groups of multiple control bodies in the large space component are combined to form a pressure solution matrix. Thus, through discrete processing, complex continuous equations can be solved on discrete control bodies, so that mature numerical methods can be used to solve them, and then the local detailed changes of the large space component can be captured, and the difficulty of solving is reduced while calculating more accurate vapor and liquid phase distribution and energy and momentum transfer conditions, thereby improving the feasibility of calculation.

[0089] Furthermore, considering that the mixing of cold and hot fluids may cause the stratification of cold and hot fluids in the large space component, after establishing the vapor and liquid phase conservation equations of the large space component, the parameter calculation method of the large space component of the nuclear reactor also includes: updating the vapor and liquid phase conservation equations according to the equivalent thermal diffusion coefficient of the mass transfer process of the large space component.

[0090] Among them, the general expression of the updated vapor and liquid phase conservation equations of large space components is:

[0091]

[0092] Where ρ is density; v is fluid velocity; t is time; z is spatial coordinate; λ is thermal diffusion coefficient; ε is equivalent thermal diffusion coefficient, which characterizes the heat transfer effect of flow mixing enhancement in large space components of the system; S φ is the generalized source term; φ is the node control parameter of the control body, and the pressure equation group is the equation when φ is the pressure.

[0093] In this embodiment, the equivalent thermal diffusion coefficient can more accurately reflect the heat transfer and material diffusion characteristics in the mass transfer process of large space components. By introducing it into the update of the vapor and liquid phase conservation equations, the updated conservation equations can better adapt to the mass transfer characteristics of large space components, facilitate the model to better capture the actual physical process, better simulate the violent mixing and heat conduction effects of large space components of nuclear reactor systems, reduce unnecessary calculations, and improve the accuracy of calculation results.

[0094] Step 103, solving the first structure matrix block in parallel to obtain the vapor phase pressure of each control volume in the large space component.

[0095] Step 104, solving the second structural matrix block in parallel to obtain the liquid phase pressure of each control volume in the large space component.

[0096] In this embodiment, parallel solving can make full use of multi-core processors or distributed computing resources, distribute computing tasks to multiple computing units and perform them simultaneously, which can greatly shorten the computing time, especially for the calculation and simulation of large-scale large-space components. Moreover, the solution of the vapor phase pressure and liquid phase pressure of each control body is an independent process, but is interrelated in the large-space component as a whole. Through parallel solving, the interaction between each control body can be processed more accurately, which helps to improve the accuracy of the calculation and reduce the risk of calculation interruption and data loss due to software errors.

[0097] Specifically, step 103 specifically includes: determining the coefficients of the first structure matrix block according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle; using the LU decomposition method, based on the coefficients of the first structure matrix block and the initial parameters of the large space components, calling the matrix solution algorithm on different first structure matrix blocks for parallel calculation to obtain the vapor phase pressure of different control bodies in the current iteration cycle.

[0098] Similarly, step 104 specifically includes: determining the coefficients of the second structure matrix block according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle;

[0099] The LU decomposition method is adopted, and according to the coefficients of the second structure matrix block and the initial parameters of the large space component, the matrix solving algorithm is called on different second structure matrix blocks for parallel calculation to obtain the liquid phase pressure of different control volumes in the current iteration cycle.

[0100] In this embodiment, on the one hand, parallel computing can make full use of multi-core processors or distributed computing resources, distribute computing tasks to multiple computing units and perform them simultaneously, significantly shortening computing time and improving computing efficiency. On the other hand, the matrix is ​​decomposed into the product of the lower triangular matrix and the upper triangular matrix by LU decomposition, and then the linear equations are solved by forward and backward substitution, thereby reducing the error accumulation in the calculation process and helping to improve the computing performance when processing large-scale matrices.

[0101] The initial parameters include the initial inlet pressure, initial outlet pressure, initial inlet flow, initial outlet flow, initial inlet temperature, initial outlet temperature, etc. of the large space component.

[0102] Step 105, iteratively updating the interface pressure parameters of the large space component according to the vapor phase pressure and liquid phase pressure corresponding to the large space component until convergence.

[0103] In this embodiment, by continuously iteratively updating the interface pressure parameters, the interface pressure parameters gradually converge, making the calculation results closer to the actual situation, thereby achieving goals such as better fluid flow control and improved energy transfer efficiency, optimizing the performance of large-space components, and ensuring stable operation and good performance of large-space components.

[0104] The parameter calculation method for the large space component of the nuclear reactor provided in the embodiment of the present application can be applied to the terminal, can also be applied to the server, and can also be software running in the terminal or the server. In some embodiments, the terminal can be a smart phone, a tablet computer, a laptop computer, a desktop computer, etc.; the server side can be configured as an independent physical server, or a server cluster or distributed system composed of multiple physical servers, and can also be configured as a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms; the software can be an application that implements the parameter calculation method for the large space component of the nuclear reactor, etc., but is not limited to the above forms.

[0105] In one embodiment, after step 103 and step 104, the parameter calculation method of the large space component of the nuclear reactor also includes: solving the mass-energy conservation equations according to the initial boundary conditions of the mass-energy solution matrix, and the vapor phase pressure, liquid phase pressure and fluid velocity corresponding to the large space component, to obtain other node control parameters of each control body in the large space component except the pressure.

[0106] Among them, other node control parameters include temperature or cavitation fraction. The initial boundary conditions include initial conditions and boundary conditions. The initial boundary conditions are used to provide a starting point and boundary constraints for the calculation, so as to provide specific constraints and restrictions for solving the problem, so that the solution is unique. The initial boundary conditions can be reasonably set according to the functional machine control requirements of the large space equipment, and the embodiments of this application are not specifically limited.

[0107] In this embodiment, in addition to the calculated pressure and the constant fluid velocity, other node control parameters such as temperature, concentration, etc. are also crucial to fully describe the physical state of the large space component. To this end, with the initial boundary conditions as constraints, the known vapor phase pressure, liquid phase pressure and fluid velocity are substituted into the mass-energy conservation equations for solution, and these other node control parameters can be obtained, thereby gaining a deeper understanding of the thermodynamic and fluid dynamic characteristics of the large space component, and facilitating a more accurate analysis of the performance and behavior of the large space component.

[0108] In one embodiment, after obtaining the node control parameters of each control body in the large space component, the parameter calculation method of the large space component of the nuclear reactor also includes: updating the physical property parameters of the large space component through state equations and physical property functions based on the node control parameters of the control body in the large space component.

[0109] The node control parameters include at least one of pressure, temperature, fluid velocity, and cavitation fraction; the pressure includes vapor phase pressure and liquid phase pressure; and the physical property parameters include at least one of flow rate, mass, and enthalpy.

[0110] In this embodiment, the node control parameters can reflect the actual physical state of different positions inside the large space component, and then combined with the state equation and the physical property function for calculation, the physical property parameters of the large space component can be determined more accurately. In order to accurately predict the key physical processes such as temperature distribution, pressure change, fluid flow, etc., the precise and complex physical property parameters of the large space component can make the simulation of the large space component closer to the actual situation, so that the operator can understand the changes of the large space component in time, assist them in selecting the optimal design parameters, and improve the performance and efficiency of the nuclear reactor.

[0111] It is worth mentioning that for other components in the nuclear reactor system, the same method as that for large space components can also be used to calculate node control parameters (pressure, temperature, cavitation fraction, etc.). Specifically, the pressure solution matrix of multiple control bodies in other components is divided into a third structure matrix block and a fourth structure matrix block, wherein the length of the pressure solution matrix is ​​A times the number of control bodies of other components, A is a multiple of 2, and the length of the third structure matrix block or the fourth structure matrix block is the number of control bodies of other components; the third structure matrix block is solved in parallel to obtain the vapor phase pressure of each control body in other components; the fourth structure matrix block is solved in parallel to obtain the liquid phase pressure of each control body in other components; according to the vapor phase pressure and liquid phase pressure corresponding to other components, the interface pressure parameters of the other components are iteratively updated until convergence. According to the initial boundary conditions of the mass-energy solution matrix obtained by discretizing the vapor and liquid phase conservation equations of other components, and the vapor phase pressure, liquid phase pressure and fluid velocity corresponding to other components, the mass-energy conservation equation group is solved to obtain other node control parameters of each control body in other components except pressure.

[0112] In one embodiment, the parameter calculation method of the large space component of a nuclear reactor also includes: exchanging the interface parameters of the connection between the large space component and the target component obtained in the current iteration cycle to update the initial boundary conditions of the pressure solution matrix of the next iteration cycle and / or the initial boundary conditions of the mass-energy solution matrix obtained by discretizing the vapor and liquid phase conservation equations.

[0113] The target component is another component connected to the large space component, and the interface parameter includes at least one of an interface pressure parameter, an interface temperature parameter, and an interface flow parameter.

[0114] In this embodiment, the interface parameters of different interconnected components in the system are exchanged, the initial values ​​of the pressure solution matrix and the mass-energy solution matrix are updated, and the calculation of the next iteration cycle is started until the iterative calculation terminates and the output parameters are output. Thus, by exchanging the interface parameters, the large space component and the target component can better handle the interaction and coupling relationship between the various parts of the system, so that the information between the two can be transmitted and fed back more accurately, so as to continuously correct the initial boundary conditions in the iteration process, make the calculation results closer to the actual situation, and improve the accuracy of the calculation.

[0115] For example, the core water supply tank under special conditions such as ocean conditions is taken as the large space component for research, and multiple different control body nodes are divided.

[0116] Step 1: Taking the inlet and outlet of the core water supply tank as the boundary, the core water supply tank is regarded as a special process component, and the entire nuclear reactor system is divided into blocks, namely, the core water supply tank and other system components;

[0117] Step 2: Establish the steam and liquid phase mass, momentum and energy conservation equations of the core water supply tank and other system components respectively; perform a certain form of numerical discretization (such as mass sum and mass difference discretization; momentum sum and momentum difference discretization, etc., depending on the discretization form of the specific nuclear reactor system analysis program) on the steam and liquid phase mass, momentum and energy conservation equations of the core water supply tank and other system components respectively, so as to divide the core water supply tank into M control bodies, and divide other system components into N control bodies;

[0118] The general form of the conservation equations for vapor and liquid mass, momentum, and energy is as follows:

[0119]

[0120] Where ρ is density; v is fluid velocity; t is time; z is spatial coordinate; λ is thermal diffusion coefficient; ε is equivalent thermal diffusion coefficient; S φ is the generalized source term; φ is the node control parameter of the control body, including cavitation fraction, velocity, temperature or enthalpy, etc. The equivalent thermal diffusion coefficient can be determined experimentally. For the core water tank, the equivalent thermal diffusion coefficient can be related to parameters such as temperature difference and volume expansion coefficient through experimental determination; for other components of the system, the equivalent thermal diffusion coefficient is 0.

[0121] Step 3: By transforming the mass-energy conservation equation and momentum conservation equation in the form of two-fluid dual pressure, respectively establish the pressure solution matrix of the control volume divided by the core water tank and other system components, forming different 2M×2M and 2N×2N matrix equation groups;

[0122] Step 4: According to the form of the pressure solution matrix equation group, the 2M×2M or 2N×2N pressure solution matrix can be divided into 4 M×M or N×N structural matrix blocks with similar structures by adopting the idea of ​​block matrix;

[0123] Step 5: First set the initial and boundary value conditions of each pressure solution matrix, as well as the known inlet and outlet pressures, inlet and outlet temperatures, and inlet and outlet flow rates of the core water supply tank; adopt a block matrix solution algorithm, based on the general matrix LU symbolic decomposition, matrix LU numerical decomposition and equation forward and backward solution steps, and solve four M×M and N×N structure matrix blocks with similar structures in parallel to obtain the vapor / liquid phase pressure of each control body node in the core water supply tank and other components of the system, update the pressure interface parameters, iterate the calculation until convergence, and then update the coefficients to start the calculation of the next iteration cycle.

[0124] In this embodiment, step 5 specifically includes the following steps:

[0125] Step 5.1: First set the initial and boundary value conditions of the pressure matrix of the core water tank and other system components, and set the inlet and outlet pressures of the core water tank (Pin , P out ), inlet and outlet temperatures (T in , T out ), inlet and outlet flow (W in , W out ); inlet and outlet pressures of other components (P s,in , P s,out ), inlet and outlet temperatures (T s,in , T s,out ), inlet and outlet flow (W s,in , W s,out );

[0126] Step 5.2: Update the equation matrix coefficients according to the initial boundary conditions set in step 5.1. Based on the general matrix LU symbolic decomposition, matrix LU numerical decomposition and equation forward and backward solution steps, call the matrix solution algorithm on different computing nodes for parallel calculation to obtain the steam and liquid phase pressures of each control volume node of different components under the current time step;

[0127] Step 5.3: Keep the core inlet and outlet temperatures and flow rates unchanged, and exchange the pressure parameters calculated by the core makeup water tank and other system components. Other system components will flow into the core makeup water tank outlet pressure P s,out Assign a value to the core water tank inlet pressure P in ; The core water tank will make the core water tank outlet pressure P out Assigned to the inlet pressure P of other components of the system s,in .

[0128] Step 5.4: Iterative calculation, according to the pressure convergence criterion Determine whether to start calculation of the next iteration cycle.

[0129] Step 6: By associating with step 5, the node pressure and velocity parameters of each control body are obtained by calculation, and the M and N mass-energy solution matrix equation groups obtained by discretizing the mass and energy conservation equations corresponding to the core water supply tank and other system components are solved in parallel to obtain the node temperature, cavitation fraction and other parameters of each control body of the core water supply tank and other system components at the new time;

[0130] Step 7: Based on the pressure, temperature, velocity, cavitation fraction and other parameters of each control body node, other physical parameters such as flow rate, mass, enthalpy value, etc. are updated through the state equation and physical property function;

[0131] Step 8: Figure 2 As shown, the interface parameters of the core water supply tank and other system components are exchanged, the initial values ​​of the pressure solution matrix and the mass-energy solution matrix are updated, and the calculation at the next moment is started until the iterative calculation terminates and the output parameters are output.

[0132] Step 8 specifically includes: other components of the system will flow into the core water tank flow rate W s,out , Temperature T s,out Assign a value to the core water tank inlet flow rate W in , inlet temperature T in ; The core water supply tank will be the core water supply tank outlet flow W out , outlet temperature T out Assign values ​​to the inlet flow W of other components in the system s,in , inlet temperature T s,in .

[0133] In this embodiment, the entire nuclear reactor system is divided into blocks with the inlet and outlet of the system large space component as the boundary, and is divided into the system large space component and other system components. For the energy conservation equation established for the system large space component, the influence of the equation heat diffusion term is considered, and the equivalent heat diffusion coefficient of the mass transfer process is added to the heat diffusion term to characterize the heat transfer effect of the flow mixing enhancement of the system large space component. Based on the idea of ​​block matrix solution, the large pressure matrix solution matrix is ​​divided into matrix blocks with similar structural forms. Based on the general matrix LU symbolic decomposition, matrix LU numerical decomposition and equation pre-substitution back-substitution solution steps, the vapor / liquid phase pressure of each control body node of different components is obtained by parallel calculation. On different computing nodes, the mass-energy solution matrix equation group obtained by discretizing the mass and energy conservation equations corresponding to the system large space component and other system components is calculated in parallel, which can effectively reduce the difficulty of large matrix discrete solution of the system analysis program and improve the calculation efficiency of the system analysis program.

[0134] It should be noted that the size of the serial numbers of the steps in the above embodiments does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0135] Furthermore, if Figure 3 As shown, as a specific implementation of the above-mentioned parameter calculation method for large-space components of a nuclear reactor, an embodiment of the present application provides a parameter calculation device 300 for large-space components of a nuclear reactor, and the parameter calculation device 300 for large-space components of a nuclear reactor includes: a component classification module 301, a calculation module 302 and a parameter update module 303.

[0136] The component classification module 301 is used to classify the nuclear reactor system and determine the large space components and other components in the nuclear reactor system;

[0137] The calculation module 302 is used to divide the pressure solution matrix of multiple control bodies in the large space component into a first structure matrix block and a second structure matrix block, wherein the length of the pressure solution matrix is ​​A times the number of control bodies in the large space component, A is a multiple of 2, and the length of the first structure matrix block or the second structure matrix block is the number of control bodies in the large space component; and, the first structure matrix block is solved in parallel to obtain the vapor phase pressure of each control body in the large space component; and, the second structure matrix block is solved in parallel to obtain the liquid phase pressure of each control body in the large space component;

[0138] The parameter updating module 303 is used to iteratively update the interface pressure parameters of the large space component according to the vapor phase pressure and liquid phase pressure corresponding to the large space component until convergence.

[0139] Furthermore, the calculation module 302 is also used to establish the vapor and liquid phase conservation equations of the large space component; discretize the vapor and liquid phase conservation equations, determine multiple control bodies in the large space component, and obtain the mass-energy conservation equations and momentum conservation equations corresponding to the large space component; transform the mass-energy conservation equations and momentum conservation equations, and establish the pressure solution matrix of multiple control bodies in the large space component.

[0140] Furthermore, the calculation module 302 is also used to solve the mass-energy conservation equations according to the initial boundary conditions of the mass-energy solution matrix, and the vapor phase pressure, liquid phase pressure and fluid velocity corresponding to the large space component, to obtain other node control parameters of each control body in the large space component except the pressure, wherein the other node control parameters include temperature or cavitation fraction.

[0141] Further, the parameter calculation device 300 for a large space component of a nuclear reactor further includes: an equivalent module (not shown in the figure);

[0142] The equivalent module is used to update the vapor and liquid phase conservation equations according to the equivalent thermal diffusion coefficient of the mass transfer process of the large space component; the vapor and liquid phase conservation equations of the large space component are expressed as follows:

[0143]

[0144] Where ρ is density; v is fluid velocity; t is time; z is spatial coordinate; λ is thermal diffusion coefficient; ε is equivalent thermal diffusion coefficient; S φ is the generalized source term; φ is the node control parameter of the control body.

[0145] Furthermore, the calculation module 302 is also used to update the physical property parameters of the large space component through the state equation and the physical property function according to the node control parameters of the control body in the large space component; wherein the node control parameters include at least one of pressure, temperature, fluid velocity, and cavitation fraction, the pressure includes vapor phase pressure and liquid phase pressure, and the physical property parameters include at least one of flow rate, mass, and enthalpy.

[0146] Furthermore, the calculation module 302 is specifically used to determine the coefficients of the first structure matrix block according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle; using the LU decomposition method, according to the coefficients of the first structure matrix block and the initial parameters of the large space components, the matrix solution algorithm is called on different first structure matrix blocks for parallel calculation to obtain the vapor phase pressure of different control bodies in the current iteration cycle; wherein the initial parameters include initial inlet pressure, initial outlet pressure, initial inlet flow, initial outlet flow, initial inlet temperature, and initial outlet temperature.

[0147] Furthermore, the calculation module 302 is specifically used to determine the coefficients of the second structure matrix block according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle; using the LU decomposition method, according to the coefficients of the second structure matrix block and the initial parameters of the large space components, the matrix solution algorithm is called on different second structure matrix blocks for parallel calculation to obtain the liquid phase pressure of different control bodies in the current iteration cycle; wherein the initial parameters include initial inlet pressure, initial outlet pressure, initial inlet flow, initial outlet flow, initial inlet temperature, and initial outlet temperature.

[0148] Furthermore, the parameter calculation device 300 for a large space component of a nuclear reactor further includes: a condition updating module (not shown in the figure);

[0149] A condition update module is used to exchange the interface parameters at the connection between the large space component and the target component obtained in the current iteration cycle to update the initial side conditions of the pressure solution matrix of the next iteration cycle and / or the initial side conditions of the mass-energy solution matrix obtained by discretizing the vapor and liquid phase conservation equations; wherein the target component is other components connected to the large space component, and the interface parameters include at least one of interface pressure parameters, interface temperature parameters, and interface flow parameters.

[0150] The specific definition of the parameter calculation device for the large space component of the nuclear reactor can be found in the definition of the parameter calculation method for the large space component of the nuclear reactor mentioned above, which will not be repeated here. Each module in the parameter calculation device for the large space component of the nuclear reactor can be implemented in whole or in part by software, hardware and a combination thereof. The above modules can be embedded in or independent of the processor in the computer device in the form of hardware, or can be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

[0151] Based on the above Figure 1 The method shown in the embodiment of the present application accordingly provides a readable storage medium on which a computer program is stored, and when the program is executed by a processor, the above-mentioned Figure 1 The parameter calculation method of the large space components of the nuclear reactor is shown.

[0152] Based on this understanding, the technical solution of the present application can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (which can be a CD-ROM, a USB flash drive, a mobile hard disk, etc.), and includes a number of instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each implementation scenario of the present application.

[0153] Based on the above Figure 1 The method shown, and Figure 3 In order to achieve the above-mentioned purpose, the embodiment of the present application further provides a computer device, which can be a personal computer, a server, a network device, etc. The computer device includes a storage medium and a processor; the storage medium is used to store a computer program; the processor is used to execute the computer program to achieve the above-mentioned Figure 1 The parameter calculation method of the large space components of the nuclear reactor is shown.

[0154] Optionally, the computer device may further include a user interface, a network interface, a camera, a radio frequency (RF) circuit, a sensor, an audio circuit, a WI-FI module, etc. The user interface may include a display, an input unit such as a keyboard, etc., and the optional user interface may also include a USB interface, a card reader interface, etc. The network interface may optionally include a standard wired interface, a wireless interface (such as a Bluetooth interface, a WI-FI interface), etc.

[0155] Those skilled in the art will appreciate that the computer device structure provided in this embodiment does not limit the computer device, and may include more or fewer components, or a combination of certain components, or different component arrangements.

[0156] The storage medium may also include an operating system and a network communication module. The operating system is a program that manages and saves the hardware and software resources of the computer device, and supports the operation of information processing programs and other software and / or programs. The network communication module is used to realize communication between the components inside the storage medium, and communication with other hardware and software in the physical device.

[0157] Through the description of the above implementation methods, those skilled in the art can clearly understand that the present application can be implemented by means of software plus a necessary general hardware platform, or by hardware.

[0158] Those skilled in the art will appreciate that the accompanying drawings are only schematic diagrams of a preferred implementation scenario, and the modules or processes in the accompanying drawings are not necessarily necessary for implementing the present application. Those skilled in the art will appreciate that the modules in the devices in the implementation scenario can be distributed in the devices of the implementation scenario according to the description of the implementation scenario, or can be changed accordingly and located in one or more devices different from the present implementation scenario. The modules of the above-mentioned implementation scenario can be combined into one module, or can be further split into multiple submodules.

[0159] The above serial numbers of this application are only for description and do not represent the advantages and disadvantages of the implementation scenarios. The above disclosure is only a few specific implementation scenarios of this application, but this application is not limited to them, and any changes that can be thought of by technicians in this field should fall within the scope of protection of this application.

Claims

1. A method for calculating parameters of large space components of a nuclear reactor, characterized in that: The method comprises: Classifying a nuclear reactor system and identifying large space components and other components in the nuclear reactor system; Divide the pressure solution matrix of the multiple control bodies in the large space component into a first structure matrix block and a second structure matrix block, wherein the length of the pressure solution matrix is ​​A times the number of control bodies of the large space component, A is a multiple of 2, and the length of the first structure matrix block or the second structure matrix block is the number of control bodies of the large space component; Solving the first structural matrix blocks in parallel to obtain the vapor phase pressure of each control volume in the large space component; Solving the second structural matrix blocks in parallel to obtain the liquid phase pressure of each control volume in the large space component; According to the vapor phase pressure and liquid phase pressure corresponding to the large space component, the interface pressure parameter of the large space component is iteratively updated until convergence.

2. The method for calculating parameters of a large space component of a nuclear reactor according to claim 1, characterized in that: The method further comprises: Establishing the vapor and liquid phase conservation equations of the large space component; Discretizing the vapor and liquid phase conservation equations, determining a plurality of control bodies in the large space component, and obtaining a group of mass-energy conservation equations and a group of momentum conservation equations corresponding to the large space component; The mass-energy conservation equations and the momentum conservation equations are transformed to establish a pressure solution matrix for multiple control volumes in the large space component.

3. The method for calculating parameters of a large space component of a nuclear reactor according to claim 2, characterized in that: The method further comprises: According to the initial boundary conditions of the mass-energy solution matrix, and the vapor phase pressure, liquid phase pressure and fluid velocity corresponding to the large space component, the mass-energy conservation equations are solved to obtain other node control parameters of each control body in the large space component except the pressure, wherein the other node control parameters include temperature or cavitation fraction.

4. The method for calculating parameters of a large space component of a nuclear reactor according to claim 1, characterized in that: The method further comprises: Updating the vapor and liquid phase conservation equations according to the equivalent thermal diffusion coefficient of the mass transfer process of the large space component; The vapor and liquid phase conservation equations of the large space component are expressed as follows: Where ρ is density; v is fluid velocity; t is time; z is spatial coordinate; λ is thermal diffusion coefficient; ε is equivalent thermal diffusion coefficient; S φ is the generalized source term; φ is the node control parameter of the control body.

5. The method for calculating parameters of a large space component of a nuclear reactor according to claim 1, characterized in that: The method further comprises: According to the node control parameters of the control body in the large space component, the physical property parameters of the large space component are updated through the state equation and the physical property function; The node control parameters include at least one of pressure, temperature, fluid velocity, and cavitation fraction; the pressure includes vapor phase pressure and liquid phase pressure; and the physical property parameters include at least one of flow rate, mass, and enthalpy.

6. The method for calculating parameters of a large space component of a nuclear reactor according to claim 1, characterized in that: The parallel solution of the first structural matrix block to obtain the vapor phase pressure of each control volume in the large space component includes: Determining coefficients of the first structural matrix block according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle; Using the LU decomposition method, according to the coefficients of the first structural matrix block and the initial parameters of the large space component, the matrix solving algorithm is called on different first structural matrix blocks to perform parallel calculations to obtain the vapor phase pressures of different control volumes in the current iteration cycle; The parallel solution of the second structural matrix block to obtain the liquid phase pressure of each control volume in the large space component includes: Determining coefficients of the second structural matrix block according to the initial boundary conditions of the pressure solution matrix in the current iteration cycle; Using the LU decomposition method, according to the coefficients of the second structural matrix block and the initial parameters of the large space component, the matrix solving algorithm is called on different second structural matrix blocks to perform parallel calculations to obtain the liquid phase pressures of different control volumes in the current iteration cycle; The initial parameters include initial inlet pressure, initial outlet pressure, initial inlet flow, initial outlet flow, initial inlet temperature, and initial outlet temperature.

7. The method for calculating parameters of a large space component of a nuclear reactor according to claim 2 or 6, characterized in that: The method further comprises: Exchange the interface parameters of the connection between the large space component and the target component obtained in the current iteration cycle to update the initial boundary conditions of the pressure solution matrix and / or the initial boundary conditions of the mass-energy solution matrix obtained by discretizing the vapor and liquid phase conservation equations in the next iteration cycle; Wherein, the target component is the other component connected to the large space component, and the interface parameter includes at least one of an interface pressure parameter, an interface temperature parameter, and an interface flow parameter.

8. A parameter calculation device for a large space component of a nuclear reactor, characterized in that: The device comprises: A component classification module, used for classifying the nuclear reactor system and determining large space components and other components in the nuclear reactor system; A calculation module, used for dividing the pressure solution matrix of multiple control bodies in the large space component into a first structure matrix block and a second structure matrix block, wherein the length of the pressure solution matrix is ​​A times the number of control bodies of the large space component, A is a multiple of 2, and the length of the first structure matrix block or the second structure matrix block is the number of control bodies of the large space component; and Solving the first structural matrix blocks in parallel to obtain the vapor phase pressure of each control volume in the large space component; and, Solving the second structural matrix blocks in parallel to obtain the liquid phase pressure of each control volume in the large space component; The parameter updating module is used to iteratively update the interface pressure parameters of the large space component according to the vapor phase pressure and liquid phase pressure corresponding to the large space component until convergence.

9. A readable storage medium having a program or instruction stored thereon, characterized in that: When the program or instruction is executed by a processor, the steps of the method for calculating parameters of a large-space component of a nuclear reactor as claimed in any one of claims 1 to 7 are implemented.

10. A computer device comprising a storage medium, a processor, and a computer program stored in the storage medium and executable on the processor, characterized in that: When the processor executes the program, the method for calculating parameters of a large-space component of a nuclear reactor according to any one of claims 1 to 7 is implemented.

Images