Analyse vectorielle

Un article de Wikipédia, l'encyclopédie libre.
Aller à : navigation, rechercher

L'analyse vectorielle est une branche des mathématiques qui étudie les champs de scalaires et de vecteurs suffisamment réguliers des espaces euclidiens, c'est-à-dire les applications différentiables d'un ouvert d'un espace euclidien E à valeurs respectivement dans \R et dans E. Du point de vue du mathématicien, l'analyse vectorielle est donc une branche de la géométrie différentielle. Cette dernière inclut l'analyse tensorielle qui apporte des outils plus puissants et une analyse plus concise entre autres des champs de vecteurs.

Mais l'importance de l'analyse vectorielle provient de son utilisation intensive en physique et dans les sciences de l'ingénieur. C'est de ce point de vue que nous la présenterons, et c'est pourquoi nous nous limiterons le plus souvent au cas où E = \R^3 est l'espace usuel à trois dimensions. Dans ce cadre, un champ de vecteurs associe à chaque point de l'espace un vecteur (à trois composantes réelles), tandis qu'un champ de scalaires y associe un réel. Imaginons par exemple l'eau d'un lac. La donnée de sa température en chaque point forme un champ de scalaires, celle de sa vitesse en chaque point, un champ de vecteurs (pour une approche plus théorique, voir géométrie différentielle).

Le calcul vectoriel et l'analyse vectorielle furent développés à la fin du XIXe siècle par J. Willard Gibbs et Oliver Heaviside à partir de la théorie des quaternions (due à Hamilton) ; la plupart des notations et de la terminologie furent établies par Gibbs et Edwin Bidwell Wilson dans leur livre de 1901, Vector Analysis (Analyse vectorielle).

Principaux opérateurs différentiels linéaires[modifier | modifier le code]

Le gradient, la divergence et le rotationnel sont les trois principaux opérateurs différentiels linéaires du premier ordre. Cela signifie qu'ils ne font intervenir que des dérivées partielles (ou différentielles) premières des champs, à la différence, par exemple, du laplacien qui fait intervenir des dérivées partielles du second ordre.

On les rencontre en particulier

Opérateur formel nabla[modifier | modifier le code]

Article détaillé : nabla.

L'opérateur nabla \nabla tire son nom d'une lyre antique qui avait la même forme de triangle pointant vers le bas. Il s'agit d'un opérateur formel de \R^3 défini en coordonnées cartésiennes par


\nabla =
\begin{pmatrix}
 \frac {\partial}{\partial x} \\[7pt]
 \frac {\partial}{\partial y} \\[7pt]
 \frac {\partial}{\partial z}
\end{pmatrix}
.

On écrit aussi \vec\nabla pour souligner que formellement, l'opérateur nabla a les caractéristiques d'un vecteur. Il ne contient certes pas de valeurs scalaires, mais on va utiliser ses éléments constitutifs (que l'on peut voir comme des opérations en attente d'argument — des opérateurs différentiels) très exactement comme on aurait utilisé les valeurs scalaires composant un vecteur.

La notation nabla fournit un moyen commode pour exprimer les opérateurs vectoriels en coordonnées cartésiennes ; dans d'autres systèmes de coordonnées, elle est encore utilisable au prix de précautions supplémentaires ; pour plus de précisions, et des interprétations plus théoriques (en particulier la relation avec la dérivée covariante), voir les articles détaillés nabla et connexion de Koszul.

Opérateur différentiel gradient[modifier | modifier le code]

Article détaillé : gradient.

Le gradient est un opérateur qui s'applique à un champ de scalaires et décrit un champ de vecteurs qui représente la variation de la valeur du champ scalaire dans l'espace . Pratiquement, le gradient indique la direction de la plus grande variation du champ scalaire, et l'intensité de cette variation. Par exemple, le gradient de l'altitude est dirigé selon la ligne de plus grande pente et sa norme augmente avec la pente.

En mathématiques, le gradient du champ f, supposé continûment différentiable, au point a, est défini par la relation

\mathrm d f(a)\cdot h = \left(\overrightarrow{\mathrm{grad}}_a f\right) \cdot h,

