Study on Leakage Model of Typical Penetration of Closed Structures Based on Porous Media Seepage Theory

. Quantitative prediction of penetration leakage rate has not been reported so far. The theoretical prediction method of leakage rate of typical penetration static seal structure of closed structure is studied. Based on porous media seepage mechanics, the relationship between leakage rate and microstructure parameters is established, and a new prediction model of the leakage rate of sealing structures is proposed. Then, the relationship between seal-speciﬁc pressure and microstructure parameters is obtained by the Hertz contact mechanics model. Finally, the leakage rate calculation, which is independent of any experimental data, is obtained innovatively. The new model is used to predict the leakage rate of penetration, and the theoretical prediction results are compared with the experimental measurement results. It is proved that the two agree well, which veriﬁes the eﬀectiveness of the model. The prediction model based on this method can well reﬂect the eﬀects of rough surface morphology, material mechanical properties, sealing load, high temperature, and high pressure on the leakage rate. The leakage rate prediction model proposed in this paper is independent of the experimental regression coeﬃcient and can realize the conversion relationship between diﬀerent sealing media, sealing materials, and working conditions.


Introduction
Nuclear power plant containment is a typical closed structure. e pressure boundary formed by the containment and much perforated equipment, components, and pipelines penetrating the containment is not only the last barrier between the internal equipment or system and the environment but also a very important barrier related to safety. Penetration is usually the main way of containment leakage in nuclear power plants. Penetration refers to the gaping hole equipment, components, pipelines, and isolation valves penetrating the closed structure. ere are a large number of penetrations. Ensuring their sealing integrity under various working conditions is an important content in the design of closed structures [1][2][3]. Penetration can be divided into electrical penetration and mechanical penetration. e principal function of the electrical penetration assembly is to connect the electrical equipment to the outside of the containment to control the equipment inside the containment. ere are two types of mechanical penetration: one is fuel transportation penetration; the other is mechanical pipe penetration.
In order to find an effective method to improve the reliability of penetration seal, researchers at home and abroad have carried out a lot of research. Argonne National Laboratory, Sandia National Laboratory, and Idaho nuclear engineering laboratory have conducted extensive research on the reliability of penetration seal structure of 22 nuclear power plants [4,5], put forward some common problems in penetration seal, and discussed the influence of high temperature and high pressure on gasket sealing performance from the perspective of flange warpage and deformation. Based on experimental research, some researchers qualitatively discussed the damage of mechanical properties of sealing materials caused by high temperature and irradiation [4][5][6][7]. In addition, many scholars have used three-dimensional finite element technology to study the mechanical deformation of penetration system under serious accidents [8,9]. In addition, many scholars have established a series of models and theories using porous media seepage theory. However, there are few reports on the quantitative prediction of the leakage rate of penetration seal structure.
Quantitative analysis of the leakage rate is the basis of the safety design of penetration seal structure. However, the quantitative prediction research on the leakage rate of penetration seal structure has not been reported so far. e main reason is that people's understanding of the leakage mechanism of seal structure is not clear enough. Only by deeply studying the microleakage mechanism of sealing structure can the quantitative relationship between leakage rate and various influencing factors be established, and the research on penetration sealing characteristics can be developed from "qualitative" to "quantitative." Existing leakage rate calculation models are semiempirical without exception; that is, there are regression coefficients in the model that has unclear physical meaning and need to be fitted with experimental results. e dependence on experimental research makes these calculation models unable to be applied to the prediction of gasket leakage rate of nonstandard structure, special working conditions, and special gas without experimental data. e key problem is that there is no appropriate leakage mechanism model of penetration seal structure. e purpose of this paper is to introduce the porous media seepage theory and build the sealing microleakage mechanism model based on the essence of gas leakage phenomenon. e focus is to establish the quantitative relationship between microfine structure and macromeasure and then put forward a prediction and calculation method of leakage rate independent of experimental regression parameters. On this basis, the accuracy of the prediction model is verified by experiments. is method can effectively analyze the coupling effect of interface microstructure and material macromechanical properties on the leakage rate. It provides a theoretical basis for the development and design of penetration seal structure and provides guarantee for environmental quality and safety.

