Estimation of Global Solar Radiation in Onitsha and Calabar Using Empirical Models

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

Review Topic 14: Relationships between two numerical variables

ANALYSIS AND MODELLING OF RAINFALL EVENTS

AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix

A Study on the Properties of Rational Triangles

QUADRATIC EQUATION. Contents

THE INFLUENCE OF MODEL RESOLUTION ON AN EXPRESSION OF THE ATMOSPHERIC BOUNDARY LAYER IN A SINGLE-COLUMN MODEL

1B40 Practical Skills

SECTION A STUDENT MATERIAL. Part 1. What and Why.?

Lecture Notes No. 10

Numbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point

Section 6: Area, Volume, and Average Value

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES

Project 6: Minigoals Towards Simplifying and Rewriting Expressions

Linear Algebra Introduction

Introduction to Olympiad Inequalities

for all x in [a,b], then the area of the region bounded by the graphs of f and g and the vertical lines x = a and x = b is b [ ( ) ( )] A= f x g x dx

Activities. 4.1 Pythagoras' Theorem 4.2 Spirals 4.3 Clinometers 4.4 Radar 4.5 Posting Parcels 4.6 Interlocking Pipes 4.7 Sine Rule Notes and Solutions

AP CALCULUS Test #6: Unit #6 Basic Integration and Applications

University of Sioux Falls. MAT204/205 Calculus I/II

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx

Section 4: Integration ECO4112F 2011

Intermediate Math Circles Wednesday 17 October 2012 Geometry II: Side Lengths

Chapter Gauss Quadrature Rule of Integration

6.5 Improper integrals

Lesson 2: The Pythagorean Theorem and Similar Triangles. A Brief Review of the Pythagorean Theorem.

H (2a, a) (u 2a) 2 (E) Show that u v 4a. Explain why this implies that u v 4a, with equality if and only u a if u v 2a.

Integration. antidifferentiation

( ) { } [ ] { } [ ) { } ( ] { }

CS 491G Combinatorial Optimization Lecture Notes

Generalization of 2-Corner Frequency Source Models Used in SMSIM

Polynomials and Division Theory

7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus

7.1 Integral as Net Change Calculus. What is the total distance traveled? What is the total displacement?

Comparing the Pre-image and Image of a Dilation

QUADRATIC EQUATION EXERCISE - 01 CHECK YOUR GRASP

TOPIC: LINEAR ALGEBRA MATRICES

Interpreting Integrals and the Fundamental Theorem

Chapter 9 Definite Integrals

2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).

MAT 403 NOTES 4. f + f =

Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.

2. There are an infinite number of possible triangles, all similar, with three given angles whose sum is 180.

Green s Theorem. (2x e y ) da. (2x e y ) dx dy. x 2 xe y. (1 e y ) dy. y=1. = y e y. y=0. = 2 e

The practical version

Linear Inequalities. Work Sheet 1

Torsion in Groups of Integral Triangles

Now we must transform the original model so we can use the new parameters. = S max. Recruits

Factorising FACTORISING.

A Lower Bound for the Length of a Partial Transversal in a Latin Square, Revised Version

50 AMC Lectures Problem Book 2 (36) Substitution Method

MATH 122, Final Exam

8 THREE PHASE A.C. CIRCUITS

, g. Exercise 1. Generator polynomials of a convolutional code, given in binary form, are g. Solution 1.

PYTHAGORAS THEOREM WHAT S IN CHAPTER 1? IN THIS CHAPTER YOU WILL:

8.3 THE HYPERBOLA OBJECTIVES

5. Every rational number have either terminating or repeating (recurring) decimal representation.

( ) as a fraction. Determine location of the highest

1 This question is about mean bond enthalpies and their use in the calculation of enthalpy changes.

CHENG Chun Chor Litwin The Hong Kong Institute of Education

NEW CIRCUITS OF HIGH-VOLTAGE PULSE GENERATORS WITH INDUCTIVE-CAPACITIVE ENERGY STORAGE

9.1 Day 1 Warm Up. Solve the equation = x x 2 = March 1, 2017 Geometry 9.1 The Pythagorean Theorem 1

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

September 13 Homework Solutions

On the Scale factor of the Universe and Redshift.

ILLUSTRATING THE EXTENSION OF A SPECIAL PROPERTY OF CUBIC POLYNOMIALS TO NTH DEGREE POLYNOMIALS

First Semester Review Calculus BC

= x x 2 = 25 2

Section 4.4. Green s Theorem

Formula for Trapezoid estimate using Left and Right estimates: Trap( n) If the graph of f is decreasing on [a, b], then f ( x ) dx

Section 7.1 Area of a Region Between Two Curves

1.3 SCALARS AND VECTORS

Polynomials. Polynomials. Curriculum Ready ACMNA:

PYTHAGORAS THEOREM,TRIGONOMETRY,BEARINGS AND THREE DIMENSIONAL PROBLEMS

TOPPER SAMPLE PAPER - 5 CLASS XI MATHEMATICS. Questions. Time Allowed : 3 Hrs Maximum Marks: 100

Probability. b a b. a b 32.

LESSON 11: TRIANGLE FORMULAE

5: The Definite Integral

APPENDIX. Precalculus Review D.1. Real Numbers and the Real Number Line

I 3 2 = I I 4 = 2A

NON-DETERMINISTIC FSA

System Validation (IN4387) November 2, 2012, 14:00-17:00

Tutorial Worksheet. 1. Find all solutions to the linear system by following the given steps. x + 2y + 3z = 2 2x + 3y + z = 4.

Iowa Training Systems Trial Snus Hill Winery Madrid, IA

April 8, 2017 Math 9. Geometry. Solving vector problems. Problem. Prove that if vectors and satisfy, then.

( ) as a fraction. Determine location of the highest

Thomas Whitham Sixth Form

Figure 1. The left-handed and right-handed trefoils

Chapter 6 Continuous Random Variables and Distributions

M344 - ADVANCED ENGINEERING MATHEMATICS

( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x).

Algebra 2 Semester 1 Practice Final

Ellipses. The second type of conic is called an ellipse.

Section 1.3 Triangles

1. Twelve less than five times a number is thirty three. What is the number

Homework Assignment 6 Solution Set

QUADRATIC EQUATIONS OBJECTIVE PROBLEMS

Review of Gaussian Quadrature method

Effects of Drought on the Performance of Two Hybrid Bluegrasses, Kentucky Bluegrass and Tall Fescue

Transcription:

Communitions in Applied Sienes ISS 0-77 Volume, umer, 0, 5-7 Estimtion of Glol Solr dition in Onitsh nd Clr Using Empiril Models M.. nuhi, J. E. Ekpe nd G. F Ieh Deprtment of Industril Physis, Eonyi Stte University, Akliki, igeri Corresponding uthor: G. F Ieh, Deprtment of Industril Physis, Eonyi Stte University, Akliki, igeri ABSAC Understnding solr rdition dt is essentil for modelling solr energy systems. he purpose of the present study ws to estimte glol solr rdition on horizontl surfe using meteorologil prmeters for period of eleven yers (996 006). Monthly verge glol solr rdition, H hs een estimted for Onitsh nd Clr, igeri using prediting models generted y simple liner nd multiple regression nlyses. he models inluded one vrile model with reltive humidity s the independent vrile nd two three-vrile models orrelting H with reltive humidity nd eh of verge temperture, loudiness index nd numer of sunshine hours. he vlues of the glol solr rdition estimted y the models nd the mesured solr rdition were tested using the men is error (MBE), root men squre error (MSE) nd men perentge error (MPE) sttisl tehniques. he vlues of the orreltion oeffiient (CC) were lso determined for eh model. he model tht indite good greement etween the mesured nd estimted vlues for Onitsh nd Clr re nd respetively. he developed models n e used for estimting glol solr rdition in Onitsh nd Clr nd other lotions with similr limti ftors. Keywords: Sunshine, verge temperture, reltive humidity, loudiness, glol solr rdition Introdution Knowledge of the glol solr rdition is of fundmentl importne for ll solr energy onversion systems. Informtion on glol solr rdition reeived t ny site (preferly gined over long period) should e useful not only to the Copyright 0 the uthors. 5

Communitions in Applied Sienes 6 lolity where the rdition dt is olleted ut lso for the wider world ommunity (Mssquoi, 988). Prtilly, solr rdition dt re esily otined using the relevnt equipment. Pyrheliometer nd pyrnometer n e used redily to otin the diffused omponent of the rdition nd the glol solr rdition respetively. Wether sttions hve een used mostly for this purpose. However, there re few of suh sttions ross the gloe nd worse still in the developing ntions. In n effort to generte rdition dt, reserhers hd extrpolted vlues from one lotion for pplition in different lotion. Hene solr rdition predition from estimtion models hs een widely utilized glolly to generte solr rdition dtse for vrious lotions of the world. he development of the solr rdition dt se for vrious igeri lotions hs een n on-going tsk for reserhers in the field for mny yers now. With the very few meteorologil sttions, the option of using estimting models hs een widely dopted in igeri for prediting solr rdition t speifi lotion nd t regionl sle (Ago, et l 007). ESEACH MEHODOLOGY he dt suh s sunshine hours, verge tempertures, loud over, reltive humidity nd glol solr rdition dt for Onitsh nd Clr used for this study were otined from the igeri meteorologil Ageny, Federl ministry of Avition Oshodi, Lgos, igeri. he dt olleted overed period of eleven yers (996-006) for Onitsh (Ltitude, Longitude, nd ltitude 56 metres ove se level) nd Clr (Ltitude, Longitude nd ltitude 58 metres ove se level). hese vlues were otined y the use of GPS (generl position stellite) equipment. DAA AALYSIS he monthly dt proessed in preprtion for the orreltion re presented in les nd for Onitsh nd Clr respetively.

7 Communitions in Applied Sienes le : Monthly men vlue of limti prmeters for Onitsh 996-006 S/ MOH v( o C) % S/ /C H(MJ/m /dy) JA.6 67.56 0.5 0.56.5 FEB 5.57 69.77 0.55 0.56 5.65 MA.809 75.6 0.56 0.600.77 AP.555 79.8 0.58 0.69.7 5 MAY.6 8.77 0.58 0.6.85 6 JU.07 8.56 0.9 0.6.6 7 JULY 9.55 87.000 0.5 0.6.65 8 AUG 9.6 88.8 0.7 0.69 0.80 9 SEP 0.6 85.66 0. 0.67.6 0 OC.6 8.88 0.8 0.60 5.8 OV.500 75.909 0.580 0.65 6.5 DEC.86 7.7 0.55 0.576 5. otl 89.55 98.9 5.758 7.60 69. le : Monthly men vlues of limti prmeters for Clr (996-006) S/ MOH v( o C) % S/ /C H(MJ/m /dy) JA.8 76.6 0.89 0.609.00 FEB.809 76.55 0.9 0.6 6.7 MA.78 8.66 0.55 0.65 5.5 AP.955 8.7 0. 0.66 6.6 5 MAY.509 8.000 0.0 0.66 5. 6 JU 0.6 86.88 0. 0.6.09 7 JULY 8.6 90 0. 0.65.6 8 AUG 8.7 9.66 0.9 0.67.9 9 SEP 9.55 89.000 0. 0.6.9 0 OC 0.08 87.000 0.9 0.60. OV.056 85.7 0.6 0.6. DEC.909 80.55 0.58 0.6.6 otl 7.7 0.9.07 7.59 69.56

Communitions in Applied Sienes 8 v = monthly verge dily temperture (%) = eltive humidity @ 9 hours H = vlue of mesured verge dily solr rdition on the horizontl surfe. = he loudiness index C S = he monthly verge dily reltive sunshine durtion. o develop the model, the glol solr rdition mesured using Gun-Bellini distillte were onverted to useful form (MJM /dy) using onversion ftor of.6 proposed y (Smo, 985). Professionl modellers hve proposed orreltion tht should e used to estimte the glol solr rdition. he first orreltion proposed for estimting the monthly men glol solr rdition on the horizontl surfe H (MJ/m /dy) using the sunshine durtion dt ws done y (Angstrom, 9). (Presott, 90) hs put the Angstrom orreltion in more onvenient form s: ( ) () Where H is the mesured monthly men dily glol solr rdition, Ho is the monthly men extrterrestril solr rdition on horizontl surfe, n is the monthly men dily right sunshine hours. is the mximum possile monthly men dily sunshine hour or the dy length, is the lerness index, is the frtion of sunshine hours, nd is regression onstnts. A numer of orreltions whih inlude more meteorologil prmeters suh s mient temperture, the totl preipittion, reltive nd speifi humidity, mount of totl loud over et. hve een developed y different workers (Attili nd Adll, 99). he ove eqution hs een found to e very onvenient to lrge numer of lotions nd most widely used orreltion. he extrterrestril solr rdition on horizontl surfe is given y (Igl, 98) s written elow: H o I s E o 80 s sin sin os os sin s ()

9 Communitions in Applied Sienes Where Is is the solr onstnt, Eo is the eentriity orretion ftor, = ltitude, =solr delintion, s = hour ngle.(igl,98) gve the expressions for Is,Eo,, s follows: s 67 x 600 I s 000000 E o 0.0 os 60 65 s 60 8.5 sin 65 os 5 os tn tn tn tn Where is the hrteristi dy numer for eh month where nd re onstnts whih n e determined, Y is eing repled s Hm nd it is dependent vrile nd is the independent vrile whih n reple ny of the meteorologil dt like verge temperte, v nd reltive humidity,. Equtions () nd () re equtions of lest squre line nd lest squre prol or first nd seond order regressions respetively (Murry, 96).o exeute the regression nlysis of the first order, oth sides of eqution () hs to e multiply y nd suessively nd summing oth sides to otin Y 5 Y 6 If we pply our vriles like (v nd ) s the independent vrile in eqn. (.), we hve these equtions: v H 7 H 8 At this point, pplying verge temperture, v s our independent vrile in equtions (5) nd (6) we hve:

Communitions in Applied Sienes 0 0 9 v v v v H H Equtions (9) nd (0) will e used to evlute the regression onstnts nd he si eqution for seond order regression nlysis is known s lest squre prol. o otin these equtions we hve to multiply eqution (.) y, nd suessively, nd dding them to otin equtions (), () nd (). Y Y Y Applying these equtions to solve our prolem, tht is using our vriles, suh s verge temperture nd reltive humidity we hve 5 hen sustituting v s our independent vrile in equtions (), () nd (), we otin the eqution outlined elow: 8 7 6 v v v v v v v v v v H H H 0 9 H H H o further uttress the nlysis, the multiple orreltion nlysis ws lso pplied to evlute for H5, H6, nd H7. In this se, the multiple vriles suh

Communitions in Applied Sienes s verge temperture, nd reltive humidity, re used whih led to the eqution written elow. Where nd re oeffiients. If we tke = H, = nd H = v, we n now employ regression eqution of on nd nd rry out orreltion nlysis on eqution () written elow s By multiplying oth sides of eqution () y, nd suessively nd dding them to otin the outlined equtions: 5 6 Sustituting our vriles in the lst two equtions we rrive s; 5 5 v H 7 5 v 5 v H v 8 Eqution 7 nd 8 were used to evlute for H5. Similrly, H6we e estimted y using -S omintion, tht is reltive humidity nd sunshine hours whih leds to: H6 = 6 + 6S 9 he sme proedure used in otining equtions 7 nd 8 were pplied to rrived t equtions 0 nd. 6 6 S H 0 6 S 6 S 6 HS

Communitions in Applied Sienes Under this we n lso estimte y using nd in omintion. Hene, in rrying out the orreltion nlysis in equtions (7) nd (8) we otin: 7 7 C H 7 C 7 C HC he ury of the predited vlues ws tested y lulting the Men Bis Error (MBE), the oot Men Squre Error (MSE), the Men Perentge Error (MPE) nd the oeffiient of orreltion CC. esults le : Mesured nd Predited Solr dition for Onitsh, (996-006) S/ Month MJ/M /dy MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy JA.5 5.5 6. FEB 5.65 6.7 5.7 MA.77 6.66.75 AP.7.8.086 5 MAY.85.958.6 6 JU.6.65.5 7 JULY.65.9.76 8 AUG 0.80.07.50 9 SEP.6.777.96 0 OC 5.8.86.5 OV 6.5.787.656 DEC 5. 7 5.6.00..80 5.05 5.95.5.95.966.800.95.8.50.5.55.95.875.90.95.95.5 0.5 0.89 0.0.65.00.5 5.75 5.5 6.895 6.675 6.95 5.900 5.95.80.56.565.5.66.90 6..657 6.5 6.50.608.60.67.5 0. 0.60.59.5.99 5.05 6.0.59 5.0.89

Communitions in Applied Sienes le : Mesured nd Predited Solr dition for Clr, (996-006) S/ MO H MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy MJ/M / dy JA.00 5.099 5. FEB 6.7 6. 5.0 MA 5.5 5.6.09 AP 6.6.8.78 5 MAY 5..509.770 6 JU.09.56.5 7 JULY.6.70.60 8 AUG.9.6.9 9 SEP.9.850.80 0 OC..59.00 OV..9.58 DEC.6.79..85 6.05 5.50 6.50.550.5.55.880.005.5.95.50 5.0.968 5.05.58 5.97 5.895 5..07 5.7.875.08.99.666 6.75.067.8. 0..76.605.5.86 5.85.85.65.77 6.0..86. 5.9.98.55.80.8.88.5.67.5..95.099.78.767.5.0 le nd show the omprison etween the mesured glol solr rdition nd the predited one. A plot of the mesured glol solr rdition with the predited vlues gve rise to figure s nd for Onitsh nd Clr respetively, whih sientifilly gree well with le nd grphilly.

Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Communitions in Applied Sienes GAPHICAL EPESEAIO OF HE MEASUED AD PEDICED VALUES OF SOLA ADIAIO When the mesured nd predited vlues of solr rdition otined from model equtions were plotted, they showed lmost the sme urve. his reveled tht the model equtions n e used to predit the solr rdition of ny prt of the world tht possessessimilr limtologil ftors like tht of Onitsh nd Clr. he following plotted grphs illustrted these points etter. Eqution. (Onitsh) Mesured S 5 6 7 8 9 0 Months of the Yer Fig : Comprison etween the mesured nd predited vlues of orreltion eqution for Onitsh (996 006) Eqution. (Onitsh) 5 6 7 8 9 0 Months of the Yer Mesured S Eqution. (Onitsh) 5 7 9 Months of the Yer Mesured S Fig : Comprison etween the mesured nd predited vlues of orreltion eqution for Onitsh (996 006) Fig : Comprison etween the mesured nd predited vlues of orreltion eqution for Onitsh (996 006)

Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) 5 Communitions in Applied Sienes Eqution.5 (Onitsh) 5 7 9 Months of the Yer Mesured S Eqution. (Onitsh) Mesured S 5 7 9 Months of the Yer Fig d: Comprison etween the mesured nd predited vlues of orreltion eqution for Onitsh (996 006) Eqution. (Onitsh) Fig e: Comprison etween the mesured nd predited vlues of orreltion eqution for Onitsh (996 006) Fig f: Comprison etween the mesured nd predited vlues of orreltion eqution for Onitsh (996 006) 5 7 9 Months of the Yer 5 7 9 Months of the Yer Mesured S Eqution.8 (Onitsh) Mesured S Fig g: Comprison etween the mesured nd predited vlues of orreltion eqution for Onitsh (996 006) Fig : Comprison etween the mesured nd predited vlues of orreltion eqution for Clr (996 006) Fig : Comprison etween the mesured nd predited vlues of orreltion eqution for Clr (996 006) Eqution. (Clr) 5 7 9 Months of the Yer Eqution. (Clr) 5 7 9 Months of the Yer Mesured S Mesured S

Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Glol Solr dition (MJ/M /dy) Communitions in Applied Sienes 6 Fig : Comprison etween the mesured nd predited vlues of orreltion eqution for Clr (996 006) Fig d: Comprison etween the mesured nd predited vlues of orreltion eqution for Clr (996 006) Fig e: Comprison etween the mesured nd predited vlues of orreltion eqution for Clr (996-006) Eqution. (Clr) 5 7 9 Months of the Yer Eqution.5 (Clr) 5 7 9 Months of the Yer Mesured S Eqution. (Clr) 5 7 9 Months of the Yer Mesured S Mesured S Fig f: Comprison etween the mesured nd predited vlues of orreltion eqution for Clr (996 006) Fig g: Comprison etween the mesured nd predited vlues of orreltion eqution for Clr (996 006) COCLUSIO Eqution.8 (Clr) 5 7 9 Months of the Yer Eqution. (Clr) 5 7 9 Months of the Yer Mesured S Mesured S he proposed model equtions for Onitsh nd Clr n e used for estimting the monthly verge solr rdition on horizontl surfe for ny lotion in the ountry with solute vlues of the MPE less thn 8%

7 Communitions in Applied Sienes nd other ples outside the ountry whih hve the sme vlues of the mximum lerness index. herefore, from the ove results, it n e onluded tht the following simple first order Angstrom type orreltion n e used for the estimtion of glol solr rdition H on horizontl surfe t the lotions under study rn. So. S. Austr. Vol. 6, pp. - 8 Smo A. S, (985). Solr rdition with meteorologil dt.igeri Journl of Solr Energy Vol.. pp 59-6 EFEECES Angstrom, A.S. 9 Solr nd terrestril rdition meteorologil soiety. Vol. 50. Pp. -7 Ago, S., F.I. Ezem, nd P.E. ugwoke 007 Solr dition estimtes from eltive humidity-bsed D-model igeri Journl of solr Energy. vol. 8,pp. -8 Attili, P. nd Adll, A. (99).Sttistil omprison of glol nd diffuse solr rdition orreltion. Meteorologil Soiety, Vol 50, pp.. Iql, M. 98 An Introdution to solr rdition Ademy press.ew York. Murry, B. 96 Correltion of solr rdition with louds Solr Energy, Vol. () pp. 07. Presott, J.A 90 Evportion from wter surfe in reltion to solr rdition