\mathrm d f(a). h désigne la valeur sur le vecteur h de la différentielle de la fonction f au point a.

C'est donc tout simplement la définition de l'application linéaire tangente du champ scalaire f(M)=f(x,y,z) en M=a. De plus, pour une surface d'équation f(x,y,z)=0, le vecteur normal à la surface au point a=(x_a,y_a,z_a) est donné par \overrightarrow{\mathrm{grad}}_a f, ce qui se déduit facilement de ce qui précède.

Il en résulte immédiatement que la dérivée de la fonction en a par rapport au vecteur v est donnée par

\overrightarrow{\mathrm{grad}}_a f \cdot v.

En dimension 3 et coordonnées cartésiennes, le champ de gradients vérifie (en base orthonormale)

\overrightarrow{\mathrm{grad}} f = \vec\nabla f =
\begin{pmatrix}
\frac{\partial f}{\partial x} \\
\frac{\partial f}{\partial y} \\
\frac{\partial f}{\partial z}
\end{pmatrix}.

Cette relation peut servir, dans le cas particulier où elle s'applique, de définition du gradient. Elle se généralise naturellement en dimension quelconque en ajoutant des composantes au nabla.

Application linéaire tangente d'un champ de vecteurs \vec{F}(M)[modifier | modifier le code]

Soit M' le point translaté de M par la translation de vecteur \vec{h} ; alors :

\vec{F(}M') - \vec{F}(M) = (\widehat{\partial \vec{F}})_M \cdot \vec{h}  + o(\|\vec{h}\|)

définit l'opérateur linéaire noté par un chapeau pour signifier que sa représentation dans une base est une matrice carrée [3-3], application linéaire tangente du champ de vecteurs F(M).

Le déterminant de cet opérateur est le Jacobien de la transformation qui à M associe F(M).

Sa trace définira la divergence du champ de vecteurs F(M).

Cela permettra de donner du rotationnel du champ de vecteurs F(M) une définition intrinsèque.

On pourra vérifier que symboliquement :

(\widehat{\partial \vec{F}})_M \cdot \vec{h} = (\vec{h}\cdot \vec{\nabla})\vec{F}

Opérateur divergence[modifier | modifier le code]

Article détaillé : divergence (analyse vectorielle).

La divergence s'applique à un champ de tenseurs d'ordre n et le transforme en un champ de tenseurs d'ordre n-1. Pratiquement, la divergence d'un champ de vecteurs exprime sa tendance à fluer localement hors d'un petit volume entourant le point M où est calculée la divergence.

En dimension 3 et en coordonnées cartésiennes, si  \vec{F} est un tenseur d'ordre 1, alors c'est un vecteur et on peut définir la divergence par la relation


\mathrm{div} \vec F = \vec \nabla \cdot \vec{F}
= \frac {\partial F_x} {\partial x} +
\frac {\partial F_y} {\partial y} +
\frac {\partial F_z} {\partial z}

\vec{F} = (F_x, F_y, F_z) désigne le champ de vecteurs auquel est appliqué l'opérateur divergence. La divergence peut être vue, formellement, comme le produit scalaire de l'opérateur nabla par le vecteur « générique » du champ auquel elle est appliquée, ce qui justifie la notation \vec\nabla\cdot. Bien entendu, cette définition se généralise naturellement en dimension quelconque.

La définition indépendante du choix de la base est :

\mathrm{div} \vec{F} = \mbox{Tr} (\widehat{\partial \vec{F}})

Une autre définition possible, plus générale mais plus difficile à formaliser, consiste à définir la divergence d'un champ de vecteurs en un point comme le flux local du champ autour de ce point.

Opérateur rotationnel[modifier | modifier le code]

Article détaillé : rotationnel.

Le rotationnel transforme un champ de vecteurs en un autre champ de vecteurs. Plus difficile à se représenter aussi précisément que le gradient et la divergence, il exprime la tendance qu'a un champ à tourner autour d'un point : sa circulation locale sur un petit lacet entourant le point M est non nulle. Par exemple :

  • dans une tornade, le vent tourne autour de l'œil du cyclone et le champ de vecteurs vitesse du vent a un rotationnel non nul autour de l'œil. Le rotationnel de ce champ de vitesse (autrement dit le champ de vorticité ou encore champ tourbillon) est d'autant plus intense que l'on est proche de l'œil.
  • le rotationnel du champ des vitesses \vec{V}(M) = \vec{\Omega}_0\wedge \vec{OM} d'un solide qui tourne à vitesse constante \vec{\Omega}_0 est constant, dirigé selon l'axe de rotation et orienté de telle sorte que la rotation ait lieu, par rapport à lui, dans le sens direct, et vaut simplement 2 \vec{\Omega}_0

Dans un espace à 3 dimensions et en coordonnées cartésiennes, on peut définir le rotationnel par la relation

{\overrightarrow{\mathrm{rot}}}\ \vec F = \vec \nabla \wedge \vec F
= \begin{pmatrix}
{\partial F_z / \partial y} - {\partial F_y / \partial z} \\ 
{\partial F_x / \partial z} - {\partial F_z / \partial x}\\ 
{\partial F_y / \partial x} - {\partial F_x / \partial y}
\end{pmatrix}

\vec{F} = (F_x, F_y, F_z) désigne le champ de vecteurs auquel est appliqué l'opérateur rotationnel. L'analogie formelle avec un produit vectoriel justifie la notation \vec\nabla\wedge.

Cela peut aussi s'écrire, par abus de notation (c'est aussi une astuce mnémotechnique), à l'aide d'un déterminant :

