Krummlinige Koordinaten
Krummlinige Koordinaten
Krummlinige Koordinaten
Erfolgreiche ePaper selbst erstellen
Machen Sie aus Ihren PDF Publikationen ein blätterbares Flipbook mit unserer einzigartigen Google optimierten e-Paper Software.
<strong>Krummlinige</strong> <strong>Koordinaten</strong><br />
Einige <strong>Koordinaten</strong>systeme im R 3 haben wir bereits kennengelernt :<br />
x 1 , x 2 , x 3<br />
. . . kartesische <strong>Koordinaten</strong><br />
r, φ, x 3 . . . Zylinderkoordinaten<br />
r, φ, ϑ . . . Kugelkoordinaten<br />
Sind andere <strong>Koordinaten</strong> u 1 , u 2 , u 3 gegeben, sodass wir die kartesischen<br />
<strong>Koordinaten</strong> x 1 , x 2 , x 3 als Funktionen von u 1 , u 2 , u 3 schreiben können,<br />
i.e.<br />
x 1 = x 1 (u 1 , u 2 , u 3 ) , x 2 = x 2 (u 1 , u 2 , u 3 ) , x 3 = x 3 (u 1 , u 2 , u 3 ) ,<br />
dann können wir umgekehrt (lokal) die u i durch die x i ausdrücken, wenn<br />
die Jacobi-Determinante verschieden von Null ist, d.h.<br />
∂x 1 ∂x 2 ∂x 3<br />
∂u<br />
∣ 1 ∂u 1 ∂u 1<br />
∣<br />
∂x 1 ∂x 2 ∂x 3<br />
∂u 2<br />
∂x 1<br />
∂u 2<br />
∂x 2<br />
∂u 2<br />
∂x 3<br />
∂u 3 ∂u 3<br />
∣ ∣∣∣∣∣∣<br />
= ∂(x 1,x 2 ,x 3 )<br />
∂(u 1 ,u 2 ,u 3 ) ≠ 0<br />
∂u 3<br />
Fixieren wir einen Punkt (u 0 1, u 0 2, u 0 3) , dann erhalten wir drei Kurven<br />
durch diesen Punkt, indem wir zwei <strong>Koordinaten</strong> festhalen und die dritte<br />
als Kurvenparameter nehmen. Dies sind die <strong>Koordinaten</strong>linien durch<br />
diesen Punkt.<br />
⃗r(u 1 ) = (u 1 , u 0 2, u 0 3) , ⃗r(u 2 ) = (u 0 1, u 2 , u 0 3) , ⃗r(u 3 ) = (u 0 1, u 0 2, u 3 )<br />
Wenn sich die <strong>Koordinaten</strong>linien in jedem Punkt jeweils paarweise im<br />
rechten Winkel schneiden, spricht man von rechtwinkeligen oder auch<br />
von orthogonalen <strong>Koordinaten</strong>.<br />
Die Tangenteneinheitsvektoren in einem Punkt ⃗r 0 = (u 0 1, u 0 2, u 0 3) ergeben<br />
sich durch<br />
∣ ∣ ∣<br />
∣∣⃗r=⃗r0<br />
⃗T i = ⃗e ui | ⃗r=⃗r0<br />
= 1 ∂⃗r ∣∣⃗r=⃗r0 h ui ∂u i<br />
mit h ui = ∣ ∂⃗r ∣∣<br />
∂u i<br />
1
Die Vektoren e u1 , e u2 , e u3 hängen von u 1 , u 2 , u 3 ab und haben daher im<br />
allgemeinen eine vom Ort abhängige Richtung. Dies ist bei kartesischen<br />
<strong>Koordinaten</strong> nicht der Fall.<br />
Für orthogonale <strong>Koordinaten</strong> gilt ⃗e ui · ⃗e uj = δ ij , i, j = 1, 2, 3. Die Vektoren<br />
⃗e u1 , ⃗e u2 , ⃗e u3 bilden eine (orthogonale) Basis des R 3 , und da diese<br />
Vektoren auch Einheitsvektoren sind, spricht man von einer Orthonormalbasis.<br />
Des weiteren spricht man von einem Rechtssystem (oder einer rechtshändigen<br />
Basis), wenn gilt : ⃗e u1 × ⃗e u2 = ⃗e u3 , ⃗e u2 × ⃗e u3 = ⃗e u1 , ⃗e u3 × ⃗e u1 = ⃗e u2<br />
Beispiel. (Zylinderkoordinaten)<br />
⎛ ⎞ ⎛ ⎞<br />
x 1 (ρ, φ) ρ cos φ<br />
⃗r = ⎝ x 2 (ρ, φ) ⎠ = ⎝ ρ sin φ ⎠<br />
x 3<br />
x 3<br />
mit<br />
0 ≤ ρ < ∞ , 0 ≤ φ < 2π , −∞ < x 3 < +∞<br />
Die Umkehrung ist durch ρ = √ ) (x<br />
x 2 1 + x2 2 , φ = arctan 2<br />
x 1<br />
gegeben.<br />
Die <strong>Koordinaten</strong>linien sind<br />
⎛ ⎞<br />
ρ cos φ 0<br />
• ⃗r 1 (ρ) = ⎝ ρ sin φ 0<br />
⎠ . . . eine von der x 3 -Achse in der Höhe x 0<br />
x 0 3<br />
3<br />
Winkel φ 0 ausgehende Halbgerade<br />
⎛ ⎞<br />
ρ 0 cos φ<br />
• ⃗r 2 (φ) = ⎝ ρ 0 sin φ ⎠ . . . ein Kreis mit Radius ρ 0 in der Höhe x 0<br />
x 0 3<br />
3<br />
⎛<br />
• ⃗r 3 (x 3 ) = ⎝<br />
ρ 0 , φ 0<br />
ρ 0 cos φ 0<br />
ρ 0 sin φ 0<br />
x 0 3<br />
⎞<br />
im<br />
⎠ . . . eine Gerade von x 3 = −∞ bis +∞ bei<br />
Die normierten Basisvektoren in einem Punkt (ρ, φ, x 3 ) sind<br />
2
⎛<br />
⃗e ρ = 1 ∂⃗r<br />
h ρ ∂ρ = ⎝<br />
⎛<br />
⃗e φ = 1 ∂⃗r<br />
h φ ∂φ = ⎝<br />
⃗e x3 = 1<br />
h x3<br />
∂⃗r<br />
∂x 3<br />
=<br />
cos φ<br />
sin φ<br />
0<br />
⎞<br />
− sin φ<br />
cos φ<br />
0<br />
⎛<br />
⎝<br />
0<br />
0<br />
1<br />
⎞<br />
⎠ (h ρ = 1)<br />
⎞<br />
⎠ (h φ = ρ)<br />
⎠ (h x3 = 1)<br />
Die Zylinderkoordinaten sind damit offenbar orthogonale <strong>Koordinaten</strong>, und<br />
die Vektoren ⃗e ρ , ⃗e φ , ⃗e x3 bilden ein Rechtssystem.<br />
Beispiel. (Kugelkoordinaten)<br />
⎛ ⎞ ⎛<br />
⎞<br />
x 1 (r, ϑ, φ) r sin ϑ cos φ<br />
⃗r = ⎝ x 2 (r, ϑ, φ) ⎠ = ⎝ r sin ϑ sin φ ⎠<br />
x 3 (r, ϑ, φ) r cos ϑ<br />
0 ≤ r < ∞ , 0 ≤ ϑ ≤ π , 0 ≤ φ < 2π<br />
mit<br />
Die Umkehrung ist durch<br />
r = √ x 2 1 + x2 2 + x2 3 , ϑ = arctan √<br />
x<br />
2<br />
1 +x 2 2<br />
x 2 3<br />
Die <strong>Koordinaten</strong>linien sind<br />
, φ = arctan<br />
( )<br />
x 2<br />
x 1<br />
gegeben.<br />
• ⃗r 1 (r) = ⃗r(r, ϑ 0 , φ 0 ) . . . eine Halbgerade vom Ursprung in Richtung<br />
ϑ 0 , φ 0 .<br />
• ⃗r 2 (r) = ⃗r(r 0 , ϑ, φ 0 ) . . . ein Längenkreis mit Radius r 0<br />
durch die x 3 -Achse im Winkel φ 0 zur x 1 x 3 -Ebene.<br />
in der Ebene<br />
• ⃗r 3 (r) = ⃗r(r 0 , ϑ 0 , φ) . . . ein Breitenkreis mit Radius r 0 sin ϑ 0<br />
zur x 1 x 2 -Ebene.<br />
parallel<br />
Die normierten Basisvektoren in einem Punkt (r, ϑ, φ) sind<br />
3
⎛<br />
⃗e r = 1 ∂⃗r<br />
h r ∂r = ⎝<br />
⎛<br />
⃗e ϑ = 1 ∂⃗r<br />
h ϑ ∂ϑ = ⎝<br />
⎛<br />
⃗e φ = 1 ∂⃗r<br />
h φ ∂φ = ⎝<br />
sin ϑ cos φ<br />
sin ϑ sin φ<br />
cos ϑ<br />
cos ϑ cos φ<br />
cos ϑ sin φ<br />
− sin ϑ<br />
− sin φ<br />
cos φ<br />
0<br />
⎞<br />
⎞<br />
⎠ (h r = 1)<br />
⎞<br />
⎠ (h ϑ = r)<br />
⎠ (h φ = r sin ϑ)<br />
Die Kugelkoordinaten sind damit ebenfalls orthogonale <strong>Koordinaten</strong>, und<br />
die Vektoren ⃗e r , ⃗e ϑ , ⃗e φ bilden ein Rechtssystem.<br />
Sobald wir eine Basis in einem Vektorraum haben, können wir die Vektoren<br />
des Vektorraums als Linearkombination der Basisvektoren darstellen.<br />
Im Falle einer orthonormalen Basis (hier im R 3 ) gilt offenbar<br />
⃗F = F u1 ⃗e u1 + F u2 ⃗e u2 + F u3 ⃗e u3 ⇒ F ui = ⃗ F · ⃗e ui = 1<br />
h ui<br />
⃗ F · ∂⃗r<br />
∂u i<br />
Beispiel. Der Ortsvektor ⃗r ist in der kartesischen Basis durch ⃗r =<br />
x 1 ⃗e 1 + x 2 ⃗e 2 + x 3 ⃗e 3 dargestellt, hat also die Komponenten x 1 , x 2 , x 3 .<br />
In der Basis der Zylinderkoordinaten ist ⃗r = r ρ ⃗e ρ + r φ ⃗e φ + r x3 ⃗e x3<br />
r ρ = ⃗r · ⃗e ρ = (ρ cos φ, ρ sin φ, x 3 ) · (cos φ, sin φ, 0) = ρ<br />
r φ = ⃗r · ⃗e φ = (ρ cos φ, ρ sin φ, x 3 ) · (− sin φ, cos φ, 0) = 0<br />
mit<br />
r x3 = ⃗r · ⃗e x3 = (ρ cos φ, ρ sin φ, x 3 ) · (0, 0, 1) = x 3<br />
Beispiel. In Kugelkoordinaten ist ⃗r = (r sin ϑ cos φ, r sin ϑ sin φ, r cos ϑ)<br />
und damit offenbar (siehe vorher) ⃗r = r⃗e r .<br />
(Die Koeffizienten von ⃗e ϑ und ⃗e φ sind Null).<br />
Beispiel. Wir bestimmen die Komponenten des Vektors ⃗ F = (x3 , 0, 0)<br />
4
im Basissystem der Kugelkoordinaten.<br />
F r = F ⃗ · ⃗e r = (r cos ϑ, 0, 0) · (sin ϑ cos φ, sin ϑ sin φ, cos ϑ) =<br />
= r sin ϑ cos ϑ cos φ<br />
F ϑ = F ⃗ · ⃗e ϑ = (r cos ϑ, 0, 0) · (cos ϑ cos φ, cos ϑ sin φ, − sin ϑ) =<br />
= r cos 2 ϑ cos φ<br />
F φ = F ⃗ · ⃗e φ = −r cos ϑ sin φ<br />
Beim Differenzieren von Vektoren (in allgemeinen krummlinigen <strong>Koordinaten</strong>)<br />
ist nun zu beachten, dass nicht nur die Komponenten, sondern auch<br />
die Basisvektoren differenziert werden müssen.<br />
⃗F = 3 ∑<br />
i=1<br />
F ui ⃗e ui ⇒ ∂ ⃗ F<br />
∂u j<br />
= 3 ∑<br />
i=1<br />
( )<br />
∂Fui<br />
∂⃗e<br />
∂u j<br />
⃗e ui + F<br />
ui<br />
ui<br />
∂u j<br />
Ähnliches gilt für die Differenziation nach anderen Parametern.<br />
Beispiel. Man bestimme aus dem Ortsvektor ⃗r(t) durch Ableitung nach<br />
der Zeit den Geschwindigkeitsvektor ⃗v(t) .<br />
Im Falle kartesischer <strong>Koordinaten</strong> sind die Basisvektoren konstante Vektoren<br />
und deren Ableitung verschwindet. Man braucht also nur die Komponenten<br />
nach t ableiten.<br />
Betrachten wir nun Zylinderkoordinaten. Hier gilt<br />
⃗v(t) = ˙⃗r(t) = d dt (ρ⃗e ρ + x 3 ⃗e x3 ) = ˙ρ⃗e ρ + ρ ˙⃗e ρ + ẋ 3 ⃗e x3 + x 3 ˙⃗ex3<br />
Die Ableitung eines Einheitsvektors steht senkrecht auf diesen und ist<br />
damit eine Linearkombination der beiden anderen Vektoren. In unserem<br />
Fall gilt<br />
˙⃗e ρ = ˙φ⃗e φ , ˙⃗ex3 = ⃗0 und folglich ⃗v(t) = ˙ρ⃗e ρ + ρ ˙φ⃗e φ + ẋ 3 ⃗e x3<br />
Beispiel. Beschreiben wir die Bewegung eines Punktes durch einen<br />
zeitabhängigen Ortsvektor ⃗r(t) in Kugelkoordinaten, so gilt<br />
5
⃗r(t) = r(t)⃗e r ⇒ ⃗v(t) = ˙⃗r(t) = r(t) d⃗e r<br />
dt + dr(t)<br />
dt ⃗e r .<br />
Da<br />
d⃗e r<br />
dt ⊥⃗e r<br />
gilt, muß d⃗e r<br />
dt<br />
durch ⃗e ϑ und ⃗e φ darstellbar sein.<br />
Tatsächlich gilt<br />
d⃗e r<br />
dt<br />
= ˙ϑ⃗e ϑ + sin ϑ ˙φ⃗e φ und damit ist<br />
⃗v(t) = dr(t)<br />
dt ⃗e r + r(t) ˙ϑ⃗e ϑ + r(t) sin ϑ ˙φ⃗e φ<br />
⃗v(t) = ṙ⃗e r + r ˙ϑ⃗e ϑ + r sin ϑ ˙φ⃗e φ .<br />
bzw.<br />
Bemerkung. Die Komponenten eines Vektors sind immer bezogen auf<br />
ein zuvor definiertes Basissystem.<br />
Der Vektor ⃗a habe in der kartesischen Basis ⃗e 1 , ⃗e 2 , ⃗e 3<br />
⃗a = 5x 1 ⃗e 1 + 5x 2 ⃗e 2 + 5x 3 ⃗e 3 .<br />
Dann sind 5x 1 , 5x 2 , 5x 3<br />
Basis.<br />
die Darstellung<br />
die Komponenten von ⃗a bzgl. der kartesischen<br />
Diese Komponenten (bzgl. der kartesischen Basis !) können auch durch<br />
Kugelkoordinaten ausgedrückt werden.<br />
5x 1 = 5r sin ϑ cos φ , 5x 2 = 5r sin ϑ sin φ , 5x 3 = 5r cos ϑ<br />
Wählen wir nun die Basis der Kugelkoordinaten ⃗e r , ⃗e ϑ , ⃗e φ , dann gilt<br />
⃗a = 5r⃗e r .<br />
Die Komponenten von ⃗a bzgl. der Basis der Kugelkoordinaten sind also<br />
5r, 0, 0 .<br />
Diese Komponenten können wiederum durch kartesische <strong>Koordinaten</strong> ausgedrückt<br />
werden, nämlich durch 5 √ x 2 1 + x2 2 + x2 3 , 0, 0 .<br />
Als nächstes untersuchen wir die Form von ∇Φ , ∇ · ⃗A , ∇ × A ⃗ , also<br />
Gradient, Divergenz und Rotation für orthonormale Basissysteme.<br />
• Gradient. ∇Φ(u 1 , u 2 , u 3 )<br />
6
∑<br />
∇Φ = 3 (∇Φ) ui ⃗e ui , wobei (∇Φ) ui zu bestimmen sind.<br />
∂Φ<br />
∂u i<br />
i=1<br />
= 3 ∑<br />
j=1<br />
∂Φ<br />
∂x j<br />
∂x j<br />
∂u i<br />
Also<br />
1<br />
h ui<br />
∂Φ<br />
∂u i<br />
= ∇Φ(u 1 , u 2 , u 3 ) ·<br />
= ∇Φ(x 1 (u 1 , u 2 , u 3 ), x 2 (u 1 , u 2 , u 3 ), x 3 (u 1 , u 2 , u 3 )) · ∂⃗r<br />
∂u i<br />
(<br />
1<br />
h ui<br />
∂⃗r<br />
∂u i<br />
)<br />
= ∇Φ(u 1 , u 2 , u 3 ) · ⃗e ui = (∇Φ) ui .<br />
Der Nabla-Operator hat also im neuen Basissystem die Form<br />
∇ = 3 ∑<br />
i=1<br />
1<br />
h ui<br />
∂<br />
∂u i<br />
⃗e ui .<br />
Im besonderen gilt damit etwa ∇u i = 1<br />
h ui<br />
⃗e ui .<br />
Bemerkung. Bilden ⃗e u1 , ⃗e u2 , ⃗e u3 ein Rechtssystem, dann ist<br />
∇u 2 × ∇u 3 = 1<br />
h u2<br />
⃗e u2 × 1<br />
h u3<br />
⃗e u3 = 1<br />
h u2 h u3<br />
⃗e u2 × ⃗e u3 = 1<br />
h u2 h u3<br />
⃗e u1 .<br />
Folglich ist ⃗e u1 = h u2 h u3 ∇u 2 × ∇u 3 . Analog gilt<br />
⃗e u2 = h u3 h u1 ∇u 3 × ∇u 1 und ⃗e u3 = h u1 h u2 ∇u 1 × ∇u 2 .<br />
Damit zeigen wir jetzt, dass mit einer Skalarfunktion ψ gilt<br />
1 ∂<br />
(a) ∇ · (ψ⃗e u1 ) =<br />
h u1 h u2 h u3 ∂u 1<br />
(ψh u2 h u3 )<br />
(b) ∇ × (ψ⃗e u1 ) = 1<br />
h u3 h u1<br />
∂<br />
∂u 3<br />
(ψh u1 )⃗e u2 − 1<br />
h u1 h u2<br />
∂<br />
∂u 2<br />
(ψh u1 )⃗e u3<br />
(c) Analoges gilt bzgl. der Vektoren ψ⃗e u2 und ψ⃗e u3 .<br />
Beweis.<br />
zu a) ∇ · (ψ⃗e u1 ) = ∇ · (ψh u2 h u3 ∇u 2 × ∇u 3 ) =<br />
= ∇(ψh u2 h u3 ) · (∇u 2 × ∇u 3 ) + ψh u2 h u3 ∇ · (∇u 2 × ∇u 3 ) =<br />
= ∇(ψh u2 h u3 ) · ( 1<br />
h u2<br />
⃗e u2 × 1<br />
h u3<br />
⃗e u3 ) + 0 = ∇(ψh u2 h u3 ) ·<br />
1<br />
h u2 h u3<br />
⃗e u1 =<br />
= ( 1<br />
h u1<br />
∂<br />
∂u 1<br />
(ψh u2 h u3 )⃗e u1 + 1<br />
h u2<br />
∂<br />
∂u 2<br />
(ψh u2 h u3 )⃗e u2 + 1<br />
h u3<br />
∂<br />
1<br />
∂u 3<br />
(ψh u2 h u3 )⃗e u3 )·<br />
h u2 h u3<br />
⃗e u1 =<br />
7
=<br />
1<br />
h u1 h u2 h u3<br />
∂<br />
∂u 1<br />
(ψh u2 h u3 )<br />
zu b) ∇ × (ψ⃗e u1 ) = ∇ × (ψh u1 ∇u 1 ) = ∇(ψh u1 ) × ∇u 1 + ψh u1 ∇ × ∇u 1 =<br />
= ∇(ψh u1 ) × 1<br />
h u1<br />
⃗e u1 + ⃗0 =<br />
= ( 1<br />
h u1<br />
= 1<br />
h u3 h u1<br />
∂<br />
∂u 1<br />
(ψh u1 )⃗e u1 + 1<br />
h u2<br />
∂<br />
∂u 3<br />
(ψh u1 )⃗e u2 − 1<br />
h u1 h u2<br />
∂<br />
∂u 2<br />
(ψh u1 )⃗e u2 + 1<br />
h u3<br />
∂<br />
∂u 2<br />
(ψh u1 )⃗e u3<br />
∂<br />
∂u 3<br />
(ψh u1 )⃗e u3 ) × 1<br />
h u1<br />
⃗e u1 =<br />
• Divergenz. ∇ · ⃗A(u 1 , u 2 , u 3 )<br />
⃗A = A u1 ⃗e u1 + A u2 ⃗e u2 + A u3 ⃗e u3 . Damit ist ∇ · ⃗A =<br />
= ∇ · (A u1 ⃗e u1 + A u2 ⃗e u2 + A u3 ⃗e u3 ) = ∇ · (A u1 ⃗e u1 ) + ∇ · (A u2 ⃗e u2 ) + ∇ · (A u3 ⃗e u3 )<br />
Mit dem vorhergehenden Ergebnis folgt sofort<br />
∇ · ⃗A(u 1 , u 2 , u 3 ) =<br />
]<br />
∂u 1<br />
(h u2 h u3 A u1 ) + ∂<br />
∂u 2<br />
(h u1 h u3 A u2 ) + ∂<br />
∂u 3<br />
(h u1 h u2 A u3 )<br />
[<br />
1 ∂<br />
=<br />
h u1 h u2 h u3<br />
• Rotation. ∇ × ⃗ A(u 1 , u 2 , u 3 )<br />
⃗A = A u1 ⃗e u1 + A u2 ⃗e u2 + A u3 ⃗e u3 . Damit ist ∇ × ⃗ A =<br />
= ∇ × (A u1 ⃗e u1 + A u2 ⃗e u2 + A u3 ⃗e u3 ) =<br />
= ∇ × (A u1 ⃗e u1 ) + ∇ × (A u2 ⃗e u2 ) + ∇ × (A u3 ⃗e u3 )<br />
Mit dem vorhergehenden Ergebnis ist dann<br />
∇ × ⃗ A(u 1 , u 2 , u 3 ) =<br />
= 1<br />
h u3 h u1<br />
+ 1<br />
h u1 h u2<br />
+ 1<br />
h u2 h u3<br />
∂<br />
∂u 3<br />
(A u1 h u1 )⃗e u2 − 1<br />
h u1 h u2<br />
∂<br />
∂u 1<br />
(A u2 h u2 )⃗e u3 − 1<br />
h u2 h u3<br />
∂<br />
∂u 2<br />
(A u3 h u3 )⃗e u1 − 1<br />
h u3 h u1<br />
∂<br />
∂u 2<br />
(A u1 h u1 )⃗e u3 +<br />
∂<br />
∂u 3<br />
(A u2 h u2 )⃗e u1 +<br />
∂<br />
∂u 1<br />
(A u3 h u3 )⃗e u2 =<br />
8
[<br />
= 1 ∂<br />
h u2 h u3<br />
[<br />
+ 1<br />
h u3 h u1<br />
[<br />
+ 1<br />
h u1 h u2<br />
]<br />
∂u 2<br />
(A u3 h u3 ) − ∂<br />
∂u 3<br />
(A u2 h u2 ) ⃗e u1 +<br />
]<br />
∂<br />
∂u 3<br />
(A u1 h u1 ) − ∂<br />
∂u 1<br />
(A u3 h u3 ) ⃗e u2 +<br />
]<br />
∂<br />
∂u 1<br />
(A u2 h u2 ) − ∂<br />
∂u 2<br />
(A u1 h u1 ) ⃗e u3<br />
Dies kann wiederum in folgender Form geschrieben werden.<br />
∣ ∣ ∣∣∣∣∣ h<br />
∇ × A(u ⃗ u1 ⃗e u1 h u2 ⃗e u2 h u3 ⃗e u3 ∣∣∣∣∣<br />
1<br />
∂ ∂ ∂<br />
1 , u 2 , u 3 ) =<br />
h u1 h u2 h u3 ∂u 1 ∂u 2 ∂u 3<br />
h u1 A u1 h u2 A u2 h u3 A u3<br />
Bemerkung. Daraus folgt etwa, dass auch in krummlinigen (orthonormalen)<br />
<strong>Koordinaten</strong>systemen der Rotor eines Gradienten verschwindet.<br />
Sei<br />
⃗ A = ∇Φ(u1 , u 2 , u 3 ) wobei Φ zweimal stetig differenzierbar ist.<br />
∣ h u1 ⃗e u1 h u2 ⃗e u2 h u3 ⃗e u3 ∣∣∣∣∣<br />
rot grad Φ(u 1 , u 2 , u 3 ) = ∇ × (∇Φ) =<br />
∂ ∂ ∂<br />
∂u 1 ∂u 2 ∂u 3<br />
= 0<br />
∣<br />
∂Φ<br />
∂u 3<br />
∂Φ<br />
∂u 1<br />
∂Φ<br />
∂u 2<br />
Bemerkung. Damit kann auch der Laplace-Operator in orthonormalen<br />
<strong>Koordinaten</strong>systemen angegeben werden.<br />
△Φ(u 1 , u 2 , u 3 ) =div grad Φ(u 1 , u 2 , u 3 ) = ∇ · ∇Φ(u 1 , u 2 , u 3 ) =<br />
) ( ) (<br />
∂Φ<br />
h u1 ∂u 1<br />
+ ∂ hu3 h u1 ∂Φ<br />
∂u 2 h u2 ∂u 2<br />
+ ∂ hu1 h u2<br />
∂u 3 h u3<br />
{ (<br />
1 ∂ hu2 h<br />
=<br />
u3<br />
h u1 h u2 h u3 ∂u 1<br />
)}<br />
∂Φ<br />
∂u 3<br />
Beispiel. Zylinderkoordinaten ρ, φ, x 3<br />
[<br />
]<br />
∂<br />
∇Φ(ρ, φ, x 3 ) = ⃗e ρ ∂ρ + 1 ρ ⃗e φ ∂<br />
∂φ + ⃗e ∂<br />
x 3 ∂x 3<br />
Φ(ρ, φ, x 3 )<br />
[<br />
∆Φ(ρ, φ, x 3 ) =<br />
1 ∂<br />
ρ<br />
∂ρ ρ ∂ ∂ρ + 1 ρ 2 ∂ 2<br />
∂φ 2 + ∂2<br />
∂x 2 3<br />
]<br />
Φ(ρ, φ, x 3 ) =<br />
9
[<br />
=<br />
∂ 2<br />
∂ρ<br />
+ 1 ∂<br />
2 ρ ∂ρ + 1 ∂ 2<br />
ρ 2 ∂φ<br />
+ ∂2<br />
2<br />
∂x 2 3<br />
]<br />
Φ(ρ, φ, x 3 )<br />
Der Laplace-Operator im R 2<br />
weggelassen wird.<br />
x 3<br />
ergibt sich daraus, indem die Ableitung nach<br />
Beispiel. Kugelkoordinaten r, ϑ, φ<br />
[<br />
]<br />
∂<br />
∇Φ(r, ϑ, φ) = ⃗e r ∂r + 1 r ⃗e ϑ ∂<br />
∂ϑ + 1<br />
r sin ϑ ⃗e φ ∂<br />
∂φ<br />
Φ(r, ϑ, φ)<br />
[<br />
∆Φ(r, ϑ, φ) =<br />
[<br />
=<br />
∂<br />
∂r 2 + 2 r<br />
1<br />
r 2<br />
( )<br />
∂<br />
∂r r<br />
2 ∂<br />
∂r +<br />
1<br />
∂<br />
∂r + 1 ∂ 2<br />
r 2 ∂ϑ<br />
+ cos ϑ<br />
2<br />
∂<br />
r 2 sin ϑ ∂ϑ<br />
∂<br />
r 2 sin ϑ ∂ϑ + 1<br />
r 2 sin 2 ϑ<br />
( )<br />
sin ϑ<br />
∂<br />
∂ϑ +<br />
1<br />
∂ 2<br />
r 2 sin 2 ϑ<br />
∂φ 2 ]<br />
Φ(r, ϑ, φ)<br />
∂ 2<br />
∂φ 2 ]<br />
Φ(r, ϑ, φ) =<br />
Zum Abschluß betrachten wir noch Bogen-, Flächen- und Volumenselement<br />
in orthonormalen <strong>Koordinaten</strong>.<br />
∂⃗r<br />
Aus<br />
∂u i<br />
= h ui ⃗e ui ergibt sich für die Differenziale in Richtung der u i -<br />
<strong>Koordinaten</strong>linie<br />
d⃗r ui = ∂⃗r<br />
∂u i<br />
du i = h ui du i ⃗e ui .<br />
Der Größe h ui du i entspricht daher die lineare Näherung eines kleinen<br />
Bogenabschnittes entlang einer u i -Linie.<br />
Das totale Differenzial ds der Bogenlänge ist weiters (für orthogonale<br />
<strong>Koordinaten</strong>) gegeben durch<br />
ds = √ (h u1 du 1 ) 2 + (h u2 du 2 ) 2 + (h u3 du 3 ) 2<br />
Damit kann auch entlang Kurven in krummlinigen <strong>Koordinaten</strong> integriert<br />
werden.<br />
Beispiel. Bestimme die Länge eines spiralförmigen Drahtes (Zylinderradius<br />
a , Ganghöhe des Drahtes 2πb ).<br />
10
Die Kurve ist durch ρ = a , 0 ≤ φ ≤ 4π , z = bφ gegeben.<br />
(ds) 2 = (dρ) 2 + ρ 2 (dφ) 2 + (dz) 2 = 0 + a 2 (dφ) 2 + b 2 (dφ) 2<br />
Also ist ds = √ a 2 + b 2 dφ und S =<br />
∫4π<br />
0<br />
√<br />
a2 + b 2 dφ = 4 √ a 2 + b 2 π<br />
Für das Flächenelement einer <strong>Koordinaten</strong>fläche ⃗r(u 1 , u 2 , u 0 3) gilt<br />
dA(u 1 , u 2 ) = |(h u1 ⃗e u1 ) × (h u2 ⃗e u2 )|du 1 du 2 = |h u1 h u2 |du 1 du 2 .<br />
Analoges gilt für die anderen <strong>Koordinaten</strong>flächen. Für einen Zylindermantel<br />
(<strong>Koordinaten</strong>fläche bei Zylinderkoordinaten) gilt also etwa<br />
dA(φ, x 3 ) = ρ 0 dφdx 3<br />
Das Volumselement wird in linearer Näherung als kleiner Quader<br />
dV = |(h u1 du 1 ⃗e u1 ) · [(h u2 du 2 ⃗e u2 ) × (h u3 du 3 ⃗e u3 )]| =<br />
= |h u1 h u2 h u3 [e u1 · (e u2 × e u3 )]|du 1 du 2 du 3 = |h u1 h u2 h u3 |du 1 du 2 du 3<br />
Bei dem auftretenden Faktor handelt es sich um eine Determinante, die<br />
Jacobi-Determinante, die schon erwähnt wurde. Im speziellen erhalten wir<br />
für<br />
• Zylinderkoordinaten: dV = ρdρdφdx 3<br />
• Kugelkoordinaten:<br />
dV = r 2 sin ϑdrdϑdφ<br />
11