Next Article in Journal
Research on the Concept of Hydrogen Supply Chains and Power Grids Powered by Renewable Energy Sources: A Scoping Review with the Use of Text Mining
Next Article in Special Issue
Investigating the Impact of Plumbing Configuration on Energy Savings for Falling-Film Drain Water Heat Recovery Systems
Previous Article in Journal
Experimental Flow Boiling Study of R245a at High Reduced Pressures in a Large Diameter Horizontal Tube
Previous Article in Special Issue
Designing for the Environment: An Example of Multi-Criteria Analysis Used for Solar Hot Water System Selection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Aerodynamics and Complicated Heat Transfer with the Mixed Motion of Air in the Flat Duct of a High-Temperature Heat Exchanger

1
Institute of Engineering Thermophysics, National Academy of Sciences of Ukraine, 03057 Kiev, Ukraine
2
Department of Building Physics and Renewable Energy, Kielce University of Technology, al. Tysiąclecia Państwa Polskiego 7, 25-314 Kielce, Poland
*
Author to whom correspondence should be addressed.
Energies 2022, 15(3), 865; https://doi.org/10.3390/en15030865
Submission received: 27 December 2021 / Revised: 18 January 2022 / Accepted: 23 January 2022 / Published: 25 January 2022
(This article belongs to the Special Issue Heat Transfer and Heat Recovery Systems)

Abstract

:
The purpose of the research is to study the aerodynamics and heat transfer in the duct of a high-temperature recuperation system. The object of the research is a flat duct with a thickness-to-height ratio a:b = 1:10, length c = 400·a with one-sided heat input to the duct surface, complicated heat transfer and mixed air movement in the duct. The objectives of the research are to determine: (a) average temperatures on the duct surfaces; (b) air temperature distribution along the length of the duct; (c) local and average integral temperature values along the length of the duct; (d) local and average integral heat flow densities on the duct surfaces; (e) local and average integral heat exchange coefficients and Nusselt numbers on the duct surfaces; and (f) pressure distribution along the length and total pressure drop in the duct. The research method is based on conducting a mathematical numerical experiment in a stationary three-dimensional Cartesian formulation and physical modeling of processes on a prototype. It is established that the dominant heat transfer between the duct walls and the air is a complicated convective heat transfer, in which forced convection is affected by free convection. There is a 1.5–2.0-fold difference in the values of the heat flow density on the hot and non-heated surfaces of the duct. Generalizing dependences of Nusselt numbers, temperature pressures and friction resistance coefficients as functions of a dimensionless coordinate are obtained.

1. Introduction

Issues of high-temperature recuperation have an attitude to the processes of transferring of substance and energy, which are complicated by the simultaneous action of various factors. These are and different kinds and regimes of substance motion, and different mechanisms of transferring heat in the objects that are investigated. Many works are devoted to studying complicated heat transfers in compact high-temperature heat exchangers of various types [1,2,3,4]. The dominant mechanisms of transferring in them are convective and radiative heat transfer [5,6,7,8,9,10]. The feature of the work consists of the applied nature of the research aimed at studying the conditions of flow and heat transfer in an air-cooled duct, which is the basis of heat recuperation devices for high-temperature industrial plants. The formulation of the research tasks most fully describes the conditions-temperature, magnitude and direction of the heat flow density, the type of coolant and its mass flow, the geometry of the duct and the main physical and mechanical characteristics of its surfaces observed in modern high-temperature processes for the production of glass and stone fibers. The research and processing of the results were carried out using classical methods [11,12]. At the same time, for each of the duct surfaces involved in heat transfer and differing in the magnitude and direction of the supplied heat flow, Nusselt numbers and temperature pressures are determined. The results of the research were compared with the case of convective heat transfer in pipes and ducts with symmetrical heat supply to the duct surfaces, which allowed us to establish the nature and magnitude of the influence of complicated heat transfer on the characteristics of the process, in order to obtain their qualitative and quantitative indicators. The generalized dependences obtained make it possible to calculate the aerodynamic drag, temperature and energy parameters of the process, with complicated heat transfer and mixed motion necessary for the design of high-temperature heat recuperation systems in order to increase the energy efficiency of production processes. Utilization and reduction of heat losses in high-temperature industrial plants in order to increase their energy efficiency is implemented on the basis of the use of innovative heat exchange devices or heat recuperation systems. Heat transfer processes in them are carried out at high temperature heads, one-sided heat input, complicated radiation-conductive-convective heat transfer and mixed motion of air as a heat carrier.
Problem definition.
The simplified scheme of heat flow densities for a fragment of the lining is shown, which is typical for high-temperature heat recuperation systems with an air heat curtain, Figure 1.
In the summary conductive heat flow through the lining q c o n d   1 , one consists of losses with convective q c o n v   2 and radiation q r a d   2 heat flows from the internal surface of the air duct q c o n d   1 = q c o n v   2 + q r a d   2 ; the radiation heat flow q r a d   2 makes up for convective heat transfer from the surface of the duct q c o n v   3 and heat loss by conductive q c o n d   3 through thermal insulation q r a d   2 = q c o n v   3 + q c o n d   3 . The heat losses through thermal insulation q c o n d   3 consists from convective q c o n v   4 and radiation q r a d   4 heat flows from external surface of the equipment to environment. As can be seen, the heat losses to the environment comprise complicated radiation-convective heat transfer. For an integral description of all heat flows, the classical expressions presented in [13,14,15,16] are used to calculate the value of the convective and radiation heat exchange coefficients. It is assumed that each of the heat transfer mechanisms is independent, and there is no interaction between them.

2. Methodology of the Research

