diff options
Diffstat (limited to 'sample_notebooks/amit kumarsaini/Chapter2.ipynb')
| -rwxr-xr-x | sample_notebooks/amit kumarsaini/Chapter2.ipynb | 202 |
1 files changed, 202 insertions, 0 deletions
diff --git a/sample_notebooks/amit kumarsaini/Chapter2.ipynb b/sample_notebooks/amit kumarsaini/Chapter2.ipynb new file mode 100755 index 00000000..9a5cce62 --- /dev/null +++ b/sample_notebooks/amit kumarsaini/Chapter2.ipynb @@ -0,0 +1,202 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:ddaa68fa50b83699f8a6f35deb059bff01c78a46802bd57850513a730bade48f"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter2-Solar Radiation and its Measurement"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.1-pg 59"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##Ex2.4.1.;Detremine local solar time and declination\n",
+ "import math\n",
+ "##The local solar time=IST-4(standard time longitude-longitude of location)+Equation of time correstion\n",
+ "##IST=12h 30min;for the purpose of calculation we are writing it as a=12h,b=29 min 60sec;\n",
+ "a=12.;\n",
+ "b=29.60;\n",
+ "##(standard time longitude-longitude of location)=82 degree 30min - 77 degree 30min;\n",
+ "##for the purpose of calculation we are writing it as\n",
+ "STL3=82.5-72.5;\n",
+ "##Equation of time correstion: 1 min 01 sec\n",
+ "##for the purpose of calculation we are writing it as\n",
+ "c=1.01;\n",
+ "##The local solar time=IST-4(standard time longitude-longitude of location)+Equation of time correstion\n",
+ "LST=b-STL3-c;\n",
+ "print'%s %.2f %s %.2f %s'%(\" The local solar time=\",a,\"\"and \"\",LST,\" in hr.min.sec\");\n",
+ "##Declination delta can be obtain by cooper's eqn : delta=23.45*sin((360/365)*(284+n))\n",
+ "n=170.;##(on June 19)\n",
+ "##let\n",
+ "a=(360./365.)*(284.+n)\n",
+ "aa=(a*math.pi)/180.\n",
+ "##therefore\n",
+ "delta=23.45*math.sin(aa);\n",
+ "print'%s %.2f %s'%(\"\\n delta=\",delta,\" degree\");\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ " The local solar time= 12.00 18.59 in hr.min.sec\n",
+ "\n",
+ " delta= 23.43 degree\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.2-pg 59"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "##Ex2.4.2.;Calculate anglr made by beam radiation with the normal to a flat collector.\n",
+ "gama=0.;##since collector is pointing due south.\n",
+ "##For this case we have equation : cos_(theta_t)=cos(fie-s)*cos(delta)*cos(w)+sin(fie-s)*sin(delta)\n",
+ "##with the help of cooper eqn on december 1,\n",
+ "n=335.;\n",
+ "##let\n",
+ "a=(360./365.)*(284.+n);\n",
+ "aa=(a*math.pi)/180;\n",
+ "##therefore\n",
+ "delta=23.45*math.sin(aa);\n",
+ "print'%s %.2f %s'%(\" delta=\",delta,\" degree\");\n",
+ "##Hour angle w corresponding to 9.00 hour=45 Degreew\n",
+ "w=45.;##degree\n",
+ "##let\n",
+ "a=math.cos(((28.58*math.pi)/180.)-((38.58*math.pi)/180.))*math.cos(delta*math.pi*180**-1)*math.cos(w*math.pi*180**-1);\n",
+ "b=math.sin(delta*math.pi*180**-1)*math.sin(((28.58*math.pi)/180.)-((38.58*math.pi)/180.));\n",
+ "##therefore\n",
+ "cos_of_theta_t=a+b;\n",
+ "theta_t=math.acos(cos_of_theta_t)*57.3;\n",
+ "print'%s %.2f %s'%(\"\\n theta_t=\",theta_t,\" Degree\");\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ " delta= -22.11 degree\n",
+ "\n",
+ " theta_t= 44.73 Degree\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.1-pg 68"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "##Ex2.7.1.;Determine the average values of radiation on a horizontal surface\n",
+ "import math\n",
+ "##Declination delta for June 22=23.5 degree, sunrice hour angle ws\n",
+ "delta=(23.5*math.pi)/180;##unit=radians\n",
+ "fie=(10*math.pi)/180;##unit=radians\n",
+ "##Sunrice hour angle ws=acosd(-tan(fie)*tan(delta))\n",
+ "ws=math.acos(-math.tan(fie)*math.tan(delta));\n",
+ "print'%s %.2f %s'%(\" Sunrice hour angle ws=\",ws,\" Degree\");\n",
+ "n=172.;##=days of the year (for June 22)\n",
+ "##We have the relation for Average insolation at the top of the atmosphere\n",
+ "##Ho=(24/%pi)*Isc*[{1+0.033*(360*n/365)}*((cos (fie)*cos(delta)*sin(ws))+(2*%pi*ws/360)*sin(fie)*sin(delta))]\n",
+ "Isc=1353.;##SI unit=W/m^2\n",
+ "ISC=1165.;##MKS unit=kcal/hr m^2\n",
+ "##let\n",
+ "a=24./math.pi;\n",
+ "aa=(360.*172.)/365.;\n",
+ "aaa=(aa*math.pi)/180.;\n",
+ "b=math.cos(aaa);\n",
+ "bb=0.033*b;\n",
+ "bbb=1+bb;\n",
+ "c=(10*math.pi)/180.;\n",
+ "c1=math.cos(c);\n",
+ "cc=(23.5*math.pi)/180;\n",
+ "cc1=math.cos(cc);\n",
+ "ccc=(94.39*math.pi)/180;\n",
+ "ccc1=math.sin(ccc);\n",
+ "c=c1*cc1*ccc1;\n",
+ "d=(2*math.pi*ws)/360.;\n",
+ "e=(10*math.pi)/180;\n",
+ "e1=math.sin(e);\n",
+ "ee=(23.5*math.pi)/180;\n",
+ "ee1=math.sin(ee);\n",
+ "e=e1*ee1;\n",
+ "##therefoe Ho in SI unit\n",
+ "Ho=a*Isc*(bbb*(c+(d*e)));\n",
+ "print'%s %.2f %s'%(\"\\n SI UNIT->Ho=:\",Ho,\" W/m^2\");\n",
+ "Hac=Ho*(0.3+(0.51*0.55))\n",
+ "print'%s %.2f %s'%(\"\\n SI UNIT->Hac=\",Hac,\" W/m^2 day\");\n",
+ "ho=a*ISC*(bbb*(c+(d*e)));\n",
+ "print'%s %.2f %s'%(\"\\n MKS UNIT->Ho=\",ho,\" kcal/m^2\");\n",
+ "hac=ho*0.58;\n",
+ "print'%s %.2f %s'%(\"\\n MKS UNIT->Hac=\",hac,\" kcal/m^2 day\");\n",
+ "\n",
+ "##The values are approximately same as in textbook\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ " Sunrice hour angle ws= 1.65 Degree\n",
+ "\n",
+ " SI UNIT->Ho=: 9025.25 W/m^2\n",
+ "\n",
+ " SI UNIT->Hac= 5239.16 W/m^2 day\n",
+ "\n",
+ " MKS UNIT->Ho= 7771.19 kcal/m^2\n",
+ "\n",
+ " MKS UNIT->Hac= 4507.29 kcal/m^2 day\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |