Next Article in Journal
Short-Term Electric Load Forecasting Based on Variational Mode Decomposition and Grey Wolf Optimization
Next Article in Special Issue
An Experimental Performance Evaluation for a Swirl-Vane Separator Using an Air-Water Test Facility
Previous Article in Journal
Investigating the Impact of Public Transport Service Disruptions upon Passenger Travel Behaviour—Results from Krakow City
Previous Article in Special Issue
Evaporation of Methanol Solution for a Methanol Steam Reforming System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics of Refrigerant Boiling Heat Transfer in Rectangular Mini-Channels during Various Flow Orientations

Faculty of Mechatronics and Mechanical Engineering, Kielce University of Technology, 25-314 Kielce, Poland
*
Author to whom correspondence should be addressed.
Energies 2021, 14(16), 4891; https://doi.org/10.3390/en14164891
Submission received: 11 July 2021 / Revised: 30 July 2021 / Accepted: 7 August 2021 / Published: 10 August 2021
(This article belongs to the Special Issue Two-Phase Flow with Boiling Heat Transfer)

Abstract

:
This paper reports the results of heat transfer during refrigerant flow in rectangular mini-channels at stationary conditions. The impacts of selected parameters on boiling are discussed, i.e., thermal and flow parameters, dimensions and orientation of the channels. Four refrigerants (FC-72, HFE-649, HFE-7000 and HFE-7100) were used as the working fluid. Research was carried out on the experimental set-up with the test section with a single rectangular mini-channel of 180 mm long and with a group of five parallel mini-channels, each 32 mm long. The temperature of the mini-channel’s heated wall was measured by infrared thermography. Local values of the heat transfer coefficient at the contact surface between the fluid and the plate were calculated using the 1D mathematical method. The results are presented as the relationship between the heat transfer coefficient and the distance along the mini-channel length and boiling curves. Two-phase flow patterns are shown. Moreover, the results concerning various refrigerants and the use of modified heater surfaces are discussed. The main factors influencing the heat transfer process were: mini-channel inclination to the horizontal pane (the highest heat transfer coefficient at 270° and 0°), using modified heater surfaces (especially electroerosion texturing and vibration-assisted laser No. 2 texturing) and working fluids (FC-72 and HFE-7000).

1. Introduction

