Method and device for determining total atmospheric pressure in containment

Abstract

The invention discloses a method and a device for determining total atmospheric pressure in a containment, and relates to the technical field of nuclear industry. The method comprises the following steps: S110, acquiring a hydrogen concentration range in a containment; and S120, determining a first atmospheric total pressure Pa2 in the containment according to the hydrogen concentration range. According to the embodiment of the invention, the first atmospheric total pressure in the containment is determined according to the hydrogen concentration range, and compared with an integrated serious accident analysis program, a typical serious accident sequence of the nuclear power plant is calculated to obtain a hydrogen release amount curve and a containment atmospheric pressure temperature curve; and then the hydrogen release amount curve and the containment atmospheric pressure and temperature curve serve as input, a conventional method of calculating the atmospheric total pressure of the containment atmosphere after the hydrogen is completely combusted in an adiabatic and isometric manner is adopted, the method breaks through conventional thinking, and the first atmospheric total pressure in the containment can be accurately determined.

Description

Technical Field

[0001] The present application belongs to the field of nuclear industry technology, and specifically relates to a method and device for determining the total atmospheric pressure in a containment vessel. Background Art

[0002] In the event of a serious accident at a nuclear power plant, the core fuel cladding is oxidized by water or water vapor to generate a large amount of hydrogen. This hydrogen is released into the containment through the reactor primary loop rupture or valve. There is a risk of combustion when hydrogen mixes with the containment atmosphere. Hydrogen combustion can be roughly divided into five types, namely ignition, slow diffusion combustion, flame acceleration, combustion-explosion transition and explosion. Among them, flame acceleration, combustion-explosion transition and explosion are not easy to occur in the containment, and hydrogen slow diffusion combustion is the most likely type of containment hydrogen combustion. The Nuclear Safety Guide HAD 102 / 06-2020 stipulates that "when the conditions of flame acceleration and high dynamic pressure load are not met, the adiabatic isochoric complete combustion pressure curve under the hydrogen slow combustion state should be used to calculate and determine the overall and local pressure boundary loads". At the same time, the Nuclear Safety Guide HAD 102 also stipulates the quality of hydrogen combustion, requiring the assessment of the risks caused by "100% hydrogen produced by zirconium-water reaction" in the core active area cladding.

[0003] How to evaluate the total atmospheric pressure of the containment atmosphere after adiabatic isochoric complete combustion of hydrogen produced by 100% zirconium-water reaction is a research hotspot in this field. Summary of the invention

[0004] The technical problem to be solved by the embodiments of the present application is to provide a method and device for determining the total atmospheric pressure in a containment vessel in view of the above-mentioned deficiencies in the prior art, which can accurately determine the first total atmospheric pressure in the containment vessel.

[0005] In a first aspect, an embodiment of the present application provides a method for determining the total atmospheric pressure in a containment vessel, comprising:

[0006] S110, obtaining a hydrogen concentration range in the containment;

[0007] S120. Determine the first atmospheric total pressure P in the containment according to the hydrogen concentration range. a2 .

[0008] Based on the same inventive concept, in a second aspect, an embodiment of the present application further provides a device for determining total atmospheric pressure, comprising:

[0009] The first acquisition module is used to obtain the hydrogen concentration range in the containment;

[0010] The first determination module is connected to the first acquisition module and is used to determine the first atmospheric total pressure P in the containment according to the hydrogen concentration range. a2 .

[0011] According to the method and device for determining the total atmospheric pressure in the containment provided in the embodiments of the present application, the first total atmospheric pressure in the containment is determined by the hydrogen concentration range. Compared with the conventional practice of using an integrated severe accident analysis program to calculate the typical severe accident sequence of a nuclear power plant to obtain a hydrogen release curve and a containment atmospheric pressure-temperature curve, and then using the hydrogen release curve and the containment atmospheric pressure-temperature curve as inputs and using a self-compiled program to calculate the total atmospheric pressure of the containment atmosphere after adiabatic isochoric complete combustion of hydrogen, the method of the present application breaks the conventional thinking and can accurately determine the first total atmospheric pressure in the containment. BRIEF DESCRIPTION OF THE DRAWINGS

[0012] Figure 1 A schematic flow chart showing a method for determining the total atmospheric pressure in a containment vessel provided in an embodiment of the present application;

[0013] Figure 2 A hydrogen flammability judgment diagram provided in an embodiment of the present application is shown;

[0014] Figure 3 Another schematic flow chart of a method for determining the total atmospheric pressure provided in an embodiment of the present application is shown;

[0015] Figure 4 A schematic structural diagram of a device for determining the total atmospheric pressure in a containment vessel provided in an embodiment of the present application is shown. DETAILED DESCRIPTION

[0016] In order to enable those skilled in the art to better understand the technical solution of the present application, the present application is further described in detail below with reference to the accompanying drawings and embodiments.

[0017] The features and exemplary embodiments of various aspects of the present application will be described in detail below. In order to make the purpose, technical solutions and advantages of the present application clearer, the present application will be further described in detail below in conjunction with the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only configured to explain the present application and are not configured to limit the present application. For those skilled in the art, the present application can be implemented without the need for some of these specific details. The following description of the embodiments is only to provide a better understanding of the present application by illustrating the examples of the present application.

[0018] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the statement "include..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0019] Example 1

[0020] Before introducing the method for determining the total atmospheric pressure in the containment vessel provided in the embodiment of the present application, the assumptions adopted in the embodiment of the present application are first explained.

[0021] The assumptions used in the embodiments of the present application include:

[0022] 1. The water vapor in the containment is always in a saturated state;

[0023] 2. All mixed gases in the containment are ideal gases and can be solved using the ideal gas state equation;

[0024] 3. Before the accident, the containment atmosphere included nitrogen N2, oxygen O2 and water vapor; after the accident, the active area fuel cladding was completely oxidized by water, and the generated hydrogen was completely released into the containment. The containment atmosphere included nitrogen N2, oxygen O2, water vapor and hydrogen H2;

[0025] 4. There was no leakage in the containment. The increase in atmospheric pressure in the containment after the accident was caused by the increase in atmospheric temperature and the release of water vapor and hydrogen.

[0026] The following describes a method for determining the total atmospheric pressure in a containment vessel provided in an embodiment of the present application.

[0027] Figure 1 A flow chart of a method for determining the total atmospheric pressure in a containment vessel provided in an embodiment of the present application is shown.

[0028] The method for determining the total atmospheric pressure in a containment provided in the embodiment of the present application can be executed by an electronic device and a device for determining the total atmospheric pressure in a containment, etc. The following description is given by taking the method for determining the total atmospheric pressure in a containment as executed by an electronic device.

[0029] like Figure 1 As shown, the method for determining the hydrogen concentration range provided in the embodiment of the present application may include steps S110 to S120.

[0030] S110. Obtain the hydrogen concentration range in the containment.

[0031] S120. Determine the first atmospheric total pressure P in the containment according to the hydrogen concentration range. a2 .

[0032] According to the method and device for determining the total atmospheric pressure in the containment provided in the embodiments of the present application, the first total atmospheric pressure in the containment is determined by the hydrogen concentration range. Compared with the conventional practice of using an integrated severe accident analysis program to calculate the typical severe accident sequence of a nuclear power plant to obtain a hydrogen release curve and a containment atmospheric pressure-temperature curve, and then using the hydrogen release curve and the containment atmospheric pressure-temperature curve as inputs and using a self-compiled program to calculate the total atmospheric pressure of the containment atmosphere after adiabatic isochoric complete combustion of hydrogen, the method of the present application breaks the conventional thinking and can accurately determine the first total atmospheric pressure in the containment.

[0033] The specific implementation of the above steps S110 to S120 is described below.

[0034] In step S110, the containment, namely the nuclear reactor containment, or reactor containment, containment building or enclosure, safety building, safety bunker, is the outermost building constituting the pressurized water reactor, and refers to the outer shell building that contains most of the systems and equipment of the nuclear steam supply system, and is used to accommodate the reactor pressure vessel and some safety systems (including the primary circuit main system and equipment, and the shutdown cooling system), completely isolating them from the external environment, in the hope of realizing the function of a safety protection barrier.

[0035] The hydrogen concentration range may include one of hydrogen being in the combustible zone, hydrogen being at the boundary between the combustible zone and the incombustible zone, and hydrogen being in the incombustible zone.

[0036] In some embodiments, the electronic device may store the hydrogen concentration range in the containment vessel, and the hydrogen concentration range in the containment vessel may be directly called later.

[0037] In some other implementations, step S110 may include:

[0038] S111. Obtain the first atmospheric temperature T1 of the containment atmosphere before hydrogen combustion after the accident, and the first nitrogen mole number N in the containment during normal operation n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 ;

[0039] S112, according to the first atmospheric temperature T1, the first nitrogen mole number N n20 、The first oxygen mole number N o20 and the first hydrogen mole number Nh20 , determine the first nitrogen partial pressure P in the containment after the accident n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 ;

[0040] S113, according to the first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 , determine the water vapor volume concentration f in the containment after the accident st and the hydrogen volume concentration f h2 ;

[0041] S114, according to the water vapor volume concentration f st and hydrogen volume concentration f h2 , determine the hydrogen concentration range in the containment.

[0042] In this embodiment, the first atmospheric temperature T1 of the containment atmosphere before the hydrogen combustion after the accident and the first nitrogen mole number N in the containment during normal operation are obtained. n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 Then, based on the above molar number and other data, determine the first nitrogen partial pressure P in the containment after the accident n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 , and then determine the water vapor volume concentration f in the containment after the accident based on the above partial pressure data st and hydrogen volume concentration f h2 , and then according to the water vapor volume concentration f st and hydrogen volume concentration f h2 In the embodiment of the present application, by determining the hydrogen concentration range based on the hydrogen flammability experimental data, it is possible to subsequently determine the total atmospheric pressure of the containment atmosphere simply and quickly.

[0043] In step S111 , the first atmospheric temperature T1 may include the atmospheric temperature of the containment atmosphere before the combustion of hydrogen after the accident.

[0044] The first nitrogen mole number N n20 May include the number of moles of nitrogen in the containment during normal operation.

[0045] The first oxygen mole number N o20 May include the number of moles of oxygen in the containment during normal operation.

[0046] The first hydrogen mole number N h20 May include the number of moles of hydrogen in the containment during normal operation.

[0047] Exemplarily, the first atmospheric temperature T1 may be pre-stored in the electronic device and may be directly called later. The value of the first atmospheric temperature T1 may be set according to actual conditions and is not limited here. For example, the value of the first atmospheric temperature T1 may be 100°C, 200°C, 350°C, etc.

[0048] In some embodiments, the first nitrogen mole number N in the containment during normal operation is obtained. n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 , that is, step S111 may include:

[0049] Obtain the second atmospheric total pressure P of the containment atmosphere during normal operation a0 、The number of moles of zirconium in the fuel cladding of the reactor active area N zr , and the average temperature T0 of the containment atmosphere when the water vapor in the containment is saturated and operating normally;

[0050] According to the average temperature T0, determine the water vapor partial pressure P in the containment during normal operation. st0 ;

[0051] According to the second atmospheric total pressure P a0 and water vapor partial pressure P st0 , determine the second nitrogen partial pressure P in the containment during normal operation n20 and the second oxygen partial pressure P o20 ;

[0052] According to the second nitrogen partial pressure P n20 , the second oxygen partial pressure P o20 and the average temperature T0, determine the first nitrogen mole number N n20 and the first oxygen mole number N o20 ;

[0053] According to the number of zirconium moles N zr , determine the first hydrogen mole number N h20 .

[0054] In this embodiment, the first nitrogen mole number N can be obtained by the above method. n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 , which can be used to determine the first nitrogen partial pressure P in the containment after the accident. n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure Ph21 and saturated water vapor partial pressure P st1 Provide a foundation.

[0055] Second atmospheric total pressure P a0 May include the total atmospheric pressure of the containment atmosphere during normal operation.

[0056] Zirconium mole number N zr The number of moles of zirconium in the reactor active fuel cladding may include the containment atmosphere during normal operation.

[0057] The average temperature T0 may include an average value of the temperature of the containment atmosphere when the water vapor in the containment is in a saturated state and operates normally.

[0058] It should be noted that, in the embodiment of the present application, the value of the first atmospheric temperature T1 may be the same as or different from the value of the average temperature T0, which is not limited here.

[0059] Water vapor partial pressure P st0 May include the partial pressure of water vapor in the containment during normal operation.

[0060] The second nitrogen partial pressure P n20 May include the nitrogen partial pressure in the containment during normal operation.

[0061] The second oxygen partial pressure P o20 May include the oxygen partial pressure in the containment during normal operation.

[0062] For example, the second atmospheric total pressure P can be obtained according to the target nuclear power plant reactor and containment design data. a0 and average temperature T0, wherein the target nuclear power plant reactor can be any nuclear power plant reactor. Alternatively, the second atmospheric total pressure P a0 The average temperature T0 can be pre-stored in the electronic device for subsequent direct retrieval.

[0063] In some examples, the number of moles of zirconium N in the reactor active region fuel cladding is obtained. zr , which may include:

[0064] Obtain the total mass M of zirconium in the fuel cladding of the reactor active area zr and zirconium molar mass M zrm ;

[0065] The total mass of zirconium M zr and zirconium molar mass M zrm Substituting into formula (1) for calculation, we can obtain the molar number of zirconium in the fuel cladding of the reactor active area: N zr ;

[0066] Formula (1) includes:

[0067] Nzr =M zr / M zrm (1).

[0068] In this example, the total mass of zirconium M is zr and zirconium molar mass M zrm Substituting into formula (1), the number of zirconium moles N can be quickly and accurately determined zr .

[0069] For example, the total mass of zirconium M can be determined by the design data of the target nuclear power plant reactor and containment. zr Or, the total mass of zirconium M zr It can be stored in the electronic device in advance for subsequent direct retrieval.

[0070] For example, the zirconium molar mass M can be determined from the periodic table: zrm Alternatively, the zirconium molar mass M zrm It can be stored in the electronic device in advance for subsequent direct retrieval.

[0071] In some examples, the water vapor partial pressure P in the containment during normal operation is determined based on the average temperature T0. st0 ,include:

[0072] The average temperature T0 is calculated using the international water and steam property calculation model IAPWS-IF97, and the water vapor partial pressure P in the containment during normal operation is obtained. st0 .

[0073] In this embodiment, since it is assumed that the water vapor in the containment is in a saturated state, the average temperature T0 can be calculated by the international water and water vapor property calculation model IAPWS-IF97, and the water vapor partial pressure P in the containment during normal operation can be obtained quickly and accurately. st0 .

[0074] The International Water and Steam Properties Calculation Model IAPWS-IF97 (International Association for the Properties of Water and Steam, Industrial Formulation 1997) is an international standard property calculation model used to calculate various thermodynamic and transport properties of water and water vapor. This model is developed and maintained by the International Association for the Properties of Water and Steam (IAPWS). IAPWS-IF97 provides an accurate method to calculate the thermodynamic properties of water and water vapor, including density, specific volume, enthalpy, entropy, specific enthalpy and specific entropy.

[0075] In some examples, according to the second atmospheric total pressure P a0 and water vapor partial pressure P st0 , determine the second nitrogen partial pressure P in the containment during normal operation n20 and the second oxygen partial pressure P o20 , which may include:

[0076] The second atmospheric total pressure P a0 and water vapor partial pressure P st0 Substitute into formula (2) to calculate and obtain the second nitrogen partial pressure P in the containment during normal operation: n20 ;

[0077] The second atmospheric total pressure P a0 and water vapor partial pressure P st0 Substitute into formula (3) to calculate and obtain the second oxygen partial pressure P in the containment during normal operation: o20 ;

[0078] Formula (2) includes:

[0079] P n20 =0.79×(P a0 -P st0 ) (2)

[0080] Formula (3) includes:

[0081] P o20 =0.21×(P a0 -P st0 ) (3).

[0082] In this example, the second nitrogen partial pressure P can be determined quickly and accurately by using the above formulas (2) and (3). n20 and the second oxygen partial pressure P o20 .

[0083] In some examples, according to the second nitrogen partial pressure P n20 , the second oxygen partial pressure P o20 and the average temperature T0, determine the first nitrogen mole number N n20 and the first oxygen mole number N o20 ,include:

[0084] Get the ideal gas constant R and the free volume V of the containment c ;

[0085] The ideal gas constant R, the free volume of the containment V c , the second nitrogen partial pressure P n20 Substitute the average temperature T0 into formula (4) to calculate and obtain the first nitrogen mole number Nn20 ;

[0086] According to the ideal gas constant R, the free volume V of the containment c , the second oxygen partial pressure P o20 Substitute the average temperature T0 into formula (5) to calculate the first oxygen mole number N o20 ;

[0087] Formula (4) includes:

[0088] N n20 =P n20 V c / (RT0) (4)

[0089] Formula (5) includes:

[0090] N o20 =P o20 V c / (RT0) (5).

[0091] In this example, the first nitrogen mole number N can be quickly and accurately determined by the above formula (4) and formula (5): n20 and the first oxygen mole number N o20 .

[0092] The free volume of the containment can include the space in the containment that is not occupied by solid matter (concrete structure, steel structure, equipment, etc.) and where air can freely enter and exit. In the initial stage of an accident, the free volume of the containment is the main carrier of the accident energy.

[0093] For example, the ideal gas constant R may be pre-stored in the electronic device and may be directly called later. The value of the ideal gas constant R may be 8.314 J / (mol·K).

[0094] For example, the free volume V of the containment can be determined by the target nuclear power plant reactor and containment design data. c Alternatively, the containment free volume V c It can be stored in the electronic device in advance and can be directly called later.

[0095] In some examples, according to the number of moles of zirconium N zr , determine the first hydrogen mole number N h20 , which may include:

[0096] The number of moles of zirconium N zr Substitute into formula (6) to calculate and obtain the first hydrogen mole number N h20 ;

[0097] Formula (6) includes:

[0098] Nh20 =2×N zr (6).

[0099] In this example, the first hydrogen mole number N can be quickly and accurately determined by formula (6): h20 .

[0100] It should be noted that formula (6) can be derived from the chemical equation Zr+2H2O→ZrO2+2H2, assuming that all the hydrogen is released into the containment after a serious accident. Wherein, Zr represents zirconium, and ZrO2 represents zirconium dioxide, commonly known as zirconium oxide.

[0101] In step S112, the electronic device obtains the first atmospheric temperature T1 of the containment atmosphere before the hydrogen combustion after the accident, the first nitrogen mole number N in the containment during normal operation, and the n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 Then, according to the first atmospheric temperature T1, the first nitrogen mole number N n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 , determine the first nitrogen partial pressure P in the containment after the accident n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 .

[0102] The first nitrogen partial pressure P n21 May include the nitrogen partial pressure in the containment after an accident.

[0103] The first oxygen partial pressure P o21 May include the oxygen partial pressure in the containment after an accident.

[0104] The first hydrogen partial pressure P h21 May include the hydrogen partial pressure in the containment after an accident.

[0105] Saturated water vapor partial pressure P st1 It may include the partial pressure of saturated water vapor in the containment after the accident.

[0106] In some embodiments, according to the first atmospheric temperature T1, the first nitrogen mole number N n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 , determine the first nitrogen partial pressure P in the containment after the accident n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure Pst1 , that is, step S112 may include:

[0107] Get the ideal gas constant R and the free volume V of the containment c ;

[0108] The first nitrogen mole number N n20 , ideal gas constant R, first atmosphere temperature T1 and containment free volume V c Substitute into formula (7) to calculate and obtain the first nitrogen partial pressure P in the containment after the accident: n21 ;

[0109] The first oxygen mole number N o20 , ideal gas constant R, first atmosphere temperature T1 and containment free volume V c Substitute into formula (8) to calculate and obtain the first oxygen partial pressure P in the containment after the accident: o21 ;

[0110] The first hydrogen mole number N h20 , ideal gas constant R, first atmosphere temperature T1 and containment free volume V c Substitute into formula (9) to calculate and obtain the first hydrogen partial pressure P in the containment after the accident: h21 ;

[0111] The first atmospheric temperature T1 was calculated using the international water and water vapor property calculation model IAPWS-IF97, and the saturated water vapor partial pressure P in the containment after the accident was obtained. st1 ;

[0112] Among them, formula (7) includes:

[0113] P n21 =N n20 RT1 / V c (7)

[0114] Formula (8) includes:

[0115] P o21 =N o20 RT1 / V c (8)

[0116] Formula (9) includes:

[0117] P h21 =N h20 RT1 / V c (9).

[0118] In this example, the first nitrogen partial pressure P can be determined quickly and accurately by using formulas (7) to (9). n21 , the first oxygen partial pressure P o21, the first hydrogen partial pressure P h21 ; The saturated water vapor partial pressure P can be determined quickly and accurately through the international water and water vapor physical property calculation model IAPWS-IF97 st1 .

[0119] The above formulas (7) to (9) can be obtained from the ideal gas state equation PV=nRT, where P represents pressure, V represents volume, n represents the number of moles, R represents the molar gas constant, that is, the ideal gas constant mentioned above, and T represents temperature.

[0120] In step S113, the electronic device generates a temperature signal according to the first atmospheric temperature T1 and the first nitrogen mole number N n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 , determine the first nitrogen partial pressure P in the containment after the accident n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 Then, according to the first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 , determine the water vapor volume concentration f in the containment after the accident st and the hydrogen volume concentration f h2 .