{\overrightarrow{\mathrm{rot}}}\ \vec F 
 = \begin{vmatrix}
 \vec{i} & \vec{j} & \vec{k} \\
 \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\
 F_x & F_y & F_z \end{vmatrix}

(\vec i, \vec j, \vec k) désigne la base canonique. Cette dernière expression est un peu plus compliquée que la précédente, mais elle se généralise facilement à d'autres systèmes de coordonnées.

  • Une définition intrinsèque (parmi d'autres) du rotationnel est la suivante :

À partir du champ \vec{F}, on peut construire le champ \vec{X_0} \wedge \vec{F} (où \vec{X_0} est un vecteur uniforme) dont la divergence est une forme linéaire de \vec{X_0} et donc exprimable par un produit scalaire \vec{K} \cdot \vec{X_0}, où \vec{K} est l'opposé du rotationnel de \vec{F} :

 \mathrm{div}(\vec{X_0} \wedge \vec{F}) = - \overrightarrow{\mathrm{rot}} \vec{F} \cdot \vec{X_0}

Une autre définition possible, plus générale mais plus difficile à formaliser, consiste à définir le rotationnel d'un champ de vecteurs en un point comme la circulation locale du champ autour de ce point (voir rotationnel en physique).

Opérateurs d'ordre supérieur[modifier | modifier le code]

Opérateur laplacien[modifier | modifier le code]

Le plus utilisé des opérateurs d'ordre 2 est le laplacien, du nom du mathématicien Pierre-Simon de Laplace. Le laplacien d'un champ est égal à la somme des dérivées secondes de ce champ par rapport à chacune des variables.

En dimension 3 et en coordonnées cartésiennes, il s'écrit :

\Delta=\nabla^2 =\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2}.

Cette définition a un sens aussi bien pour un champ de scalaires que pour un champ de vecteurs. On parle respectivement de laplacien scalaire et de laplacien vectoriel. Le laplacien scalaire d'un champ de scalaires est un champ de scalaires alors que le laplacien vectoriel d'un champ de vecteurs est un champ de vecteurs. Pour distinguer ce dernier, on le note parfois \operatorname{\vec{\Delta}} (afin que les novices n'oublient pas qu'il s'agit de l'opérateur  \overrightarrow{\mathrm{grad}}\ \mathrm{div} - \overrightarrow{\mathrm{rot}}\ \overrightarrow{\mathrm{rot}}) ; la notation \vec{\Delta} est plutôt à déconseiller.

L'autre notation du laplacien qui apparaît ci-dessus, \nabla^2, invite à le considérer, formellement, comme le carré scalaire de l'opérateur nabla « \nabla ».

Le laplacien apparaît dans l'écriture de plusieurs équations aux dérivées partielles qui jouent un rôle fondamental en physique.

Opérateur laplacien vectoriel[modifier | modifier le code]

Le laplacien d'un champ de vecteurs \vec A est un vecteur défini par le laplacien scalaire de chacune des composantes du champ vectoriel, ainsi en coordonnées cartésiennes, il est défini par :

\operatorname{\vec{\Delta}} \vec A = \operatorname{\vec \nabla^2} \vec A = \begin{bmatrix} \frac{\partial^2 A_x}{\partial x^2} + \frac{\partial^2 A_x}{\partial y^2} + \frac{\partial^2 A_x}{\partial z^2} \\ \frac{\partial^2 A_y}{\partial x^2} + \frac{\partial^2 A_y}{\partial y^2} + \frac{\partial^2 A_y}{\partial z^2} \\ \frac{\partial^2 A_z}{\partial x^2} + \frac{\partial^2 A_z}{\partial y^2} + \frac{\partial^2 A_z}{\partial z^2} \end{bmatrix} = \begin{bmatrix} \Delta A_x \\ \Delta A_y \\ \Delta A_z \end{bmatrix}

Le laplacien vectoriel est présent :

Quelques formules différentielles[modifier | modifier le code]

Attention : les formules suivantes sont valables à condition que certaines hypothèses soient vérifiées (la fonction scalaire dans la première formule doit être \mathcal{C}^2(\Omega), où \Omega \subset \mathbb{R}, par exemple. De même, si \vec f désigne la fonction vectorielle concernée dans la seconde formule, il faut vérifier \vec f \in \mathcal{C}^2(\Omega), \Omega \subset \mathbb{R}^n.)

  • \overrightarrow{\mathrm{rot}}(\overrightarrow{\mathrm{grad}})=\vec{0}
  • \mathrm{div}(\overrightarrow{\mathrm{rot}})=0
  • \Delta = \mathrm{div}(\overrightarrow{\mathrm{grad}}) (appliqué à un scalaire)

