Citat:
Ursprungligen postat av BengtZz
Men varför får man totalreflektion om det är varmare vid marken? Den kontinuerliga förändringen i brytningsindex bör ju i sådana fall göra att ljuset bryts mer och mer. Tittar vi på en infinitesmal bit av stålen i en infinitesmal del av atmossfären så kommer vi inte i någon punkt kunna visa på totalreflektion. Så varför är det en totalreflektion? Rita gärna upp en bild där ljusets bryts kontinuerligt och titta på en liten bit åt gången.
Nja, fast även när det inte är totalreflektion så reflekteras ju en del av ljuset. När ljusets infallsvinkel blir flackare och flackare så kommer mer och mer av ljuset reflekteras, så på så sätt så blir det totalreflektion.
Frågan om huruvida det blir totalreflektion handlar om huruvida ljusvinkeln "hinner" bli horisontell, innan den nått marken. Om den gör det, så blir det ju totalreflektion, i bemärkelse att inget ljus når marken. Om däremot ljusets kurva såg ut som t.ex. y = 1 - sqrt(x), så hinner den inte bli horisontell innan den når marken (oavsett hur långt upp man börjar faktiskt), så totalreflektion uppstår inte.
Så hur avgöra vad som händer? Tja, matte såklart. Antag att luften består av horisontella skikt, och att skiktet på höjd z har brytningsvinkeln n(z). Antag att vi har en infallande ljusstråle, som går mot marken, och så att ljusstrålen har vinkeln v(z) mot lodlinjen när den korsar skiktet på höjd z. Målet är att ställa upp en differentialekvation som relaterar v(z) till n(z).
Men, om vill tillåts missbruka Leibniz notation lite, så går det att resonera såhär: Låt säga att vi har två skikt, väldigt nära varandra, ett på höjd z och ett på höjd dz. Antag att dz är så litet att vi kan se det som att vi bara har ett brytningsskikt. Brytningen här ges då av
sin(v(z+dz)) / sin(v(z)) = n(z)/n(z + dz).
Logaritmera båda sidor, och vi får
log sin(v(z+dz)) - log sin(v(z)) = - (log n(z+dz) - log n(z)).
Dividerar med dz, och låt dz gå mot 0:
d/dz (log sin(v(z))) = d/dz (log n(z)).
Detta är en diffekvation av den übersnälla varianten, som bara kan lösas genom att integrera båda sidor:
log sin(v(z)) = log n(z) + C
eller med andra ord
sin(v(z)) = C / n(z).
Men låt säga att ursprungsvillkoret är att strålen hade vinkeln v0 på en höjd där brytningsindexed var n0. Sätter vi in detta, så kan vi lösa ut C som C = n0 sin(v0), och stoppar vi in detta, så får vi
sin(v(z)) = sin(v0) n0 / n(z).
Detta visar att ljusstrålens vinkel vid ett viss skikt faktiskt bara beror på brytningsindexet vid det skiktet, och inte på brytningsindex vid mellanliggande skikt. (Hur långt ljustrålen har färdats horisontellt mellan skikten beror dock på vad som händer mellan.) Hursomhelst, om nu n(z) är så litet så att sin(v0) n0 / n(z) > 1 så kan ekvationen ovan omöjligen gälla. Ljusstrålen kommer aldrig nå ett skikt så långt ner att n(z) vid det skiktet uppfyller sin(v0) n0 / n(z) > 1, och alltså måste totalreflektion ha uppstått på vägen.