A new approach to two-phase flow analysis in a rod bundle

Abstract This paper deals with the development of a mathematical and numerical technique for the steady-state subchannel analysis. In view of the importance of the drift-flux model (DFM) in two phase flow analysis, the conservation equations are based on this model. The conservation equations are expressed in terms of five field equations: mixture and momentum continuity; liquid energy equations; gas continuity; and radial cross flow equation for the adjacent subchannels. The numerical algorithm of the single-phase subchannel analysis used in the DIYANA code are extended for the present two phase flow subchannel analysis. The transfer of mass, momentum and energy between adjacent subchannels are split into diversion, turbulent mixing and void drift cross-flow components. The transfer of mass by turbulent mixing is assumed to occur in a volume-for-volume scheme reflecting experimental observations. The phenomenon of lateral vapor drift and turbulent mixing enhancement with flow regime are included in the model. In order to validate the prediction, the GE 3 × 3 rod bundle experiments with both uniform and non-uniform radial power distribution are used. Good agreement has been obtained between the present numerical prediction and the available experimental data. Simulation results showed that the turbulent mixing and void drift phenomena flatten the void fraction and equilibrium vapor quality profile in the rod bundle cross section. The mathematical formulation is considered to be a major step toward a more basic understanding of two-phase flow analysis in fuel rod bundles. Therefore, a fast, efficient and stable numerical technique is presented in this paper for thermal-hydraulic analysis of nuclear reactors.

[1]  J. Tu,et al.  On numerical modelling of low-pressure subcooled boiling flows , 2002 .

[2]  Yong-bum Lee,et al.  Analysis of three different types of blockage in a sodium flow path with the MATRA-LMR-FB code , 2009 .

[3]  N. Zuber,et al.  Average volumetric concentration in two-phase flow systems , 1965 .

[4]  Alberto Teyssedou,et al.  The effect of turbulent mixing models on the predictions of subchannel codes , 1994 .

[5]  G. Wallis.  One Dimensional Two-Phase Flow , 1969 .

[6]  M. Ishii,et al.  Thermo-Fluid Dynamics of Two-Phase Flow , 2007 .

[7]  M. Ishii,et al.  One-dimensional drift-flux model and constitutive equations for relative motion between phases in various two-phase flow regimes , 2003 .

[8]  Seizo Hirao,et al.  Diana — A fast and high capacity computer code for interchannel coolant mixing in rod arrays , 1974 .

[9]  M. Glück,et al.  Sub-channel analysis with F-COBRA-TF – Code validation and approaches to CHF prediction , 2007 .

[10]  M. Avramova,et al.  Evaluation and enhancement of COBRA-TF efficiency for LWR calculations , 2006 .

[11]  M. Glück,et al.  Validation of the sub-channel code F-COBRA-TF: Part II. Recalculation of void measurements , 2008 .

[12]  Paul Coddington,et al.  A study of the performance of void fraction correlations used in the context of drift-flux two-phase flow models , 2002 .

[13]  Kwon-Yeong Lee,et al.  DEVELOPMENT OF A WALL-TO-FLUID HEAT TRANSFER PACKAGE FOR THE SPACE CODE , 2009 .

[14]  Yasunori Bessho,et al.  Boiling transition predictions for BWR fuel bundles by the subchannel analysis program MENUETT , 1987 .

[15]  J. M. Healzer,et al.  A void fraction correlation for generalized applications , 1992 .

[16]  Hadi Davilu,et al.  A 2D numerical simulation of sub-cooled flow boiling at low-pressure and low-flow rates , 2009 .