The progressive development of miniaturization and the rapid growth in the areas of technology application have resulted in the need to dissipate increasing amounts of heat from ever smaller surfaces, and therefore it becomes necessary to look for new ways of heat dissipation. One of the effective methods of cooling, enabling the dissipation of high-density heat in a wide temperature range, is boiling heat transfer during flow in mini-channels with the change of the state of aggregation of the working medium. This process is characterized by obtaining high values of the heat transfer coefficient and high heat flux, with a small temperature difference between the heating wall and the saturated liquid. The results of investigations reported in the literature confirm the importance of various factors’ influence on boiling heat transfer in mini-channels, such as thermal and flow parameters (temperature and pressure of the fluid, mass flux) and the spatial orientation of mini-channels and their geometry. Modifying the heated wall surfaces of mini-channels could intensify the boiling heat transfer.
In the literature, the influence of various parameters on flow boiling heat transfer has been widely discussed for over a dozen years, however the results relate to different geometries and spatial orientations of the mini-gaps, and various fluids are still investigated.
In Reference [1], an experimental investigation concerning boiling heat transfer during refrigerant flow (R134a, R245fa and R600a) along small-diameter tubes (ranging from 0.38 to 2.6 mm) at a wide range of mass and heat fluxes (mass flux from 49 to 2200 kg/(m2·s) and heat flux up to 185 kW/m2) were discussed. The effect of the influence of selected experimental parameters, such as geometrical (diameter of a tube) and thermal and flow (mass flux, heat flux, saturation temperature), on the heat transfer coefficient and dry-out vapor quality was examined, for several fluids. It was concluded that the heat transfer coefficient increased with increasing mass flux, heat flux and saturation temperature, while it decreased with increasing tube diameter. The dry-out vapor quality decreased with increasing mass flux and vapor-specific volume. The results and comments were evaluated only for circular mini-channels.
The experimental research in [2] concerned flow boiling heat transfer of refrigerants R134A and R410A flowing along a smooth tube or two surface-enhanced tubes with a 11.5 mm inside diameter. The experimental tests were conducted with the following parameters: saturation temperatures of 6 and 10 °C, mass flux from 70 to 200 kg/(m2∙s), heat flux from 10 to 35 kW/m2 and inlet and outlet vapor quality of 0.2 and 0.8, respectively. The boiling heat transfer coefficient of the enhanced surface tubes was 1.15–1.66 times higher than that of the smooth tubes. Coefficients were evaluated using several widely used correlations. A good predictive ability was obtained when only two correlations were used.
Investigations on flow boiling heat transfer of R32 and R1234yf, pressure drop and flow patterns in a horizontal square mini-channel with a hydraulic diameter of 2 mm were described in [3]. The authors showed the effects of mass flux, vapor quality, heat flux and refrigerant properties on the flow boiling characteristics. The tested parameters were in the ranges: mass flux: 50–400 kg/(m2∙s), heat flux: 5–40 kW/m2 and vapor quality: 0.1–0.9, at a saturation temperature of 15 °C. Detected flow patterns were classified as plug, wavy, churn and annular flows. The data obtained using R32 and R1234yf fluids were analyzed. It was observed that the heat transfer coefficient for R32 achieved higher values in comparison with R1234yf at the same mass flux and vapor quality. Moreover, the frictional pressure drop of R32 was 1.5–2 times higher in comparison to R1234yf at the same mass flux and vapor quality. Additionally, the small effect of channel shape on the frictional pressure drop was observed. The authors also provided a theoretical explanation for these observed relationships.
The aim of the study described in [4] was to investigate the effects of parameters such as mass flux, heat flux, channel diameter and inlet pressure on boiling during FC-72 flow in a semicircular mini-channel. The ranges of the experimental parameters which were tested are as follows: mass flux: 300–2000 kg/(m2∙s), heat flux: 100–300 kW/m2 and diameter of the channel: 1–5 mm. The results indicated that an increase in the mass flux caused a change in the entropy generation when the flow changed from laminar to transitional flow in small-diameter channels. The increase of the channel diameter resulted in a lower variation in saturated temperature and the reduction of the heat transfer coefficient since the vapor quality achieved high values close to the intermittent dry-out range. The larger channel diameters produced higher entropy generation in comparison to the smaller diameters for each heat flux. During the increase in the inlet pressure, the entropy generation decreased for all tested mass and heat fluxes. Comparative results for fluids with different physical properties than Fluorinert FC-72 would be interesting and allow for the formulation of broader observations and remarks.
In Reference [5], boiling heat transfer during R141b flow in bi-porous mini-channels with dimensions of 1.5 × 1 mm was investigated. The heat transfer performance, pressure drop, two-phase flow structures and flow instability at three mass fluxes (158, 253 and 348 kg/(m2∙s)) were discussed. It was noticed that the heat transfer coefficient of bi-porous mini-channels rose faster with heat flux and was larger than that of the plain surface mini-channel. The pressure drop obtained for the enhanced surface mini-channel increased significantly at high heat flux compared to that of the plain surface mini-channel. Furthermore, it was observed that the flow instability of the bi-porous mini-channel was smaller than that of the plain surface at different mass fluxes. Novel bi-porous mini-channels sintered with copper woven tape, which is proposed to enhance the flow boiling performance, need more experimental testing. The analysis was based on the experimental research in which only one working fluid, i.e., R141b, was used. This fluid characterizes a low boiling point of 32.1 °C under normal pressure conditions.
In Reference [6], the experimental investigations of two kinds of model plate heat exchangers, namely a mini-gap exchanger and a mini-channel heat exchanger, were described. A single plate was equipped with a mini-gap of 0.1 m long and 0.2 m wide with a thickness of 0.5 and 1.0 mm or with a set of 50 rectangular mini-channels with a 1 × 1 mm cross-section. Ethanol was used as the working fluid. Experiments were conducted in the following conditions: heat flux supplied to the evaporator from 4.1 to 28 kW/m2 and mass flux from 10 to 92 kg/(m2·s). Values of the heat transfer coefficient and pressure drop were compared for both heat exchangers as a function of various parameters. It was reported that the coefficient increased with the decreasing thickness of the mini-gap. Furthermore, the coefficient was higher for a mini-channel multiport than for a mini-gap of the same cross-section area. It was indicated that the pressure drop at higher mass flow rate slightly increased for a mini-channel heat exchanger in comparison to a mini-gap heat exchanger. The modified mini-gap surface that was used did not result in a significant heat transfer coefficient increase. The results obtained for such different geometries of heat exchangers with mini-spaces, as presented in this work, in comparative analyses require further study in a wide range of thermal and flow parameters. It should be taken into consideration that the mechanisms of heat transfer processes differ due to the substantial differences in the mini-gap cross-section geometry.
Experimental investigation of water flow boiling heat transfer in a grooved-wall micro-channel was the main aim of Reference [7]. The influence of selected geometrical parameters of the test section was investigated, including the silicon base channel aspect ratios (1, 2.5 and 4), groove spacing ratios (2, 5 and 8) and groove depths (15, 30 and 45 μm) on the flow boiling. During the tests, flow patterns in the micro-channels were observed by using a high-speed camera. Experiments were conducted for mass flux of 446–963 kg/(m2∙s) and heat flux of 36.3–502.8 kW/m2. It was observed that the best improvement of heat transfer was achieved using the aspect ratio of 2.5 in the grooved-wall channel. The heat transfer coefficient value firstly increased as the aspect ratio increased from 1 to 2.5, and then decreased when the aspect ratio increased from 2.5 to 4. Research was limited to testing one modified developed surface (a grooved wall). Using other developed wall surfaces could significantly affect the results. As the tested micro-channel heat sink was fabricated by the standard silicon MEMS fabrication technique, the geometry variation could be over a wide range.
In Reference [8], the study of an experimental research on the water flow boiling heat transfer in a mini-channel with a 2.12 mm inner diameter, vertically oriented with the upward flow, was presented. The experimental parameters tested were the mass flux from 20 to 120 kg/(m2·s) and the average vapor quality in the test section inlet of 0.05 to 0.9. The two-phase flow patterns were observed at the test section and classified as three patterns: slug, annular and churn flows. The characteristics of the heat transfer and fluid flow were analyzed. It was observed that the highest mass flux yielded the highest convective heat transfer coefficient. Moreover, the pressure loss was higher at a higher mass flux. Only one channel geometry was applied in the experiments. Experimental research focused on using water as the working fluid. Using other coolers is needed to compare the results. Moreover, setting the downward flow or other spatial orientation of the test section will significantly influence the data obtained from experiments, especially forming two-phase structures.
The main aim of Reference [9] was boiling heat transfer during deionized water flow in a single mini-channel of 1.8 mm width and 0.9 mm height. The following experimental parameters were set during the tests: the wall heat flux in the range of 204.7–431.3 kW/m2 and mass fluxes (20, 40 and 60 kg/(m2·s)). It was noticed that the heat transfer coefficient showed an irregular trend with increasing heat flux or vapor quality, while it increased with increasing mass flux. Furthermore, the pressure drop increased with an increase in heat flux or in exit vapor quality. During the experiments, the bubbly, slug and annular flow patterns were recorded. Experiments were conducted only for low mass fluxes. The authors did not compare their own research results with refrigerants and other geometries of the mini-channel.
The study of flow boiling of water and n-pentane/water and freon-11/water emulsion of heat transfer in a mini-channel was presented in Reference [10]. In the test section, a mini-channel of 1.1 mm hydraulic diameter, horizontally oriented, was produced. The mean heat transfer coefficient from the channel wall to emulsion was higher than to water at bubble and slug flow conditions of a two-phase flow. The variations in the liquid temperature at the channel inlet and outlet and the excess pressure in the system during the instability of two-phase flow were examined. The experiments were conducted at mass fluxes of 42 and 75 kg/(m2·s), and with an increase in the mass flow, the segment where water boiling begins additionally shifts closer to the channel outlet. The presence of an additional vapor phase in a mini-channel, which consists of the vapor of a low-boiling liquid, contributes to the appearance of instability and boiling crisis at lower heat flux, as compared with pure water. The intensification of boiling heat transfer in a mini-channel was observed in experiments in which emulsions were used in comparison to pure water. The temperature of a two-phase flow at the outlet during the boiling of an emulsion proved to be higher than the boiling temperature of a dispersed phase. The study does not take into account the influence of the spatial orientation of the mini-channel on forming two-phase flow structures. The influence of the variable geometry of the channel on heat transfer could also be investigated.
Pool boiling heat transfer was the subject of numerous publications. Examples of articles presenting the results of such a topic can be found in References [11,12,13,14,15,16]. Research on distilled water and ethyl alcohol pool boiling on copper micro-channels modified by a laser beam was described in [11] and wire mesh coatings made of copper and phosphor bronze in [12]. It was detected that the laser texturing provided a significant improvement in pool boiling heat transfer and enabled the production of highly efficient phase-change heat exchangers [11]. Moreover, it was shown that the metal meshes sintered onto the heater surfaces can be very effective in providing enhancement of boiling heat transfer [12].
The authors of Reference [13] studied pool boiling heat transfer of water at sub-atmospheric pressure and low saturation temperature. The experiments were conducted for three vapor pressures (2.4, 3.1 and 4.1 kPa), four levels of liquid and five applied heat fluxes (in the range 3.6–7.1 kW/m2). The following boiling regimes were identified: convection or small popping bubbles, isolated bubbles, intermittent boiling and fully developed boiling. For the high level of liquid, mostly convection or small popping bubbles were shown.
Experimental results relating to pool boiling heat transfer on open micro-channel surfaces with dimensions from 0.2 to 0.4 mm wide and from 0.2 to 0.5 mm deep [14,16], and on plain micro-fins and micro-fins with a porous layer of 1.0 mm high [15], were discussed. The boiling heat transfer coefficients achieved from the surface open micro-channel surfaces were 3 times higher [14] or 2.5–4.9 times higher [16] compared to the smooth surface. The plain micro-fin surfaces helped to obtain an increase in the heat transfer coefficient of two times that of the micro-fins with perforated foil [15].
Using processes of heat transfer with the change of phase, which occur both during boiling and in the reverse process, i.e., condensation, are an effective way of cooling. In the literature, investigations of heat transfer during condensation in small channels are widely discussed [17,18].
At Kielce University of Technology, previous investigations on boiling heat transfer during cooling liquids’ flow [19,20,21] along rectangular mini-channels of changeable spatial orientation, with smooth or modified heater surfaces, were discussed in [20,22,23]. The thermoelasticity problem was the main topic of References [24,25] dealing with a two-layer plate [24] and non-contacting face seals [25].
This work provides a summary of the results collected during experimental work in a laboratory setup, in which an essential element was a model mini-channel heat exchanger with a singular long mini-channel or a set of short parallel mini-channels. It aims to deliver a summary of observations and comments concerning boiling heat transfer during refrigerant flow in rectangular mini-channels of different dimensions (depth, length and width), with changeable inclination to the horizontal plane and asymmetrically heated—with a smooth heated wall or a modified one. A wide range of experimental tests and a large number of cooling fluids used in the study are additional advantages of this work and a novelty in experimental investigation. Up to now, preliminary research has been performed on the 3 and 5 mini-channel test sections [26,27]. Previous studies were carried out using a singular mini-channel, which are summarized here, with many additional parameters of the experiment being tested.

