Sustainable Energy Technologies and Assessments 39 (2020) 100704

Contents lists available at ScienceDirect

Sustainable Energy Technologies and Assessments

journal homepage: www.elsevier.com/locate/seta

Original article

Numerical performance modelling of solar chimney power plants: Influence T

of chimney height for a pilot plant in Manzanares, Spain
Erdem Cucea,b,⁎, Harun Sena,b, Pinar Mert Cucea,c
a
Low/Zero Carbon Energy Technologies Laboratory, Faculty of Engineering, Recep Tayyip Erdogan University, Zihni Derin Campus, 53100 Rize, Turkey
b
Department of Mechanical Engineering, Faculty of Engineering, Recep Tayyip Erdogan University, Zihni Derin Campus, 53100 Rize, Turkey
c
Department of Energy Systems Engineering, Faculty of Engineering, Recep Tayyip Erdogan University, Zihni Derin Campus, 53100 Rize, Turkey

ARTICLE INFO ABSTRACT

Keywords: 3D axisymmetric CFD model is developed for a solar chimney power plant (SCPP) in Manzanares, Spain, and
Solar chimney power plant potential impacts of chimney height (H ) on main performance parameters are comprehensively analysed. Mesh-
CFD model independent solutions are achieved, and accuracy justification is done over the previous numerical and ex-
Chimney height perimental attempts prior to parametric research. Discrete ordinate (DO) non-grey radiation model with solar
Solar intensity
ray tracing approach is adopted in the research. A very good accordance is achieved between the numerical
Optimisation
findings and in-situ data. For five different H values, temperature, pressure and velocity distributions within the
pilot plant are achieved as well as maximum air velocity, mass flow rate of air, temperature rise in collector,
dynamic pressure difference at the turbine position, overall system efficiency and potential electrical power. It is
found that maximum air velocity thus mass flow rate shows an exponential growth in H . On the contrary,
temperature rise in collector notably reduces with the increasing H . Overall system efficiency is determined to
be 0.67% whenH = 500 m. Power output (P ) linearly rises with H . The system is capable of generating 55 and
134 kW electrical power, forH = 200 and 500 m, respectively.

Introduction renewable energy application. Working principle of a typical solar

chimney can be basically explained by the functions of two main parts
Solar energy is harnessed for various purposes today. The most of the whole plant, which are collector and chimney. Collector is the
common use of solar energy in practice appears in solar thermal and part which enables the penetration of solar thermal radiation into the
solar electricity applications. Solar thermal systems basically aim at air medium between the collector and the ground. The incoming solar
capturing incoming solar radiation within a blackened absorber surface radiation is absorbed by the ground material at a rate of its absorptivity
in which there is a direct or indirect contact between the absorber and coefficient, and as a consequence of the greenhouse effect, collector air
the working fluid [1]. On the other hand, solar electricity refers to the gets warmer with time. Momentum of air increases with the improve-
direct conversion of photons into electrical power through the semi- ment of thermal energy content of the air medium, and it is accelerated
conductor devices like photovoltaic (PV) modules [2]. It is well-docu- owing to the remarkable pressure difference provided by the massive
mented in literature that both thermal performance figures and PV cell chimney. Although solar energy is a driving force for solar chimneys,
parameters are dependent on environmental conditions [3]. In other the system is capable of producing energy even in night time because of
words, climatic parameters and locational characteristics play a key the aforesaid pressure difference between the inlet and outlet sections
role in the overall performance assessment of solar energy applications. of the chimney [5].
On the contrary to this common fact, there is a clean energy technology The idea of solar chimney dates back to 1903, however, first suc-
called solar chimney which is capable of generating electricity con- cessful applications in the world appear in 1980s. solar chimney power
tinuously irrespective of change in environmental conditions owing to plants (SCPPs) are highly in the centre of interest especially over the
its unique structure [4]. last three decades owing to their characteristic ability in large-scale
Solar chimneys do not require fossil fuels for operation, hence they harnessing of solar energy. SCPPs are giant structures, however, it is
do not cause any greenhouse gas emissions as a completely eco-friendly usually possible to express that a typical system consists of three main

⁎
Corresponding author at: Department of Mechanical Engineering, Faculty of Engineering, Recep Tayyip Erdogan University, Zihni Derin Campus, 53100 Rize,
Turkey.
E-mail address: [email protected] (E. Cuce).

https://doi.org/10.1016/j.seta.2020.100704
Received 8 January 2020; Received in revised form 10 March 2020; Accepted 31 March 2020
2213-1388/ © 2020 Elsevier Ltd. All rights reserved.
E. Cuce, et al. Sustainable Energy Technologies and Assessments 39 (2020) 100704

parts, which are chimney, collector and turbine [6]. Collector is usually Manzanares pilot plant. Through the MATLAB code developed, the
constructed by transparent materials like polythene or clear glass. It is impacts of collector diameter and chimney height on main performance
fixed horizontally or at a specific slope via steel frames, and mounted at figures of the plant are analysed. On the contrary to the findings of
a certain height from the ground to be able to welcome air from outside Zhou et al. [13], they report that there is no optimum value for the
to inside through air vents [7]. The level of greenhouse effect beneath chimney height in SCPPs. When chimney height is increased to 2000 m,
the collector is directly associated with the collector radius, transmis- power output is predicted to enhance about sevenfold and exceeds
sivity coefficient of the collector material and the thermophysical 400 kW according to their results. Larbi et al. [16] propose a mathe-
properties of the ground. Materials with high specific heat capacity are matical model for a typical plant with a chimney height of 250 m and a
usually considered on the ground of SCPPs. Stored energy within the collector diameter of 600 m. For the said geometric parameters, power
ground material during the daytime contributes to maintain the output of the plant is predicted to be 3.42 MW. However, it reaches
buoyant forces even after the sunset. Chimney as the core of the plant is 7 MW when the chimney height is enhanced to 500 m. They also report
mounted in the middle of collector, and provides a notable pressure that power output reduces with the increase in ambient temperature.
difference between chimney inlet and outlet. Geometric parameters Shahi et al. [17] develop a 2D model for the Manzanares pilot plant,
such as chimney height and chimney diameter play a key role in the and examine the influence of chimney height on power output. For the
said pressure difference. Maximum air velocities are usually observed at existing geometric parameters, power output is found to be 48 kW for
a certain height above the chimney inlet, and a turbine is placed at that the solar radiation of 1000 W/m2 whereas it is reported to be 50 kW in
position to convert kinetic energy of air to electrical energy. Current the experimental studies. Then, numerical research is repeated for the
turbines utilised in SCPPs can succeed this energy conversion by an chimney heights of 400, 600, 800 and 1000 m. It is achieved through
efficiency over 80% as reported in literature [8]. the findings that the power output exceeds 300 kW for the chimney
The first pilot application of SCPPs appears in 1982 in Manzanares, height of 1000 m. Dhahri et al. [18] evaluate the performance figures of
Spain. The pilot plant has a chimney height of 194.6 m, chimney dia- SCPPs with respect to changes in environmental and geometric para-
meter of 10.16 m, collector diameter of 244 m and collector height of meters. When chimney height is risen from 200 to 300 m, mass flow
1.85 m. A fixed horizontal collector design is preferred for the system rate of the system is found to improve by about 20%. On the other
covering 40,000 m2 membrane and 6000 m2 clear glass. In the early hand, this change reduces the temperature rise of air along the collector
1980s, Haaf et al. [9] conduct several experimental and numerical re- from 22 to 16 °C. Toghraie et al. [19] conduct a 3D CFD research for a
search on the pilot plant for an overall performance evaluation of the SCPP with a chimney height of 100 m, a chimney diameter of 4 m, a
entire system. Their results indicate that there is a temperature rise of collector diameter of 200 m and a collector height of 2 m. Standard k-ε
20 °C in the collector air from inlet to outlet. For the solar intensity of turbulence model is considered in the numerical methodology. For the
1000 W/m2, maximum velocity of air is reported to be about 15 m/s solar intensity of 800 W/m2, power output is determined to enhance
around the turbine inlet. The power output of the pilot plant corre- more than 100% when the chimney height is increased twofold. They
sponding to this velocity is given to be 55 kW. It is also underlined that also report that power output is expected to decrease after a certain
chimney height is a significant geometric parameter on the perfor- value of chimney height.
mance figures of a SCPP. It is concluded from the literature survey that SCPPs are in-
Chimney height is widely studied in other works as well over the vestigated by many researchers through theoretical and numerical
last three decades through theoretical and numerical analyses. For in- models, and most of the works merely focus on main performance
stance, Mullett [10] carries out a research for overall performance as- figures such as power output. It is addressed in different studies that
sessment of SCPPs with different chimney heights. It is found that the chimney height is a significant geometric parameter in SCPPs, however
efficiency of a plant with a chimney height of 1000 m is about 1%, there are still contradictions between the results of previous researchers
which is also verified in similar works. The results also reveal that the concerning with the influence of chimney height on power output. In
estimated power output of a typical plant with a chimney height of order to elucidate the abovementioned contradictions and incon-
200 m and chimney diameter of 250 m is 100 kW whereas it is sistencies, an elaborative CFD research is presented within the scope of
10000 kW with a chimney height and diameter of 400 and 1780 m, this research. A unique 3D axisymmetric CFD model is developed with a
respectively. Ngala et al. [11] numerically investigate the dependency DO oriented solar ray tracing radiation model. Following the accuracy
of power output on chimney height. It is achieved from their results that justification of the numerical results through the performance reports of
there is a linear relationship between the chimney height and the power Manzanares pilot plant, impact of chimney height is comprehensively
output in SCPPs. Choi et al. [12] develop a mathematical model to assessed. Temperature, velocity and pressure distributions as well as
evaluate how the geometric parameters and heat storage medium on mass flow rates and power output values within the plant are illustrated
the ground affect the power output of a typical SCPP. The results in- as a function of chimney height for a reliable and scientific metho-
dicate that the power output is enhanced by the rise in chimney height dology.
and diameter. On the other hand, considering heat storage medium on
the ground is found to reduce the power output in the daytime, whereas Theoretical background
it notably contributes to energy generation after sunset. Zhou et al. [13]
develop a numerical model for determining the optimum chimney The present study aims at proposing a reliable and accurate simu-
height for the pilot plant in Manzanares. They report that the power lation model for SCPPs by taking geometric and material characteristics
output rises with the chimney height until a specific value of 615 m. For of Manzanares pilot plant into consideration. Specifically, chimney
greater chimney heights, power output is observed to decrease for the height is investigated in the numerical model as an independent vari-
existing geometric parameters of the plant. They also address that for able for its potential influence on velocity, temperature and pressure
optimum operational conditions, chimney height needs to be increased distributions as well as power output. Accuracy justification is done for
in proportion to the collector diameter. Another model is developed by maximum air velocities for various solar intensity levels in the range of
Ming et al. [14] for performance investigation of SCPPs. The findings of 100–1000 W/m2. Continuity, momentum, energy and turbulence model
the 3D model show that there is a linear relationship between the equations are solved simultaneously to achieve the aforesaid perfor-
chimney height and power output. System efficiency is also found im- mance parameters. Following assumptions are considered in the present
proving with the rise in chimney height. For a solar intensity of work: Numerical analyses are conducted for the steady-state conditions.
1000 W/m2, enhancing chimney height from 200 to 500 m yields to a Changes in climatic conditions are assumed to be negligible. Boussinesq
threefold improvement in the power output. approximation is preferred to calculate the change of air density. Heat
Li et al. [15] theoretically investigate the performance of losses across the chimney walls can be ignored. There is a uniform and

2
E. Cuce, et al. Sustainable Energy Technologies and Assessments 39 (2020) 100704

constant solar radiation during each process. The radiative transfer where Mt is the Mach number in turbulent flow. Mt can be calcu-
equation is represented by the discrete ordinate (DO) non-grey radia- lated by
tion model. Following continuity, momentum and energy equations are
k
solved coupled in the numerical analyses: Mt =
a2 (11)
▪ Continuity equation: The inverse effective Prandtl numbers for and k are represented by
( u) ( v) ( w) and k , respectively. The additional term which differs RNG k-ε
+ + + =0 turbulence model from the standard model is given below:
t x y z (1)
Cµ 3 (1 / 0) 2
R =
▪ Momentum equation: 1+ 3 k (12)
( u) ( uu) ( uv ) ( uw ) p 2u 2u 2u
In the aforesaid equation, 0 = 4.38, =0.012 and η = Sk/ [20].
+ + + = +µ + +
t x y z x x2 y2 z2 Instant streamline and strain curvatures can be adapted in the RNG k-ε
(2) turbulence model. Boussinesq model helps to calculate the change in air
density in the numerical calculations. By considering the buoyancy
( v) ( vu) ( vv ) ( vw ) p 2v 2v 2v
parameter in the momentum equation, air density is defined in the
+ + + = +µ + 2 + 2
t x y z y x2 y z model as follows:
(3) ( a)g a (T Ta ) (13)
( w) ( wu) ( wv ) ( ww ) In the said equation, a is the air density at ambient temperature of
+ + +
t x y z Ta , and refers to thermal expansion coefficient. When the literature on
p 2w 2w 2w SCPPs is evaluated, it is easy to say that there are various methods to
= +µ + + + g (T Ta)
z x2 y2 z2 (4) determine the power output of the system. In the meantime, main ef-
fective parameters are listed as turbine pressure drop, chimney effi-
▪ Energy equation: ciency, turbine efficiency and collector efficiency in most cases.
Different figures are utilised in literature for the turbine efficiency ( t ).
However, t is taken to be 0.8 in the present work [21]. Collector ef-
( cT ) (pcuT ) (pcvT ) (pcwT ) 2T 2T 2T
+ + + = + +
t x y z x2 y2 z2 ficiency ( coll ) is another significant term of the numerical model. Its
(5) value is highly dependent on the incoming solar radiation, and calcu-
lated as follows:
It is well-documented in literature that the air flow in SCPPs belongs
to turbulent flow, and a wide range of turbulence models is utilised by Q
=
researchers in the numerical modelling works of these systems. When coll
Acoll G (14)
the whole turbulence models are briefly evaluated, it can be easily as-
In the said equation, Q refers to available thermal power of the
serted that RNG k-ε and standard k-ε turbulence models are much more
system, and given by
in the centre of interest as a consequence of their efficacy and simpli-
city. Further evaluations reveal that RNG k-ε turbulence model leads in Q = mcp (T1 T2) = mcp T (15)
terms of a thorough accordance and coherence owing to the turbulent
where temperature increase through the collector is expressed with
viscosity term that the model has.
T . For the determination of turbine pressure drop, there are also dif-
ferent expression in literature. Turbine pressure drop ( Pt ) is associated
▪ RNG k-ε turbulence model equation:
with the air velocity at the turbine inlet through a linear regression as
k given below:
( k) + ( kui) = k µeff + Gk + Gb + YM + Sk
t xi xj xj (6) Pt = 18, 87v 57, 59 (16)
Power output in a SCPP can be calculated through the following
( )+ ( ui )
t xi equation:
2 T
= µeff + C1 (Gk + C3 Gb) C2 Re + S Ptot = mgHc
xj xj k k (7) Ta (17)

