Vektorgradient

mathematischer Operator
(Weitergeleitet von Gradient eines Vektorfeldes)

Der Gradient eines Vektorfeldes oder kurz Vektorgradient (von lateinisch gradiens ‚schreitend‘[1]) fasst das Gefälle oder den Anstieg der Komponenten eines Vektorfeldes zu einem mathematischen Objekt zusammen. Während mit dem Gradient eines Skalarfeldes das Gefälle oder der Anstieg in einer bestimmten Richtung (sog. Richtungsableitung) als Skalar angegeben wird, stellt die Richtungsableitung mit dem Vektorgradient einen Vektor dar.

Abb. 1: Die Drehung und Streckung von materiellen Linienelementen h (rot) bei einer Deformation wird mit einem Vektorgradient beschrieben

Ein anschauliches Beispiel ist das Vektorfeld der Bewegung der Partikel eines Körpers. Der Deformationsgradient transformiert die Strecke von einem Partikel zu einem benachbarten Partikel des Körpers im undeformierten Zustand – das ist die vorgegebene Richtung h – in die entsprechende Strecke im deformierten Zustand, was die Richtungsableitung in Richtung h ist, siehe Bild. Die Strecke h kann bei der Deformation gedreht und gestreckt werden. Die Richtungsableitung mit maximalem Wert ist hier diejenige Richtung, in der der Körper die größte Dehnung erfährt; in dieser Richtung benachbarte Partikel entfernen sich im Zuge der Verformung am weitesten voneinander, siehe auch #Verformungen und die #Beispiele.

Der Gradient eines Vektorfeldes entsteht aus dem Vektorfeld durch Anwendung des Gradientenoperators grad, der eine Verallgemeinerung der Ableitung in der mehrdimensionalen Analysis ist. Zur besseren Abgrenzung zwischen Operator und Resultat seiner Anwendung bezeichnen manche Quellen[2]:353[3]:112 den Gradient vektorieller Feldgrößen als Vektorgradient.

Der Gradient hat tensorielle Eigenschaften[4]:421: der Gradient eines Skalarfeldes (Tensorfeld nullter Stufe) führt auf ein Gradientenvektorfeld, das ein Tensorfeld erster Stufe ist. Entsprechend ist der Vektorgradient ein Tensorfeld zweiter Stufe; das Ergebnis lässt sich bezüglich einer Orthonormalbasis als Matrix schreiben. Die Komponenten des Vektorgradienten sind die kovarianten Ableitungen der Komponenten des Vektorfeldes in einem Punkt; bei den Basisvektoren sind dies die Christoffelsymbole.

Der Gradient wird zusammen mit anderen Differentialoperatoren wie Divergenz und Rotation in der Tensoranalysis untersucht.

Definition

Bearbeiten

Der Gradient   eines differenzierbaren Vektorfeldes   nähert dieses in der Umgebung eines Punkts   linear an:[5]

  für  

Das Landau-Symbol 𝓞(x) steht für Terme, die langsamer als x wachsen, und   stellt eine lineare Funktion von   dar. Wenn der Gradient existiert, ist er eindeutig und kann aus dem Gâteaux-Differential

 

berechnet werden. In einem euklidischen Vektorraum mit Standardskalarprodukt „·“ ergibt sich der Vektorgradient aus der Anwendung des skalaren Operators  , der mit dem Nabla-Operator 𝜵 gebildet wird:

 

So werden auch Gradienten für Tensorfelder zweiter Stufe oder allgemein Tensorfelder n-ter Stufe definiert.[2]:358[4]:420[6]:43 Durch Gradientenbildung entsteht aus einem Tensorfeld n-ter Stufe ein Tensorfeld der Stufe n+1.

Bei einem Vektorfeld, das ein Tensorfeld erster Stufe ist, ergibt sich als Gradient ein Tensorfeld zweiter Stufe und zwar durch Nutzung des dyadischen Produkts „⊗“:

 

Das hochgestellte bedeutet eine Transponierung. In einigen Quellen[5]:4[7]:23[8]:34 wird

 

und in anderen[6]:43[2]:354

 

definiert, was wegen des nicht kommutativen dyadischen Produkts einen nicht unerheblichen Unterschied ausmacht, der beispielsweise bei der #Produktregel und der Richtungsableitung zu beachten ist.

Hier wird die erstere Form (die ohne Tilde) als Konvention benutzt.

Schreibweisen

Bearbeiten

Die vielfältigen Anwendungen haben zu variantenreichen Schreibweisen geführt.

In der Kontinuumsmechanik ist es üblich, Größen, die sich auf den undeformierten Ausgangszustand eines Körpers beziehen, groß zu schreiben, und solche, die sich auf den deformierten Zustand beziehen, klein. Entsprechend bedeuten GRAD oder Grad Gradienten im undeformierten Körper und grad einen Gradient im deformierten. Andere Notationen mit dem Nabla-Operator benutzen 𝜵X, 𝜵0 für den Operator im undeformierten Körper und 𝜵x, 𝜵t für den im deformierten.

Es wird auch

 

geschrieben.

Geometrische Interpretation

Bearbeiten

Das eingangs aufgeführte Beispiel des Deformationsgradienten soll hier vertieft werden. Dazu sei das Vektorfeld in einer nahen Umgebung eines Punkts als Bewegungsfeld der Partikel einer Gummihaut interpretierbar, was der Fall ist, wenn der Vektorgradient invertierbar ist, es also eine Eins-zu-eins-Beziehung zwischen Raumpunkten und ihren Bildern gibt. Bezeichnen im undeformierten Körper Großbuchstaben die Orte von Partikeln und Kleinbuchstaben ihre Orte im deformierten, dann stellt man folgendes fest.

Auf der Haut wird ein Kreis gezeichnet und wenn man nun die Gummihaut lang zieht, wird der Kreis zu irgendeiner geschlossenen Kurve, siehe Abb. 3. Mit kleiner werdendem Radius des Kreises nähert sich sein Bild jedoch einer Ellipse. Ein Pfeil MP vom Mittelpunkt M zu einem Partikel P auf dem Umfang des Kreises wird zu mp gedehnt und verdreht, wobei p auf der Ellipse liegt. Die Transformation von MP zu mp leistet beim kleinen Kreis die mit dem Vektorgradient gebildete Richtungsableitung. Diese lineare Annäherung des Vektorfeldes mittels der Richtungsableitung in der Umgebung des Punkts ist die definierende Eigenschaft eines Gradienten.

Die betraglich größte Änderung der Positionsdifferenzen mp zu MP tritt in der Richtung auf, in der der Körper die größte Dehnung erfährt. Auf der Gummihaut landet das Partikel P im Kreis auf der Hauptachse der Ellipse (bei p). Die Richtung der größten betraglichen Änderung erhält man hier als Lösung eines Eigenwertproblems, jedoch nicht des Deformationsgradienten, sondern des mit ihm gebildeten Strecktensors, siehe #Zusammenhang mit der Richtungsableitung und #Beispiele.

Markiert man im undeformierten Körper ein Partikel A und ein (infinitesimal) nahe benachbartes B und trägt im deformierten Körper vom Ort a des Partikels A die Richtungsableitung in Richtung AB auf, dann landet man im deformierten Körper am Ort b des Partikels B. Genauso kann man in b die Richtungsableitung in Richtung BC zu einem benachbarten Partikel C auftragen und landet im deformierten Körper an dessen Ort c. Diese Prozedur kann man beliebig oft wiederholen, die Integralrechnung gestattet sogar unendliche Wiederholungen mit infinitesimalen Schritten. So gelangt man von a aus an den Ort p eines beliebigen Partikels P und zwar unabhängig vom eingeschlagenen Weg von A nach P. Diese Wegunabhängigkeit zeichnet Gradientenfelder aus.

Koordinatendarstellung

Bearbeiten

Zwecks kompakter Darstellung bezeichnet im Folgenden ein Index hinter einem Komma die Ableitung nach einer Koordinate:

 

Hier ist die vereinbarte #Konvention   zu beachten.

Kartesische Koordinaten

Bearbeiten

In kartesischen Koordinaten   mit Standardbasis êi lautet der Vektorgradient eines Vektorfeldes   mit Koeffizienten  

 

mit

 .

In drei Dimensionen ist speziell

 

Die Vektorgradienten der Basisvektoren verschwinden hier, da sie ortsunabhängig sind.

Zylinderkoordinaten

Bearbeiten

In Zylinderkoordinaten mit radialer Koordinate ρ, Azimut φ und Höhe z über der ρφ-Ebene lauten die Basisvektoren mit dem Sinus und Cosinus

 

und der Vektorgradient einer Vektorfunktion  :

 

mit

 

Die Terme in der zweiten Zeile oben sind Ergebnis der Vektorgradienten der Basisvektoren:

 

Kugelkoordinaten

Bearbeiten

In Kugelkoordinaten mit Abstand r vom Ursprung, Zenitwinkel ϑ und Azimut φ lauten die Basisvektoren mit dem Sinus und Cosinus

 

und der Vektorgradient einer Vektorfunktion  :

 

mit

 

dem Tangens tan und dem Einheitstensor 1 = êr ⊗ êr + êϑ ⊗ êϑ + êφ ⊗ êφ.

Die Terme in der zweiten Zeile oben sind Ergebnis der Vektorgradienten der Basisvektoren:

 

Allgemein krummlinige Koordinaten

Bearbeiten

In krummlinigen Koordinaten   lauten die ko- und kontravarianten Basisvektoren

 

Das Symbol   ist das Kronecker-Delta und der Index ,i bedeutet in diesem Abschnitt eine Ableitung nach yi. Der Nabla-Operator schreibt sich in krummlinigen Koordinaten

 

Hier wie im Folgenden muss die Einsteinsche Summenkonvention angewendet werden, der gemäß über in einem Produkt doppelt vorkommende Indizes, hier nur i, von eins bis zur Dimension des Raumes zu summieren ist.

Kontravariantes Vektorfeld

Bearbeiten

Die #Produktregel angewandt auf ein kontravariantes Vektorfeld[2]:146 führt zu

 

Der Gradient des kovarianten Basisvektors kann mit den Christoffelsymbolen[2]:340[7]:58   ausgedrückt werden:

 

Mit   lautet der Gradient schließlich

 

Darin ist   die sogenannte kovariante Ableitung der Komponente  .[2]:341[7]:61

Kovariantes Vektorfeld

Bearbeiten

Bei einem kovarianten Vektorfeld[2]:146   wird die Ableitung des kontravarianten Basisvektors benötigt, eine Ableitung, die auch mit Christoffelsymbolen ausgedrückt werden kann:

 

Der Gradient eines kontravarianten Basisvektors schreibt sich damit

 

Die #Produktregel liefert analog zum kontravarianten Vektor

 

mit der kovarianten Ableitung   der Komponente  .

Eigenschaften

Bearbeiten

Es ist die hier vereinbarte #Konvention   zu beachten.

Zusammenhang mit dem totalen Differenzial

Bearbeiten

Betrachtet wird eine infinitesimale Verschiebung in einem Vektorfeld:

 

Das vollständige oder totale Differenzial eines Vektorfeldes   ist:

    bzw. in Indexschreibweise    

Das totale Differenzial eines Skalarfeldes und eines Vektorfeldes haben somit (formal) dieselbe Form. Beim totalen Differenzial eines Skalarfeldes wird der Gradient mit dem Differenzial skalar multipliziert. Beim totalen Differenzial eines Vektorfeldes ist die Multiplikation zwischen dem Gradient (Matrixform) mit dem Differenzialvektor als Matrix-Vektor-Produkt durchzuführen.