[0121] Water vapor volume concentration f st May include the volume concentration of water vapor in the containment after the accident.

[0122] Hydrogen volume concentration f h2 May include the volume concentration of hydrogen in the containment after the accident.

[0123] In some embodiments, according to the first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 , determine the water vapor volume concentration f in the containment after the accident st and the hydrogen volume concentration f h2 , which may include:

[0124] The first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1Substitute into formula (10) to calculate and obtain the third atmospheric total pressure P of the containment atmosphere after the accident: a1 ;

[0125] The saturated water vapor partial pressure P st1 and the third atmospheric total pressure P a1 Substitute into formula (11) to calculate and obtain the water vapor volume concentration f in the containment after the accident st ;

[0126] The first hydrogen partial pressure P h21 and the third atmospheric total pressure P a1 Substitute into formula (12) to calculate and obtain the hydrogen volume concentration f in the containment after the accident: h2 ;

[0127] Formula (10) includes:

[0128] P a1 =P n21 +P o21 +P st1 +P h21 (10)

[0129] Formula (11) includes:

[0130] f st =P st1 / P a1 (11)

[0131] Formula (12) includes:

[0132] f h2 =P h21 / P a1 (12).

[0133] In this embodiment, the water vapor volume concentration f can be quickly and accurately determined by the above formulas (10) to (12). st and the hydrogen volume concentration f h2 .

[0134] In step S114, the electronic device generates a pressure drop according to the first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 , determine the water vapor volume concentration f in the containment after the accident st and the hydrogen volume concentration f h2 Afterwards, the water vapor volume concentration f st and the hydrogen volume concentration f h2 , determine the hydrogen concentration range in the containment.

[0135] Figure 2 A hydrogen flammability determination diagram provided by an embodiment of the present application is shown. Figure 2 As shown, Figure 2 The horizontal axis in the figure can represent the hydrogen volume share, the vertical axis can represent the water vapor volume share, the gray area can represent the hydrogen flammable area, and the white area can represent the hydrogen non-flammable area. Figure 2 The boundary line X between the hydrogen flammable zone and the hydrogen non-flammable zone can be described by the formula (27) summarized from the experiment. Formula (27) includes:

[0136] f st '=1-(f h2 +0.373exp(-0.7f h2 ')+5.18exp(-48.8f h2 ')) (27)

[0137] Among them, f st ' can represent the volume fraction of water vapor; f h2 ' can represent the volume fraction of hydrogen.

[0138] Based on the above analysis, in some embodiments, according to the water vapor volume concentration f st and hydrogen volume concentration f h2 , determine the hydrogen concentration range in the containment, that is, step S114 may include:

[0139] In f st <1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )) In the case of hydrogen concentration in the containment, it is determined that the hydrogen is in the flammable zone;

[0140] Or, in f st =1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )) In the case of, the hydrogen concentration range in the containment is determined as hydrogen in the containment at the junction of the combustible zone and the non-combustible zone;

[0141] Or, in f st >1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )) is carried out, the hydrogen concentration range in the containment is determined to be that hydrogen is in the non-flammable zone.

[0142] In this embodiment, by fst and 1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )) can accurately and quickly determine the hydrogen concentration range.

[0143] It is understandable that in f st =1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )), the hydrogen concentration range in the containment is determined to be that the hydrogen is at the junction of the combustible zone and the non-combustible zone in the containment, and the total atmospheric pressure of the containment atmosphere is at the maximum value at which hydrogen can be ignited.

[0144] In some embodiments, after determining that the hydrogen concentration range in the containment is in the flammable zone, the method may further include:

[0145] S115, updating the first atmospheric temperature T1 according to a preset step size;

[0146] S116. Repeat steps S112 to S115 on an iterative basis until it is determined that the hydrogen concentration range in the containment is at the junction of the combustible zone and the incombustible zone, and the iterative first atmospheric temperature T1 is obtained.

[0147] In this embodiment, through the above method, it is possible to determine that the hydrogen concentration range in the containment is the first atmospheric temperature T1 when the hydrogen is at the junction of the combustible zone and the incombustible zone, that is, the iterated first atmospheric temperature T1, and also the first atmospheric temperature T1 corresponding to the limit state point of the combustible state of hydrogen can be determined, providing a basis for the subsequent determination of the first atmospheric total pressure.

[0148] The preset step size can be set according to actual conditions and is not limited here. For example, the preset step size can be 1K°C.

[0149] It should be noted that step S112 may be the beginning of a cyclic search for the containment atmosphere state before hydrogen combustion.

[0150] It is understandable that in f st <1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2)), it is determined that the hydrogen concentration range in the containment is that hydrogen is in the flammable zone, and then the first atmospheric temperature of the containment atmosphere after the accident can be further increased to increase the atmospheric pressure of the containment before the hydrogen combustion, thereby obtaining a conservative total atmospheric pressure of the containment atmosphere after the hydrogen is completely burned. The first atmospheric temperature T1 is updated according to the preset step size, and steps S112 to S115 are repeatedly executed on an iterative basis until it is determined that the hydrogen concentration range in the containment is that hydrogen is at the junction of the flammable zone and the incombustible zone, and the iterative first atmospheric temperature T1 is obtained.

[0151] In S120, after obtaining the hydrogen concentration range in the containment, the electronic device can also determine the first atmospheric total pressure P in the containment according to the hydrogen concentration range. a2 .

[0152] In some embodiments, the first atmospheric total pressure P in the containment is determined according to the hydrogen concentration range. a2 , i.e. step S120, may include:

[0153] S121, when the hydrogen concentration in the containment is at the boundary between the combustible zone and the incombustible zone, determining a first atmospheric temperature T1 corresponding to the hydrogen concentration range, and a second atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned;

[0154] S122, determine the first atmospheric total pressure P in the containment according to the first atmospheric temperature T1 and the second atmospheric temperature T2 a2 .

[0155] In this embodiment, by determining the first atmospheric temperature T1 corresponding to the hydrogen range concentration when the hydrogen is at the boundary of the combustible zone and the non-combustible zone in the containment, and the second atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned, the first atmospheric total pressure P in the containment can be accurately determined according to the first atmospheric temperature T1 and the second atmospheric temperature T2. a2 .

[0156] The second atmospheric temperature T2 may include the assumed atmospheric temperature of the containment atmosphere after the hydrogen is completely burned. The value of the second atmospheric temperature T2 may be set according to actual conditions and is not limited here.

[0157] It should be noted that, when the hydrogen concentration range in the containment is such that hydrogen is at the boundary between the combustible zone and the incombustible zone, the first atmospheric temperature T1 corresponding to the hydrogen concentration range may be the first atmospheric temperature T1 after iteration in step S116 in Example 1.

[0158] Exemplarily, the value of the second atmospheric temperature T2 may be pre-stored in the electronic device and may be directly called later, or may be obtained in response to an input operation of the user based on the electronic device. The input operation may include an operation of outputting the value of the second atmospheric temperature T2, for example, the input operation may be an operation of directly inputting the value of the second atmospheric temperature T2 in an input box displayed by the electronic device. For another example, the input operation may be an operation of selecting the value of the second atmospheric temperature T2 from a plurality of temperature options displayed by the electronic device.

[0159] In some embodiments, the first atmospheric total pressure P in the containment is determined according to the first atmospheric temperature T1 and the second atmospheric temperature T2. a2 , i.e. step S122, may include:

[0160] S1221. Obtain first mass information, second mass information, first specific heat capacity at constant volume information, second specific heat capacity at constant volume information, and first hydrogen mole number N h20 The total heat released by complete combustion is Q all ;

[0161] S1222, based on the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, the total heat Q all and the first atmospheric temperature T1, determine the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * ;

[0162] S1223, according to the second atmospheric temperature T2 and the third atmospheric temperature T2 * , determine the total pressure of the first atmosphere in the containment after the adiabatic isochoric complete combustion of hydrogen P a2 .

[0163] In this embodiment, the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, the total heat Q all and the first atmospheric temperature T1, determine the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * , and then according to the second atmospheric temperature T2 and the third atmospheric temperature T2 * , determine the total pressure of the first atmosphere in the containment after the adiabatic isochoric complete combustion of hydrogen P a2 , that is, based on the experimental data of hydrogen flammability; then search for the first atmospheric temperature T1 (that is, the limit state point of hydrogen flammability) corresponding to the junction of hydrogen in the flammable zone and the inflammable zone; then perform adiabatic isochoric complete combustion of hydrogen at this state point to obtain the first atmospheric total pressure P after combustion a2 In this way, there is no need to perform severe accident sequence calculations, which is simple and quick.

[0164] In step S1221, the total heat Qall may include a first hydrogen mole number N h20 Complete combustion releases all the heat.

[0165] The first mass information may include the first nitrogen mass M in the containment during normal operation. n20 、First oxygen mass M o20 、First hydrogen mass M h20 and the first water vapor mass M st1 .

[0166] The second mass information includes the second nitrogen mass M in the containment after the hydrogen is completely burned n22 、The second oxygen mass M o22 and the second water vapor mass M st2 .

[0167] The first constant volume specific heat capacity information includes the first constant volume specific heat capacity CV of nitrogen at the first atmospheric temperature T1 n21 、First oxygen constant volume specific heat capacity CV o21 、The first hydrogen constant volume specific heat capacity CV h21 and the first water vapor constant volume specific heat CV st1 .

[0168] The second constant volume specific heat capacity information includes the first constant volume specific heat capacity CV of nitrogen at the second atmospheric temperature T2 n22 、Second oxygen constant volume specific heat capacity CV o22 、Second specific heat capacity of hydrogen at constant volume CV h22 and the second water vapor constant volume specific heat capacity CV st2 .

[0169] In some embodiments, the first mass information may include the first nitrogen mass M in the containment during normal operation. n20 、First oxygen mass M o20 、First hydrogen mass M h20 and the first water vapor mass M st1 ;

[0170] Acquiring the first quality information may include:

[0171] Get the molar mass of nitrogen M n2m 、Oxygen molar mass M o2m 、Hydrogen molar mass M h2m , water vapor molar mass M stm 、The first nitrogen mole number N n20 、The first oxygen mole number N o20 、The first hydrogen mole number N h20 、Free volume of containment V c , and the density of saturated water vapor in the containment ρ st1 ;

[0172] The molar mass of nitrogen M n2m and the first nitrogen mole number N n20 Substitute into formula (15) to calculate and obtain the first nitrogen mass M in the containment during normal operation: n20 ;

[0173] The molar mass of oxygen M o2m and the first oxygen mole number N o20 Substitute into formula (16) to calculate and obtain the first oxygen mass M in the containment during normal operation: o20 ;

[0174] The molar mass of hydrogen M h2m and the first hydrogen mole number N h20 Substitute into formula (17) to calculate and obtain the first hydrogen mass M in the containment during normal operation: h20 ;

[0175] The free volume of the containment V c and density ρ st1 Substitute into formula (18) to calculate and obtain the first steam mass M in the containment during normal operation: st1 ;

[0176] Formula (15) includes:

[0177] M n20 =M n2m N n20 (15)

[0178] Formula (16) includes:

[0179] M o20 =M o2m N o20 (16)

[0180] Formula (17) includes:

[0181] M h20 =M h2m N h20 (17)

[0182] Formula (18) includes:

[0183] M st1 =ρ st1 V c (18).

[0184] In this embodiment, the first mass information can be determined quickly and accurately through the above formulas (15) to (18), which is used for the subsequent determination of the first atmospheric total pressure P a2 Provide the foundation.

[0185] The first nitrogen mass M n20 May include the mass of nitrogen in the containment during normal operation.

[0186] The first oxygen mass M o20 May include the mass of oxygen in the containment during normal operation.

[0187] The first hydrogen mass M h20 May include the mass of hydrogen in the containment during normal operation.

[0188] Molar mass of nitrogen M n2m The molar mass of nitrogen may be included.

[0189] Oxygen molar mass M o2m The molar mass of oxygen may be included.

[0190] Molar mass of hydrogen M h2m The molar mass of hydrogen may be included.

[0191] Molar mass of water vapor M stm The molar mass of water vapor may be included.

[0192] Density ρ st1 May include the density of saturated water vapor in the containment after an accident.

[0193] The first water vapor mass M st1 May include the mass of water vapor in the containment during normal operation.

[0194] Exemplarily, the nitrogen molar mass M can be obtained according to the periodic table: n2m 、Oxygen molar mass M o2m 、Hydrogen molar mass M h2m and the water vapor molar mass M stm Alternatively, the nitrogen molar mass M can be pre-stored in the electronic device. n2m 、Oxygen molar mass M o2m 、Hydrogen molar mass M h2m and the water vapor molar mass M stm , so that it can be called directly later.

[0195] For example, the first atmospheric temperature T1 can be calculated using the international water and water vapor property calculation model IAPWS-IF97 to obtain the density of saturated water vapor in the containment: st1 Alternatively, the density of saturated water vapor in the containment can be pre-stored in the electronic device: st1 , so that it can be called directly later.

[0196] In some embodiments, the second mass information may include the second nitrogen mass M in the containment after the hydrogen is completely burned. n22 、The second oxygen mass Mo22 and the second water vapor mass M st2 ;

[0197] Obtain second quality information, including:

[0198] Get the first nitrogen mole number N n20 、The first oxygen mole number N o20 , the first hydrogen mole number N h20 , Nitrogen molar mass M n2m 、Oxygen molar mass M o2m , water vapor molar mass M stm , and the first steam mole number N in the containment during normal operation st1 ;

[0199] The first nitrogen mole number N n20 Substitute into formula (19) to calculate and obtain the second nitrogen mole number N in the containment after the hydrogen is completely burned: n22 ;

[0200] The first oxygen mole number N o20 and the first hydrogen mole number N h20 Substitute into formula (20) to calculate and obtain the second oxygen mole number N in the containment after the hydrogen is completely burned: o22 ;

[0201] The first water vapor mole number N st1 and the first hydrogen mole number N h20 Substitute into formula (21) to calculate and obtain the second water vapor mole number N in the containment after the hydrogen is completely burned: st2 ;

[0202] The second nitrogen mole number N n22 Substitute the molar mass of nitrogen into formula (22) to calculate the second nitrogen mass M in the containment after the hydrogen is completely burned. n22 ;

[0203] The second oxygen mole number N o22 and oxygen molar mass M o2m Substitute into formula (23) to calculate and obtain the second oxygen mass M in the containment after the hydrogen is completely burned: o22 ;

[0204] The second water vapor mole number N st2 and the water vapor molar mass M stm Substituting into formula (24) for calculation, we can obtain the second steam mass M in the containment after the hydrogen is completely burned: st2 ;

[0205] Formula (19) includes:

[0206] Nn22 =N n20 (19)

[0207] Formula (20) includes:

[0208] N o22 =N o20 -N h20 / 2 (20)

[0209] Formula (21) includes:

[0210] N st2 =N st1 +N h20 (twenty one)

[0211] Formula (22) includes:

[0212] M n22 =M n2m N n22 (twenty two)

[0213] Formula (23) includes:

[0214] M o22 =M o2m N o22 (twenty three)

[0215] Formula (24) includes:

[0216] M st2 =M stm N st2 (twenty four).

[0217] In this embodiment, the second quality information can be determined quickly and accurately through the above formulas (19) to (24).

[0218] In some examples, the first water vapor mole number N in the containment during normal operation is obtained. st1 , which may include:

[0219] Get the first water vapor mass M st1 and the water vapor molar mass M stm ;

[0220] The first water vapor mass M st1 and the water vapor molar mass M stm Substitute into formula (25) to calculate and obtain the first water vapor mole number N in the containment during normal operation: st1 ;

[0221] Formula (25) includes:

[0222] N st1 =M st1 / Mstm (25).

[0223] In this embodiment, by setting the first water vapor mass M st1 and the water vapor molar mass M stm Substituting into formula (25), the first water vapor mole number N can be quickly and accurately determined st1 .

[0224] Exemplarily, the water vapor molar mass M stm It can be calculated by the above formula (18). Alternatively, the electronic device can pre-store the water vapor molar mass M stm , which can be called directly later.

[0225] In some embodiments, obtaining a first hydrogen mole number N h20 The total heat released by complete combustion is Q all , which may include:

[0226] Obtain the combustion heat Q of hydrogen and the first number of hydrogen moles N h20 ;

[0227] The heat of combustion Q and the first number of hydrogen moles N h20 Substitute into formula (26) to calculate and obtain the first hydrogen mole number N h20 The total amount of heat released by complete combustion;

[0228] Formula (26) includes:

[0229] Q all =N h20 Q (26).

[0230] In this embodiment, the heat of combustion Q and the first number of hydrogen moles N are h20 Substituting into formula (26), the total heat Q can be quickly and accurately determined all .

[0231] For example, the combustion heat Q can be pre-stored in the electronic device and can be directly called later. For example, the combustion heat Q of 1 mol of hydrogen is 285.8 kJ / mol.

[0232] It can be understood that according to the hydrogen combustion chemical equation: The above formula (26) can be derived.

[0233] For example, the first specific heat capacity at constant volume information and the second specific heat capacity at constant volume information may be determined according to the fluid property table. Alternatively, the first specific heat capacity at constant volume information and the second specific heat capacity at constant volume information may be pre-stored in the electronic device for subsequent direct call.

[0234] In step S1222, the electronic device obtains the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, and the first hydrogen mole number N h20 The total heat released by complete combustion is Q all Afterwards, the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, the total heat Q all and the first atmospheric temperature, determine the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * .

[0235] The third atmospheric temperature T2 * May include the actual atmospheric temperature of the containment atmosphere after complete combustion of the hydrogen.

[0236] In some embodiments, the first mass information may include the first nitrogen mass M in the containment during normal operation. n20 、First oxygen mass M o20 、First hydrogen mass M h20 and the first water vapor mass M st1 ;

[0237] The second mass information includes the second nitrogen mass M in the containment after the hydrogen is completely burned n22 、The second oxygen mass M o22 and the second water vapor mass M st2 ;

[0238] The first constant volume specific heat capacity information includes the first constant volume specific heat capacity CV of nitrogen at the first atmospheric temperature T1 n21 、First oxygen constant volume specific heat capacity CV o21 、The first hydrogen constant volume specific heat capacity CV h21 and the first water vapor constant volume specific heat CV st1 ;

[0239] The second constant volume specific heat capacity information includes the first constant volume specific heat capacity CV of nitrogen at the second atmospheric temperature T2 n22 、Second oxygen constant volume specific heat capacity CV o22 、Second specific heat capacity of hydrogen at constant volume CV h22 and the second water vapor constant volume specific heat capacity CV st2 ;

[0240] According to the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, the total heat Q all and the first atmospheric temperature T1, determine the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * , that is, step S1222 may include:

[0241] The first nitrogen mass M n20、First oxygen mass M o20 、First hydrogen mass M h20 、First water vapor mass M st1 、The second nitrogen quality M n22 、The second oxygen mass M o22 、The second water vapor mass M st2 、First nitrogen constant volume specific heat capacity CV n21 、First oxygen constant volume specific heat capacity CV o21 、The first hydrogen constant volume specific heat capacity CV h21 、First water vapor constant volume specific heat capacity CV st1 、First nitrogen constant volume specific heat capacity CV n22 、Second oxygen constant volume specific heat capacity CV o22 、Second specific heat capacity of hydrogen at constant volume CV h22 、Second specific heat capacity of water vapor at constant volume CV st2 Total calories Q all Substitute the first atmospheric temperature into formula (13) to calculate and obtain the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * ;

[0242] Formula (13) includes:

[0243] T2 * =(Q all +T1(M n20 CV n21 +M o20 CV o21 +M h20 CV h21 +M st1 CV st1 )) / (M n22 CV n22 +M o22 CV o22 +M st2 CV st2 ) (13).

[0244] In this embodiment, the third atmospheric temperature T2 can be determined quickly and accurately by the above formula (13): * .

[0245] It should be noted that the derivation principle of the above formula (13) is the energy conservation formula.

[0246] In step S1223, the electronic device calculates the total heat Q according to the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, and the total heat Q all and the first atmospheric temperature, determine the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned *Afterwards, the second atmospheric temperature T2 and the third atmospheric temperature T2 * , determine the total pressure P of the first atmosphere in the containment atmosphere after the adiabatic isochoric complete combustion of hydrogen a2 .

[0247] In some embodiments, according to the second atmospheric temperature T2 and the third atmospheric temperature T2 * , determine the first atmospheric total pressure P of the containment atmosphere after adiabatic isochoric complete combustion of hydrogen a2 , that is, step S1223 may include:

[0248] S1, in abs(T2-T2 * ) / T2≤N, according to the second nitrogen mole number N in the containment after the hydrogen is completely burned n22 、The second oxygen mole number N o22 、The second water vapor mole number N st2 and the second atmosphere temperature T2, determine the total pressure P of the first atmosphere in the containment atmosphere after the adiabatic isochoric complete combustion of hydrogen a2 ; or, in abs(T2-T2 * ) / T2>N, let T2=(T2+T2 * ) / 2, and repeatedly perform step S120 to step S1 on an iterative basis until the total pressure of the first atmosphere in the containment atmosphere after the adiabatic isochoric complete combustion of hydrogen is determined to be P a2 ;

[0249] Wherein, N is a positive number.

[0250] In this embodiment, by abs(T2-T2 * The relationship between ) / T2 and N can accurately determine the total pressure P of the first atmosphere in the containment atmosphere after the adiabatic isochoric complete combustion of hydrogen. a2 .