The presence of mean velocity gradient yields to the turbulent ki- where Ta is the ambient temperature, T is the temperature difference
netic energy production. This is represented in equations (6) and (7) between inlet and outlet sections of the collector, Hc is the chimney
with Gk , and the said term is expressed as follows: height, g is the gravitational acceleration and m is the mass flow rate of
air. By considering Qv is the volumetric flow rate, power output can also
ui be obtained by
Gk = ui uj
uj (8)
Ptot = t Pt Qv (18)
On the other hand, buoyant forces also contribute to turbulent ki-
netic energy production, which is given with the term Gb in the same
CFD analyses
equations. Gb can be calculated by the following equation:
µt t The numerical methodology adopted in the present research is
Gb = gi
Prt x i (9) based on a unique 3D axisymmetric CFD model developed for the SCPP
In the turbulence model equation, YM refers to how the fluctuating in Manzanares, Spain. Exact material properties and geometric para-
dilatation in compressible turbulence affects the overall dissipation meters of the pilot plant listed in Table 1 are utilised in the CFD model
rate. YM can be expressed by the following equation: for a realistic, reliable and scientific approach. It would be useful to
illustrate the outline of the solar chimney concept analysed in the CFD
YM = 2 Mt 2 (10) research with the system components and boundary conditions as

3
E. Cuce, et al. Sustainable Energy Technologies and Assessments 39 (2020) 100704

Table 1 Table 2
Geometric characteristics of the SCPP in Manzanares, Glass, ground and chimney characteristics of the plant utilised in CFD model.
Spain.
Physical property (unit) Glass Ground Chimney
Geometric parameter Value
3
Density (kg/m ) 2500 2160 2719
Mean collector radius 122.0 m Thermal conductivity (W/mK) 1.15 1.83 202.4
Mean collector height 1.85 m Specific heat (J/kgK) 750 710 871
Chimney height 194.6 m Absorption coefficient 0.03 0.9 0
Chimney radius 5.08 m Transmissivity 0.9 Opaque Opaque
Ground thickness 0.5 m Emissivity 0.1 0.9 1
Refractive index 1.526 1 1
Thickness (m) 0.004 0.5 0.00125

depicted in Fig. 1. Physical properties of the ground, chimney and

collector glass are of vital importance for the efficacy of the numerical
Table 3
approach, and thus the aforesaid data is obtained from the in-situ
Main parameters considered in CFD modelling.
measurements as shown in Table 2. Besides, there are some construc-
tional and operational parameters, climatic data and thermophysical Solar radiation (W/m2) 1000
properties shown in Table 3, which are utilised in the CFD analyses. As
Ambient pressure (Pa) 101,325
a consequence of intensive cell counts, iterations take a significant Ambient temperature (K) 293.15
period of time ranging from 24 to 36 h for each case. Since the case is Gravitational acceleration (m/s2) 9.81
axisymmetric, one-quarter part of the pilot plant is modelled in the Ambient air density (kg/m3) 1.2046
Ideal gas constant (J/kgK) 287
analyses via two plane symmetries (YZ and XZ planes) to expedite the
Air conductivity (W/mK) 0.0259
convergence. When the whole geometry is considered, it can be as- Kinematic viscosity of air (m2/s) 1.48 × 10−5
serted that curvatures outweigh the smooth surfaces, hence triangular Air heat capacity (J/kgK) 1006.43
mesh type is preferred in the CFD analyses as illustrated in Fig. 2. Stefan-Boltzmann constant (W/m2K4) 5.667 × 10−8
ANSYS FLUENT is utilised in the CFD analyses as one of the most re- Turbine pressure drop ratio 2/3

liable and well-documented commercial software in market. Finite

volume approach, which is used by ANSYS FLUENT is an ideal method
Results and discussion
for solving the governing equations in the present work. In order to
solve the momentum equation in the numerical research, RNG k-ε
Mesh-independent solution is comprehensively proposed prior to
turbulence model is considered. To be able to determine the interaction
the demonstration of the results of parametric research. In this respect,
between the velocity and pressure figures of air, SIMPLE algorithm is
over the three main performance parameters, which are mass flow rate
used. PRESTO technique is adopted for the pressure interpolation.
of air (m ) at the turbine inlet, maximum air temperature (Tm ) and
Governing equations are discretised by second order upwind method.
maximum air velocity (Vm ), dependency on cell counts is determined as
DO non-grey radiation approach with solar ray tracing method is pre-
shown in Table 4. For the cell count of 386041, % change in m and Vm is
ferred for radiation equation. Exact latitude and longitude data of
observed to be 0.24 and 1.75, respectively. As the aforesaid figures are
Manzanares, Spain is utilised in the radiation model for accurate cal-
satisfactory enough in comparison with the similar works in literature,
culation of solar beam direction. Boussinesq approximation is adopted
the whole CFD is research is carried out over this cell count. Following
to determine the variation of air density with temperature. An ideal
the satisfaction of mesh-independent solution, accuracy verification of
value of 10-6 is selected for the convergence criterion in the iterations.
the CFD model is done through the Vm values in the pilot plant acquired
Parametric research is conducted once the mesh-independent solutions
for various solar intensities in the range of 100–1000 W/m2. It can be
are obtained.
easily understood from the preliminary findings of the CFD model
shown in Fig. 3 that present work is in very good accordance with the
experimental results [22], and provides the best fitted tendency of Vm

Fig. 1. Structural details and boundary conditions considered in the CFD research.

4
E. Cuce, et al. Sustainable Energy Technologies and Assessments 39 (2020) 100704

Fig. 2. Mesh model adopted in the CFD analyses of Manzanares pilot plant.

Table 4 with respect to changes in environmental conditions. When pioneer

Mesh-independent solution results for the CFD analyses. numerical works in the scope [14] are analysed in detail, it is achieved
Cell count m Tm Vm % change in % change that notable deviations occur especially for high solar intensity levels.
Vm inm This can be ascribed to the considerable role of turbulence and radia-
tion models in the overall efficacy after a certain value of solar radia-
259,800 278.298 kg/s 318.250 K 14.420 m/s – – tion. On the contrary to the said works, present work provides an ex-
305,619 278.822 kg/s 317.864 K 13.996 m/s 2.94 0.18
386,041 279.513 kg/s 317.852 K 14.242 m/s 1.75 0.24
cellent agreement with the measurements owing to RNG k-ε turbulence
model and DO non-grey radiation model assisted by solar ray tracing
approach.Vm is found to be 14.24 m/s for the solar radiation level of

Fig. 3. Accuracy verification of the CFD model over the maximum velocity values.

5
E. Cuce, et al. Sustainable Energy Technologies and Assessments 39 (2020) 100704

beneath the collector is remarkably deteriorated by the increase in

chimney height. This can be attributed to the shorter period of time of
air while it is flowing through the collector area when the buoyant
forces are improved with H . In the CFD analyses, air is assumed to enter
the collector with a temperature of 293.15 K. At the outlet of the col-
lector, air is found to reach a maximum temperature of 304.4 K
whenH = 100 m. On the contrary, it can only rise up to 301.2 K when
the chimney height is 500 m. From this point of view, it can be con-
cluded that chimney height and collector diameter are the performance
parameters which interact each other, and for an optimum perfor-
mance, they are required to be enlarged by an appropriate ratio be-
tween the said parameters. In SCPPs there is a connection path between
the outlet of collector and the inlet of chimney. In this area, some part
of the kinetic energy of air transforms into thermal energy due to in-
tensive collisions of air molecules, and by doing so, additional increase
in air temperature is observed as shown in Fig. 9. This is expressed as
the overall temperature rise in the collector, and should not be confused
with maximum temperature increase. Following equation is developed
through the regression model between average T and H :
Fig. 4. Maximum velocity (Vm ) in the pilot plant as a function of chimney height T = 12.63 + 11.94e 0.005963H (21)
(H ).
Mass flow rate (m ) is an important parameter which needs to be
2 analysed in terms of its potential change with chimney height. Fig. 10
1000 W/m , while it is given to be 15 m/s in the experimental reports.
illustrates that m exponentially increases with H . For H = 100 m, m is
As previously mentioned, Manzanares pilot plant has a chimney
found to be 886.76 kg/s whereas it is 1496.35 kg/s for the chimney
height (H ) of 194.6 m. By taking all the other geometric parameters
height of 500 m. This result is also an expected output of improved
constant, H is considered to be as an independent variable within the
buoyant forces in the chimney with the enhanced momentum effects. A
scope of this research, and main performance figures are obtained for
similar tendency is already achieved for the Vm values as previously
five different values of H in the range of 100–500 m. Vm is the first
discussed. Regression model proposes the following correlation be-
parameter in this respect, and Fig. 4 demonstrates its variation with H .
tween m and H :
It is clear from the data that Vm exponentially rises with H . This can be
easily attributed to the improved buoyant forces with the increase in m = 686.8(1 e 0.02597H ) 1142(1 e 0.002487 )
(22)
chimney height as a consequence of enhanced pressure difference be-
Efficiency (η) of SCPPs is notably poor because of the huge collector
tween inlet and outlet sections of chimney. A regression analysis is
area. In the efficiency expression, the incoming solar power is written to
performed to the numerical data, and the following equation is devel-
the denominator which yields to very low figures. In literature, overall
oped for the first time in literature for the relationship between Vm and
efficiencies of SCPPs vary in the range of 0.2–1% depending on geo-
H:
metric parameters [19]. An efficiency analysis is also carried out within
Vm = 8.178(1 e 0.03551H ) + 14.25(1 e 0.00268H )
(19) the scope of this research to examine the influence of chimney height.
Fig. 11 indicates that there is a linear relationship between η and H .
The proposed regression model indicates that there is an expected This is also justified through the regression analysis as follows:
optimum figure for chimney height since Vm has a tendency of ex-
= 0.003458 + 0.001307H (23)
ponential rise to a maximum data. However, it is difficult to make a
generalisation regarding the exact optimum value of H since a com- Power output (P ) of the pilot plant is also analysed as a function of
prehensive life-cycle cost analysis is required for such an attempt cov- chimney height as depicted in Fig. 12. It is achieved that P rises linearly
ering the chimney cost. Detailed static contours of air velocities in the with H . For a chimney height of 200 m, the system is found to generate
pilot plant for the chimney height of 100 and 500 m are illustrated in 55 kW electrical power while it is improved to 134 kW when H is in-
Fig. 5. It is unequivocal from the findings that buoyant effects re- creased to 500 m. The relationship between the power output and
markably contribute to the momentum which yields to improvements chimney height can be expressed by the following equation:
in Vm values around the turbine inlet. P = 712.4 + 270.8H (24)
Pressure distributions through the power plant especially about the
turbine inlet are of vital importance for performance assessment of
SCPPs. In this respect, dynamic pressure differences ( Pdyn ) at 9 m Conclusions
above the chimney inlet in which the turbine is fixed are determined as
depicted in Fig. 6. It is clear from the findings that Pdyn values are Chimney height is comprehensively analysed in terms of its poten-
notably enhanced by the rise in H . Maximum Pdyn forH = 100 m is tial impacts on main performance parameters. 3D axisymmetric CFD
found to be −67.22 Pa whereas it is −147.24 and −202.82 Pa model is developed for the geometric parameters and material char-
forH = 300 and 500 m, respectively as illustrated in Fig. 7. Following acteristics of Manzanares pilot plant. Following the accuracy justifica-
regression model is obtained through the parametric research: tion of the CFD results through previous works in the scope, parametric
research is conducted for five different chimney heights in the range of
Pdyn = 328.4(1 e 0.001632H ) 19.33(1 e 0.02488)
(20) 100–500 m. Temperature, velocity and pressure contours are achieved
for each case. Mass flow rate of air, average temperature rise in col-
Exponential rise of momentum effects with the chimney height can lector, average dynamic pressure difference at the turbine position,
also be verified through the aforesaid data. Temperature rise across the system efficiency and power output are numerically determined.
collector of the pilot plant is also investigated in the present work as Following bullet points can be concluded from the research:
shown in Fig. 8. Chimney height is a significant parameter in terms of
momentum effects as previously mentioned. However, it is observed ▪ CFD model shows a very good agreement with the previous
from the temperature contours that thermal energy content of air

6
E. Cuce, et al. Sustainable Energy Technologies and Assessments 39 (2020) 100704

Fig. 5. Air velocity contours within the pilot plant for the chimney height of 100 and 500 m.

measurements in the pilot plant for certain radiation conditions. Vm

is reported to be about 15 m/s in the experimental analyses for a
solar intensity of 1000 W/m2, whereas it is found 14.24 m/s in the
CFD research.
▪ Vm thus m shows an exponential growth with the chimney height.
▪ Pressure distributions within the plant are of vital importance in
terms of accurate assessment of buoyant forces and momentum ef-
fects. Maximum Pdyn for a chimney height of 100 m is found to be
−67.22 Pa while it is −147.24 and −202.82 Pa forH = 300 and
500 m, respectively.
▪ Temperature rise in collector is negatively affected by the rise in
chimney height. Air reaches a maximum temperature of 304.4 K
whenH = 100 m. On the other hand, it can only rise up to 301.2 K
when the chimney height is 500 m.
▪ m can be enhanced from 886.76 to 1496.35 kg/s when H is in-
creased from 100 to 500 m.
▪ Efficiency of the plant shows a linearly increasing tendency with
chimney height. η is calculated to be 0.15% forH = 100 m while it
is improved to 0.67% whenH = 500 m.
▪ A linear relationship between power output of the plant and Fig. 7. Dynamic pressure difference ( Pdyn ) around the turbine as a function of
chimney height is observed. For H = 200 and 500 m, the system is chimney height (H ).
capable of providing 55 and 134 kW electrical power, respectively.
Pinar Mert Cuce: Conceptualization, Validation, Visualization.
CRediT authorship contribution statement
Declaration of Competing Interest
Erdem Cuce: Writing - original draft, Writing - review & editing,
Supervision. Harun Sen: Data curation, Formal analysis, Software. The authors declare that they have no known competing financial

Fig. 6. Pressure distributions around the turbine inlet of the pilot plant for different chimney heights.

7
E. Cuce, et al. Sustainable Energy Technologies and Assessments 39 (2020) 100704

Fig. 8. Temperature distributions through the collector of the pilot plant for different chimney heights.

Fig. 9. Average temperature rise in collector ( T ) as a function of chimney

height (H ).
Fig. 11. Efficiency of the pilot plant ( ) as a function of chimney height (H ).

Fig. 12. Power output of the pilot plant (P ) as a function of chimney height
(H ).
Fig. 10. Mass flow rate of air (m ) as a function of chimney height (H ).

8
E. Cuce, et al. Sustainable Energy Technologies and Assessments 39 (2020) 100704

interests or personal relationships that could have appeared to influ- [11] Ngala GM, Oumarou MB, Huhammad AB, Shuwa M. Evaluation of solar chimney
ence the work reported in this paper. power plant in semi-arid region of Nigeria. Arid Zone J Eng, Technol Environ
2015;11:1–12.
[12] Choi YJ, Kam DH, Park YW, Jeong YH. Development of analytical model for solar
References chimney power plant with and without water storage system. Energy
2016;112:200–7.
[13] Zhou X, Yang J, Xiao B, Hou G, Xing F. Analysis of chimney height for solar chimney
[1] Cuce E, Cuce PM, Saxena A, Guclu T, Besir AB. Performance analysis of a novel solar
power plant. Appl Therm Eng 2009;29(1):178–85.
desalination system – Part 1: The unit with sensible energy storage and booster
[14] Ming TZ, Wei L, Guoliang X. Analytical and numerical investigation of the solar
reflector without thermal insulation and cooling system. Sustainable Energy
chimney power plant systems. Int J Energy Res 2006;30(11):861–73.
Technol Assess 2020;37:100566.
[15] Li JY, Guo PH, Wang Y. Effects of collector radius and chimney height on power
[2] Cuce E, Cuce PM. Optimised performance of a thermally resistive PV glazing
output of a solar chimney power plant with turbines. Renewable Energy
technology: An experimental validation. Energy Rep 2019;5:1185–1195.
2012;47:21–8.
[3] Cuce E, Bali T. Variation of cell parameters of a p-Si PV cell with different solar
[16] Larbi S, Bouhdjar A, Chergui T. Performance analysis of a solar chimney power
irradiances and cell temperatures in humid climates. Sharjah: United Arab Emirates;
plant in the southwestern region of Algeria. Renew Sustain Energy Rev
2009.
2010;14(1):470–7.
[4] Cuce E, Cuce PM. Performance assessment of solar chimneys: Part 1 – Impact of
[17] Shahi DVV, Gupta A, Nayak VS. CFD Analysis of solar chimney wind power plant by
chimney height on power output. Energy Res J 2019;10:11–9.
Ansys Fluent. Int J Technol Res Eng 2018;5(9):3746–51.
[5] Schlaich J, Bergermann R, Schier W, Weinrebe G. Design of commercial solar up-
[18] Dhahri A, Omri A, Orfi J. Numerical study of a solar chimney power plant. Res J
draft tower systems utilization of solar induced convective flows for power gen-
Appl Sci Eng Technol 2014;8(18):1953–65.
eration. J Sol Energy Eng 2005;127(1):117–24.
[19] Toghraie D, Karami A, Afrand M, Karimipour A. Effects of geometric parameters on
[6] Cuce E, Cuce PM. Performance assessment of solar chimneys: Part 2 – Impacts of
the performance of solar chimney power plants. Energy 2018;162:1052–61.
slenderness value and collector slope on power output. Energy Res J 2019;10:20–6.
[20] Abdelmohimen MA, Algarni SA. Numerical investigation of solar chimney power
[7] Fasel HF, Meng F, Shams E, Gross A. CFD analysis for solar chimney power plants.
plants performance for Saudi Arabia weather conditions. Sustain Cit Soc
Sol Energy 2013;98:12–22.
2018;38:1–8.
[8] Ming TZ, Zheng Y, Liu C, Liu W, Pan Y. Simple analysis on thermal performance of
[21] Nizetic S, Ninic N, Klarin B. Analysis and feasibility of implementing solar chimney
solar chimney power generation systems. J Energy Inst 2010;83(1):6–11.
power plants in the Mediterranean region. Energy 2008;33(11):1680–90.
[9] Haaf W. Solar chimneys: part II: preliminary test results from the Manzanares pilot
[22] Rabehi R, Chaker A, Aouachria Z, Tingzhen M. CFD analysis on the performance of
plant. Int J Sustain Energ 1984;2(2):141–61.
a solar chimney power plant system: Case study in Algeria. Int J Green Energy
[10] Mullett LB. The solar chimney Overall efficiency, design and performance. Int J
2017;14(12):971–82.
Ambient Energy 1987;8(1):35–40.