| Computers, Materials & Continua DOI:10.32604/cmc.2022.027180 |  |
| Article |
Thermal Loss Analysis of a Flat Plate Solar Collector Using Numerical Simulation
1Institute of Information and Computational Technologies CS MES RK, Almaty, Kazakhstan
2International Information Technology University, Almaty, Kazakhstan
3Al-Farabi Kazakh National University, Almaty, Kazakhstan
4Astana IT University, Nur-Sultan, Kazakhstan
*Corresponding Author: Timur Merembayev. Email: timur.merembayev@gmail.com
Received: 13 January 2022; Accepted: 13 May 2022
Abstract: In this paper, we studied theoretically and numerically heated losses of a flat solar collector to model the solar water heating system for the Kazakhstan climate condition. For different climatic zones with a growing cost for energy or lack of central heating systems, promising is to find ways to improve the energy efficiency of the solar system. The mathematical model (based on ordinary differential equation) simulated the solar system work process under different conditions. To bridge the modeling and real values results, we studied the important physical parameters such as loss coefficient, Nu, Ra, and Pr values. They impacted the efficiency of flat solar collectors and heat losses of the system. The developed mathematical models, the design and composition of the software and hardware complex, and automated control and monitoring systems allow solar hot water heating systems to increase the energy efficiency of life support systems and heat supply of buildings by reducing energy consumption for heat supply. The simulation result showed that during the daytime, the temperature of water in the collector is 70°C; the storage of heated water since heated water is cooled at night. We defined that a work period of the system can be extended with high efficiency (April–October) for Almaty region.
Keywords: Solar heating system; heat loss coefficient; dynamic simulation; flat plate collector
Nomenclature:
| Area (m2); | |
| Temperature of output fluid from a heat exchanger (°C); | |
| Volume of fluid in a collector cycle (m3); | |
| Specific heat capacity in a collector cycle (J/kg K); | |
| Temperature of outlet fluid from collector (°C); | |
| Temperature of inlet fluid into collector (°C); | |
| Volumetric flow rate in collector cycle (m3/s); | |
| Solar irradiance (W/m2); | |
| Temperature of collector ambient (°C); | |
| Heat loss coefficient of collector (W/m2°C); | |
| Thermal conductivity coefficient (W/m2°C); | |
| Radiation coefficient from the glass coating to the absorber (W/m2°C); | |
| Wind heat transfer coefficient (10 W/m2°C); | |
| Radiation coefficient between the glass surface and the air (W/m2°C); | |
| Heat transfer coefficient (W/m2 K); | |
| Distance between the absorbent and the glass; | |
| Heat conductivity of the insulating material (W/mK); | |
| Gravitational constant ( | |
| Temperature difference between the boards; | |
| Kinematic viscosity (m2/s); | |
| Absorbent temperature; | |
| Glass surface temperature; | |
| Nusselt number | |
| Prandtl number | |
| Greek symbols: | |
| Efficiency (precentage) | |
| Density of fluid in colletor (kg/m3) | |
| Volume expansion coefficient; | |
| Thermal diffusivity (m2/s); | |
| Absorbent emission (0.88); | |
| Glass emission (0.95); | |
| Stefan-Boltzmann constant (5.670367 × 10−5 W m2 K−4) | |
Each year renewable energy becomes an important role in our daily lives. Using renewable energy can reduce energy consumption and decrease our dependence on fossil energy. The most straightforward way to use renewable energy is solar energy which allows getting clean energy without pollution. Researchers study mathematical models and designs of solar heating systems for domestic buildings.
One of the researches is devoted to estimating solar hot water heating and analyzing data in South Korea for three years [1]. A solar water heating system was installed for multi-family housing with 1179 families in 14 units. Results of the research showed that comparison of conventional boiler for heating water (efficiency = 85%) with solar hot water heating system showed a positive environmental effect and was estimated to reduce 71.9 L/year of oil and reduction of 186.3 tons CO2. The authors of the paper [2] provided the numerical experiment of a centralized solar hot water heating system developed for a high-rise residential building in Hong Kong.
The author of the research [3] considered a solar hot water heating system and simulated the solar system with the flat collector, storage tank, and circulation. It analyzed two designs of hot circulation: water–water and water–air. The result of the research showed that the design of water–air does not become effective. In the paper [4], there was a review of research on solar energy absorption by a solar system with different fluid types. The main focus of the research was a solar-based absorption cooling systems, diffusion absorption systems, ejector–based absorption systems, compression absorption systems and cogeneration/trigeneration absorption systems. The thermodynamic properties of the most common working fluids were reviewed and analyzed for a ternary mixture in solar absorption systems. The author of the paper [5] calculated optical, thermal, and thermodynamic analysis of solar collectors and described methods of an estimatied efficiency.
In the research [6], the author investigated the construction and efficiency of solar hot water heating systems. The geometry and dimension of a solar system were defined based on the material’s thermal properties. In research [7], the authors proposed to use the algorithms for calibration of small solar hot water heating systems and the application in Matlab.
Thermal performance is calculated based on the first law of thermodynamics (energy), but that does not allow estimating a flat solar collector [8]. The second law of thermodynamics (exergy) estimates the various losses in the flat solar collector and allows for estimating a solar hot water heating system [9–11]. The authors [12] used the numerical inverse Laplace homotopy technique for solving some interesting 1-D time fractional heat equations. The authors of [13] have provided the numerical arrangement of ordinary differential equations (ODE) that appeared from different problems; in the simulation of a solar hot water heating system, it can be considered to decrease the simulation time. In [14], the author provided analytical and numerical techniques for physical fields in the physical time-domain; the problems that appeared in a solar hot water heating system due to the mechanical deformation during the thermal gradient temperature are applied on the outer surface. In the paper [15], the authors presented numerically and graphically in isotherms, streamlines, local Nusselt number, global Nusselt number, and global fluid temperature.
In Kazakhstan, it is extensive attention to developing a solar hot water heating system. The system’s price is high; therefore, it is unavailable for all customers. The new hybrid energy system with mathematical methods and computer simulation, software, and hardware help optimize the price of flat solar water heating systems. One way to efficiently use energy is to use a new source of energy (renewable and environmentally friendly) in Kazakhstan’s fuel and energy system. So the development of an energy system based on the double-circuit solar system with a heat pump is an actual and urgent problem for autonomous power supply.
Kazakhstan has a high potential for the generation of solar energy potential. In the paper [16], the authors used the actual observations and theoretical calculations to generalize the most favorable period and localization of to use the solar energy in the Almaty region. The authors considered the solar collector with a storage tank, heat exchanger, and pump stations for hot water. The simulation was performed for different temperature conditions and solar irradiance for the Almaty region. The simulation process is done in Matlab and Simulink, and the tool is often used for other simulation processes and, in particular, for the simulation of solar hot water heating systems [17,18]. Using simulation, we evaluate the technical possibility of a solar water heating system and indicate system nodes to be modified and improved for Kazakhstan weather conditions.
In this paper, we investigate a solar hot water heating system in different Kazakhstan weather conditions. Previously, no studies have been devoted to simulating other Kazakhstan weather conditions. We focus on providing calculation parameters of a solar hot water heating system and the result of a simulation based on collected weather conditions. We consider a result of simulation to get knowledge about the effective heating system with parameters that will be effective in different regions.
The rest of the paper is organized as follows. In the next section, we describe the wavelet transformation, data analysis, and machine learning algorithms. The theoretical model is presented in Section 2. Energy analysis of a solar hot water heating system is presented in Section 3. The result of the simulation is presented in Section 4. Section 5 concludes the paper.
The solar system consists of several parts. To simplify the process of description and simulation, we describe each part of the system separately. Fig. 1 shows a solar hot water heating system with a heat exchanger coil in a storage tank below. The heat exchanger coil is installed inside the storage tank, and it allows heat water for consumption: heating and domestic hot water.
Figure 1: Schema of solar hot water heating system. Constructed solar hot water heating system. The solar hot water heating system. 1–flat solar collector; 2–storage tank with a heat exchanger; 3–pump
The central object of research is a flat solar collector, and we described a mathematical model of a flat solar collector separately Fig. 2 and energy analysis. The mathematical model of the collector is defined as a function:
Figure 2: Schematic diagram of flat solar collector. Constructed flat solar collector. The solar hot water heating system
Based on solar energy absorbed by the plate of the collector, the heat loss coefficient of a collector, and heat absorbed by the fluid, the energy balance equation can be given by authors in the paper [19]:
where
Input values of the flat solar collector are presented in Tab. 1. We used these parameters of the collector to analyze the heat loss coefficient.
We used common loss factors for the solar collector; similar results were performed in [20]. We consider the flat solar collector for a system with a single glass coating. The analysis of the loss coefficient for the upper surface is carried out using the formula:
The convection coefficient between absorbed and glass coating
The dependence of the parameters Nu and Ra was calculated in [21].
where positive exponent means that only positive values are used if the term is negative, then use zero.
To calculate the thermal conductivity coefficient
The dependence of the temperature difference between the receiving surface and
Figure 3: The dependence of the temperature difference and
The correlation between
Figure 4: Dependence between
Figure 5: Dependence between
Fig. 6 shows the thermal conductivity
Figure 6: Dependence between the thermal conductivity,
After the thermal conductivity coefficient has been calculated, it is necessary to calculate the radiation coefficient from the glass coating to the absorber.
The calculation of the radiation coefficient between the glass surface and the air is calculated by the equation:
The calculation of the temperature of the glass surface is calculated through the equation:
The process of calculating the loss coefficient for the top surface is iterative and depends on the temperature of the glass surface, Fig. 7 shows the flow-chart of this calculation.
Figure 7: The flow-chart of calculation a maximum loss coefficient
During calculation, we get the following result for the maximum loss coefficient in Tab. 2.
The calculation results are the maximum loss coefficient from the collector plate to the ambient
Figure 8: Dependence between the average temperature of plate and top loss coefficient
The calculation of the loss coefficient for the lower surface of the collector is calculated using the formula:
where k is the heat conductivity of the insulating material (material–0.04
where U is calculated using Eq. (11), but for lateral isolation. Collector thickness: 0.1 m; insulation thickness: 0.01 m. The total loss coefficient will be calculated by summing all the coefficients.
4 Result of Simulation and Discussion
In the simulation performed by Matlab/Simulink, the model was taken from research [23]. The Matlab codes for of simulation a solar hot heating water system can be found at https://github.com/TimurMZh/Simulation-of-solar-flat-collector. The created model has been simulated a full solar hot heating water system: flat collector, storage tank, and pumps. Based on the mathematical model of the flat solar collector Eq. (1), a block schema has been developed by Simulink; the block schema of the collector model is shown in Fig. 9. The simulation result is the temperature of outlet fluid from the collector.
Figure 9: General block schema of the solar hot water heating system
We collected data for one year. The solar hot water heating system was installed in the Almaty city, in the foothills of Alatau; the height above the mean sea level 1088 meters. Input data: solar irradiance and ambient temperature. The result of the simulation is displayed in Fig. 10. Based on the figure, we can note that the solar hot water heating system will not provide an acceptable temperature in the flat collector and storage tank from December to March. The best result (high temperature) of the system will be provided from June to the middle of October. Fig. 10 shows a high correlation between solar irradiance and the collector and storage tank temperatures. Suppose the density of solar radiation is higher. In that case, the temperature is more monotonous, and the standard deviation of the temperature in the storage is close to the average value, for example, from July to August. The working process of the pump in the collector cycle depends on the reached temperature in the collector cycle, and we used the same limit values for activation of the pump as in research [18]. Fig. 10 showed that the pump turned to operate regularly in the summertime to transfer reached temperature from the collector into a storage tank. In wintertime, the pump does not work.
Figure 10: The result of simulation for one year
We scaled simulated data to 3 days (1st July–3rd July) to analyze simulated results in detail. The simulation result showed the maximum temperature of water in the storage tank 56.9°C in the Almaty region and in the collector is 60.5°C, the result presented in Figs. 11 and 12. As previously mentioned, a pump worked when the temperature in the collected reached a high temperature; it was reached in the daytime, and in the evening and night, a pump did not work. Also, we noted that the temperature in storage at night decreased monotonously without solar irradiation. Therefore, in Kazakhstan’s north region, we recommended insulating the storage tank or locating it in a building. The problem of saving heating in a storage tank is still open for the north of Kazakhstan, and we would like to consider it in another paper.
Figure 11: Ambient temperature, the temperature in collector and storage tank
Figure 12: Temperature in collector and storage tank with pump operational mode
Before performing a water heating simulation, it is necessary to analyze the technical characteristics of the flat solar collector, such as the heat loss of the collector. This analysis takes more time but allows simulating more accurate results. After the calculated values, we can proceed to energy analysis. The main purpose of the simulation model is to determine the installation limit values and identify weak installation points for improvement. The developed simulation model makes it possible to simulate the operation of a solar hot water heating system in various climatic conditions (temperature and solar radiation) and different individual system parameters to increase the system efficiency in various climatic conditions of Kazakhstan.
The simulation result showed that during the daytime, the temperature of water in the collector is 70°C, and during that period, it can be used at home. The main problem is heated water storage since the heated water is cooled at night. Another task we can define from the simulation analysis is to find a way to extend the work period of the system with high efficiency (April–October); it will be engineering, technical, or programming of controller solutions. This point needs to investigate and improved.
The next stage of the research will be the numerical simulation of the reservoir design. It is necessary to determine the design with high and low efficiency. We will determine the optimal system parameters for different climatic conditions using machine learning algorithms to adjust the parameters. We will conduct an experiment and a comparative analysis of simulated and experimental data.
The figure displays the collector, storage tank, and pump operation mode temperatures. We can see a correlation between ambient temperature and heating water temperature when the temperature in the collector reaches 40°C. The pump switches and starts to transfer a heated fluid to the exchanger coil inside the storage tank and receives the cooled fluid in the collector cycle for heating. The second reason for cooling fluid in the cycle is weather conditions (decrease in ambient temperature, sunset).
In our future research, we intend to concentrate on machine learning and deep learning algorithm to forecast the efficiency of a solar hot water system. Another research perspective is to apply Neural ODE (normalizing flow) for computational fluid dynamics (CFD) analysis of solar hot water heaters.
Funding Statement: This work was supported by the Ministry of Education and Science of the Republic of Kazakhstan BR10965172.
Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.
References
1. J. H Yoo, “Evaluation of solar hot water heating system applications to high-rise multi-family housing complex based on three years of system operation,” Energy and Buildings, vol. 101, pp. 54–63, 2015. [Google Scholar]
2. T. Chow, K. Fong, A. Chan and Z. Lin, “Potential application of a centralized solar water-heating system for a high-rise residential building in Hong Kong,” Applied Energy, vol. 83, no. 1, pp. 42–54, 2006. [Google Scholar]
3. V. Badescu, “Simulation analysis for the active solar heating system of a passive house,” Applied Thermal Engineering, vol. 25, no. 17-18, pp. 2754–2763, 2005. [Google Scholar]
4. M. Siddiqui, and S. Said, “A review of solar powered absorption systems,” Renewable and Sustainable Energy Reviews, vol. 42, pp. 93–115, 2015. [Google Scholar]
5. S. Kalogirou, “Solar thermal collectors and applications,” Progress in Energy and Combustion Science, vol. 30, no. 3, pp. 231–295, 2004. [Google Scholar]
6. R. Maldonado, E. Huerta, J. Corona, O. Ceh, A. Leon et al., “Design and construction of a solar flat collector for social housing in México,” Energy Procedia, vol. 57, pp. 2159–2166, 2014. [Google Scholar]
7. C. Nogueira, M. Vidotto, F. Toniazzo and G. Debastiani, “Software for designing solar water heating systems,” Renewable and Sustainable Energy Reviews, vol. 58, pp. 361–375, 2016. [Google Scholar]
8. S. Farahat, F. Sarhaddi and H. Ajam, “Exergetic optimization of flat plate solar collectors,” Renewable Energy, vol. 34, no. 4, pp. 169–1174, 2009. [Google Scholar]
9. I. Luminosu and L. Fara, “Determination of the optimal operation mode of a flat solar collector by exergetic analysis and numerical simulation,” Energy, vol. 30, no. 5, pp. 731–747, 2005. [Google Scholar]
10. S. Park, A. Pandey, V. Tyagi, and S. Tyagi, “Energy and exergy analysis of typical renewable energy systems,” Renewable and Sustainable Energy Reviews, vol. 30, pp. 105–123, 2014. [Google Scholar]
11. Y. Amirgaliyev, M. Kunelbayev, B. Amirgaliyev, A. Kalizhanova, A. Kozbakova et al., “Mathematical justification of thermosyphon effect main parameters for solar heating system,” Cogent Engineering, vol. 7, no. 1, no. 5, p. 1851629, 2020. [Google Scholar]
12. M. Yavuz and N. Özdemir, “Numerical inverse Laplace homotopy technique for fractional heat equations,” Thermal Science, vol. 22, pp. 185–194, 2017. [Google Scholar]
13. A. M. S. Mahdy and E. S. M. Youssef, “Numerical solution technique for solving isoperimetric variational problems,” International Journal of Modern Physics, vol. 32, pp. 2040–2057, 2021. [Google Scholar]
14. A. M. S. Mahdy, K. Lotfy, W. Hassan and A. A. El-Bary, “Analytical solution of magneto-photothermal theory during variable thermal conductivity of a semiconductor material due to pulse heat flux and volumetric heat source,” Waves in Random and Complex Media, vol. 31, no. 6, p. 107, 2021. [Google Scholar]
15. M. Fayz-Al-Asad, M. Yavuz, M. Alam, M. Sarker, M. Alam et al., “Influence of fin length on magneto-combined convection heat transfer performance in a lid-driven wavy cavity,” Fractal and Fractional, vol. 5, no. 3, pp. 107, 2021. [Google Scholar]
16. Y. Amirgaliyev, W. Wójcik, M. Kunelbayev, T. Merembayev, D. Yedilkhan et al., “Theoretical prerequisites of electric water heating in solar collector-accumulator,” News of the National Academy of Sciences of the Republic of Kazakhstan, Series of Geology and Technical Sciences, vol. 6, pp. 54–63, 2019. [Google Scholar]
17. G. Morini and S. Piva, “The simulation of transients in thermal plant. Part I: Mathematical model,” Applied Thermal Engineering, vol. 27, no. 11–12, pp. 2138–2144, 2007. [Google Scholar]
18. G. Morini and S. Piva, “The simulation of transients in thermal plant. Part II: Applications,” Applied Thermal Engineering, vol. 28, no. 2–3, pp. 244–251, 2008. [Google Scholar]
19. W. Hao, Y. Lu, Y. Lai, H. Yu and M. Lyu, “Research on operation strategy and performance prediction of flat plate solar collector with dual-function for drying agricultural products,” Renewable Energy, vol. 127, pp. 685–696, 2018. [Google Scholar]
20. F. Jafarkazemi and E. Ahmadifard, “Energetic and exergetic evaluation of flat plate solar collectors,” Renewable Energy, vol. 56, pp. 55–63, 2013. [Google Scholar]
21. K. Hollands, T. Unny, G. Raithby and L. Konicek, “Free convective heat transfer across inclined air layers,” Free Convective Heat Transfer Across Inclined Air Layers, vol. 98, pp. 189–193, 1976. [Google Scholar]
22. H. Tabor, “Radiation, convection and conduction coefficients in solar collectors,” Bulletin of the Research Council of Israel, Sector C, vol. 6, no. 3, pp. 155–176, 1958. [Google Scholar]
23. Y. Amirgaliyev, T. Merembayev and M. Kunelbayev, “Dynamic simulation of a solar hot water heating system for Kazakhstan climate conditions,” in 2018 14th Int. Conf. on Electronics Computer and Computation (ICECCO), IEEE: Kaskelen, Kazakhstan, vol. 1, pp. 206–212, 2018. [Google Scholar]
 | This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. |