Für einen reinen Text-Ausdruck kann man die Bilder-Darstellung im Browser deaktivieren:
Internet-Explorer: Extras > Internetoptionen > Erweitert > Bilder anzeigen (Häkchen entfernen und mit OK bestätigen)
Mozilla Firefox: Extras > Einstellungen > Inhalt > Grafiken laden (Häkchen entfernen und mit OK bestätigen)
Opera: Ansicht > Bilder > Keine Bilder
Texte, die in Klappboxen stehen, werden nicht immer ausgedruckt (abhängig von der Definition). Auf jeden Fall müssen sie ausgeklappt sein, wenn sie gedruckt werden sollen.
Die Funktion „Als PDF herunterladen“ kann zu Darstellungsfehlern führen.
Die Vektoranalysis wendet die Methoden der Analysis (Differential- und Integralrechnung) auf mathematische Funktionen an, in denen Vektoren auftreten, die sich in Abhängigkeit von Ort und Zeit verändern können. Die wichtigsten Anwendungsgebiete der Vektoranalysis sind physikalische Felder, insbesondere elektromagnetische Felder.
sind Teilgebiete des Raumes , in denen jedem Punkt eindeutig ein Skalar oder ein Vektor (auch ein Tensor oder Spinor) zugeordnet ist. Je nach Art der »Feldgröße« spricht man von einem Skalarfeld oder einem Vektorfeld.
Skalare Feldgrößen sind z. B. Druck, Temperatur, Beleuchtungsstärke, Potential.
Vektorielle Feldgrößen sind z. B. elektrische und magnetische Feldstärke, magnetische Induktion, Strömungsgeschwindigkeit.
Kraftfelder sind Felder, in denen z. B. eine elektrische Ladung oder eine Masse eine Kraft erfährt.
Elektrodynamische Felder sind zeitlich veränderliche elektrische und magnetische Felder, in denen Induktionsvorgänge stattfinden.
Feldlinien sind (gedachte) Linien dergestalt, dass die Vektoren der Feldgröße ihre Tangenten sind. Bekannte Beispiele sind: Stromlinien, elektrische und magnetische Feldlinien.
Zur Schreibweise: Im Text werden – der deutschen Norm folgend – die Zeichen für Vektoren (V) kursiv und fett geschrieben. In den mit TeX gesetzten Formeln sind die Zeichen mit einem Pfeil versehen.
Für die Beschreibung eines Vektors durch seine kartesischen Komponenten sind drei Schreibweisen üblich: Mittels der Einheitsvektoren i, j, k auf der X-, Y- und Z-Achse, als einzeilige Matrix und als einspaltige Matrix.
Ich werde diese Schreibweisen je nach Zweckmäßigkeit abwechselnd verwenden.
Eine Funktion, bei der die abhängige Variable ein Vektor ist, heißt Vektorfunktion. Im einfachsten Fall sind die (skalaren) kartesischen Komponenten Vx, Vy, Vz des Vektors Funktionen einer einzigen Variablen u (einparametrige Vektorfunktion).
Wenn φ alle reellen Zahlenwerte annimmt, durchläuft der Punkt P eine Schraubenlinie mit dem Radius a und der Ganghöhe h. Die Ableitung dieser Vektorfunktion ist der Vektor
Setzen wir
wobei t die Zeit sein soll, so hat P die konstante Winkelgeschwindigkeit ω und den Geschwindigkeitsvektor
Anwendungen auf die Differentialgeometrie der Raumkurven
Die Tangente an eine Raumkurve mit dem Ortsvektor r(u) in einem ihrer Punkte P ist die Gerade durch P mit derselben Richtung wie der Vektor (dr/du)P (das bedeutet: die Vektorfunktion dr/du gebildet an der Stelle P).
Dabei ist u irgendeine Variable, durch die r beschrieben wird.
Diese Definition wird sofort plausibel, wenn wir die Variable u durch die Zeit t ersetzen. Dann ist (siehe oben):
Der Geschwindigkeitsvektor v gibt aber die momentane Bewegungsrichtung des Punktes P an, und das ist die Richtung der Kurventangente.
Führen wir nun wieder die beliebige Variable u ein, dann wird
wobei dt/du lediglich ein skalarer Faktor ist, der an der Richtung des Vektors v nichts ändert. Also hat auch der Vektor die Richtung von v und damit die Richtung der Tangente.
Für das Folgende brauchen wir den auf der Kurventangente gelegenen Einheitsvektor. Er wird mit t bezeichnet. Man findet ihn, indem man den Geschwindigkeitsvektor durch seinen Betrag dividiert:
wobei v die Bahngeschwindigkeit des Punktes P ist. Ist r als Funktion der Bogenlänge s der Kurve gegeben, wobei s von einem beliebigen Punkt der Kurve aus gemessen wird, dann kann man v durch ds/dt ersetzen:
Es ist nützlich, sich zunächst die analogen Überlegungen und Begriffe bei einer ebenen Kurve zu vergegenwärtigen. Dort liegen selbstverständlich auch alle Kurventangenten in derselben Ebene, der Ebene der Kurve. Ändert sich die Richtung der Tangente (ihr Winkel) auf der Weglänge (Bogenlänge) Δs um den Wert Δ τ so ist die »mittlere Krümmung« km auf der Strecke Δs
und die Krümmung der Kurve im betrachteten Punkt P
Unter dem Krümmungskreis der Kurve im Punkt P versteht man den Kreis durch P, der dieselbe Steigung und dieselbe Krümmung wie die Kurve in P hat. Der Radius ρ dieses Kreises heißt Krümmungsradius der Kurve in P. Es gilt
Die Tangenten einer Raumkurve liegen nicht in derselben Ebene und es gibt – im Gegensatz zu Flächen – im Punkt P auch nicht nur eine Tangentialebene, sondern unendlich viele. Unter ihnen greifen wir die Ebene heraus,in der der Tangenteneinheitsvektor t und der Vektor dt/ds liegen. Der letztgenannte Vektor gibt nämlich die Richtung an, in welcher sich der Vektor t in P dreht. Diese Ebene heißt die Schmiegungsebene der Kurve in P.
Der in der Schmiegungsebene liegende Einheitsvektor, der auf t senkrecht steht und dieselbe Richtung wie der Vektor dt/ds hat, heißt Hauptnormaleneinheitsvektor n der Kurve in P.
Hat ein Vektor v(u) eine konstante Länge v, so ist wegen v2 = v2 auch v2 = konst. Differenziert man die letzte Gleichung nach u und benutzt dabei die Regel für die Differentiation eines Skalarprodukts v·w mit w = v, so findet man
Wenn das Skalarprodukt zweier Vektoren v und dv/du null ist und keiner der beiden Vektoren selbst null ist (Nullvektor bzw. konstanter Vektor), dann müssen die beiden Vektoren aufeinander senkrecht stehen. Dies leuchtet auch unmittelbar ein: Wenn der Vektor dv/du eine Komponente in Richtung v hätte, dann würde sich die Länge von v zugleich mit u verändern.
Dieses Ergebnis wenden wir auf den Tangenteneinheitsvektor t einer Raumkurve an. Da die Länge von t konstant ist, muss seine Ableitung dt/ds auf t senkrecht stehen.
In der Abbildung liegen der Tangenten- und der Normalenvektor in der Zeichenebene, die folglich mit der Schmiegungsebene zusammenfällt. Die Kurve selbst dagegen verläuft im Allgemeinen links und rechts von P außerhalb dieser Ebene.
Unter der mittleren Krümmung einer Kurve im Bereich Δs versteht man den auf Δs bezogenen Drehwinkel Δτ der Tangente. Ihr Grenzwert für Δs gegen 0 heißt Krümmung k der Kurve im Punkt P.
Ein in der Schmiegungsebene gelegener Kreis durch P mit derselben Steigung und derselben Krümmung wie die Raumkurve, heißt Krümmungskreis der Kurve. Sein Radius heißt Krümmungsradius ρ der Kurve in P. Da für den Kreisbogen Δs (unabhängig von seiner Größe) stets gilt
gilt für seine Krümmung
Zur Berechnung der Krümmung einer Kurve aus ihrem Ortsvektor r(s) gehen wir wie folgt vor:
1. Berechnung von dt/ds:
2. Berechnung von dt/dτ:
Es ist:
und
3. Damit ergibt sich:
Da der Vektor dt/ds die Richtung des Normaleneinheitsvektors n hat, ist
In Integralen können Vektoren sowohl als Integrand (= die zu integrierende Funktion) als auch als Differential bei dem Integranden auftreten.
1. Typ: Nur der Integrand ist ein Vektor
Ein typisches Beispiel ist das Zeitintegral der Kraft, das in der Dynamik auftritt. (Dort ist es ein bestimmtes Integral; es genügt hier jedoch, nur unbestimmte Integrale zu untersuchen.)
Das Ergebnis ist also, wie zu erwarten war, ein Vektor.
Anmerkung: Dass oben die Einheitsvektoren i, j, k wie konstante Faktoren vor die Integrale gezogen werden dürfen, lässt sich wie folgt beweisen: Das Integralzeichen ist das Symbol für den Grenzwert einer Summe. Konstante Faktoren bei den Summanden können ausgeklammert werden, auch wenn sie (konstante) Vektoren sind.
2. Typ: Integrand und Differential sind Vektoren
Ein Beispiel dafür ist das Wegintegral der Kraft, mit dem die Arbeit berechnet wird.
Da F·dr ein Skalarprodukt ist, ergibt sich für das Ergebnis des Integrals erwartungsgemäß auch ein Skalar.
Ein spezielles wichtiges Beispiel hierfür ist:
Die Integration folgt hier formal derselben Regel wie bei .
Ein anderes interessantes Beispiel (unter Verwendung des erst später erklärten Operators grad, dessen Bedeutung hier erkennbar ist):
Erläuterung: Der Integrand im vorletzten Integral ist das vollständige Differential dU der Funktion U = U(x, y, z).
Ein physikalisches Feld ist – wie eingangs schon erklärt – ein Teilgebiet des Raumes, in welchem in jedem Punkt eine eindeutig bestimmte skalare oder vektorielle physikalische Größe (»Feldgröße« genannt) anzutreffen ist.
Bei Skalarfeldern ist die (skalare) Feldgröße U eine skalare Funktion des Ortsvektors r des betrachteten Punktes P:
Beispiele für skalare Feldgrößen aus der Physik sind Druck und Temperatur in der Atmosphäre, das Gravitationspotential in der Umgebung einer Masse (z. B. der Erde), das Potential in der Umgebung eines elektrisch geladenen Körpers, die Lautstärke in einem Schallfeld.
Bei Vektorfeldern ist die (vektorielle) Feldgröße V eine Vektorfunktion von r:
Beispiele für vektorielle Feldgrößen sind die elektrische und die magnetische Feldstärke, die Gravitationsfeldstärke, die Geschwindigkeit von Gasen und Flüssigkeiten in Strömungsfeldern.
Ein wichtiges Beispiel für ein Vektorfeld und ein Skalarfeld
Die elektrische Feldstärke E im Feld einer punktförmigen elektrischen Ladung vom Betrag Q, die sich in O befindet, ist
Das Potential φ eines Punktes P in einem beliebigen elektrischen Feld ist definiert als die »ladungsbezogene Arbeit« W/q, die aufzuwenden ist, um die Ladung q aus unendlicher Entfernung zu dem Punkt P zu bringen. (Ein Punkt eines Feldes besitzt nur dann ein definiertes Potential, wenn diese Arbeit vom Weg unabhängig ist, auf dem die Ladung nach P gebracht wird. – Dieses Problem wird später noch genauer untersucht.) Also:
Diese Definition gilt analog auch für das Potential eines Gravitationsfeldes, wobei lediglich q durch die Masse m des bewegten Körpers zu ersetzen ist.
Für das oben beschriebene zentralsymmetrische elektrische Feld, in dem – wie später gezeigt wird – jedem Punkt ein Potential zugeordnet werden kann, errechnet man die Arbeit durch eine Integration:
Da in diesem Feld die Arbeit vom gewählten Weg unabhängig ist, denken wir uns die Ladung q einfach radial nach innen bewegt, wobei dann Kraft- und Wegvektor gleich- oder entgegengesetzt gerichtet sind. Allerdings ist der Vektor ds dem Vektor dr entgegengesetzt gerichtet, da die Bewegung in Richtung abnehmendem r erfolgt: ds = - dr.
In einem Punkt mit der Feldstärke E erfährt die Ladung q eine Kraft vom Betrag
F = Eq, also ist
Damit erhalten wir für das Potential
Da das Potential (definitionsgemäß) ein Skalar ist, ist das Potentialfeld ein Skalarfeld.
Für r = konst. ist auch φ = konst. Die Punkte gleichen Potentials liegen also auf einer Kugelfläche um O. Die »Äquipotentialflächen« oder »Niveauflächen« dieses Feldes sind also Kugeln (siehe Abbildung). Das elektrische Potential wird in Volt (V) gemessen.
Wir begeben uns nun zu einem Punkt P(x, y, z) eines Skalarfeldes mit der Feldgröße U(r) und fragen zunächst nach dem Anstieg ΔU der Feldgröße auf der Strecke Δs und dann nach der mittleren Steigung ΔU/Δs der Feldgröße auf derselben Strecke. (Die Feldgröße könnte z. B. die Temperatur, der Luftdruck oder das Potential eines Feldes sein. )
Dazu brauchen wir zunächst die Steigung der Feldgröße in Richtung der drei Koordinatenachsen.
Die Differentialrechnung liefert für die Steigung der Feldgröße U in Richtung der drei Koordinatenachsen im Punkt P:
Der Anstieg ΔU der Feldgröße U längs einer Strecke Δ s = (Δx Δy Δz) ist dann - dies ist ebenfalls ein Ergebnis der Differentialrechnung - für hinreichend kleine Δx, Δy, Δz
und die mittlere Steigung ΔU/Δs der Feldgröße U auf der Strecke Δ s ist
Die Quotienten Δx/Δs, Δy/Δs, Δz/Δs sind die Richtungskosinus des Vektors Δ s:
Sie sind vom Betrag Δs des Vektors Δs unabhängig und bleiben auch für Δs gegen null unverändert. Sie sind die skalaren Komponenten des Einheitsvektors eΔs, der dieselbe Richtung hat wie Δ s:
Damit wird
Daraus ergibt sich für Δs gegen null die Steigung der Feldgröße U in Richtung Δs:
oder
Richtungsableitung und Gradient einer skalaren Feldgröße
Der soeben gefundene Term für die Steigung der Feldgröße U in der durch den Vektor (cos α cos β cos γ) beschriebenen Richtung kann interpretiert werden als das Skalarprodukt des Vektors
und des Vektors
Der Vektor v hat bemerkenswerte, für die Untersuchung von Feldern sehr nützliche Eigenschaften, weshalb er einen eigenen Namen erhalten hat: Gradient U (grad U). (»Gradient« ist ein aus einem lateinischen Stamm abgeleitetes Kunstwort, das man etwa mit »Steigungszeiger« übersetzen könnte.) Damit gilt für die so genannte Richtungsableitung der Feldgröße U in der Richtung des Vektors (cos α cos β cos γ)
Die besonderen Eigenschaften des Vektors grad U ergeben sich so:
Das Skalarprodukt zweier Vektoren v und w ist gleich dem Produkt ihrer Beträge v und w und dem Kosinus des Winkels δ zwischen den beiden Vektoren:
Bei gegebenen Werten von v und w ist der Wert des Skalarprodukts maximal (nämlich gleich u v), wenn δ = 0 ist. Die Richtungsableitung (= Steigung) der Feldgröße U ist also dann am größten, wenn der Vektor eΔs (oder der Vektor Δ s) dieselbe Richtung wie der Vektor grad U hat. Anders herum gesagt:
Der Vektor grad U weist in die Richtung, in der die Feldgröße U die größte Steigung hat (am stärksten steigt).
Steht dagegen der Vektor eΔs auf dem Vektor grad U senkrecht, dann ist dU/ds = 0. Das bedeutet, der Vektor eΔs liegt in der Tangentialebene der Niveaufläche U = konst. des betrachteten Punktes P. Daraus folgt:
Der Vektor grad U steht auf der Niveaufläche durch den Punkt P senkrecht.
Ferner: Der Maximalwert der Steigung (oder der Richtungsableitung) ist der Maximalwert des obigen Skalarprodukts:
Das bedeutet: Der Betrag des Vektors grad U ist gleich dem Maximalwert der Steigung der Feldgröße im betrachteten Punkt.
Beispiel: Gesucht ist der Gradient des Potentials φ einer elektrischen Punkt- (oder Kugel-)ladung Q.
Für das Potential gilt, wie früher gezeigt wurde,:
Die partiellen Ableitungen werden am einfachsten nach der Kettenregel gebildet:
Die Physik lehrt, dass elektrostatische Felder und stationäre Gravitationsfelder - unabhängig von der Anzahl und der Anordnung der Ladungen bzw. Massen, die das Feld aufbauen - so genannte Potentialfelder sind. Das bedeutet: Um eine Ladung q bzw. eine Masse m aus unendlicher Entfernung zu einem bestimmten Punkt P des Feldes zu bringen, ist eine (positive oder negative) Arbeit aufzuwenden, die unabhängig von dem Weg ist, auf dem die Ladung bzw. die Masse transportiert wird.
Da die aufzuwendende Arbeit proportional der Ladung bzw. Masse ist, erhält man eine nur von der Lage des Punktes P abhängige skalare Größe, wenn man die Arbeit durch die Ladung bzw. Masse dividiert. Diese Größe, also die »ladungs- bzw. massebezogene Arbeit«, heißt das Potentialφ des Punktes P:
Der Vektor der Feldstärke ist definiert als die ladungs- bzw. massebezogene Kraft, die eine Ladung q bzw. eine Masse m in einem Punkt des Feldes erfährt:
Wegen der formalen Übereinstimmung der entsprechenden Gleichungen für das elektrische Feld und das Gravitationsfeld und wegen der sich dadurch anbietenden Vereinfachung bezeichne ich im Folgenden die Feldstärke allgemein und neutral mit V. Die Größe q kann fortan sowohl eine elektrische Ladung als auch eine »schwere Ladung«, das heißt eine Masse, bedeuten.
In einem Punkt A des Feldes habe das Potential den Wert φA, in einem Punkt B den Wert φB. Dann ist der Potentialunterschied der beiden Punkte
Dabei ist ΔWAB die Arbeit, die aufzuwenden ist, um die Ladung Q von A nach B zu transportieren. Für sie gilt:
wobei FA die in A auf die Ladung wirkende Kraft sein soll. Damit wird
Da die »arbeitende Kraft« der vom Feld auf die Ladung ausgeübte »Feldkraft« FFeld entgegengesetzt gleich und andererseits FFeld/q gleich der Feldstärke V ist, folgt
wobei Vx usw. die Komponenten des Vektors V an der Stelle A sind.
Wählt man ΔrAB so, dass ΔyAB = ΔzAB = 0 ist, dann wird daraus
Lässt man B unbeschränkt gegen A rücken, so wird
Analog findet man
Daraus folgt weiter (jetzt ohne Indices geschrieben):
Der Term auf der linken Seite aber ist der Vektor grad φ. Daher gilt für jedes beliebige Potentialfeld
Umgekehrt gelesen:
Der Feldstärkevektor eines jeden Potentialfeldes ist gleich dem negativen Gradienten des Potentials.
In einem Potentialfeld werde eine Ladung q gegen die Kraft des Feldes von A nach B verschoben. Die dazu aufzuwendende Arbeit ist
und wegen
(Der Integrand ist das vollständige Differential dφ des nur vom Ort abhängigen Potentials φ = φ(x, y, z).)
Also:
Die Arbeit W hängt also nur vom Potential des Anfangs- und Endpunktes des Weges ab, nicht aber vom Verlauf des Weges; das entsprechende Linienintegral ist »wegunabhängig«.
Wird die Ladung q zunächst auf einem beliebigen Weg von A nach B gebracht und danach auf einem anderen Weg von B zurück nach A, so ist
und daher
Das heißt: Das Linienintegral wird null, wenn man es über einen geschlossenen Weg bildet.
Zusammenfassung:
In einem Vektorfeld,
zu dem ein Potentialfeld gehört,
oder, was dasselbe ist, dessen Feldvektor der negative Gradient eines Skalarfeldes ist,
ist das Arbeitsintegral über einen geschlossenen Weg gleich null.
Das bedeutet, dass man durch Herumführen einer Ladung auf einem geschlossenen Weg weder Arbeit gewinnen kann noch Arbeit investieren muss.
Ein solches Vektorfeld und die in ihm auf eine Ladung ausgeübte Kraft heißen konservativ.
Welche Bedingungen muss der FeldvektorVerfüllen, damit er der negative Gradient eines Skalarfeldes mit der Feldfunktion U(x, y, z) sein kann?
Wenn
sein soll, muss
sein. Diese Forderung ist keineswegs selbstverständlich oder trivial, denn Vx, Vy und Vz können im Allgemeinen drei von einander völlig unabhängige Funktionen sein.
Nach dem Satz von SCHWARZ muss
und
sein. Das heißt: Bei der Bildung der zweiten partiellen Ableitung nach verschiedenen Variablen ist die Reihenfolge beliebig.
Auf unser Problem angewendet, bedeutet das: Wenn der Feldvektor V der negative Gradient eines Skalarfeldes mit der Feldfunktion U sein soll, muss
sein. Dann und nur dann ist (Vx dx + Vy dy + Vz dz) das vollständige Differential dU einer Funktion U, und nur dann kann daraus durch Integration eine Funktion U bestimmt werden, deren negativer Gradient dann der Vektor V ist. (Und nur dann ist auch der Wert des Arbeitsintegrals vom Weg unabhängig.)
Später wird sich zeigen, dass die oben beschriebene Bedingung identisch ist mit der Forderung, dass das Feld mit dem Feldvektor V wirbelfrei ist, (d. h., dass überall rot V = 0 ist.)
Beispiel:
Der Feldvektor
erfüllt – wie man leicht durch Rechnung bestätigen kann - die oben beschriebene »Integrabilitätsbedingung«, und es ist
Gegeben sei ein »Strömungsfeld« mit dem Feldvektor v(r), wobei v die Geschwindigkeit einer Flüssigkeit ist.
Stellen wir uns ein von einem Drahtrahmen umgrenztes ebenes Flächenstück vom Größenwert A vor, das so in die Flüssigkeit eintaucht, dass es auf der zunächst als homogen angenommenen Strömung senkrecht steht.
Dann strömt in der Zeitspanne Δt das Flüssigkeitsvolumen ΔV = v·Δt·A durch den Rahmen. Der Quotient aus diesem Volumen und der Zeitspanne Δt heißt der Fluss Φ der Strömung (oder auch – nicht ganz exakt, aber gebräuchlich - der Fluss Φ des Feldvektors v) durch das Flächenstück:
Der Fluss hat demnach die Dimension Volumen/Zeit =Länge3/Zeit.
Steht das Flächenstück auf der Strömungsrichtung nicht senkrecht, dann ist
wobei φ der Winkel zwischen dem Geschwindigkeitsvektor v und dem auf der Fläche senkrecht stehenden Flächenvektor A ist.
Dann ist der Fluss durch das Flächenstück
wobei v·A das Skalarprodukt der Vektoren v und A ist.
Ist schließlich das betrachtete Flächenstück nicht eben, oder ist das Strömungsfeld nicht homogen, dann denken wir uns die Fläche in hinreichend kleine Teilstücke vom Größenwert ΔA zerlegt und den Flächenvektor ΔA in der Mitte eines jeden Teilstücks errichtet. Jeder dieser Flächenvektoren wird dann skalar mit dem Geschwindigkeitsvektor multipliziert, der dem Fußpunkt des Flächenvektors zugeordnet ist. Für den Fluss Φ durch die gesamte Fläche A gilt dann:
Denkt man sich nun die Anzahl der Teilflächen unbegrenzt wachsend, wobei ΔA gegen null geht, dann strebt diese Summe einem Grenzwert zu, welcher der Fluss der Strömung (oder, wie man etwas nachlässig sagt, der Fluss des Vektors v ) durch die Fläche A ist und durch ein Flächenintegral dargestellt wird:
Dieser Begriff des Flusses wird in der Physik auch auf andere Vektorfelder übertragen, vor allem auf elektrische und magnetische Felder. Dies mag zunächst etwas befremden, aber man kann ja – als Hilfe für die Vorstellung - jeden beliebigen Feldvektor als den Geschwindigkeitsvektor einer Flüssigkeitsströmung interpretieren. Man muss dann lediglich, wann immer vom Fluss eines Feldvektors die Rede ist, sich vergegenwärtigen, dass damit eigentlich der Fluss einer »virtuellen Flüssigkeit« gemeint ist, deren Geschwindigkeitsvektor der betrachtete Feldvektor ist. Dazu das folgende Beispiel.
Der »Fluss des Feldvektors«
des zentralsymmetrischen Feldes einer Punkt- oder Kugelladung mit dem Mittelpunkt in O durch eine konzentrische Kugelfläche mit dem Radius R ist
(Hinweis: Der Feldvektor steht überall auf der Kugelfläche senkrecht.)
Also: Wenn E der Geschwindigkeitsvektor eines Strömungsfeldes wäre, betrüge der Fluss der Flüssigkeit durch jede zur Ladung Q konzentrische Kugelfläche Q/ε0. Der Fluss ist der Ladung also proportional, und für Q = 0 wäre auch Φ = 0. Demnach könnte man die Ladung Q als die »Quelle« des Feldes der virtuellen Flüssigkeit betrachten. Für Q < 0 wäre auch Φ < 0. Dies wäre so zu interpretieren: Der Geschwindigkeitsvektor ist - wie die Feldlinien des Feldes – nach innen gerichtet und bildet mit den Flächennormalen der Kugel überall den Winkel 180°, weshalb das Skalarprodukt E·dA = - E dA ist. Die negative Ladung ist dann die »Senke« (= Gegenteil einer Quelle) des Feldes.
Das Ergebnis Φ = Q/ε0 gilt übrigens, wie sich zeigen lässt und was auch durchaus plausibel erscheint, für jede beliebige, die Ladung Q umhüllende Fläche.
Zur Vereinfachung betrachte ich im Folgenden wieder einen »echten« Geschwindigkeitsvektor v eines Strömungsfeldes, jedoch gelten die Betrachtungen und ihre Ergebnisse für jedes beliebige Vektorfeld und sein virtuelles Strömungsfeld.
Integriert man das Skalarprodukt v·dA über eine geschlossene Fläche (»Hülle«), so ist der Wert des »Hüllenintegrals« gleich dem Fluss (Volumen/Zeit), der durch die Hülle nach außen tritt. Dieser muss gleich der »Schüttung« S (= Ergiebigkeit) aller innerhalb der Hülle liegenden Quellen sein, wobei die Senken einen negativen Beitrag zur Schüttung liefern:
Betrachten wir nun ein Raumgebiet vom Volumen ΔV. Die Schüttung aller Quellen in diesem Raumgebiet sei ΔS. Der Quotient ΔS/ΔV ist dann die »mittlere Quelldichte« in diesem Gebiet:
Lässt man nun die Hüllfläche auf einen Punkt P schrumpfen und somit ΔV gegen null gehen, so ist der Grenzwert
die Quelldichte des Feldvektors (eigentlich: des Strömungsfeldes, dessen Geschwindigkeitsvektor v ist) in dem Punkt P, auf den die Hülle geschrumpft ist. Sie wird als die Divergenz des Vektors v im Punkt P bezeichnet:
Zur Berechnung der Divergenz aus dem Feldvektor v = (vx vy vz) betrachten wir einen Quader mit den Seiten Δx, Δy, Δz, dessen Mittelpunkt der Punkt P (x, y, z) ist.
Die Flächennormalen auf den Seitenflächen sind die Einheitsvektoren in Achsenrichtung: i, j, k sowie -i, -j, -k. (j und der dazu gehörige Feldvektor sind nicht eingezeichnet.)
Die Flüsse durch die einzelnen Seitenflächen sind:
Der gesamte Fluss ΔΦ durch die Flächen des Quaders ist die Summe aus diesen sechs Flüssen. Er entspricht dem Wert des Hüllenintegrals in der Definition der Divergenz. Von den sechs Summanden lassen sich je zwei wie folgt zusammenfassen:
Der Bruch auf der rechten Seite ist der »partielle Differenzenquotient« der Funktion vx für y = konst. und z = konst. Für Δx gegen 0 wird daraus die partielle Ableitung von vx nach x. Zusammen mit den übrigen vier Summanden ergibt sich dann
1. Gegeben ein Vektorfeld mit dem Feldvektor v(r) = r.
Der Feldvektor ist also radial nach außen gerichtet, seine Länge ist gleich der Länge des Ortsvektors des betreffenden Punktes. Dann ist:
Wir verifizieren an diesem Beispiel den GAUSS-Integralsatz
an einer Kugel vom Radius R um den Ursprung.
In Worten lautet der GAUSS-Integralsatz: Die Ergiebigkeit der Quellen in einem Raumgebiet V ist gleich dem Fluss durch dessen Hüllfläche.
Der Fluß Φ durch ihre Oberfläche ist Φ = 4 π R2R = 4 π R3.
Die Ergiebigkeit S aller innerhalb der Kugel liegenden Quellen ist S = V div r = 3 V = 4 π R3.
2. Es sei v(r) = r/r.
Der Feldvektor ist also radial nach außen gerichtet und hat die konstante Länge 1. Dann ist:
und
Die partielle Ableitung
berechnet man am einfachsten durch implizite Ableitung aus:
Analog findet man
Damit ergibt sich schließlich:
Test:
Der Fluss des Vektors v = r/r durch die Oberfläche einer Kugel um O mit dem Radius R ist
Das Volumen einer Kugelschale vom Radius r und der Dicke dr ist
Die Ergiebigkeit der in der Kugelschale liegenden Quellen ist
Die Ergiebigkeit S der Quellen in einer Kugel vom Radius R ist dann
3. Dieses Beispiel ist von ganz anderer Natur als die vorangegangenen. Hier ist kein Feld vorgegeben, dessen Eigenschaften untersucht werden sollen, sondern eine physikalische Anordnung, ein sehr langer, elektrisch geladener Leiter, dessen Feld gesucht ist.
Wir betrachten ein Leiterelement Δl und eine Kreisscheibe mit Radius R um dieses Leiterelement. Die im Leiterelement vorhandene elektrische Ladung sei Δq, die »Ladungsdichte« ρ also Δq/Δl. Weiter oben haben wir gesehen, dass der von einer Ladung Q erzeugte Fluss Φ des Vektors E = gleich Q/ε0 ist.
Der von der Ladung Δq erzeugte Fluss ΔΦ des Vektors E verlässt die Kreisscheibe nur an deren senkrechter Umrandung, welche die Fläche ΔA = 2πR Δl hat. Da der Fluss auf der Umrandung stets senkrecht steht, gilt für die Feldstärke am Rand
Da E radial nach außen gerichtet ist, ist
Wir berechnen nun noch div E, die außerhalb des Leiters überall null sein muss:
Vorbemerkung: Diese Einleitung ist etwas unkonventionell. Sie versucht, die Begriffe und Zusammenhänge anschaulich werden zu lassen und dem Anfänger dadurch die Chance zu bieten, sie wirklich zu verstehen.
Im Kapitel »Verschiebungsarbeit ...« (Vektoranalysis: Teil II) wurde gezeigt, dass das Linienintegral über das Skalarprodukt v·ds gleich null ist,
- wenn das Integral sich über eine geschlossene Kurve erstreckt und
- wenn ein Potentialfeld vorliegt, d. h. wenn der Vektor v der Gradient eines Skalarfeldes ist.
Letzteres ist jedoch keineswegs immer der Fall, und auch in der Physik gibt es wichtige Felder, die diese Bedingung nicht erfüllen.
Ein Beispiel dafür ist das magnetische Feld eines unendlich langen Leiters (Stromstärke I).
Der Leiter ist von konzentrischen kreisförmigen Feldlinien umgeben; der Vektor H der Feldstärke steht auf dem Radius ρ senkrecht, für seinen Betrag H gilt:
Es lohnt sich, dieses Feld etwas genauer zu betrachten und einige Überlegungen anzustellen.
1. Bewegt man einen Magnetpol (linkes Bild) aus sehr großer (unendlicher) Entfernung radial zu irgendeinem Punkt P hin, so ist dabei keine (positive oder negative) Arbeit zu erbringen. Erfolgt die Bewegung jedoch schräg, so hat der Weg eine Komponente in Richtung des Feldes, und es ist daher Arbeit aufzuwenden. Das Linienintegral
hat keinen bestimmten, vom Weg unabhängigen Wert. Daher kann man dem Punkt P kein bestimmtes Potential zusprechen.
2. Betrachten wir eine Linie, die eine viereckige Fläche umfasst, deren Seiten radial bzw. tangential verlaufen. Man kann sich leicht davon überzeugen, dass das Linienintegral über den geschlossenen Umlauf den Wert null hat. (Wäre dagegen der Betrag der Feldstärke z. B. proportional 1/ρ2, wäre das nicht so.)
3. Bei einem geschlossenen Umlauf, der den Leiter umschlingt (rechtes Bild), hat das Linienintegral dagegen – unabhängig vom Weg - den Wert I (Stromstärke). Dies zeigt, dass das Linienintegral über eine geschlossene Kurve eine besondere Bedeutung haben kann. Darum wollen wir uns genauer mit ihm befassen.
Definition: Unter der Zirkulation Γ eines Vektors v längs einer geschlossenen Kurve K versteht man das Linienintegral des Vektors längs dieser Kurve:
Beispiel: Wie oben gezeigt wurde, ist die Zirkulation des Feldstärkevektors H längs einer Feldlinie des Feldes eines unendlich langen Leiters gleich der Stromstärke I im Leiter. Umfasst dagegen die Kurve K den Leiter nicht, ist die Zirkulation null.
Im Allgemeinen wird die Zirkulation auch von der umlaufenen Fläche A abhängen. Um deren Einfluss auszuschalten, dividiert man die Zirkulation durch die Fläche und betrachtet die Größe Γ/A.
Beispiel: Bei dem oben beschriebenen Feld ist, wenn man eine Feldlinie vom Radius ρ umläuft,
Der Quotient Γ/A nimmt also mit abnehmendem Radius ρ immer mehr zu und erreicht am Umfang des Leiters (Radius a) den Wert der Stromdichte j = I/π a2 im Leiter.
Betrachtet man die Zirkulation ΔΓ des Feldvektors v längs der Umrandung einer kleinen Fläche ΔA und denkt sich diese dann auf einen Punkt P schrumpfend, dann wird der Quotient ΔΓ/ΔA dabei im Allgemeinen einem Grenzwert zustreben. Diesen Grenzwert nenne ich den Wirbel w des Vektors v in P:
Beispiel: Lässt man bei dem oben betrachteten Feld den Radius ρ des Kreises gegen null gehen, also auf den Leitermittelpunkt hin schrumpfen, so wird für den Stromfaden von punktförmigem Querschnitt der Wirbel w = j.
Es soll nun der Wirbel in einem Punkt einer gegebenen Fläche berechnet werden.
In einem Raumgebiet, in dem durch eine Funktion v = v(r) ein Vektorfeld definiert ist, befinde sich ein Flächenstück.
Durch einen Punkt P nahe an der Fläche legen wir ein kleines Koordinatensystem, dessen Achsen parallel zu den entsprechenden Achsen des großen Koordinatensystems sind. Die Ebenen des kleinen Koordinatensystems schneiden aus der Fläche ein kleines Flächenstück heraus, das wir durch ein ebenes Dreieck annähern.
Die Abschnitte auf den Achsen seien 2Δx, 2Δy, 2Δz.
Die Seitenmitten des Dreiecks sind dann
und seine Seitenvektoren
Zur (zunächst) angenäherten Berechnung des Linienintegrals über die drei Seiten multiplizieren wir den jeweiligen Seitenvektor skalar mit dem Wert, den der Feldvektor v in der Seitenmitte hat:
Für die Vektoren vA, vB und vC gilt:
Der Index PA bedeutet:
bei der Richtungsableitung dv/ds dass diese an der Stelle P und in der Richtung PA zu bilden ist,
bei Δs, dass damit die Strecke Δs = PA gemeint ist.
Zur Berechnung werden die drei »Ungefährgleichungen« in ihre Komponenten zerlegt:
Dabei sind – wie auch im Folgenden – alle partiellen Ableitungen an der Stelle P zu bilden.
Ferner ist:
Analoges gilt für die Komponenten von vB und vC.
Wenn man die zusammengehörigen Komponentengleichungen wieder zu einer Vektorgleichung zusammenfasst, erhält man:
Damit ergibt sich:
Die Summe Σ dieser drei Skalarprodukte ist
Dabei sind alle partiellen Ableitungen im Punkt P zu bilden.
Der erste Summand ist null, da die Summe der Seitenvektoren des Dreiecks null ist.
Die letzten drei Summanden können interpretiert werden als das Skalarprodukt aus einem Vektor V und einem Vektor Δ W:
Der erste Vektor erhält wegen seiner besonderen Bedeutung einen eigenen Namen: Rotation (von) v (geschrieben: rot v).
Die Komponenten des zweiten Vektors sind die Projektionen der Fläche ΔA in die Koordinatenebenen:
ΔWx = ΔAx,
ΔWy = ΔAy,
ΔWz = ΔAz.
Das heißt: Der Vektor ΔW ist identisch mit dem Flächenvektor ΔA.
Dieser Vektor kann auch geschrieben werden als
Δ A = ΔAn
wobei n der Normaleneinheitsvektor der Fläche Δ A ist.
Folglich ist
und
Lässt man nun den Punkt P unbeschränkt an die Fläche heranrücken, dann gehen ΔA und die Summe gegen null. Für den Grenzwert des Quotienten gilt:
Das heißt: Der Wirbel w des Feldvektors v im Punkt P einer Fläche ist gleich der Projektion des Vektors rot v an dieser Stelle auf die Flächennormale.
Beachten Sie: Der Vektor rot v ist nur eine Funktion des Ortsvektors r(x, y, z), während der Wirbel auch von der Richtung der Fläche abhängt, für die der Wirbel berechnet wird. Der Wirbel an einer bestimmten Stelle ist maximal, wenn die Fläche auf rot v senkrecht steht; er ist null, wenn der Vektor rot v in der Tangentialebene der Fläche liegt.
Wir betrachten ein beliebiges Flächenstück A mit dem Rand C.
Wir zerlegen die Fläche – wie in der Abbildung angedeutet – in kleine Flächenstücke ΔAi mit den nach außen gerichteten Flächenvektoren ΔAi. Die einzelnen Flächenstücke sollen alle im gleichen Sinn wie der Rand C umlaufen werden. Dabei werden alle Seiten – bis auf die, die auf dem Rand C liegen – zweimal durchlaufen. Die Richtungen der beiden Durchläufe aber sind gegensinnig.
Für jedes einzelne Flächenstück mit seinen vier Seiten gilt dann:
Summiert man über alle Flächenstücke ΔAi, so fallen alle die Summanden heraus, deren Δsi,k nicht auf dem Rand C liegt. Die übrig bleibenden werden nun mit dem Index j versehen:
Für alle ΔA gegen null (wobei natürlich auch alle Δs gegen null gehen) wird daraus
wobei die Gestalt der Fläche völlig beliebig ist. Dies ist der Integralsatz von STOKES.
1. Ein Vektorfeld habe konzentrische, kreisförmige Feldlinien um die Z-Achse. Der Größenwert v des Feldvektors sei proportional 1/ρ (ρ = Abstand des betrachteten Punktes von der Z-Achse). Gesucht die Gleichungen v = v(r) und v = v(x, y, z) sowie rot v.
Der Feldvektor v hat die Richtung des Vektors k x r und soll den Größenwert v = c/ρ haben. Ferner ist
und daher
Mit
ergibt sich
Ein Beispiel für ein solches Feld ist das magnetische Feld (Feldstärke H) eines unendlich langen Leiters. (Siehe »Einleitung – Zirkulation und Wirbel eines Vektors«)
Nach den MAXWELL-Gleichungen ist
Da die Stromdichte außerhalb des Leiters überall null ist, muss dort auch rot H = 0 sein.
Wir wenden nun zur Berechnung des Größenwerts H der Feldstärke den Integralsatz von STOKES auf eine Kreisfläche vom Radius ρ an, die mit dem Leiter konzentrisch ist und auf ihm senkrecht steht:
Nun ist einerseits
wobei q der Leiterquerschnitt und I die Stromstärke im Leiter ist. Andererseits ist
Also ist
2. Gesucht sind die Eigenschaften des Feldvektors
(Siehe dazu die Abbildung bei Beispiel 1.)
Der Vektor v steht auf k und r senkrecht und hat die konstante Länge c.
Ferner ist
und daher
Damit wird
und schließlich
3. Eine ebene Scheibe mit der Flächennormalen n rotiere mit der Winkelgeschwindigkeit ω um einen ihrer Punkte P mit n als Drehachse. Gesucht der Geschwindigkeitsvektor v eines ihrer Punkte im Abstand R von P.
Es ist v = ωR und
Einfacher ist es jedoch, ein neues Koordinatensystem einzuführen, dessen Achsen parallel zu denen des alten sind und dessen Ursprung in P liegt. Dann ist:
Der Hamiltonsche Differential-Operator Nabla (Nabla-Operator)
Bei der Berechnung von Gradient, Divergenz und Rotation werden an einer skalaren Funktion bzw. an einer Vektorfunktion bestimmte Rechenoperationen vorgenommen. Diese weisen bei aller Verschiedenheit gewisse formale Ähnlichkeiten auf: Immer werden die partiellen Ableitungen der Funktion bzw. der Vektorkomponenten gebildet.
1. Gradient:
Es werden die partiellen Ableitungen der Funktion U(x, y, z) berechnet und diese werden als Komponenten eines Vektors benutzt:
Dieser Vektor kann formal aufgefasst werden als das Produkt aus und einem »symbolischen Vektor«
und der Funktion U:
2. Divergenz:
Die Komponenten eines Vektors v = (vxvyvz) werden partiell »nach ihrem Index« abgeleitet und daraus durch Addition eine skalare Funktion gebildet:
Diese Summe kann formal aufgefasst werden als das Skalarprodukt aus Vektor Nabla und dem Vektor
:
3. Rotation:
Die Komponenten eines Vektors v werden partiell »nach den beiden anderen Indices« abgeleitet und daraus Differenzen nach Art eines Vektorprodukts gebildet. Diese werden schließlich als Komponenten eines neuen Vektors benutzt. In der Schreibweise als symbolische Determinante wird dies am anschaulichsten:
Dieser Vektor wiederum kann formal aufgefasst werden als das Vektorprodukt aus dem symbolischen Vektor Nabla und dem Vektor v:
Der symbolische Vektor Nabla wird »Hamiltonscher Differentialoperator« oder Nabla-Operator genannt. (Der Name beruht auf der Ähnlichkeit des Symbols mit der antiken Harfe Nabla.) Also:
In diesem Anhang stelle ich zunächst die Beweise der elementaren Rechengesetze für die Differentialoperatoren Gradient, Divergenz und Rotation zusammen. Dabei gehe ich auch die Anwendung des Differentialoperators Nabla ein. Danach beweise ich die Rechengesetze für die Kombinationen der Differentialoperatoren. Dabei setze ich die Rechengesetze der Analysis (skalerer Funktionen) als bekannt (und bewiesen) voraus.
Zur deutlichen und auffälligen Unterscheidung verwende ich dabei für skalare Ortsfunktionen die Buchstaben f = f(x, y, z), g = g(x, y, z) usw., für vektorielle Ortsfunktionen die Buchstaben v = v(x, y, z), w = w(x, y, z) usw.
Elementare Rechengesetze für die Differentialoperatoren