Formules dites de Leibniz pour les produits[modifier | modifier le code]

  • \overrightarrow{\mathrm{grad}}(\vec{X_0}\cdot \vec{B} ) = (\vec{X_0} \cdot \overrightarrow{\mathrm{grad}})\vec{B} + \vec{X_0} \wedge \overrightarrow{\mathrm{rot}}(\vec{B}) (où \vec{X_0} est un vecteur uniforme) et évidemment :
  • \overrightarrow{\mathrm{grad}}(\vec{A}\cdot \vec{B} ) = (\vec{A} \cdot \overrightarrow{\mathrm{grad}})\vec{B} + \vec{A} \wedge \overrightarrow{\mathrm{rot}}\vec{B} + (\vec{B} \cdot \overrightarrow{\mathrm{grad}})\vec{A} + \vec{B} \wedge \overrightarrow{\mathrm{rot}}\vec{A}
  • \overrightarrow{\mathrm{grad}}(\vec{F}\cdot \vec{F}) = 2 (\vec{F}\cdot\overrightarrow{\mathrm{grad}})\vec{F} + 2 \vec{F} \wedge (\overrightarrow{\mathrm{rot}}(\vec{F})) (dite de Bernoulli, en mécanique des fluides)
  • \mathrm{div}(\vec{X_0} \wedge \vec{B})= - \vec{X_0} \cdot \overrightarrow{\mathrm{rot}}(\vec{B}) (où \vec{X_0} est un vecteur uniforme, définition intrinsèque du rotationnel)
  • \mathrm{div}(\vec{A} \wedge \vec{B})= - \vec{A} \cdot \overrightarrow{\mathrm{rot}}\vec{B} + \vec{B} \cdot \overrightarrow{\mathrm{rot}}\vec{A}
  • \overrightarrow{\mathrm{rot}}( \vec{X_0}\wedge \vec{B}) = \vec{X_0}\cdot \mathrm{div}\vec{B} - (\vec{X_0}\cdot \overrightarrow{\mathrm{grad}})\vec{B} (où \vec{X_0} est un vecteur uniforme, par définition de l'application linéaire tangente)
  • \overrightarrow{\mathrm{rot}}( \vec{A}\wedge \vec{B}) = \vec{A} \ \mathrm{div}\vec{B} - (\vec{A}\cdot \overrightarrow{\mathrm{grad}})\vec{B} - \vec{B} \ \mathrm{div}\vec{A} + (\vec{B}\cdot \overrightarrow{\mathrm{grad}})\vec{A}


  • \overrightarrow{\mathrm{grad}}(fg) = f \cdot\overrightarrow{\mathrm{grad}}(g)  +  g\cdot\overrightarrow{\mathrm{grad}}(f) (symétrique en f et g)
  • \mathrm{div}(\rho \cdot \vec{V}) = \rho \cdot \mathrm{div} \vec{V} + \overrightarrow{\mathrm{grad}}(\rho)\cdot  \vec{V}
  • \overrightarrow{\mathrm{rot}}(\rho \cdot \vec{V}) = \rho \cdot \overrightarrow{\mathrm{rot}} \vec{V} + \overrightarrow{\mathrm{grad}}(\rho)\wedge  \vec{V}
  • \Delta (f\cdot g) = f\cdot \Delta g + 2 \overrightarrow{\mathrm{grad}}(f) \cdot \overrightarrow{\mathrm{grad}}(g)  +g\cdot \Delta f
  • \mathrm{div} ( f \cdot \overrightarrow{\mathrm{grad}}(g) -g \cdot \overrightarrow{\mathrm{grad}}(f)) = f \Delta g - g \Delta f

Quelques formules utiles[modifier | modifier le code]

  • Soient f(M) et g(M) deux champs scalaires, il existe un champ de vecteurs \vec{A}(M) tel que :
    \overrightarrow{\mathrm{rot}}\vec{A} = \overrightarrow{\mathrm{grad}}f \wedge \overrightarrow{\mathrm{grad}}\,g
  • Le champ central \vec{OM}=\vec{r} joue un rôle très important en physique. Aussi convient-il de mémoriser ces quelques évidences :
    • son application linéaire tangente est la matrice identité (cf. la définition !),
    • donc \mathrm{div}\vec{r}=3 et \overrightarrow{\mathrm{rot}}(\vec{X_0}\wedge\vec{r})=2\vec{X_0} (où \vec{X_0} est un vecteur uniforme) et \overrightarrow{\mathrm{rot}}(\vec{r})=\vec{0}
  • D'autre part \vec{X_0}=\overrightarrow{\mathrm{grad}}(\vec{X_0}\cdot\vec{r}) (où \vec{X_0} est un vecteur uniforme).
  • Et aussi : \overrightarrow{\mathrm{grad}}f(r)=f'(r)\vec{u} avec \vec{u}=\frac{\vec{r}}{r}en particulier \overrightarrow{\mathrm{grad}}(r^2)=2\vec{r} (évident car d(\vec{r}\cdot\vec{r})=d(r^2))
  • \Delta f(r)=f''(r)+\frac{2}{r}\cdot f'(r) , sauf en r=0
  • Le champ newtonien, soit \frac{\vec{r}}{r^3}, est très souvent étudié car c'est le seul champ central à divergence nulle (évident si l'on pense en termes de flux) sauf pour r=0, où elle vaut 4\pi\cdot\delta(r) ; ce résultat est le théorème de Gauss pour l'angle solide). Il en résulte que \Delta(1/r) = - 4 \pi \cdot \delta(r). Donc
    \Delta(\vec{X_0}/r) = - 4 \pi \cdot \vec{X_0}\cdot \delta(r)
    (où \vec{X_0} est un vecteur uniforme) qui se décompose en :
    \overrightarrow{\mathrm{grad}}(\mathrm{div})(\vec{X_0}/r) = - 4 \pi \cdot \vec{X_0}\cdot \delta(r)\cdot(1/3)

