Study of the viability of a glider drone for the return of experiments carried by weather balloons

Transcription

1 Bachelor s Degree in Aerospace Vehicle Engineering Study of the viability of a glider drone for the return of experiments carried by weather balloons Appendices Author: Albert Gassol Baliarda Director: Manel Soria Guerrero Co-Director: Josep Oriol Lizandra Dalmases June 12 th, 2015

2

3 Contents List of Figures List of Tables ii iii A International Standard Atmosphere 1 A.1 Temperature modeling A.2 Pressure modeling A.2.1 Pressure above the tropopause A.3 Density modeling A.4 Computational implementation and tabulated data B Jet streams study 7 B.1 Definition of the study B.2 Results B.2.1 Warm months B.2.2 Cold months B.3 Conclusions B.4 Numerical data extracted from the jet streams maps C Airfoil HQ 2.5/12 data 21 Bibliography 25 Albert Gasssol Baliarda i

4 LIST OF FIGURES List of Figures a.1 Temperature variation with altitude in the ISA a.2 Small atmosphere element b.1 Jet streams map corresponding to February 3, b.2 Relative frequencies of the jet stream speeds during the warm months 9 b.3 Relative frequencies of the jet stream speeds equals and majors than 60 kt during the warm months b.4 Relative frequencies of the jet stream speeds during the cold months 10 b.5 Relative frequencies of the jet stream speeds equals and majors than 60 kt during the cold months ii Albert Gasssol Baliarda

5 List of Tables a.1 International Standard Atmosphere, MSL conditions a.2 ISA tabulated data from the MSL up to m b.1 Results of the analysis of the warm months b.2 Global results of the analysis of the warm months b.3 Results of the analysis of the cold months b.4 Global results of the analysis of the cold months b.5 Numerical data of the wind speeds at 300 mbar over the Spanish territory c.1 Airfoil HQ 2.5/12 data corresponding to simulations at Re = Albert Gasssol Baliarda iii

6

7 Appendix A International Standard Atmosphere The International Standard Atmosphere (ISA) is a mathematical model employed to know the vertical distribution of air properties such as pressure, temperature, density and speed of sound. The ISA is published by the International Organization of Standardization (ISO) as an international standard, ISO 2533:1975. The ISA is also defined in the ICAO Document 7488/2. Since the real atmosphere does not remain constant at any particular time or place, a hypothetical model must be employed as an approximation to what may be expected, so in the ISA the air is considered as a perfect gas and assumed to be devoid of dust, moisture and water vapour. It is considered to be at rest with respect to the Earth as well (that is, no winds or turbulence). The ISA considers the mean sea level (MSL) conditions as given in Table a.1. Temperature T o = 288, 15 K Pressure P o = Pa Density ρ o = 1, 225 kg/m 3 Speed of sound a o = 340, 294 m/s Acceleration of gravity g o = 9, 807 m/s 2 Table a.1: International Standard Atmosphere, MSL conditions. A.1 Temperature modeling Figure a.1 represents the temperature variation with altitude in the ISA. As it is seen, temperature decreases with altitude at a constant rate of 6.5 K/1000 m up to the tropopause, placed at an altitude of m. Therefore, the air temperature within the troposphere is given by this expression: T = T o 6.5 h 1000 (a.1) Albert Gasssol Baliarda 1

8 Appendix A Figure a.1: Temperature variation with altitude in the ISA [1]. From the troposphere up to m the temperature remains at a constant value of K. A.2 Pressure modeling Pressure variations for the ISA are calculated using the hydrostatic equation, the perfect gas law and Equation (a.1). Figure a.2 shows a small atmosphere element and the forces acting on it. Figure a.2: Small atmosphere element [1]. Assuming that the atmosphere is in equilibrium, if one applies the first Newton s law it results: dp = ρgdh The equation of state for a perfect gas is: (a.2) P = ρrt (a.3) 2 Albert Gasssol Baliarda

9 where R = 287 J/kgK is the real gas constant for the air. Introducing (a.3) in (a.2) gives: dp P = g RT dh (a.4) Then, the relationship amongst pressure and altitude in the troposphere is obtained by integrating Equation (a.4) between h o = 0 and h: P P o dp h P = g h o RT dh = g h dh R h o T o h P = P o ( h T o ) (a.5) A.2.1 Pressure above the tropopause Temperature is constant at altitudes above the tropopause up to m, so that integrating Equation (a.4) from the tropopause to an altitude above it results in: P P 11 dp P = g RT 11 h h 11 dh P = P 11 e g RT 11 (h 11 h) (a.6) where the parameters with the subscript 11 correspond to the values at the tropopause, so T 11 = K, P 11 = Pa and h 11 = m. A.3 Density modeling Since the pressure and temperature are known for a given altitude, the density is easily calculated from the perfect gas equation of state, Equation (a.3): ρ = P RT (a.7) A.4 Computational implementation and tabulated data The ISA model can be implemented as a MATLAB function, which with the altitude as the only input, returns the atmospheric parameters T, a, P and ρ corresponding to this altitude. Following it is shown the code of this MATLAB function, which is used for the model developed in Chapter 4. 1 function [ T, a, P, rho ] = ISA ( h ) 2 % ISA r e t u r n s t h e temperature (T), speed o f sound ( a ), p r e s s u r e (P) and 3 % d e n s i t y ( rho ) o f t h e Earth s atmosphere corresponding to a g i v e n 4 % a l t i t u d e, according to t h e I n t e r n a t i o n a l Standard Atmosphere ( ISA ) 5 % model. A l l t h e magnitudes are e x p r e s s e d in SI u n i t s. 6 % The f u n c t i o n i s v a l i d f o r a l t i t u d e s between 0 m and m. For n e g a t i v e Albert Gasssol Baliarda 3

10 Appendix A 7 % a l t i t u d e s, t h e f u n c t i o n r e t u r n s t h e MSL c o n d i t i o n s. For a l t i t u d e s majors 8 % than m t h e f u n c t i o n might r e t u r n i n c o r r e c t v a l u e s % MSL Conditions 11 T0 = ; % [K] 12 a0 = ; % [m/ s ] 13 P0 = ; % [ Pa ] 14 rho0 = ; % [ kg /mˆ3] % Constants 17 g = ; % Gravity a c c e l e r a t i o n, [m/ s ˆ2] 18 Ma = e 3; % Molar mass o f air, [ kg /mol ] 19 R = ; % U n i v e r s a l gas constant, [ J /( mol*k) ] 20 Ra = R/Ma; % S p e c i f i c gas c o n s t a n t o f a i r [ J /( kg *K) ] 21 gam = 1. 4 ; % A d i a b a t i c gas c o n s t a n t o f t h e a i r 22 LR = 6.5e 3; % Lapse rate, [K/m] 23 LR2 = 1e 3; % Lapse rate, [K/m] switch l o g i c a l ( t r u e ) 26 c a s e h <= 0 27 T = T0 ; % [K] 28 a = a0 ; % [m/ s ] 29 P = P0 ; % [ Pa ] 30 rho = rho0 ; % [ kg /mˆ3] c a s e h > 0 && h <= T = T0+LR*h ; % [K] 34 a = sqrt (gam*ra*t) ; % [m/ s ] 35 P = P0*(1+LR*h/T0)ˆ( g/lr/ra ) ; % [ Pa ] 36 rho = P/Ra/T; % [ kg /mˆ3] c a s e h > && h <= T11 = T0+LR*11 e3 ; % [K] 40 P11 = P0*(1+LR*11 e3 /T0)ˆ( g/lr/ra ) ; % [ Pa ] T = T11 ; % [K] 43 a = sqrt (gam*ra*t) ; % [m/ s ] 44 P = P11*exp( g/ra/t11 *(h 11e3 ) ) ; % [ Pa ] 45 rho = P/Ra/T; % [ kg /mˆ3] o t h e r w i s e 48 T11 = T0+LR*11 e3 ; % [K] 49 T20 = T11 ; % [K] 50 P11 = P0*(1+LR*11 e3 /T0)ˆ( g/lr/ra ) ; % [ Pa ] 51 P20 = P11*exp( g/ra/t11 *(20 e3 11e3 ) ) ; % [ Pa ] T = T20+LR2*(h 20e3 ) ; % [K] 54 a = sqrt (gam*ra*t) ; % [m/ s ] 55 P = P20*(1+LR2*(h 20e3 )/ T20)ˆ( g/lr2/ra ) ; % [ Pa ] 56 rho = P/Ra/T; % [ kg /mˆ3] 4 Albert Gasssol Baliarda

11 57 58 end end In Table a.2 it has been tabulated the ISA parameters from the MSL up to an altitude of m. Altitude (m) Altitude (ft) Temperature (K) Pressure (Pa) Density (kg/m 3 ) θ = T/T o δ = P/P o σ = ρ/ρ o , , ,42 284, ,167 0,9887 0,9421 0, ,84 281, ,112 0,9774 0,8870 0, ,26 278, ,058 0,9662 0,8345 0, ,68 275, ,0063 0,9549 0,7846 0, ,10 271, ,957 0,9436 0,7371 0, ,52 268, ,909 0,9323 0,6919 0, ,94 265, ,863 0,9210 0,6490 0, ,36 262, ,819 0,9098 0,6083 0, ,78 258, ,777 0,8985 0,5697 0, ,20 255, ,736 0,8872 0,5331 0, ,62 252, ,697 0,8759 0,4985 0, ,04 249, ,660 0,8647 0,4656 0, ,46 245, ,624 0,8534 0,4346 0, ,88 242, ,590 0,8421 0,4052 0, ,30 239, ,557 0,8308 0,3775 0, ,72 236, ,525 0,8195 0,3513 0, ,14 232, ,495 0,8083 0,3267 0, ,56 229, ,466 0,7970 0,3034 0, ,98 226, ,439 0,7857 0,2815 0, ,40 223, ,413 0,7744 0,2609 0, ,82 219, ,388 0,7631 0,2415 0, ,24 216, ,364 0,7519 0,2234 0, ,66 216, ,336 0,7519 0,2064 0, ,08 216, ,311 0,7519 0,1908 0, ,50 216, ,287 0,7519 0,1763 0, ,92 216, ,265 0,7519 0,1629 0, ,34 216, ,245 0,7519 0,1506 0, ,76 216, ,227 0,7519 0,1392 0, ,18 216, ,210 0,7519 0,1286 0, ,60 216, ,194 0,7519 0,1189 0, ,02 216, ,179 0,7519 0,1099 0, ,44 216, ,165 0,7519 0,1015 0, ,86 216, ,153 0,7519 0,0938 0, ,28 216, ,141 0,7519 0,0867 0, ,70 216, ,131 0,7519 0,0801 0, ,12 216, ,121 0,7519 0,0741 0, ,54 216, ,112 0,7519 0,0685 0, ,96 216, ,103 0,7519 0,0633 0, ,38 216, ,095 0,7519 0,0585 0, ,80 216, ,088 0,7519 0,0540 0,0719 Table a.2: ISA tabulated data from the MSL up to m. Albert Gasssol Baliarda 5

12

13 Appendix B Jet streams study Jet streams are narrow fast flowing air currents, located at an altitude near the tropopause. There are mainly two jet streams in the Earth that go over the planet from West to East: the polar jet and the subtropical jet. Europe (and Spain) is mostly affected be the polar jet, which is usually located between latitudes 50 N and 60 N at altitudes around 7 12 km, and can reach speeds greater than 100 km/h. Due to the high speeds of the air at these altitudes, it is necessary to have an idea of the order of these air speeds to see the conditions at which the drone will be subjected. B.1 Definition of the study This study consists on a statistical analysis of the jet streams maps to determine how jet streams usually are in the regions of interest. To do this, it has been employed the California Regional Weather Server 1 resource from the San Francisco State University (SFSU), where there is available, amongst others, an archive with the jet streams maps of several regions of the planet and some other meteorological data of the last four years. These jet streams maps are elaborated using the Global Forecast System 2 (GFS) model and are always calculated at the altitude corresponding to 300 mbar, this means, about 9150 m. Figure b.1 shows how these maps look. Following it is explained the considerations that it have been done and how the analysis has been carried out: ˆ The analysis is done since the 1 st of March of 2014 till the 28 th of February of 2015, this means, it covers a period of one year. ˆ The data that is considered in the study belongs to the Spanish territory, concretely the north-eastern zone. 1 the resource can be found in [2]. 2 Global numerical weather prediction system containing a global computer model and variational analysis run by the US National Weather Service (NWS). Albert Gasssol Baliarda 7

14 Appendix B Figure b.1: Jet streams map corresponding to February 3 rd, Extracted from [2]. ˆ The maps from the California Regional Weather Server display the overall velocity field with vector arrows, but only gives a numerical value to the velocities exceeding 60 kt. Due to this fact, it is treated apart the days at which wind speeds majors than 60 kt were detected from the days at which wind speeds were less than 60 kt (values that can be considered as usual values). ˆ Initially, it is considered apart the data belonging to warm months (April, May, June, July, August and September) from the belonging to cold months (March, October, November, December, January and February) to see if there exists or not a significant difference amongst them. ˆ The data was extracted from the maps twice daily, at 12am and 12pm local time. It was done to see if the sun affects or not in the jet stream intensity. B.2 Results In this section there are shown the results of the numerical data extracted from the jet streams maps. The numerical data can be seen in Table b.5 of Section B.4. 8 Albert Gasssol Baliarda

15 B.2.1 Warm months Table b.1 and Figure b.2 show the results obtained from the analysis of the previously called warm months (April, May, June, July, August and September), separating the data extracted at 12am and 12pm. 12am 12pm % Speeds < 60 kt 55,74% 50,82% % Speeds > 60 kt 44,26% 49,18% Average > 60 kt(kt) 70,62 70,78 Standard deviation (kt) 11,66 12,20 Table b.1: Results of the analysis of the warm months. Figure b.2: Relative frequencies of the jet stream speeds during the warm months. It does not exist any clear difference between the values taken at 12am and the ones taken at 12pm, so it will be better considered together. Table b.2 gathers all the results, and Figure b.3 considers only the registered wind speeds equals and majors than 60 knots and represents their relative frequencies. Total % speeds < 60 kt 53,28% % speeds > 60 kt 46,72% Average > 60 kt(kt) 70,7 Standard deviation (kt) 11,91 Table b.2: Global results of the analysis of the warm months. Albert Gasssol Baliarda 9

16 Appendix B Figure b.3: Relative frequencies of the jet stream speeds equals and majors than 60 kt during the warm months. B.2.2 Cold months In the same way that it has been done previously, Table b.3 and Figure b.4 show the results of the analysis of the cold months (March, October, November, December, January and February). 12am 12pm % Speeds < 60 kt 51,10% 51,10% % Speeds > 60 kt 48,90% 48,90% Average > 60 kt(kt) 77,08 75,39 Standard deviation (kt) 19,84 18,03 Table b.3: Results of the analysis of the cold months. Figure b.4: Relative frequencies of the jet stream speeds during the cold months. 10 Albert Gasssol Baliarda

17 As it occurs in the warm months, it does not exist a significant difference between the values from 12am and 12pm, so that it will be assumed that it does not exist any dependence with the hour of the day in any of the two cases. According to this, Table b.4 gathers the results in a single table and Figure b.5 represents the relative frequencies of the registered wind speeds equals and majors than 60 kt, as Figure b.3 does. Total % Speeds < 60 kt 51,10% % Speeds > 60 kt 48,90% Average > 60 kt(kt) 76,24 Standard deviation (kt) 18,92 Table b.4: Global results of the analysis of the cold months. Figure b.5: Relative frequencies of the jet stream speeds equals and majors than 60 kt during the cold months. B.3 Conclusions According to the results, wind speeds only exceed the value of 60 knots a half of the times roughly, and when it occurs, the great majority of the registered wind speeds are between 60 kt and 70 kt. Extreme values like kt only occur in isolated days and quite randomly, so these values should not be considered as representatives. It is also remarkable that the average of the wind speeds majors than 60 kt is slightly major (5.54%) during the cold months than during the warm months, but the difference is not enough to be worth being considered. Albert Gasssol Baliarda 11

18 Appendix B B.4 Numerical data extracted from the jet streams maps Table b.5 shows the numerical data extracted from the California Regional Weather Server jet streams maps. Cells containing a hyphen mean that the registered wind speeds were less than 60 kt. Day 12am Velocity (knots) 12pm March April Albert Gasssol Baliarda

19 May Albert Gasssol Baliarda 13

20 Appendix B June July Albert Gasssol Baliarda

21 August Albert Gasssol Baliarda 15

22 Appendix B September October Albert Gasssol Baliarda

23 November Albert Gasssol Baliarda 17

24 Appendix B December Albert Gasssol Baliarda

25 January February Albert Gasssol Baliarda 19

26 Appendix B Table b.5: Numerical data of the wind speeds at 300 mbar over the Spanish territory extracted from the California Regional Weather Server maps. 20 Albert Gasssol Baliarda

27 Appendix C Airfoil HQ 2.5/12 data The following table (Table c.1) contains the numerical information of the aerodynamic airfoil HQ 2.5/12 (available in [3]) from which have been done the airfoil graphs of Chapter 3, Section 3.2. α (deg) C l C d C dp C m Top Xtr Bot Xtr -10,75-0,4130 0,0969 0,0946-0, , ,5-0,5896 0,0503 0,0479-0, , ,25-0,6196 0,0453 0,0428-0, , ,6528 0,0432 0,0407-0, ,0117-9,75-0,6799 0,0369 0,0338-0,0675 0,9950 0,0115-9,5-0,6719 0,0298 0,0261-0,0717 0,9869 0,0116-9,25-0,6471 0,0262 0,0219-0,0745 0,9832 0, ,6216 0,0236 0,0190-0,0760 0,9780 0,0123-8,75-0,5920 0,0216 0,0168-0,0777 0,9742 0,0126-8,5-0,5585 0,0204 0,0153-0,0796 0,9714 0,0130-8,25-0,5321 0,0180 0,0127-0,0808 0,9665 0, ,5029 0,0172 0,0118-0,0816 0,9604 0,0150-7,75-0,4719 0,0163 0,0107-0,0827 0,9559 0,0160-7,5-0,4468 0,0155 0,0098-0,0824 0,9475 0,0168-7,25-0,4200 0,0145 0,0087-0,0825 0,9412 0, ,3973 0,0138 0,0079-0,0817 0,9320 0,0199-6,75-0,3707 0,0133 0,0074-0,0815 0,9254 0,0221-6,5-0,3472 0,0127 0,0068-0,0808 0,9168 0,0257-6,25-0,3204 0,0125 0,0066-0,0806 0,9101 0, ,2963 0,0120 0,0060-0,0800 0,9020 0,0357-5,75-0,2694 0,0119 0,0058-0,0798 0,8955 0,0402-5,5-0,2448 0,0114 0,0053-0,0793 0,8876 0,0453-5,25-0,2181 0,0112 0,0050-0,0791 0,8813 0, ,1915 0,0111 0,0048-0,0788 0,8739 0,0531-4,75-0,1661 0,0106 0,0043-0,0784 0,8678 0,0601-4,5-0,1394 0,0105 0,0041-0,0782 0,8605 0,0653-4,25-0,1134 0,0101 0,0037-0,0779 0,8544 0, ,0871 0,0098 0,0034-0,0776 0,8474 0,0867 Albert Gasssol Baliarda 21

28 Appendix C -3,75-0,0608 0,0094 0,0031-0,0774 0,8411 0,1122-3,5-0,0356 0,0089 0,0029-0,0772 0,8345 0,1845-3,25-0,0104 0,0083 0,0026-0,0770 0,8279 0, ,0138 0,0077 0,0024-0,0768 0,8216 0,4196-2,75 0,0395 0,0073 0,0024-0,0764 0,8146 0,5077-2,5 0,0665 0,0072 0,0024-0,0761 0,8085 0,5551-2,25 0,0938 0,0072 0,0023-0,0759 0,8012 0,5841-1,75 0,1490 0,0072 0,0023-0,0755 0,7882 0,6298-1,5 0,1770 0,0072 0,0023-0,0754 0,7823 0,6450-1,25 0,2047 0,0072 0,0023-0,0752 0,7757 0, ,2320 0,0072 0,0023-0,0749 0,7685 0,6794-0,75 0,2593 0,0072 0,0023-0,0746 0,7603 0,6937-0,5 0,2869 0,0072 0,0023-0,0744 0,7525 0,7037-0,25 0,3143 0,0072 0,0023-0,0742 0,7436 0, ,3420 0,0073 0,0023-0,0740 0,7355 0,7226 0,25 0,3693 0,0073 0,0023-0,0737 0,7265 0,7321 0,5 0,3969 0,0073 0,0023-0,0736 0,7190 0,7405 0,75 0,4246 0,0073 0,0023-0,0735 0,7110 0, ,4523 0,0073 0,0023-0,0734 0,7028 0,7544 1,25 0,4797 0,0073 0,0023-0,0732 0,6943 0,7616 1,5 0,5072 0,0074 0,0023-0,0731 0,6848 0,7687 1,75 0,5345 0,0074 0,0023-0,0729 0,6754 0, ,5616 0,0074 0,0024-0,0727 0,6655 0,7838 2,25 0,5887 0,0074 0,0024-0,0725 0,6546 0,7920 2,5 0,6155 0,0075 0,0024-0,0722 0,6433 0, ,6683 0,0075 0,0026-0,0715 0,6168 0,8184 3,25 0,6942 0,0076 0,0026-0,0710 0,6004 0,8285 3,5 0,7197 0,0077 0,0027-0,0705 0,5819 0,8399 3,75 0,7444 0,0078 0,0028-0,0698 0,5613 0, ,7684 0,0079 0,0029-0,0690 0,5380 0,8681 4,25 0,7909 0,0080 0,0030-0,0679 0,5111 0,8899 4,5 0,8149 0,0081 0,0031-0,0670 0,4785 0,9327 4,75 0,8505 0,0085 0,0033-0,0690 0, ,8731 0,0089 0,0036-0,0683 0, ,25 0,8953 0,0094 0,0038-0,0675 0, ,5 0,9179 0,0098 0,0041-0,0668 0, ,75 0,9402 0,0102 0,0044-0,0660 0, ,9632 0,0106 0,0047-0,0654 0, ,25 0,9863 0,0110 0,0051-0,0647 0, ,5 1,0095 0,0114 0,0054-0,0640 0, ,75 1,0325 0,0117 0,0057-0,0633 0, ,0549 0,0121 0,0060-0,0626 0, ,25 1,0768 0,0125 0,0063-0,0617 0, ,5 1,0988 0,0129 0,0067-0,0609 0, ,75 1,1201 0,0133 0,0071-0,0600 0, ,1405 0,0138 0,0075-0,0590 0, Albert Gasssol Baliarda

29 8,25 1,1605 0,0143 0,0079-0,0579 0, ,5 1,1800 0,0148 0,0083-0,0567 0, ,75 1,1983 0,0153 0,0088-0,0553 0, ,2157 0,0159 0,0093-0,0539 0, ,25 1,2314 0,0165 0,0099-0,0521 0, ,5 1,2457 0,0171 0,0104-0,0501 0, ,75 1,2580 0,0178 0,0111-0,0478 0, ,2633 0,0190 0,0121-0,0447 0, ,25 1,2643 0,0205 0,0134-0,0412 0, ,5 1,2730 0,0216 0,0145-0,0388 0, ,75 1,2842 0,0225 0,0155-0,0369 0, ,2940 0,0235 0,0167-0,0349 0, ,25 1,3020 0,0247 0,0179-0,0328 0, ,5 1,3058 0,0263 0,0196-0,0305 0, ,75 1,3135 0,0276 0,0210-0,0288 0, ,3203 0,0290 0,0225-0,0271 0, ,25 1,3252 0,0306 0,0242-0,0255 0, ,5 1,3287 0,0324 0,0261-0,0239 0, ,75 1,3316 0,0343 0,0281-0,0226 0, ,3328 0,0364 0,0304-0,0213 0, ,25 1,3331 0,0388 0,0328-0,0203 0, ,5 1,3324 0,0413 0,0355-0,0195 0, ,75 1,3304 0,0441 0,0384-0,0189 0, ,3273 0,0471 0,0415-0,0185 0, ,25 1,3227 0,0504 0,0449-0,0184 0, ,5 1,3185 0,0539 0,0485-0,0186 0, ,75 1,3125 0,0577 0,0525-0,0191 0, ,3068 0,0617 0,0566-0,0199 0, ,25 1,2968 0,0665 0,0616-0,0212 0, ,5 1,2913 0,0709 0,0661-0,0225 0, ,75 1,2815 0,0761 0,0713-0,0243 0, ,2772 0,0806 0,0760-0,0260 0, ,25 1,2667 0,0862 0,0817-0,0282 0, ,5 1,2612 0,0912 0,0869-0,0304 0, Table c.1: Airfoil HQ 2.5/12 data corresponding to simulations at Re = Extracted from [3]. Albert Gasssol Baliarda 23

30

31 Bibliography [1] Mustafa Cavcar. The International Standard Atmosphere (ISA). In: Anadolu University, Turkey (2000), pp [2] San Francisco State University. California Regional Weather Server url: http: //squall.sfsu.edu/crws.html (visited on 03/16/2015). [3] HQ 2.5/12 AIRFOIL. url: NoTithttp://airfoiltools.com/airfoil/details? airfoil=hq2512-ille (visited on 04/30/2015). Albert Gasssol Baliarda 25