**Open Access Library Journal**

Vol.03 No.05(2016), Article ID:69330,18 pages

10.4236/oalib.1102666

Enclosed Solar Chimney Power Plants with Thermal Storage

Christos D. Papageorgiou

Electrical Engineering Department, National Technical University of Athens, Athens, Greece

Copyright © 2016 by author and OALib.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 19 April 2016; accepted 23 May 2016; published 27 May 2016

ABSTRACT

The main problem associated with intermittent operation of the wind and PV power plants is the lack of low cost electrical storage methods. In the present paper a new type of solar chimney technology that we shall heretofore call “Enclosed Solar Chimney Power Plant technology (ESCP technology) stands for an effective alternative of combining both low cost electricity generation as well as smooth, uninterrupted operation for 24 hours/day and for all 365 days/year. The ESCP technology also has the following benefits in relation to PV technologies: 1) Local manufacturing and use of common material for their construction (steel, glass, concrete and conventional axial fans as air turbines); 2) Easy and low cost maintenance; 3) Longer operating life up to 50 years (its solar chimney up to 100 years); 4) No water consumption; 5) Zero effect on living nature and ambient; 6) Their empty protected surface area beneath the glassed roofs of their greenhouses can be used for desalination with parabolic trough mirrors.

**Keywords:**

Solar Chimney, Enclosed Greenhouse

**Subject Areas:** **Electric Engineering**

1. Introduction

One of the major problems in the wide application of both solar and wind energy power plants is their intermittent operation. Although their levelized cost per kWh has already become comparable to the cost of kWh generated by fossil fuels (coal, gas, nuclear), the projected participation up to the year 2035 in overall energy demand is unsatisfactory. The main problem associated with intermittent operation is the lack of sufficiently mature electrical storage methods. According to a recent article by Bill Gates titled “We Need an Energy Miracle” http://www.theatlantic.com/magazine/archive/2015/11/we-need-an-energy-miracle/407881/, under reasonable estimations the cost of kWh for storage could be between 30 and 100 cents, an amount much higher than the levelized KWh electricity generation cost, which do not exceed 10 cents.

In the present paper we intend to show that a new type of solar chimney technology that we shall heretofore call “Enclosed Solar Chimney Power Plant” (ESCP) technology stands for an effective alternative capablity of combining both low generation cost as well as smooth, uninterrupted operation for 24 hours/day and for all 365 days/year. This is achieved with the introduction of a low cost artificial water storage facility which can always be constructed on ground near the base of the solar chimney structure and protected by external winds.

2. The Solar Chimney Power Plant as an Artificial Wind Generator

Conventional solar chimney power plants (Solar Updraft Power Plants) are made of the following cooperative systems:

・ The solar air up-drafting system that is made of two parts:

A circular solar collector, which is a circular greenhouse with a transparent roof, that is open to its perimeter. The solar irradiation through the transparent roof of the greenhouse warms the ground beneath the greenhouse, while its transparent roof blocks the thermal radiation of the ground and decreases the thermal losses due to the air convection. Thus the air temperature inside the greenhouse is increasing and this warm air becomes lighter than the ambient air.

Additionally, a cylindrical solar chimney is placed in the center of the circular solar collector, through which the warm air is forced due to its buoyancy to escape towards the upper atmosphere. The escaping warm air is constantly replaced by fresh air entering through the open periphery of the greenhouse. So the warm air stream moves constantly, so long as the greenhouse ground gets warmer than the ambient air and this phenomenon can go on for several hours after the sunset.

・ The electricity generating power system is made of:

One or more air turbines placed inside the solar chimney (of vertical axis) or around the base of the solar chimney (of horizontal axis), which are forced to rotate by the moving stream of warm air escaping through the solar chimney. The air turbines are engaged, via appropriate gear boxes, to electric generators which are also forced to rotate, generating electricity.

Extensive description and important information about solar chimney technology is posted on the website https://en.wikipedia.org/wiki/Solar_updraft_tower.

As a result we can say that the solar chimney power plants practically are systems creating artificial wind. The stream of artificial wind entering through the open periphery of their greenhouses and escaping through their solar chimneys, is forcing the air turbines and their engaged electric generators to rotate, generating electricity.

An indicative figure of a conventional solar chimney power plant with one turbine inside the solar chimney is shown in Figure 1.

3. The Enclosed Solar Chimney Power Plant (ESCP)

The enclosed solar chimney power plant (ESCP) is a simpler and lower cost solar chimney power plant, where the name “enclosed” stems from its solar collector (greenhouse) being encircled by a peripheral wall. The air turbines of the ESCPs, without gear boxes, are placed in proper openings on the peripheral wall that is enclosing the solar collector (greenhouse). An indicative model of an enclosed solar chimney power plant (ESCP) is shown in Figure 2(a) & Figure 2(b).

The overall cross section area of the air turbines is several times smaller than the solar chimney cross section area. Thus the air speed entering into the air turbines is also several higher than the speed of the up-drafting warm air inside the solar chimney. Thus the rotational speed of the air turbines is increasing consequently.

A careful choice of the air turbine diameters in relation to solar chimney diameter for optimum efficiency can define their rotational speed to be equal to the generators rotational speed (that depends upon the grid frequency and the number of their pole pairs).

Thus the gear boxes can be omitted and the ESCPs electricity generation system has a lower cost. Furthermore ordinary axial air fans, existing in abundance in the market, can be used as reciprocal machines i.e. as air turbines with ESCPs. The existing axial air fans have been tested for many years in many applications and their efficiencies are not less than 70%.

However the air speed of the rising air inside the solar chimney is varying during the 24 hour operation and day by day, due to seasonal solar irradiation variation. Also external winds can alter the chimney air speed and

Figure 1. A schematic representation of a conventional solar chimney power plant of one air turbine placed inside its solar chimney.

(a) (b)

Figure 2. (a) An ESCP made of a octagonal greenhouse enclosed by a peripheral wall and sixteen air turbines placed in proper openings on its peripheral wall. (b) A detail of the previous ESCP greenhouse with one air turbine and its opening’s air stop.

the turbines air speeds. Thus in order to match continuously the turbine air speed to chimney air speed for an optimal operation and maximum efficiency and electricity generation we should control the electro-mechanical air stop systems.

The openings, where the air turbines are placed, are equipped with electro-mechanical air stop systems. Every air stop system can be activated separately, in order to close firmly or open the respective opening and release or stop the incoming air stream that is forcing the respective air turbine to rotate.

Thus when the speed of the air stream passing through the air turbine is lower than a limit, it is advisable to close up this opening activating its air stop system and terminating the operation. By this action, that can be controlled by a microprocessor expert system, the air speed through the remaining openings will increase and consequently causing a subsequent increase of their electricity generating units power. Additionally, at the moment the chimney air speed starts increasing one or more air stop systems automatically disengage opening the operation of their respective air turbines. Finally the overall electrical power of the ESCP can be increased in a controllable way.

Increasing the operational speed of the air turbines of the ESCPs in comparison to the air speeds forcing the air turbines of a respective conventional solar chimney power plant, can multiply the annual electricity generation several times as will be shown in the next paragraph.

The ground area of the closed solar collector (greenhouse), beneath its glass roof, should be covered with a thin black plastic or elastic carpet, in order to increase the solar irradiation absorbance of the solar collector. Beneath this carpet, in a part near the center, a concrete structure filled with water can be built, in order to increase the thermal storage capacity of the ESCP. Thus the ESCP can have a smooth and continuous operation 24 h/day 365 days/year.

The wall at the perimeter of the ESCP protects the thermal storage facility by heavy thermal losses due to external winds, especially at night, increasing its electricity generation. Thus ESCP is a solar technology with smooth and non intermittent operation that can be used in proper cooperation with the rest of renewable technologies in order to achieve timely the extinction of fossil fuels from electricity generation.

Thus the main benefits of ESCPs in comparison to standard conventional solar chimney power plants can be summarized as follows:

- Avoidance of gear boxes;

- Use of existing off-the-shelve air turbines (used as axial air fans for instance);

- Ability to optimize ESCP’s power output via microprocessor controlled air openings;

- Possibility of enhancing power output via openings of overall smaller cross section area in comparison to chimney cross-section area;

- Wind protected thermal storage facility inside their enclosed greenhouse.

4. Operational Equations for the ESCP

Let us define the overall opening cross section area of the operating air turbines of the ESCP to be while the chimney cross section area is where d is chimney internal diameter. If the air speed passing through the ducted air turbines is and the chimney air speed is we may define the operational ratio

(1)

The second equality is arising by the fact that the air mass flow is the same in air turbines and the solar chimney and considering that the difference in the relative air densities is negligible.

We define the maximum air speed inside the solar chimney that can be achieved with an open greenhouse of a conventional solar chimney power plant as. We may then introduce a simple saturation relation for as a function of given by the equation

(2)

In (2) we introduce the parameter a as an approximate constant depending mainly on the geometric characteristics of the ESCP, as openings overall cross-section, chimney height and diameter, greenhouse diameter etc. The parameter is estimated to be in the range.

The electric power generated by the air turbines will be given by

(3)

An equivalent conventional solar chimney power plant with one air turbine inside the chimney, would generate a standard turbine electrical power given by

(4)

Taking the ratio of (3) and (4) will be given as

(5)

We define as the Geometric Performance Index of the ESCP. In Figure 3 and Figure 4 the

functions f(r) and GP are shown for various a values.

As we can see by the figures GP can become several times higher than 1. For example for a reasonable estimation of a around 5 and a proper choice of r = 0.2, the value of GP can be as high as 11. Thus by this theory is anticipated that this ESCP can generate up to 11 times more electrical energy in comparison with an equivalent conventional solar chimney power plant of the same dimensions.

Figure 3. The saturation relation between the chimney actual air speed and maximum air speed as function of the ratio r of cross section of operating air turbine openings to chimney cross section with various values of parameter a. Upper curve means greater α.

Figure 4. The geometric performance index GP as function of r for various values of the parameter a. Upper curve means greater α.

Such an extraordinary result begs a proper justification in terms of the true physical origin of the GP factor which we attempt to provide in the next paragraphs.

Many researchers in the field of solar chimneys technology have often questioned the reasons for the extremely low thermal efficiency of the conventional solar chimney power plants.

As a matter of fact the useful solar thermal power generated in the greenhouse entering at the bottom entrance of a solar chimney is given by, where is the moving air mass in kg/sec, C_{P} ~1004 and ΔT the temperature increase of the entering in the chimney air stream from the ambient air temperature.

It can be shown theoretically [9] that the maximum air speed inside the chimney is defined exclusively by the temperature difference ΔT and H, the height of the solar chimney, and is given by the relation

(6)

In (6), g = 9.81, T_{0} the ambient absolute temperature and k ~0.1 is a reasonable estimation for a well designed solar chimney with smooth internal surface and a ratio of height/diameter ≥8.

The maximum kinetic power of the rising air stream inside the chimney will given by. Thus the electric power, generated by the turbine placed inside the chimney, should be a part of this quantity, a fact that can be expressed with the aid of an appropriate efficiency coefficient as

(7)

The coefficient η_{T} expresses the overall efficiency of the turbine and electric generator hence we end up with the relation:

(8)

From relation (8) we see that the thermal power P_{th} generated by the greenhouse is eventually transformed to electricity, with an efficiency given by:

(9)

For g ~9.81, T_{0} ~300 K, η_{T} = 0.7 and a solar chimney height of H = 200 m we end up with an estimate of the order of. Considering furthermore that a double glassed greenhouse transforms a ~50% of the arising solar irradiation on its roof, the solar efficiency of the solar chimney technology, is in the range of ~0.2%. As a remark I mention that when a lower cost single glass greenhouse is used the efficiency drops by~33%, thus the economics of the solar chimney technology is in favor of double glass greenhouses.

The solar efficiency is defined as the electrical power generated divided by the solar irradiation on the roof of its greenhouse (or as the annual KWh of generated electricity divided by the annual solar irradiation on the roof of the greenhouse). The annual solar irradiation per m^{2} is a characteristic of the place of solar chimney power plant installation.

Thus we end up with a rather very low solar efficiency for an ordinary solar chimney power plant with a solar chimney height of 200 m.

This is the reason that in order to have an acceptable solar efficiency we should build solar chimneys of 1000 m (where the efficiency can become ~1.0%) and consequently accompanied by huge greenhouses due to the economics of solar chimney technology (see Appendix A). In this Appendix A, it is explained why the economics is forcing greenhouse dimensions to follow the huge dimensions of the solar chimney leading to a very expensive structure.

In order to evaluate the technical and operational feasibility and also the economic viability of the solar chimney technology we had to raise a huge and very expensive structure and spent for an experimental plant too much. If we compare this with the negligible funds in order to test PV technology or solar thermal technology made of a large number of similar units it becomes obvious why the solar chimney technology actually never tested after the Manzanares famous experimental solar chimney plant of Prof. Schlaich’s [1] .

This extraordinary small efficiency was theoretical supported by a Carnot cycle analysis [2] . Also the experimental power plant installed at Manzanares verified the proposed Carnot analysis. Actually the measured results were even lower due mainly to the fact that its turbine was operating without any control and thus its annual average efficiency(η_{T}) was much lower than 0.7.

The author’s objection to all these is due to the fact that because the air turbine was installed inside the solar chimney the low solar efficiency should be expected. The ducted air turbines, for low air temperatures, can transform into rotational power only a part of the kinetic power of the air stream moving inside their duct (i.e. the solar chimney). An excellent presentation of air turbines in low temperatures is given by Wright at [3] .

Thus it is the author’s opinion that the proposed Carnot cycle analysis has taken into consideration solely the internal solar chimney phenomena and not the atmospheric effects on the escaping air from the top exit of the solar chimney.

The true problem of such an approach is that in order to make the standard Carnot cycle analysis to agree with the empirical data from experimental tests, as for instance the figures shown in Prof. Schlaich’s book [1] that were obtained from a simplified analysis, one has to rely in the unjustified hypothesis that the air escaping from the chimney’s top, falls to the ground immediately after its exit.

Yet, a proper view using meteorological theory shows that this hypothesis is not true. In fact, the warm air column escaping by the solar chimney will rise up and can reach heights up to 4 km! Hence the previous Carnot model should be revised accordingly and the discrepancy should be corrected.

The hereby proposed new type of ESCP has already shown experimental indications that the efficiency can be increased say by a factor of GP ~10. In this case, an ESCP with a moderate chimney height of 200 m, could have a solar efficiency ten times higher i.e. ~2.0%. Thus with an ESCP of 200 m solar chimney we can transform 2.0% of the solar irradiation arriving on top of its greenhouse into electricity. The result, if proven, will give new prospects in solar chimney technology. We now suggest a possible physical explanation of this efficiency increase that was initially verified with a small experimental ESCP as will be presented in the next paragraph.

Several years ago an Australian company http://www.solartran.com.au/balloon_engine.htm suggested a rather unusual method for solar thermal electricity generation using uprising air balloons filled with hot air. The balloons could be filled with warm air made by a greenhouse and left to rise due to their atmospheric buoyancy. The balloons were then expected to reach a height at which the inside air density would approximately equal that of the ambient air at a certain height. The Australian concept was based in using ropes attached on the rising balloons to provide the mechanical energy in ground based equipment for electricity generation. The air temperature in the balloons of course would depend on the temperature difference achieved by the warm air generation greenhouse of this technology.

The same buoyancy phenomenon is happening also to the escaping by the solar chimney warm air column. To provide an estimate of the height reached by this warm air column we start from a standard relation for the warm air density in zero height (https://en.wikipedia.org/wiki/Density_of_air) given as

(10)

where ρ_{0} is the ambient air density at zero height. Let us assume that the thermal losses of the warm stream during its uprising are negligible thus the uprising can be considered as an adiabatic process. In this case in height H its density will become equal to.

This warm air density for “Aerostatic equilibrium” should become equal to the air density of the standard atmosphere at height H given by the relation.

Using this relation we can prove that even for a temperature difference of only 20˚C we get the “Aerostatic equilibrium” at a height around 4.3 Km. Humidity effects could force the uprising stream of warm air to climb even at higher heights. The air humidity was proposed as a possible means for a solar chimney technology without solar collectors as presented in [4] . Taking into consideration “aerostatic” and humidity effects we consider a more realistic meteorological approach for the solar chimney power plant operation.

Similar meteorological observations justify this effect for the case of uprising warm air above the oceans. Prof. Emmanuel presented a theory explaining the strong hurricanes as due to a gigantic Carnot cycle related to warm and wet air climbing above oceans at an average height of several km see for example [5] .

The same argument can be used to provide a qualitative explanation of the enhanced results that were calculated in the case of the improved ESCP due to the GP factor. The author’s present view is that by ESCPs we additionally exploit the escaping warm air column rising above the solar chimney. Thus we do not only use the kinetic energy of the air stream inside the solar chimney but also part of the dynamic energy due to the rising warm air. In order to test this hypothesis a small experimental ESCP installation was constructed in central Greece.

5. The Compotades Enclosed Solar Chimney Power Plant

Near a small village in central Greece, named Compotades, a small experimental solar chimney power plant was raised a video of its operation can be viewed at http://www.youtube.com/watch?v=RVJlM6spTjU.

Its solar chimney was a light metal skeleton structure with a height of 25 m and an internal diameter of 2.5 m. The lower 15 meters of the solar chimney were permanently wrapped by a thermally insulated plastic, while the wrapping of the top was folding fabric in order to test and receive measurements with 15 m and 24 m solar chimney.

The roof of the greenhouse was made of a clear plastic and had a surface area of 1020 m^{2}∙s depicted in the photo of Figure 5. At the entrance of the greenhouse a small turbine of 1.12 meters diameter (1 m^{2} cross section area) was placed as shown in Figure 6. This sole air turbine of the small power plant according to its specifications had a rating power of 300 W, for an air speed of 11 m/sec, in case of incoming air stream with a speed of ~4 m/sec its output electric power drops to ~10 ÷ 15 Watts.

Figure 5. A picture of the ESCP experimental plant during measurements.

Figure 6. The small air turbine used in the experiments.

As can be seen from the photographic evidence this structure was a very simple form of an enclosed solar chimney power plant (ESCP) because its greenhouse was closed firmly everywhere except of one circular opening (of 1.12 m diameter) at the side opposite to its solar chimney. At this opening the air turbine generator was placed.

The air turbine forced by an incoming stream of ambient air, was rotating during daylight continuously during the period of the experiments (on June 2013) and usually was operating and a few hours after the sunset. This stream of air was becoming warmer while moving inside the greenhouse and was updrafting through the solar chimney to upper layers of the atmosphere.

Actual implementation and measurements were completed at the small experimental ESCP for the period of June 2013.

In Table 1 we present the data of June 15 2013, at 12.00. We denote with U_{T}, the air velocity in m/sec at various openings made at the side opposite to solar chimney (where the small turbine finally was placed) and U_{Ch} the uprising air velocity inside the chimney for a chimney height (H) of 15 m. For a ratio of the opening area to the chimney cross section greater than one, we found an approximately constant chimney velocity of ~0.83 m/ sec.

In Table 2 we present similar measurements for an increased height of the chimney at ~24 m. Again for a ratio greater than one we find an approximate constant chimney air speed at ~ 1.05 m/sec.

The experimental results obtained reveal that we can definitely receive more power at the various air entrances, in stark contrast with the predictions of the Carnot previous conservative model [1] .

In fact the solar chimney maximum kinetic power for equals (ρ = 1.2) which for a height of 15 m where the speed was 0.83 m/sec, gives ~1.7 W. While for the opening of 1 m^{2} where the respective air speed was measured at 4.0 m/sec the kinetic power calculated by is ~38.4 W. The respective figures from Table 2 then for a chimney height of 24 m were calculated as ~3.3 W and ~74.0 W.

Thus in both cases of different chimney heights a GP higher than 20 was measured, that means that the air stream kinetic power in the turbine entrance before its placement was at least 20 times higher than the kinetic power of the rising air inside the solar chimney. This according to the classic Carnot cycle analysis [2] was not possible.

Furthermore measurements were taken, after the placement of the small rotating air turbine. It is obvious that the air turbine interference could lower GP, however the electric power for the 15 m solar chimney, for the estimated air speed of something lower than 4m/sec in the entrance of the air turbine, as estimated by its characteristic output curve was a figure between 10 to 15 W. Although this result was low it was much higher than the maximum kinetic power of the rising air stream inside the chimney (~1.7 W). Such a result was also unexpected from classical solar chimney theory, giving strong experimental evidence supporting the author’s theory applicable in the ESCP technology.

Thus the energy of the uprising warm air stream inside the solar chimney alone does not generate the electrical energy of the ESCP, but the plant is supplied by a much higher amount of dynamic energy due to the atmosphere. Thus the amount of the used energy by the ESCPs depends on atmospheric phenomena that are related in a more complicated way to the geometry of solar chimney, the greenhouse and the turbine openings of the power plant.

Table 1. Solar chimney height H = 15 m, Measured air speeds at turbine and solar chimney openings for varying ratio r of their surface areas.

Table 2. Solar chimney height H = 24 m, Measured air speeds at turbine and solar chimney openings for varying ratio r of their surface areas.

It is the author’s firm belief that the construction of a pilot ESCP of several KW will definitely test this new theoretical approach for the ESCP technology. This pilot ESCP could have a rating power of 18 KW and its main dimensions should be: A solar chimney H = 50 m and d = 5 m standing in the center of a single glassed hexagonal enclosed greenhouse ( as shown for example) in Figure 7, covering an area of 6000 m^{2}.

At the perimeter’s wall of the enclosed chimney six air turbines of 1.0 m diameter and 3 KW rating power each, will be placed. The turbines openings will be equipped with electronically controlled air stop systems. The estimated cost of this DEMO ESCP will not exceed 250,000 EURO.

Measurements on this pilot ESCP will evaluate the proposed theory and the new prospects and potential of the ESCPs for solar electricity generation.

6. Water Storage Modeling of an ESCP

Thermal study of a solar chimney power plant for its overall 24 h operating cycle is a complicated mathematical problem. In references [6] [7] methods of tackling this problem, taking into consideration the soil thermal storage ability were presented.

In Appendix B, a detailed method is presented, in which Fourier time series analysis is used for modeling the thermal behavior during a 24 hours cycle operation of an ESCP. With the proposed method we can tackle the presence of a water storage facility that is covering a part of the base of the ESCP greenhouse.

With the proposed code we can study the effect of partial coverage (10% of its area) of the greenhouse by the water storage facility, which is the proper thermal storage method in order to leave ground space of the greenhouse for other symbiotic activities beneath its glassed roof. However the reader can omit the detailed reading of Appendix B without missing the main arguments and results of the present paper.

The CODE was applied to a model ESCP with main dimensions as follows: Solar chimney of 200 m height (H) and 25 m diameter (d), Circular greenhouse of ~250,000 m^{2} surface area (A_{C}) and a diameter (Dc) of ~565 m, equipped with a circular concrete water storage facility of approximately ~187 m diameter and enclosed by a concrete short wall of approximately 1.1 m height and covered by a plastic cover.

The concrete structure water storage facility can be raised in the center part of the greenhouse and its water storage capacity will be ~25,000 m^{3}. In Figure 7, the estimated daily variation of the ambient temperature and achieved temperature increase ΔT in ˚C, and its daily % variations of the horizontal solar irradiation and the % calculated power output are shown (for GP ~11, a = 5 and r = 1/5), with and without the thermal storage facility. (Table 3)

The daily results are for a typical representative day (for example an April day) of average temperature of T_{0} = 23˚C and a daily variation of DT= ±4˚C in a place of an annual solar irradiation W_{Y} = 2000 kWh/m^{2}. (Figure 8)

By the Figure 8 it is evident that the ESCP, with the water thermal storage facility, operates more smoothly in comparison to the same ESCP without water thermal storage. Furthermore the rating power of the first could be 25% lower than its alternative, while both generate the same annual KWh estimated to ~11 GWh. Thus the increase in the construction cost of the first, due to its water thermal storage facility, is compensated by the cost increase of the second due to its higher power rating air turbines.

In order to test experimentally the water storage behavior of the proposed ESCP technology and its daily and annual electricity generation, a small demo ESCP could be constructed.

The proposed demo ESCP should be dimensioned as follows:

Solar Chimney:

1) Height H = 100 m.

