Analytical Determination of Nusselt Numbers for Convective Heat Transfer Coefficients in Channel Macroporous Absorbers

TL;DR

This paper presents a new analytical method to calculate heat transfer in solar absorbers with macroporous structures, improving efficiency predictions.

Contribution

A novel analytical equation for Nusselt numbers in macroporous absorbers is introduced, incorporating boundary layer effects and flow characteristics.

Findings

The new equation improves accuracy in heat transfer coefficient calculations for macroporous media.

Empirical coefficients derived using statistical methods validate the equation against numerical data.

The method enhances predictions of absorber performance in concentrated solar power plants.

Abstract

This article introduces a novel analytical equation for computing the Nusselt number within the macroporous structures of channel absorbers utilized in high-temperature solar receivers. The equation incorporates heat and mass transfer processes occurring within boundary layers as fluid flows through complex-shaped macroporous absorber channels. The importance of accounting for the length of the thermodynamic boundary layer within channel-type macroporous media when calculating heat transfer coefficients using the Nusselt equation is demonstrated. By incorporating proposed indicators of porosity and flow characteristics, this method significantly enhances the accuracy of heat transfer coefficient calculations for such media. Discrepancies observed in existing calculation relationships and experiments are attributed to the omission of certain proposed values in the Nusselt number for macroporous media. To address this, empirical coefficients for the Nusselt number are derived using statistical methods. The resulting semi-empirical equation is applied to macroporous absorbers in solar receivers. The findings enable more accurate predictions of future absorber characteristics, enhancing their efficiency. The derived equation is successfully validated against numerical data across various geometric structures of absorbers in concentrated solar power plants.

Figures

Comparative analysis of the existing dependences of the Nusselt number on Reynolds for macroporous media Re ∈ [5;150]. References to the literature: Chudnovsky A., 1962 [5], Ackermann S., 2017 [6], Wu Z.Y., 2011 [7], Wu Z.Y., 2011 [8], Yang Z., 2010 [9], Xue C., 2016 [11], Kribus A., 2014 [12], Sauerhering J, 2008 [17], Ahammad N.A., 2021 [18], Faizan, M., 2022 [24], Brendelberger S., 2012 [22].

Full text

1. Introduction

Macroporous channel media are used for various types of absorbers in high-temperature solar receivers. The optimization of such structures requires computational models of heat and mass transfer with sufficient accuracy. However, finding an analytical solution of sufficient accuracy for such structures is of scientific and practical interest.

The effective thermal conductivity $[eqn]$ and effective heat transfer coefficient for the model of heat and mass transfer in a channel-type macroporous medium are convenient values that generalize the complex processes of heat and mass transfer in a macroporous structure. The heat transfer differential equation for the model of channel homogenization in a macroporous medium is

$$[eqn]$$

where $[eqn]$ —density;

c_p_—isobaric heat capacity;m—an index referring to all macroporous material with fluids;f—an index that refers to the fluid phase;T—temperature;τ—time; $[eqn]$ —the rate of temperature change in the elementary fluid volume due to convective flows inside the macroporous structure.q_source_—thermal energy of heat sources. It may consist of radiation energy, heat of chemical reactions, heat consumption for melting and evaporation, etc.

In recent studies combining continuum and discrete models [1], the following calculation algorithm for conductive-radiative heat transfer in macroporous media is used: First, the effective thermal conductivity and the effective emissivity are approximated. Then, the Stark number is determined as an explicit function of thermal conductivity, porosity, and nominal diameter. The subsequent step involves finding the effective absorption coefficient. Finally, the effective dissipation coefficient is expressed as an explicit function of the emissivity of the solid phase, porosity, and the effective absorption coefficient. The authors recommend this model for small Stark numbers (the error was about 5%). The Rosseland model is suggested only for Stark numbers greater than 1.

For thermal evaluation, the local thermal equilibrium (LTE) method can be used, which assumes a single temperature field for both the liquid phase and the solid structure of the absorber. There are also continuum models based on the more complex Local Thermal Non-Equilibrium (LTNE) approach. LTNE modeling assumes two temperature fields for the liquid phase and the solid absorber with separate differential equations. The two phases are connected through convective heat transfer.

Quasihomogeneous approximations are one of the most common models for describing transport processes in macroporous media. In this method, the real macroporous medium is replaced by a continuous medium with effective parameters, and the equations valid in a single pore (capillary) are used as macroscopic equations. The heat transfer is modeled using the transient energy equation.

$$[eqn]$$

where $[eqn]$ —convective heat transfer coefficient W/(m^2^∙K).

s—an index that refers to the solid phase;f—an index that refers to the fluid phase;A_v_—the surface area per volume of the absorber structure.

For analytical models of heat and mass transfer, the convective heat transfer coefficient is usually calculated using the Nusselt number.

2. Methods of Contemporary Models of Convective Heat Transfer Coefficient in Macroporous Media

According to the LTNE approach, the energy equation for the macroporous material can be written as follows:

$$[eqn]$$

The energy equation for the fluid phase is defined as

$$[eqn]$$

To determine a calculation method for the convective heat transfer coefficient for macroporous channel-type media, stationary fluid flow through a single channel of an arbitrary shape is considered. Equations (3) and (4) can be rewritten as follows for a stationary process:

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ —is the effective thermal conductivity coefficient of the solid phase;

$[eqn]$ —is the effective thermal conductivity coefficient of the fluid phase.

For further consideration of the convective heat transfer coefficient, it is necessary to ensure that in the considered structures, the Knudsen number is much less than unity since it determines the conditions of heat and mass transfer in macroporous material. To meet this condition, the minimum macropore size was considered to be 0.1 mm.

In a previous publication, the authors proposed the following formulas for calculating the effective coefficient of thermal conductivity [2]:

$$[eqn]$$

where for the fluid phase $[eqn]$ is defined as follows:

$$[eqn]$$

and for the solid phase as follows:

$$[eqn]$$

where $[eqn]$ is porosity;

k_x_ and ky is the total radius vector showing the displacement of the channels relative to each other in the x and y planes, respectively (for inline array, it is equal to 1).

In this case, the coefficient $[eqn]$ is considered the displacement of elementary channels in the plane z. The z-axis is the axis perpendicular to the channel area under consideration. For repeating unit cells, $[eqn]$ . For honeycomb structures, $[eqn]$ , $[eqn]$ the offset angle of the next layer of channels.

The convective heat transfer coefficient in macroporous media is analytically found from the Nusselt number. In general form, the equation of heat transfer coefficient for the absorber channel can be written as follows:

$$[eqn]$$

where $[eqn]$ —are coefficients that consider the location of pores. In [3], it is proposed to consider c and n as a Darcy function.

Pr—Prandtl number

$$[eqn]$$

where $[eqn]$ dynamic viscosity, Pa⸱s

c_p_—specific heat, J/(kg·K)

The average method for calculating the heat transfer coefficient given in [4] is

$$[eqn]$$

where $[eqn]$ and $[eqn]$ are empirical constants depending on the Reynolds number;

d_h_—hydraulic diameter; $[eqn]$ —mass flow density.

Here is an analysis of the types of Nusselt number equations and coefficients c and n proposed in the literature. In [5], five different authors are analyzed who propose different heat transfer coefficient correlations for dispersed bulk materials. Bulk materials are considered a macroporous medium with channels formed by voids. The author concludes that for values of $[eqn]$

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ —mass flow of fluid per unit area.

For Re > 200:

$$[eqn]$$ $$[eqn]$$

In article [6], the heat transfer coefficient of gas passing through macroporous ceramics is calculated from the Nusselt number.

$$[eqn]$$

In article [7], the following models are offered:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where $[eqn]$ —specific surface area;

H = 1 is defined as the Heywood circularity factor

The author in [8] uses the following model for macroporous solar receivers of open type:

$$[eqn]$$ $$[eqn]$$

The author in [9] proposes an equation for porous media in solar thermal molten salt energy storage systems.

$$[eqn]$$

One of the seven models considered by the authors in [10] has the most adequate physical explanation and passes the test for limit transitions, namely model [11].

$$[eqn]$$

Using different models of volumetric convective heat transfer coefficients in [10] leads to a difference in efficiency of 27.7%, indicating a high sensitivity of the thermal efficiency of the absorber of the considered models to convective heat transfer coefficients. Based on a comparative analysis of data from different authors, the following formula is recommended for porous absorbers in high-temperature solar receivers [12]:

$$[eqn]$$

A formal pore diameter based on the manufacturer’s specification d = 1/PPC, a ‘passage’ diameter that includes the porosity $[eqn]$ , and an effective strut diameter we can find as follows:

$$[eqn]$$

It is worth noting that none of the proposed models considers the Mikheev correction [13] for the variability of the properties of capillary fluids $[eqn]$ . In this correction, the thermophysical properties of the fluid are calculated for the fluid temperature and wall temperature. This correction establishes the relationship between the heat transfer coefficient of the fluid and the wall temperature. Suppose that the average wall temperature can be calculated using the following formula:

$$[eqn]$$

where $[eqn]$ and $[eqn]$ are the maximum temperature and minimum temperature of the solid phase in the absorber, which lies in the area at a distance from the front of the surface of maximum radiation absorption, K;

$[eqn]$ is the relative depth of the zone of maximum absorption of solar radiation in the absorber.

Among the amendments, it should be noted that for microporous structures, the authors [14] recommend calculating the actual velocity in microchannels $[eqn]$ as the velocity is determined through the mass flow rate, considering the superficial velocity as follows:

$$[eqn]$$

and reference speed for calculations as follows:

$$[eqn]$$

where

d—hydraulic diameter, m;L—channel length, m;ρ_fw_—fluid density by a wall temperature, kg/m^3^σ_μ_—experimental value of the viscous slip coefficient = 1.1

The use of another approach, through the effective coefficient of thermal conductivity and determination of heat flow on computer models of solar receivers, can provide more accurate results, but only for a specific porous structure of the modeled object. In his dissertation, Grobbel J. [15] calculates the convective heat transfer coefficient between dispersed spherical particles as a sequence of two thermal resistances (method from [16]): contact resistance $[eqn]$ and penetration resistance $[eqn]$ , as follows:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where

A_cont_—contact area,d_p_—particle diameter,l—the modified mean free path and δ—the surface roughness of the particles.

When using the channel structure, the equation used for small-diameter pipes was also considered Equations (25.1) and (25.2). The macroporous medium is represented by staggered or corridor pipes. For them, Nu is determined from known equations, and the length of hydrodynamic and thermal stabilization in the channel is determined by

$$[eqn]$$ $$[eqn]$$

where $[eqn]$   $[eqn]$ are empirical coefficients equal to 0.07 and 0.065 for round channels.

Replacing the macroporous medium with a channel medium can be geometrically acceptable for many open-type macroporous media configurations. This method can be applied analogously to the passage of air through pipes. In the article [17], the similarity of the obtained experimental data for the Nusselt number in metal foams for the case of laminar flow inside the pipe, which gives the Nusselt numbers, is noted as Nu ≈ 3.7 (for the defining pore diameter of 0.5 mm, the value Re ≈ 50).

For air flow on top of a bundle of pipes at Re < 1000, take

$$[eqn]$$

According to the hydrodynamic theory of heat transfer for the developed turbulent continuous motion, there is a connection between heat transfer and resistance, which can be expressed as follows [13]

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ —aerodynamic resistance

Pe—Peclet number

$$[eqn]$$

E—coefficient

$$[eqn]$$

