Coupling solving method, program, equipment and storage medium suitable for performance analysis of integrated steam-water separation device of nuclear energy system

Abstract

The invention belongs to the technical field of three-dimensional refined numerical calculation of an integrated steam-water separator in a nuclear energy system, and particularly relates to a coupling solution method, program and equipment suitable for performance analysis of an integrated steam-water separator of the nuclear energy system and a storage medium. According to the method, data of an outlet of the rotary vane separator and data of an inlet of the corrugated plate dryer are coupled and serve as boundary conditions for each other, iteration is conducted at the same time, and in the iteration process, a source item elimination model is added to a hydrophobic area of the rotary vane separator to simulate the liquid seal phenomenon of a hydrophobic opening; and adding the critical Weber number model into the corrugated plate dryer, judging the liquid drop state, simulating the secondary carrying phenomenon in the corrugated plate dryer, obtaining the real flow condition of the corrugated plate dryer, and finally obtaining the flow condition in the integrated steam-water separator. According to the invention, the internal hydrophobic liquid seal phenomenon of the rotary vane separator can be finely simulated, and the inflection point of the secondary carrying speed of the corrugated plate dryer can be accurately predicted.

Description

Technical Field

[0001] The present invention belongs to the technical field of three-dimensional refined numerical calculation of integrated steam-water separators in nuclear energy systems, and specifically relates to a coupling solution method, program, equipment and storage medium suitable for performance analysis of integrated steam-water separators in nuclear energy systems. Background Art

[0002] As the core component of the dehumidification equipment in the secondary circuit, the separation performance of the steam-water separator directly affects the safety, reliability and economy of the nuclear power plant. The steam-water separator is divided into a coarse separator and a fine separator. The coarse separator is generally composed of a rotary vane separator, which is used to separate large-size droplets and separate 90% of the liquid phase in the steam; the fine separator is generally composed of a corrugated plate dryer, which is used to separate small-size droplets and separate the droplets that have not been separated in the coarse separator. The outlet humidity cannot be higher than 0.25%. The existing steam-water separation devices are designed separately for the coarse separator and the fine separator, and then combined for testing. In recent years, with the vigorous development of computer and computational fluid dynamics (CFD) technology, three-dimensional computational fluid dynamics solvers are increasingly used in two-phase flow analysis due to their many advantages such as high precision, good repeatability, low cost requirements, short implementation cycle and no safety concerns. CFD technology is also widely used in steam-water separators. However, due to the complex internal separation mechanism of rotary vane separators and corrugated plate dryers, CFD technology is difficult to reproduce some real situations, such as: the liquid seal problem of the rotary vane separator drain port and the secondary carryover phenomenon of the corrugated plate dryer.

[0003] At present, nuclear energy technology is developing in the direction of miniaturization and integration. The gas-water separation device is also developing in this direction, but the existing research on the gas-water separation device is to conduct numerical research on the rotary vane separator and the corrugated plate dryer separately, and the internal flow of the integrated gas-water separation device cannot be obtained. The difficulty lies in that the calculation models applicable to the rotary vane separator and the corrugated plate dryer cannot be used in the same fluid domain. If only one of the calculation models is used to perform numerical simulation on the integrated gas-water separator, the error is unacceptable. In view of the above problems, different applicable models are developed for the liquid seal problem of the drain port of the rotary vane separator and the secondary carryover problem of the corrugated plate dryer. Then, through data coupling, the calculations of the rotary vane separator and the corrugated plate dryer are coupled, and then the internal flow of the integrated gas-water separator is analyzed, and the three-dimensional refined numerical simulation of the integrated gas-water separator is realized. Summary of the invention

[0004] The object of the present invention is to provide a coupling solution method, program, equipment and storage medium suitable for performance analysis of an integrated steam-water separation device in a nuclear energy system.

[0005] A coupling solution method suitable for performance analysis of an integrated steam-water separation device in a nuclear energy system comprises the following steps:

[0006] Step 1: Model and mesh the rotary blade separator and corrugated plate dryer; the hydrophobic area of ​​the rotary blade separator is meshed separately, and the corrugated plate dryer is meshed by partitions to separate the baffle, uniform air hole, and corrugated plate areas;

[0007] Step 2: Use the integrated coupling calculation method to couple and iterate the rotary vane separator and the corrugated plate dryer, use the outlet data of the rotary vane separator as the inlet boundary condition of the corrugated plate dryer, and use the inlet data of the corrugated plate dryer as the outlet boundary condition of the rotary vane separator;

[0008] First, the rotary vane separator is calculated. The gas mass flow rate and liquid mass flow rate at the inlet boundary of the rotary vane separator are input, and the outlet pressure P1 of the rotary vane separator is input. The calculation is a steady-state calculation. In order to simulate the phenomenon of liquid sealing at the drain port, the source term elimination model is added to the hydrophobic area of ​​the rotary vane separator. The liquid phase is eliminated after entering the hydrophobic area, while the gas phase is not affected. The gas phase flows out from the outlet of the rotary vane separator, thereby achieving the model drain port liquid sealing effect. After the calculation converges, the gas mass flow rate and liquid mass flow rate at the outlet of the rotary vane separator are obtained.

[0009] The gas mass flow rate and liquid mass flow rate at the outlet of the rotary blade separator are used as input to calculate the corrugated plate dryer. In order to simulate the secondary carryover phenomenon inside the corrugated plate dryer, the critical Weber number model is added to the wall conditions of the corrugated plate area; if the Weber number of the droplet is less than the critical Weber number, the droplet is absorbed and no subsequent calculation is performed; if the Weber number of the droplet is greater than or equal to the critical Weber number, the droplet is given a momentum loss coefficient for rebound and continues to participate in subsequent calculations until the droplet is absorbed or leaves the calculation domain; after the calculation converges, the pressure P2 at the inlet of the corrugated plate dryer and the liquid mass flow rate at the outlet are obtained;

[0010] Step 3: Compare the pressure P2 at the inlet of the corrugated plate dryer with the outlet pressure P1 of the rotary blade separator. If the error is less than 1%, the coupling is considered successful; otherwise, take P2 as the outlet pressure of the rotary blade separator, that is, set P1=P2, and return to step 2 to recalculate.

[0011] Furthermore, in the process of model construction and meshing of the rotary vane separator in step 1, the fluid domain surrounded by the two blades is meshed as a whole, firstly a line mesh is generated for the envelope line of the blade, then a surface mesh is generated based on the line mesh of the blade envelope line, and then the two blade surface meshes are converted into volume meshes, thereby obtaining a mesh model of the blade area of ​​the rotary vane separator through the line mesh-surface mesh-volume mesh conversion generation technology.

[0012] Furthermore, in step 1, the hydrophobic region of the rotary vane separator is grid-divided separately, and the length of the hydrophobic region accounts for 20% to 30% of the length from the blade to the hydrophobic port.

[0013] Furthermore, in the step 1, partitioned grid division is adopted for the corrugated plate dryer to separate the three regions of the baffle, the uniform air hole and the corrugated plate, and then the model of each region is cut with 1 / 60 to 1 / 10 of the whole model with the Z axis as the symmetry axis as the unit model of each region, and the structured grid division technology is adopted for each unit model, and each unit grid is arranged in a periodic array, and finally assembled by using the Interface technology to obtain the grid model of the corrugated plate dryer.

[0014] Furthermore, in step 2, the source term elimination model is added to the hydrophobic area of ​​the rotary vane separator, specifically:

[0015] Source cancellation model in a rotary vane separator:

[0016]

[0017] Among them, ρ and C l is the density and volume fraction of the gas-liquid two-phase mixture in the hydrophobic area grid unit of the rotary blade separator; t is the time step;

[0018] The source cancellation model is applied to the mass equation of the hydrophobic region mesh of the rotary vane separator:

[0019]

[0020] in, represents the Hamiltonian operator, V represents the velocity vector;

