Helmholtz-Theorem
Das Helmholtz-Theorem, auch Helmholtz-Zerlegung, Stokes-Helmholtz-Zerlegung[1] oder Fundamentalsatz der Vektoranalysis besagt, dass bestimmte differenzierbare Vektorfelder als Summe eines rotationsfreien (wirbelfreien) Gradientenfelds und eines divergenzfreien (quellenfreien) Rotationsfelds geschrieben werden können.
Geschichte
Die Helmholtz-Zerlegung in drei Dimensionen wurde erstmals 1849[2] von George Gabriel Stokes für eine Theorie der Beugung beschrieben, Hermann von Helmholtz veröffentlichte 1858[3] sein Papier über die hydrodynamischen Grundgleichungen,[4] das zu seiner Forschung zu den Helmholtzschen Wirbelsätzen gehört.
Die Zerlegung hat sich zu einem wichtigen Werkzeug für viele Probleme der theoretischen Physik entwickelt,[4][5] aber auch Anwendungen in der Animation, Computervision sowie Robotik gefunden.[6] Dabei wurde die Helmholtz-Zerlegung auf höher-dimensionale Räume erweitert und als Helmholtz-Hodge-Zerlegung unter Nutzung von Differentialgeometrie und Tensorrechnung auch auf riemannschen Mannigfaltigkeiten angewandt.[4][5][6][7]
Definitionen
Das Helmholtz-Theorem besagt, dass es möglich ist, ein auf einem Gebiet definiertes, differenzierbares Vektorfeld als Superposition eines rotationsfreien (wirbelfreien) Gradientfelds für ein Skalarpotential und eines divergenzfreien (quellenfreien) Rotationsfelds darzustellen, so dass:
Diese Zerlegung ist allerdings nicht eindeutig.
Eine alternative Definition auf Basis von Räumen lautet: Für ein Gebiet wird der Raum der divergenzfreien Funktionen genannt, wobei und die -Norm bezeichnet. Die Zerlegung
mit wird Helmholtz-Zerlegung genannt, insofern die Zerlegung existiert. In diesem Fall gibt es eine Projektion mit , die sog. Helmholtz-Projektion.
Die Zerlegung existiert auf jeden Fall, falls der Halbraum, ein beschränktes Gebiet mit -Rand oder ein Außenraum mit -Rand ist. Für existiert die Zerlegung für beliebige Gebiete mit -Rand.[8] Hat einen -Rand, gilt , wobei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{\nu}} die äußere Normale ist. Für andere Fälle gibt es verschiedene Verfahren, Skalar- und Rotationsfeld zu bestimmen.
Verfahren
Im drei-dimensionalen Raum
Im drei-dimensionalen Raum kann ein divergenzfreies Feld als Rotation eines Vektorpotentials Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{A}(\vec{x})} dargestellt werden:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{r}(\vec{x}) = \operatorname{rot}(\vec{A}(\vec{x})).}
Daraus folgt:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \operatorname{rot}(\vec{g}(\vec{x})) = -\operatorname{rot}(\operatorname{grad}(\Phi(\vec{x})))\equiv 0}
und
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \operatorname{div}(\vec{r}(\vec{x})) = \operatorname{div}(\operatorname{rot}(\vec{A}(\vec{x})))\equiv 0} .
Es ist also möglich das Vektorfeld Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{f}(\vec{x})} durch Superposition (Addition) zweier unterschiedlicher Potentiale Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Phi(\vec{x})} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{A}(\vec{x})} auszudrücken.
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{f}(\vec{x}) = \vec{g}(\vec{x}) + \vec{r}(\vec{x}) = -\operatorname{grad}(\Phi(\vec{x})) + \operatorname{rot}(\vec{A}(\vec{x}))}
Die beiden einander ergänzenden Potentiale lassen sich durch die folgenden Integrale aus dem Feld Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{f}(\vec{x})} gewinnen. Es handelt sich hierbei um die Faltung der Divergenz bzw. Rotation des Vektorfelds, wobei als Integralkern die Fundamentallösung der Laplace-Gleichung genutzt wird.[9][10]
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Phi(\vec{x}) = \frac{1}{4\pi}\int_V \frac{\operatorname{div}(\vec{f}(\vec{x}'))}{|\vec{x}-\vec{x}'|}\mathrm{d}^3r'}
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{A}(\vec{x}) = \frac{1}{4\pi}\int_V \frac{\operatorname{rot}(\vec{f}(\vec{x}'))}{|\vec{x}-\vec{x}'|} \mathrm{d}^3r'}
Hierbei ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V} das die Felder enthaltende Volumen. Bei Feldern, für die Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V} unbeschränkt ist, beispielsweise der gesamte drei-dimensionale Raum, so ist die mathematische Voraussetzung für die Verwendung der Faltungsintegrale, dass Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \operatorname{div}(\vec{f}(\vec{x}))} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \operatorname{rot}(\vec{f}(\vec{x}))} für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle r \to \infty} schneller als Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \tfrac{1}{|\vec{x}|}} gegen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 0} geht, also Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \lim\nolimits_{|\vec{x}| \to \infty} \vec{f}(\vec{x}) |\vec{x}| = 0} .[11] Ansonsten divergieren die obigen Integrale, lassen sich also nicht mehr berechnen. Daher wird die Helmholtz-Zerlegung in Lehrbüchern häufig überhaupt nur für Testfunktionen definiert, die diese Eigenschaft erfüllen.[10] Otto Blumenthal zeigte allerdings bereits 1905, dass mit einem veränderten Integrationskern eine Integration für alle Felder möglich ist, die schneller als eine Potenzfunktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle |\vec{x}|^{-d}} mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle d > 0} abfallen. Hierfür ersetzt man in den Integrationsgleichungen den Kern Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K(x, x')} durch Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K'(\vec{x}, \vec{x}') = K(\vec{x}, \vec{x}') - K(0, \vec{x}')} .[12] Mit noch komplexeren Integrationskernen kann sogar für divergierende, allerdings nicht schneller als polynomial ansteigende, Funktionen eine numerische Lösung berechnet werden.[13][14]
Im n-dimensionalen Raum
Die Verallgemeinerung auf Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle n} Dimensionen kann nicht mit einem Vektorpotential erfolgen, da der Rotationsoperator und das Kreuzprodukt nur in drei Dimensionen definiert sind. Das Skalarpotential ist hingegen identisch definiert wie in drei Dimensionen als
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G(\vec{x}) = \int_{\mathbb{R}^n} \operatorname{div}(\vec{f}(\vec{x}')) K(\vec{x}, \vec{x}') dx' = \int_{\mathbb{R}^n} \sum_i \frac{\partial f_i}{\partial x_i}(\vec{x}') K(\vec{x}, \vec{x}') dx',}
wobei als Integrationskern Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K(\vec{x}, \vec{x}')} wieder die Fundamentallösung der Laplace-Gleichung im n-dimensionalen Raum eingesetzt wird:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K(\vec{x}, \vec{x}') = \begin{cases} \frac{1}{2\pi} \log{ | \vec{x}-\vec{x}' | } & n=2, \\ \frac{1}{n(2-n)V_n} | \vec{x}-\vec{x}' | ^{2-n} & \text{sonst}, \end{cases} }
mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V_n = \pi^\frac{n}{2} / \Gamma\big(\tfrac{n}{2}+1\big)} dem Volumen des n-dimensionalen Einheitsballs und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Gamma(\vec{x})} der Gamma-Funktion.
Für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle n = 3} entspricht Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle V_n} gerade Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{4 \pi}{3}} , wodurch sich derselbe Vorfaktor wie oben ergibt. Das Rotationspotential ist eine antisymmetrische Matrix mit den Elementen:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle R_{ij}(\vec{x}) = \int_{\mathbb{R}^n} \left( \frac{\partial f_i}{\partial x_j}(\vec{x}') - \frac{\partial f_j}{\partial x_i}(\vec{x}') \right) K(\vec{x}, \vec{x}') dx'. }
Oberhalb der Diagonale stehen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \textstyle\binom{n}{2}} Einträge, die an der Diagonale gespiegelt erneut auftreten, allerdings mit negativem Vorzeichen. Im drei-dimensionalen Fall entsprechen die Matrixelemente gerade die Komponenten des Vektorpotentials Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle A = [A_1, A_2, A_3] = [R_{23}, R_{31}, R_{12}]} . Ein solches Matrix-Potential lässt sich allerdings nur im dreidimensionalen Fall als Vektor schreiben, weil Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \textstyle\binom{n}{2} = n} nur für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle n = 3} gilt.
Wie im drei-dimensionalen Fall ist das Gradientenfeld als Gradient des Skalarpotentials Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G} definiert. Das Rotationsfeld ist hingegen im allgemeinen Fall definiert als:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{r}(\vec{x}) = \left[ \sum\nolimits_k \partial_{x_k} R_{ik}(\vec{x}); {1 \leq i \leq n} \right],}
also als Zeilen-Divergenz der Matrix. Im drei-dimensionalen Raum ist dies äquivalent zur Rotation des Vektorpotentials.[7][13]
Genauso wie in drei Dimensionen muss die Divergenz des Vektorfelds sowie der Term Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \tfrac{\partial f_i}{\partial x_j} - \tfrac{\partial f_j}{\partial x_i}} für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle |\vec{x}| \to \infty} schneller als Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \tfrac{1}{|\vec{x}|^2}} gegen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 0} gehen, oder bei durch Ersetzen des Integrationskerns Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K(\vec{x}, \vec{x}')} durch Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K'(\vec{x}, \vec{x}') = K(\vec{x}, \vec{x}') - K(0, \vec{x}')} schneller als Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle |\vec{x}|^{-d}} mit d > 0 abfallen.[14]
Für einige Vektorfelder, die auch nicht im Unendlichen gegen Null gehen müssen, lassen sich mit Methoden basierend auf der partiellen Integration und der Cauchy-Formel für mehrfache Integration[15] die Rotations- und Skalarpotentiale analytisch berechnen, wie bei multivariaten Polynomen, Sinus-, Cosinus- und Exponentialfunktionen.[7]
Helmholtz-Hodge-Zerlegung auf riemannschen Mannigfaltigkeiten
Im Rahmen der Hodge-Theorie und der De-Rham-Kohomologie verallgemeinert die Helmholtz-Hodge-Zerlegung das Helmholtz-Theorem von Vektorfeldern zu Differentialformen auf riemannschen Mannigfaltigkeiten.[6][16][17]
Eindeutigkeit und Eichung
Grundsätzlich ist die Helmholtz-Zerlegung nicht eindeutig bestimmt. Addiert man zum Skalarpotential Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle G(\vec{x})} eine harmonische Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle H(\vec{x})} , die also Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Delta H(\vec{x}) = 0} erfüllt, so ist
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \begin{align} \vec{g}'(\vec{x}) &= \operatorname{grad} (G(\vec{x}) + H(\vec{x})) = \vec{g}(\vec{x}) + \operatorname{grad} H(\vec{x}),\\ \vec{r}'(\vec{x}) &= \vec{r}(\vec{x}) - \operatorname{grad} H(\vec{x}) \end{align}}
ebenfalls eine Helmholtz-Zerlegung von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{f}(\vec{x})} . Für Vektorfelder Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{f}} , die im unendlichen gegen null abfallen, ist eine plausible Wahl, dass Skalar- und Rotations- bzw. Vektorpotential dies auch tun. Da Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle H(\vec{x}) = 0} die einzige harmonische Funktion mit dieser Eigenschaft ist, was sich aus einer Abwandlung des Satzes von Liouville folgern lässt, ist so die Eindeutigkeit des Gradienten- und Rotationsfelds sichergestellt.[18]
Während das ursprüngliche Vektorfeld an jedem Punkt von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Omega} durch Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle n} Komponenten zu beschreiben ist, sind im drei-dimensionalen Fall für das skalare und das Vektorpotential zusammen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle n+1} Komponenten nötig. Dadurch entsteht eine Redundanz, weil die Wahl der Potentiale nicht eindeutig ist, die physikalische Beschreibung invariant gegenüber der Wahl einer Eichung. Die Eichtheorie befasst sich damit, wie verschiedene Potentiale mittels Eichtransformationen berechnet werden können, bekannte Beispiele aus der Physik sind die Coulomb-Eichung und die Lorenz-Eichung. Die Redundanz lässt sich auch beseitigen, indem der quellfreie Anteil des Vektorfeldes der toroidal-poloidalen Zerlegung unterworfen wird, wodurch letztlich insgesamt drei Skalarpotentiale zur Beschreibung ausreichen.
Anwendungen
Elektrodynamik
Das Helmholtz-Theorem ist besonders in der Elektrodynamik von Interesse, da sich mit seiner Hilfe die Maxwell-Gleichungen im Potentialbild schreiben und einfacher lösen lassen. Mit der Helmholtz-Zerlegung lässt sich beweisen, dass bei gegebener elektrischer Stromdichte und Ladungsdichte das elektrische Feld und die magnetische Flussdichte bestimmen lassen. Sie sind eindeutig, wenn die Dichten im unendlichen verschwinden und man dasselbe für die Potentiale annimmt.[10]
Fluiddynamik
In der Fluiddynamik, speziell der Lösbarkeitstheorie der Navier-Stokes-Gleichungen, spielt die Helmholtz-Projektion eine wichtige Rolle. Wird die Helmholtz-Projektion auf die linearisierte inkompressiblen Navier-Stokes-Gleichungen angewandt, erhält man die Stokes-Gleichung. Diese ist nur noch von der Geschwindigkeit der Teilchen in der Strömung abhängig, jedoch nicht mehr vom statischen Druck, wodurch die Gleichung auf eine Unbekannte reduziert werden konnte. Beide Gleichungen, die Stokes- und die linearisierte Gleichung, sind jedoch äquivalent. Der Operator Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle P\Delta} wird Stokes-Operator genannt.[19]
Theorie dynamischer Systeme
In der Theorie dynamischer Systeme können mittels der Helmholtz-Zerlegung „Quasipotentiale“ bestimmt werden sowie in manchen Fällen Lyapunov-Funktionen berechnet werden.[20][21][22]
Für einige dynamische Systeme wie das auf Edward N. Lorenz (1963[23]) zurückgehende Lorenz-System, ein vereinfachtes Modell für atmosphärische Konvektion, lässt sich die Helmholtz-Zerlegung analytisch berechnen:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot \vec{x} = \vec{f}(\vec{x}) = \big[a (x_2-x_1), x_1 (b-x_3)-x_2, x_1 x_2-c x_3 \big].}
Die Helmholtz-Zerlegung von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \vec{f}(\vec{x})} ist, mit dem Skalarpotential Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \Phi(\vec{x}) = - \tfrac{a}{2} x_1^2 - \tfrac{1}{2} x_2^2 - \tfrac{c}{2} x_3^2} gegeben als:
- Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \begin{align} \vec{g}(\vec{x}) &= \big[-a x_1, -x_2, -c x_3 \big], \\ \vec{r}(\vec{x}) &= \big[+ a x_2, b x_1 - x_1 x_3, x_1 x_2 \big]. \end{align}}
Das quadratische Skalarpotential sorgt für eine Bewegung in Richtung des Koordinatenursprungs, was für den stabilen Fixpunkt für einige Parameterbereich verantwortlich ist. Für andere Parameter sorgt das Rotationsfeld dafür, dass ein seltsamer Attraktor entsteht, wodurch das Modell einen Schmetterlingseffekt zeigt.[7][24]
Computeranimation und Robotik
Die Helmholtz-Zerlegung wird auch im Bereich der Computertechnik verwendet. Dazu gehört die Robotik, die Bildrekonstruktion aber auch die Computeranimation, wo die Zerlegung für eine realistische Visualisierung von Fluiden oder Vektorfeldern eingesetzt wird.[6]
Einzelnachweise
- ↑ Tribikram Kundu: Ultrasonic and Electromagnetic NDE for Structure and Material Characterization. CRC Press, 2012, ISBN 1-4398-3663-9, S. 37 (eingeschränkte Vorschau in der Google-Buchsuche).
- ↑ George Gabriel Stokes: On the Dynamical Theory of Diffraction. In: Transactions of the Cambridge Philosophical Society 9, 1849, S. 1–62. doi:10.1017/cbo9780511702259.015.
- ↑ Hermann von Helmholtz: Über Integrale der hydrodynamischen Gleichungen, welche den Wirbelbewegungen entsprechen. In: Journal für die reine und angewandte Mathematik 55, 1858, S. 25–55, doi:10.1515/crll.1858.55.25 (Volltext auf sub.uni-goettingen.de und digizeitschriften.de).
- ↑ a b c Alp Kustepeli: On the Helmholtz Theorem and Its Generalization for Multi-Layers. In: Electromagnetics 36.3, 2016, S. 135–148, doi:10.1080/02726343.2016.1149755.
- ↑ a b Wolfgang Sprössig: On Helmholtz decompositions and their generalizations – An overview. In: Mathematical Methods in the Applied Sciences 33.4, 2009, S. 374–383, doi:10.1002/mma.1212.
- ↑ a b c d Harsh Bhatia, Gregory Norgard, Valerio Pascucci, Peer-Timo Bremer: The Helmholtz-Hodge Decomposition – A Survey. In: IEEE Transactions on Visualization and Computer Graphics 19.8, 2013, S. 1386–1404, doi:10.1109/tvcg.2012.316.
- ↑ a b c d Erhard Glötzl, Oliver Richters: Analytical Helmholtz Decomposition and Potential Functions for many n-dimensional unbounded vector fields. 2021, arxiv:2102.09556, doi:10.48550/arXiv.2102.09556.
- ↑ G. P. Galdi, An introduction to the mathematical theory of the Navier-Stokes equations. Vol. I, Springer Tracts in Natural Philosophy, vol. 38, Springer-Verlag, New York, 1994, ISBN 0-387-94172-X.
- ↑ Morton E. Gurtin: On Helmholtz’s theorem and the completeness of the Papkovich-Neuber stress functions for infinite domains. In: Archive for Rational Mechanics and Analysis 9.1, 1962, S. 225–233, doi:10.1007/BF00253346.
- ↑ a b c Dietmar Petrascheck: The Helmholtz decomposition revisited. In: European Journal of Physics 37.1, 2015, Artikel 015201, doi:10.1088/0143-0807/37/1/015201.
- ↑ R. Douglas Gregory: Helmholtz's Theorem when the domain is Infinite and when the field has singular points. In: The Quarterly Journal of Mechanics and Applied Mathematics 49.3, 1996, S. 439–450, doi:10.1093/qjmam/49.3.439.
- ↑ Otto Blumenthal: Über die Zerlegung unendlicher Vektorfelder. In: Mathematische Annalen 61.2, 1905, S. 235–250, doi:10.1007/BF01457564.
- ↑ a b Erhard Glötzl, Oliver Richters: Helmholtz Decomposition and Rotation Potentials in n-dimensional Cartesian Coordinates. 2020, arxiv:2012.13157, doi:10.48550/arXiv.2012.13157.
- ↑ a b Ton Tran-Cong: On Helmholtz’s Decomposition Theorem and Poissons’s Equation with an Infinite Domain. In: Quarterly of Applied Mathematics 51.1, 1993, S. 23–35, JSTOR 43637902.
- ↑ Augustin-Louis Cauchy: Trente-Cinquième Leçon. In: Résumé des leçons données à l’École royale polytechnique sur le calcul infinitésimal. Imprimerie Royale, Paris 1823, S. 133–140 (gallica.bnf.fr).
- ↑ Hersh Bhatia, Valerio Pascucci, Peer-Timo Bremer: The Natural Helmholtz-Hodge Decomposition for Open-Boundary Flow Analysis. In: IEEE Transactions on Visualization and Computer Graphics 20.11, Nov. 2014, S. 1566–1578, Nov. 2014, doi:10.1109/TVCG.2014.2312012.
- ↑ Frank W. Warner: The Hodge Theorem. In: Foundations of Differentiable Manifolds and Lie Groups. (= Graduate Texts in Mathematics 94). Springer, New York 1983, doi:10.1007/978-1-4757-1799-0_6.
- ↑ Sheldon Axler, Paul Bourdon, Wade Ramey: Bounded Harmonic Functions. In: Harmonic Function Theory (= Graduate Texts in Mathematics 137). Springer, New York 1992, S. 31–44, doi:10.1007/0-387-21527-1_2.
- ↑ Alexandre J. Chorin, Jerrold E. Marsden: A Mathematical Introduction to Fluid Mechanics (= Texts in Applied Mathematics 4). Springer US, New York 1990, doi:10.1007/978-1-4684-0364-0.
- ↑ Tomoharu Suda: Construction of Lyapunov functions using Helmholtz–Hodge decomposition. In: Discrete & Continuous Dynamical Systems – A 39.5, 2019, S. 2437–2454, doi:10.3934/dcds.2019103.
- ↑ Tomoharu Suda: Application of Helmholtz–Hodge decomposition to the study of certain vector fields. In: Journal of Physics A: Mathematical and Theoretical 53.37, 2020, S. 375703. doi:10.1088/1751-8121/aba657.
- ↑ Joseph Xu Zhou, M. D. S. Aliyu, Erik Aurell, Sui Huang: Quasi-potential landscape in complex multi-stable systems. In: Journal of The Royal Society Interface 9.77, 2012, S. 3539–3553, doi:10.1098/rsif.2012.0434.
- ↑ Edward N. Lorenz: Deterministic Nonperiodic Flow. In: Journal of the Atmospheric Sciences 20.2, 1963, S. 130–141, doi:10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2.
- ↑ Heinz-Otto Peitgen, Hartmut Jürgens, Dietmar Saupe: Strange Attractors: The Locus of Chaos. In: Chaos and Fractals. Springer, New York, S. 655–768. doi:10.1007/978-1-4757-4740-9_13.