<<

Contemporary Sciences, Vol. 11, 2018, no. 74, 3647 - 3654 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.88396

Development of a CFD Model to Predict the

Thermal Performance of Solar Collectors

Jorge Duarte1, Guillermo E. Valencia2 and Luis G. Obregón3

1 Efficient Management Research Group – Kaí, Universidad Atlántico, Carrera 30 Número 8 - 49, Puerto Colombia – Colombia

2 Mechanical Eng., Efficient Energy Management Research Group – Kaí, Universidad del Atlántico, Carrera 30 Número 8 – 49 Puerto Colombia – Colombia

3 Sustainable Chemical and Biochemical Processes Research Group Universidad del Atlántico, Carrera 30 Número 8 – 49 Puerto Colombia – Colombia

Copyright © 2018 Jorge Duarte, Guillermo E. Valencia and Luis G. Obregón. This article is distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium provided the original work is properly cited.

Abstract In this work, the thermal performance of unglazed transpired solar collector (used in preheating air), is studied in CFD (Computational Dynamics), to analyze the influence of flow parameters on their thermal performance. To do so, variations were made in flow parameters (such as ), to study the response of the system performance by using a CFD model of the fluid inside it. The results obtained by the were validated with the theoretical model, and a mesh independence study was carried out, to optimize the simulation process and computational resources; these results allowed to conclude that the developed CFD model can predict accurately the behavior of the flow inside the collector, and that the increment on the flow velocity generates a reduction in the heat loss due to recirculation of the inlet flow, which in turn improves the efficiency of the collector, making these suitable for air heating processes.

Keywords: Solar collectors, CFD, Performance, Flow parameters, Simulation

1. Introduction

In the last 50 years, the development of energy exploitation schemes has become a 3648 Jorge Duarte et al.

research topic of great interest, given the need to minimize dependence on fossil fuels. At the same time, these schemes seek to integrate alternative sources of energy within the processes that make up certain activities, whether comfort or industrial ones; therefore, they have become mechanisms to increase the efficiency of said processes, without requiring a visible increase in the amount of fuel to be used, which in turn generates money savings that can be used in other processes or applications [1]. From the energy sources available for their use in the schemes above, the solar energy is the most used, as it requires a relatively low investment in its application. Though solar energy can also be used for electricity generation, its main use is in heating appliances, such as solar thermal collectors [2], [3]. This term commonly refers to solar hot panels but may refer to installations such as solar towers or basic installations such as solar air heaters, typically used in residential and commercial buildings for space heating [4], [5]. Due to this wide range of application, many studies have been developed to analyze their performance parameters, for the development of optimization processes that allow achieving higher efficiencies in these applications [2], [3], [6]–[10]. However, the performance of solar collectors requires the understanding of the fluid flow around the collector (as the typically occurs by ); though there are some simplifications that allow studying some flow properties from an overall view, these assumptions reduce the quality and accuracy of the results, which in turn generates some issues related to performance prediction. To solve flow-related issues, CFD (Computational ) has become the most used option, as a method that involves the numerical solution of the that describe the flow [3], [11], [12], has to calculate fluid properties in less time than experimental measurements, with a comparable precision. In this article, a CFD model for an unglazed solar collector was developed, to study the influence of flow parameters, mainly the speed of the air flow, on the collector thermal performance. The results of the simulation were validated using data from the theoretical model, with the purpose of evaluating the precision of the data obtained by these models, and to make visible the importance of the CFD numerical methods in this kind of analysis.

2. Methodology

2.1. Development of the study

CAD model for an unglazed transpired solar collector, was elaborated by using the SOLIDWORKS® package. Following that, a control volume was defined in the inside of the collector, which is the region where the energy transfer takes place. The CAD model generated, was used in the ANSYS® Fluent add-on to simulate a set of fixed points, which are used to plot collector efficiency against the flow speed; regarding the boundary conditions, an inlet velocity was defined for each studied point, while for the second condition, a outlet of 101.325 kPa was used. To ensure the convergence of the solution, the outlet pressure was monitored, Development of a CFD model to predict the thermal performance 3649

with a convergence criterion of 10-4 was considered for the normalized sum of the residuals of the discretized equations. To validate the results from the simulation, the results obtained were plotted against the results given by the theoretical model, and after that, the simulation was used in a parametric study, changing the values of a perturbance applied on the collector (an external current) to compare the influence of the flow speed on the performance, where a disturbance is added.

2.2. Fundamental Equations

For the study of unglazed solar collectors, some theoretical models have been proposed and validated in the literature, to evaluate and predict their thermal performance. From these models, it can be concluded that the useful heat collected 푞푢 can be calculated by the expression [13]

푞푢 = 퐼푐 ∙ 퐴푐 ∙ 훼푠 − 푈푐 ∙ 퐴푐 ∙ (푇표푢푡 − 푇푎) (1)

Where 퐼푐 is the insolation period, 퐴푐 the collector area, and 훼푠 the absorptivity of the , while the second term denotes the losses due to simultaneous radiation and convection from the surface. Within this term, the overall heat loss coefficient 푈푐 is given as

푈푐 = ℎ푟⁄휀ℎ푥 + ℎ푐 (2)

Being 휀ℎ푥 the absorber heat-exchanger effectiveness, while ℎ푟 and ℎ푐 denote a linearized radiative heat-transfer coefficient and a convective heat-loss coefficient, respectively. Then, the heat-exchanger effectiveness for air flowing through the absorber plate is defined using the of the air that exits the collector (푇표푢푡), the collector and ambient temperatures 푇푐, 푇푎, by the [1]

푇표푢푡−푇푎 휀ℎ푥 = (3) 푇푐−푇푎

Taking into account the perturbance given by an external wind current with speed 푈∞, the convective heat transfer coefficient is calculated using the mathematical formula

푈∞∙휈∙푐푝 ℎ푐 = 0.82 ∙ (4) 퐿∙휐0

2 Where 휈 is the kinematic of air in m /s, 푐푝 is the specific heat in J/kg·K, 휐0 is the flow velocity, and 퐿 is the height of the collector. In the other hand, the heat transfer coefficient used to quantify radiative losses to the sky and ground, can be calculated through the equation

4 4 4 (푇푐 −퐹푐푠∙푇푠푘푦−퐹푐푔∙푇푔푛푑) ℎ푟 = 휀푐 ∙ 휎 ∙ (5) 푇푐−푇푎 3650 Jorge Duarte et al.

In this equation, 퐹푐푠 and 퐹푐𝑔 are the view factors between the collector and the sky and between the collector and the ground, respectively, and for a vertical wall with infinite growth und in front of it, both factors take the value of 0.5 [1].

3. Results and discussion

Below are the results of the main points of the present study:

3.1. Mesh independence study

To optimize the time and computational resources, mesh independence was done before the analysis of the simulation results. By monitoring the outlet pressure, it was verified that the convergence criterion was achieved for each point of the simulation, the convergence criterion of 10-4 was achieved. Five different meshes were used to satisfy the grid requirements for the solution approaches; as result of this study, a mesh with a total of 305462 elements or cells was utilized, as it reflects the best compromise between solution-accuracy requirements and computational resources as shown in Figure 2.

Figure 1. Mesh independence test.

3.2. Validation of the simulation

For validation purposes, Figure 2 shows the thermal efficiency calculated from the CFD simulation, as a function of the air flow velocity, along with the predictions from the theoretical equations, for a square collector of 3 m with the parameters 푇푎 2 = 10°C,푇푠푘푦=푇푎−15°C, 푇𝑔푛푑=푇푎, and 퐼푐=700 W/m . The results obtained by the simulation show a close resemblance to the theoretical model, with a maximum error of around 1-2%, due mainly to the simplifications done by the theoretical model, which generates a slight offset from the simulation results. However, even with these differences, it was verified than the simulation model is able to predict Development of a CFD model to predict the thermal performance 3651

the thermal performance of the solar collector, as function of the boundary conditions defined in the CFD package.

Figure 2. Validation of the model.

3.3. Case study: Variation of external wind velocity To study the influence of the perturbation (the external wind current), the system response plot was obtained at two values of velocity: 0 m/s (no win) and 5 m/s, maintaining constant the other parameters. As an intermediate result, the difference between the ambient air and the collector output (which have a high relevance on convective losses) was plotted as a function of velocity, as shown in Figure 3. As it can be seen in the figure, for both cases this difference is high at low velocities, but it diminishes as the flow of air is increased, reaching values close to zero at higher velocities; however, the difference is even higher, in a case where there is no perturbation, reaching values up to 90°C at low velocities, as the dominant process in this stage is the natural convection, which is slower than , increasing the convective losses of the system.

Figure 3. Temperature difference plot for the solar collector. 3652 Jorge Duarte et al.

Regarding the collector efficiency, Figure 4 made evident that the increment of the wind velocity generates a reduction of the efficiency at low velocities, as it slows the movement of the ambient air moving to the collector, which in turn increases the losses and changes the dominant mechanism to the slower natural convection, which negatively affects the heat transfer coefficient and, by extension, the efficiency. However, at higher velocities, the efficiency keeps at a constant value, as the influence of wind velocity becomes smaller when compared to the air flow through the collector.

Figure 4. Efficiency plot for the solar collector.

4. Conclusions

The study developed with the CFD package ANSYS® Fluent, allowed to evaluate the thermal performance of solar collectors and predict its behavior, without the need to do experimental tests, which involves a considerable reduction in resources, keeping the same accuracy level with a lower investment; as an example, the developed model allowed to predict the theoretical results with a offset of 1-2%, with a lower time consumption when compared to the analytical solution. Regarding the thermal performance of the solar collector, it was found that the velocity of air flow that enters the control volume, has a high influence on performance parameters, as it is closely related to the losses due to convective heat transfer. Regarding this fact, it was found that a low velocity, the heat transfer mechanism is severely affected, allowing the heat to escape from the volume control, before the air exits from the solar collector; however, it was also found that the presence of an external wind current reduces this effect, though at higher velocities its effect is negligible, and the heat losses are kept to a steady value near to zero, which also keeps the efficiency at a constant value. The results given by the present work does not take into account the losses due to radiation processes, and it is uncertain its effects on the thermal performance, which can be studied in future works. Development of a CFD model to predict the thermal performance 3653

References

[1] D. Y. Goswami, Principles of Solar Engineering, 3rd ed., Taylor & Francis, CRC Press, 2015. https://doi.org/10.1201/b18119

[2] D. Dović and M. Andrassy, Numerically assisted analysis of flat and corrugated plate solar collectors thermal performances, Sol. Energy, 86 (2012), no. 9, 2416–2431. https://doi.org/10.1016/j.solener.2012.05.016

[3] M. Badache, D. R. Rousse, S. Hallé and G. Quesada, Experimental and numerical simulation of a two-dimensional unglazed transpired solar air collector, Sol. Energy, 93 (2013), 209–219. https://doi.org/10.1016/j.solener.2013.02.036

[4] B. Stojanović, D. Hallberg and J. Akander, A steady state thermal model derived by -theory approach and applied on an unglazed solar collector, Sol. Energy, 84 (2010), no. 10, 1838–1851. https://doi.org/10.1016/j.solener.2010.06.016

[5] S. Kumar and S. C. Mullick, Wind heat transfer coefficient in solar collectors in outdoor conditions, Sol. Energy, 84 (2010), no. 6, 956–963. https://doi.org/10.1016/j.solener.2010.03.003

[6] S. Li and P. Karava, Energy modeling of photovoltaic thermal systems with corrugated unglazed transpired solar collectors - Part 2: Performance analysis, Sol. Energy, 102 (2014), 297–307. https://doi.org/10.1016/j.solener.2013.12.041

[7] S. Li, P. Karava, S. Currie, W. E. Lin and E. Savory, Energy modeling of photovoltaic thermal systems with corrugated unglazed transpired solar collectors - Part 1: Model development and validation, Sol. Energy, 102 (2014), 282–296. https://doi.org/10.1016/j.solener.2013.12.040

[8] R. Kumar and M. A. Rosen, A critical review of photovoltaic-thermal solar collectors for air heating, Appl. Energy, 88 (2011), no. 11, 3603–3614. https://doi.org/10.1016/j.apenergy.2011.04.044

[9] K. E. Amori and H. M. Taqi Al-Najjar, Analysis of thermal and electrical performance of a hybrid (PV/T) air based solar collector for Iraq, Appl. Energy, 98 (2012), 384–395. https://doi.org/10.1016/j.apenergy.2012.03.061

[10] M. R. Assari, H. Basirat Tabrizi and I. Jafari, Experimental and theoretical investigation of dual purpose solar collector, Sol. Energy, 85 (2011), no. 3, 601–608. https://doi.org/10.1016/J.SOLENER.2011.01.006

3654 Jorge Duarte et al.

[11] S. C. Chapra and R. P. Canale, MétodosNuméricos Para Ingenieros, 6th ed., México: McGraw Hill, 2011.

[12] L. A. Tagliafico, F. Scarpa and M. De Rosa, Dynamic thermal models and CFD analysis for flat-plate thermal solar collectors - A review, Renewable and Sustainable Energy Reviews, 30 (2014), 526–537. https://doi.org/10.1016/j.rser.2013.10.023

[13] M.D. Herrera, E. Castro, J. Duarte, G.E. Valencia and L.G. Obregón, Second Law Analysis of Supercritical Partial Cooling Brayton Cycle with Recompression, Contemporary Engineering Sciences, 11 (2018), no. 44, 2149 – 2156. https://doi.org/10.12988/ces.2018.85220

Received: August 20, 2018; Published: September 4, 2018