# Overview and Discussion of the OECD/NRC Benchmark Based on NUPEC PWR Subchannel and Bundle Tests

Authors: M. Avramova; A. Rubin; H. Utsuno
Journal: Science and Technology of Nuclear Installations
 (2013)

---

## Abstract

The Pennsylvania State University (PSU) under the sponsorship of the US Nuclear Regulatory Commission (NRC) has prepared, organized, conducted, and summarized the Organisation for Economic Co-operation and Development/US Nuclear Regulatory Commission (OECD/NRC) benchmark based on the Nuclear Power Engineering Corporation (NUPEC) pressurized water reactor (PWR) subchannel and bundle tests (PSBTs). The international benchmark activities have been conducted in cooperation with the Nuclear Energy Agency (NEA) of OECD and the Japan Nuclear Energy Safety Organization (JNES), Japan. The OECD/NRC PSBT benchmark was organized to provide a test bed for assessing the capabilities of various thermal-hydraulic subchannel, system, and computational fluid dynamics (CFDs) codes. The benchmark was designed to systematically assess and compare the participants’ numerical models for prediction of detailed subchannel void distribution and department from nucleate boiling (DNB), under steady-state and transient conditions, to full-scale experimental data. This paper provides an overview of the objectives of the benchmark along with a definition of the benchmark phases and exercises. The NUPEC PWR PSBT facility and the specific methods used in the void distribution measurements are discussed followed by a summary of comparative analyses of submitted final results for the exercises of the two benchmark phases.


---

## Body

## 1. Introduction

The need to refine the models for best-estimate calculations, based on good-quality experimental data, has been expressed in many recent meetings in the field of nuclear engineering. The needs arising in this respect should not be limited to the currently available macroscopic methods but should be extended to next-generation analysis techniques that focus on more microscopic processes. One of the most valuable databases identified for the thermal-hydraulics modeling was developed by the Nuclear Power Engineering Corporation (NUPEC), Japan, which includes subchannel void fraction and departure from nucleate boiling (DNB) measurements in a representative pressurized water reactor (PWR) fuel assembly [1–3]. From 1987 to 1995, NUPEC performed a series of void measurement tests using full-size mock-up tests for both BWRs and PWRs. Void fraction measurements and departure from nucleate boiling (DNB) tests were performed at NUPEC under steady-state and transient pressurized water reactor (PWR) conditions. Part of this database is made available for an international benchmark activity entitled as the OECD/NRC NUPEC PWR subchannel and bundle tests (PSBT) benchmark [4]. This international project is officially approved by the Japan Ministry of Economy, Trade, and Industry (METI) and US Nuclear Regulatory Commission (NRC) and endorsed by the OECD/NEA. The benchmark team is organized based on the collaboration between Japan and USA.The fine-mesh high-quality subchannel void fraction and departure from nucleate boiling data encourages the advancement in understanding and modeling complex flow behavior in real rod bundles. Considering that the present theoretical approach is relatively immature, the benchmark specification is designed to systematically assess and compare the participants’ analytical models on the prediction of detailed void distributions and DNB. It should be recognized that the purpose of this benchmark is not only the comparison of currently available computational approaches but, above all, the encouragement to develop novel next-generation approaches that focus on more microscopic processes.The benchmark consists of two phases with a total of seven exercises, as described below.(i)
Phase I: void distribution benchmark,
(a)
exercise 1: steady-state single subchannel benchmark,(b)
exercise 2: steady-state bundle benchmark,(c)
exercise 3: transient bundle benchmark,(d)
exercise 4: Pressure drop benchmark.(ii)
Phase II: DNB benchmark,
(a)
exercise 1: steady-state fluid temperature benchmark,(b)
exercise 2: steady-state DNB benchmark,(c)
exercise 3: transient DNB benchmark.The exercises in phase I of the benchmark are designed to test the codes’ ability to predict void distribution in a single subchannel and a bundle under both steady-state and transient conditions as well as to calculate the pressure drop across a bundle.The exercises in Phase II of the benchmark are designed to test the codes’ ability to predict DNB in a bundle assembly under both steady-state and transient conditions, as well as to predict fluid temperatures under these conditions.


## 2. Description of PSBT Benchmark

The benchmark exercises were performed at the NUPEC test facility, shown in Figure1. The facility consists of a high-pressure and high-temperature recirculation loop, a cooling loop, and instrumentation and data recording systems. The recirculation loop consists of a test section, circulation pump, preheater, steam drum (acting as a pressurizer), and a water mixer. The design pressure is 19.2 MPa and the design temperature is 362°C.Figure 1
System diagram of NUPEC PWR test facility.Three different test sections were used to perform the benchmark exercises: one for the subchannel void distribution test, one for the rod bundle void distribution test, and one for the bundle DNB measurements.The transmission method of gamma-ray was used to measure the density and converted to the void fraction of the gas-liquid two-phase flow. In the DNB measurements, thermocouples were used to determine the heat flux axially along the rod bundle. The thermocouples were attached to the inner surface of the heater rods and determine the boiling transition. The bundle power was increased gradually by fine steps (about 30 kW) to the vicinity of DNB power, which was based on preliminary analysis and experience. The occurrence of DNB was confirmed by a rod temperature rise of more than 11°C (20°F) as measured by the thermocouples.The range of operating conditions for the facility is given in Table1, and the operating conditions for the four transient scenarios are given in Table 2.Table 1
Range of NUPEC PWR test facility operating conditions.
Quantity
Range
Pressure
4.9–16.6 MPa
Mass velocity
550–4150 kg/m2s
Inlet coolant temperature
140–345°CTable 2
Transient parameters of NUPEC PWR test facility.
Transient scenario
Transient change
Depressurization
−0.03 MPa/s
Temperature increase
1°C/s
Flow reduction
−25%/s
Power increase
15%/s


### 2.1. Subchannel Test Assembly

Four different test assemblies were used to model void distribution in a single subchannel. The subchannels represented by these assemblies are described in Table3. Figure 2 shows the test section used for the central subchannel void measurement. It simulates one of the subchannel types found in a PWR assembly. Coolant flows in the pressure vessel horizontally through the coolant inlet nozzle located just below the heated section. Similar test sections (Figure 3) were used for the central with thimble, side, and corner subchannel types. The effective heated length is 1555 mm, and the void measurement section begins at 1400 mm from the bottom of the heated section [4].Table 3
Geometry and power shape for subchannel test assemblies for void distribution benchmark.
Item
Data
Assembly (subjected subchannel)
Subchannel type
Center (typical)
Center (thimble)
Side
Corner
Number of heaters
4 × 1/4
3 × 1/4
2 × 1/4
1 × 1/4
Axial heated length (mm)
1555
1555
1555
1555
Axial power shape
Uniform
Uniform
Uniform
Uniform
Black part: subjected subchannel, white circles: heated rod, and circles with dots: thimble rod.Figure 2
Test section for central subchannel void distribution measurement.Figure 3
Cross-sectional view of subchannel test assembly.


### 2.2. Bundle Test Assembly

