Rotation eines Vektorfeldes

{{#if: beschreibt den mathematischen Differentialoperator „Rotation“. Zu der Drehbewegung eines Körpers siehe Rotation (Physik), zu der geometrischen Abbildung siehe Drehung.

 | Vorlage:Hinweisbaustein 
 | {{#ifeq: 0 | 0 |}}

}}

Als Rotation oder Rotor<ref>{{#invoke:Vorlage:Literatur|f}}</ref><ref>{{#invoke:Vorlage:Literatur|f}}</ref> bezeichnet man in der Vektoranalysis, einem Teilgebiet der Mathematik, einen Differentialoperator, der einem Vektorfeld im dreidimensionalen euklidischen Raum mit Hilfe der Differentiation ein neues Vektorfeld zuordnet.

Die Rotation eines Strömungsfeldes gibt für jeden Ort das Doppelte der Winkelgeschwindigkeit an, mit der sich ein mitschwimmender Körper dreht („rotiert“). Dieser Zusammenhang ist namensgebend.

Datei:Uniform curl.svg

Das Geschwindigkeitsfeld einer rotierenden Scheibe besitzt eine konstante Rotation parallel zur Drehachse
<math>\left( \omega < 0 \right)</math>

Es muss sich aber nicht immer um ein Geschwindigkeitsfeld und eine Drehbewegung handeln; beispielsweise betrifft das Induktionsgesetz die Rotation des elektrischen Feldes.

Ein Vektorfeld, dessen Rotation in einem Gebiet überall gleich null ist, nennt man wirbelfrei oder, insbesondere bei Kraftfeldern, konservativ. Ist das Gebiet einfach zusammenhängend, so ist das Vektorfeld genau dann der Gradient einer Funktion, wenn die Rotation des Vektorfeldes im betrachteten Gebiet gleich null ist.

Die Divergenz der Rotation eines Vektorfeldes ist gleich null. Umgekehrt ist in einfach zusammenhängenden Gebieten ein Feld, dessen Divergenz gleich null ist, die Rotation eines anderen Vektorfeldes.

Beispiele:

Das Vektorfeld, das an jedem Ort die Windrichtung und -geschwindigkeit eines Wirbelsturms angibt, hat in der Umgebung des Auges (der Rotationsachse) eine von null verschiedene Rotation.
Das Vektorfeld <math>\vec{v}(x, y, z) = \omega\cdot(x\,\hat{e}_y-y\,\hat{e}_x)\,,</math> das an jedem Punkt einer rotierenden Scheibe die Geschwindigkeit angibt, hat an jedem Punkt dieselbe von null verschiedene Rotation. Die Rotation beträgt das Zweifache der Winkelgeschwindigkeit, <math>\operatorname{rot}\,\vec{v}(x,y,z) = 2\,\omega \,\hat{e}_z\,.</math> Siehe Abbildung
Das Kraftfeld, das an jedem Punkt die Gravitationskraft der Sonne auf ein Testteilchen angibt, ist wirbelfrei. Das Kraftfeld ist der negative Gradient der potentiellen Energie des Teilchens.

Definition der Rotation

Definition in kartesischen Koordinaten

Seien <math>(x,y,z)</math> die kartesischen Koordinaten des dreidimensionalen euklidischen Raumes und <math>\hat{e}_x</math>, <math>\hat{e}_y</math> und <math>\hat{e}_z</math> die auf Einheitslänge normierten, zueinander senkrechten Basisvektoren, die an jedem Punkt in Richtung der zunehmenden Koordinaten zeigen.

Die Rotation eines dreidimensionalen, differenzierbaren Vektorfeldes

ist das dreidimensionale Vektorfeld

<math>\operatorname{rot}\,\vec F(x,y,z) =

\left (\frac{\partial F_z}{\partial y} - \frac{\partial F_y}{\partial z}\right )\hat{e}_x

\left (\frac{\partial F_x}{\partial z} - \frac{\partial F_z}{\partial x}\right )\hat{e}_y

\left (\frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}\right )\hat{e}_z \,.</math> Man kann <math>\operatorname{rot}\, \vec F</math> wie das Kreuzprodukt als formale Determinante einer Matrix auffassen, deren erste Spalte die kartesischen Basisvektoren enthält, die zweite die partiellen Ableitungen nach den kartesischen Koordinaten und die dritte die zu differenzierenden Komponentenfunktionen

<math>\operatorname{rot}\,\vec F =\operatorname{det}\,

\begin{pmatrix} \hat{e}_x & \frac{\partial}{\partial x} & F_x\\ \hat{e}_y & \frac{\partial}{\partial y} & F_y\\ \hat{e}_z & \frac{\partial}{\partial z} & F_z \end{pmatrix}\,=\operatorname{det}\, \begin{pmatrix} \hat{e}_x & \hat{e}_y & \hat{e}_z\\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z}\\ F_x & F_y & F_z \end{pmatrix}\,. </math> Allerdings sind hier die verschiedenen Spalten nicht Vektoren desselben Vektorraumes.

Gibt man die Vektoren als Spaltenvektoren ihrer kartesischen Komponenten an, dann ist <math>\operatorname{rot}\,\vec F</math> das formale Kreuzprodukt des Spaltenvektors der partiellen Ableitungen nach den kartesischen Koordinaten, des Nabla-Operators <math>\nabla</math>, mit dem Spaltenvektor der kartesischen Komponentenfunktionen

<math>\operatorname{rot}\,\vec F(x,y,z) = \nabla\times\vec F

=\sum_{i=1}^3\hat e_i\times\frac{\partial\vec F}{\partial x_i} =\begin{pmatrix}

 \frac{\partial}{\partial x} \\
 \frac{\partial}{\partial y} \\
 \frac{\partial}{\partial z}

\end{pmatrix} \times \begin{pmatrix}

 F_x\\
 F_y\\
 F_z

\end{pmatrix} = \begin{pmatrix}

 \frac{\partial F_z}{\partial y} - \frac{\partial F_y}{\partial z} \\
 \frac{\partial F_x}{\partial z} - \frac{\partial F_z}{\partial x} \\
 \frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}

\end{pmatrix} </math>

wo die Koordinaten nach dem üblichen Schema x → 1, y → 2 und z → 3 durchnummeriert wurden.

Koordinatenunabhängige Definition mit dem Nabla-Operator

Der Nabla-Operator ist auch in anderen Koordinatensystemen definiert und so kann mit ihm die Rotation koordinatenunabhängig durch

<math>\mathrm{rot}\ \vec F:=\nabla\times\vec F</math>

definiert werden. Mit dem Nabla-Operator können auch der Gradient- sowie die Divergenz eines Vektorfeldes dargestellt und Produktregeln hergeleitet werden.

Definition in Kugelkoordinaten

Schreibt man das Vektorfeld in Kugelkoordinaten <math>(r, \theta, \varphi)</math> als Linearkombination

<math>\vec F(r,\theta,\varphi)=F_r(r,\theta,\varphi)\, \hat{e}_r + F_{\theta}(r,\theta,\varphi)\,\hat{e}_\theta + F_\varphi(r,\theta,\varphi)\,\hat{e}_\varphi </math>

der auf Einheitslänge normierten Vektoren

<math>

\begin{align} \hat{e}_r &=\begin{pmatrix} \sin(\vartheta)\cos(\varphi)\\ \sin(\vartheta)\sin(\varphi)\\ \cos(\vartheta) \end{pmatrix} = \frac{1}{\sqrt{x^2+y^2+z^2}} \begin{pmatrix} x\\y\\z \end{pmatrix} \,,\\ \hat{e}_\theta &=\begin{pmatrix} \cos(\vartheta)\cos(\varphi)\\ \cos(\vartheta)\sin(\varphi)\\ -\sin(\vartheta) \end{pmatrix} = \frac{1}{\sqrt{(x^2+y^2+z^2)(x^2+y^2)}} \begin{pmatrix} z\,x\\z\,y\\-x^2-y^2 \end{pmatrix}\,,\\ \hat{e}_\varphi &=\begin{pmatrix} -\sin(\varphi)\\ \cos(\varphi)\\ 0 \end{pmatrix} = \frac{1}{\sqrt{x^2+y^2}} \begin{pmatrix} -y\\x\\0 \end{pmatrix} \,, \end{align} </math> die an jedem Punkt in Richtung zunehmender <math>r,\theta,\varphi</math>-Koordinaten zeigen, so ist die Rotation

<math>

\begin{align} \operatorname{rot}\,\vec F = \,

 &\frac{1}{r \sin \theta} \left[ \frac{\partial}{\partial \theta} \left( F_\varphi \sin \theta \right) - \frac{\partial F_\theta}{\partial \varphi }\right]\hat{e}_r

+

\left [ \frac{1}{r \sin \theta}\frac{\partial F_r}{\partial \varphi} - \frac{1}{r}\frac{\partial}{\partial r} \left( r F_\varphi \right)\right]\hat{e}_\theta \,+\,
 \frac{1}{r} \left[ \frac{\partial}{\partial r} \left( r F_\theta \right) - \frac{\partial F_r}{\partial \theta} \right] \hat{e}_\varphi

\,. \end{align} </math>

Definition in Zylinderkoordinaten

Gibt man das Vektorfeld in Zylinderkoordinaten <math>(r, \varphi, z)</math> als Linearkombination

<math>\vec F(r,\varphi,z)=F_r(r,\varphi,z)\, \hat{e}_r

+F_\varphi(r,\varphi,z)\,\hat{e}_\varphi +F_{z}(r,\varphi,z)\,\hat{e}_z</math>

der Vektoren

<math>

\begin{align} \hat{e}_r &=\begin{pmatrix}\cos(\varphi)\\ \sin(\varphi)\\ 0\end{pmatrix} = \frac{1}{\sqrt{x^2+y^2}} \begin{pmatrix} x\\y\\0 \end{pmatrix} \,,\\ \hat{e}_\varphi&=\begin{pmatrix}-\sin(\varphi)\\ \cos(\varphi)\\0\end{pmatrix} = \frac{1}{\sqrt{x^2+y^2}} \begin{pmatrix} -y\\x\\0 \end{pmatrix}\,,\\ \hat{e}_z &= \begin{pmatrix} 0\\0\\1 \end{pmatrix} \,, \end{align} </math> an, die auf Einheitslänge normiert an jedem Punkt in Richtung zunehmender <math>r,\varphi,z</math>-Koordinaten zeigen, so ist die Rotation

<math>

\begin{align} \operatorname{rot}\,\vec F = \left[ \frac 1 r \frac {\partial F_z}{\partial \varphi} - \frac{\partial F_\varphi}{\partial z} \right]\hat{e}_r +\left [ \frac{\partial F_r}{\partial z} - \frac{\partial F_z}{\partial r} \right ]\hat{e}_\varphi \,+ \, \frac 1 r \left[ \frac \partial {\partial r} \left( r \cdot F_\varphi \right) - \frac{\partial F_r}{\partial \varphi} \right]\hat{e}_z\,. \end{align}</math>

Rotation in zwei Dimensionen

Ein Vektorfeld im zweidimensionalen euklidischen Raum kann als Vektorfeld

in drei Dimensionen aufgefasst werden, das nicht von der dritten Koordinate abhängt und dessen dritte Komponente verschwindet. Seine Rotation ist kein Vektorfeld dieser Art, sondern besteht gemäß

<math>\operatorname{rot}\,\vec F(x,y,z) =

\left (\frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}\right )\hat{e}_z </math> aus einer Komponente, die senkrecht zum Vektorfeld in drei Dimensionen ist. Definiert man in zwei Dimensionen die Rotation als den Differentialoperator

