Measurements of the angular distribution of diffuse irradiance

Similar documents
Solar radiation and thermal performance of solar collectors for Denmark

Evaluation of long-term global radiation measurements in Denmark and Sweden

An Equivalent Source Method for Modelling the Lithospheric Magnetic Field Using Satellite and Airborne Magnetic Data

Modelling Salmonella transfer during grinding of meat

Validation of Boundary Layer Winds from WRF Mesoscale Forecasts over Denmark

Reducing Uncertainty of Near-shore wind resource Estimates (RUNE) using wind lidars and mesoscale models

Kommerciel arbejde med WAsP

Lidar calibration What s the problem?

Numerical modeling of the conduction and radiation heating in precision glass moulding

Radioactivity in the Risø District July-December 2013

Exposure Buildup Factors for Heavy Metal Oxide Glass: A Radiation Shield

Global Wind Atlas validation and uncertainty

Sound absorption properties of a perforated plate and membrane ceiling element Nilsson, Anders C.; Rasmussen, Birgit

Geometry of power flows and convex-relaxed power flows in distribution networks with high penetration of renewables

Aalborg Universitet. Lecture 14 - Introduction to experimental work Kramer, Morten Mejlhede. Publication date: 2015

Accelerating interior point methods with GPUs for smart grid systems

Validation and comparison of numerical wind atlas methods: the South African example

Diffusion of sulfuric acid in protective organic coatings

The dipole moment of a wall-charged void in a bulk dielectric

Optimizing Reliability using BECAS - an Open-Source Cross Section Analysis Tool

Aalborg Universitet. CLIMA proceedings of the 12th REHVA World Congress Heiselberg, Per Kvols. Publication date: 2016

Radioactivity in the Risø District January-June 2014

Design of Robust AMB Controllers for Rotors Subjected to Varying and Uncertain Seal Forces

Optimizing the Photovoltaic Solar Energy Capture on Sunny and Cloudy Days Using a Solar Tracking System

A Note on Powers in Finite Fields

Vindkraftmeteorologi - metoder, målinger og data, kort, WAsP, WaSPEngineering,

Realtime 3D stress measurement in curing epoxy packaging

Optimal acid digestion for multi-element analysis of different waste matrices

A Procedure for Classification of Cup-Anemometers

Application of the wind atlas for Egypt

Aalborg Universitet. Comparison between Different Air Distribution Systems Nielsen, Peter Vilhelm. Publication date: 2006

SU solar resource measurement station: Sonbesie metadata

Application of mesoscale models with wind farm parametrisations in EERA-DTOC

Simulating wind energy resources with mesoscale models: Intercomparison of stateof-the-art

Matrix representation of a Neural Network

On numerical evaluation of two-dimensional phase integrals

Feasibility of non-linear simulation for Field II using an angular spectrum approach

Evaluation the impact tilt angle on the sun collectors

A systematic methodology for controller tuning in wastewater treatment plants

BSRN STATION DESCRIPTION

Aalborg Universitet. Published in: Proceedings of Roomvent Publication date: 2007

Generating the optimal magnetic field for magnetic refrigeration

Robustness of a cross contamination model describing transfer of pathogens during grinding of meat

Determination of packaging induced 3D stress utilizing a piezocoefficient mapping device

Impact of the interfaces for wind and wave modeling - interpretation using COAWST, SAR and point measurements

GHI CORRELATIONS WITH DHI AND DNI AND THE EFFECTS OF CLOUDINESS ON ONE-MINUTE DATA

Crack Tip Flipping: A New Phenomenon yet to be Resolved in Ductile Plate Tearing

Reflection of sound from finite-size plane and curved surfaces

Temporal and Spatial Variations of the Solar Radiation Observed in Singapore

Potential solution for rain erosion of wind turbine blades

On the Einstein-Stern model of rotational heat capacities

Wake decay constant for the infinite wind turbine array Application of asymptotic speed deficit concept to existing engineering wake model

Comparisons of winds from satellite SAR, WRF and SCADA to characterize coastal gradients and wind farm wake effects at Anholt wind farm

Steady State Comparisons HAWC2 v12.2 vs HAWCStab2 v2.12

Prediction of airfoil performance at high Reynolds numbers.

Reflection and absorption coefficients for use in room acoustic simulations

Thermal Analysis of Bending Under Tension Test

Extreme sea levels and the assessment of future coastal flood risk

Benchmarking of optimization methods for topology optimization problems

Design possibilities for impact noise insulation in lightweight floors A parameter study

Aalborg Universitet. Publication date: Document Version Publisher's PDF, also known as Version of record

Dynamic Effects of Diabatization in Distillation Columns

A comparison study of the two-bladed partial pitch turbine during normal operation and an extreme gust conditions

Application of the MCNPX-McStas interface for shielding calculations and guide design at ESS

Berøringsløs optisk måling af gassammensætning ved UV og IR spektroskopi

ME 476 Solar Energy UNIT THREE SOLAR RADIATION

Battling Bluetongue and Schmallenberg virus Local scale behavior of transmitting vectors

Quality control of airborne geophysical data from the EU Mining Sector Support Programme, Ghana

Unknown input observer based scheme for detecting faults in a wind turbine converter Odgaard, Peter Fogh; Stoustrup, Jakob

Single Shooting and ESDIRK Methods for adjoint-based optimization of an oil reservoir

Concept Testing of a Simple Floating Offshore Vertical Axis Wind Turbine

Multi (building)physics modeling

Strongly 2-connected orientations of graphs

Extreme Design Loads Calibration of Offshore Wind Turbine Blades through Real Time Measurements

Neutronics calculations for the ITER Collective Thomson Scattering Diagnostics

Higher-order spectral modelling of the diffraction force around a vertical circular cylinder

The Wien Bridge Oscillator Family

Analytical Studies of the Influence of Regional Groundwater Flow by on the Performance of Borehole Heat Exchangers

Global, direct and diffuse radiation measurements at ground by the new Environmental Station of the University of Rome Tor Vergata

Alkali Release from Typical Danish Aggregates to Potential ASR Reactive Concrete

Chapter 2 Available Solar Radiation

Characterization of the Airflow from a Bottom Hung Window under Natural Ventilation

ScienceDirect. Experimental study on the sun tracking ability of a spherical solar collector

Multiple Scenario Inversion of Reflection Seismic Prestack Data

On machine measurements of electrode wear in micro EDM milling

Plane Wave Medical Ultrasound Imaging Using Adaptive Beamforming

Available online at ScienceDirect. Energy Procedia 81 (2015 )

Investigation of a Gas-Solid Separation Process for Cement Raw Meal

Citation for published version (APA): Burcharth, H. F. (1988). Stress Analysis. Aalborg: Department of Civil Engineering, Aalborg University.

The IEA Annex 20 Two-Dimensional Benchmark Test for CFD Predictions

Development of a Dainish infrastructure for spatial information (DAISI) Brande-Lavridsen, Hanne; Jensen, Bent Hulegaard

A practical implementation of microphone free-field comparison calibration according to the standard IEC

Solar Resource Mapping in South Africa

Graphs with given diameter maximizing the spectral radius van Dam, Edwin

Radiometric Measurements of Environmental Radioactivity Beta Counting, Alpha and Gamma Spectrometry

Aalborg Universitet. Windows Heiselberg, Per Kvols; Svidt, Kjeld; Nielsen, Peter Vilhelm. Publication date: 2000

Validity of PEC Approximation for On-Body Propagation

Aalborg Universitet. Land-use modelling by cellular automata Hansen, Henning Sten. Published in: NORDGI : Nordic Geographic Information

Transcription:

Downloaded from orbit.dtu.dk on: Nov 02, 2018 Measurements of the angular distribution of diffuse irradiance Nielsen, Elsabet Nomonde Noma; Nielsen, Kristian Pagh ; Dragsted, Janne; Furbo, Simon Published in: Energy Procedia Publication date: 2015 Document Version Publisher's PDF, also known as Version of record Link back to DTU Orbit Citation (APA): Andersen, E., Nielsen, K. P., Dragsted, J., & Furbo, S. (2015). Measurements of the angular distribution of diffuse irradiance. Energy Procedia, 70, 729-736. General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

Available online at www.sciencedirect.com ScienceDirect Energy Procedia 70 (2015 ) 729 736 International Conference on Solar Heating and Cooling for Buildings and Industry, SHC 2014 Measurements of the angular distribution of diffuse irradiance Elsa Andersen a, Kristian Pagh Nielsen b, Janne Dragsted a, Simon Furbo a a Technical University of Denmark, Institute of Civil Engineering, Nordvej 119, 2800, Kgs. Lyngby, Denmark b Danish Meterological Institute, Lyngbyvej 100, 2100 Copenhagen Ø, Denmark Abstract Advanced solar resource assessment and forecasting is necessary for optimal solar energy utilization. In order to investigate the short-term resource variability, for instance caused by clouds it is necessary to investigate how clouds affect the solar irradiance, including the angular distribution of the solar irradiance. The investigation is part of the Danish contribution to the task 46 within the International Energy Agency and financed by the Danish Energy Agency. The investigation focuses on the distribution of the diffuse solar irradiance and is based on horizontal measurements of the solar irradiance from 8 different parts of the sky as well as horizontal measurements of the total beam and total diffuse irradiance. 2015 2015 The The Authors. Authors. Published Published by by Elsevier Elsevier Ltd. Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). Peer-review by the scientific conference committee of SHC 2014 under responsibility of PSE AG. Peer-review by the scientific conference committee of SHC 2014 under responsibility of PSE AG Keywords: Measurement of diffuse solar irradiance; Angular distribution of diffuse solar irradiance. 1. Introduction The International Energy Agency (IEA) Solar Heating and Cooling programme (SHC) Task 46: Solar Resource Assessment and Forcasting has two major goals: 1) To understand short-term resource variability caused by clouds and 2) To assess the application of short term forecast for planning the operation of multiple solar energy systems. In order to investigate the short-term resource variability caused by clouds it is necessary to investigate how fractional 3-dimensional cloud fields affect the solar irradiance, including the angular distribution of solar irradiance. This investigation is done by applying high temporal resolution measurements of the solar irradiance from ground based stations operated by the Technical University of Denmark, Department of Civil Engineering (DTU) in Kgs. Lyngby and in Sisimiut, Greenland, and at a number of stations across Denmark operated by the Danish Meterological Institute (DMI) in Copenhagen. 1876-6102 2015 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). Peer-review by the scientific conference committee of SHC 2014 under responsibility of PSE AG doi:10.1016/j.egypro.2015.02.182

730 Elsa Andersen et al. / Energy Procedia 70 ( 2015 ) 729 736 At the DTU climate station in Kgs. Lyngby detailed measurements of solar irradiance have been performed for 25 years, providing a long set of historical data, which are very valuable in this investigation. In particular, DTU have focused on measuring both the global irradiance and the diffuse irradiance on horizontal. The objective of the investigation described in this paper is to better understand short-term resource variability caused by clouds. The short-term resource variability caused by clouds is investigated with a set of 8 high temporal resolution pyranometers that measure the total horizontal solar irradiance from 8 different sky directions and measurements of the beam normal and global irradiances. The detailed measurements are carried out at the DTU climate station. 2. The experimental set-up 2.1. The climate station The climate station is located at the technical university of Denmark. The latitude is 55.6 N and the longitude is 12.5 E. See in Fig. 1, left. The total diffuse irradiance on horizontal is measured with a pyranometer mounted on a sun tracker that screens of the beam irradiance. The beam normal irradiance is measured with a pyrheliometer mounted on the sun tracker, see Fig 1, right. Fig. 1. (Left) The climate station. (Right) The sun tracker. The equipment is from Kipp and Zonen. The pyranometer is type CMP 11, the pyrheliometer is type CHP 1 and the sun tracker is type SOLYS 2. The sun tracker is equipped with shadow balls that screen of the beam irradiance from the pyranometer mounted on the horizontal top plate on the sun tracker and the pyranometer measures only the diffuse irradiance. The sun tracker is further equipped with a sun sensor that fine tunes the position of the sun tracker. The view angle from the pyrheliometer to the sky is 5. The view angle from the pyranometer to the shadow ball is also 5. 2.2. Fractional irradiance measurement set-up The fractional solar irradiance is measured with 8 pyranometers also type CMP 11 from Kipp and Zonen. On top of each pyranometer, domes with openings corresponding to 1/8 part of the sky are mounted. The domes point in 8 different directions in such a way that solar irradiance is measured from the entire sky. The domes are half balls with diameters of 291 mm. The domes are held in position with metal plates with multiple holes which allows for good ventilation inside the domes. The inside of the domes and the metal plate is painted with black low reflection paint. The reflection from the black paint is < 5%. The domes are shiny white on the outside. Fig. 2 shows a top view of the climate station where it can be seen how the domes point in 8 different directions, covering the whole sky and

Elsa Andersen et al. / Energy Procedia 70 ( 2015 ) 729 736 731 pictures of the fractional measurement set-up. The sun tracker is located in the North West corner of the climate station. Fig. 2. (Left) Top view of the climate station. (Right) The experimental set-up with the dome mounted on top of the pyranometer. 3. Measurements 3.1. Side-by-side calibration measurements Initially the pyranometers are tested side-by-side in order to investigate if the measured values from the different pyranomerers are identical when the pyranometers are exposed to the exact same weather conditions. Fig. 3 shows a picture of 16 pyranometers mounted on horizontal at the climate station. 8 of the pyranometers are used for the fractional horizontal measurements and 1 pyranometer is used for measuring the total horizontal diffuse irradiance with the sun tracker. Fig. 4 shows the measurements from all 16 pyranometers on a sunny day. It can be seen that the pyranometers measure exactly the same global solar irradiance, except in the morning and in the evening with low solar elevation where shading from the climate station influence the measurements. Fig. 3. The 16 pyranometers mounted on horizontal for side-by-side calibration measurements.

732 Elsa Andersen et al. / Energy Procedia 70 ( 2015 ) 729 736 Fig. 4. Side-by-side measurements of the global irradiance with 16 different pyranometers on a sunny day. The pyranometers are named P1, P2,, P16. 3.2. Measurement period The measurement period is from June 19 2014 to October 6 2014. Unfortunately there are periods without measurements, see Table 1. Table 1. Measurement period and missing data during the measurement period. Measurement period in 2014 Number of missing days The missing days June 19 June 30 1 30. July 1 July 31 11 1.-11. August 1 August 31 6 19.-20., 26.-29. September 1 September 30 0 October 1 October 6 0 3.3. Fractional horizontal diffuse irradiance measurements The fractional horizontal diffuse irradiance is calculated by subtracting the horizontal beam irradiance from the fractional measurement which is exposed to beam irradiance. The geometry of the domes result in shadows on the pyranometer close to the beginning and the end of the dome opening, see Fig 5. Measurements from the shaded areas are left out and the fractional horizontal diffuse irradiance is interpolated between the unshaded areas.

Elsa Andersen et al. / Energy Procedia 70 ( 2015 ) 729 736 733 Pyranometer Fig. 5. (Left) Picture of dome with the opening. (Right) Schematic illustration of how the dome leads to shading on the pyranometer where the dome opening begins and ends. 3.4. Investigations on measured inaccuracies and correction of the measured inaccuracies During cloudy days, the sum of the 8 fractional horizontal diffuse irradiance measurements is similar to the total horizontal diffuse irradiance. However, on sunny days, the sum of the 8 fractional horizontal diffuse irradiance measurements is larger than the total horizontal diffuse irradiance. Fig. 6 shows the fractional horizontal diffuse irradiance from the 8 pyranometers and the sum of the 8 fractional measurements and the horizontal diffuse measurement carried out with the sun tracker on a sunny day. In the figure, the areas shaded by the domes are left out. It is clear to see that the sum of the fractional horizontal diffuse irradiance measurements differ from the total horizontal diffuse irradiance when there is beam irradiance. 160 06-09-2014 140 120 100 80 60 40 20 0 00:00 01:04 02:08 03:12 04:16 05:20 06:24 07:28 08:32 09:36 10:40 11:44 12:48 13:52 14:56 16:00 17:04 18:08 19:12 20:16 21:20 22:24 23:28 Solar irradiance [W/m 2 ] N NE E SE S SW W NW Diffuse 8/8 Fig. 6. Measured fractional horizontal diffuse irradiance, the sum of the fractional irradiances and the total horizontal diffuse irradiance measured with the sun tracker on a sunny day. In order to investigate if a temperature difference between the inside of the dome where the glass dome of the pyranometer is located and outside the dome where the body of the pyranometer is located is causing the measured inaccuracies, measurement are compared between the experimental set-up with the dome and a similar set-up where the dome is removed and the glass dome of the pyranometer is covered leaving only 1/8 part open to the sky. The 1/8 openings in the two different experimental set-ups are pointing in the same direction. The measured irradiances are exactly the same when the sun is shining through the middle of the openings. Also smoke tests are carried out to

734 Elsa Andersen et al. / Energy Procedia 70 ( 2015 ) 729 736 see how fast the air is exchanged inside the dome. The smoke is ventilated away instantly. Based on these investigations it is concluded that the measured inaccuracies are due to interreflections inside the domes. The visible opening areas of the domes are projected on a plane perpendicular to the beam solar irradiance and the corrections are distributed according to the visible opening areas of the domes. Also the area of the glass dome of the pyranometer is subtracted from the visible opening area of the dome with the sun, See Fig. 7. It is clear to see that even if only one pyranometer is exposed directly to beam irradiance, the sun shines through the openings of the neighboring domes as well. In Fig. 8 the measured fractional horizontal diffuse irradiances and the sum of the fractional horizontal diffuse irradiances corrected for interreflections and interpolated between the shaded areas caused by the domes and the total horizontal diffuse irradiance measured with the sun tracker on a sunny day are shown. The dotted lines show the uncorrected measurements and the full blown lines show the corrected and interpolated measurements. The figure shows that the corrections result in a good degree of similarity between the sum of the fractional horizontal diffuse irradiance measurements and the total horizontal diffuse irradiance measured with the sun tracker. Plane perpendicular to the beam solar irradiance Fig. 7. Illustration of the applied method for correcting the measurements for interreflections inside the domes. Solar irradiance [W/m 2 ] 160 140 120 100 80 60 40 20 0 06-09-2014 N NE E SE S SW W NW Diffuse 8/8 N NE E SE S SW W NW Diffuse 8/8 Diffuse.tracker 00:00 01:06 02:12 03:18 04:24 05:30 06:36 07:42 08:48 09:54 11:00 12:06 13:12 14:18 15:24 16:30 17:36 18:42 19:48 20:54 22:00 23:06 Fig. 8. Measured fractional horizontal diffuse irradiances and the sum of the fractional irradiances corrected for interreflections and interpolated between the shaded areas and the total horizontal diffuse irradiance measured with the sun tracker on a sunny day. The dotted lines show the uncorrected measurements and the full blown lines show the corrected measurements.

Elsa Andersen et al. / Energy Procedia 70 ( 2015 ) 729 736 735 4. Results The fractional horizontal diffuse solar radiation measurements from the whole measurement period are shown in Fig. 9. The figure shows as expected that most diffuse solar radiation comes from the pyranometers with the south facing dome openings. 30 June 19 2014 - October 6 2014 Horizontal diffuse solar radiation [kwh/m 2 ] 25 20 15 10 5 0 N NE E SE S SW W NW Fig. 9. Measured fractional horizontal diffuse solar radiation from the whole measurement period. The solar diffuse fraction, SDF is defined as shown in eq. (1). Horizontaldiffusesolar radiation from thedomeopening with thesun SDF (1) Total horizontaldiffusesolar radiation Fig. 10 shows the daily solar diffuse factor as function of the daily beam solar radiation on horizontal. The figure shows that the daily horizontal diffuse solar radiation is highly influenced by the daily beam solar radiation on horizontal. From the figure it can be seen that during sunny days, most horizontal diffuse solar radiation comes from the area around the sun. On sunny days, the solar diffuse fraction can be up to 42 % and the figure also shows that there are large variations from day to day. Hence, the distribution of the horizontal diffuse irradiance on sunny days is unevenly distributed or non-isotropic distributed. The figure also shows that during days without any beam irradiance, the horizontal solar diffuse factor is around 12.5 % which exactly correspond to the opening area of a dome. Hence, it is seen that during days without any beam irradiance the horizontal diffuse irradiance is evenly distributed or isotropic distributed over the azimuth direction.

736 Elsa Andersen et al. / Energy Procedia 70 ( 2015 ) 729 736 45 40 Daily solar diffuse fraction [%] 35 30 25 20 15 10 5 0 0 1 2 3 4 5 6 7 Daily beam radiation on horizontal [kwh/m 2 ] Fig. 10. Daily horizontal solar diffuse fraction from the whole measurement period as function of the measured daily beam radiation on horizontal. 5. Conclusions The fractional diffuse solar irradiances on horizontal from 8 different parts of the sky were investigated by means of ground based measurements with high temporal resolution pyranometers. Also the total horizontal diffuse and beam irradiances were measured. The measurements were carried out at the climate station located at the Technical University of Denmark at latitude 55.6 N and longitude12.5 E. The measurement period is from June 19 2014 to October 6 2014. The fractional horizontal measurements are facilitated by experimental set-ups with domes with openings corresponding to 1/8 part of the sky mounted over the glass domes of the pyranometers. The domes are half balls with diameters of 291 mm painted black on the inside with low reflection paint. The experimental set-ups result in shading on the pyranometers at the beginning and end of the dome openings and some measured inaccuracies caused by unintended interreflections inside the domes during periods with beam irradiance. The measurements are corrected for the shades and the interreflections, resulting in a good degree of similarity between the sum of the fractional diffuse irradiance on horizontal and the total diffuse irradiance on horizontal. The results show as expected that the distribution of horizontal diffuse irradiance is highly influenced by beam irradiance. On sunny days the horizontal diffuse irradiance is much higher around the sun than in the rest of the sky. The solar diffuse fraction around the sun, that is: the sum of the horizontal diffuse solar radiation from the dome opening with the sun divided with the total horizontal diffuse radiation can be up to 42% and the variations are large from day to day. On cloudy sky days without any beam irradiance, the horizontal diffuse irradiance is evenly distributed over all azimuth directions. The measurements will be continued, analyzed further and used to improve solar radiation models for estimating beam and diffuse irradiance on horizontal and on tilted surfaces.