CN113704959B - Simulation method and system for heat exchanger type passive containment cooling system - Google Patents
Simulation method and system for heat exchanger type passive containment cooling system Download PDFInfo
- Publication number
- CN113704959B CN113704959B CN202110785696.3A CN202110785696A CN113704959B CN 113704959 B CN113704959 B CN 113704959B CN 202110785696 A CN202110785696 A CN 202110785696A CN 113704959 B CN113704959 B CN 113704959B
- Authority
- CN
- China
- Prior art keywords
- tube
- heat exchange
- heat
- exchange coefficient
- coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000001816 cooling Methods 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 title claims abstract description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 26
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000002474 experimental method Methods 0.000 claims abstract description 9
- 230000007246 mechanism Effects 0.000 claims abstract description 6
- 239000012530 fluid Substances 0.000 claims description 77
- 238000012546 transfer Methods 0.000 claims description 21
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000007789 gas Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000005507 spraying Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
The invention discloses a simulation method and a system of a heat exchanger type passive containment cooling system, wherein the simulation method comprises the following steps: s100, determining an in-tube heat exchange coefficient and an out-tube heat exchange coefficient of a heat exchanger type passive containment cooling system and a heat conduction coefficient of a pipeline based on a mechanism experiment; s200, determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the pipe, the determined heat exchange coefficient outside the pipe and the determined heat conduction coefficient and a heat balance equation; s300, correcting the heat exchange coefficient outside the pipe, the heat exchange coefficient inside the pipe and the heat conduction coefficient of the pipeline based on the determined total heat exchange coefficient; s400, embedding the corrected heat exchange coefficient outside the pipe, the corrected heat exchange coefficient inside the pipe and the corrected heat conduction coefficient into a thermal hydraulic general calculation program of the pressurized water reactor nuclear power station, and completing a simulation process of the heat exchanger type passive containment cooling system. The invention obviously improves the calculation speed on the premise of ensuring the accuracy.
Description
Technical Field
The invention relates to the field of nuclear power station safety, in particular to a simulation method and a simulation system of a heat exchanger type passive containment cooling system.
Background
The third generation nuclear power technology, for example, hualong one adopts a tube bundle type heat exchanger as an passive containment cooling system (PCS for short), as shown in fig. 1. The PCS system is used for long-term heat rejection of the containment under design expansion working conditions, including accidents related to power failure of the whole plant and fault of the spraying system. When design expansion conditions (including severe accidents) occur in the power station, the containment pressure and temperature are reduced to acceptable levels, and containment integrity is maintained.
When the design expansion working condition occurs in the power station, the pressure and the temperature in the containment vessel rise rapidly. When the containment pressure is high and the safety spraying is unavailable, a containment electric isolation valve on a system downcomer receives an opening signal from a main control room or an emergency command center, and the PCS system is put into operation. The high temperature steam-air or steam-air-hydrogen (or other non-condensable gas) mixture flushes the PCS system heat exchanger surfaces. The low-temperature water from the heat exchange water tank outside the containment vessel heats and expands in the heat exchanger, and the heat in the containment vessel is guided out to the heat exchange water tank outside the containment vessel along the rising pipe of the PCS system. The temperature difference between the high-temperature mixed gas in the containment and the heat exchange water tank and the height difference between the heat exchange water tank and the heat exchanger are driving forces for driving the PCS system to perform natural circulation and leading out heat in the containment. Along with the continuous rising of the temperature of the water tank, the temperature of the heat exchange water tank reaches the saturation temperature under the corresponding pressure, and the discharged part of water vapor finally enters the external environment.
The rapid and accurate simulation of the PCS system in the accident analysis of the prototype power station is a precondition for accurately judging the state in the containment; in a prototype power station modeling experiment, accurate simulation of the PCS system can ensure the accuracy of mass energy release modeling and the authenticity of an experimental result.
The heat exchange process of the PCS heat exchange tube is shown in fig. 5. The GOTHIC program does not have a heat exchanger type PCS system model, and a conservative method for simulating the Hualong No. one PCS system by adopting the GOTHIC program generally has a certain problem.
(1) Setting PCS power, namely selecting Sp Heat Flux by the Heat exchange type of one side of the PCS Heat exchange tube Heat member 5 connected with the containment control body 6, as shown in the method of FIG. 2; and setting the PCS fixed wall temperature, namely the method for selecting SpTemp by the heat exchange type of one side of the PCS heat exchange tube heat member 5 connected with the containment control body 6 can not accurately simulate the change process of PCS power.
(2) The method shown in fig. 3 is that a water tank control body 1, a descending section control body 2, an ascending section control body 3 and a heat exchanger control body 4 are added, one side heat exchange type of the PCS heat exchange pipe thermal component 5 connected with the containment control body 6 selects a containment side PCS heat exchange coefficient relational expression, and one side heat exchange type of the PCS heat exchange pipe thermal component 5 connected with the heat exchanger control body 4 selects a film or other pipeline side PCS heat exchange coefficient model. Compared with the method shown in the step (1), the calculation result of the method is closer to the actual working condition, but the calculation speed is obviously reduced due to the fact that more control bodies are added.
Thus, to achieve both computational accuracy and efficiency, the PCS system may be modeled as in FIG. 4. Namely, considering the water tank control body 1, selecting the PCS heat exchange coefficient relation formula of the containment side through experimental verification by the heat exchange type of the side, which is connected with the containment control body 6, of the PCS heat exchange tube heat component 5, selecting the film or the PCS heat exchange coefficient model of the pipeline side through mechanical experiment verification by the heat exchange type of the side, which is connected with the water tank control body 1, and correcting the heat exchange coefficients of the pipeline side and the containment side, which are obtained through the mechanical experiment, in order to accurately simulate the simulation of the method.
Disclosure of Invention
Aiming at the defects existing in the prior art, the invention aims to provide a simulation method and a simulation system of a heat exchanger type passive containment cooling system, which can obviously improve the calculation speed on the premise of ensuring the accuracy.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows:
a method of simulating a heat exchanger type passive containment cooling system, comprising:
s100, determining an in-pipe heat exchange coefficient and an out-pipe heat exchange coefficient of a heat exchanger type passive containment cooling system and a heat conduction coefficient of a pipeline based on a mechanism experiment, wherein the out-pipe heat exchange coefficient is a heat exchange coefficient between an outside pipe fluid and an outside pipe wall, and the in-pipe heat exchange coefficient is a heat exchange coefficient between an inside pipe fluid and an inside pipe wall;
s200, determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the pipe, the determined heat exchange coefficient outside the pipe and the determined heat conduction coefficient and a heat balance equation;
s300, correcting the heat exchange coefficient outside the pipe, the heat exchange coefficient inside the pipe and the heat conduction coefficient of the pipeline based on the determined total heat exchange coefficient;
s400, embedding the corrected heat exchange coefficient outside the pipe, the corrected heat exchange coefficient inside the pipe and the corrected heat conduction coefficient into a thermal hydraulic general calculation program of the pressurized water reactor nuclear power station, and completing a simulation process of the heat exchanger type passive containment cooling system.
Further, as described above, S200 includes:
tube outside fluid heat exchange equation:
q out ＝h cont ·(T cont,bulk -T tube,out )·S out (1)
tube internal heat conduction equation:
tube inside fluid heat exchange equation:
q in ＝h in ·(T tube,in -T in,bulk )·S in (3)
S out ＝π·d 2 ·l (4)
S in ＝π·d 1 ·l (5)
according to the heat balance equation, there are:
q in ＝q tube ＝q out (6)
at the outside fluid temperature T cont,bulk And tube inside fluid temperature T in,bulk Based on the total thermal resistance R for the outside area of the tube total The method comprises the following steps:
total heat exchange coefficient h total The method comprises the following steps:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the coefficient of thermal conductivity of the pipe, W/(mK).
Further, as described above, S300 includes:
neglecting fluid in a pipeline in the simulation process of the heat exchanger type passive containment cooling system, directly exchanging heat between the atmosphere in the containment and a heat exchange water tank through a heat component, and building in order to ensure that the heat rejection capability of a simplified scheme is in accordance with the actual requirementVertical total heat exchange coefficient h total And the corrected total heat exchange coefficient h' total The relation between:
h total ·(T cont,bulk -T in,bulk )·S out ＝h' total (T cont,bulk -T tank )·S out (9)
further, there are:
according to formula (10), the modified heat exchange coefficient h 'in the tube' in The method comprises the following steps:
the corrected heat exchange coefficient h 'outside the tube' cont The method comprises the following steps:
the corrected thermal conductivity lambda' is:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the coefficient of thermal conductivity of the pipe, W/(mK).
A simulation system of a heat exchanger type passive containment cooling system, comprising:
the first determining module is used for determining an in-pipe heat exchange coefficient and an out-pipe heat exchange coefficient of the heat exchanger type passive containment cooling system and a heat conduction coefficient of a pipeline based on a mechanism experiment, wherein the out-pipe heat exchange coefficient is a heat exchange coefficient between an outside pipe fluid and an outside pipe wall, and the in-pipe heat exchange coefficient is a heat exchange coefficient between an inside pipe fluid and an inside pipe wall;
the second determining module is used for determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the pipe, the determined heat exchange coefficient outside the pipe and the determined heat conduction coefficient and the determined heat balance equation;
the correction module is used for correcting the heat exchange coefficient outside the pipe, the heat exchange coefficient inside the pipe and the heat conduction coefficient of the pipeline based on the determined total heat exchange coefficient;
and the simulation module is used for embedding the corrected heat exchange coefficient outside the pipe, the corrected heat exchange coefficient inside the pipe and the corrected heat conduction coefficient into a thermal hydraulic general calculation program of the pressurized water reactor nuclear power station to complete the simulation process of the heat exchanger type passive containment cooling system.
Further, in the system as described above, the second determining module is specifically configured to:
tube outside fluid heat exchange equation:
q out ＝h cont ·(T cont,bulk -T tube,out )·S out (1)
tube internal heat conduction equation:
tube inside fluid heat exchange equation:
q in ＝h in ·(T tube,in -T in,bulk )·S in (3)
S out ＝π·d 2 ·l (4)
S in ＝π·d 1 ·l (5)
according to the heat balance equation, there are:
q in ＝q tube ＝q out (6)
at the outside fluid temperature T cont,bulk And tube inside fluid temperature T in,bulk Based on the total thermal resistance R for the outside area of the tube total The method comprises the following steps:
total heat exchange coefficient h total The method comprises the following steps:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the coefficient of thermal conductivity of the pipe, W/(mK).
Further, in the system as described above, the correction module is specifically configured to:
in heat exchanger type passive containment coolingNeglecting fluid in a pipeline in the simulation process of the system, directly exchanging heat between the atmosphere in the containment and a heat exchange water tank through a heat component, and establishing a total heat exchange coefficient h to ensure that the heat rejection capability of a simplified scheme is practically consistent total And the corrected total heat exchange coefficient h' total The relation between:
h total ·(T cont,bulk -T in,bulk )·S out ＝h' total (T cont,bulk -T tank )·S out (9)
further, there are:
according to formula (10), the modified heat exchange coefficient h 'in the tube' in The method comprises the following steps:
the corrected heat exchange coefficient h 'outside the tube' cont The method comprises the following steps:
the corrected thermal conductivity lambda' is:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the coefficient of thermal conductivity of the pipe, W/(mK).
The invention has the beneficial effects that: the invention is based on the current hot hydraulic general calculation program GOHTIC of the pressurized water reactor nuclear power station, adopts a simplified calculation method to calculate the temperature and pressure in the containment after serious accidents, and provides a basis for reducing the pressure and the temperature of the containment to acceptable levels under the action of a heat exchanger type passive containment cooling system after the serious accidents happen to the nuclear power station and keeping the integrity of the containment. On the premise of ensuring accuracy, the calculation speed is obviously improved.
Drawings
FIG. 1 is a schematic diagram of a heat exchanger type passive containment cooling system (PCS system) provided in an embodiment of the present invention;
FIG. 2 is a schematic diagram of a method for conservative simulation of a PCS system according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a loop-type simulation method of a PCS system according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a method for rapidly simulating a PCS system according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of a heat transfer process of a heat exchange tube of a PCS system according to an embodiment of the present invention;
FIG. 6 is a schematic flow chart of a simulation method of a heat exchanger type passive containment cooling system according to an embodiment of the present disclosure;
fig. 7 is a schematic structural diagram of a simulation system of a heat exchanger type passive containment cooling system according to an embodiment of the present disclosure.
Detailed Description
In order to make the technical problems solved, the technical scheme adopted and the technical effects achieved by the invention more clear, the technical scheme of the embodiment of the invention will be further described in detail with reference to the accompanying drawings.
The embodiment of the invention provides a simulation method of a heat exchanger type passive containment cooling system, as shown in fig. 6, comprising the following steps:
s100, determining an in-pipe heat exchange coefficient and an out-pipe heat exchange coefficient of a heat exchanger type passive containment cooling system and a heat conduction coefficient of a pipeline based on a mechanism experiment, wherein the out-pipe heat exchange coefficient is a heat exchange coefficient between an outside pipe fluid and an outside pipe wall, and the in-pipe heat exchange coefficient is a heat exchange coefficient between an inside pipe fluid and an inside pipe wall;
s200, determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the pipe, the determined heat exchange coefficient outside the pipe and the determined heat conduction coefficient and a heat balance equation;
s200 includes:
tube outside fluid heat exchange equation:
q out ＝h cont ·(T cont,bulk -T tube,out )·S out (1)
tube internal heat conduction equation:
tube inside fluid heat exchange equation:
q in ＝h in ·(T tube,in -T in,bulk )·S in (3)
S out ＝π·d 2 ·l (4)
S in ＝π·d 1 ·l (5)
according to the heat balance equation, there are:
q in ＝q tube ＝q out (6)
at the outside fluid temperature T cont,bulk And tube inside fluid temperature T in,bulk Based on the total thermal resistance R for the outside area of the tube total The method comprises the following steps:
total heat exchange coefficient h total The method comprises the following steps:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the coefficient of thermal conductivity of the pipe, W/(mK).
S300, correcting the heat exchange coefficient outside the pipe, the heat exchange coefficient inside the pipe and the heat conduction coefficient of the pipeline based on the determined total heat exchange coefficient;
s300 includes:
neglecting fluid in a pipeline in the simulation process of the heat exchanger type passive containment cooling system, directly exchanging heat between the atmosphere in the containment and a heat exchange water tank through a heat component, and establishing a total heat exchange coefficient h in order to ensure that the heat rejection capacity of a simplified scheme is in accordance with the actual heat rejection capacity total And the corrected total heat exchange coefficient h' total The relation between:
h total ·(T cont,bulk -T in,bulk )·S out ＝h' total (T cont,bulk -T tank )·S out (9)
further, there are:
according to formula (10), the modified heat exchange coefficient h 'in the tube' in The method comprises the following steps:
the corrected heat exchange coefficient h 'outside the tube' cont The method comprises the following steps:
the corrected thermal conductivity lambda' is:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the coefficient of thermal conductivity of the pipe, W/(mK).
S400, embedding the corrected heat exchange coefficient outside the pipe, the corrected heat exchange coefficient inside the pipe and the corrected heat conduction coefficient into a thermal hydraulic general calculation program of the pressurized water reactor nuclear power station, and completing a simulation process of the heat exchanger type passive containment cooling system.
By adopting the method provided by the embodiment of the invention, based on the current thermal hydraulic general calculation program GOHTIC of the pressurized water reactor nuclear power station, the temperature and pressure in the containment after serious accidents are calculated by adopting a simplified calculation method, so that the basis is provided for reducing the pressure and the temperature of the containment to acceptable levels under the PCS effect after serious accidents of the nuclear power station and keeping the integrity of the containment.
It should be noted that, for simplicity of description, the foregoing method embodiments are all described as a series of acts, but it should be understood by those skilled in the art that the present invention is not limited by the order of acts described, as some steps may be performed in other orders or concurrently in accordance with the present invention. Further, those skilled in the art will also appreciate that the embodiments described in the specification are all preferred embodiments, and that the acts and modules referred to are not necessarily required for the present invention.
According to another aspect of an embodiment of the present invention, there is also provided a simulation system of a heat exchanger type passive containment cooling system, as shown in fig. 7, including:
the first determining module 100 is configured to determine an in-tube heat exchange coefficient and an out-tube heat exchange coefficient of the heat exchanger type passive containment cooling system, and a heat conduction coefficient of a pipeline, where the out-tube heat exchange coefficient is a heat exchange coefficient between an outside-tube fluid and an outside-tube wall, and the in-tube heat exchange coefficient is a heat exchange coefficient between an inside-tube fluid and an inside-tube wall;
a second determining module 200, configured to determine a total thermal resistance and a total heat exchange coefficient based on the determined heat exchange coefficient in the pipe, the determined heat exchange coefficient outside the pipe, and the determined heat balance equation;
the correction module 300 is used for correcting the heat exchange coefficient outside the pipe, the heat exchange coefficient inside the pipe and the heat conduction coefficient of the pipeline based on the determined total heat exchange coefficient;
the simulation module 400 is used for embedding the corrected heat exchange coefficient outside the pipe, the heat exchange coefficient inside the pipe and the heat conduction coefficient into a thermal hydraulic general calculation program of the pressurized water reactor nuclear power station to complete the simulation process of the heat exchanger type passive containment cooling system.
It should be noted that, the simulation system of the heat exchanger type passive containment cooling system and the method for simulating the heat exchanger type passive containment cooling system disclosed by the invention belong to the same inventive concept, and specific embodiments are not repeated.
The system provided by the embodiment of the invention is based on the current hot hydraulic general calculation program GOHTIC of the pressurized water reactor nuclear power station, and a simplified calculation method is adopted to calculate the temperature and pressure in the containment after serious accidents, so that the basis is provided for reducing the pressure and the temperature of the containment to acceptable levels under the PCS effect after serious accidents of the nuclear power station and keeping the integrity of the containment.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (2)
1. A method of simulating a heat exchanger type passive containment cooling system, comprising:
s100, determining an in-pipe heat exchange coefficient and an out-pipe heat exchange coefficient of a heat exchanger type passive containment cooling system and a heat conduction coefficient of a pipeline based on a mechanism experiment, wherein the out-pipe heat exchange coefficient is a heat exchange coefficient between an outside pipe fluid and an outside pipe wall, and the in-pipe heat exchange coefficient is a heat exchange coefficient between an inside pipe fluid and an inside pipe wall;
s200, determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the pipe, the determined heat exchange coefficient outside the pipe and the determined heat conduction coefficient and a heat balance equation;
s300, correcting the heat exchange coefficient outside the pipe, the heat exchange coefficient inside the pipe and the heat conduction coefficient of the pipeline based on the determined total heat exchange coefficient;
s400, embedding the corrected heat exchange coefficient outside the pipe, the corrected heat exchange coefficient inside the pipe and the corrected heat conduction coefficient into a thermal hydraulic general calculation program of the pressurized water reactor nuclear power station to complete a simulation process of the heat exchanger type passive containment cooling system;
s200 includes:
tube outside fluid heat exchange equation:
q out ＝h cont ·(T cont,bulk -T tube,out )·S out (1)
tube internal heat conduction equation:
tube inside fluid heat exchange equation:
q in ＝h in ·(T tube,in -T in,bulk )·S in (3)
S out ＝π·d 2 ·l (4)
S in ＝π·d 1 ·l (5)
according to the heat balance equation, there are:
q in ＝q tube ＝q out (6)
at the outside fluid temperature T cont,bulk And tube inside fluid temperature T in,bulk Based on the total thermal resistance R for the outside area of the tube total The method comprises the following steps:
total heat exchange coefficient h total The method comprises the following steps:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the λ is the coefficient of thermal conductivity of the pipe, W/(mK);
s300 includes:
neglecting fluid in a pipeline in the simulation process of the heat exchanger type passive containment cooling system, directly exchanging heat between the atmosphere in the containment and a heat exchange water tank through a heat component, and establishing a total heat exchange coefficient h in order to ensure that the heat rejection capacity of a simplified scheme is in accordance with the actual heat rejection capacity total And the corrected total heat exchange coefficient h' total The relation between:
h total ·(T cont,bulk -T in,bulk )·S out ＝h' total (T cont,bulk -T tank )·S out (9)
further, there are:
according to formula (10), the modified heat exchange coefficient h 'in the tube' in The method comprises the following steps:
the corrected heat exchange coefficient h 'outside the tube' cont The method comprises the following steps:
the corrected thermal conductivity lambda' is:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the heat conductivity coefficient of the pipeline, W/(m.K), T tank Is the tank temperature.
2. A simulation system of a heat exchanger type passive containment cooling system, comprising:
the first determining module is used for determining an in-pipe heat exchange coefficient and an out-pipe heat exchange coefficient of the heat exchanger type passive containment cooling system and a heat conduction coefficient of a pipeline based on a mechanism experiment, wherein the out-pipe heat exchange coefficient is a heat exchange coefficient between an outside pipe fluid and an outside pipe wall, and the in-pipe heat exchange coefficient is a heat exchange coefficient between an inside pipe fluid and an inside pipe wall;
the second determining module is used for determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the pipe, the determined heat exchange coefficient outside the pipe and the determined heat conduction coefficient and the determined heat balance equation;
the correction module is used for correcting the heat exchange coefficient outside the pipe, the heat exchange coefficient inside the pipe and the heat conduction coefficient of the pipeline based on the determined total heat exchange coefficient;
the simulation module is used for embedding the corrected heat exchange coefficient outside the pipe, the corrected heat exchange coefficient inside the pipe and the corrected heat conduction coefficient into a thermal hydraulic general calculation program of the pressurized water reactor nuclear power station to complete the simulation process of the heat exchanger type passive containment cooling system;
the second determining module is specifically configured to:
tube outside fluid heat exchange equation:
q out ＝h cont ·(T cont,bulk -T tube,out )·S out (1)
tube internal heat conduction equation:
tube inside fluid heat exchange equation:
q in ＝h in ·(T tube,in -T in,bulk )·S in (3)
S out ＝π·d 2 ·l (4)
S in ＝π·d 1 ·l (5)
according to the heat balance equation, there are:
q in ＝q tube ＝q out (6)
at the outside fluid temperature T cont,bulk And tube inside fluid temperature T in,bulk Based on the total thermal resistance R for the outside area of the tube total The method comprises the following steps:
total heat exchange coefficient h total The method comprises the following steps:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the λ is the coefficient of thermal conductivity of the pipe, W/(mK);
the correction module is specifically configured to:
neglecting fluid in a pipeline in the simulation process of the heat exchanger type passive containment cooling system, directly exchanging heat between the atmosphere in the containment and a heat exchange water tank through a heat component, and establishing a total heat exchange coefficient h in order to ensure that the heat rejection capacity of a simplified scheme is in accordance with the actual heat rejection capacity total And the corrected total heat exchange coefficient h' total The relation between:
h total ·(T cont,bulk -T in,bulk )·S out ＝h' total (T cont,bulk -T tank )·S out (9)
further, there are:
according to formula (10), the modified heat exchange coefficient h 'in the tube' in The method comprises the following steps:
the corrected heat exchange coefficient h 'outside the tube' cont The method comprises the following steps:
the corrected thermal conductivity lambda' is:
wherein q out The heat exchange power of the fluid outside the tube to the tube is W; q tube The heat conduction power W of the outer wall of the tube faces the inner wall surface; q in The heat exchange power from the inner wall of the tube to the fluid in the tube is W; h is a cont W/(m) is the heat exchange coefficient outside the tube 2 ·K)；h in W/(m) is the heat exchange coefficient in the tube 2 ·K)；T cont,bulk The temperature of the fluid outside the tube, K; t (T) in,bulk Is the temperature of the fluid inside the tube, K; t (T) tube,out The temperature of the outer wall of the tube, K; t (T) tube,in The temperature of the inner wall of the tube is K; l is the length of the pipeline, m; d, d 1 Tube inside diameter, m; d, d 2 Is the outer diameter of the tube, m; s is S out For heat transfer area outside the tube, m 2 ；S in For heat transfer area in the tube, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Lambda is the heat conductivity coefficient of the pipeline, W/(m.K), T tank Is the tank temperature.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110785696.3A CN113704959B (en) | 2021-07-12 | 2021-07-12 | Simulation method and system for heat exchanger type passive containment cooling system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110785696.3A CN113704959B (en) | 2021-07-12 | 2021-07-12 | Simulation method and system for heat exchanger type passive containment cooling system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113704959A CN113704959A (en) | 2021-11-26 |
| CN113704959B true CN113704959B (en) | 2024-04-09 |
Family
ID=78648485
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110785696.3A Active CN113704959B (en) | 2021-07-12 | 2021-07-12 | Simulation method and system for heat exchanger type passive containment cooling system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113704959B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102915776A (en) * | 2012-09-27 | 2013-02-06 | 中国核电工程有限公司 | Testing method for passive safety shell heat lead-out system |
| WO2016011569A1 (en) * | 2014-07-24 | 2016-01-28 | 哈尔滨工程大学 | Containment cooling system, and containment and reactor pressure vessel joint cooling system |
| CN105788682A (en) * | 2016-03-22 | 2016-07-20 | 上海交通大学 | Nuclear reactor containment cooling simulation system and method based on external spraying |
| CN107958712A (en) * | 2017-11-23 | 2018-04-24 | 中国核电工程有限公司 | A kind of heat pipe exchanging type passive containment thermal conduction system |
| CN108010592A (en) * | 2016-10-28 | 2018-05-08 | 国核华清（北京）核电技术研发中心有限公司 | For nuclear power plant containment shell performance test and the contracting proportion composite experimental rig of research |
| CN109902433A (en) * | 2019-03-15 | 2019-06-18 | 西安交通大学 | Across the dimension coupling process of presurized water reactor passive containment residual heat removal system |
| CN112164482A (en) * | 2020-09-25 | 2021-01-01 | 中国核电工程有限公司 | Thermal hydraulic comprehensive test device for containment |
| CN112307626A (en) * | 2020-11-02 | 2021-02-02 | 西安交通大学 | Operation and influence analysis method for multiple groups of passive heat export systems in containment |
-
2021
- 2021-07-12 CN CN202110785696.3A patent/CN113704959B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102915776A (en) * | 2012-09-27 | 2013-02-06 | 中国核电工程有限公司 | Testing method for passive safety shell heat lead-out system |
| WO2016011569A1 (en) * | 2014-07-24 | 2016-01-28 | 哈尔滨工程大学 | Containment cooling system, and containment and reactor pressure vessel joint cooling system |
| CN105788682A (en) * | 2016-03-22 | 2016-07-20 | 上海交通大学 | Nuclear reactor containment cooling simulation system and method based on external spraying |
| CN108010592A (en) * | 2016-10-28 | 2018-05-08 | 国核华清（北京）核电技术研发中心有限公司 | For nuclear power plant containment shell performance test and the contracting proportion composite experimental rig of research |
| CN107958712A (en) * | 2017-11-23 | 2018-04-24 | 中国核电工程有限公司 | A kind of heat pipe exchanging type passive containment thermal conduction system |
| CN109902433A (en) * | 2019-03-15 | 2019-06-18 | 西安交通大学 | Across the dimension coupling process of presurized water reactor passive containment residual heat removal system |
| CN112164482A (en) * | 2020-09-25 | 2021-01-01 | 中国核电工程有限公司 | Thermal hydraulic comprehensive test device for containment |
| CN112307626A (en) * | 2020-11-02 | 2021-02-02 | 西安交通大学 | Operation and influence analysis method for multiple groups of passive heat export systems in containment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113704959A (en) | 2021-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107784156B (en) | Method for calculating parameters of steam discharge system of nuclear power plant | |
| CN109902433B (en) | Cross-dimension coupling method for pressurized water reactor passive containment waste heat discharge system | |
| CN106952669B (en) | Stagnation pressure external container cooling test rack in a kind of fusant heap | |
| CN111274748B (en) | Cross-dimension coupling calculation method for pool type sodium-cooled fast reactor passive waste heat removal system | |
| CN110288135B (en) | Drainage water level energy-saving optimization method for high-pressure heating system | |
| CN112307626B (en) | Operation and influence analysis method for multiple groups of passive heat export systems in containment | |
| CN110472270B (en) | Modeling optimization method of sodium cooled reactor passive waste heat removal system | |
| CN111723450A (en) | Safety analysis method and system for nuclear power plant | |
| CN113704959B (en) | Simulation method and system for heat exchanger type passive containment cooling system | |
| CN113689966B (en) | Small-size villiaumite cooling high temperature reactor passive exhaust system comprehensive experiment device | |
| CN108565034B (en) | U-tube flows backwards the system and method for thermal environment simulation under passive service condition | |
| Lv et al. | Experimental study of DRACS thermal performance in a low-temperature test facility | |
| CN113948229B (en) | Analysis method for stable operation condition of ocean nuclear power heat pipe type waste heat discharging system | |
| Lv et al. | Experimental study of DRACS steady-state and transient performance | |
| CN206247373U (en) | A kind of pipe heat exchanger residual neat recovering system | |
| CN113283189B (en) | Cross-dimension coupling analysis method for containment and passive containment air cooling system | |
| CN111735648B (en) | System and method for testing thermal performance of pressurized water reactor main evaporator proportional prototype | |
| Wang et al. | Numerical analysis of transient pressure variation in the condenser of a nuclear power station | |
| Polidori–ENEA et al. | Validazione dei Modelli RELAP5 dello Scambiatore di Calore HERO-2 | |
| CN218585188U (en) | Power station heat exchange simulation system | |
| CN115036049A (en) | Passive residual heat removal test device based on hot water boiler simulation boundary | |
| Peterson et al. | Functional Requirements Overview For a 50-MW (t) Liquid-Salt Intermediate Loop for NGNP | |
| Li et al. | Research on Characteristics of Condensation Heat Transfer for Marine Nuclear Power Platform PRHR HX Under Tilting Condition | |
| Kim et al. | A 2.5 MWT Passive Decay Heat Removal System Design for PGSFR | |
| CN115910400A (en) | Test device and method for simulating heat pipe failure transient accident in heat pipe micro-stack |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |