|
Lexikon auf Ihrer Homepage |
|
Lexikon als Lesezeichen hinzufügen |
Das Oberflächenintegral oder Flächenintegral ist eine Verallgemeinerung des Integralbegriffes auf ebenen oder gekrümmten Flächen. Integrationsgebiet <math>\mathcal F</math> ist also nicht ein eindimensionales Intervall, sondern eine zweidimensionale Menge im dreidimensionalen Raum <math> \mathbb R^3 </math>. Für eine allgemeinere Darstellung im <math> \mathbb R^n </math> mit <math> n \geq 2 </math> siehe: Integration auf Mannigfaltigkeiten.
Es wird generell zwischen einem skalaren und einem vektoriellen Oberflächenintegral unterschieden, je nach Form des Integranden und des sogenannten Oberflächenelements. Sie lauten
Inhaltsverzeichnis |
Bei der Integration über Flächen treten Parametrisierungen der Fläche an die Stelle der Integrationsvariable und Oberflächenelemente an die Stelle der infinitesimalen (unendlich kleinen) Intervallbreite <math> \mathrm d x </math>.
Als zweidimensionale Menge lässt sich eine Oberfläche als Funktion von zwei Variablen darstellen (parametrisieren). Ist <math> B \subset \mathbb R ^2 </math> eine Menge, deren Rand keine doppelten Punkte enthält, stetig differenzierbar, nicht unendlich lang und ferner <math> \varphi </math> eine Abbildung von <math> B </math> in den <math> \mathbb R^3 </math> ist, so sagt man, <math> \varphi </math> ist Parametrisierung der Fläche <math> \mathcal F </math>, wenn <math> \mathcal F = \varphi(B) </math> ist. An dieser Stelle sei darauf hingewiesen, dass ein Großteil der Schwierigkeiten im Umgang mit Oberflächenintegralen mit der Parametrisierung zusammenhängt. Es ist a priori nicht klar, dass unterschiedliche Parametrisierungen den gleichen Wert für das Integral erzeugen. Ein Koordinatenwechsel für Oberflächenintegrale ist nicht trivial und ist mithin Motivation für die Verwendung von Differentialformen.
Allgemein lässt sich eine Fläche im <math>\mathbb{R}^{3}</math> mit zwei Parametern <math>u</math> und <math>v</math> in folgender Form darstellen:
x\left( u,v \right) \\ y\left( u,v \right) \\ z\left( u,v \right) \\
\end{matrix} \right)</math>
Auf der Fläche <math>\vec{\varphi} \left( u,v \right)</math> bilden die Kurvenscharen <math>u=\mathrm{const}</math> bzw. <math>v=\mathrm{const}</math> die Koordinatenlinien. Diese überziehen die Fläche mit einem Koordinatennetz, wobei durch jeden Punkt zwei Koordinatenlinien verlaufen. Somit hat jeder Punkt auf der Fläche eindeutige Koordinaten <math>\left( u_0, v_0 \right)</math>.
Die Oberfläche einer Kugel mit Radius <math> R </math> lässt sich wie folgt parametrisieren: <math> B </math> ist das Rechteck <math> [0, \pi] \times [0, 2\pi[ </math> und
Man rechnet leicht nach, dass diese Parametrisierung die Kugelgleichung <math> x^2 + y^2 + z^2 = R^2 </math> erfüllt (siehe auch Kugelkoordinaten). <math>u</math> ist hier der Polarwinkel (meist <math>\vartheta\,</math> oder <math>\theta\,</math>) und <math>v</math> der Azimutwinkel (meist <math>\varphi\,</math> oder <math>\phi\,</math> bezeichnet).
Ist <math>f:B\to \mathbb{R},\ \left( x,y \right)\mapsto f\left( x,y \right)</math> eine Funktion und die Fläche in der Form <math>z=f(x,y)</math> angegeben so sind <math>x</math> und <math>y</math> die beiden Parameter; die Parametrisierung der Fläche sieht also wie folgt aus:
x \\ y \\ f\left( x,y \right) \\
\end{matrix} \right)</math>
Wenn im eindimensionalen Fall das <math> \mathrm{d}x </math> die Breite eines unendlich kleinen Intervalls darstellt, so liegt es nahe, es im zweidimensionalen Fall durch die Fläche eines unendlich kleinen Flächenstückes <math>\mathrm d \sigma</math> zu ersetzen. Durch die im vorhergehenden Abschnitt beschriebene Parametrisierung kann man an jeden Punkt der Oberfläche zwei Tangenten legen (siehe auch: Krummlinige Koordinaten): Einmal die Tangente, die entsteht, wenn man <math> v </math> konstant lässt und <math> u </math> minimal variiert, und einmal mit vertauschten Variablen. Das heißt also zwei Tangenten an die beiden Koordinatenlinien im betrachteten Punkt <math>\left( u_0, v_0 \right)</math>. Diese Tangenten lassen sich durch zwei infinitesimale Tangentenvektoren ausdrücken (sei <math>\vec{\varphi}\left( u,v \right)</math> die parametrische Form der Fläche):
Im Folgenden wird die kompakte Schreibweise für die partiellen Ableitungen verwendet:
Sind diese Tangenten in keinem Punkt der Fläche parallel, so spricht man von einer regulären Parametrisierung. Das Kreuzprodukt der Tangentenvektoren ist dann ein Vektor, dessen Länge ungleich Null ist.
Die beiden Tangentenvektoren liegen in der Tangentialebene der Fläche am betrachteten Punkt. Der Flächeninhalt des von beiden Tangentenvektoren aufgespannten Parallelogramms entspricht nun gerade dem Betrag ihres Kreuzproduktes.
Ist nun <math> \vec{\varphi}(u,v) </math> eine reguläre Parametrisierung der Oberfläche, so definiert man:
Gemäß den Eigenschaften des Kreuzprodukts steht das vektorielle Oberflächenelement senkrecht auf der Fläche, sein Betrag entspricht gerade der Größe des infinitesimalen Flächenstücks.
In der oben vorgestellten Form ist das vektorielle Oberflächenelement nicht wohldefiniert, da seine Richtung davon abhängt ob man <math> \vec{\varphi}_u \times \vec{\varphi}_v </math> oder <math> \vec{\varphi}_v \times \vec{\varphi}_u = - \left( \vec{\varphi}_u \times \vec{\varphi}_v \right)</math> berechnet. Die beiden Möglichkeiten sind antiparallel zueinander. Betrachtet man geschlossene Oberflächen, vereinbart man meist, dass das nach außen weisende vektorielle Oberflächenelement zu verwenden ist.
Die Oberfläche der Kugel mit Radius R kann, wie oben gezeigt, durch den Polarwinkel <math>u</math> und den Azimutwinkel <math>v</math> parametrisiert werden. Das Flächenelement ergibt sich aus folgender Rechnung:
& \vec{\varphi}=R\left( \begin{matrix}
\sin u\ \cos v \\
\sin u\ \sin v \\
\cos u \\
\end{matrix} \right),\quad \vec{\varphi}_{u}=R\left( \begin{matrix}
\cos u\ \cos v \\ \cos u\ \sin v \\ -\sin u \\
\end{matrix} \right),\quad \vec{\varphi}_{v}=R\left( \begin{matrix}
-\sin u\ \sin v \\ \sin u\ \cos v \\ 0 \\
\end{matrix} \right), \\
& \pm \left( \vec{\varphi}_{u}\times \vec{\varphi}_{v} \right)=\pm R^{2}\sin u\left( \begin{matrix}
\sin u\ \cos v \\
\sin u\ \sin v \\
\cos u \\
\end{matrix} \right), \quad \left| \left| \pm \left( \vec{\varphi}_{u}\times \vec{\varphi}_{v} \right) \right| \right|=R^{2}\sin u,\\
& \hat{n}=\pm \left( \begin{matrix}
\sin u\ \cos v \\
\sin u\ \sin v \\
\cos u \\
\end{matrix} \right), \quad \mathrm{d}\vec{\sigma }=\hat{n}\, \mathrm{d}\sigma =\hat{n}\ R^{2}\sin u\ \mathrm{d}u \,\mathrm{d}v \\ \end{align}</math>
Beim Normalenvektor sind zwei Lösungen möglich (<math>\pm </math>), abhängig von der Reihenfolge von <math>\vec{\varphi}_{u}</math> und <math>\vec{\varphi}_{v}</math> im Kreuzprodukt. Typischerweise wählt man hier die positive Lösung, bei der <math>\hat n</math> von der konvexen Kugeloberfläche weg zeigt (sog. „äußere Normale“).
Ist die Fläche in der Form <math>z=f(x,y)</math> angegeben, so drückt man das Flächenelement durch die Differentiale der Koordinaten <math>x</math>, <math>y</math> aus.
y\\ f(x,y)\end{matrix}\right),\quad\vec{\varphi}_{x}=\left(\begin{matrix}1\\ 0\\ f_{x}\end{matrix}\right),\quad\vec{\varphi}_{y}=\left(\begin{matrix}0\\ 1\\ f_{y}\end{matrix}\right)</math>
-f_{y}\\ 1\end{matrix}\right),\quad\left|\left|\pm\left(\vec{\varphi}_{x}\times\vec{\varphi}_{y}\right)\right|\right|=\sqrt{f_{x}^{2}+f_{y}^{2}+1}\ ,\quad \hat{n}=\pm\frac{1}{\sqrt{f_{x}^{2}+f_{y}^{2}+1}}\left(\begin{matrix}-f_{x}\\ -f_{y}\\ 1\end{matrix}\right)</math>
Somit sind Flächenelement und vektorielles Flächenelement gleich:
-f_{y}\\ 1\end{pmatrix}\ \mathrm{d}x\,\mathrm{d}y</math>
Wir gehen im Folgenden davon aus, dass eine Fläche <math>A</math> mit ihrem Flächenelement <math>\mathrm{d}A</math> und zugehörigem Normalenvektor <math>\hat{n}_A</math> bekannt ist. Z. B.
Für eine weitere Fläche <math>\mathcal F</math> mit Normalenvektor <math>\hat{n}_{\mathcal{F}}</math> soll das Flächenelement <math>\mathrm{d}\sigma</math> ermittelt werden. Die Fläche ist etwa durch <math>g(x,y,z)=0</math> gegeben und somit der Normalenvektor gleich <math>\hat{n}_{\mathcal{F}}=\nabla g/||\nabla g||</math>.
Wir projizieren nun <math>\mathcal F</math> entlang von <math>\hat{n}_{A}</math> auf <math>A</math>. Dann lassen sich die Flächenelemente mittels <math>\mathrm{d}A=\mathrm{d}\vec{A}\cdot\hat{n}_{A}=|\mathrm{d}\vec{\sigma}\cdot\hat{n}_{A}|=\mathrm{d}\sigma\,|\hat{n}_{\mathcal{F}}\cdot\hat{n}_{A}|</math> für <math>\hat{n}_{\mathcal{F}}\cdot\hat{n}_{A}\neq 0</math> verknüpfen:
Dabei darf jede Gerade entlang der Normalenvektoren <math>\hat{n}_{A}</math> die Fläche <math>\mathcal F</math> nur einmal schneiden. Sonst muss man die Fläche <math>\mathcal F</math> aufteilen in kleinere Flächen <math>\mathcal F_1,\,\mathcal F_2\,\ldots</math>, deren Projektion dann eindeutig ist, oder eine andere Grundfläche <math>A</math> wählen.
Das vektorielle Flächenelement ist:
Sei eine Fläche <math>\mathcal F</math> der Form <math>z=f(x,y)</math> gegeben, so gilt <math>g(x,y,z)=z-f(x,y)</math> und damit:
-f_{y}\\ 1\end{pmatrix}\ ,\quad||\nabla g||=\sqrt{f_{x}^{2}+f_{y}^{2}+1}\ ,\quad\hat{n}_{\mathcal{F}}=\frac{\nabla g}{||\nabla g||}=\frac{1}{\sqrt{f_{x}^{2}+f_{y}^{2}+1}}\begin{pmatrix}-f_{x}\\ -f_{y}\\ 1\end{pmatrix}</math>
Diese Fläche wird nun in die xy-Ebene projiziert mit <math>\mathrm{d}A=\mathrm{d}x\mathrm{d}y</math> und <math>\hat{n}_A=\hat{e}_z</math>; dabei ist
-f_{y}\\ 1\end{pmatrix}\mathrm{d}x\mathrm{d}y</math>
Gesucht ist das Flächenelement eines Rotationskörper um die z-Achse mit <math>\rho=f(z)</math>, also <math>g(\rho,\varphi,z)=\rho-f(z)</math>.
Durch Projektion auf die Mantelfläche eines Kreiszylinders mit Radius <math>\rho=f(z)</math> erhält man das Flächenelement:
Mit den Parametrisierungen und den Oberflächenelementen kann man nun die Oberflächenintegrale definieren. Diese mehrdimensionalen Integrale sind Lebesgue-Integrale, können aber in den meisten Anwendungsfällen als mehrfache Riemann-Integrale berechnet werden.
Das skalare Oberflächenintegral einer skalaren Funktion <math> f: \mathbb R^3 \rightarrow \mathbb R </math> über eine Oberfläche <math> \mathcal F </math> mit regulärer Parametrisierung <math> \varphi : B \rightarrow \mathbb R^3 </math> mit <math> B \subset \mathbb R ^2 </math> ist definiert als
Setzt man beispielsweise <math> f(\vec{x}) = 1 </math>, so ist das skalare Oberflächenintegral einfach der Flächeninhalt der Oberfläche.
Das vektorielle Oberflächenintegral einer vektorwertigen Funktion <math> f: \mathbb R^3 \rightarrow \mathbb R^3 </math> über eine Oberfläche <math> \mathcal F </math> mit regulärer Parametrisierung <math> \varphi : B \rightarrow \mathbb R^3 </math> mit <math> B \subset \mathbb R ^2 </math> ist definiert als
Eine anschauliche Vorstellung dieses Integrals geschieht über den Fluss <math>\Phi</math> eines Vektorfeldes <math>\vec{f}</math> durch die Fläche <math> \mathcal F </math>: Die Größe <math> \vec{f} \cdot \mathrm d \vec{\sigma} </math> gibt an, welchen Beitrag zum Gesamtfluss <math>\Phi_{\mathcal F}(\vec f)</math> der infinitesimal-kleine Oberflächen-Vektor <math> \mathrm d \vec{\sigma}=\hat{n}\,\mathrm d\sigma </math> liefert; nämlich wie viel von <math>\vec{f}</math> durch das Oberflächenstück <math>\mathrm d {\sigma}</math> fließt. Der Fluss ist maximal, wenn das Vektorfeld <math>\vec{f}</math> parallel zur Flächennormale <math>\hat{n}</math> steht, und null, wenn <math>\vec{f}</math> senkrecht zu <math>\hat{n}</math> steht, also tangential zur Oberfläche ist - dann "fließt" <math>\vec{f}</math> entlang der Oberfläche, aber nicht durch sie hindurch.
