Thermal Systems Oriented Two-Phase Heat Exchanger Models. Focus on Numerical Robustness.

The simulation of complex refrigeration architectures (that usually include vapor compression cycles) provides useful information for design, study and optimization purposes. Such arrangements may include several interconnected systems and a large variety and quantity of components. All components must meet two crucial requirements, namely, low CPU resolution time and high numerical robustness, in order to achieve relatively fast simulations and to prevent solver resolution issues at the architecture level. Among the usual components present, heat exchangers are the most challenging to address considering both the phenomenological and the numerical point of views. A generic heat exchanger model oriented for flexible purposes and meeting the aforementioned requirements has been developed under the Modelica programming language. The model can handle both single-phase and two-phase flows based on a simplified approach that considers three different zones for the refrigerant phase. Its numerical robustness has been extensively tested focusing on different boundary characteristics (definition, values, and signal types) and on demanding operating conditions (null mass flow rate, reversed flow, and reversed heat direction). This document presents the main characteristics of the model and a complete assessment of its numerical behaviour in terms of robustness and CPU time consumption.This project has received funding from the Clean Sky 2 Joint Undertaking (JU) under grant agreement No 886533. The JU receives support from the European Union’s Horizon 2020 research and innovation programme and the Clean Sky 2 JU members other than the Union. Carles Oliet, as a Serra Húnter lecturer, acknowledges the Catalan Government for the support through this Programme.Postprint (author's final draft