Zusammenhang mit der Richtungsableitung

Bearbeiten

Mit dem Vektorgradient kann die Richtungsableitung in Richtung eines Vektors   berechnet werden:

 

Das hochgestellte bedeutet eine Transponierung. In der Strömungsmechanik wird die linke Darstellung mit dem Nabla-Operator gegenüber der rechten bevorzugt, die in der Kontinuumsmechanik üblich ist. Die mithilfe des Vektorgradienten berechnete Richtungsableitung entspricht der Richtungsableitung, die man durch Grenzwertbildung bekommt:

  für alle  

Interessiert diejenige Richtung  , in der die Richtungsableitung maximalen Betrag hat, ergibt sich das Eigenwertproblem

 

Der Tensor   ist symmetrisch und positiv semidefinit, sodass alle Eigenwerte reell und nicht negativ sind. Der zum größten Eigenwert gehörende Eigenvektor liefert die Richtung, in der die Richtungsableitung den größten Betrag hat.

Denn die Zielgröße ist

 

Ein Extremum unter der Nebenbedingung   berechnet sich mit einem Lagrange-Multiplikator λ:

 

Im Extremum müssen die Ableitungen nach allen Variablen verschwinden. Die Ableitung nach dem Lagrange-Multiplikator

 

bedeutet, dass wie gewünscht die Nebenbedingung notwendig eingehalten wird. Die Ableitung nach   in Richtung   liefert

 

weil C symmetrisch ist. Da dies für alle   gelten soll, ist das gleichbedeutend mit dem oben angegebenen Eigenwertproblem.

Zusammenhang mit Rotation und Divergenz

Bearbeiten

Der Vektorgradient beinhaltet alle partiellen Ableitungen der Komponenten eines Vektorfeldes, die bei der Rotation und Divergenz eines Vektorfeldes gebraucht werden. Es ist zu vermuten, dass diese Operatoren aus dem Gradient eines Vektorfeldes ableitbar sind. Tatsächlich ist[9]

 

für alle konstanten Vektoren  . Der Tensor in der eckigen Klammer ist schiefsymmetrisch und dessen dualer axialer Vektor (·)× ist die Rotation. Der duale axiale Vektor ist die negative Hälfte der Vektorinvariante  , bei der das dyadische Produkt ⊗ durch das Kreuzprodukt × ersetzt ist:

 

Die Spur des Vektorgradienten liefert die Divergenz:

 

Rechenregeln

Bearbeiten

Für alle Konstanten  ,  , total differenzierbaren Skalarfelder   und Vektorfelder   gilt mit der vereinbarten #Konvention  :

Linearität
 
 
 
Produktregel
 
 
 
In drei Dimensionen ist speziell[2]:367
 
Integralsätze[6]:45
 
Dabei ist der Integrationsweg von   nach   beliebig. Diese Wegunabhängigkeit zeichnet Gradientenfelder aus[2]:433.
 
Hier ist   ein zweimal stetig differenzierbares Feld und   der nach außen gerichtete Normaleneinheitsvektor auf der geschlossenen Oberfläche A des Volumens V.

Anwendungen

Bearbeiten

Es ist die hier vereinbarte #Konvention   zu beachten.

Verformungen

Bearbeiten
 
Abb. 2: Verlagerung von materiellen Linien

Der schon angesprochene Deformationsgradient ist die grundlegende Größe zur Beschreibung von Verformungen von Körpern. Lokal stellen sich bei einer Verformung Längenänderungen und Winkeländerungen zwischen materiellen Linienelementen ein, die man sich in das Material eingeritzt denken kann, siehe Bild. Die Längenänderungen korrespondieren mit Dehnungen und die Winkeländerungen mit Scherungen im Material.

In der Kontinuumsmechanik gibt die Bewegungsfunktion   den Ort   an, an dem zur Zeit t ein Partikel ist, das zu einer definierten Zeit t0 am Ort   war. Der Deformationsgradient F kann aus

 

berechnet werden, was seine Transformationseigenschaften der Linienelemente im undeformierten Zustand ( ) in den deformierten ( ) verdeutlicht.

Substanzielle Beschleunigung

Bearbeiten

In der Fluidmechanik wird die Eulersche Betrachtungsweise eingenommen, die das Vektorfeld der Geschwindigkeit   als Funktion des Ortes   und der Zeit t benutzt. Der Impulssatz eines Kontinuums besagt, dass eine volumenverteilte Kraft, wie die Schwerkraft eine ist, die Partikel des Körpers beschleunigt. Um das darzustellen, wird die Geschwindigkeit des Partikels mittels der Bewegungsfunktion   eingeführt, die den Ort angibt, an dem sich das Partikel   zur Zeit t befindet:

 

Der Überpunkt bildet hier die Substanzielle Zeitableitung. Für den Impulssatz kann nun die Substanzielle Beschleunigung als Zeitableitung der Geschwindigkeit bei festgehaltenem Partikel   berechnet werden:

 

Der zweite Summand stellt einen konvektiven Anteil dar, der physikalisch daraus resultiert, dass das Partikel auch dadurch beschleunigt werden kann, dass es von einem schneller oder langsamer fließenden Stromfaden mitgenommen wird. Der Geschwindigkeitsgradient[10]:42 ff.[6]:230 ff.   hat eine fundamentale Bedeutung in der Fluidmechanik.

Objektive Zeitableitung

Bearbeiten

Ein Insasse eines fahrenden Zuges wird die Geschwindigkeit eines vorbeifliegenden Vogels anders beurteilen als ein in der Nähe befindlicher Fußgänger. Die Geschwindigkeit ist demnach vom Standpunkt abhängig, sie ist genauer nicht bezugssysteminvariant oder kürzer nicht objektiv.

Für die Formulierung eines Materialmodells, in dem die Raten konstitutiver Variablen auftreten, wie beispielsweise beim newtonschen Fluid, werden jedoch objektive Zeitableitungen dieser Variablen benötigt. Denn es entspricht nicht der Erfahrung, dass ein bewegter Beobachter ein anderes Materialverhalten misst als ein ruhender.

 
Geschwindigkeitsfeld mit Stromlinien (blau). Wo die roten Linien nahe beieinander liegen ist die Geschwindigkeit   hoch, andernorts gering.

Ein Newtonsches Fluid besitzt Viskosität, bei der zur Aufrechterhaltung einer Scherströmung eine Kraft erforderlich ist. Scherströmungen weisen unterschiedliche Geschwindigkeiten von benachbarten Fluidelementen auf, d. h. es treten Geschwindigkeitsgradienten[10]:12[6]:83 auf wie im unteren Teil des Bildes. Weil die Geschwindigkeit selbst nicht objektiv ist, siehe Hauptartikel, ist es ihr Gradient ebenfalls nicht[10]:37 f. und letzterer ist deshalb für die Modellierung von viskosen Fluiden ungeeignet. Der symmetrische Anteil des Geschwindigkeitsgradienten ist aber objektiv und wird bei der Modellierung von newtonschen Fluiden benutzt, was auf die Navier-Stokes-Gleichungen führt, die auf diese Weise invariant gegenüber einer Galilei-Transformation sind.

Tensorgradient

Bearbeiten

Mit dem skalaren Operator   kann auch der Gradient eines Tensorfeldes T gebildet werden, wobei ein Tensorgradient[2]:356 entsteht:

 

