CFD Numerical Investigation of a New Solar Flat Air-Collector Having Different Obstacles with Various Configurations and Arrangements

,is work deals with a numerical parametric optimization study of a new Solar Flat Air Collector (SFAC) configuration.,e CFD numerical parametric study investigates various SFAC structures inside the air cavity without obstacles and with spherical, cubic, cylindrical, and pyramidal obstacles. ,e study optimizes the most convenient configuration and arrangement that allow for the increase of the heat-transfer surface and to make the flow homogeneous in order to generate turbulence zones inside the SFAC air cavity.,e result shows that the thermal performances of the cubic form are close to those of the spherical obstacles. Another set of simulations was performed to evaluate the performances of the cubic shape baffles for three orientation angles equal to 0°, 22.5°, and 45°, respectively. Each configuration has three forms of arrangement with a relative roughness pitch (b/a) varying between 2, 4, and 6.,e results of the simulation study showed that the relative roughness pitch, the Reynolds number, as well as the angle of orientation influence the performance and the operation of the SFAC.,e results of the simulations showed that the combination of an orientation of 45° with a roughness pitch of b/a� 2 increases the SFAC thermal performances, which can reach 85%.


Introduction
e exploitation of renewable energies such as solar energy, which is a clean and inexhaustible source, presents a promising solution against environmental problems and the future exhaustion of fossil fuels. is solar energy can be converted through using solar collectors for the production of thermal energy. e range of application of this thermal energy produced is broad such as drying, heating, or desalination. Given the multiple difficulties that characterize these sensors, such as the presence of dead zones in certain places, thermal losses, and low inertia, various researchers have worked on improving thermal performance on solar collector's air plane. Among the improvements proposed by these researchers are the judicious choice of forms of design parameters (geometric, thermo-physical, and optical), such as the length, the width, and the height of the channel; the thickness of the air space; the glazing; the insulation; and the choice of the absorbing material or the intention to bring back modifications to the absorbing plate. In this context, and according to numerous research studies, the integration of obstacles in the air passage channel presents an effective method for improving the convective heat transfer process by creating turbulence, thereby minimizing dead zones. Indeed, the use of fin types, OT, OIF1, WT, WDL1, and WOL1 [1], and the integrating plus shaped baffles [2] or the planting of triangular, circular, and hexagonal obstructions [3], prove the creation of a turbulent flow, a good distribution of the fluid resulting in the reduction of dead zones.
is improves the outlet temperature as well as the efficiency of the collector. For the same purpose, an experimental study carried out by Sompol et al. [4] reveals that the combination of engraved grooves and vortex generators which correspond to fins in the delta form with holes in the middle placed on the plate of the absorber with a very precise ratio of fin area and hole area with an angle of attack equal to 45°, provides the greatest number of Nusselt, and gives a friction factor around 6 and 30 times greater in the case of smooth channel. e technique of using a multi-pass solar collector gives rather important results for the improvement of yield. Using numerical and experimental analysis, Doguş Tuncer et al. [5] showed that heat transfer can be improved by dividing the interior of the solar collector channel into four parts. ey reported that this quadruple pass sensor (QPSAC) used for drying can achieve an average thermal efficiency equal to 80.66% considered efficient compared to ordinary collector and can reduce the drying time. In this context, the form and the geometry of the solar collector can enhance its thermal performance as proposed by Singh [6].
is method consists of making a double pass model containing three in-lines, offset, and hybrid fin shapes. ese models are digitally studied (CFD) and prove that these types of fins offer better thermo-hydraulic performance which can reach 79%, 13% more than a single air passage design. e change of the internal structure by increasing the residence time of the air, considerably improves the thermal efficiency. Hu et al. [7] conducted experimental and numerical studies to highlight the influence of the installation of internal deflectors on the prolongation of the residence time of the air.
e results indicate that this technique enhances the internal air disturbance that affects convective heat transfer. e authors also assume that radiant heat losses can be reduced by using three upper panes which improve performance. Numerical and experimental studies of the tube-type solar air heater are performed by Afshari et al. [8] e effect of adding various turbulators was modeled using CFD approach, and the obtained numerical findings were used in the experimental device. e maximum instantaneous heater efficiency was ascertained to be 72.41%. In the same context, Sözen et al. [9] studied the thermal performance improvement of tube-type indirect solar dryer with iron mesh modification. e results show that the use of mesh modification improved the thermal efficiency of the tube-type heater by about 11%. e same idea is the subject of considerable work done by Ebru and Fatih [10] who presented a study dealing with the analysis of the thermal performances of four types of a solar collector. e first is a simple solar collector and the three other solar collectors having different obstacle shapes. To reduce dead zones and to create turbulence, they also added small baffles of triangular, rectangular, and sheet types. It was found that the solar collector equipped with small obstacles reached the most important efficiency. It is also found that this configuration allows the reach of the highest output temperature.
Considerable work has been done by Zheng et al. [11] who has studied thermal performance as well as developed mathematical models of a new metal-corrugated air solar collector in the heating field in cold-climate buildings. is study has focused on the effects of the key structural and operating parameters and the analysis of the results obtained prove that they have a great influence on the thermal performance of the collector. In fact, the thermal efficiency of the collector increases with the increase of the settlement height, the specific surface area, the intensity of the radiation, the air velocity, and the ambient temperature, whereas it decreases with the increase of the width of the corrugated pad and the inlet air temperature. e thermal efficiency and the net exegetical efficiency of the sensor can reach 73%. e improvement of the thermal performance of the solar collector with the increase of convective heat transfer can be effected by the deformation of the absorbent plate which has an artificial roughness that favors the increase of the exchange surface and the creation of turbulence. Among the rough surfaces that are treated in this perspective are Vshaped rib roughness [12], Discrete W-shaped rib roughness [13], Multi V-shaped rib roughness [14], Arc shaped rib roughness [15], W-shaped rib roughness [16], and Transverse wire rib roughness [17]. ese studies prove that the variation of relative roughness step affects the operating performance of the sensors, hence the need to identify the optimal values. With the development of tools and computer software, three-dimensional simulation through CFD software provides a solution to see some physical phenomena or calculate some parameters that are difficult to measure experimentally. By the using od a CFD program, Manjunath et al. [18] proposed a numerical simulation analysis of an air solar collector having a sinsoidal form absorber plate. is study was done in a range of the Reynolds numbers between 4000 and 24000. e authors established that the presence of sinusoidal ripples increases the flow disturbances and significantly increases the heat transfer and as a result the average thermal efficiency is increased by about 12.5% compared to the ordinary air solar collector where the plate is smooth. Manjunath et al. [19] studied the influence of spherical obstacles on thermal efficiency and thermohydraulic performance on the solar collector. e analysis is carried out with a three-dimensional simulation with the CFD software where the turbulent model chosen is K-epsilon. e simulation shows that the thermal efficiency increases with the increase of the diameter of the sphere and the reduction of the relative roughness step. e thermal efficiency can reach up to 23.4% in percentage and the highest number of Nusselt is 2.5 times higher compared to the ordinary collector. e two numerical simulations done by the CFD software are in good agreement with the experimental results. is proves that this tool is potentially profitable for parametric studies of an air solar heater with a saving of time and without having costs like the experimental works. Various numerical and experimental studies presented in the literature show that the use of voluminous baffles inside the solar collector cavity makes it possible to optimize heat transfer conditions better than simple baffles that do not have large surfaces.
In this sense, we present a numerical study using the CFD software comparing the different volume shapes of the most famous baffles (spheres, cubes, cylinders, and pyramids) and identify the most efficient shapes. Our study aims to perform a special air solar collector configuration and design by developing a parametric optimization and investigation. e 2 Mathematical Problems in Engineering originality of our work is that we present a solar collector configuration with cubic obstacles integrated on the plate absorber of the solar collector. In this study, many angles of orientation are tested in order to enhance the thermal efficiency of the solar collector. e investigation is conducted by CFD numerical parametric study which intended to test the effect of the forms of baffles (Obstacles) on the heat transfer inside the air cavity. e study also deals with the effect of baffles and their orientation on the global thermal efficiency of the SFAC by creating turbulence zones and by minimizing the dead zones in order to avail of the absorbent surface inertia. e study was carried out for different relative roughness pitch b/a equal to 2, 4, and 6 for a turbulent regime in the range 2500 ≤ Re ≤ 20000.

Dimensional Numbers.
e dimensionless numbers that characterize a flow and are required for our study are mainly: where D h presents the hydraulic diameter given by

ermal Efficiency.
e thermal efficiency of SFAC is defined as the ratio of the actual thermal power exchanged and the theoretical maximum power: where Q u � m . C p ΔT and Q max � I.S

Mechanical Pumping
Power. e mechanical pumping power is defined as the force required for forcing the air to circulate inside the solar collector. is power noted Pm is given in the form: 2.4. Effective Efficiency. By using two expressions of thermal efficiency and mechanical power, we can write the mathematical expression of effective efficiency as follows as defined by Cortes and Piacentini [20] and Mittal and Varshney [21]: where the value of the constant C recommended by [20,21] is 0.18. is effective efficiency is used to evaluate the thermo-hydraulic performance of the solar collector.

CFD Simulation
e simulation was performed by using the CFD Fluent software. In this numerical program, various steps precede the launch of the simulation: (i) e first step is to create the geometry to study. It is followed by the generation of a mesh that can be structured or unstructured in certain areas. (ii) After the integration of the parameters and the input variables of the model, the software can solve the equations necessary for the simulation in every element of the considered mesh until the convergence. (iii) When the convergence is obtained, the numerical model will be exploited for the evaluation of the Solar Flat Air Collector (SFAC) performances.

Description of the Solar Collector
Configurations. e solar collector used for our study is a Solar Flat Air Collector (SFAC) that consists of an inlet section, an outlet section, and a cavity that contains an absorber containing baffles. is configuration is inspired by a publication by Manjunath et al. [19] illustrated in Figure 1, which presents a validation support for our study.
(i) e absorber situated at the top of the SFAC has a dimension of 1200 * 150 mm 2 and a thickness equal to 0.5 mm. It was conceived with aluminum material whose thermal-physical properties are as follows: density 2719 kg/m 3 , specific heat Cp � 871 J/kg K, and thermal conductivity 202.4 W/m.K. (ii) e SFAC is equipped with inlet and outlet air duct sections that are evaluated and selected on the basis of the ASHRAE 93-77 [22] standard for the turbulent regime. (iii) On the SFAC absorber is fixed a series of baffles fabricated with aluminum materials and having a volume of 20 mm 3 and having different shapes (Tables 1 and 2). e different SFAC configurations obtained by these management will be studied and compared on thermal basis in order to find the best configurations characterised by the highest performances. (iv) e arrangements of these baffles inside the cavity and on the absorber plate are carried out in the following arrangement ( Figure 2). e values of a, b, and c are summarized in Table 3 for all proposed arrangements.

Meshing Treatment.
e developed FLUENT code is based on the discretization of the studied configuration by cutting it into structured or unstructured elementary volumes which will be taken as control volume (hexahedrons, tetrahedrons, and prisms) using the finite volumes method. In order to minimize the computation time and to keep the Mathematical Problems in Engineering precision of the meshes, which have an importance in the quality of the results, the following approaches were taken: (i) e use of symmetry along the solar collector for certain symmetrical configurations. (ii) e simplification of the solar collector geometry by eliminating irrelevant structures for the simulation. We mainly consider two sections of air entry and air exit. (iii) e choice of the most suitable mesh which will enhance the quality of the results and reduce the time of computation by conducting mesh independence ( Figure 3). Figure 3 and Table 4 illustrate this test for the solar collector with the simple absorbent plate with different Reynolds numbers (Re 1 � 5000, Re 2 � 7500, Re 3 � 10000, Re 4 � 12500, Re 5 � 15000, and Re 6 � 20000) and each time the ΔT (difference between inlet and outlet temperature) is calculated for each mesh. Table 1: Dimension of the baffels to be integrated inside the air cavity.

Form Dimensions Spheric
Diameter � 20 mm Cubic Length of an edge � 20 mm Cylinderical Diameter of the base � 20 mm/Height � 20 mm Pyramidal Side length of the base � 20 mm/Height � 20 mm

Turbulence Model Choice for Simulation.
For the resolution of the RANS (Navier-Stokes equations obtained by passing to the Average Reynolds), we have in Fluent several classes of turbulence models, which are equations that complete the balance equations and allow the closure of the system [23][24][25]. Considering the geometry and the flow conditions, we have chosen the models k − ε which is judged by several authors who prove that the latter is better adapted to the physical problems of flows and gives better flow conservation [26,27].

e Equations Needed for the Simulation.
e necessary equations for the flow simulation are classical NAVIER-STOKES mass and momentum conservation equations, which are written as follows for an incompressible fluid: (ii) Momentum equation (iii) Energy conservation equation (iv) Transport equations e introduction of transport equations makes it possible to follow the evolution of certain characteristic variables of turbulence. eir expressions given below were obtained after manipulation of averaged and nonaveraged Navier-Stokes equations [23]: One related to turbulent kinetic energy k: And, the other one about the rate of kinetic energy dissipation ε: For the model, we chose k − ε Realizable:

Properties of the Air Inside the Solar Collector
Cavity. e thermo-physical properties of air vary with pressure, temperature, and humidity. For the accuracy of the results obtained from the simulation, we will integrate the following thermo-physical properties in the selected turbulence model that are valid for a the temperature ranges between 280 and 470 K [22,[23][24][25][26][27][28][29][30]:   Mathematical Problems in Engineering

Boundary
Conditions. e boundary conditions chosen for the computer model are as follows: (i) At the inlet of the cavity is maintained an initial air velocity that corresponds to a Reynolds number between 2500 and 20 000 with a temperature equal to 307 K. (ii) At the outlet of the channel, applied pressure is equal to 101325 Pa. (iii) e solar collector is considered as an isolated and an adiabatic cavity. (iv) e solar radiation flux applied on the solar collector surface and on its components is equal to 1084 W/m 3 .

Parameters of the Solver and the Simulation
Convergence.
(i) For our simulation, the solver is based on double precision pressure as well as a SIMPLE algorithm for coupling pressure and velocity. (ii) e choice of a Second-Order discretization scheme leads to a rapid pressure convergence compared to other schemes. (iii) e convergence of the iterative process is determined by the residue concept. Discretized equations are considered to converge when the curves become stable and relatively flat: less than 10 −6 for energy and less than 10 −5 for momentum and continuity equations. e comparison of the variation of the Nusselt number and thermal efficiency versus the Reynolds number between the referent numerical study and our digital study is illustrated in Figures 5 and 6 , respectively. e obtained results and the results given by Manjunath et al. through the two curves followed the same trends and present very similar values. However, a little difference between the two studies is noted. It is found that the thermal efficiency for Reynold's superior to 17103 and Reynold's inferior to 7103 ranges, respectively, between 17 and 14%. It was found that the two CFD models are very close. e mismatch between the two results can be explained by the specific conditions of the CFD program proposed by Manjunath et al., compared to our numerical program. However, we can globally consider that our model with the imposed limiting conditions can reproduce the results proposed by Manjunath et al. Consequently, we can exploit our developed model in the investigation of the different arrangements and configurations in order to obtain the most appropriate and feasible Solar Flat Air Collector (SFAC). adapted configuration that can create more turbulence inside the SFAC air cavity. e results are represented in Figures 7(a)-7(c). It was found that the lowest Nusselt values are attained for the ordinary configuration (without obstacles). Besides, it was found that the Nusselt number increases with the Reynolds number for the different configurations. It was seen that the spherical shape with a relative roughness pitch of b/a � 2 has the highest Nusselt number which can exceed 120. We also noted that the results obtained with the use of the cubic form are close to those found for spherical obstacles. It is important to assert that, for the same configuration, the same relative roughness pitch and for the same Reynolds number, the present study offers a good enhancement of the Nusselt number compared with results published in a previous paper [31].

Effects of Obstacle
To highlight the SFAC with a cubic obstacles configuration, a numerical simulation was proposed. It investigates the effect of the cubic obstacles on the air flow inside the SFAC air cavity. In Figures 8 and 9 are represented the evolution of the air-flow velocity inside the SFAC cavity without ( Figure 8) and with cubic obstacles (Figure 9) (Relative roughness pitch b/a � 2) for a chosen Reynolds number equal to 15000. It was found that for an ordinary configuration (without obstacles), the air flow is almost uniform inside the cavity and the velocity becomes slower towards the ends. e results also showed that inside a cavity containing cubic obstacles, we can observe the creation of turbulence in the vicinity of the obstacles which makes the medium homogeneous and cancels the dead zones inside the air cavity of the solar collector.

Effects of Relative Roughness Pitch (b/a) on the Cubic
Obstacles SFAC Configuration Performances. In Figure 10 is represented the variation of the Nusselt numbers as a function of the Reynolds number without and with cubic baffles for different pitch; b/a � 2; b/a � 4; and b/a � 6. It is found that the Nusselt number increases as a function of the Reynolds number. Indeed, the Nusselt number reaches a value close to 120 for a Reynolds number equal to 23 10 3 and a relative roughness pitch b/a � 2. In Figure 11 is represented the influence of the relative roughness pitch (b/a) of the cubic obstacle on the turbulent intensity distribution and on the air flow inside the SFAC air cavity. It was found that the turbulence varies remarkably depending on the position inside the solar collector cavity. It becomes intense in the vicinity of obstacles. It was also clear that this turbulence increases by reducing the relative roughness pitch. It was found that for a relative roughness pitch b/a � 2, the SFAC reaches the most appropriate turbulence intensity rate. e performance of our solar collector with its four configurations was also evaluated by considering other parameters such as the thermal efficiency, the pressure drop, the mechanical pumping power, and the effective efficiency, which are studied according to the Reynolds number changes.
In Figure 12 is represented the variation of the thermal efficiency of the solar collector as a function of the Reynolds number. Figure 12 shows the variation of the thermal efficiency of solar collector with a smooth absorber and with cubic obstacles having different relative roughness pitches. It was noted that the thermal efficiency which characterizes the performance of our solar collector varies in an increasing way for the four solar collector configurations. It was also seen that the evolution of thermal efficiency depends on the roughness step. e highest values are recorded for the solar collector with a roughness pitch b/a � 2 which can reach a maximum of about 80% for a Reynolds value equal to 23 10 3 .
In Figure 13 is represented the Reynolds impact on the pressure drop for the different SFAC configurations. It was showed that during the air-flow circulation, the initial pressure increases with the Reynolds number increasing. We also found that this pressure drop increases by decreasing the space between the cubic obstacles. e results is explained by the fact the volume of the obstacles integrated inside the air cavity reduces the airflow inside the air cavity. In Figure 14 is represented the variation of the mechanical pumping power of the solar collector for different configurations. It was seen that the mechanical pumping powers of the different configurations according to the Reynolds number have the same tendency of the pressure drop variation. Due to the strong increase in pressure drop between the inlet section and the air outlet section, the model reacts in a way that ensures flow invariability from where constant air circulation is kept in the air cavity of the SFAC. It was also        Mathematical Problems in Engineering seen that the air pumping force increases by the decreasing of the distance between the cubic obstacles. In Figure 15 is illustrated the values of the effective efficiency obtained during the simulation by varying each time the Reynolds number for the four model configurations. It is noted that the effective efficiency which characterizes the performance of the solar collector and its capacity to achieve its ends increases considerably as a function of Reynolds number up to a value of 17500 with the favoring of the configuration with the relative roughness pitch of b/a � 2, which can reach 70%. e second configuration considered effective is that of relative roughness pitch b/a � 4 as long as the third efficient model is reserved for the largest roughness pitch. e least efficient configuration is that of a solar collector with the smooth absorbent plate. Beyond Reynolds value of 17500, the values of the effective efficiency undergo a decrease. Consequently, the most appropriate conditions to have an optimal thermal efficacy (about 70%) is to assure an air pumping with a Re number equal to 17 500 and with a relative roughness pitch b/a � 2.

Area Attack Orientation Effects on Flow and on ermal
Performances of the SFAC. Considering the promising results of the use of the cubic form as baffle obstacles for the performance improvement of the SFAC, another set of simulations were conducted. e numerical investigation was directed towards the study of the effect of the cube angle orientation on the thermal performances of the SFAC and to compare with spherical obstacles in order to show that the change of the orientation improves the cubic configuration. In this context, three orientations are proposed (Table 5): In Figures 16(a)-16(c) is represented the variation of the Nusselt number for the spherical and the cubic configurations as a function of Reynolds for a relative roughness pitches: b/a � 2, 4, and 6, respectively. e simulation shows that with an orientation of 45°we can reach a total attack surface equal to 800 mm 2 since the attack is made on two surfaces of 400 mm 2 each. It is shown that the Nusselt numbers increase with the increase of the Reynolds number and sense its low values for a Reynolds number equal to 2500. By analyzing these curves, we can notice that the Nusselt value for the configuration which contains cubes with an orientation of 45°with a roughness pitch b/a � 2 is greater compared to the other configurations studied, even those that contain spheres. For a Reynolds number value of 23000 and for a configuration which combines the orientation of 45°and the relative roughness pitch of b/a � 2, the     Nusselt number reaches a large value which exceeds 7% of the value of Nusselt for the configuration which contains spherical baffles and 55.46% higher than a configuration which contains a smooth plate.
To highlight the effect of the cubic obstacle angle orientation, we presented in Figure 17 the variation of the Nusselt number for the three angles of orientation 0°, 22.5°, and 45°for a relative roughness pitch equal to b/a � 2 and  configurations. e promising performance of the cubic shape with an orientation of 45°is also visible on the curves below (Figures 18-20 ) which describe the evolution of thermal efficiency as a function of the Reynolds number for the three studied relative roughness pitches. Figures 18-20 show that the thermal efficiency of the sensor is sensitive by the cube orientation angle as well as the roughness pitch. e results of the simulations show that the thermal efficiency increases with the increase of the orientation angle of the cubic baffle and with the decrease of the relative roughness pitch. We can therefore see that the combination of an orientation of 45°with a roughness pitch of b/a � 2 increases the SFAC thermal performances, which can reach 85%.
In Figure 21 is represented the variation of the pressure drop with the Reynolds number for the different orientations with respect to the air flow for a roughness step equal to b/ a � 2. It was seen that the pressure drop increases with Reynolds increases for the different configurations. e low values of the pressure difference are recorded for the configuration where the absorbent plate is smooth as long as the highest values are those of the configuration containing cubic baffles with 45°orientation. ese intense values are explained by the existence of tabulators compared to ordinary solar collectors which break the speed of air flow during flow and by the increase in the attack surface compared to other configurations.
In Figure 22 is represented the evolution of the mechanical pumping power defined as being the power necessary to ensure the uniformity of air flow along the solar collector cavity vs. the Reynolds number variation for the spherical and cubic configurations. e results show that the pumping power for the two configurations increases with Reynolds number increases. It was found the lowest power is obtained for a model without obstacles fixed on the absorber. It was also seen that for the same Reynolds value, the cube oriented 45°requires greater pumping power than the other models. is result is explained by the large surface exposed directly to the air flow compared to the other orientations.
is increase in surface increases the resistance which favors an increase in the pressure drop, hence a high power to maintain the flow uniformity. In Figure 23 is represented the    It was found that the collector with cubic baffles having an orientation angle of 45°is more efficient than the other solar collectors with cubic baffles having the orientation of 0°22.5°, respectively. It was also seen that cubic baffles having an orientation angle of 45°are more effective than the spherical baffles. It can be added that the importance of achieving effective efficiency for a lower Reynolds number is keeping the solar collector inertia for a long time.

Effects of Reynolds Number on Performance.
In order to study the impact of the Reynolds number on performance, we represent in Figure 24 the variation of thermal efficiency as well as (Pm/C * S * I) vs. the Reynolds number variation for a relative roughness pitch b/a � 2. It was found that for high Reynolds numbers, the rate of increase in mechanical pumping power is very significant, while the rate of thermal   efficiency becomes almost constant. On the other hand, for low Reynolds numbers, the rate of increase in mechanical pumping power, necessary to keep the air flow constant, is low, while the rate of thermal efficiency is important. Between these two extremes of the Reynolds number, we can reach a stage where we obtain the optimal value of the energy gain measures of 76% for a Reynolds number equal to 15 000. Figure 24 also shows that for a large Reynolds number value, the net energy gain disappears and decreases towards lower values contrary to the desired values. Consequently, as a summary of this treatment, we can limit the Reynolds value to 15 000 for the correct functioning of the solar collector.

Conclusions
In this paper, we present a numerical study using the CFD software to perform a special air solar collector configuration of the most famous baffles (spheres, cubes, cylinders, and pyramids) to fix on the absorber of the SFAC by developing a parametric optimization and investigation by defining a geometric form of an orientation angle. e study was carried out for different relative roughness pitches b/a equal to 2, 4, and 6 for a turbulent regime 2500 ≤ Re    configurations in order to obtain the most appropriate and feasible Solar Flat Air Collector (SFAC). In this paper, we have also investigated the obstacle geometries and arrangements effects on air flow inside the SFAC air cavity. In this context, a simulation based on the most adapted configuration that can create more turbulence inside the SFAC air cavity was performed. To highlight the SFAC with a cubic obstacles configuration, another numerical simulation was proposed. It studied the effect of the cubic obstacles on the air flow inside the SFAC air cavity. e results of the simulation also showed that inside a cavity containing cubic obstacles, we can observe the creation of turbulence in the vicinity of the obstacles which makes the medium homogeneous and cancels the dead zones inside the air cavity of the solar collector. We have also studied the effect of the relative roughness pitch (b/a) on the Nusselt numbers. It was found that the Nusselt number increases as a function of the Reynolds number. e simulation showed that the turbulence fluctuate remarkably depending on the position inside the solar collector cavity. We have also studied the configurations and the arrangement of the obstacles effects on the thermal efficiency variation of the SFAC. e impact on the pressure drop for the different SFAC configurations was also investigated in this paper. It was showed that during the air-flow circulation, the initial pressure increases with the Reynolds number increasing. A simulation of the Nusselt number variation for the spherical and the cubic obstacles configuration is carried out. e simulation is achieved for a relative roughness pitches equals to b/a � 2, 4, and 6, respectively. e simulation showed that the Nusselt value for the configuration which contains cubes obstacles with 45°orientation and with a roughness pitch b/ a � 2 is more important compared to the other configurations studied, even those which contain spheres. e effect of the cubic obstacle angle orientation on the Nusselt number variation for a relative roughness pitch equal to b/a � 2 and a flow rate that corresponds to a Reynolds number equal to 23000 was also studied. It was showed that with an orientation of 45°and with a roughness pitch of b/a � 2, the thermal efficiency can exceed 85%.
We conclude that SFAC contains cubes obstacles with 45°orientation and with a roughness pitch b/a � 2 confers a promising solution to enhance the solar collector performances. It is obvious that further studies on SFAC applications in Tunisian households are necessary. e results of this investigation will be accompanied with economic study to highlight the SFAC configuration and arrangement.

Data Availability
e data used to support the findings of this article were constituted by authors mentioned in the paper. Also, all the data are included within the article and no external data were used to support this study.

Conflicts of Interest
e authors declare that they have no known conflicts of interest or personal relationships that could have appeared to persuade the work reported in this paper.