CN113704959A - Simulation method and system of heat exchanger type passive containment cooling system - Google Patents

Simulation method and system of heat exchanger type passive containment cooling system Download PDF

Info

Publication number
CN113704959A
CN113704959A CN202110785696.3A CN202110785696A CN113704959A CN 113704959 A CN113704959 A CN 113704959A CN 202110785696 A CN202110785696 A CN 202110785696A CN 113704959 A CN113704959 A CN 113704959A
Authority
CN
China
Prior art keywords
tube
heat exchange
heat
coefficient
exchange 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.)
Granted
Application number
CN202110785696.3A
Other languages
Chinese (zh)
Other versions
CN113704959B (en
Inventor
李精精
王辉
黄政
孙婧
周喆
石雪垚
常愿
王贺南
雷宁博
李汉辰
蔡盟利
孙晓晖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Nuclear Power Engineering Co Ltd
Original Assignee
China Nuclear Power Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Nuclear Power Engineering Co Ltd filed Critical China Nuclear Power Engineering Co Ltd
Priority to CN202110785696.3A priority Critical patent/CN113704959B/en
Publication of CN113704959A publication Critical patent/CN113704959A/en
Application granted granted Critical
Publication of CN113704959B publication Critical patent/CN113704959B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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)
  • Monitoring And Testing Of Nuclear Reactors (AREA)
  • Heat-Exchange Devices With Radiators And Conduit Assemblies (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-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; s200, determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the tube, the heat exchange coefficient and the heat conductivity coefficient outside the tube 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 external heat exchange coefficient, internal heat exchange coefficient and heat conductivity coefficient into a thermodynamic and hydraulic general calculation program of the pressurized water reactor nuclear power station, and completing the 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

Simulation method and system of heat exchanger type passive containment cooling system
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
For the third generation nuclear power technology, for example, hualong number one adopts a tube bundle type heat exchanger as a passive containment cooling system (PCS for short), as shown in fig. 1. The PCS system is used for long-term heat removal of the containment vessel under the design expansion working condition, and comprises accidents related to station blackout and spraying system faults. When the power station has a design expansion working condition (including serious accidents), the pressure and the temperature of the containment vessel are reduced to acceptable levels, and the integrity of the containment vessel is maintained.
When the power station has a design expansion working condition, the pressure and the temperature in the containment vessel rapidly rise. When the pressure of the containment vessel is high and the safety spray is unavailable, the electric containment vessel isolation valve on the system descending pipe receives an opening signal from a main control room or an emergency command center, and the PCS system is put into operation. A mixture of high temperature steam-air or steam-air-hydrogen (or other non-condensable gas) flushes the PCS system heat exchanger surfaces. The low-temperature water from the containment heat exchange water tank is heated and expanded in the heat exchanger, and the heat in the containment is led out to the containment heat exchange water tank along the PCS system ascending pipe. 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 naturally circulate and lead out heat in the containment. With the continuous rise of the temperature of the water tank, the temperature of the heat exchange water tank reaches the saturation temperature under the corresponding pressure, and partial water vapor is discharged and 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 experiment result.
The heat exchange process of the PCS heat exchange pipe is shown in figure 5. The GOTHIC program does not have a heat exchanger type PCS system model, and a conservative method for simulating the Hualongyi PCS system by adopting the GOTHIC program has certain problems.
(1) As shown in fig. 2, PCS power is set, that is, Sp Heat Flux is selected as the Heat exchange type at one side of the PCS Heat exchange tube thermal member 5 connected to the containment control body 6; and the PCS fixed wall temperature is set, namely the method of selecting Sp Temp for the heat exchange type at one side of the PCS heat exchange tube hot member 5 connected with the containment control body 6 cannot accurately simulate the change process of PCS power.
(2) As shown in fig. 3, a water tank control body 1, a descent control body 2, an ascent control body 3 and a heat exchanger control body 4 are added, a safe shell side PCS heat exchange coefficient relation is selected for the heat exchange type of one side of a PCS heat exchange tube thermal member 5 connected with a containment control body 6, and a film or other pipeline side PCS heat exchange coefficient models are selected for the heat exchange type of one side of the PCS heat exchange tube thermal member 5 connected with the heat exchanger control body 4. The calculation result of the method is closer to the actual working condition than the method shown in (1), but the calculation speed is obviously reduced due to the fact that more control bodies are added.
Therefore, the PCS system can be simulated by the method shown in fig. 4 for both accuracy and efficiency of the calculation. The method comprises the steps that a safety shell side PCS heat exchange coefficient relation formula which is verified through experiments is selected according to the heat exchange type of one side of a water tank control body 1, a PCS heat exchange pipe thermal member 5 connected with a containment control body 6, and a film or a pipeline side PCS heat exchange coefficient model which is verified through mechanism experiments is selected according to the heat exchange type of one side of the PCS heat exchange pipe thermal member 5 connected with the water tank control body 1.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a method and a system for simulating 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 purpose, the technical scheme adopted by the invention is as follows:
a simulation method of a heat exchanger type passive containment cooling system comprises the following steps:
s100, 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 out-pipe fluid and an outer wall of a pipe, and the in-pipe heat exchange coefficient is a heat exchange coefficient between the in-pipe fluid and an inner wall of the pipe;
s200, determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the tube, the heat exchange coefficient and the heat conductivity coefficient outside the tube 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 external heat exchange coefficient, internal heat exchange coefficient and heat conductivity coefficient into a thermodynamic and hydraulic general calculation program of the pressurized water reactor nuclear power station, and completing the simulation process of the heat exchanger type passive containment cooling system.
Further, the method as described above, S200 includes:
equation of heat exchange of fluid outside the tube:
qout＝hcont·(Tcont，bulk-Ttube，out)·Sout (1)
the equation for the heat transfer inside the tube:
Figure BDA0003158653820000031
equation of heat exchange of fluid inside the tube:
qin＝hin·(Ttube，in-Tin，bulk)·Sin (3)
Sout＝π·d2·l (4)
Sin＝π·d1·l (5)
according to the heat balance equation, there are:
qin＝qtube＝qout (6)
with outside fluid temperature Tcont，bulkAnd the temperature T of the fluid inside the tubein，bulkFor reference, the total thermal resistance R for the area outside the tubetotalComprises the following steps:
Figure BDA0003158653820000032
total heat transfer coefficient htotalComprises the following steps:
Figure BDA0003158653820000033
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont，bulkThe temperature of the fluid outside the tube, K; t isin，bulkIs the temperature of the fluid inside the tube, K; t istube，outIs the tube outer wall temperature, K; t istube，inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
Further, the method as described above, S300 includes:
neglecting the fluid in the pipeline in the simulation process of the heat exchanger type passive containment cooling system, directly exchanging heat with the heat exchange water tank by the atmosphere in the containment through the thermal component, and establishing a total heat exchange coefficient h for ensuring that the heat removal capability of the simplified scheme is consistent with the actual capabilitytotalAnd the corrected total heat exchange coefficient h'totalThe relation between:
htotal·(Tcont，bulk-Tin，bulk)·Sout＝h′total(Tcont，bulk-Ttank)·Sout (9)
further, there are:
Figure BDA0003158653820000041
according to a formula (10), the corrected heat exchange coefficient h 'in the tube'inComprises the following steps:
Figure BDA0003158653820000042
corrected external heat exchange coefficient h'outComprises the following steps:
Figure BDA0003158653820000043
the corrected thermal conductivity λ' is:
Figure BDA0003158653820000044
Figure BDA0003158653820000045
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont，bulkThe temperature of the fluid outside the tube, K; t isin，bulkIs the temperature of the fluid inside the tube, K; t istube，outIs the tube outer wall temperature, K; t istube，inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
A simulation system for 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 fluid and an outer wall of a pipe, and the in-pipe heat exchange coefficient is a heat exchange coefficient between the inside fluid and an inner wall of the pipe;
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 tube, the heat exchange coefficient and the heat conductivity coefficient outside the tube and a 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 external heat exchange coefficient, internal heat exchange coefficient and heat conductivity coefficient into a general calculation program of the thermal and hydraulic power 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:
equation of heat exchange of fluid outside the tube:
qout＝hcont·(Tcont，bulk-Ttube，out)·Sout (1)
the equation for the heat transfer inside the tube:
Figure BDA0003158653820000051
equation of heat exchange of fluid inside the tube:
qin＝hin·(Ttube，in-Tin，bulk)·Sin (3)
Sout＝π·d2·l (4)
Sin＝π·d1·l (5)
according to the heat balance equation, there are:
qin＝qtube＝qout (6)
with outside fluid temperature Tcont，bulkAnd the temperature T of the fluid inside the tubein，bulkFor reference, the total thermal resistance R for the area outside the tubetotalComprises the following steps:
Figure BDA0003158653820000052
total heat transfer coefficient htotalComprises the following steps:
Figure BDA0003158653820000061
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont，bulkThe temperature of the fluid outside the tube, K; t isin，bulkIs the temperature of the fluid inside the tube, K; t istube，outIs the tube outer wall temperature, K; t istube，inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
Further, in the system as described above, the modification module is specifically configured to:
neglecting the fluid in the pipeline in the simulation process of the heat exchanger type passive containment cooling system, directly exchanging heat with the heat exchange water tank by the atmosphere in the containment through the thermal component, and establishing a total heat exchange coefficient h for ensuring that the heat removal capability of the simplified scheme is consistent with the actual capabilitytotalAnd the corrected total heat exchange coefficient h'totalThe relation between:
htotal·(Tcont，bulk-Tin，bulk)·Sout＝h′total(Tcont，bulk-Ttank)·Sout (9)
further, there are:
Figure BDA0003158653820000062
according to a formula (10), the corrected heat exchange coefficient h 'in the tube'inComprises the following steps:
Figure BDA0003158653820000063
modified outside tube exchangeThermal coefficient h'outComprises the following steps:
Figure BDA0003158653820000064
the corrected thermal conductivity λ' is:
Figure BDA0003158653820000071
Figure BDA0003158653820000072
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont，bulkThe temperature of the fluid outside the tube, K; t isin，bulkIs the temperature of the fluid inside the tube, K; t istube，outIs the tube outer wall temperature, K; t istube，inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
The invention has the beneficial effects that: the invention is based on the current general calculation program GOHTIC of the thermal hydraulic power of the pressurized water reactor nuclear power station, adopts a simplified calculation method to calculate the temperature and the pressure in the containment after a serious accident, 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 accident of the nuclear power station and keeping the integrity of the containment. On the premise of ensuring the accuracy, the calculation speed is obviously improved.
Drawings
FIG. 1 is a schematic structural 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 conservative simulation method for 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 rapid simulation method for 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 provided in an embodiment of the 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 invention;
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 invention.
Detailed Description
In order to make the technical problems solved, the technical solutions adopted, and the technical effects achieved by the present invention clearer, the technical solutions of the embodiments of the present 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, which comprises the following steps of:
s100, 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 out-pipe fluid and an outer wall of a pipe, and the in-pipe heat exchange coefficient is a heat exchange coefficient between the in-pipe fluid and an inner wall of the pipe;
s200, determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the tube, the heat exchange coefficient and the heat conductivity coefficient outside the tube and a heat balance equation;
s200 comprises the following steps:
equation of heat exchange of fluid outside the tube:
qout＝hcont·(Tcont，bulk-Ttube，out)·Sout (1)
the equation for the heat transfer inside the tube:
Figure BDA0003158653820000081
equation of heat exchange of fluid inside the tube:
qin＝hin·(Ttube，in-Tin，bulk)·Sin (3)
Sout＝π·d2·l (4)
Sin＝π·d1·l (5)
according to the heat balance equation, there are:
qin＝qtube＝qout (6)
with outside fluid temperature Tcont，bulkAnd the temperature T of the fluid inside the tubein，bulkFor reference, the total thermal resistance R for the area outside the tubetotalComprises the following steps:
Figure BDA0003158653820000082
total heat transfer coefficient htotalComprises the following steps:
Figure BDA0003158653820000083
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont，bulkThe temperature of the fluid outside the tube, K; t isin，bulkIs the temperature of the fluid inside the tube, K; t istube，outIs the tube outer wall temperature, K; t istube，inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
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 comprises the following steps:
neglecting the fluid in the pipeline in the simulation process of the heat exchanger type passive containment cooling system, directly exchanging heat with the heat exchange water tank by the atmosphere in the containment through the thermal component, and establishing a total heat exchange coefficient h for ensuring that the heat removal capability of the simplified scheme is consistent with the actual capabilitytotalAnd the corrected total heat exchange coefficient h'totalThe relation between:
htotal·(Tcont，bulk-Tin，bulk)·Sout＝h′total(Tcont，bulk-Ttank)·Sout (9)
further, there are:
Figure BDA0003158653820000091
according to a formula (10), the corrected heat exchange coefficient h 'in the tube'inComprises the following steps:
Figure BDA0003158653820000092
corrected external heat exchange coefficient h'outComprises the following steps:
Figure BDA0003158653820000093
the corrected thermal conductivity λ' is:
Figure BDA0003158653820000094
Figure BDA0003158653820000095
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont，bulkThe temperature of the fluid outside the tube, K; t isin，bulkIs the temperature of the fluid inside the tube, K; t istube，outIs the tube outer wall temperature, K; t istube，inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
S400, embedding the corrected external heat exchange coefficient, internal heat exchange coefficient and heat conductivity coefficient into a thermodynamic and hydraulic general calculation program of the pressurized water reactor nuclear power station, and completing the 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 general calculation program GOHTIC of the thermal and hydraulic power of the pressurized water reactor nuclear power station, the simplified calculation method is adopted to calculate the temperature and the pressure in the containment after a serious accident, so that the pressure and the temperature of the containment can be reduced to acceptable levels under the action of the PCS after the serious accident of the nuclear power station, and the basis is provided for keeping the integrity of the containment.
It should be noted that, for simplicity of description, the above-mentioned method embodiments are described as a series of acts or combination of acts, but those skilled in the art will recognize that the present invention is not limited by the order of acts, as some steps may occur in other orders or concurrently in accordance with the invention. Further, those skilled in the art should also appreciate that the embodiments described in the specification are preferred embodiments and that the acts and modules referred to are not necessarily required by the invention.
According to another aspect of the 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, based on a mechanism experiment, an in-tube heat exchange coefficient and an out-of-tube heat exchange coefficient of the heat exchanger type passive containment cooling system, and a heat conductivity coefficient of a pipeline, where the out-of-tube heat exchange coefficient is a heat exchange coefficient between a fluid on the outer side of the tube and the outer wall of the tube, and the in-of-tube heat exchange coefficient is a heat exchange coefficient between a fluid on the inner side of the tube and the inner wall of the tube;
a second determining module 200, configured to determine a total thermal resistance and a total heat transfer coefficient based on the determined in-tube heat transfer coefficient, the determined out-of-tube heat transfer coefficient, and the determined heat transfer coefficient, and a heat balance equation;
a correction module 300 for correcting the heat transfer coefficient outside the tube, the heat transfer coefficient inside the tube and the heat transfer coefficient of the pipe based on the determined total heat transfer coefficient;
and the simulation module 400 is used for embedding the corrected external heat exchange coefficient, internal heat exchange coefficient and heat conductivity coefficient into a general calculation program of the thermal and hydraulic power 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 of the present invention and the method for simulating the heat exchanger type passive containment cooling system belong to the same inventive concept, and detailed description is omitted.
The system provided by the embodiment of the invention is based on the current general calculation program GOHTIC of the thermal and hydraulic power of the pressurized water reactor nuclear power station, adopts a simplified calculation method to calculate the temperature and the pressure in the containment after a serious accident, and provides a basis for reducing the pressure and the temperature of the containment to acceptable levels under the action of the PCS after the serious accident of the nuclear power station and keeping the integrity of the containment.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is intended to include such modifications and variations.

Claims (6)

1. A simulation method of a heat exchanger type passive containment cooling system is characterized by comprising the following steps:
s100, 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 out-pipe fluid and an outer wall of a pipe, and the in-pipe heat exchange coefficient is a heat exchange coefficient between the in-pipe fluid and an inner wall of the pipe;
s200, determining total thermal resistance and total heat exchange coefficient based on the determined heat exchange coefficient in the tube, the heat exchange coefficient and the heat conductivity coefficient outside the tube 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 external heat exchange coefficient, internal heat exchange coefficient and heat conductivity coefficient into a thermodynamic and hydraulic general calculation program of the pressurized water reactor nuclear power station, and completing the simulation process of the heat exchanger type passive containment cooling system.
2. The method of claim 1, wherein S200 comprises:
equation of heat exchange of fluid outside the tube:
qout＝hcont·(Tcont,bulk-Ttube,out)·Sout (1)
the equation for the heat transfer inside the tube:
Figure FDA0003158653810000011
equation of heat exchange of fluid inside the tube:
qin＝hin·(Ttube,in-Tin,bulk)·Sin (3)
Sout＝π·d2·l (4)
Sin＝π·d1·l (5)
according to the heat balance equation, there are:
qin＝qtube＝qout (6)
with outside fluid temperature Tcont,bulkAnd the temperature T of the fluid inside the tubein,bulkFor reference, the total thermal resistance R for the area outside the tubetotalComprises the following steps:
Figure FDA0003158653810000012
total heat transfer coefficient htotalComprises the following steps:
Figure FDA0003158653810000021
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont,bulkThe temperature of the fluid outside the tube, K; t isin,bulkIs the temperature of the fluid inside the tube, K; t istube,outIs the tube outer wall temperature, K; t istube,inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
3. The method of claim 2, wherein S300 comprises:
simulation process of heat exchanger type passive containment cooling systemNeglecting the fluid in the pipeline, directly exchanging heat between the atmosphere in the containment vessel and the heat exchange water tank through the thermal component, and establishing a total heat exchange coefficient h for ensuring that the heat extraction capability of the simplified scheme is consistent with the realitytotalAnd the corrected total heat exchange coefficient h'totalThe relation between:
htotal·(Tcont,bulk-Tin,bulk)·Sout＝h′total(Tcont,bulk-Ttank)·Sout (9)
further, there are:
Figure FDA0003158653810000022
according to a formula (10), the corrected heat exchange coefficient h 'in the tube'inComprises the following steps:
Figure FDA0003158653810000023
corrected external heat exchange coefficient h'outComprises the following steps:
Figure FDA0003158653810000024
the corrected thermal conductivity λ' is:
Figure FDA0003158653810000031
Figure FDA0003158653810000032
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinFor changing the inner wall of the tube to the fluid in the tubeThermal power, W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont,bulkThe temperature of the fluid outside the tube, K; t isin,bulkIs the temperature of the fluid inside the tube, K; t istube,outIs the tube outer wall temperature, K; t istube,inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
4. A simulation system of a heat exchanger type passive containment cooling system is characterized by 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 fluid and an outer wall of a pipe, and the in-pipe heat exchange coefficient is a heat exchange coefficient between the inside fluid and an inner wall of the pipe;
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 tube, the heat exchange coefficient and the heat conductivity coefficient outside the tube and a 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 external heat exchange coefficient, internal heat exchange coefficient and heat conductivity coefficient into a general calculation program of the thermal and hydraulic power of the pressurized water reactor nuclear power station to complete the simulation process of the heat exchanger type passive containment cooling system.
5. The system of claim 4, wherein the second determination module is specifically configured to:
equation of heat exchange of fluid outside the tube:
qout＝hcont·(Tcont,bulk-Ttube,out)·Sout (1)
the equation for the heat transfer inside the tube:
Figure FDA0003158653810000033
equation of heat exchange of fluid inside the tube:
qin＝hin·(Ttube,in-Tin,bulk)·Sin (3)
Sout＝π·d2·l (4)
Sin＝π·d1·l (5)
according to the heat balance equation, there are:
qin＝qtube＝qout (6)
with outside fluid temperature Tcont,bulkAnd the temperature T of the fluid inside the tubein,bulkFor reference, the total thermal resistance R for the area outside the tubetotalComprises the following steps:
Figure FDA0003158653810000041
total heat transfer coefficient htotalComprises the following steps:
Figure FDA0003158653810000042
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont,bulkThe temperature of the fluid outside the tube, K; t isin,bulkIs the temperature of the fluid inside the tube, K; t istube,outIs the tube outer wall temperature, K; t istube,inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
6. The system of claim 5, wherein the modification module is specifically configured to:
neglecting the fluid in the pipeline in the simulation process of the heat exchanger type passive containment cooling system, directly exchanging heat with the heat exchange water tank by the atmosphere in the containment through the thermal component, and establishing a total heat exchange coefficient h for ensuring that the heat removal capability of the simplified scheme is consistent with the actual capabilitytotalAnd the corrected total heat exchange coefficient h'totalThe relation between:
htotal·(Tcont,bulk-Tin,bulk)·Sout＝h′total(Tcont,bulk-Ttank)·Sout (9)
further, there are:
Figure FDA0003158653810000051
according to a formula (10), the corrected heat exchange coefficient h 'in the tube'inComprises the following steps:
Figure FDA0003158653810000052
corrected external heat exchange coefficient h'outComprises the following steps:
Figure FDA0003158653810000053
the corrected thermal conductivity λ' is:
Figure FDA0003158653810000054
Figure FDA0003158653810000055
wherein q isoutThe heat exchange power from the fluid outside the tube to the inside of the tube, W; q. q.stubeThe heat conduction power, W, of the outer wall of the tube facing the inner wall surface; q. q.sinThe heat exchange power from the inner wall of the tube to the fluid in the tube is W; h iscontIs the heat transfer coefficient outside the tube, W/(m)2·K)；hinIs the heat exchange coefficient in the tube, W/(m)2·K)；Tcont,bulkThe temperature of the fluid outside the tube, K; t isin,bulkIs the temperature of the fluid inside the tube, K; t istube,outIs the tube outer wall temperature, K; t istube,inIs the tube inner wall temperature, K; l is the length of the pipeline, m; d1Pipe inner diameter, m; d2Is the outside diameter of the tube, m; soutIs the heat transfer area outside the tube, m2；SinIs the heat transfer area in the tube, m2(ii) a λ is the thermal conductivity of the pipe, W/(m.K).
CN202110785696.3A 2021-07-12 2021-07-12 Simulation method and system for heat exchanger type passive containment cooling system Active CN113704959B (en)

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 true CN113704959A (en) 2021-11-26
CN113704959B 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)

