# Experimental Investigation on Flow-Induced Vibration of Fuel Rods in Supercritical Water Loop

Authors: Licun Wu; Daogang Lu; Yu Liu
Journal: Science and Technology of Nuclear Installations
 (2014)

---

## Abstract

The supercritical water-cooled reactor (SCWR) is one of the most promising Generation IV reactors. In order to make the fuel qualification test for SCWR, a research plan is proposed to test a small scale fuel assembly in a supercritical water loop. To ensure the structure safety of fuel assembly in the loop, a flow-induced vibration experiment was carried out to investigate the vibration behavior of fuel rods, especially the vibration caused by leakage flow. From the experiment result, it can be found that: the vibration of rods is mainly caused by turbulence when flow rate is low. However, the effects of leakage flow become obvious as flow rate increases, which could changes the distribution of vibrational energy in spectrum, increasing the vibrational energy in high-frequency band. That is detrimental to the structure safety of fuel rods. Therefore, it is more reasonable to improve the design by using the spacers with blind hole, which can eliminate the leakage flow, to assemble the fuel rods in supercritical water loop. On the other hand, the experimental result could provide a benchmark for the theoretical studies to validate the applicability of boundary condition set for the leakage-flow-induced vibration.


---

## Body

## 1. Introduction

The supercritical water reactor is a Generation IV reactor concept which uses supercritical water as the working fluid. The high thermodynamic efficiency and plant simplification make SCWR attractive for consideration as a promising advanced nuclear system. In order to develop a viable design for the core, accurately estimate the heat transfer coefficient, and develop materials for the fuel and core structure that will be sufficiently corrosion-resistant to withstand SCWR conditions, a research plan, which is made by the steering committee of the SCWR system in the Generation IV International Forum, is proposed to test a small scale fuel assembly in a supercritical water loop in the LVR-15 research reactor, which is located in Nuclear Research Institute in Rez, Czech Republic. This water cooled reactor with an enough core height enables replacing one of its assemblies with a pressure tube containing four fuel rods, as Figure1 shows. These four fuel rods can reach a fissile power of more than 50 kW. A recuperator inside the pressure tube, situated right above the fuel rods, is boosting the feed-water temperature of 300°C to typical evaporator conditions. A cooler in the top section of the pressure tube acts as the heat sink to remove the fissile and gamma power. The pressure tube will be connected with coolant pumps and safety and auxiliary systems, forming the supercritical water loop to simulate a supercritical water environment at the fuel assembly while the rest of the core operates at ambient pressure (Figure 2). Before carrying out the project, a series of safety analyses on the small scale fuel assembly will be made, including the experimental investigation on the flow-induced vibration of fuel assembly, which is presented in this paper.Figure 1
Concept of the active channel for SCWR fuel qualification.Figure 2
Schematic diagram of the primary loop and of the safety system for the SCWR fuel qualification test.In general, there is a tiny clearance between the fuel rod and the spacer to make charging easier, and the wear between the fuel rods and spacer will make the clearance become larger over time. Consequently, the fuel rods would vibrate in case of coolant flowing through the narrow clearance. This phenomenon is called leakage-flow-induced vibration, which has been studied by many researchers (Blevins [1]; Hobson [2]; Chen [3]). Païdoussis [4] wrote a review to introduce some serious accidents in the nuclear reactors caused by it.Mateescu and Paidoussis [5] researched the unsteady potential flow in a narrow annular passage, which is formed by a motionless rigid duct and an oscillating rigid center-body, both of axially variable cross-sections. The fluid dynamic force exerted on the center-body was obtained on the basis of unsteady potential flow theory. It was concluded that the stability of center-body decreases as the pivot of rotation is shifted towards the downstream end of the center-body, and a divergent annular passage has a destabilizing effect. Arai and Tajima [6] studied a similar problem by using a different approach.Li et al. [7] developed a numerical method to analyze the flutter instability of the flow vibration of the inner cylinder. Meanwhile, an experiment on the annular leakage-flow-induced vibration was carried out, and a critical flow rate of the flutter instability was obtained for several annular leakage-flow systems with different passage increment ratios as well as the eccentricities. They found that a divergent annular passage may cause a negative fluid damping, so the flutter instability will happen when the flow rate is larger than the critical flow rate.Langthjem et al. [8] investigated the linear stability of a flexible, cylindrical rod exposed in annular leakage flow. It is found that when flow speed reaches a certain critical value, simply supported rods may become unstable by either flutter or divergence.Fujita et al. [9] investigated the dynamical behavior of an axisymmetric elastic beam subjected to axial leakage flow. By making complex eigenvalue analysis, the variation of the dynamic behavior during pre- and postinstability is researched with respect to increasing axial leakage-flow velocity. The relationships between the unstable phenomena and axial flow velocities are clarified for the transition from the lower mode to higher mode.Since the study on leakage-flow-induced vibration of the fuel rods has not been reported previously, and the stability associated with the leakage flow mechanism strongly depends on the detailed geometries and fixed mode of the structure, it is necessary to perform an experiment to investigate the effects of leakage flow on the vibration of the small scale fuel rods in supercritical water loop. The present research is an extension of previous studies and aims at investigating the vibration characteristics of fuel rods under the influence of leakage flow. On one hand, the experimental results could provide a reference for optimal design of the fuel assembly in supercritical water loop. On the other hand, it could be treated as the benchmark to validate the availability of mathematical model built for the leakage-flow-induced vibration.


