Draught vs. tower height for a contraction ratio 10.
1. Introduction
The thermal draft principle is currently used in exhaust chimneys to enhance combustion in domestic or industrial heating installations. An introductory level theory of gravity draught in stacks was issued by the old German research institute for heating and ventilation (HermannRietschelInstitut) in Charlottenburg, in a widely translated reference book (Raiss 1970). Technological and practical aspects of air draught management are clearly exposed in this works, but a widepredicting theory still lacks. As early as in 1931 a surprisingly advanced proposal to use thermal draught as a propelling system to generate electricity from solar energy was forwarded by another German researcher (Günter 1931). Major advancements in convective flows prediction during the last decades of the 20th century were accompanied by a series of publications and we cite first the basic book due to a work from Darmstadt (Unger 1988). The related topic of convective
However, a simple and efficient solution exists which is here demonstrated by means of energy conservation. This method provides a superior energetic efficiency with moderate costs and a high reliability through simplicity. It consists of optimally heating the freshair by means of a mirror array concentrator and an efficient solar receiver, and accelerating it further in the tall towers through gravity draught (Fig. 1, Rugescu 2005).
This genuine combination has already a history of theoretical study (Rugescu 2005) and an incipient experimental history too (Rugescu et al. 2005). First designed for air acceleration without any moving parts or drivers with application to infraturbulence aerodynamics and aeroacoustics, the project was further extended for green energy applications along a series of published studies (Rugescu et al. 2006, Rugescu 2008, Rugescu et al. 2008, Rugescu et al. 2009, Rugescu et al. 2010, Cirligeanu et al. 2010, Rugescu et al. 2011a, Rugescu et al. 2011b, Rugescu 2012, Rugescu et al. 2012a, Rugescu et al. 2012b). The demonstration of the high draught tower energetic efficiency provided below is expected to convince the skeptics and to bolster again the direct solar energy exploitation in tall tower power plants (Rugescu et al. 2012 b).
2. Gravitydraught accelerator modeling
A schematic diagram of a generic draught tower is drawn in Fig. 2. The fresh air in its ascending motion up the tower, due to the gravity draught, is first absorbed, from the immobile atmosphere (
The ideal gas behavior under the influence of a gravity field of intensity g, flowing upward with the local velocity w into a vertical duct of cross area
The air flow of the material, infinitesimal control volume
wher
The computational solution of the stack flow further depends on the initial and limit conditions that must fit the physical process of thermal draught (Bejan 1984) and may be managed in simple thermodynamic terms. In its general form, the dynamic equilibrium of the stack flow was first debated in a dedicated book (Unger 1988), with emphasize on the static pressure equilibrium within and outside the stack at the openings, the key of the entire stack problem. The onedimensional steady flow assumption with negligible friction was accounted and we add the proofs that this approach is consistent with the problem. In that regard we analyze in a new way the flow with friction losses, estimate their magnitude and add a different accounting for compressibility at entrance. Our point of view faintly modifies the foregoing results regarding the compressibility of the air during inlet and exit acceleration, still consists of a necessary improvement.
The aerostatic influence of the gravitation is then given by the pressure gradient equation inside (density ρ) and outside (density ρ_{0}) the tower,
The right hand term in these equations is nothing but the slope to the left of the vertical in each pressure diagram from Fig. 3.
This means that the inner pressure in the stack (left, doted line) is decreasing less steeply and remains closer to the vertical than the outer pressure of the atmosphere. The dynamic equilibrium is established when, following a series of transforms, the stagnation pressures inside and outside become equal (Fig. 3). While the air outside the stack preserves immobile and due to the effect of gravitation its pressure decreases with altitude from
Under the assumption of a slender tower with constant cross area
and energy for the compressible flow, with the assumption
for the entrance into the stack. In other words the air acceleration takes place at tower inlet between 01 as governed by the energy compressible equation with constant density
where
The air is warmed in the heat exchanger/solar receiver between the sections 12 with the heat
The impulse equation gives now the value of the pressure loss due to air dilatation,
where a possible pressure loss into the heat exchanger Δ
equivalent to
The gravitational effect (4) continues to decrease the value of the inner pressure up to the exit rim of the stack, where the inner pressure becomes
Either the impulse equation in the form
or the energy equation in the form
appears for the receiver, heated zone, and
for the free ascending flow above the receiver, with
At the upper exit from the stack the gas is diluting and braking into the still atmosphere, thus the compressible Bernoulli equation applies,
when constant density during this process is assumed again. The pressure variation at stack’s exit is very small and this ends in the fact that other simplifying hypotheses do not give results consistent with the physical phenomena.
Modifying eq. (14) the inner static pressure at stage
where the relative heating of the air is expressed in terms of densities,
with a given control value for
Using eq. (17) the static pressure of the exhausted air becomes
which is used in the equilibrium condition as follows.
The values of the pressures and velocities into the main sections result from the equilibrium condition of the pressures above the upper exit, where the inner
In this way (Unger 1988, Rugescu 2005, Rugescu et al. 2005), the mass flow rate through the stack mainly depends on the relative heating of the air, expressed in terms of densities, and results when the pressure difference between the interior and the exterior of the tower exit recovers by dynamic braking of the air (Fig. 3).
For negligible friction losses (
slightly higher than the predicted value of the previous models (Unger, 1988).
When the friction losses are considered, the actual value for the quadratic mass flow rate results from the second degree equation (2212) which gets the form,
where at the nominator a reference freefall mass flow rate appears,
based on the Torricelli freefall velocity
with the constants
For an example slender, tall stack with the inner channel of elongation
meaning that the difference from the frictionless flow is actually smaller than 0.5 ‰. Consequently the nonfriction result in (2313) should be considered as accurate. Its quadratic form shows the known fact that the heating of the inner air presents an optimal value and there exist an upper limit of the heating where the flow in the stack ceases.
Formula (2313) shows that the nondimensional quadratic mass flow rate
and is given by
The entrance speed exhibits a maximum at the theoretically optimal heating
namely
The optimal heating for the standard air appears at a relative density reduction
meaning an equal increase of the absolute temperature of (1+
A slightly improved model is delivered when the following conditions at the upper exit are introduced, starting from equation (19). The constant density assumption along the upper stack
is here replaced with the condition (Unger 1988) of an isobaric exit
This means that the dynamic equilibrium is reestablished when the stagnation pressure from inside the tower equals the one from outside, at the exit level,
This equation is the end element that allows determining the equilibrium value of the air mass flow rate passing through the stack. Equaling (20) and (21),
Reducing by the quotient
Depending on the construction of the heat exchanger the drag largely varies. For simple, tubular channels the pressure loss due to frictions stands negligible (Rugescu et al. 2005a, Rugescu 2005, Rugescu et al 2005b) and the reduced mass flow rate (RMF) results from the simple equation
It gives an alternative to the previous solution of Unger (Unger 1988)
or to the one from above (Rugescu et al. 2005a)
and gives optimistic values in the region of smaller values of heating (Fig. 4).
The behavior of the chimney flow for various heating intensities of the airflow, in the limit case of equal far stagnation pressures (FSP) and for the three different models described is reproduced in Fig. 4, where the limiting, linear cases of the dynamic equilibrium are drawn through straight, tangent lines. These are in fact the derivatives of the mass flow rate in respect to
Differences between the present solution and the previous ones, as given in the above diagram, are nonnegligible and show the sensible effect of the variation in modeling of the compressibility behavior at entrance and exit of the stack. This is explained by the tinny variations in pressure and density during the very small acceleration of the air at tower inlet that makes the flow highly sensible to pressure perturbations, either natural or numerical. The same applies for the tower exit. For this reason the previous solution was obtained by completely neglecting the air compressibility at tower upper exit, where the static pressure was taken into consideration instead of the dynamic one.
Numerical simulations of the ducted airflow and the experimental measurements on a scale model support of the present model. The conclusion of this very simplified but efficient modeling of the selfsustained gravity draught, with no energy extraction, is that the heating of the air must be limited to between 0.3÷0.5 in terms of the relative density reduction through heating, or to between 90÷150ºC in terms of air temperature after heating, because under the accepted assumptions the product
The accelerating potential and the expense of heat to perform this acceleration at optimal conditions result from equations (37)÷(39). In a practical manner, the velocity
In fact these formulae render identical results for the optimal values for










7  11.72  4.85  8.28  82.8 
14  16.57  6.86  11.72  117.2 
30  24.26  10.05  17.15  171.5 
70  37.05  15.35  26.20  262.0 
140  52.40  21.71  37.05  370.5 
The value of
3. Experimental results
With the existing smallscale test rig built by the team of University “Politehnica” of Bucharest, the tests that have been conducted led to the values for air velocity in the tube as given in the diagram below. The average values, measured at a distance of 1.7 m from the entrance area of the tube, were registered as 2.115 m/s air speed with the contracted area effect (simulation of a turbine) and of 6.216 m/s without turbine simulation. Air temperature at the exit section was recorded to be of 195^{o}C and 123^{o}C, respectively (Tache et al. 2006).
The turbine simulation and the image of the inner electrical heater, simulating the solar receiver, are shown in figures below.
The experimental values recorded during the measurement session and the ones obtained from numerical simulations are listed in Table 2.




1  2.101  5.703  2.714  Vmin = 2.19  Vmin = 5.90  
2  2.190  7.110  3.247  
3  2.051  5.767  2.811  
4  1.996  7.310  3.662  
5  2.127  6.920  3.253  
6  1.867  5.521  2.957  
7  2.414  7.208  2.986  
8  2.027  5.966  2.943  
9  2.051  5.703  2.780  Vmax = 3.29  Vmax = 7.07  
10  2.307  6.920  3.000  
11  2.276  5.966  2.621  
12  2.027  5.351  2.639  
13  2.076  5.767  2.778  
14  2.276  6.329  2.780  
15  1.937  5.703  2.945  







The differences between these values are small, with greater values (~29.55%) when accounting for the turbine effects and much smaller values (~4.33%) in the other case.
4. Design example
As already stated, the optimal air heating for a good draught effect (Fig. 4) stays between 50÷100º
the discharge of the stack exhibits a good value of 2/3 of the maximal one,
At half of the optimal heating, that means at 100º
Under these circumstances it is fairly reasonable to accept for the further computation a moderate rarefaction of
where the density of the air is still the normal one
Considering now a rough constant pressure specific heat of the air of
the power consumed with the heating of the air raises to
Under a global heating efficiency of 80% the required total solar irradiation is
The lunaraveraged solar irradiation in Bucharest with the daily and annual values respectively are given below (University of Massachusetts 2004),
for a local horizontal surface, under averaged turbidity conditions. From the ESRA database, the value of 3.7 results. In the same database, the optimal irradiation angle is given equal to 35º, although the local latitude is 45º. The difference is coming from the Earth inclination to the ecliptic. As far as the mirror system is optimally controlled, the radiation at the optimal angle must be accounted, as equal to:
and the mean diurnal insolation time at the same location in Bucharest equal to
The following solar irradiation intensity received during the daylight time results
The reflector area, directly facing the Sun results, with the value of
Due to different angular positions of the mirrors versus the straight direction to the Sun, due to their individual location on the positioning circle, at least 50% extra reflector area is required to collect the desired radiating power from the Sun, or
When 3
In contrast to the natural gravity air advent, when a turbine or other means of energy extraction are present, the characteristic of the tower suffers a major change however. The tower characteristic includes now the kinetic energy removal by the turbine under the form of externally delivered mechanical work.
5. Turbine effect over the gravitydraught acceleration
The turbine could be inserted after or before the air heater. For practical reasons, the turbine block is better imbedded right upwind the solar receiver (Fig. 9), forcing the raising of the position of the receiver and thus a better insolation of the heater along the whole daylight.
According to the design in Fig. 16, a turbine is introduced in the SEATTLER facility next to the solar receiver, with the role to extract at least a part of the energy recovered from the sun radiation and transmit it to the electric generator, where it is converted to electricity. The heat from the flowing air is thus transformed into mechanical energy with the payoff of a supplementary air rarefaction and cooling in the turbine. The best energy extraction will take place when the air recovers entirely the ambient temperature before the solar heating, although this desire remains for the moment rather hypothetical. To search for the possible amount of energy extraction, the quotient
To describe the model for the air draught with mechanical energy extraction we shall resume some of the formulas from above. First, the process of air acceleration at tower inlet is governed by the same incompressible energy (constant density
The air is heated in the solar receiver with the amount of heat
the continuity condition shows that the variation of the speed is given by
No global impulse conservation appears in the tower in this case, as long as the turbine is a source of impulse extraction from the airflow. Consequently the impulse equation will be written for the heating zone only, where the loss of pressure due to the air dilatation occurs,
A possible pressure loss due to friction into the lamellar solar receiver is considered through Δ
The dilatation drag is thus perfectly identified and the total pressure loss Δ
Observing the definition of the rarefaction factor in (54) and using some arrangements the equation (58) gets the simpler form
The thermal transform further into the turbine stator grid is considered as isentropic, where the amount of enthalpy of the warm air is given by
If the simplifying assumption is accepted that, under this aspect only, the heating progresses at constant pressure, then a far much simpler expression for the enthalpy fall in the stator appears,
To better describe this process a choice between a new rarefaction ratio of densities
The air pressure at stator exit follows from combining (62) and (59) to render
Considering the utilization of a Zöllytype turbine, its rotor wheel keeps thermally neutral by definition and thus no variation in pressure, temperature and density appears in the rotor channel. The only variation is in the direction of the air motion, preserving its kinetic energy as constant.
Thus the absolute velocity of the airflow decreases from the value
The air ascent in the tube is only accompanied by the gravity updraught effect due to its reduced density, although the temperature could drop to the ambient value. We call this quite strange phenomenon the
The simplification was assumed again that the air density varies insignificantly during the tower ascent. The value for
The energy equation in the form of Bernoulli is used to retrieve the stagnation pressure of the moving air above the upper exit from the tower, under incompressible condition when the density remains constant,
Value for
It is observed again that up to this point the entire motion into the tower hangs on the value of the mass flowrate, yet unknown. The mass flowrate itself will manifest the value that fulfils now the condition of outside pressure equilibrium, or
This way the air pressure at the local altitude of the outside atmosphere equals the stagnation pressure of the escaping airflow from the inner tower. Introducing the equation (68) in equation (69), after some rearrangements of the terms, the dependence of the global mass flowrate along the tower, when a turbine is inserted after the heater, is given by the developed formula:
where the notations are again recollected
All other variables are already specified in the previous chapters. It is clearly noticed that by zeroing the turbine effect (
For different and given values of the efficiency
6. Discussion on the equations
Notice must be made that the result in (70) is based on the convention (60). The exact expression of the energy
The larger term here is the ratio
7. Discussion on denominator
The expression from the denominator of the formulae (70), which gave the flow reportedly, it can be canceled (becomes 0) for the usual values of the dilatation rapport (ratio) gamma and respectively quota part from energy extracted omega. This strange behavior must be explained. The separate denominator in (72) is,
The curve of zeros and the zones with opposite signs are:
It is yet hard to accept that such a selfamplification or pure resonance of the flow can be real and in fact the formulae (71) does not allow, in its actual form, the geometrical scaling of the tunnel and of the turbine. The rigor of computational formulae is out of any discussion, this showing that the previous result outcomes from the hypotheses adopted. Among those, the hypothesis of isobaric heating before the turbine is obviously the most doubtful.
8. Improved model
Analyzing the simple draught only, observe how easily the hypothesis of isobaric heating leads to an incomplete result, by eliminating the drag produced by the thermal dilatation and the acceleration throw heating, thus reducing the problem to a linear one, without physical anchorage. It could be presumed that the acceptance of relation (57) for the cooling in the stator, relation where it was presumed that the anterior heating performed isobaric, induces an excessive rigidity in the computational model. Replacing this very simple relation between the temperatures and the heat added to the fluid through a nonisobaric relation complicates drastically the model, which becomes completely nonlinear.
It remains to be analyzed whether such an inconvenient model leads to physically acceptable results for the values of mass flowrate in the turbine tower.
The isobaric relation (60) will be replaced by the exact equation,
where the heat received in non isobaric heat exchanger is expressed, through the equation of energy, in the complete form:
to take also into account the possible pressure losses due to friction in the solar receiver ∆
The absorbed heat (74) will also be used in its complete form in the relation that supplies the pressure at stator exit:
fact that obviously induces another level of nonlinearity. Using also the equation of state, the pressure from the stator exits writes from (71),
and for the value
Taking into account the draught from the tower (66) and the fluid brake at exit (67), the equilibrium of static pressure reads
Here the notation was used:
In the followings the undimensionalised flowrate
where also naturally appears the characteristic velocity
The characteristic velocity
in connection with which the relative flowrate couls also be expressed, in the form
in other words this flowrate is proportional to the squared local Mach number of the flow.
From (78) the equation of the flowrate
where the constant coefficients are again reproducing those working conditions,
The algebraic, non linear equation (80) is now solved using a standard numerical method to obtain solutions for the mass flowrate, as depending on the different working conditions concerning the heating level applied in the solar receiver
It proves however that the above given model is not properly reproducing the StackTurbine (ST) characteristic at low heating rates (
or
A very slight change in the assumptions could therefore deeply affect the result of the simulation modeling, due to the small overall magnitudes of pressure and density gradients along the ST channel.


0  3,50 
0,1   
0,2  (1,280) 
0,3  (0,875) 
0,4  0,611331000 
0,5  0,428261326 
0.6  0,298397500 
0,7  0,199248700 
0,8  0,1098315 şi 0,012898027 
0,9  0,0634500 şi 0,055130000 
1,0  0,00 
The results are plotted in the diagram from Fig. 11. The discharge characteristic of the tunnel resulting from the given assumptions is drawn in dark red.
9. Energy output of the gravitydraught accelerator
The main concern and reluctance for the classical solar towers comes from the regular perception that the energetic efficiency of those systems is unsatisfactory. Largely correct, this perception does not further stand valid for gravity draught towers and to prove this a piece of attention must be allocated to the energy balance.
The equation of energy in its rough form (3) needs thus further attention. Pointing the values to the exit station “2” of the receiver (Fig. 2) we first observe that the gain in kinetic energy
is equal to
where the first righthand term is the total heat introduced into the stack per one kg of air, and the last term represents the heat consumed for directly heating the air to the final temperature
The quantity of kinetic energy transferred to the air is the difference that remains available. This entire amount could be used to produce energy, without any thermal or mechanical loss. Physically, the heat introduced in the air to create the updraught along the tower could entirely be extracted into useful mechanical work through a low temperature wind turbine, and the draught is maintained due to the low air density despite the energy extraction in the tower.
The process remains however greatly dependent to the optimal selection of the heating level and of the utilization of the solar radiation in an efficient manner. The problem with the cloudy weather and the energy stocking during the night are solved through heat accumulators of specific construction.
10. Conclusion
The principle of a solar energy power plant, based on a mirrortype collector, is depicted in the nearby drawing. It represents the application of the WINNDER thermal accelerator principle into the ecological and sustainable means of accelerating the air without any moving device and, consequently, with a very low noise and turbulence level, ideal for aeroacoustic applications. A multiplerows array of controllable ground mirrors are installed around. In this manner a highly efficient utilization of the solar energy is available, due to the known high release coefficient of the mirror surfaces. Means to follow the Sun along its apparent trajectory are common and available at low cost today. Problems regarding the maintenance of the system can be solved through a proper technological and economic management of the facility.
It does not seem however equally attractive for energy production, despite the clean method involved, but this represents a first sight impression, easily dismounted through an indepth analysis. The computational model depicted above shows that the resources for producing energy trough the solar gravity draught are high enough and represent an interesting resource of green energy of a new and yet unexplored type.
Although the equipment costs of the present project are much higher than the Greenhouse power plant ones, it is believed that the overall costs are still competitive and the proposed solution of reflector tower is useful. One of the explanations resides in the fact that the reflexivity of the mirrors is very high. The design example given above highlights the main factors.
This example shows that a circular ground surface of roughly 0.8
After the data in (Schleich et al. 2005) this value equals 0,94 and this adds to the very high absorbing properties of the tower walls. It serves here as a nice illustration of possible extra applications of the chimney draught effects in directly producing electrical power.
As another comparison item, the newly renovated Solar Two solar thermal electric generating station, located in California’s Mojave Desert, consists of 1,900 motorized mirrors surrounding a generating station with 10 megawatts of capacity, which began operation in early 1996. It is part of an effort to build a commercially viable 100
The molten salt (an environmentally benign combination of sodium nitrate and potassium nitrate) allows a summer capacity factor as high as 60%, compared with 25% without storage. The plant consists of 1,926 motorized mirrors focused on a 300
Equipment costs of WINNDER are higher than for the Greenhouse power plants, still the overall costs of exploitation and maintenance are competitive and the proposed combination of mirror array and draught tower is literally efficient. It remains to convince the investors of the efficiency of this exotic energy producer.
The gravitational updraught due to Archimedes’s effect does not contribute, in any way, to the balance of energy. It simply remains the driver of the air into the stack and the solar energy introduced in the system is the only source of air acceleration and further production of electric energy within a turbogenerator. Consequently it does not seem specifically attractive for energy production, although it provides the cleanest energy ever and involves the lowest levels of losses.
References
 1.
Bejan, A., Convection Heat Transfer, New York, Wiley and Sons, 1984.  2.
Cirligeanu, R., Rugescu, R. D., and AlinaBogoi (2010), TRANSIT code for turbine flows in solargravity draught power plants, Paper GT201022518, Proceedings of the ASME International Conference on Gas Turbines, June 1418, 2010, Glasgow, Skotland, UK.  3.
Gannon, A. J., T. W. Von Backström, Solar Chimney Turbine Performance, ASME Journal of Solar Energy Engineering, 125, 2003, p101106.  4.
Haaf, W., Solar Chimneys, Part II: Preliminary Test Results from the Manzanares Pilot Plant, Int. J. Sol. Energy, 2 (1984), pp. 141–161.  5.
Jaluria, Y., Natural Convection, Heat and Mass Transfer, Oxford, New York, Pergamon Press, 1980.  6.
Mancini, Th. R., SolarElectric Dish Stirling System Development, SAND972924C, 1998.  7.
Mueller R.W., Dagestad K.F., Ineichen P., Schroedter M., Cros S., Dumortier D., Kuhlemann R., Olseth J.A., Piernavieja G., Reise C., Wald L., Heinnemann D. (2004), Rethinking satellite based solar irradiance modelling  The SOLIS clear sky module. Remote Sensing of Environment, 91, 160174.  8.
Raiss, W., Heiz und Klimatechnik, Springer, Berlin, vol. 1, pp. 180188, 1970.  9.
Rugescu, R. D., ThermischeTurbomaschinen, ISBN 9733018465, E. DP. Bucharest 2005.  10.
Rugescu, R. D., T. G. Chiciudean, A. C. Toma, F. Tache, Thermal Draught Driver Concept and Theory as a Tool for Advanced InfraTurbulence Aerodynamics, in DAAAM Scientific Book 2005, ISBN 3901509437 (Ed. B. Katalinic), DAAAM International Viena, 2005.  11.
Rugescu, R. D., D. A. Tsahalis, and V. Silivestru (2006), SolarGravity Power Plant Modeling Uncertainties, WSEAS Transactions on Power Systems, ISSN 17905060, 10(2006), pp.17131720.  12.
Rugescu, R. D. (2008), “Technology of CFD in space engines and solargravitational power plants”, International Journal on Energy Technology and Policies, ISSN 14728923, Online ISSN 1741508X, 6, 12(2008), p.124142.  13.
Rugescu, R. D., Silivestru V., and Ionescu M. D. (2008),Technology of Thermal Receivers and Heat Transfer for SEATTLER towers, Chapter 59 in the DAAAM International Scientific Book 2008, pp. 697742, B. Katalinic (Ed.), published by DAAAM International Vienna, ISBN 3901509690, ISSN 17269687, Vienna, Austria.  14.
Rugescu, D. R. and R. D. Rugescu (2009), Chapter 1 “Potential of the solar energy on Mars”, pp. 0124, in Rugescu, D. R. (editor), “Solar Energy”, ISBN 9789533070520, Ed. INTECH, Rjeka, Croatia, 2010, 462 p.  15.
Rugescu, R. D., R. B. Cathcart, Dragos R. Rugescu, and S. P. Vataman (2010), Electricity and Freshwater Macroproject in the Arid African Landscape of Djibouti, ARGEOC3 Conference, Exploring and harnessing the renewable and promising geothermal energy, 2225 November, Djibouti.  16.
Rugescu, R. D., R. B. Cathcart, and Dragos R. Rugescu, Electricity and Freshwater Macroproject in an Arid, Ochre African Landscape: Lac Assal (Djibouti), Eighth International Conference on Structural Dynamics, Proceedings EURODYN2011, Leuven, Belgium, 46 July 2011 (pp NA).  17.
Rugescu, R. D., A. Bogoi, and R. Cirligeanu (2011), Intricacy of the TRANSIT manifold concept paidoff by computational accuracy, Proceedings of ICMERA2011 Conference, ISBN 9789810704209, ISSN 2010460X, Bucharest 2022 October, Romania, pp. 6165.  18.
Rugescu, R. D. (Ed.), Solar Power, ISBN 9789535100140, INTECH Press, Rjeka, Croatia, Hard cover, 378 pages, Feb. 2012.  19.
Rugescu, R. D., and C. Dumitrache (2012), Solar Draught Tower Simulation and Proof by NumericalAnalytical Methods, Paper V&V20126191, ASME 2012 Verification and Validation Symposium, Las Vegas, May 24, 2012, p. 73.  20.
Rugescu, R. D., Demos T. Tsahalis, CiprianDumitrache (2012), Extremely High Draught Tower Efficiency Proof by Conservation of Energy, Invited plenary lecture, Proceedings of the 5th ICSCCE 2012 From Scientific Computing to Computational Engineering, 47th July 2012, Athens, Greece, pp. 243248.  21.
Scharmer, K., Greif, J. (2000), The European Solar Radiation Atlas, Presses de l'Ecole des Mines, Paris, France.  22.
Schiel, W., T. Keck, J. Kern, and A. Schweitzer, Long Term Testing of Three 9 kW DishStirling Systems, ASME International Solar Energy Conference, San Francisco, CA, USA, March 1994.  23.
SchleichBergermann und Partners (2005), EuroDish System Description.  24.
Tache, F., Rugescu, R. D., Slavu, B., Chiciudean, T. G., Toma, A. C., and Galan, V. (2006), Experimental demonstrator of the draught driver for infraturbulence aerodynamics, Proceedings of the 17th International DAAAM Symposium, "Intelligent Manufacturing & Automation: Focus on Mechatronics & Robotics", Vienna, 811th November 2006.  25.
Unger, J., Konvektionsströmungen, B. G. Teubner, ISBN 3519030330, Stuttgart, 1988.  26.
Energy Information Administration, Form EIA63B, “Annual Photovoltaic Module/Cell Manufacturers Survey.”  27.
Energy Information Administration, Office of Coal, Nuclear, Electric and Alternate Fuels, Renewable Energy Annual 1996, U.S. Department of Energy, Washington, DC 20585, April 1997.  28.
University of Massachusetts (2004), Lowell Photovoltaic Program, International Solar Irradiation Database, Version 1.0, Monthly solar irradiation.