* Cited by examiner, † Cited by third party
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

Patent Citations (8)

* Cited by examiner, † Cited by third party
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
CN113704959B (en) 2024-04-09

Similar Documents

Publication Publication Date Title
CN109902433B (en) Cross-dimension coupling method for pressurized water reactor passive containment waste heat discharge system
CN108010592B (en) Shrinkage ratio comprehensive test device for nuclear power station containment performance test and research
CN106952669B (en) Stagnation pressure external container cooling test rack in a kind of fusant heap
CN112164482B (en) Comprehensive thermal hydraulic test device for containment
CN111274748B (en) Cross-dimension coupling calculation method for pool type sodium-cooled fast reactor passive waste heat removal system
Reyes et al. Testing of the multi-application small light water reactor (MASLWR) passive safety systems
CN103545000B (en) Secondary side passive residual heat removal heat exchanger analogue means and analogy method
CN111723450A (en) Safety analysis method and system for nuclear power plant
CN110472270B (en) Modeling optimization method of sodium cooled reactor passive waste heat removal system
KR101651466B1 (en) Verification test device for passive residual heat removal system of a research reactor
CN112307626B (en) Operation and influence analysis method for multiple groups of passive heat export systems in containment
CN111680417B (en) Thermal design method suitable for sodium water direct-current steam generator
CN113704959B (en) Simulation method and system for heat exchanger type passive containment cooling system
Lv et al. Experimental study of DRACS thermal performance in a low-temperature test facility
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
CN115036049A (en) Passive residual heat removal test device based on hot water boiler simulation boundary
Polidori–ENEA et al. Validazione dei Modelli RELAP5 dello Scambiatore di Calore HERO-2
CN110489796B (en) Effectiveness judgment method of heat pipe-based nuclear reactor heat exporting system
Dai et al. Research on Thermal Hydraulic Characteristics of Passive Containment Cooling System for Ship Nuclear Power Platform
CN115910400A (en) Test device and method for simulating heat pipe failure transient accident in heat pipe micro-stack
Kang et al. CFD Analysis for a Westinghouse Natural Circulation Experiment during Severe Accidents
Qin et al. The application of heat pipe discharge containment heat experimental study: new type of passive containment heat removal system concept design
CN117423486A (en) High temperature heat pipe tube bank inefficacy experimental apparatus

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