_{1}

^{*}

The results of the numerical studies of vortex formation inside short heat pipes (HP’s) with profiled vapour channel in the Laval-liked nozzle form are presented. For the first time, it was found that the vapour vortex of moist compressible vapour flow in the cooled part of vapour channel changes its rotational motion direction. The rotation direction of the toroidal vapour vortex, obtained by solving the Navier Stokes equations is dependent on the heat power value, entering to the HP’s evaporator. With low heat power loads the rotational direction of the circular toroidal vapour ring due to the Coanda effect and sticking moving vapour jets to the channel’s walls occurs from the periphery to the longitudinal axis of the vapour channel. While the heat power load increasing, the direction of the circular toroidal vapour ring rotation changes to the opposite, from the longitudinal axis to the periphery of the vapour channel. The thickness of the formed working fluid condensate film located under the toroidal vapour vortex also related to the evaporator heat power load and the associated toroidal vapour vortex rotation direction. The numerical thickness calculation of the formed working fluid condensate film located under the toroidal vapour vortex was compared with experimental values, obtained by capacitive sensors. The thickness values of the calculated condensate film thickness and experimentally measured values using capacitive sensors are close in magnitude order.

Heat pipes (HP’s) are very interesting heat transfer devices, currently widely used in almost all cases of heat-stressed cooling of spacecraft and satellites, military facilities and household devices and computers. Together with the wide practical application, the HP’s are of great interest for the scientific study of moist vapour currents, the occurrence of pulsations in the vapour channel during the intermittent vaporization operation of the capillary-porous evaporator, the formation of toroidal vortices near the cooled condensation surface and many other interesting effects [

In this study was applied stainless steel HP’s with shaped vapour channel in the converging-diverging nozzle form of nonlinear varying diameter and with flat top and bottom covers [

All details of the HP’s manufacture, profiled vapour channel formation along the entire HP’s length and hydraulic diethyl ether delivery system to the evaporator, design and calibration of capacitance measuring sensors, development of high-frequency electronic generators and measuring devices are given in the [

The dimensions are chosen in such a way that it is possible to establish capacitive sensors for measuring film thickness. Experimental results on the diethyl ether condensate film thickness measurements depending on the heat power value, entering to the HP’s evaporator are shown in

The present paper represents the numerical modeling of the vortex flows inside a profiled vapour channel of short HP’s, finite volume and pressure based CFD 10.0 code, for predicting compressible vapour flow with friction, heat transfer and converging-diverging nozzle area change. Navier-Stokes and heat-transfer equations with measured boundary conditions were solved, i.e. using fixed temperature values of heat source and heat outlet. A flat resistive heater is used as a located at the HP’s bottom heat source, and a original design flow calorimeter was used as a heat drain [

The model was studied as a longitudinal section along the axes of the two injector channels, which helps to preserve all the specific features of whirling instability under the conditions of continuous circulation motion of the working fluid during liquid and vapour phases. CFD 10.0 code employs a finite volume formulation of mass, momentum and energy conservation equations in a network consisting of nodes and branches. The mesh discretizes the nozzle space into small tetrahedral solid elements second order of accuracy and high quality, with ten nodes and enough accurate in the problem modeling. The forming vortex flow in the HP’s upper part, cooled by the flow calorimeter, are shown on

The cross-section profile of the vapour toroidal vortex ring near the cooled by the flow calorimeter upper HP’s cover 1,

On the moving in a radial direction on the cover to the capillary-porous insert 4,

If the toroidal rotation direction of the vapour vortex due to the Coanda effect occurs from the periphery to the longitudinal axis of the HP’s vapour channel as shown in _{f} in the near-axis area on the HP’s upper cover will have inflated values.

With increasing the heat power value, the direction of toroidal rotation of the vapour vortex acquires the opposite direction, from the longitudinal axis to the periphery of the vapour channel, as it shown in _{f} of the film diethyl ether on the cover will accelerate and its thickness δ_{f} will have to decrease.

_{B} = 15 K, and the temperature difference between the evaporator and the condensation surface T_{ev} − T_{cond} = 25 K.

Giving HP’s vapour channel profiled Laval-liked nozzle form at a large heat power value and the boil beginning in the capillary-porous evaporator and the formation of large vapour amount leads to the toroidal vortices formation near the cooled condensation surface.

The vortices occurrence is a decisive factor that changes the longitudinal components of vapour velocity and pressure distribution in the Laval-liked nozzle vapour channel, and explains the increase in the heat transfer coefficient K_{HP}, in comparison with standard HP’s with a cylindrical vapour channel and a simple vapour return currents in them at equal overall dimensions.

The calculated confirmation of the heat transfer coefficient K_{HP} increase for the HP’s with a Laval nozzle-liked vapour channel is shown in

upper part of the

In the center of the toroidal vortex ring there is a rarefaction (−δP_{v}), the value of which is about |−δP_{v}| ~ 2 × 10^{3} Pa, or 0.5% of the condensation vapour pressure 4.1 × 10^{5} Pa in the cooled part of the HP’s vapour channel.

The lower part of the

The counter current (negative value) velocity on the toroidal vortex outer surface forms a negative dynamic pressure about―2 × 10^{3} Pa, and the total pressure drop on the vortex ring near the cooled HP’s condensation surface reaches 4 × 10^{3} Pa, or about 1% of the condensing vapour pressure. The vapour pressure gradient inside the toroidal vortex with the vapour channel radius 1×10^{−2} m is no less than 4 × 10^{5} Pa/m. This value is more than enough to increase the condensation intensity of the diethyl ether along the vortex constant pressure gradient line (surface).

Let’s consider the interaction of a vapour jet with a HP’s flat cover. The origin of the coordinate system is placed at the junction central axis of the vapour channel with a HP’s flat cover, the z coordinate is measured in the opposite direction to the vapour jet direction.

The gas dynamics equations in the Cartesian coordinate system represent a description of a unsteady viscous compressible vapour flow interaction using the following equations:

∂ Q ∂ t + ∂ F ∂ x + ∂ G ∂ y + ∂ H ∂ z = 0. (1)

Equation (1) is supplemented by the ideal gas equation of state:

P v p = ( γ − 1 ) ρ v p [ e − ( u 2 + v 2 + w 2 ) 2 ] . (2)

Vector of conservative variables Q and vectors of vapour flows F, G, H have the following form:

Q = ( ρ v p ρ v p u ρ v p v ρ v p w ρ v p e ) ; F = ( ρ v p u ρ v p u u + p − τ x x ρ v p u v − τ x y ρ v p u w − τ x z ( ρ v p e + p ) u − u τ x x − v τ x y − w τ x z + q x ) . (3)

G = ( ρ v p v ρ v p v u − τ y x ρ v p v v + p − τ y y ρ v p v w − τ y z ( ρ v p e + p ) v − u τ y x − v τ y y − w τ y z + q y ) ; H = ( ρ v p w ρ v p w u − τ z x ρ v p w v − τ z y ρ v p w w + p − τ z z ( ρ v p e + p ) w − u τ z x − v τ z y − w τ z z + q z ) . (4)

The components of the viscous stress tensor and the components of the heat flux vector are written in a standard way:

τ i j = μ e f f ( ∂ i ∂ x j + ∂ v j ∂ x i − 2 3 ∂ v k ∂ x k δ i j ) ; q i = − λ e f f ∂ T ∂ x i . (5)

In the above equations, t―time, s; ρ_{vp}―the moist vapour density, kg/m^{3}; u, v, w, m/s, components of the vapour flow velocity in the coordinate directions x, y, z in the HP’s channel; p―pressure, Pa; e―total energy per unit mass of moist vapour, kJ/kg; T_{vp}―vapour temperature, K; γ―the specific heats ratio; λ_{eff}―effective thermal conductivity, W/m・K; μ_{eff}―effective viscosity is calculated in the standard way as the sum of molecular μ and turbulent μ_{t} viscosity:

d M d t = J v p . (6)

where M ˙ ―the increase in the mass of the condensed film in a unit of time, kg/s; J_{vp}―the condensing vapour flow, kg/s.

In the vapour flow braking area on the HP’s upper cover near the central axis (r = 0), the speed of the vapour jets significantly reduced, in the turning and vortex formation section are strongly curved. The depth (thickness) parameter of vortex is b.

In the near-wall section (vortex area formation), the shear stresses τ_{s} become dominant compared to normal ones τ_{n}. Spreading in the radial (and tangential) directions vapour jet loses momentum (amount of motion), the thickness of the boundary layer and the condensate film thickness δ_{f} increase and the static pressure increases also. The calculation results show the occurrence of a large vortex formation formed in the mixing layer of the vapour jet in confined conditions of the vapour channel with HP’s flat cover surface. The vapour pressure pulsations Reynolds number [_{vp} ~ (1.5 - 1.65) × 10^{5}, the Prandtl number Pr = 0.77.

The structure of the vortex formation on the HP’ upper cover have a defining effect on the heat exchange in this region of the vapour jets interaction with a flat end cap (flat wall). The heat transfer characteristics are non-uniform even at small values of the Reynolds number Re_{vp} ~ 5 × 10^{2}. The Nusselt number Nu has a characteristic maximum at the breaking point (r = 0), and the minimum value in the upper cover area near the walls where the inverse (counter) vortex flow is formed [

The justification for the heat exchange characteristics change in the vapour jets interaction point (r = 0) with the HP’s flat top cover and the occurrence of a Nusselt number local maximum are explained by the formation of a large-scale toroidal vortex [

To describe the flows arising from the interaction of a vapour jets with HP’s flat top cover, the Navier-Stokes equations (averaged over the Reynolds number) are usually used, in solving which the shortcomings of various turbulence models are apparent manifested. In the event of vapour toroidal vortices, an additional source of turbulence energy generation of the condensed vapour occurs due to the large curvature of the current lines near the HP’s condensation surface.

To assess the effect of the moving vapour jets curvature in the toroidal vapour vortices occurrence and a more complete assessment of the turbulent energy generation, the so-called Kato-Loner correction [

_{v}, with a temperature T_{v} and axial component velocity u_{v} and tangential component velocity v_{v}. We consider the circumferential (tangential) velocity component in the Laval-liked nozzle vapour channel is not equal to be zero. Static vapour pressure at the inlet of toroidal vapour ring vortex is equal to P_{v}, the cylindrical coordinate system is used.

As a main model variant, to simplify the analysis we consider the stationary state flow of a diethyl ether film condensate between the stationary state vapour toroidal ring vortex and the flat surface of the HP’s top cover. The vapour toroidal vortex ring and the condensate film move in the radial and tangential directions. The task is to determine the thickness δ_{f} of the diethyl ether condensate film under the toroidal ring vortex.

The proposed mathematical model uses the principle of the Karman and Pohlhausen [_{f} is considered. The peculiarity of the problem under consideration is that the longitudinal and tangential components of the velocity in vapour channel are commensurate.

Therefore, the momentum theorem is applied twice―for radial and circumferential directions of the vapour flow:

d d r ( ρ v p ( u v 2 − u f 2 ) r δ r * * ) + ρ v p r u f d ( u v − u f ) d r δ r * − ρ v p ( v v − v f ) 2 δ φ * * − 2 ρ v p v f ( v − v f ) δ φ * = r τ r s (7)

where u_{f}, v_{f} ―radial and tangential condensate film flow velocity under the vapour vortex, m/s; u_{v}, v_{v}―radial and tangential vapour velocity in the vapour vortex, m/s; δ_{f} ―the hydrodynamic thickness of the diethyl ether condensate film under the vapour vortex ring, m, in general depends on the vapour flow rate; δ r * * ―the condensate film thickness of flow momentum loss due to the friction in the radial direction, m; δ r * ―the condensate film displacement thickness, taking into account the displacement of the current lines in the condensate film due to the viscosity of working fluid, m; δ φ * * ―the condensate film thickness of flow momentum loss due to the friction in the tangential direction, m; τ_{rs} = μ(du_{v}/dz)z = δ_{f}―radial shear friction stress on the outer surface of condensate film under the vapour vortex ring, Pa.

Momentum conservation equation for tangential (circumferential) direction of the condensate film flow look like this:

d d r ( ρ v p u r 2 ( v − v f ) δ r φ * * ) = r 2 τ φ s . (8)

where τ_{φ}_{s} = μ(dv_{v}/dz)z = δ_{f}―tangential shear friction stress on the outer surface of the condensate film under the vortex ring, Pa; δ r φ * * ―the condensate film thickness of flow momentum loss due to the friction in the radial and tangential directions, m.

The thicknesses of the liquid condensate film δ^{*} and δ^{**} are determined in the standard way [

δ r * = ∫ 0 δ ( 1 − ρ l u f ρ v u v ) d z ; δ r * * = ∫ 0 δ ρ l u f ρ v u v ( 1 − ρ l u f ρ v u v ) d z .

δ φ * = ∫ 0 δ ( 1 − ρ l v f ρ v v v ) d z ; δ φ * * = ∫ 0 δ ρ l v f 2 ρ v v v 2 d z ; δ r φ * * = ∫ 0 δ ρ l u f ρ v u v ( 1 − ρ l v f ρ v v v ) d z . (9)

where ρ_{l}―the density of liquid condensate, kg/m^{3}.

Equations (7) and (8) contains more than two unknowns, therefore additional physically based linking equations are needed. The velocity vector deviates from the geometric radial direction, and the degree of deviation in the flux core and near the surface of the film are different.

Tangential (circumferential) vapour velocity v in the flux core in the in the shaped vapour channel we can define in the simplest way:

v v r v m = v r v p m

where m is a numerical parameter of unity order; r_{v}―the inner radius of the toroidal vortex, m.

The relationship between the pressure, density and temperature of the condensing vapour in the first approximation can be given by the ideal gas equation of state in the following form:

ρ v p = P R T ; T = T v ( P P v ) ( k − 1 ) / k ; ρ v p ρ v = ( P P v ) 1 / k . (10)

where R = 8.3144 J/mol・K―universal gas constant; k = 1.31―value of the adiabatic index of diethyl ether vapour; T_{v}―the temperature inside the toroidal vapour vortex, K; ρ_{v}―the vapour density inside the toroidal vortex, kg/m^{3}; P_{v}―the vapour pressure inside the toroidal vortex, Pa.

After integrating the differential Equation (7) for the vapour core flow in the Laval nozzle shaped vapour channel along the radius r and taking into account the equation of state (10), we obtain an equation to estimate the pressure distribution over the radius r in the core flow:

P P v = [ 1 + ( k − 1 ) k 1 R T v ( u 2 − u v 2 2 + v 2 − v v 2 2 m ) ] k / ( k − 1 ) . (11)

From the continuity equation of the vapour flow in the integral form and subject to the vortex size we obtain the next correlation:

ρ v u v r v 2 = ρ v p u r ( b − δ r * ) .

where b―the thickness of the vapour toroidal ring vortex, m.

From this relation, taking into account the density ρ_{vp} from (7) and (8), an equation can be obtained to estimate the radial velocity in the core of the vapor flow at the inlet to the vapour toroidal ring vortex:

u = u v r v 2 r ( b − δ r * ) [ 1 + ( k − 1 ) k 1 R T v ( u 2 − u v 2 2 + v 2 − v v 2 2 m ) ] k / ( k − 1 ) . (12)

In general, the vapour flow in the Laval nozzle shaped vapour channel is swirling flow, which is due to both the uneven jet character of operation of capillary-porous evaporator with injection channels and the variable area of the HP’s passage section. For swirled vapour flows, the spin parameter is the determining in the HP’s vapour channel parameter, characterizing the ratio between the circumferential (tangential) and longitudinal (axial) velocity components inside the channel.

To determine the local twist of the passing into the liquid phase (condensing) vapour stream through the thickness of the diethyl ether film condensate (boundary layer), the parameter of the local spin angle is used [

tg φ = v v u v = ε ( r , z ) . (13)

where tgφ―parameter of local spin angle.

Within the variable thickness of the of liquid condensate film under vapour toroidal ring vortex parameter of local twist tgφ depends on the radius on the flat surface of the HP’s condensation top cover and the linear parameter of the film thickness z. At the vapour stream core in the diffuser part of the vapour channel tgφ_{0} parameter characterizes the degree of deviation of the vapour flow from the radial direction:

tg φ 0 = v u = ε ( r ) 0 . (14)

where tgφ_{0}―parameter, defining the degree of deviation in the diffuser part of the HP’s vapour channel from the radial direction.

The liquid condensate film, formed under the vapour toroidal ring vortex, due to surface interfacial friction is fond of it and moves in the radial and tangential directions on the HP’s smooth condensation flat (wall) surface, moreover, the wall shear stress τ_{w} ~10^{2} Pa [

The parameter tgφ_{max} represents the limiting (surface) tangent of the twist angle of the formed liquid condensate film and is the ratio of surface shear stresses of friction in tangential (circumferential) and radial directions on the film surface of the film of condensate of variable thickness and can be recorded as follows:

tg φ max = ε ( r ) max = τ φ s τ r s . (15)

For evaluation the circumferential velocity component in the turbulent diethyl ether condensate film under the vapor toroidal ring vortex in accordance with the Nikuradze’s experiments results [

v f v v = ( z δ f ) n . (16)

where the exponent of the parabolas’ n depends on the Reynolds number Re_{vp} of the vapour flow and changes in the interval 1/6 … 1/10 (0.17 - 0.1).

The radial component of the liquid condensate film velocity u_{f} under the vapor toroidal vortex ring can be defined as follows:

u f u v = tg φ χ ( v v − v f ) v v = tg φ χ ⋅ ( z δ f ) n . (17)

where tgφ―twist parameter of the surface layer of the liquid condensate film, formed on the HP’s condensation surface under the vapour toroidal vortex; χ― friction parameter of the toroidal vapour vortex ring.

The twist parameter tgφ decreases within the thickness of the condensate z, varies according to the quadratic parabola law [

tg φ = ε − ( ε − χ ) ( 2 − z δ f ) z δ f . (18)

where ε―the tangent of the limiting (surface) angle of twist of the vapour flow is the ratio of surface shear stresses of friction in the circumferential (tangential) and radial directions on the condensate film surface.

In the equations for the momentum conservation of the vapour moving (7) and (8), the projection of the radial and tangential friction stress τ_{s} on the condensate film surface must be determined. To do this, we use the Prandtl hypothesis [_{s} (with components τ_{φs}, τ_{rs}) on the film surface as follows:

β = τ s ρ l δ f ν l = 11.65 ; τ s = τ φ s 2 + τ r s 2 . (19)

where β―the criterion of diethyl ether film condensate flow stability under the vapour toroidal ring vortex; ν_{l}―kinematic coefficient of diethyl ether viscosity, m^{2}/s.

The limit value of the working fluid film thickness δ_{f} under the toroidal vortex ring can be determined from the condition of velocities equality of the outside film surface and the toroidal vapour ring vortex. Given that the tangential constituent of the condensate film flow velocity provided the sticking on the inner surface of the top cover (v = 0) and can be represented as:

u f = τ r s μ l z . (20)

therefore, the calculation of the diethyl ether film flow rate can be carried out by the formula:

τ r s μ l δ f = u f + ( u v − u f ) ( δ f δ ) n . (21)

where μ_{l}―dynamic coefficient of diethyl ether viscosity, Pa・s; δ―condensate film thickness at immovable (stationary) vapour, m.

Measurements of the film flow rate were carried out earlier [

δ f = ( u v − u f δ n ⋅ μ l τ r s ) 1 / ( 1 − n ) . (22)

After replacing the limit δ_{f} in Equation (22), we can obtain an expression for the radial component of the vapour―condensate film shear stress τ_{rs} as the following:

τ r s = β 2 ( n − 1 ) n + 1 ρ v p u v 2 R e δ − 2 n n + 1 ( 1 + ε 2 ) − 0.5 ; R e δ = u v δ f ν . (23)

where Re_{δ}―is the Reynolds number, calculated on the radial velocity u_{f} and thickness δ_{f} of the diethyl ether condensate film flow.

The friction coefficient at the interfacial surface is determined in a standard way using the Blasius ratio:

τ r s = c f ρ v p u v 2 2 ; с f = 0.664 ( u v z / υ l ) 1 / 2 ~ 10 − 2 ÷ 10 − 1 . (24)

When the film condensation occurs on the oriented in the horizontal direction HP’s top cover, the projection of gravity on the film flow direction is zero, and the provision of shear motion of the film occurs only due to the dynamic effects of vapour flow moving along the film outer surface. Motion in a thin film is considered to be inertialess.

Measurement average values results of the of the diethyl ether condensate film thickness on the cooled surface of the HP’s upper cover depending on the thermal power on the HP’s evaporator are shown in ^{−3} mm [

A sharply non-linear and decreasing values of the film thickness of the diethyl ether on the cooled smooth surface of top cover in dependence of the HP’s evaporator’s overheating compared to the boiling point of working fluid δT = T_{ev} − T_{B}. At low thermal power loads and the vapour vortex flow from the walls to the longitudinal axis of the vapour channel, condensate film has a large thickness. As the thermal power on the HP’s evaporator increases, the direction of the vortex flow changes to the opposite, the direction of vapour and condensate film flow become wake (in one direction), and the film thickness begins to decrease sharply.

The standard mean-square time average error of the film thickness measurement is ~ 1 × 10^{−3} mm, but begins to increase substantially starting from the evaporator superheat value δT = T - T_{B} ~ (11 - 12) К, at which begins the boiling and a two-phase moist vapour form in the capillary porous evaporator. The condensing of the two-phase moist vapour leads to a noticeable increase in the experimental spread in the capacitive measurements of the diethyl ether film thickness on the cooled condensation surface.

The Reynolds number of the moving film of diethyl ether with thickness δ = 10^{−2} mm, ν_{l} = 0.032 × 10^{−4} m^{2}/s, Prandtl number Pr = 0.77, at equal velocities of the vortex (30 m/s) and the film is equal to:

R e δ = u v δ f ν = 30 m / s × 10 − 5 m 0.032 × 10 − 4 m 2 / s = 0.94 × 10 2 . (25)

τ r s = β 2 ( n − 1 ) n + 1 ρ v u v 2 R e δ − 2 n n + 1 ( 1 + ε 2 ) − 0.5 = 2.59 × 10 2 Pa ; u f ~ τ r s ρ l ~ 0.1 - 10 m / s (26)

The calculated values of the film thickness are in the range of 10^{−2} - 10^{−4} mm when the value of friction coefficient is c_{f} = 10^{−2} - 10^{−3}, condensing vapor Reynolds number Re_{vp} = 10^{3} - 10^{4} and the vortex velocity up to 30 m/s. That is close to the experimentally obtained values of the film thickness of the diethyl ether condensate [

1) For the first time by the calculation method, it was noted the possibility of changing the rotation direction of toroidal vortices of condensed vapour near the cooled HP’s condensation surface with Laval-liked nozzle vapour channel. At low values of the thermal power entering to the evaporator, the toroidal (axial) rotation of the vapour vortex due to the Coanda effect occurs from the periphery to the longitudinal axis of the vapour channel. A rotating toroidal vapour vortex acts on the moving condensate film and changes its thickness. With low thermal power and vortex rotation from the periphery to the center of the vapour channel, the condensate film movement due to the vortex friction about the surface is braked and slowing down, whereby the film thickness in the radial area of the upper cover effectively increases, has inflated values.

2) With increasing thermal power, the toroidal rotation direction of the vapour vortex changes to the opposite, from the z axis to the vapour channel walls and due to the friction on the outer surface of moving to the capillary-porous insert condensate film, the film radial flow velocity increases and its thickness decreases. The toroidal vortex additionally ‘blows off’ condensate film to the channel periphery along the HP’s upper flat cover surface.

3) The experimental results of the diethyl ether condensate film thickness depending on the HP’s thermal power (the evaporator overheat relative to the boiling point of the diethyl ether at atmospheric pressure) show a sharp decrease in the film thickness with an increase in thermal power. Such a sharp dependence of the film thickness on the thermal power indirectly confirms the change in the toroidal rotation direction of the HP’s vapour vortex. Increased film thickness (and HP’s high thermal resistance) at low loads and a sharp decrease in film thickness (and significantly reduced HP’s thermal resistance) with increasing thermal power can be associated with a change in the direction of the toroidal rotation of vapour vortex.

The author declares no conflicts of interest regarding the publication of this paper.

Seryakov, A.V. (2019) Numerical Modeling of the Vapour Vortex Formation in the Short Heat Pipes. Journal of High Energy Physics, Gravitation and Cosmology, 5, 218-234. https://doi.org/10.4236/jhepgc.2019.51013