Zoznam úloh

1. Janči mal tieň

Zadanie

Janči si pokojne vykračoval a zrazu ho oziabalo. O koľko klesne množstvo energie dopadajúce na zem kvôli zatmeniu Slnka v mieste, kde je zatmenie úplné?

Uvažujme dva modely:

  1. Model s konštantným vyžarovaním: Povrch Slnka vyžaruje rovnomerne do každého smeru.

  2. Lineárny model s korekciou do prvého rádu: Povrch Slnka vyžaruje ako

$$ \frac{I(\psi)}{I(0)} = (1-u) + u \cos \psi \, , $$

kde $I(\psi)$ je intenzita žiarenia z povrchu Slnka pod uhlom $\psi$ (meraný od normály) a $u$ je fixný parameter modelu.

Pre prvý model vypočítajte profil poklesu. Pre oba modely vypočítajte priemernú hodnotu poklesu počas celého zatmenia a popíšte závislosť priemerného poklesu na parametri $u$.

Predpokladajte, že:

  • Slnko je guľa.

  • Slnko a Mesiac majú kruhový obrys rovnakej veľkosti.

  • Dráhy Slnka a Mesiaca sú zarovnané a stred Mesiaca prechádza cez stred Slnka.

Vyriešme si najprv \textbf{jednoduchší model}, bez parametra $u$. Na popis polohy Mesiaca použíjeme súradnicu $t$ rovnú polovici vzdialenosti stredov kruhov: zatmenie začína na $t = R$ a úplné zatmenie je pre $t =0$. $$ S_\text{zakr} = 2\left( \pi R^2 \frac{2\varphi(t)}{2\pi} - t \sqrt{R^2 - t^2}\right) \quad \text{, kde } \cos\varphi = \frac tR. $$ Zakrytá plocha Slnka je dvojnásobok plochy kruhového odseku:

Súradnica $t$ a plocha kruhového odseku. Súradnica $t$ a plocha kruhového odseku.

Vzdialenosť $t$ sa mení lineárne s časom, a teda priemerný pokles výkonu je rovný priemernému relatívnemu zakrytiu. To stačí navyše počítať od začiatku zatmenia po úplné zatmenie, keďže odkrytie prebieha symetricky.

Priemerný percentuálny pokles výkonu je teda $$ \frac 1R \int_{0}^R \frac{S_\text{zakr}}{\pi R^2} dt = \frac{2}{\pi R} \int_0^R \left(\varphi(t) - \frac tR \sqrt{1-\frac{t^2}{R^2}}\right) dt. $$ V prvom člene použijeme vzťah medzi $t$ a $\varphi$, teda $\sin\varphi d\varphi = dt/R$, aby sme dostali1 $$\text{1. člen} = \frac 2\pi \int_{0}^{\pi/2} \underbrace{\varphi \sin\varphi}_{ = (-\varphi \cos\varphi)’ + \cos \varphi} d\varphi = \frac 2\pi \left( [-\varphi \cos\varphi]_0^{\pi/2} + \int_0^{\pi/2} \cos \varphi d \varphi\right) = \frac 2 \pi( 0 + 1 )= \frac 2 \pi. $$ Druhý člen dáva $$\text{2. člen} = -\frac 2\pi \underbrace{\int_0^R \frac tR \sqrt{1-\frac{t^2}{R^2}} \frac {dt}R}_{\text{nech } s = t^2/R^2} = -\frac 1 \pi \int_0^1 \sqrt{1-s} ds = \frac 2 {3\pi} [(1-s)^{3/2}]_0^1 = -\frac{2}{3\pi}.$$ Spolu je teda priemerný relatívny pokles výkonu za zatmenie $$ \frac{2}\pi (1-1/3) = \frac{4}{3\pi} \doteq 42\%. $$ Zatmenie, ktoré trvá dve a pol hodiny, teda spôsobí pokles výkonu na jednotku plochy o približne $5.2 \mathrm{MJ}/\mathrm{m}^2$, ak uvažujeme dopadajúci výkon Slnka rovný $1360 \mathrm{W}/\mathrm{m}^2$ (skutočný výkon je okolo $1000\text{-}1100 \mathrm{W}/\mathrm{m}^2$ kvôli rozptylu v atmosfére, a teda aj pokles výkonu bude nižší).

V \textbf{komplikovanejšom modeli} sa mení intenzita s polohou v slnečnom disku podľa $$I(\psi) = I_0 ( 1-u + u \cos\psi ),$$ kde uhol $\psi$ závisí od vzdialenosti od stredu disku ako $\sin \psi = r/R$.

Súvis uhla $\psi$ a vzdialenosti od stredu slnečného disku $r$. Súvis uhla $\psi$ a vzdialenosti od stredu slnečného disku $r$.

Uvažujme teda bod disku na polárnych súradniciach $(r, \phi)$ a určime, ako dlho bude tento bod zakrytý. Tento čas závisí lineárne od $b(r, \phi)$, šírky Mesiaca v tejto výške, lebo presne táto časť Mesiaca zakrýva tento bod2.

Šírka Mesiaca $b$ zakrývajúca daný bod Slnka. Šírka Mesiaca $b$ zakrývajúca daný bod Slnka.

Táto šírka je rovná $$ b(r, \phi) = 2 \sqrt{R^2 - r^2 \sin^2 \phi} $$ a bod $(r, \phi)$ je zakrytý $\tfrac{b}{4R}$ času zatmenia. Celkový \textbf{stratený} výkon je integrál týchto príspevkov cez slnečný disk $$ I_\text{strata}(u) = \int_\text{kruh} I(\psi; u) \frac{b(r, \phi)}{4R} dS = \frac{I_0}{2} \int_\text{kruh} \left(1-u + u \sqrt{1-\frac{r^2}{R^2}} \right) \sqrt{1- \frac{r^2}{R^2} \sin^2 \phi} \,dS. $$ Prejdime do kartézskych súradníc normovaných $R$, teda $x = r/R \cos \phi$ a $y = r/R \sin \phi$, element plochy je $R^2 dx dy$ $$ I_\text{strata}(u) = \frac{I_0 (1-u) R^2}{2} \int_\text{kruh} \sqrt{1-y^2} dx dy + \frac{I_0 u R^2}{2} \int_\text{kruh} \sqrt{1-y^2} \sqrt{1-x^2 -y^2} dx dy. $$ Prvý integrál je jednoducho $$ \int_\text{kruh} \sqrt{1-y^2} dx dy = \int_{y=-1}^{y=1} \sqrt{1-y^2}\int_{-\sqrt{1-y^2}}^{\sqrt{1-y^2}} dx dy = 2\int_{-1}^1 (1-y^2) dy = 2[y-y^3/3]_{-1}^1 = \frac 83. $$ Druhý integrál je komplikovanejší, ale dá sa vypočítať nasledovne $$ \int_\text{kruh} \sqrt{1-y^2} \sqrt{1-x^2 -y^2} dx dy = \int_{y=-1}^{y=1} \sqrt{1-y^2} \int_{-\sqrt{1-y^2}}^{\sqrt{1-y^2}} \sqrt{(1-y^2) \left(1-\frac{x^2}{1-y^2}\right)} dx dy = \dots $$ Vo vnútornom integrále použijeme substitúciu $\xi = x/\sqrt{1-y^2}$ a dostaneme $$\dots = \int_{y=-1}^{y=1} (1-y^2)^{3/2} \underbrace{\int_{-1}^1 \sqrt{1-\xi^2} d\xi}_{\text{plocha polkruhu}} dy = \frac\pi2 \underbrace{\int_{-1}^{1} (1-y^2)^{3/2} dy}_{y = \sin \alpha} = \frac\pi2 \int_{-\pi/2}^{\pi/2} \cos^4\alpha d\alpha = \frac \pi 2 \int_{-\pi/2}^{\pi/2} \frac{(e^{4i\alpha} + e^{-i\alpha})^4}{16}. $$ Rozviňme $$ (e^{i\alpha} + e^{-i\alpha})^4 = e^{4i\alpha} + 4 e^{2i\alpha} + 6 + 4e^{-2i\alpha} + e^{-4i\alpha} $$ a všimnime si, že len prostredný člen má nenulový integrál cez $[-\pi/2, \pi/2]$, keďže ostatné členy majú periódu rovnú alebo zlomok dĺžky tohoto intervalu. Pre celkový integrál, ktorý počítame, teda dostávame $$ \dots = \frac \pi 2 \frac {6}{16} \pi = \frac {3\pi^2}{16}.$$ Spočítali sme teda celkový stratený výkon $$ I_\text{strata}(u) = \frac{I_0 R^2}{2} \left( (1-u) \frac 83 + u \frac{3\pi^2}{16} \right). $$ Tento výkon musíme podeliť celkovým priemerným výkonom v prípade nezakrytého Slnka $$I_\text{ref}(u) = I_\text{tot}(u) = \int_\text{kruh} I(\psi; u) dS = I_0 \int_\text{kruh} \left(1-u + u \sqrt{1-\frac{r^2}{R^2}} \right) dS. $$ Jediný netriviálny integrál je $$ \int_\text{kruh} \sqrt{1-\frac{r^2}{R^2}} \, dS = 2\pi \int_0^R \sqrt{1-\frac{r^2}{R^2}} rdr =\frac{2\pi R^2}{3} [-(1-r^2/R^2)]_0^R = \frac{2\pi R^2}{3} $$. Spolu teda dostaneme $$I_\text{ref}(u) = I_0 \pi R^2 \left( (1-u) + u\frac{2}{3} \right). $$ Priemerný pokles výkonu je teda pomer $$ \frac{I_\text{strata}(u)}{I_\text{ref}(u)} = \frac{(1-u)\frac{4}{3\pi}+ u \frac{3\pi}{32} }{1-u + u \frac 23} ,$$ čo sa mení od hodnoty $42\%$ pre $u=0$ po $44\%$ pre $u=1$.


  1. Znamienka pri substitúciách ignorujeme a poradie hraníc integrálu vždy volíme tak, aby sme dostali kladný výsledok. 

  2. Takúto šikovnú parametrizáciu si všimol Richard Dudek. 

Pre odovzdávanie sa musíš prihlásiť.
Trojsten

FX zastrešuje občianske združenie Trojsten.

Kontakt
Ďalšie projekty