1. Introduction
Increasing renewable energy demands are leading to the enhancement of the performance of solar thermal systems, which are applied in various domains [
1]. A lower thermal efficiency in terms of the solar thermal collector (
) due to an inferior convective heat transfer coefficient is a major obstacle to harvesting thermal energy. The deficiency can be overcome by unsettling the boundary layer formed on the heat transfer surface by creating artificial roughness of different forms such as transverse and angled ribs, V, W and Z-shaped ribs, a dimple and protrusion-shaped roughness, baffles and blockages [
2]. The impinging jet is reported to be an effective method for heat transfer enhancement [
3]. The integration of an artificially roughened heated plate with the impinging of air to enhance the convective heat transfer coefficient [
4,
5] is an effective technique when used to harvest good amount of thermal energy. A lower thermal efficiency in terms of the solar thermal collector (
) due to inferior convective heat transfer coefficient is a major obstacle to harvesting thermal energy.
Nadda et al. [
6] analyzed a heated surface roughened with a protrusion in an arc shape by impinging the air jets on a solar air heater (SAH) absorber plate. A maximum thermal hydrodynamic performance of 1.5 was achieved. Kercher and Tabakoff [
7] analyzed streamwise and spanwise pitches in a low velocity jet impingement process. Nadda et al. [
8] analyzed an SAH duct for thermal performance by impinging air jets on the heated absorber plate. The investigation reported that the heat transfer enhancement was significant, but the corresponding pressure drop penalty decreased the overall thermohydraulic performance of the SAH, as determined by the preference selection index (
) methodology. Metzger et al. [
9] studied the effect of the staggered impingement of air jets on a heated surface. The results reported that the incorporation of jets in staggered arrangements are not better than the in-line arrangement of jets in terms of thermal performance for the selected parameters. Chauhan et al. [
10] analyzed the effect of employing the air jets impinging on a heated absorber plate surface with an aim of analyzing the thermal performance of the
duct. The various geometric parameters considered were the pitches in the flow direction and transverse direction, the diameter of the jet and the flow (
). The outcome of the analysis revealed a substantial improvement in the
ratio of 2.67. Brevet et al. [
11] investigated the array of air jet impinging on a heated surface for a range of jet impingement distances and spanwise jet-to-jet spacing. The air leaving the heated section after impingement was directed in a particular direction. The analysis resulted in the determination of the parametric values for spanwise jet-to-jet spacing and impingement distance to be 3 to 6 and 4 to 5, respectively. Nadda et al. [
12] experimentally analyzed an
duct for its heat transfer and friction characteristics by impinging circular air jets on the heated surface. The analysis was carried out by selecting various geometric parameters, i.e., relative height ratio, relative width ratio, relative pitch ratio and the angle of the arc. The
ratios, when compared with the smooth duct of an
, were found to be 6.29 and 9.25, respectively, and the extreme thermohydraulic performance parameter achieved was 3.64. Mishra et al. [
13] experimentally studied the effect of various parameters on a solar air passage with air jet impingement. An analytical hierarchy preference selection process using the (AHP-TOPSIS) technique was used to determine the set of geometric parameters that deliver the highest thermal performance with minimum friction losses. Sedighi et al. [
14] analyzed the cooling effect by impinging four air jets in a turbulent regime on a heated plate with the parameters selected as follows: the distance between the jets and the heated plate, the distance between the jets, the angle of the jets and the distance between the outer and inner jets. The results reported that the
is enhanced when the impinging jet angles are fixed at a higher
, when the jet-to-jet and jet-to-heated surface spaces are smaller and when there is a broader opening in terms of the outer jets. The swirling air jet impinging on the heated
absorber plate was analysed numerically by Afroz and Sharif [
15], who reported an 8% enhancement in the
by incorporating a swirling motion, as compared to a non-swirling motion.
The analysis of thermal performance using jet impingement is supported by the numerical techniques of various investigators. Issac et al. [
16] experimentally and numerically analyzed the performance of a round jet for various nozzle to plate spacing. The numerical analysis was carried out using a different RANS turbulence model. The results reported that the magnitudes of the inlet turbulent intensity and eddy viscosity play a major role in performance enhancement. N. Celik [
17] applied the design of experiment (DoE) and analysis of variance (
) to an experimental investigation of impinging jets. The major factors considered for investigation were the roughness of the heated surface and the jet geometry. It was reported that the
highly was affected by the radial distance, whereas the effect in terms of surface roughness was much lower. Singh et al. [
18] examined the influence of air jet impinging on a rectangular duct of double pass
with a corrugated wavy shape jet plate. A bed porosity of 98% and impinging plate perforation of 0.48% yielded a maximum thermal efficiency of 94%. Kumar et al. [
19] investigated an
duct with air jet impinging through a circular inside conical ring and found that the maximum thermos-hydraulic performance parameter was 2.16. Different configuration in terms of inner conical ring obstacles were analyzed by Kumar et al. [
20] using an
rectangular duct with impinging air jets and an of
from 5000 to 23,000. The maximum rate of heat transfer obtained was of the order of 2.33. Statistical correlations using the regression analysis were developed for
. Erasmus et al. [
21] experimentally and numerically investigated the impinging of an air jet on a concave hemispherical surface using computational tools. The correlation in terms of Nusselt number and the total pressure loss was developed in terms of Reynolds number and the Prandtl number.
The above literature study reveals that impinging jets are a significant heat enhancement technique. Artificial roughness generates a secondary flow because of the generation of vortices, which enhances the convective heat transfer. A staggered fluid jet impinging on the heated surface with a staggered roughness will boost the heat transfer, and it was this objective that the present study investigated. The parameters selected for investigation are novel, as there is no literature review present for the analysis of staggered jets on a staggered spherical protrusion location. The jets were impinging on the heated surface at the location where the surface was roughened in form of spherical protrusions. The location of the jet was exactly below the hemispherical protrusions, and the jet and protrusion diameters were the same. The location of the perpendicular striking staggered jets was exactly below the spherical protrusions to diverge the span of the jet to a larger area. A schematic of the jet position on the plate and the protruded plate is shown in
Figure 1. The range of flow and geometrical parameters selected for the experimental investigation were: The relative jet diameter to hydraulic diameter ratio (
) ranged from 0.043–0.086, the relative streamwise pitch (
) was in the range of 0.869–2.173 and the relative spanwise pitch (
) was in the range of 0.434–1.08; the flow Reynolds number (
) was selected in the range of 4000–18,000.
The methodology followed in the present investigation is displayed in
Figure 2. The process of data collection was followed by thermos-hydraulic performance analysis, correlation development, exergy calculation, a sustainability assessment and finally a cot analysis of the thermal system.
3. Data Reduction
The temperature and pressure drop data collected experimentally is presented in the dimensionless form as Nusselt number (
) and friction factor (
), respectively. The computation starts with the calculation of the air mass flow rate (
) from the pressure drop (
) determined by the calibrated orifice plate with area of
. Air is the cooling medium in the rectangular channel, with physical properties (
) calculated on the basis of the temperature of the air at the experimental location. The
is determined using the following formula.
where
and
is the coefficient of discharge.
The flow Reynolds number (
) through the setup of the hydraulic diameter (
) is found by:
The pressure drop (
) generated across the duct determines the friction factor (
) by applying the Darcy equation as
The heat transfer coefficient (
) of the heated plate is used to determine the Nusselt number (
)
Using the uncertainty method proposed by Kline [
23], the highest uncertainties in terms of the mass flow rate: ±2%, Reynolds number: ±2.5%, Nusselt number: ±4.9% and friction factor: ±4.7% was determined.
6. Results and Discussion
Figure 5 shows the variations in the
with
for different ranges of
with a constant
of 1.739 and
of 1.08. The
increases with an increase in
and its values are highest for a
of 0.086. As the
increases, the air flow through the staggered jet increases, which affects more of the heated surface by jet impingement, and this increased flow of air leads to heat transfer augmentation. The intensity with which the air jet impinges on the heated surface increases until a ratio of 0.086, beyond which the turbulent intensity of jet impinging on the surface decreases due to the higher diameter of the jet and thus the lower intensity of turbulence caused by a higher mass air flow rate, with a low striking velocity causing a downfall in the heat transfer rate. The
data for a range of
in terms of roughened ducts at different relative streamwise pitches (
) for a relative spanwise pitch (
) of 1.08 and a relative jet diameter to hydraulic diameter ratio (
) of 0.086 is shown in
Figure 6. The introduction of a protrusion on the absorber plate allows for an increment in heat transfer. As the air mass flow and the number of air impinging jets increases on the heated absorber plate, the heat transfer increases by raising the value of parameter
. The highest
was obtained at an
of 1.739, and beyond this the
decreased because the number of rows in a streamwise direction decreased, which leads to lower reattachment points.
The air jet impinging on the spherical protrusion surface is seen in
Figure 7. The air jet striking the spherical protrusion spreads over the protruded surface as well as the flat heated surface and disrupts the larger laminar sub-layer, thus eliminating larger hot zone areas, which leads to augmented heat transfer.
Figure 8 reveals the effect of
on
at various
values for an
of 1.739 and a
of 0.086. It is seen that the
continuously rises with an escalation in the
up to a value of 1.08, with the
subsequently decreasing. The probable reason for this enhancement is that as the spanwise spacing increases, it provides an area for the jet stream to spread-out on and merge with neighboring jet streams. This amalgamation of the jet streams enhances the turbulence and hence the heat transfer is improved. Whereas an increase in
beyond 1.08 results in larger area between air jet streams, the intensity with which they interact with the neighboring jet stream is lower. Thus, it can be concluded that the
plays an important role in enhancing the
by impinging the spanwise air to eliminate larger hot zones.
Figure 9 shows the effect of
on
at various
values for fixed
value of 1.739 and a
value of 1.08. The results show that the
for a
of 0.043 is in the maximum range of
values investigated. As the size jet holes increase from a
ratio of 0.043 to one of 0.130, more air should be pass through it, which results in a lower
. The variation in
with
for various
values is shown in
Figure 10 for a
of 1.08 and a
of 0.086. It is seen that the value of
goes down with an increasing
. The protrusion-facing surface experiences higher values of
than that of the smooth duct. Since the pressure drop increases because of the presence of protrusions, a secondary flow is generated on either side of the protrusion ducts. Any increase in the
value beyond 1.739 reduces the value of
, as number of jets striking the heated plate reduces, which leads to a lower intensity in terms of the turbulence. The plots show that the maximum
for the duct occur in the case of an
of 1.739. Strength and intensity of secondary flow deteriorated in case of an
of 0.869, as compared to an
of 1.739 and 2.173; hence, the
is lower in such cases. The variation in
with an
of a roughened duct at different
values, while keeping the other parameters values constant, is displayed in
Figure 11, with corresponding values for the smooth duct. The
value decreases with an increase in
. As the
value changes from 0.434 to 1.08, the number of jets impinging on the heated plate reduces, which leads to smaller area for air transfer, which in turn lowers the intensity of turbulence and causes a higher
.
The thermo-hydraulic performance parameter (
) is a standard parameter that determines the thermal and friction performance of heat
a transferring surface, as described by Han et al. [
26] and Webb and Eckert [
27]. The thermo-hydraulic performance parameter (
) is given as:
The
determined for the
,
and
parameters are presented in
Figure 12 for a range of
values. It can be seen from
Figure 12a that the highest value in terms of
is found at a
0.086, and the physics behind this is that at
0.043, the diameter that provides a higher pressure drop penalty is small, and at
0.13, a larger diameter with a lower turbulent intensity produces a small
. Similarly, the effect of
and
on
is presented in
Figure 12b,c. The maximum
is achieved at parameters values for
and
of 1.739 and 0.869, respectively. The best
considering the entire set of geometrical parameters was found to be 3.18 at an
of 15,000, and the staggered pattern also plays a major role in enhancing the heat transfer.
The use of impinging jets on the heated surface leads to a higher heat transfer, and augmentation is seen to be further enhanced by the use of roughness on the heated surface. A comparative analysis was carried out on the basis of the thermohydraulic performance of the previously investigated heated surfaces with air jet impinging on them. The comparative graph plotted in
Figure 13 represents the present investigation results and the results of various investigations.
7. Development of Correlations
The extensive experimentation yielded data in dimensionless terms for the Nusselt number (), friction factor () and thermohydraulic performance parameter () for the selected operating and geometric parameters. The data acquired for the and was observed to be a strong function of the geometric and operating parameters .
The data obtained was used to develop the correlation of
and
by regression analysis to bring about an efficient way of calculating the
and
data for the selected range of parameters by researchers in future. The regression analysis of the data resulted in the correlation for
being:
Correspondingly, the correlation for
developed is written as
Figure 14a–c displays the deviation in the predicted and experimental values for
and
, and the deviation was found to be ±11%, ±12%, and ±9%, respectively.
8. Exergetic Performance
The energetic performance of a staggered impinging jet on an
rectangular duct was evaluated in terms of the temperature of the absorber plate, which is a major parameter for measuring performance. It is evident that when the
increases, there is a higher degree of turbulence intensity, and this leads to the breakdown of the boundary layer formed in the vicinity of the heated surface. This disruption in the boundary layer increases the heat transfer rate from the heated surface to the air, and thus the absorber plate temperature (
) decreases. The lower plate temperature reveals that the convective heat transfer coefficient is higher in the case of air jet impinging on a plate, with method resulting in a higher heat transfer coefficient compared to heat transfer without jet impingement on a smooth plate. Values of 0.869, 1.739 and 1.5 for the investigated geometric parameters
,
and
, respectively, yields a minimum plate temperature in terms of the
duct. The extraction of thermal energy from the heated surface is accompanied by number of exergy losses in the form of heat transfer loss (
), optical exergy loss (
), irradiation loss ((
), friction loss (
) and fluid heat transfer loss (
). The selected geometrical parameters values that deliver the highest
results in the exergy losses encountered in
Figure 15. An analysis of the plot reveals that the impact of
on the optical exergy loss (
) is not prominent. The optical exergy loss (
) is the incident radiation reflected back into the atmosphere from the glass cover. Likewise, it is seen that the other exergy losses follow an asymptotic curve, and these exergy losses does not have a dependency on
. A computation of the total exergy inlet to the system (
), total exergy losses (
) and the exergy efficiency (
) was carried out and the plot is shown in
Figure 16. The plot reveals a maximum exergetic efficiency of 3.87 at a set of geometrical parameters values for
,
and
and an
of 4000.
9. Sustainability Index
The sustainability index (
), waste energy ratio (
) and improvement potential (
) are the suggestive parameters for evaluating exergy losses and corresponding irreversibility in thermal systems [
28,
29,
30]. These three variables are based on an exergy study that ensures long-term viability. The
and
increase with an increment in exergy losses, in which case, correspondingly, the
decreases. The higher the exergy losses, the higher the improvement potential and waste heat that needs to be recovered. These three components are mathematically expressed as:
Improvement potential (W)
The results for the
were calculated for all the studied parameters and are represented in
Figure 17,
Figure 18 and
Figure 19. The
results imply the life time of the system and
describes the exergy loss per unit exergy inlet. The
variation as a function of
for different values in terms of relative streamwise pitch
are shown in
Figure 17a, while the relative spanwise pitch (
) and relative jet diameter to hydraulic diameter ratio (
) are kept constant. The
values vary between 1.0171–1.0347. A higher
value signifies a higher exergy efficiency and, hence, the system is sustainable on these parameters. An
, its optimum value, results in the highest
. Beyond this
value, the
decreases because of a higher pressure drop and lower levels of heat transfer. Similar trends were observed when the (
) and (
) varied, as can be seen from
Figure 17b,c, respectively. The range in terms of
was 1.0219–1.0349 when the (
) varied from 0.434–1.08. Analogous to previous results, the
variation was 1.0218–1.0349 for a variation in (
) from 0.043 to 0.13.
The waste energy ratio (
) represents the exergy loss per unit of exergy available at an inlet. The exergy loss decreases significantly with an increase in the Reynolds number. The decrement in exergy loss signifies that the energy available at the inlet is used for the heating process. As the usage of exergy increases, the
values decrease, as can be seen in the
Figure 18. The
varies between 0.9665–0.9832, 0.962–0.9832 and 0.97–0.9832 according to varying
, (
) and (
) values, respectively. At higher values in terms of
, a higher heat transfer is received and, hence, the lower the
.
Figure 18a, represents the variation in
as a function of
by varying the
from 0.869 to 2.17. The lowest
value was obtained at an
, and the maximum at an
= 0.869. These trends signify that, among the studied parameters, the exergy losses reach a maximum at lower
value and a minimum at higher values until the
reaches 1.739. Beyond an
= 1.739 and for all
values, the
was higher. The apparent reason for this is that at an
= 2.17, the heat transfer decreases and, hence, the
increases. Similar trends in terms of
variations can be seen in
Figure 18b, where
and (
) are kept constant and the
results as function of
by varying the (
) are shown. The results testify that the optimum
value is 0.869, where the
is at a minimum. The results in terms of
by varying the (
) are shown in
Figure 18c, which shows that a
= 0.086 results in a minimum
.
In the preceding section, the findings in terms of
and
variations with geometric and operational parameters were reviewed. The preceding explanation suggests that the system consumption can be sustained even if exergy losses occur during system operation within the stated range. The third parameter employed in this study was improvement potential (
), which is used to demonstrate the maximum improvement in the exergy efficiency of the process. The results on
as a function of
by varying
, (
) and (
) values are appended in
Figure 19a–c. The results signify that the maximum improvement of the system can be achieved by considering an
= 1.739 and corresponding parameter values of
= 0.869 and
= 0.086. The apparent reason for this is that at these values the exergy losses are at a minimum and the exergy efficiency is at its maximum, which ultimately results in a maximal improvement in system processes. The
range lies between 2.69–10.84 W.
Figure 19a shows the variation in terms of
at different
values. The maximum improvement in system performance achieved was 10.84 W at an
1.739, while the minimum vape in terms of IP was 5.52 W at an
= 0.869. The range in terms of
lies between 6.56–10.84 W for a variation in
from 0.434 to 1.08 while
and
are kept constant, as shown in
Figure 19b. The IP results when varying the
and
are shown in
Figure 19c, which shows the range in terms of
as 6.65–10.84 W.