2) Internal diameter d = 10 m.

Greenhouse:

1) Octagonal of D = ~185 m diameter.

2) Area Ac = ~24,000 m^{2}.

3) Double glass roof of > 2% inclination, entrance height ~3 m central height ~5 m.

Internal Water Storage Facility of ~60 m diameter, enclosed by a wall of 1.2 m with a water capacity of ~3000 m^{3}, covered by a plastic cup.

Air Turbines:

1) 16 air turbines, 1.25 m diameter, 7.5 KW rating electric power output.

2) 16 symmetric openings on the external eight sides of greenhouse wall.

3) 16 electro-mechanically controlled air barriers of the openings.

The estimated construction cost of this demo ESCP reaches approximately 1.6 million EURO analyzed as shown in Table 3. This demo ESCP can be installed in any place with an obvious preference for a sunny place of annual solar irradiation ~2000 KWh/m^{2}.

The electric generators of this demo ESPC can be connected to a nearby electric grid thus proper induction generators should be used.

Figure 7. An initial figure of the proposed small demo of 18 KW.

Figure 8. Daily operation in an April day of an ESCP with thermal storage covering ~10% of its ground with 1 m water and without water thermal storage.

Table 3. Dimensions and construction cost figures of the Model ESCP.

If no electric grid is present the demo ESCP can be operated as a local autonomous power plant in which case its electric generators should be synchronous and their output loads could be for instance ohmic resistances giving their generated thermal energy to nearby water tanks

7. Conclusions

The presented ESCP technology already tested in a small experimental power plant in Greece could be a major player in renewable technology market after its successful testing with relatively low cost demonstration plants.

The ESCP technology can combine low levelized cost of electricity with smooth 24 hours/day operation for 365 days/year.

The ESCP technology also has the following benefits in relation to its major competitive technologies (wind and PV and CSP):

・ No water consumption;

・ Zero effect on living nature and ambient;

・ Local manufacturing and use of common material (steel, glass, concrete and conventional axial fans as air turbines);

・ Easy and low cost maintenance;

・ Longer operating life up to 50 years (its solar chimney up to 100 years);

・ Their empty protected surface area beneath the glassed roofs of their greenhouses can be used partially for agriculture and or placement of parabolic trough mirrors for desalination.

Cite this paper

Christos D. Papageorgiou, (2016) Enclosed Solar Chimney Power Plants with Thermal Storage. *Open Access Library Journal*,**03**,1-18. doi: 10.4236/oalib.1102666

References

- 1. Schlaich, J. (1995) The Solar Chimney: Electricity from the Sun. Axel Mengers Edition, Stuttgart.
- 2. Gannon, A. and Von Backstrom, T. (2000) Solar Chimney Cycle Analysis with System loss and solar Collector Performance. Journal of Solar Energy Engineering, 122, 133-137.
- 3. Wright, T. (1999) Fluid Machinery, Performance Analysis and Design. CRC Press.
- 4. Papageorgiou, C.D., Psalidas, M. and Katopodis, P. (2010) Solar Chimney Technology without Solar Collectors. SCPT 2010 Solar Chimney Power Technology International Conference, Bochum, Germany, 28-30 September 2010, 275-282.
- 5. Emanuel, K.A. (1986) An Air-Sea Interaction Theory for Tropical Cyclones. Part I: Steady State Maintenance. Journal of the Atmospheric Sciences, 43, 585-604.

http://dx.doi.org/10.1175/1520-0469(1986)043<0585:aasitf>2.0.co;2 - 6. dos S. Bernades, M.A., Vob, A. and Weinrebe, G. (2003) Thermal and Technical Analyses of Solar Chimneys. Solar Energy, 75, 511-524.

http://dx.doi.org/10.1016/j.solener.2003.09.012 - 7. Pretorius, J.P. and Kroger, D.G. (2006) Solar Chimney Power Plant Performance. Journal of Solar Energy Engineering, 128, 302-311.
- 8. Gengel, Y.A. (1998) Heat Transfer. McGraw Hill, New York.
- 9. Papageorgiou, C.D. (2010) Floating Solar Chimney technology. In: Rugescu, R.D., Ed., Solar Energy, Intech, 187-222.

Appendix A: Cost Optimization of the Solar Chimney Technology

Let us consider a Solar Chimney Power Plant (SCPP) generating annually a certain amount of electricity in KWh defined by the symbol W, and made of three major parts:

・ A Solar Chimney of height H and an average internal diameter d

For stability reasons we assume that

・ A circular Greenhouse of area Ac

・ An electricity generating system of air turbines of rating power P = W/h (h = 8760 X the capacity factor of the technology)

The basic construction costs of the three main parts are approximately given by the relations:

SCC = Solar Chimney Cost =

GHC = Greenhouse cost =

ATC = Air Turbine cost =

Thus the overall SCPP cost (TC) is equal to the following sum

where ε represents the overhead and unforeseen costs. According to previous presented theory We may take a constant W and ask under what condition TC is minimized. In order TC to be minimum the sum should also be minimized. Replacing the following function of H should be minimized, are independent of H

From the above one obtains

or

i.e.

Conclusively, in order to minimize the construction cost for constant annual KWh generation of a solar chimney power plant, its greenhouse cost should be approximately twice the cost of its solar chimney.

Appendix B: ESCP Thermal Heat Transfer Analysis

1. The ESCP Model

The ESCP model with a circular collector is shown in the indicative diagram in Figure B1. The circular solar collector of this ESCP is divided in a series of M circular sectors of equal width Δr as shown in Figure B1.

Figure B1. The cut of a circular sector of the solar collector of the SCPP.

The m^{th} circular sector (m = 1 to M) will have a width, an average radius and an average height. For a linear variation of the roof height, where: D_{c} = solar collector diameter and D_{in} = Final internal diameter of the solar collector. If not available data are given D_{in} ≈ 2∙d.

These consecutive circular sectors, for the moving stream of air of mass flow, are special tubes of almost parallel flat surfaces and an equivalent diameters. The ambient air as moving towards the entrance of the first circular sector is assumed that it is increasing its temperature T_{0} to T_{0} + dT .The exit temperature of the first sector is the inlet temperature for the second etc. and finally the exit temperature of the final M^{th} sector is the T_{03}, i.e. the inlet temperature to the air turbines.

The m^{th} circular sector (m = 1 to M) will have a width, an average radius and an average height.

For a linear variation of the roof height, where: D_{c} = solar collector diameter and D_{in} = Final internal diameter of the solar collector. If not available data are given D_{in} ≈ 2∙d.

These consecutive circular sectors, for the moving stream of air of mass flow, are special tubes of almost parallel flat surfaces and an equivalent diameters.

The ambient air as moving towards the entrance of the first circular sector is assumed that it is increasing its temperature T_{0} to T_{0} + dT.

The exit temperature of the first sector is the inlet temperature for the second etc. and finally the exit temperature of the final M^{th} sector is the T_{03}, i.e. the inlet temperature to the air turbines.

2. Daily Variations of Ambient and Sky Temperature and Solar Irradiation

The ambient (T_{0}) and sky (T_{sk}) temperatures are in general functions of time. If there are not sufficient data for a place of ESCP installation the following formulae in ˚C and t in hours from midnight can be used for the average annual representation day:

For the average day if not available data, T_{0} = 20˚C.

as given by Bernades M.A. dos

S., Vob A., Weinrebe G. in [6] .

where T_{dp}(t) is the due point temperature in ˚C given approximately by the formula:

, where RH is the relative humidity of the ambient air %. Without any avail-

able data, RH can be assumed to be equal in average to 50%, i.e. in average T_{dp}(t) = T_{0}(t) − 10 in ˚C.

The solar irradiation on horizontal surface is a very important characteristic of the place of SCPP’s installation. The annual solar irradiation per m^{2} on horizontal surface W_{y} is a characteristic figure of the place and should be given or estimated approximately by existing solar maps. The average horizontal solar irradiance on the sector roof is equal to the annual solar irradiation divided by 8760 hours. However for the average representation day of the year, if other data are not available, the horizontal solar irradiance G(t) can be assumed as a sinusoidal function of time between the 6^{th} and the 18^{th} hour and zero outside, hence:

where: G_{M} = π∙W_{y}/8760. If data are available T_{0}(t), T_{dp}(t), T_{sk}(t), G(t) can be described by their real equivalent daily functions for the representative average day of the year or any other day. As representation day usually is received the day of which its calculated daily electrical energy production is equal in approximation to the average daily electrical energy. Thus the annual electricity production by the ESCP will be 365 times this average daily electricity production. As a general rule an April the day’s duration is very close to the average representation day. The solar irradiation energy per m^{2} arriving on the ground of the sector is given approximately by t_{a}∙G(t), where t_{a} is the mixed coefficient {overall covers transmissivity ground absorptivity} for short wavelength solar light irradiation. If there are not available data for the cover and curtain for one cover the overall transmissivity is estimated to 0.9 while for two covers is estimated to 0.81. The absorptivity for treated ground if data are not available is estimated to 0.9.

Without any available data, RH can be received as equal in average to 50%, i.e. in average T_{dp}(t) = T_{0}(t) − 10 in ˚C.

The solar irradiation on horizontal surface is a very important characteristic of the place of SCPP’s installation. The annual solar irradiation per m^{2} on horizontal surface W_{y} is a characteristic figure of the place and should be given or estimated approximately by existing solar maps. The average horizontal solar irradiance on the sector roof is equal to the annual solar irradiation divided by 8760 hours. However for the average representation day of the year, if other data are not available, the horizontal solar irradiance G(t) can be assumed as a sinusoidal function of time between the 6^{th} and the18^{th} hour and zero outside. Thus:

where: G_{M} = π∙W_{y}/8760.

If data are available T_{0}(t), T_{dp}(t), T_{sk}(t), G(t) can be described by their real equivalent daily functions for the representative average day of the year or any other day. As representation day usually is received the day of which its calculated daily electrical energy production is equal in approximation to the average daily electrical energy. Thus the annual electricity production by the ESCP will be 365 times this average daily electricity production. As a general rule an April day is very close to the average representation day. The solar irradiation energy per m^{2} arriving on the ground of the sector is given approximately by, where t_{a} is the mixed coefficient {overall covers transmissivity ground absorptivity} for short wavelength solar light irradiation. If data are not available for the cover and curtain for one cover, the overall transmissivity is estimated to 0.9 while for two covers is estimated to 0.81. The absorptivity for treated ground if data are not available is estimated to 0.9.

3. Heat Transfer Analysis of the Basic Two Glazing Model

The mass flow of the sector could be varying during the daily cycle or constant. In the present analysis the mass flow is considered constant, as the value that gives to the ESCP its maximum electricity generation. An initial estimation for is, where air speed is υ is estimated to 7 - 8 m/sec and the air density given by ρ = p_{0}/(287 × 307.15). The internal solar chimney diameter is constant and equal to d.

Following the circular sector diagram for a double glazing solar collector, as shown in Figure 2, the following equations are derived:

・ For outer cover (w), where:

・ For inner cover or curtain (c),

・ For the moving air, where:

・ For ground surface (s) where is the heat power going down to the ground.

The first three equations form a linear system as follows:

Through of which the temperatures or outer cover T_{w}, of the curtain T_{c} and the air T can be expressed as linear functions of the ground temperatures T_{s}. and the irradiance G are functions of time given in the previous paragraph. These linear forms are in general in the form:

,

In general c_{1}, c_{2}, c_{3}, c_{4}, c_{5} and c_{6} are functions of time.

Replacing T_{c} and T_{in} the ground equation:

The following relation is derived

,

where

We can always write for each cover and the moving air a similar number of linear relations between the temperatures with heat transfer parameters as coefficients and always we can solve the linear system in order to derive a relation for the ground surface (s), connecting the ground temperature T_{s} the power heat flow down to the ground surface q_{gr}. In general all the functions appearing on the ground relation are functions of time (t). In order to calculate the time functions T_{s}(t) and q_{gr}(t) it is necessary to have one more relation between them. This relation it will be derived by the ground heat flow partial differential equations as explained in the next paragraph.

4. Ground Equivalent Equation in Fourier Space

The ground bellow the sector is considered as infinitely wide thus the surface temperature it depends in the depth distance (z) and time (t). It is also considered homogeneous and characterized by three parameters, its density ρ_{gr}, its specific heat capacity c_{gr} and its thermal conductivity k_{gr}. The ground partial differential equations, relating the surface temperature T_{s}(z,t) and the down-moving heat power q_{g}(z,t), starting from its surface (z = 0) and going down (see Figure 2), can be written as the following set of two equations:

where q_{gr}(z) for z → ∞ is zero.

Let us consider that T_{s}(z,t) and q_{g}(z,t) are sinusoidal functions of time as follows:

Thus the set of two differential equations becomes:

This set represents an infinite homogeneous transmission line with “equivalent voltage” the function T_{s}(z) and “equivalent current” the function q_{gr}(z).

According to transmission line theory this infinite transmission line is equivalent to its admittance defined by the relation:.

Thus the Fourier transform of heat flow q_{gr} for z = 0, is related to the Fourier transform of the Ground temperature T_{s} in the surface z = 0 by the linear relation:

For the daily variation the ground temperature in z = 0, T_{s}(t) can be analysed in a Fourier series with a fundamental harmonic ω_{0} = 2 × π/(24 × 60 × 60) as follows:

For each harmonic the previous is valid for ω = n∙ω_{0} thus:

In order to take into consideration the water thermal storage effect, let us assume that the equivalent water height dh is very small, with an assumed very high thermal conductivity thus its temperatures in both of its surfaces are almost equal to the ground temperature T_{s}. Under these reasonable approximations the water thermal storage sheet of thickness dh is equivalent in Fourier space to an admittance

, where c_{w} = water specific heat = 4,186,800.

Thus in case of ground plus artificial water thermal storage of equivalent height dh the linear relation for the Fourier transform harmonics of temperature T_{s}_{ }and heat flow q_{gr} are given by the relation: Where:._{ }

5. Fourier Analysis of the Sector

On the ground surface (s) the time equation relating the ground temperature T_{s}(t) and the power heat transfer to the ground q_{gr}(t) is given by the equation:

In the equation two known functions q(t) and h(t) appear.

These functions can be expressed by the following Fourier series:

,

While the Fourier series of the function q_{gr}(t) is given below where the Fourier harmonics of the function q_{gr}(t) have been replaced by the harmonics of the T_{s}(t).

We can replace the functions with their Fourier series from −M_{2} up to M_{2} harmonic, in the ground relation. The bilinear form of the produced equation is equivalent to a set of 2 × M_{2} + 1 linear relations. The equations are derived equating the coefficients of the time functions exp(j∙n∙ω_{0}∙t) for −Μ_{2} ≤ n ≤ +M_{2} and considering the coefficients of these functions negligible for n > M_{2} and n < −M_{2}. This is a very reasonable estimation due to the fact that the harmonics of all time functions are rapidly decreasing.

For example for M_{2} = 2 the following matrix equation it is the set of relevant equations:

The matrix equation can be solved and T_{s}_{,n} can be calculated for −Μ_{2} ≤ n ≤ +M_{2}.

Thus T_{s}(t) can be calculated by its Fourier time series as follows:

It is suggested to use the Real part operator, before the Fourier sum, in order to get rid of the negligible imaginary parts due to computer and method approximations.

Usually it is enough to calculate N values of the previous function for: i.e. for 1 ≤ ν < N. Using the linear relations connecting T_{w}(t), T_{c}(t) and T(t) with T_{s}(t) these temperature time functions can be for the same time sequence.

The temperature dependent heat transfer coefficients can be re estimated and the same program can be run again defining the temperatures with a better accuracy etc. (usually two to three times running is enough).

The equation defines the exit air temperature by the sector that is equal to the inlet air temperature to the next sector etc. The exit temperature of the final sector is equal to the inlet air temperature to the air turbine T_{03}(t) for the given mass flow. The number of sectors can be between 10 and 50 for acceptable accuracy in T_{03}(t)calculation.

Thus using an approximate proportion coefficient relating the heat Power output (defined by) of the greenhouse as function of time for 24 hours daily period the Electric Power output can be calculated. The code can be used to calculate the best daily average value for the average daily mass flow in order to maximize the power output of the ESCP. Using control methods through the openings we can vary this average mass flow during the 24 hours daily period in order to achieve an optimal performance of the ESCP.

M_{2} defines the number of harmonics, usually M_{2} can received in the range of 5 to 10, for all the sectors, while the time points ν can be a reasonable multiple of 24 (96 is enough).