<math>\operatorname{rot}\colon \vec F \mapsto \frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}\,,</math>

dann ist das Ergebnis ein Skalarfeld und kein Vektorfeld.

Eigenschaften

Koordinatenfreie Darstellung der Rotation als Volumenableitung

Mit Hilfe des Satzes von Stokes kann die Rotation, ähnlich wie die Divergenz (Quellendichte), als Volumenableitung dargestellt werden. Diese Darstellung hat den Vorteil, dass sie koordinatenunabhängig ist. Aus diesem Grund wird die Rotation im Bereich der Ingenieurwissenschaften oftmals direkt so definiert.

Ist <math>\mathcal{V}</math> ein Raumgebiet mit stückweise glattem Rand <math>\partial \mathcal{V}</math> und dem Volumeninhalt <math>V</math>, dann kann die Rotation des Vektorfelds <math>\vec{F} \colon \mathcal{V} \to \R^3</math> im Punkt <math>p \in \mathcal{V}</math> mittels der Volumenableitung durch

<math>

\mathrm{rot}\,\vec F(p)=\lim_{V\rightarrow 0}\frac{\oint_{\partial \mathcal{V}} \hat n\times\vec F\,\mathrm{d}A}{V} </math>

berechnet werden. Dabei bezeichnet <math>\hat n\mathrm{d}A=\mathrm{d}\vec{A}</math> das äußere vektorielle Flächenelement von <math>\partial \mathcal{V},</math> wobei <math>\hat n</math> der nach außen zeigende Normaleneinheitsvektor und <math>\mathrm{d}A</math> das skalare Flächenelement ist. Zur Grenzwertbildung wird das Raumgebiet <math>\mathcal{V}</math> auf den Punkt <math>p</math> zusammengezogen, sodass sein Inhalt <math>V</math> gegen null geht, siehe auch #Integralsatz von Stokes weiter unten.<ref>{{#invoke:Vorlage:Literatur|f}}</ref>

Ersetzt man <math>\vec F</math> durch eine Strömungsgeschwindigkeit, erscheint die Rotation als Wirbeldichte. Ähnlich gebildete Synonyme existieren auch für die Divergenz (Quellendichte) und den Gradienten (Kraftdichte). Die Koordinatendarstellungen des vorigen Abschnitts ergeben sich aus der Volumenableitung, wenn man das jeweilige Volumenelement als Raumgebiet <math>\mathcal{V}</math> wählt.

Axialvektorfeld

Die Rotation eines Vektorfeldes ist ein Pseudovektorfeld. Ein Vektorfeld geht bei Spiegelung am Ursprung in sein Negatives am gespiegelten Ort über, die Rotation des Vektorfeldes ändert bei dieser Spiegelung ihr Vorzeichen nicht,

<math>\begin{align}

\vec F^\prime(\vec{x}) & = - \vec F(-\vec{x})\,, \\ \bigl(\operatorname{rot}\,\vec F^\prime\bigr)(\vec{x}) & = \bigl(\operatorname{rot}\,\vec F\bigr)(-\vec{x})\,. \end{align}</math>

Rechenregeln

Die Rotation ist linear. Für alle Konstanten <math>c\in\R</math> und differenzierbaren Vektorfelder <math>\vec{F}</math> und <math>\vec{G}</math> gilt

<math>

\operatorname{rot}\,(c \,\vec{F}+\vec G)
 = c\,\operatorname{rot}\,\vec{F} + \operatorname{rot}\,\vec{G}\,.</math>

Die Rotation eines Vektorfeldes verschwindet genau dann, wenn es lokal ein Gradientenfeld ist und die Divergenz eines Vektorfeldes verschwindet genau dann, wenn es lokal die Rotation eines anderen Feldes ist:

<math>

\operatorname{rot~grad}\,f= 0\,,\ \operatorname{div~rot}\,\vec F= 0 </math>

Für differenzierbare Funktionen <math>f\,</math> und Vektorfelder <math>\vec{F}</math> und <math>\vec{G}</math> gelten die Produktregeln

<math>\begin{align}

\operatorname{rot}\,(f\,\vec{F}) =& f\,\operatorname{rot}\,\vec{F} + (\operatorname{grad}\,f)\,\times \vec{F} \\ \operatorname{rot}\,(\vec{F}\times\vec{G}) =& \left(\vec{G}\cdot\nabla\right)\vec{F} - \left(\vec{F}\cdot\nabla\right)\vec{G} + \vec{F}\,(\nabla\cdot\vec{G}) - \vec{G}\,(\nabla\cdot\vec{F}) \\=& \left(\operatorname{grad}\vec{F}\right)\cdot\vec{G} - \left(\operatorname{grad}\vec{G}\right)\cdot\vec{F} + \vec{F}\,(\operatorname{div}\,\vec{G}) - \vec{G}\,(\operatorname{div}\,\vec{F}) \,.\end{align}</math>

Darin ist <math>\nabla</math> der Nabla-Operator und in der letzten Formel bildet grad den Vektorgradient. Für die zweifache Anwendung der Rotation gilt

<math>

\operatorname{rot}\,\operatorname{rot}\,\vec{F}= \operatorname{grad\,div}\,\vec{F} - \operatorname{div\,grad}\,\vec{F}= \operatorname{grad\,div}\,\vec{F} - \Delta\,\vec{F} </math>

wo <math>\Delta</math> der Laplace-Operator ist. Für einen Vektor <math> \vec{v} </math>, der von einem Skalar <math>s</math> abhängt, und dieser in 3D vom Ort, gilt die Kettenregel

<math>

\operatorname{rot}\, \vec{v} ( s ) = \operatorname{grad}\,s \times \frac{\mathrm{d}\vec{v}}{\mathrm{d}s}. </math>

Anwendungen

Zusammenhang zur Winkelgeschwindigkeit

Bei der Drehung eines starren Körpers um die <math>z</math>-Achse mit konstanter Winkelgeschwindigkeit <math>\omega</math> wächst der Drehwinkel <math>\varphi</math> gleichmäßig mit der Zeit an, <math>\varphi = \omega\,t\,,</math> und jeder Punkt durchläuft eine Bahn

<math>\begin{pmatrix}

x(t)\\y(t)\\z(t) \end{pmatrix} = \begin{pmatrix} \cos(\omega\,t)\,x(0) -\sin(\omega\,t)\,y(0)\\ \sin(\omega\,t)\,x(0) +\cos(\omega\,t)\,y(0)\\ z(0) \end{pmatrix}\,. </math> Die Geschwindigkeit beträgt

<math>

\frac{\mathrm d}{\mathrm d\,t} \begin{pmatrix} x(t)\\y(t)\\z(t) \end{pmatrix} =\omega\, \begin{pmatrix} -\sin(\omega\,t)\,x(0) -\cos(\omega\,t)\,y(0)\\ \ \cos(\omega\,t)\,x(0) -\sin(\omega\,t)\,y(0)\\ 0 \end{pmatrix}= \omega\, \begin{pmatrix} -y(t)\\x(t)\\0 \end{pmatrix} \,. </math> Das Geschwindigkeitsfeld einer starren Drehung um die <math>z</math>-Achse ist also, wie oben im Beispiel angegeben,

<math>\vec{v}(x,y,z) = \omega\,(-y\,\hat{e}_x + x \, \hat{e}_y)\,.</math>

Seine Rotation ist die doppelte Winkelgeschwindigkeit

<math>\operatorname{rot}\,\vec{v} = 2\,\omega\,\hat{e}_z\,.</math>

Veranschaulichung durch Drehmoment

