THE STUDY OF THE FLOW IN THE VAPOUR CHANNEL OF SHORT LINEAR HEAT PIPES
- Autores: Seryakov A.V.1
-
Afiliações:
- LLC Research and Development Company “Rudetransservice”
- Edição: Volume 18, Nº 3 (2017)
- Páginas: 592-603
- Seção: Articles
- URL: https://journals.eco-vector.com/2712-8970/article/view/503369
- ID: 503369
Citar
Texto integral
Resumo
E-mail: seryakovav@yandex.ru The results of studies of flow of moist vapour in Laval-like vapour channels of short linear heat pipes (HPs) are pre- sented. The increase in heat transfer coefficient of short linear HPs, intended for creation of the cooling systems of heat-stressed designs of spacecraft, is carried out by making the HPs vapour channel forms similar to the shape of the Laval-like nozzle. Comparison of the heat transfer coefficients of short HPs with the standard cylindrical vapour channel and the channel, made in the form of the Laval-like nozzle with the equality of all dimensions, flat evaporator and the same amount of the working fluid, shows that the HPs with the vapour channel in the form of the Laval-like nozzle exceeds the heat transfer characteristics of the standard HPs with a cylindrical vapour channel under high ther- mal loads. The study of the flow and condensation in such shaped vapour channels of the short HPs at high thermal loads gives an opportunity to analyze in detail the advantages of using such HPs and make the conclusion about the necessity of wide introduction of such HPs for the cooling systems of the spacecraft. Capacitive sensors were addition- ally installed in cooled top covers of the HPs, and electromagnetic pulses with a frequency of 100 kHz were supplied to them from the external generator. At heating the HPs evaporator, starting from a certain thermal power threshold value, electromagnetic pulses became modulated. It is related to the formations of the boiling process in the capillary- porous evaporator and large amount of vapour over it and its discontinuous distribution. It was discovered that the frequencies of the pulsations are more and they occur at lower values of heat load at the evaporator in the HPs with a vapour channel of a Laval-like nozzle, compared to the HPs with the standard cylindrical vapour channel with equal overall dimensions.
Palavras-chave
Texto integral
Introduction. The increase of heat transfer coefficient of short linear HPs is a very important task when design- ing a spacecraft with a rigid regulation of take-off weight. The application of the linear HPs is justified in case of structural impossibility of placing a loop (contoured) HPs for it can improve long-term reliability of the cooling sys- tems, which are distributed inlet and outlet closed lines of the loop (contour) HPs. Pulsating (or oscillating) loop HPs - siphons, was invented in 1990 by Akachi [1]. When the temperature difference between evaporator and condenser exceeds a certain threshold, Tong et al. [2], R. Naik, et al. [3], the vapour bubbles and condensate liquid plugs begin to oscillate spontaneously back and forth. The amplitude of oscillations is quite strong and the liquid plugs penetrate into both condenser and evaporator. The heat is thus transferred not only by the latent heat transfer liquid- vapour-liquid like in other types of heat pipes, but also by conductive heat transfer of the hot oscillating liquid plugs to cold walls in the condensation zone of loop HPs. This phenomenon is the reason of high heat transfer coefficient of loop HPs in comparison with other types of HPs. The increase of the heat transfer coefficient in short linear open HPs with capillary porous insert is achieved using inner vapour channel of Laval-liked nozzle form. There is an availability to use a jet vapour nozzle, similar to the Laval nozzle, surrounded by a capillary-porous insert along the full length of the HPs, in case if the heat flow is axially directed that is typical for short HPs. The shape of this vapour Laval-like channel, is defined by the inner capillary porous insert. The implementation of the nozzle increases flow velocity [4], pulse rate of two-phase vapour flow and heat-transfer coefficient in the HPs. Pulsation of the flow in short linear HPs with the capillary porous insert represents more complex pheno- mena, associated with boiling of the working fluid in the capillary porous evaporator. Heat-transfer coefficient and operational effectiveness of short linear HPs with capillary-porous inserts and Laval-like vapour channel are defined by closed circulation motion of the working fluid, which undergoes liquid- vapour transition and heat absorption in the evaporation region of the capillary porous insert, vapour phase trans- fer through a convergent-divergent vapour channel, vapour-liquid transition with heat liberation in the con- densation region, and liquid return to the evaporation region through capillary-porous inserts of the HP. Improvement of the heat-transfer coefficient of heat pipes is a complex task, concerning the device’s detailed mechanical design and gas dynamics. Flowing vapour, with microdrops of condensate, appears as a nonlinear flow, with attendant internal processes, such as interpha- sic heat and mass exchange, and power dissipation, Gupta et all [5]. Static regain of supersaturated vapour flow in the diffuser condensation region of the HPs with turbula- tor, causes flow stagnation, strong vorticity production near the condensation surface and reverse vapour flow. Thermal load-related processes of vorticity formation and density and pressure fluctuation, proceeding in the convergent-divergent vapour channel of HPs, are of the utmost interest. For a field research of these processes, stainless steel heat pipes have been manufactured with diethyl ether as the working fluid [6-9]. Materials and methods. The diethyl ether C4H10O is selected as the main working fluid, which has the boiling temperature TB = 308.55 K (35.4 °C) under the atmos- pheric pressure, freezing temperature TF = 156.95 K (-116.2 °C) and critical parameters TC = 466.55 K (193.4 °C), PC = 3.61 MPa. Successful use of diethyl ether as the working fluid for Wilson chambers, on long liquid phase at 413.15 K (140 °C), shows its heat resistance and allows using it as the HP working fluid. Fig. 1. Layout of HP equipped with capacitance sensors: 1 - flat upper lid with a smooth surface; 2 - cylinder body of HP; 3 - cone-shaped turbulence agitator; 4 - capillary-porous insert defining the vapour channel’s form; 5 - flat bottom lid; 6 - injector channels; 7 - capillary-porous evaporator; 8, 9 - capacitance sensors, one of which is intended for a condensate film thickness measurement, while another one has a sensing element of CT3-19 thermistor mounted on its electrodes to measure the film temperature [11] Рис. 1. Схема тепловых труб: 1 - верхняя крышка; 2 - цилиндрический корпус ТТ; 3 - конический турбулизатор; 4 - капиллярно-пористая вставка; 5 - нижняя крышка; 6 - инжекторные каналы; 7 - капиллярно- пористый испаритель; 8, 9 - емкостные датчики конденса- ции, измеряющие толщину слоя жидкого конденсата. Изме- рительная поверхность емкостных датчиков совпадает с внутренней поверхностью плоской верхней крышки 593 The developed sensor enables us to perform measure- ments of local characteristics of the film flows primarily, the film thickness and temperature, without making any major disturbances in the flow. To measure the fluid con- densate of film thickness, short HPs were used with vapour channel made in the form of a Laval-like nozzle [10-12]. The HPs are described in details below. The capillary-porous evaporator 7 is constructed of 43 circular layers of thin stainless steel net (fig. 1). Each layer is 0.07 mm thick with the mesh size of 0.04 mm and the aggregate thickness of 3 mm; all layers are spot- welded to the flat bottom lid 5. Outer diameters of the circular layers gradually decrease in such a manner that the lateral side of the assembled evaporator forms a pointed cone with a 45° total point angle that aligns with To ensure accurate measuring of thermal power and heat removal augmentation in the HPs, jet flow of input water is swirled, values of flow velocity and vorticity due to air bubbles are recorded. The HPs evaporators, also equipped with thermocouples, is heated using a resistance heater, and the temperature is maintained at δТ °C higher than the diethyl ether boiling temperature of 308.55 K (35.4 °C) under atmospheric pressure. The heater temperature is stabilized and HPs evaporator’s overheat value is set in the range of δТ = 0-20 K, herewith thermal power of single HP does not exceed 120 W. Measurement were performed to disclose variability of the heat transfer coefficients in Laval-like HPs and with standard cylindrical vapor channel. The heat transfer coefficients in HPs, Faghri [13] is defined by the formula: the angle of the end part of the capillary-porous insert 4. The evaporator as a whole is further equipped with 4 vapour KHP = E F (z) (T - T ) , (1) injection channels 6 with increasing diameter of 1 mm, 1.1 mm, 1.2 mm and 1.3 mm diameter that are inclined to the longitudinal axis of HP at a specific angle and create a swirl- ing jet stream of vapour. The capillary-porous insert 4 is assembled on a dedi- cated profiled frame made of layers of thin stainless steel net same as the evaporator. Each layer is 0.07 mm thick with mesh size of 0.04 mm, the total thickness of the insert is 1.5 mm at the edges and 7.5 mm at the critical diameter of the inner nozzle. All layers of the insert are radially stitched together along the forming insert with a ev cond where KHP is heat transfer coefficient of the HPs, W/K; E is thermal power, supplied in the evaporator of the HPs, W; F(z) is surface area of the evaporator inside the vapour channel of the HPs, m2; Tev is temperature of the evapora- tor, K; Tcond is temperature of the condensation surface, K. Temperature of evaporator surface and condensation surface of HPs were measured using differential thermo- couples copper-constantan, heat power E, W; transferred by condensations zones of HPs to the calorimeter, was defined by the formula: thin wire 0.05 mm in diameter and stitch spacing of 5-7 mm. æ E = G × C - CK ö V × T (t) - T (t) , (2) r Eight lines of the wire stitches that radially fix the net layers together are tilted by 45° towards each other and çç H2O è ÷÷ [ K K 0 ] H2O 0 ø form a rigid structure of the capillary-porous insert with an inner vapour channel made in the form of a nozzle, similar to the Laval nozzle. At the edge of insert 4 that is adjacent to the evaporator 7, an end portion is placed; its total angle is somewhat wider than that of the diffuser portion of the vapour channel and equals 45°. The longi- tudinal length of the end portion is 3 mm and matches the thickness of the evaporator 7 exactly. The 45° angle of the end portion of insert 4 is formed by the gradually decreas- ing net layers’ lengths, from the outer layer to the inner one, and matches the angle of the evaporator 7 exactly. The multilayered evaporator with a tapering lateral side is tightly inserted into the rigid capillary-porous insert 4 at the full length of the end portion and is spot-welded to it in 8 spots along the perimeter of the insert. After being cooled down at the boiling temperature of fluid nitrogen 77 K (-196 °С), the assembled capillary-porous insert together with the evaporator on the flat bottom lid 5 is tightly inserted in the stainless steel cylindrical shell 2 of the HP. Porosity of the insert and evaporator is 72 % and to- gether they form one hydraulic system designed to deliver the service fluid to the evaporator when the HP is operat- ing. The HP length is 100 mm, its diameter is 20 mm, the max. diameter of the vapour nozzle in the convergent and divergent regions is 16 mm, the critical nozzle diameter is 4 mm, the length of the nozzle convergent region is 13 mm, the total angle of the convergent region is 41°, the length of the nozzle divergent region is 81 mm, the total angle of the divergent region is 8.5°, and the length of the cylindrical region in the nozzle throat section is 1 mm. Measurement of the heat transfer coefficient. Con- densation zones of the HPs are provided with insulated thermocouples and set into the vortex continuous-flow calorimeter, shown in fig. 2, with stabilized water flow. where G is water flow in a vortex flow calorimeter, kg/s; CH2O is the specific heat of water, J/kg·K; CK - the heat capacity of the calorimeter, J/K; ρH2O is the density of water in the calorimeter, kg /m3; V0 is the internal volume of the calorimeter, m3; TK0(τ), TK(τ) is the water tempera- ture at the inlet and outlet of the calorimeter, K. Input temperature of water at the calorimeter, tem- perature differences (heat value) of the flowing water at entrance and exit of the calorimeter were also measured using differential thermocouples copper-constantan, com- parator Р3003 and digital voltmeter V7-34A. Water flow rate was measured using ultrasonic flowmeter. Total measuring inaccuracy of the thermal power E of the HPs does not exceed 1.7 %. Comparison of heat transfer coefficients in Laval-like HPs and HPs with standard cylindrical vapour channel with the equality of all dimensions showed following results, presented in fig. 3. The increase of the heat transfer coefficient in the HP with a Laval-like vapour channel appears to be related to the emergence of a vortex ring near the condensation surface of the HP. Numerical model details. Numerical simulations of the vortex flows inside a vapour channel of the Laval-like HPs have been performed in finite element modeling in CFD 10.0 code Fluent 6.3.26 2D under double preci- sion axis-symmetric conditions. Navier-Stokes equations with measured boundary conditions were solved, i. e. using fixed temperature values of heat source and heat outlet. In the construction of the design model we used about 457233 finite elements with increased meshing at injection capillary channels sections, nozzle throat section and turbulence element. Fig. 4 shows a diagram of the vortex obtained by numerical simulation of CFD 10.0 code. Fig. 2. Vortical calorimeter: 1 - vortical continuous-flow calorimeter; 2 - heat pipes bolting flange; 3 - glass cover; 4 - cover fastening; 5 - heat pipes; 6 - resistance heater; 7 - outlet stub tube for water flow; 8 - inlet stub tube for water flow; 9 - silicone sealant of the sensing wire; 10 - capacitive sensors for measuring the thickness of the condensed layer of the working fluid; 11 - the measuring and reference generators of the capacitive transducer; 12 - external digital generator; 13 - power assist element; 14 - digital oscilloscope; 15 - computer; 16 - commutation switch; 17 - digital voltmeter; 18 - container for constant water head; 19 - source of air bubble; 20 - water flow meter; 21 - vacuum-jacketed zero temperature container Рис. 2. Вихревой калориметр: 1 - вихревой проточный калориметр; 2 - фланец крепления тепловых труб; 3 - стеклянная крышка; 4 - крепление крышки; 5 - тепловые трубы; 6 - резистивный нагреватель; 7 - выходной штуцер для воды; 8 - входной штуцер для воды; 9 - уплотнение измерительных прово- дов; 10 - конденсаторные датчики измерения толщины слоя конденсата; 11 - высокочастотные гене- раторы; 12 - генератор импульсов Г5-56; 13 - усилитель; 14 - осциллограф АКИП-4116/2; 15 - ком- пьютер; 16 - коммутатор; 17 - вольтметр В7-34А; 18 - сосуд постоянного напора воды; 19 - генератор пузырьков воздуха; 20 - расходомер воды; 21 - сосуд Дьюара Fig. 3. Comparison of heat transfer coefficients in Laval-like HP and HP with standard cylindrical vapour channel: 1 - Laval- like HPs; 2 - HPs with standard cylindrical vapour channel 595 Рис. 3. Сравнение коэффициентов теплопередачи ТТ с паро- вым каналом в виде сопла, близкого к соплу Лаваля, и ТТ с цилиндрическим паровым каналом при равных габаритных размерах: 1 - ТТ с паровым каналом в виде сопла; 2 - ТТ со стандартным паровым каналом Fig. 4. The occurrence of a vortex ring near the condensation surface inside of Laval-like HPs. The numerical analysis of the flow in the Laval-like HPs condensation zone shows that the vortex structure has a spatial nature, at the same time the flow asymmetry becomes apparent being determined by non-linear friction against the underlying surface and the two-dimensional compressibility of the wet condensing vapour as well. Overheating of the evaporator in reference to the boiling temperature of diethyl ether at atmos- pheric pressure is 15 K, the temperature difference between the evaporator and the condensation surface equals 25 K Рис. 4. Возникновение вихревого кольцевого течения вблизи поверхности конденсации вокруг турбулизирующего элемента внутри ТТ с паровым каналом в виде сопла, близкого к соплу Лаваля. Перегрев испарителя относительно температуры кипения диэтилового эфира 15 К, разница температур между испарителем и поверхностью конденсации ТТ в кало- риметре 25 К Fig. 5 shows the results of numerical analysis of the distribution of the longitudinal component of the velocity of compressible moist vapour within the Laval-like HP. The formation of a wet condensing vapour vortex ring near the condensation surface inside of Laval-like HPs leads to quite interesting results. Vortex ring is a highly gradient zone of the condensable moist vapour velocity with the opposite directions inside and outside the vortex ring. The numerical analysis of the velocity distribution along the centerline of the vapour channel in the Laval- like HPs including the vortex ring shows the occurrence of two positive peaks of the velocity, one of them is in the critical section of the nozzle and the other is near the con- densation surface. In a critical section of the nozzle the axial component of velocity reaches 85 m/s, and near the condensation surface the counter-flow reaches the veloc- ity of 33 m/s. The distribution of the axial component of velocity of the moist vapour shows the presence of a counter current due to the formation of a ring vortex within the vapour channel near the HPs condensation surface. This figure of the velocity distribution, shown in fig. 5, confirms the fact that the vortex ring is a zone of sharp velocity gradients, and consequently is a pressure gradient. In the central part the vortex ring has a notice- able positive dynamic pressure, and it has a negative dy- namic pressure in the peripheral part of the vortex ring. This means that lower static pressure occurs in the central part of the vortex ring and this leads to additional absorption of moist vapor in the condensation zone of Laval-like HPs. Apparently, this means that the static pressure in the vortex center in the Laval-like HP is less than the pressure in the axial part of the condensation zone of the HP with standard cylindrical vapor channel. The appearance of high static pressure in the peripheral part of the vortex ring, where there is a negative dynamic pressure leads to intensification of the condensation process near the top cap of the Laval-like HPs that leads, apparently, to an increased value of heat transfer coefficient in the Laval- like HPs. This explains the reason why the heat transfer coefficient of Laval-like HPs becomes greater than that of HPs with the cylindrical vapor channel at boiling point in evaporator HP and there is formation of large amounts of vapor and increase in the efficiency of the nozzle in comparison with a cylindrical channel. Measurement of pulse characteristics of HPs. The pulsation characteristics of HPs were measured and discussed previously [6; 7]. Starting from the specific overheating δTS = (Т - ТB) of HPs evaporators, diethyl ether, which is the working fluid in the capillary-porous evaporators, starts to boil and produce large amount of vapour. The vapor pressure increases to values of P*, at which the boiling point of diethyl ether is less than the evaporator temperature Tev, and the process of evapora- tion stops. The pressure pulsations with amplitude up to 3·10-4 Pa spread through a vapour channel and condense in the cooled part of HPs. The vapor pressure decreases, the rarefaction wave comes to the evaporator and the next cycle of vapour formation and its moving to the HPs con- densation zone begins. The frequency of pulsations in the Laval-like vapour channel is larger than the pulsations in the HPs with the standard cylindrical vapour channel. The obtained values of the pulsations frequencies are given in table. 597 Fig. 5. The distribution of the axial component of velocity of the moist vapour inside the vapour channel of Laval-like HPs along the HPs centerline. Upper pattern with a positive veloc- ity peak near the condensation surface of Laval-like HPs represents the velocity distribution along the longitudinal axis of vapour channel and inside the vortex ring. The figure on the bottom with a negative velocity value peak near the surface of condensation represents the velocity distribution along the longitudinal axis and outside of a vortex ring and shows the presence of a counter current due to the formation of a ring vortex near the HPs condensation surface. Overheating of the evaporator in reference to the boiling temperature of diethyl ether at atmospheric pressure is 15 K, the temperature differ- ence between the evaporator and the condensation surface equals 25 K Рис. 5. Распределение расчетных значений осевой состав- ляющей скорости влажного пара внутри выполненного в виде сопла, близкого к соплу Лаваля, парового канала ТТ. Перегрев испарителя относительно температуры кипения диэтилового эфира 15 К, разница температур между испа- рителем и поверхностью конденсации ТТ в калориметре 25 К Pulsation rate values Overheating of the evaporator, δT, K Pulsation frequency f1, Hz Pulsation frequency f2, Hz 9.05 386 ± 5 - 10.1 396 ± 5 - 11.03 426 ± 5 406 ± 5 12.15 450 ± 5 420 ± 5 13.0 456 ± 5 437 ± 5 14.07 474 ± 5 440 ± 5 15.03 474 ± 5 454 ± 5 16.0 478 ± 5 453 ± 5 17.1 490 ± 5 460 ± 5 18.06 491 ± 5 472 ± 5 19.02 495 ± 5 473 ± 5 20.12 502 ± 5 474 ± 5 Table shows the pulsation rate values (modulation rate), obtained in the Laval-like HP, f1 [Hz], and in the HP with standard cylindrical vapour channel, f2 [Hz], depend- ing on overheating δT, [K], of the evaporators. Both HPs have equal outer diameters 20 mm and equal sectional areas of capillary-porous inserts near the condensation region. As the overheating of the evaporator δT ~ 20 K increases, pulsation frequency in the Laval-like nozzle goes up to ~ 502 Hz, the derivative of the relation- ship between pulsation frequency and temperature is approximately 10.5 Hz/K. While investigating HPs with standard cylindrical vapour channel and equal outer diameter of 20 mm, length of 100 mm and thickness of the evaporator and capillary-porous insert of 3 mm, initial pulsed flows occur in event of overheating of the evapora- tor δT ~ 11 К, frequency f ~ 406 Hz. As the overheating of the evaporator δT ~ 20 K increases, pulsation frequency in the cylindrical vapour channel goes up to 474 Hz, the derivative of the relationship between pulsation frequency and temperature is approximately 7.5 Hz/K. The insensitivity zone of the capacity sensors in the cylindrical channel, defined by the initial convective nature of the vapour flow, is greater than in the vapour The figure presents test values of the vapour flow pulsa- tion in the vapour channel in the Laval-like form of the HPs. There are seven steps of calculation of the moist vapour flow velocity, from up to down. When boiling in the capillary-porous insert and the initial velocity of the flow of the two-phase vapour directly over the surface of the evaporator is 0.5-1 m/s in the throat section of the nozzle the axial flow velocity reaches 5 m/s and more. Obvious occurrence of pulsations, pulsation frequency rises, what is consistent with the experimental results [6; 7]. Validation of numerical scheme. Evident conversion from convection mode to convection-vortex mode and then to pulsation mode of the vapour flow inside of the HP is recorded at evaporator overheating value δT ~ 9 К. Maximum value of velocity of diethyl ether moist vapour flow, obtained using calculation method in the throat section of the vapour channel using the colour indication of the software CFD Design 10.0, reaches 100-110 m/s when pulsation occurs, fig. 5. This fact gives the opportu- nity to evaluate the Reynolds number Re of the vapour channel during pulsation, which is defined by the for- mula: rmixu D channel of Laval-like HPs. Measuring inaccuracy does not exceeds 3-5 Hz. Re = vp vp C . hmix (3) Numerical model details. Numerical simulations of We substitute moist vapour density (vapour and the vortex pulsation flows inside a vapour channel of the Laval-like HP have been performed in finite element drops) [14; 15] mix r np ~ 3 kg/m3, dynamic coefficient of modeling in CFD 10.0 2D Fluent 6.3.26 under double precision axis-symmetric conditions. Navier-Stokes equa- tions with measured boundary conditions were solved, using fixed temperature values of heat source and heat outlet. The model was studied as a longitudinal section along the axes of the two injector channels. It 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. In the construction of the design model about 457233 finite viscosity of the moist vapour ηmix ~ 8·10-6 Pa·s; maximum value of the vapour flow velocity (colour indication) near the throat section of the vapor channel uvp ~ (100-110) m/s; critical diameter of the vapour channel DC ~ 4·10-3m, and obtain the value Re ~ (1.5-1.65)·105, the Prandtl number Pr = 0.77. Duration of the pulsation period Δτ0 inside the diver- gent part of the HPS vapour channel can be estimated using the formula: Dl elements were used, with increased meshing at injection capillary channels sections, nozzle throat section and Dt0 ~ . u vp (4) turbulence element. The model size is a compromise between available computer resources and computational investigation error. Transition from a stationary convective flow regime in the vapour channel to a pulsatile flow regime is clearly visible. At high heat capacity, received by the HPs, and boiling in the grid evaporator, excessive vapour occurs in the convergent part of the nozzle, leading to pressure increase up to the value, at which the average temperature of layers of the flat grid evaporator becomes lower than the boiling temperature of the working fluid, and boiling in the evaporator stops. Excessive vapour spreads through the divergent part of the vapour channel to the cooled area of the HP and is partially condensed. Due to the conden- sation effect, the pressure in the vapour channel decreases Then substitute values of distance between pulsation crests (maximums) obtained in fig. 6 Δl ~ (2-3)·10-2 m, moist vapour flow velocity in divergent part of the vapour channel obtained using colour indication by the software CFD Design 10.0 uvp ~ (20-30) m/s, and you will obtain numerical value of duration of the pulsation Δτ0 ~ (0.75- -1.5)·10-3 s. Analytical evaluation. Evaporation cycle in the HPs exist at low heat load of the evaporator, up to 4-6 W/cm2, and is characterized by convective flow in the vapour channel. Evaporation cycle of short HPs, when evaporator heating power is constant and maximum value is limited to prevent the development of bubble boiling in the flat grid evaporator, is defined in the following way: DQ and the boiling process in the flat evaporator proceeds. Pressure increase period, vapour wave spreading in the E = Dt < EB , (5) condensation area of the HP and return expansion wave define the pulsation period in the vapour channel. The results of flow simulation of compressible super- saturated vapour environment inside a vapour channel. where ΔQ is heat energy, absorbed in the evaporator over a period of Δτ, J; Δτ is time, sec; EB is thermal power, wherein the process of bubble boiling begins in the grid evaporator, W. Fig. 6. The results of flow simulation of compressible supersaturated vapour environment inside a vapour channel. The figure presents test values of the vapour flow pulsation in the vapour channel in the Laval-like form of the HPs. There are seven steps of calculation of the moist vapour flow velocity, from up to down. When boiling in the capillary-porous insert and the initial velocity of the flow of the two-phase vapour directly over the surface of the evaporator is 0.5-1 m/s in the throat section of the nozzle the axial flow velocity reaches 5 m/s and more. Obvious occurrence of pulsations, pulsation frequency rises, what is consistent with the experimental results [6-7] Rate of evaporation of dry monomolecular vapour over the evaporator is defined by the following equation: the evaporator and mass flow of vapour in HPs, is calcu- lated from the equation: · · E · E EN A M = nvp mvp = r (T ) , (6) nvp @ r(T )m = , r(T )m (8) · where M B is the amount of the dry vapour, generated over B vp B vp where NA is the Avogadro constant, mol-1; μvp is molar mass of diethyl ether vapour, kg/mol. · the evaporator per time unit, kg/s; nvp is the velocity of growth of the number of vapour molecules over the evaporator per time unit, s-1; mvp is the mass of molecule of diethyl ether, kg; r(TB) is specific vaporization heat of the diethyl ether in HPs, depending upon temperature and Assuming approximately equal velocities and without considering the jet type nature of the flow pattern, the linear velocity of the hydrodynamic flotation of the vapour flow over the evaporator surface is calculated from the equation: · pressure, J/kg; TB(P) is the boiling temperature of the diethyl ether, K. uvp @ nvp F (z)n (T = ) F (z)r(T E )m n = (T ) vp ev B vp vp ev Mass flow of saturated dry monomolecular vapour over the evaporator is defined by the following equation: · = EN A . F (z)r(TB )mvpnvp (Tev ) (9) Gvp = M = F (z)nvp (Tev ) mvp uvp = E Insert expressions (9) and (8) in the equation (7) to obtain the formula for calculation of molecular flow = , = F (z)rvp (Tev )uvp r (TB ) (7) of the dry vapour over the evaporator: EN where Gvp is a mass flow of dry saturated vapour over the Gvp @ rvp (Tev ) A . r(TB )mvpnvp (Tev ) (10) evaporator, kg/s; nvp(Tev) is an average number of dry vapour molecules in the unit volume over the evaporator, m-3; uvp is an average velocity of the vapour phase, m/s; ρvp(Tev) is vapour density of diethyl ether over the evaporator, kg/m3; Tev is the temperature of the evaporator, K. Hydrodynamic flow of saturated dry vapour is defined by pressure difference value between the evaporator and condensation region in HP, according to the formula: rvp (Tcond ) F (z)2 (P(Tev ) - P (Tcond )) Growth rate of the number of molecules of vapour over the evaporator in the convergent region of the nozzle of the vapour channel, which defines excess pressure over Gvp @ A hL F (z)2 DPvp = A hL , = (11) where A is nondimensional constant about a unit; η is coefficient of dynamic viscosity of dry vapour, Pa·s; L is the length of the vapour channel of the HPs, m; P(Tev) is vapour pressure under the HPs evaporator, Pa; P(Tcond) rvp (Tcond ) F (z)2 r(Tev )2 mvpnvp (Tcond ) E @ ´ NAhLTcond ´(Tev -Tcond ) = KHPDT , (17) is vapour pressure near the condensation surface of the HPs, Pa; ΔPvp is vapour pressure difference over the evaporator and near the condensation surface in the vapour channel of the HPs, Pa; Tcond is temperature of the condensation surface, K. Equate mass and hydrodynamic flow of saturated dry where nvp(Tcond) is an average number of dry vapour molecules in the unit volume near the condenser, m-3. Heat transfer coefficient KHP at cross-section of the vapour channel of the HP is defined according to the fol- lowing expression: rvp (Tcond ) F (z)2 r(T )2 mvpnvp (Tcond ) vapour, and obtain the following equation: r (T ) F (z)2 (P(T ) - P (T )) KHP @ NAhLTcond . (18) Gvp @ A vp cond ev cond = hL When the rate of heat supply into the thin evaporator is high, its average temperature is greater than the boiling temperature of the working fluid during bubble boiling or = rvp (Tev ) EN A . r(T )m n (T ) (12) vaporization. The hydrodynamic vapour flow at the con- verging region of the nozzle has no time to carry out heat- B vp vp ev Excess pressure over the evaporator defines transfer of the vapour flow in the vapour channel of the HPs, and is calculated in linear approximation according to the equation: ev cond ev cond P(T ) @ P (T ) + dP (T -T ). (13) dT With respect to the Clapeyron-Clausius equation, the pressure derivative of the vapour by temperature is calcu- lated in a conventional manner, nevertheless considering the fact that for liquid the specific volume ratio is small, νL/ νvp < 10-2-10-3, hence in the Clapeyron-Clausius equa- tion the value of the specific volume of fluid νL is ignored, and in ideal gas state the following equation is obtained: dP = 1 r(T ) @ 1 r (Tcond ) @ ing capacity generated through boiling in the evaporator. Vapour density and pressure increase and boiling temperature of the working fluid goes up until it exceeds average temperature of the evaporator. Upon pressure increase the boiling in the evaporator stops (slows), and an overpressure wave spreads in the vapour channel to the condensation region of the HPs, where the vapour becomes supersaturated and condenses. The condensation process is not instant, and when the evaporator’s activity is slowed (stopped), condensation lasts until the pressure decreases to the vapour saturation pressure at the conden- sation temperature, thereafter the condensation process stops. The slow process of saturated vapour pressure decrease, due to condensation, feeds back through the vapour chan- nel of the HP back to the evaporator, and the boiling dT T (nvp -nL ) Tcond nvp process resumes there. Pressure pulsations in the vapour channel of the HPs @ rvp (Tcond ) r (Tcond ) , Tcond (14) result from boiling and intensive vapour generation in the evaporator, non-instantaneous mass-transfer through the vapour channel to the condensation region of the HP, and where νvp is the specific volume of the saturated vapour, m3/kg; νL is the specific volume of the diethyl ether on the liquid-vapour coexistence line, m3/kg. Substitute expression (10) in (8), and obtain the equa- tion to calculate vapour quantity in the HP: T - T @ F (z)2 r (Tcond )rvp (Tcond ) ( ) the slow condensation process, which provides pressure decrease in the condensation region at first and then in the evaporator of the HPs. Following all described processes the next pulsed evaporation cycle resumes. The heating capacity, entering to the flat grid evaporator of the short HP, when the evaporator’s temperature exceeds the boiling temperature hL @ EN A ev cond TB(P) of the working fluid, is defined according to the next equation: r (T )m n . (T ) (15) [Tev -TB (P) ] F (z) = E . (19) cond vp vp ev Vapour temperature over the evaporator’s surface with low evaporation and without boiling is defined by the equation (16): Tev @ Tcond ´ Rev (T ) It is considered that moist vapour stream is formed of two subsystems: microdrops stream system Gdr and dry vapour stream system, Gvp. Rate of vaporization of the moist droplet vapour flow is defined in the standard way: æ EN hL ö · ´ç1 + A ÷ £ TB . (16) M = Gmix = Gvp + Gdr . (20) ç r (T ) F (z)2 r (T )2 m n (T ) ÷ è vp cond ev vp vp ev ø Steady-state evaporation conditions in HPs mean that the temperature in the evaporator does not exceed boiling temperature of the working fluid. Heating capacity of the HPs, W, is defined according to the following equation: For the purpose of simplification of the analytical model construction, the real droplet vapour flow over the evaporator, having microdrops dimensioned by compli- cated double-humped distribution function, Lee, Reges, Almenas, [16], is reported in terms of a mono-dispersal system of spherical microdrops with arithmetic middling radius ra, which is frequently used while analyzing two- phase droplet vapour flows: G = G + G rmix F (z)DP @ vp vp . (27) = 1 ¥ r år n , Подпись: a (21) mix vp dr hmix L ndr i=0 dri dri As a result of the boiling process in the evaporator, the where ndr is the total number of microdrops of all sizes per unit volume of the droplet vapour flow over the evaporator, 1/m3 ; ndri is the number of microdrops of the diethyl ether with radius rdri per unit volume of the droplet vapour medium, 1/m3. Considering the accepted assumption over spherical shape of microdrops, the expression for absolute moisture pressure of the vapour over the evaporator increases up to P*, whereby the boiling process in the surface layers and further in the whole of the thin evaporator (3 mm thick) is slowed (stopped) due to the fact that average temperature of the evaporator Tev becomes lower than the boiling tem- perature TB of the working fluid in the evaporator under the increased pressure and the confined space: γ of the droplet vapour flow is defined in the standard way: Tev < TB (P* ). (28) M é r æ 3 1 öù-1 Furthermore, pulsing of the overpressured vapour g = dr = ê1 + vp ç - 1÷ú , (22) begins to spread to the condensation region through the è a dr ø M dr + Mvp êë rdr ç 4p r3n ÷úû where Мdr is total mass of microdrops per unit volume of the droplet vapour flow over the evaporator, kg; Мvp is vapour mass per unit volume of the dry vapour flow over the evaporator, kg; ndr is the total number of microdrops of all sizes per unit volume of the droplet vapour flow over the evaporator, 1/m3; ndr is the total number of microdrops of all sizes per unit volume of the droplet vapour flow over the evaporator, 1/m3; ρdr is density of r (T @ @ ) microdrops of the diethyl ether, kg/m3. vapour channel. Cessation of boiling and retarding of va- porization in the evaporator of the HP lead to significant reduction (cessation) of the heat release and decrease in moist vapour transportation along the vapour channel of the HP to the condensation region of the HP. The time period Δτev of pressure increase up to P* and cessation of boiling in the capillary-porous evaporator is estimated in linear approximation in ideal gas state and laminar heat transfer inside of the vapour channel of the HP, according to the formula: The density of moisture vapour mix vp ev without Dtev éëP* - P (Tcond )ùû F ( z ) L taking into account the relative motion phases is: kBTevnvp r mix vp (T ) = rdrrvp . ev gr + (1 - g r) vp dr (23) éëP* - P (Tcond )ùû F ( z ) Lr (TB ) mvp EkBTev (29) @ . Synergies between consumable and absolute mass concentrations of microdrops γG or consumable and abso- lute moisture levels are defined according to the following relation: Surplus energy (increased pressure) release time ΔτHP in the evaporating region, by means of vapour flow trans- fer to the condensation region of the HP, is estimated according to the following formula: g = gy ; y = udr . (24) Dt @ EDtev . (30) G (1-g) + gy uvp HP r (T )rmix (T )u (T ) F ( z ) The value γG represents the relation between the quantity of the condensed droplet phase and the total quantity of the two-phase droplet vapour, and, considering the expression (22), is as follows: B vp ev ev Surplus pressure release time ΔτHP partially determines the time value of oscillations of the vapour flow in the vapour channel of the HP, during which the vapour pressure pulse, arisen over the evaporator, reaches the HP é r æ 3 1 öù-1 G condensation surface, and condenses partially. The time Подпись: G g = ê1 + vp ç - 1÷ú = dr . (25) period Δτcond of pressure decrease up to P(Tcond) is esti- ê r y ç 4p r3n ÷ú G ë dr è a dr øû mix mated according to the formula: Mass flow rate of microdrops on the evaporator’s surface is considered to be proportional to vaporization velocity and mass flow rate of the vapour: Dtcond @ EDtev vp r (Tcond )rmix (Tcond )u (Tcond ) F ( z ) . (31) æ E öa æ G öb As a result of the liquid phase formation, pressure near dr G @ B ç ÷ ç vp ÷ , (26) the cooling condensation surface reduces to: è r(Tev ) ø è F ø P* @ P (Tcond ). (32) where a and b are numerical coefficients; B is a coeffi- cient with account for thermophysical properties of the diethyl ether and structural parameters of the evaporator (porosity, typical dimension of channels and pores). To evaluate the vapour quantity transported from the evaporator to the condensation region of the HP, equation This leads to the slowing of heat transfer through the vapour channel, the rarefaction wave propagation from the condensation zone to the evaporator and the start of the next cycle of pulsed increase of the pressure near the evaporator’s surface. Thus, cycle pulse duration Δτ0 in the vapour channel of the HP is: (31) is used, where the microdrops subsystem contribu- Dt0 @ Dtev + DtHP + Dtcond + Dtsound . (33) tion to the generation of overpressure of the two-phase droplet mixture over the evaporator P(Tev) in the vapour channel is neglected: The analysis of the obtained measured data of heat transfer over the evaporator shows that the designed HPs work in the boiling mode. This fact is confirmed by influ- ence of operating parameters of the vaporization process (heat flow density q and pressure P(Tev)) on the heat transfer coefficient α. The influence of these parameters is related to the similar influence of heat flow density and pressure over the bubble boiling in the substantial volume. Performed numerical estimation of parts of pul- sation duration Δτ0 from the expression (33) show the following. Pressure increase duration Δτev over the HPs evapora- tor leads to increased values Δτev ~ 10-1-10-2 c, which means incomplete adequateness of flat model of estima- tion of evaporability in the evaporator. It is necessary to con- sider nonequilibrium volumetric vaporization with regard to the structure and porosity of the evaporator, tempera- ture and moisture gradient in the boiling two-phase work- ing fluid. Moving duration ΔτHP of the vapour cluster along the convergent-divergent HPs vapour channel at vapour flow velocity in the channel equal ~ 100-300 m/s, reaches ΔτHP ~ (4.8-1.6)·10-3 s, which is acceptable only in the case when sonic velocity in dry vapour is ~300 m/s and correspondingly large temperature gradients in the vapour channel. Maximum duration of condensation period of moist vapour is also excessive and reaches Δτcond ~ 10-1-10-2 s. Duration of expansion wave motion from condensa- tion region to the HP evaporator at sonic speed in moist vapour about 100 m/s does not exceed Δτsound ~ 1·10-3 s. Comparison of experimental measured results of pul- sation period duration in the convergent-divergent vapour channel of short low temperature HPs Δτ0 = (2-2.5)·10-3 s with the valued Δτ0 ~ (0.75-1.5)·10-3 s obtained using numerical method and the software CFD Design 10.0 shows that they have good fit values. Comparison of experimental measured results of pul- sation period duration in the convergent-divergent vapour channel of short low temperature HPs Δτ0 = (2-2.5)·10-3 with estimated values obtained using analytical method shows excessive values of analytical results. Suggested simple theoretical model of the flat surface evaporation gives excessive values of pulsation period duration and does not fully correspond to acutely nonequilibrium evaporation processes in the grid capillary porous evapo- rator with boiling working fluid. Conclusions. 1. Heat transfer coefficient of the short linear HPs with the Laval-like vapour channel reaches the value of (22 ± 1)·103 W/m2K; heat transfer coefficient of the HPs with the cylindrical vapour channel does not exceed (15 ± 1)·103 W/m2K, if thickness of the capillary-porous insert layer is 3 mm, diameter is 20 mm and length is 100 mm. Comparison of heat transfer coefficients in short HPs shows sufficient advantages of HP with the Laval-like vapour channel in contrast with standard cylindrical channel. 2. The pulsatile flow regime in the vapour channel of short linear HPs occurs when the working fluid in the evaporator starts to boil. Vapour density and pressure increase and boiling temperature of the working fluid goes up until it exceeds average temperature of the evapo- rator. Upon pressure increase, the boiling in the evaporator stops (slows), and an overpressure wave spreads in the convergent-divergent vapour channel to the condensation zone of the HPs, where the vapour becomes supersatu- rated, and condenses. Condensation of moist vapour results in pressure decrease in condensation region, and after arrival of the expansion wave at the HPs evaporator, pressure over the evaporator also decreases, the boiling process in the evaporator proceeds and the pulsation cycle in the vapour channel repeats. 3. Flow stagnation leads to vortex formation, and their interaction causes pulse-coupled vortex decay, static pres- sure boost and complex reverse flows. Calculation results show that flow stagnation during the pulsation cycle leads to an enlargement of the recirculating region and augmen- tation of condensation, and this effect is significant. With all the operating parameters of heat pipes, two effects were observed: pulsating flow regime of two- phase vapour flow and film-type condensation.×
Sobre autores
A. Seryakov
LLC Research and Development Company “Rudetransservice”
Email: seryakovav@yandex.ru
55, Nekhinskaya Str., Veliky Novgorod, 173025, Russian Federation
Bibliografia
- Pat. 4921041 US. Structure of Heat Pipe / Akachi H. 1990.
- Tong B. Y., Wong T. N., Ooi K. T. Closed-loop pulsating heat pipe // Applied Thermal Engineering. 2001. Vol. 21, № 18. P. 1845-1862.
- Experimental Investigation and Performance Evaluation of a Closed Loop Pulsating Heat Pipe / R. Naik [et al.] // Journal of Applied Fluid Mechanics. 2013. Vol. 6, № 2. P. 267-275.
- Gupta A. K., Lilley D., Sayred N. Swirling flow. New York : Wiley, 1987. 588 p.
- Seryakov A. V. Pulsation flow in the vapour chan- nel of low temperature range heat pipes // Direct Research Journal of Engineering and Information Technology. 2014. Vol. 2(1). P. 1-10.
- Seryakov A. V. Pulsation flow in the vapour chan- nel of short low temperature range heat pipes // Interna- tional Journal on Heat and Mass Transfer Theory and Application. 2014. Vol. 2, № 2. Pp. 40-49.
- Seryakov A. V., Ananiev V. I., Orlov A. V. Con- densation research in the short low-temperature range heat pipes // Proceedings of the 9th Minsk International Seminar of Heat Pipes, Heat Pumps, Refrigerators, Power Sources. (7-10 September 2015, Minsk). Vol. 2. P. 168-176.
- Seryakov A. V., Ananiev V. I. Condensation research in the short low-temperature range heat pipes // Proceedings of the VIII Intern. Symposium on Turbu- lence, Heat and Mass Transfer. (15-18 September 2015, Sarajevo, Bosnia and Herzegovina). Begell House Inc. P. 693-696.
- Пат. 2431101 RF, F 28 D 15/00. Способ заполне- ния тепловых труб / Серяков А. В. Опубл. 10.10.2011, Бюл. 28.
- Полезная модель 152108. Емкостный датчик для определения толщины слоя жидкости / Серяков А. В. Опубл. 27.06.2015, Бюл. № 18/2015.
- Серяков А. В., Конькин А. В., Белоусов В. К. Применение струйного парового сопла в тепловых трубах среднетемпературного диапазона // Вестник СибГУ. 2012. Вып. 1(41). C. 142-147.
- Seryakov A. V. Velocity measurements in the va- pour channel of low temperature range heat pipes // Inter- national Journal of Engineering Research & Technology. 2013. Vol. 2, № 8. P. 1595-1603.
- Faghri A. Heat Pipe Science and Technology. Washington : Taylor and Francis, 1995.
- Варгафтик Н. Б. Справочник по теплофизиче- ским свойствам газов и жидкостей. М. : Физматгиз. 1963. 708 с.
- Таблицы физических величин : cправочник / под ред. акад. И. К. Кикоина. Атомиздат, 1976. 1008 с.
- Lee R., Reges J., Almenas K. Size and number density change of droplet populations above front during reflood // International Journal of Heat and Mass Trans- fer. 1984. Vol. 27, № 4. Р. 573-585.
Arquivos suplementares
![](/img/style/loading.gif)