1. Introduction
The issues of heat and mass flow are of great importance in current technological development [
1,
2]. Two-phase flows with heat exchange play an important role here. The flow structures formed during the condensation process have a significant impact on the intensity of the condensation process. It is influenced by void fraction
φ, vapor quality
x, the mass flux density
G, and the shape of the interface. For this reason, in addition to the study of heat transfer or flow resistance during the condensation process in minichannels, the observation of two-phase flow structures formed in a nonadiabatic flow is of great importance. In a nonadiabatic flow, it is very important to define the range of the above-mentioned parameters in which individual flow structures occur. The structures of two-phase adiabatic flow; e.g., water and air, are the most widely described in the literature. However, in heating and refrigeration devices, two-phase nonadiabatic flows, such as boiling or condensation, are most often used. Due to the difficulty level of the condensation process and its observation in minichannels, studies of the flow structures in the boiling process are conducted and described much more frequently. The investigations of Elazhary and Soliman [
3] are an example of research on adiabatic flow structures. The authors conducted experimental investigations of two-phase flow regimes and pressure drop in horizontal minichannels with a rectangular cross-section and dimensions of 1.87 × 20 mm
2. They investigated the flow of air–water mixtures at 200 kPa pressure and standard temperature. The authors identified bubbly, plug, churn, and annular flow structures. The studies were conducted in gas and liquid phases with superficial velocities ranging from 0.04 ≤
jv ≤ 10 m/s to 0.02 ≤
jl ≤ 0.7 m/s. The most extensive studies regarding condensation processes in minichannels were carried out by Coleman and Garimella [
4]. They investigated the condensation of R134a refrigerant in nine minichannels of various cross-sections (round, square, and rectangular) with a hydraulic diameter
dh = 1–4.91 mm and mass flux density in the range
G = 150–50 kg/m
2s. The channels were cooled by air. The flow was visualized during the condensation process, and based on these findings, the flow structures and substructures were classified. Four structures were predominantly observed: annular, wave, intermittent, and dispersed. The influence of gravity, interfacial interactions, and inertial force was also studied. It has been shown that the changes in vapor quality and mass flow influence the flow structure. The impact of the channel size on the existing structures has also been demonstrated [
5]. On the other hand, Jige et al. [
6] investigated the condensation of R32 refrigerant in rectangular multiports with a hydraulic diameter
dh = 0.5–1 mm. The results demonstrated that the plug structure occurred at low vapor quality, while the annular structure occurred at high vapor quality and high mass flux density. Xiao and Hrnjak [
7] conducted research on the condensation of R134a, R1234ze, R32, R245fa, and R1233zd refrigerants in channels with a hydraulic diameter
dh = 1.4–6 mm. The transition between various structures formed during the process of condensation was described. The annular flow is mainly characterized by the presence of surface tension and shear forces. As the thickness of the condensate film increases, the gravity force obtains a value higher than the surface tension and shear forces. At this point, the flow changes to stratified flow. When the thickness of the liquid film approaches the radius of the channel, the ratio of gravity force to the surface tension can be given as the Bond number Bo:
where
ρv and
ρl are the densities of the gas and liquid phases, respectively;
dh is the hydraulic diameter;
g is the gravity acceleration; and
σ is the surface tension factor. The relationship between inertia forces and surface tension can be described by the Weber number (
We):
The Weber number is used for calculations in two-phase flows, with particular emphasis on disturbances at the interface. This number describes the ratio of inertia forces to surface tension. The transition between stratified and wave structures is caused by the formation of Kelvin–Helmholtz instability, and the height of the waves depends on the differences in the velocity of the liquid and gas phases. For a transition to an intermittent flow, the wave amplitude should be close to the internal diameter of the channel (or otherwise should reach the highest value). At low values of the mass flux density
G, the amplitude of the waves is low, which prevents the formation of an intermittent flow despite an increase in the thickness of the liquid film. Enoki et al. [
8] observed the flow structures of two-phase R410A refrigerant in circular, rectangular, and triangular cross-sectional channels with a hydraulic diameter
dh = 1 mm. The study was carried out in the field with a mass flux density
G = 30–400 kg/m
2s, vapor quality
x = 0.05–0.9, and saturation temperature
ts = 10 °C. The results revealed the formation of slug, stratified, wave, froth, and annular flow structures. On the basis of these results, the authors made modifications to the flow structure map developed by Chen et al. [
9], which was created on the basis of a two-phase flow test of the R134a refrigerant in pipe channels with hydraulic diameter
dh = 1.10, 2.01, 2.88, and 4.26 mm [
6]. Nema et al. [
10] proposed new transition boundaries for individual flow structures, using the database for the condensation of R134a refrigerant in the channels with a hydraulic diameter
dh = 1–4.91 mm and mass flow density range
G = 150–750 kg/(m
2s). Thus, they proved the influence of the channel diameter and heat transfer on the occurrence of some flow structures, by comparing it with adiabatic, air–water, two-phase flow under similar conditions. In the case of adiabatic flow, the occurrence of “discrete wave” and “dispersed wave” structures was not observed, and the boundaries of individual structures differ significantly in both cases. It is similar to the influence of the channel diameter on the range of occurrence of particular structures. The decrease in the diameter of the channel causes the shifting of most of the boundaries towards a higher vapor quality
x. This shift is in turn affected by an increase in the influence of viscosity
μ and surface tension
σ on the condensation process, following a decrease in the diameter of the channel and the influence of gravitational forces. Nema et al. [
11] used the Bond number (
Bo) to determine the formation of the wave structure during the flow. When the Bond number exceeds the critical value (
Bocr), there is a possibility for the occurrence of a wave structure. Otherwise, there is a direct transition from the annular to an intermittent structure. The critical Bond number is determined from the following equation:
If the Bo number is lower than
Bocr, the minimum amount of liquid phase required for the formation of gas slugs is described by the Lockhart–Martinelli parameter
χtt,slug:
The transition between a discrete wave and dispersed wave structure is described by the modified Froude number
Fr*, which determines the mutual relationship between inertia force and gravity in the flow:
Zhuang et al. [
12] performed a study on the two-phase flow of R170 refrigerant in a horizontal, smooth glass channel with a hydraulic diameter
dh = 4 mm. During the investigation, structures such as slug, transitional, annular wave, and annular were observed. The experimental results were compared with the maps constructed by Breber, Tandon, Cavallini, El Hajal, and Barbieri, and then the authors modified the map developed by Kim et al. [
13]. They also determined the limits of transition between individual flow structures using the modified Weber number (
We*) and the Lockhart–Martinelli parameter:
Annular/annular wave:
Annular wave/transitional: ;
Transitional/slug: ;
Slug/plug: .
In the process of identifying two-phase flow structures, flow structure maps are a very useful tool. Nasrfard et al. [
14] analyzed the efficacy of flow maps proposed by El Hajal et al., Kim et al., Tandon et al., and Cavalini et al. by comparing them with the results of their own experimental studies on the condensation of R141b refrigerant in a horizontal, smooth circular channel with a hydraulic diameter
dh = 8 mm. It was found that the previously described maps of flow structures did not match with the results of visualization. Therefore, the authors proposed their own map, which was a modification of the El Hajal et al. map that functions in the low range of the vapor quality values. Xiao and Hrnjak [
7] proposed a map of the flow structures formed during the condensation of R134a, R1234ze, R32, R245fa, and R1233zd refrigerants in the channels with hydraulic diameter
dh = 1.4 and 6 mm. Specific enthalpy of the refrigerant was adopted as one of the parameters that define the limits of the occurrence of individual flow structures. This procedure took into account the heat transfer that occurred during the process of condensation. The map thus constructed was also a modification of El Hajal et al.’s map.
In the literature, there are many adiabatic two-phase flow studies. Some authors published a flow structures map for adiabatic and nonadiabatic two-phase flows. There is a lack of visualization studies of condensation of new, environmentally friendly refrigerants with low ODP (Ozone Depletion Potential) and GWP (Global Warming Potential). The condensation process of these refrigerants in minichannels has not been tested in a wide range. Flow structure maps for the condensation process of new refrigerants are lacking in the literature. The novelty of this work is represented by the investigation results of the condensation flow structures of new low-pressure refrigerants from 3M’s Novec group in minichannels. The study was conducted over a very wide range of condensation processes in minichannels. So far, only Al-Zaidi et al. [
15] investigated the flow structures formed during the condensation of HFE7100 refrigerant in a multiport made of rectangular minichannels with
dh = 0.57 mm,
ts = 60 °C, and
G = 48–126 kg/m
2s. Mikielewicz et al. [
16] carried out investigations on HFE 7000 and HFE 7100 refrigerant condensation in a minichannel with an internal diameter
dh = 2.3 mm for the mass flux density range
G = 240–850 kg/m
2s, but the authors did not investigate the flow structures during condensation of these refrigerants. The Novec refrigerants are used in wide range of heating devices. Woodoc et al. [
11] investigated the use of HFE7000 refrigerant in a silicone mini heat exchanger (MECH-X), which is a heat sink reactor with a thickness of 800 μm. The exchanger was used to cool electronic components that were heated to 90 °C. The boiling studies of Novec refrigerants (HFE7000, HFE7100, FC-72, and HFE-649) were conducted by Eraghubia et al. [
17], Piasecka et al. [
18,
19,
20], and Cao et al. [
21], who created a visualization of pool boiling of HFE-649. Mohamadi et al. [
22], on the other hand, investigated heat transfer with the use of a nanofluid that consisted of an HFE 7000 refrigerant and Al
2O
3 as well as SiO
2 nanoparticles. Adebayo et al. [
23] proposed the use of an HFE 7000 refrigerant in combination with CO
2 in a cascade refrigeration device. Other scientists proposed the use of Novec refrigerant in the organic Rankine cycle [
24] and a three-stage refrigeration system [
25]. Currently, two-phase flow modeling also is very often used to reduce costly experimental studies in the field of flow structures [
26].
The next novelty of this article is represented by flow structure maps, developed based on visualization and flow investigations. So far, no results of flow structures during the process of condensation of this type refrigerants in minichannels have been published. The proposed flow structure maps are also new in this field.
3. Image Analysis
The term “image analysis” refers to the automatic processing and analysis of images of selected objects or surroundings, and the primary objective is to obtain useful mathematical information that can be used for their interpretation that may further affect the results of the observed process. The image analysis algorithms process the images in three stages. The first stage involves preprocessing of the video signal. Its purpose is to eliminate disturbances, extract the examined object from the background, determine the level of disturbances, and balance the histogram. The second stage includes image segmentation, location of the examined objects, recognition of their shape, distinguishing the characteristic features of the object, etc. The final stage comprises analysis of the movement of the examined object, if any; interpretation of the object control; assigning processing parameters; and analyzing the image [
30].
Visualization studies of adiabatic, two-phase, gas–liquid flow in vertical channels have been carried out by many authors [
31,
32]. Quantitative stereology technique was used in this work. Stereology is a scientific discipline in which information about the three-dimensional structure of objects can be obtained from their two-dimensional counterparts (mainly images). It provides information about volume (
V), surface area of the object (
A), and length of the line (
L). The application of stereology methods requires the inclusion of probability calculations and mathematical statistics because the studied structures are sometimes stochastic, as their sizes and shapes change over time and space (especially during boiling and condensation processes). In stereology, one dimension is reduced at the moment it is mapped on the solid plane. Reduced quantities indicate the volume that the “real structure” occupies. Measurements made on the basis of the structural image can be related to the surface (
AA,
LA,
PA) or a line drawn on the image (
LL,
PL). The latter case is called the Cavaliere–Hacquert linear method. Among other applications, this method is used for the determination of the void fraction (
φ) of the channel, which is defined as the volume fraction of the dispersed phase per unit volume of a two-phase mixture. In this case, the void fraction (
φ) is equal to the ratio of the length of the section passing through the gas phase structure (
Lki) in the two-dimensional image to the total length of the section outlined in the image. The section
Lk is a straight line that connects opposite walls of the channel in a two-dimensional image. It is influenced by the most characteristic flow elements inside the channel. The operation can be repeated several times to obtain the average value of the void fraction for a given image.
In investigations involving two-phase flow structures, different methods of image analysis are used to determine the geometric parameters of the flow structure, velocity of individual phases, vapor quality (
x), and void fraction (
φ) [
33]. The void fraction (
φ) is defined as the volume fraction of the gas phase with respect to the total volume. According to the rules of stereology, it can be determined by reducing one dimension of the three-dimensional image; that is, by reducing the volume of individual phases to the cross-sectional area. In this case, the static void fraction can be calculated as follows:
where
Av is the area of the image occupied by the vapor phase and
Al is the area of the image occupied by the liquid phase. If the process takes place in a circular minichannel, then the two-phase flow structures exhibit circular symmetry. Therefore, it can be assumed that the fraction of the gas phase in the total volume, determined for the two-dimensional image, is close to the real fraction [
34]. By knowing the static void fraction, it is possible to determine the static vapor quality (for example, in the case of condensation or boiling process):
where
ρl and
ρv are the densities of liquid and vapor phases, respectively; and
φ is the static void fraction. This formula can be used during flow without interphase slip [
35,
36].
This method was used for our investigations. In image analyses, frames were generated from the recorded images and processed. The image analysis used in this study was based on the algorithm presented by Michalska-Pożoga et al. [
37], which was written in MATLAB (version, manufacturer, city, country) and modified by the author and adapted to determine the void fraction
φ. The protocol followed for image processing and analysis is presented in
Figure 4. Based on the results of this analysis, vapor quality x was determined using the dependencies defined by Equation (8).
To apply the above algorithm, it was necessary to capture the correct image of the tested flow structure to enable the algorithm to binarize the image.
Figure 4 shows the individual phases of image processing using the described algorithm [
36]. The first step was to cut an image of the inside of the minichannel from the test frame (that is, remove the background). This was necessary because the visualization studies were carried out inside the minichannel. The next step was the possible closing of the phase separation line and binarization of the image (
Figure 4c). The black pixels now occupied the fields filled with the liquid phase, and the white pixels occupied the gas phase. Then, the algorithm inverted the image (
Figure 4d) and reduced the so-called noise (
Figure 4e). During inversion, the color of the pixels changed from white to black and vice versa, which meant that after inversion, the liquid phase was represented by white pixels and the gas phase by black pixels. Only such a prepared image could be analyzed using the MATLAB algorithm. The algorithm also allowed the author to align shapes and change colors, contrast, and brightness, as well as the filter. Based on the given dimensions of the image and the number of black and white pixels counted by the algorithm, it was possible to determine the area of the image occupied by individual phases. The ratio of the black pixel area (after image conversion and filtering by the algorithm) to the total area of the image described the void fraction
ϕ. Based on this value, the vapor quality
x was determined. A very similar method in the study of the flow structures during boiling was used by Płaczkowski et al. [
38].
The accuracy of the determination of the void fraction and vapor quality using stereology technique depends on observed flow structures, the size of minichannel, the flow rate of the refrigerant, and the image recording rate. Accuracy varies from about 2 to 20%. The best accuracy of stereology methods is for plug and annular flows; the mist and transition structures, where froth is observed, have an accuracy of about 20%. The larger diameter of the minichannel gives a higher accuracy when determining the area occupied by individual phases. Accuracy also depends on the image quality.
4. Experimental Visualization Results
The visualization studies of the two-phase flow structures formed during the condensation of HFE7000, HFE7100, and Novec649 refrigerants were performed on the test stand shown in
Figure 1. This article presents the results characteristic of this study. Vapor quality and void fraction were determined using a picture analysis system described by Bohdal and Sikora [
36].
Figure 5 and
Figure 6 present the structures formed during the visualization of the condensation process of HFE700 refrigerant in minichannels with an internal diameter
dh = 2.0 and 0.8 mm. For this refrigerant, structures such as mist, annular, annular-wavy, slug, plug, and bubble were produced.
Figure 7 shows the results of the visualization of the condensation structures of HFE7100 refrigerant in a minichannel with a diameter
dh = 2.0 mm. In this case, the following structures were observed: mist, annular, annular-wavy, slug, plug, and bubble. Exemplary flow structures, observed during the study of the condensation of the HFE7100 refrigerant in a minichannel with a diameter of 0.8 mm, are shown in
Figure 8.
Figure 9 and
Figure 10 show pictures of flow structures observed during Novec649 refrigerant condensation in minichannels with internal diameter
dh = 1.2 and 0.5 mm. Similar two-phase flow structures were observed with two diameters; however, during condensation in the channel with
dh = 1.2 mm, no plug structure was observed. The descriptions given in
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9 and
Figure 10 specify the parameters used during the two-phase condensation flow, which allowed for drawing conclusions regarding the occurrence of individual structures.
In
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9 and
Figure 10, it can be seen that the type of flow structure depended on the flow parameters, such as mass flux density
G, velocity of individual phases, vapor quality
x, and void fraction
ϕ. With the increase of the mass flux density, the range of occurrence of stratified structures decreased in favor of discontinuous ones. The plug and slug flows were included in the group of slug structures. They differed from each other in shape and thermal flow parameters. The slug flow was particularly observed for high mass flux densities
G, the value of which depended on the type of refrigerant tested and vapor quality. It was related to the disturbances on the interface. Below
G = 500 kg/m
2s, there was practically no slug flow in the tested conditions, but above this value, plug flow disappeared in favor of slug flow. The plug flow occurred in the form of a regularly shaped slug, and was mainly observed in small-diameter channels, especially at
dh ˂ 1 mm. Under these conditions, an increase in the influence of surface tension on the condensation process became apparent. At
G > 3000 kg/m
2s, no bubbly flow was observed.
6. Maps of Two-Phase Flow during Condensation of Novec Refrigerants
Based on the results obtained after visualization studies on the condensation of HFE7000, HFE7100, and Novec649 refrigerants in a minichannel, the author constructed maps of flow structures. Owing to the fact that the HFE7100 refrigerant had a much higher vapor phase density compared to the HFE7000 and Novec649 refrigerants, significant difficulties were encountered during the construction of a common map of flow structures in
G–x coordinates for the three tested refrigerants. In an attempt to generalize the results of experimental studies obtained for these three refrigerants, we proposed to present an experimental map of two-phase flow structures formed during condensation in minichannels with an internal diameter
dh = 2.0–0.5 mm, in the coordinates
G–
ϕ.
Figure 17 and
Figure 18 present original maps of the flow structures created as a result of visualization investigations, and the construction of maps was supported by image analysis.
Figure 18 shows the influence of the channel diameter on the occurrence of individual flow structures. This was directly related to the mass flux density, which increased as the channel diameter decreased at a constant flow rate. Based on the presented flow structure map, it was concluded that the annular, plug, and bubbly substructures occurred mainly at lower mass flux density (
G) values. As can be seen in
Figure 17 and
Figure 18, at the mass flux density
G = 2500 kg/m
2s, mist flow occurred with void fraction
φ > 0.85, annular flow at 0.72 <
φ < 0.85, annular wave flow at 0.55 <
φ < 0.72, slug flow at 0.35 <
φ < 0.55, and bubbly flow at
φ < 0.35. In this range of mass flux densities, no occurrence of plug flow was observed. At high values of mass flux density (both phases at high velocities), the probability for the formation of mist, annular wave, and slug flow structures was widened, and the remaining observed substructures disappeared. With an increase in the velocity of the gas phase and the amount of condensate (with a decrease in the void fraction
φ), the annular substructure transformed to the annular wave substructure. It was similar to intermittent structures. With the increase in the difference of velocity between the vapor and liquid phases, the plug substructure disappeared and underwent transition to the slug substructure. The plug substructure exhibited regular shapes, which resulted from the low value of the mass flux density (
G), which in turn could be attributed to the low velocity of both phases and the relatively small difference in their velocities (slip).
The results of the present study were compared with the flow structure maps proposed by other authors. It was observed that the maps presented in
Figure 19,
Figure 20 and
Figure 21 had different ranges of mass flux density (
G). Unfortunately, the ranges of the presented maps of flow structures only partially corresponded to the scope of the presented research. The map proposed by El Hajal [
40] (
Figure 19) has been verified for its possible application in two-phase flows in minichannels during the boiling process and in the presence of high-pressure refrigerants. In the case of condensation of the tested low-pressure refrigerants in the minichannels, the boundaries of the occurrence of individual structures were shifted towards lower values of the vapor quality.
Figure 20 depicts the map of two-phase flow structures proposed by Soliman [
41]. It was created for the two-phase, nonadiabatic flow of R12 and R113 refrigerants in conventional channels. The lack of compliance also occurred when comparing this study’s experimental results with the map of two-phase flow structures developed by Olivier et al. [
42] (
Figure 21). They created the map for the condensation of high- and medium-pressure refrigerants in finned minichannels. It can be seen that this map does not take into account the occurrence of intermittent structures at low values of mass flux density (up to 300 kg/m
2s). This was related to an increase in surface roughness due to microribbing. Such a phenomenon was not observed in the current study. Thus, the discrepancy in the results of experimental studies can be attributed to differences in both physicochemical properties of the tested factors and the dimensions of the channel.