2. The Experimental Rig and Methodology

2.1. Main Circuits

The main systems of the experimental setup are illustrated in Figure 1. It consists of: a FC-72 circulating loop including the test section with a singular mini-channel or five parallel mini-channels, the current supply and control system with the heat source and the data and image acquisition system comprising the lighting system. Infrared cameras (FLIR, model E60 or A655SC, USA) were used to determine changes in temperature on the outer surface of the mini-channel’s heated wall. Two-phase flow patterns were observed through the mini-channel’s glass wall using a high-speed camera (JAI, model SP-50 0 0M-CXP2, Denmark/Japan) and illuminated by high-power LEDs.
The mass flow rate was measured due to the Coriolis mass flow meter (Endress+Hauser, model Proline Promass A 100, Switzerland). The pressure of the fluid at the inlet and outlet of the mini-channels was measured by pressure sensors (Endress+Hauser, model PMP71, Cerabar S, Switzerland).
The test section was set at four inclination angles in relation to the horizontal plane (0°, 90°, 180° and 270°). Four refrigerants (FC-72, HFE-649, HFE-7000 and HFE-7100) were used as the working fluid flowing along mini-channels, and asymmetrically heated.

2.2. Test Sections

Two constructions of the test section were applied during experiments: the first one with a singular mini-channel of 1.7 mm depth and 180 mm length (named as “the 1st”) and the second with five parallel mini-channels of 1 mm depth and 32 mm length (named as “the 2nd”). Geometrical characteristics of the test sections are specified in Table 1, while their views and the longitudinal section are shown in Figure 2.
Figure 2a presents a view of the 1st test section components (the 2nd test section is similarly constructed). The mini-channel is positioned between the glass plate (7) and the heated plate (6), which is made of Haynes-230 alloy. The channel’s side walls are formed by Teflon plates (2). At the inlet and outlet of the test section, a K-type thermocouple and pressure sensors are mounted (9). A silicone O-ring (8) is placed in the groove to avoid the leakage of the fluid, whereas a graphite plate (3) is applied for electrical insulation between the heated plate and other metal elements of the test section. In Figure 2b, a longitudinal view of the test section is presented. In Figure 2c,d, the views of the test sections from the side of the heated plate are shown, of the 1st test section (Figure 2c) and of the 2nd one (Figure 2d).

2.3. Experimental Parameters and Accuracies

The main features of the experiments are specified in Table 2, including: basic experimental parameters, spatial positions of the test section, working fluids and heated plate specification. Boiling regions are also listed.
Absolute errors of the main experimental parameters, estimated on the basis of the data of manufacturers of measurement devices, are provided in Table 3. The details of the conducted error analysis can be found in [21].
The temperature of the fluid at the inlet and outlet of the mini-channel was measured due to K-type thermocouples manufactured by Czaki Thermo-Product (Poland). Thermocouples were connected via a thermocouple module (DBK84) to the data acquisition station (DaqLab/2005). Such temperature measurement system has been used in the laboratory for many years. All thermocouples were additionally calibrated by using a calibrator of thermocouple type 9102 HDRC, manufactured by Hart Scientific (USA). The calibrating procedure, the detailed information concerning the temperature measurement system with K-type thermocouples and the absolute error of temperature measurement were described in [28].
The basic physical properties of the refrigerants under experimental investigation, i.e., FC-72, HFE-649, HFE-7000 and HFE-7100 (3M), are shown in Table 4.
The heated plates were made of Haynes-230 alloy or Hastelloy X alloy, and were 0.45 or 0.65 mm thick. The arithmetic mean deviations of the roughness profile of the smooth plate surface are as follows: Ra = 0.175 μm (Haynes-230) and Ra = 0.133 μm (Hastelloy X). Modified surfaces of the heated plates were investigated in terms of intensification of heat transfer. Views of the investigated modified heated plate surfaces are revealed in Figure 3. Several types of modified surfaces were tested, as follows: vibration-assisted laser surface texturing (three textures dependent on the applied parameters of the production process, named as No. 1, No. 2 and No. 3), electroerosion surface texturing and sintering or soldering iron powder to the alloy plate base (porous surfaces).
During experiments, the test section was set at various inclinations to the horizontal plane. Four spatial orientations of the test section were tested: 0°, 90°, 180° and 270°. Figure 4 illustrates the positioning of the heated wall and the mini-channel in the test section. The flow direction in each case is also indicated.

3. HTC Determination

Heat transfer coefficient (HTC) was determined from the 1D approach of heat flow across the test section layers. In this simple mathematical approach, only one direction of heat flow was assumed, i.e., perpendicular to the fluid flow direction and the mini-channel width. The 2D approach using Trefftz functions for solving the inverse problem and the commercial software ADINA were also applied in the computations. The results were discussed in the authors’ previous works [19,21,26,27,29]. It allowed to verify the correctness of the 1D approach. According to the results obtained from the 1D and 2D approaches at the saturated boiling region presented in [30], it was concluded that the distribution of HTC values was similar for both methods at the same heat flux.
Local HTCs (HTC(x)) were calculated from the following dependence [20]:
H T C x = q w / T P x , δ P T f x q w · δ P λ P
where:
x—direction perpendicular to the fluid flow direction and the mini-channel width,
TP(x)—local temperature of the heated plate outer surface measured by infrared thermography,
Tf(x)—local fluid temperature based on the assumption of the linear distribution from the channel inlet to the outlet: (i) of the liquid bulk temperature (Tl) at the subcooled boiling region, and (ii) of the liquid saturated temperature (Tsat) at the saturated boiling region, both based on experimental measurements (of fluid temperature or pressure at the inlet and outlet, respectively),
λP—thermal conductivity of the heated plate,
δP—thickness of the heated plate,
qw—heat flux transferred from the heated plate to the fluid, defined as:
q w = I · Δ U / A
where:
I—current supplied to the heated plate,
U—voltage drop across the heated plate,
A—area of the heated plate.
The total heat flux transferred from the heated plate to the fluid was reduced by heat loss to the ambient air from the heated plate (calculated similarly as in [22]).
It is assumed that the average uncertainties of the HTC and qw achieve similar values, as determined in previous experiments conducted by the authors, discussed in works such as References [21,22]. Average uncertainties were calculated in accordance with the principles of measurement accuracy analysis. The uncertainties of HTC were in the range 11.0–19.8% at the subcooled boiling region [22]. At the saturated boiling region, the errors of the HTC were noticeably higher, from 14.2% to 47.9% [21]. The average uncertainties of qw for both boiling regions were in the range from 5.5% to 10.5% [22]. When analyzing the tested experiments, the highest value of the uncertainty of qw equals 4.7% at the lowest heat flux.

4. Results and Discussion

4.1. Database of Selected Experiments

The database of experiments due to the study of the influence of the parameter/factor on flow boiling was large. To clearly compile the data adopted for analysis from a wide experimental database source according to the basic criteria of their selection, Table 5 illustrates their summary. In the second column of this table, the criterion of the analysis is provided, referring to the number of figures illustrating the selected results from experiments. In other columns, more detailed information concerning the geometrical parameters of the test section, working fluid and thermal-flow parameters are shown with specifications of the boiling region.

4.2. The Results—Analysis and Observations

As in most scientific works, the results of heat transfer are presented in the form of temperature graphs, HTC, flow structures and boiling curves.
Selected data from the experiments are compiled as the relationship between:
Selected two-phase flow structures showing characteristic data are also illustrated graphically in Figure 12, Figure 13, Figure 14 and Figure 15.
Comprehensive characteristics of the influence of selected factors on boiling heat transfer during flow in mini-channels are shown in Table 6. A database of selected experiments showing the main information concerning the investigations was listed in Table 5. It is worth underlining that at the subcooled boiling region, HTC was relatively low in comparison to HTC obtained at the saturated boiling region.
Boiling curves are essential to recognize heat transfer mechanisms in flow boiling. The boiling curves shown in Figure 16 were generated on the basis of the experimental data during increasing heat flux supplied to the mini-channel’s heated plate.
Boiling curves show the relationship between heat flux, qw, and the difference in temperature, TP – Tf. They were generated at selected distances from the mini-channel inlet. The curves presented on one graph were constructed for: various mass fluxes (Figure 16a,c), two vertical orientations of the mini-channel (Figure 16b), two horizontal orientations of the mini-channel and vertical mini-channel with fluid up-flow (Figure 16d), four refrigerants (Figure 16e) and mini-channels with smooth and modified heated plate surfaces (Figure 16f). It is known that at the subcooled boiling region, an increase in heat flux may cause a sudden decrease in the heated-wall temperature during the incipience of boiling and detachment of the first bubbles. This phenomena is known as “nucleation hysteresis”. Points marked as ‘ONB’ (onset of nucleate boiling) on the boiling curves’ courses were selected due to the highest difference in temperature, TP – Tf. When analyzing the boiling curves presented in Figure 16i, it can be pointed out that:
  • With the increase of mass flux, the temperature difference corresponding to ONB increased—see Figure 16a,c.
  • The temperature difference at ONB was higher for O—270° in comparison to O—90°, see Figure 16b.
  • The highest temperature difference at ONB was observed for the vertical channel with fluid up-flow (O—90°), in comparison to both horizontal orientations (O—0° and O—180°)—see Figure 16d.
  • The highest temperature difference at ONB occurred for HFE 7100 in comparison to other tested fluids—see Figure 16e.
  • The lowest temperature difference at ONB was gained for electroerosion and vibration-assisted laser No. 2 surface texturing, and the highest one for porous surfaces produced by soldering iron powder—see Figure 16f.

5. Conclusions

This work presented the results of boiling heat transfer during flow in rectangular mini-channels. The influence of selected factors on heat transfer were discussed, i.e., main experimental parameters (inlet pressure and mass flux), mini-channel dimensions and spatial orientation, using modified heated plate surfaces. Four refrigerants (FC-72, HFE-649, HFE-7000, HFE-7100) were applied as the working fluid. The test sections with a single rectangular mini-channel (180 mm long) or with a group of five parallel mini-channels (each 32 mm long) were used during experiments performed at stationary condition.
The results were presented as local heated-wall temperatures and heat transfer coefficients and boiling curves. Two-phase flow patterns were also illustrated and discussed.
The major remarks and conclusions can be drawn concerning the influence of selected factors on heat transfer during flow boiling in mini-channels, as follows:
  • Inlet pressure (subcooled boiling region):
    Higher HPT was achieved at higher pin, and higher HTC at lower pin (180 mm long channel).
  • Mass flux (subcooled boiling region):
    Higher HPT was reached for lower G (both test sections), the highest HTC was obtained for the lowest G at the highest qw (180 mm long channel), and for the higher G at lower qw (32 mm long channel), and larger bubbles were observed near channel side walls for the lowest G.
  • Spatial orientation of a mini-channel (subcooled boiling region):
    The highest HPT was noticed for the lowest G and at O—270°, for lower qw. The highest HTC was achieved for higher G and all qw, for O—270°, and more gaseous phase with locally elongated agglomerates was observed for O—270° (180 mm long channel).
    The highest HPT was detected for all qw, the highest HTC was obtained for the highest qw and the lowest and the least fragmented bubbles (32 mm long channel) for O—0°.
  • Working fluids:
    The highest HTC was achieved for HFE-7000 for all qw, at the subcooled boiling region, and for FC-72 and HFE-7000, for all qw, at the saturated boiling region. More developed vapor structures as elongated agglomerates were observed for FC-72 in comparison to other refrigerants, and for HFE-7000, bubbly, bubbly-slug and slug structures were recorded.
  • Heated plate surface of a mini-channel:
    The highest HTC was obtained for electroerosion and vibration-assisted laser No. 2 surface texturing at both boiling regions. The highest vapor share in two-phase flow structures with elongated agglomerates in comparison to other tested plate surfaces was observed for electroerosion surface texturing, and bubbles of larger sizes and elongated agglomerates were recorded for vibration-assisted laser No. 2 surface texturing.
When analyzing boiling curves, it was observed that with the increase of mass flux, the temperature difference corresponding to ONB increased, the temperature difference at ONB was higher for O—270° in comparison to O—90°, the highest temperature difference at ONB was observed for vertical channel orientation with fluid up-flow (O—90°), the highest temperature difference at ONB was detected for HFE 7100 in comparison to other tested fluids and the lowest temperature differences at ONB were gained for electroerosion and vibration-assisted laser No. 2 surface texturing, and the highest for porous surfaces produced by soldering iron powder.

Author Contributions

Conceptualization, methodology, M.P. and K.S.; investigations, calculations, visualization, K.S.; validation, formal analysis, writing—original draft preparation, review and editing, M.P. and K.S.; funding acquisition, M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Centre, Poland, grant number UMO-2018/31/B/ST8/01199.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

The research reported herein was supported by a grant from the National Science Centre, Poland, No. UMO-2018/31/B/ST8/01199.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Aarea of the heated plate, m2
ffluid
Gmass flux, kg/(m2·s)
HPTheated plate temperature, K
HTCheat transfer coefficient, W/(m2·K)
Icurrent, A
Llength of the mini-channel, m
Ochannel orientation
ONBonset of nucleate boiling
pinpressure at the mini-channel inlet, Pa
qwheat flux density, W/m2
Raarithmetic mean deviation of the roughness profile of the plate surface, μm
Sheated plate surface
Ttemperature, K
the 1stthe test section with a singular mini-channel of LIRT = 180 mm
the 2ndthe test section with five parallel mini-channels of LIRT = 32 mm
xdirection perpendicular to the fluid flow direction and the channel width, m
Greek symbols
δdepth, thickness, m
Uvoltage drop, V
λthermal conductivity, W/(m·K)
Subscripts
ffluid
IRTinfrared thermography
lliquid
Pplate
satsaturated

