Computer Programs

Computer Programs This thesis presents a novel package of computer programs based on Matlab software. This package presents a complete design and control strategy of interconnecting PV system, WES and PV/WES HEPS accompanied with or without battery storage, BS with utility grid, UG. It can be applied for design any hybrid electric power system consists of PV, WES or PV/WES to feed a load in any site in the world as stand-alone or interconnected with UG. This package includes:1 Design of PV System To design PV system, set the variable of penetration ration in the program equal to 1. The format of this variable as shown pent=1. The flowchart of this program is shown in Fig. 2-7 see chapter 2. The input data of this program is: 1- Site latitude, Degrees. Zafarâna site is located on the western coast of the Seuz Gulf, latitude 29.07o N and longitude 31.36o E, Egypt. Input Site latitude in degrees as shown below under the variable:phia=29.07*pi/180 2- Hourly radiation, kW/m2. The data can be obtained from the Egyptian Metrological Authority to Zafarâna site or any site in the world from 7 A.M to 16 P.M. Format of the input data of radiation for twelve month as shown below:a1=[0,0,0,0,0,0,0.003,0.242,0.785,1.396,1.712,2.,2.1,2.,1.772,1.158,0.384,0.004,0,0,0,0,0,0]; a2=[0,0,0,0,0,0,0.138,0.675,1.02,1.403,2.116,2.622,2.553,2.591,2.561,1.886,1.196,0.406,0.008,0,0,0,0,0]; a3=[0,0,0,0,0,0,0.069,0.797,1.687,2.285,2.944,3.143,3.113,3.044,1.963,1.449,0.59,0.113,0,0,0,0,0,0]; a4=[0,0,0,0,0,0,0.123,0.399,0.659,1.066,2.131,2.545,2.959,3.097,2.568,1.894,0.92,0.118,0,0,0,0,0,0]; a5=[0,0,0,0,0,0,0.138,0.805,1.503,1.671,2.062,2.806,3.381,3.251,3.051,2.415,2.131,1.395,0.537,0.054,0,0,0] a6=[0,0,0,0,0,0,0,0.184,0.843,1.671,2.492,3.021,3.281,3.358,3.373,3.182,2.752,2.17,1.449,0.667,0.107,0,0,0] a7=[0,0,0,0,0,0,0.038,0.337,0.989,2.024,2.852,3.373,3.634,3.626,3.35,2.791,2.024,1.357,0.652,0.0115,0,0,0] a8=[0,0,0,0,0,0,0.038,0.56,1.319,2.101,2.599,3.22,3.45,3.45,3.12,2.775,2.139,1.311,0.498,0.031,0,0,0,0];

234

a9=[0,0,0,0,0,0,0.008,0.36,1.173,1.963,2.461,3.044,3.373,3.289,2.975,2.415,1.687,0.859,0.119,0,0,0,0,0]; a10=[0,0,0,0,0,0,0.176,0.621,1.058,1.51,2.461,2.07,2.17,2.07,1.579,0.82,0.406,0.015,0,0,0,0,0,0]; a11=[0,0,0,0,0,0,0.038,0.498,1.119,1.802,2.246,2.384,2.269,1.986,1.441,0.705,0.115,0,0,0,0,0,0,0]; a12=[0,0,0,0,0,0,0,0.222,0.966,1.434,1.518,1.449,1.664,1.035,0.36,0.238,0.054,0,0,0,0,0,0,0];

3- Hourly load demand, kW. It is assumed here that the load demand varies monthly. This means that each month has daily load curve different from other months. Therefore, there are twelve daily load curves through the year. Format of the input data of hourly load demand for twelve month as shown below:pl11=[30,34,37,41,46,57,62,70,74,68,61,53,42,31,28,26,26,24,26,27,28,28,29,30; 32,35,38,43,48,56,64,72,78,71,63,55,44,38,31,28,26,25,27,27,29,30,30,31; 33,37,40,44,47,57,66,74,79,81,78,65,54,41,36,32,27,26,28,28,29,30,32,33; 32,36,39,41,45,51,68,76,81,86,80,73,64,53,46,37,30,27,26,25,27,31,31,33; 34,38,40,45,51,58,67,75,83,87,85,78,68,60,51,43,37,32,28,29,30,30,32,36; 37,42,51,58,62,71,78,83,87,90,84,79,71,68,62,51,47,42,36,32,28,33,35,38; 40,43,54,61,68,75,81,84,91,93,87,81,76,71,64,56,51,48,39,32,30,33,37,40; 41,43,56,64,67,73,82,87,93,95,90,86,81,74,70,63,54,47,41,35,32,34,37,38; 38,41,48,62,63,68,78,81,87,90,84,78,69,63,58,54,42,38,35,34,36,37,37,34; 35,38,45,51,58,63,69,74,83,89,82,74,65,58,55,48,41,32,30,28,29,31,32,32; 31,36,41,50,57,61,65,71,79,84,80,71,63,56,49,44,37,28,28,29,30,30,31,33 ; 31,34,36,42,48,56,61,70,75,69,60,54,47,40,32,27,28,28,27,28,30,31,31,30 ]*1e6;

4- Ambient Temperature The fourth input of the program is the ambient temperature to the site. Format of the input data of ambient temperature as shown below:Tmin=[9.7,12.8,14,27.8,19.5,24.9,27.6,25.8,24.9,23.2,17.8,14.8]; Tmax=[18.9,21,21.9,40.8,27.9,35.4,39,35.4,37.3,30.2,25.3,21.6]; At the command prompt denoted by the symbol >> Type the name solarwindphd8battery53. Output Files The output file with name subsys.txt contains the data in the Table (2-3) and (2-5). 235

Data in Table (2-6) will display if you type surdef'/1e6 after the symbol >>. Data in Table (2-7) will display if you type [Egenh1'/1e6,Eyload'/1e6] after the symbol >> 1-2 Design of the DC/DC Boost Converter Input data for this program as shown below: vo=600;po=500e3;fs=50e3; Type boostMohan1 after the symbol >> The output figure of this program is shown in Figs. 2-13 to 2-16. 2- Operation Control Strategy and Simulation of PV System Connected With UG The Matlab simulink Program novlegridpv600.mdl is used to modeling and simulation of PV system connected with EU. 3- Design of WES To design WES, set penetration ration in the program equal to 0. The format of this variable as shown pent=0. Input wind speed. Format of the input data of wind speed as shown below:% u is the wind speed uo1=[6,6,5,3,0,0,0,10,8,12,13,12,10,8,11,12,13,10,8,6,6,7,5,5]; uo2=[12,6,5,4,3,0,0,0,0,3,5,11,15,14,12,14,13,12,15,16,14,15,15,14]; uo3=[12,10,12,10,10,11,8,10,14,12,12,13,14,15,13,12,10,11,12,14,16,15,17,16]; uo4=[18,17,20,18,19,17,12,13,8,12,15,14,18,20,23,22,19,20,16,15,16,14,18,20]; uo5=[11,13,7,2,5,7,6,3,0,7,12,8,8,10,8,8,5,6,8,10,10,8,9,10]; uo6=[12,10,10,11,10,12,10,10,11,12,11,11,12,12,12,10,12,16,14,15,14,14,14,16]; uo7=[10,8,3,5,8,10,9,10,8,13,9,8,7,10,12,10,10,10,13,14,8,6,7,8]; uo8=[10,8,2,0,0,5,6,7,3,4,4,5,6,8,8,7,8,7,6,10,13,12,13,11]; uo9=[9,10,10,8,9,7,8,7,10,10,10,8,9,10,7,13,12,12,13,15,11,10,8,9]; uo10=[6,7,0,0,0,0,0,0,0,10,12,11,14,10,8,9,8,8,13,10,12,10,8,9]; uo11=[5,8,7,9,6,8,5,7,8,12,11,13,10,12,15,13,10,11,13,10,8,9,7,8]; uo12=[12,15,12,11,10,12,8,6,10,8,6,5,10,0,0,0,0,5,10,12,13,10,8,5];

236

3- Design of PV/ WES To design PV/WES set penetration ration in the program equal to 0,.1,.2,…….1. 4- Operation Control Strategy and Simulation of WES System Connected With UG The Matlab simulink Program novlegridwtg600.mdl is used to modeling and simulation of WTG connected with EU. 5- Application And Results Of FL Configuration of PV/UG The program phdpv.m is used to estimate LOLP for this configuration. Configuration of WES/UG The program phdwind.m is used to estimate LOLP for this configuration. Configuration of PV/WES/UG The program phdpvwes.m is used to estimate LOLP for this configuration.

237