In krummlinigen Koordinaten   und dem Nabla-Operator   (Notation siehe #Allgemein krummlinige Koordinaten) wird daraus:

 

Soll das Argument wie beim Vektorgradient rechts vom Operator stehen, siehe #Konvention, dann lautet der Tensorgradient

 

Für einen Tensor zweiter Stufe gibt es in krummlinigen Koordinaten vier Darstellungen:

 

Mit den Ableitungen der Basisvektoren

 

ergibt sich in der ersten Ausführung:

 

mit der kovarianten Ableitung der Tensorkomponente

 

Analog ergibt sich in den anderen Darstellungen:[2]:348, 356 f.

 

siehe auch die Anwendung der Christoffelsymbole bei Tensorfeldern.

Beispiele

Bearbeiten

Sei   der Ortsvektor und r sein Betrag. Dann ist mit dem Einheitstensor 1:

 
 
 

siehe Gradient (Mathematik)#Nützliche Formeln. Mit der #Produktregel berechnet sich damit

 

Die beiden letzten Formeln werden z. B. bei der kartesischen Multipolentwicklung verwendet.

Als weiteres Beispiel wird das Vektorfeld

 

angeführt, wo a eine beliebige Konstante ist. Der Gradient wird mit der vereinbarten #Konvention wie folgt aus der Richtungsableitung berechnet:

 

 
Abb. 3: Bei a=3/2 auf gleich großes Urbild skalierte Abbildungen des Kreises, die sich der Ellipse (blau) bei Radien 0,1, 0,01 und 0,001 annähern, sowie Eigenvektoren 1 (grün) und 2 (rot) mit ihren Bildern (türkis bzw. magenta).

Im Ursprung nimmt der Gradient die Form

 

an. Die maximale Richtungsableitung ergibt sich aus dem Eigensystem des Tensors

 

Er hat bei   die Eigenwerte und Eigenvektoren

 
 

Die größte Richtungsableitung ist in Richtung  , die durch den Gradient auf   abgebildet wird.

Siehe auch

Bearbeiten

Literatur

Bearbeiten
  1. Bedeutungsübersicht: Gradient. Duden online, abgerufen am 28. Oktober 2020.
  2. a b c d e f g h i j k Wolfgang Werner: Vektoren und Tensoren als universelle Sprache in Physik und Technik. Tensoralgebra und Tensoranalysis. Band 1. Springer Vieweg Verlag, Wiesbaden 2019, ISBN 978-3-658-25271-7, doi:10.1007/978-3-658-25272-4.
  3. Hugo Sirk: Einführung in die Vektorrechnung: Für Naturwissenschaftler, Chemiker und Ingenieure. Springer-Verlag, 2013, ISBN 3-642-72313-6, Kap. 5.4 "Das Vektorfeld und der Vektorgradient".
  4. a b C. B. Lang, N. Pucker: Mathematische Methoden in der Physik. Springer Spektrum, Berlin, Heidelberg 2016, ISBN 978-3-662-49312-0.
  5. a b M. E. Gurtin: The Linear Theory of Elasticity. In: S. Flügge (Hrsg.): Handbuch der Physik. Band VI2/a, Bandherausgeber C. Truesdell. Springer, 1972, ISBN 3-540-05535-5, S. 10.
  6. a b c d e Holm Altenbach: Kontinuumsmechanik. Einführung in die materialunabhängigen und materialabhängigen Gleichungen. Springer-Verlag, Berlin, Heidelberg 2012, ISBN 978-3-642-24118-5, S. 43, doi:10.1007/978-3-642-24119-2.
  7. a b c P. Haupt: Continuum Mechanics and Theory of Materials. Springer, 2002, ISBN 3-540-43111-X.
  8. J. Betten: Kontinuumsmechanik. Elastisches und inelastisches Verhalten isotroper und anisotroper Stoffe. 2. erw. Auflage. Springer, Berlin, Heidelberg u. a. 2001, ISBN 978-3-642-62645-6, doi:10.1007/978-3-642-56562-5.
  9. Johannes Wandinger: Gradient, Divergenz und Rotation. (Pdf) 13. November 2017, abgerufen am 2. November 2020.
  10. a b c Ralf Greve: Kontinuumsmechanik. Springer, 2003, ISBN 978-3-642-62463-6, S. 4, doi:10.1007/978-3-642-55485-8 (eingeschränkte Vorschau in der Google-Buchsuche).