The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials and Industrial Machines

Yakov B. Gorelik and Artur H. Khabitov

## Abstract

Refrigeration systems based on free convection (two-phase thermosyphons) are used for cooling equipment units in chemical, nuclear power, and steel-making industries, as well as for thermal stabilization of natural materials with temperature-dependent properties, such as permafrost. Results of laboratory testing are reported for two types of thermosyphons applied mainly to the thermal stabilization of frozen ground: (a) vertical tubes with finning and (b) systems with horizontal evaporation tubes (HET systems). Their uses are currently restricted to relatively small thermal loads, but the effect of the loads on the cooling performance remains poorly investigated. Theoretical analysis of internal and external heat transfer in a vertical thermosyphon provides constraints on the boundary conditions at the evaporator wall, to be used in formulating and solving problems on the temperature regime of frozen ground stabilized with thermosyphons. Comparison of measured and calculated parameters that characterize the operation of a model HET system allows improving the calculation quality by applying the respective corrections.

Keywords: thermosyphons, refrigeration systems, free convection, two-phase flow, hydraulic resistance, cooling effect, experimental modeling, theoretical modeling, permafrost

## 1. Introduction

Refrigeration systems based on free convection (two-phase thermosyphons) are used for cooling equipment units in chemical, nuclear power, and steel-making industries, as well as for the thermal stabilization of natural materials with temperature-dependent properties, such as permafrost. Thermosyphons (also called "heat pipes" [1]) have different design features and sizes depending on application and operation conditions [1, 2]. Generally, two-phase thermosyphons consist of an evaporator and a condenser. The evaporator has a tightly closed case partly filled with liquid working fluid (coolant) and is placed in immediate vicinity of the cooled equipment or inside the cooled ground. The condenser contacts with a colder ambience (e.g., atmospheric air) and gives up heat. Heat flux to the evaporator boils up the fluid, and the released vapor moves to the condenser driven by the saturation vapor pressure gradient. The vapor condenses in the condenser, releasing heat into the environment, while the liquid phase of the fluid returns into the evaporator driven by gravity, capillary, or other forces, and the cycle repeats. The capillary liquid flow is provided by special microcellular lining on the inner walls of the thermosyphon. In the presence of constant sources of heat and cold, the system can operate successfully without other energy costs. The operation of thermosyphons has a large literature (see, for instance, an overview in [3]). Most publications deal with the conditions of high thermal loads (100 W/m or higher per tube unit length), when the heat transfer is especially active [4]. However, the conditions of low thermal loads (≤10 W/m) applied to the stabilization of permafrost, when the low kinetic energy of the two-phase flow hardly overcomes hydraulic resistance, have received much less attention. The action of the thermosyphons with the simplest construction (as vertical heat pipe) on permafrost is investigated mostly [5, 6]. The operation of some thermosyphon types, such as systems with horizontal evaporator tubes, has been poorly studied theoretically and experimentally. Nevertheless, such devices apply successfully on many constructions of northern areas [7–9].

To bridge the gap, model thermosyphons operated at low thermal loads, which simulate typical permafrost conditions, have been tested in the laboratory [10, 11]. Most thermosyphons used for thermal stabilization of permafrost run only during the cold season and are either vertical tubes or systems with horizontal evaporation tubes (HET). Fluid flow in such systems is commonly driven by gravity, as the use of capillary mechanisms has economic limitations. Where possible, the testing results are formulated as a boundary condition on the evaporator wall at the contact with the ambience, which allows exhaustive characterization of heat exchange between the system and the cooled natural or industrial materials and ensures correct formulation for the respective theoretical problem of external cooling.

## 2. Vertical thermosyphons

A vertical two-phase free convection system (thermosyphon) consists of relatively short (not over 10–15 m in length) tubes with an inner diameter of 30–50 mm [7, 12–14], based on Long's thermo-valve piling design [15]. The thermosyphons are inserted into the ground to most of their height and have a 1.5–2.0 m stick-up exposed to lower temperatures (Figure 1). Their operation commonly stops in summer, and special engineering solutions proceeding from thermal design are required to keep them running all year round or to increase the cooling effect (heat insulation, a certain number of tubes, etc.).

Two-phase free convention thermosyphons have been used to stabilize permafrost under kilometer-long aboveground pipelines (Trans-Alaska system or pipelines in northern, West, and Eastern Siberia), as well as buildings and utility structures in Russian Arctic oil and gas fields [13, 14, 16]. They are intended to cool down the ground to a designed temperature and keep this temperature stable for the whole lifetime of buildings and structures. Thermosyphons are most often installed prior to main construction, when foundations are being prepared. They are almost never used for the energy-consuming freezing of unfrozen parts of the ground, except for few cases of deeply buried systems for thermal stabilization of dam cores or for repairing thaw-related failures in the course of operation [8, 9].

The condenser can be equipped with heat absorption fins which are presumed to enhance the cooling effect, as it has been observed in practice. On the other hand, it follows from the theory [17] that fins should be mounted on the side of the medium with worse thermal properties. Special research is needed, however, to assess the potential efficiency of finning, given that heat transfer in thermosyphons is mainly

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

#### Figure 1.

A sketch of a two-phase thermosyphon: (1) tubular case, (2) fins, and (3 and 4) liquid and vapor phases of cooling fluid.

by convection above the ground and by conduction below the ground which has quite a low thermal conductivity. We have seen no evidence that would prove the thermal performance gain due specifically to fins.

#### 2.1 Condenser

Heat transfer in thermosyphons is modeled jointly for processes inside (fluid phase change and flow) and outside (heat exchange with ground and air) such systems.

The two phases of the working fluid in the tube move as follows [1, 2]. The vapor pressure can be assumed constant along the tube length, which is short in our case, and the hydraulic resistance to vapor flow toward the condenser is very small. Correspondingly, the saturation temperature ts (related to pressure via the saturation curve) is constant along both evaporator and condenser tubes. The condensed liquid flows along the inner tube as a thin film which thickens up progressively on its way to the evaporator. The film is the thinnest (zero) at the top of the condenser, is the thickest at the evaporator input, and then thins down upon evaporation while moving along the evaporator tube. In the simplest theoretical case, the hydraulic interaction of the two counterflowing phases can be assumed vanishing.

Heat exchange between the condenser and the air can be estimated easily [6, 18]. The total amount of heat Qk dissipated from the condenser to air per unit time is related to the temperature difference Δtka = tk–ta between the condenser wall (tk) and the incoming air (ta) as

$$Q\_k = a\_{tot} \cdot \mathbb{S}\_k \cdot \Delta t\_{st} \tag{1}$$

$$a\_{\rm tot} = \begin{cases} \frac{1}{\overline{\alpha\_f}} \cdot \frac{S\_k}{S\_{ik}} + \frac{1}{\alpha\_k} \end{cases}^{-1} \tag{2}$$

$$
\overline{a}\_{\overline{f}} = \frac{\lambda\_{\overline{c}}}{\overline{\delta}}.\tag{3}
$$

Note that equations for flat walls are used here and below, because the condensate film is very thin (δ <<R, where δ is the film thickness and R is the tube radius), according to observations on transparent physical models. In Eqs. (1)–(3), αtot is the total effective heat loss; Sk is the total area of the condenser outer surface; Sik is the total area of the condenser inner surface; α<sup>k</sup> is the heat exchange between the condenser and the air, without the film thermal resistance (found using known relationships via heat exchange of the condenser wall with air [6]); α<sup>f</sup> is the average heat transfer across the film; δ is the film thickness averaged over the condenser height; and λ<sup>c</sup> is the condensate thermal conductivity. The thermal resistance of the condenser wall is neglected as it has quite a high thermal conductivity while the wall is thin.

The average heat transfer across the sinking condensate film α<sup>f</sup> can be found from known depth dependence of the film thickness δ(z). At a constant temperature of the condenser wall tk, the depth-dependent film thickness δ(z) follows the Nusselt theory [19]. The equation for depth dependence of film thickness, corrected for heat exchange with air, neglecting vapor density (much smaller than liquid density) is given below for a condenser with the axis Oz directed down from the coordinates origin at the condenser top (Figure 1):

$$\int\_0^x q(z)dz = \frac{\kappa \cdot \mathbf{g} \cdot \rho\_c}{\mathfrak{Z} \cdot \nu\_c} \cdot \delta^3(z). \tag{4}$$

The heat flux density q(z) at the condenser refers to its inner surface, of the area Sik, and is

$$q(\mathbf{z}) = \frac{\Delta t\_{\rm sat}}{\delta(\mathbf{z})/\lambda\_{\rm c} + \mathbf{S}\_{ik}/\alpha\_{k} \cdot \mathbf{S}\_{k}}.\tag{5}$$

In Eqs. (4) and (5), ν<sup>c</sup> and ρ<sup>c</sup> are the kinematic viscosity and density of the condensate, respectively; κ is the heat of vapor-to-liquid phase transition; and g is the gravity acceleration.

Differentiation of both sides in expression (4) along z leads to the equation for δ(z), at δ(0) = 0:

$$\frac{\nu\_c \cdot \Delta t\_{sa}}{\kappa \cdot \mathbf{g} \cdot \rho\_c} = \left(\frac{\delta(z)}{\lambda\_c} + \frac{\mathbf{S}\_{ik}}{a\_k \cdot \mathbf{S}\_k}\right) \left(\delta^2(z) \cdot \frac{d\delta}{dz}.\tag{6}$$

Integration of this equation leads to a biquadratic algebraic equation with respect to δ(z):

$$\delta^4 + \frac{4 \cdot \delta^3}{3} \cdot \frac{\lambda\_c \cdot \mathbf{S}\_{ik}}{\alpha\_k \cdot \mathbf{S}\_k} - \frac{4 \cdot \lambda\_c \cdot \nu\_c \cdot \Delta t\_{\text{sat}} \cdot \mathbf{z}}{\kappa \cdot \mathbf{g} \cdot \rho\_c} = \mathbf{0}. \tag{7}$$

The solution of Eq. (7) can be written in elementary functions, but it would be cumbersome and inconvenient for further analysis. A more simple relation can be obtained for the small film thickness, with the first term in Eq. (7) much less than the second one:

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

$$
\delta \le \frac{4\cdot}{3} \cdot \frac{\lambda\_{\rm c} \cdot \mathbf{S}\_{ik}}{\alpha\_{k} \cdot \mathbf{S}\_{k}}.\tag{8}
$$

As follows from the comparison of relations (8) and (5) (factor 4/3 being insignificant), the film thermal resistance becomes negligible when inequality (8) fulfills. It means that the performance of highly efficient condensers (with large Sc and high αc, with a large dot product in the denominator of the right-hand side) is limited by the thermal resistance of the film, even commensurate with the resistance of the condenser. Therefore, the fulfillment of inequality (8) places constraints on the design of finned condensers in thermosyphons and in this respect is an economically and technologically important criterion. The respective requirement imposed on the condenser design has to be taken into account in further consideration.

With regard to the formulated requirement, we obtain the equation for the film thickness from Eq. (7), neglecting the first term in the left-hand side:

$$\delta(z) = \left(\frac{\mathbf{3} \cdot \nu\_{\rm c} \cdot \Delta t\_{\rm st} \cdot z}{\kappa \cdot \mathbf{g} \cdot \rho\_{\rm c}} \cdot \frac{a\_{\rm k} \cdot \mathbf{S}\_{\rm k}}{\mathbf{S}\_{\rm ik}}\right)^{\frac{1}{\mathbf{\tilde{s}}}}.\tag{9}$$

The film is the thickest on the condenser bottom, at z = l<sup>1</sup> (Figure 1). If inequality (8) fulfills at the maximum film thickness, it obviously fulfills for any thinner film. From relations (8) and (9), it follows that

$$\frac{\alpha\_k \cdot \mathbf{S}\_k}{\mathbf{S}\_{ik}} \cdot \left(\frac{\mathbf{S}^4 \cdot \nu\_c \cdot \Delta t\_{st} \cdot l\_1}{\mathbf{4}^3 \cdot \lambda^3 \boldsymbol{\iota} \cdot \mathbf{\kappa} \cdot \mathbf{g} \cdot \rho\_c}\right)^{\frac{1}{4}} \le \le 1. \tag{10}$$

The left-hand side of inequality (10) is dimensionless and depends on the fluid (coolant) type, as well as on the condenser design. The inner and outer surface areas of an unfinned condenser are equal (the wall thickness being neglected), and the first factor in the left-hand side is equal to the heat loss α<sup>a</sup> for unfinned (smooth-walled) tubes in the case of cross flow around tubes defined by the known relationship from the specified wind speed [20]. Let the left-hand side of inequality (10) be denoted as Φ with subscripts that refer to different coolants. For R = 16 mm, l<sup>1</sup> = 1.5 m, and the wind speed 3 and 8 m/s, the α<sup>a</sup> values are 26.4 and 58.3 W/(m2 �deg), while <sup>Φ</sup> values are <sup>Φ</sup>am <sup>=</sup> 8.26 � <sup>10</sup>�<sup>3</sup> , 1.83 � <sup>10</sup>�<sup>2</sup> for ammonia, <sup>Φ</sup>fr-12 <sup>=</sup> 1.50 � <sup>10</sup>�<sup>1</sup> , 3.33 � <sup>10</sup>�<sup>1</sup> for Freon-12, <sup>Φ</sup>fr22 <sup>=</sup> 1.04 � <sup>10</sup>�<sup>1</sup> , 2.31 � <sup>10</sup>�<sup>1</sup> for Freon-22, <sup>Φ</sup>CO2 <sup>=</sup> 2.63 � <sup>10</sup>�<sup>2</sup> , and 5.79 � <sup>10</sup>�<sup>2</sup> for CO2 (at the temperature difference about the maximum Δtsa = 20°C). Thus, inequality (10) fulfills in all cases in the absence of fins (though not perfectly for Freons at high wind speed), and the contribution of film thickness can be neglected for all fluid considered above. On the other hand, fins on the outer surface of the condenser improve external heat exchange, but their enhancement loses sense when the thermal resistance of the film becomes significant.

For a finned condenser, the parameter Φ depends on thickness, spacing, and other parameters of fins (see the fin thickness dependence of Φ in Figure 2 for different fluids but the same fin parameters and a wind speed of 5 m/s). For the chosen condenser design, inequality (10) fulfills well for all fin thicknesses with ammonia and slightly worse with carbon dioxide. With Freons, only very thin fins can be efficient, but they may be problematic to fabricate. The effect of the thermal resistance of a condensate film is controlled by its thickness and the thermal conductivity of the liquid phase in different coolants. Variations of film thickness at the condenser output as a function of fin thickness for ammonia and Freon-12

Figure 2.

Behavior of parameter Φ as a function of fin thickness (h) for ammonia (1), CO2 (2), HCFC (Freon)-22 (3), and HCFC (Freon)-12 (4) fluids, at R = 16 mm and l<sup>1</sup> = 1.5 m.

(Figure 3) show that ammonia films are at least 5 times thinner, have about 8 times higher thermal conductivity, and, hence, have about 40 times lower thermal resistance.

As relation (8) fulfills, Eq. (2) for the total heat loss from the condenser becomes simpler:

$$a\_{tot} = a\_k.\tag{11}$$

#### 2.2 Evaporator

The thickness of the sinking condensate film continuously decreases from the maximum δ<sup>0</sup> at the evaporator input, which coincides with that at the condenser output, and can be estimated by expression (9) at z = l1. The depth-dependent film thickness variations in the evaporator likewise follow the Nusselt theory, applied to evaporation instead of condensation; at a constant temperature on the evaporator wall, tg is equal to the current temperature of the ambient ground. The previously estimated [21] fluid flow rates inside the system largely exceed the rate of ground

#### Figure 3.

Film thickness (δ) inside a condenser as a function of fin thickness (h) for ammonia (1) and HCFC (Freon)-12 (2) fluids.

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

temperature redistribution. Therefore, the formation and flow of the condensate film can be considered as a quasi-stationary process at the above conditions. In this case, the heat flux to the evaporator q(z) is

$$q(\mathbf{z}) = \frac{\lambda\_{\mathbf{c}} \cdot \Delta t\_{\mathbf{g}}}{\delta(\mathbf{z})}.\tag{12}$$

The equation of the evaporator heat budget (similar to Eq. (4)) makes basis for the differential equation with respect to film thickness. With the initial condition δ = δ<sup>0</sup> at z = l1, the equation for δ(z) becomes

$$\delta^4(\mathbf{z}) = \delta\_0^4 + \frac{4 \cdot \lambda\_\varepsilon \cdot \Delta t\_\text{g} \cdot \nu\_\varepsilon \cdot (l\_1 - z)}{\kappa \cdot \rho\_\varepsilon \cdot \mathbf{g}}; l \succeq \mathbf{z} \succeq l\_1,\tag{13}$$

where l is the total tube length. According to expression (13), the nonzero film thickness is limited by some natural depth reasonably corresponding to the maximum evaporator length (see below).

Eq. (13), with averaging along z, can be used to calculate the average heat loss over the evaporator αif . Note that the parameters ts and tg remain unconstrained in all above equations. The saturation temperature ts is estimated using the total heat budget for the system as a whole; heat coming to the evaporator equals that dissipating to the air from the condenser; Qk = Qf, where Qk is defined by relation (1), while the heat from the evaporator is Qf ¼ αif � Sf � Δtgs (where the evaporator surface area Sf = 2�π�R�l2 and its length is l2). Solving the heat budget equation with respect to ts gives

$$t\_s = \frac{\mathcal{S}\_f \overline{\alpha}\_{\circ f} \cdot t\_\mathcal{g} + \mathcal{S}\_k \cdot \alpha\_k \cdot t\_\mathcal{a}}{\mathcal{S}\_f \overline{\alpha}\_{\circ f} + \mathcal{S}\_k \cdot \alpha\_k}. \tag{14}$$

With Eqs. (12) and (13), the average heat loss from the evaporator αif likewise depends on the temperature ts and is given by

$$
\overline{a}\_{\circ \circ} = \left(\frac{\lambda\_{\circ}^{\circ} \cdot \kappa \cdot \rho\_{\circ} \cdot \mathbf{g}}{4 \cdot \nu\_{\circ} \cdot l\_{2} \cdot \Delta t\_{\mathfrak{g}^{\circ}}}\right)^{0.25}.\tag{15}
$$

To write the equation for ts, we first define the difference Δtgs ¼ tg � ts using expression (14), then pick the terms with the difference Δtgs, and obtain the transcendent algebraic equation with respect to this value:

$$
\Delta t\_{\rm gs} + \frac{3}{4} \cdot \frac{\text{S}\_{\rm f}}{a\_k \cdot \text{S}\_k} \cdot \left(\frac{\lambda\_\text{c}^3 \cdot \kappa \cdot \rho\_\text{c} \cdot \text{g}}{4 \cdot \nu\_\text{c} \cdot l\_2}\right)^{0.75} \cdot \Delta t\_{\rm gs}^{0.75} = \Delta t\_{\rm ga}, \tag{16}
$$

which is solved numerically. In the range of ground temperatures tg from 0°C to ta, the temperature ts is close to tg to thousandth fractions for ammonia and 0.6°C for Freon-12.

The difference Δtgs becomes important for estimating the maximum length l ∗ <sup>2</sup> of the evaporator tube and cannot be neglected. At δ = 0, the maximum tube length is obtained from expression (13) as

$$d\_2^\* = \delta\_0^4 \cdot \frac{\boldsymbol{\kappa} \cdot \boldsymbol{\rho}\_c \cdot \mathbf{g}}{\mathbf{4} \cdot \boldsymbol{\lambda}\_c \cdot \Delta t\_{\text{g}^s} \cdot \nu\_c}. \tag{17}$$

Let the heat loss from the condenser in Eq. (9) be α<sup>k</sup> � Sk=Sik � F (see its behavior as a function of fin thickness in Figure 4) with the values 40, 200, and 400 W/ (m<sup>2</sup> �deg). Then, l ∗ 2= 1.7, 6.0, and 11.9 m, respectively, for ammonia and can reach 200 m for Freon-12. This length may be a major advantage of Freon evaporators in thermosyphons used for cooling deep ground.

In principle, evaporator tubes in thermosyphons used in deep ground can be longer than the estimated maximum value, but they should be filled with liquid coolant till the upper limit. In this case, heat transfer should be considered with regard to a shift in the temperature of boiling (in the lower tube section) under the effect of hydrostatic pressure. However, we have no reliable evidence of thermosyphon operation in these conditions.

The reported estimates for the maximum evaporator length require further experimental checks. The issue is especially important because one of us repeatedly observed sporadic ebullition of fluid at the evaporator bottom during laboratory testing of vertical thermosyphons in physical models with transparent walls. The vapor-liquid flow from such ebullition wets the evaporator walls and reaches the condenser, while the evaporator temperature falls abruptly, and moisture condenses on its outer wall surfaces. This effect, which may play a significant role, was demonstrated in a video at TICOP [22], but it is neglected in this consideration. Laboratory testing of thermosyphons has to be continued. The sinking film moves wavelike upon interaction with the down-going flow of vapor [23], but the wavelike flow is stable and even improves the heat exchange to some extent [19].

#### 2.3 External problem for thermal stabilization

The distribution of ground temperatures in the zone of thermal stabilization is estimated as follows. The key issue is to constrain the boundary conditions on the evaporator outer wall which contacts the ground. For this, the saturation temperature (ts) is expressed from the equation of total heat budget. The heat loss from the condenser Qk is defined by the same Eq. (11), while the heat coming to the evaporator Qf from the ground is

$$
\Delta a\_k \cdot \mathbf{S}\_k \cdot \Delta t\_{\rm sat} = \mathbf{S}\_f \cdot \lambda\_f \cdot \left(\frac{\partial t}{\partial r}\right)\_R,\tag{18}
$$

Figure 4. Behavior of parameter F as a function of fin thickness (h); R = 16 mm.

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

where λ<sup>f</sup> is the thermal conductivity of frozen ground and t is the ground temperature depending on the radial coordinate r (in the axisymmetric case) and the time τ. With the Eq. (14) for ts and inequality (8), the sought condition is

$$\frac{\alpha\_{\text{k}} \cdot \text{S}\_{\text{k}}}{\text{S}\_{\text{f}}} \cdot \left(t\_{\text{a}} - t\_{\text{g}}\right) = -\lambda\_{\text{f}} \cdot \left(\frac{\partial t}{\partial r}\right)\_{\text{R}}.\tag{19}$$

The same equation was obtained earlier for an unfinned condenser [21], with the only difference in the factor F in the left-hand side of condition (19), before the temperature difference (in the absence of fins, Ftot = αal1/l2). Eq. (19) is a thirdorder boundary condition on the evaporator outer wall (tg � t(R)), which fully characterizes the cooling effect in problems for temperature variation in permafrost stabilized by vertical thermosyphons. The condition allows formulating and solving a large scope of problems concerning the temperature field maintained with any number of vertical thermosyphons in construction. This is the most important result of this part.

## 3. HET systems

Thermosyphons with horizontal evaporation tubes allow using pad foundations instead of piling in some types of structures, which saves labor costs. These systems have been largely used for permafrost stabilization in terms of works designed mainly by the R&D companies Fundamentstroiarkos (direct-flow HET systems) and Fundamentproekt (counter-flow systems with inclined evaporator tubes) [8, 24], but no methods are still available for estimating their efficiency. The known theoretical studies characterize the operation of HET systems at high thermal loads of tens or hundreds W/m [4], while the loads in thermal stabilization of permafrost do not exceed 3–5 W/m [11]. As shown by laboratory testing at small thermal loads [10, 11], the two phases of the working fluid have commensurate volumetric contents along the tube length and interact mechanically. This fact does not allow simple description of fluid dynamics and heat exchange in the two-phase flow (e.g., [25–28]). A pure theoretical description of these processes is hardly possible at the time being and requires updating with reference to testing results for specific systems. The original design of HET systems implied unidirectional flow of both phases all along the tubes, which reduces hydraulic resistance to the coolant motion and improves the system performance.

#### 3.1 Methods and results of laboratory testing

The operation of a laboratory model of an HET system is discussed below for a general experiment layout as in Figures 5 and 6. The laboratory model was described in detail earlier [10, 11]. The test equipment includes temperature sensors (thermistors), an electronic vacuum meter, a level gauge, and an automatic recorder (Figure 5). The condenser is made of metal and placed in a large cooling chamber, while the evaporator consists of glass tubes connected by chemically inert rubber hoses and laid horizontally on wooden pads on the room floor. The transparent glass tubes make visible the flow behavior in different segments of the system, which is important for documenting its operation. Acetone was used as a cooling fluid as its saturation vapor pressure is below the ambient pressure within the applied temperature range (unlike most coolants used in the industry). This ensures cheap, straightforward, and safe handling of the model but requires

#### Figure 5.

A sketch of a laboratory test equipment. 1–7 are thermistors.

#### Figure 6.

Structural scheme of the model system (symbols are explained in the text).

pumping out atmospheric gases which are present inside the system and interfere with its work.

The model of the HET system (Figure 6) consists of a condenser (1) with a narrower bottom part which accommodates the condensed liquid (2) and an evaporator (3). The vertical evaporator segment leading to the condenser input (segment CК in Figure 6) is made of a nontransparent insulating material, which masks the fluid flow. The flow behavior within this segment can be judged from examination of the neighbor glass horizontal segment C0C inside the cooling chamber. When the system operates normally, both fluid phases flow in the same direction, toward the condenser input К, while the share of vapor is quite large. The flow of each phase toward point C within the main tube length undergoes successively the plug, wave, and projectile regimes.

The system operates normally as long as the tube length-to-diameter (L/D) ratio remains ˜500 (in industrial systems, it may reach 5000 or more). As this ratio increases, the operation becomes progressively less stable, while the cooling effect (temperature difference between the evaporator and the room air) reduces. On the other hand, some amount of condensate liquid always flows back into the

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

evaporator during the normal operation. This return flow is visible from outside the cooling chamber [11], except for the nontransparent segment C0C, and thus forms apparently all along C0CK (Figure 7). The evaporator cools down till the lowest temperature during normal operation, which is recorded by thermistors (Figure 8) as distinct from the unstable one (Figure 9). The flow of phases within segment C0C at normal (stable) regime can be either laminar or turbulent (Figures 10 and 11).

During the unstable operation, till failure, the liquid phase moves very slowly (visually seeming immobile) within the segment C0C, while the vapor phase is restricted to bubbles of different sizes (Figures 12 and 13) which move relatively slowly toward the condenser input. The cooling effect in this case is often oscillatory and fails to reach the designed temperature at the normal regime (Figure 9).

#### Figure 7.

Two-phase flow near the condenser input (segment C0CK): Dark and light shades are condensate and vapor; heavy and thin arrows show vapor and liquid flows, respectively; dash line is the boundary between the two flows.

#### Figure 9.

Thermistor readings along the evaporator tube at unstable operation of the system.

Figure 10. Laminar flows within segment C0C at stable operation.

Figure 11. Turbulent flows within segment C0C at stable operation.

Cooling is the weakest at the smallest amount and slowest flow of vapor within C0C, that is, the system actually fails. In this case, the condensate level H in the condenser (Figure 6) reaches its maximum [11], and constant return flow is visible from outside the evaporator, like during the normal operation. Thus, the return

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

Figure 12. Small vapor bubbles at unstable operation.

Figure 13. Big vapor bubbles at unstable operation.

flow along C0CK is a typical feature of fluid behavior (at least at small thermal loads on the evaporator). Other features (operation instability, slow start, poor cooling effect) become more prominent at longer evaporator tubes (greater L/D ratios) but are minor at short tubes.

The evaporator tubes of full-size HET systems are commonly fully filled with the liquid phase, while the condenser is filled partly, till the specified level H (relative to the evaporator position in the horizontal plane). However, the effect of the fluid volume on the system performance at small thermal loads remains poorly investigated. On the other hand, systems with a lower fluid volume are cheaper and pose lower hazard to environment in case of emergency spilling. The respective laboratory tests were performed also with a partly filled evaporator and an empty condenser (β = L0/L < 1, where L0 is the tube length filled with condensate and L is the total tube length), in order to initially increase the vapor share in the vapor– liquid mixture and thus reduce the hydraulic resistance and improve the system performance. The measured stationary temperature distribution along the tube for different H levels (Figures 14–17) demonstrates that the cooling effect reduces notably as the fluid volume increases.

Main testing results can be summarized as follows: (i) there is always a return flow of liquid sinking into the evaporator opposite to the main flow of the twophase mixture in the segment near the condenser input (segment C0CК in Figure 6) [11]. The presence of this return flow differs the real fluid behavior from the theoretical idea of unidirectional flow of the fluid phases along the whole circulation path (at least for small thermal loads); (ii) at low thermal loads, the systems with return flow develop a slug and plug regime that leads to failure [25]. Specifically, the condensate plugs increase the hydraulic resistance and obstruct the flow to the

#### Figure 14.

Temperatures at thermistor points (1–7) and excess pressure (P) during the operation of the model at β = 0.5.

#### Figure 15. Measured (1) and calculated (2, 3) temperatures along the tube at β = 0.5, ω = 1 (2), and ω = 50 (3).

Figure 16. Measured (1) and calculated (2, 3) temperatures along the tube at β = 1.0, ω = 1 (2), and ω = 70 (3).

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

Figure 17. Measured (1) and calculated (2) temperatures along the tube at standard H level; ω = 1 (2).

condenser within this segment. The flow behavior may change frequently within a single test, which leads to oscillations of the evaporator temperatures (Figure 9). This fluid flow behavior at small thermal loads has to be taken into account in calculations for estimating the system performance.

The experimental uncertainty analysis of the experimental data was not done yet. It is planned for the next publications.

#### 3.2 Calculations for thermosyphon design

According to the existing methods of fluid-dynamic calculations [25–28], the state of a two-phase flow in each tube cross section depends on the flow rates of vapor Gw(z) and condensate Gl(z), in kg/s and volumetric vapor content αw(z), where z is the coordinate along the evaporator tube, with the origin at point O (Figure 6). The condensate flow rate is expressed via the total flow rate G: Gl(z) = G-Gw(z). In the stationary process, G is constant along the flow, but it is unknown a priori and has to be calculated.

The function αw(z) in the simplest model of Lockhart-Martinelli [26, 27] is defined by Gw(z) and Gl(z). The fluid flow along the tube is split into two parts: heating and two-phase flow. In the heating part, from the point O to the boiling point B at zb (Figure 6), αw(z) = 0, and the coordinate zb is found from the known solutions for temperature and pressure for the single-phase flow at a given lateral heat transfer and from the condition at which the flow parameters reach the fluid saturation at the given point. In the simplest case of constant heat flux q along the evaporator wall (W/m), the temperature change between O and zb is

$$\mathbf{t}(\mathbf{z}) = \mathbf{t}\_t + \mathbf{q} \cdot \mathbf{z}/\mathbf{G} \cdot \mathbf{c}\_l \tag{20}$$

where ts and cl are the fluid temperatures at the output of the cooling chamber (point O, Figure 6) and the specific heat of the liquid phase, respectively. At the liquid–vapor equilibrium, the saturation pressure Ps correlates with the temperature of the phases t, and the Ps(t) function can be determined empirically as a threeparameter curve, to a sufficient accuracy. Here we use the quadratic equation:

$$P\_s(t) = b\_1 + b\_2 \cdot t + b\_3 \cdot t^2,\tag{21}$$

with the coefficients b1 = 10<sup>4</sup> Pa, b2 = 600 Pa/°C and b3 = 12.5 Pa/(°C)2 for an acetone fluid. The z-dependent saturation pressure can be expressed as in Eq. 20: Ps(z) = Ps(t(z)). The liquid flow within the OB segment is laminar, and expression for hydrodynamic pressure Ph(z) is given by

$$P\_h(\mathbf{z}) = P\_k + \rho\_l \cdot \mathbf{g} \cdot H - \frac{128 \cdot \mu\_l \cdot \mathbf{G} \cdot (H + z)}{\pi \cdot \rho\_l \cdot D^4},\tag{22}$$

where Pk = Ps(tk); tk and Pk are, respectively, the pressure and temperature of saturation vapor inside the condenser; ρ<sup>l</sup> and μ<sup>l</sup> are the density and dynamic viscosity of the condensate; and D and g are the inner diameter of the tube and the acceleration due to gravity. The position of the boiling point B (zb) is found from the equation Ps(z) = Ph(z) (the respective equation is omitted). Both tk and Pk inside the condenser are a minimal value of temperature and pressure for all two-phase segment BCK.

In the two-phase flow part, between zb and the condenser input K (Figure 6), the functions Gw(z) and αw(z) increase monotonically from zero at zb, while the z dependence is calculated according to the specified conditions of heat transfer on the evaporator wall. In the current model, all incoming heat is spent on fluid evaporation, and Gw(z) is

$$G\_w(z) = q(z - z\_b)/\kappa. \tag{23}$$

Note that all variables depend on G, H, H0, and other constructive and external parameters (e.g., on the temperature in the cooling chamber ta). The pressure gradient within the two-phase flow segment dP/dz is a sum of three components (see [27, 28] for their expressions), dP/dz = (dP/dz)f + (dP/dz)a + (dP/dz)g, which refer to friction, flow acceleration, and gravity, respectively. The gravity component disappears within the horizontal tube, and the acceleration component disappears within the vertical tube segment CK (with the assumption of ideal insulation of the inlet tube). Integration of the dP/dz equation along the segment between the condenser top (where the pressure is Pk) and condenser input K leads to zero total pressure, and the sought variables can be found by equation

$$
\oint \frac{dP(z, G, H)}{dz} dz = 0.\tag{24}
$$

In this equation, P depends explicitly on the parameters G and H, which is relevant to the further consideration. The variables tk and ts in Eqs. (20) and (22) can be found from the equation for condenser-air heat exchange (with regard to the condenser design and processes inside it) and the integral balance equation for the stationary operation of the system (similar to Eq. (14). In this consideration though, focused on the applicability of the existing calculation methods to small thermal loads, the temperatures are found in a simpler empirical way: from thermistor readings (#1 and #7 in Figure 5 for ts and tk, respectively).

If the condenser design ensures invariable H at any change inside the evaporator during the system operation (i.e., H is a specified and known value in the calculations), Eq. (24), together with known ts and tk, is sufficient to find all process parameters. In this case, integration in relation (24) leads to an algebraic equation with respect to G which has a simple numerical solution. Once G has been found, all flow parameters can be found as well: zb, Gw(z), αw(z), Ps(z), and the temperature distribution within the two-phase flow part which is related with the saturation pressure Ps(z) as in expression (21). The temperature distribution along the tube

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

should be a continuous function of z, given the way of zb determination. It is convenient in our case to divide the friction contribution to the integral of Eq. (24) into two components (IfL from the origin to C0 and IfH0 along C0CK), as high hydraulic resistance within C0CK is especially important for this consideration.

If the evaporator tube is filled to a part of its length (β ≤ 1), the initial fluid level is H = 0, while H is set unprompted and is unknown (to be calculated) in operant condition. For this, an additional equation is used, which is a mathematical form of the mass conservation law for fluid in the system, likewise including two unknowns (G and H), assuming equal tube diameters:

$$(\mathbf{1} - \boldsymbol{\beta}) \cdot \mathbf{L} + (\mathbf{1} - a\_w(\mathbf{L})) \cdot H\_0 + H = \int\_{x\_b}^{L} a\_w(z) dz,\tag{25}$$

where zb and αw(z) parametrically depend on G and H. Thus, Eqs. (24) and (25) make up a system of algebraic equations with respect to G and H. Solving them jointly yields the unknowns and the temperature along the tube at arbitrary external parameters (temperature in the cooling chamber, heat transfer at the evaporator wall, fluid properties, and system specifications). Cooling of the evaporator relative to its ambience (room air in the case of laboratory testing) is the main parameter of its heat exchange with the host material.

#### 3.3 Calculated and measured data compared

The calculated and measured data can be compared with an example of three laboratory tests at β = 0.5, β = 1, and at a standard fluid level (H = 1.6 m). Other experimental conditions are identical: cooling chamber temperature tc = �12.5°C and room air temperature ta = +23°C, L = 17.6 m, and D = 0.01 m (L/D�2000). The thermal load on the evaporator is assumed to be q = 2.8 W/m proceeding from free convection as the basic mechanism of heat exchange between the evaporator and the ambient air. Time-dependent temperature variations recorded by all thermistors and saturation pressure for the run with β = 0.5 (Figure 14) allow estimating the values tk = 1.5 and ts =5(°C) in the case of a quasi-stationary process (similar records for the other runs are not shown). The measured evaporator temperatures for the quasi-stationary process approximated by curve 1 in Figure 15 differ markedly from those calculated (curve 2) with Eqs. (24) and (25). To improve the fit, the calculated values require correction by a factor of ω (at the integral IfH0), which increases hydraulic resistance within the segment C0CK. A good fit between measured and calculated temperatures can be achieved with a 50 times higher hydraulic resistance or at ω = 50 (curve 3). In the second run (β = 1.0, Figure 16), the fit improves at ω = 70. This correction for higher hydraulic resistance within the segment C0CK is valid for partly filled evaporators with the filled-to-total tube length ratios (β) from 0.3 to 1.0. According to preliminary estimates, the correction factor ω in this β range can vary from a few to hundreds of times, depending on variations of other parameters. In the case of standard H level (run 3, Figure 17), the fit is good without correction (ω = 1).

### 4. Conclusions

Improving the performance of vertical thermosyphons by increasing the surface area of fins is limited by the internal thermal resistance of the condensate film that sinks down the condenser inner walls. The thermosyphon-ground heat exchange

can be described by the third-order boundary condition, with the respective problem formulation.

Laboratory testing shows that the model of an HET system in the tested design can operate with both partly and fully filled evaporator tubes. The calculations in the former case should additionally include the mass conservation equation for the working fluid. Calculated and measured data fit well with a correction for higher hydraulic resistance in the flow along the vertical evaporator segment leading to the condenser input. The modified calculations remain semiempirical through and require further updating.

## Acknowledgements

The authors wish to thank G.M. Dolgikh, director of the Fundamentstroiarkos R&D company, and S.N. Okunev, chief engineer of the same company, for the discussions and constructive criticism of this subject matter.

## Notes

The work was made in accordance with Program of Fundamental Researches RAN IX. 135.2, Project IX. 135.2.4.

## Author details

Yakov B. Gorelik\* and Artur H. Khabitov Earth Cryosphere Institute, Siberian Branch of the Russian Academy of Sciences, Russia, Tyumen

\*Address all correspondence to: gorelik@ikz.ru

© 2018 The Author(s). Licensee IntechOpen. This chapteris distributed underthe terms oftheCreative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The Laboratorial Research of Two-Phase Free Convection Devices for Cooling of Materials… DOI: http://dx.doi.org/10.5772/intechopen.80789

## References

[1] Dunn PD, Reay DA. Heat Pipes. Oxford: Pergamon Press; 1978. p. 334

[2] Pioro IL, Antonenko VA, Pioro LS. Efficient Heat Exchange Systems with Two-Phase Thermosyphons. Kiev: Naukova Dumka; 1991. p. 246 (in Russian)

[3] Vasil'ev LL, Konev SV. Heat pipes. Heat Transfer (Soviet Research). 1974; 6(1):1-105

[4] Kutepov AM, Sterman AS, Stiushin AG. Fluid Dynamics and Heat Transfer in Vapor Formation. Moscow: Vysshaya Shkola; 1986. p. 448 (in Russian)

[5] Khrustalev LN. Principles of the Geotechnics at Cryolithozone. Moscow: MSU; 2005. p. 554 (in Russian)

[6] Gorelik JB, Seleznev AA. Efficiency of finning in short vertical two-phase thermosyphons for construction on permafrost. Earth's Cryosphere. 2016; XX(2):78-89

[7] Dolgikh GM, Okunev SN. Construction and stabilization of cooled foundations in permafrost by Fundamentstroiarkos company. In: Proceedings of the International Conference on Theory and Practice of Permafrost: Assessment and Predictions; 29-31 May 2006; Russia. Tyumen: TNGGU; 2006. pp. 228-232

[8] Dolgikh GM, editor. Systems for Thermal Stabilization of Permafrost. A Collection of Papers by Specialists from Fundamentstroiarkos R&D Company for 2010–2014. Novosibirsk: Geo; 2014. 218 p (in Russian)

[9] Kutvitskaya NB, Minkin MA. Design of foundations for infrastructure of oil, gas, and condensate field in complex permafrost conditions. Osnovaniya, Fundamenty, Mekhanika Gruntov. 2014;1:21-25

[10] Gorelik JB, Gorelik RJ. Laboratory research of work of the natural convectional biphase cooling device with horizontal evaporative part. Kriosfera Zemli. 2011;XV(2):34-43

[11] Gorelik JB. Operation of two-phase thermosyphons with horizontal evaporator tubes: Causes of instability. Earth's Cryosphere. 2015;XIX(4):81-93

[12] Vasil'ev LL, Vaaz SL. Freezing and Heating of Ground with Thermosyphons. Minsk: Nauka i Tekhnika; 1986. 192 p (in Russian)

[13] Vyalov SS, editor. Ground Temperature Stabilization with Thermosyphons. Yakutsk: IM SO RAN; 1983. 124 p (in Russian)

[14] Makarov VI. Thermosyphons in High-Latitude Construction Engineering. Novosibirsk: Nauka; 1985. 169 p (in Russian)

[15] Long EL. Means for Maintaining Permafrost Foundations. Patents USA, Cl. 165–45; 1964. No. 3, 217, 791

[16] Dolgikh GM, Dolgikh DG, Okunev SN. Freezing of soils: Engineering solutions of Fundamentstroiarkos company. In: Proceedings of the International Conference on the Cryosphere of Petroleum Provinces; 26-28 May 2004; Russia. Tyumen: TNGGU; 2004. pp. 56-57

[17] Arnold LV, Mikhailovskiy GA, Seliverstov VM. Engineering Thermodynamics and Heat Transfer. Moscow: Vysshaya Shkola; 1979. p. 445 (in Russian)

[18] Roizen LI, Dulkin IN. Thermal Design of Finned Surfaces. Moscow: Energiya; 1977. p. 255 (in Russian)

[19] Mikheev MA, Mikheeva IM. Fundamentals of Heat Transfer.

Moscow: Energiya; 1973. 320 p (in Russian)

[20] Wong HY. Handbook of Essential Formulae and Data on Heat Transfer for Engineers. London, New York: Longman Group Ltd; 1977. p. 236

[21] Gorelik JB. Modeling ground temperature field around vapor-liquid thermal piles. Problemy Nefti i Gaza Tyumeni. 1980;47:58-61

[22] Gorelik JB. Laboratory study of the devices for thermostabilization of frozen ground. Initiative oral reporting. In: Program of 10th International Conference on Permafrost (TICOP); 25-29 June 2012; Russia. Salekhard: Pechatnik; 2012. p. 18

[23] Kapitsa PL. Wave flow of thin liquid layers. Zhurnal Teoreticheskoy Fiziki. 1948;18(1):3-18

[24] Minkin MA. Basements and Foundations in Permafrost. Moscow: GASIS; 2005. 213 p (in Russian)

[25] Wallis G. One Dimensional Two-Phase Flow. New York: McGraw-Hill; 1969. p. 440

[26] Butterworth D, Hewitt GF. Two-Phase Flow and Heat Transfer. Oxford University Press (Harwell Series); 1977. 514 p

[27] Crowe CT, editor. Multiphase Flow Handbook. London—New York: Taylor & Francis; 2006. 1128 p

[28] Labuntsov DA, Yagov VV. Mechanics of Two-Phase Systems. Moscow: MEI; 2000. 374 p (in Russian)

## Chapter 5
