• Nuclear Engineering and Technology
• /
• v.55 no.3
• /
• pp.1152-1166
• /
• 2023

The bubbly flow of air or steam in subcooled water are investigated here in several test cases, characterized by different pipe sizes, bubble dimensions and flow rates, by means of CFD using a Eulerian-Eulerian approach. The performance of models that differ for the turbulence closure in the continuous phase, as well as for the description of the lift force on the dispersed phase, are compared in detail. When air is considered, the space of the experimental parameters leading to a reasonable performance for the selected models are identified and discussed, while the issues left in the modelling of the concurrent condensation are highlighted for the cases where steam is used.

• Nuclear Engineering and Technology
• /
• v.54 no.6
• /
• pp.2276-2296
• /
• 2022

This study validates the applicability of the CUPID code for simulating subcooled wall boiling under high-pressure conditions against number of DEBORA tests. In addition, a new numerical technique in which the interfacial momentum non-drag forces are calculated at the cell faces rather than the center is presented. This method reduced the numerical instability often triggered by calculating these terms at the cell center. Simulation results showed good agreement against the experimental data except for the bubble sizes in the bulk. Thus, a new model to calculate the Sauter mean diameter is proposed. Next, the effect of the relationship between the bubble departure diameter (D_dep) and the nucleation site density (N) on the performance of the Wall Heat Flux Partitioning (WHFP) model is investigated. Three correlations for D_dep and two for N are grouped into six combinations. Results by the different combinations show that despite the significant difference in the calculated D_dep, most combinations reasonably predict vapor distribution and liquid temperature. Analysis of the axial propagations of wall boiling parameters shows that the N term stabilizes the inconsistences in D_dep values by following a behavior reflective of D_dep to keep the total energy balance. Moreover, ratio of the heat flux components vary widely along the flow depending on the combinations. These results suggest that separate validation of D_dep correlations may be insufficient since its performance relies on the accompanying N correlations.