Abstract
Solar air heating thermal systems have found extensive utilization in a broad array of industrial and residential settings, playing a pivotal role in the conversion and reclamation of solar energy. Implementing repeated artificial roughness in the surfaces has the potential to augment thermal performance in solar air heaters (SAHs). This study presents a numerical investigation of SAHs with artificial rough surfaces, consisting of polygonal-shaped ribs and grooves located at different places inside the rectangular duct, that improve thermal efficiency. ANSYS Fluent software was employed to simulate the SAH with different relative pitch distances of p = 10 mm and 20 mm and relative rib heights e/d = 0.09–0.045. The working fluid air flows at different Reynolds numbers (Re), ranging from 3,800 to 18,000. Nusselt number (Nu), friction factor (f), and Thermal Hydraulic Performance (THP) are parameters to evaluate the performance of the SAH. The renormalized group (RNG) k-ϵ turbulent model was implemented in this simulation. The study outcomes indicate that increasing the rib height improves the heat transfer rate and nonetheless increases pressure drop while increasing the pitch distance. The higher Nusselt number (Nu) is 3.762 attained at p = 10 mm and 3.420 at p = 20 mm in the center-positioned rib at Re 3,800. The lower friction factor (ƒ) obtained in p = 20 mm is 1.681 and 0.785 in p = 10 mm in the staggered positioned rib at higher Re 15,000. The optimal THP was achieved at 2.813 at a staggered rib height at a pitch distance of p = 10 mm at Re 8,000. The study’s findings suggest that the incorporation of artificial rough surfaces has the potential to enhance the THP of an SAH.
1 Introduction
A solar air heater (SAH) is a device that harnesses the abundant solar energy available worldwide for various applications, including space heating, textile, and chemical industries. However, the conventional SAH suffers from lower thermal performances due to poor thermal convection between the working fluid air and the surface of an absorber plate. In recent decades, researchers have explored the use of artificial roughness in the absorber plate surfaces in augmenting heat transfer in the SAH. The implementation of rough surfaces in absorber plates is achieved by introducing ribs or protrusions of different shapes and sizes. Some of the significant outcomes of these numerical and experimental investigations include an increase in heat transfer rates, a decrease in pressure drops, and improved thermal performance of the SAH (; ). Various rib and groove shapes have been selected in the investigation to attain better heat transfer rates, including rectangular (), circular (), triangular (), trapezoidal, and serrated ribs () at Reynolds numbers ranging from Re 3,000 to Re 18,000. Among them, rectangular and circular rib shapes are the most commonly investigated due to their ease of manufacturing and the availability of experimental data. The absorber plate efficiency of the SAH is enhanced using various rib shapes, such as a transverse wire rib producing a maximum THP of 1.85 at Reynolds number ranging from Re 3,800 to Re 18,000 (), an arc-shaped wire rib producing a THP of 1.70 at Reynolds number ranging from Re 3,800 to Re 18,000 (), an inverted U-shape rib producing a THP of 1.82 at Reynolds number ranging from Re 3,800 to Re 18,000 (), a dimple-shaped rib attaining a THP of 1.77 at Reynolds number ranging from Re 3,800 to Re 18,000 (), and a W-shaped rib producing a maximum THP of 1.87 at Reynolds number ranging from Re 3,800 to Re 18,000 (). These ribs are incorporated to advance the THP of the SAH. Research conducted by the authors demonstrates that the THP ranging from 1.72 to 1.87 has been achieved, with the maximum heat transfer rate obtained at an attack angle of 90o.
Several rib shapes have been explored by researchers to enhance the THP. used an L-shaped rib that yielded a maximum THP of 1.90, and Yadaba Mahanand () used an inverse T-shaped rib resulting in a higher THP of 1.87. A novel shape was articulated () through a detailed examination of circular ribs through analytical means. Variations in height and pitch distance were performed, which resulted in an exceptional thermal enhancement, surpassing that of the smooth duct by a factor of 1.65. Yadaba Mahanand () used a quarter circular rib in the SAH and observed that the Nu increased with increasing Reynolds number, while the ƒ increased with increasing rib pitch and a maximum value of THP was attained (1.88).
studied multiple arc-shaped rib configurations with varying gap settings and rib arc angles, ranging between 30o and 75o. Their study revealed that the Nu was 5.85 times greater than that of the smooth duct. Similarly, studied the tapered shape of ribs in rectangular ducts, and their results showed that the use of ribs enhanced heat transfer significantly. Ribs with angles ranging between 1.6o and 3.2o and varying pitch distances from 10 mm to 25 mm were examined. The report unveiled that, at a Reynolds number of 12,000, an inclination of 1.6° resulted in the most significant thermal enhancement. In another study, investigated the influence of ribs in the SAH, and they found the thermal efficiency in the SAH to be augmented significantly, utilizing diverse shapes and placements. It was found that the polygonal-shaped rib with minimum pitch distance yields higher thermal performances than ribs of other shapes and pitch distances. took a semi-circular rough surface for enhancement of THP at a relative pitch distance of 5–25 at Re 3,800 to Re 18,000. Researchers revealed that a maximum THP value of 1.98 had been obtained at Re 15,000. presented a detailed overview of various rough surfaces with different pitch distances, rib height in Reynolds number ranging from 2000 to Re 20,000, and their thermal enhancement in different applications. investigated V, arc, and transverse broken shape of the rib at various attack angles of 45°–75°, pitch distance of 15–25, relative roughness height 0.024–0.036, and relative wire length of 0.14–0.2. The analytical report presented that attack angles of 45° produced higher thermal performances at Re 8,750, and the THP is 1.63. Rajesh conducted an exergy analysis in the SHC using delta winglets at various attack angles in Re ranging from 4,000 to 18,000. The researchers reported that an angle of 60° achieved 29.88% exergy efficiency at a pitch distance of 5. proposed a new design NACA 0020 to enhance the THP at p = 13 and roughness height 0.06 at Re ranging from 6,000 to 18,000. The result revealed that a maximum THP of 2.21 was attained at Re 6,000. investigated the conical protrusion rib at a relative pitch distance of 6–12 and relative height of 0.020–0.044 at Re 4,000 to 16,000. The researchers reported extreme thermal efficiency achieved at 69.8% and recorded an efficiency enhancement factor of 1.346%. took a semi-elliptical rough surface with an attack angle of 30°–90° and positioned it in an inline and staggered technique at Re 6,000 to Re 18,000 to achieve higher THP. The results show that an attack angle of 75° produced higher Nu and ƒ of 2.05 and 6.93 for inline and 1.73 and 6.12 for staggered arrangement, respectively. investigated an S-shaped dimpled rib surface to augment the THP in the SAH. The result revealed that Nu was 3.55 times and ƒ 2.81 to 4.54 times greater than that of a smooth surface and the maximum THP achieved was 2.15.
The literature demonstrates a wide range of strategies employed to improve the thermal efficiency of the SAH, with a particular focus on the use of various rib and groove shapes. The studies have used different rib shapes, such as rectangular, inclined trapezoidal, serrated, semicircular, and helical, and have reported significant improvements in the heat transfer characteristics in the SAH by varying the rib pitch distance, rib height, and rib positioning angle. However, despite the extensive research in this area, there is still a research gap that needs to be addressed. Specifically, there is a lack of comparing rib positioning on both sides of the wall in the inline and staggered arrangement. Moreover, ribs placed in the center of the duct were not compared with those observed in earlier research. Therefore, further research is needed to compare the effectiveness of ribs placed in different parts of SAHs and their ability to enhance heat transfer and minimize pressure drop in SAHs. Such studies can help identify the most suitable rib shape for a given application and optimize the design of SAHs for extreme performance.
The main aim of the research work is as follows:
1. Explore the impact of relative roughness pitch and position of ribs on the average heat transfer and flow friction characteristics.
2. Optimize artificial roughness and flow parameters to improve the thermal performance of an artificially roughened SAH while maintaining a constant pumping power requirement.
The investigation methodology comprised three steps. First, the smooth surface of the SAH was studied by varying the entry flow velocity based on the Re to investigate the augmentation of the Nusselt number (Nu) and friction factor (f). Second, the proposed test section employed rough surfaces individually, by changing pitch distances, and THP and pressure drop are measured. Third, the results obtained were compared to determine the heat transfer enhancement achieved by the application of rough surfaces.
2 Numerical procedure
illustrates a schematic layout of the SAH, while
Figure 1. (ii) displays the proposed geometries for the b) polygonal rib in an inclined position, c) polygonal rib in a staggered position, d) polygonal rib and groove, and e) polygonal rib at the center of the rectangular duct. The rectangular duct of the SAH shown in
Figure 1(i) has dimensions of 640 mm in length (245 mm for the entry section—L1, 280 mm for the test section—L2, and 115 mm for the exit section—L3) with a height (H) and width (W) of 20 mm and 100 mm, respectively, as described in
. The proposed geometry with their various positions are presented in detail in
Figure 2. Furthermore, the rib arrangement is outlined in
Table 1.
Figure 3demonstrates the position of the polygonal rib and groove in a rectangular shape at varying pitch distances of
p= 10 mm and
p= 20 mm, and sharp edges with tilted profiles of the rib, which were selected to disrupt the fluid flow and maximize the turbulence in the flow region of laminar sub-layer to enhance heat transfer. To avoid unnecessary complexity, the present computational fluid dynamics (CFD) investigation followed certain assumptions.
a. The flow considered in this study is steady state, turbulent, and incompressible ().
b. The flow fluid is air, which enters the duct at atmospheric conditions. The properties of air are listed in Table 2.
c. A no-slip condition is applied to the bottom and sidewall of the duct.
FIGURE 1
FIGURE 2
TABLE 1
| Type | Shape | Pitch distance P (mm) | Hydraulic diameter D (mm) | Rib height, e (mm) | Relative roughness pitch, p/e | Relative roughness height, e/D |
|---|---|---|---|---|---|---|
| 1 | Polygonal rib | 10 | 33.33 | e1 = 3 | 3.33 | 0.09 |
| 2 | 20 | e3 = 1.5 | 13.33 | 0.045 |
Arrangement of artificial roughness.
FIGURE 3
TABLE 2
| S. No | Property | Air |
|---|---|---|
| 1 | Density, ρ (kg m−3) | 1.117 |
| 2 | Specific heat, Cp (J Kg-1 K−1) | 1,007 |
| 3 | Thermal conductivity, K (W m−1 K−1) | 0.0262 |
| 4 | Viscosity, µ (N m−2) | 1.857e−05 |
| 5 | Prandtl number, Pr | 0.71 |
Present study—(air) working fluid—(thermo-physical property) ().
2.1 Grid-independent study
To ensure accurate results for the proposed geometry, a grid-independent study is needed to be conducted using ANSYS Fluent version 13.0. In this study, the element sizes were varied between 0.4 mm and 0.2 mm, and four different values were used to identify the appropriate element cell size for the analysis. The computed values for Nu and ƒ were closer to acceptable results when using an element cell size of 0.2 mm than the previously used element cell size (). Figure 4A,B demonstrate the variation in the grid element size and node count on a smooth surface of an SAH and shows how the change in grid size affects a smooth surface in the SAH, while Figure 5 shows the mesh generated for the polygonal rib positioned at a) inline, b) staggered, c) rib and groove, and d) center placed at a pitch distance of p = 20 mm. Creation of a non-uniform mesh in the shape of a tetrahedron grid and convergence was taken between 114,565 and 1 52,145 element cells, depending on the proposed element size (), and the details are presented in Table 3. In addition to CFD simulations, the y + values are crucial for assessing result accuracy, particularly in the context of near-wall flow behavior. In the current simulation, the y + value achieved was 0.1993, which is a strong indicator that the simulation is accurately capturing the details of flow near surfaces or walls. This low y + value suggests that the mesh resolution near the walls is sufficient to model the boundary layer and viscous effects with precision. In essence, a y + value of 0.1993 signifies that the simulation’s near-wall predictions are reliable and trustworthy. The maximum deviation in Nu was observed in the 152,145 element cells, while the average Nu ranges showed the least deviation of +1.50%. Therefore, 152,145 element cells were used for the analytical investigation.
FIGURE 4
FIGURE 5
TABLE 3
| S. No | Element size (mm) | Number of elements | Nusselt number | % Deviation |
|---|---|---|---|---|
| 1 | 0.4 | 114,565 | 77.9 | 4.1 |
| 2 | 0.36 | 127,489 | 82.3 | 3.4 |
| 3 | 0.32 | 149,981 | 83.7 | 1.6 |
| 4 | 0.3 | 152,145 | 84.5 | 0.8 |
Grid-independent study.
2.2 Governing equations
The examination of heat transfer within a working fluid in an SAH incorporating a rough surface in simulations involves the application of the following formulated steady-state 2D continuity, momentum, and energy equations ():
Mass conservation or continuity equation:
Momentum conservation or Navier–Stokes equation:
Energy conservation equation:
and t are molecular thermal diffusivity and turbulent thermal diffusivity, respectively, represented by
Additionally, non-dimensional flow parameters are assessed using the following equation:
i)
ii) The friction factor was computed using the following equation…
It is essential to study the heat transfer and pressure drop characteristics of the SAH through THP. The THP with constant pumping, as referenced in , is determined by the following equation:
THP parameter:
i) Nur represents the Nu in a rough surface
ii) Nus represents an Nu for a smooth surface and was taken from the Dittus–Boelter equation (; ):
The Reynolds number varies from 3,800 to 18,000 (
), and the Prandtl number (Pr) of the working fluid air = 0.71
Table 2.
i) fr represents the friction factor for rough surfaces.
ii) fs represents a friction factor for the smooth surface of an SAH, obtainable from the Blasius equation ()
2.3 Identification of the turbulence model
Validating turbulent models is crucial for accurately predicting thermal performance enhancement in computational exercises. This research validates the average Nusselt number for a smooth surface by employing empirical equations such as Dittus–Boelter (9) and Taler (10), in conjunction with various turbulence models:
a) Standard k-ε epsilon,
b) RNG k-ε epsilon,
c) Realizable k-ε epsilon,
d) Standard k-ω omega, and
e) SST k-ω omega.
The study outcomes are compared and visually presented in Figure 6. The renormalization k-ε epsilon model is found to have the closest Nu values to the empirical correlation, with only a + 0.35% deviation. As a result, the RNG k-ε epsilon model is preferred for simulating the SAH in this research, consistent with similar studies by other researchers (; ).
FIGURE 6
2.4 Transport equation for the RNG k-ϵ model
The renormalized group (RNG) k-ϵ model incorporates the following equation to investigate the impact of small-scale turbulence adjacent to the rough surface (rib) and wall surface in the direction of the fluid stream ().
and
The term Pk denotes the production of turbulent kinematic energy resulting from the velocity gradient and serves as an indicator of turbulence intensity.i) represents the effective turbulent viscosity
In this context, the turbulent viscosity represented by μt is calculated by combining k and ε in the following manner:
Turbulent model RNG k- constant values have default values as follows:
2.5 Solution methods
In the present study, ANSYS Fluent version 13.0 was used to employ the finite volume method (FVM). The governing equations of continuity, momentum, and energy were solved by designing artificial roughness, generating the grid, and ensuring solution convergence. The specified Reynolds numbers, as listed in Table 4, were used to set the boundary conditions for the inlet velocity magnitude, while a constant heat flux of 1000 W/m2 was applied to the upper wall surface in the test section and the outer surface of the polygonal rib as solar radiation. Atmospheric conditions were set as the outlet pressure. The pressure–velocity coupling in the simulation utilized a SIMPLE algorithm, and a second-order double-precision was opted for the computations. Second-order discretization was applied to all transport equations. Criteria for continuity, x-velocity, and y-velocity were set at 10−6 to ensure accurate results during residual monitoring (). The standard initialization method was employed for computing all zones. In the calculation part, the number of iterations and reporting intervals were adjusted depending on the rib placement within the duct. After achieving convergence, the average Nusselt number (Nu) and average friction factor (f) values were compared with empirical Equations (9, 10, and 11) to assess thermal enhancement.
TABLE 4
| shape | Reynolds Number | Calculating u | |||||
|---|---|---|---|---|---|---|---|
| 3,800 | 5,000 | 8,000 | 12,000 | 15,000 | 18,000 | from Re no. | |
| Polygonal rib | 1.91 m/s | 2.51 m/s | 4.03 m/s | 6.04 m/s | 7.55 m/s | 9.06 m/s | |
Inlet velocity taken from Reynolds number ().
3 Validation of numerical outcomes through comparison with previously published findings
In our current study, we have conducted a comparative analysis of the smooth surface analytical data concerning average Nu and ƒ characteristics. We compared our findings with significant results from both numerical and experimental investigations. Those of , , and are a few investigations that are similar to ours that are presented in Figure 7.
FIGURE 7
Our observations indicate that the Nu values exhibit the least variation in the data presented by ; , with the maximum deviation being limited to +5% in the case of . Additionally, the friction factor (f) data display a minimal variance in the values obtained from . This outcome instills confidence in our research, encouraging further exploration of the proposed geometric parameters.
4 Result and discussion
This numerical analysis investigates the mean Nu, ƒ, and THP of the projected rough surface rib by varying the relative roughness pitch p = 10 mm and p = 20 mm, and the relative roughness height e/D = 0.9 performance of the rib shapes is evaluated by contrasting them with established outcomes from published studies to determine the extent of thermal augmentation.
4.1 Nusselt number enhancement
Figure 8 shows the Nusselt number for the proposed rib shapes in the SAH, indicating a continuous increase in the Nu with the Reynolds number for all rough ribs due to forced convection. Additionally, the rib configuration induces a vortex that disrupts the secondary layer of the flow region, leading to improved heat transfer. As the Reynolds number increases, the length of recirculation in the rectangular duct decreases, thereby enhancing thermal performance. Among the rib shapes, the staggered polygonal rib at Re = 18,000 for p = 20 mm and center rib for p = 10 mm significantly increase the heat transfer rate primarily due to the formation of eddies near the ribs positioned in opposite directions, which were closer to the laminar sub-layer on both sides of the wall, and closer pitch distances in the center rib produced greater heat transfer than others.
FIGURE 8
It was observed that the vortex size near the initial set of ribs closest to the boundary-wall surface was gradually expanding while boosting the Reynolds number. This expansion results in enhanced local convective heat transfer, as indicated by the contour plots of TKE in Figure 9 and Figure 10. These plots are based on relative pitch distances of p = 20 mm and p = 10 mm, with the Reynolds number ranging from 3,000 to 18,000. In the first scenario (case 1), ribs positioned inline lead to the development of strong turbulence originating from the first pair of ribs. These turbulent effects intensify as they propagate downward. For the second scenario (case 2), which involves rib and groove configurations, a high level of circulation is generated when p = 10 mm, while a robust vortex forms at p = 20 mm. In the third case (case 3), which features staggered rib positioning, vortices are generated on both sides of the laminar surface close to the wall. A substantial circulation is observed at higher Reynolds numbers (Re = 18,000) in the upper wall region, while lower Reynolds numbers (Re = 3,800) result in circulation on the bottom wall. Both the entry and exit regions are depicted in Figure 9 and Figure 10. In the fourth case (case 4), where ribs are positioned at the center, strong turbulence is generated at lower Reynolds number, and minimal turbulence occurs at higher Reynolds number.
FIGURE 9
FIGURE 10
The fluid flow velocity and their contour plot results in Figure 11 and Figure 12 demonstrate that increasing inlet velocity results in higher turbulent intensity, causing stronger vortices near the staggered and center-positioned rib at a relative roughness pitch of 10 mm compared to 20 mm shown in Figure 11 and Figure 12 (c-entry and exit of the test section). This results in more robust turbulence and increased heat energy transfer as the Reynolds number increases in the SAH.
FIGURE 11
FIGURE 12
The pronounced interaction between the initial pair of ribs amplifies the average Nusselt number along the streamwise distribution. Figure 8 provides the Nu enhancement ratio for the suggested rib configuration, with the highest Nu/Nus ratio of 3.76 observed for the polygonal-shaped rib located at the center with a relative pitch distance of p = 10 mm and at a Re 3,800.
The fluid flow velocity contour plot in Figure 13 and Figure 14 at p = 20 mm and p = 10 mm clearly illustrates the influence of various polygonal-shaped ribs on the absorber plate of the SAH. It indicates that the size of the circulation generated in the flow direction is larger for the envisaged rib shapes at closer relative pitch distances of p = 10 mm at a Reynolds number of 18,000. This phenomenon is linked to the impact of the rib’s apex edge shape, which is then accompanied by a subsequent change in the rib’s orientation. This combination creates a strong vortex in the vicinity of the initial rib. This induces strong turbulence near the wall surface in the flow direction, leading to a higher heat transfer rate. Furthermore, it was observed from all rib shapes that the center polygonal rib had a more significant impact than other ribs at both pitch distances, as shown in Figure 13 and Figure 14. (ii) This observation suggests that the average Nu increases with increasing Reynolds number, primarily because higher Reynolds numbers lead to increased recirculation in the fluid flow direction, thereby resulting in a higher convective heat transfer coefficient.
FIGURE 13
FIGURE 14
4.2 Friction factor characterization
Implementing a rough surface in the SAH may lead to extra friction losses compared to a smooth surface, but it is crucial to examine the hydraulic (heat transfer) performances in the SAH to achieve a better heat transfer rate. For maintaining a balance in pressure drop within a duct across different Re values, a constant pumping power is necessary. Figures 15A,B illustrate the ƒ characteristics for dissimilar polygonal rib shapes with relative roughness pitch p = 20 mm and p = 10 mm in varying Reynolds numbers. It can be inferred that the staggered positioned rib has a record least average ƒ, as shown in Figure 15. The findings suggest that the average ƒ decreases as the Reynolds number increases, which is attributed to the conquest of the viscous region of the sub-layer in a fully developed turbulent flow. The staggered polygonal-shaped rib, due to the intermission in the flow path caused by the enhanced relative roughness height of the rib, produces lower ƒ. In addition, it is observed that greater pitch distances lead to a reduction in the average friction factor, especially at the maximum Reynolds number of 18,000. From the depicted data, it can be inferred that the average ƒ decreases as the pitch distance increases, as well as when the relative roughness height of the rib increases.
FIGURE 15
The pressure contour plots for relative pitch distances of p = 20 mm and p = 10 mm are illustrated in Figure 16 and Figure 17. These plots unveil that, at a higher Reynolds number Re = 18,000, there is a sudden increase in pressure near the entry region of the test section due to the presence of the rough surface. Notably, a significant impact is observed, with a larger circulation created at the smaller pitch distance of p = 10 mm compared to p = 20 mm in the case of ribs positioned in-line. A similar pattern is observed in the rib and groove configurations. In contrast, staggered rib positioning and center rib positioning result in the generation of substantial turbulence that is symmetrically directed toward the wall region. It is noteworthy that at a maximum Reynolds number Re = 15,000 and a pitch distance of p = 20 mm, a lower friction factor is obtained, whereas at Re = 18,000 and p = 10 mm, the ƒ is higher.
FIGURE 16
FIGURE 17
4.3 Thermo-hydraulic performances
The aim of the present research work is to investigate the impact of the position of polygonal-shaped ribs with different relative pitch distances on heat transfer and friction in an SAH (). Eq. 8 is used to evaluate the thermal and hydraulic performance of the SAH, and the resulting average Nu and friction factor values are presented in Figure 18 for relative pitch distances of p = 20 mm and p = 10 mm. The analysis illustrates that the thermo-hydraulic performance parameter (THPP) values range from 0.43 to 2.81 for all proposed rib shapes, as presented in Table 5. It is noted that decreasing the relative pitch distances of the ribs leads to an increase in the average Nu due to boundary layer interruption, while increasing the pitch distances leads to a decrease in the average ƒ. The staggered polygonal-shaped rib with a relative pitch distance of p = 10 mm at Reynolds number 8000 achieved the highest THP value of 2.81. Table 6 provides a comparative evaluation of rough surface shapes that have been previously published, alongside their associated THP values. This comparison showcases a noteworthy improvement associated with the projected rib, as validated by the thermal enhancement factor. The table further illustrates the diverse configurations of polygonal-shaped ribs, each characterized by varying relative pitch distance.
FIGURE 18
TABLE 5
| Position of ribs | Re = 18,000 | Re = 15,000 | Re = 12,000 | Re = 8,000 | Re = 5,000 | Re = 3,800 | ||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| p = 20 mm | p = 10 mm | p = 20 mm | p = 10 mm | p = 20 mm | p = 10 mm | p = 20 mm | p = 10 mm | p = 20 mm | p = 10 mm | p = 20 mm | p = 10 mm | |
| Inline | 0.589 | 1.917 | 0.848 | 2.591 | 0.811 | 2.584 | 0.848 | 2.387 | 0.989 | 2.700 | 0.912 | 2.485 |
| Staggered | 0.714 | 2.709 | 1.305 | 2.731 | 1.299 | 2.743 | 1.379 | 2.813 | 1.185 | 2.651 | 1.094 | 2.172 |
| Rib and groove | 0.476 | 2.276 | 0.822 | 2.337 | 0.918 | 1.932 | 0.984 | 1.868 | 1.240 | 1.875 | 1.000 | 1.648 |
| Centre | 0.435 | 1.968 | 0.739 | 2.741 | 0.786 | 2.418 | 0.985 | 1.523 | 1.256 | 1.522 | 1.140 | 1.447 |
Thermo-hydraulic performance parameter.
Bold values represents that the maximum value among the attained values.
TABLE 6
| S. No | Reference | Geometry | Parameter | Maximum THP |
|---|---|---|---|---|
| 1 | Square ribs | Re: 3,800–18,000 | 1.82 | |
| 2 | Circular ribs | Re: 3,800–18,000 | 1.65 | |
| 3 | Discrete W-shaped | Re: 3,800–18,000 | 1.72 | |
| 4 | Circular wire in an arc shape | Re: 6,000–18,000 | 1.70 | |
| 5 | Inverted U-shaped rib | Re: 3,800–18,000 | 1.82 | |
| 6 | Quarter circular rib | Re: 3,800–18,000 | 1.88 | |
| 7 | Present study | Polygonal–staggered rib | Re: 3,800–18,000 | 2.81 |
THP values compared with earlier published work ().
Bold values represents that the maximum value among the attained values.
It assesses heat transfer improvement through thermal performance and evaluates the additional frictional resistance through hydraulic performance (). Eq. 8 is used to predict the THP of the SAH, and the resulting average Nu and ƒ values are plotted in Figure 18 for pitch distances of p = 20 mm and p = 10 mm. By analyzing the present proposed shapes with varying relative pitch distances, it was observed that the THPP value varied between 0.43 and 2.81, as shown in Table 5. This variation in THPP helps identify the effectiveness of the proposed ribs in enhancing the THP of the SAH. The observation indicates that reducing the pitch distances between the ribs significantly disrupts the boundary layer, leading to an increase in the average Nu.
Conversely, increasing pitch distances results in a decrease in the average friction factor. Notably, the staggered polygonal rib shape with a relative pitch distance of p = 10 mm at an Re of 8000 achieved a maximum THP of 2.81. The THP values obtained were compared with those of previously published rib shapes and their respective THP outcomes, as outlined in Table 6. This comparison underscores a significant enhancement in the proposed rib configuration, a conclusion further validated by the thermal enhancement factor.
5 Conclusion
The study analyzed the numerical results of a 2D rectangular duct in an SAH with artificial rough surfaces of polygonal-shaped ribs placed at various positions inside the duct. The study investigated the impact of relative roughness pitch distances on the average Nu, ƒ, and THP parameters. The following conclusions were made
1. The study’s analysis unequivocally validates that the incorporation of artificial roughness in a rectangular duct significantly elevates the heat transfer coefficient compared to a smooth surface.
2. Diverse turbulent models were assessed, and it was found that the RNG k-ϵ model consistently produces results closely aligned with the empirical values of the Dittus–Boelter equation, bolstering confidence in its applicability.
3. The investigation identifies that the most effective heat transfer enhancement is achieved by using polygonal-shaped ribs positioned at the center of the duct with a relative pitch distance of p = 10 mm, particularly in cases of lower Reynolds numbers (Re 3800). This configuration attains a maximum average Nusselt number of 3.76, surpassing the 3.42 Nusselt number recorded at p = 20 mm and Re 3,800.
4. Notably, the study reveals a substantial increase in the friction factor (ƒ) with the introduction of artificial roughness. At a relative pitch distance of p = 20 mm and Reynolds number of 18000, the maximum ƒ is 3.66 times that of a smooth surface, and at p = 10 mm, it is 2.6 times higher.
5. In the context of staggered rib positioning at a relative pitch distance of p = 10 mm and Reynolds number Re 8000, the research identifies the attainment of the highest peak THP, which reaches 2.813.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
BK: conceptualization, data curation, formal analysis, investigation, validation, and writing–original draft. CS: data curation, investigation, software, and writing–original draft. PK: conceptualization, formal analysis, methodology, writing–original draft, and writing–review and editing. DT: data curation, formal analysis, supervision, and writing–review and editing. MS: data curation, validation, and writing–review and editing. JT: conceptualization, supervision, validation, and writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Nomenclature
| L | Length of the duct, mm |
| L1 | Inlet length of the duct, mm |
| L2 | Test length of the duct, mm |
| L3 | Outlet length of the duct, mm |
| Cp | Specific heat of air, J/kg k |
| D | Equivalent or hydraulic diameter of the duct, mm |
| e | Rib height, mm |
| h | Heat transfer coefficient, W/m2 k |
| H | Depth of the duct, mm |
| I | Turbulence intensity/intensity of solar radiation, W/m2 |
| k | Thermal conductivity of air, W/m k |
| m | Mass flow rate, kg/s |
| dp | Pressure drop, Pa |
| p | Pitch, mm |
| To | Fluid outlet temperature, K |
| Ti | Fluid inlet temperature, K |
| Ta | Ambient temperature, K |
| v | Velocity of air in the duct, m/s |
| W | Width of the duct, mm |
| TKE | Turbulent kinetic Energy |
| I | Turbulent intensity |
| Dimensionless parameters | |
| Nu | Nusselt number |
| Nur | Nusselt number for the rough duct |
| Nus | Nusselt number for the smooth duct |
| f | Friction factor |
| fr | Friction factor for a rough surface |
| fs | Friction factor for a smooth surface |
| Pr | Prandtl number |
| Re | Reynolds number |
| St | Stanton number |
| W/H | Duct aspect ratio |
| Greek symbols | |
| ρ | Density of air, kg/m3 |
| ɛ | Dissipation rate |
| ω | Specific dissipation rate |
| δ | Transition sub-layer thickness, m |
| µ | Dynamic viscosity, Ns/m2 |
| µt | Turbulent viscosity, Ns/m2 |
| k | Turbulence kinetic energy, m |
| Γ | Molecular thermal diffusivity |
| Γt | Turbulent thermal diffusivity |
References
1
AlamT.KimM.-H. (2016). Numerical study on thermal hydraulic performance improvement in solar air heater duct with semi ellipse shaped obstacles. Energy, V.112, 588–598. 10.1016/j.energy.2016.06.105
2
AlamT.KimM.-H. (2017). Heat transfer enhancement in solar air heater duct with conical protrusion roughness ribs. Appl. Therm. Eng.126, 458–469. 10.1016/j.applthermaleng.2017.07.181
3
AryaN.GoelV.SundenB. (2023). Solar air heater performance enhancement with differently shaped miniature combined with dimple shaped roughness: CFD and experimental analysis. Sol. Energy250, 33–50. 10.1016/j.solener.2022.12.024
4
ASHRAE (2003). Method of testing to determine the thermal performance of solar collectors. Atlanta, GA: American Society of Heating, Refrigeration and Air Conditioning Engineers, 42.
5
BopcheS. B.TandaleM. S. (2009). Experimental investigations on heat transfer and frictional characteristics of a turbulator roughened solar air heater duct. Int. Commun. Heat Mass Transf.52, 2834–2848. 10.1016/j.ijheatmasstransfer.2008.09.039
6
DuffieJ. A.BeckmanW. A. (1980). Solar engineering of thermal processes. New York: Wiley.
7
Eiamsa-ardS.PromvongeP. (2008). Numerical study on heat transfer of turbulent channel flow over periodic grooves. Int. Commun. Heat. Mass Transf.35, 844–852. 10.1016/j.icheatmasstransfer.2008.03.008
8
Eiamsa-ardS.PromvongeP. (2009). Thermal characteristics of turbulent rib-grooved channel flows. Int. Commun. Heat. Mass Transf.36, 705–711. 10.1016/j.icheatmasstransfer.2009.03.025
9
GawandeV. B.DhobleA. S.ZodpeD. B.ChamoliS. (2016). Experimental and CFD investigation of convection heat transfer in solar air heater with reverse L-shaped ribs. Sol. Energy131, 275–295. 10.1016/j.solener.2016.02.040
10
KumarA.BhagoriaJ. L.SarviyaR. M. (2009). Heat transfer and friction correlations for artificially roughened solar air heater duct with discrete W-shaped ribs. Energy Conversat. Manag.50, 2106–2117. 10.1016/j.enconman.2009.01.025
11
KumarS.DasR. K.KulkarniK.AlamT.EldinS. M. (2023). Designing of low-cost solar air heater equipped with roughness of streamlined cross-section. Case Stud. Therm. Eng.45, 102915, 102915. 10.1016/j.csite.2023.102915
12
MahanandY.Ranjan SenapatiJ. (2020). Thermal enhancement study of a transverse inverted-T shaped ribbed solar air heater. Int. Commun. Heat Mass Transf.119, 104922. 10.1016/j.icheatmasstransfer.2020.104922
13
MahanandY.SenapatiJ. R. (2021). Thermo-hydraulic performance analysis of a solar air heater (SAH) with quarter-circular ribs on the absorber plate: a comparative study. Int. J. Therm. Sci.161, 106747. 10.1016/j.ijthermalsci.2020.106747
14
MaithaniR.KumarA.AliM. A.GuptaN. K.SharmaS.AlamT.et al (2023). Exergy and sustainability analysis of a solar heat collector with wavy delta winglets as turbulent promoters: a numerical analysis. Case Stud. Therm. Eng.49, 103293. 10.1016/j.csite.2023.103293
15
McAdamsW. H. (1942). Heat transmission. New York: McGraw-Hill.
16
PandeyN. K.BajpaiV. K.Varun (2016). Experimental investigation of heat transfer augmentation using multiple arcs with gap on absorber plate of solar air heater. Sol. Energy134, 314–326. 10.1016/j.solener.2016.05.007
17
PromvongeP.SkullongS.KwankaomengS.ThiangpongC. (2012). Heat transfer in square duct fitted diagonally with angle-finned tape–Part 1: experimental study. Int. Commun. Heat. Mass Transf.39 (5), 617–624. 10.1016/j.icheatmasstransfer.2012.03.007
18
PromvongeP.ThianpongC. (2008). Thermal performance assessment of turbulent channel flows over different shaped ribs. Int. Commun. Heat. Mass Transf.35, 1327–1334. 10.1016/j.icheatmasstransfer.2008.07.016
19
SainiR. P.SainiS. K. (2008). Development of correlations for Nusselt number and friction factor for solar air heater with roughened duct having arc-shaped wire as artificial roughness. Sol. Energy82, 1118–1130. 10.1016/j.solener.2008.05.010
20
SainiR. P.VermaJ. (2008). Heat transfer and friction factor correlations for a duct having dimple shape artificial roughness for solar air heaters. Energy33, 1277–1287. 10.1016/j.energy.2008.02.017
21
SinghJ.Singh BishtV.BhandariP.KumarK.SinghJ.AlamT.et al (2023). Computational parametric investigation of solar air heater with dimple roughness in S-shaped pattern. Int. J. Interact. Des. Manuf. (IJIDeM. 10.1007/s12008-023-01392-8
22
Singh YadavA.AgrawalA.SharmaA.GuptaA. (2022b). Revisiting the effect of ribs on performance of solar air heater using CFD approach. Mater. today Proc.63, 240–252. 10.1016/j.matpr.2022.02.549
23
Singh YadavA.AlamT.GuptaG.SaxenaR.GuptaN. K.Viswanath AllamrajK.et al (2022a). A numerical investigation of an artificially roughened solar air heater. Energies15 (21), 8045. 10.3390/en15218045
24
Singh YadavA.BhagoriaJ. L. (2013). A CFD (computational fluid dynamics) based heat transfer and fluid flow analysis of a solar air heater provided with circular transverse wire rib roughness on the absorber plate. Energy55, 1127–1142. 10.1016/j.energy.2013.03.066
25
Singh YadavA.BhagoriaJ. L. (2014). A numerical investigation of square-sectioned transverse rib roughened solar air heater. Int. J. Therm. Sci.79, 111–131. 10.1016/j.ijthermalsci.2014.01.008
26
SripattanapipatS.PromvongeP. (2009). Numerical analysis of laminar heat transfer in a channel with diamond-shaped baffles. Int. Commun. Heat. Mass Transf.36, 32–38. 10.1016/j.icheatmasstransfer.2008.09.008
27
TalerD.JanT. (2017). Simple heat transfer correlations for turbulent tube flow. E3S Web Conf.13, 02008. 10.1051/e3sconf/20171302008
28
VarshneyL.GuptaA. D. (2017). Performance prediction for solar air heater having rectangular sectioned tapered rib roughness using CFD. Therm. Sci. Eng. Prog.4, 122–132. 10.1016/j.tsep.2017.09.005
29
Varun Kumar.B.ManikandanG.Rajesh KannaP.TalerD.TalerJ.Nowak-OcłońM.et al (2018). A performance evaluation of a solar air heater using different shaped ribs mounted on the absorber plate—a review. Energies11, 3104. 10.3390/en11113104
30
Varun KumarBManikandanGRajesh KannaG. P. (2021). Enhancement of heat transfer in SAH with polygonal and trapezoidal shape of the rib using CFD. Energy, V.234, 121154. 10.1016/j.energy.2021.121154
31
VermaS.PrasadB. (2000). Investigation for the optimal thermohydraulic performance of artificially roughened solar air heaters. Renew. Energy20, 19–36. 10.1016/S0960-1481(99)00081-6
32
WebbR. L.EckertE. R. G. (1972). Application of rough surfaces to heat exchanger design. Int. J. Heat. Mass Transf.15 (9), 1647–1658. 10.1016/0017-9310(72)90095-6
Summary
Keywords
rough surface, polygonal rib with groove, THP, solar air heater, Nusselt number, numerical simulation
Citation
Kumar BV, Selvan CP, Kanna PR, Taler D, Szymkiewicz M and Taler J (2023) Numerical investigation of heat transfer enhancement in solar air heaters using polygonal-shaped ribs and grooves. Front. Energy Res. 11:1279225. doi: 10.3389/fenrg.2023.1279225
Received
17 August 2023
Accepted
30 October 2023
Published
28 November 2023
Volume
11 - 2023
Edited by
Mahendran Samykano, Universiti Malaysia Pahang, Malaysia
Reviewed by
Tabish Alam, Central Building Research Institute (CSIR), India
Subbarama Kousik Suraparaju, Sri Vasavi Engineering College, India
Updates
Copyright
© 2023 Kumar, Selvan, Kanna, Taler, Szymkiewicz and Taler.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jan Taler, jan.taler@pk.edu.pl; P. Rajesh Kanna, prkanna@gmail.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.