Diferenciální operátory ve sférických souřadnicích
Úloha číslo: 262
Sférické souřadnice (r, θ, φ) bodu P jsou definovány následovně: r je vzdálenost od počátku vztažné soustavy. Souřadnice θ je orientovaný úhel, který svírá spojnice bodu P a počátku s osou z, přičemž jej měříme od osy z. A souřadnice φ je úhel, který měříme od kladné poloosy x do přímky, kterou určuje počátek a průmět bodu P do roviny xy.
Vztah sférických a kartézských souřadnic popisují rovnice
\[x = r\,\sin\theta\,\cos\varphi,\] \[y = r\,\sin\theta\,\sin\varphi,\] \[z = r\,\cos\theta.\]Ve sférických souřadnicích definujeme tři jednotkové vektory
\[\vec r_0, \vec \theta_0, \vec \varphi_0\]mířící ve „směru růstu jednotlivých souřadnic“. Jestliže označíme jednotkové vektory v kartézské soustavě mířící po pořadě ve směru kladných poloos x, y, z postupně
\[\vec x_0, \vec y_0, \vec z_0\]potom formálně můžeme jednotkové vektory ve sférických souřadnicích definovat vztahy
\[\vec r_0 = \sin\theta\,\cos\varphi\,\vec x_0 + \sin\theta\,\sin\varphi\,\vec y_0 + \cos\theta\,\vec z_0,\] \[\vec \theta_0 = \cos\theta\,\cos\varphi\,\vec x_0 + \cos\theta\,\sin\varphi\,\vec y_0 - \sin\theta\,\vec z_0,\] \[\vec \varphi_0 = -\sin\varphi\,\vec x_0 + \cos\varphi\,\vec y_0.\]
(a) Ověřte, že vektory \[\vec r_0,\ \vec\theta_0,\ \vec\varphi_0\] mají jednotkovou velikost a jsou navzájem kolmé.
(b) Ukažte, že gradient skalárního pole T má ve sférických souřadnicích tvar
\[\nabla T = \frac{\partial T}{\partial r}\vec r_0 + \frac{1}{r}\frac{\partial T}{\partial \theta}\vec \theta_0 + \frac{1}{r\sin\theta}\frac{\partial T}{\partial \varphi}\vec \varphi_0.\](c) Ukažte, že divergence vektorového pole \(\vec v\) má tvar
\[\nabla\cdot\vec v = \frac{1}{r^2}\frac{\partial}{\partial r}(r^2v_r)+\frac{1}{r\sin\theta}\frac{\partial}{\partial\theta}(\sin\theta\,v_\theta)+\frac{1}{r\sin\theta}\frac{\partial v_\varphi}{\partial \varphi}.\](d) Ukažte, že rotace vektorového pole \(\vec v\) má tvar
\[\nabla\times\vec v = \frac{1}{r\sin\theta}\left[\frac{\partial}{\partial \theta}(\sin\theta\,v_\varphi)-\frac{\partial v_\theta}{\partial \varphi}\right]\,\vec r_0 + \frac{1}{r}\left[\frac{1}{\sin\theta}\frac{\partial v_r}{\partial \varphi\varphi}-\frac{\partial}{\partial r}(rv_\varphi)\right]\,\vec \theta_0 + \frac{1}{r}\left[\frac{\partial}{\partial r}(rv_\theta)-\frac{\partial v_r}{\partial \theta}\right]\,\vec \varphi_0.\]Rozbor
Úlohu budeme řešit přímým výpočtem.
S výjimkou části (a), kterou lze vyřešit přímočarým počítáním, máme v ostatních částech úlohy převádět různé diferenciální operátory. K tomu potřebujeme zvládnout:
- vyjádření parciální derivace podle kartézských souřadnic pomocí parciálních derivací podle sférických souřadnic,
- výpočet složek vektoru ve sférických souřadnicích z jeho složek v souřadnicích kartézských.
Komentář – k vyjádření „sférických“ jednotkových vektorů
Asi nejnázornější cestou, jak tyto vektory vyjádřit, je využít geometrie. Tady poukážeme na způsob, jak k nim dojít na základě ryze početních úvah.
V jakýchkoli souřadnicích definujeme základní jednotkové vektory tak, aby mířili ve směru přírůstku jednotlivých souřadnic. V kartézské soustavě tedy míří ve směru kladných poloos x, y, z.
Má-li být stejná konvence dodržena také u sférických souřadnic, pak musí platit, že jednotkové vektory \(\vec r_0,\,\vec\theta_0,\,\vec\varphi_0\) musí mít směr přírůstku souřadnic r, θ, φ. To jest, popořadě musí mít směr vektorů
\[\frac{\partial\vec r}{\partial r}, \quad \frac{\partial\vec r}{\partial \theta}, \quad \frac{\partial\vec r}{\partial \varphi},\]kde \(\vec r\) je polohový vektor. V kartézských souřadnicích samozřejmě platí, že
\[\vec r = (x,y,z),\]neboli v jiném zápisu
\[\vec r = x\vec x_0 + y\vec y_0 + z\vec z_0.\]Nyní využijeme definičních vztahů sférických souřadnic a dosazením za x, y, z dostaneme, že
\[\vec r = r\,\sin\,\theta\,\cos\,\varphi\, \vec x_0 + r\,\sin\,\theta\,\sin\,\varphi\,\vec y_0 + r\,\cos\,\theta\,\vec z_0.\]Nyní vypočteme hledané derivace (uvědomte si, že vektory \(\vec x_0, \, \vec y_0, \, \vec z_0\) jsou konstantní)
\[\frac{\partial\vec r}{\partial r} = \sin\,\theta\,\cos\,\varphi\, \vec x_0 + \sin\,\theta\,\sin\,\varphi\,\vec y_0 + \cos\,\theta\,\vec z_0,\] \[\frac{\partial\vec r}{\partial \theta} = r\,\cos\,\theta\,\cos\,\varphi\, \vec x_0 + r\,\cos\,\theta\,\sin\,\varphi\,\vec y_0 - r\,\sin\,\theta\,\vec z_0,\] \[\frac{\partial\vec r}{\partial \varphi} = -r\,\sin\,\theta\,\sin\,\varphi\, \vec x_0 + r\,\sin\,\theta\,\cos\,\varphi\,\vec y_0.\]Nyní stačí tyto tři vektory podělit jejich velikostí. Protože platí (podle Pythagorovy věty pro vektorovou algebru, neboť „kartézské“ vektory jsou jednotkové a na sebe kolmé)
\[\left|\frac{\partial\vec r}{\partial r}\right|^2 = \sin^2\theta\,\cos^2\varphi + \sin^2\theta\,\sin^2\varphi + \cos^2\theta = 1,\] \[\left|\frac{\partial\vec r}{\partial \theta}\right|^2 = r^2\,\cos^2\theta\,\cos^2\varphi + r^2\,\cos^2\theta\,\sin^2\varphi + r^2\,\sin^2\theta = r^2,\] \[\left|\frac{\partial\vec r}{\partial \varphi}\right|^2 = r^2\,\sin^2\theta\,\sin^2\varphi + r^2\,\sin^2\theta\,\cos^2\varphi = r^2\sin^2\theta,\]dostáváme, že
\[\vec r_0 = \frac{\frac{\partial\vec r}{\partial r}}{|\frac{\partial\vec r}{\partial r}|} = \sin\,\theta\,\cos\,\varphi\, \vec x_0 + \sin\,\theta\,\sin\,\varphi\,\vec y_0 + \cos\,\theta\,\vec z_0,\] \[\vec \theta_0 = \frac{\frac{\partial\vec r}{\partial \theta}}{|\frac{\partial\vec r}{\partial \theta}|} = \cos\,\theta\,\cos\,\varphi\, \vec x_0 + \cos\,\theta\,\sin\,\varphi\,\vec y_0 - \sin\,\theta\,\vec z_0,\] \[\vec \varphi_0 = \frac{\frac{\partial\vec r}{\partial \varphi}}{|\frac{\partial\vec r}{\partial \varphi}|} = -\sin\,\varphi\, \vec x_0 + \cos\,\varphi\,\vec y_0.\]Řešení části (a)
Nejprve ověříme, že vektory \(\vec r_0, \ \vec\theta_0, \ \vec\varphi_0\) jsou jednotkové. K tomu využijeme, že „kartézské“ jednotkové vektory \(\vec x_0,\ \vec y_0,\ \vec z_0\) jsou jednotkové a na sebe kolmé. To znamená, že
\[\vec x_0\cdot \vec x_0 = 1, \quad \vec y_0\cdot\vec y_0 = 1, \quad \vec z_0\cdot \vec z_0 = 1\] \[\vec x_0\cdot \vec y_0 = 0, \quad \vec x_0\cdot\vec z_0 = 0, \quad \vec y_0\cdot \vec z_0 = 0.\]Díky těmto vztahům dostáváme (podle Pythagorovy věty pro vektorovou algebru), že
\[|\vec r_0|^2 = \vec r_0\cdot\vec r_0 =\] \[=(\sin\,\theta\,\cos\,\varphi)^2 + (\sin\,\theta\,\sin\,\varphi)^2 + (\cos\,\theta)^2 = \sin^2\theta\,(\sin^2\varphi+\cos^2\varphi) + \cos^2\theta = \sin^2\theta + \cos^2\theta = 1.\]Obdobně platí
\[|\vec \theta_0|^2 = \vec \theta_0\cdot\vec \theta_0=\] \[= (\cos\,\theta\,\cos\,\varphi)^2 + (\cos\,\theta\,\sin\,\varphi)^2 + (\sin\,\theta)^2 = \cos^2\theta\,(\sin^2\varphi+\cos^2\varphi) + \sin^2\theta = \cos^2\theta + \sin^2\theta = 1.\]A zcela analogicky také
\[|\vec \varphi_0|^2 = \vec \varphi_0\cdot\vec \varphi_0 = (-\sin\,\theta)^2 + (\cos\,\theta)^2 = 1.\]Nyní ověříme kolmost dvojic různých vektorů \(\vec r_0,\ \vec \theta_0, \ \vec\varphi_0\). K tomu stačí spočítat jejich skalární součin
\[\vec r_0 \cdot \vec \theta_0 = (\sin\,\theta\,\cos\,\theta\,\cos^2\varphi)+(\sin\,\theta\,\cos\,\theta\,\sin^2\varphi)-(\sin\,\theta\,\cos\,\theta) = \sin\,\theta\,\cos\,\theta\,(\cos^2\varphi+\sin^2\varphi-1) = 0.\]Obdobně je
\[\vec r_0 \cdot \vec \varphi_0 = -(\sin\,\theta\,\cos\,\varphi\,\sin\,\varphi)+(\sin\,\theta\,\sin\,\varphi\,\cos\,\varphi) + 0 = 0\]a také platí
\[\vec \theta_0 \cdot \vec \varphi_0 = -(\cos\,\theta\,\cos\,\varphi\,\sin\,\varphi)+(\cos\,\theta\,\sin\,\varphi\,\cos\,\varphi) + 0 = 0.\]Vyjádření kartézských jednotkových vektorů pomocí sférických
Ukažte, že platí
\[\vec x_0 = \sin\,\theta\,\cos\,\varphi\,\vec r_0 + \cos\,\theta\,\cos\,\varphi\,\vec \theta_0 - \sin\,\varphi\,\vec\varphi_0\] \[\vec y_0 = \sin\,\theta\,\sin\,\varphi\,\vec r_0 + \cos\,\theta\,\sin\,\varphi\,\vec \theta_0 + \cos\,\varphi\,\vec\varphi_0\] \[\vec z_0 = \cos\,\theta\,\vec r_0 -\sin\,\theta\,\vec\theta_0.\]Sférické souřadnice vyjádřené pomocí kartézských souřadnic
Ukažte, že platí
\[r = \sqrt{x^2+y^2+z^2}\] \[\theta = \arccos\,\frac{z}{r}\] \[\varphi = \arctan\,\frac{y}{x}.\]Výpočet parciálních derivací — shrnutí vztahů
Pro další účely vypočtěte postupně parciální derivace proměnných r, θ, φ podle proměnných x, y, z (celkem devět derivací). K tomu použijte vyjádření r, θ, φ jako funkcí kartézských souřadnic.
Nakonec tyto parciální derivace vyjádřete pomocí souřadnic sférických, to jest pomocí proměnných r,θ,φ.
To jest, ukažte, že platí
\[\frac{\partial r}{\partial x} = \sin\,\theta\,\cos\,\varphi\tag{a1}\] \[\frac{\partial r}{\partial y} = \sin\,\theta\,\sin\,\varphi\tag{a2}\] \[\frac{\partial r}{\partial z} = \cos\,\theta\tag{a3}\] \[ \frac{\partial \theta}{\partial x} = \frac{1}{r}\cos\varphi\cos\theta\tag{b1}\] \[ \frac{\partial \theta}{\partial y} = \frac{1}{r}\sin\varphi\cos\theta\tag{b2}\] \[ \frac{\partial \theta}{\partial z} = -\frac{1}{r}\sin\,\theta\tag{b3}\] \[\frac{\partial \varphi}{\partial x} = -\frac{1}{r\sin\theta}\sin\varphi\tag{c1}\] \[\frac{\partial \varphi}{\partial y} = \frac{1}{r\sin\theta}\cos\varphi\tag{c2}\] \[\frac{\partial \varphi}{\partial z} = 0.\tag{c3}\]Nápověda k části (b)
Pro přepočítání parciálních derivací
\[\frac{\partial T}{\partial x}, \quad \frac{\partial T}{\partial y}, \quad \frac{\partial T}{\partial z},\]použijte řetízkové pravidlo:
\[\frac{\partial T}{\partial x} = \frac{\partial T}{\partial r}\frac{\partial r}{\partial x} + \frac{\partial T}{\partial \theta}\frac{\partial \theta}{\partial x} + \frac{\partial T}{\partial \varphi}\frac{\partial \varphi}{\partial x}\] \[\frac{\partial T}{\partial y} = \frac{\partial T}{\partial r}\frac{\partial r}{\partial y} + \frac{\partial T}{\partial \theta}\frac{\partial \theta}{\partial y} + \frac{\partial T}{\partial \varphi}\frac{\partial \varphi}{\partial y}\] \[\frac{\partial T}{\partial z} = \frac{\partial T}{\partial r}\frac{\partial r}{\partial z} + \frac{\partial T}{\partial \theta}\frac{\partial \theta}{\partial z} + \frac{\partial T}{\partial \varphi}\frac{\partial \varphi}{\partial z}\]a předchozích výsledků.
Řešení části (b)
Nechť T je skalární pole. V kartézských souřadnicích je jeho gradient definován jako vektorové pole
\[\nabla T = \left(\frac{\partial T}{\partial x},\ \frac{\partial T}{\partial y},\ \frac{\partial T}{\partial z}\right),\]nebo můžeme také psát, že
\[\nabla T = \frac{\partial T}{\partial x}\vec x_0 + \frac{\partial T}{\partial y}\vec y_0 + \frac{\partial T}{\partial z}\vec z_0.\]V tomto vyjádření použijeme řetízkové pravidlo
\[\nabla T = \frac{\partial T}{\partial x}\vec x_0 + \frac{\partial T}{\partial y}\vec y_0 + \frac{\partial T}{\partial z}\vec z_0 =\left(\frac{\partial T}{\partial r}\frac{\partial r}{\partial x} + \frac{\partial T}{\partial \theta}\frac{\partial \theta}{\partial x} + \frac{\partial T}{\partial \varphi}\frac{\partial \varphi}{\partial x}\right)\,\vec x_0 +\] \[+ \left(\frac{\partial T}{\partial r}\frac{\partial r}{\partial y} + \frac{\partial T}{\partial \theta}\frac{\partial \theta}{\partial y} + \frac{\partial T}{\partial \varphi}\frac{\partial \varphi}{\partial y}\right)\, \vec y_0 + \left(\frac{\partial T}{\partial r}\frac{\partial r}{\partial z} + \frac{\partial T}{\partial \theta}\frac{\partial \theta}{\partial z} + \frac{\partial T}{\partial \varphi}\frac{\partial \varphi}{\partial z}\right)\, \vec z_0 =\]a posléze dosadíme za parciální derivace a jednotkové vektory z nápověd výše
\[= \left(\frac{\partial T}{\partial r}\,\sin\,\theta\,\cos\,\varphi + \frac{\partial T}{\partial \theta}\,\frac{1}{r}\,\cos\,\theta\,\cos\,\varphi - \frac{\partial T}{\partial \varphi}\,\frac{1}{r\sin\theta}\,\sin\,\varphi\right) \,\cdot\,\left(\sin\,\theta\,\cos\,\varphi\,\vec r_0 + \cos\,\theta\,\cos\,\varphi\,\vec \theta_0 - \sin\,\varphi\,\vec\varphi_0\right) +\] \[+ \left(\frac{\partial T}{\partial r}\,\sin\,\theta\,\sin\,\varphi + \frac{\partial T}{\partial \theta}\,\frac{1}{r}\,\cos\,\theta\,\sin\,\varphi + \frac{\partial T}{\partial \varphi}\,\frac{1}{r\sin\theta}\,\cos\,\varphi\right) \,\cdot\,\left(\sin\,\theta\,\sin\,\varphi\,\vec r_0 + \cos\,\theta\,\sin\,\varphi\,\vec \theta_0 + \cos\,\varphi\,\vec\varphi_0\right) +\] \[+ \left(\frac{\partial T}{\partial r}\,\cos\,\theta - \frac{\partial T}{\partial \theta}\,\frac{1}{r}\,\sin\,\theta\right)\,\cdot\,\left(\cos\,\theta\,\vec r_0 -\sin\,\theta\,\vec\theta_0\right)\]a po úpravách dostaneme hledaný vztah
\[\nabla T = \frac{\partial T}{\partial r}\vec r_0 + \frac{1}{r}\frac{\partial T}{\partial \theta}\vec \theta_0 + \frac{1}{r\sin\theta}\frac{\partial T}{\partial \varphi}\vec \varphi_0.\]Nápověda k části (c)
Divergence vektorového pole \(\vec v\) je v kartézských souřadnicích definována vztahem
\[\nabla\cdot\vec v = \frac{\partial v_x}{\partial x}+\frac{\partial v_y}{\partial y}+\frac{\partial v_z}{\partial z}.\]Nejprve použijeme řetízkového pravidla. Platí, že
\[\frac{\partial v_x}{\partial x} = \frac{\partial v_x}{\partial r}\frac{\partial r}{\partial x} + \frac{\partial v_x}{\partial \theta}\frac{\partial \theta}{\partial x} + \frac{\partial v_x}{\partial \varphi}\frac{\partial \varphi}{\partial x}\] \[\frac{\partial v_y}{\partial y} = \frac{\partial v_y}{\partial r}\frac{\partial r}{\partial y} + \frac{\partial v_y}{\partial \theta}\frac{\partial \theta}{\partial y} + \frac{\partial v_y}{\partial \varphi}\frac{\partial \varphi}{\partial y}\] \[\frac{\partial v_z}{\partial z} = \frac{\partial v_z}{\partial r}\frac{\partial r}{\partial z} + \frac{\partial v_z}{\partial \theta}\frac{\partial \theta}{\partial z} + \frac{\partial v_z}{\partial \varphi}\frac{\partial \varphi}{\partial z}.\]A poté ještě musíme kartézské složky vektoru \(\vec v\) pomocí „sférických složek“ a sférických souřadnic.
To provedeme tak, že do vyjádření
\[\vec v = v_r\vec r_0 + v_\theta\vec \theta_0 + v_\varphi\vec \varphi_0\]nahradíme vektory \(\vec r_0, \, \vec \theta_0, \, \vec \varphi_0\) jejich vyjádřením pomocí kartézských souřadnic.
Řešení části (c)
Divergence vektorového pole \(\vec v\) je v kartézských souřadnicích definována vztahem
\[\nabla\cdot\vec v = \frac{\partial v_x}{\partial x}+\frac{\partial v_y}{\partial y}+\frac{\partial v_z}{\partial z}.\]Nejprve použijeme řetízkového pravidla. Platí, že
\[\frac{\partial v_x}{\partial x} = \frac{\partial v_x}{\partial r}\frac{\partial r}{\partial x} + \frac{\partial v_x}{\partial \theta}\frac{\partial \theta}{\partial x} + \frac{\partial v_x}{\partial \varphi}\frac{\partial \varphi}{\partial x}\] \[\frac{\partial v_y}{\partial y} = \frac{\partial v_y}{\partial r}\frac{\partial r}{\partial y} + \frac{\partial v_y}{\partial \theta}\frac{\partial \theta}{\partial y} + \frac{\partial v_y}{\partial \varphi}\frac{\partial \varphi}{\partial y}\] \[\frac{\partial v_z}{\partial z} = \frac{\partial v_z}{\partial r}\frac{\partial r}{\partial z} + \frac{\partial v_z}{\partial \theta}\frac{\partial \theta}{\partial z} + \frac{\partial v_z}{\partial \varphi}\frac{\partial \varphi}{\partial z}.\]A poté ještě musíme kartézské složky vektoru \(\vec v\) vyjádřit pomocí „složek sférických“ – tedy vx, vy a vz vyjádřit pomocí sférických souřadnic a složek vr, vθ a vφ.
Nejprve ve vyjádření
\[\vec v = v_r\vec r_0 + v_\theta\vec \theta_0 + v_\varphi\vec \varphi_0\]nahradíme sférické jednotkové vektory kartézskými podle vztahů uvedených v zadání úlohy
\[\vec r_0 = \sin\theta\,\cos\varphi\,\vec x_0 + \sin\theta\,\sin\varphi\,\vec y_0 + \cos\theta\,\vec z_0\] \[\vec \theta_0 = \cos\theta\,\cos\varphi\,\vec x_0 + \cos\theta\,\sin\varphi\,\vec y_0 - \sin\theta\,\vec z_0\] \[\vec \varphi_0 = -\sin\varphi\,\vec x_0 + \cos\varphi\,\vec y_0.\]Dostaneme tak, že
\[\vec v = \left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right)\vec x_0 +\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi + v_\varphi\cos\varphi\right)\vec y_0 + \] \[+\left(v_r\cos\theta -v_\theta\sin\theta\right)\vec z_0.\]Porovnáním se vztahem
\[\vec v = v_x\vec x_0 + v_y\vec y_0 + v_z\vec z_0\]dostaneme, že
\[v_x = v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\] \[v_y = v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi + v_\varphi\cos\varphi\] \[v_z = v_r\cos\theta -v_\theta\sin\theta.\]Kombinací řetízkového pravidla, parciálních derivací vypočtených v nápovědě výše a právě odvozených vztahů dostaneme:
\[\frac{\partial v_x}{\partial x} = \frac{\partial v_x}{\partial r}\frac{\partial r}{\partial x} + \frac{\partial v_x}{\partial \theta}\frac{\partial \theta}{\partial x} + \frac{\partial v_x}{\partial \varphi}\frac{\partial \varphi}{\partial x} =\] \[ = \frac{\partial}{\partial r}\left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right) \cdot (\sin\theta\cos\varphi) +\] \[+ \frac{\partial}{\partial \theta}\left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right) \cdot \frac{1}{r}\cos\theta\cos\varphi +\] \[+ \frac{\partial}{\partial \varphi}\left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right)\cdot \left(-\frac{1}{r\sin\theta}\sin\varphi\right) =\]a nyní je potřeba provést derivování
\[= \frac{\partial v_r}{\partial r}\sin^2\theta\cos^2\varphi + \frac{\partial v_\theta}{\partial r}\sin\theta\,\cos\theta\,\cos^2\varphi - \frac{\partial v_\varphi}{\partial r}\sin\theta\,\sin\varphi\,\cos\varphi +\] \[+ \frac{1}{r}\frac{\partial v_r}{\partial\theta}\cos\theta\,\sin\theta\,\cos^2\varphi + \frac{1}{r}v_r\cos^2\theta\cos^2\varphi +\] \[+ \frac{1}{r}\frac{\partial v_\theta}{\partial \theta}\cos^2\theta\cos^2\varphi - \frac{1}{r}v_\theta\sin\theta\,\cos\theta\,\cos^2\varphi -\frac{1}{r}\frac{\partial v_\varphi}{\partial\theta}\cos\theta\,\sin\varphi\,\cos\varphi - \] \[- \frac{1}{r}\frac{\partial v_r}{\partial\varphi}\cos\varphi\sin\varphi + \frac{1}{r}v_r\sin^2\varphi - \frac{1}{r\sin\theta}\frac{\partial v_\theta}{\partial\varphi}\cos\theta\,\cos\varphi\,\sin\varphi +\] \[+ \frac{1}{r\sin\theta}v_\theta\cos\theta\,\sin^2\varphi + \frac{1}{r\sin\theta}\frac{\partial v_\varphi}{\partial\varphi}\sin^2\varphi + \frac{1}{r\sin\theta}\,\cos\varphi\,\sin\varphi.\]Obdobně vypočteme y-ový člen:
\[\frac{\partial v_y}{\partial y}= \frac{\partial v_y}{\partial r}\frac{\partial r}{\partial y} + \frac{\partial v_y}{\partial \theta}\frac{\partial \theta}{\partial y} + \frac{\partial v_y}{\partial \varphi}\frac{\partial \varphi}{\partial y} =\] \[ = \frac{\partial}{\partial r}\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi+v_\varphi\cos\varphi\right) \cdot (\sin\theta\sin\varphi) +\] \[+ \frac{\partial}{\partial \theta}\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi+v_\varphi\cos\varphi\right) \cdot \frac{1}{r}\cos\theta\sin\varphi +\] \[+ \frac{\partial}{\partial \varphi}\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi+v_\varphi\cos\varphi\right)\cdot \left(\frac{1}{r\sin\theta}\cos\varphi\right) =\]a nyní je potřeba provést derivování
\[= \frac{\partial v_r}{\partial r}\sin^2\theta\sin^2\varphi + \frac{\partial v_\theta}{\partial r}\sin\theta\,\cos\theta\,\sin^2\varphi + \frac{\partial v_\varphi}{\partial r}\sin\theta\,\sin\varphi\,\cos\varphi +\] \[+ \frac{1}{r}\frac{\partial v_r}{\partial\theta}\cos\theta\,\sin\theta\,\sin^2\varphi + \frac{1}{r}v_r\cos^2\theta\sin^2\varphi +\] \[+ \frac{1}{r}\frac{\partial v_\theta}{\partial \theta}\cos^2\theta\sin^2\varphi - \frac{1}{r}v_\theta\sin\theta\,\cos\theta\,\sin^2\varphi + \frac{1}{r}\frac{\partial v_\varphi}{\partial\theta}\cos\theta\,\sin\varphi\,\cos\varphi +\] \[+ \frac{1}{r}\frac{\partial v_r}{\partial\varphi}\sin\varphi\,\cos\varphi + \frac{1}{r}v_r\cos^2\varphi + \frac{1}{r\sin\theta}\frac{\partial v_\theta}{\partial\varphi}\cos\theta\,\sin\varphi\,\cos\varphi +\] \[+ \frac{1}{r\sin\theta}v_\theta\cos\theta\,\cos^2\varphi + \frac{1}{r\sin\theta}\frac{\partial v_\varphi}{\partial\varphi}\cos^2\varphi - \frac{1}{r\sin\theta}v_\varphi\,\sin\varphi\,\cos\varphi .\]A nakonec vypočteme z-ový člen:
\[\frac{\partial v_z}{\partial z}= \frac{\partial v_z}{\partial r}\frac{\partial r}{\partial z} + \frac{\partial v_z}{\partial \theta}\frac{\partial \theta}{\partial z} + \frac{\partial v_z}{\partial \varphi}\frac{\partial \varphi}{\partial z} = \frac{\partial}{\partial r}\left(v_r\cos\theta -v_\theta\sin\theta\right) \cdot \cos\theta + \frac{\partial}{\partial \theta}\left(v_r\cos\theta -v_\theta\sin\theta\right) \cdot \left(-\frac{1}{r}\sin\theta\right) =\]a nyní je potřeba provést derivování
\[= \frac{\partial v_r}{\partial r}\cos^2\theta - \frac{\partial v_\theta}{\partial r}\sin\theta\,\cos\theta + \frac{1}{r}\frac{\partial v_r}{\partial\theta}\cos\theta\,\sin\theta + \frac{1}{r}v_r\sin^2\theta + \frac{1}{r}\frac{\partial v_\theta}{\partial \theta}\sin^2\theta + \frac{1}{r}v_\theta\sin\theta\,\cos\theta.\]Sečtením všech předchozích vztahů dostaneme
\[\nabla \cdot \vec v = \frac{\partial v_r}{\partial r} + \frac{2}{r}v_r + \frac{1}{r}\frac{\partial v_\theta}{\partial\theta} + \frac{1}{r\sin\theta}v_\theta\cos\theta + \frac{1}{r\sin\theta}\frac{\partial v_\varphi}{\partial\varphi}.\]S přihlédnutím ke vztahům
\[\frac{1}{r^2}\frac{\partial}{\partial r}(r^2v_r) = \frac{2}{r}v_r + \frac{\partial v_r}{\partial r}\]a
\[\frac{\partial}{\partial\theta}(v_\theta\sin\theta) = v_\theta\cos\theta + \frac{\partial v_\theta}{\partial\theta}\sin\theta\]ihned dostáváme hledaný vztah.
Nápověda k části (d)
Rotace v kartézských souřadnicích je vektor, který značíme \(\nabla\times\vec v\) a jehož složky jsou definovány následovně:
\[(\nabla\times\vec v)_x = \frac{\partial v_z}{\partial y} - \frac{\partial v_y}{\partial z}\] \[(\nabla\times\vec v)_y = \frac{\partial v_x}{\partial z} - \frac{\partial v_z}{\partial x}\] \[(\nabla\times\vec v)_z = \frac{\partial v_y}{\partial x} - \frac{\partial v_x}{\partial y}.\]Naším úkolem je najít složky tohoto vektoru ve sférických souřadnicích a jednotlivé členy vyjádřit pomocí složek vektoru \(\vec v\) ve sférických souřadnicích – ot jest pomocí vr, vθ a vφ – a jejich derivací podle sférických souřadnic.
Začneme tak, že složky \((\nabla\times\vec v)_{x,y,z}\) vyjádříme pomocí sférických souřadnic vektoru \(\vec v\) a jejich derivací podle sférických souřadnic. K tomu použijeme řetízkové pravidlo a již výše spočtené parciální derivace.
Poté použijeme následující postup. Jestliže libovolný vektor \(\vec u\) má souřadnice (ux, uy, uz), pak jej lze psát ve tvaru
\[\vec u = u_x\vec x_0 + u_y\vec y_0 + u_z\vec z_0\]a dosazením za jednotkové kartézské vektory jejich vyjádřením pomocí sférických jednotkových vektorů umíme získat souřadnice vektoru \(\vec u\) ve sférických souřadnicích.
Řešení části (d)
Složky vektoru rotace mají tvar
\[(\nabla\times\vec v)_x = \frac{\partial v_z}{\partial y} - \frac{\partial v_y}{\partial z}\] \[(\nabla\times\vec v)_y = \frac{\partial v_x}{\partial z} - \frac{\partial v_z}{\partial x}\] \[(\nabla\times\vec v)_z = \frac{\partial v_y}{\partial x} - \frac{\partial v_x}{\partial y}.\]Tyto parciální derivace vyjádříme pomocí sférických souřadnic.
Z předchozí části (c) víme, že platí
\[ v_x = v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\] \[v_y = v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi + v_\varphi\cos\varphi\] \[v_z = v_r\cos\theta -v_\theta\sin\theta.\]z-ová složka rotace
Podle řetízkového pravidla pak máme
\[\frac{\partial v_x}{\partial y} = \frac{\partial v_x}{\partial r}\frac{\partial r}{\partial y} + \frac{\partial v_x}{\partial \theta}\frac{\partial \theta}{\partial y} + \frac{\partial v_x}{\partial \varphi}\frac{\partial \varphi}{\partial y} = ,\]nyní dosadíme z předchozích vztahů pro složku vx a příslušné parciální derivace (viz nápověda výše)
\[= \frac{\partial}{\partial r}\left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right)\cdot(\sin\theta\,\sin\varphi) +\] \[+ \frac{\partial}{\partial \theta}\left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right)\cdot\left(\frac{1}{r}\sin\varphi\,\cos\theta\right) +\] \[+ \frac{\partial}{\partial \varphi}\left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right)\cdot\left(\frac{1}{r\sin\theta}\cos\varphi\right).\]Znovu podle řetízkového pravidla můžeme psát
\[\frac{\partial v_y}{\partial x} = \frac{\partial v_y}{\partial r}\frac{\partial r}{\partial x} + \frac{\partial v_y}{\partial \theta}\frac{\partial \theta}{\partial x} + \frac{\partial v_y}{\partial \varphi}\frac{\partial \varphi}{\partial x} = ,\]nyní dosadíme z předchozích vztahů pro složku vy a příslušné parciální derivace (viz nápověda výše)
\[= \frac{\partial}{\partial r}\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi+v_\varphi\cos\varphi\right)\cdot(\sin\theta\,\cos\varphi) +\] \[+ \frac{\partial}{\partial \theta}\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi+v_\varphi\cos\varphi\right)\cdot\left(\frac{1}{r}\cos\varphi\,\cos\theta\right) +\] \[+ \frac{\partial}{\partial \varphi}\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi+v_\varphi\cos\varphi\right)\cdot\left(-\frac{1}{r\sin\theta}\sin\varphi\right).\]Využitím předchozím vztahů dostaneme po úpravách, že
\[(\nabla\times \vec v)_z = \frac{\partial v_y}{\partial x}-\frac{\partial v_x}{\partial y} =\] \[= \frac{\partial v_\varphi}{\partial r}\sin\theta + \frac{1}{r}\frac{\partial v_\varphi}{\partial\theta}\cos\theta - \frac{1}{r}\frac{\partial v_r}{\partial\varphi} - \frac{1}{r\sin\theta}\frac{\partial v_r}{\partial\varphi}\cos\theta - \frac{1}{r\sin\theta}\frac{\partial v_\theta}{\partial\varphi}\cos\theta + \frac{1}{r\sin\theta}v_\varphi.\tag{R1}\]y-ová složka rotace
Podle řetízkového pravidla pak máme
\[\frac{\partial v_z}{\partial x} = \frac{\partial v_z}{\partial x}\frac{\partial r}{\partial x} + \frac{\partial v_z}{\partial \theta}\frac{\partial \theta}{\partial x} + \frac{\partial v_z}{\partial \varphi}\frac{\partial \varphi}{\partial x} = ,\]nyní dosadíme z předchozích vztahů pro složku vz a příslušné parciální derivace (viz nápověda výše)
\[= \frac{\partial}{\partial r}\left(v_r\cos\theta - v_\theta\sin\theta\right)\cdot(\sin\theta\,\cos\varphi) +\] \[+ \frac{\partial}{\partial \theta}\left(v_r\cos\theta - v_\theta\sin\theta\right)\cdot\left(\frac{1}{r}\cos\varphi\,\cos\theta\right) +\] \[+ \frac{\partial}{\partial \varphi}\left(v_r\cos\theta - v_\theta\sin\theta\right)\cdot\left(-\frac{1}{r\sin\theta}\sin\varphi\right).\]Znovu podle řetízkového pravidla můžeme psát
\[\frac{\partial v_x}{\partial z} = \frac{\partial v_x}{\partial r}\frac{\partial r}{\partial z} + \frac{\partial v_x}{\partial \theta}\frac{\partial \theta}{\partial z} + \frac{\partial v_x}{\partial \varphi}\frac{\partial \varphi}{\partial z} =, \]nyní dosadíme z předchozích vztahů pro složku vx a příslušné parciální derivace (viz nápověda výše)
\[= \frac{\partial}{\partial r}\left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right)\cdot(\cos\theta) +\] \[+ \frac{\partial}{\partial \theta}\left(v_r\sin\theta\,\cos\varphi + v_\theta\cos\theta\,\cos\varphi-v_\varphi\sin\varphi\right)\cdot\left(-\frac{1}{r}\,\sin\theta\right) + 0.\]Využitím obou předchozích vztahů dostaneme, že
\[(\nabla\times \vec v)_y = \frac{\partial v_x}{\partial z}-\frac{\partial v_z}{\partial x} =\] \[ =\frac{\partial v_\theta}{\partial r}\cos\varphi - \frac{\partial v_\varphi}{\partial r}\sin\varphi\cos\theta - \frac{1}{r}\frac{\partial v_r}{\partial \theta}\cos\varphi + \frac{1}{r}v_\theta\cos\varphi +\tag{R2}\] \[+ \frac{1}{r}\frac{\partial v_\varphi}{\partial \theta}\sin\theta\,\sin\varphi + \frac{1}{r\sin\theta}\frac{\partial v_r}{\partial \varphi}\cos\theta\sin\varphi - \frac{1}{r}\frac{\partial v_\theta}{\partial \varphi}\sin\varphi.\]x-ová složka rotace
Podle řetízkového pravidla pak máme
\[\frac{\partial v_y}{\partial z} = \frac{\partial v_y}{\partial r}\frac{\partial r}{\partial z} + \frac{\partial v_y}{\partial \theta}\frac{\partial \theta}{\partial z} + \frac{\partial v_y}{\partial \varphi}\frac{\partial \varphi}{\partial z} = ,\]nyní dosadíme z předchozích vztahů pro složku vz a příslušné parciální derivace (viz nápověda výše)
\[= \frac{\partial}{\partial r}\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi + v_\varphi\cos\varphi\right)\cdot(\cos\theta) +\] \[+ \frac{\partial}{\partial \theta}\left(v_r\sin\theta\,\sin\varphi + v_\theta\cos\theta\,\sin\varphi + v_\varphi\cos\varphi\right)\cdot\left(-\frac{1}{r}\sin\theta\right) + 0.\]Znovu podle řetízkového pravidla můžeme psát
\[\frac{\partial v_z}{\partial y} = \frac{\partial v_z}{\partial r}\frac{\partial r}{\partial y} + \frac{\partial v_z}{\partial \theta}\frac{\partial \theta}{\partial y} + \frac{\partial v_z}{\partial \varphi}\frac{\partial \varphi}{\partial y} = ,\]nyní dosadíme z předchozích vztahů pro složku vx a příslušné parciální derivace (viz nápověda výše)
\[= \frac{\partial}{\partial r}\left(v_r\cos\theta -v_\theta\sin\theta\right)\cdot(\sin\theta\,\sin\varphi) +\] \[+ \frac{\partial}{\partial \theta}\left(v_r\cos\theta -v_\theta\sin\theta\right)\cdot\left(\frac{1}{r}\sin\varphi\,\cos\theta\right) + \] \[+ \frac{\partial}{\partial \varphi}\left(v_r\cos\theta -v_\theta\sin\theta\right)\cdot\left(\frac{1}{r\sin\theta}\cos\varphi\right).\]Využitím obou předchozích vztahů dostaneme, že
\[(\nabla\times \vec v)_x = \frac{\partial v_z}{\partial y}-\frac{\partial v_y}{\partial z} =\] \[ = -\frac{\partial v_\theta}{\partial r}\sin\varphi -\frac{\partial v_\varphi}{\partial r}\cos\varphi\,\cos\theta + \frac{1}{r}\frac{\partial v_r}{\partial \theta}\sin\varphi -\frac{1}{r}v_\theta\sin\varphi +\tag{R3}\] \[+ \frac{1}{r}\frac{\partial v_\varphi}{\partial \theta}\sin\theta\,\cos\varphi + \frac{1}{r\sin\theta}\frac{\partial v_r}{\partial \varphi}\cos\theta\,\cos\varphi - \frac{1}{r}\frac{\partial v_\theta}{\partial \varphi}\cos\varphi.\]Přepočet do sférických souřadnic
Nyní ještě musíme přepočítat souřadnice vektoru rotace do sférických souřadnic. Obecně platí, že pokud má vektor \(\vec u\) souřadnice (ux, uy, uz), potom lze psát
\[\vec u = u_x\vec x_0 + u_y\vec y_0 + u_z\vec z_0\]a souřadnice (ur, uθ, uφ) lze získat tak, že za vektory \(\vec x_0, \vec y_0, \vec z_0\) dosadíme jejich vyjádření pomocí vektorů \(\vec r_0, \vec \theta_0, \vec \varphi_0\) (viz výše). Platí, že
\[\vec x_0 = \sin\,\theta\,\cos\,\varphi\,\vec r_0 + \cos\,\theta\,\cos\,\varphi\,\vec \theta_0 - \sin\,\varphi\,\vec\varphi_0\] \[\vec y_0 = \sin\,\theta\,\sin\,\varphi\,\vec r_0 + \cos\,\theta\,\sin\,\varphi\,\vec \theta_0 + \cos\,\varphi\,\vec\varphi_0\] \[\vec z_0 = \cos\,\theta\,\vec r_0 -\sin\,\theta\,\vec\theta_0.\]Dosazením do vztahu výše dostaneme, že
\[\vec u = (u_x\sin\theta\,\cos\varphi + u_y\,sin\theta\,\sin\varphi + u_z\cos\theta)\vec r_0 + (u_x\cos\theta\,\cos\varphi + u_y\,cos\theta\,\sin\varphi - u_z\sin\theta)\vec \theta_0 +\] \[+ (-u_x\sin\varphi + u_y\,\cos\varphi)\vec \varphi_0,\]odkud ihned vyplývá, že
\[u_r = u_x\sin\theta\,\cos\varphi + u_y\,sin\theta\,\sin\varphi + u_z\cos\theta\] \[u_\theta = u_x\cos\theta\,\cos\varphi + u_y\,cos\theta\,\sin\varphi - u_z\sin\theta\] \[u_\varphi = -u_x\sin\varphi + u_y\,\cos\varphi.\]Aplikací tohoto obecného vztahu na \(\vec u = (\nabla\times\vec v)\) dostaneme, že
\[(\nabla\times\vec v)_r = (\nabla\times\vec v)_x\sin\theta\,\cos\varphi + (\nabla\times\vec v)_y\,sin\theta\,\sin\varphi + (\nabla\times\vec v)_z\cos\theta ,\]což s využitím vztahů (R1), (R2), (R3) výše dává po úpravě
\[(\nabla\times\vec v)_r = \frac{1}{r}\frac{\partial v_\varphi}{\partial \theta} -\frac{1}{r\sin\theta}\frac{\partial v_\theta}{\partial \varphi} + \frac{1}{r\sin\theta}v_\varphi\cos\theta .\]Porovnejme dosažený vztah se vztahem uvedeným v zadání úlohy, z něhož plyne
\[(\nabla\times\vec v)_r = \frac{1}{r\sin\theta}\left[\frac{\partial}{\partial \theta}(\sin\theta\,v_\varphi)-\frac{\partial v_\theta}{\partial\varphi}\right].\]Zjišťujeme, že po provedení derivace v hranaté závorce a roznásobení dostáváme tentýž vztah.
Kýžené rovnosti pro ostatní složky vektoru rotace dostaneme analogicky dosazením (R1), (R2), (R3) do vztahů
\[(\nabla\times\vec v)_\theta = (\nabla\times\vec v)_x\cos\theta\,\cos\varphi + (\nabla\times\vec v)_y\,cos\theta\,\sin\varphi - (\nabla\times\vec v)_z\sin\theta\] \[(\nabla\times\vec v)_\varphi = -(\nabla\times\vec v)_x\sin\varphi + (\nabla\times\vec v)_y\,\cos\varphi.\]Odpověď
Ověřili jsme vztahy v zadání úlohy.