Trace:

rainbow_project

A PCRE internal error occured. This might be caused by a faulty plugin

====== Differences ====== This shows you the differences between two versions of the page.

rainbow_project [2014/05/29 06:09] wikimanager [Questions about theory and coding] |
rainbow_project [2014/06/03 19:18] (current) delakins [Coding tasks] |
||
---|---|---|---|

Line 117: | Line 117: | ||

* $P_\text{sc}$ $\,=\,I_\text{sc}\!\cdot\!{\text{Area}}_\text{sc}$ $\,\approx\,I_\text{sc}\!\cdot\!2\pi\,(D\sin\gamma)\,D\Delta \gamma\;$, where $\,D$ is the distance from the droplet to the observer. | * $P_\text{sc}$ $\,=\,I_\text{sc}\!\cdot\!{\text{Area}}_\text{sc}$ $\,\approx\,I_\text{sc}\!\cdot\!2\pi\,(D\sin\gamma)\,D\Delta \gamma\;$, where $\,D$ is the distance from the droplet to the observer. | ||

* $P_\text{sc}$ $=s\!\cdot\!P_\text{inp}$ | * $P_\text{sc}$ $=s\!\cdot\!P_\text{inp}$ | ||

- | * $I_\text{sc}$ $=s\!\cdot\!I_\text{inp}\frac{d}{D^2\sin\gamma}\frac{\Delta d}{\Delta\gamma}$ $=s \cdot I_\text{inp}\frac{R^2\delta}{D^2\sin\gamma}\frac{\Delta \delta}{\Delta\gamma}$ $\sim I_\text{inp}\frac{\delta}{\sin\gamma}\frac{\Delta \delta}{\Delta\gamma}\;$: here the factors //s//, //R//, and //D// are pretty constant (or too complicated to calculate) and do not contribute to the relative intensity at various angles and wavelengths | + | * $I_\text{sc}$ $=s\!\cdot\!I_\text{inp}\frac{d}{D^2\sin\gamma}\left|\frac{\Delta d}{\Delta\gamma}\right|$ $=s \cdot I_\text{inp}\frac{R^2\delta}{D^2\sin\gamma}\left|\frac{\Delta \delta}{\Delta\gamma}\right|$ $\sim I_\text{inp}\frac{\delta}{\sin\gamma}\left|\frac{\Delta \delta}{\Delta\gamma}\right|\;$: here the factors //s//, //R//, and //D// are pretty constant (or too complicated to calculate) and do not contribute to the relative intensity at various angles and wavelengths |

| {{ :projects:rainbow:rainbowdrawing2.png?nolink |}} | | | {{ :projects:rainbow:rainbowdrawing2.png?nolink |}} | | ||

Line 124: | Line 124: | ||

====Questions about theory and coding==== | ====Questions about theory and coding==== | ||

- I was wondering if the relative scattered intensity was suppose to be calculated with the Rayleigh forumula? | - I was wondering if the relative scattered intensity was suppose to be calculated with the Rayleigh forumula? | ||

- | * **Answer**: No, Rayleigh formula refers to scattering of light off of objects whose size is much smaller than the wavelength of light (see Wikipedia article on [[wp>Rayleigh scattering]]). The rain droplets are 3 mm on average (as somebody posted in the data section below), much larger than the wavelengths of light. instead, you should use the ray diagrams and Snell's law and the law of reflection, see the Figures above. You need to figure out the angles in Fig. 4 using these laws, then put the formulas into excel. | + | * **Answer**: No, Rayleigh formula refers to scattering of light off of objects whose size is much smaller than the wavelength of light (see Wikipedia article on [[wp>Rayleigh scattering]]). The rain droplets are 3 mm on average (as somebody posted in the data section below), much larger than the wavelengths of light. Instead, you should use the ray diagrams and Snell's law and the law of reflection, see the Figures above. You need to figure out the angles in Fig. 4 using these laws, then put the formulas into excel. |

- I figured out the trigonometry/geometry of the raindrop and light rays. I am having trouble with Excel, however. Please view the attached spreadsheet (highlighted cells). The value that Excel is calculating for Sin(-4....) (''CELL AJ06'') is different from what I get on my calculator, Wolframalpha, and Google. This is preventing me from formatting any further as I continue to get an ERROR message (''CELL AP06''). | - I figured out the trigonometry/geometry of the raindrop and light rays. I am having trouble with Excel, however. Please view the attached spreadsheet (highlighted cells). The value that Excel is calculating for Sin(-4....) (''CELL AJ06'') is different from what I get on my calculator, Wolframalpha, and Google. This is preventing me from formatting any further as I continue to get an ERROR message (''CELL AP06''). | ||

* **Answer**: I think I know what your problem is. Excel always used radians for all trig functions. That's why you have to use ''DEGREES()'' to convert angles to degrees. But then, in columns ''X:AC'', you use ''SIN(R5)'' and so on, but cells like ''R5'' contain the angle in degrees. The ''SIN()'' always expects the radians though, and calculates the wrong value. Then it is all corrupt after that. | * **Answer**: I think I know what your problem is. Excel always used radians for all trig functions. That's why you have to use ''DEGREES()'' to convert angles to degrees. But then, in columns ''X:AC'', you use ''SIN(R5)'' and so on, but cells like ''R5'' contain the angle in degrees. The ''SIN()'' always expects the radians though, and calculates the wrong value. Then it is all corrupt after that. | ||

* Remember, radians and degrees are related through the factor of $\frac{\pi}{180}$. You can either keep everything in degrees up to the last final plot, or use a conversion factor every time you take a ''SIN()'' or ''TAN()''. | * Remember, radians and degrees are related through the factor of $\frac{\pi}{180}$. You can either keep everything in degrees up to the last final plot, or use a conversion factor every time you take a ''SIN()'' or ''TAN()''. | ||

- I think my problem is that I am apparently not understanding what needs to be done... | - I think my problem is that I am apparently not understanding what needs to be done... | ||

- | * **Answer**: You've done the first step -- got your hand on the actual data for $n(\lambda)$ dependence for pure water. Now you need to do the rest of the analysis. It is all explained on this wiki. Briefly, you need to | + | * **Answer**: You've done the first step -- got your hands on the actual data for $n(\lambda)$ dependence for pure water. Now you need to do the rest of the analysis. It is all explained on this wiki. Briefly, you need to |

- Look at Fig. 4 above and, using the parameter $d$ as a given, as well as the index of refraction $n$, calculate $\beta_1+\beta_2+\beta_3'$ (this is the total deflection angle for the 1<sup>st</sup> rainbow), and, optionally, $\beta_1+\beta_2+\beta_3+\beta_4'$ (this is the total deflection angle for the 2<sup>nd</sup> rainbow). Basically, you need to obtain these two formulas (a bunch of sin() functions, really) that depend on $n$ and on $d/R$ ratio. | - Look at Fig. 4 above and, using the parameter $d$ as a given, as well as the index of refraction $n$, calculate $\beta_1+\beta_2+\beta_3'$ (this is the total deflection angle for the 1<sup>st</sup> rainbow), and, optionally, $\beta_1+\beta_2+\beta_3+\beta_4'$ (this is the total deflection angle for the 2<sup>nd</sup> rainbow). Basically, you need to obtain these two formulas (a bunch of sin() functions, really) that depend on $n$ and on $d/R$ ratio. | ||

- Set up a grid of possible $d/R$ ratios from 0 (normal incidence) to 1 (grazing beam), with some small step, e.g. 0.001 | - Set up a grid of possible $d/R$ ratios from 0 (normal incidence) to 1 (grazing beam), with some small step, e.g. 0.001 | ||

- | - On that grid, for a particular lambda, calculate the total deflection angle as a function of d/R (i.e. your grid) | + | - On that grid, for a particular $\lambda$, calculate the total deflection angle as a function of $d/R$ (i.e. your grid) |

- Then, on the same grid, calculate the intensity of the deflected light, assuming some (any, really) intensity of the sunlight | - Then, on the same grid, calculate the intensity of the deflected light, assuming some (any, really) intensity of the sunlight | ||

- (//optional, for maximum grade//) - convert this intensity as a function of the $d/R$ grid to the total deflection angle grid. That would be the model of what you observe when you look at the actual rainbow, for each color. | - (//optional, for maximum grade//) - convert this intensity as a function of the $d/R$ grid to the total deflection angle grid. That would be the model of what you observe when you look at the actual rainbow, for each color. | ||

- | | + | - I am having difficulty substituting $\beta_1$ into Snell's law. Is it possible I could have a hint based on the geometry of the diagram for the rainbow project to lead me in the right direction? |

+ | * **Answer**: If you look at a point where the beam enters the droplet, and notice that $\theta_i$ is equal to $\theta_i$ immediately on the other side of the water surface, then you can just see that $\beta_1$ is the difference $\theta_i-\theta_r$. See that on the drawing (Fig. 4)? The next reflection point is symmetric with this one, but the meaning of angles is different. Also, if you spot a right-angled triangle with $R$ (radius of the sphere) and $d$ as its sides, and the angle $\theta_i$ as one of its angles, that will allow you to relate $\sin\theta_i$ to the $d/R$ ratio. The flow of math is: | ||

+ | - from $\;d/R\;$ find $\;\sin\theta_i$ | ||

+ | - from that find $\theta_i$ | ||

+ | - from that, using Snell's law, find $\theta_r$ | ||

+ | - from that, find $\beta_1$ | ||

+ | - How would I start solving for beta $\beta_2$, $\beta_3$, and $\beta_4$? | ||

+ | * **Answer**: If you look at the next reflection (Fig. 4), you can see that $\beta_2+\theta_r+\theta_r$ add up to 180$^\circ$ (or, $\pi$ radians). There, from symmetry, $\theta_r$ is the same as $\theta_r$ in the first refraction. All the other points are also symmetric to the original, so same logic applies. $\beta_3'$ is the difference between two thetas already found earlier. And so on. | ||

+ | | ||

=====Data===== | =====Data===== | ||

Give yourself 20 minutes to google or search the literature for the index of refraction of water as a function of wavelength. Paste references and data (formulas and raw numbers if compact) below. | Give yourself 20 minutes to google or search the literature for the index of refraction of water as a function of wavelength. Paste references and data (formulas and raw numbers if compact) below. | ||

- **Index of refraction for water** = 1.33. This is the average value. We need its dependence on the wavelength $\lambda$! | - **Index of refraction for water** = 1.33. This is the average value. We need its dependence on the wavelength $\lambda$! | ||

- **Average radius of water droplet** = 3 mm ([[rainbow project#references|Ref. [4]]]) | - **Average radius of water droplet** = 3 mm ([[rainbow project#references|Ref. [4]]]) | ||

+ | - **Average [[wp>Sunlight]] Intensity on Earth** = 1367 $\frac{W}{m^2}$ | ||

{{visible_spectrum.jpg}} | {{visible_spectrum.jpg}} | ||

| **Table 1**. Refraction index of pure water. ||||| | | **Table 1**. Refraction index of pure water. ||||| | ||

- | | | at $\rho=999.7\,\frac{\text{kg}}{ {\text m}^3}$ \\ //Calculated by N.K. from Ref. 2//<color white>... ... ... ... ... ... ... ...</color> || at $T=20\,^\circ$C,\\ //adapted from Refs. 1, 5 by B.V.// | at $T=25\,^\circ$C,\\ //adapted from Refs. 1, 6 by B.V.// | | + | | | at $T=10\,^\circ$C, $\rho=999.7\,\frac{\text{kg}}{ {\text m}^3}$ \\ //Calculated by N.K. from Ref. 2//<color white>... ... ... ... ... ... ... ...</color> || at $T=20\,^\circ$C,\\ //adapted from Refs. 1, 5 by B.V.// | at $T=25\,^\circ$C,\\ //adapted from Refs. 1, 6 by B.V.// | |

^ $\lambda$ (nm) ^ $(n^2−1)/(n^2\!+2)$ $\text{at}$ $10\,^\circ{\text C}$ ^ $n\big(10\,^\circ{\text C}\big)$ ^ $n\big(20\,^\circ{\text C}\big)$ ^ $n\big(25\,^\circ{\text C}\big)$ ^ | ^ $\lambda$ (nm) ^ $(n^2−1)/(n^2\!+2)$ $\text{at}$ $10\,^\circ{\text C}$ ^ $n\big(10\,^\circ{\text C}\big)$ ^ $n\big(20\,^\circ{\text C}\big)$ ^ $n\big(25\,^\circ{\text C}\big)$ ^ | ||

| 400 | 0.211798 | 1.343924521 | 1.34356 | 1.339 | | | 400 | 0.211798 | 1.343924521 | 1.34356 | 1.339 | |

rainbow_project.1401343752.txt.gz · Last modified: 2014/05/29 06:09 by wikimanager

Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Noncommercial-Share Alike 3.0 Unported