Kugelkoordinaten

In räumlichen Polarkoordinaten wird ein Punkt durch seinen Abstand vom Ursprung und durch zwei Winkel angegeben.

Bei Punkten auf einer Kugeloberfläche (Sphäre) ist der Abstand vom Ursprung (Kugelmittelpunkt) konstant. Die Zahl der Variablen reduziert sich auf die beiden Winkel, die als sphärische Koordinaten bezeichnet werden.

Der Begriff Kugelkoordinaten kann als Oberbegriff für den allgemeinen Fall und die sphärischen Koordinaten angesehen werden.

Für Polarkoordinaten in der Ebene (ein Abstand, ein Winkel) und Zylinderkoordinaten (zwei Abstände, ein Winkel) siehe den Artikel Polarkoordinaten.

Übliche Konvention

Kugelkoordinaten und kartesisches Koordinatensystem

Die nebenstehende Abbildung zeigt einen Punkt $P$ mit den Kugelkoordinaten $(r, \varphi, \theta)$. Die Transformationsgleichungen von kartesischen Koordinaten $(x,y,z)$ in Kugelkoordinaten lauten:

${r}=\sqrt{x^2+y^2+z^2}$

$\varphi = \operatorname{atan2}(y, x) = \begin{cases} \arctan(\frac y x) & \text{, wenn } x > 0, \\ \sgn(y)\frac \pi 2 & \text{, wenn } x = 0, \\ \arctan(\frac y x) + \pi & \text{, wenn } x < 0 \land y \geq 0, \\ \arctan(\frac y x) - \pi & \text{, wenn } x < 0 \land y < 0. \end{cases}$

${\theta}=\arccos\frac z{\sqrt{x^2+y^2+z^2}} \ = \arccos\frac z{r} \ = \ \frac{\pi}{2} - \arctan \frac{z}{\sqrt{x^2+y^2}}$

In der Analysis und ihren Anwendungen werden Kugelkoordinaten meist im Bogenmaß angegeben. Ist $\vec r$ der Ortsvektor von $P$ (also der Vektor, der den Koordinatenursprung $O$ mit $P$ verbindet) und $\vec r_{xy}$ die senkrechte Projektion von $\vec r$ in die xy-Ebene, dann haben die Kugelkoordinaten von $P$ folgende anschauliche Bedeutung:

• $r$ (Radius) ist der Abstand des Punktes $P$ vom Koordinatenursprung $O$, also die Länge des Vektors $\vec r$
• $\phi$ bzw. $\varphi$ (Azimutwinkel) ist der Winkel zwischen der positiven x-Achse und $\vec r_{xy}$, gezählt von $-\pi$ bis $\pi$ (-180° bis 180°) gegen den Uhrzeigersinn
• $\theta$ bzw. $\vartheta$ (Polarwinkel) ist der Winkel zwischen der positiven z-Achse und $\vec r$, gezählt von $0$ bis $\pi$ (0° bis 180°)

Die Rücktransformation erfolgt nach den Gleichungen (Parametrisierung):

$x = r \cdot \sin \theta \cdot \cos \varphi$

$y = r \cdot \sin \theta \cdot \sin \varphi$

$z = r \cdot \cos \theta$

Die Einheitsvektoren der sphärischen Koordinaten lauten damit:

$\mathbf{e}_r=\begin{pmatrix} \sin\theta \cdot \cos\varphi\\ \sin\theta \cdot \sin\varphi\\ \cos\theta \end{pmatrix} \qquad \mathbf{e}_\theta= \begin{pmatrix} \cos\theta \cdot \cos\varphi\\ \cos\theta \cdot \sin\varphi\\ -\sin\theta \end{pmatrix} \qquad \mathbf{e}_\varphi = \begin{pmatrix} -\sin\varphi\\ {\color{White}+}\cos\varphi\\ 0 \end{pmatrix}$

Um diese Einheitsvektoren zu erhalten, muss die Parametrisierung nach der jeweiligen Koordinate (hier $r, \theta, \varphi$) abgeleitet und auf $1$ normiert werden. Sie bilden in der Reihenfolge $\mathbf{e}_r, \mathbf{e}_\theta, \mathbf{e}_\varphi$ ein Rechtssystem.

Andere Konventionen

Die obige Koordinatenwahl ist internationaler Konsens in der theoretischen Physik. Manchmal werden die Zeichen θ und φ gerade im umgekehrten Sinne verwendet, insbesondere in amerikanischer Literatur.

Der Polarwinkel θ ist nicht die geographische Breite. Diese ist vielmehr als Winkel zwischen der Äquatorialebene und dem Ortsvektor definiert und nimmt Werte zwischen −90° und 90° an. Wird sie mit φ bezeichnet, so ist φ = 90° − θ, θ = 90° − φ. Hingegen kann man das oben benutzte φ ohne weiteres mit der geographischen Länge λ östlich von Greenwich gleichsetzen. Siehe dazu den Artikel: geographische Koordinaten.

Des Weiteren ist die obige Konstruktion in gewisser Hinsicht inkonsistent zum Aufbau der ebenen Polarkoordinaten. Für manche Probleme ist es praktischer, die Darstellung

$x = r \cos \theta \, \cos \varphi$

$y = r \cos \theta \, \sin \varphi$

$z = r \sin \theta \quad$

zu benutzen. In dieser Darstellung entspricht $\theta$ der geographischen Breite.

Die Rücktransformation des Punktes bzw. Vektors $\vec{p}$ in die Winkelbestandteile erfolgt dann mit:

$\theta = \arcsin(z / r)$

$\varphi = \mathrm{atan2}(y, x)$

wobei $r = |\vec{p}\,|$.

Verallgemeinerung auf n-dimensionale Kugelkoordinaten

Eine Verallgemeinerung der Kugelkoordinaten auf $n$ Dimensionen:

$\begin{align} x_1 & =r\cos(\phi_1)\\ x_2 & =r\sin(\phi_1)\cos(\phi_2)\\ x_3 & =r\sin(\phi_1)\sin(\phi_2)\cos(\phi_3)\\ & {}\,\,\, \vdots\\ x_{n-1} & =r\sin(\phi_1)\cdots\sin(\phi_{n-2})\cos(\phi_{n-1})\\ x_n & =r\sin(\phi_1)\cdots\sin(\phi_{n-2})\sin(\phi_{n-1}) \end{align}$

Die Winkel entwickeln sich nach:

$\begin{align} \tan(\phi_{n-1}) & =\frac{x_{n}}{x_{n-1}} \\ \tan(\phi_{n-2})& =\frac{\sqrt{{x_n}^2+{x_{n-1}}^2}}{x_{n-2}} \\ & {}\,\,\,\vdots\\ \tan(\phi_1) & =\frac{\sqrt{{x_n}^2+{x_{n-1}}^2+\cdots+{x_2}^2}}{x_{1}} \end{align}$

Durch Umnummerierung erhält man eine Rekursionsformel für die Winkel:

$\begin{align} x_n & =r\cos(\phi_{n-1})\\ x_{n-1} & =r\sin(\phi_{n-1})\cos(\phi_{n-2})\\ x_{n-2} & =r\sin(\phi_{n-1})\sin(\phi_{n-2})\cos(\phi_{n-3})\\ & {}\,\,\, \vdots\\ x_{2} & =r\sin(\phi_{n-1})\cdots\sin(\phi_{2})\cos(\phi_{1})\\ x_1 & =r\sin(\phi_{n-1})\cdots\sin(\phi_{2})\sin(\phi_{1}) \end{align}$

Woraus sich die folgenden Winkel ergeben:

$||\vec L_k||=\sgn(x_k)\sqrt{x_k^2+||\vec L_{k-1}||^2}=\frac{x_k}{||x_k||}\sqrt{x_k^2+||\vec L_{k-1}||^2}$

mit $||\vec L_0||=0$ und

$\tan(\phi_k)=\frac{\sqrt{x_k^2+||\vec L_{k-1}||^2}}{x_{k+1}}=\frac{||\vec L_k||}{x_{k+1}}$

Der Radius ist:

$r=||\vec L_n||$

Eine Fallunterscheidung liefert mittels Arkustangens den passenden Winkel zur gegebenen kartesischen Koordinate, wobei $\arctan(\pm \, \infty) = \pm \, \tfrac{\pi}{2}$:

$\begin{align} \phi_k=\begin{cases} \arctan\left(\frac{||\vec L_k||}{x_{k+1}}\right)+\pi, & \text{(1) wenn: }x_{k+1}<0 \; \land \; k = n-1 \\ \arctan\left(\frac{||\vec L_k||}{x_{k+1}}\right), & \text{(2) wenn: } \text{nicht (1)} \land \; \text{nicht (3)} \\ 0 , & \text{(3) wenn: }x_{k+1}=||\vec L_k||=0\\ \end{cases} \end{align}$

Dabei fällt auf, dass $\begin{align} \vec L_k \end{align}$ immer ein zweidimensionaler Vektor ist für $\begin{align} k>0 \end{align}$ .

Beispiel

Zuordnung am Beispiel $n=3$ mit den geläufigen Koordinatenachsen $x,y, z$:

$\begin{align} x_3&=z= r\cos(\phi_{2})\\ x_2&=x= r\sin(\phi_{2})\cos(\phi_{1})\\ x_1&=y= r\sin(\phi_{2})\sin(\phi_{1})\\ \end{align}$

Die Winkel sind dann:

$\begin{align} \tan(\phi_{2})=\frac{||\vec L_2||}{x_{3}}&= \frac{\sqrt{x_2^2 + x_1^2} }{x_3} =\frac{\sqrt{x^2 + y^2} }{z} \\ \tan(\phi_{1})=\frac{||\vec L_1||}{x_{2}}&= \frac{\sqrt{x_1^2} }{x_2} =\frac{y }{x} \end{align}$

Transformation von Differentialen

Jacobi-Matrix

Die lokalen Eigenschaften der Koordinatentransformation werden durch die Jacobi-Matrix beschrieben. Für die Transformation von Kugelkoordinaten in kartesische Koordinaten lautet diese

$J =\frac{\partial(x,y,z)}{\partial(r,\theta,\varphi)} =\begin{pmatrix} \sin\theta\cos\varphi&r\cos\theta\cos\varphi&-r\sin\theta\sin\varphi\\ \sin\theta \sin\varphi&r\cos\theta\sin\varphi&r\sin\theta\cos\varphi\\ \cos\theta&-r\sin\theta&0 \end{pmatrix};$

bei sphärischen Polarkoordinaten (nur θ, φ) fällt die erste Spalte (und Zeile) weg.

Die zugehörige Funktionaldeterminante lautet:

$|J|=r^2\sin\theta$

Die Jacobi-Matrix der entgegengesetzten Transformation ist nur für räumliche, nicht für sphärische Polarkoordinaten definiert; man berechnet sie am einfachsten als Inverse von J:

$J^{-1} =\frac{\partial(r,\theta,\varphi)}{\partial(x,y,z)} =\begin{pmatrix} \sin\theta\cos\varphi & \sin\theta\sin\varphi & \cos\theta\\ \frac{1}{r}\cos\theta\cos\varphi & \frac{1}{r}\cos\theta\sin\varphi & -\frac{1}{r}\sin\theta \\ -\frac{1}{r}\frac{\sin\varphi}{\sin\theta} & \frac{1}{r}\frac{\cos\varphi}{\sin\theta} & 0 \end{pmatrix}.$

Einige Komponenten dieser Matrix sind Brüche, an deren Nennern man die Uneindeutigkeit der Polarkoordinaten bei r=0 und bei sin θ=0 (also θ=0 oder π) erkennt. Ungebräuchlicher ist die Darstellung in kartesischen Koordinaten:

$J^{-1} =\begin{pmatrix} \frac{x}{r}&\frac{y}{r}&\frac{z}{r}\\\\ \frac{xz}{r^2\sqrt{x^2+y^2}}&\frac{yz}{r^2\sqrt{x^2+y^2}}&\frac{-(x^2+y^2)}{r^2\sqrt{x^2+y^2}}\\\\ \frac{-y}{x^2+y^2}&\frac{x}{x^2+y^2}&0 \end{pmatrix}.$

Differentiale, Volumenelement, Linienelement

Die Jacobi-Matrix erlaubt es, die Umrechnung von Differentialen übersichtlich als lineare Abbildung zu schreiben:

$(\mathrm{d}x,\mathrm{d}y,\mathrm{d}z)^T=J\cdot(\mathrm{d}r,\mathrm{d}\theta,\mathrm{d}\varphi)^T$

beziehungsweise

$(\mathrm{d}r,\mathrm{d}\theta,\mathrm{d}\varphi)^T=J^{-1}\cdot(\mathrm{d}x,\mathrm{d}y,\mathrm{d}z)^T$.

Das bedeutet eine Koordinatentransformation $\,K\, \rightarrow\, K'$ gegeben durch die Jacobi-Matrix, und zwar hier von kartesischer zu sphärischer Geometrie. Dabei sind die Komponenten $\textstyle \frac{\partial}{\partial x_i}$ in $J$ Elemente des zu $dx_j$ dualen (Tangential-)Raums (d. h. im Prinzip $dx^i|_p \left(\left. \tfrac{\partial}{\partial x^j}\right|_p \right)=\delta^i\,_j$). Dieser Koordinatenwechsel spielt auch eine Rolle bei der Parametrisierung einer Kurve für Kurvenintegrale.

Das Volumenelement $\mathrm{d}V=\mathrm{d}x \mathrm{d}y \mathrm{d}z$ lässt sich besonders einfach mit Hilfe der Funktionaldeterminante

$|J|=r^2\sin\theta$

umrechnen:

$\, \mathrm{d}V=r^2 \sin \theta \, \mathrm{d}\varphi \, \mathrm{d}\theta \, \mathrm{d}r$.

Hier gilt $|J|=\frac{\sqrt{|g|}}{\sqrt{|g'|}}=\sqrt{|g|}$ (mit $\, |g|$ als Determinante des metrischen Tensors) sodass man gleich sieht, dass:

$\mathrm{d}V=\sqrt{|g|} \, \mathrm{d}\varphi \, \mathrm{d}\theta \, \mathrm{d}r$.

Das spielt eine Rolle im Falle einer Verallgemeinerung im Sinne allgemeiner Skalarprodukte und z. B. des Laplace-Beltrami-Operators (g aus dem allgemeinen Skalarproduktes beinhaltend).

Durch Differentiation dV/dr erhält man für ein Flächenelement dA auf einer Sphäre mit Radius r

$\mathrm{d}A=r^2 \sin \theta \, \mathrm{d}\varphi \, \mathrm{d}\theta$.

Ein Linienelement ds rechnet man gemäß

$\mathrm{d}s^2=\mathrm{d}x^2 + \mathrm{d}y^2 + \mathrm{d}z^2 =\mathrm{d}r^2 + r^2\mathrm{d}\theta^2 + r^2\sin^2\theta\mathrm{d}\varphi^2=\, \Sigma_i\, dx_i dx^i=\Sigma_{ij}\, dx_i g^{ij} dx_j$

um (mit i,j=1 bis 3). Dabei sind $\, g_{ij}$ die Komponenten des metrischen Tensors, der die Multiplikation in nichttrivialen Räumen definiert (kartesisch gilt $\, g_{ij}=\delta_{ij}$), beziehungsweise benutzt wird, um Indizes zu heben und zu senken, d. h. kovariante Tensoren in Kontravarianten und umgekehrt zu transformieren, sodass $\, g$, bzw. $\, ds^2$ eine besonders große Rolle in der allgemeinen und in der speziellen Relativitätstheorie, d.h im Standardmodell der Kosmologie und der Elementarteilchen, spielten.

Metrik und Rotationsmatrix

Im Fehlen gemischter Glieder im Linienelement ds spiegelt sich wider, dass der metrische Tensor

$g=J^T J=\begin{pmatrix}1&0&0\\ 0&r^2&0\\ 0&0&r^2\sin^2\theta\end{pmatrix}$

auch in Kugelkoordinaten keine Außerdiagonalelemente hat.

Der metrische Tensor ist offensichtlich das Quadrat der Diagonalmatrix

$h=\operatorname{diag}(1,r,r\sin\theta)$.

Mit Hilfe dieser Matrix lässt sich die Jacobi-Matrix als J=Sh schreiben, wobei S die Rotationsmatrix

$S =\begin{pmatrix} \sin\theta\cos\varphi&\cos\theta\cos\varphi&-\sin\varphi\\ \sin\theta\sin\varphi&\cos\theta\sin\varphi&\cos\varphi\\ \cos\theta&-\sin\theta&0 \end{pmatrix}$

ist.

Transformation von Vektorfeldern und -Operatoren

Im Folgenden soll die Transformation von Vektoren und Differentialoperatoren exemplarisch hergeleitet werden. Die Ergebnisse werden bevorzugt in kompakter Form unter Benutzung von Transformationsmatrizen geschrieben.

Transformation der Vektorraumbasis

Der Basisvektor e_φ zur Koordinate φ gibt an, in welche Richtung sich ein Punkt P(r, θ, φ) (oft auch Event genannt) bewegt, wenn die Koordinate φ um einen infinitesimalen Betrag dφ verändert wird:

$\mathbf{e}_\varphi \sim \frac{\partial \mathrm{P}}{\partial\varphi}$.

Daraus erhält man

$\mathbf{e}_\varphi \sim \frac{\partial \mathrm{P}}{\partial\varphi} = \frac{\partial x}{\partial \varphi}\frac{\partial \mathrm{P}}{\partial x} +\frac{\partial y}{\partial \varphi}\frac{\partial \mathrm{P}}{\partial y} +\frac{\partial z}{\partial \varphi}\frac{\partial \mathrm{P}}{\partial z} =-r\sin\theta\sin\varphi\mathbf{e}_x +r\sin\theta\cos\varphi\mathbf{e}_y$.

Um eine orthonormale Basis zu erhalten, muss e_φ noch auf die Länge 1 normiert werden:

$\mathbf{e}_\varphi = -\sin\varphi\mathbf{e}_x + \cos\varphi\mathbf{e}_y$.

In ähnlicher Weise erhält man die Basisvektoren e_r und e_θ. Um die folgenden Transformationen in kompakter Form zu schreiben, verwenden wir die oben eingeführte Rotationsmatrix S. Diese Matrix ist orthogonal, das heißt, $S^{-1}=S^T$. Die normierten Basisvektoren des Kugelkoordinatensystems kann man dann zusammengefasst so mitteilen:

$(\mathbf{e}_r,\mathbf{e}_\theta,\mathbf{e}_\varphi)^T =S^T \cdot (\mathbf{e}_x,\mathbf{e}_y,\mathbf{e}_z)^T$.

Entsprechend lautet die Transformation in die Gegenrichtung

$(\mathbf{e}_x,\mathbf{e}_y,\mathbf{e}_z)^T =S\cdot (\mathbf{e}_r,\mathbf{e}_\theta,\mathbf{e}_\varphi)^T$.

Transformation eines Vektorfeldes

Ein Vektor, als ein geometrisches Objekt, muss vom Koordinatensystem unabhängig sein:

$A_x\mathbf{e}_x + A_y\mathbf{e}_y + A_z\mathbf{e}_z = \mathbf{A} =A_r\mathbf{e}_r + A_\theta\mathbf{e}_\theta + A_\varphi\mathbf{e}_\varphi.$

Diese Bedingung wird erfüllt durch

$(A_r,A_\theta,A_\varphi)^T =S^T \cdot (A_x,A_y,A_z)^T$

beziehungsweise

$(A_x,A_y,A_z)^T =S \cdot (A_r,A_\theta,A_\varphi)^T$.

Transformation der partiellen Ableitungen

Die partiellen Ableitungen transformieren sich wie die Basisvektoren, aber ohne Normierung. Man kann genau wie oben rechnen, nur lässt man den Punkt P im Zähler weg (tatsächlich werden in der modernen Formulierung der Differentialgeometrie die Koordinatenbasisvektoren des Tangentialraums und die partiellen Ableitungen gleichgesetzt) und verwendet die Jacobi-Matrix J=Sh anstelle der Rotationsmatrix S. Die Transformation lautet also:

$\left(\frac{\partial}{\partial r},\frac{\partial}{\partial\theta},\frac{\partial}{\partial\varphi} \right) = \left(\frac{\partial}{\partial x},\frac{\partial}{\partial y},\frac{\partial}{\partial z} \right) \cdot J$,

und in die Gegenrichtung

$\left(\frac{\partial}{\partial x},\frac{\partial}{\partial y},\frac{\partial}{\partial z} \right) = \left(\frac{\partial}{\partial r},\frac{\partial}{\partial\theta},\frac{\partial}{\partial\varphi} \right) \cdot J^{-1}$.

Transformation des Nabla-Operators

Der Nabla-Operator $\nabla$ hat nur in kartesischen Koordinaten die einfache Form

$\mathbf{\nabla} = \mathbf{e}_x\frac{\partial}{\partial x} +\mathbf{e}_y\frac{\partial}{\partial y} +\mathbf{e}_z\frac{\partial}{\partial z}$.

Sowohl die partiellen Ableitungen als auch die Einheitsvektoren muss man in der oben hergeleiteten Weise transformieren. Man findet:

$\mathbf{\nabla} =\mathbf{e}_r\frac{\partial}{\partial r} + \mathbf{e}_\theta\frac{1}{r}\frac{\partial}{\partial\theta} + \mathbf{e}_\varphi\frac{1}{r\sin\theta}\frac{\partial}{\partial\varphi}$.

In dieser Form kann der transformierte Nabla-Operator unmittelbar angewandt werden, um den Gradienten eines in Kugelkoordinaten gegebenen Skalarfeldes zu berechnen.

Um die Divergenz eines in Kugelkoordinaten gegebenen Vektorfeldes A zu berechnen, ist hingegen zu berücksichtigen, dass $\nabla$ nicht nur auf die Koeffizienten A_r, … wirkt, sondern auch auf die in A implizit enthaltenen Basisvektoren e_r

$\mathbf{\nabla}\cdot\mathbf{A} = \frac{1}{r^2}\frac{\partial}{\partial r}(r^2 A_r) + \frac{1}{r\sin\theta}\frac{\partial}{\partial\theta}(\sin\theta A_\theta) +\frac{1}{r\sin\theta}\frac{\partial}{\partial\varphi}A_\varphi.$

Um die Rotation eines in Kugelkoordinaten gegebenen Vektorfeldes A zu berechnen, ist selbiges zu berücksichtigen:

$\mathbf{\nabla} \times \mathbf{A}= {1 \over r\sin\theta}\left({\partial \over \partial \theta} ( A_\varphi\sin\theta ) - {\partial A_\theta \over \partial \varphi}\right) \mathbf{e}_r + {1 \over r}\left({1 \over \sin\theta}{\partial A_r \over \partial \varphi} - {\partial \over \partial r} ( r A_\varphi ) \right) \mathbf{e}_\theta + {1 \over r}\left({\partial \over \partial r} ( r A_\theta ) - {\partial A_r \over \partial \theta}\right) \mathbf{e}_\varphi$

Transformation des Laplace-Operators

Wenn man in der Divergenzformel als Vektorfeld A den Gradientenoperator $\nabla$ einsetzt, findet man den Laplace-Operator

$\mathbf{\Delta}=\mathbf{\nabla}^2 = \frac{1}{r^2}\frac{\partial}{\partial r} \left( r^2 \frac{\partial}{\partial r} \right) + \frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta} \left( \sin\theta \frac{\partial}{\partial\theta} \right) +\frac{1}{r^2\sin^2\theta}\frac{\partial^2}{\partial\varphi^2}$.

bzw.

$\mathbf{\Delta} = \frac{\partial^2}{\partial r^2} + \frac{2}{r}\frac{\partial}{\partial r} + \frac{1}{r^2}\frac{\partial^2}{\partial \theta^2} + \frac{1}{r^2}\frac{\cos \theta}{\sin \theta}\frac{\partial}{\partial \theta} +\frac{1}{r^2\sin^2\theta}\frac{\partial^2}{\partial\varphi^2}$.

Eine besonders einfache Form um diese Form des Laplace-Operators herzuleiten ergibt sich mit Hilfe der Formel für den Laplace-Beltrami-Operator. Dabei gilt:

$\operatorname{div}\, X=\sum_i \frac{1}{\sqrt{|g|}}\partial _i \sqrt{|g|}X^i= \sum_{ij}\frac{1}{\sqrt{|g|}}\partial _i \sqrt{|g|} g^{ij} X_j$

$(\operatorname{grad}\, f)^i=\partial^i f=\sum_j g^{ij}\partial_j f$.

Dabei gilt $\, g_{ij} g^{jk}=\delta_i^k$, also $\, g^{ij}=(g_{ij})^{-1}$.

Der Laplace-Beltrami-Operator ist dann genauso definiert wie der gewöhnliche Laplace-Operator, gilt aber für allgemeinere riemannsche Mannigfaltigkeiten:

$\Delta\, f= \operatorname{div}\,\operatorname{grad}\, f= \sum_i\frac{1}{\sqrt{|g|}}\partial _i \sqrt{|g|}\partial^i f= \sum_{ij} \frac{1}{\sqrt{|g|}}\partial_i \sqrt{|g|} g^{ij} \partial_j f$.

Der metrische Tensor in Kugelkoordinaten ist diagonal und von der Form

$g= \operatorname{diag}(1,r^2,(r \sin\theta)^2)$.

Somit entfällt die Summe über $i$ und die $g^{ij}$ sind die Kehrwerte der Komponenten des metrischen Tensors, also $1,\ r^{-2}$ und $(r\, \sin\, \theta)^{-2}.$

Die Wurzel des Betrags der Determinante ist $\sqrt{|g|}=r^2 \sin\theta$, also gleich der Jacobi-Determinante.

Setzt man die verschiedenen Größen in die Formel des Laplace-Beltrami-Operators ein, dann ergibt sich

$(\Delta f)_{(1)} = \frac{1}{r^2 \sin\theta}\,\partial_r \left(r^2\sin\theta\,\partial _r f\right)=\frac{1}{r^2}\frac{\partial}{\partial r} r^2 \frac{\partial}{\partial r} f$

$(\Delta f)_{(2)} = \frac{1}{r^2 \sin\theta}\,\partial_\theta\, \left(r^2\sin\theta\,\partial _\theta (r^{-2}f)\right) = \frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\sin\theta \frac{\partial}{\partial\theta}f$

$(\Delta f)_{(3)} = \frac{1}{r^2 \sin\theta}\,\partial_\varphi\, \left(r^2\sin\theta\,\partial _\varphi \left((r^2\sin^2\theta)^{-1} f\right)\right)= \frac{1}{r^2\sin^2\theta}\frac{\partial^2}{\partial\varphi^2}f$.

Addiert man zu $\Delta f= (\Delta f)_{(1)}+ (\Delta f)_{(2)}+ (\Delta f)_{(3)}$, so ergibt sich die obige Formel für den Laplace-Operator in Kugelkoordinaten.

Weblinks

• Matroids Matheplanet: Einführung in die Vektoranalysis (als PDF) von Eckard Specht