Next Article in Journal
Scientific Approaches to Solving the Problem of Joint Processes of Bubble Boiling of Refrigerant and Its Movement in a Heat Pump Heat Exchanger
Previous Article in Journal
Short-Term Load Forecasting Based on Outlier Correction, Decomposition, and Ensemble Reinforcement Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of Thermo-Hydraulic Performances of Artificial Ribs Mounted in a Rectangular Duct

1
Velammal College of Engineering and Technology, Madurai 625107, India
2
CO2 Research and Green Technologies Centre, Vellore Institute of Technology, Vellore 632014, India
3
Department of Thermal Processes, Air Protection and Waste Management, Cracow University of Technology, 31-155 Cracow, Poland
4
Department of Energy, Cracow University of Technology, 31-864 Cracow, Poland
*
Author to whom correspondence should be addressed.
Energies 2023, 16(11), 4404; https://doi.org/10.3390/en16114404
Submission received: 27 April 2023 / Revised: 26 May 2023 / Accepted: 27 May 2023 / Published: 30 May 2023
(This article belongs to the Section J1: Heat and Mass Transfer)

Abstract

:
This research investigates the fluid flow characterization and thermohydraulic performances (THP) of rib surfaces, using computational and experimental methods. ANSYS computational fluid dynamics (CFD) software was used to predict and validate the findings in the experimental setup. Artificial 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 flow rate of air (working fluid) 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 thermohydraulic 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.

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

2. 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,
(iii)
digital thermocouples,
(iv)
a U-tube manometer,
(v)
an inclined manometer.
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.

2.1. 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  W H and the exit Section 2.5  W H .

2.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 nickel-chromium heated plate standing on it to provide heat energy equivalent to solar radiation at a constant heat flux of 1000 W/m2.
Figure 4 depicts the placement of electric resistances with 60-volt energy. A voltage regulator manages and regulates the energy supply to the heating plate. For heat protection, glass wool and thermo-form of 4 cm thickness were employed to cover the entire experimental apparatus.

2.3. 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/m3 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.

2.4. 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/m2. 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.

2.5. 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 X1, X2, X3,…, Xn. In that result, R is the function of independent variables (X1, X2, X3,…, Xn), where R = R (X1; X2; X3;…; Xn).
The uncertainty in the result denoted as WR, can be calculated using the equation provided by Holman [30]. Additionally, W1, W2, W3,…, Wn represent the uncertainties associated with the independent variables X1, X2, X3,…, Xn, respectively.
The following equation was used to calculate the values of WR:
W R = W 1 R X 1 2 + W 2 R X 21 2 + . . . + W n R X n 2 + 0.5
The relative error can be calculated using the following equation:
ER = WR/R%
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.

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

3.1. 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:
i.
Mass conservation or continuity equation:
x i ρ u i = 0
ii.
Momentum conservation or Navier–Stokes equation:
x i ρ u i u j = ρ x i + μ x j u i x j + u j x i
iii.
Energy conservation equation:
ρ u j T x i = x j Γ + Γ t T x j
Further analysis of non-dimensional flow parameters includes:
iv.
Reynolds number
Re = ρ u D μ
v.
Nusselt number
Nu = h D ν
vi.
Friction factor
fr = p / l D 2 ρ v 2
vii.
Thermo-hydraulic performances parameter
THP = N u r / N u s f r / f s 1 / 3
An empirical equation is applied to validate the smooth surface Nusselt number by:
viii.
Dittus–Boelter equation [24]:
Nu = 0.023 Re0.8 Pr0.4
ix.
Blasius equation:
f = 0.00791 Re−0.25

3.2. 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:
(a)
The standard k-ε model,
(b)
The renormalization k-ε model,
(c)
The realizable k-ε model,
(d)
The standard k-ω model,
(e)
The SST k-ω model.
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].

3.3. 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 106 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.

4. Result and Discussion

4.1. Experimental 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. Figure 9 and Figure 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.

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

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

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

4.2. Numerical Results

4.2.1. Nusselt Number

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.

4.2.2. Friction Factor

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

5. 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]:
Nu = 0.0251 Re0.833 (p/e)0.218·exp [0.01318 Ln (p/e)2]
Furthermore, another correlation was developed from the observed effects of the parameters on the friction factor:
  • Friction factor decreases with increasing Reynolds number Re.
  • 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]

6. 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:
  • Implementing artificial roughness in the absorber plate significantly improved the Nusselt number (Nu) compared to a smooth surface, at all Reynolds numbers.
  • The RNG renormalized group k-ε model predicts very close results (±3%) to the Dittus–Boelter empirical values in the numerical and experimental work.
  • 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.
  • 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.
  • The thermo-hydraulic performance (THP) of the proposed shapes achieved the highest values of 2.95 in the polygonal rib at p/e = 10 and Reynolds number 4000.
  • An empirical correlation was developed for Nusselt number Nu and friction factor ƒ from the experimental values in (Equation (10) and (11)), which showed maximum accuracy of ±8% between experimental and calculated values.
Overall, the use of artificial roughness in SAHs can be an effective way to enhance heat transfer and improve system efficiency. The polygonal shape of the rib with p/e = 6.7/p = 20 mm (numerical) and p/e = 10/p = 30 mm (experimental) showed the most significant improvement in the Nusselt number, making it a promising choice for future SAH design.

Author Contributions

Conceptualization, methodology, investigation: B.V.K., P.R.K., Software, validation writing—review and editing: P.R.K., G.M., J.T. Formal analysis, resources, data curation: D.T., M.N.-O. Writing—original draft preparation, supervision: B.V.K., J.T., T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data available on request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Nomenclature

Across-sectional area of the duct, [A = WH], m2
Acsurface area of collector plate, m2
Cdcoefficient of discharge for orifice meter
Cpspecific heat of air at constant pressure, kJ/(kg K)
D h, dequivalent diameter of the air passage, [D = 4A/[2(W + H)], m
Iturbulence intensity/intensity of solar radiation, W/m2
eroughness height, mm
e/drelative roughness height
Hheight of air channel, m
hheat transfer coefficient, W m−2 K−1
kthermal conductivity of air, W m−1 K−1
Lduct length, m
L1inlet length of duct, mm
L2test length of duct, mm
L3outlet length of duct, mm
mmass flow rate of air, kg/s
proughness pitch, m
Dimensionless parameters
ffriction factor
frfriction factor for rough surface
fsfriction factor for a smooth surface
fexpfriction factor experimental
NuNusselt number
NusNusselt number for rough duct
NurNusselt number for smooth duct
Nu expNusselt number experimental
PrPrandtl number
P/erelative roughness pitch
ReReynolds number
W/Hduct aspect ratio
Greek symbols
µdynamic viscosity, Ns/m2
µtturbulent viscosity, Ns/m2
ρdensity of air, kg/m3
αangle of attack, degree
εdissipation rate
ωspecific dissipation rate
δtransition sub-layer thickness, m
kturbulence kinetic energy, m
Γmolecular thermal diffusivity
Γtturbulent thermal diffusivity

References

  1. Saha, S.K.; Islam, M.S.; Saha, S.K.; Alam, M.M. Experimental study of heat transfer and friction characteristics of rectangular duct with ribbed walls. Int. J. Heat Mass Transf. 2013, 60, 559–567. [Google Scholar]
  2. Kim, K.H.; Kim, K.Y.; Ha, M.Y. Heat transfer and pressure drop characteristics of a rectangular duct with perforated ribs. Int. J. Heat Mass Transf. 2014, 71, 283–291. [Google Scholar]
  3. Kefayati, G.R.; Safaei, M.R.; Esfahani, J.A. Numerical study on the heat transfer and pressure drop in a rectangular duct with longitudinal rib turbulators. Appl. Therm. Eng. 2018, 138, 621–632. [Google Scholar]
  4. Ji, H.; Luo, L.; Yan, X.; Zhang, B. Numerical investigation on heat transfer and fluid flow characteristics in rectangular ducts with twisted tape inserts. Appl. Therm. Eng. 2020, 169, 114919. [Google Scholar]
  5. Balakrishnan, A.R.; Khader, M.A. Numerical investigation of heat transfer augmentation in a rectangular duct using discrete multiple V-ribs. Int. J. Heat Mass Transf. 2015, 91, 421–430. [Google Scholar]
  6. Narayanan, R.; Bhaskaran, R. Numerical analysis of turbulent flow and heat transfer in a rectangular duct with multiple ribs of varying geometry. J. Therm. Anal. Calorim. 2019, 137, 431–446. [Google Scholar]
  7. Li, H.; Li, J.; Li, J.; Li, Q. Experimental study on turbulent heat transfer enhancement in a rectangular duct with staggered sawtooth ribs. Exp. Therm. Fluid Sci. 2021, 124, 111210. [Google Scholar]
  8. Wen, D.; He, Y.L.; Tang, G.H. Numerical investigation of heat transfer and flow characteristics in a rectangular duct with elliptical dimples. Int. J. Heat Mass Transf. 2019, 132, 9–19. [Google Scholar]
  9. Wang, W.; He, Y.L.; Tao, W.Q. Numerical investigation of heat transfer and fluid flow characteristics in a rectangular duct with multiple dimple-shaped ribs. Int. J. Heat Mass Transf. 2021, 181, 121458. [Google Scholar]
  10. Karwa, R.; Solanki, S.C.; Saini, J.S. Thermo-hydraulic performance of solar air heaters having integral chamfered rib roughness on absorber plates. Energy 2001, 26, 161–176. [Google Scholar] [CrossRef]
  11. Ahmed, A.E.; Ahmed, N.M.; Ahmed, M.S. Experimental and numerical investigation of heat transfer enhancement in a rectangular duct using wavy and perforated ribs. J. Therm. Anal. Calorim. 2020, 139, 465–483. [Google Scholar]
  12. Sahiti, N.; Cao, B.Y. Experimental and numerical investigation of heat transfer enhancement in a rectangular duct with a novel serrated twisted tape. Exp. Therm. Fluid Sci. 2017, 81, 77–89. [Google Scholar]
  13. Siddiqui, M.F.; Kim, M.H. Experimental and numerical investigation of flow and heat transfer in a rectangular duct with multiple V-shaped ribs. Int. J. Heat Mass Transf. 2018, 117, 689–702. [Google Scholar]
  14. Islam, M.S.; Alim, M.A. Numerical study on thermal-hydraulic characteristics of turbulent flow in a rectangular duct with novel converging/diverging rib arrangements. Appl. Therm. Eng. 2021, 188, 116680. [Google Scholar]
  15. Lee, H.J.; Kang, B.H. Heat transfer and fluid flow characteristics in a rectangular duct with multi-vortex generator. Int. J. Heat Mass Transf. 2020, 160, 120221. [Google Scholar]
  16. Seo, J.; Kim, Y. Experimental study of heat transfer and friction factor characteristics in a rectangular duct with novel compound turbulator. Appl. Therm. Eng. 2021, 194, 113098. [Google Scholar]
  17. Karthikeyan, K.; Kannan, S. Numerical study of fluid flow and heat transfer in a rectangular duct with discrete vortex generator. Appl. Therm. Eng. 2021, 187, 116438. [Google Scholar]
  18. Verma, S.; Prasad, B. Investigation for the optimal thermo-hydraulic performance of artificially roughened solar air heaters, Renewable. Energy 2000, 20, 19–36. [Google Scholar]
  19. Saini, R.P.; Saini, S.K. Development of correlations for Nusselt number and friction factor for solar air heater with roughened duct having arc-shaped wire as artificial roughness. Sol. Energy 2008, 82, 1118–1130. [Google Scholar] [CrossRef]
  20. Bopche, S.B.; Tandale, M.S. Experimental investigations on heat transfer and frictional characteristics of a turbulator roughened solar air heater duct. Int. Commun. Heat Mass Transf. 2009, 52, 2834–2848. [Google Scholar] [CrossRef]
  21. Saini, R.P.; Verma, J. Heat transfer and friction factor correlations for a duct having dimple shape artificial roughness for solar air heaters. Energy 2008, 33, 1277–1287. [Google Scholar] [CrossRef]
  22. Kumar, A.; Bhagoria, J.L.; Sarviya, R.M. Heat transfer and friction correlations for artificially roughened solar air heater duct with discrete W-shaped ribs. Energy Conversat. Manag. 2009, 50, 2106–2117. [Google Scholar] [CrossRef]
  23. Kumar, B.V.; Manikandan, G.; Kanna, P.R. Enhancement of heat transfer in SAH with polygonal and trapezoidal shape of the rib using CFD. Energy 2021, 234, 121154. [Google Scholar] [CrossRef]
  24. Kumar, B.V.; Manikandan, G.; Kanna, P.R. Performances of solar air heater using polygonal and trapezoidal rib absorber plate for augmentation of heat transfer. Tierärztliche Prax. 2020, 40, 888–903. [Google Scholar]
  25. Varun Kumar, B.; Manikandan, G.; Kanna, P.R.; Taler, D.; Taler, J.; Nowak-Ocłoń, M.; Mzyk, K.; Toh, H.T. A Performance Evaluation of a Solar Air Heater Using Different Shaped Ribs Mounted on the Absorber Plate—A Review. Energies 2018, 11, 3104. [Google Scholar] [CrossRef]
  26. ASHRAE Standard 93; Method of Testing to Determine the Thermal Performance of Solar Collectors. 30329. American Society of Heating, Refrigeration and Air Conditioning Engineers: Atlanta, GA, USA, 2003.
  27. Holman, J.P. Experimental Methods for Engineers; Tata McGraw-Hill Publishing Company Limited: New Delhi, India, 2007. [Google Scholar]
  28. Joint Committee for Guides in Metrology. JCGM 100: Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement; JCGM. 2008. (Also Known as GUM). Available online: https://www.bipm.org/documents/20126/2071204/JCGM_100_2008_E.pdf/cb0ef43f-baa5-11cf-3f85-4dcd86f77bd6 (accessed on 26 May 2023).
  29. Kline, S.J.; McClintock, F.A. Policy on reporting uncertainties in experimental measurements and results. ASME J. Heat Transf. 2000, 122, 411–413. [Google Scholar]
  30. Duffie, J.A.; Beckman, W.A. Solar Engineering of Thermal Processes; Wiley: New York, NY, USA, 1980. [Google Scholar]
  31. McAdams, W.H. Heat Transmission; McGraw-Hil: New York, NY, USA, 1942. [Google Scholar]
  32. Webb, R.L.; Eckert, E.R.G. Application of rough surface to heat exchanger design. Int. J. Heat Mass Transf. 1972, 15, 1647–1658. [Google Scholar] [CrossRef]
  33. Taler, D.; Taler, J. Simple heat transfer correlations for turbulent tube flow. E3S Web Conf. 2017, 13, 02008. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of Solar Air Heater (SAH) [24].
Figure 1. Schematic diagram of Solar Air Heater (SAH) [24].
Energies 16 04404 g001
Figure 2. Pictorial view of the experimental setup with the rib surface [24]. (a) Side view; (b) front view and shape of rib surface; (c) top; (d) front view.
Figure 2. Pictorial view of the experimental setup with the rib surface [24]. (a) Side view; (b) front view and shape of rib surface; (c) top; (d) front view.
Energies 16 04404 g002
Figure 3. Proposed rib surfaces (a) Trapezoidal rib, (b) Polygonal rib at various pitch distances [24].
Figure 3. Proposed rib surfaces (a) Trapezoidal rib, (b) Polygonal rib at various pitch distances [24].
Energies 16 04404 g003
Figure 4. Electric heater plate position in the test section [24].
Figure 4. Electric heater plate position in the test section [24].
Energies 16 04404 g004
Figure 5. Location of thermocouple in absorber plate [24].
Figure 5. Location of thermocouple in absorber plate [24].
Energies 16 04404 g005
Figure 6. Schematic layout of the (a) rectangular duct and (b) zoomed polygonal rib at p = 20 mm.
Figure 6. Schematic layout of the (a) rectangular duct and (b) zoomed polygonal rib at p = 20 mm.
Energies 16 04404 g006
Figure 7. Grid-independent study of proposed geometry at various element cells (a,b).
Figure 7. Grid-independent study of proposed geometry at various element cells (a,b).
Energies 16 04404 g007
Figure 8. Comparison of Nu for a turbulent model with empirical Equation (8).
Figure 8. Comparison of Nu for a turbulent model with empirical Equation (8).
Energies 16 04404 g008
Figure 9. Comparison of experimental values with a Dittus–Boelter equation (Nu vs. Re) [24].
Figure 9. Comparison of experimental values with a Dittus–Boelter equation (Nu vs. Re) [24].
Energies 16 04404 g009
Figure 10. Comparison of experimental values with a Blasius equation (f vs. Re) [24].
Figure 10. Comparison of experimental values with a Blasius equation (f vs. Re) [24].
Energies 16 04404 g010
Figure 11. Nusselt number vs. Reynolds number for smooth and proposed roughness shapes.
Figure 11. Nusselt number vs. Reynolds number for smooth and proposed roughness shapes.
Energies 16 04404 g011
Figure 12. Friction factor vs. Reynolds number for smooth shape and proposed roughness.
Figure 12. Friction factor vs. Reynolds number for smooth shape and proposed roughness.
Energies 16 04404 g012
Figure 13. Thermo hydraulic performances of proposed ribs with various Reynolds numbers.
Figure 13. Thermo hydraulic performances of proposed ribs with various Reynolds numbers.
Energies 16 04404 g013
Figure 14. Numerical and experimental Nu number vs. Re number at various pitch distances.
Figure 14. Numerical and experimental Nu number vs. Re number at various pitch distances.
Energies 16 04404 g014
Figure 15. Contour plot of velocity at various Reynolds numbers Re = 5000 to Re = 20,000 at p/e = 6.7, p/e = 10, p/e = 13.4.
Figure 15. Contour plot of velocity at various Reynolds numbers Re = 5000 to Re = 20,000 at p/e = 6.7, p/e = 10, p/e = 13.4.
Energies 16 04404 g015aEnergies 16 04404 g015b
Figure 16. Contour velocity plot at various Reynolds numbers Re = 5000 to Re = 20,000 at p/e = 6.7, p/e = 10, p/e = 13.4.
Figure 16. Contour velocity plot at various Reynolds numbers Re = 5000 to Re = 20,000 at p/e = 6.7, p/e = 10, p/e = 13.4.
Energies 16 04404 g016aEnergies 16 04404 g016b
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.
Energies 16 04404 g017
Figure 18. Numerical and experimental f vs. Re number at various pitch distances.
Figure 18. Numerical and experimental f vs. Re number at various pitch distances.
Energies 16 04404 g018
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.
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.
Energies 16 04404 g019aEnergies 16 04404 g019b
Table 1. Uncertainties and relative error values.
Table 1. Uncertainties and relative error values.
S. NoParameterAttained Values in SAHUncertaintiesRelative Error (%)
1Temperature differences (°C)16±0.170.67
2Heat flux (W/m2)1000±70.71
3Thermal efficiency (%)77.9±0.020.027
4Connective heat transfer (W/m2 °C)47.9±0.711.71
5Useful heat gained by SAH (W)872.8±8.11.57
Table 2. Thermo-Physical Properties of Air and Absorber plate (Aluminium) for analysis [23].
Table 2. Thermo-Physical Properties of Air and Absorber plate (Aluminium) for analysis [23].
S. NoPropertyAirAluminum
1Density, ρ (kg m−3)1.1172719
2Specific heat, Cp (J kg−1 K−1)1007871
3Thermal conductivity k (W m−1 K−1)0.0262202.4
4Viscosity µ (N m−2)1.857 × 10−5-
5Prandl number, Pr.0.71-
Table 3. Experimental values of Nusselt number vs. Reynolds number.
Table 3. Experimental values of Nusselt number vs. Reynolds number.
S. NoRe NoSmoothPolygonal p/e = 6.7Polygonal p/e = 10Polygonal p/e = 13.4Trapezoidal
p/e = 6.7
Trapezoidal
p/e = 10
Trapezoidal
p/e = 13.4
1200010.116.120.519.215.218.517.9
2400014.238.448.24636.747.545.9
3600023.260.770.86957.668.167
4800026.570.980.179.467.177.675.7
510,00032.276.986.48570.881.680.2
612,00039.188.998.697.885.595.794.2
714,00044.5107.2117.5115103.2113.9112.5
817,00052.7115.0125.9123.5109.7120.1119
920,00059.1134.5145.8144.1124.1139.9137
Table 4. Experimental values of friction factor vs. Reynolds number.
Table 4. Experimental values of friction factor vs. Reynolds number.
S. NoRe NoSmoothPolygonal
p/e = 6.7
Polygonal
p/e = 10
Polygonal
p/e = 13.4
Trapezoidal
p/e = 6.7
Trapezoidal
p/e = 10
Trapezoidal
p/e = 13.4
120000.01610.0180.0210.020.01740.020.019
240000.01240.01710.0190.0180.01670.01850.0179
360000.01030.01670.0180.01750.01610.01710.0165
480000.00960.01610.01740.01690.01530.01630.0157
510,0000.009010.01450.01650.01590.01410.01550.015
612,0000.00870.01340.01550.0150.01290.01410.0135
714,0000.008410.01280.01450.01370.01220.01340.013
817,0000.00790.01240.01340.01250.01190.01250.0115
920,0000.00710.01140.01270.0120.0110.01150.01
Table 5. Thermo hydraulic performances of proposed ribs.
Table 5. Thermo hydraulic performances of proposed ribs.
S. NoShapep/e200040006000800010,00012,00014,00017,00020,000
1Polygonal6.71.542.432.232.262.041.972.101.881.95
2101.862.952.542.482.202.082.212.012.04
313.41.772.862.502.492.192.092.202.012.05
4Trapezoidal6.71.472.342.142.171.901.922.051.821.82
5101.712.932.482.462.122.092.191.962.02
613.41.682.862.472.432.112.082.191.992.07
Table 6. Experimental and numerical values of Nusselt number vs. Reynolds number.
Table 6. Experimental and numerical values of Nusselt number vs. Reynolds number.
S. NoReNumericalExperimental
p = 20p = 30p = 40p = 20p = 30p = 40
120,000128123.5116.8134.5145.8144.1
217,000110.3114.5108.2115125.9123.5
314,000107.8103.798.2107.2117.5115
412,00094.892.287.188.998.697.8
5800078.2375.972.170.980.179.4
6600064.0163.660.860.770.869
Table 7. Experimental and numerical values of friction factor vs. Reynolds number.
Table 7. Experimental and numerical values of friction factor vs. Reynolds number.
S. NoReNumericalExperimental
p = 20p = 30p = 40p = 20p = 30p = 40
120,0000.00950.00940.00930.01140.01270.012
217,0000.01130.01110.010.01240.01340.0125
314,0000.01250.0120.0110.01280.01450.0137
412,0000.0140.01390.01380.01340.01550.015
580000.01620.01610.0160.01610.01740.0159
660000.01790.01780.01770.01670.0180.0161
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kumar, B.V.; Kanna, P.R.; Manikandan, G.; Taler, D.; Taler, J.; Sobota, T.; Nowak-Ocłoń, M. Investigation of Thermo-Hydraulic Performances of Artificial Ribs Mounted in a Rectangular Duct. Energies 2023, 16, 4404. https://doi.org/10.3390/en16114404

AMA Style

Kumar BV, Kanna PR, Manikandan G, Taler D, Taler J, Sobota T, Nowak-Ocłoń M. Investigation of Thermo-Hydraulic Performances of Artificial Ribs Mounted in a Rectangular Duct. Energies. 2023; 16(11):4404. https://doi.org/10.3390/en16114404

Chicago/Turabian Style

Kumar, B. Varun, P. Rajesh Kanna, G. Manikandan, Dawid Taler, Jan Taler, Tomasz Sobota, and Marzena Nowak-Ocłoń. 2023. "Investigation of Thermo-Hydraulic Performances of Artificial Ribs Mounted in a Rectangular Duct" Energies 16, no. 11: 4404. https://doi.org/10.3390/en16114404

APA Style

Kumar, B. V., Kanna, P. R., Manikandan, G., Taler, D., Taler, J., Sobota, T., & Nowak-Ocłoń, M. (2023). Investigation of Thermo-Hydraulic Performances of Artificial Ribs Mounted in a Rectangular Duct. Energies, 16(11), 4404. https://doi.org/10.3390/en16114404

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

Article Metrics

Back to TopTop