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)
Spotlight: Akron Brass FireCalc Pocket Computer
Spotlight: Akron Brass FireCalc Pocket Computer Welcome to a special Monday Edition of Eddie’s Math and Calculator blog. Thi...