[0251] That is, if abs(T2-T2 * ) / T2≤N, the iteration converges, and according to the second atmospheric temperature T2 after the hydrogen is completely burned and the ideal gas state equation, the first atmospheric total pressure P in the containment atmosphere after the adiabatic isochoric complete combustion of hydrogen is determined. a2 ; If abs(T2-T2 * ) / T2>N, then update the value of the second atmospheric temperature T2, let T2=(T2+T2 * ) / 2, and repeatedly perform step S120 to step S1 on an iterative basis until the total pressure of the first atmosphere in the containment atmosphere after the adiabatic isochoric complete combustion of hydrogen is determined to be P a2The acquisition of the second constant volume specific heat capacity information in step S1221 can be used as the start of searching for the second atmospheric temperature T2 in the containment vessel after the hydrogen is completely burned, and the second atmospheric temperature T2 can be updated with the iterative process.

[0252] It should be noted that the value of N can be set according to actual conditions and is not limited here. For example, N can be 0.01.

[0253] In some examples, the second nitrogen mole number N in the containment after the hydrogen is completely burned is n22 、The second oxygen mole number N o22 、The second water vapor mole number N st2 and the second atmosphere temperature T2, determine the total pressure P of the first atmosphere in the containment atmosphere after the adiabatic isochoric complete combustion of hydrogen a2 ,include:

[0254] Get the ideal gas constant R and the free volume V of the containment c ;

[0255] The second nitrogen mole number N n22 、The second oxygen mole number N o22 、The second water vapor mole number N st2 , the second atmosphere temperature T2, the ideal gas constant R and the free volume V of the containment c Substituting into formula (14) for calculation, we can obtain the first atmospheric total pressure P of the containment atmosphere after adiabatic isochoric complete combustion of hydrogen: a2 ;

[0256] Formula (14) includes:

[0257] P a2 =(N n22 +N o22 +N st2 )RT2 / V c (14).

[0258] In this example, the first atmospheric total pressure P can be quickly and accurately determined by the above formula (14): a2 .

[0259] In order to better understand the method for determining the total atmospheric pressure provided in the embodiment of the present application, a specific embodiment is described below.

[0260] like Figure 3 As shown, the method for determining the total atmospheric pressure provided in the embodiment of the present application may include steps S300 to S310.

[0261] S300, consult the nuclear power plant design data to determine the zirconium mass (i.e. the total zirconium mass M zr ), the atmospheric pressure of the containment during normal operation (i.e. the second atmospheric total pressure Pa0 ) and temperature (i.e. average temperature T0), free volume (i.e. containment free volume V c ).

[0262] S301, calculate the hydrogen production after the active area fuel cladding is completely oxidized (i.e. the total number of hydrogen moles N h20 and hydrogen mass M h20 ).

[0263] S302, calculate the content of each component of the containment atmosphere during normal operation (i.e. the first nitrogen mole number N n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 wait.

[0264] S303. Assume / update the atmospheric temperature T1 before hydrogen combustion after the accident (i.e., the first atmospheric temperature T1), and calculate the content of each component of the containment atmosphere.

[0265] S304, determine whether the hydrogen is at the flammable boundary (i.e. determine that the hydrogen concentration range in the containment is such that the hydrogen is at the boundary between the flammable zone and the non-flammable zone).

[0266] S305: If no, update the atmospheric temperature T1 before hydrogen combustion.

[0267] S306: If yes, assume / update the atmospheric temperature T1 after hydrogen combustion.

[0268] S307. According to the law of energy conservation, calculate the atmospheric temperature T2 after hydrogen combustion * (i.e. the third atmospheric temperature T2 * ).

[0269] S308, determine abs(T2-T2 * ) / T2≤N holds true.

[0270] S309: If not established, update the atmospheric temperature T2 after hydrogen combustion.

[0271] S310, if established, calculate the atmospheric pressure after hydrogen combustion (i.e., the first atmospheric total pressure P) according to the ideal gas state equation a2 ).

[0272] The method for determining the total atmospheric pressure in the containment provided in the embodiment of the present application is based on the experimental data of hydrogen flammability, searches for the limit state point of the hydrogen flammability state, performs adiabatic isochoric complete combustion of hydrogen at this state point, and obtains the atmospheric pressure after combustion. This method does not require severe accident sequence calculation and is simple and fast.

[0273] The method for determining the total atmospheric pressure in the containment provided in the embodiment of the present application can be used as a rapid assessment method for the atmospheric pressure in the containment after adiabatic isochoric complete combustion of hydrogen after a severe accident, providing support for hydrogen risk assessment in nuclear power plants.

[0274] Based on the same inventive concept as the method for determining the total atmospheric pressure in a containment vessel, an embodiment of the present application further provides a device for determining the total atmospheric pressure in a containment vessel. The method for determining the total atmospheric pressure in a containment vessel is introduced below in conjunction with the accompanying drawings.

[0275] Example 2

[0276] Figure 4 It is a structural schematic diagram of a device for determining the total atmospheric pressure in a containment vessel provided in an embodiment of the present application.

[0277] like Figure 4 As shown, the device for determining the total atmospheric pressure in the containment provided in the embodiment of the present application may include a first acquisition module 410 and a first determination module 420 .

[0278] The first acquisition module 410 is used to obtain the concentration range of hydrogen in the containment;

[0279] The first determination module 420 is connected to the first acquisition module 410 and is used to determine the first atmospheric total pressure P in the containment according to the hydrogen concentration range. a2 .

[0280] According to the device for determining the total atmospheric pressure in the containment provided in the embodiment of the present application, the first total atmospheric pressure in the containment is determined by the hydrogen concentration range. Compared with the conventional practice of using an integrated severe accident analysis program to calculate the typical severe accident sequence of a nuclear power plant to obtain a hydrogen release curve and a containment atmospheric pressure-temperature curve, and then using the hydrogen release curve and the containment atmospheric pressure-temperature curve as inputs and using a self-compiled program to calculate the total atmospheric pressure of the containment atmosphere after adiabatic isochoric complete combustion of hydrogen, the method of the present application breaks the conventional thinking and can accurately determine the first total atmospheric pressure in the containment.

[0281] In some implementations, the first determining module 420 may include:

[0282] The first determination submodule is used to determine the first atmospheric temperature T1 corresponding to the hydrogen range concentration when the hydrogen range concentration in the containment is at the junction of the combustible zone and the non-combustible zone, and the second atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned;

[0283] The second determination submodule is connected to the first determination submodule and is used to determine the first atmospheric total pressure P in the containment according to the first atmospheric temperature T1 and the second atmospheric temperature T2. a2 .

[0284] In some implementations, the second determining submodule includes:

[0285] An acquisition unit, used to acquire the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, and the first hydrogen mole number N h20 The total heat released by complete combustion is Q all ;

[0286] The first determining unit is connected to the acquiring unit and is used to determine the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, the total heat Q all and the first atmospheric temperature T1, determine the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * ;

[0287] The second determining unit is connected to the first determining unit and is used to determine the temperature of the second atmosphere T2 and the third atmosphere T2 according to the second atmosphere temperature T2 * , determine the total pressure of the first atmosphere in the containment after the adiabatic isochoric complete combustion of hydrogen P a2 .

[0288] In some implementations, the first acquisition module 410 may further include:

[0289] The first acquisition submodule is used to obtain the first atmospheric temperature T1 of the containment atmosphere before the hydrogen combustion after the accident, and the first nitrogen mole number N in the containment during normal operation. n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 ;

[0290] The third determination submodule is connected to the first acquisition submodule and is used to determine the first atmospheric temperature T1 and the first nitrogen mole number N n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 , determine the first nitrogen partial pressure P in the containment after the accident n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 ;

[0291] The fourth determination submodule is connected to the third determination submodule and is used to determine the nitrogen partial pressure P according to the first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 , determine the water vapor volume concentration f in the containment after the accident st and the hydrogen volume concentration fh2 ;

[0292] The fifth determination submodule is connected to the fourth determination submodule and is used to determine the water vapor volume concentration f st and the hydrogen volume concentration f h2 , determine the range of hydrogen concentration in the containment.

[0293] In some embodiments, the hydrogen concentration range includes one of hydrogen being in the flammable zone, hydrogen being at the boundary between the flammable zone and the non-flammable zone, and hydrogen being in the non-flammable zone;

[0294] The fifth determination submodule may be specifically used for:

[0295] In f st <1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )) In the case of hydrogen concentration in the containment, it is determined that the hydrogen is in the flammable zone;

[0296] Or, in f st =1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )) In the case of, the concentration range of hydrogen in the containment is determined as hydrogen in the containment at the junction of the combustible zone and the non-combustible zone;

[0297] Or, in f st >1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )), determine that the hydrogen concentration range in the containment is that hydrogen is in the non-flammable zone.

[0298] In some embodiments, the device may further include:

[0299] An updating module, used for updating the first atmospheric temperature T1 according to a preset step size;

[0300] The fifth determination module is connected to the update module and is used to repeatedly execute steps S112 to S115 on an iterative basis until it is determined that the hydrogen concentration range in the containment is at the junction of the combustible zone and the incombustible zone, and the iterative first atmospheric temperature T1 is obtained.

[0301] In some implementations, the first acquisition submodule may include:

[0302] The acquisition unit is used to obtain the second atmospheric total pressure P of the containment atmosphere during normal operation.a0 、The number of moles of zirconium in the fuel cladding of the reactor active area N zr , and the average temperature T0 of the containment atmosphere when the water vapor in the containment is saturated and operating normally;

[0303] The third determination unit is used to determine the water vapor partial pressure P in the containment during normal operation according to the average temperature T0. st0 ;

[0304] The fourth determining unit is used to determine the total atmospheric pressure P according to the second atmospheric total pressure P a0 and water vapor partial pressure P st0 , determine the second nitrogen partial pressure P in the containment during normal operation n20 and the second oxygen partial pressure P o20 ;

[0305] A fifth determining unit is configured to determine the nitrogen partial pressure P according to the second nitrogen partial pressure P n20 , the second oxygen partial pressure P o20 and the average temperature T0, determine the first nitrogen mole number N n20 and the first oxygen mole number N o20 ;

[0306] The sixth determining unit is used to determine the molar number N of zirconium according to zr , determine the first hydrogen mole number N h20 .

[0307] In some embodiments, the acquisition unit may be specifically configured to:

[0308] Obtain the total mass M of zirconium in the fuel cladding of the reactor active area zr and zirconium molar mass M zrm ;

[0309] The total mass of zirconium M zr and zirconium molar mass M zrm Substituting into formula (1) for calculation, we can obtain the molar number of zirconium in the fuel cladding of the reactor active area: N zr ;

[0310] Formula (1) includes:

[0311] N zr =M zr / M zrm (1).

[0312] In some implementations, the third determining unit may be specifically configured to:

[0313] The average temperature T0 is calculated using the international water and steam property calculation model IAPWS-IF97, and the water vapor partial pressure P in the containment during normal operation is obtained. st0 .

[0314] In some implementations, the fourth determining unit may be specifically configured to:

[0315] The second atmospheric total pressure P a0 and water vapor partial pressure P st0 Substitute into formula (2) to calculate and obtain the second nitrogen partial pressure P in the containment during normal operation: n20 ;

[0316] The second atmospheric total pressure P a0 and water vapor partial pressure P st0 Substitute into formula (3) to calculate and obtain the second oxygen partial pressure P in the containment during normal operation: o20 ;

[0317] Formula (2) includes:

[0318] P n20 =0.79×(P a0 -P st0 ) (2)

[0319] Formula (3) includes:

[0320] P o20 =0.21×(P a0 -P st0 ) (3).

[0321] In some implementations, the fifth determining unit may be specifically configured to:

[0322] Get the ideal gas constant R and the free volume V of the containment c ;

[0323] The ideal gas constant R, the free volume of the containment V c , the second nitrogen partial pressure P n20 Substitute the average temperature T0 into formula (4) to calculate and determine the first nitrogen mole number N n20 ;

[0324] According to the ideal gas constant R, the free volume V of the containment c , the second oxygen partial pressure P o20 Substitute the average temperature T0 into formula (5) to calculate and determine the first oxygen mole number N o20 ;

[0325] Formula (4) includes:

[0326] N n20 =P n20 V c / (RT0) (4)

[0327] Formula (5) includes:

[0328] N o20 =P o20 V c / (RT0) (5).

[0329] In some implementations, the sixth determining unit may be specifically configured to:

[0330] The number of moles of zirconium N zr Substitute into formula (6) to calculate and obtain the first hydrogen mole number N h20 ;

[0331] Formula (6) includes:

[0332] N h20 =2×N zr (6).

[0333] In some implementations, the third determining submodule may be specifically used to:

[0334] Get the ideal gas constant R and the free volume V of the containment c ;

[0335] The first nitrogen mole number N n20 , ideal gas constant R, first atmosphere temperature T1 and containment free volume V c Substitute into formula (7) to calculate and obtain the first nitrogen partial pressure P in the containment after the accident: n21 ;

[0336] The first oxygen mole number N o20 , ideal gas constant R, first atmosphere temperature T1 and containment free volume V c Substitute into formula (8) to calculate and obtain the first oxygen partial pressure P in the containment after the accident: o21 ;

[0337] The first hydrogen mole number N h20 , ideal gas constant R, first atmosphere temperature T1 and containment free volume V c Substitute into formula (9) to calculate and obtain the first hydrogen partial pressure P in the containment after the accident: h21 ;

[0338] The first atmospheric temperature T1 was calculated using the international water and water vapor property calculation model IAPWS-IF97, and the saturated water vapor partial pressure P in the containment after the accident was obtained. st1 ;

[0339] Among them, formula (7) includes:

[0340] P n21 =N n20 RT1 / V c (7)

[0341] Formula (8) includes:

[0342] P o21 =N o20 RT1 / V c (8)

[0343] Formula (9) includes:

[0344] P h21 =N h20 RT1 / V c (9).

[0345] In some implementations, the fourth determining submodule may be specifically configured to:

[0346] The first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 Substitute into formula (10) to calculate and obtain the third atmospheric total pressure P of the containment atmosphere after the accident: a1 ;

[0347] The saturated water vapor partial pressure P st1 and the third atmospheric total pressure P a1 Substitute into formula (11) to calculate and obtain the water vapor volume concentration f in the containment after the accident st ;

[0348] The first hydrogen partial pressure P h21 and the third atmospheric total pressure P a1 Substitute into formula (12) to calculate and obtain the hydrogen volume concentration f in the containment after the accident: h2 ;

[0349] Formula (10) includes:

[0350] P a1 =P n21 +P o21 +P st1 +P h21 (10)

[0351] Formula (11) includes:

[0352] f st =P st1 / P a1 (11)

[0353] Formula (12) includes:

[0354] f h2 =P h21 / Pa1 (12).

[0355] In some embodiments, the first quality information includes the first nitrogen quality M in the containment during normal operation. n20 、First oxygen mass M o20 、First hydrogen mass M h20 and the first water vapor mass M st1 ;

[0356] The second mass information includes the second nitrogen mass M in the containment after the hydrogen is completely burned n22 、The second oxygen mass M o22 and the second water vapor mass M st2 ;

[0357] The first constant volume specific heat capacity information includes the first constant volume specific heat capacity CV of nitrogen at the first atmospheric temperature T1 n21 、First oxygen constant volume specific heat capacity CV o21 、The first hydrogen constant volume specific heat capacity CV h21 and the first water vapor constant volume specific heat CV st1 ;

[0358] The second constant volume specific heat capacity information includes the first constant volume specific heat capacity CV of nitrogen at the second atmospheric temperature T2 n22 、Second oxygen constant volume specific heat capacity CV o22 、Second specific heat capacity of hydrogen at constant volume CV h22 and the second water vapor constant volume specific heat capacity CV st2 ;

[0359] The first determining unit may be specifically used for:

[0360] The first nitrogen mass M n20 、First oxygen mass M o20 、First hydrogen mass M h20 、First water vapor mass M st1 、The second nitrogen quality M n22 、The second oxygen mass M o22 、The second water vapor mass M st2 、First nitrogen constant volume specific heat capacity CV n21 、First oxygen constant volume specific heat capacity CV o21 、The first hydrogen constant volume specific heat capacity CV h21 、First water vapor constant volume specific heat capacity CV st1 、First nitrogen constant volume specific heat capacity CV n22 、Second oxygen constant volume specific heat capacity CV o22 、Second specific heat capacity of hydrogen at constant volume CV h22 、Second specific heat capacity of water vapor at constant volume CV st2 Total calories Q allSubstitute the first atmospheric temperature into formula (13) to calculate and obtain the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * ;

[0361] Formula (13) includes:

[0362] T2 * =(Q all +T1(M n20 CV n21 +M o20 CV o21 +M h20 CV h21 +M st1 CV st1 )) / (M n22 CV n22

[0363] +M o22 CV o22 +M st2 CV st2 )(13).

[0364] In some implementations, the first determining unit may be specifically configured to:

[0365] S1, in abs(T2-T2 * ) / T2≤N, according to the second nitrogen mole number N in the containment after the hydrogen is completely burned n22 、The second oxygen mole number N o22 、The second water vapor mole number N st2 and the second atmosphere temperature T2, determine the total pressure P of the first atmosphere in the containment after the hydrogen is completely burned adiabatically and at the same volume a2 ; or, in abs(T2-T2 * ) / T2>N, let T2=(T2+T2 * ) / 2, and repeatedly perform step S120 to step S1 on an iterative basis until the total pressure of the first atmosphere in the containment vessel after the adiabatic isochoric complete combustion of hydrogen is determined to be P a2 ;

[0366] Wherein, N is a positive number.

[0367] In some implementations, the first determining unit may be specifically configured to:

[0368] Get the ideal gas constant R and the free volume V of the containment c ;

[0369] The second nitrogen mole number N n22 、The second oxygen mole number N o22 、The second water vapor mole number N st2, the second atmosphere temperature T2, the ideal gas constant R and the free volume V of the containment c Substituting into formula (14) for calculation, we can obtain the total pressure of the first atmosphere in the containment after adiabatic isochoric complete combustion of hydrogen: a2 ;

[0370] Formula (14) includes:

[0371] P a2 =(N n22 +N o22 +N st2 )RT2 / V c (14).

[0372] In some embodiments, the first quality information includes the first nitrogen quality M in the containment during normal operation. n20 、First oxygen mass M o20 、First hydrogen mass M h20 and the first water vapor mass M st1 ;

[0373] The acquisition unit can be used specifically for:

[0374] Get the molar mass of nitrogen M n2m 、Oxygen molar mass M o2m 、Hydrogen molar mass M h2m , water vapor molar mass M stm 、The first nitrogen mole number N n20 、The first oxygen mole number N o20 、The first hydrogen mole number N h20 、Free volume of containment V c , and the density of saturated water vapor in the containment ρ st1 ;

[0375] The molar mass of nitrogen M n2m and the first nitrogen mole number N n20 Substitute into formula (15) to calculate and obtain the first nitrogen mass M in the containment during normal operation: n20 ;

[0376] The molar mass of oxygen M o2m and the first oxygen mole number N o20 Substitute into formula (16) to calculate and obtain the first oxygen mass M in the containment during normal operation: o20 ;

[0377] The molar mass of hydrogen M h2m and the first hydrogen mole number N h20 Substitute into formula (17) to calculate and obtain the first hydrogen mass M in the containment during normal operation: h20 ;

[0378] The free volume of the containment V c and density ρ st1 Substitute into formula (18) to calculate and obtain the first steam mass M in the containment during normal operation: st1 ;

[0379] Formula (15) includes:

[0380] M n20 =M n2m N n20 (15)

[0381] Formula (16) includes:

[0382] M o20 =M o2m N o20 (16)

[0383] Formula (17) includes:

[0384] M h20 =M h2m N h20 (17)

[0385] Formula (18) includes:

[0386] M st1 =ρ st1 V c (18).

[0387] In some embodiments, the second mass information includes the second nitrogen mass M in the containment after the hydrogen is completely burned. n22 、The second oxygen mass M o22 and the second water vapor mass M st2 ;

[0388] The acquisition unit can be used specifically for:

[0389] Get the first nitrogen mole number N n20 、The first oxygen mole number N o20 , the first hydrogen mole number N h20 , Nitrogen molar mass M n2m 、Oxygen molar mass M o2m , water vapor molar mass M stm , and the first steam mole number N in the containment during normal operation st1 ;

[0390] The first nitrogen mole number N n20 Substitute into formula (19) to calculate and obtain the second nitrogen mole number N in the containment after the hydrogen is completely burned: n22 ;

[0391] The first oxygen mole number N o20 and the first hydrogen mole number N h20 Substitute into formula (20) to calculate and obtain the second oxygen mole number N in the containment after the hydrogen is completely burned: o22 ;

[0392] The first water vapor mole number N st1 and the first hydrogen mole number N h20 Substitute into formula (21) to calculate and obtain the second water vapor mole number N in the containment after the hydrogen is completely burned: st2 ;

[0393] The second nitrogen mole number N n22 Substitute the molar mass of nitrogen into formula (22) to calculate the second nitrogen mass M in the containment after the hydrogen is completely burned. n22 ;

[0394] The second oxygen mole number N o22 and oxygen molar mass M o2m Substitute into formula (23) to calculate and obtain the second oxygen mass M in the containment after the hydrogen is completely burned: o22 ;

[0395] The second water vapor mole number N st2 and the water vapor molar mass M stm Substituting into formula (24) for calculation, we can obtain the second steam mass M in the containment after the hydrogen is completely burned: st2 ;