Having a set of calculation equations for the Nusselt number and several different methods for determining the convective heat transfer coefficient in a macroporous medium (numerical calculation method [18,19], analytical calculation method based on the Nusselt and Reynolds number, Nusselt and Peclet number, analytical method of heat transfer calculation as thermal resistance [15], analytical method of Fourier solution or homotropy analysis method [20], etc.), it can be concluded that an integrated approach is required to determine the convective heat transfer coefficient of macroporous media. At the same time, it is not possible to compare the methods and their contribution to the model of heat transfer in a macroporous medium of channel type for different absorber structures with different coolants in order to determine the best one for all open-type solar absorbers. This is primarily due to the lack of clear mathematical numbers for describing the macroporous medium of channel type. Numerical calculation methods are set only by our idea of heat transfer in a macroporous medium, considering only the general, averaged distribution of thermal energy over the material without considering the physical essence of the heat transfer itself and neglecting microprocesses. However, at the same time, they strictly describe a given geometric characteristic of the macroporous medium of channel type and allow for the determination of a result based on the current level of knowledge about heat transfer processes. Therefore, it is necessary to find a new computational Nusselt equation that will include indicators of heat and mass transfer in the porous medium. To research these indicators, the existing dependencies were analyzed.

To construct a graphical dependence of the Nusselt number on Re for various equations and experimental data (Figure 1), it was assumed that the characteristics of a macroporous medium are constant. The Reynolds number was calculated for the temperature at which most types of absorbers are currently tested, according to [16] $[eqn]$ K; Pr = 0.7023; $[eqn]$ W/(m∙K); and the desired temperature level, $[eqn]$ K; Pr(des) = 0.708; $[eqn]$ W/(m∙K). The fluid or void cross-section Af defines the porosity; therefore, it can be directly recalculated from the porosity and the unit cell cross-section: $[eqn]$ . The solid cross-section is then simply given as $[eqn]$ .

It was considered that the coolant velocity is 1.5 m/s, the equivalent channel diameter is 2 mm, the porosity of the absorber is 90%, and the value of Re = 92.5 was obtained. Therefore, in Figure 1, the range Re ∈ [5;150] is chosen. The data obtained using a numerical model for high-temperature solar receivers [21] (based on the numerical LTNE method) are listed as the “numerical DLR method”. “Numerical DLR method front” is a Nusselt calculation for the flow perpendicular to the porous medium (used for the surface of a porous absorber). The numerical results are compared against experimental data taken from [14,22,23,24].

Comparative analysis of the existing dependences of the Nusselt number on Reynolds for macroporous media Re ∈ [5;150]. References to the literature: Chudnovsky A., 1962 [5], Ackermann S., 2017 [6], Wu Z.Y., 2011 [7], Wu Z.Y., 2011 [8], Yang Z., 2010 [9], Xue C., 2016 [11], Kribus A., 2014 [12], Sauerhering J, 2008 [17], Ahammad N.A., 2021 [18], Faizan, M., 2022 [24], Brendelberger S., 2012 [22].

From Figure 1, it can be concluded that for macroporous absorbers of open type, it is most appropriate to use the equations given by [5,11,12,17]. For microporosity with Re < 3, the equation [11,17] can be used. However, it should be noted that the experimental data on various macroporous structures of open-type absorbers are not enough to give priority to a specific equation.

According to [25], the following equation holds as an approximate solution for a hydrodynamically developed laminar flow in a pipe:

$$[eqn]$$ $$[eqn]$$

There is a proposal to find Nu in [25] as follows:

$$[eqn]$$

However, in our case, this value did not coincide with the results obtained.

Figure 1 shows the uncertainty of experimental data and existing dependencies. Moreover, by analyzing the dependency of the Nusselt number on the Reynolds number, it can be concluded that as the geometric scale of the porosity decreases, the effect of the Reynolds number becomes smaller, i.e., the exponential coefficient of the Reynolds number is close to one (cf., e.g., Equation (10)).

An attempt is made to calculate the Nusselt numbers in the channel of a porous absorber by the classical method of determining the coefficients of dependence of the Nusselt number on the Reynolds number using data from the numerical method [21]. The channel of the porous absorber StepRec2 is shown in Figure 2. This type of absorber was chosen as an example for two reasons: according to the authors, it is the most promising absorber for modern solar towers, and it has a rather complex shape that can be divided into simple sections. Other types of absorbers were also considered by the authors and are briefly described in the following sections:

