Kegelkoordinaten

aus Wikipedia, der freien Enzyklopädie
Zur Navigation springen Zur Suche springen
Koordinatenflächen der Kegel Koordinaten mit b=1 und c=2. Auf der roten Kugel ist r=2, auf dem blauen Kegel mit senkrechter Achse ist θ=cosh(1) und auf dem gelben Kegel, dessen Achse grün ist, ist λ2 = 2/3. Der schwarze Kreis markiert einen Punkt mit den kartesischen Koordinaten (1.26, -0.78, 1.34).

Kegelkoordinaten sind orthogonale Koordinaten, in denen ein Punkt des dreidimensionalen Raums durch Angabe der Lage auf einer Kugel und zwei elliptischen Kegeln bestimmt wird, siehe Bild.

Kegelkoordinaten (englisch conical coordinates[1] :37[2] :659) erlauben immer eine Trennung der Veränderlichen in der Laplace- und Helmholtz-Gleichung,[1] :8 was deren Lösung vereinfacht. Kegelkoordinaten bieten sich zur Lösung von Randwertaufgaben dort an, wo die Ränder des Gebiets kugel- oder kegelförmig sind. Zur Anpassung an diese Ränder dienen zwei Parameter, im Bild b und c, die die Form der #Koordinatenflächen beeinflussen. Kegelkoordinaten wurden auf mehrere verschiedene Arten definiert.[3]

Sie sind nicht zu verwechseln mit den nicht-orthogonalen Polarkoordinaten gleichen Namens.

Koordinatenflächen

[Bearbeiten | Quelltext bearbeiten ]

Die kartesischen Koordinaten ( x , y , z ) {\displaystyle (x,y,z)} {\displaystyle (x,y,z)} berechnen sich aus den Kegelkoordinaten r , θ , λ R , r , θ > 0 {\displaystyle r,\theta ,\lambda \in \mathbb {R} ,r,\theta >0} {\displaystyle r,\theta ,\lambda \in \mathbb {R} ,r,\theta >0} bei gegebenem b , c R {\displaystyle b,c\in \mathbb {R} } {\displaystyle b,c\in \mathbb {R} } gemäß:[1] :37

r := ( x y z ) = r b c c 2 b 2 ( θ λ c 2 b 2 c ( θ 2 b 2 ) ( b 2 λ 2 ) b ( c 2 θ 2 ) ( c 2 λ 2 ) ) {\displaystyle {\vec {r}}:={\begin{pmatrix}x\\y\\z\end{pmatrix}}={\frac {r}{bc{\sqrt {c^{2}-b^{2}}}}}{\begin{pmatrix}\theta \lambda {\sqrt {c^{2}-b^{2}}}\\c{\sqrt {(\theta ^{2}-b^{2})(b^{2}-\lambda ^{2})}}\\b{\sqrt {(c^{2}-\theta ^{2})(c^{2}-\lambda ^{2})}}\end{pmatrix}}} {\displaystyle {\vec {r}}:={\begin{pmatrix}x\\y\\z\end{pmatrix}}={\frac {r}{bc{\sqrt {c^{2}-b^{2}}}}}{\begin{pmatrix}\theta \lambda {\sqrt {c^{2}-b^{2}}}\\c{\sqrt {(\theta ^{2}-b^{2})(b^{2}-\lambda ^{2})}}\\b{\sqrt {(c^{2}-\theta ^{2})(c^{2}-\lambda ^{2})}}\end{pmatrix}}}

Damit das möglich ist, muss

0 λ 2 < b 2 < θ 2 < c 2 {\displaystyle 0\leq \lambda ^{2}<b^{2}<\theta ^{2}<c^{2}} {\displaystyle 0\leq \lambda ^{2}<b^{2}<\theta ^{2}<c^{2}}

sein. Die Koordinatenflächen bestehen aus einer Kugel (r=const., rot im Bild oben),

x 2 + y 2 + z 2 = r 2 {\displaystyle x^{2}+y^{2}+z^{2}=r^{2}} {\displaystyle x^{2}+y^{2}+z^{2}=r^{2}}

einem elliptischen Kegel um die z-Achse (θ=const., blau im Bild oben),

x 2 θ 2 + y 2 θ 2 b 2 z 2 c 2 θ 2 = 0 {\displaystyle {\frac {x^{2}}{\theta ^{2}}}+{\frac {y^{2}}{\theta ^{2}-b^{2}}}-{\frac {z^{2}}{c^{2}-\theta ^{2}}}=0} {\displaystyle {\frac {x^{2}}{\theta ^{2}}}+{\frac {y^{2}}{\theta ^{2}-b^{2}}}-{\frac {z^{2}}{c^{2}-\theta ^{2}}}=0}

und einem elliptischen Kegel um die x-Achse (λ=const., gelb im Bild oben),

x 2 λ 2 y 2 b 2 λ 2 z 2 c 2 λ 2 = 0 {\displaystyle {\frac {x^{2}}{\lambda ^{2}}}-{\frac {y^{2}}{b^{2}-\lambda ^{2}}}-{\frac {z^{2}}{c^{2}-\lambda ^{2}}}=0} {\displaystyle {\frac {x^{2}}{\lambda ^{2}}}-{\frac {y^{2}}{b^{2}-\lambda ^{2}}}-{\frac {z^{2}}{c^{2}-\lambda ^{2}}}=0}

Aus obigen drei Gleichungen können die Koordinatenquadrate bestimmt werden:

r 2 = x 2 + y 2 + z 2 θ 2 = b 2 ( r 2 y 2 ) + c 2 ( r 2 z 2 ) + [ b 2 ( r 2 y 2 ) + c 2 ( r 2 z 2 ) ] 2 ( 2 b c r x ) 2 2 r 2 λ 2 = b 2 ( r 2 y 2 ) + c 2 ( r 2 z 2 ) [ b 2 ( r 2 y 2 ) + c 2 ( r 2 z 2 ) ] 2 ( 2 b c r x ) 2 2 r 2 {\displaystyle {\begin{aligned}r^{2}=&x^{2}+y^{2}+z^{2}\\\theta ^{2}=&{\frac {b^{2}(r^{2}-y^{2})+c^{2}(r^{2}-z^{2})+{\sqrt {[b^{2}(r^{2}-y^{2})+c^{2}(r^{2}-z^{2})]^{2}-(2bcrx)^{2}}}}{2r^{2}}}\\\lambda ^{2}=&{\frac {b^{2}(r^{2}-y^{2})+c^{2}(r^{2}-z^{2})-{\sqrt {[b^{2}(r^{2}-y^{2})+c^{2}(r^{2}-z^{2})]^{2}-(2bcrx)^{2}}}}{2r^{2}}}\end{aligned}}} {\displaystyle {\begin{aligned}r^{2}=&x^{2}+y^{2}+z^{2}\\\theta ^{2}=&{\frac {b^{2}(r^{2}-y^{2})+c^{2}(r^{2}-z^{2})+{\sqrt {[b^{2}(r^{2}-y^{2})+c^{2}(r^{2}-z^{2})]^{2}-(2bcrx)^{2}}}}{2r^{2}}}\\\lambda ^{2}=&{\frac {b^{2}(r^{2}-y^{2})+c^{2}(r^{2}-z^{2})-{\sqrt {[b^{2}(r^{2}-y^{2})+c^{2}(r^{2}-z^{2})]^{2}-(2bcrx)^{2}}}}{2r^{2}}}\end{aligned}}}

Alternative Formulierungen

[Bearbeiten | Quelltext bearbeiten ]

Eine alternative Formulierung[2] :659 benutzt die Koordinaten

ξ 1 = r , ξ 2 = b 2 λ 2 c = α c n ( μ , α ) , ξ 3 = θ 2 b 2 c = β c n ( ν , β ) x = ξ 1 α ( α 2 ξ 2 2 ) ( α 2 + ξ 3 2 ) = ξ 1 s n ( μ , α ) d n ( ν , β ) y = ξ 1 β ( β 2 + ξ 2 2 ) ( β 2 ξ 3 2 ) = ξ 1 d n ( μ , α ) s n ( ν , β ) z = ξ 1 ξ 2 ξ 3 α β = ξ 1 c n ( μ , α ) c n ( ν , β ) {\displaystyle {\begin{aligned}\xi _{1}=&r,\;\xi _{2}={\frac {\sqrt {b^{2}-\lambda ^{2}}}{c}}=\alpha ,円{\rm {cn}}(\mu ,\alpha ),\;\xi _{3}={\frac {\sqrt {\theta ^{2}-b^{2}}}{c}}=\beta ,円{\rm {cn}}(\nu ,\beta )\\x=&{\frac {\xi _{1}}{\alpha }}{\sqrt {(\alpha ^{2}-\xi _{2}^{2})(\alpha ^{2}+\xi _{3}^{2})}}=\xi _{1},円{\rm {sn}}(\mu ,\alpha ),円{\rm {dn}}(\nu ,\beta )\\y=&{\frac {\xi _{1}}{\beta }}{\sqrt {(\beta ^{2}+\xi _{2}^{2})(\beta ^{2}-\xi _{3}^{2})}}=\xi _{1},円{\rm {dn}}(\mu ,\alpha ),円{\rm {sn}}(\nu ,\beta )\\z=&{\frac {\xi _{1}\xi _{2}\xi _{3}}{\alpha \beta }}=\xi _{1},円{\rm {cn}}(\mu ,\alpha ),円{\rm {cn}}(\nu ,\beta )\end{aligned}}} {\displaystyle {\begin{aligned}\xi _{1}=&r,\;\xi _{2}={\frac {\sqrt {b^{2}-\lambda ^{2}}}{c}}=\alpha ,円{\rm {cn}}(\mu ,\alpha ),\;\xi _{3}={\frac {\sqrt {\theta ^{2}-b^{2}}}{c}}=\beta ,円{\rm {cn}}(\nu ,\beta )\\x=&{\frac {\xi _{1}}{\alpha }}{\sqrt {(\alpha ^{2}-\xi _{2}^{2})(\alpha ^{2}+\xi _{3}^{2})}}=\xi _{1},円{\rm {sn}}(\mu ,\alpha ),円{\rm {dn}}(\nu ,\beta )\\y=&{\frac {\xi _{1}}{\beta }}{\sqrt {(\beta ^{2}+\xi _{2}^{2})(\beta ^{2}-\xi _{3}^{2})}}=\xi _{1},円{\rm {dn}}(\mu ,\alpha ),円{\rm {sn}}(\nu ,\beta )\\z=&{\frac {\xi _{1}\xi _{2}\xi _{3}}{\alpha \beta }}=\xi _{1},円{\rm {cn}}(\mu ,\alpha ),円{\rm {cn}}(\nu ,\beta )\end{aligned}}}

mit den drei grundlegenden Jacobischen Funktionen sinus– sn, cosinus– cn bzw. delta amplitudinis dn, dem elliptischen Modul α = b / c {\displaystyle \alpha ={\sqrt {b/c}}} {\displaystyle \alpha ={\sqrt {b/c}}}, dem komplementären Parameter β = 1 α 2 {\displaystyle \beta ={\sqrt {1-\alpha ^{2}}}} {\displaystyle \beta ={\sqrt {1-\alpha ^{2}}}} und beliebigen Variablen μ,ν∈R.

Die Koordinatenflächen sind hier eine Kugel und elliptische Kegel um die x- und y-Achse:

x 2 + y 2 + z 2 ξ 1 2 = 1 x 2 α 2 ξ 2 2 y 2 β 2 + ξ 2 2 z 2 ξ 2 2 = 0 x 2 α 2 + ξ 3 2 y 2 β 2 ξ 3 2 + z 2 ξ 3 2 = 0 {\displaystyle {\begin{aligned}{\frac {x^{2}+y^{2}+z^{2}}{\xi _{1}^{2}}}=&1\\{\frac {x^{2}}{\alpha ^{2}-\xi _{2}^{2}}}-{\frac {y^{2}}{\beta ^{2}+\xi _{2}^{2}}}-{\frac {z^{2}}{\xi _{2}^{2}}}=&0\\{\frac {x^{2}}{\alpha ^{2}+\xi _{3}^{2}}}-{\frac {y^{2}}{\beta ^{2}-\xi _{3}^{2}}}+{\frac {z^{2}}{\xi _{3}^{2}}}=&0\end{aligned}}} {\displaystyle {\begin{aligned}{\frac {x^{2}+y^{2}+z^{2}}{\xi _{1}^{2}}}=&1\\{\frac {x^{2}}{\alpha ^{2}-\xi _{2}^{2}}}-{\frac {y^{2}}{\beta ^{2}+\xi _{2}^{2}}}-{\frac {z^{2}}{\xi _{2}^{2}}}=&0\\{\frac {x^{2}}{\alpha ^{2}+\xi _{3}^{2}}}-{\frac {y^{2}}{\beta ^{2}-\xi _{3}^{2}}}+{\frac {z^{2}}{\xi _{3}^{2}}}=&0\end{aligned}}}

Hieraus lässt sich für die in diesem Artikel benutzte Formulierung

θ = c d n ( ν , β ) , λ = b s n ( μ , α ) x = r s n ( μ , α ) d n ( ν , β ) , y = r c n ( μ , α ) c n ( ν , β ) , z = r d n ( μ , α ) s n ( ν , β ) {\displaystyle {\begin{aligned}\theta =&c,円{\rm {dn}}(\nu ,\beta ),\;\lambda =b,円{\rm {sn}}(\mu ,\alpha )\\x=&r,円{\rm {sn}}(\mu ,\alpha ),円{\rm {dn}}(\nu ,\beta ),\;y=r,円{\rm {cn}}(\mu ,\alpha ),円{\rm {cn}}(\nu ,\beta ),\;z=r,円{\rm {dn}}(\mu ,\alpha ),円{\rm {sn}}(\nu ,\beta )\end{aligned}}} {\displaystyle {\begin{aligned}\theta =&c,円{\rm {dn}}(\nu ,\beta ),\;\lambda =b,円{\rm {sn}}(\mu ,\alpha )\\x=&r,円{\rm {sn}}(\mu ,\alpha ),円{\rm {dn}}(\nu ,\beta ),\;y=r,円{\rm {cn}}(\mu ,\alpha ),円{\rm {cn}}(\nu ,\beta ),\;z=r,円{\rm {dn}}(\mu ,\alpha ),円{\rm {sn}}(\nu ,\beta )\end{aligned}}}

ableiten.

Metrische Faktoren, Weg- und Volumenelemente

[Bearbeiten | Quelltext bearbeiten ]

Die kovarianten Basisvektoren sind mit r = ( x , y , z ) {\displaystyle {\vec {r}}=(x,y,z)^{\top }} {\displaystyle {\vec {r}}=(x,y,z)^{\top }}:

g r := r r = 1 r ( x y z ) , g θ := r θ = ( x θ y θ θ 2 b 2 z θ c 2 θ 2 ) , g λ := r λ = ( x λ y λ λ 2 b 2 z λ c 2 λ 2 ) {\displaystyle {\vec {g}}_{r}:={\frac {\partial {\vec {r}}}{\partial r}}={\frac {1}{r}}{\begin{pmatrix}x\\y\\z\end{pmatrix}},\;{\vec {g}}_{\theta }:={\frac {\partial {\vec {r}}}{\partial \theta }}={\begin{pmatrix}{\frac {x}{\theta }}\\{\frac {y\theta }{\theta ^{2}-b^{2}}}\\{\frac {-z\theta }{c^{2}-\theta ^{2}}}\end{pmatrix}},\;{\vec {g}}_{\lambda }:={\frac {\partial {\vec {r}}}{\partial \lambda }}={\begin{pmatrix}{\frac {x}{\lambda }}\\{\frac {y\lambda }{\lambda ^{2}-b^{2}}}\\{\frac {-z\lambda }{c^{2}-\lambda ^{2}}}\end{pmatrix}}} {\displaystyle {\vec {g}}_{r}:={\frac {\partial {\vec {r}}}{\partial r}}={\frac {1}{r}}{\begin{pmatrix}x\\y\\z\end{pmatrix}},\;{\vec {g}}_{\theta }:={\frac {\partial {\vec {r}}}{\partial \theta }}={\begin{pmatrix}{\frac {x}{\theta }}\\{\frac {y\theta }{\theta ^{2}-b^{2}}}\\{\frac {-z\theta }{c^{2}-\theta ^{2}}}\end{pmatrix}},\;{\vec {g}}_{\lambda }:={\frac {\partial {\vec {r}}}{\partial \lambda }}={\begin{pmatrix}{\frac {x}{\lambda }}\\{\frac {y\lambda }{\lambda ^{2}-b^{2}}}\\{\frac {-z\lambda }{c^{2}-\lambda ^{2}}}\end{pmatrix}}}

die, wie es sein muss, senkrecht zueinander sind, und deren Beträge die metrischen Faktoren sind:

h r := | g r | = 1 , h θ := | g θ | = r θ 2 λ 2 ( θ 2 b 2 ) ( c 2 θ 2 ) , h λ := | g λ | = r θ 2 λ 2 ( b 2 λ 2 ) ( c 2 λ 2 ) {\displaystyle {\begin{aligned}h_{r}:=|{\vec {g}}_{r}|=1,\quad h_{\theta }:=|{\vec {g}}_{\theta }|=r{\sqrt {\frac {\theta ^{2}-\lambda ^{2}}{(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})}}},\\h_{\lambda }:=|{\vec {g}}_{\lambda }|=r{\sqrt {\frac {\theta ^{2}-\lambda ^{2}}{(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2})}}}\end{aligned}}} {\displaystyle {\begin{aligned}h_{r}:=|{\vec {g}}_{r}|=1,\quad h_{\theta }:=|{\vec {g}}_{\theta }|=r{\sqrt {\frac {\theta ^{2}-\lambda ^{2}}{(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})}}},\\h_{\lambda }:=|{\vec {g}}_{\lambda }|=r{\sqrt {\frac {\theta ^{2}-\lambda ^{2}}{(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2})}}}\end{aligned}}}

Das Orthonormalsystem ist dann

c ^ 1 = 1 b c c 2 b 2 ( θ λ c 2 b 2 c ( θ 2 b 2 ) ( b 2 λ 2 ) b ( c 2 θ 2 ) ( c 2 λ 2 ) ) / ( ) c ^ 2 = 1 b c ( c 2 b 2 ) ( θ 2 λ 2 ) ( λ c 2 b 2 ( θ 2 b 2 ) ( c 2 θ 2 ) θ c ( c 2 θ 2 ) ( b 2 λ 2 ) θ b ( θ 2 b 2 ) ( c 2 λ 2 ) ) c ^ 3 = 1 b c ( c 2 b 2 ) ( θ 2 λ 2 ) ( θ c 2 b 2 ( c 2 λ 2 ) ( b 2 λ 2 ) λ c ( θ 2 b 2 ) ( c 2 λ 2 ) λ b ( c 2 θ 2 ) ( b 2 λ 2 ) ) {\displaystyle {\begin{aligned}{\hat {c}}_{1}=&{\frac {1}{bc{\sqrt {c^{2}-b^{2}}}}}{\begin{pmatrix}\theta \lambda {\sqrt {c^{2}-b^{2}}}\\c{\sqrt {(\theta ^{2}-b^{2})(b^{2}-\lambda ^{2})}}\\b{\sqrt {(c^{2}-\theta ^{2})(c^{2}-\lambda ^{2})}}\end{pmatrix}}/()\\{\hat {c}}_{2}=&{\frac {1}{bc{\sqrt {(c^{2}-b^{2})(\theta ^{2}-\lambda ^{2})}}}}{\begin{pmatrix}\lambda {\sqrt {c^{2}-b^{2}}}{\sqrt {(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})}}\\\theta c{\sqrt {(c^{2}-\theta ^{2})(b^{2}-\lambda ^{2})}}\\-\theta b{\sqrt {(\theta ^{2}-b^{2})(c^{2}-\lambda ^{2})}}\end{pmatrix}}\\{\hat {c}}_{3}=&{\frac {1}{bc{\sqrt {(c^{2}-b^{2})(\theta ^{2}-\lambda ^{2})}}}}{\begin{pmatrix}\theta {\sqrt {c^{2}-b^{2}}}{\sqrt {(c^{2}-\lambda ^{2})(b^{2}-\lambda ^{2})}}\\-\lambda c{\sqrt {(\theta ^{2}-b^{2})(c^{2}-\lambda ^{2})}}\\-\lambda b{\sqrt {(c^{2}-\theta ^{2})(b^{2}-\lambda ^{2})}}\end{pmatrix}}\end{aligned}}} {\displaystyle {\begin{aligned}{\hat {c}}_{1}=&{\frac {1}{bc{\sqrt {c^{2}-b^{2}}}}}{\begin{pmatrix}\theta \lambda {\sqrt {c^{2}-b^{2}}}\\c{\sqrt {(\theta ^{2}-b^{2})(b^{2}-\lambda ^{2})}}\\b{\sqrt {(c^{2}-\theta ^{2})(c^{2}-\lambda ^{2})}}\end{pmatrix}}/()\\{\hat {c}}_{2}=&{\frac {1}{bc{\sqrt {(c^{2}-b^{2})(\theta ^{2}-\lambda ^{2})}}}}{\begin{pmatrix}\lambda {\sqrt {c^{2}-b^{2}}}{\sqrt {(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})}}\\\theta c{\sqrt {(c^{2}-\theta ^{2})(b^{2}-\lambda ^{2})}}\\-\theta b{\sqrt {(\theta ^{2}-b^{2})(c^{2}-\lambda ^{2})}}\end{pmatrix}}\\{\hat {c}}_{3}=&{\frac {1}{bc{\sqrt {(c^{2}-b^{2})(\theta ^{2}-\lambda ^{2})}}}}{\begin{pmatrix}\theta {\sqrt {c^{2}-b^{2}}}{\sqrt {(c^{2}-\lambda ^{2})(b^{2}-\lambda ^{2})}}\\-\lambda c{\sqrt {(\theta ^{2}-b^{2})(c^{2}-\lambda ^{2})}}\\-\lambda b{\sqrt {(c^{2}-\theta ^{2})(b^{2}-\lambda ^{2})}}\end{pmatrix}}\end{aligned}}}

Das Linien- und Volumenelement lauten[1] :38

d r = g r d r + g θ d θ + g λ d λ d s 2 := | d r | 2 = d r 2 + r 2 ( θ 2 λ 2 ) ( θ 2 b 2 ) ( c 2 θ 2 ) d θ 2 + r 2 ( θ 2 λ 2 ) ( c 2 λ 2 ) ( b 2 λ 2 ) d λ 2 d V = r 4 ( θ 2 λ 2 ) ( c 2 b 2 ) b c y z d r d θ d λ {\displaystyle {\begin{aligned}{\rm {d}}{\vec {r}}=&{\vec {g}}_{r},円{\rm {d}}r+{\vec {g}}_{\theta },円{\rm {d}}\theta +{\vec {g}}_{\lambda },円{\rm {d}}\lambda \\{\rm {d}}s^{2}:=&|{\rm {d}}{\vec {r}}|^{2}={\rm {d}}r^{2}+{\frac {r^{2}(\theta ^{2}-\lambda ^{2})}{(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})}}{\rm {d}}\theta ^{2}+{\frac {r^{2}(\theta ^{2}-\lambda ^{2})}{(c^{2}-\lambda ^{2})(b^{2}-\lambda ^{2})}}{\rm {d}}\lambda ^{2}\\{\rm {d}}V=&{\frac {-r^{4}(\theta ^{2}-\lambda ^{2})}{(c^{2}-b^{2})bcyz}}{\rm {d}}r,円{\rm {d}}\theta ,円{\rm {d}}\lambda \end{aligned}}} {\displaystyle {\begin{aligned}{\rm {d}}{\vec {r}}=&{\vec {g}}_{r},円{\rm {d}}r+{\vec {g}}_{\theta },円{\rm {d}}\theta +{\vec {g}}_{\lambda },円{\rm {d}}\lambda \\{\rm {d}}s^{2}:=&|{\rm {d}}{\vec {r}}|^{2}={\rm {d}}r^{2}+{\frac {r^{2}(\theta ^{2}-\lambda ^{2})}{(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})}}{\rm {d}}\theta ^{2}+{\frac {r^{2}(\theta ^{2}-\lambda ^{2})}{(c^{2}-\lambda ^{2})(b^{2}-\lambda ^{2})}}{\rm {d}}\lambda ^{2}\\{\rm {d}}V=&{\frac {-r^{4}(\theta ^{2}-\lambda ^{2})}{(c^{2}-b^{2})bcyz}}{\rm {d}}r,円{\rm {d}}\theta ,円{\rm {d}}\lambda \end{aligned}}}

Die Basisvektoren bilden demnach ein Rechtssystem, wo das Produkt bcyz im Nenner des Volumenelements negativ ist.

Operatoren in Kegelkoordinaten

[Bearbeiten | Quelltext bearbeiten ]

Wegen der länglichen Ausdrücke für die metrischen Faktoren, wird auf die allgemeine Darstellung der Operatoren Gradient, Divergenz und Rotation eines Vektorfeldes im Hauptartikel verwiesen.

Der Laplace-Operator ist:[1] :38

Δ f = 2 f r 2 + 2 r f r + + 1 r 2 ( θ 2 λ 2 ) { ( θ 2 b 2 ) ( c 2 θ 2 ) 2 f θ 2 θ ( 2 θ 2 b 2 c 2 ) f θ + + ( b 2 λ 2 ) ( c 2 λ 2 ) 2 f λ 2 + λ ( 2 λ 2 b 2 c 2 ) f λ } {\displaystyle {\begin{aligned}\Delta f=&{\frac {\partial ^{2}f}{\partial r^{2}}}+{\frac {2}{r}}{\frac {\partial f}{\partial r}}+\dots \\&\dots +{\frac {1}{r^{2}(\theta ^{2}-\lambda ^{2})}}{\Bigg \{}(\theta ^{2}-b^{2})(c^{2}-\theta ^{2}){\frac {\partial ^{2}f}{\partial \theta ^{2}}}-\theta (2\theta ^{2}-b^{2}-c^{2}){\frac {\partial f}{\partial \theta }}+\dots \\&\qquad \qquad \qquad \qquad \;\dots +(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2}){\frac {\partial ^{2}f}{\partial \lambda ^{2}}}+\lambda (2\lambda ^{2}-b^{2}-c^{2}){\frac {\partial f}{\partial \lambda }}{\Bigg \}}\end{aligned}}} {\displaystyle {\begin{aligned}\Delta f=&{\frac {\partial ^{2}f}{\partial r^{2}}}+{\frac {2}{r}}{\frac {\partial f}{\partial r}}+\dots \\&\dots +{\frac {1}{r^{2}(\theta ^{2}-\lambda ^{2})}}{\Bigg \{}(\theta ^{2}-b^{2})(c^{2}-\theta ^{2}){\frac {\partial ^{2}f}{\partial \theta ^{2}}}-\theta (2\theta ^{2}-b^{2}-c^{2}){\frac {\partial f}{\partial \theta }}+\dots \\&\qquad \qquad \qquad \qquad \;\dots +(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2}){\frac {\partial ^{2}f}{\partial \lambda ^{2}}}+\lambda (2\lambda ^{2}-b^{2}-c^{2}){\frac {\partial f}{\partial \lambda }}{\Bigg \}}\end{aligned}}}

Lösung der Laplace- und Helmholtz-Gleichung in Kegelkoordinaten

[Bearbeiten | Quelltext bearbeiten ]

Kegelkoordinaten bieten sich bei der Lösung von Randwertaufgaben an, in denen die Ränder kugel- oder kegelförmig sind.[1] :1 Die Lösung wird erleichtert, wenn eine Trennung der Variablen gelingt, was in Kegelkoordinaten immer möglich ist[1] :7[2] :511 Dazu wird der Separationsansatz [1] :39

ϕ ( r , θ , λ ) = R ( r ) Θ ( θ ) Λ ( λ ) {\displaystyle \phi (r,\theta ,\lambda )=R(r)\cdot \Theta (\theta )\cdot \Lambda (\lambda )} {\displaystyle \phi (r,\theta ,\lambda )=R(r)\cdot \Theta (\theta )\cdot \Lambda (\lambda )}

in die Helmholtz-Gleichung Δ ϕ + κ 2 ϕ = 0 {\displaystyle \Delta \phi +\kappa ^{2}\phi =0} {\displaystyle \Delta \phi +\kappa ^{2}\phi =0} eingesetzt. Die Faktoren bestimmen sich dann aus den drei gewöhnlichen Differenzialgleichungen [1] :39

2 R r 2 + 2 r R r + ( κ 2 α 2 r 2 ) R = 0 ( θ 2 b 2 ) ( c 2 θ 2 ) 2 Θ θ 2 θ ( 2 θ 2 b 2 c 2 ) Θ θ + ( α 2 θ 2 α 3 ) Θ = 0 ( b 2 λ 2 ) ( c 2 λ 2 ) 2 Λ λ 2 + λ ( 2 λ 2 b 2 c 2 ) Λ λ ( α 2 λ 2 α 3 ) Λ = 0 {\displaystyle {\begin{aligned}{\frac {\partial ^{2}R}{\partial r^{2}}}+{\frac {2}{r}}{\frac {\partial R}{\partial r}}+\left(\kappa ^{2}-{\frac {\alpha _{2}}{r^{2}}}\right)R=&0\\(\theta ^{2}-b^{2})(c^{2}-\theta ^{2}){\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}-\theta (2\theta ^{2}-b^{2}-c^{2}){\frac {\partial \Theta }{\partial \theta }}+(\alpha _{2}\theta ^{2}-\alpha _{3})\Theta =&0\\(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2}){\frac {\partial ^{2}\Lambda }{\partial \lambda ^{2}}}+\lambda (2\lambda ^{2}-b^{2}-c^{2}){\frac {\partial \Lambda }{\partial \lambda }}-(\alpha _{2}\lambda ^{2}-\alpha _{3})\Lambda =&0\end{aligned}}} {\displaystyle {\begin{aligned}{\frac {\partial ^{2}R}{\partial r^{2}}}+{\frac {2}{r}}{\frac {\partial R}{\partial r}}+\left(\kappa ^{2}-{\frac {\alpha _{2}}{r^{2}}}\right)R=&0\\(\theta ^{2}-b^{2})(c^{2}-\theta ^{2}){\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}-\theta (2\theta ^{2}-b^{2}-c^{2}){\frac {\partial \Theta }{\partial \theta }}+(\alpha _{2}\theta ^{2}-\alpha _{3})\Theta =&0\\(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2}){\frac {\partial ^{2}\Lambda }{\partial \lambda ^{2}}}+\lambda (2\lambda ^{2}-b^{2}-c^{2}){\frac {\partial \Lambda }{\partial \lambda }}-(\alpha _{2}\lambda ^{2}-\alpha _{3})\Lambda =&0\end{aligned}}}

Denn die Helmholtz-Gleichung lautet mit dem Ansatz:

Δ ϕ = 2 R r 2 Θ Λ + 2 r R r Θ Λ + + 1 r 2 ( θ 2 λ 2 ) { ( θ 2 b 2 ) ( c 2 θ 2 ) R 2 Θ θ 2 Λ θ ( 2 θ 2 b 2 c 2 ) R Θ θ Λ + + ( b 2 λ 2 ) ( c 2 λ 2 ) R Θ 2 Λ λ 2 + λ ( 2 λ 2 b 2 c 2 ) R Θ Λ λ } + + κ 2 R Θ Λ = 0 {\displaystyle {\begin{aligned}\Delta \phi =&{\frac {\partial ^{2}R}{\partial r^{2}}}\Theta \Lambda +{\frac {2}{r}}{\frac {\partial R}{\partial r}}\Theta \Lambda +\dots \\&\dots +{\frac {1}{r^{2}(\theta ^{2}-\lambda ^{2})}}{\Bigg \{}(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})R{\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}\Lambda -\theta (2\theta ^{2}-b^{2}-c^{2})R{\frac {\partial \Theta }{\partial \theta }}\Lambda +\dots \\&\qquad \qquad \qquad \quad \dots +(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2})R\Theta {\frac {\partial ^{2}\Lambda }{\partial \lambda ^{2}}}+\lambda (2\lambda ^{2}-b^{2}-c^{2})R\Theta {\frac {\partial \Lambda }{\partial \lambda }}{\Bigg \}}+\dots \\&\dots +\kappa ^{2}R\Theta \Lambda =0\end{aligned}}} {\displaystyle {\begin{aligned}\Delta \phi =&{\frac {\partial ^{2}R}{\partial r^{2}}}\Theta \Lambda +{\frac {2}{r}}{\frac {\partial R}{\partial r}}\Theta \Lambda +\dots \\&\dots +{\frac {1}{r^{2}(\theta ^{2}-\lambda ^{2})}}{\Bigg \{}(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})R{\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}\Lambda -\theta (2\theta ^{2}-b^{2}-c^{2})R{\frac {\partial \Theta }{\partial \theta }}\Lambda +\dots \\&\qquad \qquad \qquad \quad \dots +(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2})R\Theta {\frac {\partial ^{2}\Lambda }{\partial \lambda ^{2}}}+\lambda (2\lambda ^{2}-b^{2}-c^{2})R\Theta {\frac {\partial \Lambda }{\partial \lambda }}{\Bigg \}}+\dots \\&\dots +\kappa ^{2}R\Theta \Lambda =0\end{aligned}}}

Multiplikation mit r 2 R Θ Λ {\displaystyle {\tfrac {r^{2}}{R\Theta \Lambda }}} {\displaystyle {\tfrac {r^{2}}{R\Theta \Lambda }}} liefert umgestellt:

r 2 R 2 R r 2 + 2 r R R r + κ 2 r 2 + + 1 θ 2 λ 2 { ( θ 2 b 2 ) ( c 2 θ 2 ) 2 Θ θ 2 Θ θ ( 2 θ 2 b 2 c 2 ) Θ θ Θ + + ( b 2 λ 2 ) ( c 2 λ 2 ) 2 Λ λ 2 Λ + λ ( 2 λ 2 b 2 c 2 ) Λ λ Λ } = 0 {\displaystyle {\begin{aligned}&{\frac {r^{2}}{R}}{\frac {\partial ^{2}R}{\partial r^{2}}}+{\frac {2r}{R}}{\frac {\partial R}{\partial r}}+\kappa ^{2}r^{2}+\dots \\&\dots +{\frac {1}{\theta ^{2}-\lambda ^{2}}}{\Bigg \{}(\theta ^{2}-b^{2})(c^{2}-\theta ^{2}){\frac {\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}{\Theta }}-\theta (2\theta ^{2}-b^{2}-c^{2}){\frac {\frac {\partial \Theta }{\partial \theta }}{\Theta }}+\dots \\&\qquad \qquad \qquad \dots +(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2}){\frac {\frac {\partial ^{2}\Lambda }{\partial \lambda ^{2}}}{\Lambda }}+\lambda (2\lambda ^{2}-b^{2}-c^{2}){\frac {\frac {\partial \Lambda }{\partial \lambda }}{\Lambda }}{\Bigg \}}=0\end{aligned}}} {\displaystyle {\begin{aligned}&{\frac {r^{2}}{R}}{\frac {\partial ^{2}R}{\partial r^{2}}}+{\frac {2r}{R}}{\frac {\partial R}{\partial r}}+\kappa ^{2}r^{2}+\dots \\&\dots +{\frac {1}{\theta ^{2}-\lambda ^{2}}}{\Bigg \{}(\theta ^{2}-b^{2})(c^{2}-\theta ^{2}){\frac {\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}{\Theta }}-\theta (2\theta ^{2}-b^{2}-c^{2}){\frac {\frac {\partial \Theta }{\partial \theta }}{\Theta }}+\dots \\&\qquad \qquad \qquad \dots +(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2}){\frac {\frac {\partial ^{2}\Lambda }{\partial \lambda ^{2}}}{\Lambda }}+\lambda (2\lambda ^{2}-b^{2}-c^{2}){\frac {\frac {\partial \Lambda }{\partial \lambda }}{\Lambda }}{\Bigg \}}=0\end{aligned}}}

Weil nur die Terme in der ersten Zeile vom Radius r abhängen, ergänzen sie sich in Summe zu einer Konstanten α2:

r 2 R 2 R r 2 + 2 r R R r + κ 2 r 2 = α 2 {\displaystyle {\frac {r^{2}}{R}}{\frac {\partial ^{2}R}{\partial r^{2}}}+{\frac {2r}{R}}{\frac {\partial R}{\partial r}}+\kappa ^{2}r^{2}=\alpha _{2}} {\displaystyle {\frac {r^{2}}{R}}{\frac {\partial ^{2}R}{\partial r^{2}}}+{\frac {2r}{R}}{\frac {\partial R}{\partial r}}+\kappa ^{2}r^{2}=\alpha _{2}}

Diese eingesetzt erlaubt auch θ und λ zu trennen:

( θ 2 b 2 ) ( c 2 θ 2 ) 2 Θ θ 2 Θ θ ( 2 θ 2 b 2 c 2 ) Θ θ Θ + α 2 θ 2 = = ( b 2 λ 2 ) ( c 2 λ 2 ) 2 Λ λ 2 Λ + λ ( 2 λ 2 b 2 c 2 ) Λ λ Λ α 2 λ 2 {\displaystyle {\begin{aligned}&(\theta ^{2}-b^{2})(c^{2}-\theta ^{2}){\frac {\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}{\Theta }}-\theta (2\theta ^{2}-b^{2}-c^{2}){\frac {\frac {\partial \Theta }{\partial \theta }}{\Theta }}+\alpha _{2}\theta ^{2}=\dots \\&\qquad \qquad \qquad \;\dots =(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2}){\frac {\frac {\partial ^{2}\Lambda }{\partial \lambda ^{2}}}{\Lambda }}+\lambda (2\lambda ^{2}-b^{2}-c^{2}){\frac {\frac {\partial \Lambda }{\partial \lambda }}{\Lambda }}-\alpha _{2}\lambda ^{2}\end{aligned}}} {\displaystyle {\begin{aligned}&(\theta ^{2}-b^{2})(c^{2}-\theta ^{2}){\frac {\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}{\Theta }}-\theta (2\theta ^{2}-b^{2}-c^{2}){\frac {\frac {\partial \Theta }{\partial \theta }}{\Theta }}+\alpha _{2}\theta ^{2}=\dots \\&\qquad \qquad \qquad \;\dots =(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2}){\frac {\frac {\partial ^{2}\Lambda }{\partial \lambda ^{2}}}{\Lambda }}+\lambda (2\lambda ^{2}-b^{2}-c^{2}){\frac {\frac {\partial \Lambda }{\partial \lambda }}{\Lambda }}-\alpha _{2}\lambda ^{2}\end{aligned}}}

Weil die Summe auf der linken Seite nur von θ und die auf der rechten nur von λ abhängt, ergeben beide eine Konstante α3, was auf die drei oben angegebenen gewöhnlichen Differenzialgleichungen zur Bestimmung der Faktoren R,Θ und Λ führt.

Die im Hauptartikel angegebene Methode zur Separation der Helmholtz-Gleichung führt mit der Stäckel-Matrix[1] :37

S = ( 1 1 r 2 0 0 ϑ 2 ( c 2 ϑ 2 ) ( ϑ 2 b 2 ) 1 ( c 2 ϑ 2 ) ( ϑ 2 b 2 ) 0 λ 2 ( b 2 λ 2 ) ( c 2 λ 2 ) 1 ( b 2 λ 2 ) ( c 2 λ 2 ) ) {\displaystyle \mathbf {S} ={\begin{pmatrix}1&-{\frac {1}{r^{2}}}&0\0円&{\frac {\vartheta ^{2}}{(c^{2}-\vartheta ^{2})(\vartheta ^{2}-b^{2})}}&{\frac {-1}{(c^{2}-\vartheta ^{2})(\vartheta ^{2}-b^{2})}}\0円&{\frac {-\lambda ^{2}}{(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2})}}&{\frac {1}{(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2})}}\end{pmatrix}}} {\displaystyle \mathbf {S} ={\begin{pmatrix}1&-{\frac {1}{r^{2}}}&0\0円&{\frac {\vartheta ^{2}}{(c^{2}-\vartheta ^{2})(\vartheta ^{2}-b^{2})}}&{\frac {-1}{(c^{2}-\vartheta ^{2})(\vartheta ^{2}-b^{2})}}\0円&{\frac {-\lambda ^{2}}{(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2})}}&{\frac {1}{(b^{2}-\lambda ^{2})(c^{2}-\lambda ^{2})}}\end{pmatrix}}}

und den Funktionen

f 1 = r 2 , f 2 = ( θ 2 b 2 ) ( c 2 θ 2 ) , f 3 = ( b 2 λ 2 ) ( c 2 λ 2 ) {\displaystyle f_{1}=r^{2},\;f_{2}={\sqrt {(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})}},\;f_{3}={\sqrt {(b^{2}-\lambda ^{2}),円(c^{2}-\lambda ^{2})}}} {\displaystyle f_{1}=r^{2},\;f_{2}={\sqrt {(\theta ^{2}-b^{2})(c^{2}-\theta ^{2})}},\;f_{3}={\sqrt {(b^{2}-\lambda ^{2}),円(c^{2}-\lambda ^{2})}}}

auf ein vergleichbares Ergebnis.

Einzelnachweise

[Bearbeiten | Quelltext bearbeiten ]
  1. a b c d e f g h i j P. Moon, D.E. Spencer: Field Theory Handbook. Including Coordinate Systems, Differential Equations and Their Solutions. 2. Auflage. Springer Verlag, Berlin, Heidelberg, New York 1971, ISBN 3-540-02732-7. 
  2. a b c P. M. Morse, H. Feshbach: Methods of Theoretical Physics, Part I. McGraw-Hill, New York 1953 (archive.org). 
  3. Eric Weisstein: Conical Coordinates. MathWorld, 22. Juni 2024, abgerufen am 23. Juni 2024 (englisch). 
Abgerufen von „https://de.wikipedia.org/w/index.php?title=Kegelkoordinaten&oldid=246236764"