[0021] The liquid volume fraction in the hydrophobic area is obtained by solving the mass equation and momentum equation of the hydrophobic area grid of the rotary vane separator. The gas mass flow rate and liquid mass flow rate at the outlet of the rotary vane separator are obtained by further solving the mass equation and momentum equation of the entire grid of the rotary vane separator.

[0022] Furthermore, in the calculation process of the corrugated plate dryer in step 2, the corrugated plate dryer adopts a discrete term model, and the diameter distribution of the droplet particles at the inlet of the corrugated plate dryer adopts the Rosin-Rammler distribution function:

[0023]

[0024] in, is the droplet size constant, n is the size distribution parameter; d is the diameter of the droplet particle, and each grid at the inlet plane of the corrugated plate dryer contains a diameter distribution that follows Yd droplet particles;

[0025] Set the volume of group N to V 粒子 The particles are unified into a volume V 包 The package is calculated, that is, NV 粒子 =V 包 ;in, m is the mass of the droplet particles; Ml2 is the liquid mass flow rate at the outlet of the rotary blade separator; t is the time step.

[0026] Furthermore, the calculation method of the critical Weber number of the droplet in step 2 is:

[0027]

[0028] Among them, u l represents the average velocity of the liquid film; u g represents the gas phase velocity, ρ l represents the liquid density; ρ g represents the gas phase density, l represents the plate spacing, σ represents the surface tension of the liquid film; u l represents the dynamic viscosity of the liquid phase, u g It represents the dynamic viscosity of the gas phase.

[0029] A computer device / equipment / system includes a memory, a processor and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the above-mentioned coupling solution method suitable for performance analysis of an integrated steam-water separation device of a nuclear energy system.

[0030] A computer-readable storage medium stores a computer program / instruction, which, when executed by a processor, implements the steps of the coupling solution method for performance analysis of an integrated steam-water separation device in a nuclear energy system.

[0031] A computer program product includes a computer program / instruction, which, when executed by a processor, implements the steps of the coupling solution method for performance analysis of an integrated steam-water separation device in a nuclear energy system.

[0032] The beneficial effects of the present invention are:

[0033] The present invention couples the data of the outlet of the rotary vane separator with the data of the inlet of the corrugated plate dryer. The two serve as boundary conditions for each other and are iterated at the same time. During the iteration process, a source term elimination model is added to the hydrophobic area of ​​the rotary vane separator to simulate the hydrophobic port liquid seal phenomenon; a critical Weber number model is added to the corrugated plate dryer to judge the droplet state, simulate the secondary carryover phenomenon inside the corrugated plate dryer, obtain the real flow condition of the corrugated plate dryer, and finally obtain the flow condition inside the integrated steam-water separator. The present invention can perform a refined simulation of the hydrophobic liquid seal phenomenon inside the rotary vane separator and accurately predict the inflection point of the secondary carryover velocity of the corrugated plate dryer. When simulating the integrated steam-water separator, the two calculation models are calculated in two fluid domains by data coupling, realizing the refined simulation of the integrated steam-water separator, greatly improving the calculation efficiency and calculation accuracy of the three-dimensional refined CFD numerical simulation in the application fields of rotary vane separators, corrugated plate dryers and integrated steam-water separators. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] Figure 1 Schematic diagram of the rotary blade separator grid.

[0035] Figure 2 Schematic diagram of grid partitioning of corrugated plate dryer.

[0036] Figure 3 Schematic diagram of the grid of each area of ​​the corrugated plate dryer.

[0037] Figure 4 It is the microelement equilibrium diagram of the curved liquid film.

[0038] Figure 5 It is the overall flow chart of the present invention. DETAILED DESCRIPTION

[0039] The present invention is further described below in conjunction with the accompanying drawings.

[0040] In order to solve the problems in the prior art and take into account the accuracy and rapidity of three-dimensional refined CFD numerical calculations, the present invention provides a coupling solution method suitable for the performance analysis of an integrated steam-water separation device in a nuclear energy system, which divides the integrated steam-water separator into two calculation domains, namely a rotary vane separator and a corrugated plate dryer. The two calculation domains are meshed using high-quality meshing technology, and then the source term elimination model is loaded into the rotary vane separator to simulate the real phenomenon of the hydrophobic liquid seal of the rotary vane separator. At the same time, the critical Weber number model is loaded into the corrugated plate dryer to simulate the secondary carryover phenomenon of the corrugated plate dryer to achieve three-dimensional refined simulation. Finally, the data coupling of the outlet boundary of the rotary vane separator and the inlet boundary of the corrugated plate dryer is achieved through an integrated coupling calculation method, and a converged solution is obtained in the CFD numerical iterative solution process.