To show the impossibility of defining the coefficients c and n as constants (Equation (10)), Figure 3 shows the variation in the Nusselt number from Reynolds for the channel front of the StepRec2 absorber (absorber with complex shape channel, Figure 2) calculated using the DLR numerical method.

According to the classical method of experimental determination of the constants included in Equation (10), the dependence of Nu on Re is determined at a constant fluid Prandtl number. It has been previously proven that the correction of capillary fluid $[eqn]$ has a positive effect on the accuracy of THE power law dependence of the Nusselt number on the Reynolds numbers for the absorber channel. Therefore, Figure 2 shows an approximation of the Nusselt–Reynolds dependency according to classical theory to determine the constants of the equation.

In Figure 3, the Nusselt coefficients can be grouped by absorber section. For a conventional infinite channel, this plot is a continuous step dependence with a constant power number from 0.25 to 0.75 (Equation (10)). In Figure 3, however, a discrete function is observed, which can be divided into separate functions by channel section. In addition, the function for the first segment is also discrete. The obvious difference between the sections is the different hydraulic diameters. However, even if the Nusselt coefficients are calculated separately for each section according to the classical theory, the deviation of the obtained Nusselt numbers from the average one is too large for the first section. Therefore, an innovative approach is needed to construct a single equation that describes all sections of a given absorber simultaneously. Also, sections 4 and 5 of StepRec2 are identical in terms of effective properties such as porosity, specific surface area, etc. They have the same square channel geometry but are offset by half the cell width. The data obtained show that as the velocity increases, the convective heat transfer decreases in short sections. This can be explained by the formation of a hydrodynamic and thermal boundary layer as the fluid flows around the front surface of the macroporous channel. Therefore, it is concluded that it is necessary to recalculate the obtained data considering the thickness of the hydraulic boundary layer.

It is obvious that porosity is not a unimodal criterion for choosing one or another type of dependence for the calculation of thermal processes of convective heat transfer in a macroporous medium. Therefore, in addition to porosity, it is necessary to introduce additional indicators of the porous structure into the calculation equations.

3. Finding of Nusselt Number for Convective Heat Transfer Coefficient in Macroporous Absorbers of the High-Temperature Solar Receiver

The function for the Nusselt number in the general case of the macroporous absorber of the air receiver, which most fully describes the ratio of convective to conductive heat transfer in the absorber of the solar receiver, will be

$$[eqn]$$

Suggest that $[eqn]$ considered as an analog of a homogeneous void fraction.

$$[eqn]$$

where

L—length of the channel;l_t_—length of the thermal boundary layer (Figure 3);d_e_—hydraulic diameter.

By determining the significant coefficients, it was found that considering the ratio $[eqn]$ increases the accuracy of the approximation of the Nusselt number from the Reynolds numbers for different channel geometries. Also, it was proposed to calculate the constant c as a function of the ratio of the hydraulic diameter to the length of the hydrodynamic layer and the ratio of the section length to the length of the thermodynamic layer.

$$[eqn]$$

This proposal is an innovative solution and significantly increases the accuracy of the calculation of the heat transfer coefficient in channel macroporous structures. So, the following equation is proposed for macroporous channel structures to calculate the Nusselt number:

$$[eqn]$$

where $[eqn]$ —coefficient changes the shape of the channel;

for the first section, $[eqn]$ ;for the last section, $[eqn]$ ;

$[eqn]$   $[eqn]$   $[eqn]$ are structural constants. These coefficients are the same for the entire geometry of the complex channel shape (for all sections of the macroporous channel).

The $[eqn]$ coefficients were selected based on the results of approximation and averaging between absorber sections. The values of the averaged coefficients c and n and the method of their calculation are shown in Table 1. The statistic deviations of the invented formula are shown in Table 2.

The approximation graphs are shown below (Figure 4):

The obtained coefficients are $[eqn]$   $[eqn]$   $[eqn]$ . The statistical errors for the sections of the StepRec2 model for different section lengths are shown in the table. As can be seen, the accuracy of the universal calculation for the StepRec2 model was 18.3%.

