A blog is that is all about mathematics and calculators, two of my passions in life.
Sunday, December 31, 2017
Friday, December 22, 2017
Optics: Transmitted Angle, Critical Angle
Optics: Transmitted Angle, Critical Angle
Formulas Used
Reflected Angle: θr =
θ
Transmitted Angel: θt
= asin((n1 * sin θ)/n2)
Critical Angle: θc =
asin(n1/n2), n2 > n1
Perpendicular Frensel Parameters (t, r):
t = (2 * n1 * cos θ) / (n1 * cos θ + n2 * cos θt)
r = (n1 * cos θ – n2 * cos θt) / (n1 * cos θ + n2 * cos θt)
Parallel Fresnel Parameters (t, r):
t = (2 * n1 * cos θ) / (n1 * cos θt + n2 * cos θ)
r = (n1 * cos θt – n2 * cos θ) / (n1 * cos θt + n2 * cos θ)
Table of Indices of Reflection (n)
Air
|
1.00
|
Moissanite (SiC)
|
2.65
|
Carbon Dioxide (CO2)
|
1.0045
|
Pyrex
|
1.47
|
Cubic Zirconia
|
2.15
|
Salt (NaCl)
|
1.54
|
Diamond
|
2.42
|
Sapphire
|
1.76
|
Glass
|
1.52
|
Silicon
|
3.45
|
Ice
|
1.31
|
Water
|
1.33
|
HP Prime Program FLATOPTIC
EXPORT
FLATOPTIC()
BEGIN
// 2017-12-17 EWS
// Solar Energy
// Change to degree
HAngle:=1;
// index of refraction
LOCAL L0:={1,1.0045,2.15,2.42,
1.52,1.31,2.65,1.47,
1.54,1.76,3.45,1.33};
LOCAL L1:={"Air","CO2",
"Cubic
Zirconia","Diamond",
"Glass","Ice",
"Moissanite
(SiC)","Pyrex",
"Salt (NaCl)","Sapphire",
"Silicon","Water"};
LOCAL L2:={"Perpendicular",
"Parallel",
"None"};
LOCAL n1,n2,k1,k2,θ,t,r;
LOCAL θc,θt,θr,x;
INPUT({{k1,L1},{k2,L1},θ,{x,L2}},
"Optic - Flat Interface",
{"Medium 1: ","Medium 2:
",
"Angle (°):
","Polarized?"}
);
// Calculation
PRINT();
θr:=θ;
PRINT("Reflect Angle:
"+θr+"°");
n1:=L0(k1); n2:=L0(k2);
IF n2>n1 THEN
θc:=ASIN(n1/n2);
PRINT("Critical Angle:
"+θc+"°");
END;
θt:=ASIN(n1*SIN(θ)/n2);
PRINT("Transmitted Angle: "+
θt+"°");
IF x==1 THEN
t:=(2*n1*COS(θ))/
(n1*COS(θ)+n2*COS(θt));
r:=(n1*COS(θ)-n2*COS(θt))/
(n1*COS(θ)+n2*COS(θt));
PRINT("t: "+t);
PRINT("r: "+r);
END;
IF x==2 THEN
t:=(2*n1*COS(θ))/
(n1*COS(θt)+n2*COS(θ));
r:=(n1*COS(θt)-n2*COS(θ))/
(n1*COS(θt)+n2*COS(θ));
PRINT("t: "+t);
PRINT("r: "+r);
END;
END;
TI-84 Plus CE Program
"EWS
2017-12-17"
{1,1.0045,2.15,2.42,1.52,1.31,2.65,1.47,1.54,1.76,3.45,1.33}→L6
Degree
For(K,1,2)
Disp "1. AIR 2.
CO2 3. CUB.ZIC.","4. DIAMOND 5. GLASS 6. ICE","7.
MOISANITE","8. PYREX 9. NACL","10. SAPPHIRE 11.
SIL.","12. WATER"
Input M
If K=1:M→N
End
L6(N)→N:L6(M)→M
Input
"INCIDENT: °",θ
sin^-1(N*sin(θ)/M)→A
Disp
"TRANSMITED: °",A
If M>N:Then
sin^-1(N/M)→C
Disp "CRITICAL:
°",C
End
Pause
Menu("POLARISE?","PERPENDICULAR",1,"PARALLEL",2,"NONE",3)
Lbl 1
(2Ncos(θ))/(Ncos(θ)+Mcos(A))→T
(Ncos(θ)-Mcos(T))/(Ncos(θ)+Mcos(T))→R
Disp
"T",T,"R",R
Goto 3
Lbl 2
(2Ncos(θ))/(Ncos(T)+Mcos(θ))→T
(Ncos(T)-Mcos(θ))/(Ncos(T)+Mcos(θ))→R
Disp
"T",T,"R",R
Lbl 3
Source:
Klaus Jäger,
Olindo Isabella, Arno H.M. Smets, René
A.C.M.M. van Swaaij, Miro Zeeman Solar
Energy: Fundamental, Technology, and Systems Delft University of Technology, 2014
(no money is made from this blog entry)
Hopefully you find this useful.
Eddie
This blog is property of Edward Shore, 2017
Subscribe to:
Posts (Atom)
RPN with HP 15C & DM32: Swapping Items in the Stack
RPN with HP 15C & DM32: Swapping Items in the Stack A New Series Welcome to a new series: RPN with the HP 15C & DM32. ...