A Homogeneous Relaxation Model Algorithm in Density-Based Formulation with Novel Tabulated Method for the Modeling of CO2 Flashing Nozzles

Accurate prediction of the flow physics within nozzles and ejectors working with two-phase carbon dioxide remains particularly challenging. Indeed, the flashing phenomenon in the device often comes along with the formation of the metastable phase, which can severely affect both the flow physics, and the flow rate. Due to the particular challenges associated with building and operating two-phase carbon dioxide experimental setups, computational fluid dynamics constitutes a compelling alternative to study and predict nozzle and ejector flows. In this work, the metastability effects are modeled through a Homogeneous Relaxation Model for the first time expressed in density-based formulation, within the SU2 solver. The thermodynamic relations for the equation of state are tabulated using a novel quad-tree algorithm and directly coupled to the solver, ensuring fast and accurate predictions. The simulation tool is then validated against experimental nozzle pressure profiles and used to compare the commonly used Homogeneous Equilibrium Model to the Homogeneous Relaxation Model in terms of local flow topology and mass flow rates. Most notably, it was found that the Homogeneous Relaxation Model was able to better fit the experimental pressure profiles within the first part of the expansion.