(où \vec{X_0} est un vecteur uniforme), et

\overrightarrow{\mathrm{rot}}(\overrightarrow{\mathrm{rot}})(\vec{X_0}/r) = + 4 \pi \cdot \vec{X_0}\cdot \delta(r)\cdot(2/3)

(où \vec{X_0} est un vecteur uniforme) ce qui est moins évident (cf. moment magnétique).

Expressions des opérateurs en différentes coordonnées[modifier | modifier le code]

Coordonnées cylindriques[modifier | modifier le code]

\overrightarrow{\mathrm{grad}}f=\frac{\partial f}{\partial r}\vec{u_r}+\frac{1}{r}\frac{\partial f}{\partial \theta}\vec{u_\theta}+\frac{\partial f}{\partial z}\vec{u_z}
\mathrm{div}\vec{A}=\frac{1}{r}\frac{\partial}{\partial r}\left(rA_r \right)+\frac{1}{r}\frac{\partial A_\theta}{\partial \theta}+\frac{\partial A_z}{\partial z}
\overrightarrow{\mathrm{rot}}(\vec{A})=\left(\frac{1}{r}\frac{\partial A_z}{\partial \theta}-\frac{\partial A_\theta}{\partial z}\right)\vec{u_r} + \left(\frac{\partial A_r}{\partial z}-\frac{\partial A_z}{\partial r}\right)\vec{u_\theta} + \frac{1}{r}\left(\frac{\partial}{\partial r}(rA_\theta)-\frac{\partial A_r}{\partial \theta}\right)\vec{u_z}
\Delta f=\frac{1}{r}\frac{\partial}{\partial r}\left(r\frac{\partial f}{\partial r}\right) + \frac{1}{r^2}\frac{\partial^2 f}{\partial \theta^2} + \frac{\partial^2 f}{\partial z^2}

Coordonnées sphériques[modifier | modifier le code]

\overrightarrow{\mathrm{grad}}f
 =   \frac{\partial f}{\partial r}\vec{u_r}
   + \frac{1}{r}\frac{\partial f}{\partial \theta}\vec{u_\theta}
   + \frac{1}{r \sin\theta}\frac{\partial f}{\partial \varphi} \vec{u_\varphi}
\mathrm{div}\vec{A}
 =   \frac{1}{r^2}\frac{\partial}{\partial r}(r^2A_r)
   + \frac{1}{r\sin\theta}\frac{\partial} {\partial \theta}(\sin\theta A_\theta)
   + \frac{1}{r\sin\theta}\frac{\partial A_\varphi}{\partial \varphi}
\overrightarrow{\mathrm{rot}}(\vec{A})
 =   \frac{1}{r\sin\theta}\left(\frac{\partial}{\partial \theta}(\sin\theta A_\varphi)-\frac{\partial A_\theta}{\partial \varphi}\right)\vec{u_r}
   + \left(\frac{1}{r\sin\theta}\frac{\partial A_r}{\partial \varphi}-\frac{1}{r}\frac{\partial}{\partial r}(rA_\varphi)\right)\vec{u_\theta}
   + \frac{1}{r}\left(\frac{\partial}{\partial r}(rA_\theta)-\frac{\partial A_r}{\partial \theta}\right)\vec{u_\varphi}
\Delta f
 =   \frac{1}{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial f}{\partial r}\right)
   + \frac{1}{r^2 \sin\theta}\frac{\partial}{\partial \theta}\left(\sin \theta\frac{\partial f}{\partial \theta}\right)
   + \frac{1}{r^2\sin^2\theta}\frac{\partial^2 f}{\partial \varphi^2}

Annexes[modifier | modifier le code]

Sur les autres projets Wikimedia :

Bibliographie[modifier | modifier le code]

  • Max Abraham et Paul Langevin Analyse Vectorielle dans Encyclopédie des sciences mathématiques pures et appliquées. Tome IV. Cinquième volume, fascicule 1 Jules Molk (éd.) p. 12, (Gauthier-Villars, Paris, 1912-1914)

Articles connexes[modifier | modifier le code]