STNI Science and Technology of Nuclear Installations 1687-6083 1687-6075 Hindawi Publishing Corporation 10.1155/2016/1613989 1613989 Research Article A Calculation Method for the Sloshing Impact Pressure Imposed on the Roof of a Passive Water Storage Tank of AP1000 Lu Daogang 1,2 2 http://orcid.org/0000-0001-9293-6225 Zeng Xiaojia 1,2 2 Dang Junjie 3 http://orcid.org/0000-0001-6758-2373 Liu Yu 1,2 2 Tiselj Iztok 1 School of Nuclear Science and Engineering North China Electric Power University Beijing 102206 China ncepu.edu.cn 2 Beijing Key Laboratory of Passive Safety Technology for Nuclear Energy North China Electric Power University Beijing 102206 China ncepu.edu.cn 3 China Nuclear Power Engineering Co., Ltd. Beijing 100840 China 2016 1262016 2016 18 01 2016 06 04 2016 19 05 2016 1262016 2016 Copyright © 2016 Daogang Lu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

There is a large water storage tank installed at the top of containment of AP1000, which can supply the passive cooling. In the extreme condition, sloshing of the free surface in the tank may impact on the roof under long-period earthquake. For the safety assessment of structure, it is necessary to calculate the impact pressure caused by water sloshing. Since the behavior of sloshing impacted on the roof is involved into a strong nonlinear phenomenon, it is a little difficult to calculate such pressure by theoretical or numerical method currently. But it is applicable to calculate the height of sloshing in a tank without roof. In the present paper, a simplified method was proposed to calculate the impact pressure using the sloshing wave height, in which we first marked the position of the height of roof, then produced sloshing in the tank without roof and recorded the maximum wave height, and finally regarded approximately the difference between maximum wave height and roof height as the impact pressure head. We also designed an experiment to verify this method. The experimental result showed that this method overpredicted the impact pressure with a certain error of no more than 35%. By the experiment, we conclude that this method is conservative and applicable for the engineering design.

1. Introduction

As the ultimate heat sink of AP1000 reactor, passive cooling system (PCS) is the key equipment to ensure the safety of nuclear power plant. The large water storage tank, installed at the top of the containment of AP1000, can supply plenty of water for the passive cooling. In the extreme condition, sloshing of the free surface in the tank may impact on the roof and jeopardize structural integrity under long-period earthquake. For the safety assessment of structure, it is necessary to calculate the impact pressure caused by water sloshing.

Since the behavior of sloshing impacted on the roof is involved into a strongly nonlinear phenomenon, calculations of the impact pressure with theoretical or numerical method currently are of difficulties. Ibrahim  focused on a 2D tank with simple geometry to solve the linear sloshing problems using analytical methods. A numerical model using finite element technique was presented by Pal et al.  to study the linear behavior of cylindrical tanks. Choun and Yun  used the velocity potential and the linear water wave theory to decompose the surface wave into multiple forms. There are some complex methods to solve the nonlinear sloshing problems. Li et al.  used an improved material point method (MPM) to predict the liquid impact force by a contact algorithm. Liquid sloshing experiments in a partially watered square tank were proposed to validate the results of simulation. Eswaran et al.  proposed a numerical method based on volume of fluid (VOF) techniques with arbitrary-Lagrangian-Eulerian (ALE) formulation to analyze baffled and unbaffled tanks with a nonlinear sloshing behavior. However, these researches are usually valid for simple cases with linear or weakly nonlinear liquid sloshing dynamics.

But it is applicable to calculate the height of sloshing in a tank without roof. Fujita et al.  utilized the velocity potential theory to analyze the liquid sloshing in the annular region of more intricate coaxial circular cylinders. Formulas about the maximum wave height ( η m a x ) at shell wall and the maximum pressure ( p m a x ) at the free surface were obtained. More interestingly, the correlation between maximum wave height and maximum pressure was p m a x = ρ g η m a x from equation (34) in their research. Virella et al.  used the finite element package ABAQUS to investigate the free surface wave amplitude and pressure distribution of tank wall by both linear wave theory and nonlinear wave theory models. Nayak and Biswal  used the Galerkin-weighted-residual based finite element method (FEM) to solve Laplace equation with nonlinear boundary conditions. The wave height of nonlinear sloshing was verified to be accurate.

Besides, the impact pressure is an important parameter in the assessment of safety of engineering design. Researchers had conducted large-scale experiments to investigate impact pressure .

As is mentioned above, the available studies mostly focused on the sloshing characteristics of rectangular tanks with simple geometry. However, considering the special structure of PCCWST, which is a coaxial circular cylinder tank with an inclined bottom, it is difficult in obtaining analytical expressions for the prediction of the natural modes and the liquid motion. Moreover, numerical and analytical methods to precisely describe the sloshing impact pressure are complicated because of the significant nonlinearity phenomena. In the present paper, a simplified method was proposed to calculate the impact pressure using the sloshing wave height. Moreover, an experiment was designed to verify this method.

2. Calculation Method

Sloshing of the water surface in the tank may impact on the roof under long-period earthquake and the prediction of the impact pressure is necessary. Due to the apparent nonlinear behavior caused by complicated sloshing phenomena, a simplified method was proposed to calculate the impact pressure using the sloshing wave height, in which we first marked the position of the height of roof, the solid wide line showed in Figure 1, then produced sloshing in the tank without roof and recorded the maximum wave height which can be shown as the oblique line in Figure 1, and finally regarded approximately the difference between maximum wave height and roof height as the impact pressure head. In this way, the impact pressure can be calculated by the following equation: (1) P = ρ g h m a x - h , where P is the maximum impact pressure, ρ is the density of water, g is the acceleration due to gravity, h m a x is the maximum wave height, and h is the distance between the static water surface and the roof.

Schematic of sloshing.

Besides, the impact force of water in a tank at any time is given by the following equation: (2) F w = s 1 s 2 P d s , where F w is impact force, P is the impact pressure, and s 1 , s 2 , respectively, mean the lower and upper boundary of tank roof.

By utilizing the maximum sloshing wave height, the maximum impact pressure, which appeared near the outer wall, can be calculated by (1).

The impact area is divided into n parts and the impact pressures in each part can be obtained by linear interpolation. The i th pressure is given by the following equation: (3) p i = i n ρ g h m a x - h .

According to Figure 1, L , which is related to the impact area, can be calculated by (4) L R = h m a x - h h m a x .

The total impact force can be written as (5) F t = i = 0 n - 1 - R 2 - y 2 R 2 - y 2 d x R - L + i L / n R - L + i + 1 L / n p i + p i + 1 2 d y i = 0,1 , 2,3 , , n - 1 .

Substituting (3) and (4) to (5), we can get the total impact force.

3. Experiment

In order to verify the calculation method, an experiment was designed.

3.1. Scaling Model

Since the passive storage tank of AP1000 is large, experiment with the large size of tank is difficult. In order to keep the similar properties of fluid, a scaling model is proposed. For the case of the present experiment, some main governing dimensionless numbers that should be respected are as follows [6, 13].

Froude number: F r = V 2 / g L ; Reynolds number: R e = V L / ν ; Strouhal number: S t = f L / V ; Euler number: E u = P / ρ V 2 ; Bond number: B o = ρ g L 2 / σ .

With the analysis, the Froude number, Strouhal number, Euler number, and Bond number are respected and the Re number is relaxed when ignoring viscosity effects. The scaling factors for the different parameters in sloshing experiments are summarized in Table 1 and these expressions are called the modified linear scaling method.

Scaling factors for sloshing experiments.

Parameter Scaling ratio
Length: L [m] C l
Gravity: g [m/s2] 1
Velocity: V [m/s] C l 1 / 2
Frequency: f [Hz] C l - 1 / 2
Time: t [s] C l 1 / 2
Kinematic viscosity: ν [m2/s] C l 3 / 2
Pressure: P [Pa] C ρ C l
Wave height: H [m] C l
3.2. Experimental Equipment

According to the similarity criteria, an experimental model (Figure 2) of PCCWST was built with the scale of 1 : 40. The experimental tank was mainly cylindrical (0.122 m internal radius) with an inclined bottom and the vessel was made of acrylic sheet to observe the sloshing phenomena clearly. The tank was considered to be rigid due to the 15 mm thickness of outer wall. This assumption was made to simplify the study of sloshing. In order to measure the maximum sloshing wave height, the length scaling for the height of the tank was 3 : 40. A summary of the main sizes of experimental tank is presented in Figure 3.

Experimental setup.

Sensor layout.

The internal structures, showed as the red parts in Figure 3, were added suspended from the top cover, which was fixed on the outer structure. Such internal structures contained a supporting configuration and an in-vessel circle roof mounted on the supporting configuration by 4 connecting rods with knurled caps. The whole experimental tank was glued to a square shell which was screwed in 28 threaded holes on the shaking table.

The experimental setup included a computer, a laser displacement sensor, 3 pressure sensors, an acceleration sensor, a data acquisition tool, a set of commercial acquisition software COINV DASP, and some wires. The layouts and numbers of sensors are illustrated in Figure 3.

3.3. Experimental Procedure and Results

Tabulating in Table 2, a total of 6 experiment cases involving 4 different filling depths were carried out on the experimental setup with two configurations, first without internal components (a supporting configuration and an in-vessel roof) and then with internal components. For each case, sloshing experiments under multiple excitation amplitudes were carried out and several kinds of transient sloshing responses like sloshing wave height and slosh-induced impact pressure were recorded by data acquisition system. Measurements of the wave height are presented in Section 3.3.2, in cases without internal components. Impact pressure of sloshing is discussed in Section 3.3.3 in cases with internals.

Experiment sequences.

Case number Filling depth H [mm] Distance between static water surface and in-vessel roof h [mm] Excitation amplitude A [mm]
1 217 43 12, 24, 36
2 177 83 12, 24, 36
3 197 63 12, 24, 36
4 237 23 12, 24, 36

: H and h are based on the length scale.

3.3.1. Natural Frequency of Sloshing in First Mode

Liquid motion inside a tank in periodic movement is dominated by the liquid filling ratio, the amplitude and frequency of the tank motion, and the geometry of the tank. Considering the complex geometry of the experimental tank, the fundamental natural frequency of liquid was measured through sweep frequency response analysis (SFRA) experiment. The tank was partially filled with water, which was mixed with gouache paint and covered by pearl cotton on the surface. It was aimed to make the laser displacement sensor work normally by detecting object with surface that exhibited a diffuse reflection. From the research of Lu et al. , the influence of water level on natural frequency of sloshing water was quite small. Hence, the subtle variations of natural frequency with different filling depths can be ignored. The tank was filled up to a water level of 217 mm, which was based on the length scale, for the measuring of liquid frequency. For the purpose of looking in detail at the resonance frequency of liquid, a logarithmic sweep with a range of 0.1~20 Hz was adopted. In addition, there are two types of control mode for driving the shaking table: Acceleration Control (AC) and Displacement Control (DC). At low frequency, about 0.1~4 Hz, AC is used and the displacement was controlled at the value of 2.0 mm, since the displacement sensor achieved higher accuracy than acceleration sensor at low frequency range. However, DC is invoked in high frequency of 4~20 Hz and the acceleration was set with 0.650 m/s2. The time history of sloshing displacement was transformed into frequency spectrum by a fast Fourier transform (FFT). By this analysis, the fundamental natural frequency of the liquid was evaluated to be 0.75 Hz.

For an annular region formed by coaxial circular cylinder tank in periodic motion, the resonant frequencies of the liquid can be calculated from the linear potential theory as [6, 13] (6) ω m n 2 = g a ξ m n tanh ξ m n H a m = 0,1 , 2 , , n = 1,2 , 3 , , , where g is the gravitational acceleration, a is the outer radius of the tank, ξ m n are the positive roots related to the ratios of inner radius to outer radius of the tank, and H is the still filling depth.

The calculating fundamental natural frequency from (6) is 0.887 Hz, a little larger than the experimental result. It can be explained by the inclined bottom configuration of the tank.

3.3.2. Sloshing without the Roof

Due to the fact that the periods of sloshing water ( 1 / f = 1 / 0.75 = 1.33 s ) were far more longer than the site predominant period (0.34 s), the resonant three-cycle sine wave (RTCS) was assumed to be the input excitation according to the HAF0102 of nuclear safety guide in China, which is also a widely used evaluation method of seismic response in Japan. The load function can be written as X = A sin ( 2 π f t ) , where X is the displacement of shaking table, A is the amplitude of the sine wave, f is the water natural frequency, and t is the time which lasts for three periods of sine waves. In our experiments, the frequency of RTCS was set as 0.75 Hz. When the frequency of excitation closed to the natural frequency, the liquid inside the experimental tank can exhibit strong oscillations, thus exerting a greater extent of sloshing and violent impact loads on the tank.

The experimental setup without internals was carried out for the experiment of seismic response about sloshing wave height. The sloshing wave height was measured through laser displacement sensor which was installed at point 1 near the outer tank wall which is showed in Figure 3. The acceleration transducer was installed at the bottom of the tank, in order to monitor the initial and final points of RTCS. Figure 4 reveals the vertical displacement of the wave in case 1 with 24 mm amplitude. It can be seen that the crest was almost 2 times size of the trough, which embodied an apparent feature of nonlinearity. Table 3 lists the maximum wave heights of case 1 to case 4. Moreover, through the analysis of these data, for each case, the maximum wave heights increased as the excitation amplitudes grew and the increasing tendency gradually reduced, which was also a nonlinear behavior.

The maximum wave height.

Amplitude [mm] Case  1 Case  2 Case  3 Case  4
12 79.191 65.272 67.103 98.9172
24 227.355 163.953 190.414 264.142
36 316.327 217.348 255.602 355.023

The vertical displacement time history of sloshing wave of case 1 with 24 mm amplitude.

3.3.3. Sloshing with Roof

At the second part of the experiment, experiments with internals were carried out to measure the impact pressure of in-vessel roof. CYY8 dynamic pressure sensors with high frequency response were installed at measuring points 2, 3, and 4 on the roof, as showed in Figure 3, to present the radial pressure distribution.

Pressure variations were measured by dynamic piezoresistive pressure transducers in cases 1 to 4. Figure 5 illustrates the impact pressures time history results in case 1 of measuring point 4 with 24 mm amplitude. According to the curve, the maximum impact pressure appeared when it came to the third wave of RTCS and then 2 small peaks came after it. These two small peaks can be explained by the inertia effect of water, which made water impact the inner roof despite the end of excitation. In addition, the gap between in-vessel roof and the tank wall would lead to a small portion of water splashing upward to the top surface of cap. Therefore, some negative values appeared intermittently in the measuring data showed in Figure 5.

Time history of impact pressure about case 1 of point 4 with 24 mm amplitude.

As is shown in Figure 6, impact pressure distributions along the radial direction are presented in cases with different RTCSs and the same filling depth, 217 mm. The impact pressure of point 2 was the minimum and the impact pressure of point 4 was the maximum, which showed that the impact pressure increases from inner wall to outer wall.

The distribution of impact pressure of case 1 with 24 mm amplitude.

During the operational process of experiments, there were some errors in experimental data owing to erroneous human input and the accuracy of measuring instruments to a certain extent. The use of laser displacement sensor, pressure sensors, and the indigenous vibration of shaking table caused significant errors in measurements. Besides, there also existed electromagnetic interferences during the operational process of shaking table, which adversely affected the accuracy of sensors especially pressure sensors, also contributed to the errors’ list.

By computing the area integral of pressure, the experimental force data can be derived from the measuring pressure data with some conservative assumptions.

In our research, as illustrated in Figure 7, the area from inner wall to outer wall along y axis is divided into 3 parts and the pressure of each part is assumed to be the measuring value of the sensor which is installed at the corresponding part. Therefore, the impact force in each part can be written as in (7). The total impact force is the sum of that in each part which is showed in (8). Consider the following: (7) F i = - R 2 - y 2 R 2 - y 2 d x y i y i + 1 p i d y , (8) F t = i = 1 n F i i = 1,2 , 3 , where F i and p i are the impact force and impact pressure in each part, R is the outer radius of cylinder tank, and F t is the total impact force.

Schematic of area dividing.

In our cases, the total impact forces of different conditions are listed in Table 4.

Total impact force from measuring data.

Case no. Filling depth [mm] Amplitude [mm] Impact force [N]
1 217 mm 12 11.31001
24 60.37129
36 91.28681

2 177 mm 12 4.370317
24 25.62775
36 41.8155

3 197 mm 12 7.423492
24 36.49382
36 56.19478

4 237 mm 12 23.1682
24 84.23761
36 117.9033
4. Analysis of Experimental Results

In Section 2, we proposed a simplified method to estimate the impact pressure. So we list the results of simplified method and measured data. We used the calculated impact pressure (CIP) indirectly to represent the results of simplified method. The measured data were called the measured impact pressure (MIP) directly. The comparison between calculated pressure and measured pressure is listed in Table 5.

The comparison between calculated pressure and measured pressure.

Case no. Amplitude [mm] Measured pressure (MIP) [Pa] Calculated pressure (CIP) [Pa] Error (%)
1 12 500.62 354.6718 −29.1536
24 1549.30 1806.679 16.6126
36 2074.05 2678.605 29.14851

2 12 263.27 0 −100
24 913.30 793.3394 −13.1349
36 1243.44 1316.61 5.884514

3 12 397.514 40.209 −89.8856
24 1221.34 1248.657 2.236658
36 1579.48 1887.5 19.50133

4 12 757.21 743.9886 −1.74607
24 1753.86 2363.192 34.74232
36 2445.99 3253.825 33.02691

According to Table 5, CIP are lower than MIP under small amplitudes of RTCS. A possible reason is the influence of air motion. In our experiment, the water tank cannot be thoroughly sealed. Therefore, air in the ring cavity can move with the vigorous sloshing of water. Air motion is the primary influential factor under low level of excitation. However, it is considered that small excitation has really subtle effect in engineering practice. As for the higher excitation level, CIP are larger than MIP under medium and high level of excitations. With a certain error of no more than 35%, this method of impact pressure calculation is considered to be overpredictive. Therefore, similar to the results of Lu’s  research, it can be concluded that the calculation method is reasonable and relatively conservative.

According to the simplified calculation equation for impact force proposed in Section 2, we obtained the impact force of simplified method which was called calculated force here. Showed in Section 3.3.3, the experimental force data was called experimental force here. Figure 8 compares the calculated force and experimental force. The calculated force is lower than, though generally coinciding with, the experimental force in smaller wave height cases. However, calculated force becomes close to the experimental force and then exceeds it due to the increase of maximum sloshing wave height.

Comparison of calculated force and experimental force.

By the comparison of calculated results and experimental results, some future works are left. First, there should be some holes on the in-vessel roof to reduce the influence of air motion during experimental procedure. Second, longer duration time experiments with more pressure sensors and displacement sensors are needed to investigate the impact pressure distribution and computation of impact force.

5. Conclusions

In order to estimate the impact pressure imposed on the roof of water storage tank of AP1000, we proposed a simplified method, which regarded approximately the difference between maximum wave height and roof height as the impact pressure head.

In order to verify the above method, we built a scaled experimental facility.

Using the experimental facility, we measured the maximum sloshing wave height in the case of not having roof and the impact pressure in the case of having roof.

By the above experiment, the simplified method was considered to be conservative with a certain error of no more than 35%.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The project was sponsored by National Science and Technology Major Project of the Ministry of Science and Technology of China (2015ZX06004002-003).

Ibrahim R. A. Liquid Sloshing Dynamics: Theory and Applications 2005 New York, NY, USA Cambridge University Press Pal N. C. Sinha P. K. Bhattacharyya S. K. Finite element coupled slosh analysis of rectangular liquid filled laminated composite tanks Journal of Reinforced Plastics and Composites 1999 18 15 1375 1407 2-s2.0-0033313269 Choun Y.-S. Yun C.-B. Sloshing analysis of rectangular tanks with a submerged structure by using small-amplitude water wave theory Earthquake Engineering and Structural Dynamics 1999 28 7 763 783 10.1002/(SICI)1096-9845(199907)28:7<763::AID-EQE841>3.0.CO;2-W 2-s2.0-0032839563 Li J. G. Hamamoto Y. Liu Y. Zhang X. Sloshing impact simulation with material point method and its experimental validations Computers & Fluids 2014 103 86 99 10.1016/j.compfluid.2014.07.025 2-s2.0-84906064525 Eswaran M. Saha U. K. Maity D. Effect of baffles on a partially filled cubic tank: numerical simulation and experimental validation Computers & Structures 2009 87 3-4 198 205 10.1016/j.compstruc.2008.10.008 2-s2.0-58149492530 Fujita K. Ito T. Okada K. Seismic response of liquid sloshing in the annular region formed by coaxial circular cylinders Engineering Computations 1985 2 4 299 306 10.1108/eb023628 2-s2.0-0022201021 Virella J. C. Prato C. A. Godoy L. A. Linear and nonlinear 2D finite element analysis of sloshing modes and pressures in rectangular tanks subject to horizontal harmonic motions Journal of Sound and Vibration 2008 312 3 442 460 10.1016/j.jsv.2007.07.088 2-s2.0-40249090507 Nayak S. K. Biswal K. C. Nonlinear seismic response of a partially-filled rectangular liquid tank with a submerged block Journal of Sound and Vibration 2016 368 148 173 10.1016/j.jsv.2016.01.010 Akyildiz H. Ünal E. Experimental investigation of pressure distribution on a rectangular tank due to the liquid sloshing Ocean Engineering 2005 32 11-12 1503 1516 10.1016/j.oceaneng.2004.11.006 2-s2.0-18744362160 Song Y. K. Chang K.-A. Ryu Y. Kwon S. H. Experimental study on flow kinematics and impact pressure in liquid sloshing Experiments in Fluids 2013 54 9, article 1592 10.1007/s00348-013-1592-5 2-s2.0-84883267664 Souto-Iglesias A. Bulian G. Botia-Vera E. A set of canonical problems in sloshing. Part 2: influence of tank width on impact pressure statistics in regular forced angular motion Ocean Engineering 2015 105 136 159 10.1016/j.oceaneng.2015.06.016 2-s2.0-84936103580 Pistani F. Thiagarajan K. Experimental measurements and data analysis of the impact pressures in a sloshing experiment Ocean Engineering 2012 52 60 74 10.1016/j.oceaneng.2012.06.002 2-s2.0-84864444238 Abramson H. N. The dynamic behavior of liquids in moving containers, with applications to space vehicle technology 1966 NASA-SP-106 Washington, DC, USA NASA Lu D. Liu Y. Zeng X. Experimental and numerical study of dynamic response of elevated water tank of AP1000 PCCWST considering FSI effect Annals of Nuclear Energy 2015 81 73 83 10.1016/j.anucene.2015.03.019 2-s2.0-84926226113 Lu D. Sloshing response of the free surface in the main vessel of CEFR excited by 3 sine waves Chinese Journal of Nuclear Science and Engineering 2003 23 306 310