The results of studies by solving the inverse thermal conductivity problem of the heat capacity of evaporator of the short linear heat pipes (HP’s) with a Laval nozzle-liked vapour channel and intended for cooling spacecraft and satellites with strict take-off mass regulation are presented. Mathematical for mulation of the inverse problem for the HP’s thermal conductivity in one-dimensional coordinate system is accompanied by the measurement results using the monotonic heating method in a vacuum adiabatic calorimeter the HP’s surface temperatures along the longitudinal axis over the entire temperature load range, thermal resistance, and arrays of thermal power data on the evaporator Qev and vortex flow calorimeter Qcond for the condensation surface allow us to estimate the average value of the evaporator heat capacity Cev by solving the inverse thermal conductivity problem in the HP’s evaporator region. Since at the beginning of working fluid boiling for a certain time interval, the temperature of the capillary-porous evaporator remains close to constant, and with the continuation of heating and by solving the inverse thermal conductivity problem, it becomes possible to calculate the heat capacity of the working evaporator and the evaporation specific heat of the boiling working fluid and compare it with the table values.
The current development level of the computer technology and computing technologies allows us to significantly expand the class of applied problems to be solved. Those scientific and technical problems of the HP’s heat transfer that have traditionally been tried to be considered analytically are increasingly being analysed and solved with the help of numerical methods using specialized software for engineering calculations. We are talking about the application of computational programs for solving inverse problems of thermal conductivity of solids in one-dimensional formulation and the use of the obtained methods for analysing the operation of short HP’s with a Laval nozzle-liked vapour channel and a large amount of working fluid in the capillary-porous insert and evaporator. Numerical methods have received the greatest application due to a number of specific advantages, the main of which is their relatively simple implementation on computers [
Inverse problems are characterized by instability of their solutions, which manifests itself in the occurrence of large numerical changes in the solution with small changes in the initial data.
Tasks of this type are called incorrectly assigned tasks or ill-posed problems [
The general nonlinear problem of non-stationary HP’s thermal conductivity can be considered as a set of three nonlinear problems: a problem in which the nonlinearity arises due to the temperature dependence of the coefficients of the main equation; a problem in which the boundary conditions are nonlinear, for example, due to evaporation or boiling of the working fluid in the evaporator; and a problem in which the nonlinearity arises due to the temperature dependence of internal heat sources (sinks). The boundaries between these problems are very conditional, since some transformations allow you to switch from one type of nonlinearity to another, which sometimes simplifies the use of certain methods for solving nonlinear problems.
The main difficulty in the approximate solution of ill-posed problems is the choice of regularization parameters to achieve the stability of the solution. To determine it, the following approaches are most widely used [
When solving inverse problems for equations of mathematical physics, gradient iterative methods are widely used in the variational formulation of the inverse problem [
In addition the heat source in the HP’s evaporator depends only on the time τ, and the temperature field T(z) in the entire geometric space of the HP inside the adiabatic calorimeter is experimentally determined using copper-constantan surface thermocouples. These measurements guarantee that the inverse problem has an unambiguous solution, but this solution is unstable; therefore, additional regularization and the use of the variational method [
A simple algorithm for the numerical solution of a one-dimensional inverse problem by the coefficient of thermal conductivity for calculating the heat capacity of a short HP’s evaporator with a boiling working fluid has been developed and implemented in the FORTRAN system for a PC. In this case, an important factor is the temperature of the external surface of the capillary-porous HP’s evaporator, which is close to a constant value. With further heating, the level of boiling working fluid in the HP’s evaporator slowly decreases, and the temperature of the external surface also slowly decreases. And this allows us to calculate the extreme behavior of the evaporator heat capacity and estimate the specific heat of working fluid boiling in them.
In addition, our simple algorithm is convenient for its implementation in the measuring scheme of an adiabatic vacuum calorimeter for the study of short HP’s. The capabilities of this algorithm are illustrated by solving several specific problems for calculating the heat capacity of short solids under monotonic heating at moderate temperatures, including a model problem with the evaporation of diethyl ether and, consequently, with variable initial conditions.
Thus, the main purpose of this work is to expand the research capabilities of HP heat pipes, solve inverse problems of thermal conductivity with their help and clarify the thermophysical characteristics of working fluids when boiling occurs in an HP capillary-porous evaporator. The proposed new application of the developed mathematical ideas makes it possible to significantly increase the scope of HP and raise the depth and vastness of the results obtained.
HP’s Thermophysical AnalysisWhen analyzing short HP’s, the thermophysical characteristics are represented by the so-called effective properties, and the process of heat transfer along the vertically oriented z coordinate is described by the thermal conductivity equation:
ρ H P С e f f ( T ) L H P F H P ∂ T ∂ τ = ∇ ( λ e f f ( T ) ∇ T ) . (1)
The specific heat capacity of the HP’s is determined additively by the formula (2) when the components and materials are characterized by the absence of chemical interaction between them.
C ¯ H P ( T ) = ∑ j ψ j C j ( T ) . (2)
where ψj is the mass fractions of the HP’s components, Cj is their heat capacity’s. The effective thermal conductivity λeff as an additive sum thermophysical characteristics of the components can be presented only in simple cases, in particular for HP’s in a non-working state.
During the phase transformation of the fluid component, for example, during the working fluid boiling process in the HP’s evaporator, Equation (1) must be supplemented with an additional source function:
Y ( T ) = ρ H P r j ( T ) ∂ ψ j ∂ τ = ρ H P r j ( T ) ∂ ψ j ∂ T ∂ T ∂ τ . (3)
where rj(T) is the specific heat of the phase transformation of the j – th component.
Taking into account (2), Equation (1) can also be used in the presence of phase transformations in the HP with a new average heat capacity:
C ¯ H P ( T ) = С e f f ( T ) − r s p ( T ) .
The heat conduction equation in one-dimensional Cartesian coordinate system for calculating the heat propagation in a vertically oriented short linear HP’s with a cylindrical body along the longitudinal z-axis, during monotonic heating in an adiabatic calorimeter without heat losses through the side walls is written as follows:
ρ H P С v ( T ) L H P F H P T ˙ ( z , τ ) = L H P F H P 1 z ∂ ∂ z 1 R H P ( T ) ∂ t ( z , τ ) ∂ z ; R H P ( T ) = L H P λ H P ( T ) F H P . (4)
The process of heating and heat transfer using HP’s is endothermic process, parameter Y(T) has a negative value, Y(T) ≤ 0, which implies the principle of maximum and the uniqueness of the thermal conductivity equation solution [
The measurement part of our method consists in a system for measuring the thermal power entering to the evaporator, the thermal power entering in the vortex flow calorimeter and the temperature distribution of the HP’s surface inside the adiabatic calorimeter. In the case of using a vacuum adiabatic calorimeter with slow monotonous heating, the temperature distribution along the HP’s body makes it possible to distinguish characteristic areas (zones) of the HP’s, for example, the zone of the evaporator, the temperature in which, under a high temperature load, is close to constant for a fairly long boiling time (several tens of minutes), the time when the level of the working fluid in the capillary-porous evaporator decreases from the initial value h = 3.5 mm to the final value h = 0.02 – 0.01 mm and even less than this value.
The monotonic heating mode is extremely informative for the study of the short HP’s with the phase transformation (boiling) of the working fluid in the evaporator [
At the working fluid boiling beginning during a certain time interval, the temperature of the capillary-porous evaporator remains close to constant and with the continuation of heating and by solving the inverse problem of thermal conductivity, it becomes possible to calculate the heat capacity of the working evaporator and the heat of evaporation.
The mathematical description of the monotone heating mode contains and implies several simplifying assumptions, which as a result lead to analytical relations for the direct calculation of the IPTC.
Heating mode contains and implies several simplifying assumptions, which as a result lead to analytical relations for the direct calculation of the IPTC. However, it is necessary to introduce some specific adjustments to take into account several interfering factors. It is necessary to take into account the final longitudinal dimensions of the HP, which lead to a distortion of the HP’s temperature field due to the heat exchange of the condensation end surface and the nearby side section with the cooling running water in the calorimeter, while analytical expressions imply the one-dimensionality of the temperature field in the HP.
The relatively large noticeable temperature differences within the HP during the boiling of the working fluid in the evaporator also require certain adjustments to the calculations. The values of the corrections will (can) be estimated in a special way.
The direct problem of the HP thermal conductivity is to find the HP’s temperature (surface temperature for the thin HP’s in a one-dimensional model of heat propagation in an adiabatic calorimeter) that satisfies the differential equation of thermal conductivity and has the specified boundary and initial conditions of unambiguity [
The solution of the direct problem for the HP’s evaporator in the case of large temperature load and the boiling process beginning in the temperature region of the working fluid phase transition is accompanied by a nonlinear dependence of the intrinsic properties (heat capacity CHP and thermal resistance RHP) to temperature and leads to large error as the averaged results and the ambiguity of the resulting solutions.
At the same time, the average values of the obtained heat capacity results CHP are not of great interest and make it absolutely necessary to isolate the zone of the HP’s evaporator at the working fluid boiling beginning, thereby allocating the heat capacity Cev of the HP’s evaporative fragment.
The inverse problem of thermal conductivity (IPTC) in this case is to determine the parameters of internal heat transfer also in a one-dimensional model of heat propagation in an adiabatic calorimeter, the thermal resistance RHP (or thermal conductivity λHP) and the heat capacity Cev of a short linear HP’s evaporator at high thermal loads and the working fluid boiling beginning.
IPTC is an incorrect problem in the sense of Hadamard [
Formulation and solution of coefficient inverse problems of heat conduction in a one-dimensional coordinate system with one-dimensional temperature field are the theoretical basis of the main experimental methods for the evaporator’s heat capacity study of short linear HP’s (HP’s evaporative fragment), using the results of thermal resistance measurement and the transferred heat power qev along the longitudinal axis to reflect a decrease in transmitted power due to losses on internal friction and heat transfer losses due to the nonabsolute adiabatic mode in the calorimeter [
Thus, we have a coefficient inverse problem of thermal conductivity, when the known experimental data on the thermal resistance RHP and the side surface temperature THP in adiabatic conditions, as well as the measured heat flow on the surface of the evaporator Qev (and condensation Qcond), it is possible to restore the values of the heat capacity Cev of the evaporative fragment depending on the temperature load δt.
To conduct the experimental studies of the linear HP’s thermal characteristics, we used previously developed and used in many previous studies short stainless steel HP’s with a Laval nozzle-shaped vapour channel, the detailed description of which was repeatedly given in previous publications [
The capillary-porous evaporator 7 together with insert 4 forming a unified hydraulic working fluid delivery system, are made of layers of a thin stainless steel mesh with 0.07 mm thick each layer, with a cell size of 0.04 mm. The diethyl ether C4H10O is used as the working fluid, which has the boiling temperature under the atmospheric pressure of TB = 308.65 K (35.5˚C), freezing temperature TF = 156.95 K (−116.2˚C) and critical parameters TC = 466.55 K (193.4˚C), PC = 3.61 MPa. The volume of the insertion pores is determined during manufacture, and in our case it is equal to 16.62 × 10−6 m3, mass of diethyl ether with the density of 713.5 kg/m3 (20˚С) in the insert pores is equal to 11.858 × 10−3 kg. The restriction in the HP’s filling amount is made so that the evaporator is not flooded with ether at a high temperature load and the film boiling beginning on the lower cover surface of the HP (the evaporator surface).
The charging ratio of the HP with Laval-liked vapour channel (ratio of the diethyl ether volume to the total volume of the HP) is equal to 16.62 × 10−6 m3/3.14 × 10−5 m3 = 0.529.
The length of our HP’s is LHP = 0.1 m, diameter d = 2 × 10−2 m, when filling it with diethyl ether the average density is ρ ¯ H P = 1.871 × 10 3 kg / m 3 , ρ ¯ H P ~ ρ e v , the average isochoric heat capacity C ¯ H P = 1.15 kJ / kg ⋅ K is a slightly varying temperature function due to the uniform distribution of the working fluid diethyl ether in the condensation, transport and evaporation areas of the HP and close to the total specific heat capacity Ceff(T) of diethyl ether and stainless steel with the mass fractions of the components (2).
The defining parameter of the entire capillary-porous insert and the evaporator is a small longitudinal interlayer gap (clearance) between layers of the metal mesh, which does not exceed 0.007 mm, and developed by a porous system the capillary pressure, which is sufficient for the HP’s work when the external inertial effects.
The last row of
The orientation of the longitudinal pores (clearances) with the liquid diethyl ether coincides with the skeleton wires of the capillary-porous insert and the direction of the heat flux propagation, therefore, the effective value of the thermal conductivity coefficient of a diethyl ether-saturated evaporator and insert can be calculated as follows:
λ e v = Π ⋅ λ l + ( 1 − Π ) ⋅ λ s c . (5)
At the porosity value П = 0.72 and values λl = 0.136 W/m K, λsc = 24 W/m∙K, λev = 6.82 W/m∙K. The evaporator temperature Tev (maximum value) is equal to the Tev = 329.75 ± 0.01 K (56.6˚C).
The measuring calorimeter for the HP’s studies is shown in
The HPs evaporators, also equipped with 0.1 mm diameter copper constantan wire thermocouples, is heated using a flat resistance heater, and the temperature is maintained at δТ, K higher than the diethyl ether boiling temperature of 308.65 K under atmospheric pressure. The heater temperature is stabilized and HP’s evaporators overheat value is set in the range of δТ = 0 ÷ 20 К, herewith heat power of single HP does not exceed 200 W.
| Average pore diameter, μ | Cell size, μ | The average diameter of the wire, μ | Porosity, % | The height of capillary rise of diethyl ether, mm | Permeability, Kwick, 10−10 m2 |
|---|---|---|---|---|---|
| 120 | 400 | 15 | 78 | 8 | 7.25 |
| 30 | 30 | 10 | 72 | 81 | 0.2 |
| 7 | 30 | 10 | 72 | 400 | 0.15 |
To reduce heat loss when working with HPs, they are placed in a stainless steel vacuum chamber 22 (10−3 torr), where they are further surrounded by thin-walled copper adiabatic screen 23, the inner surface of which is covered with a layer of nickel, and on the outer surface placed 4 sections guard heaters 24. The value of nonadiabaticity near the middle of the HP’s does not exceed ~2 × 10−2 K, in the field of the resistive heater H2 and HP’s evaporators nonadiabaticity is about ~10−1 K.
The main HP, called measuring, is filled with diethyl ether and the reference one, which is completely identical to the main HP, is filled with dehumidified air at a pressure of 1 bar with dew point temperature lower than 233.15 К (−40˚С). The heat transfer coefficient KHP2 of the second HP does not exceed 0.15% from the first one (measuring one) KHP1 and is not taken into account. The second HP, completely identical to the first one, performing the reference function in measurements the condensate film thickness in the first HP [
Conducted measurements of the condensate film thickness on the internal surface of the HP upper cover, using capacitive sensors and developed high frequency generators [
The most general and informative non-stationary energy equation for the work description of the HP’s flat evaporator in the monotonous heating mode is the Fourier-Kirchhoff equation [
C e v ρ ¯ e v L e v F ( z ¯ ) ( ∂ T ¯ e v ∂ τ + v ⋅ ∇ T ) = L e v F ( z ¯ ) ( ∂ ∂ τ + d i v ( v ) ) P e v + d i v ( λ H P ∇ T ) + ∑ i = 1 i = 2 μ i Φ d f i . (6)
The dimensionless length z ¯ and the average heat capacity Cev of the HP’s evaporator calculate in the usual way:
C e v = C s c + x e v C v p + ( 1 − x e v ) C f e v + C w ; z ¯ = z L H P . (7)
With an increase in the temperature load in the initial heating period of a capillary-porous HP’s mesh evaporator the wetting of the grid frame with diethyl ether deteriorates sharply, since the wetting angle increases with an increase in temperature. This effect is thermosensitive and already with a relatively small increase in temperature δT = Tev – TB ~ (3 - 5) K wetting deteriorates, the level of diethyl ether in the grid evaporator decreases and further HP’s heating δT = Tev – TB ≥ 10 K through the flat surface of the bottom cover occurs in a thin layer of boiling diethyl ether.
The thickness of the diethyl ether layer decreases from ~3.5 mm at the heating beginning to the most intense boiling, while the vapour flow becomes stationary, and the thickness of the ether micro-layer does not exceed (1 - 2) × 10−2 mm. This fact greatly simplifies the conduct of assessments and the maximum achievable result in our conditions is as follows:
Q e v = − λ H 2 ( ∂ T H 2 ∂ z ¯ ) z ¯ = 0 = E H 2 − C H 2 T ˙ e v = ( T H 2 − T ¯ f e v ) F e v ( z ) / ( δ c λ c + δ w e v λ w + δ f e v λ f e v + 1 α e v ) = 118.4 W . (8)
Δ T c = Q e v δ c λ c F e v = 5.4 K , Δ T w e v = Q e v δ w λ w F e v = 18.8 K , Δ T f e v = Q e v δ f e v λ f e v F e v = ( 35 … 0.55 ) K . (9)
The equation of the thermal energy of a working HP’s evaporator in a monotonic heating mode in the approximation of small heat losses k H P s h ( T ¯ H P − T s h e l ) in the adiabatic calorimeter can be represented using enthalpy equation as follows:
Q e v ( τ ) = H e v ( τ ) + ∑ i = 1 i = 2 μ i Φ d f i + k H P s h ( T ¯ H P − T s h e l ) . (10)
and the actual enthalpy equation of an evaporator in a monotonous heating mode also can be represented using this thermodynamic equation:
H e v ( τ ) = G v p ( τ ) r ( T B ) + G v p ( τ ) С v p ( T e v − T s c ) + G l ( τ ) C p l ( T e v − T ¯ f e v ) . (11)
The expression for the thermal energy released in the HP’s condensation region is written similarly using the heat balance equation:
Q c o n d = G v p r ( T B ) d x e v d z ¯ + G v p ( 1 − x e v ) C v p d T f e v d z ¯ + λ s c ( 1 − Π ) d 2 T s c d z ¯ 2 F e v ( z ) L e v − ∑ i = 1 i = 2 μ i Φ d f i − k H P s h ( T ¯ H P − T s h e l ) . (12)
The performed measurements and calculations make it possible to estimate the heat transfer coefficient in the vortex flow calorimeter αcal = 2.3 × 103 W/m2∙K, αcond ≥ 5 × 103 W/m2∙K, and the final expression for the maximum value of the condensation heat in our HP is the follows:
Q c o n d = ( T ¯ f c o n d − T c a l ) F c o n d / ( δ f c o n d λ l + δ w λ w + 1 α c a l + 1 α c o n d ) = 115.4 W .(13)
The thickness parameters in parentheses are equal to the following values:
δ f c o n d = 10 − 1 - 10 − 3 mm , δ w = 1 mm , α c a l = 2.4 × 10 3 W / m 2 ⋅ K , α c o n d ≥ ( 5 - 10 ) × 10 3 W / m 2 ⋅ K .
To solve the complex Fourier-Kirchhoff Equation (6), it is necessary first to determine the velocity field of the liquid and vapour phases in a capillary-porous evaporator using the Navier-Stokes equations. The results of this study will be presented later in a separate report.
To simplify all subsequent calculations, the average cross-section HP’s temperature at any point at altitude z ¯ i = z i / L H P inside the vapour channel can be considered equal to the surface temperature of the HP, measured at this height T ( z ¯ i , τ ) s u r in the vacuum chamber of the adiabatic calorimeter:
1 F ( z ¯ ) ∫ T H P ( z ¯ ) d F i = T H P ( z ¯ i ) s u r = T H P ( z ¯ ) . (14)
The vapour flow can be estimated according to the interphase mass transfer gives in case of small departures from equilibrium by the Hertz-Knudsen equation [
G v p = 2 ξ 2 − ξ ( M 2 π R ) 1 / 2 ( P e v T ¯ f e v − P f c o n d T ¯ f c o n d ) . (15)
where the diethyl ether C4H10O, selected as the working fluid has a molar mass M = 74.1216 g/mol, R is the universal gas constant and ξ is the condensation coefficient, ξ ≤ 1, defined as the ratio of the number of vapour molecules condensing over the total number of molecules which strikes the liquid film surface.
There is a redistribution of the liquid mass of diethyl ether inside the HP, while the heat capacity of the evaporator increases due to the liquid film boiling. In this case, the temperature Tfev of the diethyl ether layer and its dependence on time are determined using external thermocouples that control the temperature of the HP external surface inside the vacuum chamber. The results of the thermocouple measurements are shown in
Earlier calculations of the flow velocity and vapour density of diethyl ether using the Navier Stokes equation system in a short HP’s vapour channel similar to the Laval nozzle [
ρ v p m i x = P R T ; ρ v p m i x ρ v = ( P P v ) 1 / k ; P v ~ P f c o n d ; 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 ) . (16)
where k = 1.31 is the value of the adiabatic index of diethyl ether vapour.
Calculations of moving vapour temperature in a short HP’s vapour channel, made in the form of Laval nozzle, were performed using the program code ANSYS\CFD Fluent 6.3.26 -20090623 [
T = T v ( P P v ) ( k − 1 ) / k . (17)
The calculated values of the vapour temperature along the longitudinal axis of the HP vapour channel at the temperature of lower layer of the diethyl ether in the evaporator T ¯ f e v = 329.75 K ( 56.6 ˚ C ) . The calculation error is 0.3%, the measurement error of the HP’s surface temperature is less than 0.1 K. Heating and all calculations were carried out in the temperature range of the evaporator (298 - 348) K (25˚C - 75˚C).
The outer surface of the HP’s body in the region of a 3.5 mm high multilayer mesh evaporator, filled with boiling diethyl ether is characterized by a close to constant temperature under monotonous heating, which is clearly seen in
Inside the vapour channel, the temperature decreases sharply, the temperature drop reaches ΔT = 35 K, due to the cooling in the condensation region and the influence of the Laval nozzle shape of the vapour channel. The results of the experimental analysis of the HP’s surface temperature distribution along the generatrix also confirmed the nonlinear nature of the temperature distribution as a function of the channel length in the Laval-liked vapour channel.
Experimental data shows a close to constant temperature of the evaporator when diethyl ether boils in it, a weak nonlinearity in the confuser part of the nozzle, and a strong nonlinearity in the behavior of the temperature near the HP’s condensation surface. With a further increase in the heat load, the length of the constant temperature region in the evaporator is reduced.
The results of numerical analysis of the vapour axial temperature distribution inside the Laval-liked HP’s vapour channel of compressible moist vapour at high heat loads confirm the nonlinear character of the temperature distribution as function of the channel length, which is due to the joint influence of the channel shape in the form of a Laval-liked nozzle and a sharp gradient vortex formation near the HP’s condensation surface, taking into account an additional temperature decrease inside the vortex ring of the condensing vapour [
Experimental values in the stationary state of the thermal resistance RHP of a short HP’s are defined as the thermophysical characteristic of the heat transfer device as a whole [
The full temperature load on the evaporator δt = Tev – TB = (0 - 20) K [
Experimental values of thermal resistance RHP obtained during monotonic heating of the HP’s evaporator with a speed close to linear in time 3 × 10−3 K/s using a special precision temperature controller and a system of computer-controlled switches, surface thermocouples, heater power measurements form arrays of data on HP surface temperatures in a vacuum adiabatic calorimeter, condensation surface, heater power EH2 and heat flows Qev and Qcond for solving thermal conductivity equations and calculating thermal resistance according to formula (18) from [
R H P ( t ) = T e v − T c o n d Q e v . (18)
The small size and the low heating rates ensured that the measurements of the RHP were carried out under homogeneity conditions and local thermal equilibrium. The obtained values of the thermal resistance RHP, K/W in the continuous heating mode refer to the direction of the heat flow that coincides with heat flow direction in the stationary states measurements.
The polynomial Equation (19) describing in dimensionless form the experimental values of the thermal resistance RHPk for a time moment τk of the short HP with a Laval nozzle-liked vapour channel, depending on the overheating value of the evaporator δt: 0 ≤ δt ≤ 20 (and also in dimensionless form) looks like this:
R H P k ( δ t ) = ∑ i = 1 n R R H P k i ( δ t ) i − 1 = − 2.0795621 × 10 − 8 ( δ t ) 7 + 1.6029662 × 10 − 6 ( δ t ) 6 − 4.9921411 × 10 − 5 ( δ t ) 5 + 8.0489929 × 10 − 4 ( δ t ) 4 − 0.0071936 ( δ t ) 3 + 0.03633406 ( δ t ) 2 − 0.1113127 δ t + 0.2702057 , n R ≤ 8 (19)
Standard deviation σ = 0. 0024929, Fisher criterion R2 = 0.9980243.
The obtained experimental results of thermal resistance at temperature load δt > 10 K allow us to estimate the characteristic total thickness of the diethyl ether film on the evaporator surface (boiling in the evaporator) and the condensation surface of short HP’s with a Laval nozzle-liked vapour channel and flat upper and lower covers as follows (without the thermal resistance of covers):
δ ¯ e v + δ ¯ c o n d ~ R H P λ l F = 4 × 10 − 2 K / W × 0.136 W / m 2 ⋅ K × 3.14 × 10 − 4 m 2 ≤ 1.7 × 10 − 6 m .
From the analysis of Figure(6), it follows that the film thickness on the evaporator is close to a constant value in the range of thermal loads δT = (10 - 20) K, and with a further increase in the temperature load δT, the evaporator begins to dry out and there is a sharp increase in the thermal resistance RHP.
In addition, when the HP’s evaporator is monotonically heated, the diethyl ether layer thickness inside the evaporator δfev(τ) also decreases monotonically and within the range of the coordinate z:
0 < z < Lev = 0.035∙LHP limited by the longitudinal dimensions of the evaporator Lev, the functional dependence of the thermal resistance RHP on the coordinate z: RHP = RHP(z) is of practical interest, and has a positive sign:
∂ R H P ∂ z = ∂ R H P ∂ t ( ∂ t ∂ z ) = ∂ R H P ∂ t ( ∂ z ∂ t ) − 1 > 0 ; 0 ≤ z < L e v = 0.035 ⋅ L H P . (20)
The distributions of the experimental values of the HP’s surface temperature shown in
Let’s formulate a mathematical statement of the inverse problem of restoring the heat capacity of working HP’s. Let there be a vertically oriented HP with the length LHP = 100 mm and diameter of 20 mm, DHP/LHP ≤ 0.2 and the temperature field inside can be viewed in one-dimensional axisymmetric mode and located at the bottom of HP the capillary-porous evaporator with a thickness of Lev = 0.035∙LHP. At the lower end of the HP the evaporator is supplied with a flow of heat Q e v ( z ¯ = 0 ) from a flat heater H2,
The classical equation of spatial thermal conductivity for the linear HP’s as a solid rod (solid body) is looks like this [
ρ H P с v ( t ) L H P F ( z ¯ ) ∂ T ∂ τ = d i v ( λ H P g r a d T ) . (21)
In a working HP with monotonous heating the heat capacity Cev of the evaporator can be considered and analyzed only as the heat capacity of the HP’s lower fragment С e v = С H P ( z ¯ = 0.035 ) with a thickness of Lev = 0.035∙LHP and thermal resistance R e v = R H P ( z ¯ = 0.035 ) .
The distribution of the one-dimensional axisymmetric temperature field t H P ( z ¯ i ) s u r and t ˙ ( z ¯ , τ k ) s u r along the longitudinal dimensionless z ¯ -axis of a short linear HP’s is used to solve the heat conduction Equation (22) for the evaporator fragment without internal sources [
1 z ¯ ∂ ∂ z ¯ 1 R e v ( t ) ∂ t ( z ¯ , τ ) ∂ z ¯ = С e v ( t ) t ˙ ( z ¯ , τ ) ; z ¯ = z L H P ≤ 0.035 . (22)
The thermal resistance of the evaporator (evaporative HP fragment) is an integral part of the thermal resistance RHP of the entire HP, and the value of the thermal resistance of the evaporator during the diethyl ether boiling at temperature load (pressure) > 11 K can be estimated using the experimental results presented in
Energy losses due to friction in the vapour channel and non-absolute calorimeter adiabaticity lead to the fact that the thermal power of Qev and Qcond differ from each other:
( Q e v − Q c o n d ) / Q e v ≤ 2.6 % .
And the heat power becomes a function of the dimensionless vertical coordinate z ¯ and time with monotonous heating:
Q e v = Q e v ( z ¯ , τ ) . (23)
One of the possible formulations of the coefficient inverse problem (CIP) for Equation (23) consists in setting redefined boundary conditions, with the help of which the definition unicity of three functions – temperature t s u r ( z ¯ , τ k ) , heat capacity C e v ( z ¯ , τ k ) and thermal resistance R e v ( z ¯ , τ k ) (or thermal conductivity λ e v ( z ¯ , τ k ) ), moreover all functions must be analytical functions, for example, represented as polynomials or piecewise smooth spline functions. For experimentally determined boundary conditions, the inverse heat conduction problem for Equation (22) is usually incorrect in the classical sense [
The temperature distribution of the cylindrical HP’s housing inside the vacuum chamber with an adiabatic shell near the flat resistive heater is symmetrical:
[ ∂ t ( z ¯ , τ ) ∂ x ] | z ¯ = 0 = [ ∂ t ( z ¯ , τ ) ∂ y ] | z ¯ = 0 = 0 . (24)
The initial boundary conditions [
t ( 0 , τ ) ≡ t e v = T e v ( τ ) − ( R c + R w + R f e v ) Q e v ( τ ) . (25)
t ( 0 , 0 ) = t 0 ; t ( 0 , τ k ) = t 0 + τ k × 3 × 10 − 3 K / s . (26)
t ( L H P , τ k ) ≡ t c o n d = T c o n d ( τ k ) + ( R c a l + R w + R f c o n d + R c o n d ) Q c o n d ( τ k ) . (27)
Q e v ( z ¯ = 0 , τ k ) = E H 2 − C H 2 T ˙ e v ( z ¯ , τ k ) . (28)
The Qcond energy transmitted by HP to the vortex flow calorimeter during vapour condensation is determined using the experimental stand shown in
We write down the dimensional equation of thermal conductivity to estimate the heat capacity of a capillary-porous evaporator with a finite height of 0.035∙LHPi in the following standard form:
z С e v k ( t ) t ˙ ( z , τ k ) = ∂ ∂ z [ L H P R H P ( t ) F ( z ) ∂ t ( z , τ k ) ∂ z ] ; z L H P ≤ 0.035 . (29)
And we introduce a general expression for the heat flow q ( z , τ k ) inside the short HP along the longitudinal z axis:
q ( z , τ k ) = z λ H P ∂ T H P ( z , τ k ) ∂ z = z L H P R H P F H P ∂ T H P ( z , τ k ) ∂ z = E H 2 − C H 2 T ˙ e v ( z , τ k ) .(30)
In the first approximation, in which the heat flow is determined only by the HP’s capillary-porous evaporator without taking into account heat losses, the heat propagation equation in the vapour channel using the thermal resistance RHP of a short HP’s as a whole can be written in the following form:
z L H P R H P ( t ) F e v ∂ T H P ( z ¯ , τ k ) ∂ z ¯ = G v p r ( T B ) d x e v d z ¯ + G v p C v p d T f e v d z ¯ + λ s c ( 1 − Π ) d 2 T s c d z ¯ 2 F e v ( z ) L e v . (31)
The solution of the Equations (29) and (31) is possible only with positive values of time τk and heat capacity C ¯ e v k ( T ) :
τ k > 0 ; C ¯ e v k ( T ) > 0
To simplify all subsequent calculations, it is possible to take the initial temperature value equal to zero:
t ( z ¯ , 0 ) = 0
Equations (30) and (29) of the heat flow propagation q ( z , τ k ) can be rewritten in the dimensionless form q e v ( z ¯ , τ k ) as the heat flow of the evaporator (without taking into account heat losses) along the z-axis at time moment τk and presented as a system of two energy equations for calculating the heat capacity Cevk(t) and the heat flow of the evaporator q e v ( z ¯ , τ k ) , the value of which is related to the thermal resistance RHP of a short HP’s as a whole and for the timeτk can be represented as follows:
z ¯ C e v k ( t ) t ˙ + ∂ q e v ( z ¯ , τ k ) ∂ z ¯ = 0 . (32)
z L H P R H P F ( z ¯ ) ∂ t ( z ¯ , τ ) ∂ z ¯ + q e v ( z ¯ , τ k ) = 0 . (33)
The values of the variable thermal power q e v ( z ¯ , τ k ) on the lower HP’s flat borders, on the upper boundary of the flat evaporator and on the upper flat borders are equal to:
q e v ( 0 , τ k ) = Q e v ; q e v ( 0.035 , τ k ) = Q 0.035 ; q e v ( 1 , τ k ) = Q c o n d . (34)
and Equation (33) for determination the transmitted thermal power by the HP’s capillary-porous evaporator at a height of z = 0.035∙LHP is made using the HP’s surface temperature distribution in the adiabatic calorimeter:
q e v ( z ¯ = 0.035 , τ k ) = − z L H P R H P F ( z ¯ ) [ ∂ t ( z ¯ , τ ) ∂ z ¯ ] | z ¯ = 0.035 = − 0.035 L H P 2 R H P F ( z ¯ ) [ ∂ t ( z ¯ , τ ) ∂ z ¯ ] | z ¯ = 0.035 . (35)
The heat capacity of the evaporator allows us to evaluate the work of the HP. To calculate the
HP’s evaporator heat capacity Cevk, the heat output in the evaporator volume can be considered as reduced by the temperature reduce from the Equation (31) and (35), and an illustration of the experimental temperature distribution data from
0.035 C e v k ( t ) t ˙ = − [ ∂ q ( z ¯ , τ k ) ∂ z ¯ ] | z ¯ = 0.035 . (36)
The inverse thermal conductivity problem (ITCP) solution for the HP’s evaporator is a method of stepwise continuation of the known solution, for example equal to zero at the initial moment, for the next time interval Δτk = (τk, τk+1) and the corresponding temperature interval tk = (tk, tk+1) to calculate the heat capacity of a capillary-porous evaporator with a height (thickness) of 0.03·LHP.
The solution of the equations system (32) and (33) is carried out using experimentally determined surface temperature values t ( z ¯ , τ k ) in the range HP0.03 = (τk, τk+1) × (0, 0.035∙LHP), experimental thermal resistance values R H P ( z ¯ , τ k ) in the range HP1 = (τk, τk+1) × (0, 1∙LHP), and calculated values of the evaporator volumetric heat capacity Cevk(δt) in the range (τk, τk+1) × (0, 0.035∙LHP) and satisfying the thermal conductivity equations (22) and additional conditions (25) - (28), provided that the thermal resistance defined in stationary states is known as a polynomial (19).
Crank-Nicholson calculation scheme,
ρ e v i C e v i T i n + 1 − T i n Δ τ = λ e v i ( T i n + 1 + T i n 2 ) ( T i − 1 n + 1 − 2 T i n + 1 + T i + 1 n + 1 2 h 2 + T i − 1 n − 2 T i n + T i + 1 n 2 h 2 ) .(37)
The coefficientλevi of the evaporator thermal conductivity can be estimated by the Equation (5).
Experimental temperature tk at time τk inside the HP’s vapour channel with monotonous heating relative to the boiling point tB of diethyl ether at atmospheric pressure is determined as follows:
t k = t ( z ¯ , τ k ) = δ t ( z ¯ , τ k ) + t B ; δ t k = δ t ( z ¯ , τ k ) = t k − t B . (38)
At the initial time moment τk and at the initial temperature value tk, we should consider that the HP’s evaporator heat capacity Cevk(δt) with the initial temperature load value δ t ( z ¯ , τ k ) are known and is equal zero, so the heat capacity from now on further heating can be represented as a polynomial with a variable (δt):
C e v k ( t ) = ϕ ( ξ ) ⋅ ∑ i = 1 n c C e v k i ( δ t ) i − 1 ; n C ≤ 10. (39)
Approximation function Φ(ξ) is a mandatory regularizing parameter in the form of monotone increasing function, Φ(ξ) = ξ, or maybe exp(ξ). The value of numerical parameter Φ(ξ) determines the stability and the range of acceptable heat capacity properties, and smoothes out non-linearity’s in a blurred thermal transformation in the HP’s evaporator. With the change from linear to an exponential dependence, the error in calculating the heat capacity decreases markedly, especially when the polynomial heat capacity function leaves the peak, see
The temperature time derivative is a constant value for time values τ > τk and is a coordinating (synchronizing) parameter of all measurements in a time interval Δτk:
t ˙ i > 0 , i = 0 , ⋯ , 2 , τ > τ k ; t k = t 0 + τ k × 3 × 10 − 3 K / s ; Δ t k = Δ τ k × 3 × 10 − 3 K / s . (40)
The experimental distribution of the HP’s surface temperature tk at time τk along the HP’s body is known, are positive and can be written as follows:
t e v ( τ k ) = t k 0 ( z ¯ = 0 ) > 0 ; t 0.035 ( τ k ) = t k 0 ( z ¯ = 0.035 ) > 0 ; t 1 ( τ k ) = t k 0 ( z ¯ = 1 ) > 0. (41)
The accounting and synchronization of the thermal resistance RHPk equilibrium values (33) in a linear slow heating process to the moments of time measurements τk is performed using the value of the temperature load δtk from (38) and (40)-(41). The thermal power from the Equation (30) is calculated as follows:
τ > τ k : Q z ¯ = 0 ( τ k ) = Q e v ; Q z ¯ = 0.035 ( τ k ) = q z ¯ = 0.035 ; q z ¯ = 1 ( τ k ) = Q c o n d . (42)
When the ITCP problem at time τk is solved and the temperature at time τk+1 of the evaporator z ¯ = 0 : t e v ( τ k ) = t k 0 ( z ¯ = 0 ) , at a height z ¯ = 0.035 : t 0.035 ( τ k ) = t k 0 ( z ¯ = 0.035 ) and the temperature of the HP condensation surface z ¯ = 1 : t 1 ( τ k ) = t k 0 ( z ¯ = 1 ) becomes measured and known, this solution becomes the temperature and thermal power initial distribution at transition to the next τk+1 time interval, Δτk+1 = (τk+1, τk+2). From the solution of the problem on the time interval Δτk+1, the function Ck+1(t) is determined on the temperature interval tk+1 = (tk+1, tk+2).
To reduce the number of Ck+1 parameters, defined at each time interval Δτk+1, the crosslink should be performed with the polynomial (39) from the previous calculation at the boundary point of the temperature tk by the function and its derivatives, this greatly simplifies all calculations.
The total number of equations for crosslinking Nc = 2 (heat capacity and temperature derivative of heat capacity), and the maximum degree of the polynomial nc = 10 are primarily the main regularizing parameters of calculation, and were determined in accordance with the estimated accuracy of the heat capacity and experimental error of the temperature measuring.
In addition, the method of numerical solution of ITCP for short linear HP’s using a polynomial expansion of the heat capacity, the size of each time Δτk step is an important regularizing parameter [
The program that controls the operation of the measuring stand,
In this mode, the parameters of the control program are specified such as the duration of the sensors readings measurement, the measurement times of digital voltmeters, oscilloscopes, frequency meters, etc. After the HP evaporator and the vortex flow calorimeter reach stationary isothermal states, the program switches to the “measurement” mode.
The main resistive heater H2 and control system are switched on, and the linear time heating of the HP’s evaporators begins, the adiabatic system and all measuring sensors start working. In this mode, the temperatures distribution on the outer surface of the measuring HP, the evaporator heat power, the condensation surface temperature, the heat output of the adiabatic system and all other thermal characteristics, including the heat transfer characteristics of the HP, are measured. The obtained experimental data arrays are saved and a measurement library is formed.
We analyze the case when the thermal resistance RHP, Equation (19) of the entire HP as a whole is known, have a positive value and it is necessary to calculate the heat capacity of a capillary-porous evaporator with a height of 3.5 mm and saturated with diethyl ether. We integrate Equation (22) with respect to z ¯ and, taking into account the experimental boundary conditions (25), (26), (27) and (28) we obtain a nonlinear integral equation for calculating the evaporator heat capacity at the time moment τk:
∫ 0 0.035 z ¯ C e v k ( t ) t ˙ ( z ¯ , τ k ) d z ¯ = 0.035 q e v ( z ¯ , τ k ) . (43)
To solve such integral equations, the finite difference method is used, based on replacing the derivatives with their approximate values expressed in terms of the functions differences at individual discrete points (nodes) of the finite-difference grid, Equation (37).
We will solve the functional Equation (43) by the iterative method, for which it is necessary to proceed to the variational formulation of this task. We define the target functional as the discrepancy, corresponding to the difference between the left and right sides of Equation (43), and the calculation problem for each time interval Δτk will look like the task of minimizing the discrepancy functional:
∑ k δ C e v k ( t ) = 1 2 ∑ k ∫ τ k τ k + 1 [ ∫ 0 0.035 z ¯ C e v k ( t ) t ˙ ( z ¯ , τ k ) d z ¯ − 0.035 q e v ( τ k ) ] 2 d τ = inf С e v k i ( t ) . (44)
The most suitable way to minimize the functional is the conjugate gradient method [
The Lagrange functional L for minimization the problem (44) with additional conditions (22), (25), (26) and with the values of the heat flow at the boundaries of the evaporator (32), (36) is as follows:
L = [ ∂ ∑ k δ C e v k ( t ) ∂ c i ] + ∫ τ k τ k + 1 [ ∫ 0 0.035 η ( z ¯ C e v k ( t ) t ˙ + ∂ q ( z ¯ , τ k ) ∂ z ¯ ) ] d z ¯ d τ k + ∫ τ k τ k + 1 [ ∫ 0 0.035 χ ( q e v ( z ¯ , τ k ) + z L H P R H P F ( z ¯ ) ∂ t ( z ¯ , τ ) ∂ z ¯ ) ] d z ¯ d τ k . (45)
where the domain of variables τk and q e v ( z ¯ , τ k ) is kept the same as HP0.035 and HP1, and η and χ are Lagrange multipliers corresponding to the first and second Equations (32) and (33).
To obtain the gradient expression of the functional L [
δ L = ∑ i ∫ τ k τ k + 1 [ ∫ 0 0.035 χ С e v k i ( t ) g i δ ( z ¯ L H P 2 R H P F ( z ¯ ) ) ] d z ¯ d τ k − ∫ τ k τ k + 1 [ ∫ 0 0.035 ( z ¯ С e v k i ( t ) η ˙ + z ¯ L H P 2 R H P F ( z ¯ ) ∂ b ∂ z ¯ ) ] δ t d z ¯ d τ k − ∫ τ k τ k + 1 [ ∫ 0 0.035 ( ∂ η ∂ z ¯ − z ¯ χ ) ] δ q e v ( z ¯ , τ k ) d z ¯ d τ k + ∫ 0 0.035 z ¯ η С e v k i ( t ) δ t d z ¯ + ∫ τ k τ k + 1 [ z ¯ L H P 2 R H P F ( z ¯ ) − ∫ 0 0.035 z ¯ η С e v k i ( t ) t ˙ ( z ¯ , τ k ) d z ¯ − 0.03 q e v ( τ k ) ] δ t d τ k . (46)
where the value gk denotes the derivative of the heat capacity components in the following form:
g i = ∂ C e v k ( t , τ k ) ∂ С e v k i ( t , τ k ) ; ∂ ∑ k δ C e v k ( t ) ∂ С i ( t , τ k ) = ∫ τ k τ k + 1 [ ∫ 0 0.035 η t ˙ ( z ¯ , τ k ) g i ] d z ¯ d τ k . (47)
where η, is the solution of the conjugate problem completely corresponding to the original problem (22), (25) and (27). In the reverse time, τ = τk+1 – τk the conjugate problem is written using Equations (32) and (33) and the HP’s thermal conductivity coefficient ak:
z ¯ η k = a k ∂ ( z ¯ ∂ η k / ∂ z ¯ ) ∂ z ¯ ; a k = λ k C e v k ( t , τ k ) = 1 C e v k ( t , τ k ) L H P R H P k F ( z ¯ ) . (48)
Both the initial and boundary conditions of the conjugate task are written in the usual way, and at the upper boundary of the evaporator, the conjugate value is equal to the caloric residual from the square brackets in the Equation (44):
η ( z ¯ , 0 ) = 0 ; η ( 0 , τ k ) = 0 .
η ( z ¯ = 0.035 , τ k ) = [ ∫ 0 0.035 z ¯ C e v k ( t ) t ˙ ( z ¯ , τ k ) d z ¯ − 0.035 q e v ( τ k ) ] 2 ( τ k + 1 − τ ) . (49)
the inverse problem for the heat capacity C e v k ( t , τ k ) at each k-th interval is reduced to finding the minimum of the functional by the method of conjugate gradients.
The solution of the HP’ evaporator thermal conductivity inverse problem (32), (37) and (46) and the solution of the integral Equation (44) for determining the evaporator heat capacity was carried out using the developed program in Fortran [
In the conjugate gradient method, one-dimensional minimization in the permissible direction was carried out using the interpolation-extrapolation method [
For a measuring sample with conductive heat transfer, all details are from Equation (22). For a model sample, the boundary conditions are given below in Equations (50) and (51):
z ¯ = L e v : λ ∂ T ∂ z ¯ = Q c o n d . (50)
| sample type | Φ(ξ) | C/C0 |
|---|---|---|
| measuring sample with heat transfer by the thermal conductivity | Φ(ξ) = ξ, | 3.3 |
| Φ(ξ) = exp(ξ), | 3.3 | |
| model sample with heat transfer by the liquid evaporating and vapour condensing | Φ(ξ) = ξ, | 9.5 |
| Φ(ξ) = exp(ξ), | 11.8 |
z ¯ = 0 : λ ∂ T ∂ z ¯ = Q e v − G e v ⋅ Q v p . (51)
where the evaporation rate and the saturated vapor pressure are estimated as follows:
G e v = c o n s t P v p − P ¯ 2 π R T / M ; P v p = P ¯ exp ( − Q v p R T )
To calculate the heat capacity of the evaporator, a rather complex program was developed in FORTRAN, using computational equations based on the Simpson method [
δ C e v C e v = ( δ t ˙ t ˙ ) 2 + ( δ q e v q e v ) 2 ≤ 0.3 J / K . (52)
The results of the calculated solution of equation (44) in the range of temperature load on the HP’s evaporator δt = Tev − TB = (−5 … +25) K are shown in
The evaporation enthalpy of the working evaporator with boiling diethyl ether was estimated earlier in Equation (11), therefore we can substitute the known table values for diethyl ether [
С e v = H e v t ˙ = ( 345 kJ / kg + 1.75 kJ / kg ⋅ K × 11 K + 2.34 kJ / kg ⋅ K × 16 K ) × 1.57 × 10 − 4 kg / s 3 × 10 − 3 K / s = 21.02 kJ / K . (53)
The ratio of the heat capacity of the evaporator with diethyl ether to the dry evaporator Cev/Cev0 taking into account the heat capacity of the resistive heater CH2 = 1.2 kJ/K at the diethyl ether boiling beginning is as follows:
С e v С e v 0 = С e v 0 + H e v / t ˙ С e v 0 = 1 + H e v / t ˙ С e v 0 = 21.02 + 1.2 1.2 = 18.5 . (54)
which is very close to the maximum numerical value of the heat capacity of the HP’s evaporator.
Equations (53) and (54) confirm the correspondence of the experimental and calculated thermophysical characteristics of diethyl ether located in the evaporator at a high temperature load.
1) The measurements of the surface temperature of short linear HP’s with linear heating of the evaporator in time make it possible to unambiguously record an uneven drop in surface temperature. The HP’s surface temperature values around the evaporator, which are close to constant, allow us to estimate a small drop in the transmitted thermal power in the HP’s vapour channel. The constancy of the surface temperature around the evaporator also allows us to apply the method of solving the inverse problem of thermal conductivity (IPTC), and numerically obtain the extreme value of the heat capacity of the area (fragment) of the HP’s evaporator. This calculated result was obtained for the first time.
2) The HP’s external surface temperature of the capillary-porous evaporator, which is close to a constant value with monotonous heating, confirms the boiling process beginning and allows us to estimate the complex thermophysical parameters of diethyl ether in the evaporator.
3) With step-by-step regularization of the inverse problem solution of the HP’s thermal conductivity and piecewise smooth polynomial approximation of unknown coefficients, it is possible to minimize the error in restoring the evaporator heat capacity using a polynomial-exponential approximation function that restricts the range of acceptable heat capacity values by only positive values.
The author expresses gratitude to Tsvetkov Oleg Borisovich, Professor of ITMO, St. Petersburg, for his active participation in the discussion of this article at the Institute’s conferences. The author expresses special gratitude to the management of “RUDETRANSSERVICE”, Veliky Novgorod, for providing instruments and equipment for conducting experimental determination of thermal characteristics of heat pipes with heating close to linear in time.
The author declares no conflicts of interest regarding the publication of this paper.
Seryakov, A.V. (2022) The Solving of the Inverse Thermal Conductivity Problem for Study the Short Linear Heat Pipes. Engineering, 14, 185-216. https://doi.org/10.4236/eng.2022.146018
a = λ/CP∙RHP—coefficient of HP’s thermal diffusivity;
ak = k-th calculated component of the thermal diffusivity coefficient, m2/s;
Ceff: effective heat capacity of the HP, J/K;
CHP(t)—HP’s heat capacity as a whole, J/K;
cv(t)—HP’s volume heat capacity, J/m3∙K;
cv(T)—HP’s specific isochoric heat capacity, J/kg∙K;
Cev0—heat capacity of the dry evaporator, J/K;
Cev—heat capacity of the evaporator with diethyl ether, J/K;
Cev—calculated heat capacity of the evaporator, J/K;
Cl—diethyl ether heat capacity, J/kg∙K;
Cvp—heat capacity of diethyl ether moist vapour, J/kg∙K;
Cfev—heat capacity of diethyl ether on the evaporator surface, J/kg∙K;
EH2—thermal power of the resistive heater, W;
FHP(sz)—HP’s cross-sectional area, m2;
Gl—mass flow rate of diethyl ether liquid phase in the evaporator, kg/s;
Gmix—mass flow of the moist saturated vapour over the evaporator, kg/s;
Gvp—mass flow rate of dry saturated vapour over the evaporator, kg/s;
(1 – xev) = G’/(G’ + G) —the degree of moisture vapour in the evaporator;
gk—the derivative of the heat capacity components;
Hev—enthalpy of the HP’s evaporator, W;
LHP—HP’s length, m;
Lev—the thickness of the capillary porous evaporator, m;
P ¯ —the average pressure in the vapour channel, Pa;
Pev—pressure above the evaporator, Pa;
Pfcond—pressure on the condensation surface, Pa;
Pv—the vapour pressure inside the toroidal vortex, close to the condensation pressure Pfcond;
Pr—the Prandtl number of the diethyl ether vapour;
Qev—heat input to the evaporator, W;
Qcond—heat of condensation, W;
Rc—thermal resistance of the lubricant layer, K/W;
Rev—thermal resistance of the evaporator, K/W;
Rw—thermal resistance of the metal cover, K/W;
RHP(T)—HP’s thermal resistance, K/W;
r(TB) —specific heat of evaporation of diethyl ether, kJ/kg;
TB—boiling point of the diethyl ether, K;
Tcal—the temperature of the flowing water in the vortex calorimeter;
Tcond—condensation surface temperature T, K;
T ¯ f c o n d —the average temperature of the diethyl ether film on the HP’s condensation surface, K;
T ¯ f e v —the average temperature of diethyl ether boiling microlayer on the HP’s lower cover surface, K;
T H P ( z ¯ ) —the temperature inside the HP’s vapour channel at a height of z ¯ above the evaporator, K;
Tev—stationary temperature of the grid evaporator saturated with the diethyl ether, K;
Tsc—stationary temperature of the evaporator metal mesh, K;
Tv—he vapour temperature inside the toroidal vapour vortex, K;
t ˙ —temperature growth rate, K/s;
v—vapour (and liquid) phases flow velocity inside the capillary-porous evaporator, m/s;
xev—the degree of vapour moisture in the evaporator;
Y—parameter of the additional source function, J;
z ¯ = z / L H P —dimensionless longitudinal coordinate;
uv, vv—axial and tangential components of the vortex velocity in the HP’s vapour channel, m/s;
Greek
αcal—heat transfer coefficient in the vortex flow calorimeter, W/m2∙K;
αcond—heat transfer coefficient in the HP’s condensation region, W/m2∙K;
δfcond—the average thickness of the diethyl ether film on the condensation surface, m;
δw—the thickness of the flat metal bottom cover, m;
η—Lagrange multiplier;
λeff—effective thermal conductivity, W/m·K;
λ ¯ H P —effective thermal conductivity of the HP, W/m∙K;
λl—thermal conductivity of diethyl ether, W/m∙K;
λ(t)—HP thermal conductivity coefficient, W/m∙K;
λsc—thermal conductivity coefficient of metal mesh evaporator frame, W/m∙K;
λw—thermal conductivity coefficient of metal top and bottom covers, W/m∙K;
μi—friction coefficients;
ξ—condensation coefficient defined as the ratio of the number of vapour molecules condensing over the total number of molecules which strikes the surface;
Π—the porosity of the evaporator and insert (wick) in the HP’s vapour channel;
ρHP—HP’s density, kg/m3;
ρ v p m i x —the density of moist vapour inside the vapour channel, kg/m3;
ρv—the vapour density inside the toroidal vortex, kg/m3;
ρ ¯ e v , kg/m3—average evaporator density, kg/m3;
τ—time, s;
Δτk—time step for calculations and measurements, s;
Φdf—dissipative functions;
χ—Lagrange multiplier;
ψ—the mass fractions of the HP’s components;