Proposed porosity indicators $[eqn]$ are such that they will significantly improve the accuracy of the calculation of the Nusselt number and make it possible to predict its change in macroporous channel structures.

For complex geometries of macroporous channel structures, having only four total coefficients instead of two local ones for each section, as suggested by classical theory (10 coefficients for the absorber used as an example), simplifies the engineering task of finding an equation to solve the Nu number. As will be shown below, the proposed coefficients are also close to universal and can be used with some accuracy for all macroporous channel structures.

4. Result and Discussion

The proposed formula was validated using a numerical method on the DLR computer model. The computational model was built in the OpenFOAM Version 6 open-source software environment and included 4164 Nusselt number calculations for various configurations of the StepRec2 absorber, each with different input parameters [21]. Below are validation graphs (Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9) comparing calculations using the proposed formula with those using the standard engineering method outlined in HeatAtlas [16]. The graphs depict results for five different section lengths of the StepRec2 model.

In all cases, the calculation based on the proposed formula provides greater accuracy than the standard method. Furthermore, it demonstrates the potential to reduce the maximum deviation to 10.9% by recalculating coefficients for each channel section based on the hydraulic diameter.

For the first section of the model with manually adjusted coefficients from StepRec2, the maximum deviation is 5.5%. This accuracy is considered sufficient to prevent overfitting of the coefficients.

For the second section of the model with manually adjusted coefficients of StepRec2, the maximum deviation is 7.7%. This accuracy is considered sufficient compared to the model with universal coefficients.

For the third section of the model with manually adjusted coefficients from StepRec2, the maximum deviation is 17%. The values of the coefficients could not be chosen for greater accuracy.

For the fourth section of the model with manually adjusted coefficients of StepRec2, the maximum deviation is 14.5%. The values of the coefficients could not be chosen for higher accuracy.

For the fifth section of the model with manually adjusted coefficients of StepRec2, the maximal deviation is 11.2%. However, the values of these coefficients could not be further optimized for greater accuracy.

The increased accuracy of the proposed formula, as demonstrated in previous graphs, does not alter the functional dependence of the Nusselt number. Therefore, it is proposed to use the newly derived formula with constant coefficients for engineering calculations. This will facilitate the construction of an optimization model for heat and mass transfer in structured channel macroporous media, enabling the identification of a preliminary range of geometric characteristics that maximize thermal efficiency.

Further refinement of the optimal geometric characteristics of the channel-type macroporous medium for more precise calculations will be achieved through the use of numerical models. This approach is expected to significantly reduce both the time and computational resources required to identify the optimal macroporous structure of the absorber.

An additional potential application of the developed formula is in adjusting coefficients to enhance the calculation accuracy for new types of absorbers. The future objective will be to minimize the number of numerical models required for different geometries while maximizing the accuracy of the developed formula.

To provide insight into the variety of porous structures discussed in this article, Figure 10 illustrates the geometric structure of porous absorbers utilized in a solar tower.

The proposed method has also been tested for different types of macroporous absorber structures (existing solar absorbers of optimal configuration) for different input mass fluxes and different temperatures. The results are presented in Figure 11, Figure 12, Figure 13 and Figure 14. The proposed equation with the same coefficients was used for validation with different channel shapes and diameters. The validation results are presented below. The model of Hoffschmidt B. for calculating Nusselt numbers [26] was also used for the validation of HitRec and the first section of EmiRec.

In the figure, the deviation in the first section is from −6.6% to 14.3%. The deviation in the second section is from 2% to 33.7%.

The presented validation for completely different geometries and parameters of fluid confirms the possibility of using the proposed formula for high-temperature solar receivers and shows the importance of using the proposed porous structure indicators for macroporous channel structures.

