Next Article in Journal
Analysis of Electricity Supply and Demand Balance in Residential Microgrids Integrated with Micro-CAES in Northern Portugal
Previous Article in Journal
Heat Transfer Performance Factors in a Vertical Ground Heat Exchanger for a Geothermal Heat Pump System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solar Chimney Operation Variant

1
Institute of Environmental Engineering, University of Zielona Góra, 65-417 Zielona Góra, Poland
2
Institute of Mathematics, University of Zielona Góra, 65-417 Zielona Góra, Poland
*
Author to whom correspondence should be addressed.
Energies 2024, 17(19), 5004; https://doi.org/10.3390/en17195004
Submission received: 23 August 2024 / Revised: 23 September 2024 / Accepted: 6 October 2024 / Published: 8 October 2024
(This article belongs to the Section A2: Solar Energy and Photovoltaic Systems)

Abstract

:
This paper presents a solar chimney that acts as a heat accumulator. It is based on its alternating charging (melting of the phase change material—PCM) and discharging (solidification), which helps to save energy and ensures stable operation of the solar chimney. In this paper, special attention has been paid to the heat dissipation process (solidification of the PCM). The theoretical model of solidification has been solved in an original way. This paper presents a new simple theoretical model for the solidification of the PCM on a flat plate and presents the results of numerical tests. The theoretical model presents a method for determining the heat transfer coefficient at the solidification front of the PCM. In addition, the heat transfer coefficient from the flowing air to the outer surface of the solidifying front plate was determined experimentally in an original way. The heat transfer coefficient values resulting from the experiments may be employed in order to calculate the heat transfer coefficient for air flowing through the slot of the collector in the solar chimney. The calculated value of the heat transfer coefficient was 18.55 W/m²K.

1. Introduction

The airflow in a solar chimney, which consists of a flat collector with an air gap and a tall central cylinder, is a significant research subject among scientists in the field of heat transfer. A solar chimney can have multiple uses, which include inducing ventilation in buildings, powering an air turbine inside the chimney for electricity generation, and facilitating drying processes. Numerous articles have been published on solar chimneys, mainly discussing experimental research on particular designs and their theoretical analysis. Heat flow theory is typically applied to describe and study actual systems. Some articles also feature analytical models of the movement of air and heat transfer in the solar collector-chimney system.
A solar chimney is a device that functions as a solar energy power plant. Initially, it transforms solar energy to thermal energy through its solar collector. Afterward, it converts this energy into kinetic energy to propel the flowing air up the chimney. Finally, it harnesses the power generated by the wind turbine and generator to produce electricity.
The use of solar energy in solar chimneys helps to reduce greenhouse gas emissions.
The operation of a solar chimney relies solely on solar power. If the sun is not powerful enough, the solar chimney may not work effectively enough to make electricity. Adding a phase change material (PCM) as part of the structure can solve this issue by storing extra heat during periods with more sun and releasing it during periods with less sun, making continuous operation possible. Therefore, suggesting various ways of operating a solar chimney is necessary and appropriate.
One of the earliest works on solar collector and chimney systems is Isidoro Cabanyes’ 1903 proposal [1]. His work looked at engine design and proposed the introduction of a solar air heater attached to a house with a chimney. The inside of the chimney was fitted with an impeller to generate electricity. In 1970, the first combined solar collector and chimney system was designed by Schlaich in Manzanares, Spain [2]. The installation cost is considerable due to the necessity of incorporating a chimney and collector. Furthermore, the power conversion efficiency is notably low. In order to overcome these difficulties, researchers have put forth a number of innovative proposals with the aim of enhancing performance. One potential solution is the selection of an appropriate geometry.
Solar chimney systems have a significant impact on energy supply in nations with abundant solar radiation. Many experts have researched this technology in recent years, specifically exploring solar thermal technology. This type of technology has potential uses globally [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24].
The European Parliament’s climate and energy policy, which aims to reduce greenhouse gas emissions by 2030, has set a target of at least 27% of energy consumed in the EU coming from renewable energy by 2030. In order to reach this target, it is necessary to look for new technical solutions that enable the use of renewable energy sources [25,26]. A good example of this is energy from the sun, which can be used for drying, for example, with a solar chimney.
The article [11] examines the potential of solar chimney systems in combination with other renewable or conventional energy systems. It presents an overview of compact solar chimneys, solar panels, solar ponds, and geothermal energy, as well as an analysis of solar chimney systems, integrated with power plants.
In work [12], a three-dimensional hybrid solar chimney with an integrated external heat source was developed. Solar energy is used for continuous power generation using flue gas channels as a power source in the collector. Extensive numerical studies are also presented. The aim of this article is to theoretically analyze the solidification model of the PCM material on the slab, to analyze the operating conditions of the solar chimney, and to calculate the air heat transfer coefficient.
The origins of solar energy can be traced back 3000 years. Especially deserving of attention are structures from the Middle East. For instance, solar ventilation systems were constructed as early as three thousand years ago in Iran. The main energy driving the movement within these systems is wind power. If wind power is inadequate, solar energy is utilized to heat the walls. This results in a decrease in the density of indoor air, leading to free convection.
Additionally, the use of phase transformation in these constructions is noteworthy. In the building, there were watercourses or reservoirs on the lower floors. Water from the surface of these evaporated, lowering the inside temperature.
Isidoro Cabanyes [4] made an early effort to use solar collector and chimney systems in 1903. He proposed a solar air heater attached to a house with a chimney as part of a motor design. An impeller was installed inside the chimney to produce electricity [4]. Schlaich conceptualized the solar collector and chimney arrangement in 1970, at Manzanares, Spain. The chimney had a height of 195 m and a diameter of 10 m, and the collector’s radius was 122 m. The prototype was a miniature experimental version. In its final form, it produced approximately 50 kW [2]. In 1995, Schlaich suggested an innovative method to create a solar collector stack, where the gadget obtains power from direct and diffused radiation, as well as from the soil’s heat accumulation. The collector and stack were fashioned utilizing basic elements and techniques [15]. In the paper [27], a numerical study was presented, which examined the impact of varying chimney shapes, including round, convergent, divergent, abrupt contraction, and abrupt expansion, on the performance of a solar chimney power plant (SCPP). Moreover, the parametric effect of varying chimney divergence angles and ground absorber inclination angles on the performance of a solar chimney power plant (SCPP) was also examined. The optimal divergence angle was identified. The findings of their study were benchmarked against those of the classical Manzanares power plant, which demonstrated an increase in energy production of up to 80% (92 kW) relative to the Manzanares plant.
The objective of the paper [28] is to review SCPPs in order to gain insight into the underlying physics of the flow and to identify the key parameters and their values that are most pertinent to the optimal design of an SCPP. Additionally, hybrid models with further modifications that are suitable for commercial production of large-scale SCPPs are included.
In sunnier countries, solar chimneys have a vital role in supplying energy. Therefore, in recent decades, loads of researchers have delved into solar chimney systems, exploring them as a way of producing solar thermal energy that could be used globally. Bernardes and colleagues (2012) conducted a thorough analysis of different models to determine the performance of a solar collector stack. The aim was to estimate the output of solar chimneys and to examine how varying environmental conditions and geometric dimensions affect the output of solar chimneys. The mathematical model was verified through experiments and used to predict how commercial solar chimneys will perform. Stack height, turbine pressure drop and diameter features were discovered to be very important [4].
The paper [29] reviews the thermodynamic principles underlying the energy balance equations of the plant and those needed for the various components, as these principles are useful for design calculations. Furthermore, exergy analysis and cost estimation are also discussed. A review of the global status of power plants over the past two decades reveals that the system capacity potential ranges from 0.053 to 27 MW. Any mathematical analysis and research gaps were summarized, as they are useful for informing future research. Additionally, the factors affecting system performance, including area transfer coefficients, pressure profile, and turbine pressure drop, were provided. The main challenges encountered at the plant were identified as lower thermal efficiency, site requirements for a commercial plant, and the construction of a tall chimney. In conclusion, recommendations and prospects for the future of the plant were presented.
An intriguing paper is that of Saleh et al. [30], which presents a review of current literature examining the performance of the primary components of solar chimney power plants and the most effective approaches to enhance the system’s performance through diverse methodologies. These include the incorporation of the system with heat storage materials and phase change materials (PCMs), as well as its integration with other technologies. Additionally, the paper discusses the integration of solar chimney power plants with external heat sources and engineering solutions for tower chimney design, with the objective of optimizing the performance of solar chimney power plants.
In the article by Lipnicki et al. [7], the authors suggest a new way to solve the problem of air flow through solar collectors and chimneys. They use the classical system of conservation equations (momentum, mass, and energy) and similarity theory. This method allows the making of a solar chimney. The theoretical analysis was compared to experiments conducted on current solar chimney and solar collector systems. The experimental and theoretical studies were found to be sufficiently consistent. The flow of heat in a solar chimney was outlined at length, using dimensionless numbers, including Reynolds, Grashof, Galilean, Biot, and Prantdl numbers.
The paper deals with the solidification of a PCM material in the form of a flat disc with radius R and height H. The paper presents heat accumulation by melting the PCM material and heat release by the solidification of the PCM material, along with a new theoretical model for the solidification of PCM on a flat plate, presenting the results of numerical tests. The theoretical model presents a method for determining the heat transfer coefficient at the solidification front of the PCM. In addition, the heat transfer coefficient from the flowing air to the outer surface of the solidifying front plate was determined experimentally in an original way. Chapter two presents a theoretical model for the solidification of a PCM material in the form of a flat plate, introduces a mini solar chimney, and shows its construction and the variation possibilities of the operation of this chimney. Chapter three presents the results of numerical calculations of the solidifying slab of PCM material, and the heat transfer coefficient is determined and presented.
This work solves a new problem about solar chimney operation. It looks at heat accumulation in PCM using phase transformation. This problem has not been considered in the scientific literature.

2. Materials and Methods

2.1. Theoretical Model of Solidification of PCM Plate

Figure 1 shows the one-dimensional solidification model in the z direction, directed opposite to the acceleration vector of gravity g , for a flat liquid layer of PCM material of thickness H , thermally insulated from the top and cooled from the bottom by air flowing at a temperature T a and at a velocity u . A very small gap, called the contact layer, is formed between the lower plate made at the temperature T W of a very good conductor (copper) and the solidified PCM material, which creates additional resistance to heat flow q . The resulting gap in contact with the solidified layer has a temperature of T C . A solidified layer of thickness δ develops on the lower plate. The non-solidified liquid layer of thickness H δ decreases with time t and is assumed to be immovable (no free convection).
The lack of free convection in the liquid layer justifies the expected temperature distribution. The lower liquid layers have a lower temperature and therefore a higher density, which does not favor the formation of free convection.
Such a simplified description of the solidification phenomenon justifies its relatively long duration [5,31].
h F T ¯ L T F + ρ S L d δ d t = k S T F T C δ = h C O N T C T W = h T W T a
The first term of the above equation represents the sum of the heat fluxes: transferred from the liquid with the transfer coefficient h F , where T ¯ L is the average temperature of the liquid, and T F represents the solidification temperature and the heat of phase change, where ρ S is the density of the solidified layer, and L denotes the heat of solidification. The second term of the equation represents the heat flux flowing through the solidified layer. The third term represents the heat flux flowing through the contact layer, with the transfer coefficient h C O N . The fourth term represents the heat flux transferred by the cooling air flow with the transfer coefficient h .
The heat stream absorbed from the superheat liquid PCM is equal to the capacitive heat change in the plate of the liquid in the duration phenomena
h F δ T ¯ L T F = d d t H δ ρ L c L T ¯ L T F
The above equation is derived from the heat balance of a PCM board.
The subject of this work is the solidification of a liquid placed in a plate space. The space is shown in Figure 1. In the space, the solidification flat front with thickness δ t moves from the flat cooling surface to the outer surface of the thickness H . The inner surface with z = 0 is cooled by the flowing coolant air. The heat stream of solidification q ˙ flows into the air.
The quasi-stationary equation describing in flat space the temperature fields both in the liquid space δ t z < H and in the solidified layer space 0 z δ t describes the same quasi-stationary differential equation:
d 2 T d z 2 = 0 .
The boundary conditions are equal:
-
For the solidified layer in the space 0 z δ t :
T = T F f o r z = δ ,   and   T = T C   f o r   z = 0 .
-
For the stand still liquid in the space δ t z < H
T = T F   and   T ¯ L t = δ H T z d z H δ   f o r   z = δ ,
where the average temperature of the liquid is given by
T ¯ L T F = δ H T z T F d z H δ .
In the solidified layer 0 z δ t , the temperature distribution and the stream heat are equal:
T = T C + T F T C z δ
q ˙ S = k S T z ^ z ^ = δ = k S T F T C δ .
In the liquid space δ t < z < H , the temperature distribution and the stream heat on the solidification front are equal:
T = 2 T ¯ L T F H δ z + T F H + δ H δ 2 δ H δ T ¯ L ,
T δ = T F H δ H δ = T F
q ˙ = k L T z z = δ = k L 2 T ¯ L T F H δ .
The stream of heat absorbed from the liquid to the solidification front is equal:
q ˙ = h F δ T ¯ L T F .
Comparing heat streams, the heat transfer coefficient on a cylindrical surface with a thickness δ on the liquid side is equal:
h F δ = 2 k L H δ .
After the appropriate transformation of equations (1), one equation of heat flux balance was obtained:
h F T ¯ L T F + ρ S L d δ d t = k S h C O N h h C O N k S + h k S + h C O N h δ T F T a .
The system of differential Equations (1)–(3) can be replaced by the following dimensionless system:
2 k ~ B 1 δ ~ θ ¯ L + d δ ~ d τ = B i a B i C O N B i a + B i C O N + B i a B i C O N δ ~ ,
2 a ~ 1 δ ~ θ ¯ L = S t e d d τ 1 δ ~ θ ¯ L ,
by introducing dimensionless quantities into the above equation system: thickness of the solid layer δ ~ = δ / H , time τ = F o   S t e , the average overheat temperature of the liquid θ ¯ L = T ¯ L T F / T L 0 T F , where T L 0 is the initial temperature of liquid PCM material, F o = a S t / H 2 is the Fourier number, S t e = c S T F T a / L is the Stefan number,   a S = k S / ρ S c S is the thermal diffusion coefficient of the solidified layer, a L = k L / ρ L c L is the liquid thermal diffusion coefficient, B i C O N = h C O N H / k S is the Biot number in the contact layer, B i a = h H / k s is the Biot number on the surface of the plate, B = T L 0 T F / T F T a is the liquid overheat parameter, a ~ = a L / a S is the ratio of liquid and solid heat diffusion coefficients, k ~ = k L / k S is the ratio of the heat conduction coefficients of the liquid and solid.
The system of differential equations above satisfies the initial conditions:
δ ~ = 0 ,               θ ¯ L = 1   f o r   τ = 0 .
These conditions can be found numerically.
For small overheating parameters B 0 (see the figure in Section 3.1), which occur in reality, there is no significant difference in the change of the solidified layer thickness; therefore, the simple equations below can be used in the solidified layer thickness and in the solidification velocity calculations.
Hence, Equation (14) leads to simple equations of the solidified layer thickness [5]:
δ ~ = 1 B i a + 1 B i C O N + 1 B i a + 1 B i C O N 2 + 2 τ
and the solidification velocity is given as follows:
d δ ~ d τ = 1 1 B i a + 1 B i C O N 2 + 2 τ .
Equations (17) and (18), after performing the transformations, allow us to determine the thickness of the solidified layer ( δ ) and the heating power ( Q ˙ ) :
δ = H   δ ~   and   Q ˙ = ρ S L a S H S t e   F   d δ ~ d τ .

2.2. The Mini Solar Chimney

The mini solar chimney that inspired this work was created as part of the project titled “Cooperation of scientific partners in education and knowledge exchange in the field of energy storage technology and energy efficiency in the SNB region”, co-funded by the INTERREG VA Brandenburg—Poland 2014—2020 program. The mini solar chimney is located in the premises of the Centre for Education in Energy Storage Technology in the field of thermal energy storage and energy efficiency, located in Nowy Kisielin, in the Science and Technology Park of the University of Zielona Góra. Figure 2 shows a photo of the mini solar chimney.
Figure 3 shows a schematic diagram of a mini solar chimney, which consists of two main parts: a solar collector and a vertical thermally insulated pipe placed in the center of the panel.
Due to its size, the mini-chimney’s performance is not significant. However, obtaining practical research results using thermodynamic similarity principles may be useful for real, large-scale solar chimneys. The theoretical and practical outcomes of the small chimney can aid in the development of real chimneys, and it is clear that the cost of experimental studies on a small chimney is not high compared to the cost of working on large solar chimneys. Solar chimneys are not profitable in Poland, but they can be used as an example for other countries with more sunlight.
The solar panel consists of two flat disks horizontally positioned at a height of H and parallel to each other, with a space in between, called an air gap h. The top plate is made of metal and has a specific diameter of 2R = 2000 mm, an air gap width of h = 80 mm, and a distance between two parallel flat disks of H = 470   m m . In the middle of the concrete surface sat a vertically positioned, heat-resistant tube with an inner diameter of d = 100 mm and a height of Lc = 2500 mm (Figure 3). Depending on the research requirements, the pipe could be substituted with a longer one. Moreover, the illustration displays further measurement and testing components: a water heater with manifold and metal coils, spheres containing PCM material, and a flat PCM material plate. The standpipe may have a turbine power generator, for example, attached within, and an axial fan installed at the very top to reinforce the airflow. Sensors are appropriately placed within the pipe to gauge the airflow and temperature.
Due to the heating of the air from the sun, radiant heat, a water heater, or the solidification of PCM material within the chimney, an air flow known as the “chimney effect” emerges in the solar collector—chimney (standpipe) system. An instrument called a wing anemometer was utilized to measure the air velocity within the chimney, while a thermometer was employed to measure the temperature. For example, the air velocity measurements using the wing anemometer were compared to those made with the Pitot tube in the chimney. In addition, we measured the pressure drop at the junction between the collector and standpipe with a Recknagel differential pressure gauge. The pressure drop at the junction between the solar collector and the standpipe represents the local loss of airflow, as defined in the context of fluid mechanics. These results were the subject of a previous paper by the authors [7]. We used a radiation sensor at a nearby meteorological station to measure solar radiation.
All the measurements of the solar chimney’s functioning were saved using the Building Management System (BMS), as shown in Figure 4. The BMS tracks the temperature at two points of the chimney, the air flow in the chimney, and the temperature at the top and bottom of the collector container, as well as the power that is consumed and the active power.
In addition, external parameters can be read from a nearby meteorological station, including the temperature, the total and diffuse solar radiation, humidity, as well as wind speed and direction.
The mini solar chimney has been designed and manufactured in such a way that it can be used to test different configurations and heat flows. The mini solar chimney allows the multi-dimensional operation of such a system to be represented, i.e., the use of the solar heat flow, the use of the heat exchanger, and the use of the heat flow of the PCM material. This is possible thanks to a system of three exchangers (coils) installed in the mini solar chimney. The following are the heat exchangers:
-
The first one is placed on the bottom plate of the solar collector, within the space of the flowing air. In this case, the air in the gap of the solar collector is heated by the heat of the solar radiation.
-
The second one is placed in the tank space under the solar collector, where the PCM material is placed, which works in two cycles:
-
Melting of the PCM material by the warm air flowing in the collector gap (charging of the storage tank);
-
Solidification of the PCM material by the cold air flowing in the collector gap (discharging the heat storage tank).
-
The third is at the bottom of the tank, under the solar collector. In this case, the air in the collector gap is heated by the heat from the heat exchanger.
-
The solar chimney co-operates with a round flat plate of phase change material (PCM), which solidifies when air flows through the collector gap. The solidification process of the liquid PCM material is shown in Figure 5.
-
The air flowing in the collector gap has a velocity u and an average temperature T ¯ p lower than the solidification temperature T F of the PCM material. As a result of the process of absorbing heat q   ˙ from the liquid, it is in the form of a cylindrical plate, with a radius R and a thickness H , where the solidification process of the PCM material takes place.

3. Results and Discussion

This section is divided into three subsections: 3.1. Solidification of PCM Plate, 3.2. Determination of the Heat Transfer Coefficient, 3.3. Examples of Mini-Chimney Operation Processes.

3.1. Solidification of PCM Plate

Numerical solutions of the integral in Equations (14) and (15) for different parameters are presented in Figure 6a,b.
Figure 6a shows the thickness of the solid layer over time, Figure 6b depicts the average overheat temperature of the liquid over time. The presented figures indicate that the thickness of the slid layer increased, and the average overheat temperature of the liquid decreased over time. The influence of the overheating parameter B on the temperature distribution is insignificant.
Figure 7 shows the thickness of the solid layer over time for different Biot numbers ( B i a ) on the air side for B = 0 .
The solidified layer became thicker over time, with more the Biot numbers on the surface of the plate B i a .

3.2. Determination of the Heat Transfer Coefficient

A very important thermodynamic parameter for the air flowing in the solar collector gap is the heat transfer coefficient ( h ) . In order to determine it, original experimental research was carried out on a test bench constructed by the co-author of the paper and shown in Figure 8a. The diagram of the test bench is shown in Figure 8b. The test stand consists of a tank with a radius of R = 63.5   m m and a height of H = 104   m m , and it is made of acid-resistant steel, with a wall thickness of 1.5   m m and a mass of m s = 0.695   k g . The mass of the water in the tank is m = 1   k g , and the temperature is T w . The results of the measurements are recorded for different initial temperatures T w p . The water tank is installed in a 247   m m diameter duct (Figure 8), in which the air flows at a temperature of T p . There is a heater under the duct. The heater heats the air at different flow rates ( u ) .
The temporary heat balance of the tank is described by the formula
d Q = m c w + m s c s d T w = h T p T w F d t
where h —the heat transfer coefficient, c w —the specific heat of the water, c s —the specific heat of the solidified layer, F —the side surface of the tank ( F = 2 π R H ) .
By separating the variables, an ordinary differential equation was obtained:
m c w + m s c s d T p T w T p T w = h F d t ,
whose solution is given when the initial condition is satisfied: t = 0 ,   a n d   T w = T w p is the equation
T w = T p T p T w p · e x p A t ,
where A = h F m c w + m s c s .
The average simplified heat balance for a container of water over time t is equal to
m w c w T w k T w p + m s c s T w k T w p = h T p T ¯ w F t ,
where T ¯ w = T w p + T w k 2 is the average water temperature, T w p is the initial water temperature, T w k is the final water temperature, F = 2 π R H is the side surface of the tank, and t is the time.
The heat transfer coefficient, which also takes into account the heat absorbed by the tank wall, is
h = m w c w T w k T w p + m s c s T w k T w p T p T ¯ w F t = 18.55 W m 2 K .
It is assumed that the temperature of the water is equal to the temperature of the vessel wall.
By substituting the values of the h parameters from Equation (24) into parameter A , the time variation of the water temperature is described by the equation:
T w = T p T p T w p e A t = 62.9 62.9 31.94 e 0.000184   t
Equation (25) shows the variation of the water temperature with time as a function of the ambient air temperature ( T p ) and the initial temperature of the water ( T w p ) . The constant A represents the ratio of the heat transfer coefficient and the external surface area of the tank to the sum of the products of the mass of water and its specific heat and the mass of the solidified layer and its specific heat.
As can be seen in Figure 9, the temperature profile as a function of time is linear (straight line). Therefore, the simplified method presented in this paper can be used to calculate the heat transfer coefficient on the side surface of the tank.
The values of the heat transfer coefficient obtained from the experiments carried out can be used to calculate the heat transfer coefficient for the air flowing in the collector slot in the solar chimney. An additional assumption must then be made that the curvature of the cylindrical surface 1 / R is small. Heat transfer coefficient values h (Figure 10). They increase as the velocity of the flowing air increases. The test results obtained can be used to determine the conditions of heat transfer in the flat gap of a solar collector.
The results of the experimental and theoretical studies are presented in a graph (Figure 10).

3.3. Examples of Mini-Chimney Operation Processes

In a research study [6], a solar chimney was examined both theoretically and experimentally, and the optimum operating conditions were established based on the authors’ original research (refer to Figure 11 and Figure 12).
The data obtained from the operation of the mini solar collector, comprising measurements of air flow speed in the chimney, wind speed, sunlight, and air temperature, are presented in tabular form in Table 1.
The test results presented in Table 1 demonstrate that the air velocity within the chimney is contingent upon solar radiation, air temperature, and external wind speed. In general, the air velocity in the chimney increases with increasing solar radiation. The presence of dust in the chimney is a significant parameter affecting the speed of the air flow. The chimney contains dust on the top plate of the solar collector; therefore, it is essential to ensure that the collector plate is kept clean. The mini solar chimney, as the name suggests, achieves a low output; however, the results obtained from its operation can be useful for chimneys with large external dimensions.

3.3.1. Heating the Air in the Collector Gap with the Heat of the Solidification Stream

The influence of the heat transfer coefficient ( h ) on the thickness of the solidified layer is shown in Figure 13. The curves shown in the figure are the graphical form of Equation (16).
As the value of the heat transfer coefficient increases, the thickness of the solidified layer increases more rapidly.

3.3.2. The Heat Flow from the Flat Lower Circular Plate in the Collector Gap

The PCM material was used to calculate the heat flow from the flat-bottomed circular plate in the collector slot.
The PCM used was Rubitherm®RT64HC, with the following thermophysical parameters: the solidification heat, L = 250   k J / k g ; the thermal conductivity coefficient, k s = 0.2   W / m K ; the density of the solid, ρ S = 880   k g / m 3 ; the density of the liquid ρ L = 780   k g / m 3 ; the specific heat capacity c p = 2.0   k J / k g K ; the heat diffusion coefficient, a S = 1.25 · 10 7   m 2 / s ; density, ρ s = 880   k g / m 3 [24].
The parameters of the PCM plate were the following: the radius of the plate, R = 1   m ; the plate temperature, T w = 64   ° C ; the surface area of the plate, F = 1   m 2 .
The air parameters in the gap were as follows: the air velocity, u = 1.5   m / s ; the average air temperature, T a ¯ = 20   ° C ; the coefficient of kinematic viscosity, ν = 15.06 · 10 6 m2/s; the heat conductivity of the air, k a = 2.59 · 10 2   W / ( m K ) ; the Prandtl number, P r = 0.703 .
The growth rate of the solidified layer d δ t / d t (indirectly representing the heat flux density) and the increase in thickness of the solidified layer δ t over time according to Equation (16) are shown in Figure 14 and Figure 15.
As can be seen in Figure 14, the thickness of the solidified layer increases during the solidification process, while the solidification rate decreases. The highest solidification rate occurs at the beginning of the process.
The heat released as the phase of the PCM material changes over time is shown in Figure 16. This process can be thought of as discharging the heat storage. Figure 16 also shows the increase in thickness of the solidified layer over time.
As depicted in the figure above, both the total amount of heat released and the thickness of the solidified layer increase with time.
The total heat flux released in the solidification process, calculated according to equation 12 for a collector plate with the surface of F = 1 m 2 , is shown in Figure 16.
As can be seen in Figure 16, the total heat flux (power) decreases during the process. This is due to the increase in thickness of the solidified layer and therefore, the increase in thermal resistance of the solidified layer.

4. Conclusions

The multi-variant operation of the mini solar chimney, as presented in this paper, demonstrates in a straightforward manner that external heat can be converted into air movement in a number of different ways. Solar radiation, or radiation from a radiator, can be employed directly to heat the air flowing in the gap between the solar collectors. This particular configuration of a solar chimney represents a classic example of such a device, and it is comprehensively documented in the technical literature. In reference [6], a comprehensive and insightful theoretical analysis of the solar chimney operation was conducted.
The process presented in this work on solar chimney operation is comprised of two stages. The initial stage entails utilizing solar radiation to impart kinetic energy to the airflow within the solar chimney, while simultaneously facilitating the melting of the PCM material. The melting of the PCM material is associated with the charging of excess energy. In the second stage, in the absence of solar radiation, the energy accumulated in the PCM material is returned to the air flowing in the chimney as a result of the PCM material solidification process. This work provides a detailed analysis of the solidification process. The two stages are repeated cyclically. The proposed method of chimney operation is continuous and stable. During the nocturnal period, when solar radiation is absent, the solar chimney operates continuously, maintaining a constant temperature throughout the heat exchange process.
The solution presented in this work is highly beneficial for environmental protection. Furthermore, the work presents an original method of determining the heat transfer coefficient between the flowing air and the flat plate of the solar collector.
The achievements of this work include experimental studies on the heat transfer coefficient between the air flowing in the gap and the surface of the collector plate. There is a lack of theoretical studies in the scientific literature on the heat transfer coefficient for axially symmetric air flows over flat plates, which should be the focus of further theoretical and experimental research.
It is best to build solar chimneys in countries with a warm climate and plenty of sunlight. In our country, the climatic conditions are rather unfavorable for building large chimney construction investments. However, the theoretical research carried out on a mini-chimney and the results obtained can be transferred to large chimneys, using a similarity scale.

Author Contributions

Conceptualization, M.G., Z.L., T.M. and P.G.; methodology, M.G., Z.L., T.M. and P.G.; validation, M.G. and Z.L.; formal analysis, M.G., T.M. and Z.L.; investigation, M.G., P.G. and Z.L.; resources, M.G., P.G. and Z.L.; data curation, M.G., Z.L., T.M. and P.G.; writing—original draft preparation, M.G. and Z.L.; writing—review and editing, M.G. and Z.L.; visualization, M.G. and Z.L.; supervision, M.G. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of this study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Nomenclature

a S thermal diffusivity of the solidified layer, m 2 / s
a L thermal diffusivity of the liquid, m 2 / s
c S specific heat of the solidified layer, J / k g K
c L specific heat of the liquid, J / k g K
g gravity of acceleration, m / s 2
h heat transfer coefficient at the plate, W / m 2 K
h F heat transfer coefficient at the solidification front, W / m 2 K
h C O N heat transfer coefficient at the contact layer, W / m 2 K
k L heat conductivity of liquid, W / m K
k S heat conductivity of the solidified layer, W / m K
k C heat conductivity of cooling liquid, W / m K
H plate PCM material height, m
L latent heat of liquid, J / k g
k heat conductivity, W / ( m K )
L C chimney length, m
q ˙ heat flux, W / m 2
Q ˙ heat transfer rate, W
t time, s
T temperature, K
u velocity, m / s
z coordinate, m
Greek symbols
δ thickness of the solid layer, m
ρ density, k g / m 3
θ ¯ L overheat temperature of the liquid
τ time

References

  1. Asnaghi, A.; Ladjevardi, S.M. Solar chimney power plant performance in Iran. Renew. Sustain. Energy Rev. 2012, 16, 3383–3390. [Google Scholar] [CrossRef]
  2. Haaf, W.; Friedrich, K.; Mayr, G.; Schlaich, J. Solar Chimneys—Part I: Principle and Construction of the Pilot Plant in Manzanares; Taylor & Francis: Abingdon, UK, 1983; pp. 3–20. [Google Scholar]
  3. Bernardes, M.A.d.S.; Voß, A.; Weinrebe, G. Thermal and technical analyses of solar chimneys. Sol. Energy 2003, 75, 511–524. [Google Scholar] [CrossRef]
  4. Cabanyes, I. Proyecto de Motor Solar. Energ. Eléctrica—Rev. Gen. Electr. Sus Apl. 1903, 8, 61–65. [Google Scholar]
  5. Lipnicki, Z. Dynamics of Liquid Solidification, Thermal Resistance of Contact Layer; Mathematical Engineeringm; Springer International Publishing: Cham, Switzerland, 2017. [Google Scholar]
  6. Lipnicki, Z.; Gortych, M.; Staszczuk, A.; Kuczyński, T. The theoretical analysis of mass and energy flow through solar collector—Chimney system. Civ. Environ. Eng. Rep. 2017, 24, 117–131. [Google Scholar] [CrossRef]
  7. Lipnicki, Z.; Gortych, M.; Staszczuk, A.; Kuczyński, T.; Grabas, P. Analytical and Experimental Investigation of the Solar Chimney System. Energies 2019, 12, 2060. [Google Scholar] [CrossRef]
  8. Schlaich, J. The Solar Chimney. Electricity from the Sun; Edition Axel Menges: Stuttgart, Germany, 1995. [Google Scholar]
  9. Mehling, H.; Cobeza, L. Heat and Cold Storage with PCM; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  10. Lipnicki, Z.; Rubaszewski, A. Współpraca akumulatora ciepła z układem chłodzenia silnika spalinowego. Arch. Motoryz. 2000, 3, 177–186. [Google Scholar]
  11. Ahmed, O.K.; Algburi, S.; Ali, Z.H.; Ahmed, A.K.; Shubat, H.N. Hybrid solar chimneys: A comprehensive review. Energy Rep. 2022, 8, 438–460. [Google Scholar] [CrossRef]
  12. Yazdi, M.H.; Solomin, E.; Fudholi, A.; Sopian, K.; Chong, P.L. Numerical analysis of the performance of a hybrid solar chimney system with an integrated external thermal source. Therm. Sci. Eng. Prog. 2021, 26, 101127. [Google Scholar] [CrossRef]
  13. Behar, O.; Sbarbaro, D.; Morán, L. A Practical Methodology for the Design and Cost Estimation of Solar Tower Power Plants. Sustainability 2020, 12, 8708. [Google Scholar] [CrossRef]
  14. Albarbar, A.; Arar, A. Performance Assessment and Improvement of Central Receivers Used for Solar Thermal Plants. Energies 2019, 12, 3079. [Google Scholar] [CrossRef]
  15. Khosravi, A.; Malekan, M.; Pabon, J.J.G.; Zhao, X.; Assad, M.E.H. Design parameter modelling of solar power tower system using adaptive neuro-fuzzy inference system optimized with a combination of genetic algorithm and teaching learning-based optimization algorithm. J. Clean. Prod. 2020, 244, 118904. [Google Scholar] [CrossRef]
  16. Ahmed, O.K.; Hussein, A.S. New design of solar chimney (case study). Case Stud. Therm. Eng. 2018, 11, 105–112. [Google Scholar] [CrossRef]
  17. Arzpeyma, M.; Mekhilef, S.; Kazi, S.N.; Horan, B.; Seyedmahmoudian, M.; Akram, N.; Stojcevski, A. Solar chimney power plant and its correlation with ambient wind effect. J. Therm. Anal. Calorim. 2019, 141, 649–668. [Google Scholar] [CrossRef]
  18. Arzpeyma, M.; Mekhilef, S.; Kazi, S.N.; Horan, B.; Seyedmahmoudian, M.; Akram, N.; Stojcevski, A.; Alnaimat, F.; Mathew, B. Numerical simulation of the effect of chimney configuration on the performance of a solar chimney power plant. J. Therm. Anal. Calorim. 2021, 147, 2549–2563. [Google Scholar] [CrossRef]
  19. Guo, P.; Li, T.; Xu, B.; Xu, X.; Li, J. Questions and current understanding about solar chimney power plant: A review. Energy Convers. Manag. 2018, 182, 21–33. [Google Scholar] [CrossRef]
  20. Guo, P.; Wang, Y.; Li, J.; Wang, Y. Thermodynamic analysis of a solar chimney power plant system with soil heat storage. Appl. Therm. Eng. 2016, 100, 1076–1084. [Google Scholar] [CrossRef]
  21. Habibollahzade, A.; Houshfar, E.; Ashjaee, M.; Ekradi, K. Continuous power generation through a novel solar/geothermal chimney system: Technical/cost analyses and multi-objective particle swarm optimization. J. Clean. Prod. 2021, 283, 124666. [Google Scholar] [CrossRef]
  22. Tajemnicze Aqda i Nocleg w Świątyni Zoroastrian. Available online: https://kolemsietoczy.pl/tajemnicza-aqda-nocleg-w-swiatyni-zoroastrian/ (accessed on 8 July 2022).
  23. Modele Krzepnięcia Materiałów PCM w Różnych Układach Geometrycznych. Available online: https://projekt-energie.eu/pl/platforma-edukacyjna/article-modele-krzepniecia-materialow-pcm-w-roznych-ukladach-geometrycznych.html (accessed on 11 June 2023).
  24. RUBITHERM. RT64HC. Available online: https://www.rubitherm.eu/media/products/datasheets/Techdata_-RT64HC_EN_09102020.PDF (accessed on 12 July 2022).
  25. Płuciennik-Koropczuk, E.; Myszograj, S.; Mąkowski, M. Reducing CO2 Emissions from Wastewater Treatment Plants by Utilising Renewable Energy Sources? Case Study. Energies 2022, 15, 8446. [Google Scholar] [CrossRef]
  26. Myszograj, S.; Bocheński, D.; Mąkowski, M.; Płuciennik-Koropczuk, E. Biogas, Solar and Geothermal Energy—The Way to a Net-Zero Energy Wastewater Treatment Plant—A Case Study. Energies 2021, 14, 6898. [Google Scholar] [CrossRef]
  27. Kumar, D.; Biswas, N.; Manna, N.K.; Cemal, A. Impact of chimney divergence and sloped absorber on energy efficacy of a solar chimney power plant (SCPP). Ain Shams Eng. J. 2024, 15, 102390. [Google Scholar] [CrossRef]
  28. Pradhan, S.; Chakraborty, R.; Mandal, D.K.; Barman, A.; Bose, P. Design and performance analysis of solar chimney power plant (SCPP): A review. Sustain. Energy Technol. Assess. 2021, 47, 101411. [Google Scholar] [CrossRef]
  29. Das, P.; Chandramohan, V.P. A review on solar updraft tower plant technology: Thermodynamic analysis, worldwide status, recent advances, major challenges and opportunities. Sustain. Energy Technol. Assess. 2022, 52, 102091. [Google Scholar] [CrossRef]
  30. Saleh, M.J.; Atallah, F.S.; Algburi, S.; Ahmed, O.K. Enhancement methods of the performance of a solar chimney power plant: Review. Results Eng. 2023, 19, 101375. [Google Scholar] [CrossRef]
  31. Taler, J.; Duda, P. Rozwiązania Prostych i Odwrotnych Zagadnień Przewodzenia Ciepła; Wydawnictwo Naukowo-Techniczne: Warszawa, Poland, 2003. [Google Scholar]
Figure 1. Solidification model of the PCM plate.
Figure 1. Solidification model of the PCM plate.
Energies 17 05004 g001
Figure 2. Chimney test bench.
Figure 2. Chimney test bench.
Energies 17 05004 g002
Figure 3. Construction diagram for a mini solar chimney.
Figure 3. Construction diagram for a mini solar chimney.
Energies 17 05004 g003
Figure 4. BMS visualization of solar chimney operation.
Figure 4. BMS visualization of solar chimney operation.
Energies 17 05004 g004
Figure 5. Integration of a mini solar chimney with a solidifying plate.
Figure 5. Integration of a mini solar chimney with a solidifying plate.
Energies 17 05004 g005
Figure 6. Numerical results for the system defined by Equations (13) and (14) with parameters ( B i C O N = 160 , B i a = 12.5 , S t e = 0.112 ,   a ~ = 1.13 , k ~ = 1 ) for various values of B , (a) the thickness of the solid layer over time, (b) the dependence of the average temperature over time.
Figure 6. Numerical results for the system defined by Equations (13) and (14) with parameters ( B i C O N = 160 , B i a = 12.5 , S t e = 0.112 ,   a ~ = 1.13 , k ~ = 1 ) for various values of B , (a) the thickness of the solid layer over time, (b) the dependence of the average temperature over time.
Energies 17 05004 g006
Figure 7. The thickness of the solid layer over time for B = 0 .
Figure 7. The thickness of the solid layer over time for B = 0 .
Energies 17 05004 g007
Figure 8. Scheme (a) and the photo (b) of the research unit.
Figure 8. Scheme (a) and the photo (b) of the research unit.
Energies 17 05004 g008
Figure 9. Comparison of the results of experimental and theoretical studies of temperature versus time measurements.
Figure 9. Comparison of the results of experimental and theoretical studies of temperature versus time measurements.
Energies 17 05004 g009
Figure 10. Dependence of the coefficient of heat transfer on the velocity of the air.
Figure 10. Dependence of the coefficient of heat transfer on the velocity of the air.
Energies 17 05004 g010
Figure 11. Dependence of velocity u ¯ C and temperature T C of air in chimney on chimney stack height H C for a heat flow of q ˙ = 800   W / m 2 [6].
Figure 11. Dependence of velocity u ¯ C and temperature T C of air in chimney on chimney stack height H C for a heat flow of q ˙ = 800   W / m 2 [6].
Energies 17 05004 g011
Figure 12. Comparison of air velocity in the theoretical model chimney with experimental tests for a mini-chimney [6].
Figure 12. Comparison of air velocity in the theoretical model chimney with experimental tests for a mini-chimney [6].
Energies 17 05004 g012
Figure 13. The influence of the heat transfer coefficient on the thickness of the solidified layer for different values at the plate, with B i C O N = 160 .
Figure 13. The influence of the heat transfer coefficient on the thickness of the solidified layer for different values at the plate, with B i C O N = 160 .
Energies 17 05004 g013
Figure 14. The increase in the thickness of the solidified layer and the rate of solidification as a function of time for h = 20   W / m 2 K and B i C O N = 160 .
Figure 14. The increase in the thickness of the solidified layer and the rate of solidification as a function of time for h = 20   W / m 2 K and B i C O N = 160 .
Energies 17 05004 g014
Figure 15. Heat released from a layer that has solidified over time, for h = 20   W / m 2 K and B i C O N = 160 .
Figure 15. Heat released from a layer that has solidified over time, for h = 20   W / m 2 K and B i C O N = 160 .
Energies 17 05004 g015
Figure 16. Total heat flow as a function of time for h = 20   W / m 2 K and B i C O N = 160 .
Figure 16. Total heat flow as a function of time for h = 20   W / m 2 K and B i C O N = 160 .
Energies 17 05004 g016
Table 1. Summary of operating parameters of the mini solar chimney.
Table 1. Summary of operating parameters of the mini solar chimney.
Date/TimeAir Flow Velocity in the Chimney [m/s]External Wind Speed [m/s]Insolation [W/m2]Air Temperature
[°C]
Air Flow Velocity in the Chimney [m/s]
8 February 20200.725279.24.10.7
11 February 20201.17.2164.95.81.1
2 September 20222.11.0846.621.22.1
3 September 20220.60214.10.6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gortych, M.; Lipnicki, Z.; Małolepszy, T.; Grabas, P. Solar Chimney Operation Variant. Energies 2024, 17, 5004. https://doi.org/10.3390/en17195004

AMA Style

Gortych M, Lipnicki Z, Małolepszy T, Grabas P. Solar Chimney Operation Variant. Energies. 2024; 17(19):5004. https://doi.org/10.3390/en17195004

Chicago/Turabian Style

Gortych, Marta, Zygmunt Lipnicki, Tomasz Małolepszy, and Piotr Grabas. 2024. "Solar Chimney Operation Variant" Energies 17, no. 19: 5004. https://doi.org/10.3390/en17195004

APA Style

Gortych, M., Lipnicki, Z., Małolepszy, T., & Grabas, P. (2024). Solar Chimney Operation Variant. Energies, 17(19), 5004. https://doi.org/10.3390/en17195004

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

Article Metrics

Back to TopTop