The complicated heat transfer in the air duct is determined by the action radiation and convective mechanisms of the transfer. Cases of such heat transfer do not have universal generalizing correlations, so they should be studied individually for each variant of technological equipment [13,17,18]. In particular, for our problem, the conditions that determine the features of complicated heat transfer in the duct are as follows:
-
the Boltzmann criterion is much larger than 1, B o ( 10 2 ; 3 · 10 3 ) , so the
-
mechanism of heat transfer between the duct surfaces and the air flow is
-
convection;
-
the high temperature of the hot surface of duct-up to 700 °C;
-
presence of a temperature difference between the duct surfaces;
-
radiation transfer of heat from the hot wall of the duct to the cold one;
-
turbulent flow regime, R e ( 4 · 10 3 ; 10 4 ) ;
-
the combined influence of free and forced convection on heat transfer in the
-
duct, R a ( 10 5 ; 6 · 10 6 ) .
Research objective. The research objective is to determine: (a) average temperatures on the duct surfaces; (b) distribution of air temperature along the duct length; (c) local and average integral temperature heads along the duct length; (d) local and average integral heat flow density on the duct surfaces; (e) local and average integral heat exchange coefficients and Nusselt numbers on the duct surfaces; and (f) pressure distribution along the length and total pressure drop in the duct.
The object of the research is a flat duct with the ratio of sides of thickness to height a:b = 1:10, length c = 400·a with one-sided heat input (boundary conditions of the 2nd kind) and complicated heat transfer with mixed air motion, Figure 2.
The subject of the research is the aerodynamics of air flow and heat transfer in the duct. The research is based on mathematical and physical modeling of processes. The mathematical numerical experiment is performed on a stationary three-dimensional Cartesian model. Physical modeling of the operation of the air-cooling system is carried out during pilot tests.

3. Physical Model of Complicated Heat Transfer and Mixed Air Motion in the Flat Duct

Features of the problem are: (a) one-sided heat input, as a result of which traditional convective heat transfer in the duct is complicated by radiation heat transfer between the hot surface (source) and other surfaces of the duct that are not heated; (b) air is considered as a diathermal medium; and (c) multidirectional action of vectors of volume forces and pressure forces. The research is conducted using a CFD package for turbulent flow. In particular, to describe radiation heat transfer in the duct, the surface-to-surface (S2S) model of radiation heat transfer in space bounded by gray diffuse surfaces is used. The main assumption in the S2S model-absorption, radiation or scattering by a gas medium-is ignored, and only surface radiation is considered in the analysis. A non-isothermal flow is modeled, turbulent flow regime is considered, the air flow is taken as a Newtonian liquid and the dependence of the thermophysical properties of air on temperature is considered.
Mathematical model of flow and complicated heat transfer in the duct with mixed motion of air. In a general view, the expressions for the S2S model are written as follows [15,16]
energy flow transmitted by other surfaces— A k · q i n , k = j = 1 N A j · q o u t , j · F j k ;
mathematical formulation of the principle of reciprocity— A j · F j k = A k · F k j ;
heat flow that is gotten on the surface k from surfaces N q i n , k = j = 1 N q o u t , j · F j k ;
factor of mutual influence F i j between the two surfaces with final sizes i and j
F i j = 1 A i A i A j cos θ i · cos θ j π · r 2 δ i j d A i d A j ;
radiation energy flow that transfers from surface k
q o u t , k = ε k · σ · T k 4 + ρ k · j = 1 N q o u t , j · F j k ,
where A k —area of the surface k ; F j k —factor which consider geometrical features position of the surfaces k and j and influence on the part of the energy that transmitted from surface j to surface k ; q i n , k —heat flow that get in the surface k from environmental; δ i j —defines the visibility d A j to d A i ( δ i j = 1 if d A j is visible from d A i or 0 otherwise); ε k —radiant emittance of the surface k ; ρ k —radiant reflectance of the surface k ; σ —Stefan-Boltzmann constant. Indexes: i n ,   o u t —indexes refer to data that get in or get out, correspondingly; i ,   j ,   k ,   N to refer to the surfaces i ,   j ,   k ,   N , correspondingly.
In convective heat transfer, the equations of continuity, momentum transfer, energy, kinetic energy dissipation and kinetic energy dissipation rate were used.
The following boundary conditions were set
x = a 2 ;   y   ϵ   [ b 2 ;   b 2 ] ; z   ϵ   [ 0 ; c ] , q 0 = λ a i r ( T x ) x = a 2 + q r a d , a 2 ;
x = a 2 ;   y   ϵ   [ b 2 ;   b 2 ] ; z   ϵ   [ 0 ; c ] , 0 = λ a i r ( T x ) x = a 2 q r a d , a 2 ;
x   ϵ [ a 2 ; a 2 ] ; y = b 2 ; z   ϵ   [ 0 ; c ]   0 = λ a i r ( T y ) y = b 2 q r a d , b 2 ;
x   ϵ [ a 2 ; a 2 ] ; y = b 2 ; z   ϵ   [ 0 ; c ] 0 = λ a i r ( T y ) y = b 2 q r a d , b 2 ;
x   ϵ [ a 2 ; a 2 ] ; y   ϵ [ b 2 ; b 2 ] ; z = 0 ; G i n l e t = G 0 ; T = T i n l e t ; u = v = 0 ;
x   ϵ [ a 2 ; a 2 ] ; y   ϵ [ b 2 ; b 2 ] ; z = c ; p = 0 ; G o u t l e t = G i n l e t ;       u = v = 0 ;   ( T z ) z = c = 0 ,
where q 0   ϵ   [ 5000 ;   6000 ] W/m2, G i n l e t   ϵ   [ 100 ; 200 ] kg/h, t i n l e t = 30 °C, c = 13.0 m—the parameters of the pilot plant; radiant emittance of duct surfaces ε = 0.93 .

4. Data Processing and Displaying Technique