[0041] The present invention adopts an integrated coupling calculation method to couple and iterate the rotary vane separator and the corrugated plate dryer, taking the outlet data of the rotary vane separator as the inlet boundary condition of the corrugated plate dryer, and taking the inlet data of the corrugated plate dryer as the outlet boundary condition of the rotary vane separator. When the calculation converges, the final calculation result of the integrated steam-water separator is obtained. The iterative calculation process is shown in the attached figure. Figure 5 In the coupled calculation process, it is necessary to first construct the model and mesh the rotary blade separator and corrugated plate dryer, which specifically includes the following steps:

[0042] Step 1: Model and mesh the rotary blade separator and corrugated plate dryer; the hydrophobic area of ​​the rotary blade separator is meshed separately, and the corrugated plate dryer is meshed by partitions to separate the baffle, uniform air hole, and corrugated plate areas;

[0043] Since the torsion angle of the blade area of ​​the rotary vane separator is large, the general mesh generation technology cannot be applied to the rotary vane separator. Therefore, the fluid domain surrounded by the two blades is meshed as a whole. First, a line mesh is generated for the envelope of the blade, and then a surface mesh is generated based on the line mesh of the blade envelope. Then, the two blade surface meshes are converted into volume meshes. In this way, the mesh model of the blade area of ​​the rotary vane separator is obtained through the line mesh-surface mesh-volume mesh conversion generation technology;

[0044] To simulate the hydrophobic liquid seal of the rotary blade separator, a source elimination model needs to be added. The source elimination model needs to be added to the hydrophobic area, so the hydrophobic area needs to be meshed separately. Since the size of the hydrophobic area has a great influence on the calculation results, it is calculated that the length of the hydrophobic area is 20% to 30% of the length from the blade to the drain port, which is the best. The mesh model of the above operation is shown in the attached figure. Figure 1 As shown;

[0045] Due to the complex structure of the corrugated plate dryer, there are baffles, uniform air holes, corrugated plates and other areas inside. The general grid generation technology has the disadvantages of low calculation efficiency and large grid volume. Therefore, the above three areas are separated by partitioning, as shown in the attached figure. Figure 2 Then, the model of each region is cut with 1 / 60 to 1 / 10 of the whole model with the Z axis as the symmetry axis as the unit model of each region. The structured grid division technology is used for each unit model. The grid model of each unit region is as shown in the attached figure. Figure 3 , and then arrange each unit grid into a periodic array, and finally assemble them using Interface technology to obtain the grid model of the corrugated plate dryer;

[0046] Step 2: Use the integrated coupling calculation method to couple and iterate the rotary vane separator and the corrugated plate dryer, use the outlet data of the rotary vane separator as the inlet boundary condition of the corrugated plate dryer, and use the inlet data of the corrugated plate dryer as the outlet boundary condition of the rotary vane separator;

[0047] First, the rotary vane separator is calculated. The inlet boundary gas mass flow rate Mg1 and liquid mass flow rate Ml1 of the rotary vane separator are input, and the outlet pressure P1 of the rotary vane separator is input. The calculation is a steady-state calculation. In order to simulate the liquid sealing phenomenon of the drain port, the source term elimination model is added to the hydrophobic area of ​​the rotary vane separator through the user self-compiled function (UDF). The liquid phase is eliminated after entering the hydrophobic area, while the gas phase is not affected. The gas phase flows out from the outlet of the rotary vane separator, thereby achieving the model drain port liquid sealing effect; after the calculation converges, the gas mass flow rate and liquid mass flow rate at the outlet of the rotary vane separator are obtained;

[0048] The linearized expression of the source term is:

[0049] S=A+BΦ (1)

[0050] Where: S is the source term, A and B are constants, and Φ is the independent variable of the relevant source term;

[0051] For the rotary vane separator, Φ in formula (1) is the volume fraction of the liquid phase C l , the liquid volume fraction in the drain port area should be 0, so A is equal to 0. As the liquid phase continues to move to the drain port, the liquid content continues to decrease, so B is related to the time step t. In addition, for a grid unit, the liquid content is equal to the product of the density ρ of the mixture and the volume fraction C l , the formula (1) is rewritten as:

[0052]

[0053] Formula (2) is the source elimination model in the rotary blade separator. Since the calculation is a steady-state calculation, the size of the time step t in Formula (2) needs to be given. By comparing with the experimental data, t can be taken as 0.01s.

[0054] The source elimination model is loaded onto the mass equation of the mesh in the drain outlet area, as shown in equation (3):

[0055]

[0056] in, represents the Hamiltonian operator, V represents the velocity vector;

[0057] By solving the mass equation and momentum equation of the hydrophobic area grid of the rotary vane separator, the liquid volume fraction of the hydrophobic area is obtained. The mass equation and momentum equation of the entire grid of the rotary vane separator are further solved to obtain the gas mass flow rate Mg2 and liquid mass flow rate Ml2 at the outlet of the rotary vane separator.

[0058] Then the calculation of the corrugated plate dryer is carried out, and the gas mass flow rate Mg2 and liquid mass flow rate Ml2 at the outlet of the rotary blade separator are taken as input. In order to simulate the secondary carryover phenomenon inside the corrugated plate dryer, the critical Weber number model is added to the wall conditions of the corrugated plate area.

[0059] Since the corrugated plate dryer adopts a discrete term model (DPM), the liquid mass flow rate Ml2 cannot be directly input into the DPM model of the corrugated plate dryer, but needs to be input into the inlet of the corrugated plate dryer using UDF.

[0060] The particle size distribution at the inlet adopts the Rosin-Rammler distribution function. The Rosin-Rammler distribution function is based on the assumption that there is an exponential relationship between the droplet diameter d and the mass fraction of droplets with a diameter greater than d. The mathematical expression is:

[0061]

[0062] in, is the droplet size constant, n is the size distribution parameter; d is the diameter of the droplet particle, and each grid at the inlet plane of the corrugated plate dryer contains a diameter distribution that follows Y d droplet particles; through experimental measurement, we can get The size of n, and Y d The values ​​can be 0.95, 0.7, 0.5, 0.4, 0.2, and 0.05 respectively, and the diameter d of the particle can be obtained by formula (4).

[0063] When using integrated coupling calculation, in order to use UDF in the DPM model expression (4), the volume of N groups is V 粒子 The particles are unified into a volume V 包 The package is calculated, that is, NV 粒子 =V 包 ;in, m is the mass of the droplet particles; Ml2 is the liquid mass flow rate at the outlet of the rotary blade separator; t is the time step.

[0064] In order to simulate the secondary carryover phenomenon inside the corrugated plate dryer, the critical Weber number model is added to the corrugated plate wall condition through UDF. When the droplet hits the wall, it will be judged whether the droplet is absorbed according to the size of the Weber number. The mathematical expression of the critical Weber number model is as follows:

[0065] The critical Weber number is defined as:

[0066]

[0067] Where: l is the liquid density, u l is the average velocity of the liquid film, h is the thickness of the liquid film, and σ is the surface tension of the liquid film;

[0068] Step 4.2. Since the first place where the liquid film breaks is in contact with the airflow, when it breaks, It can be considered to be equal to the air flow velocity. In formula (5), the thickness of the liquid film must also be known, as shown in the following figure. Figure 4 As shown, the force balance analysis of the curved liquid film microelement is performed:

[0069]

[0070] Where: f is the external force on the fluid element, β is the angle corresponding to the water film, and u is the flow velocity of the liquid film.

[0071] Ignore the small terms in the above formula and integrate in the direction of liquid film thickness to get:

[0072]

