1. Introduction
Flow patterns in gas–liquid two-phase flows significantly affect heat transfer coefficients and pressure drops, as the mechanisms governing these phenomena are intrinsically associated with the spatial distribution of liquid and vapour phases [
1]. Transitions in flow patterns, which represent a form of instability in two-phase flow, have been extensively studied across various applications, including nuclear reactors, pipelines, boiling and condensation heat exchangers, and microchannels, among others. The determination of flow patterns typically relies on visual assessments, photographic imaging methods, or sudden variations in the average pressure drop of the system. This may also involve the use of conductivity measurements or high-frequency impedance probes for localised holdup sampling, as well as monitoring local pressure fluctuations and average holdup values. Recently, advancements in neural network methodologies [
2] and machine learning algorithms, which are trained on high-speed, high-resolution void fraction data [
3], have been introduced to facilitate the automatic recognition of flow regimes in two-phase systems. These innovative approaches hold the potential to enhance the precision of flow pattern identification and classification compared with traditional methods, offering new data-driven strategies for defining and analysing flow regimes in a more systematic and objective manner.
The existing body of literature on two-phase flow patterns is extensive; however, the flow dynamics and the resulting flow patterns within pulsating heat pipes (PHPs) remain inadequately explored. Typically, four distinct flow patterns can be identified within pulsating heat pipes: (i) bubbly or dispersed flow; (ii) slug flow, also referred to as slug/plug flow; (iii) semi-annular flow; and (iv) annular flow. It is common for two or more of these flow patterns to coexist in various sections of the same channel. In addition, the flow is characterised by inherent unsteadiness, which arises from alternating regions where liquid slugs occupy the entire cross-section of the pipe and areas where the flow comprises a liquid layer alongside a gas layer. The complexity of the flow structure is further increased by the simultaneous occurrence of regions exhibiting co-current flow, where both phases move in the same direction, and regions where the two phases are in counterflow, as observed when vapour bubbles advance toward the condenser while the surrounding liquid film is propelled back toward the evaporator [
4,
5].
An additional complication arises from the heterogeneity of the literature data due to the variety of heat transfer fluids utilised in the functioning of pulsating heat pipes. Besides water, two of the most common fluids used in PHPs are ethanol and perfluorohexane (FC-72), which exhibit significant differences in their properties, making them suitable for distinct heat transfer applications. Ethanol, a polar solvent with a relatively low boiling point of 78.4 °C, has high thermal conductivity and specific heat capacity, which allows it to efficiently absorb and transfer heat in various systems, particularly in cooling applications and heat exchangers. Its ability to mix with water also enhances its versatility in thermal management. In contrast, perfluorohexane, a non-polar, fluorinated hydrocarbon with a higher boiling point of 56.5 °C, offers exceptional thermal stability and low surface tension, making it ideal for applications requiring inertness and chemical resistance, such as in electronics cooling and specialised heat transfer fluids [
6]. The non-polar nature of perfluorohexane also results in lower heat capacity compared with ethanol, which can limit its effectiveness in certain heat transfer scenarios. Overall, the choice between ethanol and perfluorohexane in heat transfer applications depends on the specific thermal properties required and the environmental conditions of the application.
The inherent complexity of the flow structure, coupled with the challenges associated with flow visualisation and intrusive measurements in pulsating heat pipes, has hindered a comprehensive characterisation of flow regimes. This limitation has consequently affected the development of reliable flow pattern maps that could serve as predictive instruments in the design of pulsating heat pipes (PHPs) and assist in conducting accurate numerical simulations. While several flow pattern maps have been established for mini- and micro-channels, there exists a scarcity of flow pattern maps specific for pulsating heat pipes. Initial efforts to systematically characterise flow patterns in pulsating heat pipes [
7,
8] have simply reported empirical findings in tabular format, correlating these observations with variables such as the heat input and the filling ratio. Typically, the heat input is utilised in the formulation of flow pattern maps, in conjunction with either the inclination angle or the number of turns [
9,
10].
Thongdaeng et al. [
11] proposed a flow pattern map for a top heat, closed-loop pulsating heat pipe using as coordinates the momentum fluxes of the liquid and the vapour, respectively, calculated based on flow rate and temperature measurements. They identified six flow patterns, resulting from different combinations of elementary flow patterns: (I) slug, dispersed bubble, and annular flow; (II) slug, dispersed bubble, and churn flow; (III) slug, bubble, and annular flow; (IV) slug and dispersed bubble flow; (V) bubble and annular flows; and (VI) slug and annular flows. More recently, Bhardwaj and Das [
12] developed an empirical flow pattern map for a two-phase thermosyphon using the condenser temperature and the wall superheat in the evaporator as coordinates. The obvious limitation of these approaches is they require the measurement of various quantities during the PHP operation, which is not simple, even in a laboratory environment; therefore, the relative flow pattern maps cannot be used as predictive tools. Moreover, the high spatial and temporal unsteadiness of the flow makes the choice of where and when such measurements should be carried out quite arbitrary, as well as the identification of suitable spatial and temporal averaging. Recent attempts to generate flow pattern maps for PHPs with the assistance of machine learning do not appear particularly successful to date [
13], as opposed to their promising use in the automatic flow regime recognition of standard two-phase flows [
3].
Historically, two-phase flow pattern maps have been constructed using two independent variables representative of the two phases, respectively, such as the superficial velocity of the liquid and the superficial velocity of the gas, or the corresponding Reynolds numbers. This approach is very common in the existing literature because flow patterns arise from the reciprocal distribution of the two phases within the channel, which in turn depend on the flow rates. However, due to the dynamic nature of pulsating heat pipes, which undergo rapid changes in the mass flow rate, pressure drop, rates of bubble growth and collapse, flow direction, and liquid film thickness, the superficial velocities of the liquid and vapour phases are unknown; therefore, they cannot be used as independent variables.
A less empirical approach to the construction of flow pattern maps consists of the characterisation of flow pattern transitions using suitable dimensionless numbers, such as the Bond, Weber, and Froude numbers, or modified versions of them, to account for the fluid acceleration and the bubble length [
14]. Specifically, the modified dimensionless numbers used in these PHP flow pattern maps were originally introduced to characterise the breakup of isolated vapour slugs during saturated flow boiling; therefore, they may not be suitable to describe other flow pattern transition mechanisms. In particular, these numbers require the calculation of the acceleration of the liquid–vapour interface in the direction of the flow, based on the measurement of the meniscus velocity as well as the measurement of the bubble length. Thus, their calculation is possible only when a well-defined plug or slug flow is observed, and it becomes elusive in other cases, such as annular or churn flow. An additional limitation of this approach is that, at present, there is no consensus on which dimensionless numbers should be used as independent variables in the construction of flow pattern maps, which makes it difficult if not impossible to compare the maps produced by different authors with one another. Finally, it must be observed that the remarkably large number of experimental data points reported in these works do not correspond to distinct sets of design and operation parameters because they are generated by taking several measurements at different times during the pseudo-steady state operation of the PHP.
In the present work, it is proposed to overcome the limitations of existing ex post approaches to construct empirical flow pattern maps for pulsating heat pipes by introducing two fictitious quantities, having the dimensions of a velocity, which are representative of the superficial velocities of the two phases inside a pulsating heat pipe. These quantities are defined only based on the design and operation parameters of the PHP and are used as the independent variables in the construction of a flow pattern map, which therefore can be used as a predictive tool in designing PHPs. The proposed approach is then validated using flow visualisation experiments that record the flow patterns observed in the adiabatic section of a flat-plate polypropylene pulsating heat pipe under a wide range of design and operation parameters.
2. Methodology
2.1. Concept
The main idea underpinning the proposed approach to the construction of flow pattern maps is to identify two independent variables, to be used as map coordinates, that are expressed in terms of geometric and operation parameters and in terms of the fluid thermophysical properties so that the map can be used as a truly predictive tool in PHP design and simulation. In fact, with few exceptions, the construction of flow pattern maps available in the open literature requires the direct or indirect measurement of one or more flow parameters, which poses a severe limitation to the possibility of using these maps as predictive tools.
The distribution of the liquid and vapour phases is mainly driven by their volumetric fluxes, i.e., their superficial velocities, which are unknown a priori. However, because the velocity of each phase is determined by the balance of driving forces and resistances, it is possible to define two virtual (or “ideal”) superficial velocities based essentially on dimensional arguments. In particular, in PHPs, the liquid phase is driven by capillary forces and/or gravity, and the main resistance to flow is viscus friction; meanwhile, the vapour phase is driven by the heat input and/or buoyancy, and the flow is opposed by viscous friction as well as the latent heat of the fluid. Thus, the virtual superficial velocities of the liquid and the gas phase can be defined as:
where
is the superficial velocity of the liquid,
is the superficial velocity of the vapour,
F is the filling ratio,
is the surface tension at the saturation temperature,
is the static contact angle,
is the hydraulic diameter,
and
are, respectively, the liquid and vapour densities at the saturation temperature,
is the PHP inclination angle,
is the dynamic viscosity of the liquid at the saturation temperature,
L is the length of the adiabatic section,
is the heat input,
N is the number of turns, and
is the latent heat of the fluid. The first additive term in Equation (
1) represents the relative importance of the capillary force with respect to the viscous friction, while the second term compares gravity and viscous friction. In Equation (
2), the first additive term is representative of the relative importance of the heat input in comparison with the latent heat, while the second term expresses the relative importance of buoyancy and viscous friction for the vapour phase. Note that in Equation (
1), the filling ratio is used as a correction factor for the viscous friction, while in Equation (
2), it also accounts for the sensible heat of the liquid.
The virtual superficial velocities defined by Equations (
1) and (
2) are therefore derived through a carefully constructed dimensional analysis that captures the fundamental driving forces and resistances encountered by liquid and vapour phases within the pulsating heat pipe. In particular, for the liquid phase, the velocity encapsulates the interplay between capillary forces and viscous friction, incorporating additional gravitational effects. Similarly, for the vapour phase, they construct a velocity that reflects the relative significance of heat input against latent heat while also accounting for buoyancy and viscous friction.
By using these virtual superficial velocities as coordinate axes, one can therefore create a flow pattern map that can predict flow regimes before actual PHP operation. In other words, the approach provides a semi-theoretical framework for understanding flow pattern transitions and establishing a start-up criterion for pulsating heat pipes. The strength of the proposed methodology lies in its ability to generate a predictive tool using only design and operational parameters, thereby circumventing the challenges of intrusive measurements and the high spatial and temporal unsteadiness characteristic of the flow inside pulsating heat pipes.
2.2. Experiments
The proposed flow pattern map was validated using two independent datasets generated from flow visualisation experiments in the adiabatic section of a flat-plate polypropylene PHP measuring 250 mm in length, 100 mm in width, and 1.5 mm in thickness, which was fabricated via selective transmission laser welding.
The first dataset was obtained from flow visualisation experiments on three prototypes, each with a hydraulic diameter of the meandering channel equal to 1.23 mm and featuring three turns, five turns, and seven turns, respectively. For each prototype, experiments were carried out with inclinations of 0°, 45°, and 90° with respect to the horizontal axis and heat inputs ranging from 10 W to 25 W. Both ethanol and perfluorohexane (FC-72) were used as heat transfer fluids. High-speed video recordings of the flow pattern in the adiabatic section of the PHPs were captured at 500 fps using a CMOS camera (MC1362, Mikrotron GmbH, Eching, Germany) equipped with an 18-108/2.5 zoom lens (Zoom 7000, Navitar Inc., Rochester, NY, USA). To enhance the contrast of the fluid interface, a diffuse, uniform LED back-to-front illumination was used. The liquid–vapour interfaces were further enhanced by a proprietary digital image processing algorithm running in the Matlab environment, which performed background subtraction and adjusted the digital gain and gamma correction of the resultant images.
The second dataset was obtained from flow visualisation experiments carried out during the 77th ESA Parabolic Flight Campaign on two prototypes featuring seven turns and eleven turns in the meandering channel, with hydraulic diameters in the meandering channel equal to 1.23 mm and 1.13 mm, respectively. For each PHP prototype, experiments were carried out under standard (1 g) and zero-gravity (0 g) conditions, featuring heat inputs ranging from 20 W to 35 W and using ethanol and FC-72 as heat transfer fluids. High-speed videos of the flow pattern in the adiabatic section of the PHPs were captured at 100 fps using both a CMOS camera (XiQ, Ximea GmbH, Münster, Germany) and an FLIR high-speed camera (FAST-IR M100K, Telops Inc, Québec City, Canada). The geometries of the four PHP prototypes considered in the present work are displayed in
Figure 1, while
Figure 2 shows the schematic layouts of the two experimental kits.
4. Conclusions
The innovative methodology proposed in this study represents a significant advancement in understanding and predicting flow patterns within pulsating heat pipes (PHPs). By introducing the novel concept of virtual superficial velocities for liquid and vapour phases, a predictive tool that transcends the limitations of existing approaches has been developed. Unlike previous methodologies that relied on intrusive measurements during PHP operation, this approach uniquely derives flow pattern characterisations solely from design and operational parameters.
The proposed approach has been validated on flat-plate polypropylene PHP prototypes with varying channel geometries, different working fluids (ethanol and perfluorohexane), and different operational conditions. The resulting flow pattern map revealed well-distinguished areas for the different flow patterns observed, unlike most of the existing maps.
An additional benefit of the proposed method is the formulation of an empirical start-up criterion, expressed as a linear relationship between virtual superficial velocities, which could prove useful at the stage of preliminary PHP design and optimisation.
Although further validation in a broader range of operational parameters and PHP designs is necessary, the proposed approach has a potential that extends beyond its current implementation, offering a promising foundation for predictive modelling in thermal management systems.