Results of the research are viewed as integral characteristics, namely: average heat exchange coefficient from the duct surface in the i section α = q w ¯ ( t w ¯ t a ¯ ) ; average heat flow density at the wall of the duct q w ¯ = 1 b b 2 b 2 q w ( y ) d y ; average temperature of the wall of the duct t w ¯ = 1 b b 2 b 2 t w ( y ) d y ; average mass air temperature t a i r ¯ = f ρ w t a d f f ρ w d f ; average mass air flow rate w ¯ = f ρ w d f f ρ d f ; average integral heat exchange coefficient over the length of the duct α ¯ = 1 c 0 c α ( z ) d z ; average integral temperature head Δ t ¯ = 1 c 0 c ( t w ( z ) ¯ t a ( z ) ¯ ) d z ; average integral temperature of air t a i r ¯ = 1 c 0 c t a ( z ) ¯ d z ; dimensionless average integral temperature head Θ ¯ = Δ t ¯ t a i r ¯ = 0 c ( t w ( z ) ¯ t a ( z ) ¯ ) d z 0 c t a ( z ) ¯ d z ; Nusselt average number N u ¯ = α ¯ · X λ ; pressure losses in the duct caused by friction, Darcy-Weisbach formula, Δ p = ξ ¯ · c X · ρ · w ¯ 2 2 ; average coefficient of friction resistance ξ ¯ = 8 · τ ¯ w ρ · w ¯ 2 and Nusselt local number N u = α · X λ , where a ,   b ,   c —linear dimensions of the duct, respectively: width, height and length; q w ( y ) —local heat flow density on the wall in the i section, is determined by the numerical experiment; t a —temperature of air, is determined by the numerical experiment; t w ( y ) —local temperature of the wall in the i section, is determined by the numerical experiment; w —air velocity, is determined by the numerical experiment; x ,   y ,   z —coordinates along the axes, respectively, OX, OY, OZ; X —characteristic linear duct size; α ( z ) —average heat exchange coefficient from the duct surface; λ , ρ —the coefficient of thermal conductivity and density, respectively, depend on the air temperature; τ ¯ w —average tangential stress on the duct wall is determined by the experiment for the hot surface and the unheated surface. Indexes: a i r —value is determined for the air flow; and w —value is determined for the duct wall.
For the analysis of the obtained results of a numerical experiment and the convenience of comparing them with the known results of convective heat exchange in pipes under turbulent flow [19,20,21,22], the data are presented in a criterion form. The thermophysical properties that are part of dimensionless complexes are determined for the average mass temperature of the flow in this section. The hydraulic diameter is used as the characteristic linear size of the duct cross-section. The effect of free and forced convection on the heat transfer intensity is determined from the Rayleigh ( R a ) and Peckle ( P e ) complexes, respectively. The intensity of heat transfer in the dimensionless view is presented by Nusselt number.
Generalization of the results of the numerical research was carried out by: dimensionless local coordinate— Z = R e · P r   · X z ; average dimensionless coordinate Z ¯ = R e ¯ · P r ¯   · X c ; average Reynolds integral number along the length of the calculated region R e ¯ = 1 c 0 c R e ( z ) d z ; reduced heat transfer intensity from the hot surface ( N u ¯ 1 N u ¯ ) ; dimensionless temperature head from the hot surface Θ ¯ 1 ; dimensionless temperature head from the unheated surface Θ ¯ 2 ; average coefficient of friction resistance, where is: R e ( z ) —local Reynolds number based on the results of a numerical experiment; P r ¯ —the average Prandtl integral number along the length of the calculated region, determined by analogy with the number R e ¯ . Indexes: 1, 2—the value refers to the hot surface and the non-heated surface, respectively.
Verification of the adequacy of the mathematical model is carried out according to the test research, in which the intensity of convective heat transfer is determined for the described duct geometry under turbulent flow, double-sided heat input and boundary conditions of the 2nd kind. Radiation heat transfer in the duct is not considered. Other characteristics of the test research model correspond to the values taken for the main research. The heat flux density for double-sided heat input is chosen by the condition— q 1 = 0.5 q 0 , where q 0 is the heat flux density for double-sided heat input (main research).
The results of the test research are presented in a generalized form in the Z ¯ N u ¯ coordinate system and compared with the existing results of the research of convective heat transfer in pipes with turbulent flow at small heat flow densities, when there is still no effect on heat transfer of non-isothermicity through natural convection and there is no change in the thermophysical properties of the heat carrying agent [22]
N u ¯ = ( ξ / 8 ) R e ¯ · P r ¯ K 1 + K 2 ξ 8 ( P r ¯ 2 3 1 ) ,
where is K 1 = 1 + 900 R e ¯ , K 2 = 12.7 .
The coefficient of hydraulic friction resistance is calculated as follows:
ξ = ( 1.82 lg ( R e ¯ ) 1.64 ) 2 .
Limits of application of expression (1): R e ¯ = 4 · 10 3 ÷ 5 · 10 6 ; P r ¯ = 0.5 ÷ 5.0 [22]. The discrepancy between the experimental data and the results of calculations using expression (1) does not exceed 8%.

5. Results of Research Aerodynamics and Heat Transfer for One-Sided Heat Input

The initial data and algorithm for performing a numerical experiment are selected from real experiments of testing a pilot plant, and are shown in Table 1.
The designations of experiments are given in Table 1, and are used as captions in the figures (Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18). Figure 3, Figure 4, Figure 5 and Figure 6 show the velocity and temperature distributions in the air flow, in the vertical ZOY and horizontal ZOX central-section of the duct. The profiles are constructed for cross-sections located at a distance of 10, 100 and 200 calibers from the entry plane. As the caliber, the value of the hydraulic diameter of the duct X is taken.
For the results obtained for the mass flow rates of air of 100, 150 and 200 kg/h, an asymmetry of the temperature profiles relative to the central axes (cross-sections) of the duct symmetry is observed. The influence of free convection in the central part of the duct on vertical temperature profiles (Figure 5) is observed—the profiles are more filled than in the lower part, below the plane of symmetry. It is observed a higher temperature at the end of duct with a lower air flow rate, than with a higher air flow rate. It is observed decreasing the temperature difference between the duct surfaces as the air temperature increases. The surface temperature of the source (with one-sided heat input) is higher than the temperature of the non-heated surface along the entire length of the duct (Figure 6). Along the perimeter of the duct cross-section, its surface temperature does not have the same value, just as local heat exchange coefficients differ along the duct perimeter. Processing of the results of the research consisted in determining local, average heat exchange coefficients over the duct cross-section. At the same time, considered that the average duct surface temperatures have different values, the heat exchange coefficients are determined separately for the hot surface and the non-heated surface. In Figure 7 shows graphs of changes in the average temperatures (main axis on the left) of the hot surface and the non-heated surface and air along the length of the duct. In addition, this figure shows the temperature heads (auxiliary axis on the right) to the hot surface and the non-heated surface along the length of the duct.
The qualitative valuation, which is obtained for one-sided heat input at q 0 = 5000 W/m2, is also observed for cases of heat transfer at q 0 = 5500 and 6000 W/m2.
It is the lower the mass flow rate of air, than is the higher the overall temperature level of the process: (a) higher temperatures of the surfaces of the duct; (b) higher air temperature at the outlet of the duct; (c) smaller temperature difference between the hot surface and the non-heated surface at the outlet of the duct; and (d) smaller difference in temperature heads on the hot surface and on the non-heated surface. It should be expected that at the lower mass flow rate of air, at the end of the duct, the heat transfer rate on the hot surface of the duct and on the non-heated surface will approach a certain equal value. The different pattern of change in the value of the temperature head along the length of the duct on its hot surface and the non-heated surface indicates a qualitative difference in the heat transfer processes on these surfaces. It should also expect different patterns of changes in local heat exchange coefficients on these surfaces along the length of the duct. The distribution of the heat flow density along the length of the duct on the hot surface and the non-heated surface for different conditions of the experiment is shown in Figure 8. The results obtained have well-defined sections with a non-monotonous change in the value of the heat flow density—sections at the inlet/outlet of the duct. The non-monotonicity of the q value at the beginning of the duct- z X 10 -is due to the influence of the initial hydrodynamic and thermal regions, as well as the conditions for formulating a numerical experiment. The non-monotonicity of q at the end of the duct- 218 z X < 236 -is due only to the conditions of the numerical experiment. When processing the results of the research, empirical correlations are approximated by power dependences of the species- q ( z ) = A z B (approximation reliability R 2 > 0.96 ). When constructing approximation dependencies, the plots z X 10 and 218 z X < 236 are excluded from consideration. Further processing of the research results is carried out using approximation dependencies.
In Figure 9 shows the results of calculating local heat exchange coefficients ( α ) along the duct length ( z ), on the hot surface of the duct and on the non-heated surface of the duct. The present results demonstrate that: (a) the behavior of changes in heat exchange coefficients is different, if for the hot surface it qualitatively approaches the process of forced convection, then for the non-heated surface, such an analogy is not observed; (b) the heat exchange coefficients on the non-heated surface are greater than the heat exchange coefficients on the hot surface; and (c) as the air heats up, the difference between the heat exchange coefficients on opposite duct surfaces decreases.
Further generalization of the research results is carried out in the form of dimensionless complexes. Evaluation of the influence of gravitational and inertial forces on the intensity of heat transfer in the duct begins with determining the value of the Grashof G r and Reynolds R e numbers along the length of the duct. As noted above, the thermophysical properties that make up dimensionless complexes are determined for the average mass temperature of the flow. The behavior of the Grashof numbers along the length of the duct, on the hot surface of the duct and on the non-heated surface for different regimes of the experiment is shown in Figure 10. The behavior of the Grashof numbers remains the same for all cases provided for in the experiment plan (Table 1). By analogy with the Grashof number, dependencies are also obtained for Reynolds number. The Reynolds numbers dynamics is shown in Figure 11, and is similar for all experiments performed.
The greatest influence of gravitational forces on the heat transfer intensity is observed in one-third of the duct ( 10 z X < 100 ). At the hot surface, the effect of gravitational forces is 1.5–4.0 times higher than at the non-heated surface. The absolute values of the Grashof numbers allow us to assume about the undisputed effect of free convection on the intensity of convective heat transfer in the duct and about the expected differences in the values of heat exchange coefficients for mixed convection from their values for purely forced convection. Absolute values and the behavior of changes in Reynolds numbers indicate a developed turbulent flow regime along the entire length of the duct. The Reynolds number weakly depends on the value of the heat flux density on the hot surface. The difference in local values of the Reynolds numbers increases along the length of the duct, but does not exceed 8% of the initial value.
The dependence of the local Nusselt numbers along the length of the duct on the hot surface and on the non-heated surface is shown in Figure 12. The summary graph of the results of the research of complicated heat transfer in one-sided heat input and mixed motion of air is presented in Figure 13. The results are presented as the dependence of the local Nusselt numbers on the value of the dimensionless coordinate Z, which is determined by the expression Z = R e · P r   · X z .
The intensity of convective heat transfer on the hot surface is less than the intensity of ones on the non-heated surface. The intensity of convective heat transfer on the hot surface is more affected by the mass flow rate of air in the duct than by the value of the heat flow density on the surface. On the summary graph Figure 13, it can be seen that the local Nusselt numbers on the hot surface for different values of the q 0 value almost coincide. With sufficient accuracy for technical calculations, it can be argued that in the range of values q 0 [ 5000 ;   6000 ] W/m2, the intensities of convective heat transfer are the same. For the non-heated surface, the local Nusselt numbers depend on both the mass flow rate of air in the duct and the value of the heat flow density. In this case, the dependence from the heat flow density q 0 is more significant than for the hot surface. As the q 0 value increases, the local Nusselt number curves on the non-heated surface approach similar curves on a hot surface. As the Z coordinate decreases, the difference in the values of local Nusselt numbers decreases for different heat flow densities q 0 and different air mass flows in duct G .
Before plotting generalizing dependences of the average Nusselt numbers ( N u ¯ ) as functions on the dimensionless coordinate ( R e ¯ · P r ¯   · X L ), it is verificated the mathematical model and calculation methods used on the traditional problem of convective heat transfer in the duct with symmetric two-sided heat input.
The research of aerodynamics and heat transfer in symmetric two-sided heat input and turbulent flow regime in the duct (Table 2) has the following goals:
-
verification of the mathematical model and computational methods;
-
comparison of the results of the research with existing solutions for convective heat transfer in pipes and ducts;
-
comparison and generalization of the results of the research of flow and heat transfer in the duct under turbulent flow regime and symmetric two-sided heat input with the results of the research of one-sided heat input, complicated heat transfer and mixed motion.
The following boundary conditions were set
x = a 2 ; y ϵ [ b 2 ; b 2 ] ; z   ϵ   [ 0 ; c ] 0 = λ a i r ( T x ) x = a 2 q x = a 2 ; q x = a 2 = q 0 / 2 ;
x = a 2 ; y ϵ [ b 2 ; b 2 ] ; z   ϵ   [ 0 ; c ] 0 = λ a i r ( T x ) x = a 2 q x = a 2 ; q x = a 2 = q 0 / 2 ;
x   ϵ [ a 2 ; a 2 ] ; y = b 2 ; z   [ 0 ; c ] 0 = λ a i r ( T y ) y = b 2 q y = b 2 ; q y = b 2 = 0 ;
x   ϵ [ a 2 ; a 2 ] ; y = b 2 ; z   ϵ   [ 0 ; c ] 0 = λ a i r ( T y ) y = b 2 q y = b 2 ; q y = b 2 = 0 ;
x   ϵ [ a 2 ; a 2 ] ; y ϵ [ b 2 ; b 2 ] ; z = 0 ; G i n l e t = G 0 ;   T = T i n l e t ;   u = v = 0 ;
x   ϵ [ a 2 ; a 2 ] ; y ϵ [ b 2 ; b 2 ] ; z = c ; p = 0 ;   G o u t l e t = G i n l e t = G 0 ;   u = v = 0 ; ( T z ) z = c = 0 .
where is q 0 = 5000 W/m2; G i n l e t ϵ [ 100 ;   200 ] kg/h; t i n l e t = 30 °C; c = 13.0 m.
The results of the test computational experiment are presented in a dimensionless form. They are compared with the results of the computational experiment on the complicated heat transfer and mixed motion of air (Figure 14, Figure 15 and Figure 16). Figure 14 shows local Nusselt numbers on the hot surface, the non-heated surface and surfaces with symmetric two-sided heat input, respectively. Figure 15 shows the reduced Nusselt numbers on the hot surface and the non-heated surface, respectively. Figure 16 shows dimensionless average integral temperature heads for the three cases described above. Figure 14 and Figure 16 obtained the qualitatively identical situation—the objective function under symmetric two-sided heat input occupies the intermediate position between the values of the corresponding objective functions on the hot surface and the surface that is non-heated.
Further final processing of the research results consisted in plotting the generalizing dependence of the average Nusselt numbers ( N u ¯ ) on the value of the dimensionless coordinate ( R e ¯ · P r ¯   · X L ). On the graph, Figure 14 presents the results of researches of complicated heat transfer with the mixed motion of air on the hot surface (curve 2) and the non-heated surface (curve 3); the results of the test research and the results of convective heat transfer in pipes and ducts under turbulent flow regime, published in [22] (curve 1).
The average Nusselt numbers ( N u ¯ 1 ) and ( N u ¯ 2 ) are calculated using the known dependence for the average Nusselt number ( N u ¯ ), which describes convective heat transfer in pipes and ducts under turbulent flow regime, but using correction factors
( N u ¯ 1 N u ¯ ) = 5 3.9 · ( R e ¯ · P r ¯   · X L ) 1 8 ,
( N u ¯ 2 N u ¯ ) = 0.99 · e [ 4.26 · 10 4 · ( R e ¯ · P r ¯   · X L ) 1.95 ] ,
where is N u ¯ 1 —average Nusselt number on the hot surface, curve 2; N u ¯ 2 —average Nusselt number on the non-heated surface, curve 3; N u ¯ —average Nusselt number on the symmetrical two-sided heat input, curve 1, Formula (1). The accuracy approximation is (R2 ≥ 0.98).
The change in the dimensionless average integral temperature head along the length of the dimensionless coordinate for the hot surface (curve 2) and the non-heated surface (curve 3) is shown in Figure 16. For comparison, the same graph shows the dependence of the average integral temperature head for the case of symmetric two-sided heated input (curve 1).
Θ 1 = 1 19 · ( R e ¯ · P r ¯   · X L ) 1 1.8 ,
the accuracy approximation is (R2 > 0.99).
Θ 2 = 1 5.28 · e   [ 1 55.87 · ( R e ¯ · P r ¯   · X L ) ] ,
the accuracy approximation is (R2 = 0.90).
The distribution of tangential stresses on the walls along length of the duct, for different experimental conditions ( q 0 = 5000 W/m2, G = 100 ;   150 ;   200 kg/h), is shown in Figure 17.
The distribution plotting of the value τ w obtained for the value of the heat flux density q 0 = 5000 W/m2 remains relevant for another values q 0 . It can be seen that there is the minimal difference in the values of τ w obtained for different duct surfaces. The observed difference may not be considered, and the values of tangential stresses on the walls may be considered the same.
The results of the research are presented in the form of the generalizing dependence of the average coefficient of friction resistance on the value of the dimensionless coordinate (Figure 18). The generalizing graphical dependence is satisfactorily approximated by the power dependence (6), the accuracy approximation is ( R 2 > 0.99 ).
ξ ¯ = ( 3.633 R e ¯ 3.14 ) 0.1 .

6. Conclusions

It is established that in the considered range of values of heat flow density and mass flow rate of air, the dominant heat transfer between the walls and air is complicated convective heat transfer, in which forced convection ( R e > 4000 ) is affected by free convection ( 10 3 < G r < 10 7 ).
There is the difference in the intensity of convective heat transfer on the hot surface and on the non-heated surface. The difference in the density of heat flux transmitted by convection to air on the hot surface, and the non-heated surface is proportional to the expression ( ( N u ¯ 1 · Θ ¯ 1 ) ( N u ¯ 2 · Θ ¯ 2 ) ) and is 1.5–2.0.
Generalizing dependences are obtained that allow calculating the average Nusselt numbers and average temperature heads on the duct surfaces. Limits for applying dependencies (2)–(5): R e = 4 · 10 3 ÷ 10 4 ; P r = 0.68 ÷ 0.70 ; 10 3 < G r < 10 7 .
The generalizing dependence of the value of the coefficient of friction resistance as a function of the value of the dimensionless coordinate (6) is obtained. The obtained values of the resistance value allow us to assume that in the hardware design of high-temperature heat exchangers and heat recuperation systems, the value of pressure losses to overcome the resistance of the boundary layer will not exceed 10% of the total pressure losses in the system.
The research results became the basis for the creation of the engineering methodology for calculating the heat engineering indicators of high-temperature heat exchangers and heat recuperation systems, which are widely used in the production of basalt and glass fibers at the enterprises “Chernivtsi plant of thermal insulation materials” (Chernivtsi, Ukraine) and “Innovation Basalt Technology Sp. z.o.o.”(Gdansk, Poland). Their use in heat technologies allows you to return part of the heat with heated air to the technological cycle, provide controlled cooling of the equipment lining and obtain secondary energy resources in a form that is convenient for further use for one’s own technological needs or sanitary needs of production.

