Modifications to COBRA-TF to model dispersed flow film boiling with two flow, four field Eulerian-Eulerian approach - Part 1

Ergun S., Hochreiter L. E., Mahaffy J. H.

ANNALS OF NUCLEAR ENERGY, vol.35, no.9, pp.1663-1670, 2008 (SCI-Expanded) identifier identifier

  • Publication Type: Article / Article
  • Volume: 35 Issue: 9
  • Publication Date: 2008
  • Doi Number: 10.1016/j.anucene.2008.02.013
  • Journal Indexes: Science Citation Index Expanded (SCI-EXPANDED), Scopus
  • Page Numbers: pp.1663-1670
  • Hacettepe University Affiliated: Yes


In the case of a postulated loss of coolant accident (LOCA) in a nuclear reactor, an accurate prediction of clad temperature is needed to determine the safety margins. During the reflood phase of the LOCA, when the local void fraction is greater than 80% with the wall temperature above minimum film boiling temperature (T-min). the heat transfer process is dispersed flow film boiling (DFFB). This study has been performed to model DFFB in the reflood phase of a LOCA in a pressurized water reactor (PWR) rod bundle. The COBRA-TF computer code is utilized, since it has a detailed reflood package which takes into account the effect of spacer grids on the local heat transfer. The COBRA-TF code has also been improved to include a four field Eulerian-Eulerian modeling for the two-phase dispersed flow film boiling heat transfer regime. The modifications include adding a small droplet field to COBRA-TF as the fourth field. In addition, the spacer grid models of COBRA-TF have been revised and modified. in the first part of the paper, the results of the code predictions are presented by comparing the experimental data from rod bundle heat transfer (RBHT) experiments with the results of code simulations performed with original and modified code. Measurements and calculations for the heater rod, vapor temperatures and quench front progression have been compared and the results are described in detail. The results of the analysis performed with the modified code indicate the improvement in code predictions for the rod surface temperature, vapor temperature and quench front behavior. The results also indicate the need for improvement in the entrainment and interfacial drag models for the drop fields. The effects of spacer grids on the heat transfer, the models improved and developed for spacer grids and the results of the code calculations with these models are described in the part 2 of the paper. (c) 2008 Elsevier Ltd. All rights reserved.