**World Journal of Nuclear Science and Technology**

Vol.04 No.04(2014), Article ID:50339,12 pages

10.4236/wjnst.2014.44025

Assessment and Development of Two Phase Turbulent Mixing Models for Subchannel Analysis Relevant to BWR

Mohit P. Sharma^{1}, Arun K. Nayak^{2}

^{1}Homi Bhabha National Institute, Mumbai, India

^{2}Bhabha Atomic Research Centre, Reactor Engineering Division, Mumbai, India

Email: mohit.lucky07@gmail.com

Copyright © 2014 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 4 July 2014; revised 11 August 2014; accepted 26 August 2014

ABSTRACT

Determination of turbulent mixing rate of two phase flow between neighboring subchannels is an important aspect of sub channel analysis in reactor rod bundles. Various models have been developed for two phase turbulent mixing rate between subchannels. These models show that turbulent mixing rate is strongly dependent on flow regimes; their validity was examined against specific or limited experiments. It is vital to evaluate these models by comparing the predicted two phase turbulent mixing rate with available experimental data conducted for various subchannel geometries and operating conditions. This paper describes evaluation of different models for two phase turbulent mixing rate for both gas and liquid phase against large range of experimental data which are obtained from various subchannel geometries. The results indicate that there is large discrepancy between the predicted and experimental data for turbulent mixing rate. This paper provides important shortcoming of the previous work and need for the development of a new model. In the view of this, a two phase flow model is presented, which predicts both liquid and gas phase turbulent mixing rate between adjacent sub channels of reactor rod bundles. The model presented here is for slug churn flow regime, which is dominant as compared to the other regimes like bubbly flow and annular flow regimes, since turbulent mixing rate is the highest in slug churn flow regime. The present model has been tested against low pressure and temperature air-water and high pressure and temperature steam-water experimental data found that it shows good agreement with available experimental data.

**Keywords:**

Subchannel Analysis, Two Phase Turbulent Mixing Model, Turbulent Mixing Rate

1. Introduction

The fluid transfer among subchannels is explained by three mechanisms i.e. turbulent mixing, void drift and diversion cross flow (Lahey and Moody [1] ). The fluid exchange due to turbulent mixing is because of turbulent fluctuation. In this mechanism of fluid exchange, neither net mass transfer nor net volume transfer between adjacent subchannels occurs. Second mechanism is void drift which occurs even in absence of pressure difference. Void drift is due to redistribution of non-equilibrium flow to attain equilibrium flow. Third mechanism is diversion cross flow which occurs due to lateral pressure difference between adjacent subchannels.

1.1. Review of Experiments for Subchannel Analysis

The available experiment on two phase turbulent mixing rate between adjacent subchannels is listed in Table 1.

Insights from Previous Experiments

The two phase turbulent mixing experiments are performed by Walton [2] , Rudzinski [3] , Singh K. S. [4] , Kawahara et al. [5] , Sadatomi et al. [6] and Kawahara et al. [7] . These experiments provide important insights as given below:

a) The total two phase turbulent mixing rate is sum of liquid and gas phase turbulent mixing rate and it is strongly related to flow regimes. The liquid phase turbulent mixing rate starts at zero quality, which increases in bubbly flow reaches maximum value in slug churn and then decreases beyond churn-annular flow transition. The gas phase mixing rate starts near zero value of quality, which reaches maximum in slug churn flow and then decreases with increase in quality. It is thus rational to consider the turbulent mixing separately in each flow pattern.

b) Two phase mixing rate depends on gap between the subchannels. On increasing the gap between the subchannels, mixing rate increases.

c) The two phase turbulent mixing rate increases with increase in mass flux.

d) The two phase turbulent mixing rate decreases with increase in pressure.

1.2. Review of Models for Subchannel Analysis

The available models on two phase turbulent mixing rate between adjacent subchannels are listed in Table 2.

2. Evaluation of Turbulent Mixing Model

