Navigaatiofysiikkaa III – Etäisyyden määrityksestä

Navigaatiokirjallisuudessa, esim. [1], löytyvä laskukaava etäisyyden laskemiseksi majakkaan, silloin, kun majakan valo tulee juuri näkyviin horisontin takaa, on

 $s=2{,}08\ \left(\sqrt{\frac{H}{\mathrm{m}}}\ +\ \sqrt{\frac{h}{\mathrm{m}}}\right)\ \mathrm{M}$(1)

jossa $H$ ja $h$ ovat majakan (valon) ja silmän korkeudet merenpinnasta. M = meripeninkulma (Nautical mile) = 1852 m. Tämän kirjoitussarjan ensimmäisessä osassa [2] saimme geometrisellä tarkastelulla tuloksen

 $s=1{,}927\ \left(\sqrt{\frac{H}{\mathrm{m}}}\ +\ \sqrt{\frac{h}{\mathrm{m}}}\right)\ \mathrm{M}$(2)

Yllä oleva geometrinen tulos poikkeaa yhtälöstä (1) noin 8

 $\cos\theta=\ \frac{R}{r}\ \cdot\ \frac{1+\ a\rho_0}{1+\ a\rho_0e^{\left(-\ \frac{\rho_0gT_0}{p_0\alpha}+1\right)\ln\frac{1}{1-\ \alpha\ (r-R)/T_0}}}$(3)

Tämä yhtälö siis vastasi tilannetta, jossa valonsäde sivuaa yhdessä kohtaa meren pintaa, eli $\theta = 0$, kun $r=R$.

Oletetaan tässä normaali ilmanpaine $p_0=101{,}325\ \mathrm{kPa}$  ja lämpötilaa $T_0=+25\ °\mathrm{C}$ vastaava (ideaalikaasun) tiheys $\rho_0=1{,}185\ \mathrm{kg/m^3}$, saadaan $b=\frac{\rho_0g}{p_0}=\ 0{,}000115\ 1/\mathrm{m}$. Vakiolle $\alpha$ käytämme arvoa $\alpha=0{,}0065\ \mathrm{K/m\ }$ [4, 5]. Edelleen ilman taitekertoimen tiheysriippuvuuden vakiolle a käytämme lukuarvoa $a=0{,}000228\ \mathrm{m}^3/\mathrm{kg}$ [4].

Hyödynnämme tätä yhtälöä (3) etäisyyden laskemiseksi majakkaan. Etäisyyden määrittämiseksi tarvitaan korkeus ($h = r-R$) valon maan keskipisteen suhteen kiertyneen kulman  funktiona.  Nyt valon kulkusuunnan (vaakatason kanssa muodostaman) kulman $\theta$ tangentti on kulkukäyrän kulmakerroin

 $\tan\theta=\ \frac{dy}{dx}$(4)
Kuva 1.

Siirryttäessä karteesisesta xy-koordinaatistosta napakoordinaatteihin ($r$, $\varphi$)  pätee (lähtöpisteen ollessa y-akselilla) infinitesimaalisella tarkkuudella $dy\ =\ dr$ ja $dx=rd\varphi$ (kts. kuva 1.). Joten

 $\frac{dr}{d\varphi}=r\tan\theta=r\ \frac{\sqrt{1-\mathrm{\cos}^2\theta\ }}{\cos\theta}=\ \frac{r\ \sqrt{r^2-\ R^2\ \left(f(r-R)\right)^2}}{R\ f(r-R)}$(5)

jossa tangentti on lausuttu kosinin avulla ja funktio $f(x)$ on yhtälön (3) perusteella

 $f\left(x\right)=\ \frac{1+\ a\rho_0}{1+\ a\rho_0e^{\left(-\ \frac{b\ T_0}{\alpha}+1\right)\ln\frac{1}{1-\ \alpha\ x/T_0}}}$(6)

Yhtälöstä (5) saadaan integroitua etäisyys isoympyrää pitkin

 $s=R\varphi=\int_R^r\frac{R^2\ f\left(r-R\right)dr}{r\ \sqrt{r^2-\ R^2\ \left(f(r-R)\right)^2}}=\ \int_0^H\frac{R^2\ f\left(x\right)dx}{(R+x)\ \sqrt{(R+x)^2-\ R^2\ \left(f(x)\right)^2}}$(7)

Ottaen tällöin huomioon, että $H\ll\ R$ ja $bH\ll1$, kehitetään integroitava muuttujan $x$ johtavaan kertalukuun (vrt. McLaurinin sarja). Siten saadaan

 $s\cong\ R\ \int_0^H\frac{dx}{\sqrt{2Rx-2R^2a\rho_0\left(b-\ \alpha/T_0\right)x/(1+a\rho_0)}}=\ \frac{\sqrt{2RH}}{\sqrt{1-Ra\rho_0\left(b-\ \alpha/T_0\right)/\left(1+a\rho_0\right)}}$(8)

Kirjoitetaan yhtälö (8) meripeninkulman M avulla kuten tämän kirjoitussarjan ensimmäisessä osassa [2]

 $s=\ \sqrt{\frac{360\cdot60}{1852\ \pi}\ }\cdot\ \frac{1}{\sqrt{1-Ra\rho_0\left(b-\ \alpha/T_0\right)/\left(1+a\rho_0\right)}}\ \cdot\ \sqrt{\frac{H}{\mathrm{m}}}\ \mathrm{M}$(9)

Yhtälössä (9) ensimmäinen kerrointermi on arvoltaan 1,927, toinen termi puolestaan on 1,091, jolloin yhteensä

 $s=2{,}10\ \mathrm{M}\sqrt{\frac{H}{\mathrm{m}}}$(10)

Valon taittumisen vaikutus etäisyysmittaukseen on huomattavan suuri, noin 9

 $90°-θ≅90°$(11)

Matemaattisesti tämän taitekertoimen pienen muutoksen aiheuttaman vaikutuksen valon kulku­kulmaan näkee kirjoitussarjan aiemmassa osassa [3] johdetusta valon kulkusuunnan ja taitekertoimen välisestä relaatiosta:

 $d\theta=\ \frac{dn}{n\tan\theta}\ \cong\ \frac{dn}{n \theta}$(12)

jolloin hyvin pienikin taitekertoimen $dn$ muutos aiheuttaa suhteellisen ison kulmamuutoksen, koska kulma $\theta $ itsessään on pieni.

Kerroin yhtälössä (10) on prosentin suurempi kuin, mitä kirjallisuudessa [1] on annettu. Tässä on kuitenkin huomattava, että kertoimen tarkka arvo riippuu lämpötilasta ja paineesta ja lisäksi maan gravitaatiokentästä purjehduspaikalla. Edellä on oletettu lämpötilan olevan $+25\ °\mathrm{C}$, ilmanpaineen normaali ja painovoimakiihtyvyys ns. normipaikkakunnalla, $g=9{,}80665\frac{\mathrm{m}}{\mathrm{s}^{\mathrm{2}}}$.

Etäisyyden voisikin siten kirjoittaa lopuksi (huomioiden myös silmän korkeuden $h$):

 $s=\ \frac{1{,}927}{\sqrt{1-a\ R\cdot\frac{pM_{ilma}\ }{R_mT^2}\left(\frac{gM_{ilma}}{R_m}-\ \alpha\right)/\left(1+a\ \frac{pM_{ilma}\ }{R_mT}\right)\ }}\ \cdot\ \left(\sqrt{\frac{H}{\mathrm{m}}}\ +\ \sqrt{\frac{h}{\mathrm{m}}}\right)\ \mathrm{M}$(13)

Tämän etäisyyskaavan kerroin on isompi korkeapaineen vallitessa. Isompiin kertoimiin päädytään myös, jos tyydytään purjehtimaan vain korkeilla pohjoisilla tai eteläisillä latitudeilla, jolloin lämpötila $T$ on yleensä alempi ja myös painovoimakiihtyvyys $g$ on keskimääräistä suurempi.

Majakan etäisyyden laskemiseksi tarkasti pitäisi siis lukea merikortista majakan valon korkeus ($H$), määrittää silmän korkeus meren pinnasta ($h$), selvittää painovoimakiihtyvyys purjehdus­paikalla ($g$) ja lukea ilmapuntari ($p$) sekä lämpömittari ($T$).

Lähdeviitteet

[1] Suomen Navigaatioliitto: Veneilijän merenkulkuoppi II, Rannikko-navigointi, Unigrafia Oy, Helsinki, 2013

[2] T. Toimela, Navigaatiofysiikkaa I, Dimensio, 5.2. 2020 https://www.dimensiolehti.fi/navigaatiofysiikkaa-i-valon-kulun-geometriaa/

[3] T. Toimela, Navigaatiofysiikkaa II, Dimensio, 12.2. 2020 https://www.dimensiolehti.fi/navigaatiofysiikkaa-ii-valon-taittuminen-ilmakehassa/

[4] David R. Lide (ed.), Handbook of Chemistry and Physics, 87th edition, Taylor Francis Group, 2006

[5] Inkinen-Tuohi, Momentti 1, Insinöörikoulutuksen fysiikka, Otava, 1999

Aloituskuva Phot Daniel Kuruvilla on Unsplash

Kirjoittaja