Calcul tensoriel/Notions élémentaires

Un livre de Wikibooks.

Sections

[modifier] Introduction

Pour les besoins du calcul tensoriel, nous serons amenés a utiliser une notation indicielle, ces indices pouvant êtres aussi bien en indice qu'en exposant. Il faudra particulièrement faire attention au fait que des indices en exposants ne représentent pas une puissance (en cas de réelle utilisation de puissance on peut toujours écrire (xi)n).

Il faudra faire aussi faire attention au fait que l'on va utiliser des conventions de simplification pour alléger l'écriture. En particulier, la convention d'Einstein: si on a un vecteur:

\vec{V} = \sum_{i=1}^n V^i \vec{e_i}

On omettra le signe "somme" et on écrira juste:

\vec{V} = V^i \vec{e_i}

Remarquez la convention d'écrire les composantes avec un indice en exposant, et les vecteurs de bases avec un indice inférieurs (lorsque l'on utilise les composantes contravariantes - nous verrons ce que c'est précisément plus tard).

Exemples:

  • le vecteur: \vec{V} = a^1 \vec{e_1} + a^2 \vec{e_2} + a^3 \vec{e_3} s'écrira \vec{V} = a^i \vec{e_i} (en spécifiant évidemment que n = 3)
  • le système:

a11x1 + a12x2 = b1

a21x1 + a22x2 = b2

s'écrira aijxj = bi (n = 2: attention que l'indice de sommation est j dans ce cas et non pas i.)

Pour les dérivées partielles par rapport aux coordonnées, on utilise indifféremment les notations

\frac{\partial f}{\partial x^i} = \partial_i f = f_{,i}

[modifier] Transformations entre systèmes de coordonnées

On s'intéresse maintenant aux changements de coordonnées entre deux systèmes de coordonnées différents.

Soient deux systèmes de coordonnées a et b (on va donc noter la ième composante d'un point respectivement \left[x_a\right]^i et \left[x_b\right]^i dans l'un ou l'autre des système). Considérons les fonctions de transformation :

\left[x_b\right]^i = \left[\Theta_{ba}\right]^i \left(\left[x_a\right]^0, \left[x_a\right]^1, \ldots\right)
\left[x_a\right]^i = \left[\Theta_{ab}\right]^i \left(\left[x_b\right]^0, \left[x_b\right]^1, \ldots\right)

Le jacobien de la transformation relie les variations infinitésimales des coordonnées :

d\left[x_a\right]^i = \sum_j \left[\theta_{ab}\right]^i{}_j d\left[x_b\right]^j

Ainsi on a:

\left[\theta_{ab}\right]^i{}_j = \frac{\partial\left[\Theta_{ab}\right]^i}{\partial\left[x_b\right]^j}

Le jacobien de la transformation inverse est l'inverse du jacobien. Introduisant \delta^i_j le symbole de Kronecker, on a

\sum_j \left[\theta_{ab}\right]^i{}_j \left[\theta_{ba}\right]^j{}_k = \delta^i_k

Et le symbole de Kronecker est défini comme suit:

\delta_{ij} = \delta_i^j = \delta^{ij} = \begin{cases} 1 & \mbox{si } i=j \\ 0 & \mbox{si } i \ne j \end{cases}

[modifier] Base naturelle

Nous allons maintenant introduire la notion de "base naturelle".

Par définition, la base naturelle au point M est la base formée par les vecteurs:

\mathbf{e}_i = \frac{\partial \mathbf{OM}}{\partial x^i}

ou O représente évidemment l'origine.

\mathbf{e}_i est donc la base naturelle de vecteurs locaux associés à un système de coordonnées quelconque xi. Un élément infinitésimal s'écrit dans cette base:

d\mathbf{l} = \sum_i dx^i \mathbf{e}_i

On a par définition:

\mathbf{e}_i = \frac{\partial \mathbf{l}}{\partial x^i}

  • Remarques
    1. La lettre \mathbf{l} étant muette, on voit parfois écrit \mathbf{e}_i = \frac{\partial}{\partial x^i}, avec le risque de confusion avec l'opérateur dérivation.
    2. Sauf dans le cas d'un repère cartésien, la base naturelle varie d'un point à un autre. Chaque symbole \mathbf{e}_1, \mathbf{e}_2, etc. représente en fait un champ de vecteurs.
    3. Les vecteurs de base se transforment selon la formule \left[\mathbf{e}_a\right]_i = \sum \left[\mathbf{e}_b\right]_i \left[\theta_{ba}\right]^i{}_j. Démonstration

[modifier] composantes covariantes et contravariantes

Soit V un espace vectoriel, et soit \vec{V}\in V et {ei} une base de V.

On note, avec la convention d'Einstein:

\vec{V}=V^i e_i

Tout simplement, les nombre Vi sont les composantes contravariantes du vecteur \vec{V}. Ce sont les composantes que l'on utilise habituellement.

Les composantes covariantes d'un vecteur sont les composantes d'un vecteur sur la base duale. On note:

\vec{V}=V_i e^i

ou la base {ei} est la base duale de {ei}, définie par:

e^i e_j = \delta^i_j

Remarques:

  • Tous les vecteurs de la base duale sont orthogonaux à tous les vecteurs de la base de départ d'indices différents (produit scalaire nul).
  • Le produit scalaire entre un vecteur de la base ordinaire et un vecteur de la base duale mais de même indice cette fois, vaut 1.
  • On peut déduire qu'une base orthogonale est identique à sa base duale.
  • La base duale de la base duale est la base de départ.

[modifier] Tenseurs

Dans un espace de dimension N, un tableau à P indices T_{i_1 i_2 \ldots i_P}, chaque indice pouvant prendre N valeurs représente les composantes covariantes d'un tenseur d'ordre P dans un système de coordonnées données, si ce tableau se transforme comme les vecteurs de base locale lors d'un changement de système de coordonnées.

\left[T_a\right]_{i_1 i_2 \ldots i_P} = \left[T_b\right]_{j_1 j_2 \ldots j_P} \left[\theta_{ba}\right]^{j_1}{}_{i_1} \left[\theta_{ba}\right]^{j_2}{}_{i_2} \ldots \left[\theta_{ba}\right]^{j_P}{}_{i_P}

Un tableau T^{i_1 i_2 \ldots i_P} représente les composantes contravariantes d'un tenseur d'ordre P dans un système de coordonnées données, si ce tableau se transforme comme les variations infinitésimales des coordonnées lors d'un changement de système de coordonnées.

\left[T_a\right]^{i_1 i_2 \ldots i_P} = \left[T_b\right]^{j_1 j_2 \ldots j_P} \left[\theta_{ab}\right]^{i_1}{}_{j_1} \left[\theta_{ab}\right]^{i_2}{}_{j_2} \ldots \left[\theta_{ab}\right]^{i_P}{}_{j_P}

On verra plus bas que des composantes contravariantes et covariantes peuvent correspondre au même tenseur, la transformation étant obtenue au moyen du tenseur métrique.

De la même manière qu'on définit un champ vectoriel comme un vecteur fonction de la position dans l'espace, on définit un champ tensoriel comme un tenseur fonction de la position dans l'espace. Dans la suite, on utilisera le simple mot tenseur, alors qu'on considérera en général un champ tensoriel.

Il n'y a pas de difficulté à définir un tenseur en composantes mixtes. On aura par exemple

\left[T_a\right]^{i}{}_j{}^k = \left[T_b\right]^{l}{}_m{}^n \left[\theta_{ab}\right]^{i}{}_{l} \left[\theta_{ba}\right]^{m}{}_{j} \left[\theta_{ab}\right]^{k}{}_{n}

[modifier] Produit tensoriel

Deux tenseurs A et B d'ordre P et Q étant donnés par leurs NP et NQ composantes covariantes, contravariantes ou mixtes, le produit des composantes définit un tableau de NP + Q composantes. Ce tableau se transforme évidemment comme les composantes d'un tenseur C d'ordre P + Q, appelé produit tensoriel de A et B.

On écrira par exemple Cijk = AijBk.

[modifier] Tenseur métrique

Le tenseur métrique, noté gij, est le tenseur produit scalaire des vecteurs de la base naturelle. Il est symétrique :

g_{ij} = g_{ji} = \mathbf{e}_i \cdot \mathbf{e}_j

Le carré d'un élément de longueur est la forme quadratique

\left(dl\right)^2 = \sum_{ij} g_{ij} dx^i dx^j

À RÉDIGER

[modifier] Transformation contraco

On transforme les composantes contravariantes d'un tenseur en composantes covariantes au moyen du tenseur métrique :

xi = ∑ gijxj
j

Étant donné que gij est la matrice inverse de gij, on a

xi = ∑ gijxj
j
  • Remarques
    1. Avec la convention d'Einstein, on sous-entend le symbole \sum dès qu'un indice figure à la fois en haut et en bas. On écrit ainsi xi = gijxj, et xi = gijxj.
    2. Il y a une convention plus radicale utilisée par certains auteurs suivant laquelle on omet également le tenseur métrique dès lors qu'il intervient dans une contraction. Ainsi, à la place de aibi = gijaibj, on écrit aibi.
    3. La transformation contraco peut être appliquée à plusieurs indices d'un tenseur d'ordre quelconque. Par exemple aij = gikgjlakl.

[modifier] Contraction

Un tenseur d'ordre P étant donné, on obtient un tenseur d'ordre P − 2 en sommant toutes les termes correspondant aux mêmes valeurs de deux indices donnés, l'un covariant, l'autre contravariant.

Aik = ∑ Bijkj
j

Avec la convention d'Einstein, on écrit

Aik = Bijkj

Puisque l'on a B^{i}{}_j{}^{kn} = B^{i}{}_{}^{mkn} g_{jm}, on peut effectuer la contraction sur deux composantes contravariantes ou bien deux composantes covariantes au moyen du tenseur métrique :

Aik = Bimkngmn = Bimkngmn

[modifier] Produit scalaire

Le produit scalaire de deux vecteurs ai et bj est la contraction de leur produit tensoriel. Il s'exprime donc au moyen du tenseur métrique gij :

\mathbf{a} \cdot \mathbf{b} = a^i b^j g_{ij}

C'est un tenseur d'ordre 0, indépendant du choix du système de coordonnées.

[modifier] Tenseur dualiseur

[modifier] Introduction

Étant donné un espace de dimension N, le symbole de Levi-Civita d'ordre N \epsilon^{i_1 i_2 \ldots i_N}, aussi appelé pseudo-tenseur unité complètement antisymétrique n'est pas un tenseur. Par exemple, ses composantes devraient être multipliées par 2N lorsque le système de coordonnées est réduit d'un facteur 2.

[modifier] Définition du tenseur dualiseur

[modifier] Tenseur dualiseur en coordonnées contravariantes

La formule

*^{i_1 i_2 \ldots i_N} = \frac{1}{\sqrt{\det g}}\epsilon^{i_1 i_2 \ldots i_N}

où detg est le déterminant du tenseur métrique, définit bien un tenseur à partir du symbole de Levi-Civita d'ordre N. Ce tenseur, aussi appelé tenseur de Levi-Civita et noté \tilde{\epsilon} est ici appelé tenseur dualiseur. Il est intéressant de noter que pour qu'un scalaire soit un tenseur (d'ordre 0), il faut qu'il soit indépendant du choix du système de coordonnées. Le déterminant du tenseur métrique n'est donc pas un tenseur mais son produit avec le symbole ε est bien un tenseur, indépendant du système de coordonnées.

[modifier] Tenseur dualiseur en coordonnées covariantes

On passe aux coordonnées covariantes en mettant en jeu N fois le tenseur métrique. Le produit de ces N termes avec le symbole de Levi-Civita d'ordre N est égal au déterminant du tenseur métrique et l'on obtient

*_{i_1 i_2 \ldots i_N} = g_{i_1 j_1} g_{i_2 j_2} \ldots g_{i_N j_N} \frac{1}{\sqrt{\det g}}\epsilon^{i_1 i_2 \ldots i_N} = \sqrt{\det g} \epsilon_{i_1 i_2 \ldots i_N}
[modifier] Détermination du signe

Le tenseur dualiseur garde néanmoins un caractère particulier. C'est le choix de la détermination de la racine carrée qui définit l'orientation du système de coordonnées. Retourner une coordonnée doit s'accompagner du changement de signe des composantes du tenseur.

Dans le cas de l'espace-temps quadridimensionnel, le déterminant du tenseur métrique est négatif. Deux possibilités s'offrent alors, considérer que le tenseur dualiseur est imaginaire pur, ou remplacer g par -g. Cette seconde possibilité facilite les calculs mais n'est pas satisfaisante d'un point de vue théorique [connecter avec la mesure de la distance spatiale ou temporelle comme racine carrée de l'intervalle d'espace-temps, dont le signe donne la nature, temps ou espace]

[modifier] Propriétés du tenseur dualiseur

[modifier] Produit de tenseurs dualiseurs

La formule suivante découle directement de la formule obtenue avec le symbole de Levi-Civita d'ordre N. Le symbole \delta^i_j est le symbole de Kronecker, représentant la matrice unité.

*^{i_1 i_2 \ldots i_P k_1 \ldots k_Q} *_{j_1 j_2 \ldots j_P k_1 \ldots k_Q} = Q! \times \left| \begin{matrix} \delta^{i_1}_{j_1} & \delta^{i_1}_{j_2} & \ldots & \delta^{i_1}_{j_P} \\ \delta^{i_2}_{j_1} & \delta^{i_2}_{j_2} & \ldots & \delta^{i_2}_{j_P} \\ \ldots & \ldots & \ldots & \ldots \\ \delta^{i_P}_{j_1} & \delta^{i_P}_{j_2} & \ldots & \delta^{i_P}_{j_P} \end{matrix}\right|
[modifier] résultat d'ordre 2
*^{k \; i_2\ldots i_N}*_{l \; i_2 \ldots i_N} = \left(N-1\right) ! \times \delta^k_l
[modifier] résultat d'ordre 4
*^{k m \; i_3\ldots i_N}*_{l n \; i_3 \ldots i_N} = \left(N-2\right) ! \times \left(\delta^k_l \delta^m_n - \delta^k_n \delta^m_l\right)

[modifier] Définition du tenseur dual

Le produit du tenseur dualiseur avec un tenseur d'ordre M dans un espace de dimension N définit un tenseur d'ordre N-M, son dual. \left[* a\right]^{i_1 i_2 \ldots i_{N-M}} = \frac{1}{M !} *^{i_1 i_2 \ldots i_N} a_{i_{N-M+1} \ldots i_N}

Le choix de faire la contraction sur les M derniers indices du tenseur dualiseur est tout à fait arbitraire. La contraction sur une autre famille d'indices fournit la même expression de * a, avec un signe éventuellement opposé. Il est aussi à noter que si l'on change l'orientation de l'espace en échangeant deux coordonnées, le tenseur dual change de signe.

[modifier] Tenseurs duaux complètement antisymétriques

Si dans un espace à N dimensions, A est un tenseur à M indices, avec 0 < = M < = N, alors le dual * A est un tenseur à N − M indices, lui-aussi complètement antisymétrique, et l'on a

* * A = A

On va le démontrer en dimension 3. La généralisation est facile.

[modifier] dimension 3

Le dual d'un vecteur bk est un tenseur antisymétrique à deux indices \left[* b\right]_{ij} = *_{ijk} b^k. Réciproquement, le dual d'un tenseur antisymétrique aij est un vecteur \left[*a\right]^i = \frac{1}{2} \times *^{ijk} a_{jk}.. Tous deux ont 3 composantes indépendantes.

Le dual d'un scalaire f est un tenseur complètement antisymétrique à trois indices \left[*f\right]^{ijk} = *^{ijk} f. Réciproquement, le dual d'un tenseur complètement antisymétrique à trois indices aijk est un scalaire \left[*a\right] = \frac{1}{6} *_{ijk} a^{ijk}. Le tenseur complètement antisymétrique à trois indices dans un espace de dimension 3 n'a qu'une composante indépendante. On l'appelle aussi pseudo-scalaire.

Le dual du dual d'un tenseur complètement antisymétrique à 0, 1, 2 ou 3 indices est le tenseur lui-même. Autrement dit, le dual du dual est l'opérateur identité pour les scalaires, les vecteurs ai, les tenseurs antisymétriques à deux indices aij = − aji et les tenseurs complètements antisymétriques à trois indices aijk = − aikj = − akji = − ajik.

En effet, pour un scalaire on a\left[**f\right] = \frac{1}{6} *_{ijk} \left[*f\right]^{ijk} = \frac{1}{6} *^{ijk} *_{jkl} f = f ; pour un vecteur on a \left[**b\right]^i = \frac{1}{2} *^{ijk} \left[*b\right]_{jk} = \frac{1}{2} *^{ijk} *_{jkl} b^l = \delta^i_l b^l = b^i ; pour un tenseur antisymétrique à deux indices on a \left[**a\right]^{ij} = \frac{1}{2} *^{ijk} \left[*a\right]_{k} = \frac{1}{2} *^{ijk} *_{klm} a^{lm} = \frac{1}{2}\left(\delta^{ij}_{lm} - \delta^{ij}_{ml}\right) a^{lm} = a^{ij} ; pour un tenseur compètement antisymétrique à trois indices on a \left[**a\right]^{ijk} = *^{ijk} \left[*a\right] = *^{ijk} \frac{1}{6} *_{lmn} a^{lmn} = \frac{1}{6}\left( \delta^{ijk}_{lmn} + \delta^{ijk}_{mnl} + \delta^{ijk}_{nlm} - \delta^{ijk}_{nml} - \delta^{ijk}_{mln} - \delta^{ijk}_{lnm}\right) a^{lmn} = a^{ijk}

[modifier] dimension 4

Le dual d'un scalaire a est le tenseur complètement antisymétrique pseudo-scalaire \left[*a\right]^{ijkl} = *^{ijkl} a

Le dual d'un vecteur al est le tenseur complètement antisymétrique \left[*a\right]^{ijk} = *^{ijkl} a_l. Tous deux ont 4 composantes indépendantes.

Le dual d'un tenseur antisymétrique akl est le tenseur antisymétrique \left[*a\right]^{ij} = *^{ijkl} a_{kl}. Tous deux ont 6 composantes indépendantes.

Le dual d'un tenseur complètement antisymétrique ajkl est le vecteur \left[*a\right]^{i} = *^{ijkl} a_{jkl}

Le dual d'un tenseur complètement antisymétrique pseudo-scalaire aijkl est le scalaire \left[*a\right] = *^{ijkl} a_{ijkl}

De la même manière qu'en dimension 3, il est facile de démontrer que l'opérateur dual du dual est l'opérateur unité pour les scalaires, les vecteurs, les tenseurs antisymétriques à deux indices et les tenseurs complètement antisymétriques à trois ou quatre indices.

[modifier] Produit vectoriel

Soient aj et bk deux vecteurs dans un espace de dimension N. Le produit vectoriel de ces vecteurs est le tenseur défini par

\left(\mathrm{a} \times \mathrm{b}\right)_{i_1 i_2 \ldots i_{N-2}} = *_{i_1 i_2 \ldots i_{N-2} j k} a^j b^k

  • en dimension 2, le produit vectoriel est le scalaire défini par \mathrm{a} \times \mathrm{b} = *_{jk} a^j b^k
  • En dimension 3, le produit vectoriel est le vecteur défini par \left(\mathrm{a} \times \mathrm{b}\right)_i = *_{ijk} a^j b^k, ou, en composantes contravariantes \left(\mathrm{a} \times \mathrm{b}\right)^i = g^{il} *_{ljk} a^j b^k.

[modifier] Dérivée covariante

[modifier] Dérivée partielle d'un vecteur

Soit un champ vectoriel \mathbf{v} de coordonnées \left[v_a\right]_i dans la base naturelle associée au système de coordonnées \left[x_a\right]^\square et de coordonnées \left[v_b\right]_k dans la base naturelle associée au système de coordonnées \left[x_b\right]^\square.

La dérivée partielle ou dérivée virgule

\left[v_a\right]_{i,j} = \left[\partial_a\right]_j \left[v_a\right]_i

n'est pas un tenseur. En effet, lorsqu'on dérive la formule de changement de coordonnées

\left[v_a\right]_i ( \left[x_a\right]^\Box ) = \sum_k{ \left[\theta_{ba}\right]^k{}_i \left[v_b\right]_k \left( \left[\Theta_{ba}\right]^\Box ( \left[x_b\right]^\Box ) \right) }

on obtient

\left[v_a\right]_{i_,j} = \sum_{kl}{\left[\theta_{ba}\right]^k{}_i (\left[v_b\right]_{k_,l}) \left[\theta_{ba}\right]^l{}_j} + \sum_{k}{ (\left[\theta_{ba}\right]^k{}_{i_,j}) \left[v_b\right]_k}

formule de transformation d'un tenseur deux fois covariant, troublée par la présence d'un second terme, contenant le jacobien de la matrice de changement de base.

[modifier] Symbole de Christoffel

Le symbole de Christoffel est défini à partir de la dérivée partielle des vecteurs de la base naturelle :

\Gamma^{k}_{ij} = \left(\partial_i \mathbf{e}_j\right) \mathbf{e}^k = \mathbf{e}_{j,i} \mathbf{e}^k

Étant donné la définition de la base naturelle, on peut écrire \Gamma^{k}_{ij} = \left(\partial_i \partial_j \mathrm{l}\right) \mathbf{e}^k pour mettre en évidence la symétrie du symbole de Christoffel par échange des indices bas :

\Gamma^{k}_{ij} = \Gamma^{k}_{ji}

  • Remarques
    1. Le symbole de Christoffel est aussi appelé connexion, avec un signe parfois différent.
    2. Ce symbole n'est pas un tenseur à cause du second terme de la formule de transformation. On définit néanmoins le symbole \Gamma_{l|ij} = g_{lk} \Gamma^{k}_{ij}
    3. Ce symbole permet de calculer le tenseur dérivée covariante d'un tenseur.

[modifier] Symbole de Christoffel fonction du tenseur métrique

Le symbole de Christoffel s'écrit en fonction de la dérivée partielle g_{ij,k} = \frac{\partial g_{ij}}{\partial x^k} du tenseur métrique:

\Gamma^k_{ij} = \frac{1}{2} g^{kl} \left(g_{lj,i} + g_{li,j} - g_{ij,l}\right)

Démonstration

[modifier] Contraction du symbole de Christoffel

La contraction du symbole de Christoffel s'exprime à partir de la dérivée partielle du déterminant du tenseur métrique.

\Gamma^i_{ij} = \frac{1}{2 \det{g}} \partial_j \left(\det g\right) = \frac{1}{\sqrt{\det{g}}} \partial_j \sqrt{\det{g}} = \partial_j \left(\log \sqrt{\det{g}}\right)

Calcul tensoriel/Notions élémentaires/Symbole de Christoffel/Contraction/Démonstration

  • Remarques
    1. Le symbole de Christoffel étant symétrique, le résultat ne dépend de l'indice de contraction : \Gamma^i_{ij} = \Gamma^i_{ji}
    2. Ni le symbole de Christoffel ni la dérivée partielle ne représentent des tenseurs. Néanmoins cette formule peut figurer dans des expressions qui représentent des tenseurs, par exemple dans la formule de la

divergence.

[modifier] Transformation du symbole de Christoffel lors d'un changement de coordonnées

\left[\Gamma_a\right]^i_{kl} = \left[\Gamma_b\right]^m_{np} \frac{\partial\left[x_a\right]^i}{\partial\left[x_b\right]^m} \frac{\partial\left[x_b\right]^n}{\partial\left[x_a\right]^p} \frac{\partial\left[x_b\right]^p}{\partial\left[x_a\right]^l} + \frac{\partial^2\left[x_b\right]^m}{\partial\left[x_a\right]^k\partial\left[x_a\right]^l} \frac{\partial\left[x_a\right]^i}{\partial\left[x_b\right]^m}

[modifier] Pseudo-contraction de la dérivée partielle du tenseur métrique

Le produit contracté du tenseur métrique et de sa dérivée partielle change de signe lorsqu'on remonte les indices d'un terme du produit et que l'on descend les indices de l'autre terme : gijgij,k = − gijgij,k. Démonstration.

[modifier] Dérivée covariante

La dérivée covariante, définie par

f;i = f,i
a^j{}_{;i} = a^j{}_{,i} + \Gamma^j_{mi} a^m
a_{k;i} = a_{k,i} - \Gamma^n_{ki} a_n
a^j{}_{k;i} = a^j{}_{k,i} + \Gamma^j_{mi} a^m{}_k - \Gamma^n_{ki} a^j{}_n

etc. est un tenseur. La présence des symboles de Christoffel permet d'anihiler le jacobien dans la formule de transformation de la dérivée virgule.

[modifier] Nullité de la dérivée covariante du tenseur métrique

Calcul tensoriel/Notions élémentaires/Tenseur métrique/Nullité de la dérivée covariante

[modifier] Nullité de la dérivée covariante du tenseur dualiseur

La dérivée covariante du tenseur dualiseur est nulle :

*^{i_1 i_2 \ldots i_N}{}_{;j} = 0^{i_1 i_2 \ldots i_N}{}_{j}.

Démonstration.

[modifier] Système de coordonnées localement géodésique

Il est possible de construire un système de coordonnées qui annule les dérivées partielles du tenseur métrique, et donc le symbole de Christoffel en un point donné, sans modifier le tenseur métrique en ce point. Démonstration.

Il est possible de construire un système de coordonnées qui annule les dérivées partielles du tenseur métrique, et donc le symbole de Christoffel sur une ligne donnée.

[modifier] Opérateurs différentiels

[modifier] Gradient

Si f est un champ scalaire indépendant du système de coordonnées, c'est un tenseur d'ordre 0, et sa dérivée partielle est égale à sa dérivée covariante : \left(\nabla f\right)_i = \partial_i f = f_{,i} = f_{;i}, aussi appelée gradient de f. Ce vecteur est habituellement exprimé en composantes contravariantes

\left(\nabla f\right)^i = g^{ij} f_{;j}

Cette formule permet, une fois établi le tenseur métrique, de calculer facilement le gradient dans un système de coordonnées quelconque. Voir par exemple le gradient en coordonnées cylindriques et le gradient en coordonnées sphériques.

  • Remarques
    1. Le gradient généralisé d'un tenseur quelconque peut être défini simplement comme sa dérivée covariante. Cette opération ajoute un indice au tenseur.

[modifier] Divergence

La divergence d'un tenseur est le tenseur obtenu en contractant un des indices de la dérivée covariante avec l'indice de la dérivation.

[modifier] Divergence d'un vecteur

Pour un champ vectoriel \mathbf{v}, on a

\nabla \mathbf{v} = v^i{}_{;i} = v^i{}_{,i} + \Gamma^i_{ij} v^j

Mettant a profit la formule de contraction \Gamma^i_{ij} = \frac{1}{\sqrt{\det g}}\partial_j \sqrt{\det g}, on a

\nabla \mathbf{v} = \frac{1}{\sqrt{\det g}} \partial_j \left(\sqrt{\det g} \; v^j \right)

Cette formule permet, une fois établi le tenseur métrique, de calculer facilement la divergence dans un système de coordonnées quelconque. Voir par exemple la divergence en coordonnées cylindriques et la divergence en coordonnées sphériques.

[modifier] Divergence d'un tenseur d'ordre 2

Suivant le même chemin que pour la divergence d'un champ de vecteurs, on écrit

a^{ij}{}_{;j} = a^{ij}{}_{,j} + \Gamma^i_{lm} a^{lm} + \Gamma^l_{lm} a^{im} = \frac{1}{\sqrt{\det g}} \partial_k \left(\sqrt{\det g} \; a^{ik}\right) + \Gamma^i_{lm} a^{lm}

[modifier] Divergence d'un tenseur antisymétrique d'ordre 2

Dans le cas d'un tenseur antisymétrique, on a

a^{ij}{}_{;j} = -a^{ji}{}_{;j} = \frac{1}{\sqrt{\det g}} \partial_k \left(\sqrt{\det g} \; a^{ik}\right)

En effet, le terme \Gamma^i_{lm} a^{lm} est nul puisque \Gamma^i_{lm} a^{lm} = -\Gamma^i_{lm} a^{ml} = -\Gamma^i_{ml} a^{ml} = -\Gamma^i_{lm} a^{lm}.

[modifier] Remarques

  • En écriture quadridimensionnelle, les équations de Maxwell mettent en jeu la divergence du tenseur de champ électromagnétique et de son dual, qui sont des tenseurs antisymétriques d'ordre 2.

[modifier] Laplacien

Le laplacien est la divergence du gradient, la divergence étant prise sur l'indice tensoriel créé par le gradient.

Cette définition est valable pour un scalaire ou un tenseur quelconque a. La laplacien \nabla^2 a = a^{;i}_{;i} = g^{ij} a_{;i;j} a le même nombre d'indices que a.

  • Pour un champ scalaire

\nabla^2 a = \frac{1}{\sqrt{\det g}} \partial_i\left( \sqrt{\det g} \; g^{ij} \partial_j a \right)

  • Pour un champ vectoriel

...À RÉDIGER...

Ces formules permettent, une fois établi le tenseur métrique, de calculer facilement le laplacien dans un système de coordonnées quelconque. Voir par exemple le laplacien en coordonnées cylindriques et le laplacien en coordonnées sphériques.

[modifier] Rotationnel

[modifier] Tenseur rotationnel

Étant donné un champ de vecteurs covariants ai dans un espace de dimension quelconque, la dérivée covariante D_j a_i = a_{i;j} = a_{i,j} - \Gamma^k_{ij} a_k est un tenseur. Le tenseur rotationnel, défini comme \left[\mathrm{rot} \; a\right]_{ij} = a_{i;j} - a_{j;i} est par construction un tenseur antisymétrique.

[modifier] Expression à partir de la dérivée simple

La symétrie \Gamma^k_{ij} = \Gamma^k_{ji} du symbole de Christoffel permet d'écrire le tenseur rotationnel à partir de la dérivée simple : \left[\mathrm{rot} \; a\right]_{ij} = a_{i,j} - a_{j,i}.

[modifier] Rotationnel en dimension 3

En dimension 3, le tenseur dualiseur permet de construire le vecteur dual d'un tenseur antisymétrique d'ordre 2. Le rotationnel d'un champ de vecteurs \mathbf{a}tridimensionnel est défini comme le dual du tenseur rotationnel : \left[\mathrm{rot} \; \mathbf{a}\right]^{\alpha} = \frac{1}{2} *^{\alpha\lambda\mu} \left[\mathrm{rot} \; a\right]_{\lambda\mu}.

Partant d'un champ de vecteurs en coordonnées contravariantes aν, mettant à profit l'antisymétrie du tenseur dualiseur, la nullité de la dérivée covariante du tenseur métrique γμν ainsi que sa symétrie, on trouve \left[\mathrm{rot} \; \mathbf{a}\right]^{\alpha} = *^{\alpha\lambda\mu} \gamma_{\mu\nu} a^{\nu}{}_{;\lambda}.

[modifier] Tenseur de courbure

La dérivée covariante seconde d'un champ scalaire f est un tenseur d'ordre 2

f_{;i;j} = f_{,i,j} - \Gamma^k_{ij} f_{;k}

Elle est symétrique parce que le symbole de Christoffel \Gamma^k_{ij} est invariant par échange des indices bas (espace sans torsion).

En revanche pour un champ de vecteurs aj, les dérivations covariantes ne commutent pas. Le tenseur de courbure de Riemann Rijkl permet de calculer leur différence

ai;k;l − ai;l;k = Rijklaj

[modifier] Expression I

R_i{}^j{}_{kl} = - \Gamma^j_{ik,l} + \Gamma^j_{il,k} + \Gamma^m_{il} \Gamma^j_{mk} - \Gamma^m_{ik} \Gamma^j_{ml}

Démonstration

[modifier] Expression II

R_{ijkl} = g^{pq} \left(\Gamma_{p|ik}\Gamma_{q|jl} - \Gamma_{p|il}\Gamma_{q|jk}\right) + \frac{1}{2} \left(g_{ik,j,l} + g_{jl,i,k} - g_{il,j,k} - g_{jk,i,l}\right)

avec

\Gamma_{p|ik} = g_{pq} \Gamma^q_{ik} = \frac{1}{2} \left(g_{ip,k} + g_{kp,i} - g_{ik,p}\right)


[modifier] Propriétés

Calcul tensoriel/Notions élémentaires/Tenseur de courbure/Propriétés

[modifier] Symétries

Calcul tensoriel/Notions élémentaires/Tenseur de courbure/Symétries

[modifier] Relation cyclique

Calcul tensoriel/Notions élémentaires/Tenseur de courbure/Relation cyclique

[modifier] Identité de Bianchi

Calcul tensoriel/Notions élémentaires/Tenseur de courbure/Identité de Bianchi

En contractant le tenseur de courbure entre un indice de la première paire et un indice de la seconde paire, on obtient le tenseur de Ricci, clé des équations d'Einstein.


[modifier] Tenseur de Ricci

Le tenseur de Ricci est obtenu en contractant le tenseur de courbure entre un indice de la première paire et un indice de la seconde paire :

Rik = Rimkm = gjlRijkl

Grâce à la symétrie par paires du tenseur de courbure, le tenseur de Ricci est symétrique.

Le tenseur de Ricci complètement contracté est un scalaire :

R = gijRij

La divergence du tenseur d'Einstein R^{ij} - \frac{1}{2} g^{ij} R est nulle :

\left[R^{ij} - \frac{1}{2} g^{ij} R\right]_{;j} = 0

Cette équation fondamentale se démontre en mettant en jeu la nullité de la dérivée covariante du tenseur métrique.

C'est en identifiant le tenseur d'Einstein et le tenseur d'énergie-impulsion que l'on obtient l'équation d'Einstein qui fonde la relativité générale.

Récupérée de « http://fr.wikibooks.org/wiki/Calcul_tensoriel/Notions_%C3%A9l%C3%A9mentaires »