In this section, we evaluate the turbulent mixing models like Bues [8] model, Kazimi and Kelly [9] model, Kawahara et al. [10] model and Carlucci et al. [11] against the data obtained from various subchannel experiments of two phase turbulent mixing as discussed in section 1.1. For evaluation, we have compared the measured (experimental) liquid phase turbulent mixing rate with predicted liquid phase turbulent mixing rate and measured (experimental) gas phase turbulent mixing rate with predicted gas phase turbulent mixing rate in two phase flow. The error analysis has been done to find out maximum, minimum and average error between measured and predicted value of both liquid and gas phase turbulent mixing rate. The error analysis shows how predicted value by turbulent mixing models differs from measured experimental values.

Max. Error (1)

Min. Error (2)

The average error is calculated as

(3)

where no. of data points and

Table 1. Description of available data on two phase turbulent mixing rate.

Table 2. Description of available models on two phase turbulent mixing rate.

2.1. Evaluation of Model of Bues [8]

Bues [8] model, the calculated liquid and gas turbulent mixing rate shows large discrepancy, when compared against measured liquid and gas phase mixing rate as seen in Figure 1(a) and Figure 1(b) respectively.

2.2. Evaluation of Model of Kazimi and Kelly [9]

In Kazimi and Kelly’s [9] , the calculated liquid and gas turbulent mixing rate shows large discrepancy, when compared against measured liquid and gas phase mixing rate as seen in Figure 2(a) and Figure 2(b) respectively.

2.3. Evaluation of Model of Kawahara et al. [10]

Kawahara et al. [10] , the calculated liquid and gas turbulent mixing rate shows large discrepancy, when compared against measured liquid and gas phase mixing rate as seen in Figure 3(a) and Figure 3(b) respectively.

2.4. Evaluation of Model of Carlucci et al. [11]

In Carlucci et al. [11] model, the calculated liquid and gas turbulent mixing rate shows large discrepancy, when compared against measured liquid and gas phase mixing rate as seen in Figure 4(a) and Figure 4(b) respectively.

2.5. Insights from Previous Models

In this paper the liquid and gas phase mixing rate in two phase flow have been predicted using the models of Bues [8] , Kazimi and Kelly [9] and Carlucci et al. [11] . All these models consider all the flow regimes. Only Kawahara et al. [10] model considered models for different flow regimes. The assessment of these models shows that there is large discrepancy between models and experimental data which is shown in Table 3.

Assessment of these models provide important shortcoming which are as follows.

a) Array effect: In all these models except Kawahara et al. [10] model, the array effect like Square-Square,

(a) (b)

Figure 1. Comparison of the predictability of Bues [8] model model against subchannels experiments for liquid and gas phase turbulent mixing rate in two phase flow.

(a) (b)

Figure 2. Comparison of the predictability of Kazimi and Kelly’s [9] model model against subchannels experiments for liquid and gas phase turbulent mixing rate in two phase flow.

(a) (b)

Figure 3. Comparison of the predictability of Kawahara et al. [10] model model against subchannels experiments for liquid and gas phase turbulent mixing rate in two phase flow.

Rectangular-Rectangular, and Triangular-Triangular subchannel array has not been considered.

b) Subchannel size effect: Carlucci [11] model doesn’t predict well, when area of subchannel is very small, which is shown by large error in Triangular-Triangular subchannel experiment of Kawahara [7] where area of subchannel is very less (~16.6 mm^{2}).

(a) (b)

Figure 4. Comparison of the predictability of Carlucci et al. [11] model against subchannels experiments for liquid and gas phase turbulent mixing rate in two phase flow.

Table 3. Error analysis between calculated liquid turbulent mixing and measured liquid turbulent mixing rate in two phase flow.

c) Gap size effect: Models of Bues [8] , Kazimi and Kelly [9] and Carlucci et al. [11] , doesn’t predict well, when the gap between subchannels is more than 2.1 mm.

d) Pressure effect: In Kawahara et al. [10] model, the effect of pressure has not been considered.

3. Model Developments

Since previous models have large errors, there is need to develop a new turbulent mixing model which can predict well for various subchannel geometries. A slug-churn flow model is proposed to predict liquid and gas phase mixing rate. The model has been tested against low pressure and temperature (ambient) air-water and high pressure and temperature steam-water experimental data found that it shows good agreement with available experimental data.

From the insights of previous experiments and models of two phase turbulent mixing rate, it can be inferred that the two phase turbulent mixing depends strongly on the subchannel quality, mass flux, pressure and subchannels geometry. Thus the liquid turbulent mixing number for two phase flow in subchannels can be expressed as

(5)

where

(6)

(7)

(8)

The coefficient and exponent were obtained by fitting the test data of Rudzinski [3] , Kawahara et al. [5] , and Kawahara et al. [7] plotted on dimensionless liquid mixing number against mixture Reynolds number as shown in Figure 5.

The equation so obtained is given by relationship

(9)

The equation for gas phase turbulent mixing rate in two phase flow can be written as follows

(10)

where

(11)

(12)

(13)

The coefficient and exponent were obtained from the test data of subchannel experiments of Rudzinski [3] , Kawahara et al. [5] , and Kawahara et al. [7] plotted on dimensionless gas mixing number against combined volumetric gas fraction and mixture Reynolds number as shown in Figure 6.

The equation so obtained is given by relationship

(14)

3.1. Modeling of Geometrical Influence of Subchannel

Incorporation of Gap and Centroidal Distance between Subchannels

The two phase turbulent mixing is affected by various parameters such as subchannels geometry, spacer and gap

Figure 5. Liquid mixing rate.

spacing between subchannels. Previous models like Carlucci [11] model considers gap to rod diameter ratio, Kawahara et al. [10] model considers pitch to rod diameter ratio, Bues [8] model and Kazimi and Kelly [9] model considers gap to hydraulic diameter ratio. In the present model, the gap spacing to centroidal distance ratio of subchannels (Figure 7) is considered. The spacer effect is not considered in present model.

The equation for liquid mixing number in two phase flow can be represented by

(15)

where

(16)

The equation of gap to centroid factor for liquid mixing rate can be expressed as best fit by

(17)

where

(18)

Figure 6. Gas mixing rate.

Figure 7. Representation of geometrical parameter in R-R, S-S and T-T subchannel array.

The coefficient and exponent were obtained by the test data of subchannel experiments of Rudzinski [3] , Kawahara et al. [5] , and Kawahara et al. [7] plotted on dimensionless gap to centroid factor against combined gap to centroidal distance ratio and volumetric liquid fraction against combined gap to centroidal distance ratio and volumetric liquid fraction of individual subchannel geometries (R-R, T-T, and S-S) as shown in Figures 8(a)-(c).

The equation for gas phase can be represented by

(19)

where

(20)

The equation of gap to centroid factor for gas mixing rate can be expressed as best fit by

(21)

The coefficient and exponent were obtained by the test data of subchannel experiments of Rudzinski [3] , Kawahara et al. [5] , and Kawahara et al. [7] plotted on dimensionless gap to centroid factor against combined gap to centroidal distance ratio and volumetric liquid fraction of individual subchannel geometries (R-R, T-T, and S-S) as shown in Figures 9(a)-(c).

3.2. Modeling of Pressure Effect

Carlucci [11] is the only model which considers pressure effect in terms of bubble diameter, which changes with change in pressure. However, bubble diameter is difficult to predict in two phase flow since the size of bubble

(a) (b) (c)

Figure 8. The coefficient and exponent for various subchannel geometry in liquid phase mixing rate.

(a) (b) (c)

Figure 9. The coefficient and exponent for various subchannel geometry in gas phase mixing rate.

does not have a single value for a particular operating condition. However, the bubble diameter which has strong effect on void fraction depends on the surface tension of fluid. Hence present model considers the surface tension of fluid to model the effect of pressure.

Thus the effect of pressure is represented by the following expression

(22)

where is the best fit for pressure correction factor

surface tension at high temperature,

surface tension at reference temperature i.e. ambient temperature

The correlation so obtained by Equation (9) and Equation (14) are modified by introducing gap to centroid factor and pressure dependent factor. The modified equation for liquid and gas phase are as follows

(23)

(24)

4. Model Evaluation

In this chapter, the model proposed 3 is evaluated by comparing the prediction from present model with experimental data in a two phase slug churn flow regime

4.1. Test against Low Pressure and Temperature (Ambient) Air-Water Experimental Data

We performed error analysis for liquid and gas phase turbulent mixing rate and found that max error, min error and gross mean error for liquid phase mixing rate considering all subchannel geometry is about +91.7%, −54.3% and −4.27% respectively. The error analysis for individual geometry is shown in Table 4. Comparison between calculated and measured liquid turbulent mixing rate in two phase flow is shown in Figure 10.

The max error, min error and gross average error for gas phase mixing rate considering all subchannels geometry is about +66.2%, −55.7% and −3.29% respectively. The error analysis for individual subchannels geometry is shown in Table 5. Comparison between calculated and measured gas phase turbulent mixing rate in two phase flow is shown in Figure 11.

However considering the case of triangular-triangular subchannel experiment (Kawahara [7] ), some of data points (yellow symbol in Figure 11) in gas mixing rate showing more error between calculated and measured mixing rate. The gross mean error is about 1380%. The reason behind showing more error is not yet perfectly understood.

4.2. Test against High Pressure and Temperature Steam-Water Experimental Data

The error analysis has been performed against high pressure and temperature steam-water experiment (Rowe and angel [12] ) of 52 bar, 255˚C and 28 bar, 215˚C with mass flux 1356.8 and 2712.5 kg/m^{2}s for total turbulent mixing rate and found that max error, min error and gross average error for total mixing rate considering square- square subchannel geometry (mm and mm) is about +79.3%, −45.4% and +9.94%. Comparison between calculated and measured liquid turbulent mixing rate in two phase flow is shown in Figure 12.

4.3. Limitations of Proposed Model

In the present model, the geometrical and pressure effect has considered. However proposed model have following limitation

a) The spacer effect is not considered in present model.

b) This model is only valid for slug churn flow regime

Table 4. Error analysis between calculated liquid turbulent mixing and measured liquid turbulent mixing rate in two phase flow.

Figure 10. Comparison of the predictability of present model against subchannels experiment for liquid phase turbulent mixing rate in two phase flow.

Figure 11. Comparison of the predictability of present model against subchannels experiment for Gas phase turbulent mixing rate in two phase flow.

Table 5. Error analysis between calculated gasturbulent mixing and measured gas turbulent mixing rate in two phase flow.

Figure 12. Comparison of the predictability of present model against subchannel experiment for total turbulent mixing rate in two phase flow.

5. Conclusions

An assessment has been done for turbulent mixing models against the experimental data available in literature and found that there are large discrepancies between predicted turbulent mixing rate models and experimental data. These models are semi-empirical in nature and applicable only for a particular geometry and operating condition. Hence this requires to development of a new model, which predicts well for different subchannel geometries and operating conditions.

In the view if this, a model for slug churn flow regime is proposed in this paper to predict the liquid and gas phase turbulent mixing rate between adjacent subchannels. In this paper, we have defined new dimensionless parameters i.e. liquid mixing number and gas mixing number for two phase turbulent mixing The liquid phase mixing number is a function of mixture Reynolds number whereas the gas phase mixing number is a function of both mixture Reynolds number and volumetric fraction of gas. The effect of gap to centroid spacing between subchannels and subchannel array i.e. square, triangular and rectangular subchannel is also included in present model. The pressure effect is modeled by considering surface tension of the fluid.

To evaluate present model, we tested present model against low pressure and temperature (ambient) air-water experiment of Rudzinski [3] , Kawahara et al. [5] and Kawahara et al. [7] and high pressure and temperature steam-water experiment of Rowe and angel [12] . In case of low pressure and temperature air-water experiment, the max error, min error and gross average error for liquid phase mixing rate is about +91.7%, −54.3% and −4.27% respectively whereas the max error, min error and gross average error for gas phase mixing rate is about +66.2%, −55.7% and −3.29% respectively. In case high pressure and temperature steam-water experiment, the max error, min error and gross mean error for total mixing rate is about +79.3%, −45.4% and +9.94%. Present model showed good agreement with measured mixing rate under low pressure and high pressure condition as compared to earlier model.

References

- Lahey Jr., R.T. andMoody, F.J. (1993) The Thermal Hydraulics of Boiling Water Nuclear Reactor. 2nd Edition, ANS, La Grange Park.
- Walton, F.B. (1969) Turbulent Mixing Measurements for Single-Phase Air, Single-Phase Water and Two-Phase Air-Water Flows in Adjacent Triangular Subchannels. M.Sc. Thesis, Department of Chemical Engineering, University of Windsor, Windsor.
- Rudzinski, K.F. (1970) Two-Phase Turbulent Mixing for Air-Water Flows in Adjacent Triangular Subchannels. M.Sc. Thesis, Department of Chemical Engineering, University of Windsor, Windsor.
- Singh, K.S. (1972) Air-Water Turbulent Mixing in Simulated Rod Bundle Geometries. Ph.D. Thesis, Department of Chemical Engineering, University of Windsor, Windsor.
- Kawahara, A., Sadatomi, M. and Sato, Y. (1997) Two-Phase Turbulent Mixing between Subchannels in a Simulated Rod Bundle Geometry—The Effect of the Number of Gaps between Subchannels. Proceedings of the 2nd Japanese- German Symposium on Multi-Phase Flow, Tokyo, 25-27 September 1997, 55-64.
- Sadatomi, M., Kawahara, A., Kano, K. and Sumi, Y. (2004) Single and Two-Phase of Turbulent Mixing Rate between Adjacent Subchannels in a Vertical 2 × 3 Rod Array Channel. Nuclear Engineering and Design, 202, 27-38.
- Kawahara, A., Sadatomi, M., Kudo, H. and Kano, K. (2006) Single and Two Phase Turbulent Mixing Rate between Subchannels in Triangle Tight Lattice Rod Bundle. JSME International Journal, Series B, 49, 287-295. http://dx.doi.org/10.1299/jsmeb.49.287
- Bues, S.G. (1972) Two Phase Turbulent Mixing Model for Flow in Rod Bundle. Report WAPD-T-2438. Bettis Atomic Power Laboratory. Pittsburgh.
- Kazimi, M.S. and Kelly, J.E. (1983) Formulation of Two Fluid Model for Mixing in LWR Bundle. In: Merilo, M., Ed., Thermal Hydraulics of Nuclear Reactor, ANS, Lagrange Park, Vol. 1, 433-439.
- Kawahara, A., Sadatomi, M. and Sato, Y. (2000) Prediction of Turbulent Mixing Rate of Both Gas and Iquid Phases between Adjacent Subchannels in a Two Phase Slug Churn Flow. Nuclear Engineering and Design, 202, 27-38. http://dx.doi.org/10.1016/S0029-5493(00)00300-9
- Carlucci, L.N., Hammouda, N. and Rowe, D.S. (2003) Two Phase Turbulent Mixing and Buoyancy Drift In Rod Bun- dle. Nuclear Engineering and Design, 227, 65-84. http://dx.doi.org/10.1016/j.nucengdes.2003.08.003
- Rowe, D.S. and Angel, C.W. (1969) Cross Flow Mixing between Parallel Flow Channel during Boiling. Part III: Effect of Spacer on Mixing between Two Channels. Bettelle Nortwest, BNWL-371. http://dx.doi.org/10.2172/4823782