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

: 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 ﬂat 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 ﬂow densities on the duct surfaces; (e) local and average integral heat exchange coefﬁcients 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 ﬂow density on the hot and non-heated surfaces of the duct. Generalizing dependences of Nusselt numbers, temperature pressures and friction resistance coefﬁcients as functions of a dimensionless coordinate are obtained.


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 conditionstemperature, 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.
Energies 2022, , x FOR PEER REVIEW 2 of 18 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 , one consists of losses with convective and radiation heat flows from the internal surface of the air duct = + ; the radiation heat flow makes up for convective heat transfer from the surface of the duct and heat loss by conductive through thermal insulation = + . The heat losses through thermal insulation consists from convective and radiation 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.

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 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 cond 1 , q cond 3 -conductive heat flow; q conv 2 , q conv 3 , q conv 4 -convective heat flow; q rad 2 , q rad 4 -radiation heat flow.
In the summary conductive heat flow through the lining q cond 1 , one consists of losses with convective q conv 2 and radiation q rad 2 heat flows from the internal surface of the air duct q cond 1 = q conv 2 + q rad 2 ; the radiation heat flow q rad 2 makes up for convective heat transfer from the surface of the duct q conv 3 and heat loss by conductive q cond 3 through thermal insulation q rad 2 = q conv 3 + q cond 3 . The heat losses through thermal insulation q cond 3 consists from convective q conv 4 and radiation q rad 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.

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, Bo ∈ 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, Re ∈ 4·10 3 ; 10 4 ; -the combined influence of free and forced convection on heat transfer in the duct, Ra ∈ 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. conditions that determine the features of complicated heat transfer in the duct are as follows: the Boltzmann criterion is much larger than 1, ∈ (10 ; 3 • 10 ), 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, ∈ (4 • 10 ; 10 ); the combined influence of free and forced convection on heat transfer in the duct, ∈ (10 ; 6 • 10 ).
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.

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) 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.

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 mediumis 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 in,k = ∑ N j=1 A j ·q out,j ·F jk ; mathematical formulation of the principle of reciprocity-A j ·F jk = A k ·F kj ; heat flow that is gotten on the surface k from surfaces N − q in,k = ∑ N j=1 q out,j ·F jk ; factor of mutual influence F ij between the two surfaces with final sizes i and j- radiation energy flow that transfers from surface k where A k -area of the surface k; F jk -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 in,k -heat flow that get in the surface k from environmental; δ ij -defines the visibility dA j to dA i (δ ij = 1 if dA j is visible from dA i or 0 otherwise); ε k -radiant emittance of the surface k; ρ k -radiant reflectance of the surface k; σ-Stefan-Boltzmann constant. Indexes: in, out-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 .0 m-the parameters of the pilot plant; radiant emittance of duct surfaces ε = 0.93.

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 (tw−ta) ; Nusselt average number Nu = α· 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 Nu = α· 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 atemperature 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: air-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 (Ra) and Peckle (Pe) 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 = Re·Pr · X z ; average dimensionless coordinate Z = Re·Pr · X c ; average Reynolds integral number along the length of the calculated region Re = 1 c c 0 Re(z)dz; reduced heat transfer intensity from the hot surface (Nu1/Nu); dimensionless temperature head from the hot surface Θ 1 ; dimensionless temperature head from the unheated surface Θ 2 ; average coefficient of friction resistance, where is: Re(z)-local Reynolds number based on the results of a numerical experiment; Pr-the average Prandtl integral number along the length of the calculated region, determined by analogy with the number Re. 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 conditionq 1 = 0.5q 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 − Nu 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] Nu = (ξ/8)Re·Pr where is K 1 = 1 + 900 Re , K 2 = 12.7. The coefficient of hydraulic friction resistance is calculated as follows: Limits of application of expression (1): Re = 4·10 3 ÷ 5·10 6 ; Pr = 0.5 ÷ 5.0 [22]. The discrepancy between the experimental data and the results of calculations using expression (1) does not exceed 8%.

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          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 = 5000 W/m 2 , is also observed for cases of heat transfer at = 5500 and 6000 W/m 2 . 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 tem-   The qualitative valuation, which is obtained for one-sided heat input at = 5000 W/m 2 , is also observed for cases of heat transfer at = 5500 and 6000 W/m 2 . 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 tem- The qualitative valuation, which is obtained for one-sided heat input at q 0 = 5000 W/m 2 , is also observed for cases of heat transfer at q 0 = 5500 and 6000 W/m 2 .
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) = Az 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. 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 value at the beginning of the duct-≤ 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 at the end of the duct-218 ≤ < 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-( ) = (approximation reliability > 0.96). When constructing approximation dependencies, the plots ≤ 10 and 218 ≤ < 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 ( ), 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 and Reynolds 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.  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 Gr and Reynolds Re 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. 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 and Reynolds 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 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 = • • . Figure 13. Dependence of the local Nusselt numbers on the dimensionless coordinate. 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 = Re·Pr · X z .
Energies 2022, , x FOR PEER REVIEW 12 of 18 The greatest influence of gravitational forces on the heat transfer intensity is observed in one-third of the duct (10 ≤ < 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 = • • . Figure 13. Dependence of the local Nusselt numbers on the dimensionless coordinate. Figure 13. Dependence of the local Nusselt numbers on the dimensionless coordinate. 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 value almost coincide. With sufficient accuracy for technical calculations, it can be argued that in the range of values ∈ [5000; 6000] W/m 2 , 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 is more significant than for the hot surface. As the value increases, the local Nusselt number curves on the non-heated surface approach similar curves on a hot surface. As the coordinate decreases, the difference in the values of local Nusselt numbers decreases for different heat flow densities and different air mass flows in duct .
Before plotting generalizing dependences of the average Nusselt numbers ( ) as functions on the dimensionless coordinate ( • • ), 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)  The following boundary conditions were set 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/m 2 , 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 (Nu) as functions on the dimensionless coordinate (Re·Pr · 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)  The following boundary conditions were set where is q 0 = 5000 W/m 2 ; G inlet [100; 200] kg/h; t inlet = 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 (Figures 14-16). Figure 14 shows local Nusselt numbers on the hot surface, the non-heated surface and surfaces with symmetric twosided 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. Figures 14 and 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.
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. Figures 14 and 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 (Nu) on the value of the dimensionless coordinate (Re·Pr · 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).
Further final processing of the research results consisted in plotting the generalizing dependence of the average Nusselt numbers ( ) on the value of the dimensionless coordinate ( • • ). 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 nonheated 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 ( ) and ( ) are calculated using the known dependence for the average Nusselt number ( ), which describes convective heat transfer in pipes and ducts under turbulent flow regime, but using correction factors where is -average Nusselt number on the hot surface, curve 2; -average Nusselt number on the non-heated surface, curve 3; -average Nusselt number on the symmetrical two-sided heat input, curve 1, Formula (1). The accuracy approximation is (R 2 ≥ 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). The average Nusselt numbers (Nu 1 ) and (Nu 2 ) are calculated using the known dependence for the average Nusselt number (Nu), which describes convective heat transfer in pipes and ducts under turbulent flow regime, but using correction factors (Nu1/Nu) = 5 3.9 · Re·Pr · X L where is Nu 1 -average Nusselt number on the hot surface, curve 2; Nu 2 -average Nusselt number on the non-heated surface, curve 3; Nu-average Nusselt number on the symmetrical two-sided heat input, curve 1, Formula (1). The accuracy approximation is (R 2 ≥ 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).
The distribution of tangential stresses on the walls along length of the duct, for different experimental conditions (q 0 = 5000 W/m 2 , 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/m 2 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).

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 (Re > 4000) is affected by free convection (10 3 < Gr < 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 ( Nu 1 ·Θ 1 / Nu 2 ·Θ 2 ) and is 1.5-2.0.
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 hightemperature 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.