## 2. Theoretical Models

Consider a circular fuel rod immersed in the coolant flowing at velocityU parallel to the z axis (Figure 3). The fuel rod has flexural rigidity EI, linear density (mass per unit length) m, and total length l. All motions of the fuel rod are assumed to be confined in the y-z plane. According to the studies of Chen [3], the equation of motion is written as
(1)EI∂4u∂z4+ηI∂5u∂t∂z4+maU2∂2u∂z2-γT0∂2u∂z2-12CTmaU2D[(1-12γ)l-z]∂2u∂z2-12(1-γ)CT′maU2∂2u∂z2+2maU∂2u∂z∂t+12CNmaU2D∂u∂z+12CNmaUD∂u∂t+CV∂u∂t+(m+ma)∂2u∂t2=g(x,t),
where u is the displacement of fuel rod, E is the modulus of elasticity, ma is the added mass, T0 is the initial axial tension, CT′ is the form drag coefficient at the free end, CV is the effective viscous damping coefficient, γ=1 if the downstream end is supported so the displacement is zero, and γ=0 if it is unsupported or elastically supported. CT is equal to CN and defined as
(2)CN=CT=4CfπCm,
where Cf is the drag coefficient due to shear forces and Cm is added mass coefficient.Figure 3
Schematic diagram of the fuel rod in axial flow.The appropriate boundary conditions associated with the equation of motion are as follows.Forz=0,
(3)k1u+EI∂3u∂z3=0,c1∂u∂z-EI∂2u∂z2=0
and for z=l,
(4)k2u-EI∂3u∂z3=0,c2∂u∂z+EI∂2u∂z2=0,
where k1 and k2 are displacement spring constants and c1 and c2 are torsional spring constants (see Figure 3).Unfortunately, (1) associated with boundary conditions (3) and (4) can only be used for the studies about linear vibration of the fuel rod. Since the boundary condition with leakage flow is nonlinear and complex, (3) and (4) are inappropriate as the boundary conditions for the calculation about leakage-flow-induced vibration of the fuel rod. On the other hand, since the boundary conditions are not taken into account during (1) derivation process, (1) is still applicable for calculation of leakage-flow-induced vibration. Therefore, only a series of appropriate boundary conditions are needed for (1) to accurately calculate the vibration behavior of fuel rod caused by leakage flow, and the experimental result in this paper could provide a benchmark to validate the availability of boundary conditions.


## 3. Details of Experiment

### 3.1. Experimental Setup

The experimental facility is shown in Figure4. The flow-induced vibration test section consists of two parts: (1) a visualized channel made of acrylic with water flowing axially through it and (2) a small scale fuel assembly. The fuel assembly in the experiment has the same geometry parameters as that in supercritical water loop. The cladding is made of stainless steel 316L, which is the same as the material used in supercritical water loop. The material of pellet in the experiment is lead, whose density (11.3 g/cm3) is similar to the density of uranium dioxide (10.96 g/cm3). In order to demonstrate the rationality of using lead as a replacement for uranium dioxide, the modal analysis was made for the fuel rod. From the results of modal analysis listed in Table 1, it is easy to see there is little difference between the natural frequencies of fuel rods with two kinds of pellets made of uranium dioxide and lead, respectively. Since the vibration behavior of fuel rods is strongly dependent on their natural frequencies, it is proper to use lead as a replacement for uranium dioxide.Table 1
Natural frequency of the fuel rod.
Pellet material
First order
Second order
Third order
Fourth order
Uranium dioxide
38.1436 Hz
38.1579 Hz
111.947 Hz
111.966 Hz
Lead
37.738 Hz
37.752 Hz
110.76 Hz
110.78 HzFigure 4
Schematic diagram of experimental loop.The bottom ends of the rods are fully constrained. Meanwhile, there is a 0.2 mm clearance between the top end of the rod and the locating hole. The displacements at midpoint of fuel rod are detected by a laser displacement sensor. The geometry parameters of the test section and experimental conditions are listed in Table2. In addition, the assembly drawing and photo of the test section are shown in Figure 5.Table 2
Experimental conditions and geometry parameters of the test section.
Temperature of water
25°C
Operational pressure
0.1 Mpa
Size of channel
20 mm × 20 mm × 725 mm
Diameter of fuel rod
8 mm
Length of fuel rod
710 mmFigure 5
The collocation and geometry of the fuel assembly and test section.In order to investigate the effects of leakage flow on the characteristics of flow-induced vibration, the effect of two different spacers, which is used to constrain the top end of the rod, is researched, respectively. Figure6(a) illustrates the top end of fuel rod is constrained by the spacer with through holes, meaning there is a leakage flow across the clearance. This case is called “through hole case.” In Figure 6(b), since the top end of fuel rod is constrained by the spacer with blind holes, there is no occurrence of leakage flow. This case is called “blind hole case.”Schematic diagram of two types of spacers.
(a)
(b)


