Üblicherweise wählt man die zwei Brennpunkte an den Stellen
−
c
{\displaystyle -c}
und
+
c
{\displaystyle +c}
auf der
x
{\displaystyle x}
-Achse eines kartesischen Koordinatensystems . Der Punkt mit den elliptischen Koordinaten
u
,
ψ
∈
R
≥
0
,
ψ
<
2
π
{\displaystyle u,\psi \in \mathbb {R} ^{\geq 0},\psi <2\pi }
[2] :17 hat dann die kartesischen Koordinaten
(
x
y
)
=
(
c
⋅
cosh
u
⋅
cos
ψ
c
⋅
sinh
u
⋅
sin
ψ
)
,
(
u
ψ
)
=
(
a
c
o
s
h
(
w
2
c
)
a
t
a
n
2
(
y
w
2
,
x
w
1
)
)
{\displaystyle {\begin{pmatrix}x\\y\end{pmatrix}}={\begin{pmatrix}c\cdot \cosh u\cdot \cos \psi \\c\cdot \sinh u\cdot \sin \psi \end{pmatrix}},\quad {\begin{pmatrix}u\\\psi \end{pmatrix}}={\begin{pmatrix}{\rm {acosh}}\left({\frac {w_{2}}{\sqrt {c}}}\right)\\{\rm {atan2}}(y\,w_{2},x\,w_{1})\end{pmatrix}}}
mit
• sin, cos:
Sinus und Cosinus
• sinh, cosh:
Sinus hyperbolicus und Kosinus hyperbolicus
• acosh:
Areakosinus hyperbolicus ,
• atan2 :
eine Umkehrfunktion des Tangens
•
w
1
{\displaystyle w_{1}}
:
=
m
1
+
m
1
2
+
y
2
,
m
1
:=
x
2
+
y
2
−
c
2
2
c
{\displaystyle ={\sqrt {m_{1}+{\sqrt {m_{1}^{2}+y^{2}}}}},\;m_{1}:={\frac {x^{2}+y^{2}-c^{2}}{2c}}}
, und
•
w
2
{\displaystyle w_{2}}
:
=
m
2
+
m
2
2
−
x
2
,
m
2
:=
x
2
+
y
2
+
c
2
2
c
{\displaystyle ={\sqrt {m_{2}+{\sqrt {m_{2}^{2}-x^{2}}}}},\;m_{2}:={\frac {x^{2}+y^{2}+c^{2}}{2c}}}
Fasst man die Ebene als komplexe Ebene auf mit imaginärer Einheit i2 =-1, so gilt
x
+
i
y
=
c
⋅
cosh
(
u
+
i
ψ
)
.
{\displaystyle x+{\rm {i}}y=c\cdot \cosh(u+{\rm {i}}\psi ).}
Der Kosinus hyperbolicus ist eine Holomorphe Funktion , was die Orthogonalität der elliptischen Koordinaten in der Ebene begründet.
Die Kurven in der xy-Ebene , auf denen u konstant ist (was die Niveaulinien von u in der xy-Ebene sind,) bilden die Ellipsen [2] :17
x
2
(
c
cosh
u
)
2
+
y
2
(
c
sinh
u
)
2
=
1
{\displaystyle {\frac {x^{2}}{(c\cosh u)^{2}}}+{\frac {y^{2}}{(c\sinh u)^{2}}}=1}
Die Niveaulinien von ψ sind die konfokalen Hyperbeln
x
2
(
c
cos
ψ
)
2
−
y
2
(
c
sin
ψ
)
2
=
1
{\displaystyle {\frac {x^{2}}{(c\cos \psi )^{2}}}-{\frac {y^{2}}{(c\sin \psi )^{2}}}=1}
die nur in Vielfachen von π ⁄2 bzw. 90°, wie in Polarkoordinaten radiale Geraden sind: Für
u
=
0
{\displaystyle u=0}
ist die ψ-Koordinatenlinie zur Verbindungsstrecke der beiden Brennpunkte entartet. Für
ψ
=
0
{\displaystyle \psi =0}
ist die
u
{\displaystyle u}
-Koordinatenlinie zur Halbgeraden
[
c
,
∞
[
{\displaystyle \left[c,\infty \right[}
auf der
x
{\displaystyle x}
-Achse entartet, für
ψ
=
π
{\displaystyle \psi =\pi }
zur dazu spiegelsymmetrischen Halbgerade auf der negativen
x
{\displaystyle x}
-Achse. Für
ψ
=
π
2
{\displaystyle \psi ={\tfrac {\pi }{2}}}
und
ψ
=
3
π
2
{\displaystyle \psi ={\tfrac {3\pi }{2}}}
ist die
u
{\displaystyle u}
-Koordinatenlinie die positive bzw. die negative
y
{\displaystyle y}
-Achse.
Alle Ellipsen und Hyperbeln haben die gleiche lineare Exzentrizität
e
=
c
{\displaystyle e=c}
. Die Ellipsen, auf denen
u
{\displaystyle u}
konstant ist, haben die große Halbachse
a
=
c
cosh
u
{\displaystyle a=c\cosh u}
, die kleine Halbachse
b
=
c
sinh
u
{\displaystyle b=c\sinh u}
und numerische Exzentrizität
ε
=
1
cosh
u
{\displaystyle \varepsilon ={\tfrac {1}{\cosh u}}}
. Die Hyperbeln, auf denen
ψ
{\displaystyle \psi }
konstant ist, haben die waagerechte Halbachse
a
=
c
cos
ψ
{\displaystyle a=c\cos \psi }
, die senkrechte Halbachse
b
=
c
sin
ψ
{\displaystyle b=c\sin \psi }
und numerische Exzentrizität
ε
=
1
cos
ψ
{\displaystyle \varepsilon ={\tfrac {1}{\cos \psi }}}
.
Die Darstellung in dieser Koordinatenform ist nur möglich, weil Kosinus hyperbolicus und Sinus hyperbolicus bzw. Kosinus und Sinus die Beziehungen zwischen großer und kleiner Halbachse (
a
2
=
e
2
+
b
2
{\displaystyle a^{2}=e^{2}+b^{2}}
) bei Ellipsen bzw. reeller und imaginärer Halbachse bei Hyperbeln (
a
2
=
e
2
−
b
2
{\displaystyle a^{2}=e^{2}-b^{2}}
) trivial erfüllen.
Metrische Faktoren, Weg- und Flächenelemente in der Ebene
Bearbeiten
Die kovarianten Basisvektoren sind
g
→
u
=
∂
∂
u
(
x
y
)
=
c
(
sinh
(
u
)
cos
(
ψ
)
cosh
(
u
)
sin
(
ψ
)
)
,
g
→
ψ
=
∂
∂
ψ
(
x
y
)
=
c
(
−
cosh
(
u
)
sin
(
ψ
)
sinh
(
u
)
cos
(
ψ
)
)
{\displaystyle {\vec {g}}_{u}={\frac {\partial }{\partial u}}{\begin{pmatrix}x\\y\end{pmatrix}}=c{\begin{pmatrix}\sinh(u)\cos(\psi )\\\cosh(u)\sin(\psi )\end{pmatrix}},\quad {\vec {g}}_{\psi }={\frac {\partial }{\partial \psi }}{\begin{pmatrix}x\\y\end{pmatrix}}=c{\begin{pmatrix}-\cosh(u)\sin(\psi )\\\sinh(u)\cos(\psi )\end{pmatrix}}}
die, wie es sein muss, senkrecht zueinander sind, und deren Beträge die metrischen Faktoren sind:[2] :18
h
u
:=
|
g
→
u
|
=
c
cosh
(
u
)
2
−
cos
(
ψ
)
2
,
h
ψ
:=
|
g
→
ψ
|
=
h
u
:=
h
{\displaystyle h_{u}:=|{\vec {g}}_{u}|=c{\sqrt {\cosh(u)^{2}-\cos(\psi )^{2}}},\quad h_{\psi }:=|{\vec {g}}_{\psi }|=h_{u}:=h}
Das elliptische Orthonormalsystem ist dementsprechend
c
^
u
=
c
h
(
sinh
(
u
)
cos
(
ψ
)
cosh
(
u
)
sin
(
ψ
)
)
,
c
^
ψ
=
c
h
(
−
cosh
(
u
)
sin
(
ψ
)
sinh
(
u
)
cos
(
ψ
)
)
{\displaystyle {\hat {c}}_{u}={\frac {c}{h}}{\begin{pmatrix}\sinh(u)\cos(\psi )\\\cosh(u)\sin(\psi )\end{pmatrix}},\quad {\hat {c}}_{\psi }={\frac {c}{h}}{\begin{pmatrix}-\cosh(u)\sin(\psi )\\\sinh(u)\cos(\psi )\end{pmatrix}}}
Das Linien- und Flächenelement ergibt sich zu[2] :18
d
r
→
=
g
→
u
d
u
+
g
→
ψ
d
ψ
d
s
2
:=
|
d
r
→
|
2
=
c
2
(
cosh
(
u
)
2
−
cos
(
ψ
)
2
)
(
d
u
2
+
d
ψ
2
)
d
A
:=
c
2
(
cosh
(
u
)
2
−
cos
(
ψ
)
2
)
d
u
d
ψ
{\displaystyle {\begin{aligned}{\rm {d}}{\vec {r}}=&{\vec {g}}_{u}{\rm {d}}u+{\vec {g}}_{\psi }{\rm {d}}\psi \\{\rm {d}}s^{2}:=&|{\rm {d}}{\vec {r}}|^{2}=c^{2}(\cosh(u)^{2}-\cos(\psi )^{2})({\rm {d}}u^{2}+{\rm {d}}\psi ^{2})\\{\rm {d}}A:=&c^{2}(\cosh(u)^{2}-\cos(\psi )^{2}){\rm {d}}u\,{\rm {d}}\psi \end{aligned}}}
Die üblichen Differentialoperatoren führt die Tabelle auf[2] :18
(
h
=
c
cosh
(
u
)
2
−
cos
(
ψ
)
2
,
v
→
=
v
u
c
^
u
+
v
ψ
c
^
ψ
)
{\displaystyle (h=c{\sqrt {\cosh(u)^{2}-\cos(\psi )^{2}}},{\vec {v}}=v_{u}{\hat {c}}_{u}+v_{\psi }{\hat {c}}_{\psi })}
Gradient :
g
r
a
d
f
=
1
h
(
c
^
u
∂
f
∂
u
+
c
^
ψ
∂
f
∂
ψ
)
{\displaystyle {\rm {grad}}\,f={\frac {1}{h}}\left({\hat {c}}_{u}{\frac {\partial f}{\partial u}}+{\hat {c}}_{\psi }{\frac {\partial f}{\partial \psi }}\right)}
Divergenz :
d
i
v
v
→
=
1
h
2
(
∂
(
h
v
u
)
∂
u
+
∂
(
h
v
ψ
)
∂
ψ
)
{\displaystyle {\rm {div}}\,{\vec {v}}={\frac {1}{h^{2}}}\left({\frac {\partial (hv_{u})}{\partial u}}+{\frac {\partial (hv_{\psi })}{\partial \psi }}\right)}
Rotation :
r
o
t
v
→
=
1
h
2
(
∂
(
h
v
ψ
)
∂
u
−
∂
(
h
v
u
)
∂
ψ
)
{\displaystyle {\rm {rot}}\,{\vec {v}}={\frac {1}{h^{2}}}\left({\frac {\partial (hv_{\psi })}{\partial u}}-{\frac {\partial (hv_{u})}{\partial \psi }}\right)}
Laplace-Operator :
Δ
f
=
1
h
2
(
∂
2
f
∂
u
2
+
∂
2
f
∂
ψ
2
)
{\displaystyle \Delta f={\frac {1}{h^{2}}}\left({\frac {\partial ^{2}f}{\partial u^{2}}}+{\frac {\partial ^{2}f}{\partial \psi ^{2}}}\right)}
Lösung der Laplace- und Helmholtz-Gleichung in der Ebene
Bearbeiten
Die besondere Form des Laplace-Operators erlaubt eine Lösung der Helmholtz-Gleichung durch multiplikative Trennung der Veränderlichen gemäß dem Separationsansatz [2] :20
ϕ
(
u
,
ψ
)
=
U
(
u
)
⋅
Ψ
(
ψ
)
{\displaystyle \phi (u,\psi )=U(u)\cdot \Psi (\psi )}
Mit obigem Laplace-Operator entsteht die Helmholtz-Gleichung:
Δ
ϕ
(
u
,
ψ
)
=
1
c
2
[
cosh
(
u
)
2
−
cos
(
ψ
)
2
]
(
d
2
U
d
u
2
Ψ
+
U
d
2
Ψ
d
ψ
2
)
=
λ
⋅
U
⋅
Ψ
{\displaystyle \Delta \phi (u,\psi )={\frac {1}{c^{2}[\cosh(u)^{2}-\cos(\psi )^{2}]}}\left({\frac {{\rm {d}}^{2}U}{{\rm {d}}u^{2}}}\Psi +U{\frac {{\rm {d}}^{2}\Psi }{{\rm {d}}\psi ^{2}}}\right)=\lambda \cdot U\cdot \Psi }
Multiplikation beider Seiten mit
c
2
[
cosh
(
u
)
2
−
cos
(
ψ
)
2
]
U
⋅
Ψ
{\displaystyle {\tfrac {c^{2}[\cosh(u)^{2}-\cos(\psi )^{2}]}{U\cdot \Psi }}}
liefert umgestellt
λ
c
2
cosh
(
u
)
2
−
d
2
U
d
u
2
U
=
λ
c
2
cos
(
ψ
)
2
+
d
2
Ψ
d
ψ
2
Ψ
{\displaystyle \lambda c^{2}\cosh(u)^{2}-{\frac {\frac {{\rm {d}}^{2}U}{{\rm {d}}u^{2}}}{U}}=\lambda c^{2}\cos(\psi )^{2}+{\frac {\frac {{\rm {d}}^{2}\Psi }{{\rm {d}}\psi ^{2}}}{\Psi }}}
Weil die linke Seite nur von u und die rechte nur von ψ abhängt, stehen auf beiden Seiten Konstanten:
λ
c
2
cosh
(
u
)
2
−
d
2
U
d
u
2
U
=
κ
2
→
d
2
U
d
u
2
+
[
κ
2
−
λ
c
2
cosh
(
u
)
2
]
U
=
0
λ
c
2
cos
(
ψ
)
2
+
d
2
Ψ
d
ψ
2
Ψ
=
κ
2
→
d
2
Ψ
d
ψ
2
−
[
κ
2
−
λ
c
2
cos
(
ψ
)
2
]
Ψ
=
0
{\displaystyle {\begin{aligned}\lambda c^{2}\cosh(u)^{2}-{\frac {\frac {{\rm {d}}^{2}U}{{\rm {d}}u^{2}}}{U}}=\kappa ^{2}\rightarrow {\frac {{\rm {d}}^{2}U}{{\rm {d}}u^{2}}}+[\kappa ^{2}-\lambda c^{2}\cosh(u)^{2}]U=0\\\lambda c^{2}\cos(\psi )^{2}+{\frac {\frac {{\rm {d}}^{2}\Psi }{{\rm {d}}\psi ^{2}}}{\Psi }}=\kappa ^{2}\rightarrow {\frac {{\rm {d}}^{2}\Psi }{{\rm {d}}\psi ^{2}}}-[\kappa ^{2}-\lambda c^{2}\cos(\psi )^{2}]\Psi =0\end{aligned}}}
Im Fall der Laplace-Gleichung ist λ=0 und die Lösungsfunktion kann mit dem Sinus und Cosinus sowie dem Sinus hyperbolicus und Kosinus hyperbolicus ausgedrückt werden:
ϕ
(
u
,
ψ
)
=
[
A
sin
(
κ
u
)
+
B
cos
(
κ
u
)
]
[
C
sinh
(
κ
ψ
)
+
D
cosh
(
κ
ψ
)
]
{\displaystyle \phi (u,\psi )=[A\sin(\kappa u)+B\cos(\kappa u)][C\sinh(\kappa \psi )+D\cosh(\kappa \psi )]}
Die Konstanten A, B, C, D und κ dienen der Anpassung an Randbedingungen . Wenn die Separationskonstante κ2 mit negativem Vorzeichen angesetzt wird, vertauschen sich in der Lösungsfunktion die Winkelfunktionen durch die Hyperbelfunktionen und umgekehrt.
Koordinatenflächen in abgeplatteten Sphäroidkoordinaten.
Sphäroidkoordinaten sind orthogonale Koordinaten , in denen ein Punkt im Raum durch Angabe der Lage auf konfokalen Rotationsellipsoiden (rot im Bild), Rotationshyperboloiden (blau) und einer Halbebene (gelb) bestimmt wird, siehe Bild. Diese Koordinaten gibt es in zwei Varianten:
Gestreckte Sphäroidkoordinaten
Hier wird die Ellipse um ihre große Halbachse rotiert und das Rotationshyperboloid ist zweischalig, siehe unten.
Abgeplattete Sphäroidkoordinaten
Hier wird die Ellipse wie im Bild um ihre kleine Halbachse rotiert und das Rotationshyperboloid ist einschalig.
Beide Formen werden in zwei verschiedenen Parametrisierungen der Rotationsflächen benutzt.
Diese Koordinaten bieten sich zur Lösung von Randwertaufgaben dort an, wo die Ränder des Gebiets Rotationsflächen von Ellipsen oder Hyperbeln sind.
Gestreckte Sphäroidkoordinaten, Variante 1
Bearbeiten
Koordinatenflächen in gestreckten Sphäroidkoordinaten.
Die kartesischen Koordinaten
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
berechnen sich aus den gestreckten Sphäroidkoordinaten
η
,
θ
,
ψ
∈
R
≥
0
,
θ
≤
π
,
ψ
≤
2
π
{\displaystyle \eta ,\theta ,\psi \in \mathbb {R} ^{\geq 0},\,\theta \leq \pi ,\,\psi \leq 2\pi }
gemäß:[2] :28
r
→
:=
(
x
y
z
)
=
a
(
sinh
(
η
)
sin
(
θ
)
cos
(
ψ
)
sinh
(
η
)
sin
(
θ
)
sin
(
ψ
)
cosh
(
η
)
cos
(
θ
)
)
{\displaystyle {\vec {r}}:={\begin{pmatrix}x\\y\\z\end{pmatrix}}=a{\begin{pmatrix}\sinh(\eta )\sin(\theta )\cos(\psi )\\\sinh(\eta )\sin(\theta )\sin(\psi )\\\cosh(\eta )\cos(\theta )\end{pmatrix}}}
Koordinatenflächen in gestreckten Sphäroidkoordinaten, Variante 1
Bearbeiten
In gestreckten Sphäroidkoordinaten[2] :28 (η,θ,ψ ) bestehen die Koordinatenflächen aus einem Ellipsoid (η =const., rot im Bild),
z
2
(
a
cosh
η
)
2
+
x
2
+
y
2
(
a
sinh
η
)
2
=
1
{\displaystyle {\frac {z^{2}}{(a\cosh \eta )^{2}}}+{\frac {x^{2}+y^{2}}{(a\sinh \eta )^{2}}}=1}
einem zweischaligen Rotationshyperboloid (θ =const., blau)
z
2
(
a
cos
θ
)
2
−
x
2
+
y
2
(
a
sin
θ
)
2
=
1
{\displaystyle {\frac {z^{2}}{(a\cos \theta )^{2}}}-{\frac {x^{2}+y^{2}}{(a\sin \theta )^{2}}}=1}
und einer Halbebene (ψ =const., gelb) mit
y
=
x
tan
ψ
{\displaystyle y=x\tan \psi }
Hieraus ergibt sich andererseits
cosh
(
η
)
2
=
x
2
+
y
2
+
z
2
+
a
2
+
(
x
2
+
y
2
+
z
2
+
a
2
)
2
−
4
a
2
z
2
2
a
2
cos
(
θ
)
2
=
x
2
+
y
2
+
z
2
+
a
2
−
(
x
2
+
y
2
+
z
2
+
a
2
)
2
−
4
a
2
z
2
2
a
2
tan
(
ψ
)
=
y
x
{\displaystyle {\begin{aligned}\cosh(\eta )^{2}=&{\frac {x^{2}+y^{2}+z^{2}+a^{2}+{\sqrt {(x^{2}+y^{2}+z^{2}+a^{2})^{2}-4a^{2}z^{2}}}}{2a^{2}}}\\\cos(\theta )^{2}=&{\frac {x^{2}+y^{2}+z^{2}+a^{2}-{\sqrt {(x^{2}+y^{2}+z^{2}+a^{2})^{2}-4a^{2}z^{2}}}}{2a^{2}}}\\\tan(\psi )=&{\frac {y}{x}}\end{aligned}}}
Metrische Faktoren, Weg- und Flächen- und Volumenelemente in gestreckten Sphäroidkoordinaten, Variante 1
Bearbeiten
Die kovarianten Basisvektoren sind mit
r
→
=
(
x
,
y
,
z
)
⊤
{\displaystyle {\vec {r}}=(x,y,z)^{\top }}
:
g
→
η
:=
∂
r
→
∂
η
=
(
x
coth
(
η
)
y
coth
(
η
)
z
tanh
(
η
)
)
,
g
→
θ
:=
∂
r
→
∂
θ
=
(
x
cot
(
θ
)
y
cot
(
θ
)
−
z
tan
(
θ
)
)
,
g
→
ψ
:=
∂
r
→
∂
ψ
=
(
−
y
x
0
)
{\displaystyle {\vec {g}}_{\eta }:={\frac {\partial {\vec {r}}}{\partial \eta }}={\begin{pmatrix}x\coth(\eta )\\y\coth(\eta )\\z\tanh(\eta )\end{pmatrix}},\;{\vec {g}}_{\theta }:={\frac {\partial {\vec {r}}}{\partial \theta }}={\begin{pmatrix}x\cot(\theta )\\y\cot(\theta )\\-z\tan(\theta )\end{pmatrix}},\;{\vec {g}}_{\psi }:={\frac {\partial {\vec {r}}}{\partial \psi }}={\begin{pmatrix}-y\\x\\0\end{pmatrix}}}
worin cot und coth die Kehrwerte des Tangens bzw. Tangens Hyperbolikus sind. Die Basisvektoren sind senkrecht zueinander und bilden in dieser Reihenfolge im gesamten Wertebereich ein Rechtssystem . Die metrischen Faktoren sind die Beträge der kovarianten Basisvektoren und lauten:[2] :28
h
:=
h
η
=
h
θ
=
a
sinh
(
η
)
2
+
sin
(
θ
)
2
,
h
ψ
=
a
sinh
(
η
)
sin
(
θ
)
{\displaystyle h:=h_{\eta }=h_{\theta }=a{\sqrt {\sinh(\eta )^{2}+\sin(\theta )^{2}}},\,h_{\psi }=a\sinh(\eta )\sin(\theta )}
Das gestreckt sphäroidische Orthonormalsystem ist dementsprechend
c
^
η
=
1
h
(
x
coth
(
η
)
y
coth
(
η
)
z
tanh
(
η
)
)
,
c
^
θ
=
1
h
(
x
cot
(
θ
)
y
cot
(
θ
)
−
z
tan
(
θ
)
)
,
c
^
ψ
=
(
−
sin
(
ψ
)
cos
(
ψ
)
0
)
{\displaystyle {\begin{aligned}{\hat {c}}_{\eta }={\frac {1}{h}}{\begin{pmatrix}x\coth(\eta )\\y\coth(\eta )\\z\tanh(\eta )\end{pmatrix}},\;{\hat {c}}_{\theta }={\frac {1}{h}}{\begin{pmatrix}x\cot(\theta )\\y\cot(\theta )\\-z\tan(\theta )\end{pmatrix}},\;{\hat {c}}_{\psi }={\begin{pmatrix}-\sin(\psi )\\\cos(\psi )\\0\end{pmatrix}}\end{aligned}}}
Das Linien-, Flächen- und Volumenelement ergibt sich zu[2] :28
d
r
→
=
g
→
η
d
η
+
g
→
θ
d
θ
+
g
→
ψ
d
ψ
d
s
2
:=
|
d
r
→
|
2
=
a
2
[
sinh
(
η
)
2
+
sin
(
θ
)
2
]
(
d
η
2
+
d
θ
2
)
+
a
2
sinh
(
η
)
2
sin
(
θ
)
2
d
ψ
2
d
A
=
a
2
{
[
sinh
(
η
)
2
+
sin
(
θ
)
2
]
c
^
ψ
d
η
d
θ
+
sinh
(
η
)
2
+
sin
(
θ
)
2
sinh
(
η
)
sin
(
θ
)
(
c
^
η
d
θ
+
c
^
θ
d
η
)
d
ψ
}
d
V
=
a
3
[
sinh
(
η
)
2
+
sin
(
θ
)
2
]
sinh
(
η
)
sin
(
θ
)
d
η
d
θ
d
ψ
{\displaystyle {\begin{aligned}{\rm {d}}{\vec {r}}=&{\vec {g}}_{\eta }{\rm {d}}\eta +{\vec {g}}_{\theta }{\rm {d}}\theta +{\vec {g}}_{\psi }{\rm {d}}\psi \\{\rm {d}}s^{2}:=&|{\rm {d}}{\vec {r}}|^{2}=a^{2}[\sinh(\eta )^{2}+\sin(\theta )^{2}]\,({\rm {d}}\eta ^{2}+{\rm {d}}\theta ^{2})+a^{2}\sinh(\eta )^{2}\sin(\theta )^{2}\,{\rm {d}}\psi ^{2}\\{\rm {d}}A=&a^{2}\left\{[\sinh(\eta )^{2}+\sin(\theta )^{2}]{\hat {c}}_{\psi }\,{\rm {d}}\eta \,{\rm {d}}\theta +{\sqrt {\sinh(\eta )^{2}+\sin(\theta )^{2}}}\sinh(\eta )\sin(\theta )({\hat {c}}_{\eta }\,{\rm {d}}\theta +{\hat {c}}_{\theta }\,{\rm {d}}\eta )\,{\rm {d}}\psi \right\}\\{\rm {d}}V=&a^{3}[\sinh(\eta )^{2}+\sin(\theta )^{2}]\sinh(\eta )\sin(\theta )\,{\rm {d}}\eta \,{\rm {d}}\theta \,{\rm {d}}\psi \end{aligned}}}
Differentialoperatoren in gestreckten Sphäroidkoordinaten, Variante 1
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:[2] :29
Δ
f
=
1
a
2
[
sin
(
θ
)
2
+
sinh
(
η
)
2
]
(
∂
2
f
∂
η
2
+
coth
(
η
)
∂
f
∂
η
+
∂
2
f
∂
θ
2
+
cot
(
θ
)
∂
f
∂
θ
)
+
…
⋯
+
1
a
2
sinh
(
η
)
2
sin
(
θ
)
2
∂
2
f
∂
ψ
2
{\displaystyle {\begin{aligned}\Delta f=&{\frac {1}{a^{2}[\sin(\theta )^{2}+\sinh(\eta )^{2}]}}\left({\frac {\partial ^{2}f}{\partial \eta ^{2}}}+\coth(\eta ){\frac {\partial f}{\partial \eta }}+{\frac {\partial ^{2}f}{\partial \theta ^{2}}}+\cot(\theta ){\frac {\partial f}{\partial \theta }}\right)+\dots \\&\dots +{\frac {1}{a^{2}\sinh(\eta )^{2}\sin(\theta )^{2}}}{\frac {\partial ^{2}f}{\partial \psi ^{2}}}\end{aligned}}}
worin cot und coth die Kehrwerte des Tangens bzw. Tangens Hyperbolikus sind.
Lösung der Laplace- und Helmholtz-Gleichung in gestreckten Sphäroidkoordinaten, Variante 1
Bearbeiten
In Ellipsoid-Koordinaten gelingt immer eine Trennung der Variablen in der Laplace- und Helmholtz-Gleichung . Mit dem Separationsansatz
ϕ
(
η
,
θ
,
ψ
)
=
H
(
η
)
⋅
Θ
(
θ
)
⋅
Ψ
(
ψ
)
{\displaystyle \phi (\eta ,\theta ,\psi )=H(\eta )\cdot \Theta (\theta )\cdot \Psi (\psi )}
lautet die Helmholtz-Gleichung
Δ
ϕ
+
κ
2
ϕ
=
0
{\displaystyle \Delta \phi +\kappa ^{2}\phi =0}
:
1
a
2
[
sin
(
θ
)
2
+
sinh
(
η
)
2
]
(
∂
2
H
∂
η
2
Θ
Ψ
+
coth
(
η
)
∂
H
∂
η
Θ
Ψ
+
H
∂
2
Θ
∂
θ
2
Ψ
+
cot
(
θ
)
H
∂
Θ
∂
θ
Ψ
)
+
…
⋯
+
1
a
2
sinh
(
η
)
2
sin
(
θ
)
2
H
Θ
∂
2
Ψ
∂
ψ
2
+
κ
2
H
Θ
Ψ
=
0
{\displaystyle {\begin{aligned}{\frac {1}{a^{2}[\sin(\theta )^{2}+\sinh(\eta )^{2}]}}\left({\frac {\partial ^{2}H}{\partial \eta ^{2}}}\Theta \Psi +\coth(\eta ){\frac {\partial H}{\partial \eta }}\Theta \Psi +H{\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}\Psi +\cot(\theta )H{\frac {\partial \Theta }{\partial \theta }}\Psi \right)+\dots &\\\dots +{\frac {1}{a^{2}\sinh(\eta )^{2}\sin(\theta )^{2}}}H\Theta {\frac {\partial ^{2}\Psi }{\partial \psi ^{2}}}+\kappa ^{2}H\Theta \Psi &=0\end{aligned}}}
Die drei Faktoren in der separierten Lösungsfunktion bestimmen sich aus den unabhängigen gewöhnlichen Differenzialgleichungen[2] :30
∂
2
H
∂
η
2
+
coth
(
η
)
∂
H
∂
η
+
(
κ
2
a
2
sinh
2
η
−
α
2
−
α
3
sinh
2
η
)
H
=
0
∂
2
Θ
∂
θ
2
+
cot
(
θ
)
∂
Θ
∂
θ
+
(
κ
2
a
2
sin
2
θ
+
α
2
−
α
3
sin
2
θ
)
Θ
=
0
∂
2
Ψ
∂
ψ
2
+
α
3
Ψ
=
0
{\displaystyle {\begin{aligned}{\frac {\partial ^{2}H}{\partial \eta ^{2}}}+\coth(\eta ){\frac {\partial H}{\partial \eta }}+\left(\kappa ^{2}a^{2}\sinh ^{2}\eta -\alpha _{2}-{\frac {\alpha _{3}}{\sinh ^{2}\eta }}\right)H=&0\\{\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}+\cot(\theta ){\frac {\partial \Theta }{\partial \theta }}+\left(\kappa ^{2}a^{2}\sin ^{2}\theta +\alpha _{2}-{\frac {\alpha _{3}}{\sin ^{2}\theta }}\right)\Theta =&0\\{\frac {\partial ^{2}\Psi }{\partial \psi ^{2}}}+\alpha _{3}\Psi =&0\end{aligned}}}
und den Randbedingungen . Bei der Laplace-Gleichung ist
κ
=
0
{\displaystyle \kappa =0}
.
Denn Multiplikation der Helmholtz-Gleichung mit
a
2
sinh
(
η
)
2
sin
(
θ
)
2
H
Θ
Ψ
{\displaystyle {\tfrac {a^{2}\sinh(\eta )^{2}\sin(\theta )^{2}}{H\Theta \Psi }}}
liefert:
sinh
(
η
)
2
sin
(
θ
)
2
sin
(
θ
)
2
+
sinh
(
η
)
2
(
∂
2
H
∂
η
2
H
+
coth
(
η
)
∂
H
∂
η
H
+
∂
2
Θ
∂
θ
2
Θ
+
cot
(
θ
)
∂
Θ
∂
θ
Θ
)
+
…
⋯
+
κ
2
a
2
sinh
(
η
)
2
sin
(
θ
)
2
+
∂
2
Ψ
∂
ψ
2
Ψ
=
0
{\displaystyle {\begin{aligned}{\frac {\sinh(\eta )^{2}\sin(\theta )^{2}}{\sin(\theta )^{2}+\sinh(\eta )^{2}}}\left({\frac {\frac {\partial ^{2}H}{\partial \eta ^{2}}}{H}}+\coth(\eta ){\frac {\frac {\partial H}{\partial \eta }}{H}}+{\frac {\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}{\Theta }}+\cot(\theta ){\frac {\frac {\partial \Theta }{\partial \theta }}{\Theta }}\right)+\dots &\\\dots +\kappa ^{2}a^{2}\sinh(\eta )^{2}\sin(\theta )^{2}+{\frac {\frac {\partial ^{2}\Psi }{\partial \psi ^{2}}}{\Psi }}&=0\end{aligned}}}
Weil auf der rechten Seite eine Konstante (null) steht und nur der letzte Term auf der linken Seite von ψ abhängt, ist letzterer ebenfalls konstant:
∂
2
Ψ
∂
ψ
2
Ψ
=
−
α
3
{\displaystyle {\frac {\frac {\partial ^{2}\Psi }{\partial \psi ^{2}}}{\Psi }}=-\alpha _{3}}
Einsetzen dieser Konstanten erlaubt auch η und θ zu trennen:
−
[
∂
2
Θ
∂
θ
2
Θ
+
cot
(
θ
)
∂
Θ
∂
θ
Θ
+
κ
2
a
2
sin
(
θ
)
2
−
α
3
sin
(
θ
)
2
]
=
…
⋯
=
∂
2
H
∂
η
2
H
+
coth
(
η
)
∂
H
∂
η
H
+
κ
2
a
2
sinh
(
η
)
2
−
α
3
sinh
(
η
)
2
{\displaystyle {\begin{aligned}&-\left[{\frac {\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}{\Theta }}+\cot(\theta ){\frac {\frac {\partial \Theta }{\partial \theta }}{\Theta }}+\kappa ^{2}a^{2}\sin(\theta )^{2}-{\frac {\alpha _{3}}{\sin(\theta )^{2}}}\right]=\dots \\&\qquad \qquad \qquad \dots ={\frac {\frac {\partial ^{2}H}{\partial \eta ^{2}}}{H}}+\coth(\eta ){\frac {\frac {\partial H}{\partial \eta }}{H}}+\kappa ^{2}a^{2}\sinh(\eta )^{2}-{\frac {\alpha _{3}}{\sinh(\eta )^{2}}}\end{aligned}}}
Weil die linke Seite nur von θ abhängt und die rechte Seite nur von η , sind beide Seiten gleich einer Konstanten α 2 . So entstehen die oben angegebene Differenzialgleichungen[2] :30 Ein gleichbedeutendes Ergebnis wird mit dem im Hauptartikel beschriebenen Verfahren und der Stäckel-Matrix
S
=
(
a
2
sinh
(
η
)
2
−
1
−
1
sinh
(
η
)
2
a
2
sin
(
θ
)
2
1
−
1
sin
(
θ
)
2
0
0
1
)
{\displaystyle \mathbf {S} ={\begin{pmatrix}a^{2}\sinh(\eta )^{2}&-1&{\frac {-1}{\sinh(\eta )^{2}}}\\a^{2}\sin(\theta )^{2}&1&{\frac {-1}{\sin(\theta )^{2}}}\\0&0&1\end{pmatrix}}}
erzielt.
Gestreckte Sphäroidkoordinaten, Variante 2
Bearbeiten
Die gestreckten Sphäroidkoordinaten der Variante 2 benutzen nicht die Variablen η,θ,ψ der ersten Variante, sondern deren Funktionswerte:[3] :661
ξ
1
=
a
cosh
(
η
)
,
ξ
2
=
cos
(
θ
)
,
ξ
3
=
cos
(
ψ
)
{\displaystyle \xi _{1}=a\cosh(\eta ),\;\xi _{2}=\cos(\theta ),\;\xi _{3}=\cos(\psi )}
Die kartesischen Koordinaten
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
berechnen sich aus den gestreckten Sphäroidkoordinaten
ξ
1
,
2
,
3
∈
R
,
ξ
1
≥
0
,
−
1
≤
ξ
2
,
3
≤
1
{\displaystyle \xi _{1,2,3}\in \mathbb {R} ,\,\xi _{1}\geq 0,-1\leq \xi _{2,3}\leq 1}
gemäß:
r
→
:=
(
x
y
z
)
=
(
ξ
3
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
(
1
−
ξ
3
2
)
ξ
1
ξ
2
)
{\displaystyle {\vec {r}}:={\begin{pmatrix}x\\y\\z\end{pmatrix}}={\begin{pmatrix}\xi _{3}{\sqrt {(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})}}\\{\sqrt {(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})(1-\xi _{3}^{2})}}\\\xi _{1}\xi _{2}\end{pmatrix}}}
Koordinatenflächen in gestreckten Sphäroidkoordinaten, Variante 2
Bearbeiten
Koordinatenflächen in gestreckten Sphäroidkoordinaten.
In der Variante 2 ist auf dem Ellipsoid (ξ 1 =const.,rot)
z
2
ξ
1
2
+
x
2
+
y
2
ξ
1
2
−
a
2
=
1
{\displaystyle {\frac {z^{2}}{\xi _{1}^{2}}}+{\frac {x^{2}+y^{2}}{\xi _{1}^{2}-a^{2}}}=1}
auf dem zweischaligen Rotationshyperboloid (ξ 2 =const., blau)
z
2
a
2
ξ
2
2
−
x
2
+
y
2
a
2
(
1
−
ξ
2
2
)
=
1
{\displaystyle {\frac {z^{2}}{a^{2}\xi _{2}^{2}}}-{\frac {x^{2}+y^{2}}{a^{2}(1-\xi _{2}^{2})}}=1}
und in der Halbebene (ξ 3 =const., gelb)
ξ
3
1
+
(
x
y
)
2
−
x
y
=
0
{\displaystyle \xi _{3}{\sqrt {1+\left({\frac {x}{y}}\right)^{2}}}-{\frac {x}{y}}=0}
Hieraus ergibt sich andererseits
ξ
1
2
=
x
2
+
y
2
+
z
2
+
a
2
+
(
x
2
+
y
2
+
z
2
+
a
2
)
2
−
4
a
2
z
2
2
ξ
2
2
=
x
2
+
y
2
+
z
2
+
a
2
−
(
x
2
+
y
2
+
z
2
+
a
2
)
2
−
4
a
2
z
2
2
a
2
ξ
3
2
=
x
2
x
2
+
y
2
{\displaystyle {\begin{aligned}\xi _{1}^{2}=&{\frac {x^{2}+y^{2}+z^{2}+a^{2}+{\sqrt {(x^{2}+y^{2}+z^{2}+a^{2})^{2}-4a^{2}z^{2}}}}{2}}\\\xi _{2}^{2}=&{\frac {x^{2}+y^{2}+z^{2}+a^{2}-{\sqrt {(x^{2}+y^{2}+z^{2}+a^{2})^{2}-4a^{2}z^{2}}}}{2a^{2}}}\\\xi _{3}^{2}=&{\frac {x^{2}}{x^{2}+y^{2}}}\end{aligned}}}
Metrische Faktoren, Weg- und Flächen- und Volumenelemente in gestreckten Sphäroidkoordinaten, Variante 2
Bearbeiten
Die kovarianten Basisvektoren sind mit
r
→
=
(
x
,
y
,
z
)
⊤
{\displaystyle {\vec {r}}=(x,y,z)^{\top }}
:
g
→
1
:=
∂
r
→
∂
ξ
1
=
(
ξ
1
x
ξ
1
2
−
a
2
ξ
1
y
ξ
1
2
−
a
2
z
ξ
1
)
,
g
→
2
:=
∂
r
→
∂
ξ
2
=
(
−
ξ
2
x
1
−
ξ
2
2
−
ξ
2
y
1
−
ξ
2
2
z
ξ
2
)
,
g
→
3
:=
∂
r
→
∂
ξ
3
=
(
x
ξ
3
−
ξ
3
y
1
−
ξ
3
2
0
)
{\displaystyle {\vec {g}}_{1}:={\frac {\partial {\vec {r}}}{\partial \xi _{1}}}={\begin{pmatrix}{\frac {\xi _{1}x}{\xi _{1}^{2}-a^{2}}}\\{\frac {\xi _{1}y}{\xi _{1}^{2}-a^{2}}}\\{\frac {z}{\xi _{1}}}\end{pmatrix}},\;{\vec {g}}_{2}:={\frac {\partial {\vec {r}}}{\partial \xi _{2}}}={\begin{pmatrix}-{\frac {\xi _{2}x}{1-\xi _{2}^{2}}}\\-{\frac {\xi _{2}y}{1-\xi _{2}^{2}}}\\{\frac {z}{\xi _{2}}}\end{pmatrix}},\;{\vec {g}}_{3}:={\frac {\partial {\vec {r}}}{\partial \xi _{3}}}={\begin{pmatrix}{\frac {x}{\xi _{3}}}\\-{\frac {\xi _{3}y}{1-\xi _{3}^{2}}}\\0\end{pmatrix}}}
Die Basisvektoren sind senkrecht zueinander und bilden in dieser Reihenfolge ein Rechtssystem . Die metrischen Faktoren sind die Beträge der kovarianten Basisvektoren und lauten:
h
1
=
ξ
1
2
−
a
2
ξ
2
2
ξ
1
2
−
a
2
,
h
2
=
ξ
1
2
−
a
2
ξ
2
2
1
−
ξ
2
2
,
h
3
=
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
1
−
ξ
3
2
{\displaystyle h_{1}={\sqrt {\frac {\xi _{1}^{2}-a^{2}\xi _{2}^{2}}{\xi _{1}^{2}-a^{2}}}},\;h_{2}={\sqrt {\frac {\xi _{1}^{2}-a^{2}\xi _{2}^{2}}{1-\xi _{2}^{2}}}},\;h_{3}={\sqrt {\frac {(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})}{1-\xi _{3}^{2}}}}}
Das gestreckt sphäroidische Orthonormalsystem schreibt sich in dieser Formulierung
c
^
1
=
1
ξ
1
2
−
a
2
ξ
2
2
(
ξ
1
ξ
3
1
−
ξ
2
2
ξ
1
(
1
−
ξ
2
2
)
(
1
−
ξ
3
2
)
ξ
2
ξ
1
2
−
a
2
)
c
^
2
=
1
ξ
1
2
−
a
2
ξ
2
2
(
−
ξ
2
ξ
3
ξ
1
2
−
a
2
−
ξ
2
(
ξ
1
2
−
a
2
)
(
1
−
ξ
3
2
)
ξ
1
1
−
ξ
2
2
)
,
c
^
3
=
(
1
−
ξ
3
2
−
ξ
3
0
)
{\displaystyle {\begin{aligned}{\hat {c}}_{1}=&{\frac {1}{\sqrt {\xi _{1}^{2}-a^{2}\xi _{2}^{2}}}}{\begin{pmatrix}\xi _{1}\xi _{3}{\sqrt {1-\xi _{2}^{2}}}\\\xi _{1}{\sqrt {(1-\xi _{2}^{2})(1-\xi _{3}^{2})}}\\\xi _{2}{\sqrt {\xi _{1}^{2}-a^{2}}}\end{pmatrix}}\\{\hat {c}}_{2}=&{\frac {1}{\sqrt {\xi _{1}^{2}-a^{2}\xi _{2}^{2}}}}{\begin{pmatrix}-\xi _{2}\xi _{3}{\sqrt {\xi _{1}^{2}-a^{2}}}\\-\xi _{2}{\sqrt {(\xi _{1}^{2}-a^{2})(1-\xi _{3}^{2})}}\\\xi _{1}{\sqrt {1-\xi _{2}^{2}}}\end{pmatrix}},\;{\hat {c}}_{3}={\begin{pmatrix}{\sqrt {1-\xi _{3}^{2}}}\\-\xi _{3}\\0\end{pmatrix}}\end{aligned}}}
Das Linien- und Volumenelement ergibt sich zu
d
r
→
=
g
→
1
d
ξ
1
+
g
→
2
d
ξ
2
+
g
→
3
d
ξ
3
d
s
2
:=
|
d
r
→
|
2
=
ξ
1
2
−
a
2
ξ
2
2
ξ
1
2
−
a
2
d
ξ
1
2
+
ξ
1
2
−
a
2
ξ
2
2
1
−
ξ
2
2
d
ξ
2
2
+
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
1
−
ξ
3
2
d
ξ
3
2
d
V
=
ξ
1
2
−
a
2
ξ
2
2
1
−
ξ
3
2
d
ξ
1
d
ξ
2
d
ξ
3
{\displaystyle {\begin{aligned}{\rm {d}}{\vec {r}}=&{\vec {g}}_{1}{\rm {d}}\xi _{1}+{\vec {g}}_{2}{\rm {d}}\xi _{2}+{\vec {g}}_{3}{\rm {d}}\xi _{3}\\{\rm {d}}s^{2}:=&|{\rm {d}}{\vec {r}}|^{2}={\frac {\xi _{1}^{2}-a^{2}\xi _{2}^{2}}{\xi _{1}^{2}-a^{2}}}\,{\rm {d}}\xi _{1}^{2}+{\frac {\xi _{1}^{2}-a^{2}\xi _{2}^{2}}{1-\xi _{2}^{2}}}\,{\rm {d}}\xi _{2}^{2}+{\frac {(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})}{1-\xi _{3}^{2}}}\,{\rm {d}}\xi _{3}^{2}\\{\rm {d}}V=&{\frac {\xi _{1}^{2}-a^{2}\xi _{2}^{2}}{\sqrt {1-\xi _{3}^{2}}}}\,{\rm {d}}\xi _{1}\,{\rm {d}}\xi _{2}\,{\rm {d}}\xi _{3}\end{aligned}}}
Differentialoperatoren in gestreckten Sphäroidkoordinaten, Variante 2
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:
Δ
f
=
(
ξ
1
2
−
a
2
)
∂
2
f
∂
ξ
1
2
+
2
ξ
1
∂
f
∂
ξ
1
+
(
1
−
ξ
2
2
)
∂
2
f
∂
ξ
2
2
−
2
ξ
2
∂
f
∂
ξ
2
ξ
1
2
−
a
2
ξ
2
2
+
(
1
−
ξ
3
2
)
∂
2
f
∂
ξ
3
2
−
ξ
3
∂
f
∂
ξ
3
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
{\displaystyle \Delta f={\frac {(\xi _{1}^{2}-a^{2}){\frac {\partial ^{2}f}{\partial \xi _{1}^{2}}}+2\xi _{1}{\frac {\partial f}{\partial \xi _{1}}}+(1-\xi _{2}^{2}){\frac {\partial ^{2}f}{\partial \xi _{2}^{2}}}-2\xi _{2}{\frac {\partial f}{\partial \xi _{2}}}}{\xi _{1}^{2}-a^{2}\xi _{2}^{2}}}+{\frac {(1-\xi _{3}^{2}){\frac {\partial ^{2}f}{\partial \xi _{3}^{2}}}-\xi _{3}{\frac {\partial f}{\partial \xi _{3}}}}{(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})}}}
Lösung der Laplace- und Helmholtz-Gleichung in gestreckten Sphäroidkoordinaten, Variante 2
Bearbeiten
In Ellipsoid-Koordinaten gelingt immer eine Trennung der Variablen in der Laplace- und Helmholtz-Gleichung . Mit dem Separationsansatz
ϕ
(
ξ
1
,
ξ
2
,
ξ
3
)
=
X
(
ξ
1
)
⋅
Y
(
ξ
2
)
⋅
Z
(
ξ
3
)
{\displaystyle \phi (\xi _{1},\xi _{2},\xi _{3})=X(\xi _{1})\cdot Y(\xi _{2})\cdot Z(\xi _{3})}
lautet die Helmholtz-Gleichung
Δ
ϕ
+
κ
2
ϕ
=
0
{\displaystyle \Delta \phi +\kappa ^{2}\phi =0}
:
(
ξ
1
2
−
a
2
)
∂
2
ϕ
∂
ξ
1
2
+
2
ξ
1
∂
ϕ
∂
ξ
1
+
(
1
−
ξ
2
2
)
∂
2
ϕ
∂
ξ
2
2
−
2
ξ
2
∂
ϕ
∂
ξ
2
ξ
1
2
−
a
2
ξ
2
2
+
(
1
−
ξ
3
2
)
∂
2
ϕ
∂
ξ
3
2
−
ξ
3
∂
ϕ
∂
ξ
3
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
+
κ
2
ϕ
=
(
ξ
1
2
−
a
2
)
∂
2
X
∂
ξ
1
2
Y
Z
+
2
ξ
1
∂
X
∂
ξ
1
Y
Z
+
(
1
−
ξ
2
2
)
X
∂
2
Y
∂
ξ
2
2
Z
−
2
ξ
2
X
∂
Y
∂
ξ
2
Z
ξ
1
2
−
a
2
ξ
2
2
+
…
⋯
+
(
1
−
ξ
3
2
)
X
Y
∂
2
Z
∂
ξ
3
2
−
ξ
3
X
Y
∂
Z
∂
ξ
3
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
+
κ
2
X
Y
Z
=
0
{\displaystyle {\begin{aligned}{\frac {(\xi _{1}^{2}-a^{2}){\frac {\partial ^{2}\phi }{\partial \xi _{1}^{2}}}+2\xi _{1}{\frac {\partial \phi }{\partial \xi _{1}}}+(1-\xi _{2}^{2}){\frac {\partial ^{2}\phi }{\partial \xi _{2}^{2}}}-2\xi _{2}{\frac {\partial \phi }{\partial \xi _{2}}}}{\xi _{1}^{2}-a^{2}\xi _{2}^{2}}}+{\frac {(1-\xi _{3}^{2}){\frac {\partial ^{2}\phi }{\partial \xi _{3}^{2}}}-\xi _{3}{\frac {\partial \phi }{\partial \xi _{3}}}}{(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})}}+\kappa ^{2}\phi &=\\{\frac {(\xi _{1}^{2}-a^{2}){\frac {\partial ^{2}X}{\partial \xi _{1}^{2}}}YZ+2\xi _{1}{\frac {\partial X}{\partial \xi _{1}}}YZ+(1-\xi _{2}^{2})X{\frac {\partial ^{2}Y}{\partial \xi _{2}^{2}}}Z-2\xi _{2}X{\frac {\partial Y}{\partial \xi _{2}}}Z}{\xi _{1}^{2}-a^{2}\xi _{2}^{2}}}+\dots \qquad &\\\dots +{\frac {(1-\xi _{3}^{2})XY{\frac {\partial ^{2}Z}{\partial \xi _{3}^{2}}}-\xi _{3}XY{\frac {\partial Z}{\partial \xi _{3}}}}{(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})}}+\kappa ^{2}XYZ&=0\end{aligned}}}
Die drei Faktoren in der separierten Lösungsfunktion bestimmen sich aus den unabhängigen gewöhnlichen Differenzialgleichungen
2
ξ
1
∂
X
∂
ξ
1
+
(
ξ
1
2
−
a
2
)
∂
2
X
∂
ξ
1
2
+
(
κ
2
ξ
1
2
−
α
2
+
α
3
a
2
ξ
1
2
−
a
2
)
X
=
0
2
ξ
2
∂
Y
∂
ξ
2
−
(
1
−
ξ
2
2
)
∂
2
Y
∂
ξ
2
2
+
(
κ
2
a
2
ξ
2
2
−
α
2
−
α
3
1
−
ξ
2
2
)
Y
=
0
(
1
−
ξ
3
2
)
∂
2
Z
∂
ξ
3
2
−
ξ
3
∂
Z
∂
ξ
3
−
α
3
Z
=
0
{\displaystyle {\begin{aligned}2\xi _{1}{\frac {\partial X}{\partial \xi _{1}}}+(\xi _{1}^{2}-a^{2}){\frac {\partial ^{2}X}{\partial \xi _{1}^{2}}}+\left(\kappa ^{2}\xi _{1}^{2}-\alpha _{2}+{\frac {\alpha _{3}a^{2}}{\xi _{1}^{2}-a^{2}}}\right)X=&0\\2\xi _{2}{\frac {\partial Y}{\partial \xi _{2}}}-(1-\xi _{2}^{2}){\frac {\partial ^{2}Y}{\partial \xi _{2}^{2}}}+\left(\kappa ^{2}a^{2}\xi _{2}^{2}-\alpha _{2}-{\frac {\alpha _{3}}{1-\xi _{2}^{2}}}\right)Y=&0\\(1-\xi _{3}^{2}){\frac {\partial ^{2}Z}{\partial \xi _{3}^{2}}}-\xi _{3}{\frac {\partial Z}{\partial \xi _{3}}}-\alpha _{3}Z=&0\end{aligned}}}
und den Randbedingungen . Bei der Laplace-Gleichung ist
κ
=
0
{\displaystyle \kappa =0}
.
Denn Multiplikation der Helmholtz-Gleichung mit
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
X
Y
Z
{\displaystyle {\tfrac {(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})}{XYZ}}}
liefert:
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
ξ
1
2
−
a
2
ξ
2
2
[
(
ξ
1
2
−
a
2
)
∂
2
X
∂
ξ
1
2
X
+
2
ξ
1
∂
X
∂
ξ
1
X
+
(
1
−
ξ
2
2
)
∂
2
Y
∂
ξ
2
2
Y
−
2
ξ
2
∂
Y
∂
ξ
2
Y
]
+
…
⋯
+
{
(
1
−
ξ
3
2
)
∂
2
Z
∂
ξ
3
2
Z
−
ξ
3
∂
Z
∂
ξ
3
Z
}
+
κ
2
(
ξ
1
2
−
a
2
)
(
1
−
ξ
2
2
)
=
0
{\displaystyle {\begin{aligned}{\frac {(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})}{\xi _{1}^{2}-a^{2}\xi _{2}^{2}}}\left[(\xi _{1}^{2}-a^{2}){\frac {\frac {\partial ^{2}X}{\partial \xi _{1}^{2}}}{X}}+2\xi _{1}{\frac {\frac {\partial X}{\partial \xi _{1}}}{X}}+(1-\xi _{2}^{2}){\frac {\frac {\partial ^{2}Y}{\partial \xi _{2}^{2}}}{Y}}-2\xi _{2}{\frac {\frac {\partial Y}{\partial \xi _{2}}}{Y}}\right]+\dots &\\\dots +\left\{(1-\xi _{3}^{2}){\frac {\frac {\partial ^{2}Z}{\partial \xi _{3}^{2}}}{Z}}-\xi _{3}{\frac {\frac {\partial Z}{\partial \xi _{3}}}{Z}}\right\}+\kappa ^{2}(\xi _{1}^{2}-a^{2})(1-\xi _{2}^{2})&=0\end{aligned}}}
Weil auf der rechten Seite eine Konstante (null) steht und nur der Term in der geschweiften Klammer von ξ 3 abhängt, ist dieser ebenfalls eine Konstante α 3 :
(
1
−
ξ
3
2
)
∂
2
Z
∂
ξ
3
2
Z
−
ξ
3
∂
Z
∂
ξ
3
Z
=
α
3
{\displaystyle (1-\xi _{3}^{2}){\frac {\frac {\partial ^{2}Z}{\partial \xi _{3}^{2}}}{Z}}-\xi _{3}{\frac {\frac {\partial Z}{\partial \xi _{3}}}{Z}}=\alpha _{3}}
Einsetzen dieser Konstanten erlaubt auch ξ 1 und ξ 2 zu trennen:
2
ξ
1
∂
X
∂
ξ
1
X
+
(
ξ
1
2
−
a
2
)
∂
2
X
∂
ξ
1
2
X
+
κ
2
ξ
1
2
+
α
3
a
2
ξ
1
2
−
a
2
=
…
⋯
=
2
ξ
2
∂
Y
∂
ξ
2
Y
−
(
1
−
ξ
2
2
)
∂
2
Y
∂
ξ
2
2
Y
+
κ
2
a
2
ξ
2
2
−
α
3
1
−
ξ
2
2
{\displaystyle {\begin{aligned}&2\xi _{1}{\frac {\frac {\partial X}{\partial \xi _{1}}}{X}}+(\xi _{1}^{2}-a^{2}){\frac {\frac {\partial ^{2}X}{\partial \xi _{1}^{2}}}{X}}+\kappa ^{2}\xi _{1}^{2}+{\frac {\alpha _{3}a^{2}}{\xi _{1}^{2}-a^{2}}}=\dots \\&\qquad \qquad \qquad \dots =2\xi _{2}{\frac {\frac {\partial Y}{\partial \xi _{2}}}{Y}}-(1-\xi _{2}^{2}){\frac {\frac {\partial ^{2}Y}{\partial \xi _{2}^{2}}}{Y}}+\kappa ^{2}a^{2}\xi _{2}^{2}-{\frac {\alpha _{3}}{1-\xi _{2}^{2}}}\end{aligned}}}
Weil die linke Seite nur von ξ 1 abhängt und die rechte Seite nur von ξ 2 , sind beide Seiten gleich einer Konstanten α 2 . Damit resultieren die oben angegebene Differenzialgleichungen. Ein gleichbedeutendes Ergebnis entsteht mit der Stäckel-Matrix[3] :661
S
=
(
1
1
ξ
1
2
−
a
2
a
2
(
ξ
1
2
−
a
2
)
2
a
2
−
1
1
−
ξ
2
2
1
(
1
−
ξ
2
2
)
2
0
0
−
1
1
−
ξ
3
2
)
{\displaystyle \mathbf {S} ={\begin{pmatrix}1&{\frac {1}{\xi _{1}^{2}-a^{2}}}&{\frac {a^{2}}{(\xi _{1}^{2}-a^{2})^{2}}}\\a^{2}&{\frac {-1}{1-\xi _{2}^{2}}}&{\frac {1}{(1-\xi _{2}^{2})^{2}}}\\0&0&{\frac {-1}{1-\xi _{3}^{2}}}\end{pmatrix}}}
und der im Hauptartikel beschriebenen Methode.
Abgeplattete Sphäroidkoordinaten, Variante 1
Bearbeiten
Koordinatenflächen in abgeplatteten Sphäroidkoordinaten.
Die kartesischen Koordinaten
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
berechnen sich aus den abgeplatteten Sphäroidkoordinaten im selben Wertebereich wie die gestreckten Sphäroidkoordinaten der Variante 1
η
,
θ
,
ψ
∈
R
≥
0
,
θ
≤
π
,
ψ
≤
2
π
{\displaystyle \eta ,\theta ,\psi \in \mathbb {R} ^{\geq 0},\,\theta \leq \pi ,\,\psi \leq 2\pi }
gemäß:[2] :31
r
→
:=
(
x
y
z
)
=
a
(
cosh
(
η
)
sin
(
θ
)
cos
(
ψ
)
cosh
(
η
)
sin
(
θ
)
sin
(
ψ
)
sinh
(
η
)
cos
(
θ
)
)
{\displaystyle {\vec {r}}:={\begin{pmatrix}x\\y\\z\end{pmatrix}}=a{\begin{pmatrix}\cosh(\eta )\sin(\theta )\cos(\psi )\\\cosh(\eta )\sin(\theta )\sin(\psi )\\\sinh(\eta )\cos(\theta )\end{pmatrix}}}
#Gestreckte Sphäroidkoordinaten, Variante 1 , benutzen im Vergleich hierzu cosh statt sinh und umgekehrt.
Koordinatenflächen in abgeplatteten Sphäroidkoordinaten, Variante 1
Bearbeiten
In abgeplatteten Sphäroidkoordinaten[2] :31 (η,θ,ψ ) bestehen die Koordinatenflächen aus einem Ellipsoid (η =const., rot im Bild),
x
2
+
y
2
(
a
cosh
η
)
2
+
z
2
(
a
sinh
η
)
2
=
1
{\displaystyle {\frac {x^{2}+y^{2}}{(a\cosh \eta )^{2}}}+{\frac {z^{2}}{(a\sinh \eta )^{2}}}=1}
einem einschaligen Rotationshyperboloid (θ =const., blau)
x
2
+
y
2
(
a
sin
θ
)
2
−
z
2
(
a
cos
θ
)
2
=
1
{\displaystyle {\frac {x^{2}+y^{2}}{(a\sin \theta )^{2}}}-{\frac {z^{2}}{(a\cos \theta )^{2}}}=1}
und einer Halbebene (ψ =const., gelb) mit
y
=
x
tan
ψ
{\displaystyle y=x\tan \psi }
Hieraus ergibt sich andererseits
cosh
(
η
)
2
=
x
2
+
y
2
+
z
2
+
a
2
+
(
x
2
+
y
2
+
z
2
−
a
2
)
2
+
4
a
2
z
2
2
a
2
sin
(
θ
)
2
=
x
2
+
y
2
+
z
2
+
a
2
−
(
x
2
+
y
2
+
z
2
−
a
2
)
2
+
4
a
2
z
2
2
a
2
tan
(
ψ
)
=
y
x
{\displaystyle {\begin{aligned}\cosh(\eta )^{2}=&{\frac {x^{2}+y^{2}+z^{2}+a^{2}+{\sqrt {(x^{2}+y^{2}+z^{2}-a^{2})^{2}+4a^{2}z^{2}}}}{2a^{2}}}\\\sin(\theta )^{2}=&{\frac {x^{2}+y^{2}+z^{2}+a^{2}-{\sqrt {(x^{2}+y^{2}+z^{2}-a^{2})^{2}+4a^{2}z^{2}}}}{2a^{2}}}\\\tan(\psi )=&{\frac {y}{x}}\end{aligned}}}
Metrische Faktoren, Weg- und Flächen- und Volumenelemente in abgeplatteten Sphäroidkoordinaten, Variante 1
Bearbeiten
Die kovarianten Basisvektoren sind mit
r
→
=
(
x
,
y
,
z
)
⊤
{\displaystyle {\vec {r}}=(x,y,z)^{\top }}
:
g
→
η
:=
∂
r
→
∂
η
=
(
x
tanh
(
η
)
y
tanh
(
η
)
z
coth
(
η
)
)
,
g
→
θ
:=
∂
r
→
∂
θ
=
(
x
cot
(
θ
)
y
cot
(
θ
)
−
z
tan
(
θ
)
)
,
g
→
ψ
:=
∂
r
→
∂
ψ
=
(
−
y
x
0
)
{\displaystyle {\vec {g}}_{\eta }:={\frac {\partial {\vec {r}}}{\partial \eta }}={\begin{pmatrix}x\tanh(\eta )\\y\tanh(\eta )\\z\coth(\eta )\end{pmatrix}},\;{\vec {g}}_{\theta }:={\frac {\partial {\vec {r}}}{\partial \theta }}={\begin{pmatrix}x\cot(\theta )\\y\cot(\theta )\\-z\tan(\theta )\end{pmatrix}},\;{\vec {g}}_{\psi }:={\frac {\partial {\vec {r}}}{\partial \psi }}={\begin{pmatrix}-y\\x\\0\end{pmatrix}}}
worin cot und coth die Kehrwerte des Tangens bzw. Tangens Hyperbolikus sind. Die Basisvektoren sind senkrecht zueinander und bilden in dieser Reihenfolge im gesamten Wertebereich ein Rechtssystem . Die metrischen Faktoren sind die Beträge der kovarianten Basisvektoren und lauten:[2] :31
h
:=
h
η
=
h
θ
=
a
cosh
(
η
)
2
−
sin
(
θ
)
2
,
h
ψ
=
a
cosh
(
η
)
sin
(
θ
)
{\displaystyle h:=h_{\eta }=h_{\theta }=a{\sqrt {\cosh(\eta )^{2}-\sin(\theta )^{2}}},\,h_{\psi }=a\cosh(\eta )\sin(\theta )}
Das gestreckt sphäroidische Orthonormalsystem ist dementsprechend
c
^
η
=
1
h
(
x
tanh
(
η
)
y
tanh
(
η
)
z
coth
(
η
)
)
,
c
^
θ
=
1
h
(
x
cot
(
θ
)
y
cot
(
θ
)
−
z
tan
(
θ
)
)
,
c
^
ψ
=
(
−
sin
(
ψ
)
cos
(
ψ
)
0
)
{\displaystyle {\begin{aligned}{\hat {c}}_{\eta }={\frac {1}{h}}{\begin{pmatrix}x\tanh(\eta )\\y\tanh(\eta )\\z\coth(\eta )\end{pmatrix}},\;{\hat {c}}_{\theta }={\frac {1}{h}}{\begin{pmatrix}x\cot(\theta )\\y\cot(\theta )\\-z\tan(\theta )\end{pmatrix}},\;{\hat {c}}_{\psi }={\begin{pmatrix}-\sin(\psi )\\\cos(\psi )\\0\end{pmatrix}}\end{aligned}}}
Das Linien- und Volumenelement ergibt sich zu[2] :31
d
r
→
=
g
→
η
d
η
+
g
→
θ
d
θ
+
g
→
ψ
d
ψ
d
s
2
:=
|
d
r
→
|
2
=
a
2
[
cosh
(
η
)
2
−
sin
(
θ
)
2
]
(
d
η
2
+
d
θ
2
)
+
a
2
cosh
(
η
)
2
sin
(
θ
)
2
d
ψ
2
d
A
=
a
2
{
[
cosh
(
η
)
2
−
sin
(
θ
)
2
]
c
^
ψ
d
η
d
θ
+
cosh
(
η
)
2
−
sin
(
θ
)
2
cosh
(
η
)
sin
(
θ
)
(
c
^
η
d
θ
+
c
^
θ
d
η
)
d
ψ
}
d
V
=
a
3
[
cosh
(
η
)
2
−
sin
(
θ
)
2
]
cosh
(
η
)
sin
(
θ
)
d
η
d
θ
d
ψ
{\displaystyle {\begin{aligned}{\rm {d}}{\vec {r}}=&{\vec {g}}_{\eta }{\rm {d}}\eta +{\vec {g}}_{\theta }{\rm {d}}\theta +{\vec {g}}_{\psi }{\rm {d}}\psi \\{\rm {d}}s^{2}:=&|{\rm {d}}{\vec {r}}|^{2}=a^{2}[\cosh(\eta )^{2}-\sin(\theta )^{2}]\,({\rm {d}}\eta ^{2}+{\rm {d}}\theta ^{2})+a^{2}\cosh(\eta )^{2}\sin(\theta )^{2}\,{\rm {d}}\psi ^{2}\\{\rm {d}}A=&a^{2}\left\{[\cosh(\eta )^{2}-\sin(\theta )^{2}]{\hat {c}}_{\psi }\,{\rm {d}}\eta \,{\rm {d}}\theta +{\sqrt {\cosh(\eta )^{2}-\sin(\theta )^{2}}}\cosh(\eta )\sin(\theta )({\hat {c}}_{\eta }\,{\rm {d}}\theta +{\hat {c}}_{\theta }\,{\rm {d}}\eta )\,{\rm {d}}\psi \right\}\\{\rm {d}}V=&a^{3}[\cosh(\eta )^{2}-\sin(\theta )^{2}]\cosh(\eta )\sin(\theta )\,{\rm {d}}\eta \,{\rm {d}}\theta \,{\rm {d}}\psi \end{aligned}}}
Differentialoperatoren in abgeplatteten Sphäroidkoordinaten, Variante 1
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:[2] :32
Δ
f
=
1
a
2
[
cosh
(
η
)
2
−
sin
(
θ
)
2
]
(
∂
2
f
∂
η
2
+
tanh
(
η
)
∂
f
∂
η
+
∂
2
f
∂
θ
2
+
cot
(
θ
)
∂
f
∂
θ
)
+
…
⋯
+
1
a
2
cosh
(
η
)
2
sin
(
θ
)
2
∂
2
f
∂
ψ
2
{\displaystyle {\begin{aligned}\Delta f=&{\frac {1}{a^{2}[\cosh(\eta )^{2}-\sin(\theta )^{2}]}}\left({\frac {\partial ^{2}f}{\partial \eta ^{2}}}+\tanh(\eta ){\frac {\partial f}{\partial \eta }}+{\frac {\partial ^{2}f}{\partial \theta ^{2}}}+\cot(\theta ){\frac {\partial f}{\partial \theta }}\right)+\dots \\&\dots +{\frac {1}{a^{2}\cosh(\eta )^{2}\sin(\theta )^{2}}}{\frac {\partial ^{2}f}{\partial \psi ^{2}}}\end{aligned}}}
worin cot der Kehrwert des Tangens ist.
Lösung der Laplace- und Helmholtz-Gleichung in abgeplatteten Sphäroidkoordinaten, Variante 1
Bearbeiten
In Ellipsoid-Koordinaten gelingt immer eine Trennung der Variablen in der Laplace- und Helmholtz-Gleichung . Mit dem Separationsansatz
ϕ
(
η
,
θ
,
ψ
)
=
H
(
η
)
⋅
Θ
(
θ
)
⋅
Ψ
(
ψ
)
{\displaystyle \phi (\eta ,\theta ,\psi )=H(\eta )\cdot \Theta (\theta )\cdot \Psi (\psi )}
lautet die Helmholtz-Gleichung
Δ
ϕ
+
κ
2
ϕ
=
0
{\displaystyle \Delta \phi +\kappa ^{2}\phi =0}
:
1
a
2
[
cosh
(
η
)
2
−
sin
(
θ
)
2
]
(
∂
2
H
∂
η
2
Θ
Ψ
+
tanh
(
η
)
∂
H
∂
η
Θ
Ψ
+
H
∂
2
Θ
∂
θ
2
Ψ
+
cot
(
θ
)
H
∂
Θ
∂
θ
Ψ
)
+
…
⋯
+
1
a
2
cosh
(
η
)
2
sin
(
θ
)
2
H
Θ
∂
2
Ψ
∂
ψ
2
+
κ
2
H
Θ
Ψ
=
0
{\displaystyle {\begin{aligned}{\frac {1}{a^{2}[\cosh(\eta )^{2}-\sin(\theta )^{2}]}}\left({\frac {\partial ^{2}H}{\partial \eta ^{2}}}\Theta \Psi +\tanh(\eta ){\frac {\partial H}{\partial \eta }}\Theta \Psi +H{\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}\Psi +\cot(\theta )H{\frac {\partial \Theta }{\partial \theta }}\Psi \right)+\dots &\\\dots +{\frac {1}{a^{2}\cosh(\eta )^{2}\sin(\theta )^{2}}}H\Theta {\frac {\partial ^{2}\Psi }{\partial \psi ^{2}}}+\kappa ^{2}H\Theta \Psi &=0\end{aligned}}}
Die drei Faktoren in der separierten Lösungsfunktion bestimmen sich aus den unabhängigen gewöhnlichen Differenzialgleichungen[2] :33
∂
2
H
∂
η
2
+
tanh
(
η
)
∂
H
∂
η
+
[
κ
2
a
2
cosh
(
η
)
2
−
α
2
+
α
3
cosh
(
η
)
2
]
H
=
0
∂
2
Θ
∂
θ
2
+
cot
(
θ
)
∂
Θ
∂
θ
+
[
−
κ
2
a
2
sin
(
θ
)
2
+
α
2
−
α
3
sin
(
θ
)
2
]
Θ
=
0
∂
2
Ψ
∂
ψ
2
+
α
3
Ψ
=
0
{\displaystyle {\begin{aligned}{\frac {\partial ^{2}H}{\partial \eta ^{2}}}+\tanh(\eta ){\frac {\partial H}{\partial \eta }}+\left[\kappa ^{2}a^{2}\cosh(\eta )^{2}-\alpha _{2}+{\frac {\alpha _{3}}{\cosh(\eta )^{2}}}\right]H=&0\\{\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}+\cot(\theta ){\frac {\partial \Theta }{\partial \theta }}+\left[-\kappa ^{2}a^{2}\sin(\theta )^{2}+\alpha _{2}-{\frac {\alpha _{3}}{\sin(\theta )^{2}}}\right]\Theta =&0\\{\frac {\partial ^{2}\Psi }{\partial \psi ^{2}}}+\alpha _{3}\Psi =&0\end{aligned}}}
und den Randbedingungen . Bei der Laplace-Gleichung ist
κ
=
0
{\displaystyle \kappa =0}
.
Denn Multiplikation der Helmholtz-Gleichung mit
a
2
cosh
(
η
)
2
sin
(
θ
)
2
H
Θ
Ψ
{\displaystyle {\tfrac {a^{2}\cosh(\eta )^{2}\sin(\theta )^{2}}{H\Theta \Psi }}}
liefert:
cosh
(
η
)
2
sin
(
θ
)
2
cosh
(
η
)
2
−
sin
(
θ
)
2
(
∂
2
H
∂
η
2
H
+
tanh
(
η
)
∂
H
∂
η
H
+
∂
2
Θ
∂
θ
2
Θ
+
∂
Θ
∂
θ
Θ
tan
(
θ
)
)
+
κ
2
a
2
cosh
(
η
)
2
sin
(
θ
)
2
+
∂
2
Ψ
∂
ψ
2
Ψ
=
0
{\displaystyle {\begin{aligned}{\frac {\cosh(\eta )^{2}\sin(\theta )^{2}}{\cosh(\eta )^{2}-\sin(\theta )^{2}}}\left({\frac {\frac {\partial ^{2}H}{\partial \eta ^{2}}}{H}}+\tanh(\eta ){\frac {\frac {\partial H}{\partial \eta }}{H}}+{\frac {\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}{\Theta }}+{\frac {\frac {\partial \Theta }{\partial \theta }}{\Theta \tan(\theta )}}\right)+\kappa ^{2}a^{2}\cosh(\eta )^{2}\sin(\theta )^{2}+{\frac {\frac {\partial ^{2}\Psi }{\partial \psi ^{2}}}{\Psi }}=0\end{aligned}}}
Weil auf der rechten Seite eine Konstante (null) steht und nur der letzte Term auf der linken Seite von ψ abhängt, ist letzterer ebenfalls konstant:
∂
2
Ψ
∂
ψ
2
Ψ
=
−
α
3
{\displaystyle {\frac {\frac {\partial ^{2}\Psi }{\partial \psi ^{2}}}{\Psi }}=-\alpha _{3}}
Einsetzen dieser Konstanten erlaubt auch η und θ zu trennen:
∂
2
H
∂
η
2
H
+
tanh
(
η
)
∂
H
∂
η
H
+
κ
2
a
2
cosh
(
η
)
2
+
α
3
cosh
(
η
)
2
=
…
⋯
=
−
[
∂
2
Θ
∂
θ
2
Θ
+
∂
Θ
∂
θ
Θ
tan
(
θ
)
−
κ
2
a
2
sin
(
θ
)
2
−
α
3
sin
(
θ
)
2
]
{\displaystyle {\begin{aligned}&{\frac {\frac {\partial ^{2}H}{\partial \eta ^{2}}}{H}}+\tanh(\eta ){\frac {\frac {\partial H}{\partial \eta }}{H}}+\kappa ^{2}a^{2}\cosh(\eta )^{2}+{\frac {\alpha _{3}}{\cosh(\eta )^{2}}}=\dots \\&\qquad \qquad \qquad \dots =-\left[{\frac {\frac {\partial ^{2}\Theta }{\partial \theta ^{2}}}{\Theta }}+{\frac {\frac {\partial \Theta }{\partial \theta }}{\Theta \tan(\theta )}}-\kappa ^{2}a^{2}\sin(\theta )^{2}-{\frac {\alpha _{3}}{\sin(\theta )^{2}}}\right]\end{aligned}}}
Weil die linke Seite nur von η abhängt und die rechte Seite nur von θ , sind beide Seiten gleich einer Konstanten α 2 . So entstehen die oben angegebene Differenzialgleichungen[2] :30 Ein gleichbedeutendes Ergebnis wird mit dem im Hauptartikel beschriebenen Verfahren und der Stäckel-Matrix[2] :31
S
=
(
a
2
cosh
(
η
)
2
−
1
1
cosh
(
η
)
2
−
a
2
sin
(
θ
)
2
1
−
1
sin
(
θ
)
2
0
0
1
)
{\displaystyle \mathbf {S} ={\begin{pmatrix}a^{2}\cosh(\eta )^{2}&-1&{\frac {1}{\cosh(\eta )^{2}}}\\-a^{2}\sin(\theta )^{2}&1&{\frac {-1}{\sin(\theta )^{2}}}\\0&0&1\end{pmatrix}}}
erzielt.
Abgeplattete Sphäroidkoordinaten, Variante 2
Bearbeiten
Die abgeplatteten Sphäroidkoordinaten der Variante 2 benutzen nicht die Variablen η,θ,ψ der ersten Variante, sondern deren Funktionswerte:[3] :662
ξ
1
=
a
sinh
(
η
)
,
ξ
2
=
cos
(
θ
)
,
ξ
3
=
cos
(
ψ
)
{\displaystyle \xi _{1}=a\sinh(\eta ),\;\xi _{2}=\cos(\theta ),\;\xi _{3}=\cos(\psi )}
#Gestreckte Sphäroidkoordinaten, Variante 2 , benutzen im Vergleich hierzu cosh statt sinh. Die kartesischen Koordinaten
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
berechnen sich aus den abgeplatteten Sphäroidkoordinaten
ξ
1
,
2
,
3
∈
R
,
ξ
1
≥
0
,
−
1
≤
ξ
2
,
3
≤
1
{\displaystyle \xi _{1,2,3}\in \mathbb {R} ,\,\xi _{1}\geq 0,-1\leq \xi _{2,3}\leq 1}
gemäß:
r
→
:=
(
x
y
z
)
=
(
ξ
3
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
(
1
−
ξ
3
2
)
ξ
1
ξ
2
)
{\displaystyle {\vec {r}}:={\begin{pmatrix}x\\y\\z\end{pmatrix}}={\begin{pmatrix}\xi _{3}{\sqrt {(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})}}\\{\sqrt {(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})(1-\xi _{3}^{2})}}\\\xi _{1}\xi _{2}\end{pmatrix}}}
Koordinatenflächen in abgeplatteten Sphäroidkoordinaten, Variante 2
Bearbeiten
Koordinatenflächen in abgeplatteten Sphäroidkoordinaten.
In der Variante 2 ist auf dem Ellipsoid (ξ 1 =const.,rot)
x
2
+
y
2
ξ
1
2
+
a
2
+
z
2
ξ
1
2
=
1
{\displaystyle {\frac {x^{2}+y^{2}}{\xi _{1}^{2}+a^{2}}}+{\frac {z^{2}}{\xi _{1}^{2}}}=1}
auf dem einschaligen Rotationshyperboloid (ξ 2 =const., blau)
x
2
+
y
2
a
2
(
1
−
ξ
2
2
)
−
z
2
a
2
ξ
2
2
=
1
{\displaystyle {\frac {x^{2}+y^{2}}{a^{2}(1-\xi _{2}^{2})}}-{\frac {z^{2}}{a^{2}\xi _{2}^{2}}}=1}
und in der Halbebene (ξ 3 =const., gelb)
y
x
=
1
−
ξ
3
2
ξ
3
{\displaystyle {\frac {y}{x}}={\frac {\sqrt {1-\xi _{3}^{2}}}{\xi _{3}}}}
Hieraus ergibt sich andererseits
ξ
1
2
=
(
x
2
+
y
2
+
z
2
−
a
2
)
2
+
4
a
2
z
2
−
a
2
+
x
2
+
y
2
+
z
2
2
ξ
2
2
=
(
x
2
+
y
2
+
z
2
−
a
2
)
2
+
4
a
2
z
2
+
a
2
−
x
2
−
y
2
−
z
2
2
a
2
ξ
3
2
=
x
2
x
2
+
y
2
{\displaystyle {\begin{aligned}\xi _{1}^{2}=&{\frac {{\sqrt {(x^{2}+y^{2}+z^{2}-a^{2})^{2}+4a^{2}z^{2}}}-a^{2}+x^{2}+y^{2}+z^{2}}{2}}\\\xi _{2}^{2}=&{\frac {{\sqrt {(x^{2}+y^{2}+z^{2}-a^{2})^{2}+4a^{2}z^{2}}}+a^{2}-x^{2}-y^{2}-z^{2}}{2a^{2}}}\\\xi _{3}^{2}=&{\frac {x^{2}}{x^{2}+y^{2}}}\end{aligned}}}
Metrische Faktoren, Weg- und Flächen- und Volumenelemente in abgeplatteten Sphäroidkoordinaten, Variante 2
Bearbeiten
Die kovarianten Basisvektoren sind mit
r
→
=
(
x
,
y
,
z
)
⊤
{\displaystyle {\vec {r}}=(x,y,z)^{\top }}
:
g
→
1
:=
∂
r
→
∂
ξ
1
=
(
ξ
1
x
ξ
1
2
+
a
2
ξ
1
y
ξ
1
2
+
a
2
z
ξ
1
)
,
g
→
2
:=
∂
r
→
∂
ξ
2
=
(
−
ξ
2
x
1
−
ξ
2
2
−
ξ
2
y
1
−
ξ
2
2
z
ξ
2
)
,
g
→
3
:=
∂
r
→
∂
ξ
3
=
(
x
ξ
3
−
ξ
3
y
1
−
ξ
3
2
0
)
{\displaystyle {\vec {g}}_{1}:={\frac {\partial {\vec {r}}}{\partial \xi _{1}}}={\begin{pmatrix}{\frac {\xi _{1}x}{\xi _{1}^{2}+a^{2}}}\\{\frac {\xi _{1}y}{\xi _{1}^{2}+a^{2}}}\\{\frac {z}{\xi _{1}}}\end{pmatrix}},\;{\vec {g}}_{2}:={\frac {\partial {\vec {r}}}{\partial \xi _{2}}}={\begin{pmatrix}-{\frac {\xi _{2}x}{1-\xi _{2}^{2}}}\\-{\frac {\xi _{2}y}{1-\xi _{2}^{2}}}\\{\frac {z}{\xi _{2}}}\end{pmatrix}},\;{\vec {g}}_{3}:={\frac {\partial {\vec {r}}}{\partial \xi _{3}}}={\begin{pmatrix}{\frac {x}{\xi _{3}}}\\-{\frac {\xi _{3}y}{1-\xi _{3}^{2}}}\\0\end{pmatrix}}}
Die Basisvektoren sind senkrecht zueinander und bilden in dieser Reihenfolge ein Rechtssystem . Die metrischen Faktoren sind die Beträge der kovarianten Basisvektoren und lauten:
h
1
=
ξ
1
2
+
a
2
ξ
2
2
ξ
1
2
+
a
2
,
h
2
=
ξ
1
2
+
a
2
ξ
2
2
1
−
ξ
2
2
,
h
3
=
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
1
−
ξ
3
2
{\displaystyle h_{1}={\sqrt {\frac {\xi _{1}^{2}+a^{2}\xi _{2}^{2}}{\xi _{1}^{2}+a^{2}}}},\;h_{2}={\sqrt {\frac {\xi _{1}^{2}+a^{2}\xi _{2}^{2}}{1-\xi _{2}^{2}}}},\;h_{3}={\sqrt {\frac {(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})}{1-\xi _{3}^{2}}}}}
Das gestreckt sphäroidische Orthonormalsystem schreibt sich in dieser Formulierung
c
^
1
=
1
ξ
1
2
+
a
2
ξ
2
2
(
ξ
1
ξ
3
1
−
ξ
2
2
ξ
1
(
1
−
ξ
2
2
)
(
1
−
ξ
3
2
)
ξ
2
ξ
1
2
+
a
2
)
c
^
2
=
1
ξ
1
2
+
a
2
ξ
2
2
(
−
ξ
2
ξ
3
ξ
1
2
+
a
2
−
ξ
2
(
ξ
1
2
+
a
2
)
(
1
−
ξ
3
2
)
ξ
1
1
−
ξ
2
2
)
,
c
^
3
=
(
1
−
ξ
3
2
−
ξ
3
0
)
{\displaystyle {\begin{aligned}{\hat {c}}_{1}=&{\frac {1}{\sqrt {\xi _{1}^{2}+a^{2}\xi _{2}^{2}}}}{\begin{pmatrix}\xi _{1}\xi _{3}{\sqrt {1-\xi _{2}^{2}}}\\\xi _{1}{\sqrt {(1-\xi _{2}^{2})(1-\xi _{3}^{2})}}\\\xi _{2}{\sqrt {\xi _{1}^{2}+a^{2}}}\end{pmatrix}}\\{\hat {c}}_{2}=&{\frac {1}{\sqrt {\xi _{1}^{2}+a^{2}\xi _{2}^{2}}}}{\begin{pmatrix}-\xi _{2}\xi _{3}{\sqrt {\xi _{1}^{2}+a^{2}}}\\-\xi _{2}{\sqrt {(\xi _{1}^{2}+a^{2})(1-\xi _{3}^{2})}}\\\xi _{1}{\sqrt {1-\xi _{2}^{2}}}\end{pmatrix}},\;{\hat {c}}_{3}={\begin{pmatrix}{\sqrt {1-\xi _{3}^{2}}}\\-\xi _{3}\\0\end{pmatrix}}\end{aligned}}}
Das Linien- und Volumenelement ergibt sich zu
d
r
→
=
g
→
1
d
ξ
1
+
g
→
2
d
ξ
2
+
g
→
3
d
ξ
3
d
s
2
:=
|
d
r
→
|
2
=
ξ
1
2
+
a
2
ξ
2
2
ξ
1
2
+
a
2
d
ξ
1
2
+
ξ
1
2
+
a
2
ξ
2
2
1
−
ξ
2
2
d
ξ
2
2
+
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
1
−
ξ
3
2
d
ξ
3
2
d
V
=
ξ
1
2
+
a
2
ξ
2
2
1
−
ξ
3
2
d
ξ
1
d
ξ
2
d
ξ
3
{\displaystyle {\begin{aligned}{\rm {d}}{\vec {r}}=&{\vec {g}}_{1}{\rm {d}}\xi _{1}+{\vec {g}}_{2}{\rm {d}}\xi _{2}+{\vec {g}}_{3}{\rm {d}}\xi _{3}\\{\rm {d}}s^{2}:=&|{\rm {d}}{\vec {r}}|^{2}={\frac {\xi _{1}^{2}+a^{2}\xi _{2}^{2}}{\xi _{1}^{2}+a^{2}}}\,{\rm {d}}\xi _{1}^{2}+{\frac {\xi _{1}^{2}+a^{2}\xi _{2}^{2}}{1-\xi _{2}^{2}}}\,{\rm {d}}\xi _{2}^{2}+{\frac {(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})}{1-\xi _{3}^{2}}}\,{\rm {d}}\xi _{3}^{2}\\{\rm {d}}V=&{\frac {\xi _{1}^{2}+a^{2}\xi _{2}^{2}}{\sqrt {1-\xi _{3}^{2}}}}\,{\rm {d}}\xi _{1}\,{\rm {d}}\xi _{2}\,{\rm {d}}\xi _{3}\end{aligned}}}
Differentialoperatoren in abgeplatteten Sphäroidkoordinaten, Variante 2
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:
Δ
f
=
(
ξ
1
2
+
a
2
)
∂
2
f
∂
ξ
1
2
+
2
ξ
1
∂
f
∂
ξ
1
+
(
1
−
ξ
2
2
)
∂
2
f
∂
ξ
2
2
−
2
ξ
2
∂
f
∂
ξ
2
ξ
1
2
+
a
2
ξ
2
2
+
(
1
−
ξ
3
2
)
∂
2
f
∂
ξ
3
2
−
ξ
3
∂
f
∂
ξ
3
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
{\displaystyle \Delta f={\frac {(\xi _{1}^{2}+a^{2}){\frac {\partial ^{2}f}{\partial \xi _{1}^{2}}}+2\xi _{1}{\frac {\partial f}{\partial \xi _{1}}}+(1-\xi _{2}^{2}){\frac {\partial ^{2}f}{\partial \xi _{2}^{2}}}-2\xi _{2}{\frac {\partial f}{\partial \xi _{2}}}}{\xi _{1}^{2}+a^{2}\xi _{2}^{2}}}+{\frac {(1-\xi _{3}^{2}){\frac {\partial ^{2}f}{\partial \xi _{3}^{2}}}-\xi _{3}{\frac {\partial f}{\partial \xi _{3}}}}{(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})}}}
Lösung der Laplace- und Helmholtz-Gleichung in abgeplatteten Sphäroidkoordinaten, Variante 2
Bearbeiten
In Ellipsoid-Koordinaten gelingt immer eine Trennung der Variablen in der Laplace- und Helmholtz-Gleichung . Mit dem Separationsansatz
ϕ
(
ξ
1
,
ξ
2
,
ξ
3
)
=
X
(
ξ
1
)
⋅
Y
(
ξ
2
)
⋅
Z
(
ξ
3
)
{\displaystyle \phi (\xi _{1},\xi _{2},\xi _{3})=X(\xi _{1})\cdot Y(\xi _{2})\cdot Z(\xi _{3})}
lautet die Helmholtz-Gleichung
Δ
ϕ
+
κ
2
ϕ
=
0
{\displaystyle \Delta \phi +\kappa ^{2}\phi =0}
:
(
ξ
1
2
+
a
2
)
∂
2
ϕ
∂
ξ
1
2
+
2
ξ
1
∂
ϕ
∂
ξ
1
+
(
1
−
ξ
2
2
)
∂
2
ϕ
∂
ξ
2
2
−
2
ξ
2
∂
ϕ
∂
ξ
2
ξ
1
2
+
a
2
ξ
2
2
+
(
1
−
ξ
3
2
)
∂
2
ϕ
∂
ξ
3
2
−
ξ
3
∂
ϕ
∂
ξ
3
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
+
κ
2
ϕ
=
(
ξ
1
2
+
a
2
)
∂
2
X
∂
ξ
1
2
Y
Z
+
2
ξ
1
∂
X
∂
ξ
1
Y
Z
+
(
1
−
ξ
2
2
)
X
∂
2
Y
∂
ξ
2
2
Z
−
2
ξ
2
X
∂
Y
∂
ξ
2
Z
ξ
1
2
+
a
2
ξ
2
2
+
…
⋯
+
(
1
−
ξ
3
2
)
X
Y
∂
2
Z
∂
ξ
3
2
−
ξ
3
X
Y
∂
Z
∂
ξ
3
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
+
κ
2
X
Y
Z
=
0
{\displaystyle {\begin{aligned}{\frac {(\xi _{1}^{2}+a^{2}){\frac {\partial ^{2}\phi }{\partial \xi _{1}^{2}}}+2\xi _{1}{\frac {\partial \phi }{\partial \xi _{1}}}+(1-\xi _{2}^{2}){\frac {\partial ^{2}\phi }{\partial \xi _{2}^{2}}}-2\xi _{2}{\frac {\partial \phi }{\partial \xi _{2}}}}{\xi _{1}^{2}+a^{2}\xi _{2}^{2}}}+{\frac {(1-\xi _{3}^{2}){\frac {\partial ^{2}\phi }{\partial \xi _{3}^{2}}}-\xi _{3}{\frac {\partial \phi }{\partial \xi _{3}}}}{(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})}}+\kappa ^{2}\phi &=\\{\frac {(\xi _{1}^{2}+a^{2}){\frac {\partial ^{2}X}{\partial \xi _{1}^{2}}}YZ+2\xi _{1}{\frac {\partial X}{\partial \xi _{1}}}YZ+(1-\xi _{2}^{2})X{\frac {\partial ^{2}Y}{\partial \xi _{2}^{2}}}Z-2\xi _{2}X{\frac {\partial Y}{\partial \xi _{2}}}Z}{\xi _{1}^{2}+a^{2}\xi _{2}^{2}}}+\dots \qquad &\\\dots +{\frac {(1-\xi _{3}^{2})XY{\frac {\partial ^{2}Z}{\partial \xi _{3}^{2}}}-\xi _{3}XY{\frac {\partial Z}{\partial \xi _{3}}}}{(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})}}+\kappa ^{2}XYZ&=0\end{aligned}}}
Die drei Faktoren in der separierten Lösungsfunktion bestimmen sich aus den unabhängigen gewöhnlichen Differenzialgleichungen
2
ξ
1
∂
X
∂
ξ
1
+
(
ξ
1
2
+
a
2
)
∂
2
X
∂
ξ
1
2
+
(
κ
2
ξ
1
2
−
α
2
−
α
3
a
2
ξ
1
2
+
a
2
)
X
=
0
2
ξ
2
∂
Y
∂
ξ
2
−
(
1
−
ξ
2
2
)
∂
2
Y
∂
ξ
2
2
−
(
κ
2
a
2
ξ
2
2
+
α
2
+
α
3
1
−
ξ
2
2
)
Y
=
0
(
1
−
ξ
3
2
)
∂
2
Z
∂
ξ
3
2
−
ξ
3
∂
Z
∂
ξ
3
−
α
3
Z
=
0
{\displaystyle {\begin{aligned}2\xi _{1}{\frac {\partial X}{\partial \xi _{1}}}+(\xi _{1}^{2}+a^{2}){\frac {\partial ^{2}X}{\partial \xi _{1}^{2}}}+\left(\kappa ^{2}\xi _{1}^{2}-\alpha _{2}-{\frac {\alpha _{3}a^{2}}{\xi _{1}^{2}+a^{2}}}\right)X=&0\\2\xi _{2}{\frac {\partial Y}{\partial \xi _{2}}}-(1-\xi _{2}^{2}){\frac {\partial ^{2}Y}{\partial \xi _{2}^{2}}}-\left(\kappa ^{2}a^{2}\xi _{2}^{2}+\alpha _{2}+{\frac {\alpha _{3}}{1-\xi _{2}^{2}}}\right)Y=&0\\(1-\xi _{3}^{2}){\frac {\partial ^{2}Z}{\partial \xi _{3}^{2}}}-\xi _{3}{\frac {\partial Z}{\partial \xi _{3}}}-\alpha _{3}Z=&0\end{aligned}}}
und den Randbedingungen . Bei der Laplace-Gleichung ist
κ
=
0
{\displaystyle \kappa =0}
.
Denn Multiplikation der Helmholtz-Gleichung mit
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
X
Y
Z
{\displaystyle {\tfrac {(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})}{XYZ}}}
liefert:
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
ξ
1
2
+
a
2
ξ
2
2
[
(
ξ
1
2
+
a
2
)
∂
2
X
∂
ξ
1
2
X
+
2
ξ
1
∂
X
∂
ξ
1
X
+
(
1
−
ξ
2
2
)
∂
2
Y
∂
ξ
2
2
Y
−
2
ξ
2
∂
Y
∂
ξ
2
Y
]
+
…
⋯
+
{
(
1
−
ξ
3
2
)
∂
2
Z
∂
ξ
3
2
Z
−
ξ
3
∂
Z
∂
ξ
3
Z
}
+
κ
2
(
ξ
1
2
+
a
2
)
(
1
−
ξ
2
2
)
=
0
{\displaystyle {\begin{aligned}{\frac {(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})}{\xi _{1}^{2}+a^{2}\xi _{2}^{2}}}\left[(\xi _{1}^{2}+a^{2}){\frac {\frac {\partial ^{2}X}{\partial \xi _{1}^{2}}}{X}}+2\xi _{1}{\frac {\frac {\partial X}{\partial \xi _{1}}}{X}}+(1-\xi _{2}^{2}){\frac {\frac {\partial ^{2}Y}{\partial \xi _{2}^{2}}}{Y}}-2\xi _{2}{\frac {\frac {\partial Y}{\partial \xi _{2}}}{Y}}\right]+\dots &\\\dots +\left\{(1-\xi _{3}^{2}){\frac {\frac {\partial ^{2}Z}{\partial \xi _{3}^{2}}}{Z}}-\xi _{3}{\frac {\frac {\partial Z}{\partial \xi _{3}}}{Z}}\right\}+\kappa ^{2}(\xi _{1}^{2}+a^{2})(1-\xi _{2}^{2})&=0\end{aligned}}}
Weil auf der rechten Seite eine Konstante (null) steht und nur der Term in der geschweiften Klammer von ξ 3 abhängt, ist dieser ebenfalls eine Konstante α 3 :
(
1
−
ξ
3
2
)
∂
2
Z
∂
ξ
3
2
Z
−
ξ
3
∂
Z
∂
ξ
3
Z
=
α
3
{\displaystyle (1-\xi _{3}^{2}){\frac {\frac {\partial ^{2}Z}{\partial \xi _{3}^{2}}}{Z}}-\xi _{3}{\frac {\frac {\partial Z}{\partial \xi _{3}}}{Z}}=\alpha _{3}}
Einsetzen dieser Konstanten erlaubt auch ξ 1 und ξ 2 zu trennen:
2
ξ
1
∂
X
∂
ξ
1
X
+
(
ξ
1
2
+
a
2
)
∂
2
X
∂
ξ
1
2
X
+
κ
2
ξ
1
2
−
α
3
a
2
ξ
1
2
+
a
2
=
…
⋯
=
2
ξ
2
∂
Y
∂
ξ
2
Y
−
(
1
−
ξ
2
2
)
∂
2
Y
∂
ξ
2
2
Y
−
κ
2
a
2
ξ
2
2
−
α
3
1
−
ξ
2
2
{\displaystyle {\begin{aligned}&2\xi _{1}{\frac {\frac {\partial X}{\partial \xi _{1}}}{X}}+(\xi _{1}^{2}+a^{2}){\frac {\frac {\partial ^{2}X}{\partial \xi _{1}^{2}}}{X}}+\kappa ^{2}\xi _{1}^{2}-{\frac {\alpha _{3}a^{2}}{\xi _{1}^{2}+a^{2}}}=\dots \\&\qquad \qquad \qquad \dots =2\xi _{2}{\frac {\frac {\partial Y}{\partial \xi _{2}}}{Y}}-(1-\xi _{2}^{2}){\frac {\frac {\partial ^{2}Y}{\partial \xi _{2}^{2}}}{Y}}-\kappa ^{2}a^{2}\xi _{2}^{2}-{\frac {\alpha _{3}}{1-\xi _{2}^{2}}}\end{aligned}}}
Weil die linke Seite nur von ξ 1 abhängt und die rechte Seite nur von ξ 2 , sind beide Seiten gleich einer Konstanten α 2 . Damit resultieren die oben angegebene Differenzialgleichungen. Ein gleichbedeutendes Ergebnis entsteht mit der Stäckel-Matrix[3] :662
S
=
(
1
1
ξ
1
2
+
a
2
−
a
2
(
ξ
1
2
+
a
2
)
2
−
a
2
−
1
1
−
ξ
2
2
1
(
1
−
ξ
2
2
)
2
0
0
−
1
1
−
ξ
3
2
)
{\displaystyle \mathbf {S} ={\begin{pmatrix}1&{\frac {1}{\xi _{1}^{2}+a^{2}}}&{\frac {-a^{2}}{(\xi _{1}^{2}+a^{2})^{2}}}\\-a^{2}&{\frac {-1}{1-\xi _{2}^{2}}}&{\frac {1}{(1-\xi _{2}^{2})^{2}}}\\0&0&{\frac {-1}{1-\xi _{3}^{2}}}\end{pmatrix}}}
und der im Hauptartikel beschriebenen Methode.