Three different test assemblies were used to model void distribution in a rod bundle. The bundles represented by these assemblies are described in Table4. An electrically heated rod bundle was used to simulate a partial section and full length of a PWR fuel assembly. Figure 4 shows the test section used for the rod bundle void measurements. The effective heated length is 3658 mm. The measurements were performed at three axial elevations: upper—3177 mm, middle—2669 mm, and lower—2216 mm, respectively (please note that the measurement elevations are taken from the bottom of heating section). Coolant flows into the pressure vessel horizontally through the coolant inlet nozzle and down through the section between the flow channel and the pressure vessel. The coolant continues into the flow channel, flowing from the bottom of the pressure vessel up through the test assembly, where the bottom of the heated section is located 630 mm above the bottom of the pressure vessel.Table 4
Geometry and power shape for bundle test assemblies for void distribution benchmark.
Item
Data
Assembly
Rods array
5 × 5
5 × 5
5 × 5
Number of heated rods
25
25
24
Number of thimble rods
0
0
1
Heated rod outer diameter (mm)
9.50
9.50
9.50
Thimble rod outer diameter (mm)
—
—
12.24
Heated rods pitch (mm)
12.60
12.60
12.60
Axial heated length (mm)
3658
3658
3658
Flow channel inner width (mm)
64.9
64.9
64.9
Radial power shape
A
A
B
Axial power shape
Uniform
Cosine
Cosine
Number of MV spacers
7
7
7
Number of NMV spacers
2
2
2
Number of simple spacers
8
8
8
MV spacer location (mm)
471, 925, 1378, 1832, 2285, 2739, 3247
NMV spacer location (mm)
2.5, 3755
Simple spacer location (mm)
237, 698, 1151, 1605, 2059, 2512, 2993, 3501
White circles: heated rod and circles with dots: thimble rod.
MV: mixing vane and NMV: non-mixing vane.
Spacer location is the distance from the bottom of heated length to spacer bottom face.Figure 4
Test section for rod bundle void distribution measurement.Figures5 and 6 show the two different radial power distributions, named A and B, respectively. All powers shown are relative powers.Figure 5
Radial power distribution A.Figure 6
Radial power distribution B.Table5 shows the axial power distribution (based on a cosine shape) that was used in the bundle tests.Table 5
Cosine axial power distribution.
Node
Relative power
(Bottom)
1
0.42
2
0.47
3
0.56
4
0.67
5
0.80
6
0.94
7
1.08
8
1.22
9
1.34
10
1.44
11
1.51
12
1.55
13
1.55
14
1.51
15
1.44
16
1.34
17
1.22
18
1.08
19
0.94
20
0.80
21
0.67
22
0.56
23
0.47
24
0.42
(Top)The test assemblies used for the DNB measurement benchmark are described in Table6. The DNB measurements were performed for full-length partial 5×5 and 6×6 array rod bundles, which simulate 17×17 PWR fuel assemblies. Measurements were performed for both steady state and transients. The heater rods used in these bundles are of the same type as those used in the bundles in the void distribution measurements. Three different types of spacer grids were used in these assemblies: simple spacers, nonmixing vane grids, and mixing vane grids [4].Table 6
Test assemblies for DNB measurements.
Assembly
Reference fuel type
Rods array
Type of cell
Power distribution
Radial
Axial
A0
17 × 17 M
Typical cell
A
Uniform
A1
5 × 5
Typical cell
C
Uniform
A2
Typical cell
A
Uniform
A3
6 × 6
Typical cell
D
Uniform
A4
Typical cell
A
Cosine
A8
5 × 5
Thimble cell
B
Cosine
A11
Typical cell
A
Cosine
A12
Thimble cell
B
CosineThe available experimental data for the steady-state single subchannel benchmark includes CT scanner measurements of the void fraction (subchannel averaged) of four representative subchannel types: typical central, central with a guide tube, side, and corner subchannels and images of the void distribution within two representative subchannel types: typical central and central with guide tube. The data can be used to assess and improve the current models for void generation (subchannel/system and CFD codes) and void distribution within subchannels (CFD codes).The steady-state bundle benchmark utilizes experimental data that include X-ray densitometer measurements of void fraction (chordal averaged) at three axial elevations. The averaging is over the four central subchannels. The transient bundle benchmark is based on data collected for four anticipated transient scenarios: power increase, flow reduction, depressurization, and temperature increase.The experimental data for the steady-state DNB benchmark include the power at which departure from nucleate boiling occurs and the corresponding location in the bundle. Measured DNB power for four anticipated transient scenarios—power increase, flow reduction, depressurization, and temperature increase—was made available for the transient DNB benchmark. In addition to the void distribution and DNB power measurements, data is available for the subchannel exit fluid temperature.


### 2.3. Spacer Grid Data

There were three types of spacers instrumented along the axial length: simple spacer (SS), spacer with non-mixing vanes (NMVs), and spacer with mixing vanes (MVs). The simple spacer has only dimples while NMV and MV have dimples and springs. The grids straps are made out of Inconel 600 alloy. Detailed geometrical description of the grids used in the experiment was not available to the benchmark. As a result, the benchmark team, with the assistance of a benchmark participant, was forced to develop a grid model based on the understanding that the grids used in the experiments were similar to grids for which data was readily available. Table7 summarizes the grid data that was available as part of the benchmark. Figures 7, 8, and 9 provide three-dimensional views of the simple spacer, nonmixing vane, and mixing vane grids. The simple spacer has dimples while the mixing vane and nonmixing vane grids have both dimples and springs. These dimples provide a gap (~0.1 mm) around each heating rod, which prevents bowing of these rods when they linearly expand at high temperatures.Table 7
Bundle average spacer pressure loss coefficients.
Spacer type
Loss coefficient
Simple spacer (SS)
0.4
Nonmixing vanes (NMVs) spacer
0.7
Mixing vanes (MVs) spacer
1.0Figure 7
View of simple spacer grid.Figure 8
View of nonmixing vane spacer grid.Figure 9
View of mixing vane spacer grid.Table7 provides the bundle average spacer pressure loss coefficients for the three types of grids. Depending on the participants’ computer code and using the provided spacer data, each participant may choose the subchannel grids loss coefficients or other required input values.


### 2.4. Heater Rod Data

Table8 summarizes the material and geometrical properties of the heater rods used in the rod bundle tests. Figure 10 provides a cross-sectional view of the heater rods and gives dimensions.Table 8
Properties of heater rods.
Item
Data
Heater
Outer diameter (mm)
9.5
Thickness (mm)
0.65
Material
Inconel 600
Heating method
Direct heating
Insulator
Outer diameter (mm)
8.2
Inner diameter (mm)
5.8
Material
AluminaFigure 10
Cross-sectional View of Heater Rod.


### 2.5. Measurement Techniques

#### 2.5.1. Void Distribution

A gamma-ray transmission method was used to measure the density of the flow, which was converted to the void fraction of the gas-liquid two-phase flow. Figure11 shows the procedure used to perform the void fraction measurements for the entire rod bundle. The top half of the figure shows the procedure used to perform the void fraction measurements for a single subchannel. In the subchannel experiments, a narrow gamma-ray beam CT scanner was used to measure the subchannel averaged void fraction, and a wide gamma-ray beam was used to measure the chordal averaged void fraction. For each subchannel type—corner, side, or center—a relationship between the subchannel averaged and the chordal averaged void fractions was individually derived. These relationships were then used to correct the subchannel averaged void fraction measured with the wide beam in the bundle tests.Figure 11
Void fraction measurement procedure.The void measurement systems shown in Figure11 consist of gamma-ray sources (Cs137), detectors, collimators, and signal processing units. The attenuation of the gamma-rays, which depends on the void fraction, was detected. The intensity of the gamma-ray source was determined to obtain the count rate of the signal processing (30×104 cps). The CT scanner system was used to determine the distribution of density/void fraction over the subchannel at steady-state flow and to define the subchannel averaged void fraction. The system was operated by translate/rotate method. At each translation/rotation location, the intensity of gamma-ray attenuated by the object, the so-called “projection data”, was detected. An image reconstruction was then performed by a filtered back-projection algorithm to obtain the distribution of the linear attenuation coefficient. A sufficient measuring time was given in order to avoid the effect of the flow motion. Two densitometer systems—in x- and y-directions—were used in the chordal averaged void fraction measurements in the single subchannel tests. Each of them consisted of a gamma source and a detector. They were located at the same tables (the same elevation) as CT, which was fixed during the measurements. A multibeam system was used to measure each subchannel void fraction of the rod bundle. Six transmission data of x-direction and six transmission data of y-direction between the rod and rod/channel wall were used to reconstruct the void fraction of the 36 subchannels by an iterative method. These subchannel void fractions corresponded to the chordal measurements of the single subchannel tests. The relationships between the subchannel averaged void fraction and the chordal averaged void fraction obtained in the single subchannel tests were used to determine the subchannel averaged void fractions in the rod bundle tests. Such measurements were performed simultaneously at three axial elevations. Table 9 shows the accuracy of the various parameters involved in the experiment. Table 10 summarizes the sources of error in the experiment.Table 9
Accuracy of process parameters in void distribution measurement.
Quantity
Accuracy
Process parameters
Pressure
1%
Flow
1.5%
Power
1%
Fluid temperature
1 Celsius
Void fraction measurement
CT measurement
Gamma-ray beam width
1 mm
Subchannel averaged (steady state)
3% void
Spatial resolution of one pixel
0.5 mm
Chordal measurement
Gamma-ray beam width (center)
3 mm
Gamma-ray beam width (side)
2 mm
Subchannel averaged (steady state)
4% void
Subchannel averaged (transient)
5% voidTable 10
Sources of error for void distribution measurement.
Error source
Chordal averaged
CT averaged
Steady state
Transient
γ-ray measurement
Effect of surrounding condition (magnetic field and temperature) on measurement system
0.1%
0.1%
0.1%
Randomness ofγ-ray source decay
0.02%
0.2%
0.1%
Correction error due to background
0.0%
0.0%
0.0%
Correction error due to counting loss
<0.5%
<0.5%
<0.1%
Calibration error
0.1%
0.1%
0.1%
Correction error due to attenuation  by surrounding water
0.0%
0.0%
—
Correction error due to scattering from multi-γ-rays
<0.2%
<0.2%
—
Total
<0.55%
<0.6%
<0.2%
Subchannel density
Transfer to density
<9 kg/m3
<10 kg/m3
<15 kg/m3
Distribution error to subchannel
<5 kg/m3
<5 kg/m3
—
Correlation error from chordal averaged to CT averaged
<6 kg/m3
<6 kg/m3
—
Subchannel density
<20 kg/m3
<21 kg/m3
<15 kg/m3
Subchannel void*
0.040
0.042
0.030
Uncertainty (1σ)
4%
5%
3%
*Reference averaged density is 500 kg/m3.Figure12 illustrates the relationship between chordal and CT averaged densities as a function of pressure. It was used to determine the uncertainties inherent in the void measurements and to correct the measured values based on the pressure of the test case. The correlation between the chordal and CT averaged values is given by the best-fit curves and was introduced for the high- and low-pressure conditions, respectively. The reference averaged density was 500 kg/m3. The uncertainty of the correlation was determined to be less than 18 kg/m3, which was regarded as three standard deviations (3σ). Therefore, one standard deviation (1σ) is 6 kg/m3. Table 11 shows the number of gamma-ray beams used in the fluid density measurement for both subchannel and rod bundle exercises. Table 12 shows the amount of time required to perform the density measurement.Table 11
Number of gamma-ray beams.
Test assembly
CT measurement
Chordal measurement
Subchannel
2(X and Y direction)
2(X and Y direction)
Rod bundle
—
6 beams × 2 × 3 section(total 36 beams)Table 12
Time required to perform void fraction measurements.
Item
CT measurement
Chordal measurement
Steady state
Time needed
5 s/step ×T33 ×  R17 step(it takes 2 h)
100 ssampling cycle 0.1 s
Measurement
2 times
3 times
Transient
Time needed
—
200 s
Measurement
—
1 timeFigure 12
Relation Between Chordal and CT Averaged Densities (for S1).


#### 2.5.2. Departure from Nucleate Boiling

The bundle power was gradually increased in fine steps to the expected vicinity of DNB, which was based on previous analysis operator experience. The onset of DNB was confirmed by a rod temperature rise greater than 11°C (20°F) as measured by the thermocouples seen in Figure13. The DNB power is defined as the power corresponding to the step immediately preceding the step in which this temperature rise is seen. The accuracy of the process parameters involved in this process is seen in Table 13. The exit fluid temperatures were measured by the thermocouples shown in Figure 14.Table 13
Accuracy of process parameters in DNB measurement.
Quantity
Accuracy
Process parameters
Pressure
1%
Flow
1.5%
Power
1%
Fluid temperature
1 CelsiusFigure 13
Locations of thermocouples for test assemblies.Figure 14
Diagram of fluid temperature measurement setup.


### 2.6. Experimental Data Analysis

Based on the participants feedback, several studies were performed to determine the validity of the benchmark data. These included (1) study performed on recalculation of the void fraction, (2) study performed on calculation of the quality, (3) study performed on the subchannel CT masking, and (4) study performed on bundle test section downcomer region. In this paper, the studies (1), (2), and (4) performed by the benchmark team are presented.


#### 2.6.1. Recalculation of Void Fraction

At the first PSBT workshop, it was noted by several participants that the “measured” void fractions (which were actually calculated from measured density data) were not consistent with void fractions calculated using the measured densities. As a result, a study was performed to recalculate the void fraction and the quality for each test case in the void distribution benchmark.Solving the standard equation for mixture density (ρ-=αρg+(1-α)ρf) for the void fraction gives α=(ρ--ρf)/(ρg-ρf), where the liquid and vapor densities (ρf and ρg, resp.) are assumed to be at saturation and the mixture density is taken from the benchmark data. After the recalculation (see Figure 15), it was noted that the measured void fraction was consistently higher than the recalculated void fraction. This recalculation was only performed for the subchannel test cases since those are the only test cases for which fluid density data was available.Figure 15
Deviation of measured void fraction from recalculated void fraction.


#### 2.6.2. Calculation of Quality

Upon completion of the study performed on void distribution, the benchmark team began a study of the calculation of quality based on the experimentally determined densities.It is recalled that the quality can be expressed using mixture enthalpy asx=(hmix-hf)/(hg-hf), where hf and hg are the liquid and vapor enthalpies, respectively. A number of different expressions were derived to determine the mixture enthalpy in the test sections assuming conservation of energy. After verifying that the axial power distribution was normalized for both the uniform and cosine power shapes, the following equations were obtained.All four subchannel test sections utilized a uniform axial power distribution. Thus, for all subchannel assemblies, the mixture enthalpy at the measurement section can be given by(1)hmix=hin+(1400mm1555mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h]).Figure16 shows the resulting deviation of the experimental quality from the recalculated quality.Figure 16
Deviation of measured quality from recalculated quality for test series 1 and 2.Assembly B5 utilized a uniform axial power distribution, so the mixture enthalpies at the three measurement locations can be given as(2)hmix=hin+(2216mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(lower elevation),hmix=hin+(2669mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(middle elevation),hmix=hin+(3177mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(upperelevation).Assemblies B6 and B7 utilized a cosine axial power shape. Recalling that the power shape is normalized, it is possible to determine what fraction of the total power has been imparted to the fluid between the flow inlet and the measurement sections. The mixture enthalpies for these two assemblies can be given as(3)hmix=hin+0.6598×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(Lower Elevation),hmix=hin+0.8172×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(MiddleElevation),hmix=hin+0.9353×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(UpperElevation).After applying these equations and calculating the qualities for each case, the following results were obtained. Figures17, 18, 19, and 20 show the deviation of experimental quality from recalculated quality for test series 5, 6, 7, and 8 (resp.).Figure 17
Deviation of measured quality from recalculated quality for test series 5.Figure 18
Deviation of measured quality from recalculated quality for test series 6.Figure 19
Deviation of measured quality from recalculated quality for test series 7.Figure 20
Deviation of measured quality from recalculated quality for test series 8.


## 2.1. Subchannel Test Assembly

Four different test assemblies were used to model void distribution in a single subchannel. The subchannels represented by these assemblies are described in Table3. Figure 2 shows the test section used for the central subchannel void measurement. It simulates one of the subchannel types found in a PWR assembly. Coolant flows in the pressure vessel horizontally through the coolant inlet nozzle located just below the heated section. Similar test sections (Figure 3) were used for the central with thimble, side, and corner subchannel types. The effective heated length is 1555 mm, and the void measurement section begins at 1400 mm from the bottom of the heated section [4].Table 3
Geometry and power shape for subchannel test assemblies for void distribution benchmark.
Item
Data
Assembly (subjected subchannel)
Subchannel type
Center (typical)
Center (thimble)
Side
Corner
Number of heaters
4 × 1/4
3 × 1/4
2 × 1/4
1 × 1/4
Axial heated length (mm)
1555
1555
1555
1555
Axial power shape
Uniform
Uniform
Uniform
Uniform
Black part: subjected subchannel, white circles: heated rod, and circles with dots: thimble rod.Figure 2
Test section for central subchannel void distribution measurement.Figure 3
Cross-sectional view of subchannel test assembly.


## 2.2. Bundle Test Assembly

Three different test assemblies were used to model void distribution in a rod bundle. The bundles represented by these assemblies are described in Table4. An electrically heated rod bundle was used to simulate a partial section and full length of a PWR fuel assembly. Figure 4 shows the test section used for the rod bundle void measurements. The effective heated length is 3658 mm. The measurements were performed at three axial elevations: upper—3177 mm, middle—2669 mm, and lower—2216 mm, respectively (please note that the measurement elevations are taken from the bottom of heating section). Coolant flows into the pressure vessel horizontally through the coolant inlet nozzle and down through the section between the flow channel and the pressure vessel. The coolant continues into the flow channel, flowing from the bottom of the pressure vessel up through the test assembly, where the bottom of the heated section is located 630 mm above the bottom of the pressure vessel.Table 4
Geometry and power shape for bundle test assemblies for void distribution benchmark.
Item
Data
Assembly
Rods array
5 × 5
5 × 5
5 × 5
Number of heated rods
25
25
24
Number of thimble rods
0
0
1
Heated rod outer diameter (mm)
9.50
9.50
9.50
Thimble rod outer diameter (mm)
—
—
12.24
Heated rods pitch (mm)
12.60
12.60
12.60
Axial heated length (mm)
3658
3658
3658
Flow channel inner width (mm)
64.9
64.9
64.9
Radial power shape
A
A
B
Axial power shape
Uniform
Cosine
Cosine
Number of MV spacers
7
7
7
Number of NMV spacers
2
2
2
Number of simple spacers
8
8
8
MV spacer location (mm)
471, 925, 1378, 1832, 2285, 2739, 3247
NMV spacer location (mm)
2.5, 3755
Simple spacer location (mm)
237, 698, 1151, 1605, 2059, 2512, 2993, 3501
White circles: heated rod and circles with dots: thimble rod.
MV: mixing vane and NMV: non-mixing vane.
Spacer location is the distance from the bottom of heated length to spacer bottom face.Figure 4
Test section for rod bundle void distribution measurement.Figures5 and 6 show the two different radial power distributions, named A and B, respectively. All powers shown are relative powers.Figure 5
Radial power distribution A.Figure 6
Radial power distribution B.Table5 shows the axial power distribution (based on a cosine shape) that was used in the bundle tests.Table 5
Cosine axial power distribution.
Node
Relative power
(Bottom)
1
0.42
2
0.47
3
0.56
4
0.67
5
0.80
6
0.94
7
1.08
8
1.22
9
1.34
10
1.44
11
1.51
12
1.55
13
1.55
14
1.51
15
1.44
16
1.34
17
1.22
18
1.08
19
0.94
20
0.80
21
0.67
22
0.56
23
0.47
24
0.42
(Top)The test assemblies used for the DNB measurement benchmark are described in Table6. The DNB measurements were performed for full-length partial 5×5 and 6×6 array rod bundles, which simulate 17×17 PWR fuel assemblies. Measurements were performed for both steady state and transients. The heater rods used in these bundles are of the same type as those used in the bundles in the void distribution measurements. Three different types of spacer grids were used in these assemblies: simple spacers, nonmixing vane grids, and mixing vane grids [4].Table 6
Test assemblies for DNB measurements.
Assembly
Reference fuel type
Rods array
Type of cell
Power distribution
Radial
Axial
A0
17 × 17 M
Typical cell
A
Uniform
A1
5 × 5
Typical cell
C
Uniform
A2
Typical cell
A
Uniform
A3
6 × 6
Typical cell
D
Uniform
A4
Typical cell
A
Cosine
A8
5 × 5
Thimble cell
B
Cosine
A11
Typical cell
A
Cosine
A12
Thimble cell
B
CosineThe available experimental data for the steady-state single subchannel benchmark includes CT scanner measurements of the void fraction (subchannel averaged) of four representative subchannel types: typical central, central with a guide tube, side, and corner subchannels and images of the void distribution within two representative subchannel types: typical central and central with guide tube. The data can be used to assess and improve the current models for void generation (subchannel/system and CFD codes) and void distribution within subchannels (CFD codes).The steady-state bundle benchmark utilizes experimental data that include X-ray densitometer measurements of void fraction (chordal averaged) at three axial elevations. The averaging is over the four central subchannels. The transient bundle benchmark is based on data collected for four anticipated transient scenarios: power increase, flow reduction, depressurization, and temperature increase.The experimental data for the steady-state DNB benchmark include the power at which departure from nucleate boiling occurs and the corresponding location in the bundle. Measured DNB power for four anticipated transient scenarios—power increase, flow reduction, depressurization, and temperature increase—was made available for the transient DNB benchmark. In addition to the void distribution and DNB power measurements, data is available for the subchannel exit fluid temperature.


## 2.3. Spacer Grid Data

There were three types of spacers instrumented along the axial length: simple spacer (SS), spacer with non-mixing vanes (NMVs), and spacer with mixing vanes (MVs). The simple spacer has only dimples while NMV and MV have dimples and springs. The grids straps are made out of Inconel 600 alloy. Detailed geometrical description of the grids used in the experiment was not available to the benchmark. As a result, the benchmark team, with the assistance of a benchmark participant, was forced to develop a grid model based on the understanding that the grids used in the experiments were similar to grids for which data was readily available. Table7 summarizes the grid data that was available as part of the benchmark. Figures 7, 8, and 9 provide three-dimensional views of the simple spacer, nonmixing vane, and mixing vane grids. The simple spacer has dimples while the mixing vane and nonmixing vane grids have both dimples and springs. These dimples provide a gap (~0.1 mm) around each heating rod, which prevents bowing of these rods when they linearly expand at high temperatures.Table 7
Bundle average spacer pressure loss coefficients.
Spacer type
Loss coefficient
Simple spacer (SS)
0.4
Nonmixing vanes (NMVs) spacer
0.7
Mixing vanes (MVs) spacer
1.0Figure 7
View of simple spacer grid.Figure 8
View of nonmixing vane spacer grid.Figure 9
View of mixing vane spacer grid.Table7 provides the bundle average spacer pressure loss coefficients for the three types of grids. Depending on the participants’ computer code and using the provided spacer data, each participant may choose the subchannel grids loss coefficients or other required input values.


## 2.4. Heater Rod Data

Table8 summarizes the material and geometrical properties of the heater rods used in the rod bundle tests. Figure 10 provides a cross-sectional view of the heater rods and gives dimensions.Table 8
Properties of heater rods.
Item
Data
Heater
Outer diameter (mm)
9.5
Thickness (mm)
0.65
Material
Inconel 600
Heating method
Direct heating
Insulator
Outer diameter (mm)
8.2
Inner diameter (mm)
5.8
Material
AluminaFigure 10
Cross-sectional View of Heater Rod.


## 2.5. Measurement Techniques

### 2.5.1. Void Distribution

A gamma-ray transmission method was used to measure the density of the flow, which was converted to the void fraction of the gas-liquid two-phase flow. Figure11 shows the procedure used to perform the void fraction measurements for the entire rod bundle. The top half of the figure shows the procedure used to perform the void fraction measurements for a single subchannel. In the subchannel experiments, a narrow gamma-ray beam CT scanner was used to measure the subchannel averaged void fraction, and a wide gamma-ray beam was used to measure the chordal averaged void fraction. For each subchannel type—corner, side, or center—a relationship between the subchannel averaged and the chordal averaged void fractions was individually derived. These relationships were then used to correct the subchannel averaged void fraction measured with the wide beam in the bundle tests.Figure 11
Void fraction measurement procedure.The void measurement systems shown in Figure11 consist of gamma-ray sources (Cs137), detectors, collimators, and signal processing units. The attenuation of the gamma-rays, which depends on the void fraction, was detected. The intensity of the gamma-ray source was determined to obtain the count rate of the signal processing (30×104 cps). The CT scanner system was used to determine the distribution of density/void fraction over the subchannel at steady-state flow and to define the subchannel averaged void fraction. The system was operated by translate/rotate method. At each translation/rotation location, the intensity of gamma-ray attenuated by the object, the so-called “projection data”, was detected. An image reconstruction was then performed by a filtered back-projection algorithm to obtain the distribution of the linear attenuation coefficient. A sufficient measuring time was given in order to avoid the effect of the flow motion. Two densitometer systems—in x- and y-directions—were used in the chordal averaged void fraction measurements in the single subchannel tests. Each of them consisted of a gamma source and a detector. They were located at the same tables (the same elevation) as CT, which was fixed during the measurements. A multibeam system was used to measure each subchannel void fraction of the rod bundle. Six transmission data of x-direction and six transmission data of y-direction between the rod and rod/channel wall were used to reconstruct the void fraction of the 36 subchannels by an iterative method. These subchannel void fractions corresponded to the chordal measurements of the single subchannel tests. The relationships between the subchannel averaged void fraction and the chordal averaged void fraction obtained in the single subchannel tests were used to determine the subchannel averaged void fractions in the rod bundle tests. Such measurements were performed simultaneously at three axial elevations. Table 9 shows the accuracy of the various parameters involved in the experiment. Table 10 summarizes the sources of error in the experiment.Table 9
Accuracy of process parameters in void distribution measurement.
Quantity
Accuracy
Process parameters
Pressure
1%
Flow
1.5%
Power
1%
Fluid temperature
1 Celsius
Void fraction measurement
CT measurement
Gamma-ray beam width
1 mm
Subchannel averaged (steady state)
3% void
Spatial resolution of one pixel
0.5 mm
Chordal measurement
Gamma-ray beam width (center)
3 mm
Gamma-ray beam width (side)
2 mm
Subchannel averaged (steady state)
4% void
Subchannel averaged (transient)
5% voidTable 10
Sources of error for void distribution measurement.
Error source
Chordal averaged
CT averaged
Steady state
Transient
γ-ray measurement
Effect of surrounding condition (magnetic field and temperature) on measurement system
0.1%
0.1%
0.1%
Randomness ofγ-ray source decay
0.02%
0.2%
0.1%
Correction error due to background
0.0%
0.0%
0.0%
Correction error due to counting loss
<0.5%
<0.5%
<0.1%
Calibration error
0.1%
0.1%
0.1%
Correction error due to attenuation  by surrounding water
0.0%
0.0%
—
Correction error due to scattering from multi-γ-rays
<0.2%
<0.2%
—
Total
<0.55%
<0.6%
<0.2%
Subchannel density
Transfer to density
<9 kg/m3
<10 kg/m3
<15 kg/m3
Distribution error to subchannel
<5 kg/m3
<5 kg/m3
—
Correlation error from chordal averaged to CT averaged
<6 kg/m3
<6 kg/m3
—
Subchannel density
<20 kg/m3
<21 kg/m3
<15 kg/m3
Subchannel void*
0.040
0.042
0.030
Uncertainty (1σ)
4%
5%
3%
*Reference averaged density is 500 kg/m3.Figure12 illustrates the relationship between chordal and CT averaged densities as a function of pressure. It was used to determine the uncertainties inherent in the void measurements and to correct the measured values based on the pressure of the test case. The correlation between the chordal and CT averaged values is given by the best-fit curves and was introduced for the high- and low-pressure conditions, respectively. The reference averaged density was 500 kg/m3. The uncertainty of the correlation was determined to be less than 18 kg/m3, which was regarded as three standard deviations (3σ). Therefore, one standard deviation (1σ) is 6 kg/m3. Table 11 shows the number of gamma-ray beams used in the fluid density measurement for both subchannel and rod bundle exercises. Table 12 shows the amount of time required to perform the density measurement.Table 11
Number of gamma-ray beams.
Test assembly
CT measurement
Chordal measurement
Subchannel
2(X and Y direction)
2(X and Y direction)
Rod bundle
—
6 beams × 2 × 3 section(total 36 beams)Table 12
Time required to perform void fraction measurements.
Item
CT measurement
Chordal measurement
Steady state
Time needed
5 s/step ×T33 ×  R17 step(it takes 2 h)
100 ssampling cycle 0.1 s
Measurement
2 times
3 times
Transient
Time needed
—
200 s
Measurement
—
1 timeFigure 12
Relation Between Chordal and CT Averaged Densities (for S1).


### 2.5.2. Departure from Nucleate Boiling

The bundle power was gradually increased in fine steps to the expected vicinity of DNB, which was based on previous analysis operator experience. The onset of DNB was confirmed by a rod temperature rise greater than 11°C (20°F) as measured by the thermocouples seen in Figure13. The DNB power is defined as the power corresponding to the step immediately preceding the step in which this temperature rise is seen. The accuracy of the process parameters involved in this process is seen in Table 13. The exit fluid temperatures were measured by the thermocouples shown in Figure 14.Table 13
Accuracy of process parameters in DNB measurement.
Quantity
Accuracy
Process parameters
Pressure
1%
Flow
1.5%
Power
1%
Fluid temperature
1 CelsiusFigure 13
Locations of thermocouples for test assemblies.Figure 14
Diagram of fluid temperature measurement setup.


## 2.5.1. Void Distribution

A gamma-ray transmission method was used to measure the density of the flow, which was converted to the void fraction of the gas-liquid two-phase flow. Figure11 shows the procedure used to perform the void fraction measurements for the entire rod bundle. The top half of the figure shows the procedure used to perform the void fraction measurements for a single subchannel. In the subchannel experiments, a narrow gamma-ray beam CT scanner was used to measure the subchannel averaged void fraction, and a wide gamma-ray beam was used to measure the chordal averaged void fraction. For each subchannel type—corner, side, or center—a relationship between the subchannel averaged and the chordal averaged void fractions was individually derived. These relationships were then used to correct the subchannel averaged void fraction measured with the wide beam in the bundle tests.Figure 11
Void fraction measurement procedure.The void measurement systems shown in Figure11 consist of gamma-ray sources (Cs137), detectors, collimators, and signal processing units. The attenuation of the gamma-rays, which depends on the void fraction, was detected. The intensity of the gamma-ray source was determined to obtain the count rate of the signal processing (30×104 cps). The CT scanner system was used to determine the distribution of density/void fraction over the subchannel at steady-state flow and to define the subchannel averaged void fraction. The system was operated by translate/rotate method. At each translation/rotation location, the intensity of gamma-ray attenuated by the object, the so-called “projection data”, was detected. An image reconstruction was then performed by a filtered back-projection algorithm to obtain the distribution of the linear attenuation coefficient. A sufficient measuring time was given in order to avoid the effect of the flow motion. Two densitometer systems—in x- and y-directions—were used in the chordal averaged void fraction measurements in the single subchannel tests. Each of them consisted of a gamma source and a detector. They were located at the same tables (the same elevation) as CT, which was fixed during the measurements. A multibeam system was used to measure each subchannel void fraction of the rod bundle. Six transmission data of x-direction and six transmission data of y-direction between the rod and rod/channel wall were used to reconstruct the void fraction of the 36 subchannels by an iterative method. These subchannel void fractions corresponded to the chordal measurements of the single subchannel tests. The relationships between the subchannel averaged void fraction and the chordal averaged void fraction obtained in the single subchannel tests were used to determine the subchannel averaged void fractions in the rod bundle tests. Such measurements were performed simultaneously at three axial elevations. Table 9 shows the accuracy of the various parameters involved in the experiment. Table 10 summarizes the sources of error in the experiment.Table 9
Accuracy of process parameters in void distribution measurement.
Quantity
Accuracy
Process parameters
Pressure
1%
Flow
1.5%
Power
1%
Fluid temperature
1 Celsius
Void fraction measurement
CT measurement
Gamma-ray beam width
1 mm
Subchannel averaged (steady state)
3% void
Spatial resolution of one pixel
0.5 mm
Chordal measurement
Gamma-ray beam width (center)
3 mm
Gamma-ray beam width (side)
2 mm
Subchannel averaged (steady state)
4% void
Subchannel averaged (transient)
5% voidTable 10
Sources of error for void distribution measurement.
Error source
Chordal averaged
CT averaged
Steady state
Transient
γ-ray measurement
Effect of surrounding condition (magnetic field and temperature) on measurement system
0.1%
0.1%
0.1%
Randomness ofγ-ray source decay
0.02%
0.2%
0.1%
Correction error due to background
0.0%
0.0%
0.0%
Correction error due to counting loss
<0.5%
<0.5%
<0.1%
Calibration error
0.1%
0.1%
0.1%
Correction error due to attenuation  by surrounding water
0.0%
0.0%
—
Correction error due to scattering from multi-γ-rays
<0.2%
<0.2%
—
Total
<0.55%
<0.6%
<0.2%
Subchannel density
Transfer to density
<9 kg/m3
<10 kg/m3
<15 kg/m3
Distribution error to subchannel
<5 kg/m3
<5 kg/m3
—
Correlation error from chordal averaged to CT averaged
<6 kg/m3
<6 kg/m3
—
Subchannel density
<20 kg/m3
<21 kg/m3
<15 kg/m3
Subchannel void*
0.040
0.042
0.030
Uncertainty (1σ)
4%
5%
3%
*Reference averaged density is 500 kg/m3.Figure12 illustrates the relationship between chordal and CT averaged densities as a function of pressure. It was used to determine the uncertainties inherent in the void measurements and to correct the measured values based on the pressure of the test case. The correlation between the chordal and CT averaged values is given by the best-fit curves and was introduced for the high- and low-pressure conditions, respectively. The reference averaged density was 500 kg/m3. The uncertainty of the correlation was determined to be less than 18 kg/m3, which was regarded as three standard deviations (3σ). Therefore, one standard deviation (1σ) is 6 kg/m3. Table 11 shows the number of gamma-ray beams used in the fluid density measurement for both subchannel and rod bundle exercises. Table 12 shows the amount of time required to perform the density measurement.Table 11
Number of gamma-ray beams.
Test assembly
CT measurement
Chordal measurement
Subchannel
2(X and Y direction)
2(X and Y direction)
Rod bundle
—
6 beams × 2 × 3 section(total 36 beams)Table 12
Time required to perform void fraction measurements.
Item
CT measurement
Chordal measurement
Steady state
Time needed
5 s/step ×T33 ×  R17 step(it takes 2 h)
100 ssampling cycle 0.1 s
Measurement
2 times
3 times
Transient
Time needed
—
200 s
Measurement
—
1 timeFigure 12
Relation Between Chordal and CT Averaged Densities (for S1).


## 2.5.2. Departure from Nucleate Boiling

The bundle power was gradually increased in fine steps to the expected vicinity of DNB, which was based on previous analysis operator experience. The onset of DNB was confirmed by a rod temperature rise greater than 11°C (20°F) as measured by the thermocouples seen in Figure13. The DNB power is defined as the power corresponding to the step immediately preceding the step in which this temperature rise is seen. The accuracy of the process parameters involved in this process is seen in Table 13. The exit fluid temperatures were measured by the thermocouples shown in Figure 14.Table 13
Accuracy of process parameters in DNB measurement.
Quantity
Accuracy
Process parameters
Pressure
1%
Flow
1.5%
Power
1%
Fluid temperature
1 CelsiusFigure 13
Locations of thermocouples for test assemblies.Figure 14
Diagram of fluid temperature measurement setup.


## 2.6. Experimental Data Analysis

Based on the participants feedback, several studies were performed to determine the validity of the benchmark data. These included (1) study performed on recalculation of the void fraction, (2) study performed on calculation of the quality, (3) study performed on the subchannel CT masking, and (4) study performed on bundle test section downcomer region. In this paper, the studies (1), (2), and (4) performed by the benchmark team are presented.


### 2.6.1. Recalculation of Void Fraction

At the first PSBT workshop, it was noted by several participants that the “measured” void fractions (which were actually calculated from measured density data) were not consistent with void fractions calculated using the measured densities. As a result, a study was performed to recalculate the void fraction and the quality for each test case in the void distribution benchmark.Solving the standard equation for mixture density (ρ-=αρg+(1-α)ρf) for the void fraction gives α=(ρ--ρf)/(ρg-ρf), where the liquid and vapor densities (ρf and ρg, resp.) are assumed to be at saturation and the mixture density is taken from the benchmark data. After the recalculation (see Figure 15), it was noted that the measured void fraction was consistently higher than the recalculated void fraction. This recalculation was only performed for the subchannel test cases since those are the only test cases for which fluid density data was available.Figure 15
Deviation of measured void fraction from recalculated void fraction.


### 2.6.2. Calculation of Quality

Upon completion of the study performed on void distribution, the benchmark team began a study of the calculation of quality based on the experimentally determined densities.It is recalled that the quality can be expressed using mixture enthalpy asx=(hmix-hf)/(hg-hf), where hf and hg are the liquid and vapor enthalpies, respectively. A number of different expressions were derived to determine the mixture enthalpy in the test sections assuming conservation of energy. After verifying that the axial power distribution was normalized for both the uniform and cosine power shapes, the following equations were obtained.All four subchannel test sections utilized a uniform axial power distribution. Thus, for all subchannel assemblies, the mixture enthalpy at the measurement section can be given by(1)hmix=hin+(1400mm1555mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h]).Figure16 shows the resulting deviation of the experimental quality from the recalculated quality.Figure 16
Deviation of measured quality from recalculated quality for test series 1 and 2.Assembly B5 utilized a uniform axial power distribution, so the mixture enthalpies at the three measurement locations can be given as(2)hmix=hin+(2216mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(lower elevation),hmix=hin+(2669mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(middle elevation),hmix=hin+(3177mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(upperelevation).Assemblies B6 and B7 utilized a cosine axial power shape. Recalling that the power shape is normalized, it is possible to determine what fraction of the total power has been imparted to the fluid between the flow inlet and the measurement sections. The mixture enthalpies for these two assemblies can be given as(3)hmix=hin+0.6598×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(Lower Elevation),hmix=hin+0.8172×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(MiddleElevation),hmix=hin+0.9353×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(UpperElevation).After applying these equations and calculating the qualities for each case, the following results were obtained. Figures17, 18, 19, and 20 show the deviation of experimental quality from recalculated quality for test series 5, 6, 7, and 8 (resp.).Figure 17
Deviation of measured quality from recalculated quality for test series 5.Figure 18
Deviation of measured quality from recalculated quality for test series 6.Figure 19
Deviation of measured quality from recalculated quality for test series 7.Figure 20
Deviation of measured quality from recalculated quality for test series 8.


## 2.6.1. Recalculation of Void Fraction

At the first PSBT workshop, it was noted by several participants that the “measured” void fractions (which were actually calculated from measured density data) were not consistent with void fractions calculated using the measured densities. As a result, a study was performed to recalculate the void fraction and the quality for each test case in the void distribution benchmark.Solving the standard equation for mixture density (ρ-=αρg+(1-α)ρf) for the void fraction gives α=(ρ--ρf)/(ρg-ρf), where the liquid and vapor densities (ρf and ρg, resp.) are assumed to be at saturation and the mixture density is taken from the benchmark data. After the recalculation (see Figure 15), it was noted that the measured void fraction was consistently higher than the recalculated void fraction. This recalculation was only performed for the subchannel test cases since those are the only test cases for which fluid density data was available.Figure 15
Deviation of measured void fraction from recalculated void fraction.


## 2.6.2. Calculation of Quality

Upon completion of the study performed on void distribution, the benchmark team began a study of the calculation of quality based on the experimentally determined densities.It is recalled that the quality can be expressed using mixture enthalpy asx=(hmix-hf)/(hg-hf), where hf and hg are the liquid and vapor enthalpies, respectively. A number of different expressions were derived to determine the mixture enthalpy in the test sections assuming conservation of energy. After verifying that the axial power distribution was normalized for both the uniform and cosine power shapes, the following equations were obtained.All four subchannel test sections utilized a uniform axial power distribution. Thus, for all subchannel assemblies, the mixture enthalpy at the measurement section can be given by(1)hmix=hin+(1400mm1555mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h]).Figure16 shows the resulting deviation of the experimental quality from the recalculated quality.Figure 16
Deviation of measured quality from recalculated quality for test series 1 and 2.Assembly B5 utilized a uniform axial power distribution, so the mixture enthalpies at the three measurement locations can be given as(2)hmix=hin+(2216mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(lower elevation),hmix=hin+(2669mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(middle elevation),hmix=hin+(3177mm3658mm)×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(upperelevation).Assemblies B6 and B7 utilized a cosine axial power shape. Recalling that the power shape is normalized, it is possible to determine what fraction of the total power has been imparted to the fluid between the flow inlet and the measurement sections. The mixture enthalpies for these two assemblies can be given as(3)hmix=hin+0.6598×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(Lower Elevation),hmix=hin+0.8172×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(MiddleElevation),hmix=hin+0.9353×Q[kW]×3600[s/h](A[m2]G[kg/m2h])(UpperElevation).After applying these equations and calculating the qualities for each case, the following results were obtained. Figures17, 18, 19, and 20 show the deviation of experimental quality from recalculated quality for test series 5, 6, 7, and 8 (resp.).Figure 17
Deviation of measured quality from recalculated quality for test series 5.Figure 18
Deviation of measured quality from recalculated quality for test series 6.Figure 19
Deviation of measured quality from recalculated quality for test series 7.Figure 20
Deviation of measured quality from recalculated quality for test series 8.


## 3. Selected Examples of Comparative Analysis of Participant Results

Detailed comparative analysis of submitted participants’ results for different exercises of the two phases of the benchmark has been reported in [5, 6]. In this paper, selected examples are shown to illustrate the methodology of comparisons and indicate the observed tendencies.


### 3.1. Phase I—Void Distribution

There were a total of twenty (20) participants for exercise 1 of phase I—Steady-state single subchannel benchmark [5]. Sixteen (16) of these participants submitted void fraction results, and eight (8) submitted axial void distribution results. Table 14 lists these participants as well as the codes for which results were submitted.Table 14
Phase I, exercise 1 participants and code list.
Participant
Code
Code type
ANL
STAR-CD
CFD
ANSYS
ANSYS
CFD
EDF
NEPTUNE
CFD
GRS
ANSYS-CFX
CFD
HZDR
ANSYS-CFX
CFD
PSI
STAR-CD
CFD
EDF
THYC
Porous media
AREVA
F-COBRA-TF
Subchannel
CEA-Saclay
FLICA-OVAP
Subchannel
CSA
VIPRE
Subchannel
JNES
CHAMP-ITA
Subchannel
KAERI
MATRA
Subchannel
KIT
SUBCHANFLOW
Subchannel
McMaster
ASSERT-PV
Subchannel
PSI
FLICA
Subchannel
CEA-Grenoble
CATHARE 3
System
IRSN
CATHARE 2
System
KTH
TRACE
System
PSI
TRACE
System
UNIPI
CATHARE 2
SystemThe axial nodalization of the problem varied widely depending on the code type. CFD codes, as expected, had significantly more nodes than the system and subchannel codes: both NEPTUNE [7] and ANSYS (ANSYS) [8] used nodalizations of ~100,000–1,000,000 nodes, while both KTH and PSI used about 30 axial nodes for TRACE. CATHARE 2 (UNIPI) [9] and STAR-CD (PSI) both used approximately 100 axial nodes.Finally, a variety of wall drag and heat transfer coefficients were also used. THYC, for example, used Gautier for monophasic and Chen for diphasic heat transfer, while FLICA (PSI) [10] used the Blasius correlation for single-phase friction losses with the Friedel correlation as a two-phase multiplier. A number of codes used Dittus-Boelter for single-phase heat transfer and a variety of correlations (Jens-Lottes for FLICA (PSI)) were used.Statistical methodology has been utilized for comparative analysis. Mean error and standard deviation were calculated for each code for each test series. The mean error (whereN is the total number of test cases) was represented as α-=∑n=1Nαn/N, where the void fraction error for test case “n” is represented as αn=(αcoden-αexpn).The standard deviation was given asσ=±∑n=1N(αn-α-)2/(N-1).Figure21 summarizes the mean error of participant void fraction calculations as compared to measurements, while Figure 22 summarizes the standard deviation of those results.Figure 21
Void fraction mean error by test series.Figure 22
Void fraction standard deviation by test series.There was no clear bias in the calculation of void fraction for any of the four subchannels. Although some of the codes consistently predicted the correct thermal equilibrium quality, there was a tendency to overpredict it at the low elevation and underpredict it at the high elevation.There were a total of seventeen (17) participants for exercise 2 of phase I—steady-state bundle benchmark [5]. All sixteen submitted void fraction results, while two submitted axial void distribution results. Table 15 lists these participants as well as the codes for which results were submitted.Table 15
Phase I, exercise 2 participants and code list.
Participant
Code
Code type
EDF
THYC
Porous media
AREVA
F-COBRA-TF
Subchannel
CEA-Saclay
FLICA-OVAP
Subchannel
CSA
VIPRE
Subchannel
JNES
CHAMP-ITA
Subchannel
KAERI
MATRA
Subchannel
KIT
SUBCHANFLOW
Subchannel
McMaster
ASSERT-PV
Subchannel
NRI
VIPRE
Subchannel
PSI
FLICA
Subchannel
WEC/INVAP
VIPRE
Subchannel
CEA-Grenoble
CATHARE 3
System
Chalmers
RELAP-5
System
KTH
TRACE
System
PSI
TRACE
System
UNIPI
CATHARE 2
System
USNRC
TRACE
SystemSome modeling considerations specific to this exercise include increasing the number of axial nodes and modeling spacer grids. A number of the codes (such as THYC) did not model spacer grid effects, while others (CATHARE 3 (CEA-Grenoble) [11]) applied a pressure loss one cell in front of the spacer grid. A number of cross-flow models were also used. THYC used Chexal-Lellouche for drift flux as well as a no-slip condition and Cheng-Todreas for turbulent diffusivity, while CATHARE 3 (CEA-Grenoble) used a turbulent dispersion term. Yet, other codes did not consider cross-flow effects. Statistical methodology also was utilized to calculate the mean error and standard deviation for each measurement location for each test series. Figures 23, 25, 27, and 29 summarize the error of participants’ void fraction predictions, while Figures 24, 26, 28, and 30 summarize the standard deviation of these results.Figure 23
Test series 5 void fraction mean error.Figure 24
Test series 5 void fraction standard deviation.Figure 25
Test series 6 void fraction mean error.Figure 26
Test series 6 void fraction standard deviation.Figure 27
Test series 7 void fraction mean error.Figure 28
Test series 7 void fraction standard deviation.Figure 29
Test series 8 void fraction mean error.Figure 30
Test series 8 void fraction standard deviation.It was noted that the codes consistently overpredicted the void fraction at the lower elevation in the bundle. However, the results were generally improved at higher elevations, although some underprediction could be seen. Both of these observations are clearly represented in the plots of mean error for each test series.The majority of the codes also consistently predicted the correct thermal equilibrium quality at the lower elevations, with the only exceptions being KTH’s TRACE and UNIPI’s CATHARE, which overpredicted the quality. All of the codes tended to underpredict the quality at the upper bundle elevations.


### 3.2. Phase II—Departure of Nucleate Boiling

The list of the participants as well as the codes for which results were submitted for exercise 2 of phase II (steady-state DNB) benchmark is shown in Table16.Table 16
Phase II, exercise 2 participants.
Participant
Code
Code type
EDF
THYC
Porous media
CEA-Saclay
FLICA-OVAP
Subchannel
KAERI
MATRA
Subchannel
KIT
SUBCHANFLOW
Subchannel
McMaster
ASSERT-PV
Subchannel
PSI
FLICA
Subchannel
NRI
VIPRE
Subchannel
CEA-Grenoble
CATHARE 3
System
KTH
TRACE
SystemThe codes used in this exercise varied widely in the DNB modeling approach used. ASSERT-PV, TRACE (KTH), CATHARE 3, and FLICA (PSI) all used Groeneveld lookup tables, while MATRA used the EPRI CHF correlation and FLICA-OVAP [12] used the Shah correlation.Figure31 shows the mean error of each code for DNB power for each test series. Figure 32 shows the standard deviation of each code for DNB power for each test series.Figure 31
Mean error of predicted DNB power.Figure 32
Standard deviation of predicted DNB power.The codes were generally able to calculate the DNB power satisfactorily, and there was no observable bias across test series. The DNB power was consistently overpredicted in Test series 0, while it was underpredicted in test series 2, 3, 4, and 13. There was also considerable variation in the predictions of axial elevation of first detected DNB. It should be noted that the measured data represents the first thermocouple at which DNB was detected. Therefore, it is the latest (axially speaking) that the onset of DNB would have occurred and is not an exact value as DNB could have occurred lower on the bundle.There were a total of eight (8) participants for exercise 3 of phase II—transient DNB benchmark [6]. Table 17 lists these participants as well as the codes for which results were submitted.Table 17
Phase II, exercise 3 participants.
Participant
Code
Code type
EDF
THYC
Porous media
CEA-Saclay
FLICA-OVAP
Subchannel
KAERI
MATRA
Subchannel
KIT
SUBCHANFLOW
Subchannel
McMaster
ASSERT-PV
Subchannel
PSI
FLICA
Subchannel
NRI
VIPRE
Subchannel
KTH
TRACE
SystemThe same statistical methodology was utilized also for these comparisons. The mean error and standard deviation were calculated for each measurement location for each test series. Figure33 shows the mean error of each code for time of the detected DNB for each test series. Figure 34 shows the mean error of each code for DNB power for each test series.Figure 33
Mean error of detected time of DNB.Figure 34
Mean error of predicted DNB power.In general, codes were not able to predict the time of DNB during the simulated transients. Most of the codes calculate earlier DNB occurrence for both bundle types (with, test series 11, and without, test series 12, thimble rod in the central region). THYC always predicts later DNB. FLICA-OVAP and SUBCHANFLOW [13] predictions show dependence on the radial power shape for both DNB power and DNB time. MATRA [14] seems to give reasonable predictions of the DNB time but significantly overpredicts the DNB power.


## 3.1. Phase I—Void Distribution

There were a total of twenty (20) participants for exercise 1 of phase I—Steady-state single subchannel benchmark [5]. Sixteen (16) of these participants submitted void fraction results, and eight (8) submitted axial void distribution results. Table 14 lists these participants as well as the codes for which results were submitted.Table 14
Phase I, exercise 1 participants and code list.
Participant
Code
Code type
ANL
STAR-CD
CFD
ANSYS
ANSYS
CFD
EDF
NEPTUNE
CFD
GRS
ANSYS-CFX
CFD
HZDR
ANSYS-CFX
CFD
PSI
STAR-CD
CFD
EDF
THYC
Porous media
AREVA
F-COBRA-TF
Subchannel
CEA-Saclay
FLICA-OVAP
Subchannel
CSA
VIPRE
Subchannel
JNES
CHAMP-ITA
Subchannel
KAERI
MATRA
Subchannel
KIT
SUBCHANFLOW
Subchannel
McMaster
ASSERT-PV
Subchannel
PSI
FLICA
Subchannel
CEA-Grenoble
CATHARE 3
System
IRSN
CATHARE 2
System
KTH
TRACE
System
PSI
TRACE
System
UNIPI
CATHARE 2
SystemThe axial nodalization of the problem varied widely depending on the code type. CFD codes, as expected, had significantly more nodes than the system and subchannel codes: both NEPTUNE [7] and ANSYS (ANSYS) [8] used nodalizations of ~100,000–1,000,000 nodes, while both KTH and PSI used about 30 axial nodes for TRACE. CATHARE 2 (UNIPI) [9] and STAR-CD (PSI) both used approximately 100 axial nodes.Finally, a variety of wall drag and heat transfer coefficients were also used. THYC, for example, used Gautier for monophasic and Chen for diphasic heat transfer, while FLICA (PSI) [10] used the Blasius correlation for single-phase friction losses with the Friedel correlation as a two-phase multiplier. A number of codes used Dittus-Boelter for single-phase heat transfer and a variety of correlations (Jens-Lottes for FLICA (PSI)) were used.Statistical methodology has been utilized for comparative analysis. Mean error and standard deviation were calculated for each code for each test series. The mean error (whereN is the total number of test cases) was represented as α-=∑n=1Nαn/N, where the void fraction error for test case “n” is represented as αn=(αcoden-αexpn).The standard deviation was given asσ=±∑n=1N(αn-α-)2/(N-1).Figure21 summarizes the mean error of participant void fraction calculations as compared to measurements, while Figure 22 summarizes the standard deviation of those results.Figure 21
Void fraction mean error by test series.Figure 22
Void fraction standard deviation by test series.There was no clear bias in the calculation of void fraction for any of the four subchannels. Although some of the codes consistently predicted the correct thermal equilibrium quality, there was a tendency to overpredict it at the low elevation and underpredict it at the high elevation.There were a total of seventeen (17) participants for exercise 2 of phase I—steady-state bundle benchmark [5]. All sixteen submitted void fraction results, while two submitted axial void distribution results. Table 15 lists these participants as well as the codes for which results were submitted.Table 15
Phase I, exercise 2 participants and code list.
Participant
Code
Code type
EDF
THYC
Porous media
AREVA
F-COBRA-TF
Subchannel
CEA-Saclay
FLICA-OVAP
Subchannel
CSA
VIPRE
Subchannel
JNES
CHAMP-ITA
Subchannel
KAERI
MATRA
Subchannel
KIT
SUBCHANFLOW
Subchannel
McMaster
ASSERT-PV
Subchannel
NRI
VIPRE
Subchannel
PSI
FLICA
Subchannel
WEC/INVAP
VIPRE
Subchannel
CEA-Grenoble
CATHARE 3
System
Chalmers
RELAP-5
System
KTH
TRACE
System
PSI
TRACE
System
UNIPI
CATHARE 2
System
USNRC
TRACE
SystemSome modeling considerations specific to this exercise include increasing the number of axial nodes and modeling spacer grids. A number of the codes (such as THYC) did not model spacer grid effects, while others (CATHARE 3 (CEA-Grenoble) [11]) applied a pressure loss one cell in front of the spacer grid. A number of cross-flow models were also used. THYC used Chexal-Lellouche for drift flux as well as a no-slip condition and Cheng-Todreas for turbulent diffusivity, while CATHARE 3 (CEA-Grenoble) used a turbulent dispersion term. Yet, other codes did not consider cross-flow effects. Statistical methodology also was utilized to calculate the mean error and standard deviation for each measurement location for each test series. Figures 23, 25, 27, and 29 summarize the error of participants’ void fraction predictions, while Figures 24, 26, 28, and 30 summarize the standard deviation of these results.Figure 23
Test series 5 void fraction mean error.Figure 24
Test series 5 void fraction standard deviation.Figure 25
Test series 6 void fraction mean error.Figure 26
Test series 6 void fraction standard deviation.Figure 27
Test series 7 void fraction mean error.Figure 28
Test series 7 void fraction standard deviation.Figure 29
Test series 8 void fraction mean error.Figure 30
Test series 8 void fraction standard deviation.It was noted that the codes consistently overpredicted the void fraction at the lower elevation in the bundle. However, the results were generally improved at higher elevations, although some underprediction could be seen. Both of these observations are clearly represented in the plots of mean error for each test series.The majority of the codes also consistently predicted the correct thermal equilibrium quality at the lower elevations, with the only exceptions being KTH’s TRACE and UNIPI’s CATHARE, which overpredicted the quality. All of the codes tended to underpredict the quality at the upper bundle elevations.


## 3.2. Phase II—Departure of Nucleate Boiling

The list of the participants as well as the codes for which results were submitted for exercise 2 of phase II (steady-state DNB) benchmark is shown in Table16.Table 16
Phase II, exercise 2 participants.
Participant
Code
Code type
EDF
THYC
Porous media
CEA-Saclay
FLICA-OVAP
Subchannel
KAERI
MATRA
Subchannel
KIT
SUBCHANFLOW
Subchannel
McMaster
ASSERT-PV
Subchannel
PSI
FLICA
Subchannel
NRI
VIPRE
Subchannel
CEA-Grenoble
CATHARE 3
System
KTH
TRACE
SystemThe codes used in this exercise varied widely in the DNB modeling approach used. ASSERT-PV, TRACE (KTH), CATHARE 3, and FLICA (PSI) all used Groeneveld lookup tables, while MATRA used the EPRI CHF correlation and FLICA-OVAP [12] used the Shah correlation.Figure31 shows the mean error of each code for DNB power for each test series. Figure 32 shows the standard deviation of each code for DNB power for each test series.Figure 31
Mean error of predicted DNB power.Figure 32
Standard deviation of predicted DNB power.The codes were generally able to calculate the DNB power satisfactorily, and there was no observable bias across test series. The DNB power was consistently overpredicted in Test series 0, while it was underpredicted in test series 2, 3, 4, and 13. There was also considerable variation in the predictions of axial elevation of first detected DNB. It should be noted that the measured data represents the first thermocouple at which DNB was detected. Therefore, it is the latest (axially speaking) that the onset of DNB would have occurred and is not an exact value as DNB could have occurred lower on the bundle.There were a total of eight (8) participants for exercise 3 of phase II—transient DNB benchmark [6]. Table 17 lists these participants as well as the codes for which results were submitted.Table 17
Phase II, exercise 3 participants.
Participant
Code
Code type
EDF
THYC
Porous media
CEA-Saclay
FLICA-OVAP
Subchannel
KAERI
MATRA
Subchannel
KIT
SUBCHANFLOW
Subchannel
McMaster
ASSERT-PV
Subchannel
PSI
FLICA
Subchannel
NRI
VIPRE
Subchannel
KTH
TRACE
SystemThe same statistical methodology was utilized also for these comparisons. The mean error and standard deviation were calculated for each measurement location for each test series. Figure33 shows the mean error of each code for time of the detected DNB for each test series. Figure 34 shows the mean error of each code for DNB power for each test series.Figure 33
Mean error of detected time of DNB.Figure 34
Mean error of predicted DNB power.In general, codes were not able to predict the time of DNB during the simulated transients. Most of the codes calculate earlier DNB occurrence for both bundle types (with, test series 11, and without, test series 12, thimble rod in the central region). THYC always predicts later DNB. FLICA-OVAP and SUBCHANFLOW [13] predictions show dependence on the radial power shape for both DNB power and DNB time. MATRA [14] seems to give reasonable predictions of the DNB time but significantly overpredicts the DNB power.


## 4. Conclusions

The OECD/NRC PSBT benchmark was designed to provide a set of data for the development and validation of the next generation of thermal-hydraulic codes. It consisted of two phases: a void fraction benchmark and a departure from nucleate boiling benchmark. Data regarding the test sections and conditions was provided to participants for use in calculations. The code results from all participants were then compiled and analyzed. In the development of the benchmark specification, a number of support studies were performed. The experimental void fraction and quality were recalculated using the experimentally determined fluid density for each of the benchmark test cases, and a deviation between these recalculated values and the measured values was observed. Studies were also conducted to determine the effect of the downcomer region on the fluid temperature and the effect of the CT mask on void fraction at subcooled boiling conditions. The participants’ results for each benchmark exercise were analyzed, and conclusions were drawn. In the results for the first phase, it was observed that the codes tended to overpredict the thermal equilibrium quality at lower elevations and underpredict it at higher elevations. There was also a tendency to overpredict void fraction at lower elevations and underpredict it at high elevations for the bundle test cases. The overprediction of void fraction at low elevations is likely caused by the X-ray densitometer measurement method used. Under subcooled boiling conditions, the voids accumulate at heated surfaces (and are therefore not seen in the center of the subchannel, where the measurements are being taken), so the experimentally determined void fractions will be lower than the actual void fraction. The underprediction of void fraction at high elevations is smaller as compared to overprediction at low elevations, and the code results were generally improved at higher elevations. Some of the best results were achieved by codes that used either turbulent mixing or dispersion terms for modeling cross-flow. It was also noted that, for the bundle cases, some of the codes were not correctly calculating the bundle-averaged thermal equilibrium quality, and this may indicate an inability to predict the correct void fraction. A time shift was noted in the void fraction results for the temperature increase transient cases, indicating that the test apparatus may have experienced unexpected heat transfer between the downcomer and test section. This heat transfer is only expected to be of significance in the transient test cases as the steady-state cases allow the system to reach thermal equilibrium.The participants’ results for each benchmark exercise of phase II were analyzed, and conclusions were drawn. Based on the final results for the second phase, a number of conclusions were drawn. The codes were generally able to replicate the results of the mixing test, staying within 5% of the experimental fluid exit temperature. THYC showed the largest errors, while all codes had difficulty with the test case 01–1237. The uncertainty related to the spacer grid orientation in the mixing tests may explain the asymmetries seen in the experimental data, but no sensitivity study was performed to determine the exact effect on the flow, so this uncertainty cannot be quantified. The improvement between FLICA-OVAP and FLICA was noticeable, with the main difference between the two codes being that FLICA-OVAP implicitly assigned flow regime based on drift flux, while FLICA assumes single phase flows. In exercises 2 and 3, the codes were generally able to predict the DNB power as well as the axial location of the onset of DNB (for the steady-state cases) and the time of DNB (for the transient cases). It was noted that the codes that used the EPRI CHF correlation (such as MATRA and SUBCHANFLOW) had the lowest mean error in exercise 2 for predicted DNB power.



---

*Source: 946173-2013-03-23.xml*