[0396] Formula (19) includes:

[0397] N n22 =N n20 (19)

[0398] Formula (20) includes:

[0399] N o22 =N o20 -N h20 / 2 (20)

[0400] Formula (21) includes:

[0401] N st2 =N st1 +N h20 (twenty one)

[0402] Formula (22) includes:

[0403] M n22 =M n2m N n22 (twenty two)

[0404] Formula (23) includes:

[0405] M o22 =M o2m N o22 (twenty three)

[0406] Formula (24) includes:

[0407] M st2 =M stm N st2 (twenty four).

[0408] In some embodiments, the acquisition unit may be specifically configured to:

[0409] Get the first water vapor mass M st1 and the water vapor molar mass M stm ;

[0410] The first water vapor mass M st1 and the water vapor molar mass M stm Substitute into formula (25) to calculate and obtain the first water vapor mole number N in the containment during normal operation: st1 ;

[0411] Formula (25) includes:

[0412] N st1 =M st1 / M stm (25).

[0413] In some embodiments, the acquisition unit may be specifically configured to:

[0414] Obtain the combustion heat Q of hydrogen and the first number of hydrogen moles N h20 ;

[0415] The heat of combustion Q and the first number of hydrogen moles N h20 Substitute into formula (26) to calculate and obtain the first hydrogen mole number N h20 The total amount of heat released by complete combustion;

[0416] Formula (26) includes:

[0417] Q all =N h20 Q (26).

[0418] The device for determining the total atmospheric pressure in a containment shell provided in the embodiment of the present application can be used to execute the method in Example 1, that is, it has the beneficial effects and implementation methods of the method for determining the total atmospheric pressure in a containment shell provided in Example 1 of the present application. For details, please refer to the specific description of the method for determining the total atmospheric pressure in a containment shell in the above-mentioned Example 1, which will not be repeated in this embodiment.

[0419] It is to be understood that the above embodiments are merely exemplary embodiments used to illustrate the principles of the present application, but the present application is not limited thereto. For those skilled in the art, various modifications and improvements can be made without departing from the spirit and substance of the present application, and these modifications and improvements are also considered to be within the scope of protection of the present application.

Claims

1. A method for determining the total atmospheric pressure in a containment vessel, characterized in that: include: S110, obtaining the hydrogen concentration range in the containment; S120, determining the first atmospheric total pressure P in the containment according to the hydrogen concentration range a2 .

2. The method according to claim 1, characterized in that Step S120 includes: S121, when the hydrogen concentration range in the containment is such that the hydrogen is at the boundary of the combustible zone and the incombustible zone, determining a first atmospheric temperature T1 corresponding to the hydrogen concentration range, and a second atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned; S122: Determine the first atmospheric total pressure P in the containment according to the first atmospheric temperature T1 and the second atmospheric temperature T2. a2 .

3. The method according to claim 2, characterized in that Step S122 includes: S1221. Obtain first mass information, second mass information, first specific heat capacity at constant volume information, second specific heat capacity at constant volume information, and first hydrogen mole number N h20 The total heat released by complete combustion is Q all ; S1222, according to the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, the total heat Q all and the first atmospheric temperature T1, determine the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * ; S1223, according to the second atmospheric temperature T2 and the third atmospheric temperature T2 * , determine the total pressure of the first atmosphere in the containment after the adiabatic isochoric complete combustion of hydrogen P a2 .

4. The method according to claim 1, characterized in that: Step S110 includes: S111. Obtain the first atmospheric temperature T1 of the containment atmosphere before hydrogen combustion after the accident, and the first nitrogen mole number N in the containment during normal operation n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 ; S112, according to the first atmospheric temperature T1, the first nitrogen molar number N n20 , the first oxygen mole number N o20 and the first hydrogen mole number N h20 , determine the first nitrogen partial pressure P in the containment after the accident n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 ; S113, according to the first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and the saturated water vapor partial pressure P st1 , determine the water vapor volume concentration f in the containment after the accident st and the hydrogen volume concentration f h2 ; S114, according to the water vapor volume concentration f st and the hydrogen volume concentration f h2 , determine the hydrogen concentration range in the containment.

5. The method according to claim 4, characterized in that The hydrogen concentration range includes one of hydrogen being in the combustible zone, hydrogen being at the junction of the combustible zone and the non-combustible zone, and hydrogen being in the non-combustible zone; Step S114 includes: In f st <1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )), determining that the hydrogen concentration range in the containment is such that hydrogen is in the flammable zone; Or, in f st =1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )), determining that the hydrogen concentration range in the containment is that the hydrogen is at the junction of the combustible zone and the non-combustible zone in the containment; Or, in f st >1-(f h2 +0.373exp(-0.7f h2 )+5.18exp(-48.8f h2 )), it is determined that the hydrogen concentration range in the containment is that hydrogen is in the non-flammable zone.

6. The method according to claim 5, characterized in that After determining that the hydrogen concentration range in the containment is that the hydrogen is in a flammable zone, the method further includes: S115, updating the first atmospheric temperature T1 according to a preset step size; S116. Repeat steps S112 to S115 on an iterative basis until it is determined that the hydrogen concentration range in the containment is at the junction of the combustible zone and the incombustible zone, and the iterative first atmospheric temperature T1 is obtained.

7. The method according to claim 4, characterized in that Get the first nitrogen mole number N in the containment during normal operation n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 ,include: Obtain the second atmospheric total pressure P of the containment atmosphere during normal operation a0 、The number of moles of zirconium in the fuel cladding of the reactor active area N zr , and the average temperature T0 of the atmosphere in the containment when the water vapor in the containment is in a saturated state and operating normally; According to the average temperature T0, the water vapor partial pressure P in the containment during normal operation is determined. st0 ; According to the second atmospheric total pressure P a0 and the water vapor partial pressure P st0 , determine the second nitrogen partial pressure P in the containment during normal operation n20 and the second oxygen partial pressure P o20 ; According to the second nitrogen partial pressure P n20 , the second oxygen partial pressure P o20 and the average temperature T0, determine the first nitrogen mole number N n20 and the first oxygen mole number N o20 ; According to the molar number N of zirconium zr , determine the first hydrogen mole number N h20 .

8. The method according to claim 7, characterized in that Get the number of zirconium moles N in the fuel cladding of the reactor active area zr ,include: Obtain the total mass M of zirconium in the fuel cladding of the reactor active area zr and zirconium molar mass M zrm ; The total mass of zirconium M zr and zirconium molar mass M zrm Substituting into formula (1) for calculation, we can obtain the molar number of zirconium in the fuel cladding of the reactor active area: N zr ; The formula (1) includes: N zr =M zr / M zrm (1)。 9. The method according to claim 7, characterized in that: The water vapor partial pressure P in the containment during normal operation is determined according to the average temperature T0. st0 ,include: The average temperature T0 is calculated using the international water and steam property calculation model IAPWS-IF97 to obtain the water vapor partial pressure P in the containment during normal operation. st0 .

10. The method according to claim 7, characterized in that The second atmospheric total pressure P a0 and the water vapor partial pressure P st0 , determine the second nitrogen partial pressure P in the containment during normal operation n20 and the second oxygen partial pressure P o20 ,include: The second atmospheric total pressure P a0 and the water vapor partial pressure P st0 Substitute into formula (2) to calculate and obtain the second nitrogen partial pressure P in the containment during normal operation: n20 ; The second atmospheric total pressure P a0 and the water vapor partial pressure P st0 Substitute into formula (3) to calculate and obtain the second oxygen partial pressure P in the containment during normal operation: o20 ; The formula (2) includes: P n20 =0.79×(P a0 -P st0 ) (2) The formula (3) includes: P o20 =0.21×(P a0 -P st0 ) (3)。 11. The method according to claim 7, characterized in that According to the second nitrogen partial pressure P n20 , the second oxygen partial pressure P o20 and the average temperature T0, determine the first nitrogen mole number N n20 and the first oxygen mole number N o20 ,include: Get the ideal gas constant R and the free volume V of the containment c ; The ideal gas constant R, the free volume V of the containment c , the second nitrogen partial pressure P n20 Substitute the average temperature T0 into formula (4) to calculate and obtain the first nitrogen mole number N n20 ; According to the ideal gas constant R, the free volume V of the containment c , the second oxygen partial pressure P o20 Substitute the average temperature T0 into formula (5) to calculate the first oxygen mole number N o20 ; The formula (4) includes: N n20 =P n20 V c / (RT0) (4) The formula (5) includes: N o20 =P o20 V c / (RT0) (5)。 12. The method according to claim 7, characterized in that According to the molar number N of zirconium zr , determine the first hydrogen mole number N h20 ,include: The molar number N of zirconium zr Substitute into formula (6) to calculate and obtain the first hydrogen mole number N h20 ; The formula (6) includes: N h20 =2×N zr (6)。 13. The method according to claim 4, characterized in that Step S112 includes: Get the ideal gas constant R and the free volume V of the containment c ; The first nitrogen mole number N n20 , the ideal gas constant R, the first atmospheric temperature T1 and the free volume V of the containment vessel c Substitute into formula (7) to calculate and obtain the first nitrogen partial pressure P in the containment after the accident: n21 ; The first oxygen mole number N o20 , the ideal gas constant R, the first atmospheric temperature T1 and the free volume V of the containment vessel c Substitute into formula (8) to calculate and obtain the first oxygen partial pressure P in the containment after the accident: o21 ; The first hydrogen mole number N h20 , the ideal gas constant R, the first atmospheric temperature T1 and the free volume V of the containment vessel c Substitute into formula (9) to calculate and obtain the first hydrogen partial pressure P in the containment after the accident: h21 ; The first atmospheric temperature T1 is calculated using the international water and water vapor property calculation model IAPWS-IF97 to obtain the saturated water vapor partial pressure P in the containment after the accident. st1 ; Wherein, the formula (7) includes: P n21 =N n20 RT1 / V c (7) The formula (8) includes: P o21 =N o20 RT1 / V c (8) The formula (9) includes: P h21 =N h20 RT1 / V c (9)。 14. The method according to claim 4, characterized in that Step S113 includes: The first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and the saturated water vapor partial pressure P st1 Substitute into formula (10) to calculate and obtain the third atmospheric total pressure P of the containment atmosphere after the accident: a1 ; The saturated water vapor partial pressure P st1 and the third atmospheric total pressure P a1 Substitute into formula (11) to calculate and obtain the water vapor volume concentration f in the containment after the accident: st ; The first hydrogen partial pressure P h21 and the third atmospheric total pressure P a1 Substitute into formula (12) to calculate and obtain the hydrogen volume concentration f in the containment after the accident: h2 ; The formula (10) includes: P a1 =P n21 +P o21 +P st1 +P h21 (10) The formula (11) includes: f st =P st1 / P a1 (11) The formula (12) includes: f h2 =P h21 / P a1 (12)。 15. The method according to claim 3, characterized in that The first quality information includes the first nitrogen quality M in the containment during normal operation n20 、First oxygen mass M o20 , the first hydrogen mass M h20 and the first water vapor mass M st1 ; The second mass information includes the second nitrogen mass M in the containment after the hydrogen is completely burned. n22 、The second oxygen mass M o22 and the second water vapor mass M st2 ; The first constant volume specific heat capacity information includes the first constant volume specific heat capacity CV of nitrogen at the first atmospheric temperature T1 n21 、First oxygen constant volume specific heat capacity CV o21 、The first hydrogen constant volume specific heat capacity CV h21 and the first water vapor constant volume specific heat CV st1 ; The second constant volume specific heat capacity information includes the first constant volume specific heat capacity CV of nitrogen at the second atmospheric temperature T2 n22 、Second oxygen constant volume specific heat capacity CV o22 、Second specific heat capacity of hydrogen at constant volume CV h22 and the second water vapor constant volume specific heat capacity CV st2 ; Step S1222 includes: The first nitrogen mass M n20 、The first oxygen mass M o20 , the first hydrogen mass M h20 , the first water vapor mass M st1 , the second nitrogen mass M n22 , the second oxygen mass M o22 , the second water vapor mass M st2 、The first nitrogen constant volume specific heat capacity CV n21 、The first oxygen constant volume specific heat capacity CV o21 、The first hydrogen constant volume specific heat capacity CV h21 、The first water vapor constant volume specific heat capacity CV st1 、The first nitrogen constant volume specific heat capacity CV n22 The second oxygen constant volume specific heat capacity CV o22 、The second hydrogen constant volume specific heat capacity CV h22 、The second water vapor constant volume specific heat capacity CV st2 The total heat Q all Substitute the first atmospheric temperature into formula (13) to calculate and obtain the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * ; The formula (13) includes: T2 * =(Q all +T1(M n20 CV n21 +M o20 CV o21 +M h20 CV h21 +M st1 CV st1 )) / (M n22 CV n22 +M o22 CV o22 +M st2 CV st2 )(13)。 16. The method according to claim 15, characterized in that The second atmospheric temperature T2 and the third atmospheric temperature T2 * , determine the total pressure of the first atmosphere in the containment vessel after the adiabatic isochoric complete combustion of hydrogen P a2 ,include: S1, in abs(T2-T2 * ) / T2≤N, according to the second nitrogen mole number N in the containment after the hydrogen is completely burned n22 、The second oxygen mole number N o22 、The second water vapor mole number N st2 and the second atmospheric temperature T2, determine the first atmospheric total pressure P in the containment after adiabatic isochoric complete combustion of hydrogen a2 ; or, in abs(T2-T2 * ) / T2>N, let T2=(T2+T2 * ) / 2, and repeatedly perform step S120 to step S1 on an iterative basis until it is determined that the first atmospheric total pressure P in the containment vessel after adiabatic isochoric complete combustion of hydrogen a2 ; Wherein, N is a positive number.

17. The method according to claim 16, characterized in that The second nitrogen mole number N in the containment after the hydrogen is completely burned n22 、The second oxygen mole number N o22 、The second water vapor mole number N st2 and the second atmospheric temperature T2, determine the first atmospheric total pressure P in the containment after adiabatic isochoric complete combustion of hydrogen a2 ,include: Get the ideal gas constant R and the free volume V of the containment c ; The second nitrogen mole number N n22 、The second oxygen mole number N o22 、The second water vapor mole number N st2 , the second atmospheric temperature T2, the ideal gas constant R and the containment free volume V c Substituting into formula (14) for calculation, the total pressure of the first atmosphere in the containment after adiabatic isochoric complete combustion of hydrogen is obtained as: a2 ; The formula (14) includes: P a2 =(N n22 +N o22 +N st2 )RT2 / V c (14).

18. The method according to claim 3, characterized in that The first quality information includes the first nitrogen quality M in the containment during normal operation n20 、First oxygen mass M o20 , the first hydrogen mass M h20 and the first water vapor mass M st1 ; Get first quality information, including: Get the molar mass of nitrogen M n2m 、Oxygen molar mass M o2m , hydrogen molar mass M h2m , water vapor molar mass M stm 、The first nitrogen mole number N n20 、The first oxygen mole number N o20 、The first hydrogen mole number N h20 、Free volume of containment V c , and the density of saturated water vapor in the containment vessel ρ st1 ; The nitrogen molar mass M n2m and the first nitrogen mole number N n20 Substitute into formula (15) to calculate and obtain the first nitrogen mass M in the containment during normal operation: n20 ; The oxygen molar mass M o2m and the first oxygen mole number N o20 Substitute into formula (16) to calculate and obtain the first oxygen mass M in the containment during normal operation: o20 ; The hydrogen molar mass M h2m and the first hydrogen mole number N h20 Substitute into formula (17) to calculate and obtain the first hydrogen mass M in the containment during normal operation: h20 ; The free volume V of the containment c and the density ρ st1 Substitute into formula (18) to calculate and obtain the first steam mass M in the containment during normal operation: st1 ; The formula (15) includes: M n20 =M n2m N n20 (15) The formula (16) includes: M o20 =M o2m N o20 (16) The formula (17) includes: M h20 =M h2m N h20 (17) The formula (18) includes: M st1 =ρ st1 V c (18)。 19. The method according to claim 3, characterized in that The second mass information includes the second nitrogen mass M in the containment after the hydrogen is completely burned. n22 、The second oxygen mass M o22 and the second water vapor mass M st2 ; Obtain second quality information, including: Get the first nitrogen mole number N n20 、The first oxygen mole number N o20 , the first hydrogen mole number N h20 , Nitrogen molar mass M n2m 、Oxygen molar mass M o2m , water vapor molar mass M stm , and the first water vapor mole number N in the containment during normal operation st1 ; The first nitrogen mole number N n20 Substitute into formula (19) to calculate and obtain the second nitrogen mole number N in the containment after the hydrogen is completely burned: n22 ; The first oxygen mole number N o20 and the first hydrogen mole number N h20 Substitute into formula (20) to calculate and obtain the second oxygen mole number N in the containment after the hydrogen is completely burned: o22 ; The first water vapor mole number N st1 and the first hydrogen mole number N h20 Substitute into formula (21) to calculate and obtain the second water vapor mole number N in the containment after the hydrogen is completely burned: st2 ; The second nitrogen mole number N n22 Substitute the molar mass of nitrogen into formula (22) to calculate and obtain the second nitrogen mass M in the containment after the hydrogen is completely burned. n22 ; The second oxygen mole number N o22 and the oxygen molar mass M o2m Substitute into formula (23) to calculate and obtain the second oxygen mass M in the containment after the hydrogen is completely burned: o22 ; The second water vapor mole number N st2 and the water vapor molar mass M stm Substitute into formula (24) to calculate and obtain the second water vapor mass M in the containment after the hydrogen is completely burned: st2 ; The formula (19) includes: N n22 =N n20 (19) The formula (20) includes: N o22 =N o20 -N h20 / 2 (20) The formula (21) includes: N st2 =N st1 +N h20 (21) The formula (22) includes: M n22 =M n2m N n22 (22) The formula (23) includes: M o22 =M o2m N o22 (23) The formula (24) includes: M st2 =M stm N st2 (24)。 20. The method according to claim 19, characterized in that Get the first water vapor mole number N in the containment during normal operation st1 ,include: Get the first water vapor mass M st1 and the water vapor molar mass M stm ; The first water vapor mass M st1 and the water vapor molar mass M stm Substitute into formula (25) to calculate and obtain the first water vapor mole number N in the containment during normal operation: st1 ; Formula (25) includes: N st1 =M st1 / M stm (25)。 21. The method according to claim 3, characterized in that Get the first hydrogen mole number N h20 The total heat released by complete combustion is Q all ,include: Obtain the combustion heat Q of hydrogen and the first number of hydrogen moles N h20 ; The combustion heat Q and the first hydrogen mole number N h20 Substitute into formula (26) to calculate and obtain the first hydrogen mole number N h20 The total amount of heat released by complete combustion; The formula (26) includes: Q all =N h20 Q (26)。 22. A device for determining the total atmospheric pressure in a containment vessel, characterized in that: include: The first acquisition module is used to obtain the hydrogen concentration range in the containment; A first determination module, connected to the first acquisition module, is used to determine the first atmospheric total pressure P in the containment according to the hydrogen concentration range. a2 .

23. The device according to claim 22, characterized in that The first determining module includes: A first determination submodule is used to determine, when the hydrogen concentration range in the containment is at the junction of the combustible zone and the incombustible zone, a first atmospheric temperature T1 corresponding to the hydrogen concentration range, and a second atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned; The second determination submodule is connected to the first determination submodule and is used to determine the first atmospheric total pressure P in the containment according to the first atmospheric temperature T1 and the second atmospheric temperature T2. a2 .

24. The device according to claim 23, characterized in that The second determining submodule includes: An acquisition unit, used to acquire the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, and the first hydrogen mole number N h20 The total heat released by complete combustion is Q all ; A first determining unit is connected to the acquiring unit and is used for determining the first mass information, the second mass information, the first constant volume specific heat capacity information, the second constant volume specific heat capacity information, the total heat Q all and the first atmospheric temperature T1, determine the third atmospheric temperature T2 of the containment atmosphere after the hydrogen is completely burned * ; A second determining unit is connected to the first determining unit and is used to determine the temperature of the second atmosphere T2 and the third atmosphere T2 according to the second atmosphere temperature T2 * , determine the total pressure of the first atmosphere in the containment after the adiabatic isochoric complete combustion of hydrogen P a2 .

25. The device according to claim 22, characterized in that The first acquisition module includes: The first acquisition submodule is used to obtain the first atmospheric temperature T1 of the containment atmosphere before the hydrogen combustion after the accident, and the first nitrogen mole number N in the containment during normal operation. n20 、The first oxygen mole number N o20 and the first hydrogen mole number N h20 ; The third determination submodule is connected to the first acquisition submodule and is used to determine the first atmospheric temperature T1, the first nitrogen mole number N n20 , the first oxygen mole number N o20 and the first hydrogen mole number N h20 , determine the first nitrogen partial pressure P in the containment after the accident n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and saturated water vapor partial pressure P st1 ; A fourth determination submodule is connected to the third determination submodule and is used to determine the nitrogen partial pressure P according to the first nitrogen partial pressure P n21 , the first oxygen partial pressure P o21 , the first hydrogen partial pressure P h21 and the saturated water vapor partial pressure P st1 , determine the water vapor volume concentration f in the containment after the accident st and the hydrogen volume concentration f h2 ; The fifth determination submodule is connected to the fourth determination submodule and is used to determine the water vapor volume concentration f st and the hydrogen volume concentration f h2 , determine the hydrogen concentration range in the containment.

Images