[0073] When the liquid film breaks in this microelement, F = 0. So: Right now

[0074] Considering the steady state of the liquid film, the shear stress of the airflow on the liquid film surface is equal to the shear stress of the wall on the liquid film. Considering that the thickness of the liquid film accounts for a small proportion of the channel width, the shear stress of the airflow on the wall of the corrugated plate is regarded as the shear stress of the airflow on the liquid film, so:

[0075]

[0076] From equations (5)-(8), we can deduce:

[0077]

[0078] Among them, u g represents the gas phase velocity, ρ g represents the gas phase density, l represents the plate spacing; μ l Indicates the dynamic viscosity of the liquid phase, μ g It represents the dynamic viscosity of the gas phase.

[0079] Substitute the relevant physical property parameters in the calculation condition into equation (9) to obtain the critical Weber number. The critical Weber number is incorporated into the wall condition through a user-defined function (UDF) and the following judgment is made: if the Weber number of the droplet is less than the critical Weber number (We <We c), then the droplet is absorbed and no subsequent calculations are performed on the droplet; if the Weber number of the droplet is greater than or equal to the critical Weber number (We ≥ We c ), the droplet is given a momentum loss coefficient ε to rebound, and the droplet continues to participate in subsequent calculations until the droplet is absorbed or leaves the calculation domain; after the calculation converges, the pressure P2 at the inlet of the corrugated plate dryer and the liquid mass flow rate Ml3 at the outlet are obtained;

[0080] Step 3: Compare the pressure P2 at the inlet of the corrugated plate dryer with the outlet pressure P1 of the rotary blade separator. If the error is less than 1%, the coupling is considered successful; otherwise, take P2 as the outlet pressure of the rotary blade separator, that is, set P1=P2, and return to step 2 to recalculate.

[0081] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A coupling solution method suitable for performance analysis of integrated steam-water separation device in nuclear energy system, characterized in that: The following steps are involved: Step 1: Model and mesh the rotary blade separator and corrugated plate dryer; the hydrophobic area of ​​the rotary blade separator is meshed separately, and the corrugated plate dryer is meshed by partitions to separate the baffle, uniform air hole, and corrugated plate areas; Step 2: Use the integrated coupling calculation method to couple and iterate the rotary vane separator and the corrugated plate dryer, use the outlet data of the rotary vane separator as the inlet boundary condition of the corrugated plate dryer, and use the inlet data of the corrugated plate dryer as the outlet boundary condition of the rotary vane separator; First, the rotary vane separator is calculated. The gas mass flow rate and liquid mass flow rate at the inlet boundary of the rotary vane separator are input, and the outlet pressure P1 of the rotary vane separator is input. The calculation is a steady-state calculation. In order to simulate the liquid sealing phenomenon of the drain port, the source term elimination model is added to the hydrophobic area of ​​the rotary vane separator. The liquid phase is eliminated after entering the hydrophobic area, while the gas phase is not affected. The gas phase flows out from the outlet of the rotary vane separator, thereby achieving the model drain port liquid sealing effect. After the calculation converges, the gas mass flow rate and liquid mass flow rate at the outlet of the rotary blade separator are obtained; The gas mass flow rate and liquid mass flow rate at the outlet of the rotary blade separator are used as input to calculate the corrugated plate dryer. In order to simulate the secondary carryover phenomenon inside the corrugated plate dryer, the critical Weber number model is added to the wall conditions of the corrugated plate area; if the Weber number of the droplet is less than the critical Weber number, the droplet is absorbed and no subsequent calculation is performed; if the Weber number of the droplet is greater than or equal to the critical Weber number, the droplet is given a momentum loss coefficient for rebound and continues to participate in subsequent calculations until the droplet is absorbed or leaves the calculation domain; after the calculation converges, the pressure P2 at the inlet of the corrugated plate dryer and the liquid mass flow rate at the outlet are obtained; Step 3: Compare the pressure P2 at the inlet of the corrugated plate dryer with the outlet pressure P1 of the rotary blade separator. If the error is less than 1%, the coupling is considered successful; otherwise, take P2 as the outlet pressure of the rotary blade separator, that is, set P1=P2, and return to step 2 to recalculate.

2. A coupled solution method for performance analysis of an integrated steam-water separation device for a nuclear energy system according to claim 1, characterized in that: In the process of model building and meshing the rotary vane separator in step 1, the fluid domain surrounded by the two blades is meshed as a whole, firstly a line mesh is generated for the envelope line of the blade, then a surface mesh is generated based on the line mesh of the blade envelope line, and then the two blade surface meshes are converted into volume meshes, thereby obtaining a mesh model of the blade area of ​​the rotary vane separator through the line mesh-surface mesh-volume mesh conversion generation technology.

3. A coupled solution method for performance analysis of an integrated steam-water separation device in a nuclear energy system according to claim 1, characterized in that: In the step 1, the hydrophobic region of the rotary vane separator is grid-divided separately, and the length of the hydrophobic region accounts for 20% to 30% of the length from the blade to the hydrophobic port.

4. A coupled solution method for performance analysis of an integrated steam-water separation device in a nuclear energy system according to claim 1, characterized in that: In the step 1, the corrugated plate dryer is divided into partitioned grids to separate the three regions of the baffle, the uniform air hole, and the corrugated plate. Then, the model of each region is cut with 1 / 60 to 1 / 10 of the whole model with the Z axis as the symmetry axis as the unit model of each region. The structured grid division technology is used for each unit model, and each unit grid is arranged in a periodic array. Finally, the interface technology is used to assemble them to obtain the grid model of the corrugated plate dryer.

5. A coupled solution method for performance analysis of an integrated steam-water separation device in a nuclear energy system according to claim 1, characterized in that: In step 2, the source term elimination model is added to the hydrophobic area of ​​the rotary vane separator, specifically: Source cancellation model in a rotary vane separator: Among them, ρ and C l is the density and volume fraction of the gas-liquid two-phase mixture in the hydrophobic area grid unit of the rotary blade separator; t is the time step; The source cancellation model is applied to the mass equation of the hydrophobic region mesh of the rotary vane separator: in, represents the Hamiltonian operator, V represents the velocity vector; The liquid volume fraction in the hydrophobic area is obtained by solving the mass equation and momentum equation of the hydrophobic area grid of the rotary vane separator. The gas mass flow rate and liquid mass flow rate at the outlet of the rotary vane separator are obtained by further solving the mass equation and momentum equation of the entire grid of the rotary vane separator.

6. A coupled solution method for performance analysis of an integrated steam-water separation device in a nuclear energy system according to claim 1, characterized in that: In the calculation process of the corrugated plate dryer in step 2, the corrugated plate dryer adopts a discrete term model, and the diameter distribution of the droplet particles at the inlet of the corrugated plate dryer adopts the Rosin-Rammler distribution function: in, is the droplet size constant, n is the size distribution parameter; d is the diameter of the droplet particle, and each grid at the inlet plane of the corrugated plate dryer contains a diameter distribution that follows Y d droplet particles; Set the volume of group N to V 粒子 The particles are unified into a volume V 包 The package is calculated, that is, NV 粒子 =V 包 ;in, m is the mass of the droplet particles; Ml2 is the liquid mass flow rate at the outlet of the rotary blade separator; t is the time step.

7. A coupled solution method for performance analysis of an integrated steam-water separation device in a nuclear energy system according to claim 1, characterized in that: The calculation method of the critical Weber number of the droplet in step 2 is: Among them, u l represents the average velocity of the liquid film; u g represents the gas phase velocity, ρ l represents the liquid density; ρ g represents the gas phase density, l represents the plate spacing, σ represents the surface tension of the liquid film; μ l Indicates the dynamic viscosity of the liquid phase, μ g It represents the dynamic viscosity of the gas phase.

8. A computer device / equipment / system comprising a memory, a processor and a computer program stored in the memory, characterized in that: The processor executes the computer program to implement the steps of the method according to any one of claims 1 to 7.

9. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instructions are executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer program product comprising a computer program / instructions, characterized in that: When the computer program / instructions are executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images