02.11.2013 Aufrufe

Krummlinige Koordinaten

Krummlinige Koordinaten

Krummlinige Koordinaten

MEHR ANZEIGEN
WENIGER ANZEIGEN

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

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!