Investigation of Thermo-Hydraulic Performances of Artiﬁcial Ribs Mounted in a Rectangular Duct

: This research investigates the ﬂuid ﬂow characterization and thermohydraulic performances (THP) of rib surfaces, using computational and experimental methods. ANSYS computational ﬂuid dynamics (CFD) software was used to predict and validate the ﬁndings in the experimental setup. Artiﬁcial rib surfaces, including polygonal and forward trapezoidal-shaped ribs, were placed in the absorber plate at different relative pitch distances (p/e) = 6.7, 10, 13.4 and relative height (e/d) = 20, and the mass ﬂow rate of air (working ﬂuid) varied at Reynolds numbers ranging from 2000 to 20,000. According to the validation results, the RNG renormalization k- ε model was selected for the investigation. The results show that strong turbulence occured closer to the wall surface and behind the rib surface, enhancing thermal performances due to the sharp edge shape of the rib. A polygonal rib with a pitch distance of p/e = 6.7 achieved a higher Nusselt number (Nu) and thermo-hydraulic performance of 2.95 at Re 4000. An empirical correlation between the Nusselt number (Nu) and friction factor ( f ) was developed using linear regression analysis and was compared with the predicted values. The comparison results show a close range of ± 8% between the experimental and predicted values.


Introduction
The availability of solar energy is abundant in nature and can be used for drying agro-based food products such as moringa leaf, dhal, banana, etc, to convert value addition products by removing moisture levels. A solar air heater (SAH) compact device (space heating) plays a significant role in removing moisture content and is suitable for drying vegetables, fruits, grains, etc. Traditional SAH models receive heat energy directly from solar radiation and produce minimum thermal performances inside the rectangular duct. Many researchers have conducted investigations via numerical and experimental methods to augment heat efficiency. Implementing artificial rib surfaces, such as ribs, grooves, etc in the absorber plate of the SAH can produce an enhancement of the heat transfer rate. Saha et al. [1] and Kim et al. [2] conducted an experimental study on a rectangular duct with ribbed walls to investigate the heat transfer and friction characteristics of artificial rib turbulators [1] and perforated ribs [2]. The results showed that using perforated ribs and turbulators increased heat transfer by 18% and pressure drop by 12% in comparison with the smooth surface. the thermal performances of longitudinal rib turbulators were numerically investigated by Kefayati et al. [3], twisted tape inserts by Ji et al. [4], and discrete multiple V-ribs by Balakrishnan and Khader [5]. The results indicated that the Nusselt number increased by 21% and there was a significant pressure drop compared with a smooth surface at a higher Reynolds number. Narayanan and Bhaskaran [6] and Li et al. [7] conducted their research numerically, using turbulent flow for augmentation of heat transfer in a rectangular duct with multiple ribs of varying geometry [6] and staggered saw tooth ribs [7]. They reported that using multiple ribs in the absorber plate increased heat transfer in comparison with other models and produced nominal pressure drop owing to optimum rib spacing.
Numerical simulation was carried out to investigate heat transfer and flow characteristics a rectangular duct's, for elliptical dimples rib by Wen et al. [8] and multiple dimple-shaped ribs by Wang et al. [9] The results showed that the use of multiple elliptical dimples significantly enhanced heat transfer and pressure drop compared with single elliptical dimples. Rajendra Karwa et al. [10] conducted an experimental investigation using integral chamfered ribs. The research results revealed enhancement of maximum thermal efficiency by 10% to 40%.
Ahmed et al. [11] and Sahiti and Cao [12] conducted both experimental and numerical investigations of the augmentation of heat transfer rate in a rectangular duct using wavy and perforated ribs [11] and serrated twisted tape [12]. The results showed that the use of wavy and perforated ribs increased heat transfer compared with others and pressure drop of 14% was attained. Siddiqui and Kim [13] performed experimental and numerical investigation of flow and heat transfer in a rectangular duct by using multiple V-shaped ribs. They reported that V-shaped ribs significantly increased heat transfer and pressure drop compared with a smooth surface and produced higher thermal performances. Islam and Alim [14] presented a numerical study on the thermal-hydraulic (THP) characteristics of turbulent flow in a rectangular duct with novel converging/diverging rib arrangements. The results showed that the use of these rib arrangements significantly enhanced heat transfer and pressure drop, and the optimum rib geometry was identified. Lee and Kang [15] and Seo and Kim [16] investigated the heat transfer and fluid flow characteristics in a rectangular duct with a multi-vortex generator and novel compound turbulator. The study showed that the multi-vortex generator significantly enhanced heat transfer and pressure drop, and the optimum rib geometry was identified. Karthikeyan and Kannan [17] presented a numerical study of fluid flow and heat transfer in a rectangular duct with a discrete vortex generator. The results showed that the use of the vortex generator significantly enhanced heat transfer and pressure drop, and the optimum vortex generator spacing was identified. Furthermore, various rib shapes such as transverse wire ribs [18], arc-shaped wire ribs [19], inverted U-shape ribs [20], dimple-shaped ribs [21], and W-shaped ribs [22] have been positioned in the absorber plate to augment THP by varying pitch distance, rib height, and attack angles. It was observed that higher THP values ranged between 1.72 to 1.85 and attack angle 90 • yielded a better heat transfer rate than other angles.
From the detailed literature review, we observed that implementing a rib surface in the absorber plate of SAH enhances thermal performance compared with a smooth surface. However, many studies focused only single side inclination in rib surface, and comparison between numerical and experimental study was limited. Moreover, implementing a Sharpe edge on both sides of the rib surface produces a higher heat transfer rate closer to the wall surface. In this study, we attempted to implement a proven polygonal rib surface [23][24][25] in an absorber plate at various pitch distances, to identify the maximum convective heat transfer rate and THP from the rib at various mass flow rates, both numerically and experimentally.

Experimental Procedure
The experimental setup for the SAH is presented in schematic form in Figure 1. The setup consists of several main components, including: (i) a rectangular duct divided into three sections, (ii) a heating plate,  Based on the Reynolds number, a centrifugal blower with a control valve regulates the mass flow rate/airflow through the connecting pipes. The entire rectangular duct is covered with insulating materials to prevent heat loss.

Test Section
The cross-sectional measurement for a rectangular wooden duct is 300 mm in width and 25 mm in height, and its entire length (L) is 2390 mm (H). The aspect ratio (W/H) was maintained at 12. Figure 2 shows the schematic views of the SAH and rib shapes [25] that were taken for the study. Rectangular duct lengths are (i) 890 mm for the entry portion, (ii) 500 mm for the exit section with plenum, and (iii) 1000 mm for the test section. As per the ASHRAE guidelines [26], the rectangular duct with entry and exit sections has been designed for the entry Section 5 √ and the exit Section 2.5 √ . Based on the Reynolds number, a centrifugal blower with a control valve regulates the mass flow rate/airflow through the connecting pipes. The entire rectangular duct is covered with insulating materials to prevent heat loss.

Test Section
The cross-sectional measurement for a rectangular wooden duct is 300 mm in width and 25 mm in height, and its entire length (L) is 2390 mm (H). The aspect ratio (W/H) was maintained at 12. Figure 2 shows the schematic views of the SAH and rib shapes [25] that were taken for the study. Rectangular duct lengths are (i) 890 mm for the entry portion, (ii) 500 mm for the exit section with plenum, and (iii) 1000 mm for the test section. As per the ASHRAE guidelines [26], the rectangular duct with entry and exit sections has been designed for the entry Section 5 √ WH and the exit Section 2.5 √ WH.  Based on the Reynolds number, a centrifugal blower with a control valve regulates the mass flow rate/airflow through the connecting pipes. The entire rectangular duct is covered with insulating materials to prevent heat loss.

Test Section
The cross-sectional measurement for a rectangular wooden duct is 300 mm in width and 25 mm in height, and its entire length (L) is 2390 mm (H). The aspect ratio (W/H) was maintained at 12. Figure 2 shows the schematic views of the SAH and rib shapes [25] that were taken for the study. Rectangular duct lengths are (i) 890 mm for the entry portion, (ii) 500 mm for the exit section with plenum, and (iii) 1000 mm for the test section. As per the ASHRAE guidelines [26], the rectangular duct with entry and exit sections has been designed for the entry Section 5 √ and the exit Section 2.5 √ .

Absorber Plate
An aluminum sheet of 0.8 mm thickness was used as an absorber plate for the experimental research. A polygonal and trapezoidal rib with a pitch distance of 20 mm (p/e = 6.7) was positioned. Ribs were placed one by one independently at a pitch distance of 30 mm and 40 mm (p/e = 10 and p/e = 13.5). Figure 3 illustrates the rough shapes of ribs with various pitch distances. The top surface of the absorber plate has a 6 mm thick nickel-chromium heated plate standing on it to provide heat energy equivalent to solar radiation at a constant heat flux of 1000 W/m 2 .

Absorber Plate
An aluminum sheet of 0.8 mm thickness was used as an absorber plate for the experimental research. A polygonal and trapezoidal rib with a pitch distance of 20 mm (p/e = 6.7) was positioned. Ribs were placed one by one independently at a pitch distance of 30 mm and 40 mm (p/e = 10 and p/e = 13.5). Figure 3 illustrates the rough shapes of ribs with various pitch distances. The top surface of the absorber plate has a 6 mm thick nickelchromium heated plate standing on it to provide heat energy equivalent to solar radiation at a constant heat flux of 1000 W/m 2 .

Measuring Instruments
The measurement of the mass flow rate and pressure drop across the orifice was conducted across the test section; pressure was measured using an inclined manometer. To ensure accuracy, butyl alcohol with a density of 800 kg/m 3 was used. Sixteen thermocouples with calibrated copper constants were placed in the SAH, as illustrated in Figure 5, to measure the orifice temperature, with a U-tube manometer between the plenum and the centrifugal blower to measure the pressure drop at different locations.

Absorber Plate
An aluminum sheet of 0.8 mm thickness was used as an absorber plate for the exper imental research. A polygonal and trapezoidal rib with a pitch distance of 20 mm (p/e = 6.7) was positioned. Ribs were placed one by one independently at a pitch distance of 30 mm and 40 mm (p/e = 10 and p/e = 13.5). Figure 3 illustrates the rough shapes of ribs with various pitch distances. The top surface of the absorber plate has a 6 mm thick nickel chromium heated plate standing on it to provide heat energy equivalent to solar radiation at a constant heat flux of 1000 W/m 2 .

Measuring Instruments
The measurement of the mass flow rate and pressure drop across the orifice was con ducted across the test section; pressure was measured using an inclined manometer. To ensure accuracy, butyl alcohol with a density of 800 kg/m 3 was used. Sixteen thermocou ples with calibrated copper constants were placed in the SAH, as illustrated in Figure 5 to measure the orifice temperature, with a U-tube manometer between the plenum and the centrifugal blower to measure the pressure drop at different locations.

Measuring Instruments
The measurement of the mass flow rate and pressure drop across the orifice was conducted across the test section; pressure was measured using an inclined manometer. To ensure accuracy, butyl alcohol with a density of 800 kg/m 3 was used. Sixteen thermocouples with calibrated copper constants were placed in the SAH, as illustrated in Figure 5, to measure the orifice temperature, with a U-tube manometer between the plenum and the centrifugal blower to measure the pressure drop at different locations.

Working Procedure
The experimental setup was designed according to ASHRAE guidelines [26]. To achieve steady-state conditions, the centrifugal blower was initially turned on and oper ated for 30 min to establish the fluid flow inside the channel. The mass flow rate was de termined using a U-tube manometer at a constant velocity. All sensors were calibrated and carefully mounted to ensure accurate data collection, while the mass flow rate within the conduit was verified to remain consistent. The upper surface of the absorber plate wa utilized as the heating plate, receiving a constant heat flux of 1000 W/m 2 . The test section consisted of a flat absorber plate, followed by a rough plate, tested sequentially to enhanc thermal performance. The data obtained were analyzed and compared on a plot in orde to determine the most effective approach. Various parameters were measured in the in vestigative work.
(a) In total, 12 different temperatures recorded in the test section (tp), inlet (ti), outle (to), and ambient temperature (ta) were measured for analysis (b) Pressure drop readings were observed at the exit section of the duct (∆p). (c) ƒ was measured by the pressure drop across the test section(∆po).
To record readings in the SAH's smooth and rib surface, the experiment was run constantly. The empirical equations developed by Dittus and Boelter and by Blausis wer applied to validate the experimental findings for the smooth surface of Nu, and f wa subsequently expanded to account for surface roughness. This led to the development o the empirical correlation for the experimental Nu and f from the regressive analysis fo further research.
The following equation was used to calculate the values of WR: (1 The relative error can be calculated using the following equation:

Working Procedure
The experimental setup was designed according to ASHRAE guidelines [26]. To achieve steady-state conditions, the centrifugal blower was initially turned on and operated for 30 min to establish the fluid flow inside the channel. The mass flow rate was determined using a U-tube manometer at a constant velocity. All sensors were calibrated and carefully mounted to ensure accurate data collection, while the mass flow rate within the conduit was verified to remain consistent. The upper surface of the absorber plate was utilized as the heating plate, receiving a constant heat flux of 1000 W/m 2 . The test section consisted of a flat absorber plate, followed by a rough plate, tested sequentially to enhance thermal performance. The data obtained were analyzed and compared on a plot in order to determine the most effective approach. Various parameters were measured in the investigative work.
(a) In total, 12 different temperatures recorded in the test section (tp), inlet (ti), outlet (to), and ambient temperature (ta) were measured for analysis (b) Pressure drop readings were observed at the exit section of the duct (∆p). (c) ƒ was measured by the pressure drop across the test section(∆po).
To record readings in the SAH's smooth and rib surface, the experiment was run constantly. The empirical equations developed by Dittus and Boelter and by Blausis were applied to validate the experimental findings for the smooth surface of Nu, and f was subsequently expanded to account for surface roughness. This led to the development of the empirical correlation for the experimental Nu and f from the regressive analysis for further research.

Uncertainties Analysis
The primary measurements contained errors that contributed to the uncertainties observed in the experimental results. To conduct the uncertainty analysis, numerous methods are available ( [27][28][29]). The analysis followed here is the methodology outlined by Holman [30]. Consider a result R, determined by a specific function involving independent variables X 1 , X 2 , X 3 , . . . , X n . In that result, R is the function of independent variables (X 1 , X 2 , X 3 , . . . , X n ), where R = R (X 1 ; X 2 ; X 3 ; . . . ; X n ).
The uncertainty in the result denoted as W R , can be calculated using the equation provided by Holman [30]. Additionally, W 1 , W 2 , W 3 , . . . , W n represent the uncertainties associated with the independent variables X 1 , X 2 , X 3 , . . . , X n , respectively.
The following equation was used to calculate the values of W R: The relative error can be calculated using the following equation: The uncertainties and relative error equations were used to evaluate the temperature difference across the heater, air mass flow rate, thermal efficiency, useful heat gained, and the convective heat transfer coefficient of the SAH. The corresponding values are provided in Table 1.

Numerical Setup
The schematic layout of the proposed experimental setup with the proposed rib surface, as depicted in Figure 6, was simulated in the computational fluid dynamics (CFD) software. The boundary conditions included an inlet section (890 mm) where air entered the rectangular duct from the left side, a test section at 1000 mm, and an exit section at 1000 mm. The 2-D schematic diagram and the proposed rib surface are illustrated in Figure 6a,b, respectively. Based on the findings from the experimental work, a sudden increase in heat transfer was observed starting from Reynolds number (Re) 5000. To simulate this phenomenon numerically, we assumed a mass flow rate range of Re 6000 to Re 20,000, encompassing the Reynolds numbers where significant heat transfer enhancements were observed in the experimental results. Thermo-physical properties of the working fluid and absorber plate are presented in Table 2. Furthermore, a grid-independent study was conducted at four different element sizes to verify Nusselt number variation. It was found that 168,314 node and element size 0.22 mm produced the least difference compared with the earlier element size 0.24 mm, as shown in Figure 7a,b. Therefore, an element size of 0.22 mm was considered for further study.
The uncertainties and relative error equations were used to evaluate the temperature difference across the heater, air mass flow rate, thermal efficiency, useful heat gained, and the convective heat transfer coefficient of the SAH. The corresponding values are provided in Table 1.

Numerical Setup
The schematic layout of the proposed experimental setup with the proposed rib surface, as depicted in Figure 6, was simulated in the computational fluid dynamics (CFD) software. The boundary conditions included an inlet section (890 mm) where air entered the rectangular duct from the left side, a test section at 1000 mm, and an exit section at 1000 mm. The 2-D schematic diagram and the proposed rib surface are illustrated in Figure 6a,b, respectively. Based on the findings from the experimental work, a sudden increase in heat transfer was observed starting from Reynolds number (Re) 5000. To simulate this phenomenon numerically, we assumed a mass flow rate range of Re 6000 to Re 20,000, encompassing the Reynolds numbers where significant heat transfer enhancements were observed in the experimental results. Thermo-physical properties of the working fluid and absorber plate are presented in Table 2. Furthermore, a grid-independent study was conducted at four different element sizes to verify Nusselt number variation. It was found that 168,314 node and element size 0.22 mm produced the least difference compared with the earlier element size 0.24 mm, as shown in Figure 7a,b. Therefore, an element size of 0.22 mm was considered for further study.

Governing Equation
The continuity, momentum, and energy equations can be utilized to express the governing equations for steady-state, two-dimensional airflow in a solar air heater with artificial roughness. These are the equations: Mass conservation or continuity equation: ii. Momentum conservation or Navier-Stokes equation: iii. Energy conservation equation: Further analysis of non-dimensional flow parameters includes: iv. Reynolds number vi. Friction factor vii. Thermo-hydraulic performances parameter An empirical equation is applied to validate the smooth surface Nusselt number by: viii. Dittus-Boelter equation [24]:

Validation of Turbulent Model
The selection and validation of a turbulent model play crucial roles in accurately predicting the thermal coefficient augmentation for a given geometry. In order to investigate this, a literature survey was conducted. The Dittus-Boelter equation (Equation 10) was used for validation purposes [31,32]. For the specific case of a smooth duct in a solar air heater, the Nusselt number (Nu) was compared among various turbulent models, including: The detailed results of this comparison are presented in Figure 8. It was observed that the RNG renormalization k-ε model significantly disrupted the streamline flow and provided Nusselt number values closer to those predicted by the empirical equation. In contrast, the other models either overestimated or underestimated the Nusselt number. As a result, for simulating heat transfer in the current study, the RNG renormalization k-ε model was selected, in agreement with reference [23]. In addition, Figure 8 shows a new form of the correlation formula for the Nusselt number Nu, taking into account 3 different ranges of the Prandtl number Pr [33].

Validation of Turbulent Model
The selection and validation of a turbulent model play crucial roles in accurately predicting the thermal coefficient augmentation for a given geometry. In order to investigate this, a literature survey was conducted. The Dittus-Boelter equation (Equation 10) was used for validation purposes [31,32]. For the specific case of a smooth duct in a solar air heater, the Nusselt number (Nu) was compared among various turbulent models, including: The detailed results of this comparison are presented in Figure 8. It was observed that the RNG renormalization k-ε model significantly disrupted the streamline flow and provided Nusselt number values closer to those predicted by the empirical equation. In contrast, the other models either overestimated or underestimated the Nusselt number. As a result, for simulating heat transfer in the current study, the RNG renormalization k-ε model was selected, in agreement with reference [23]. In addition, Figure 8 shows a new form of the correlation formula for the Nusselt number Nu, taking into account 3 different ranges of the Prandtl number Pr [33].

Solution Procedure
The current simulation solved the governing equations of continuity, momentum, and energy using ANSYS-Fluent version 13.0. The geometry of suggested artificial roughness, grid generation study, and solution convergent process was completed in this portion. From the model validation, the RNG renormalization k-ε model was selected for the numerical analysis. The SIMPLE (semi-implicit methods for pressure linked equation) pressure-velocity model was chosen with a second-order upwind discretization scheme. To achieve better results, a residual monitor was fixed at 10 −6 in the x, y, T-k model. A

Solution Procedure
The current simulation solved the governing equations of continuity, momentum, and energy using ANSYS-Fluent version 13.0. The geometry of suggested artificial roughness, grid generation study, and solution convergent process was completed in this portion. From the model validation, the RNG renormalization k-ε model was selected for the numerical analysis. The SIMPLE (semi-implicit methods for pressure linked equation) pressure-velocity model was chosen with a second-order upwind discretization scheme. To achieve better results, a residual monitor was fixed at 10 −6 in the x, y, T-k model. A standard initialization was carried out in all zones. Finally, a convergent calculation was activated individually for all geometry to identify flow characterization and friction factors. When the solution convergence was completed, an average Nusselt number and average friction factor were calculated from the results. Figure 9 shows the Nusselt number and friction factor values from an experimental investigation of a smooth surface absorber plate, compared to the Dittus-Boelter and Blausis empirical equations (Equations (8) and (9)) for varying Reynolds numbers. The results of the surface-averaged heat transfer (SAH) analysis reveal that as the Reynolds number increases, the Nusselt number (Nu) demonstrates improved performance, indicating enhanced convective heat transfer. Conversely, the friction factor (f) decreases with higher Reynolds numbers, implying reduced flow resistance. Notably, at a Reynolds number of 20,000, the Nusselt number was found to be higher, while the friction factor was lower, reinforcing the observed trends. Figures 9 and 10 indicate that the experimental results are in good agreement with the empirical equations, providing confidence in the validity of the research and the possibility of further modification.

Experimental Results
standard initialization was carried out in all zones. Finally, a convergent calculation wa activated individually for all geometry to identify flow characterization and friction fac tors. When the solution convergence was completed, an average Nusselt number and av erage friction factor were calculated from the results. Figure 9 shows the Nusselt number and friction factor values from an experimenta investigation of a smooth surface absorber plate, compared to the Dittus-Boelter and Blau sis empirical equations (Equations (8) and (9)) for varying Reynolds numbers. The result of the surface-averaged heat transfer (SAH) analysis reveal that as the Reynolds numbe increases, the Nusselt number (Nu) demonstrates improved performance, indicating en hanced convective heat transfer. Conversely, the friction factor (f) decreases with highe Reynolds numbers, implying reduced flow resistance. Notably, at a Reynolds number o 20,000, the Nusselt number was found to be higher, while the friction factor was lower reinforcing the observed trends. Figures 9 and 10 indicate that the experimental result are in good agreement with the empirical equations, providing confidence in the validity of the research and the possibility of further modification.

Nusselt Number Characterization
The performance of the Nusselt number (Nu) for different relative pitch distances o polygonal rib, forward trapezoidal rib, and smooth absorber plates is presented in Tabl 3, revealing that the Nusselt number increased with an increased Reynolds number. Th standard initialization was carried out in all zones. Finally, a convergent calculation wa activated individually for all geometry to identify flow characterization and friction fac tors. When the solution convergence was completed, an average Nusselt number and av erage friction factor were calculated from the results. Figure 9 shows the Nusselt number and friction factor values from an experimenta investigation of a smooth surface absorber plate, compared to the Dittus-Boelter and Blau sis empirical equations (Equations (8) and (9)) for varying Reynolds numbers. The result of the surface-averaged heat transfer (SAH) analysis reveal that as the Reynolds numbe increases, the Nusselt number (Nu) demonstrates improved performance, indicating en hanced convective heat transfer. Conversely, the friction factor (f) decreases with highe Reynolds numbers, implying reduced flow resistance. Notably, at a Reynolds number o 20,000, the Nusselt number was found to be higher, while the friction factor was lower reinforcing the observed trends. Figures 9 and 10 indicate that the experimental result are in good agreement with the empirical equations, providing confidence in the validity of the research and the possibility of further modification.

Nusselt Number Characterization
The performance of the Nusselt number (Nu) for different relative pitch distances o polygonal rib, forward trapezoidal rib, and smooth absorber plates is presented in Tabl 3, revealing that the Nusselt number increased with an increased Reynolds number. Th

Nusselt Number Characterization
The performance of the Nusselt number (Nu) for different relative pitch distances of polygonal rib, forward trapezoidal rib, and smooth absorber plates is presented in Table 3, revealing that the Nusselt number increased with an increased Reynolds number. The implemented rib surfaces attained a higher Nusselt number than the smooth surface, as shown in Figure 11. It is evident that incorporating a rib surface in the solar air heater (SAH) improves the heat transfer coefficient. The study suggests that implementing polygonal ribs at some relative pitch distances (p/e) can enhance heat transfer and improve the SAH's efficiency more than others. Among them, the polygonal rib shape at a relative pitch distance (p/e) = 10 (p = 30 mm) provides the highest heat transfer rate due to the presence of sharp edges and an inclination on both sides of the rib, inducing significant turbulence in the fluid flow direction at appropriate pitch distances. implemented rib surfaces attained a higher Nusselt number than the smooth surface, a shown in Figure 11. It is evident that incorporating a rib surface in the solar air heate (SAH) improves the heat transfer coefficient. The study suggests that implementing po lygonal ribs at some relative pitch distances (p/e) can enhance heat transfer and improv the SAH's efficiency more than others. Among them, the polygonal rib shape at a relativ pitch distance (p/e) = 10 (p = 30 mm) provides the highest heat transfer rate due to th presence of sharp edges and an inclination on both sides of the rib, inducing significan turbulence in the fluid flow direction at appropriate pitch distances.

Friction Factor Characterization
In the study, Figure 12 illustrates the evaluation of friction factor for smooth and postulated surfaces (referring to the polygonal rib plate). It can be observed that the friction factor on the trapezoidal rib plate p/e = 13.4 is lower than that of the other surface at higher Reynolds numbers Re 20,000. This phenomenon can be attributed to the larger shear stresses occurring closer to the wall surface on the ribbed surface, which results in a greater pressure drop in the fluid flow direction. These findings suggest that incorporating trapezoidal rib plates can be advantageous in reducing the friction factor and slightly enhancing the system's overall efficiency compared with the polygonal ribbed surface. The detailed findings are presented in Table 4.  The study compared the thermal and hydraulic performance (THP) of a rib surfac with a smooth surface in a solar air heater. In the experimental investigation, the research ers were interested in understanding how the proposed surface roughness affects the hea transfer and pressure drop in the SAH compared with the smooth surface. Equation (7 was used to calculate the THP of the SAH, which is an essential factor in evaluating it effectiveness and optimizing its design. The analysis report is presented in Figure 13. I reveals that using polygonal ribs with p/e = 10 at Reynolds number Re = 4000 resulted i THP that was 2.95 times greater than that of a smooth surface. The detailed findings ar presented in Table 5. The use of polygonal ribs with a certain p/e ratio can be an effectiv way to enhance the performance of solar air heaters and improve their efficiency.

Thermo Hydraulic Performance (THP)
The study compared the thermal and hydraulic performance (THP) of a rib surface with a smooth surface in a solar air heater. In the experimental investigation, the researchers were interested in understanding how the proposed surface roughness affects the heat transfer and pressure drop in the SAH compared with the smooth surface. Equation (7) was used to calculate the THP of the SAH, which is an essential factor in evaluating its effectiveness and optimizing its design. The analysis report is presented in Figure 13. It reveals that using polygonal ribs with p/e = 10 at Reynolds number Re = 4000 resulted in THP that was 2.95 times greater than that of a smooth surface. The detailed findings are presented in Table 5. The use of polygonal ribs with a certain p/e ratio can be an effective way to enhance the performance of solar air heaters and improve their efficiency.    Figure 14 illustrates the heat transfer and fluid flow characteristics of a polygonal rib surface at various pitch distances and Reynolds numbers. The results indicate that lower Reynolds numbers produced lower Nusselt numbers, whereas higher Reynolds numbers produced higher Nusselt numbers for both experimental and numerical procedures. The experimental study showed lower Nusselt numbers at p/e = 6.7 (p = 20 mm) and p/e = 13.4 (p = 40 mm) for Re = 6000, and maximum Nusselt numbers at p/e = 10 (p = 30 mm) for higher Reynolds numbers (Re = 20,000) compared with the numerical procedure. Additionally, it was found that the minimum pitch distance of p/e = 6.7 yielded a higher heat transfer rate than p/e = 10 and p/e = 13.4. The velocity contour plot results in Figure 15 demonstrate the heat flow distribution at various pitch distances, where a strong vortex is generated near the first pair of ribs, followed by its distribution in the fluid flow direction. Moreover, increasing the Reynolds number leads to a larger vortex size generated inside the duct, resulting in higher heat transfer at higher Reynolds numbers and normal disturbances occurring in the fluid at lower Reynolds numbers for all pitch distances. Among them, the p/e = 6.7 pitch distance produced stronger turbulence in the fluid flow than other pitch distances.  Figure 14 illustrates the heat transfer and fluid flow characteristics of a polygonal rib surface at various pitch distances and Reynolds numbers. The results indicate that lower Reynolds numbers produced lower Nusselt numbers, whereas higher Reynolds numbers produced higher Nusselt numbers for both experimental and numerical procedures. The experimental study showed lower Nusselt numbers at p/e = 6.7 (p = 20 mm) and p/e = 13.4 (p = 40 mm) for Re = 6000, and maximum Nusselt numbers at p/e = 10 (p = 30 mm) for higher Reynolds numbers (Re = 20,000) compared with the numerical procedure. Additionally, it was found that the minimum pitch distance of p/e = 6.7 yielded a higher heat transfer rate than p/e = 10 and p/e = 13.4. The velocity contour plot results in Figure 15 demonstrate the heat flow distribution at various pitch distances, where a strong vortex is generated near the first pair of ribs, followed by its distribution in the fluid flow direction. Moreover, increasing the Reynolds number leads to a larger vortex size generated inside the duct, resulting in higher heat transfer at higher Reynolds numbers and normal disturbances occurring in the fluid at lower Reynolds numbers for all pitch distances. Among them, the p/e = 6.7 pitch distance produced stronger turbulence in the fluid flow than other pitch distances.   The heat dissipation between rib surfaces in the fluid flow direction was investigated using a turbulent kinematic energy (T-k) model, and the results are presented in Figure  16 for all pitch distances with Reynolds numbers ranging from Re 20,000 to Re 5000. The findings reveal that the rib surface with a closer pitch distance produced a higher heat transfer rate than the other surfaces as the Reynolds number increased. The Nusselt number values obtained from the investigation are presented in Table 6. Furthermore, the maximum heat dissipation between the wall surface and rib surface was clearly traced in the turbulence intensity, as shown in Figure 17. It was observed that turbulence occurred from the first pair of ribs and gradually increased with each subsequent pair of ribs, and the size of the vortex also increased in the fluid flow direction. Among them, the vortex with the highest heat transfer rate at a higher Reynolds number was recorded for p/e = 6.4. The heat dissipation between rib surfaces in the fluid flow direction was investigated using a turbulent kinematic energy (T-k) model, and the results are presented in Figure 16 for all pitch distances with Reynolds numbers ranging from Re 20,000 to Re 5000. The findings reveal that the rib surface with a closer pitch distance produced a higher heat transfer rate than the other surfaces as the Reynolds number increased. The Nusselt  Table 6. Furthermore, the maximum heat dissipation between the wall surface and rib surface was clearly traced in the turbulence intensity, as shown in Figure 17. It was observed that turbulence occurred from the first pair of ribs and gradually increased with each subsequent pair of ribs, and the size of the vortex also increased in the fluid flow direction. Among them, the vortex with the highest heat transfer rate at a higher Reynolds number was recorded for p/e = 6.4.
The heat dissipation between rib surfaces in the fluid flow direction was investigated using a turbulent kinematic energy (T-k) model, and the results are presented in Figure  16 for all pitch distances with Reynolds numbers ranging from Re 20,000 to Re 5000. The findings reveal that the rib surface with a closer pitch distance produced a higher heat transfer rate than the other surfaces as the Reynolds number increased. The Nusselt number values obtained from the investigation are presented in Table 6. Furthermore, the maximum heat dissipation between the wall surface and rib surface was clearly traced in the turbulence intensity, as shown in Figure 17. It was observed that turbulence occurred from the first pair of ribs and gradually increased with each subsequent pair of ribs, and the size of the vortex also increased in the fluid flow direction. Among them, the vortex with the highest heat transfer rate at a higher Reynolds number was recorded for p/e = 6.4.
T-k p = 6.7 T-k p = 10 T-k p = 13.6 Re 20,000 Re 20,000 Re 20,000 Re 5000 Re 5000 Re 5000   Re 5000 Re 5000 Re 5000 Figure 17. Contour plot of turbulence intensity at various Reynolds numbers Re = 5000 to Re 0 20,000 at p/e = 6.7, p/e = 10, p/e = 13.4. Figure 17. Contour plot of turbulence intensity at various Reynolds numbers Re = 5000 to Re 0 20,000 at p/e = 6.7, p/e = 10, p/e = 13.4. The pressure drop due to the rib surface at varying Reynolds numbers was recorded in all case studies in this investigation, as shown in Figure 18. The friction factor values in relation to pitch distance and varying Reynolds numbers are presented in Table 6. It was observed that the friction factor decreased with increasing Reynolds number and increased with decreasing Reynolds number. Additionally, the experimental results produced more intense pressure drops than the numerical study. The pitch distance of p = 30 mm (p/e = 10) resulted in the highest friction losses at lower Reynolds numbers of Re 5000. The minimum pressure drops were recorded for pitch distances of p = 20 mm (p/e = 6.7) in the experimental study and p = 40 mm (p/e = 13.7) in the numerical analysis. The investigation results are summarized in Table 7, which shows that the minimum friction factor was recorded for a closer pitch distance of p = 20 mm (p/e = 6.7), owing to the higher Reynolds number causing higher heat transfer rates than the other pitch distances. The pressure drop due to the rib surface at varying Reynolds numbers was recorded in all case studies in this investigation, as shown in Figure 18. The friction factor values i relation to pitch distance and varying Reynolds numbers are presented in Table 6. It wa observed that the friction factor decreased with increasing Reynolds number and in creased with decreasing Reynolds number. Additionally, the experimental results pro duced more intense pressure drops than the numerical study. The pitch distance of p = 3 mm (p/e = 10) resulted in the highest friction losses at lower Reynolds numbers of Re 5000 The minimum pressure drops were recorded for pitch distances of p = 20 mm (p/e = 6.7 in the experimental study and p = 40 mm (p/e = 13.7) in the numerical analysis. The inves tigation results are summarized in Table 7, which shows that the minimum friction facto was recorded for a closer pitch distance of p = 20 mm (p/e = 6.7), owing to the highe Reynolds number causing higher heat transfer rates than the other pitch distances.  The pressure drop between the polygonal rib surfaces was clearly observed and i presented in the pressure contour plot results shown in Figure 19. It was observed that a higher Reynolds numbers, a strong vortex was generated with the closer pitch distance o  The pressure drop between the polygonal rib surfaces was clearly observed and is presented in the pressure contour plot results shown in Figure 19. It was observed that at higher Reynolds numbers, a strong vortex was generated with the closer pitch distance of p/e = 6.7 for all Reynolds numbers ranging from 20,000 to 5000. Other relative pitch distances produced gradual vortices at all Reynolds numbers. It was found that lower friction occurred at relative pitch distances of p/e = 13.7 in numerical investigations and p/e = 10 in experimental investigations. Re 5000 Re 5000 Re 5000 Figure 19. Contour plot of pressure at various Reynolds numbers Re = 5000 to Re = 20,000 at p/e = 20 mm, p/e = 30 mm, p/e = 40 mm.

Correlations Developed for Nusselt Number and Friction Factor
The development of an empirical correlation based on experimental investigations is important for predicting data and extending the study for further work. In this study, a correlation for the Nusselt number was developed by examining the effects of various parameters such as (a) Nusselt number increase with increasing Reynolds number, Re and (b) Nusselt number increase with relative pitch distance, p/e. Developing an accurate empirical correlation can be an essential tool for predicting the performance of the solar air heater and optimizing its design for maximum efficiency.
A second-order polynomial equation was used in straight-line regression analysis for the proposed parameters Nusselt number, Reynolds number, and relative pitch distance, and a correlation was developed [24]: Furthermore, another correlation was developed from the observed effects of the parameters on the friction factor: a. Friction factor decreases with increasing Reynolds number Re. b. Friction factor increases with relative pitch distance, p/e. This correlation can be used to predict the friction factor for a given combination of Reynolds number and relative pitch distance in a solar air heater. It is a valuable tool for optimizing the design and operation of the system, as it allows for the prediction of pressure drop under different conditions without the need for extensive experimentation [24]. ƒ = 4.06·10 −6 Re −0.225 (p/e)·exp [−0.0117 Ln (p/e) 2 ] (13)

Conclusions
The experimental and 2D numerical analysis focused on investigating the heat transfer and friction factor of polygonal and trapezoidal ribs with varying pitch distances in an SAH. The calibrated experimental and numerical results showed an enhancement of the average Nusselt number (Nu) with an increasing Reynolds number, and a decrease in average friction factor (ƒ) with an increasing Reynolds number. Based on the experimental and numerical analysis, it was found that: a. Implementing artificial roughness in the absorber plate significantly improved the Nusselt number (Nu) compared to a smooth surface, at all Reynolds numbers. b. The RNG renormalized group k-ε model predicts very close results (±3%) to the Dittus-Boelter empirical values in the numerical and experimental work. c. The polygonal-shaped rib, with p = 20 mm (numerical) and p = 30 mm (experimental), recorded the highest average Nusselt number at Reynolds number 20,000 compared with other rib shapes, indicating that it can effectively enhance heat transfer in the SAH.

Correlations Developed for Nusselt Number and Friction Factor
The development of an empirical correlation based on experimental investigations is important for predicting data and extending the study for further work. In this study, a correlation for the Nusselt number was developed by examining the effects of various parameters such as (a) Nusselt number increase with increasing Reynolds number, Re and (b) Nusselt number increase with relative pitch distance, p/e. Developing an accurate empirical correlation can be an essential tool for predicting the performance of the solar air heater and optimizing its design for maximum efficiency.
A second-order polynomial equation was used in straight-line regression analysis for the proposed parameters Nusselt number, Reynolds number, and relative pitch distance, and a correlation was developed [24]: Furthermore, another correlation was developed from the observed effects of the parameters on the friction factor: a.
Friction factor decreases with increasing Reynolds number Re. b.
Friction factor increases with relative pitch distance, p/e. This correlation can be used to predict the friction factor for a given combination of Reynolds number and relative pitch distance in a solar air heater. It is a valuable tool for optimizing the design and operation of the system, as it allows for the prediction of pressure drop under different conditions without the need for extensive experimentation [24]. ƒ = 4.06·10 −6 Re −0.225 (p/e)·exp [−0.0117 Ln (p/e) 2 ] (13)

Conclusions
The experimental and 2D numerical analysis focused on investigating the heat transfer and friction factor of polygonal and trapezoidal ribs with varying pitch distances in an SAH. The calibrated experimental and numerical results showed an enhancement of the average Nusselt number (Nu) with an increasing Reynolds number, and a decrease in average friction factor (ƒ) with an increasing Reynolds number. Based on the experimental and numerical analysis, it was found that: a.
Implementing artificial roughness in the absorber plate significantly improved the Nusselt number (Nu) compared to a smooth surface, at all Reynolds numbers. b.
The RNG renormalized group k-ε model predicts very close results (±3%) to the Dittus-Boelter empirical values in the numerical and experimental work. c.
The polygonal-shaped rib, with p = 20 mm (numerical) and p = 30 mm (experimental), recorded the highest average Nusselt number at Reynolds number 20,000 compared with other rib shapes, indicating that it can effectively enhance heat transfer in the SAH. d.
The minimum friction factor (f) of 0.01 at pitch distance of p = 20 in experimental work and 0.0093 at p = 40 mm in numerical work was attained at Reynolds number Re 20,000.