Parameters and dimension of the water-cooled condenser.
Refrigeration for commercial purposes is one of the industrial sectors with the largest energy consumption in the global market. Therefore, research and development of more efficient components such as compressors, condensers, and refrigerants continue to render promising results in terms of GWP and operational costs. However, Due to the urgency typically found in industry to develop prototypes, finding scalable solutions can be challenging. Arguably, this is the case for condenser and evaporators that are designed and assembled under the assumption that refrigeration systems operate at steady condition, where in real circumstances such systems operate under transients based on ambient temperatures or unfavorable startup conditions. The aim of this study is to characterize the thermal and fluid dynamics behavior of refrigerant R404a in a water-cooled condenser at startup conditions. The boundary conditions to solve the CFD simulations are taken from experimental values and set as user defined functions in a commercial software. The results displayed the time dependent oscillatory phase-transition details of the refrigerant throughout the domain.
- Heat transfer
- fluid mechanics
- multi-phase flow
Water-cooled condensers are commonly used in vapor compression refrigeration for industrial applications due to their high rate of heat transfer , simple configuration and the relatively low cost of manufacturability compared to other types of condensers . On the other hand, population growth, and the exponential increase in energy consumption throughout the world keep stressing freshwater resources . Typical components used in refrigeration systems such as compressors, evaporators, condensers, and refrigerants are constantly iterating to provide more energy efficient solutions  in order to minimize the environmental impact due to toxic molecules in refrigerants and reduce the production of residual heat . Unfortunately, the steps imposed by regulatory agencies in United States and Europe sometimes lead to production ready solutions that have not been thoroughly tested for industrial applications. Perhaps, one common example is the race in banning HFC refrigerants for compressed vapor refrigeration cycles . Ideally, the components in the refrigeration system can remain intact by changing refrigerants, but this fails to be the case for most systems since the oil in the compressor and even the compressor specifications need to be tailored to the new refrigerant .
Virtual experimentation using numerical simulations offer great insights at the micro and macroscopic level to support design and research activities of heat exchangers, condensers and evaporators  given that experimental studies can reach elevated costs and access to sensorial points to capture significant physical details are difficult to incorporate in pressurized systems [9, 10]. In fact, the analysis of startup condition for physical components in thermodynamic cycles is particularly challenging given the small-time scales and the large variations detected in the variables of interest . This alone calls for special attention in the sampling rate and physical location of feedback sensors given that dynamic and unstable conditions can be reached in a thermodynamic cycle by unpredicted changes in loads or variations in energetic states, and false positive readings can be measured depending on the flow characteristics through the system.
Concomitantly, the design and fabrication of condensers and heat exchangers require the knowledge and understanding of the physical and chemical phenomenon and other factors such as the properties of the working fluids or multi-phase flows , their regimes and states, and the geometrical characteristics of the system . But rigorous experimentation in thermodynamic systems can be particularly costly because environmental chambers are often needed to obtain reproducible results . Evidently, the lack of accessibility to an environmental chamber means that the ambient temperature is a random variable and the uncertainty of the experimental results can increase significantly  since the thermodynamic states of the thermodynamic cycle tend to fluctuate as a function of the heat rejected and the heat absorbed to and from the environment .
A less rigorous alternative to this approach is to perform a design of experiments at ambient conditions and execute the test multiple times to assess the error and response variations . Unfortunately, this route can also lead to an expensive and time-consuming solution on top of developing variable startup conditions that can affect the performance of the system. Consequently, Numerical simulations and theoretical analysis of thermodynamic cycles should help alleviate the cost of materials, manufacturing, and redundant testing.
Heat and time-dependent multi-phase computational solutions targeting problems with complicated geometries that cannot be simplified by axisymmetric conditions or reduction of dimensions tend to be challenging to develop, and the results are always suspicious if no experimental validation is provided. Therefore, a balance between experiments and simulations seems to be the ideal recipe to obtain cost-effective solutions with solid scientific ground.
In the last few decades, significant efforts have been focused on developing more rigorous models that account for turbulent fluctuations on mufti-phase flows [18, 19, 20] evidencing the complicated nature of such problems even for simple geometries  such as tubes and channels . However, no three-dimensional studies have been reported on the CFD multi-phase modeling of R404a condensation in a medium temperature water-cooled condenser at startup conditions. This research aims to provide a study of the unsteady condensing and evaporating characteristics of R404a in a water-cooled condenser by modeling the raising temperature and pressure in the system induced by the compressed vapor entering the condenser based on experimental data.
1.1 Numerical analysis of turbulent phase-transition in the condensation process of refrigerant flow
Multi-phase flow refers to the flow of a mixture of phases or species such as gases in liquids or liquids in gases of different densities  where the liquids, and gases are considered distinct phases. There are different mathematical characterizations targeted to study multi-phase flows. For instance, for dispersed flows, one phase consists of discrete, noncontinuous elements (such as bubbles) and the second phase is considered as a continuum. This approach is widely accepted for the analysis and characterization of oil and refrigerant relationships in vapor compressed cycles, given that at certain conditions, oil can scape the compressor and affect the system performance due to the variations in the thermal properties of the mixture induced by the oil .
In bubbly flows or other types of applications where dispersed multiphase flows are applicable, the properties of the inhomogeneous mixture are defined relative to the continuous phase . The mathematical manipulation for such flows can happen under a Lagrangian representation for the dispersed phase , and an Eulerian reference for the continuous phase. For the continuous phase, the material properties are defined at every point within a control volume , whereas in a separated flow under the Euler–Euler approach the phases are continuous but divided by a single contact surface .
This research focuses on the incompressible volume of fluid model where the mass transfer between phases is given by:
for evaporation and
for a condensation process. In the study of phase-transition of refrigerants, the saturation temperature is found based on the operating conditions of the system . In this study, the inlet and outlet temperatures and pressures were measured experimentally to compute
From a modeling perspective, the saturating temperature is a dominant factor in the prediction of the evaporating and condensing process. Therefore, the measured temperatures, pressures and the saturation temperature were loaded in the simulation as dynamics boundary conditions to drive the phase-transition process given the dynamic value of .
In order to find the quality of the mixture the volume fraction is defined as:
where the sum of volume occupied by the phases is 100%. Furthermore, the microscopic dynamics of the phase-transition process of refrigerants in heat exchangers and condensers can be studied by considering both the vapor and liquid as continuous and interpenetrating fluids under the Eulerian frame of reference . Even though this approach is computationally expensive given that the governing equations for each phase are solved together with the momentum equation, the mass diffusion can be monitored by a straightforward calculation. For this case, the microscopic behavior of the phases can be monitored by:
for which and are the densities of the liquid and vapor phases. Then. the mixture density and mixture viscosity are:
For each phase, the momentum conservation is given by:
where and are the momentum sources and the stress tensor is given by:
given that and are the shear and bulk viscosities of the phase Then, the two equation turbulent model for multi-phase flow is:
for the turbulent kinetic energy and
for turbulent dissipation where and
The thermal energy in the process is conserved and computed for each phase by:
Where the energy is and the enthalpy and effective thermal conductivity are:
2. Experimental procedure
Since the energetic state of a refrigerant far from its critical point can only exist for stable conditions at a single pressure and temperature pair , it is crucial to properly identify the coupled operating pressures and temperatures to obtain the saturated temperature of the refrigerant. For this reason, an experimental setup was designed and constructed so that the inlet and outlet conditions of the water and the refrigerant under transient conditions can be detected and used as dynamic boundary conditions for the numerical simulation.
The water-cooled condenser was fabricated using copper tubing. The dimensions for the line of water, refrigerant line and copper wall separation are shown in Table 1. To enhance the heat transfer, thermal insulation (Polyethylene foam) was placed around the water copper jacket and the rest of the tubbing connecting the components in the refrigeration system.
Figure 1 shows the configuration adopted to characterize the condensing details of R404a. The figure shows the standard cycle corresponding to a vapor compression cycle including the compressor (Copeland ZB15-KCE), the water-cooled condenser, the evaporator, and the thermostatic expansion valve (Danfoss T2). The temperature and pressure of the vapor and liquid refrigerant were measured at points (1) and (2) by attaching thermocouples to the refrigeration line and connecting pressure transducers to measure the pressure drop of the refrigerant across the helix. The water inlet and outlet temperatures are monitored (3) and (4) in a non-invasive fashion.
The test methodology consisted of 6 consecutive runs for 4.75 minutes each in order to stabilize the system and reach thermal equilibrium given that an environmental chamber was not used for this study. A flow meter (5) and pressure regulator are located at the water inlet to measure the flow rate and sustain the water pressure.
Table 2 shows the magnitude of the variables monitored to control the steady performance of the system.
Table 3 shows the sensors and instruments used for the experimental study as well as the accuracy reported by the manufacturers. The sampling rate was fixed to 100 ms using LabView 18 software and the corresponding hardware with the data acquisition system to gather the temperatures and pressures values of both water and refrigerant.
|Transducer (PX-309)||Full Scale|
The water and refrigerant temperature at the inlet and outlet of the condenser are reported in Figure 2 with the corresponding variation to document the test repeatability. The figure also shows the measured refrigerant pressure and the computed liquid saturation temperature using Coolproops 6.4.1 .
The experimental data demonstrated that the system reached maximum operating pressure close to 300 ms after starting the compressor, and a steady-state condition for the vapor refrigerant pressure was reached around 1200 ms after initialization. In order to characterize the startup conditions, the first 2000 ms of the average between the different experiments corresponding to each variable displayed in the Figure 2 were programed in ANSYS Fluent as dynamic boundary conditions.
3. Computational domain and boundary conditions
A time-dependent numerical model in three-dimension was developed to simulate the phase-transition and heat transfer of R404a in a water-cooled condenser. The geometry accounts for two fluid regions as well as one solid domain (copper) that separates both flows and provides the thermal mechanism for heat transfer. The computational domain ( contains three distant regions belonging to the space assigned as follows:
Fluid region for water .
Solid copper region for physical division of water and refrigerant .
Fluid region for refrigerant R404a .
The properties of the materials and working fluids are show in Table 4. The temperature dependent material properties were programmed in Fluent using transient data tables and user defined routines. A schematic representation of the domains describing the computational domain is shown in Figure 3. The figure shows the inlet and outlet faces assigned for the water and refrigerant flows. The 3D geometry was modeled in SpaceClaim from the real dimensions and the numerical solution was obtained using Ansys Fluent 19.2. The contacting surfaces were discretized with conformal elements and the computational grid was created using tetrahedral objects with a maximum Skewness of 0.74. The mesh contains 36040897 nodes and 119573669 elements with a mean size of 3.12e-2[m].
The residual target (RMS) for the momentum and energy equations were set to and respectively. The time step was fixed to for a total simulation time of 2 seconds.
Model validation and verification against experimental data is an important aspect of any numerical solution. For turbulent multi-phase flows, the model validation for the phase-transition process is particularly important because the solution of the volume fraction equation depends on the saturation temperature . However, the model validation and identification of the saturation temperature can be challenging if the geometry is complex enough to obstruct the path of light for particle image velocimetry and particle image thermometry or block the access of probes for electric tomography techniques . Unfortunately, the phase-transition process in the water-cooled condenser complies with the restrictions mentioned above because the refrigerant phase-transition takes place in a region where standard field experimental techniques offer limited information due to the water and copper jacket around the refrigerant line. Therefore, in order to develop a model that resembles reality, the saturation temperature has been programmed as a boundary condition based on the inlet and outlet temperatures and pressures measured experimentally.
3.1 Boundary conditions
The inlet boundary conditions for water and refrigerant are set to the mean mass flows. The flow rate of water measured experimentally going into the system ranges from to and the inlet refrigerant mass flow corresponds to the compressor mass flow provided by the manufacturer at the operating conditions.
The inlet and outlet time dependent temperature, pressures, and saturation temperature are loaded as user functions in Fluent from the experimental values shown in Figure 2. A no slip velocity condition was assigned to the boundary walls and a thermal insulation condition was assigned on the wall between the water and the environment. The outlet quality was set to 100% liquid given that the actual system was designed and fabricated with the charge and specifications to hold a constant sub-cool temperature after the condenser.
4. Results and discussion
The condensing behavior and evaporating dynamics of refrigerant R404a have been simulated using the VOF approach. Due to the temperature difference, and heat flux, the refrigerant condenses throughout the spiral tubing array, whereas sporadic evaporation occurs given the inertial and body forces present in the flow. The simulation was initialized by filling the total volume of the condenser section with both, the primary fluid (liquid water) and the secondary fluid domain (refrigerant) in a vapor state.
Figure 4 shows the domain initialization at t = 0.0 s. After this instant, the water starts to flow while the refrigerant is still in a full vapor state. The magnitude and direction of the water velocity is presented by the streak line in the internal flow section, and the quality of the secondary flow is indicated by the contour interface.
At t = 0.15 s, the refrigerant leaving the condenser has fully condensed near the outlet section based on the sub-cooled temperature boundary condition. Upstream the copper tube, the refrigerant begins a phase-transitioning throughout the lower section of the spiral geometry (shown in Figure 5) induced by the fluid temperature distribution.
The two-phase flow in the copper tube shows complicated details, with various phase-transition patterns and irregular oscillation manifesting at the periphery of the geometry due to the centrifugal  and inertial forces. The mean vapor quality found in the refrigerant region is 60%.
The water around the copper elbow carrying vapor refrigerant near the inlet develops a high velocity zone due to the reduction of area enhancing the heat flux in the system . On the other hand, the vapor refrigerant develops a vortex pair at the beginning of the first spiral. This is probably induced by the tube curvature close to the inlet elbow. Interestingly, no other vortices develop in subsequent curved regions which suggests that the vortex develops due to the thermal and inertial characteristics of the vapor refrigerant entering the tubular array and the proximity to the elbow located at the inlet.
Figure 6 shows the refrigerant phase transition details at 0.5 s. At this instant, the refrigerant has completely transitioned to liquid at the last quarter of the spiral leading to the outlet of the condenser. In the middle section of the condenser, the saturated vapor condenses and evaporates along the outer wall in repeated patters of saturated vapor and multi-phase flow swirling due to gravity and rotational speed.
Since the wall temperature across the copper wall was assumed constant, the phase-transition region extends across a large region of the tubular array. This boundary condition prevented a flash point transition between phases . Furthermore, the outlet temperature was specified to remain in a sub-cool state allowing the saturated liquid to develop further into the inner regions of the condenser.
Figure 7 shows the condensing and evaporating characteristics of R404a at t = 1.25 s. At the tube periphery and after the first revolution, the saturated vapor begins to condense, transitioning into a two-phase flow with a 90% quality and further decreasing into a quality of 82% while displaying a parabolic profile along the outer surface wall. Then, the mixture evaporates back to a saturated vapor state before it reaches its full condensation state after the third revolution.
Since the heat transfer can be enhanced or penalized by geometrical factors , the widely used Dittus-Boelter (Eq. (17)) was employed to find the Nusselt number in the internal tubing section for the vapor refrigerant flow. Based on the fundamental restrictions of the geometry, parametric modifications to the equation determining the Nusselt number (Eq. (18)) need to be considered based on constant wall temperature assumptions , curvature ration , and fluid regime  such that:
where the friction factor is given by:
given that the Prandtl and Reynolds numbers are:
Figure 8 shows the Nusselt number with respect to vapor refrigerant mass flow as a function different tube diameter. For , the Nusselt number calculated ranges between 900 and 1500 at different vapor mass flows. On the other hand, the Nusselt number at , ranges between 2250 and 4200, suggesting that the heat transfer enhancement occurs by increasing the diameter of the refrigerant line and reducing the number of spirals.
A time and temperature dependent three-dimensional multi-phase simulation with dynamics boundary conditions determined experimentally has been developed to characterize the thermal, inertial, and phase-transition details of refrigerant R404a in a water-cooled condenser at startup conditions. Given the complicated geometrical configuration of the system, where the secondary fluid (refrigerant) is concentric and flows inside the copper jacket in contact with the primary fluid (water), it becomes a challenging task to capture experimental information about the phase-transition details of the refrigerant using common optical or fluid field experimental methods such as optical tagging, particle image velocimetry, or tomography techniques. Therefore, indirect methods have been applied to obtain temperature and pressure values at the inlet and outlet sections of the condenser.
The numerical solution showed the condensing and evaporating oscillatory nature of the phase-transition throughout the spiral tubbing as well as the progression of saturated liquid at the outlet of the refrigerant line. Condensation occurs due to the temperature differences between the water and the refrigerant. Evaporation takes place given the body and inertial forces present in the system that prevent further bubble nucleation to sustain the condensation process.
The velocity profiles for both water and refrigerant show vortices developing at the inlet and outlet zones due to the curvatures and elbows present in the geometry. Nusselt numbers show that heat transfer can be optimized by slightly increasing the diameter of the refrigeration line without inducing a disruption in the flow regime or affecting the manufacturability of the part.
The author is grateful to FBD L.P. for the financial support to perform the experimental analysis.
Conflict of interest
The author declares no conflict of interest.
Vapor volume fraction
Mass transfer intensity
J. Bustamante, A. Rattner and S. Garimella, “Achieving near-water-cooled power plant performance with air-cooled condensers,” Applied Thermal Energy, vol. 105, pp. 362–371, 2016.
G. Heidinger, S. Nascimento, P. Gaspar and P. Silva, “Experimental evaluation of the thermal performance at different environmental conditions of a low temperature display case with built-in compressor and water-cooled condenser,” Applied Thermal Engineering, vol. 144, pp. 825–835, 2018.
T. Veldkamp, Y. Wada, J. Aerts, P. Doll, S. Gosling, L. Ja, Y. Masaki, T. Oki, S. Ostberg, Y. Pokhrel, Y. Satoh, H. Kim and P. Ward, “Water scarcity hotspots travel downstream due to human interventions in the 20th and 21st century,” Nat Commun, vol. 8, pp. 1–12, 2017.
C. Aprea, R. Mastrullo and C. Renno, “Experimental analysis of the scroll compressor performances varying its speed,” Applied Thermal Engineering, vol. 26, no. 10, pp. 983–992, 2006.
P. Donnellan, K. Cronin and E. Byrne, “Recycling waste heat energy using vapour absorption heat transformers: A review,” Renew. Sustain. Energy Rev, vol. 42, pp. 1290–1304, 2015.
A. Babiloni, J. Esbrí, C. Cervera, F. Molés and B. Peris, “Analysis based on EU Regulation No 517/2014 of new HFC/HFO mixtures as alternatives of high GWP refrigerants in refrigeration and HVAC systems,” International Journal of Refrigeration, vol. 52, pp. 21–31, 2015.
“System Drop-in Testing of R-410A Replacements in Split System Heat Pump,” Air-Conditioning, Heating, and Refrigeration Institute (AHRI) Low-GWP Alternative Refrigerants Evaluation Program (Low-GWP AREP),vol. Test report #22., 2013.
C. Acosta, D. Yanes, A. Bhalla, R. Guo, E. Finol and J. Frank, “Numerical and experimental study of the glass-transition temperature of a non-Newtonian fluid in a dynamic scraped surface heat exchanger,” International Journal of Heat and Mass Transfer, vol. 152, pp. 1–8, 2020.
C. Acosta, A. Bhalla and R. Guo, “Empirical and numerical determination of the freezing point depression of an unsteady flow in a scraped surface crystallizer,” Applied Thermal Engineering, vol. 179, pp. 1–8, 2020.
Y. Shang, . Y. Zhang, . Y. Hou, B. Bai and X. Zhong, “Effects of surface subcooling on the spreading dynamics of an impact water droplet,” Physics of Fluids, vol. 32, pp. 1–10, 2020.
S. Richter, S. Fleischer, M. Aritomi and R. Hamp, “Transient two-phase flow in arbitrary inclined tubes caused by depressurization of liquid with dissolved gasses,” International Journal of Heat and Mass Transfer, vol. 44, no. 1, pp. 1–15, 2001.
R. Akasaka, “Thermodynamic property models for the difluoromethane (R-32) + trans-1,3,3,3-tetrafluoropropene (R-1234ze(E)) and difluoromethane + 2,3,3,3-tetrafluoropropene (R-1234yf) mixtures,” Fluid Phase Equilibria, vol. 358, pp. 98–104, 2013.
D. Frenkel and B. Smit, Understanding molecular simulation: from algorithms to applications, San Diego: Academic Press, 2001.
T. He, C. Mei and J. Longtin, “Thermosyphon-assisted cooling system for refrigeration applications,” International Journal of Refrigeration, vol. 74, pp. 165–176, 2017.
J. Kim, C. Han and S. Jeong, “Disturbance observer-based robust control against model uncertainty and disturbance for a variable speed refrigeration system,” International Journal of Refrigeration, vol. 116, pp. 49–58, 2020.
C. Truesdell, The Tragicomedy of Classical Thermodynamics, Vienna: 1971, Springer.
Z. Wu and R. Du, “Design and experimental study of a miniature vapor compression refrigeration system for electronics cooling,” Applied Thermal Engineering, vol. 31, no. 2, pp. 385–390, 2011.
S. Mimouni, A. Foissac and J. Lavieville, “CFD modelling of wall steam condensation by a two-phase flow approach,” Nuclear Engineering and Design, vol. 241, no. 11, pp. 4445–4455, 2011.
D. Sun, J. Xu and Q. Chen,, “Modeling of the evaporation and condensation phasechange problems with FLUENT,” Heat Transf., Part B: Fundamentals, vol. 66, no. 4, pp. 326–342, 2014.
S. Mimouni, “CFD modeling of wall steam condensation: two-phase flow approach versus homogeneous flow approach,” Sci. Technol. Nucl. Install, 2011.
H. Mohammed, D. Giddings and G. Walker, “CFD multiphase modelling of the acetone condensation and evaporation process in a horizontal circular tube,” International Journal of Heat and Mass Transfer, vol. 134, pp. 1159–1170, 2019.
N. Padoin and C. Soares,, “CFD modeling of steam condensation in industrial pipes,,” Blucher Chem. Eng. Proc, vol. 1, no. 2, p. 12904–12911, 2015.
M. Konopacki, M. Kordas, K. Fijalkowski and R. Rakocz, “Computational fluid dynamics and experimental studies of a new mixing element in a static mixer as a heat exchanger,” Chem. Process Eng, vol. 36, no. 1, pp. 59–72, 2015.
G. Anjos, G. Oliveira and M. Norberto, “Arbitrary Lagrangian–Eulerian Method for Two-Phase Flows: Applications,” in Encyclopedia of Two-Phase Heat Transfer and Flow III, World Scientific, 2018, pp. 141–183.
J. Parekh and R. Rzehak, “Euler–Euler multiphase CFD-simulation with full Reynolds stress model and anisotropic bubble-induced turbulence,” International Journal of Multiphase Flow, vol. 99, pp. 231–245, 2018.
M. Göz, S. Lain and M. Sommerfeld, “Study of the numerical instabilities in Lagrangian tracking of bubbles and particles in two-phase flow,” Comput. Chem. Eng, vol. 28, p. 2727–2733, 2004.
E. E. Michaelides, “Review—The Transient Equation of Motion for Particles, Bubbles, and Droplets,” J. Fluids Eng., vol. 119, no. 2, pp. 1–16, 1997.
H. Hu, “Direct simulation of flows of solid-liquid mixtures,” International Journal of Multiphase Flow, vol. 22, no. 2, pp. 335–352, 1996.
A. Arora and S. Kaushik, “Theoretical analysis of a vapour compression refrigeration system with R502, R404A and R507A,” International Journal of Refrigeration, vol. 31, no. 6, pp. 998–1005, 2008.
M. Zareh, H. Shokouhmand, M. Salimpour and M. Taeibi, “Numerical simulation and experimental analysis of refrigerants flow through adiabatic helical capillary tube,” International Journal of Refrigeration, vol. 38, pp. 299–309, 2014.
S. Lee, M. Kim and S. Ro, “Thermal Conductivity of 1,1,1-Trifluoroethane (R143a) and R404A in the Liquid Phase,” Journal of Chemical and Engineering Data, vol. 46, no. 5, pp. 1013–1015, 2011.
“Pure and Pseudo-pure Fluid Thermophysical Property Evaluation and the Open-Source Thermophysical Property Library CoolProp,” Industrial and Engineering Chemistry Research,vol. 53, no. 6, pp. 2498--2508, 2014.
D. Green and R. Perry, Perry’s Chemical Engineers’ Handbook, London: McGraw-Hill, 2007.
B. Fadhl, L. Wrobel and H. Jouhara, “CFD modelling of a two-phase closed thermosyphon charged with R134a and R404a,” Applied Thermal Engineering, vol. 78, no. 5, pp. 482–490, 2015.
W. Callister and D. Rethwisch, Materials sicnce and engineering an introduction, Wiley, 2010.
V. Geller, D. Bivens and A. Yokozeki, “Viscosity of Mixed Refrigerants, R404A, R407C, R410A, and R507C,” in Int. Ref. and Air cond., West lafayette, 2000.
E. Lemmon, M. Huber and M. McLinden, NIST standard reference database 23: reference fluid thermodynamic and transport properties-REFPROP, version 8.0, Gaithersburg: Standard Reference Data, 2007.
A. Gadelha, S. Neto, R. Swarnakar and A. Lima, “Thermo-Hydrodynamics of Core-Annular Flow of Water, Heavy Oil and Air Using CFX,” Advances in Chemical Engineering and Science, vol. 3, no. 4A, pp. 37–45, 2013.
R. Heide, “The surface tension of HFC refrigerants and mixtures,” nt J. Refrig, vol. 20, pp. 496–503, 1997.
W. Lee, A pressure iteration scheme for two-phase flow modeling. Multiphase Transport Fundamentals, Reactor Safety, Applications, Washington, DC: Hemisphere Publishing, 1980.
C. Acosta, A. Bhalla and R. Guo, “Phase-transition temperature determination using optical spectroscopy in a rotating flow inside a scrape surface crystallizer,” in Photonic Fiber and Crystal Devices: Advances in Materials and Innovations in Device Applications XIV, San Diego, 2020.
H. Wu, X. Peng, P. Ye and Y. Gong, “Simulation of refrigerant flow boiling in serpentine tubes,” Transfer, International Journal of Heat and Mass, vol. 50, no. 5, pp. 1186–1195, 2007.
H. Ito, “Flow in Curved Pipes,” JSME international journal, vol. 30, pp. 543–551, 1987.
A. Seixlack and M. Barbazelli, “Numerical analysis of refrigerant flow along non-adiabatic capillary tubes using a two-fluid model,” Applied Thermal Engineering, vol. 29, no. 2, pp. 523–531, 2009.
J. Abraham, E. Sparrow and W. Minkowycz, “Internal-flow Nusselt numbers for the low-Reynolds-number end of the laminar-to-turbulent transition regime,” International Journal of Heat and Mass Transfer, vol. 54, no. 1, pp. 584–588, 2011.
R. Seban and E. McLaughlin, “Heat transfer in tube coils with laminar and turbulent flow,” International Journal of Heat and Mass Transfer, vol. 6, no. 5, pp. 387–395, 1963.
H. Zhao, X. Li, Y. Wu and X. Wu, “Friction factor and Nusselt number correlations for forced convection in helical tubes,” International Journal of Heat and Mass Transfer, vol. 155, p. 119759, 2020.
L. Janssen and C. Hoogendoorn, “Laminar convective heat transfer in helical coiled tubes,” International Journal of Heat and Mass Transfer, vol. 21, no. 9, pp. 1197–1206, 1978.