Supplementary information<\/a><\/h5>\n\n\n\n CSV Data: <\/strong><\/p>\n\n\n\n\n- Click hyperlink to access CSV Data: Fig2212_data<\/a> <\/li>\n\n\n\n
- Click hyperlink to access CSV Data: Change_in_EnthalpyVSDelta<\/a><\/li>\n\n\n\n
- Click hyperlink to access CSV Data: H2O_Equilibrium_Constant<\/a> <\/li>\n<\/ul>\n\n\n\n
Book: <\/strong><\/p>\n\n\n\n\n- Example Question 3 Solution: Solar Syngas Production_Problem 3<\/a> <\/li>\n\n\n\n
- Example Questions 4-7 Solutions: Book Chapter Problems_Gordon-v2<\/a><\/li>\n<\/ul>\n\n\n\n
Exercises Inspired by a Solar Energy Utilization Class<\/h5>\n\n\n\n
For your reference, below is a list of of the exercises in this series:<\/p>\n\n\n\n
\n- Exercise 1: Getting Familiar with Python: an Extraterrestrial Solar Radiation exercise. <\/a><\/li>\n\n\n\n
- Exercise 2: Solar Time and Solar Time Python Calculator. <\/a><\/li>\n<\/ul>\n\n\n\n
<\/p>\n\n\n\n
Exercise 1 Instructions<\/a><\/h5>\n\n\n\nExtraterrestrial Solar Radiation<\/strong><\/p>\n\n\n\nSolar radiation flux, or irradiance (G) in W\/m2<\/sup>, outside the earth\u2019s atmosphere and normal to earth\u2019s surface (Gon<\/sub>, where o stands for outside and n for normal direction) may be calculated as a function of the day in the year according to the following equation below.<\/p>\n\n\n\n![\"\"](\"https:\/\/faculty.eng.ufl.edu\/jonathan-scheffe\/wp-content\/uploads\/sites\/765\/2020\/07\/Exer1Equ.png\")
<\/figure>\n<\/div>\n\n\n<\/p>\n\n\n\n
Gsc<\/sub> is the solar constant and is equal to 1367 W\/m2 and n is the day in the year. The term contained within cosine is in units of degrees.<\/p>\n\n\n\nBrief Python Introduction<\/strong><\/p>\n\n\n\nPackages and Calling Functions<\/em>: In Python, there are several powerful computing packages commonly used that should be declared at the top of your .py, or python file. For example, some of the more popular computational packages useful for engineering include numpy, scipy, matplotlib and pandas, among others. These must be installed on your computer, either manually using a package installer (e.g. pip installer<\/a>), or are installed by default with many Python distributions, including the Anaconda Distribution that we will use in this exercise. For example, if we want to utilize the numpy package, we can use the import function and say:<\/p>\n\n\n\n\nimport numpy as np<\/p>\n<\/div>\n\n\n\n
which allows us to call the numpy package simply by typing np. If we wish to use numpy package to calculate the cosine of an angle, we would first need to call it by typing np followed by a period and then the function contained within the package we wish to utilize, in this case cos, and the value within it (in this case the numerical value pi which is also called using the numpy package as np.pi) contained within parenthesis. For example,<\/p>\n\n\n\n
\nnp.cos(np.pi)<\/p>\n<\/div>\n\n\n\n
The input contained within parentheses should have units radians and the output will be in units radians. For example, to correctly use the equation for Gon, the cosine term will look like the following.<\/p>\n\n\n\n
\nnp.cos(360*n\/365*np.pi\/180)<\/p>\n<\/div>\n\n\n\n
because we need to convert from degrees to radians by multiplying everything by np.pi\/180.<\/p>\n\n\n\n
Declaring Functions:<\/em> Functions are declared in Python by defining them using def followed by the function name, parentheses that include any variables the function is dependent on, followed by a colon. The body of the function is then started on a new, indented line and the value that is returned by the function must be declared at the end using return<\/em>. For example, if we wish to define a function y that depends on x and is a linear line with slope m and intercept y and later call it, we would write<\/p>\n\n\n\n\ndef y(x):<\/p>\n
return m*x+b<\/p>\n
y(x)<\/p>\n<\/div>\n\n\n\n
<\/p>\n\n\n\n
Exercise<\/strong>
Complete the following questions that deal with determining Gon<\/sub> using the above equation and a more accurate form of the equation (not shown, Gon_accurate function), as a way to get comfortable computing in Python using some basic techniques, learning to call and declare functions, plotting data using the package matplotlib, and using two different methods to evaluate python code, namely a traditional editor named Spyder and a web based notebook style editor called Jupyter.<\/p>\n\n\n\n\n- Download Python<\/a> (\u201cAnaconda Distribution\u201d). This includes many useful packages but most commonly we will use Spyder and Jupyter, two different ways to evaluate Python code.<\/li>\n\n\n\n
- Open and run \u201cExtraterrestrial Radiation Calculator.py\u201d using Spyder. It can be copied and pasted from the appendix of <\/span>this document<\/a><\/strong> and saved to your computer as a .py file. After running, you should see a plot with two curves.
![\"\"](\"https:\/\/faculty.eng.ufl.edu\/jonathan-scheffe\/wp-content\/uploads\/sites\/765\/2020\/07\/rb_dayyear-300x199.png\")
<\/p>
<\/li>\n\n\n\n
- Replot the data using the same code but change the plotted lines from red \u2018r\u2019 and blue \u2018b\u2019in the plt.plot lines of code, to other colors of your choice. Here are examples and more documentation<\/a>.
![\"\"](\"https:\/\/faculty.eng.ufl.edu\/jonathan-scheffe\/wp-content\/uploads\/sites\/765\/2020\/07\/DAyyear-300x179.png\")
<\/figure><\/li>\n\n\n\n- We have plotted data every day using the function np.arrange with the linedays=np.arange(1,365,1). This creates an array from 1 to 365 with a spacing of 1. However, let’s imagine we want to see individual data points plotted \u2013 it would be difficult to distinguish between them with so many. Therefore, lets plot a point every 10 days by changing the spacing from 1 to 10. Show the results using blue circles for the blue curve by inserting \u2018bo\u2019 in place of \u2018b\u2019.
![\"\"](\"https:\/\/faculty.eng.ufl.edu\/jonathan-scheffe\/wp-content\/uploads\/sites\/765\/2020\/07\/Dayyear2-300x192.png\")
<\/figure><\/li>\n\n\n\n- Use Python Jupyter Notebook to the plot the extraterrestrial solar radiation (Gon) between March 1st and July 1st. Copy and paste lines from the prior \u201cExtraterrestrial Radiation Calculator.py.py file\u201d. You can group common lines of code in their own cells that can be evaluated independently from the remaining cells. At the end of that file are a few lines of code that you can use to determine the day in the year given a date. Adding a # sign at the beginning of each line is used to comment and is not evaluated. Copy and paste these lines into a new Jupyter cell to determine what days the required dates correspond to. For example, to determine the day in the year on March 1st see below. Thus the day on March 1st is 61 and July 1st is 183.<\/li>\n<\/ol>\n\n\n\n
![\"\"](\"https:\/\/faculty.eng.ufl.edu\/jonathan-scheffe\/wp-content\/uploads\/sites\/765\/2020\/07\/GCS-small-300x40.png\")
<\/figure><\/p>\n\n\n\nNow we will plot between only these days by changing days to equal days=np.arange(61,183,1). Below is a screenshot of the cells and output.<\/p>\n\n\n\n
![\"\"](\"https:\/\/faculty.eng.ufl.edu\/jonathan-scheffe\/wp-content\/uploads\/sites\/765\/2020\/07\/GSC1367-mitchell-300x186.png\")
<\/figure><\/p>\n\n\n\n Back to top <\/a><\/p>\n\n\n\n <\/p>\n\n\n\n
Exercise 2 Instructions<\/a><\/h5>\n\n\n\n Solar Time<\/strong><\/p>\n\n\n\nSolar time is used in all sun-angle relationships. It is based on the apparent angular motion of the sun across the sky, with solar noon the time the sun crosses the meridian of the observer. Two corrections are needed to convert from standard time.<\/p>\n\n\n\n
\n- Correction for difference in longitude between observer\u2019s meridian and the meridian at which local standard time is based.<\/em><\/li>\n\n\n\n
- Correction from the equation of time which accounts for perturbations in the earth\u2019s rate of rotation.<\/em><\/li>\n<\/ol>\n\n\n\n
The equation used to calculate solar time is below.<\/p>\n\n\n\n
\nsolar time – standard time = 4(Lst<\/sub> – Lloc<\/sub>) + E<\/p>\n<\/div>\n\n\n\nLst<\/sub> is the standard meridian for the local time zone, Lloc is the longitude of the location in question and E is the equation of time in minutes. E can be determined graphically from the figure below (from Duffie and Beckmann, Solar Engineering of Thermal Processes, 4th Edition) or from the equations below from Duffie and Beckmann (Duffie and Beckmann, Solar Engineering of Thermal Processes, 4th Edition).<\/p>\n\n\n\n\nE = 229.2(0.000075 + 0.001868 \u2212 0.032077 \u2212 0.014615 2 \u2212 0.04089 2)<\/p>\n<\/div>\n\n\n
\n![\"\"](\"https:\/\/faculty.eng.ufl.edu\/jonathan-scheffe\/wp-content\/uploads\/sites\/765\/2020\/08\/Exercise-2-Equation-300x69.png\")
<\/figure>\n<\/div>\n\n\n<\/p>\n\n\n\n
Here, n is the day in the year and B has units degrees. The units for the right hand side of the equation used to calculate solar time are minutes.<\/p>\n\n\n
\n![\"\"](\"https:\/\/faculty.eng.ufl.edu\/jonathan-scheffe\/wp-content\/uploads\/sites\/765\/2020\/08\/Equation-of-time-vs-month_Ex2-300x243.png\")
<\/figure>\n<\/div>\n\n\nTo determine Lst, multiply the difference in time between local standard clock time and Greenwich Mean Time (GMT) by 15\u00b0. This relationship comes from thefact that the sun takes 4 minutes to traverse 1\u00b0 of longitude. Thus, if your local standard clock is 1 hour behind GMT then LST is 15\u00b0. You must account for Daylight Savings Time (DST): for example, DST in the United States takes place between the dates of March 8th and November 1st. During this time, the clocks are advanced by 1 hour and thus the difference between the clocks in the United States and GMT changes by 1 hour. For example, if DST is 14:00 then Local Standard Time (LST) which should be used to determine solar time and Lst, is 13:00.<\/p>\n\n\n\n
<\/p>\n\n\n\n
Figure 1. Equation of time verses month in the year. From Duffie and Beckmann, Solar Engineering of Thermal Processes, 4th Edition<\/p>\n\n\n\n
Exercise<\/strong>
Find the instructions and full Python code here<\/a>. <\/i><\/p>\n\n\n\n
- \n
- Click hyperlink to access CSV Data: Fig2212_data<\/a> <\/li>\n\n\n\n
- Click hyperlink to access CSV Data: Change_in_EnthalpyVSDelta<\/a><\/li>\n\n\n\n
- Click hyperlink to access CSV Data: H2O_Equilibrium_Constant<\/a> <\/li>\n<\/ul>\n\n\n\n
Book: <\/strong><\/p>\n\n\n\n
- \n
- Example Question 3 Solution: Solar Syngas Production_Problem 3<\/a> <\/li>\n\n\n\n
- Example Questions 4-7 Solutions: Book Chapter Problems_Gordon-v2<\/a><\/li>\n<\/ul>\n\n\n\n
Exercises Inspired by a Solar Energy Utilization Class<\/h5>\n\n\n\n
For your reference, below is a list of of the exercises in this series:<\/p>\n\n\n\n
- \n
- Exercise 1: Getting Familiar with Python: an Extraterrestrial Solar Radiation exercise. <\/a><\/li>\n\n\n\n
- Exercise 2: Solar Time and Solar Time Python Calculator. <\/a><\/li>\n<\/ul>\n\n\n\n
<\/p>\n\n\n\n
Exercise 1 Instructions<\/a><\/h5>\n\n\n\n
Extraterrestrial Solar Radiation<\/strong><\/p>\n\n\n\n
Solar radiation flux, or irradiance (G) in W\/m2<\/sup>, outside the earth\u2019s atmosphere and normal to earth\u2019s surface (Gon<\/sub>, where o stands for outside and n for normal direction) may be calculated as a function of the day in the year according to the following equation below.<\/p>\n\n\n
\n<\/figure>\n<\/div>\n\n\n<\/p>\n\n\n\n
Gsc<\/sub> is the solar constant and is equal to 1367 W\/m2 and n is the day in the year. The term contained within cosine is in units of degrees.<\/p>\n\n\n\n
Brief Python Introduction<\/strong><\/p>\n\n\n\n
Packages and Calling Functions<\/em>: In Python, there are several powerful computing packages commonly used that should be declared at the top of your .py, or python file. For example, some of the more popular computational packages useful for engineering include numpy, scipy, matplotlib and pandas, among others. These must be installed on your computer, either manually using a package installer (e.g. pip installer<\/a>), or are installed by default with many Python distributions, including the Anaconda Distribution that we will use in this exercise. For example, if we want to utilize the numpy package, we can use the import function and say:<\/p>\n\n\n\n
\nimport numpy as np<\/p>\n<\/div>\n\n\n\n
which allows us to call the numpy package simply by typing np. If we wish to use numpy package to calculate the cosine of an angle, we would first need to call it by typing np followed by a period and then the function contained within the package we wish to utilize, in this case cos, and the value within it (in this case the numerical value pi which is also called using the numpy package as np.pi) contained within parenthesis. For example,<\/p>\n\n\n\n
\nnp.cos(np.pi)<\/p>\n<\/div>\n\n\n\n
The input contained within parentheses should have units radians and the output will be in units radians. For example, to correctly use the equation for Gon, the cosine term will look like the following.<\/p>\n\n\n\n
\nnp.cos(360*n\/365*np.pi\/180)<\/p>\n<\/div>\n\n\n\n
because we need to convert from degrees to radians by multiplying everything by np.pi\/180.<\/p>\n\n\n\n
Declaring Functions:<\/em> Functions are declared in Python by defining them using def followed by the function name, parentheses that include any variables the function is dependent on, followed by a colon. The body of the function is then started on a new, indented line and the value that is returned by the function must be declared at the end using return<\/em>. For example, if we wish to define a function y that depends on x and is a linear line with slope m and intercept y and later call it, we would write<\/p>\n\n\n\n
\ndef y(x):<\/p>\n
return m*x+b<\/p>\n
y(x)<\/p>\n<\/div>\n\n\n\n
<\/p>\n\n\n\n
Exercise<\/strong>
Complete the following questions that deal with determining Gon<\/sub> using the above equation and a more accurate form of the equation (not shown, Gon_accurate function), as a way to get comfortable computing in Python using some basic techniques, learning to call and declare functions, plotting data using the package matplotlib, and using two different methods to evaluate python code, namely a traditional editor named Spyder and a web based notebook style editor called Jupyter.<\/p>\n\n\n\n- \n
- Download Python<\/a> (\u201cAnaconda Distribution\u201d). This includes many useful packages but most commonly we will use Spyder and Jupyter, two different ways to evaluate Python code.<\/li>\n\n\n\n
- Open and run \u201cExtraterrestrial Radiation Calculator.py\u201d using Spyder. It can be copied and pasted from the appendix of <\/span>this document<\/a><\/strong> and saved to your computer as a .py file. After running, you should see a plot with two curves.
<\/p>
<\/li>\n\n\n\n- Replot the data using the same code but change the plotted lines from red \u2018r\u2019 and blue \u2018b\u2019in the plt.plot lines of code, to other colors of your choice. Here are examples and more documentation<\/a>.
<\/figure><\/li>\n\n\n\n- We have plotted data every day using the function np.arrange with the linedays=np.arange(1,365,1). This creates an array from 1 to 365 with a spacing of 1. However, let’s imagine we want to see individual data points plotted \u2013 it would be difficult to distinguish between them with so many. Therefore, lets plot a point every 10 days by changing the spacing from 1 to 10. Show the results using blue circles for the blue curve by inserting \u2018bo\u2019 in place of \u2018b\u2019.
<\/figure><\/li>\n\n\n\n- Use Python Jupyter Notebook to the plot the extraterrestrial solar radiation (Gon) between March 1st and July 1st. Copy and paste lines from the prior \u201cExtraterrestrial Radiation Calculator.py.py file\u201d. You can group common lines of code in their own cells that can be evaluated independently from the remaining cells. At the end of that file are a few lines of code that you can use to determine the day in the year given a date. Adding a # sign at the beginning of each line is used to comment and is not evaluated. Copy and paste these lines into a new Jupyter cell to determine what days the required dates correspond to. For example, to determine the day in the year on March 1st see below. Thus the day on March 1st is 61 and July 1st is 183.<\/li>\n<\/ol>\n\n\n\n
<\/figure><\/p>\n\n\n\nNow we will plot between only these days by changing days to equal days=np.arange(61,183,1). Below is a screenshot of the cells and output.<\/p>\n\n\n\n
<\/figure><\/p>\n\n\n\nBack to top <\/a><\/p>\n\n\n\n
<\/p>\n\n\n\n
Exercise 2 Instructions<\/a><\/h5>\n\n\n\n
Solar Time<\/strong><\/p>\n\n\n\n
Solar time is used in all sun-angle relationships. It is based on the apparent angular motion of the sun across the sky, with solar noon the time the sun crosses the meridian of the observer. Two corrections are needed to convert from standard time.<\/p>\n\n\n\n
- \n
- Correction for difference in longitude between observer\u2019s meridian and the meridian at which local standard time is based.<\/em><\/li>\n\n\n\n
- Correction from the equation of time which accounts for perturbations in the earth\u2019s rate of rotation.<\/em><\/li>\n<\/ol>\n\n\n\n
The equation used to calculate solar time is below.<\/p>\n\n\n\n
\nsolar time – standard time = 4(Lst<\/sub> – Lloc<\/sub>) + E<\/p>\n<\/div>\n\n\n\n
Lst<\/sub> is the standard meridian for the local time zone, Lloc is the longitude of the location in question and E is the equation of time in minutes. E can be determined graphically from the figure below (from Duffie and Beckmann, Solar Engineering of Thermal Processes, 4th Edition) or from the equations below from Duffie and Beckmann (Duffie and Beckmann, Solar Engineering of Thermal Processes, 4th Edition).<\/p>\n\n\n\n
\nE = 229.2(0.000075 + 0.001868 \u2212 0.032077 \u2212 0.014615 2 \u2212 0.04089 2)<\/p>\n<\/div>\n\n\n
\n<\/figure>\n<\/div>\n\n\n<\/p>\n\n\n\n
Here, n is the day in the year and B has units degrees. The units for the right hand side of the equation used to calculate solar time are minutes.<\/p>\n\n\n
\n<\/figure>\n<\/div>\n\n\nTo determine Lst, multiply the difference in time between local standard clock time and Greenwich Mean Time (GMT) by 15\u00b0. This relationship comes from thefact that the sun takes 4 minutes to traverse 1\u00b0 of longitude. Thus, if your local standard clock is 1 hour behind GMT then LST is 15\u00b0. You must account for Daylight Savings Time (DST): for example, DST in the United States takes place between the dates of March 8th and November 1st. During this time, the clocks are advanced by 1 hour and thus the difference between the clocks in the United States and GMT changes by 1 hour. For example, if DST is 14:00 then Local Standard Time (LST) which should be used to determine solar time and Lst, is 13:00.<\/p>\n\n\n\n
<\/p>\n\n\n\n
Figure 1. Equation of time verses month in the year. From Duffie and Beckmann, Solar Engineering of Thermal Processes, 4th Edition<\/p>\n\n\n\n
Exercise<\/strong>
Find the instructions and full Python code here<\/a>. <\/i><\/p>\n\n\n\n - Correction from the equation of time which accounts for perturbations in the earth\u2019s rate of rotation.<\/em><\/li>\n<\/ol>\n\n\n\n
- Open and run \u201cExtraterrestrial Radiation Calculator.py\u201d using Spyder. It can be copied and pasted from the appendix of <\/span>this document<\/a><\/strong> and saved to your computer as a .py file. After running, you should see a plot with two curves.
- Exercise 2: Solar Time and Solar Time Python Calculator. <\/a><\/li>\n<\/ul>\n\n\n\n
- Example Questions 4-7 Solutions: Book Chapter Problems_Gordon-v2<\/a><\/li>\n<\/ul>\n\n\n\n
- Click hyperlink to access CSV Data: Change_in_EnthalpyVSDelta<\/a><\/li>\n\n\n\n