Leakage Mechanism Model of the Static
Seal Interface e interface leakage of static seal generally refers to the leakage of fluid through the rough contact surface between the flange and sealing material. e workpiece surface is always rough in varying degrees. e contact surface between the flange and the sealing material is not particularly close, and the sealing medium can leak out from these gaps. e close contact between the two sealing surfaces is realized by applying a compressive load, to reduce the pore space in the sealing gap and increase the flow resistance. Under the action of internal and external pressure differences, the medium will leak from a high-pressure environment to a low-pressure environment. According to different flow paths, leakage can be structured in interface leakage (seepage through different contact material surfaces) and seepage leakage (seepage from the inside of sealing material). e density of materials used for sealing is generally bulky, and the proportion of seepage leakage is very small. erefore, this paper is only focused on the study of interface leakage. e gas flow characteristics in a rough gap are studied using the mesoscopic lattice Boltzmann method (LBM). First, the virtual interface leakage channel is established based on the three-dimensional numerical reconstruction technology of rough surface, and then, the influence of gas rough interface on gas flow is analyzed by LBM simulation. Figures 1 and 2 show the rough numerical surface constructed based on Gauss distribution and the LBM model used.
e simulation results show that when the channel height remains unchanged, the fluid flow in the interface still retains the characteristics of plane Poiseuille flow. erefore, a rough flow factor can be defined Φ σ : In the formula, Q 0 represents the volume flow rate of the rough surface channel, and Q p represents the volume flow rate of the smooth plate channel with the same gap height. According to the model of parallel plate flow, In the formula, h represents the clearance height of the smooth flat plate channel, B is the flow channel width, l is the flow channel length, p 2 is the internal pressure, and p 1 is the external ambient pressure. Φ σ is only related to the characteristics of rough surface, which can be defined as dimensionless roughness σ * function of (σ * � σ/T), where s and T is two statistical parameters of a rough surface, namely, roughness(μm) and autocorrelation scale(μm). According to the numerical results Φσ − σ * , the curve is shown in Figure 3 (1), and its function expression is Based on the numerical results, the height flow factor is further defined Φ h to characterize the effect of gap height change on gap flow at the rough interface, and its expression is In the formula, Q 0 is the initial clearance height h 0 obtained in (2). Q represents the volume flow corresponding to the actual height h. When the characteristics of a rough surface are certain, Φ h is only the dimensionless gap height h * (h * � h/σ). Figure 3 (2) shows the curve fitted according to the numerical calculation results, and its functional relationship is In this way, when the characteristic parameters of a rough surface and the change of interface gap height are   Mathematical Problems in Engineering known, the gas volume flow through the interface microporous structure can be calculated using the following formula: Formula (6) effectively reveals the flow characteristics of gas in the sealing interface, and it offers an effective way to build the interface leakage mechanism model based on the numerical calculation method.

Construction Method of a Theoretical Model of the Penetration Leakage Rate of Typical Structures
Taking the manhole pressure cover as an example, the penetration leakage theoretical model of typical structures is constructed [4][5][6][7].

Finite Element Analysis
Model. e pressure cover of the manhole is one of the main penetrations installed on the closure. Structurally, it can classify as a large-diameter nonstandard bolt flange gasket structure. e sealing principle is to use the pretightening force of the bolt to make the cover compress the gasket and from contact pressure on the sealing surface through the elastic deformation of the gasket, to achieve the sealing effect. Figure 4 is the assembly diagram of the pressure cover of the manhole and the established three-dimensional schematic diagram [8].
Because of the characteristics of nonstandard and large diameter, some theoretical analysis methods only apply to standard bolt flange gasket and cannot be directly applied to the analyses of mechanical deformation characteristics of pressure cover, and the finite element analysis method is still necessary. To control the leakage rate of the static seal structure of the manhole, it is necessary to clarify the main influencing factors and the variation law of the leakage rate for these factors. erefore, it is the key to constructing the theoretical prediction model of the leakage rate of the pressure cover of the manhole [8,9]. e method of constructing the theoretical model of pressure covers leakage rate is divided into three steps. First, according to the microcosmic mechanical analysis of the rough element of sealing contact surface, the microcosmic leakage mechanism model is constructed. Second, the influence of various factors on the contact seal of sealing surface is analyzed according to the macrofinite element method.
ird, using the coupling relationship between sealing surface contact seal and the microcosmic leakage mechanism, the leakage rate prediction model of pressure cover sealing structure is established. is method realizes the effective coupling of microdetailed structural analysis and macromechanical characteristic analysis and does not include any experimental regression coefficient. It can effectively analyze the change law and influencing factors of pressure cover leakage rate, and it is not difficult to control its leakage rate [5,[10][11][12]. e cover plate is in the form of flat plate reinforcement, and it is not rotational axisymmetry. erefore, the 1/4 model is utilized to simplify the analyses. e physical parameters used in the calculation are as follows. e elastic modulus of the cover plate, lower flange, and cylinder is 210000 MPa and Poisson's ratio is 0.3. e elastic modulus of bolts and nuts is set to 203000 MPa, and Poisson's ratio is 0.3; e gasket is made of neoprene with a thickness of 5 mm, and its stress-strain curve is obtained by consulting the literature (refer to Figure 5) [13,14].
e flange joint has axisymmetric properties in structure, load, and constraint. erefore, boundary symmetry constraints are imposed on the two symmetry planes of the whole finite element model to constrain the circumferential displacement of the two plane nodes. At the same time, to limit the overall rigid body displacement of the joint, fixed boundary conditions are imposed on the truncated surface of the cylinder.
Calculate the bolt design load of flange joint during assembly and operation, with gasket pretension ratio pressure Y and gasket factors M.
(1) Determine the minimum bolt load required under operating conditions W m1 where D G is the iameter at the center circle of the gasket compression stress, mm, b is the effective sealing width of gasket, mm, p is the pressure of medium in the vessel, MPa, and m is the gasket coefficient. (2) Determine the minimum bolt load required for the pretension gasket W m2 where y is the gasket pretension ratio pressure, MPa. e design load of the bolt can be obtained by bringing in the gasket parameters: ODB: Job-719-F=70KN.odb Abaqus/Standard 6.14-4 Wed Nov 18 02:27:13 GMT+08:00 2020 Step: Step-pressure1 Increment 1:  Figure 6: Mathematical Problems in Engineering 5 e design load of a single bolt can be obtained by bringing in the gasket parameters: e calculation process is roughly as follows. To make the analysis easier to converge, a tiny bolt force is applied in the first analysis step to establish the contact relationship between various surfaces smoothly. In the second analysis step, the bolt force is brought to the pretension. In the third analysis step, the pressure under various working conditions is applied to the flange and the inner wall of the cover [15,16].

Finite Element Analysis Results.
e gasket material is neoprene. Considering the stress distribution and the leakage rate of the gasket when the bolt pretension f is 70 kN, 80 kN, 90 kN, 100 kN, and 110 kN, the medium working pressure is, respectively, Δ P � 0.45 MPa, 1 MPa, 1.5 MPa, and 2 MPa. When the pretension is 70 kN and considering the change of working pressure, the overall stress distribution and gasket stress distribution are shown in Figure 6 respectively. When the medium working pressure is 0.45 MPa and considering the change of pretension, the overall stress distribution and gasket stress distribution are shown in Figures 8 and 9, respectively [17].
It can be perceived from Figures 6 and 7 that, under the same pretension, the stress distribution of the overall structure of the bolt flange is similar, the stress is concentrated at the connection between the bolt and nut and the connection between the stiffener and the cover, and the maximum stress value increases significantly with the increase of pressure, while the gasket stress value decreases with the increase of pressure. According to Figures 8 and 9, under the same medium pressure, with the increase of pretension, the overall maximum stress value increases, and the gasket stress value also increases. Considering that largediameter flange is prone to warpage under large internal pressure, the following figure and table list the stress value of the specific structure and the change of flange angle under pressure change and pretension change, respectively.
It can be seen from Tables 1 and 2 that, under the given pretension and pressure, the flange angle is less than 0.3°s pecified by ASME. With the increase of pressure and pretension, the flange angle increases. Considering the maximum stress strength, it can be seen from Figure 10 that when the medium pressure is not high (e.g., △P � 0.45 MPa), the stress of the cover plate is less than the maximum value of the material given in the project, and when the pressure increases, the maximum stress may exceed the material yield strength (e.g., △P � 2.0 MPa). e maximum stress in the whole seal occurs at the welding between the stiffener and the cover, and the maximum stress in the bolt is also large. It is necessary to consider whether the bolt strength meets the requirements [13][14][15].  As the sealing performance is related to the stress at the contact surface, Figure 11 shows a schematic diagram of the stress distribution on the gasket contact surface varying with pressure under four pretensions. It can be seen from the figure that, under the following conditions, all actual sealing width is gasket width, and the contact stress of gaskets increases with the increase of radial distance. When the pretension is constant, the contact stresses decreases with the increase of medium pressure. When the pressure is constant, the contact stress increases with the increase of pretension. For further comparison, Figure 12 is a schematic diagram of the variation of the average contact stress with pressure at different pretensions. It can be seen that the average contact stress decreases linearly with the medium pressure. e greater the pretension, the greater the average contact stresses under the same pressure. Under different preloads, when the medium pressure increases from 0.45 MPa to 2.0 MPa, the average contact stress decreases by about 26%.

Fitting Formula of Leakage Rate.
e calculation data and details fitted by the leakage rate formula of the pressure cover of the manhole are shown in Table 3. ODB: Job-719-F=80KN.odb Abaqus/Standard 6.14-4 Wed Nov 18 04:29:50 GMT+08:00 2020

Quantitative Prediction of Leakage
Step: Step-pressure1 bolt pretension, but the structural stress increases at the same time.
erefore, the material damage and failure caused by excessive pretension should be avoided. According to the theoretical analysis of leakage rate, the actual roughness is Ra3.2 [18,19]. Under different pretensions, the intermediate parameters and final leakage rate are calculated, respectively, as shown in Tables 4-7. Figure 13 is a comparison diagram of gasket leakage rate under two different pretensions. It can be seen that, regardless of the pretension, the gasket leakage rate increases with the increase of medium pressure, and the gap of leakage rate between different pretensions increases with the increase of pressure [20][21][22].

Leakage Prediction and Test Verification of
Typical Penetration e leakage rate prediction method of complex sealing structure of containment penetration is shown in Figure 13. In addition to the theoretical prediction of typical penetration, the accuracy of the prediction model is verified by experiments.
ere are five types of typical penetration tested, including a pressure gate, lifting hole, isolation valve, butterfly valve, and ball valve. Figures 14 and 15         concave face outward and convex face inward. e inside of the gate is connected to a fully welded sealed sealing chamber, and the test gas is supplied to the sealing chamber through the air inlet hole. e flow supplement method is adopted in the test, and its schematic diagram is given in Figures 1-3. e gas cylinder fills the test space, and the test pressure is automatically maintained constant through the pressure regulator. e measured supplementary flow is the     Figure 16 displays the test results and theoretical calculation results of sealing rings with two hardness. In the calculation, the flange is 06Cr19Nil0 stainless steel, and surface roughness is 3.2 μm. e diameter is 800 mm. e pretightening compression ratio of the sealing ring is 12.5%. e sealing medium is dry air. e pressure of the sealing chamber gauge increases gradually from 0.05Mpa to 0.6 MPa, which simulates the process of gradual pressure rise under accident conditions. It can be seen from Figure 16, under normal pressure (△P � 0.05 Mpa), the leakage rate of the sealing ring is only 10 −8 kg/s, and the corresponding volume flow is less than 1 ml/min. When the pressure increases to 0.4∼0.5 MPa, the leakage rate rises to the order of 10 −6 kg/s, and the corresponding volume flow is 10∼14 ml/min. ese are meeting with the engineering design requirements. e data    Mathematical Problems in Engineering 13 in Figure 17 show that the theoretically predicted values of the two sealing rings are in good agreement with the experimentally measured values, with the same trend and similar values, which well verifies the accuracy of the prediction method proposed in this study.

Conclusion
(1) In this paper, the constitutive relationship between leakage rate and microstructure parameters is established based on porous media seepage mechanics, and a new prediction model of leakage rate of sealing structure is proposed. en, the relationship between seal specific pressure and microstructure parameters is obtained by Hertz contact mechanics model. e prediction model of penetration leakage rate established by this method is in good agreement with the experimental results, and it can better reflect the effects of rough surface morphology, material mechanical properties, sealing load, high temperature, and high pressure on the leakage rate.
(2) In the calculation, the roughness of the flange is Ra3.2, and the minimum leakage rate can be 4.9E −3 m 3 /day. When other conditions remain unchanged, the leakage rate will decrease significantly with the decrease of roughness. According to the analysis of leakage rate mechanism model, the main factors affecting the leakage rate of pressure cover include medium pressure, bolt preload, flange roughness, and so on. Within the material strength limit, try to select larger preload and lower roughness. (3) e sealing principle of penetration can be classified as contact static seal, where interface leakage is the main leakage path of this kind of sealing structure under normal working condition. e main factors affecting the leakage rate of various penetration parts can be divided into several categories: the first is the micromorphology characterization of the contact surface of the sealing pair, that is, the composite surface roughness. e increase of roughness often leads to the magnitude change of leakage rate. e second type is the mechanical properties of sealing pair materials, especially the mechanical deformation properties of seals. When the working conditions are bad, its aging characteristics also need to be considered. e third type is the macrostructure of the sealing pair, which is often the key to the leakage rate control. e contact stresses distribution of the sealing pair produced by different structures will be quite different. (4) It should be pointed out that there are some uncertain factors in the model calculation: ① preload; ② composite roughness of sealing pair; ③ effects of high temperature, high pressure, and high humidity on mechanical properties of rubber seals. Each of these factors will have a great impact on the leakage rate. erefore, the theoretical research of leakage rate must be based on certain measurement results and experimental data, and then, more detailed prediction model research is performed on various penetrations.

Future Directions
e current research still has certain limitations, and the following aspects need to be further strengthened and improved.  (1) e numerical simulation calculation in this paper is carried out based on making assumptions and simplifies the sealing structure. In order to make the simulation results more accurate, it is necessary to establish a model more in line with the actual situation. (2) Limited by the experimental conditions, this paper selects several typical penetrators for experiments, and more penetration experiments can be added to improve the prediction model in the future. (3) In the future, we will further study the influencing factors affecting the leakage rate, further improve the accuracy and scope of application of this prediction model, refer to the research results of other scholars, and provide a theoretical basis for the development and design of penetration seal structure.

Data Availability
e data used to support the findings of this study are included within the article.