References

  1. Kanizawa, F.T.; Tibiriçá, C.B.; Ribatski, G. Heat transfer during convective boiling inside microchannels. Int. J. Heat Mass Transf. 2016, 93, 566–583. [Google Scholar] [CrossRef]
  2. Sun, Z.-C.; Ma, X.; Ma, L.-X.; Li, W.; Kukulka, D. Flow Boiling Heat Transfer Characteristics in Horizontal, Three-Dimensional Enhanced Tubes. Energies 2019, 12, 927. [Google Scholar] [CrossRef] [Green Version]
  3. Jige, D.; Inoue, N. Boiling heat transfer, pressure drop, and flow pattern in a horizontal square minichannel. Int. J. Heat Fluid Flow 2019, 78, 108433. [Google Scholar] [CrossRef]
  4. See, Y.S.; Leong, K.C. Entropy generation for flow boiling on a single semi-circular minichannel. Int. J. Heat Mass Transf. 2020, 154, 119689. [Google Scholar] [CrossRef]
  5. He, B.; Luo, X.; Yu, F.; Zhou, J.; Zhang, J. Flow boiling characteristics in bi-porous minichannel heat sink sintered with copper woven tape. Int. J. Heat Mass Transf. 2020, 158, 119988. [Google Scholar] [CrossRef]
  6. Klugmann, M.; Dąbrowski, P.; Mikielewicz, D. Flow distribution and heat transfer in minigap and minichannel heat exchangers during flow boiling. Appl. Therm. Eng. 2020, 181, 116034. [Google Scholar] [CrossRef]
  7. Cheng, X.; Yao, Y.; Wu, H. An experimental investigation of flow boiling characteristics in silicon-based groove-wall microchannels with different structural parameters. Int. J. Heat Mass Transf. 2021, 168, 120843. [Google Scholar] [CrossRef]
  8. Enoki, K.; Ono, M.; Okawa, T.; Kristiawan, B.; Wijayanta, A.T. Water flow boiling heat transfer in vertical minichannel. Exp. Therm. Fluid Sci. 2020, 117, 110147. [Google Scholar] [CrossRef]
  9. Markal, B.; Candan, A.; Aydin, O.; Avci, M. Experimental investigation of flow boiling in single minichannels with low mass velocities. Int. Commun. Heat Mass Transf. 2018, 98, 22–30. [Google Scholar] [CrossRef]
  10. Gasanov, B.M. Flow boiling of water and emulsions with a low-boiling disperse phase in minichannels. Int. J. Heat Mass Transf. 2018, 126, 9–14. [Google Scholar] [CrossRef]
  11. Orman, Ł.J.; Radek, N.; Pietraszek, J.; Szczepaniak, M. Analysis of enhanced pool boiling heat transfer on laser-textured surfaces. Energies 2020, 13, 2700. [Google Scholar] [CrossRef]
  12. Dąbek, L.; Kapjor, A.; Orman, Ł.J. Distilled water and ethyl alcohol boiling heat transfer on selected meshed surfaces. Mech. Ind. 2019, 20, 701. [Google Scholar] [CrossRef]
  13. Wojtasik, K.; Rullière, R.; Krolicki, Z.; Zajaczkowski, B.; Bonjour, J. Subcooled boiling regime map for water at low saturation temperature and subatmospheric pressure. Exp. Therm. Fluid Sci. 2020, 118, 110150. [Google Scholar] [CrossRef]
  14. Hożejowska, S.; Kaniowski, R.; Pastuszko, R. Application of the Trefftz Method for Pool Boiling Heat Transfer on Open Microchannel Surfaces. Heat Transf. Eng. 2021, 1–13. [Google Scholar] [CrossRef]
  15. Pastuszko, R.; Kaniowski, R.; Wójcik, T.M. Comparison of pool boiling performance for plain micro-fins and micro-fins with a porous layer. Appl. Therm. Eng. 2020, 166, 114658. [Google Scholar] [CrossRef]
  16. Kaniowski, R.; Pastuszko, R. Pool Boiling of Water on Surfaces with Open Microchannels. Energies 2021, 14, 3062. [Google Scholar] [CrossRef]
  17. Bohdal, T.; Kruzel, M. Refrigerant condensation in vertical pipe minichannels under various heat flux density level. Int. J. Heat Mass Transf. 2020, 146, 118849. [Google Scholar] [CrossRef]
  18. Kruzel, M.; Bohdal, T.; Sikora, M. Heat transfer and pressure drop during refrigerants condensation in compact heat exchangers. Int. J. Heat Mass Transf. 2020, 161, 120283. [Google Scholar] [CrossRef]
  19. Strąk, K.; Piasecka, M.; Maciejewska, B. Spatial orientation as a factor in flow boiling heat transfer of cooling liquids in enhanced surface minichannels. Int. J. Heat Mass Transf. 2018, 117, 375–387. [Google Scholar] [CrossRef]
  20. Strąk, K.; Piasecka, M. The applicability of heat transfer correlations to flows in minichannels and new correlation for subcooled flow boiling. Int. J. Heat Mass Transf. 2020, 158, 119933. [Google Scholar] [CrossRef]
  21. Piasecka, M.; Strąk, K.; Maciejewska, B. Heat transfer characteristics during flow along horizontal and vertical minichannels. Int. J. Multiph. Flow 2021, 137, 103559. [Google Scholar] [CrossRef]
  22. Piasecka, M.; Strąk, K. Influence of the Surface Enhancement on the Flow Boiling Heat Transfer in a Minichannel. Heat Transf. Eng. 2019, 40, 1162–1175. [Google Scholar] [CrossRef]
  23. Depczyński, W.; Piasecki, A.; Piasecka, M.; Strąk, K. Impact of Fe powder sintering and soldering in production of porous heating surface on flow boiling heat transfer in minichannels. E3S Web Conf. 2017, 19, 03012. [Google Scholar] [CrossRef] [Green Version]
  24. Maciąg, A.; Grysa, K. Trefftz Method of Solving a 1D Coupled Thermoelasticity Problem for One- and Two-Layered Media. Energies 2021, 14, 3637. [Google Scholar] [CrossRef]
  25. Blasiak, S. Influence of Thermoelastic Phenomena on the Energy Conservation in Non-Contacting Face Seals. Energies 2020, 13, 5283. [Google Scholar] [CrossRef]
  26. Piasecka, M.; Maciejewska, B.; Łabędzki, P. Heat Transfer Coefficient Determination during FC-72 Flow in a Minichannel Heat Sink Using the Trefftz Functions and ADINA Software. Energies 2020, 13, 6647. [Google Scholar] [CrossRef]
  27. Piasecka, M.; Hożejowska, S.; Maciejewska, B.; Pawińska, A. Time-Dependent Heat Transfer Calculations with Trefftz and Picard Methods for Flow Boiling in a Mini-Channel Heat Sink. Energies 2021, 14, 1832. [Google Scholar] [CrossRef]
  28. Michalski, D.; Strąk, K.; Piasecka, M. Estimating uncertainty of temperature measurements for studies of flow boiling heat transfer in minichannels. EPJ Web Conf. 2019, 213, 02059. [Google Scholar] [CrossRef]
  29. Maciejewska, B.; Łabędzki, P.; Piasecki, A.; Piasecka, M. Comparison of FEM calculated heat transfer coefficient in a minichannel using two approaches: Trefftz base functions and ADINA software. EPJ Web Conf. 2017, 143, 02070. [Google Scholar] [CrossRef] [Green Version]
  30. Piasecka, M.; Strąk, K.; Maciejewska, B. Calculations of Flow Boiling Heat Transfer in a Minichannel Based on Liquid Crystal and Infrared Thermography Data. Heat Transf. Eng. 2017, 38, 332–346. [Google Scholar] [CrossRef]
  31. Strąk, K.; Piasecka, M. Research on flow boiling heat transfer for vertical upward and downward flows along a minichannel with a smooth heated surface. E3S Web Conf. 2018, 70, 02014. [Google Scholar] [CrossRef]
  32. Strąk, K.; Piasecka, M.; Strąk, D. Enhanced Surfaces Used in Research on Flow Boiling Heat Transfer in Minichannels. Mater. Res. Proc. 2018, 5, 154–159. [Google Scholar]
Figure 1. The main systems of the experimental setup: 1—test section, 2—gear pump, 3—compensating tank/pressure regulator, 4—tube-type heat exchanger, 5—filter, 6—mass flow meter, 7—deaerator, 8—pressure sensor, 9—high-speed camera, 10—LEDs (400 W), 11—data acquisition station, 12—PC, 13—infrared camera, 14—power supply system.
Figure 1. The main systems of the experimental setup: 1—test section, 2—gear pump, 3—compensating tank/pressure regulator, 4—tube-type heat exchanger, 5—filter, 6—mass flow meter, 7—deaerator, 8—pressure sensor, 9—high-speed camera, 10—LEDs (400 W), 11—data acquisition station, 12—PC, 13—infrared camera, 14—power supply system.
Energies 14 04891 g001
Figure 2. General characteristics of test sections: (a) view of the 1st test section components, (b) longitudinal section; (c) view of the 1st test section; (d) view of the 2nd test section; 1—channel body, 2, 4—Teflon plates, 3—graphite plate, 5—front cover, 6—heated plate, 7—glass panel, 8—O-ring, 9—temperature and pressure sensors.
Figure 2. General characteristics of test sections: (a) view of the 1st test section components, (b) longitudinal section; (c) view of the 1st test section; (d) view of the 2nd test section; 1—channel body, 2, 4—Teflon plates, 3—graphite plate, 5—front cover, 6—heated plate, 7—glass panel, 8—O-ring, 9—temperature and pressure sensors.
Energies 14 04891 g002
Figure 3. View of modified heated plate surfaces produced by: (ac) vibration-assisted laser texturing, No. 1 (a), No. 2 (b), No. 3 (c); (d) electroerosion texturing, and (e,f) porous metal surfaces produced by soldering (e) and (f) sintering iron powder to the alloy base. Ra—arithmetic mean deviation of the roughness profile of the plate surface.
Figure 3. View of modified heated plate surfaces produced by: (ac) vibration-assisted laser texturing, No. 1 (a), No. 2 (b), No. 3 (c); (d) electroerosion texturing, and (e,f) porous metal surfaces produced by soldering (e) and (f) sintering iron powder to the alloy base. Ra—arithmetic mean deviation of the roughness profile of the plate surface.
Energies 14 04891 g003
Figure 4. Spatial orientations of the test section.
Figure 4. Spatial orientations of the test section.
Energies 14 04891 g004
Figure 5. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; various inlet pressures; vertical channel orientation (90°); mass flux of 400 kg/(m2∙s); modified heated plate surface (vibration-assisted laser surface texturing No. 2); HFE-649; LIRT = 180 mm.
Figure 5. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; various inlet pressures; vertical channel orientation (90°); mass flux of 400 kg/(m2∙s); modified heated plate surface (vibration-assisted laser surface texturing No. 2); HFE-649; LIRT = 180 mm.
Energies 14 04891 g005
Figure 6. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; various mass fluxes; vertical channel orientation (90°); smooth heated plate; FC-72; LIRT = 180 mm.
Figure 6. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; various mass fluxes; vertical channel orientation (90°); smooth heated plate; FC-72; LIRT = 180 mm.
Energies 14 04891 g006
Figure 7. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; various mass fluxes; two vertical channel orientations (90° and 270°); smooth heated plate; HFE-649; LIRT = 180 mm.
Figure 7. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; various mass fluxes; two vertical channel orientations (90° and 270°); smooth heated plate; HFE-649; LIRT = 180 mm.
Energies 14 04891 g007
Figure 8. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; various mass fluxes; vertical channel orientation (90°); smooth heated plate; FC-72; LIRT = 32 mm.
Figure 8. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; various mass fluxes; vertical channel orientation (90°); smooth heated plate; FC-72; LIRT = 32 mm.
Energies 14 04891 g008
Figure 9. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; three spatial channel orientations (0°, 90° and 180°); mass flux of 440 kg/(m2∙s), smooth heated plate; FC-72; LIRT = 32 mm.
Figure 9. Plate temperature measurements using the IR camera (a,c,e) and heat transfer coefficient (b,d,f), both vs. distance from the mini-channel inlet; subcooled boiling region; three spatial channel orientations (0°, 90° and 180°); mass flux of 440 kg/(m2∙s), smooth heated plate; FC-72; LIRT = 32 mm.
Energies 14 04891 g009
Figure 10. Heat transfer coefficient vs. distance from the mini-channel inlet: subcooled boiling region (a,b) and saturated boiling region (c,d); various refrigerants; vertical channel orientation (90°); average mass flux of 435 kg/(m2∙s), modified heated plate surface (vibration-assisted laser surface texturing No. 2); FC-72; LIRT = 180 mm.
Figure 10. Heat transfer coefficient vs. distance from the mini-channel inlet: subcooled boiling region (a,b) and saturated boiling region (c,d); various refrigerants; vertical channel orientation (90°); average mass flux of 435 kg/(m2∙s), modified heated plate surface (vibration-assisted laser surface texturing No. 2); FC-72; LIRT = 180 mm.
Energies 14 04891 g010
Figure 11. Heat transfer coefficient vs. distance from the mini-channel inlet: subcooled boiling region (a,b) and saturated boiling region (c,d); smooth and several modified heated plate surfaces; vertical channel orientation (90°); average mass flux of 353 kg/(m2∙s); FC-72; LIRT = 180 mm.
Figure 11. Heat transfer coefficient vs. distance from the mini-channel inlet: subcooled boiling region (a,b) and saturated boiling region (c,d); smooth and several modified heated plate surfaces; vertical channel orientation (90°); average mass flux of 353 kg/(m2∙s); FC-72; LIRT = 180 mm.
Energies 14 04891 g011
Figure 12. Images of two-phase flow structures recorded for: various mass fluxes and two vertical channel orientations (90° and 270°), HFE-649, smooth heated plate, LIRT = 180 mm.
Figure 12. Images of two-phase flow structures recorded for: various mass fluxes and two vertical channel orientations (90° and 270°), HFE-649, smooth heated plate, LIRT = 180 mm.
Energies 14 04891 g012
Figure 13. Images of two-phase flow structures recorded for: various mass fluxes and three channel orientations (0°, 90° and 180°), FC-72, smooth heated plate, LIRT = 32 mm.
Figure 13. Images of two-phase flow structures recorded for: various mass fluxes and three channel orientations (0°, 90° and 180°), FC-72, smooth heated plate, LIRT = 32 mm.
Energies 14 04891 g013
Figure 14. Images of two-phase flow structures recorded for: various refrigerants, constant mass flux, vertical channel orientation (90°), modified heated plate surface (vibration-assisted laser surface texturing No. 2), LIRT = 180 mm.
Figure 14. Images of two-phase flow structures recorded for: various refrigerants, constant mass flux, vertical channel orientation (90°), modified heated plate surface (vibration-assisted laser surface texturing No. 2), LIRT = 180 mm.
Energies 14 04891 g014
Figure 15. Images of two-phase flow structures recorded for: various modified heated plates, constant mass flux, vertical channel orientation (90°), FC-72, LIRT = 180 mm.
Figure 15. Images of two-phase flow structures recorded for: various modified heated plates, constant mass flux, vertical channel orientation (90°), FC-72, LIRT = 180 mm.
Energies 14 04891 g015
Figure 16. Boiling curves generated for two distances from the mini-channel inlet: (a,c,d) 0.030 m (LIRT = 180 mm), (b) 0.011 m (LIRT = 32 mm). Experimental parameters are shown in Table 5, Nos. #2 to #9. ONB - onset of nucleate boiling.
Figure 16. Boiling curves generated for two distances from the mini-channel inlet: (a,c,d) 0.030 m (LIRT = 180 mm), (b) 0.011 m (LIRT = 32 mm). Experimental parameters are shown in Table 5, Nos. #2 to #9. ONB - onset of nucleate boiling.
Energies 14 04891 g016
Table 1. Geometrical characteristics of mini-channels in the two test sections.
Table 1. Geometrical characteristics of mini-channels in the two test sections.
Geometrical Characteristics of Mini-ChannelsTest Section
the 1stthe 2nd
Number of mini-channels15
Length of a mini-channel, mm18032
Width of a mini-channel, mm166
Depth of a mini-channel, mm1.71
Table 2. Main features of the experiments.
Table 2. Main features of the experiments.
Features of ExperimentsSpecific Data
Mass flux, G89 ÷ 616 (kg/(m2∙s))
Inlet pressure, pin110 ÷ 370 (kPa)
Inlet liquid subcooling35 ÷ 58 (K)
Heat flux (density), qw2.5 ÷ 102 (kW/m2)
Spatial orientation of the test section0°, 90°, 180°, 270°
RefrigerantFC-72, HFE-649, HFE-7000, HFE-7100
Modified heated wall surface produced by:
  • vibration-assisted laser surface texturing (Nos. 1, 2 and 3)
  • electroerosion surface texturing
  • sintering iron powder to the plate base
  • soldering iron powder to the plate base
Boiling regionsubcooled/saturated
Table 3. Absolute errors of the main experimental parameters.
Table 3. Absolute errors of the main experimental parameters.
Experimental ParameterAbsolute ErrorMore Specific Data
Mass flow rate1.67 ×∙10−5 (kg/s)maximum error of reading
Absolute Pressure1 285 (Pa)absolute error of the measurement
Temperature of the heated plate, Tp1 (K)absolute error of the measurement (IRT)
Temperature of the fluid, Tf, at the inlet/outlet0.34 (K)absolute error of the measurement [21]
Table 4. Basic physical properties of refrigerants under experimental investigation.
Table 4. Basic physical properties of refrigerants under experimental investigation.
RefrigerantBoiling Point Temperature (K)Density (kg/m3)Kinematic Viscosity (m2/s)Vapor Pressure (Pa)Heat of Evaporation (J/kg)
FC-7232916800.38 ×∙10−630∙10388 ×∙103
HFE-64932216000.40 ×∙10−640∙10388 ×∙103
HFE-700030714000.32 ×∙10−665∙103142 ×∙103
HFE-710033415100.38 ×∙10−627∙103112 ×∙103
Table 5. The database of research on the influence of selected factors on flow boiling in mini-channels.
Table 5. The database of research on the influence of selected factors on flow boiling in mini-channels.
No.Tested
Factor
(Figure)
Boiling RegionChannel Length, LIRT,
Heated Plate Surface, S
Fluid, fSpatial
Orientation of
Mini-Channels, O
Heat Flux, qw
(kW/m2)
Mass Flux, G
(kg/(m2s))
Inlet
Pressure, pin
(kPa)
Inlet
Liquid
Subcooling (K)
#1Inlet pressure
(Figure 5)
subcooled180 mm, v.a. las. 2HFE-64990°50
57
67
400160, 257,
165, 322,
215, 370
53
#2Mass flux
(Figure 6)
subcooled180 mm, smoothFC-7290°48, 62, 78209, 414, 61616058
#3Mass flux,
channel orientation (Figure 7)
subcooled180 mm, smoothHFE-64990°, 270°42, 56, 72209, 414, 61617048
#4Mass flux
(Figure 8)
subcooled32 mm, smoothFC-7290°2.5, 9.8, 21.389, 36811138
#5Channel orientation
(Figure 9)
subcooled32 mm, smoothFC-72
90°
180°
2.5, 9.8, 21.344011035
#6Working fluid
(Figure 10a,b)
subcooled180 mm, v.a. las. 2FC-72
HFE-649
HFE-7000
HFE-7100
90°55, 8142317039
#7Working fluid
(Figure 10c,d)
saturated180 mm, v.a. las. 2FC-72
HFE-7000
HFE-7100
90°91, 10242820041
#8Heated plate surface
(Figure 11a,b)
subcooled180 mm, v.a. las. 2FC-7290°47, 7335612142
#9Heated plate surface
(Figure 11c,d)
saturated180 mm, smooth, v.a. las. 1,2,3, electr., sold. and sin. i.p.FC-7290°84, 9334914043
v.a. las.—vibration-assisted laser surface texturing, electr.—electroerosion texturing (spark erosion), sold/sin. i.p.—porous metal surfaces produced by soldering and sintering iron powder.
Table 6. Characteristics of the influence of selected factors on flow boiling in mini-channels.
Table 6. Characteristics of the influence of selected factors on flow boiling in mini-channels.
Tested
Factor
Main Experimental DataHPTHTCFlow Structures
Inlet
pressure
(pin)
LIRT = 180 mm
S—v.a. las. 2
f—HFE-649
O—90°
Figure 5a,c,e
  • Higher HPT at higher pin, for all qw
Figure 5b,d,f (subcooled boiling region)
  • The highest HTC at lower pin for all qw
  • Slightly higher HTC at lower pin for higher qw (Figure 5d,f)
  • HTC increase with the distance from the channel inlet at all pin for all qw
  • Similar HTC distribution
Remarks based on experiments
  • More vapor bubbles in the two-phase pattern observed at lower pin
  • Tiny bubbles observed only at higher pin for higher qw
Mass
flux (G)
LIRT = 180 mm
S—smooth
f—FC-72
f—HFE-649
O—90°
Figure 6a,c,e (FC-72)
  • Higher HPT for the lowest G, for all qw
  • The lowest HPT for the highest G at the channel inlet, for lower qw (Figure 6a,c)
Figure 6b,d,f (subcooled boiling region, FC-72)
  • The highest HTC for the lowest G and higher qw at the channel outlet (Figure 6d,f)
  • Slight impact G on HTC for the lowest qw, (Figure 6b)
  • HTC increase with the distance from the channel inlet for higher qw, (Figure 6b,d)
Figure 12a,b (HFE-649)
  • Only tiny vapor bubbles observed near channel side walls at the outlet
  • Larger bubbles observed near channel side walls for the lowest G (Figure 12a)
Mass
flux (G)
orientation (O)
LIRT = 180 mm
S—smooth
f—HFE-649
O—90°, 270°
Figure 7a,c,e
  • The highest HPT for the lowest G and lower qw, for O—270° (Figure 7a,c)
  • Higher HPT for the lowest G and all qw, for O—90°
  • The lowest HTP for higher G and all qw, for O—270°
Figure 7b,d,f (subcooled boiling region)
  • The highest HTC for higher G and all qw, for O—270°
  • The lowest HTC for lower G and lower qw, for O—270°
  • Slight impact qw on HTC, for O—90°
Figure 12a–d
  • Only tiny, single vapor bubbles observed near channel side walls at the outlet, for O—90°
  • More gaseous phase, locally elongated agglomerates, for O—270°
  • Developed annular structures formed near the channel inlet, for O—270° [31]
  • Elongated vapor agglomerates causing sudden changes in HTP and HTC, for O—270°
Mass flux
(G)
LIRT = 32 mm
S—smooth
f—FC-72
O—90°
Figure 8a,c,e
  • The highest HPT for lower G and lower qw (Figure 8a,c)
  • The above remark confirms the results for a longer mini-channel LIR = 180 mm (Figure 6a,c,e)
Figure 8b,d,f (subcooled boiling region)
  • The highest HTC for higher G and lower qw (Figure 8b,d)
  • The highest HTC for lower G and the highest qw ( Figure 8f)
  • HTC decrease with the distance from the channel inlet for lower qw (Figure 8b,d)
Figure 13b,d
  • Increasing amount of vapor bubbles and their dynamic formation with qw increase, for lower G
  • More spherical bubbles of larger size for the lowest G (Figure 13b)
  • Only tiny single bubbles observed mainly at the channel outlet, for higher G and higher qw (Figure 13d)
Orientation (O) LIRT = 32 mm
S—smooth
f—FC-72
O—0°, 90°, 180°
Figure 9a,c,e
  • The highest HPT for O—0° for all qw
  • The lowest HTP for O—90° and the lowest qw (Figure 9a)
  • Similar HPT for the highest qw for O—90° and O—180°
Figure 9b,d,f (subcooled boiling region)
  • HTC increase with the distance from the channel inlet (Figure 9f)
  • The highest HTC for the highest qw, for O - 0º (Figure 9f)
  • Relatively minor differences in HTC with the channel distance for lower qw, for all orientations (Figure 9b,d)
Figure 13a,b,c
  • The greatest size of bubbles, for O—180° (Figure 13c)
  • The lowest and the least fragmented bubbles, for O—0° (Figure 13a)
Fluid (f)LIRT = 180 mm
S—v.a. las. 2
f—FC-72
HFE-649
HFE-7000
HFE-7100
O—90°
Figure 10a,b (subcooled boiling region)
  • HTC increase with the distance from the channel inlet for all fluids, lower HTC in comparison to those obtained at saturated boiling region [19,22]
  • The highest HTC for HFE-7000, for all qw
  • The lowest HTC for HFE-7100, for lower qw despite the channel inlet (Figure 10a)
Figure 14a–d
  • Higher share of vapor in two-phase structures and more developed vapor structures as elongated agglomerates, observed for FC-72 in comparison to other refrigerants (Figure 14a)
  • Bubbly, bubbly-slug or slug structures observed for HFE-7000 [19] (Figure 14c)
Fluid (f)LIRT = 180 mm
S—v.a. las. 2
f—FC-72
HFE-7000,
HFE-7100
O—90°
Figure 10c,d (saturated boiling region)
  • The highest HTC at saturated boiling region, confirming previous results [19,22]
  • The highest HTC for FC-72 and HFE-7000, for all qw, which confirms the results shown in [19]
  • The lowest HTC for HFE-7100 for higher qw (Figure 10d)
Figure 14a–d
  • Only single bubbles visible for HFE-649 and HFE-7100, mainly near the side walls, at the channel outlet (Figure 14b,d)

The least numerous bubbles observed for the HFE-649 (Figure 14b)
Heated
plate
surface (S)
LIRT = 180 mm
S—smooth
v.a. las. 1,2,3, electr., sold.
and sin. i.p.
f—FC-72
O—90°
Figure 11a,b (subcooled boiling region)
  • The highest HTC for electroerosion and vibration-assisted laser No. 2 surface texturing for higher qw, confirming the results discussed in [32]
  • The lowest HTC for vibration-assisted laser No. 3 surface texturing for lower qw (Figure 11a) and porous surfaces produced by sintering iron powder for higher qw (Figure 11b)
  • Higher HTC for higher qw, for the surface produced by soldering iron powder [23]
Figure 15a–d
  • Spherical and tiny bubbles observed for a smooth heated plate (Figure 15a)
  • Small spherical bubbles observed for vibration-assisted laser No. 1 surface texturing(Figure 15b)
  • Bubbles of larger size and elongated agglomerates observed for the vibration-assisted laser No. 2 surface texturing (Figure 15c)
  • The highest vapor share in two-phase flow structures with elongated agglomerates in comparison to other tested plate surfaces, observed for electroerosion surface texturing(Figure 15d)
Heated
plate
surface (S)
LIRT = 180 mm
S—smooth
v.a. las. 1,2,3, electr., sold. and sin. i.p.
f—FC-72
O—90°
Figure 11c,d (saturated boiling region)
  • The highest HTC for electroerosion and vibration-assisted laser No. 2 surface texturing, similarly as in [32]
  • Lower HTC for porous metal surfaces produced by soldering and sintering iron powder in comparison to other tested surfaces, for all qw According to [32]
  • Very low HTC for porous metal surfaces with capillary-fibrous structures
LIRT—long channel measured by IR camera, S—heated surface, f—fluid, O—orientation, v.a. las.—vibration-assisted laser surface texturing, electr.—electroerosion texturing (spark erosion), sold/sin. i.p.—porous metal surfaces produced by soldering and sintering iron powder.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Piasecka, M.; Strąk, K. Characteristics of Refrigerant Boiling Heat Transfer in Rectangular Mini-Channels during Various Flow Orientations. Energies 2021, 14, 4891. https://doi.org/10.3390/en14164891

AMA Style

Piasecka M, Strąk K. Characteristics of Refrigerant Boiling Heat Transfer in Rectangular Mini-Channels during Various Flow Orientations. Energies. 2021; 14(16):4891. https://doi.org/10.3390/en14164891

Chicago/Turabian Style

Piasecka, Magdalena, and Kinga Strąk. 2021. "Characteristics of Refrigerant Boiling Heat Transfer in Rectangular Mini-Channels during Various Flow Orientations" Energies 14, no. 16: 4891. https://doi.org/10.3390/en14164891

APA Style

Piasecka, M., & Strąk, K. (2021). Characteristics of Refrigerant Boiling Heat Transfer in Rectangular Mini-Channels during Various Flow Orientations. Energies, 14(16), 4891. https://doi.org/10.3390/en14164891

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop