Ciencias naturales e ingenierías

 

Simulation of direct solar radiation transmission in a building well

 

Simulación de la transmisión de radiación solar directa en un hueco aire-luz de edificio

 

Juan José Eliçabe Urriol^{1, 2} y Horacio A. Belluccia²

 

¹ Consejo de Investigaciones, Universidad Nacional de Rosario, Argentina.

² Facultad de Ciencias Exactas, Ingeniería y Agrimensura, Universidad Nacional de Rosario, Argentina.

 

J. J. Eliçabe Urriol. Dpto. de Física. Facultad de Ingeniería. Av. Pellegrini 250 - 2000 Rosario. E-Mail: elicabe@fceia.unr.edu.ar.

 

Recepción: 03-03-10
Aceptación: 16-08-10

 

Abstract

Propagation of direct solar radiation in a simple design light-well of quasi-specular inner surface is simulated using a stochastic model. Irradiances on inner windows and also on the well bottom were calculated for typical situations. The numerical results are of significance from the daylighting point of view.

Keywords: solar radiation, reflection, dayligthing, light-wells, stochastic, building.

 

Resumen

Con un modelo estocástico se simula la propagación de radiación solar directa en un hueco aire-luz de superficie interior cuasiespecular. Para situaciones típicas fueron calculadas las irradiancias sobre las ventanas internas y también la base del hueco. Los resultados numéricos, desde el punto de vista de la iluminación natural (de edificios), son significativos.

Palabras claves: radiación solar, reflexión, iluminación natural, hueco aire-luz, estocástico, edificio.

 

Introduction

When designing multi-storey buildings, it is sound practice to place windows facing a common vertical duct as a way of illuminating and ventilating central rooms.

These ordinary air/light wells supply a certain amount of solar radiation through the windows (those placed on the well walls). Wall surfaces are usually rough, and consequently, they generate a rather diffuse reflection of the incoming natural light.

Using light colors, it is possible to transmit radiation in depth; but at the sight of the well interior, the observer behind any inner window is exposed to unpleasant glare effects due to horizontal, or almost horizontal, reflected rays. To avoid this phenomenon one sort of high specular covering must be applied to the well surface. In this way the apparent solar altitude angle will remain almost the same after several inner reflections.

The daylight performance of a simple design light-well, but with a quasi-specular inner surface is here numerically simulated, for typical situations. This modeling field has been little explored up to now.

Similar problems related with the transmission of daylight by means of successive reflections on the inner surface of hollow ducts (or pipes) have received some attention since the eighties [1, 2, 3, 4]. In [5, 6], building light-wells of rather complex design has been experimentally studied.

A review of the scientific literature on the subject shows that in general, the basic purpose is to determine, either the radiation levels, or the photometric conditions at the end of the ducts, as effects of a known solar radiation input.

In the present contribution, instead of a unique outgoing flux (radiant or luminous), two quite different ones are computed: a) the flux which reaches the bottom of the well; and b), the lateral fluxes through the windows placed facing the well interior.

The model was built on the basis of the Snell´s law. Stochasticity is introduced into it by using a Monte Carlo scheme [7]. As outputs of the model, the inner irradiances relative to the outer solar irradiance on a horizontal plane were chosen.

To characterize the single light reflection phenomenon of a one direction incident radiation onto a low roughness reflecting surface, it was assumed that the radiation scattering is given by a narrow Gaussian distribution which is centered around the specular reflection direction.

 

THE MODEL

Solar radiation flux entering into a light well can be written as

where I is the vector defined by the direction, sense and intensity of direct solar radiation, k is the z-axe versor ( see figure 1 ) and S is the area of the well cross section.

The radiation incidence area (at the top of the well ) is divided into a given number of equal cells, being the incident energy flux on each cell: N ΔE/Δt, where ΔE is an elemental quantity of radiant energy ( or "photon"), Δt an arbitrary lapse time and N, the number of incident photons on the cell in a Δt time (here, "elemental" means that there are not any partition of ΔE in the model.

A vectorial algorithm simulates, step by step, the photons trajectories and photon-escape events through control surfaces (inner windows and the top and bottom of the well).

When a photon is intercepted by a reflecting surface the probability of absorption is evaluated by means of a uniform distribution. In the case of no absorption, specular reflection direction is calculated by the well known vectorial formula:

Where n is the normal at the reflecting surface, i_o is the versor which defines the sense and direction of the incident photon and ρ_o is the versor corresponding to the reflected photon.

To simulate the quasi-especular reflection phenomena it is used a Gaussian perturbation of the perfect specular reflection direction given by (2) [8].

The absolute energy gain for the i window , in a time Δt is

In a similar way the absolute energy at the bottom (neglecting any reflection from it) can be expressed as

In order to quantify the radiation energy distribution it is convenient to define:

a) the energy gain of i window as a fraction of the total incident direct radiation: G_wi= Φ_wi/Φ , where Φ_wi = Δ E_wi / Δt

b) the relation between the flux which reaches the well bottom Φ_b = Δ E_b / Δt, and the incident one:

Multiplying G_wi by the factor f = S/S_w (S_w : window area), the irradiance on i window is obtained as a fraction of the external irradiation on a horizontal plane , without considering diffuse radiation. This new quantity, I_wi = G_wi.f, supplies relevant information about the daylighting method in study.

Also it is calculated the global fraction of absorbed radiant energy by means of the expression α = Φ_a/Φ, where Φ_a is the total absorbed flux by the reflecting material.

The energy lost by retroreflection, that is , the radiant flux from the top of the well, is negligible for the cases here shown. So, from the energy conservation law, the following relation holds:

A reasonable validation of the computational code was obtained by performing a series of numerical experiments regarding energy conservation, convergence and fitting of exact solutions. It can be added that the solar intensity is automatically generated by a well known algorithm that process initial data of time (date, solar hour) and latitude.

 

Calculations Results

The dimensions and orientation of the light well that were chosen to perform the numerical examples are displayed in fig. 1. There are three stories and four windows in each storey.

It was selected an absortivity value of the reflecting material equal to 0.2. The 3σ value of the normal distribution fit a 10 degree angle of deviation respect of the specular reflection direction. For basic Montecarlo program running it was used a total photons number equal to 3. 10⁴.

Diffuse radiation was not taken into account.

Most of the simulations correspond to local latitude of 33 ^o S which, in the country, approximates the main cities latitudes (Table I). Due to the problem symmetry with respect to the North–South direction, all the results are also valid for the North hemisphere.

Additional results, for higher and lower latitudes are presented too (Tables II and III).

To calculate absolute irradiances, Tables I to III (II) can be used in connection with standard models which determine the atmosphere transmittance for direct solar radiation, as for instance the Hottel one [9]. Using this latter model it was built Fig. 2, where absolute irradiances values are shown. In all cases, 0 solar hour corresponds to solar noon.

An inspection of the irradiances curves in Fig. 2 shows that all of them decrease in a monotone way for solar time greater than 3hr. This common trend is due to the great number of inner reflections at low values of the solar altitude angle. But, when considering the complete solar time range, other factors must be taken into account in order to explain the mean irradiance evolution, as the variation of the solar azimuth angle. At solar noon, for instance, the east and west walls of the well are almost not irradiated; the solar irradiation on these walls increases from 0hr to 1hr, and consequently, the mean irradiance also increases.

 

Conclusions and Perspectives

At the present research stage the gaussian model seems a suitable approach to treat quasi-specular successive reflections. Experimental bidirectional spectral reflectivity [7] data for specific reflecting surfaces (of practical utility) could be used in the future to improve the calculation code.

Numerical values displayed in tables I to III (II) can be used with any solar direct radiation data to predict the corresponding internal irradiances values.

These irradiances, assuming that the solar spectrum is little modified by the reflections, can be considered as "solar irradiances", and then be used as inputs to determine illumination distributions into the rooms which are illuminated from the well; provided that all the complementary data are known: geometry, windows glass transmittance (spectral), reflectivities of the room surfaces, optical properties of diffusers (if needed), etc. Undoubtedly this is a complex task, but not so different to the one involved in the dayligthing calculations for a room which is illuminated through a window of the building facade. That is to say, in both situations, primary data of solar radiation can be used to obtain inner illuminances.

It can be expected significant solar radiation levels for either illumination or growing plant (photosynthetical active radiation) purposes.

In particular, to visualize the availability of natural illumination, a luminous efficiency of solar radiation equal to 100 lumen/watt [ 5 ] can be assumed.

As an example, it can be calculated the mean illuminances onto the windows at the lowest level , for 43 degree latitude, during the equinoxes. Using this latter efficiency figure in combination with Table II irradiances values and the Hottel model [ 9 ], it is found that mean illuminances are over 2200 lux between 4am to 4pm , with a maximum near 11000 lux at solar noon. On the contrary, when the well inner surface is smooth but of a dark color (absortivity around 0.5) , the model predicts that the corresponding illuminances values are much smaller.

It is worth to observe the considerable high values of the solar transmittance. This supports the idea of using the light- well as a sort of atrium, with a garden.

To appreciate the potential of this new passive daylighting strategy the results here discussed are enough; but, it is clear that more numerical-experimental research is required to gain a better insight on the subject.

Regarding the reasons to do additional research efforts, it must be said:

On one hand, the method can be implemented in existing buildings with an air-light well of similar characteristics to the one here studied .This can be done by attaching a reflecting film to the well surface (comments on specific techniques are beyond the scope of the this work).

On the other hand, the development of low cost high reflecting films would permit the method to be applied in deeper wells.

Also that, in any case, a reduction of the electrical energy consumption for diurnal illumination implies more sustainable buildings.

 

Acknowledgements

In memoriam Professor L.C.A. Fernandez.

 

References

1.- Courret G., Scartezzini J. L, Francioli D., Meyer J. J. Design and assessment of an anidolic light-duct. Energy and Buildings, V. 28, n. 1, pp. 79-99, (1998).         [ Links ]

2.- Oakley G., Riffat S. B. and Shao L. Daylight performance of lightpipes. Solar Energy, V. 69, n. 2, pp. 89-98, (2000).         [ Links ]

3.- Eliçabe Urriol J. J., Lara M. A. and Piacentini R. D. Passive solar light-ducts, Applied Energy, V. 28, pp. 95-195, (1987).         [ Links ]

4.- Zhang X and Muneer T. Mathematical model for the performance of light pipes. Lighting Research and Technology, V. 32, n. 3, pp. 141-146, (2000).         [ Links ]

5.- Zastrow A. y Vittwer V. Daylighting with fluorescent and highly reflective silver coated films. New applications for new materials. 3nd. International Symposium of Optical and Optoelectronic Applied Sciences and Engineering - Optical Materials Technology for Energy Efficiency and Solar Energy Conversion V, 14-18 April, Innsbruck, (1986).         [ Links ]

6.-Kristl Z, Krainer A. Light wells in residential buildings as a complementary daylight source. Solar Energy, V.65, pp. 197-206, (1999).         [ Links ]

7.- Siegel R. and Howell J. Thermal Radiation Heat Transfer. Mc Graw Hill, New York, (1965).         [ Links ]

8.- Eliçabe Urriol J. J. , Belluccia H. A., Vital M. V,, Navone H. D. Solar radiation transmission by means of rectangular cross section ducts with quasi-specular inner surface. Energías Renovables y Medio Ambiente, V. 16, pp. 15-20, (2005).         [ Links ]

9.- Hottel H.C. A simple model for estimating the transmittance of direct solar radiation through clear atmospheres. Solar Energy, V. 18, pp.129-134, (1976).         [ Links ]