In einem Flächenkraftdichte-Feld<ref><templatestyles src="Webarchiv/styles.css" />{{#if:20160810043839

      | {{#ifeq: 20160810043839 | *
    | Vorlage:Webarchiv/Wartung/Stern{{#if: Formelsammlung Mechanik | {{#invoke:WLink|getEscapedTitle|Formelsammlung Mechanik}} | {{#invoke:Webarchiv|getdomain|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf}} }} (Archivversionen)
    | {{#iferror: {{#time: j. F Y|20160810043839}}
         | {{#if:  || }}Vorlage:Webarchiv/Wartung/DatumDer Wert des Parameters {{#if: wayback | wayback | Datum }} muss ein gültiger Zeitstempel der Form YYYYMMDDHHMMSS sein!
         | {{#if: Formelsammlung Mechanik | {{#invoke:WLink|getEscapedTitle|Formelsammlung Mechanik}} | {{#invoke:Webarchiv|getdomain|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf}} }} {{#ifeq:  | [] | [ | ( }}Memento{{#if: {{#if:  | {{{archiv-bot}}} |  }} |  des Vorlage:Referrer }} vom {{#time: j. F Y|20160810043839}} im Internet Archive{{#if:  | ;  }}{{#ifeq:  | [] | ] | ) }}
      }}
  }}
      | {{#if:
          | {{#iferror: {{#time: j. F Y|{{{webciteID}}}}}
    | {{#switch: {{#invoke:Str|len|{{{webciteID}}}}}
       | 16= {{#if: Formelsammlung Mechanik | {{#invoke:WLink|getEscapedTitle|Formelsammlung Mechanik}} | {{#invoke:Webarchiv|getdomain|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf}} }} {{#ifeq:  | [] | [ | ( }}Memento{{#if: {{#if:  | {{{archiv-bot}}} |  }} |  des Vorlage:Referrer }} vom {{#time: j. F Y| 19700101000000 + {{#expr: floor {{#expr: {{#invoke:Str|sub|{{{webciteID}}}|1|10}}/86400}} }} days}} auf WebCite{{#if:  | ;  }}{{#ifeq:  | [] | ] | ) }}
       | 9 = {{#if: Formelsammlung Mechanik | {{#invoke:WLink|getEscapedTitle|Formelsammlung Mechanik}} | {{#invoke:Webarchiv|getdomain|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf}} }} {{#ifeq:  | [] | [ | ( }}Memento{{#if: {{#if:  | {{{archiv-bot}}} |  }} |  des Vorlage:Referrer}} vom {{#time: j. F Y| 19700101000000 + {{#expr: floor {{#expr: {{#invoke:Str|sub|{{#invoke:Expr|base62|{{{webciteID}}}}}|1|10}}/86400}} }} days}} auf WebCite{{#if:  | ;  }}{{#ifeq:  | [] | ] | ) }}
       | #default= Der Wert des Parameters {{#if: webciteID | webciteID | ID }} muss entweder ein Zeitstempel der Form YYYYMMDDHHMMSS oder ein Schüsselwert mit 9 Zeichen oder eine 16-stellige Zahl sein!Vorlage:Webarchiv/Wartung/webcitation{{#if:  || }}
      }}
    | c|{{{webciteID}}}}} {{#if: Formelsammlung Mechanik | {{#invoke:WLink|getEscapedTitle|Formelsammlung Mechanik}} | {{#invoke:Webarchiv|getdomain|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf}} }} (Memento{{#if: {{#if:  | {{{archiv-bot}}} |  }} |  des Vorlage:Referrer}} vom {{#time: j. F Y|{{{webciteID}}}}} auf WebCite{{#if:  | ;  }}{{#ifeq:  | [] | ] | ) }}
  }}
          | {{#if: 
              | Vorlage:Webarchiv/Today
              | {{#if:
                      | Vorlage:Webarchiv/Generisch
                      | {{#if: Formelsammlung Mechanik | {{#invoke:WLink|getEscapedTitle|Formelsammlung Mechanik}} | {{#invoke:Webarchiv|getdomain|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf}} }}  
                 }}}}}}}}{{#if:
    | Vorlage:Webarchiv/archiv-bot
  }}{{#invoke:TemplatePar|check
     |all      = url=
     |opt      = text= wayback= webciteID= archive-is= archive-today= archiv-url= archiv-datum= ()= archiv-bot= format= original=
     |cat      = Wikipedia:Vorlagenfehler/Vorlage:Webarchiv
     |errNS    = 0
     |template = Vorlage:Webarchiv
     |format   = *
     |preview  = 1
  }}{{#ifexpr: {{#if:20160810043839|1|0}}{{#if:|+1}}{{#if:|+1}}{{#if:|+1}}{{#if:|+1}} <> 1
    | {{#if:  || }}Vorlage:Webarchiv/Wartung/Parameter{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: Genau einer der Parameter 'wayback', 'webciteID', 'archive-today', 'archive-is' oder 'archiv-url' muss angegeben werden.|1}}
  }}{{#if: 
    | {{#switch: {{#invoke:Webarchiv|getdomain|{{{archiv-url}}}}}
        | web.archive.org = 
          {{#if:  || }}{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: Im Parameter 'archiv-url' wurde URL von Internet Archive erkannt, bitte Parameter 'wayback' benutzen.|1}} 
        | webcitation.org = 
          {{#if:  || }}{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: Im Parameter 'archiv-url' wurde URL von WebCite erkannt, bitte Parameter 'webciteID' benutzen.|1}} 
        | archive.today |archive.is |archive.ph |archive.fo |archive.li |archive.md |archive.vn = 
          {{#if:  || }}{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: Im Parameter 'archiv-url' wurde URL von archive.today erkannt, bitte Parameter 'archive-today' benutzen.|1}}
      }}{{#if: 
         | {{#iferror: {{#iferror:{{#invoke:Vorlage:FormatDate|Execute}}|Vorlage:FormatDate/Wartung/Error}}
             | {{#if:  || }}Vorlage:Webarchiv/Wartung/Parameter{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: Der Wert des Parameter 'archiv-datum' ist ungültig oder hat ein ungültiges Format.|1}}
          |  }} 
         | {{#if:  || }}Vorlage:Webarchiv/Wartung/Parameter{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: Der Pflichtparameter 'archiv-datum' wurde nicht angegeben.|1}}
      }}
    | {{#if: 
         | {{#if:  || }}Vorlage:Webarchiv/Wartung/Parameter{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: Der Parameter 'archiv-datum' ist nur in Verbindung mit 'archiv-url' angebbar.|1}}
      }}
  }}{{#if:{{#invoke:URLutil|isHostPathResource|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf}}
    || {{#if:  || }} URL/Pfad ungültig.
  }}{{#if: Formelsammlung Mechanik
    | {{#if: {{#invoke:WLink|isBracketedLink|Formelsammlung Mechanik}}
        | {{#if:  || }} Linktext ungültig.
      }}
    | {{#if:  || }}Vorlage:Webarchiv/Wartung/Linktext_fehlt Linktext fehlt.
  }}{{#switch: 
    |addlarchives|addlpages= {{#if:  || }}{{#if: 1 |Vorlage:Webarchiv/Wartung/Parameter}}{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: enWP-Wert im Parameter 'format'.|1}}
  }}{{#ifeq: {{#invoke:Str|find|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf%7Carchiv}} |-1
    || {{#ifeq: {{#invoke:Str|find|{{#invoke:Str|cropleft|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf%7C4}}%7Chttp}} |-1
         || {{#switch: {{#invoke:Webarchiv|getdomain|http://n.ethz.ch/student/nuetzig/zwickers/documents/basis/mech12_norman.pdf }}
              | abendblatt.de | daserste.ndr.de | inarchive.com | webcitation.org = 
              | #default = {{#if:  || }}{{#if: 1 |Vorlage:Webarchiv/Wartung/URL}}{{#invoke:TemplUtl|failure| Fehler bei Vorlage:Webarchiv: Archiv-URL im Parameter 'url' anstatt URL der Originalquelle. Entferne den vor der Original-URL stehenden Mementobestandteil und setze den Archivierungszeitstempel in den Parameter 'wayback', 'webciteID', 'archive.today' oder 'archive-is' ein, sofern nicht bereits befüllt.|1}}
            }} 
       }}
  }}</ref> <math> \vec f</math>, das jedem Körperoberflächenelement mit dem Inhalt <math>\mathrm{d}A</math> unabhängig von seiner Ausrichtung die Kraft <math>\vec f \mathrm{d}A</math> einprägt, erfährt eine Kugel mit dem Radius <math> R</math> (und dem zugehörigen Volumeninhalt <math>V</math>) das Drehmoment

<math>\vec M= R\, V\,\mathrm{rot}\,\vec f.</math>

Vorausgesetzt ist, dass <math> \mathrm{rot}\,\vec f</math> im Bereich der Kugel konstant ist. Die Gleichung folgt aus dem #Integralsatz von Stokes

<math>\int_\mathcal{V}\mathrm{rot}\,\vec f\mathrm{d}V = \oint_\mathcal{A}\hat n\times\vec f\,\mathrm{d}A</math>

mit <math>\textstyle \int_\mathcal{V}\mathrm{rot}\,\vec f\,\mathrm{d}V = V\mathrm{rot}\,\vec f</math> und <math>\textstyle \vec M:=\oint_\mathcal{A}R\hat n\times\vec f\,\mathrm{d}A</math>.

Sätze, in denen die Rotation eine Rolle spielt

Zerlegung in quellen- und wirbelfreien Teil

Zweifach stetig differenzierbare Vektorfelder <math>\vec{v}(\vec r)</math>, die mit ihren Ableitungen für große Abstände hinreichend rasch gegen null gehen, kann man eindeutig in einen wirbelfreien Teil <math>\vec{E}\,,\ \textsf{mit}\ \operatorname{rot}\,\vec{E} = \vec{0}\,,</math> und einen quellenfreien Teil <math>\vec{B}\,,\ \textsf{mit}\ \operatorname{div}\,\vec{B} = 0\,,</math> zerlegen,

<math>\begin{array}{lll}

\vec{v} &= \vec{E} +\vec{B}\,,& \vec{E} = -\operatorname{grad}\,\phi\,,\ \vec{B} = \operatorname{rot}\,\vec{A}\,, \\ \phi (\vec{x}) &= \frac{1}{4\,\pi}\int\!\mathrm{d}^{3}y\,\, \frac{\operatorname{div}\,\vec{v}(\vec{y})}{|\vec{x} -\vec{y}|}\,,\ &\vec{A}(\vec{x}) =\frac{1}{4\,\pi}\int\!{\mathrm d}^{3}y\,\, \frac{\operatorname{rot}\,\vec{v}(\vec{y})}{|\vec{x} -\vec{y}|} \,. \end{array} </math>

Dabei bezeichnen <math>\operatorname{div}</math> und <math>\operatorname{grad}</math> den Divergenz- bzw. Gradient-Operator, wobei die Definition <math> E = -\operatorname{grad}\,\phi\,</math> die in der Physik übliche Konvention ist. Mathematisch ist: <math> E = \operatorname{grad}\,\phi\,</math>

Diese Zerlegung ist Bestandteil des Helmholtz-Theorems.

Integralsatz von Stokes

Datei:Stokesrot.png

Fläche <math>\mathcal{F}</math> mit Berandung <math>\partial\mathcal{F}</math>

{{#if: Satz von Stokes|{{#ifexist:Satz von Stokes|

|{{#if: |{{#ifexist:{{{2}}}|

→ Haupt{{#if:|seite|artikel}}: [[{{{2}}}{{#if: ||{{{titel2}}}}}]]{{#if: |{{#ifexist:{{{3}}}| und [[{{{3}}}{{#if: ||{{{titel3}}}}}]]|}}|}}

|{{#if: |{{#ifexist:{{{3}}}|

→ Haupt{{#if:|seite|artikel}}: [[{{{3}}}{{#if: ||{{{titel3}}}}}]]

|}}|}}|}}|}}|}}|Einbindungsfehler: Die Vorlage Hauptartikel benötigt immer mindestens ein Argument.}}

Das Integral über eine Fläche <math>\mathcal{F}</math> über die Rotation eines Vektorfeldes <math>\vec{A}</math> ist nach dem klassischen Integralsatz von Stokes gleich dem Kurvenintegral über die Randkurve <math>\partial \mathcal{F}</math> über <math>\vec{A}\,,</math>

<math> \iint_{\mathcal{F}}\!\!(\operatorname{rot}\,\vec {A})\cdot\mathrm{d}\vec f= \oint_{\partial \mathcal{F}}\!\!\vec {A}\cdot\mathrm{d}\vec{x}\,.</math>

Durch das Doppelintegral wird links betont, dass man von einer zweidimensionalen Integration ausgeht. Auf der rechten Seite soll das Kreissymbol im Integralzeichen unterstreichen, dass es sich um ein Integral über einen geschlossenen Weg handelt. Die Orientierung entspricht dabei der Drei-Finger-Regel, siehe Abbildung rechts: die folgenden drei Vektoren, nämlich erstens der Vektor <math>\mathrm d\vec f</math> in Richtung der Flächennormalen, zweitens der Vektor <math>\mathrm d\vec{x}</math> in Tangentialrichtung der Kurve und drittens der vom Rand in die Fläche zeigenden Vektor, entsprechen Daumen, Zeigefinger und Mittelfinger der rechten Hand, das heißt, sie bilden ein Rechtssystem. Oft schreibt man <math>\mathrm{d}\vec{f}=\vec{n}\,\mathrm{d}f\,</math> indem man mit dem Normalenvektor <math>\vec n</math> die Richtung der Größe hervorhebt.

Der allgemeinere Satz von Stokes beinhaltet auch das Rotations-Theorem<ref>Altenbach (2012), S. 46.</ref>

<math>\int_V\mathrm{rot}\vec F\,\mathrm{d}V

=\int_A\vec n\times\vec F\,\mathrm{d}A</math>

Darin ist <math>\vec F</math> ein stetig differenzierbares Vektorfeld, <math>\vec n</math> der nach außen gerichtete Normaleneinheitsvektor auf der geschlossenen Oberfläche <math>A</math> des Volumens <math>V</math>. Wenn das Volumen so klein wird, dass die Rotation in ihm näherungsweise konstant wird, folgt hieraus die #Koordinatenfreie Darstellung der Rotation als Volumenableitung.

Rotation von Tensoren zweiter Stufe

Die Rotation von Tensorfeldern zweiter Stufe wird mit der Identität<ref name="hbphys">{{#invoke:Vorlage:Literatur|f}}</ref>

<math>\mathrm{rot}(\mathbf{T})\cdot\vec c

=\mathrm{rot}\left(\mathbf{T}^\top\cdot\vec c\right) \quad\forall\vec c</math>

definiert. Aus ihr ergibt sich

<math>\mathrm{rot}(\mathbf{T})=\nabla\times\left(\mathbf{T}^\top\right)</math>.

In kartesischen Koordinaten <math>x_{1,2,3}</math> bezüglich der Standardbasis ê_1,2,3 schreibt sich die Rotation für das Tensorfeld <math>\textstyle\mathbf{T}=\sum_{i,j=1}^3T_{ij}\hat e_i\otimes\hat e_j</math>:

<math>\mathrm{rot}(\mathbf{T})

=\sum_{i,j,k=1}^3\hat e_k\times\frac{\partial}{\partial x_k} T_{ij}\hat e_j\otimes\hat e_i =\sum_{i,j,k=1}^3\frac{\partial T_{ij}}{\partial x_k} (\hat e_k\times\hat e_j)\otimes\hat e_i </math>

Darin ist ⊗ das dyadische Produkt. Es wird aber auch die transponierte Version <math>\nabla\times\mathbf{T}</math> benutzt<ref>Altenbach (2012), S. 43.</ref>, die hieraus hervorgeht, indem die Komponenten gemäß <math>T_{ij}\rightleftarrows T_{ji}</math> vertauscht werden.

Im Zusammenhang mit Tensoren sind Klammern ein wichtiges Hilfsmittel, um die Reihenfolge der Anwendung und die Argumente der verschiedenen Operatoren klarzustellen, was auf das Ergebnis einen entscheidenden Einfluss hat. Meistens ist beispielsweise

<math>\operatorname{rot}(\mathbf{T}^\top)

=\nabla\times\mathbf{T} \ne \left[\nabla\times\left(\mathbf{T}^\top\right)\right]^\top = \operatorname{rot}(\mathbf{T})^\top </math>

weswegen die Ausdrücke <math>\operatorname{rot}\mathbf{T}^\top</math> und <math>\nabla\times\mathbf{T}^\top</math> mehrdeutig sind.

Symmetrische Tensoren

Wenn der Tensor symmetrisch ist, <math>\textstyle \mathbf{T}=\mathbf{T}^\top=\sum_{i,j=1}^3 T_{ij}\hat e_i\otimes\hat e_j</math> mit <math>T_{ij}=T_{ji}</math>, dann ist seine Rotation spurfrei:

<math>\operatorname{Sp\big(rot}(\mathbf{T})\big)

= \sum_{i,j,k=1}^3\frac{\partial T_{ij}}{\partial x_k} (\hat e_k\times\hat e_j)\cdot\hat e_i = \sum_{i,j,k=1}^3\frac{\partial T_{ij}}{\partial x_k} (\hat e_j\times\hat e_i)\cdot\hat e_k = 0 </math>

denn Terme mit vertauschten Indizes <math>i</math> und <math>j</math> sind gleich groß, besitzen aber umgekehrtes Vorzeichen und heben sich daher in der Summe gegenseitig auf, oder verschwinden bei <math>i=j</math>, siehe auch Spatprodukt.

Ableitungsregeln

Die Produktregel führt im Produkt mit einem Skalar <math>f</math>, Vektoren <math>\vec{f}, \vec{g}</math> und dem Tensor <math>\mathbf{T}</math> auf:

<math>\begin{align}

\mathrm{rot}(\vec{f}\otimes\vec{g}) =& \mathrm{rot}(\vec{g})\otimes\vec{f}-\vec{g}\times\mathrm{grad}(\vec{f})^\top \\ \mathrm{rot}(f\mathbf{T}) =& \mathrm{grad}(f)\times(\mathbf{T}^\top) + f\mathrm{rot}(\mathbf{T}) \\ \mathrm{rot}(\mathbf{T}\cdot\vec{f}) =& \mathrm{rot}(\mathbf{T}^\top)\cdot\vec{f} -\vec{\mathrm i}\left(\mathbf{T}\cdot\mathrm{grad}(\vec{f})\right) \\ \mathrm{rot}(\vec{f}\times\mathbf{T}) =& -\mathrm{rot}(\mathbf{T})\times\vec{f} +\left(\mathbf{T}\#\mathrm{grad}(\vec{f})\right)^\top \end{align}</math>

Darin bildet <math>\vec{\mathrm i}</math> die Vektorinvariante, # das äußere Tensorprodukt und grad den Gradienten. Ist T der Einheitstensor 1, dann liefert das bemerkenswerte Zusammenhänge:

<math>\begin{align}

\mathrm{rot}(f\mathbf1) =& \mathrm{grad}(f)\times\mathbf1 \\ \mathrm{rot}(\vec{f}) =& -\vec{\mathrm i}\left(\mathrm{grad}(\vec{f})\right) \\ \mathrm{rot}(\vec{f}\times\mathbf1) =& \left(\mathbf1\#\mathrm{grad}(\vec{f})\right)^\top = \left(\mathrm{Sp}\big(\mathrm{grad}(\vec{f})\big)\mathbf1 -\mathrm{grad}(\vec{f})^\top\right)^\top \\=& \mathrm{div}(\vec{f})\mathbf1-\mathrm{grad}(\vec{f}) \end{align}</math>

In divergenzfreien Feldern ist also <math>\mathrm{rot}(\vec{f}\times\mathbf1)=-\mathrm{grad}(\vec{f})</math>, was beim Poincaré-Lemma ausgenutzt wird.

Bei der Verknüpfung der Rotation mit anderen Differentialoperatoren entstehen unter Beteiligung eines Tensors teilweise ähnliche Formeln wie sie aus der Vektoranalysis bekannt sind:

<math>\begin{array}{rclcl}

\mathrm{div(rot}(\vec{f})) &=&\nabla\cdot(\nabla\times\vec{f}) &=& 0 \\ \mathrm{rot(grad}(f)) &=&\nabla\times\nabla f &=&\vec0 \\ \mathrm{div\big(rot}(\mathbf{T})^\top\big) &=&\nabla\cdot(\nabla\times\mathbf{T}) &=&\vec{0} \\ \mathrm{rot\big(grad}(\vec{f})\big) &=&\nabla\times(\nabla\otimes\vec{f}) &=&\mathbf{0} \end{array}</math>

<math>\begin{align}

\mathrm{rot\big(rot}(\vec{f})\big) =& \mathrm{grad\big(div}(\vec{f})\big) -\Delta\vec{f} \\ \mathrm{rot\big(rot}(\mathbf{T})^\top\big)^\top =&\mathrm{grad\big(div}(\mathbf{T})\big)-\Delta\mathbf{T} \end{align}</math>

oder mit den Nabla-Operator

<math>\begin{align}

\nabla\times(\nabla\times\vec f) =& \nabla(\nabla\cdot\vec{f})-\Delta\vec{f} \\ \left[\nabla\times\big(\nabla\times(\mathbf{T}^\top)\big)\right]^\top =& \big(\nabla\otimes\nabla\cdot\mathbf{T}^\top\big)^\top -\Delta\mathbf{T} \end{align}</math>

Darin ist Δ = 𝜵² der Laplace-Operator.

Siehe auch

Einzelnachweise

Literatur

{{#invoke:Vorlage:Literatur|f}}
{{#invoke:Vorlage:Literatur|f}}
{{#invoke:Vorlage:Literatur|f}}

Weblinks

Wie „krümme“ ich Nabla und Delta? Herleitung des Nablaoperators für orthonormal krummlinige Koordinaten, auf matheplanet.com.