### 3.2. Measuring Instrument Calibration

For the small size of channel, the contact measurement would have a significant influence on the flow field in the channel. Thus, the laser displacement sensor is used in the experiment to perform the noncontact measuring without influence on the flow field. It employs triangulation measurement principles whereby it projects a beam of visible laser light that creates a spot on a target surface. Reflected light from the surface is viewed from an angle by a digital camera inside the sensor, and the target’s displacement is computed from the image pixel data.However, in the process of measuring the vibration of fuel rod fixed in a channel with room-temperature water flowing inside, the refraction caused by acrylic, water, and air would lead to the measurement deviation. For accurately measuring the vibration of the fuel rods, it is necessary to correct the deviation of measurement caused by refraction before carrying out the flow-induced vibration (FIV) experiment.Generally, the measuring instrument calibration can be performed by refraction theoretical analysis for laser displacement sensor. However, since the channel in practice cannot be made evenly by the acrylic material, the nonuniformity of channel could lead to the refraction occurring not only on the surface but also inside the channel, and this phenomenon is difficult to predict in theory. Thus, a device is employed to correct the measurement deviation caused by refraction. It is composed of a micrometer and a container made of acrylic, as Figure7 shows.Figure 7
Measurement calibration device.The thickness of the containerd1 is the same as the thickness of channel d1′ in FIV experiment; meanwhile, the distance from the micrometer to the inner face of the container d2 is equal to that from fuel rod to the inner wall of the channel d2′. Therefore, the measured deviation caused by refraction is the same when the sensor is used to measure the movements of the fuel rods and micrometer, respectively. When the top end of micrometer moves forward to a little distance Δx, a corresponding measured value Δy can be obtained by the laser displacement sensor. The difference between Δx and Δy is the measured deviation caused by laser refraction. Several movements would be made to get two sets of data, Δxn(n=1,2,3,4…) and Δyn(n=1,2,3,4…). The calibration curve for laser displacement sensor could be obtained by a linear fitting with the two sets of data, as Figure 8 shows. The equation of the curve is written as
(5)y=0.72582x-17.90196,
where x is the actual position of the micrometer and y is the measured position.Figure 8
Measurement calibration curve.Assuming the actual position isx1, the measured position is y1; then (5) becomes
(6)y1=0.72582x1-17.90196.Moreover, if the actual position isx2 and the measured position is y2, (5) becomes
(7)y2=0.72582x2-17.90196.Subtracting (7) to (6) gives
(8)y2-y1=0.72582(x2-x1).Noting(x2-x1) is the actual displacement of micrometer Δx and (y2-y1) is the corresponding measured displacement Δy, so (8) becomes
(9)Δy=0.72582·Δx.From (9), it is clear that when the micrometer moves forward 1 μm, the corresponding measured value by displacement sensor is 0.72582 μm, with a deviation of 0.27418 μm, which is caused by the refraction. For eliminating the deviation, it is necessary to correct the measured value of laser displacement sensor according to (9).


### 3.3. Structure Improvement

Generally, the vibration behavior of fuel rods is random in FIV experiment, including forward motion and lateral motion. Since the fuel rods are centrosymmetric structure, statistical characteristics of the vibration in forward direction and lateral direction are the same. Thus, it is only needed to measure the forward motion during the experiment to know the vibration characteristic of fuel rods. However, when a fuel rod vibrates in both directions, as shown in Figure9(a), the position of the laser spot projected by the sensor will change with the vibration of fuel rod, moving along the surface of cylindrical rod and causing the measured value d1 to be larger than the actual displacement d2. That has a substantial influence on the result. In order to solve this problem, the structure of fuel rod has been improved by fixing a tiny metal sheet on the surface of fuel rod. The size of metal sheet is shown in Figure 9(c), and one side of the sheet is flat. In this way, even though the position of laser spot still changes with the movement of fuel rod, the laser spot could move on a flat surface instead of curved surface, as shown in Figure 9(b). As a result, in this situation the measured value d1′ is the same as the actual displacement d2′.Schematic of improved structure.
(a)
Original structured1>d2
(b)
Improved structured1′=d2′
(c)
Size of metal sheet


### 3.4. Measuring System

The hardware of the measuring system mainly consists of a noncontact and high-speed laser displacement sensor (LTC-025-04-SA, MTI), a high-speed data acquisition card (Daq3000, Iotech), and data acquisition soft system (DASYLab, National Instruments). The measuring system is shown in Figure10. In order to reduce the sampling signal bias and frequency alias, a low-pass filter is used to process the detected signal for preventing the high-frequency signal. In addition, a correction module is added to correct the measurement deviation caused by refraction according to (9).Figure 10
Schematic diagram of data acquisition system.


## 3.1. Experimental Setup

The experimental facility is shown in Figure4. The flow-induced vibration test section consists of two parts: (1) a visualized channel made of acrylic with water flowing axially through it and (2) a small scale fuel assembly. The fuel assembly in the experiment has the same geometry parameters as that in supercritical water loop. The cladding is made of stainless steel 316L, which is the same as the material used in supercritical water loop. The material of pellet in the experiment is lead, whose density (11.3 g/cm3) is similar to the density of uranium dioxide (10.96 g/cm3). In order to demonstrate the rationality of using lead as a replacement for uranium dioxide, the modal analysis was made for the fuel rod. From the results of modal analysis listed in Table 1, it is easy to see there is little difference between the natural frequencies of fuel rods with two kinds of pellets made of uranium dioxide and lead, respectively. Since the vibration behavior of fuel rods is strongly dependent on their natural frequencies, it is proper to use lead as a replacement for uranium dioxide.Table 1
Natural frequency of the fuel rod.
Pellet material
First order
Second order
Third order
Fourth order
Uranium dioxide
38.1436 Hz
38.1579 Hz
111.947 Hz
111.966 Hz
Lead
37.738 Hz
37.752 Hz
110.76 Hz
110.78 HzFigure 4
Schematic diagram of experimental loop.The bottom ends of the rods are fully constrained. Meanwhile, there is a 0.2 mm clearance between the top end of the rod and the locating hole. The displacements at midpoint of fuel rod are detected by a laser displacement sensor. The geometry parameters of the test section and experimental conditions are listed in Table2. In addition, the assembly drawing and photo of the test section are shown in Figure 5.Table 2
Experimental conditions and geometry parameters of the test section.
Temperature of water
25°C
Operational pressure
0.1 Mpa
Size of channel
20 mm × 20 mm × 725 mm
Diameter of fuel rod
8 mm
Length of fuel rod
710 mmFigure 5
The collocation and geometry of the fuel assembly and test section.In order to investigate the effects of leakage flow on the characteristics of flow-induced vibration, the effect of two different spacers, which is used to constrain the top end of the rod, is researched, respectively. Figure6(a) illustrates the top end of fuel rod is constrained by the spacer with through holes, meaning there is a leakage flow across the clearance. This case is called “through hole case.” In Figure 6(b), since the top end of fuel rod is constrained by the spacer with blind holes, there is no occurrence of leakage flow. This case is called “blind hole case.”Schematic diagram of two types of spacers.
(a)
(b)


## 3.2. Measuring Instrument Calibration

For the small size of channel, the contact measurement would have a significant influence on the flow field in the channel. Thus, the laser displacement sensor is used in the experiment to perform the noncontact measuring without influence on the flow field. It employs triangulation measurement principles whereby it projects a beam of visible laser light that creates a spot on a target surface. Reflected light from the surface is viewed from an angle by a digital camera inside the sensor, and the target’s displacement is computed from the image pixel data.However, in the process of measuring the vibration of fuel rod fixed in a channel with room-temperature water flowing inside, the refraction caused by acrylic, water, and air would lead to the measurement deviation. For accurately measuring the vibration of the fuel rods, it is necessary to correct the deviation of measurement caused by refraction before carrying out the flow-induced vibration (FIV) experiment.Generally, the measuring instrument calibration can be performed by refraction theoretical analysis for laser displacement sensor. However, since the channel in practice cannot be made evenly by the acrylic material, the nonuniformity of channel could lead to the refraction occurring not only on the surface but also inside the channel, and this phenomenon is difficult to predict in theory. Thus, a device is employed to correct the measurement deviation caused by refraction. It is composed of a micrometer and a container made of acrylic, as Figure7 shows.Figure 7
Measurement calibration device.The thickness of the containerd1 is the same as the thickness of channel d1′ in FIV experiment; meanwhile, the distance from the micrometer to the inner face of the container d2 is equal to that from fuel rod to the inner wall of the channel d2′. Therefore, the measured deviation caused by refraction is the same when the sensor is used to measure the movements of the fuel rods and micrometer, respectively. When the top end of micrometer moves forward to a little distance Δx, a corresponding measured value Δy can be obtained by the laser displacement sensor. The difference between Δx and Δy is the measured deviation caused by laser refraction. Several movements would be made to get two sets of data, Δxn(n=1,2,3,4…) and Δyn(n=1,2,3,4…). The calibration curve for laser displacement sensor could be obtained by a linear fitting with the two sets of data, as Figure 8 shows. The equation of the curve is written as
(5)y=0.72582x-17.90196,
where x is the actual position of the micrometer and y is the measured position.Figure 8
Measurement calibration curve.Assuming the actual position isx1, the measured position is y1; then (5) becomes
(6)y1=0.72582x1-17.90196.Moreover, if the actual position isx2 and the measured position is y2, (5) becomes
(7)y2=0.72582x2-17.90196.Subtracting (7) to (6) gives
(8)y2-y1=0.72582(x2-x1).Noting(x2-x1) is the actual displacement of micrometer Δx and (y2-y1) is the corresponding measured displacement Δy, so (8) becomes
(9)Δy=0.72582·Δx.From (9), it is clear that when the micrometer moves forward 1 μm, the corresponding measured value by displacement sensor is 0.72582 μm, with a deviation of 0.27418 μm, which is caused by the refraction. For eliminating the deviation, it is necessary to correct the measured value of laser displacement sensor according to (9).


## 3.3. Structure Improvement

Generally, the vibration behavior of fuel rods is random in FIV experiment, including forward motion and lateral motion. Since the fuel rods are centrosymmetric structure, statistical characteristics of the vibration in forward direction and lateral direction are the same. Thus, it is only needed to measure the forward motion during the experiment to know the vibration characteristic of fuel rods. However, when a fuel rod vibrates in both directions, as shown in Figure9(a), the position of the laser spot projected by the sensor will change with the vibration of fuel rod, moving along the surface of cylindrical rod and causing the measured value d1 to be larger than the actual displacement d2. That has a substantial influence on the result. In order to solve this problem, the structure of fuel rod has been improved by fixing a tiny metal sheet on the surface of fuel rod. The size of metal sheet is shown in Figure 9(c), and one side of the sheet is flat. In this way, even though the position of laser spot still changes with the movement of fuel rod, the laser spot could move on a flat surface instead of curved surface, as shown in Figure 9(b). As a result, in this situation the measured value d1′ is the same as the actual displacement d2′.Schematic of improved structure.
(a)
Original structured1>d2
(b)
Improved structured1′=d2′
(c)
Size of metal sheet


## 3.4. Measuring System

The hardware of the measuring system mainly consists of a noncontact and high-speed laser displacement sensor (LTC-025-04-SA, MTI), a high-speed data acquisition card (Daq3000, Iotech), and data acquisition soft system (DASYLab, National Instruments). The measuring system is shown in Figure10. In order to reduce the sampling signal bias and frequency alias, a low-pass filter is used to process the detected signal for preventing the high-frequency signal. In addition, a correction module is added to correct the measurement deviation caused by refraction according to (9).Figure 10
Schematic diagram of data acquisition system.


## 4. The Experimental Results and Analysis

### 4.1. Reliability Analysis

In practice, for reducing the difficulty of FIV experiment, the experiment is always carried out at the low temperature and pressure conditions, such as FIV experiment of heat exchange tube in steam generator. In the same way, the FIV experiment of small scale fuel assembly in supercritical water loop is also carried out at the low temperature and pressure conditions. Even though the natural frequencies of fuel rods in FIV experiment are almost the same with those in the supercritical water loop, the vibration behaviors are different because of the difference between the supercritical fluid and subcritical fluid in two cases. Thus, a reliability analysis is needed to make sure the experimental result is reliable. Because the bottom ends of rods are fully constrained, the fuel rods have no displacement in the axial direction even if there is a density gradient along the rod induced by heating. Therefore, it is only needed to analyze the difference of radial displacement in two cases.Actually, the average value of fluid exciting forceF^, which could cause radial displacement of fuel rod, is proportional to the density of fluid ρ and the square of velocity V (10). Moreover, when the size of channel is fixed, F^ is proportional to the square of flow rate Q (11). Consider
(10)F^∝ρ·V2(11)F^∝ρ·(Qs)2.Since the fluid used in the experiment is normal temperature water (25°C), whose density is larger than that of supercritical water, according to (10), it can be noted that the average value of fluid exciting force F^ in the experiment is larger than that in the supercritical water loop under the same flow rate. In consequence, the average value of axial displacement in the experiment is much larger, meaning the experimental results are conservative. Therefore, if the fuel rod can meet the mechanical performance requirements in the FIV experiment, it can definitely meet the requirements of safety and reliability in the supercritical water loop.


### 4.2. Data Validation

For the vibration induced by turbulence and leakage flow which is random [3], an infinite record should theoretically be used to describe this process. However, only a finite length record and data can be obtained in the experiment. So as to use the characteristics of a finite length time history to characterize the whole process of random vibration, this vibration has to be ergodic [10]. So before making further analysis of the experimental data, it is necessary to make sure whether the random vibration of fuel rod is ergodic or not. Therefore, the experiment has been repeated several times under the same conditions, and the comparison among the statistical values of these experimental results is made. Table 3 shows the root mean square (rms) displacement of the fuel rod in through hole case. It can easily be found from Table 3 that the rms displacement is almost the same in these experiments under the condition of the same flow rate. That means the vibration of fuel rod is a stationary ergodic process.Table 3
Root mean square displacement of the fuel rod in through hole case.
Flow rate (m3/h)
Root mean square displacement (μm)
1st
2nd
3th
Average value
0.63
10.71
11.33
11.08
11.04
1.39
13.17
13.62
13.74
13.51
2.12
15.08
14.97
15.34
15.13
2.88
17.04
17.22
16.86
17.03
3.62
19.69
19.13
19.42
19.41


### 4.3. Spectral Analysis

Since there is no Fourier transform for a random signal, its spectral characteristics are always characterized by the power spectral density (PSD), which is the measurement for the distribution of vibrational energy in frequency domain and contains the information about amplitude of vibration in time domain [11]. Figures 11 and 12 separately show the displacement PSD of the fuel rod in two different cases with different flow rate. It is shown in Figure 11 that, with the increase in flow rate, there are obvious changes in power spectrum curves at the frequency range of 6–20 Hz. According to the meaning of PSD mentioned above, it is certain that with the increase of flow rate, the proportion of vibrational energy distribution, and the amplitude of vibration increase at the frequency range of 6–20 Hz. In comparison, it can be easily concluded from Figure 12 that there is no significant change in power spectrum curve under the same conditions, except that the magnitude of spectrum tends to become larger. The reason for that is that the increase of flow velocity will lead to the enhancement of turbulent-boundary-layer pressure fluctuations, which is the most important near-field noise to induce the vibration of fuel rods [12].Figure 11
Displacement power spectral density of the fuel rod in through hole case.Figure 12
Displacement power spectral density of the fuel rod in blind hole case.As shown in Figure13, both the shape of curves and the magnitude of spectrums in two cases are essentially the same in low-frequency bands below 5 Hz at qv=3.62m3/h, but the trends of two curves become much different after 5 Hz. The curve declines gradually with frequency in blind hole case. On the contrary, in the through hole case, the curve shows an ascending trend in the frequency ranging from 6 Hz to 13 Hz and reaches the peak value at 13 Hz, and then the curve decreases after 13 Hz. As a consequence, it can be concluded that the differences in PSD between these two cases are mainly caused by leakage flow occurring on the top of fuel rods. Furthermore, the influence of leakage flow includes two aspects: one is the change in the distribution of vibrational energy in spectrum, and the other is the increase in the magnitude of spectrum at the frequency range of 6–20 Hz.Figure 13
Displacement power spectral density of the fuel rod in two cases atqv=3.62m3/h.Theoretically, since (1) is applicable for the studies of leakage-flow-induced vibration of the fuel rod, it can be used to calculate the vibration of fuel rod caused by leakage flow only if associated with appropriate boundary conditions. On the other hand, if the boundary conditions are not set properly, the computational result of (1) could not show the influence of leakage flow mentioned above. Therefore, the experimental results provide a benchmark for the theoretical calculation, and comparing the experimental result and numerical result could determine whether the boundary condition is proper for the studies related to leakage-flow-induced vibration.It is shown in Figure14 that the power spectrum curves in two cases essentially have the same changing tendency at qv=0.62m3/h, and the reason is that the influence of leakage flow weakens as the flow rate decreases. Compared with turbulent-boundary-layer pressure fluctuations, leakage flow has little influence on the vibration of fuel rods when the flow rate reduces to a certain value.Figure 14
Displacement power spectral density of the fuel rod in two cases atqv=0.62m3/h.


### 4.4. Statistical Analysis

After making a statistical analysis on the vibration behavior of fuel rods atqv=3.62m3/h, the probability density distributions of displacements at midpoint in both cases are obtained; they all could be fitted well by a normal distribution function (see Figures 15 and 16). From the comparison, it can be found that the distribution of vibration displacements in blind hole case, in contrast to the one in through hole case, is more concentrated at 0 μm. On the basis of probability statistics theory, it can be determined that the variance of random vibration, which represents the dynamic component of vibration amplitude, is smaller in blind hole case than that in through hole case. The reason is the existence of leakage flow in through hole case as well.Figure 15
Probability density distributions of displacements in through hole case atqv=3.62m3/h.Figure 16
Probability density distributions of displacements in blind hole case atqv=3.62m3/h.


### 4.5. Analysis of the Root Mean Square Displacement

The absolute value of amplitude at certain time is meaningless because the vibration of fuel rod is random. Therefore, we calculate the root mean square (rms) displacement, which is a statistical measurement of the magnitude of the varying displacement and contains the information about vibration energy. The rms displacements at midpoint changing with flow rate in two cases are compared in Figure17. It is found that the uptrend of rms displacement in through hole case is more obvious than that in blind hole case. In through hole case, if the rod has a small rightward velocity V caused by turbulence, it will lead to the reduction of right side clearance (see Figure 18). As a result, the flow across the right side clearance will decrease with time. The opposite situation exists with respect to the flow in the left side clearance. Because of flow velocity redistribution in the clearance, there are changes in the static pressure on both sides of clearance consequently. It gives rise to a net fluid force acting on the top end of fuel rod in the same direction as the motion. Hence, this results in amplified motion. Under the same other conditions, the difference between rms displacements in two different cases is strongly dependent on leakage flow. Therefore, in order to reduce the flow-induced vibration, which may lead to fatigue of the fuel rods, it is more reasonable to use the spacers with blind hole to assemble the small scale fuel rods in the supercritical water loop.Figure 17
Root-mean-square displacement changing with flow rate in two cases.Figure 18
Diagram to explain the mechanism associated with leakage-flow-induced instability.


## 4.1. Reliability Analysis

In practice, for reducing the difficulty of FIV experiment, the experiment is always carried out at the low temperature and pressure conditions, such as FIV experiment of heat exchange tube in steam generator. In the same way, the FIV experiment of small scale fuel assembly in supercritical water loop is also carried out at the low temperature and pressure conditions. Even though the natural frequencies of fuel rods in FIV experiment are almost the same with those in the supercritical water loop, the vibration behaviors are different because of the difference between the supercritical fluid and subcritical fluid in two cases. Thus, a reliability analysis is needed to make sure the experimental result is reliable. Because the bottom ends of rods are fully constrained, the fuel rods have no displacement in the axial direction even if there is a density gradient along the rod induced by heating. Therefore, it is only needed to analyze the difference of radial displacement in two cases.Actually, the average value of fluid exciting forceF^, which could cause radial displacement of fuel rod, is proportional to the density of fluid ρ and the square of velocity V (10). Moreover, when the size of channel is fixed, F^ is proportional to the square of flow rate Q (11). Consider
(10)F^∝ρ·V2(11)F^∝ρ·(Qs)2.Since the fluid used in the experiment is normal temperature water (25°C), whose density is larger than that of supercritical water, according to (10), it can be noted that the average value of fluid exciting force F^ in the experiment is larger than that in the supercritical water loop under the same flow rate. In consequence, the average value of axial displacement in the experiment is much larger, meaning the experimental results are conservative. Therefore, if the fuel rod can meet the mechanical performance requirements in the FIV experiment, it can definitely meet the requirements of safety and reliability in the supercritical water loop.


## 4.2. Data Validation

For the vibration induced by turbulence and leakage flow which is random [3], an infinite record should theoretically be used to describe this process. However, only a finite length record and data can be obtained in the experiment. So as to use the characteristics of a finite length time history to characterize the whole process of random vibration, this vibration has to be ergodic [10]. So before making further analysis of the experimental data, it is necessary to make sure whether the random vibration of fuel rod is ergodic or not. Therefore, the experiment has been repeated several times under the same conditions, and the comparison among the statistical values of these experimental results is made. Table 3 shows the root mean square (rms) displacement of the fuel rod in through hole case. It can easily be found from Table 3 that the rms displacement is almost the same in these experiments under the condition of the same flow rate. That means the vibration of fuel rod is a stationary ergodic process.Table 3
Root mean square displacement of the fuel rod in through hole case.
Flow rate (m3/h)
Root mean square displacement (μm)
1st
2nd
3th
Average value
0.63
10.71
11.33
11.08
11.04
1.39
13.17
13.62
13.74
13.51
2.12
15.08
14.97
15.34
15.13
2.88
17.04
17.22
16.86
17.03
3.62
19.69
19.13
19.42
19.41


## 4.3. Spectral Analysis

Since there is no Fourier transform for a random signal, its spectral characteristics are always characterized by the power spectral density (PSD), which is the measurement for the distribution of vibrational energy in frequency domain and contains the information about amplitude of vibration in time domain [11]. Figures 11 and 12 separately show the displacement PSD of the fuel rod in two different cases with different flow rate. It is shown in Figure 11 that, with the increase in flow rate, there are obvious changes in power spectrum curves at the frequency range of 6–20 Hz. According to the meaning of PSD mentioned above, it is certain that with the increase of flow rate, the proportion of vibrational energy distribution, and the amplitude of vibration increase at the frequency range of 6–20 Hz. In comparison, it can be easily concluded from Figure 12 that there is no significant change in power spectrum curve under the same conditions, except that the magnitude of spectrum tends to become larger. The reason for that is that the increase of flow velocity will lead to the enhancement of turbulent-boundary-layer pressure fluctuations, which is the most important near-field noise to induce the vibration of fuel rods [12].Figure 11
Displacement power spectral density of the fuel rod in through hole case.Figure 12
Displacement power spectral density of the fuel rod in blind hole case.As shown in Figure13, both the shape of curves and the magnitude of spectrums in two cases are essentially the same in low-frequency bands below 5 Hz at qv=3.62m3/h, but the trends of two curves become much different after 5 Hz. The curve declines gradually with frequency in blind hole case. On the contrary, in the through hole case, the curve shows an ascending trend in the frequency ranging from 6 Hz to 13 Hz and reaches the peak value at 13 Hz, and then the curve decreases after 13 Hz. As a consequence, it can be concluded that the differences in PSD between these two cases are mainly caused by leakage flow occurring on the top of fuel rods. Furthermore, the influence of leakage flow includes two aspects: one is the change in the distribution of vibrational energy in spectrum, and the other is the increase in the magnitude of spectrum at the frequency range of 6–20 Hz.Figure 13
Displacement power spectral density of the fuel rod in two cases atqv=3.62m3/h.Theoretically, since (1) is applicable for the studies of leakage-flow-induced vibration of the fuel rod, it can be used to calculate the vibration of fuel rod caused by leakage flow only if associated with appropriate boundary conditions. On the other hand, if the boundary conditions are not set properly, the computational result of (1) could not show the influence of leakage flow mentioned above. Therefore, the experimental results provide a benchmark for the theoretical calculation, and comparing the experimental result and numerical result could determine whether the boundary condition is proper for the studies related to leakage-flow-induced vibration.It is shown in Figure14 that the power spectrum curves in two cases essentially have the same changing tendency at qv=0.62m3/h, and the reason is that the influence of leakage flow weakens as the flow rate decreases. Compared with turbulent-boundary-layer pressure fluctuations, leakage flow has little influence on the vibration of fuel rods when the flow rate reduces to a certain value.Figure 14
Displacement power spectral density of the fuel rod in two cases atqv=0.62m3/h.


## 4.4. Statistical Analysis

After making a statistical analysis on the vibration behavior of fuel rods atqv=3.62m3/h, the probability density distributions of displacements at midpoint in both cases are obtained; they all could be fitted well by a normal distribution function (see Figures 15 and 16). From the comparison, it can be found that the distribution of vibration displacements in blind hole case, in contrast to the one in through hole case, is more concentrated at 0 μm. On the basis of probability statistics theory, it can be determined that the variance of random vibration, which represents the dynamic component of vibration amplitude, is smaller in blind hole case than that in through hole case. The reason is the existence of leakage flow in through hole case as well.Figure 15
Probability density distributions of displacements in through hole case atqv=3.62m3/h.Figure 16
Probability density distributions of displacements in blind hole case atqv=3.62m3/h.


## 4.5. Analysis of the Root Mean Square Displacement

The absolute value of amplitude at certain time is meaningless because the vibration of fuel rod is random. Therefore, we calculate the root mean square (rms) displacement, which is a statistical measurement of the magnitude of the varying displacement and contains the information about vibration energy. The rms displacements at midpoint changing with flow rate in two cases are compared in Figure17. It is found that the uptrend of rms displacement in through hole case is more obvious than that in blind hole case. In through hole case, if the rod has a small rightward velocity V caused by turbulence, it will lead to the reduction of right side clearance (see Figure 18). As a result, the flow across the right side clearance will decrease with time. The opposite situation exists with respect to the flow in the left side clearance. Because of flow velocity redistribution in the clearance, there are changes in the static pressure on both sides of clearance consequently. It gives rise to a net fluid force acting on the top end of fuel rod in the same direction as the motion. Hence, this results in amplified motion. Under the same other conditions, the difference between rms displacements in two different cases is strongly dependent on leakage flow. Therefore, in order to reduce the flow-induced vibration, which may lead to fatigue of the fuel rods, it is more reasonable to use the spacers with blind hole to assemble the small scale fuel rods in the supercritical water loop.Figure 17
Root-mean-square displacement changing with flow rate in two cases.Figure 18
Diagram to explain the mechanism associated with leakage-flow-induced instability.


## 5. Conclusions

An experiment is made to investigate flow-induced vibration of rods exposed to an axial flow within a rigid square channel, especially the vibration caused by leakage flow. In order to improve the experimental accuracy, measurement calibration is performed to correct the measurement deviation caused by refraction before the FIV experiment. According to the experimental data and analytical results, three useful conclusions are obtained as follows.(1)
When flow rate is low, the vibration of rods is mainly caused by turbulence flow, and leakage flow has little influence on that.(2)
As flow rate increases, the influence of leakage flow becomes obvious, which could change the distribution of vibrational energy in spectrum, increasing the vibrational energy in high-frequency band. Furthermore, it could be treated as a benchmark for the theoretical studies to validate the applicability of boundary condition set for leakage flow. If a boundary condition is set properly, the result of theoretical calculation would show the influence of leakage flow found in the experiment; otherwise this boundary condition is not set properly.(3)
Under the same other conditions, the difference between rms displacements in two different cases is strongly dependent on leakage flow. Therefore, from the point of view of structure safety and reliability, it is more reasonable to improve the design by using the spacers with blind hole to assemble the small scale fuel rods in the supercritical water loop.



---

*Source: 769432-2014-02-24.xml*