Implementation of CIE General Sky Model Approach in Ukraine and Effects on Room Illuminance Mode

Similar documents
International Journal of Engineering Research and Development e

Yoshida-honmati, Sakyou-ku, Kyoto, Dept.of Urban and Environmental Eng., Grad.of Eng., Kyoto Univ.

Potential daylight resources between tropical and temperate cities a case study of Ho Chi Minh City and Moscow

Simulation of daylight and artificial lighting integration and energy savings

COMPUTER PROGRAM FOR THE ANGLES DESCRIBING THE SUN S APPARENT MOVEMENT IN THE SKY

Classification of CIE standard skies using probabilistic neural networks

Chapter 1 Solar Radiation

A Critical Review of the Preetham Skylight Model

Observer-Sun Angles. ), Solar altitude angle (α s. ) and solar azimuth angle (γ s )). θ z. = 90 o α s

Chapter 2 Available Solar Radiation

ESTIMATION OF HORIZONTAL ILLUMINANCE BY MEASURED SKY LUMINANCE DISTRIBUTION

HORIZONTAL AND VERTICAL ILLUMINANCE/IRRADIANCE FROM THE IDMP STATION IN GENEVA

Perez All-Weather Sky Model Analysis

DAYLIGHT FACTOR DISTRIBUTION IN INTERIOR BASED UPON THE MEAN SKY

Daylighting on the working plane in oriented attic rooms under overcast and clear sky

SATELLITE DERIVED AND GROUND TRUTH HORIZONTAL AND VERTICAL ILLUMINANCE/IRRADIANCE FROM THE IDMP STATIONS AT GENEVA AND GÄVLE

THE AVERAGE TOTAL DAYLIGHT FACTOR

Calculating equation coefficients

1/55. Solar energy. solar radiation definitions incident solar energy

Analysis of Energy Savings and Visual Comfort Produced by the Proper Use of Windows

ft ILLUMINANCE! DAY I -i I i I

Chapter 5. Daylighting

Direct Normal Radiation from Global Radiation for Indian Stations

AVAILABLE SOLAR RADIATION THEORETICAL BACKGROUND

Motion of the Sun. View Comments

THE SOLAR RESOURCE: PART I MINES ParisTech Center Observation, Impacts, Energy (Tel.: +33 (0) )

1/71 AES-LECTURE 1. Solar energy. solar radiation definitions incident solar energy

OPTIMISING THE TEMPORAL AVERAGING PERIOD OF POINT SURFACE SOLAR RESOURCE MEASUREMENTS FOR CORRELATION WITH AREAL SATELLITE ESTIMATES

Hottel s Clear Day Model for a typical arid city - Jeddah

Regularities of Angular Distribution of Near-Horizon Sky Brightness in the Cloudless Atmosphere

Sunlight and its Properties II. EE 446/646 Y. Baghzouz

many partly cloudy skies (approximately 4-7 oktas).

Chapter Seven. Solar Energy

How to evaluate daylight. Initiated by the VELUX Group

Page 1. Name:

Chapter 6. Solar Geometry. Contents

Stochastic modeling of Extinction coefficients for solar power applications

Introduction to Photovoltaics

What is it good for? RT is a key part of remote sensing and climate modeling.

Standard Table for Reference Solar Spectral Distributions: Direct and Diffuse on 20 Tilted and Vertical Surfaces 1

Vertical Illuminance Measurement for Clear Skies in Tehran

Evaluation of ultraviolet radiation protection of a membrane structure using a UV Shade Chart

On Skylight and Aerial Perspective. A.J. Preetham ATI Research

Earth s Orbit. Sun Earth Relationships Ridha Hamidi, Ph.D. ESCI-61 Introduction to Photovoltaic Technology

DAYLIGHTING MEASURES FOR ECOHOUSING

XI. DIFFUSE GLOBAL CORRELATIONS: SEASONAL VARIATIONS

Estimation of Hourly Global Solar Radiation for Composite Climate

Daylight climatology in Athens, Greece: types of diurnal variation of illuminance levels

CHARACTERIZATION OF SKY CONDITIONS USING CLEARNESS INDEX AND RELATIVE SUNSHINE DURATION FOR ISEYIN, NIGERIA

Solar Time, Angles, and Irradiance Calculator: User Manual

Sunlight and its Properties Part I. EE 446/646 Y. Baghzouz

Mathieu Hébert, Thierry Lépine

Solar Insolation and Earth Radiation Budget Measurements

LONG -TERM OBSERVATIONS OF CLOUDNESS IN HURBANOVO

Daylighting design with climate. Francesco Anselmo

IMPROVED MODEL FOR FORECASTING GLOBAL SOLAR IRRADIANCE DURING SUNNY AND CLOUDY DAYS. Bogdan-Gabriel Burduhos, Mircea Neagoe *

Lecture 3: Atmospheric Radiative Transfer and Climate

C) the seasonal changes in constellations viewed in the night sky D) The duration of insolation will increase and the temperature will increase.

TOTAL COLUMN OZONE AND SOLAR UV-B ERYTHEMAL IRRADIANCE OVER KISHINEV, MOLDOVA

GLOBAL AND DIFFUSE ILLUMINANCE DATA IN MAKASSAR-INDONESIA SEN 166

Solar radiation in Onitsha: A correlation with average temperature

ME 430 Fundamentals of Solar Energy Conversion for heating and Cooling Applications

ESTIMATION OF HOURLY GLOBAL SOLAR RADIATION INCIDENT ON INCLINED SURFACES IN IRAQ AT DIFFERENT SKY CONDITION

Table 1-2. TMY3 data header (line 2) 1-68 Data field name and units (abbreviation or mnemonic)

COMPARISON OF GUNN BELLANI RADIOMETER DATA WITH GLOBAL SOLAR RADIATION SENSOR (PYRANOMETER CM6B) Author. Mungai Peter N.

EGEE 437: HWK #2. Brownson. Yaqdaan Alkayyoomi, Mackenzie Ailes, Nicholas Minutillo, Sheel Vora. Group 17. Due on Thursday, Feb.

Science & Technologies MODELING ULTRAVIOLET RADIATION FOR BURGAS MUNICIPALITY

HIGH PERFORMANCE MSG SATELLITE MODEL FOR OPERATIONAL SOLAR ENERGY APPLICATIONS

Exercise 6. Solar Panel Orientation EXERCISE OBJECTIVE DISCUSSION OUTLINE. Introduction to the importance of solar panel orientation DISCUSSION

EAS 535 Laboratory Exercise Solar Radiation

FUNDAMENTALS OF REMOTE SENSING FOR RISKS ASSESSMENT. 1. Introduction

Section 22. Radiative Transfer

Radiometry and Photometry

Solar Radiation Measurements and Model Calculations at Inclined Surfaces

Spectrum of Radiation. Importance of Radiation Transfer. Radiation Intensity and Wavelength. Lecture 3: Atmospheric Radiative Transfer and Climate

Satellite-to-Irradiance Modeling A New Version of the SUNY Model

Sky brightness measurements for different environmental conditions by the example of St. Petersburg

Introducing NREL s Gridded National Solar Radiation Data Base (NSRDB)

Radiation in the atmosphere

Numerical and experimental analysis on Double Light Pipe, a new system for daylight distribution in interior spaces.

Section 10. Radiative Transfer

OPAC 101 Introduction to Optics

Interaction between building services and environment.

HIGH PERFORMANCE MSG SATELLITE MODEL FOR OPERATIONAL SOLAR ENERGY APPLICATIONS

Spatiotemporal Analysis of Solar Radiation for Sustainable Research in the Presence of Uncertain Measurements

DYNAMIC DAYLIGHT AND SOLAR CONTROL IN TROPICAL CLIMATE

Determination of Photometric Properties of Materials for Energy Purposes Through the Experimental Study of a Two-Axis Goniophotometer

Estimation of solar radiation at Uturu, Nigeria

White Paper Luminance & Illuminance. Brief explanation of photometry for the application of tunnel lighting control

Donghyun Seo. 1. Introduction

Integrating the Solar Spectrum

SURFACE ORIENTATIONS AND ENERGY POLICY FOR SOLAR MODULE APPLICATIONS IN DHAKA, BANGLADESH

Purdue University Meteorological Tool (PUMET)

CHAPTER 3. The sun and the seasons. Locating the position of the sun

LECTURE 13: RUE (Radiation Use Efficiency)

Solar Energy Basics II - Estimation of Solar Radiation to a Solar Panel

Evaluation of cloudiness/haziness factor for composite climate

Research Article A Satellite-Based Sky Luminance Model for the Tropics

COST action 718 METEOROLOGICAL APPLICATIONS FOR AGRICULTURE. Spatialisation of Solar Radiation - draft report on possibilities and limitations

Simultaneous spectral albedo measurements near the Atmospheric Radiation Measurement Southern Great Plains (ARM SGP) central facility

Transcription:

, pp. 57-70 http://dx.doi.org/10.14257/ijsh.2016.10.1.07 Implementation of CIE General Sky Model Approach in Ukraine and Effects on Room Illuminance Mode Oleg Sergeychuk 1* and Dmitriy Radomtsev 2 Department of architectural constructions Kiev national university of construction and architecture, Kiev, Ukraine 1 ovsergeich@mail.ru, 2 graphichart65@gmail.com Abstract New standard CIE S 011/E:2003 CIE standard general sky was accepted in Ukraine in 2010. With its adoption the question of its appropriate usage with national documents became an important event. Feasible way of standard s utilization with Ukrainian daylighting standards and effects of its implementation are given in the paper. Keywords: Building climatology; CIE general sky; systems of room daylighting 1. Introduction In 2003 an international standard CIE S 011/E:2003 [1] was accepted, which provide 15 main types of standard skies for daylight calculation. Document is basing on the comprehensive study of R.Kittler, S.Darula and R.Perez [2]. This primary standard was identical translated and applied in Ukraine in 2010 [3]. Mentioned documents describing only mathematical modeling of daylight spatial distribution but not providing method of usage in architectural practice for estimation and designing systems of room daylighting. Ukrainian designers have some questions: which sky type is required under designing systems of room daylighting in specific geographical locality? how to take into consideration fact that standard sky types are varying in the course of year? From the moment of new standard enactment a numerous foreign investigations of standard sky types gradations for different specific localities have been taken. PR Tregenza proposed the method for defining of sky types occurrence [4]. It is basing on the best fit comparison of the observed sky spatial luminance data retrieving by the 145 patches camera and summary horizontal illuminance with the results of luminance modeling under CIE Standard General sky. Basing on the provided method W.Li, L.Tang, M.Lee and T.Muneer take study in which they perform measurements of sky luminance in Hong Kong from 1999 to 2005. On the minimal RMSE approach using neural networks they defined 3 standardized sky types with the greatest repetition frequency in the course of year [5]. Afterwards other group of researchers made more accurate derivation and proposed a set of 4 general skies for Hong Kong [6]. In addition to existing efforts in [2] R.Kittler and S.Darula proposed an approach for definition adjustment of sky types by vertical illuminance [7]. Following M. Fabian made exploration of types occurrence in Bratislava [8]. Also well-known investigations of general sky gradation in Chile [9] and Singapore [10]. Listed papers is not enough responding to the formulated questions. * Corresponding Author ISSN: 1975-4094 IJSH Copyright c 2016 SERSC

In 2010 in Ukraine standard DSTU-N B V.1.1-27:2010 «Building climatology» was accepted [11], in which presented behaviors of cloudiness, irradiance and climatology indicators for Ukraine regional centers. Quantities of zenith luminance can be derive from the Satellight program data [12]. The main aim of paper is the development of offer for utilization mentioned above standards for designing of room daylighting systems for Ukrainian localities basing on the existing climatology data [11, 12]. Effects of CIE general sky usage on reference room s illuminance mode are given in the paper. 2. CIE General Sky Application in Ukraine 2.1. Geometric Apparatus of General Sky Model In [3] all sky types deriving by 3 groups of cloudiness: overcast, intermediate and clear. Each of sky conditions is describing mathematical by single formula of spatial luminance distribution. Geometric apparatus consist of following components (Figure 1): sky dome which is representing by hemisphere with relative radius in spherical coordinate system; sun point on the hemisphere surface which position is defining by azimuth and elevation in specified time and day; arbitrary sky element calculation with specified azimuth and elevation for which relative and absolute luminance point is determining. The main parameters of system are two functions that combine relative position of elements: φ(z) luminance gradation function that incorporate vertical position of sky element and zenith describing vertical distribution form horizon to zenith; f(χ) relative scattering indicatrix that combine horizontal position of element and sun describing horizontal pattern of luminance distribution. Figure 1. Geometric Apparatus of General Sky Model Each of sky types is deriving by six light-climatic parameters and individual set of coefficients by meaning of two mentioned functions [3, Table 1] that is placing into formula of relative luminance distribution calculation: 58 Copyright c 2016 SERSC

, (1) where L γ luminance of arbitrary sky element, cd/m 2 ; L z zenith luminance, cd/m 2 ; γ angle of elevation of arbitrary sky element, rad; Z angular distance between a sky element and the zenith, rad; χ shortest angular distance between sky element and sun, rad; Z s angular distance between the sun and zenith, rad; a,b,c,d,e coefficients that deriving group of gradation and scattering indicatrix of sky type [3, table 1]. 2.2. Defining of General Sky Types for the Territory of Ukraine On the assumption of [2] basic parameters that describing sky types are: sky dome group by cloudiness; D v /E v ratio of diffuse illuminance to extraterrestrial horizontal illuminance; T v luminous turbidity factor; L z /D v ratio of absolute zenith luminance to horizontal diffuse illuminance. 1. By cloudiness behavior sky dome can be derive by the repetition frequency of days with certain cloudiness value basing on the available Satellight data [12] with subsequent attitude to specific group. Table 1. Kiev CIE General Sky Types Month І ІІ ІІІ IV V VI VII VIII IX X XI XII Cloudy 22 28 25 11 15 11 11 23 30 27 39 36 24 Intermediate Cloudless 55 51 45 23 29 37 31 31 47 38 43 43 36 23 21 30 66 56 52 46 46 23 35 18 21 40 Group II II II III III III III III II II II II II 2. For the reason that solar data in DSTU-N «Building climatology» [Table 8, 11] are radiant quantity for the comparison it has to be convert into luminous quantities. Generally luminous flux of lighting source with continuous spectrum is: 780 683 e ( ) V ( ) d, (2) 380 where Ф luminous flux, Lm; Φ eλ (λ) spectral density of radiation flux; V(λ) spectral luminous efficiency function [13]. Thereby spectral density of diffuse and direct radiation flux for horizontal surface has been got basing on the climatic data [11] for 12 00 of June, 15 and December, 15 for Kiev and other 24 regional centers. Calculations are provided with usage of bundled software SMARTS [14-16]. Atmosphere profile, volume fractions of atmosphere elements, model of atmospheric aerosol and optical thickness received from international documents and satellite data [17-19]. Basing on the diagrams of solar radiation s spectral distribution tables of the spectral density for the June, 15 and December, 15 has been got. Values of the relative spectral density for other days ca be obtain with the interpolation. Thereby Copyright c 2016 SERSC 59

using formula (2) with coefficients of the relative spectral density and knowing quantities of spectral luminous efficiency function (ISO 11664:2007(E) [13]) for the CIE 1931 standard observer, the values of the direct and diffuse illuminance were obtain. Figure 2. Solar Radiation s Spectral Distribution, Kiev, Ukraine Figure 3. Solar Radiation s Relative Spectral Distribution in Visible Range For the determining of D v /E v ratio diffuse irradiation data D e_mon from DSTU-N «Building climatology» need to be lead to instant value D e and then to instant luminous quantity D v basing on formula (2). D e mon D e n t _ (3) where D e instant horizontal diffuse radiation in real cloudiness conditions, W/m 2 ; D e _ mon summary monthly horizontal diffuse radiation in real 60 Copyright c 2016 SERSC

cloudiness conditions, MJ/m 2 [table 8, 11]; n day s amount in specific month; t duration of daylight hours, sec. As mentioned above E v horizontal illuminance reaching outer limit of Earth s atmosphere and can be define by: E v E v 0 sin, (4) γ s where E v0 luminous solar constant, 133800, lx [20]; γ s solar altitude, rad. Table 2. Defining of D v / E v Data 15.01 15.02 15.03 15.04 15.05 15.06 15.07 15.08 15.09 15.10 15.11 15.12 D e_mon, MJ/m 2 60 96 170 215 280 290 288 230 168 110 58 44 167.42 D e, W/m 2 48.90 86.50 117.1 132.6 147.5 152.4 156.8 137.3 119.6 89.53 48.85 43.78 106.75 E v, 38.24 47.79 67.92 89.36 107.5 117.4 119.5 114.3 101.4 81.17 60.78 43.20 82.426 D v, 6.25 10.33 13.57 15.09 16.60 17.02 17.65 15.36 13.85 10.69 6.25 5.94 12.386 D v/e v 0.137 0.102 0.131 0.136 0.136 0.147 0.145 0.154 0.168 0.199 0.216 0.163 0.153 3. T v luminous turbidity factor that in this mathematical model can be determine by following formulas: Tv ln( Pv / E v ), (5) a v m where P v instant direct horizontal illuminance, lx; a v luminous extinction coefficient; m optical air mass. 1 a v 9.9 0. 043, (6) m sin γ s 0.050572 1 m (7) (γ s 6.07995 1.6364 ) Instant direct horizontal illuminance can be calculated similarly to diffuse value. Table 3. Defining of T v Data 15.01 15.02 15.03 15.04 15.05 15.06 15.07 15.08 15.09 15.10 15.11 15.12 D e_mon, MJ/m 2 25 46 100 170 285 300 300 265 170 80 23 15 148.25 P e, W/m 2 20.3 41.4 68.2 104 150 157 166 158 121 65.1 19.3 14.2 90.60 P v, 8.28 14.6 19.8 22.4 24.9 25.8 26.5 23.2 20.2 15.1 8.27 2.34 17.65 m 3.07 2.19 1.64 1.31 1.16 1.11 1.13 1.23 1.49 1.94 2.78 3.46 1.88 a v 0.09 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.09 0.09 0.100 T v 5.00 5.39 7.45 10.4 12.4 13.4 13.1 12.7 10.7 8.52 7.18 8.45 9.58 Copyright c 2016 SERSC 61

4. For the derivation of L z /D v ratio zenith luminance satellite data [12] can be taken. In this work is used information of the 1996 year as the closest data to the basic CIE S 011/E 2003 measurements. Absent information in east regional centers was defined with the Satellight approach in use of ARSC-CIE sky model [21], according which zenith luminance is: / z v i i s i s i L D [ a c c o s Z c e x p ( 3 Z ) d ], (8) where a i, c i, c / i, d i coefficients that describe light-climatic parameters by e sky s clearness, [table 4, 21] which can be determine by: e [( D v E v _ n ) / D v 3 3 kz s ] /[1 kz s ], (9) E v_n instant normal direct radiation in real sky condition, Lx; Δ sky s brightness given by:, (10) D v m / E v 0 m relative optical air mass. Table 4. Defining of T v Data 15.01 15.02 15.03 15.04 15.05 15.06 15.07 15.08 15.09 15.10 15.11 15.12 L z, kcd/ m 2 3.21 4.19 6.43 7.18 7.09 11.18 12.10 8.69 7.31 4.50 3.78 2.72 6.535 L z /D v 0.392 0.330 0.268 0.212 0.168 0.143 0.144 0.175 0.245 0.359 0.479 0.432 0.279 Basing on information about cloudiness group, D v /E v, T v, L z /D v ratios monthly CIE general sky types for Kiev have been defined (Table 5). Table 5. Defining of CIE General Sky Types for Kiev, Ukraine Month І ІІ ІІІ IV V VI VII VIII IX X XI XII Sky type ІV.2 9 ІV.3 10 ІV.4 11 ІV.4 11 ІV.3 10 ІІІ.2 6 For the comparison of the defining results we can use the horizontal diffuse illumination data from the Ukrainian climatology standard [11], absolute zenith luminance quantities and calculated values of horizontal illuminance that have been obtained during modeling of the sky luminance surface according to the formula 1. Surface modeling is processing with the bundle software Diffused solar radiation. Illuminance can be calculated as the integral behavior of light field by formula: E L (, ) c o s d, (11) 62 Copyright c 2016 SERSC

where L (, ) - luminance of arbitrary sky element in specified direction by angles and, cd/m 2 ; cosβ function of radiation value; ω solid angle in which radiating surface is defining, ster; Table 6. Comparison of Meteorological and CIE General Sky Calculated Horizontal Diffused Illuminance Data for Kiev, Ukraine Data 15.01 15.02 15.03 15.04 15.05 15.06 15.07 15.08 15.09 15.10 15.11 15.12 D v, Sky type 6.25 10.3 13.5 15.1 16.6 17.0 17.6 15.3 13.8 10.6 6.25 5.94 12.386 ІV.2 9 ІV.3 10 ІV.4 11 L z, kcd/m 2 1.66 2.22 3.31 3.53 3.61 3.74 4.63 3.76 3.14 2.63 1.84 1.59 2.97 E d_calc, 6.48 9.81 13.9 15.2 15.9 16.3 19.5 16.8 12.7 10.2 6.82 5.55 12.46 E,% 3.61 5.05 2.87 1.07 3.95 3.68 10.5 9.72 7.66 4.39 8.88 6.66 5.67 ІV.4 11 ІV.3 10 ІІІ.2 6 On basis of proposed method was analyzed the remaining part of Ukraine (Table 7) and proposed a map of light-climatic zoning by D v (Figure 4). Figure 4. Map of Light-climatic Zoning of the Territory of Ukraine by D v Copyright c 2016 SERSC 63

Table 7. Source Data, Intermediate Calculations and Defining of CIE General Sky Types for the Territory Of Ukraine Regional center D e, W/m 2 D v_year_mid, Lx L z_year_mid cd/m 2 L z/d v D v_calc, Lx E,% 1 Vinnytsia 114,33 18221 4450 0,24 18251 0,16 2 Dnipropetrovsk 116,01 18461 4456 0,24 18668 1,11 3 Donetsk 113,22 18022 4430 0,24 18496 2,56 4 Zhytomir 115,44 18252 4483 0,24 18386 0,72 5 Zaporizhia 114,98 18230 4461 0,24 18625 2,12 6 Ivano-Frankivsk 112,42 18223 4460 0,24 18615 2,10 7 Kiev 117,67 18166 4315 0,23 18242 0,41 8 Kirovohrad 118,03 18695 4434 0,23 18576 0,64 9 Luhansk 114,04 18104 4418 0,24 18512 2,20 10 Lutsk 114,09 18035 4097 0,23 17635 2,26 11 Lviv 112,02 18206 4395 0,24 18517 1,68 12 Mykolaiv 119,17 18812 4671 0,24 19403 3,04 13 Odesa 120,03 19026 4617 0,23 19150 0,64 14 Poltava 115,49 18303 4378 0,24 18433 0,70 15 Rivne 114,15 18035 4175 0,23 17657 2,14 16 Simferopol 121,78 19330 4742 0,24 19490 0,82 17 Sumy 113,21 18025 4090 0,23 17297 4,20 18 Ternopil 111,90 18254 4461 0,24 18780 2,80 19 Uzhhorod 115,43 18226 4416 0,24 18504 1,50 20 Kharkiv 114,92 18222 4323 0,23 18254 0,17 21 Kherson 120,14 19014 4672 0,24 19387 1,92 22 Khmelnytskyi 111,78 18217 4445 0,24 18704 2,60 23 Cherkasy 117,17 18246 4432 0,24 18652 2,17 24 Chernivtsi 113,63 18771 4514 0,24 18901 0,68 25 Chernihiv 117,16 17734 4085 0,23 17347 2,23 sky type 2.3. Defining of Calculation Sky Types for the Territory of Ukraine As seen from the Table 5 annual sky condition of Ukraine response to the single sky type -. But over a year period light-climatic behaviors have a wide general skies range [22] that is why more appropriate approach is a data calculation of luminance distribution by per month defining of sky type with subsequent derivation of annual year data of relative luminance distribution. Evaluation of luminance distribution have been held for 15-th of every month. For the defining of daily annual distribution for each whole hour modeling luminance data by the appointed CIE general sky type (Figure 5, a) and on basis of received information annual daily/monthly values are defining (Figure 5, b). Modeling occur in software environment «MatLAB» with utilization of bundled software «Atmospheric radiation» [23]. Figure 5. Luminance Distribution Surfaces for Kiev. a Hourly Distribution, March,15; b Monthly Distribution, March 64 Copyright c 2016 SERSC

Definitive annual year data of spatial luminance distribution deriving as value from all-month annual monthly information (Figure 6). As the result we obtain array of annual year data and table of relative quantities of luminance distribution arbitrary sky elements in spherical coordinate system with 10 step. Figure 6. Surface of Year Spatial Luminance Distribution for Kiev Table 8. Year Coefficients of Luminance Distribution for Kiev Elevation/ Azimuth 0 10 20 30 40 50 60 70 80 90 0 0.44 0.43 0.37 0.32 0.28 0.26 0.25 0.25 0.26 0.28 10 0.44 0.42 0.37 0.31 0.28 0.26 0.26 0.25 0.26 0.28 20 0.45 0.43 0.38 0.32 0.29 0.26 0.26 0.25 0.26 0.28 30 0.48 0.46 0.40 0.34 0.30 0.27 0.26 0.26 0.26 0.28 340 0.51 0.49 0.42 0.35 0.31 0.28 0.26 0.26 0.26 0.28 350 0.46 0.45 0.38 0.33 0.33 0.27 0.26 0.25 0.26 0.28 360 0.44 0.43 0.37 0.32 0.32 0.26 0.25 0.25 0.26 0.28 Application in national building standard can be obtain by the replacement of calculation formula of q і, coefficient that considering irregular luminance of arbitrary sky element, with the obtained table of annual year relative luminance distribution for every regional centre of Ukraine. At the meantime coefficient is calculating for the standard overcast CIE sky: 3 q (1 2 sin θ ), i (12) 7 where θ angular height of centre of the arbitrary sky element relatively to the calculation point. Intermediate quantities can be obtain using interpolation. Subsequent calculations execute according to the existing method of daylight computation. More detailed information is given in [24]. 2.4. Influence of CIE General Sky Model on Room Luminance Conditions Analysis of house illuminance mode is basing on comparison and analytic review of illuminance data of working plane in conditions of CIE standard overcast sky and determined CIE general sky types. Study provide for Kiev (50 34 N, 30 31 E) over whole year for 09 00 AM and 15 th date of each month. For analysis typical living room Copyright c 2016 SERSC 65

with dimensions 3.0x6.0m is select. Window size is 1.7x2.1m with double glazing and transmitting efficiency 78%. Calculation points are choosing according to the national sanitation standard with step 1 meter and height 0.8m above floor. Study results of month data for CIE standard overcast sky (Figure 7, Table 9) and determined CIE general sky types (Figure 8, Table 10) are shown below. Figure 7. Illuminance Mode of Reference Room from January to December in Conditions of CIE Standard Overcast Sky Table 9. Illuminance Mode of Reference Room December in Conditions of CIE Standard Overcast Sky for Kiev, Ukraine Data 15.01 15.02 15.03 15.04 15.05 15.06 15.07 15.08 15.09 15.10 15.11 15.12 Sky type CIE standard overcast sky E d_calc, E d_ext, 42.3 75.7 121 174 212 293 232 211 171 119 68.7 40.7 155.5 1303 2355 3763 5382 6600 7121 7102 6479 5238 3655 2136 1232 4363 DF,% 3.24 3.24 3.23 3.24 3.22 4.12 3.27 3.26 3.28 3.27 3.22 3.30 3.35 66 Copyright c 2016 SERSC

Figure 8. Illuminance Mode of Reference Room from January to December in Conditions of Determined CIE General Sky Types Table 10. Illuminance Mode of Reference Room December in Conditions of Determined CIE General Sky Types for Kiev, Ukraine Data 15.01 15.02 15.03 15.04 15.05 15.06 15.07 15.08 15.09 15.10 15.11 15.12 Sky type E d_calc, E d_ext, ІV.2 9 ІV.3 10 ІV.4 11 ІV.4 11 ІV.3 10 ІІІ.2 6 78.7 150 212 297 313 291 293 316 241 191 122 65.1 214.5 1304 2331 3768 5418 6584 7161 7144 6526 5303 3656 2116 1234 4378 DF,% 6.03 6.43 5.63 5.48 4.76 4.07 4.11 4.84 4.56 5.23 5.78 5.27 5.18 As seen from the mentioned above tables usage of CIE general sky model have a great influence on room illuminance mode with identical window s dimensions in both cases annual horizontal illuminance on working plane and daylight factor is 31% and respectively 35% higher than with CIE standard overcast sky model. Detailed analysis of both cases is presented in Table 11. Table 11. Comparison of Average Illuminance Mode Data Data E d_calc, E d_ext, CIE standard overcast sky CIE standard general sky E, % 155.5 214.5 31.38 4363 4378 0.34 DF,% 3.35 5.18 35.87 Copyright c 2016 SERSC 67

3. Conclusions and Future Prospect Research Adoption of new calculation model brings an opportunity for more accurate evaluation of daylighting over a period of a whole year. In case of annual year lighting behaviors rising relative to CIE overcast sky, calculation level of daylighting will increase what will bring for more rationally utilization of artificial lighting and increment of whole energy efficiency level. For the future research derivation of annual year luminance distribution tables and detailed investigation of new sky model impact onto the lighting behaviors of buildings is planning. Acknowledgment The authors thank to DSc R.Kittler for the valuable remarks during the research and organizing committee of ACE2015 conference for the possibility of publication. References [1] J. Clerk Maxwell, A Treatise on Electricity and Magnetism, vol. 2, (1892), pp. 68 73. [2] R. Kittler, R. Perez and S. Darula, A set of standard skies, (1998). [3] Kiev: State Enterprise, "Ukrainian Scientific-Research Center of Standardization", National standard of Ukraine, (2013), p. 7. [4] P. R. Tregenza, Analysing sky luminance scans to obtain frequency distributions of CIE Standard General Skies, Light. res. and tech., vol. 36, (2004), pp. 271-279. [5] W. Li, L. Tang, M. Lee and T. Muneer, Classification of CIE standard skies using probabilistic neural networks, International journal of climatology, vol. 30, (2009), pp. 305-315. [6] E. Ng, V. Cheng, A. Gadi, J. Mu, M. Lee and A. Gadi, Defining standard skies for Hong Kong, Building and Environment, vol. 42, (2007), pp. 866-876. [7] S. Darula, R. Kittler and L. Komar, Sky type determination using vertical illuminance, Przegląd elektrotechniczny, vol. 6, (2013), pp. 315-319. [8] M. Fabian and S. Darula, Occurrences of standard skies and luminous turbidity factor in Bratislava, Advanced Materials Research, vol. 855, (2014), pp. 259-265. [9] M.. B. Piderit, C. Cauwerts and M. Diaz, Definition of the CIE standard skies and application of the HDRI technique to characterize the spatial distribution of daylight in Chile, Journal of Construction, vol.13, (2014), pp. 22-30. [10] S. Wittkopf and L. Soon., Analysing sky luminance scans and predicting frequent sky patterns in Singapore, Light. res. and tech, vol. 39, (2007), pp. 31-51. [11] Кiev, Building climatology, National standard of Ukraine, (2011), p. 127. [12] Satellight, The European database of daylight and solar radiation, http://www.satel-light.com/core.htm. [13] CIE Colorimetry - Part 1: Standard Colorimetric Observers, CIE Central Bureau Vienna, (2007). [14] Standard Tables for Reference Solar Spectral Irradiances: Direct Normal and Hemispherical on 37 Tilted Surface, http://www.astm.org/standards/g173.htm. [15] Simple Model of the Atmospheric Radiative Transfer of Sunshine (SMARTS), http://www.nrel.gov/ rredc/smarts/. [16] Standard solar constant and zero air mass solar spectral irradiance tables, http://www.astm.org/ Standards/E490.htm. [17] AFGL Atmospheric constituent profiles (0-120km), Optical physics division, Air force geophysics laboratory, http://www.dtic.mil/cgi-bin/gettrdoc?ad=ada175173. [18] Guide to reference and standard atmosphere models, http://www.spacewx.com/docs/aiaa_g_003c _2010_9-10.pdf, (2010). [19] Aerosol optical thickness, NASA database, http://neo.sci.gsfc.nasa.gov/view.php?dtasetid= MODAL2_M_AER_OD&date=2014-10-15. [20] Guide to Recommended Practice of Daylight Measurement, CIE Tech. Rep., Commission Internationale de l Eclairage, Vienna: CIE Central Bureau Vienna, (1994). [21] R. Perez, P. Ineichen, R. Seals, J. Michalsky and R. Stewart, Modeling daylight availability and irradiance components from direct and global irradiance, Solar energy, vol. 44, (1990), pp. 271-289. [22] D. Radomtsev, Defining of sky types for Kiev basing on DSTU ISO 15469:2008, Energy efficiency in construction and architecture, vol. 7, (2015), pp. 248-261. [23] V. Bazhenov, P. Lizunov, O. Pidgorny, V. Plosky and O. Sergeychyk, Applied Software «Atmospheric Radiation» for an Energy Efficient Building, 14th International Conference on Computing in Civil and Building Engineering, (2012); Moscow. 68 Copyright c 2016 SERSC

[24] O. Sergeychuk and D. Radomtsev, Employment features of CIE S 011/E:2003 (ISO 15469:2004) CIE standard general sky under designing systems of room daylighting, 9 th International Conference of Future generation communication and networking, (2015); Jeju Island, South Korea. Authors Sergeychuk Oleg, he is a Doctor of technical sciences, professor of the Department of architectural constructions. Radomtsev Dmitriy, he is a Post-graduate student at Department of architectural constructions. Copyright c 2016 SERSC 69

70 Copyright c 2016 SERSC