Understanding the type of flow in macroporous channels is important for using standard similarity models and for understanding the heat and mass transfer processes in the interchannel macroporous medium. Most authors consider the motion in the macroporous medium of the channel type of the absorber of solar tower stations laminar due to low flow velocities [27,28,29]. Indeed, based on the analysis (Figure 3: Comparative analysis of the existing dependences of the Nusselt number on Reynolds for macroporous media), the fluid Reynolds values are too small to consider the flow as turbulent, but the study shows [30] vortices and micro convective reverse currents appear at the non-linear structure of the channel. Therefore, it is proposed to improve the formula for calculating the heat transfer coefficient according to the Nusselt number by introducing an additional variable, the dislocation vector k_y_, which describes the unevenness of the structure [2]. This factor is entered as a correction to the equivalent diameter and takes the value of one for equal channels as follows:

$$[eqn]$$

With an analytical solution to the energy transfer equation, it is possible to create an optimization model for maximizing the thermal efficiency of a macroporous absorber used in open-type receivers for solar tower power plants. The obtained equation also shows the possibility of improving the existing calculation dependences of the Nusselt number for macroporous channel-type media by including the stabilization coefficient $[eqn]$ , $[eqn]$ .

To determine the influence of input factors on heat and mass transfer in porous media, we will conduct the following analysis:

$$[eqn]$$ $$[eqn]$$

From this relationship, it can be imagined that there is a decrease in the porosity of the medium while maintaining all other geometric parameters. However, the function of the coefficient c is a family of descending curves, and therefore, the porosity is uncertain, and its optimal value depends on other factors.

Consider the factor of flow stabilization

$$[eqn]$$ $$[eqn]$$

Consider the influence of other factors included in the coefficient n

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$ $$[eqn]$$ $$[eqn]$$

for the first section, $[eqn]$ ; The porosity of the first section should be maximal with $[eqn]$ for the last section, $[eqn]$ . The porosity of the last section should be minimal with $[eqn]$

5. Conclusions

The proposed porosity indicators $[eqn]$ offer significant potential for enhancing the accuracy of Nusselt number calculations and enabling the prediction of convective heat transfer coefficients in macroporous channel structures. Moreover, the article introduces a novel analytical method for computing heat transfer coefficients using the Nusselt number. This method stands out for its incorporation of additional macroporous medium indicators, such as the stabilization coefficient $[eqn]$ and the pore shape $[eqn]$ . Additionally, it includes calculations for the lengths of the hydrodynamic and thermodynamic layers within macroporous media with channel-type pores.

The universal formula for calculating the Nu number in the absorber channel is as follows:

$$[eqn]$$

In this study, constants for the proposed dependencies for modern types of solar receiver absorbers were determined. The constants of the proposed Nusselt equation were found to be $[eqn]$   $[eqn]$   $[eqn]$ .

It was determined that stabilization of the flow in the channel has a positive effect on heat transfer, and for this purpose, the desired length of each section of the channel (without changing the shape) should be greater than 0.07∙Re∙d∙Pr. However, as the channel size increases, the degree of the Reynolds criterion decreases slightly. Therefore, the rational choice of channel length is 0.07∙Re∙d∙Pr.

Increasing the hydraulic diameter of the channel with the same overall porosity of the medium increases the heat transfer of the fluid.

Verification of the results on the DLR numerical model demonstrated satisfactory reliability and a significant improvement in accuracy compared to modern analytical models.

The root mean square error for the determined geometrical and physical parameter variations were as follows: for the StepRec2 type absorber (20,820 variations), it was 0.296627; for the StepRec type absorber (35 physical parameter variations), it was 0.3597; and for the HiTRec type absorber (30 physical parameter variations), it was 0.4184.

For complex geometries of macroporous duct structures, having only four coefficients instead of the traditional ten simplifies the engineering task of solving the Nusselt number equation. This reduction in coefficients streamlines the process without sacrificing accuracy.

The potential applications of the proposed dependence are extensive. It can be used for mathematical optimization of heat and mass transfer processes in macroporous open channel media, allowing the analytical determination of preliminary design parameters without relying solely on numerical methods. Furthermore, it can improve numerical methods for optimizing the structure of macroporous open channel media and reduce iterations in machine learning algorithms.

Overall, the results of this study provide a promising way to improve the efficiency and accuracy of thermal analysis in solar receiver absorbers, paving the way for advances in solar energy technology.