Investigation of Wall Boiling Closure, Momentum Closure and Population Balance Models for Refrigerant Gas–Liquid Subcooled Boiling Flow in a Vertical Pipe Using a Two-Fluid Eulerian CFD Model

The precise design of heat exchangers in automobile air conditioning systems for more sustainable electric vehicles requires an enhanced assessment of CFD mechanistic models for the subcooled boiling flow of pure eco-friendly refrigerant. Computational Multiphase Flow Dynamics (CMFDs) relies on two-phase closure models to accurately depict the complex physical phenomena involved in flow boiling. This paper thoroughly examines two-phase CMFD flow boiling, incorporating sensitivity analyses of critical parameters such as boiling closures, momentum closures, and population balance models. Three datasets from the DEBORA experiment, involving vertical pipes with subcooled boiling flow of refrigerant at three different pressures and varying levels of inlet liquid subcooling, are used for comparison with CFD simulations. This study integrates nucleate site density and bubble departure diameter models to enhance wall boiling model accuracy. It aims to investigate various interfacial forces and examines the S-Gamma and Adaptive Multiple Size-Group (A-MuSiG) size distribution methods for their roles in bubble break up and coalescence. These proposed approaches demonstrate their efficacy, contributing to a deeper understanding of flow boiling phenomena and the development of more accurate models. This investigation offers valuable insights into selecting the most appropriate sub-closure models for both boiling closure and momentum closure in simulating boiling flows.