Author Contributions

Conceptualization, B.B. and A.P.; methodology, V.K.; software, A.T.; validation, B.B., V.K. and A.P.; formal analysis, B.B.; data curation, A.P.; writing—original draft preparation, B.B.; writing—review and editing, A.P.; visualization, A.T.; supervision, A.T.; project administration B.B.; funding acquisition, A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

References

  1. Govern, J.C.; Herman, C.V.; Nagle, D.C. High Temperature Heat Exchangers (HTHX) for Nuclear Applications. In Proceedings of the ASME 2006 International Mechanical Engineering Congress and Exposition, Chicago, IL, USA, 5–10 November 2006; Volume 1, pp. 713–718. [Google Scholar]
  2. Alkasassbeh, M.; Omar, Z.; Mebarek-Oudina, F.; Raza, J.; Chamkha, A. Heat transfer study of convective fin with temperature-dependent internal heat generation by hybrid block method. Heat Transfer-Asian Res. 2019, 48, 1225–1244. [Google Scholar] [CrossRef]
  3. Shah, P.K. Heat Exchangers for Fuel Cell Systems, in Compact Heat Exchangers and Enhancement Technology for the Process Industries; Begell House: New York, NY, USA, 2003; pp. 205–215. [Google Scholar]
  4. Marzougui, S.; Mebarek-Oudina, F.; Magherbi, M.; Mchirgui, A. Entropy generation and heat transport of Cu–water nanoliquid in porous lid-driven cavity through magnetic field. Int. J. Numer. Methods Heat Fluid Flow 2021. ahead-of-print. [Google Scholar] [CrossRef]
  5. Djebali, R.; Mebarek-Oudina, F.; Rajashekhar, C. Similarity solution analysis of dynamic and thermal boundary layers: Further formulation along a vertical flat plate. Phys. Scr. 2021, 96, 085206. [Google Scholar] [CrossRef]
  6. Kumada, M. A Study on the High Performance Ceramic Heat Exchanger for Ultra High Temperatures. Heat Transfer Enhanc. Heat Exch. 1999, IV, 301–324. [Google Scholar] [CrossRef]
  7. McDonald, C.F. Recuperator considerations for future higher efficiency microturbines. Appl. Therm. Eng. 2003, 23, 1463–1487. [Google Scholar] [CrossRef]
  8. Rafidi, N.; Blasiak, W. Thermal performance analysis on a two composite material honeycomb heat regenerators used for HiTAC burners. Appl. Therm. Eng. 2005, 25, 2966–2982. [Google Scholar] [CrossRef]
  9. Tsuji, H.; Gupta, A.; Hasegawa, T.; Katsuki, M.; Kishimoto, K.; Morita, M. High Temperature Air Combustion; from Energy Conservation to Pollution Reduction; CRC Press LLC: New York, NY, USA, 2003. [Google Scholar] [CrossRef]
  10. Dadheech, P.K.; Agrawal, P.; Mebarek-Oudina, F.; Abu-Hamdeh, N.H.; Sharma, A. Comparative Heat Transfer Analysis of MoS2/C2H6O2 and SiO2-MoS2/C2H6O2 Nanofluids with Natural Convection and Inclined Magnetic Field. J. Nanofluids 2020, 9, 161–167. [Google Scholar] [CrossRef]
  11. C Ranganayakulu, K.N.S. Compact Heat Exchangers: Analysis, Design and Optimization Using FEM and CFD Approach; ASME Press: New York, NY, USA, 2018; 546p, ISBN 9781119424185. [Google Scholar] [CrossRef]
  12. Wong, H. Basic Formulas and Data on Heat Transfer for Engineers; Addison-Wesley Longman Ltd.: New York, NY, USA, 1977; 236p, ISBN-10:0582460506. [Google Scholar]
  13. Kakaç, S.; Liu, H.; Pramuanjaroenkij, A. Heat Exchangers: Selection, Rating, and Thermal Design, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2002. [Google Scholar] [CrossRef]
  14. Minkowycz, W.J.; Sparrow, E.M.; Abraham, J.P.; Gorman, J.M. (Eds.) Advances in Numerical Heat Transfer: Numerical Simulation of Heat Exchangers, 1st ed.; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar] [CrossRef]
  15. Howell, J.R.; Mengüç, M.P.; Daun, K.; Siegel, R. Thermal Radiation Heat Transfer, 7th ed.; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar] [CrossRef]
  16. Basok, B.; Davydenko, B.; Pavlenko, A. Numerical Network Modeling of Heat and Moisture Transfer through Capillary-Porous Building Materials. Materials 2021, 14, 1819. [Google Scholar] [CrossRef] [PubMed]
  17. Asli, K.H.; Aliyev, S.A.O. Applied Research in Hydraulics and Heat Flow, 1st ed.; Apple Academic Press: New York, NY, USA, 2014. [Google Scholar] [CrossRef]
  18. Pavlenko, A.M.; Koshlak, H. Application of Thermal and Cavitation Effects for Heat and Mass Transfer Process Intensification in Multicomponent Liquid Media. Energies 2021, 14, 7996. [Google Scholar] [CrossRef]
  19. Kremnev, V.; Basok, B.; Timoshchenko, A.; Tymchyshyn, S. Energy Saving in Construction by Wide Application of High-Quality Insulation Based on Basalt Fibers. J. Mod. Phys. 2018, 09, 1724–1734. [Google Scholar] [CrossRef] [Green Version]
  20. Kremnev, V.; Basok, B.; Davydenko, B.; Timoshchenko, A.; Timoshchenko, A. Flow and Heat Transfer of Basalt Melt in the Feeder of the Smelter Furnace. J. Appl. Math. Phys. 2019, 07, 2555–2563. [Google Scholar] [CrossRef] [Green Version]
  21. Oleinik, O.A.; Samokhin, V.N. Mathematical Models in Boundary Layer Theory, 1st ed.; Chapman & Hall/CRC: London, UK, 1999. [Google Scholar] [CrossRef]
  22. Merkin, J.H. Mixed convection in a Falkner–Skan system. J. Eng. Math. 2016, 100, 167–185. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Scheme of heat transfer: 1—internal space of the technological equipment; 2—air duct for recuperation of the heat loss, type the heat curtain; 3—lining; 4—thermal insulation of the lining to reduce dissipative heat loss; 5—environment; q c o n d   1 , q c o n d   3 —conductive heat flow; q c o n v   2 , q c o n v   3 , q c o n v   4 —convective heat flow; q r a d   2 , q r a d   4 —radiation heat flow.
Figure 1. Scheme of heat transfer: 1—internal space of the technological equipment; 2—air duct for recuperation of the heat loss, type the heat curtain; 3—lining; 4—thermal insulation of the lining to reduce dissipative heat loss; 5—environment; q c o n d   1 , q c o n d   3 —conductive heat flow; q c o n v   2 , q c o n v   3 , q c o n v   4 —convective heat flow; q r a d   2 , q r a d   4 —radiation heat flow.
Energies 15 00865 g001
Figure 2. Scheme of the modelling area.
Figure 2. Scheme of the modelling area.
Energies 15 00865 g002
Figure 3. Air velocity profiles in the vertical central-section ZOY at 100 kg/h (y [m], u [m/s]), curves from left to right, respectively: 10×; 100×; 200×.
Figure 3. Air velocity profiles in the vertical central-section ZOY at 100 kg/h (y [m], u [m/s]), curves from left to right, respectively: 10×; 100×; 200×.
Energies 15 00865 g003
Figure 4. Air velocity profiles in the horizontal central-section of ZOX at 100 kg/h (x [m], u [m/s]), bottom-up curves, respectively: 10×; 100×; 200×.
Figure 4. Air velocity profiles in the horizontal central-section of ZOX at 100 kg/h (x [m], u [m/s]), bottom-up curves, respectively: 10×; 100×; 200×.
Energies 15 00865 g004
Figure 5. Air temperature profiles in the vertical central-section ZOY, 100 kg/h (t [°C], y [m]), curves from left to right, respectively: 10×; 100×; 200×.
Figure 5. Air temperature profiles in the vertical central-section ZOY, 100 kg/h (t [°C], y [m]), curves from left to right, respectively: 10×; 100×; 200×.
Energies 15 00865 g005
Figure 6. Air temperature profiles in the horizontal central-section ZOX, 100 kg/h (x [m], t [°C]): bottom-up curves, respectively: 10×; 100×; 200×.
Figure 6. Air temperature profiles in the horizontal central-section ZOX, 100 kg/h (x [m], t [°C]): bottom-up curves, respectively: 10×; 100×; 200×.
Energies 15 00865 g006
Figure 7. Change of the temperature t [˚C] and the temperature head Δt [°C] by the length of the duct z [m] ( G = 100 kg/h, q 0 = 5000 W/m2): 100-1—the average temperature of the hot surface; 100-2—the average temperature of the non-heated surface; 100-3—the average temperature of air; 100-4—the temperature head to the hot surface; and 100-5—the temperature head to the non-heated surface.
Figure 7. Change of the temperature t [˚C] and the temperature head Δt [°C] by the length of the duct z [m] ( G = 100 kg/h, q 0 = 5000 W/m2): 100-1—the average temperature of the hot surface; 100-2—the average temperature of the non-heated surface; 100-3—the average temperature of air; 100-4—the temperature head to the hot surface; and 100-5—the temperature head to the non-heated surface.
Energies 15 00865 g007
Figure 8. Distribution of the heat flow density q [W/m2] by the length of the duct z [m]: q 0 C o n s t ;   G V a r i a .
Figure 8. Distribution of the heat flow density q [W/m2] by the length of the duct z [m]: q 0 C o n s t ;   G V a r i a .
Energies 15 00865 g008
Figure 9. Changing the local heat exchange coefficients α [W/(m2·°C)] along the length of the duct z [m] q 0 = 5000 Bт/м2, G V a r i a .
Figure 9. Changing the local heat exchange coefficients α [W/(m2·°C)] along the length of the duct z [m] q 0 = 5000 Bт/м2, G V a r i a .
Energies 15 00865 g009
Figure 10. Changing the Grashof numbers along the length of the duct z [m] q 0 V a r i a ;   G C o n s t .
Figure 10. Changing the Grashof numbers along the length of the duct z [m] q 0 V a r i a ;   G C o n s t .
Energies 15 00865 g010
Figure 11. Changing the Reynolds numbers along the length of the duct z [m] q 0 V a r i a ;   G C o n s t .
Figure 11. Changing the Reynolds numbers along the length of the duct z [m] q 0 V a r i a ;   G C o n s t .
Energies 15 00865 g011
Figure 12. Changing the Nusselt numbers along the length of the duct z [m] q 0 C o n s t ;   G V a r i a .
Figure 12. Changing the Nusselt numbers along the length of the duct z [m] q 0 C o n s t ;   G V a r i a .
Energies 15 00865 g012
Figure 13. Dependence of the local Nusselt numbers on the dimensionless coordinate.
Figure 13. Dependence of the local Nusselt numbers on the dimensionless coordinate.
Energies 15 00865 g013
Figure 14. Generalizing dependences of average Nusselt numbers on the dimensionless coordinate. Curve 1—symmetrical two-sided heat input, curve 2—one-sided heat input, the hot surface; curve 3—one-sided heat input, the non-heated surface.
Figure 14. Generalizing dependences of average Nusselt numbers on the dimensionless coordinate. Curve 1—symmetrical two-sided heat input, curve 2—one-sided heat input, the hot surface; curve 3—one-sided heat input, the non-heated surface.
Energies 15 00865 g014
Figure 15. Dependence of relative Nusselt numbers on the dimensionless coordinate. Curve 1—one–sided heat input, hot surface; curve 2—one-sided heat input, non-heated surface.
Figure 15. Dependence of relative Nusselt numbers on the dimensionless coordinate. Curve 1—one–sided heat input, hot surface; curve 2—one-sided heat input, non-heated surface.
Energies 15 00865 g015
Figure 16. Generalizing dependence of dimensionless average integral temperature heads on the dimensionless coordinate. The designation of curves is similar to Figure 15.
Figure 16. Generalizing dependence of dimensionless average integral temperature heads on the dimensionless coordinate. The designation of curves is similar to Figure 15.
Energies 15 00865 g016
Figure 17. Distribution of tangential stresses τ w [Pa] on the duct walls along its length z [m] ( q 0 C o n s t ;   G V a r i a ).
Figure 17. Distribution of tangential stresses τ w [Pa] on the duct walls along its length z [m] ( q 0 C o n s t ;   G V a r i a ).
Energies 15 00865 g017
Figure 18. Generalized dependence of average coefficient of resistance of friction ξ ¯ on the dimensionless coordinate R e ¯ .
Figure 18. Generalized dependence of average coefficient of resistance of friction ξ ¯ on the dimensionless coordinate R e ¯ .
Energies 15 00865 g018
Table 1. Plan for performing the numerical experiment with one-sided heat input.
Table 1. Plan for performing the numerical experiment with one-sided heat input.
Heat Flow Density, W/m2Mass Flow Rate of Air, kg/hNotes
100150200
50005000_100_15000_150_15000_200_1Hot surface
50005000_100_25000_150_25000_200_2Non-heated surface
55005500_100_15500_150_15500_200_1Hot surface
55005500_100_25500_150_25500_200_2Non-heated surface
60006000_100_16000_150_16000_200_1Hot surface
60006000_100_26000_150_26000_200_2Non-heated surface
Table 2. Plan for performing the numerical experiment with two-sided heat input.
Table 2. Plan for performing the numerical experiment with two-sided heat input.
Heat Flow Density, W/m2Mass Flow Rate of Air, kg/hNotes
100150200
Two-Sided Heat Input
25002500_1002500_1502500_200Symmetrical heating
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Basok, B.; Kremnev, V.; Pavlenko, A.; Timoshchenko, A. Aerodynamics and Complicated Heat Transfer with the Mixed Motion of Air in the Flat Duct of a High-Temperature Heat Exchanger. Energies 2022, 15, 865. https://doi.org/10.3390/en15030865

AMA Style

Basok B, Kremnev V, Pavlenko A, Timoshchenko A. Aerodynamics and Complicated Heat Transfer with the Mixed Motion of Air in the Flat Duct of a High-Temperature Heat Exchanger. Energies. 2022; 15(3):865. https://doi.org/10.3390/en15030865

Chicago/Turabian Style

Basok, Borys, Vyacheslav Kremnev, Anatoliy Pavlenko, and Andriy Timoshchenko. 2022. "Aerodynamics and Complicated Heat Transfer with the Mixed Motion of Air in the Flat Duct of a High-Temperature Heat Exchanger" Energies 15, no. 3: 865. https://doi.org/10.3390/en15030865

APA Style

Basok, B., Kremnev, V., Pavlenko, A., & Timoshchenko, A. (2022). Aerodynamics and Complicated Heat Transfer with the Mixed Motion of Air in the Flat Duct of a High-Temperature Heat Exchanger. Energies, 15(3), 865. https://doi.org/10.3390/en15030865

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop