The impact of wind on air temperature distribution in Athens and in Santorini

Similar documents
Estimation of Global Solar Radiation at Onitsha with Regression Analysis and Artificial Neural Network Models

Tests for the Ratio of Two Poisson Rates

Redesigning Urban Open Spaces Based on Bioclimatic Criteria: Two squares in Thessaloniki, Greece

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

Measuring Electron Work Function in Metal

Entropy ISSN

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d

New data structures to reduce data size and search time

Predict Global Earth Temperature using Linier Regression

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Vyacheslav Telnin. Search for New Numbers.

APPROPRIATE ATTENUATION MODEL FOR CHIANG MAI, THAILAND FROM FIELD MEASUREMENT TO MODEL EQUATION

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

The Properties of Stars

LECTURE 14. Dr. Teresa D. Golden University of North Texas Department of Chemistry

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

CBE 291b - Computation And Optimization For Engineers

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

A Brief Review on Akkar, Sandikkaya and Bommer (ASB13) GMPE

Recitation 3: More Applications of the Derivative

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

Estimation of the particle concentration in hydraulic liquid by the in-line automatic particle counter based on the CMOS image sensor

CAPACITORS AND DIELECTRICS

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

The steps of the hypothesis test

Section 6: Area, Volume, and Average Value

Acceptance Sampling by Attributes

NUMERICAL ANALISYS OF HUMAN THERMAL COMFORT INSIDE OCCUPIED SPACES. C. O. R. Negrão(*) C. O. Carvalho F o.(**) and C. Melo (**)

8.2 ESTIMATING DETENTION VOLUMES

1 Which of the following summarises the change in wave characteristics on going from infra-red to ultraviolet in the electromagnetic spectrum?

Energy Consideration

Thermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report

1. a) Describe the principle characteristics and uses of the following types of PV cell: Single Crystal Silicon Poly Crystal Silicon

List all of the possible rational roots of each equation. Then find all solutions (both real and imaginary) of the equation. 1.

New Expansion and Infinite Series

The Velocity Factor of an Insulated Two-Wire Transmission Line

Part I: Basic Concepts of Thermodynamics

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

8 Laplace s Method and Local Limit Theorems

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

15. Quantisation Noise and Nonuniform Quantisation

MATH SS124 Sec 39 Concepts summary with examples

1 Module for Year 10 Secondary School Student Logarithm

A Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System

Non-Linear & Logistic Regression

Math 1B, lecture 4: Error bounds for numerical methods

5.7 Improper Integrals

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009

Validation of APHRODITE Precipitation Data for Humid and Sub Humid Regions of Pakistan

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Temperature influence compensation in microbolometer detector for image quality enhancement

Design Data 1M. Highway Live Loads on Concrete Pipe

Probabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

DECAMETER RADIO EMISSION OF THE SUN: RECENT OBSERVATIONS


Simulation of Eclipsing Binary Star Systems. Abstract

Student Activity 3: Single Factor ANOVA

The heat budget of the atmosphere and the greenhouse effect

Applications of Bernoulli s theorem. Lecture - 7

Review of Calculus, cont d

Data Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading

Freely propagating jet

Effects of peripheral drilling moment on delamination using special drill bits

Synoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences?

P62 In-mine Geoelectric Methods for Detection of Tectonic Disturbances of Coal Seams

A New Statistic Feature of the Short-Time Amplitude Spectrum Values for Human s Unvoiced Pronunciation

CONTRIBUTION TO THE EXTENDED DYNAMIC PLANE SOURCE METHOD

Math 113 Exam 2 Practice

Shear and torsion interaction of hollow core slabs

Summary Information and Formulae MTH109 College Algebra

Strategy: Use the Gibbs phase rule (Equation 5.3). How many components are present?

Physics 9 Fall 2011 Homework 2 - Solutions Friday September 2, 2011

User's Guide of the Inter-Batch Physico-Chemical Variability Calculator

Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I

Lecture 1. Functional series. Pointwise and uniform convergence.

Math 8 Winter 2015 Applications of Integration

Name Solutions to Test 3 November 8, 2017

Math 116 Final Exam April 26, 2013

Comparison Procedures

Vorticity. curvature: shear: fluid elements moving in a straight line but at different speeds. t 1 t 2. ATM60, Shu-Hua Chen

MATH 144: Business Calculus Final Review

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE

INVESTIGATION ON THE MODEL OF VORTEX-INDUCED

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data

Module 2: Rate Law & Stoichiomtery (Chapter 3, Fogler)

Solutions to Supplementary Problems

Computation and Modeling of the Air-Sea Heat and Momentum Fluxes. By Dr. Nasser H. Zaker

We divide the interval [a, b] into subintervals of equal length x = b a n

Key for Chem 130 Second Exam

The Form of Hanging Slinky

Quantum Mechanics Qualifying Exam - August 2016 Notes and Instructions

Riemann Sums and Riemann Integrals

Lesson 5.3 Graph General Rational Functions

A signalling model of school grades: centralized versus decentralized examinations

Transcription:

Interntionl Conference Pssive nd Low Energy Cooling 133 The impct of wind on ir temperture distribution in Athens nd in Sntorini I. Livd, K. Nichou nd M. Sntmouris Lbortory of Meteorology, Section of Applied Physics, Physics Deprtment, University of Athens, Greece ABSTRACT Reserch studies of therml comfort in buildings im mostly t mintining ir temperture, especilly during summer period, within the cceptble comfort conditions nd with the minimum energy consumption. The contribution of wind is known to be n importnt prmeter for nturl ventiltion nd pssive cooling of buildings. In the present work the correltion between the vrition of wind nd ir temperture is studied for two different meteorologicl sttions in mjor Athens re, where the het islnd effect exists. The sme study is lso performed for nother sttion in Sntorini islnd nd comprison nlysis between the two regions is mde, since Sntorini is not exposed t the urbn environment nd furthermore it is directly ffected by the etesin winds nd the vicinity of the se. 1. INTRODUCTION The increse of living stndrds in buildings nd the persistence of high mbient temperture vlues during summer especilly in big cities hve incresed the instlltion of irconditioning systems nd consequently the energy requirements for cooling. Anlysis of summer ir tempertures nd wind velocities is thus necessry in order to provide pproprite informtion to building designers. Clcultion of cooling lod, pproprite sizing of AC systems nd interction of pssive nd hybrid cooling techniques nd components within the building re directly relted with the vilbility of locl climtic conditions, of which mbient temperture nd wind velocity re the most importnt prmeters. Tselepidki nd Sntmouris (1991) reported tht the frequency of high summer ir temperture vlues exceeding o C in Athens during June to September rnges from 0.32 to 0.65. Besides the persistence of existing consecutive dys with high ir tempertures ws found sttisticlly significnt during this period. Especilly for July nd August it ws estimted tht cooling is necessry for more thn five hours per dy. Blrs et l. (1992, 1993) estimted the cooling power index for Athens, which ws proposed by Vinje. This index is bsed on ir temperture nd wind speed vlues nd corresponds to humn senstion. It ws found tht the occurrence of hot conditions over dy is sttisticlly significnt during summer period. The previling high ir tempertures together with low wind speeds ssist the persistence of hot conditions of high number of consecutive hours. Tzikopoulos et l. (04) ttempted to model the energy efficiency of bioclimtic buildings bsed on the estimted degree dys of Greece nd other Mediterrnen res nd the rest of Europe. The purpose of the present study is twofold. The cooling degree hours, DH(26), re estimted for two sttions in the extended Athens re nd for n insulr sttion in Sntorini islnd. Besides the impct of wind on the observed ir temperture distribution is nlyzed during June to September for five yer period (00-04). 2. METEOROLOGICAL DATA The study of hourly ir temperture nd wind speed vlues is of concern not only on tempo-

134 Interntionl Conference Pssive nd Low Energy Cooling rl but lso on sptil bsis. Tselepidki et l. (1992, 1994) observed importnt differentitions in neighbor regions in mjor Athens re. However, comprison nlysis between summer conditions in Athens nd in islnds is lso very interesting. Livd et l. (1998) found tht the men return period of dily mximum ir tempertures, exceeding 26 o C is bout two dys in Cycldes for the period My to September. In the present work the hourly ir temperture nd wind speed dt were used for the period from June to September 00-04 of the Ntionl Observtoty of Athens (NOA), s well s, of Ellinikon sttion tht it is locted in rurl re ner Sronikos gulf. Also nother islnd sttion is studied which is plced in Sntorini, in the southern prt of Cycldes (Fig. 1), where strong winds nd the bsence of het islnd effect re the most predominnt chrcteristics. The selection of the exmined period, 00-04, is bsed on the fct tht in ll studied sttions the men monthly ir tempertures were 0.5 o C to 3.5 o C greter thn the corresponding verge climtologicl vlues (Tselepidki et l., 1993). 3. ANALYSIS OF DEGREE HOURS The number of cooling degree hours with reference 26 o C, DH(26), is given in Tble 1 where it Figure1: Position of the three meteorologicl sttions in Athens nd in Sntorini islnd. is shown tht NOA is chrcterized by greter vlues from % up to 62% in comprison with Ellinikon sttion tht is plced in n open re djcent to the se t distnce of Km from NOA. In spite of the fct tht NOA is plced in n re with reltive vegettion, on the top of 417 meter high hill, its position ner the centre of Athens, mkes it exposed to the het islnd effect (Livd et l., 02). The estimted DH(26) vlues for NOA sttion nd for the period 00-04 re 1.9 to 2.7 greter thn the corresponding vlues for the period 1977-1986 (Tselepidki et l., 1993). This cn be explined by the higher ir tempertures nd the increse of the het islnd effect. In ddition, compring the DH(26) in Ellinikon sttion with the corresponding vlues in Izmir, in Turkey (Stmn nd Ylcinky, 00), which is plced in by of Aegen se t lmost the sme ltitude with Ellinikon, it ws found tht the DH (26) re of the sme order for the two sttions. The observed smll number of DH(26) vlues in Sntorini shows tht the bsence of the het islnd effect fvors the predominnce of lower ir tempertures during summer. Moreover the vicinity to the se, where the surfce tempertures rnge between 22 o C nd o C (Theohrtos nd Tsepidki-Livd, 1990) results in lower mbient temperture vlues in the whole region. Seeing tht during summer the locl Etesin winds (dily strong winds) re predominnt in Cycldes, it is rised question whether these winds fvor the mintennce of lower ir tempertures. Thus, the reltion between the mbient tempertures nd wind velocities is nlyzed in the following prgrphs. 4. IMPACT OF WIND ON AIR TEMPERA- TURE DISTRIBUTION The reltion between wind speed nd ir temperture ws studied for the totl dt set nd thereinfter for ech wind direction. Nmely, in Athens existing dt were clssified in three different dt sets bsed on wind direction reltive to the rel North (293 o - o, 46 o -1 o nd 136 o - 292 o ), while clms (wind speed lower thn 0.3m/s) were nlyzed seprtely. The purpose ws to study the impct of ) northern winds (293 o - o ) which correspond to locl Etesin winds during summer, b) SSE-WNW (136 o - 292 o ) winds which blow from the se, c) ENE-

Interntionl Conference Pssive nd Low Energy Cooling 1 Tble 1: Number of Dh(26) for NOA, Ellinikon nd Sntorini sttions during the period 00-04. NOA ELLINIKON SANTORINI J J A S Sum J J A S Sum J J A S Sum 00 1786 38 65 838 8647 1338 2219 2287 543 6387 8 1275 395 290 2268 01 14 32 50 9 8586 814 2478 32 6 6469 131 10 995 419 29 02 55 2898 1684 198 65 1219 2132 1863 280 5494 187 13 11 111 2824 03 2164 2742 89 621 8616 1051 1763 1953 279 5046 380 959 728 233 04 1270 48 43 600 6461 580 1797 1829 418 4624 216 548 808 90 1662 Men 1616 2990 2486 637 7729 1000 78 93 433 5604 244 1101 8 229 Tble 2: Slopes of the liner correltion between wind nd ir temperture vlues for NOA, Ellinikon nd Sntorini sttions during dy nd night (vlues in prenthesis) period. Sttions Wind Direction June July August September 136-292 o 0.26(-0.08) 0.05(-0.001) 0.62(0.) 0.60(0.41) ΝΟΑ 293- o 0.42(0.17) 0.42(0.17) 0.53(0.26) 0.29(-0.13) 46-1 o 0.71(0.49) 0.94(0.21) 0.57(0.59) 0.33(-0.07) 136-292 o -0.37(-1.19) -0.49(-0.32) -0.09(0.94) 0.06(0.) Ellinikon 293- o 0.12(0.37) 0.36(0.26) 0.42(0.) 0.09(0.32) 46-1 o 0.56(0.09) 1.03(0.14) 0.74(0.04) -0.08(-0.) Sntorini 0-179 ο -0.22(-0.03) -0.14(-0.28) -0.16(-0.) -0.26(-0.27) 180-360 ο -0.19(0.001) -0.(0.07) -0.12(-0.02) -0.22(0.02) SE (46 o -1 o ) winds tht re ffected by Ymittos, 1027m high mountin plced est of Athens nd d) clms on the mbient tempertures. In Sntorini the clssifiction of wind directions ws mde in two sectors (0 o -179 o nd 180 o - 9 o ) depending on whether or not the sttion is ffected directly by the se. Figure 2 shows tht in Sntorini, ir temperture is reduced when wind speed is incresed, independently of wind direction nd especilly for wind speeds greter thn 5m/s, either during dy or night period for the whole summer period. Not to mention tht the reduction of mbient temperture ws greter when the wind ws blowing from 0 o to 179 o, from the seside, rther thn when it ws blowing from 180 o to 9 o, where the ir strem followed continentl route (Tble 2). In Ellinikon sttion (Fig. 3), it is observed tht in ll cses (except for the dy period during September) ir temperture is incresed with wind speed, independently of wind direction. On first sight this seems prdox, however, when the sme nlysis ws performed for winds blowing from the se (SSE-WNW directions) then the mbient tempertures were reduced for incresing wind speeds during the 24 hour period (except for September). On the contrry, when the wind ws from ENE to SE directions where the effect of topogrphy ws existing then higher mbient temperture corresponded to higher wind speeds (Tble 2). Finlly, the existing northern wind speeds cused higher mbient tempertures, lthough this ws not expected. This my be ttributed to the fct the northern irstrem is mixed with the hot ir from the urbn cnopy lyer nd thus high ir tempertures re observed independently of wind speed. In cse of NOA sttion (Fig. 4), it ws found tht for ll wind directions nd during the 24- hour period, the mbient temperture is proportionl to the logrithm of wind speed (except for the night period during September). The increse of ir temperture is greter when the wind speed is from ENE to SE directions when the ctbtic winds from Ymittos mountin become wrmer. During the night period the increse is smller wheres, for winds from SSE to WNW positive correltion is found for August-September (Tble 2). During September when the se is chrcterized by higher surfce tempertures then the ir mss is heted underneth nd this is strengthen by the het islnd effect nd thus higher ir tempertures re observed t NOA sttion for high wind speeds either during the dy or night period. On the contrry during the night period for northern or estern wind directions, lower ir tempertures correspond to higher wind speeds. In generl, it hs been found tht for lower wind speeds there is n brupt increse of ir temperture, while when the wind speed is in-

136 Interntionl Conference Pssive nd Low Energy Cooling Sntorini - June (v>4.9m/s) (v>4.9m/s) Sntorini - July (v>4.9m/s) (v>4.9m/s) 18 18 b Sntorini - August (v>4.9m/s) (v>4.9) Sntorini - September (v>4.9m/s) (v>4.9m/s) c 18 d 18 Figure 2: Correltion of wind nd ir temperture vlues for Sntorini sttion (the bigger tringles nd squres correspond to V>4.9m/s). Ellinikon - June (v>4.9m/s) (v>4.9) Ellinikon - June (v>4.9m/s) (v>4.9) Ellinikon - August Ellinikon - September (v>4.9m/s) (v>4.9m/s) (v>4.9m/s) (v>4.9m/s) Tir( o C) c cresed bove threshold vlues then the ir temperture tends to stbilize presenting tendency to cncel the islnd effect. Figure 3: Correltion of wind nd ir temperture vlues for Ellinikon sttion (the bigger tringles nd squres correspond to V>4.9m/s). d 5. CONCLUSIONS From the study of Dh(26) in three meteorologicl sttions in Athens nd in Sntorini, the het islnd effect is evident in the centre of Athens with the totl number of Dh(26) vlues in NOA rnging from 65 to 8647 during the cooling

Interntionl Conference Pssive nd Low Energy Cooling 137 NOA - June dy night NOA - July dy night Tir ( o C) b NOA - August Dy NOA - September Dy Tir ( o C) c d Figure 4: Correltion of wind nd ir temperture vlues for NOA sttion (the bigger tringles nd squres correspond to V>4.9m/s). period of the yer. In n open re djcent to the seside t distnce of Km from the city centre the totl number Dh(26) is reduced vrying from 4624 up to 6469, however these vlues re considered reltive high when compred with the corresponding vlues for Sntorini (1662-29). This revels tht the het islnd effect is predominnt in the extended Athens re, even with smller intensity t rurl plces. The impct of wind speed on ir temperture distribution in Sntorini showed tht higher wind speed vlues correspond to lower mbient tempertures, either during the dy or night period nd especilly when winds flow from the seside. In Ellinikon sttion negtive correltion is found when the wind blows from the se nd especilly during June, July nd August. For ll other wind directions the ir msses re directly ffected by the het islnd effect, presenting its greter intensity in the centre of Athens nd thus resulting in higher ir tempertures independently of wind speed. Finlly, for NOA sttion ner the centre of Athens, the het islnd effect overrides the influence of wind speed during the 24-hour period nd this is mostly observed for lower wind speeds, while in no cse negtive correltion ws observed. REFERENCES Blrs, C., I. Tselepidki Livd, M. Sntmouris nd D. Asimkopoulos, 1992. Clcultions nd sttisticl nlysis of the environmentl cooling power index for Athens,Greece. J. Energy Conv. nd Mngement, 34, 2, pp. 139-146 Blrs, C., I. Tselepidki-Livd, M. Sntmouris nd D. Asimkopoulos, 1993. Anlysis of therml comfort conditions in Athens, Greece. Energy Conv. nd Mngement, 34, 4, pp. 281-285. Livd, I., S. Lykoudis, M. Sntmouris nd C. Blrs, 1998. High mbient ir temperture frequency distribution t Hellenic islnds. Energy nd Buildings, 28, pp. 119-126. Livd, I., M. Sntmouris, K. Nichou, N. Ppnikolou nd G. Mihlkkou, 02. Determintion of plces in the gret Athens re where the het islnd effect is observed, Theoreticl nd. Applied Climtology, 71, pp. 219-2. Stmn, A. nd N. Ylcinky, 1999. Heting nd cooling degree-hours for Turkey, Energy, 24, 10, pp. 833-8. Theohrtos, G. nd I. Tsepidki-Livd, 1990.The temperture of the se-surfce in the region of the Aegen nd its reltion to ir-temperture, Theoreticl nd. Applied Climtology, 42, pp. 117-119. Tselepidki-Livd, I. nd M. Sntmouris, 1991. Sttis-

138 Interntionl Conference Pssive nd Low Energy Cooling ticl nd persistence nlysis of high summer mbient tempertures in Athens for cooling purposes, Energy nd Buildings, 17, pp. 243-1. Tselepidki-Livd, I., M. Sntmouris, C. Moustris nd G. Poulopoulou, 1992. Anlysis of the summer discomfort index in Athens, Greece for cooling purposes, Energy nd Buildings, 18, pp. 51-56. Tselepidki-Livd, I., M. Sntmouris nd D. Melitsiotis, 1993. Anlysis of the summer mbient tempertures for cooling purposes, Solr Energy, 50, 3, pp. 197-4. Tselepidki-Livd, I. nd M. Sntmouris, 1993. Anlysis of summer climtic conditions in Athens, Greece for cooling purposes, Int.Conf. Solr Energy in Architecture nd Urbn Plnning, Florence, Itly, 17-21/5/1993, Proc., pp. 122-1. Tselepidki-Livd, I., M. Sntmouris, D. Asimkopoulos nd S. Kontoyinnidis, 1994. On the vribility of cooling degree-dys in n urbn environment: Appliction to Athens, Greece, Energy nd Buildings, 21, pp. 93-99. Tzikopoulos, A.F., M.C. Krtz nd J.A. Prvntis, 05. Modeling energy efficiency of bioclimtic buildings, Energy nd Buildings, 37, 5, pp. 529-544.