Divergence (analyse vectorielle)

Un article de Wikipédia, l'encyclopédie libre.
Aller à : navigation, rechercher
Page d'aide sur l'homonymie Pour les articles homonymes, voir Divergence.

En géométrie, la divergence d'un champ de vecteurs mesure le défaut de conservation du volume sous l'action du flot de ce champ.

Plus précisément, soit  t\mapsto \phi_t(x) la trajectoire du champ X issue de x. Ces trajectoires s'organisent en une famille de transformations  x\mapsto \phi_t(x) (le flot de X), et pour tout domaine D on a


\frac{\mathrm{d}}{\mathrm{d}t} vol(\phi_t(D))_{\vert t=0}  = \int_D  \operatorname{div} X \mathrm{d}\mathbf{x}

Un champ à divergence nulle est un champ qui conserve le volume, tel le champ des vecteurs vitesse d'un fluide incompressible.

Ainsi, div X est une fonction à valeurs réelles qui mesure la variation première du volume le long des trajectoires dudit champ.


En raison de son utilisation dans les calculs de flux de champ de vecteurs, la divergence intervient en physique pour exprimer des lois de conservation ainsi que pour la formulation locale des lois physiques faisant intervenir un champ suivant une loi en carré inverse de la distance. La divergence est notamment utilisée dans les équations de la mécanique des fluides ou les équations de Maxwell.

Des définitions plus précises sont données dans le corps de l'article.

Divergence d'un champ de vecteurs[modifier | modifier le code]

L'opérateur divergence est un opérateur différentiel linéaire de degré 1, défini sur les champs de vecteurs et à valeurs dans les fonctions.

Définition en dimension 3[modifier | modifier le code]

En dimension 3 et en coordonnées cartésiennes, la divergence d'un champ de vecteurs \vec{A}=\begin{pmatrix}A_x\\A_y\\A_z\end{pmatrix} a pour expression[1]

\operatorname{div}\vec A = \frac{\part A_x }{\part x }+\frac{\part A_y }{\part y }+\frac{\part A_z }{\part z }\qquad(0).

Formellement, l'opérateur divergence appliqué à un champ vectoriel \vec A peut s'interpréter comme produit scalaire du vecteur nabla \vec\nabla par le vecteur \vec A.

\vec\nabla \cdot \vec A = \operatorname{div}\vec A = \frac{\part{A_x}}{\part{x}}+\frac{\part{A_y}}{\part{y}}+\frac{\part{A_z}}{\part{z}}

Cette définition a le désavantage d'être dépendante du choix d'une base orthonormée.

Interprétation heuristique en termes de variation du volume[modifier | modifier le code]

Le flot du champ A est approximativement (pour t petit) donné par


(x,y,z)\mapsto F_t(x,y,z) =(x,y,z)+ tA(x,y,z)

Le volume de l'image par  F_t d'un "petit" pavé de centre (x,y,z) est multiplié par le déterminant de la matrice jacobienne de  F_t. Mais la dérivée par rapport à t de ce déterminant, pour t=0 est précisément


\frac{\part A_x }{\part x }+\frac{\part A_y }{\part y }+\frac{\part A_z }{\part z }.

Justification rigoureuse de l'interprétation[modifier | modifier le code]

Elle utilise les formes différentielles. Dans \mathbb{R}^3 le volume d'un domaine s'obtient en intégrant sur ce domaine la forme différentielle dx\wedge dy \wedge dz . Soit \phi_t le flot du champ A. La variation infinitésimale du volume est alors


\frac{d}{dt}\phi_t^\ast (dx\wedge dy \wedge dz)_{\mid t=0}

C'est par définition la dérivée de Lie L_A(dx\wedge dy \wedge dz) de dx\wedge dy \wedge dz . On vérifie que


L_A(dx\wedge dy \wedge dz)= \left(\frac{\part A_x }{\part x }+\frac{\part A_y }{\part y }+\frac{\part A_z }{\part z}\right)( dx\wedge dy \wedge dz)

(voir plus bas la démonstration -en dimension n'- de cette identité).

En particulier, le flot de A conserve le volume (c’est-à-dire \mathrm{vol}(\varphi_t(D)) =\mathrm{vol}(D) pour tout domaine D de  \mathbb{R}^3  ) si et seulement si la divergence est partout nulle. Le volume augmente si la divergence est positive, diminue si elle est négative.

Ce calcul montre aussi que la divergence ne dépend pas de la structure euclidienne de l'espace, mais seulement de l'élément de volume.

Exemples[modifier | modifier le code]

La divergence du champ x\frac{\part}{\part x }+y\frac{\part}{\part y }+z\frac{\part} {\part z} (champ radial) est constante (égale à 3) (le flot de ce champ est formé d'homothéties, qui multiplient les volumes par une constante).

La divergence du champ y\frac{\part}{\part x }-x \frac{\part}{\part y } est nulle. (son flot est formé de rotations)


Interprétation en termes de flux[modifier | modifier le code]

La divergence peut être définie en termes de flux d'un champ de vecteur. Si D est un domaine relativement compact de U, dont le bord est une surface lisse S, le flux de X à travers S est égal à l'intégrale sur D de la divergence, d'après le théorème de Stokes. Explicitement,

\int_D(\operatorname{div} X)dxdydz=\int_{S} i(X)(dx\wedge dy\wedge dz)=\int_{S} g(X,\nu) \Omega_S

Dans la dernière intégrale, \nu est le vecteur unitaire normal sortant de S, et \Omega_S est l'élément d'aire de la surface S. Cette égalité est valable en toute dimension, et s'etend aux domaines à bord des variétés riemanniennes orientées[2]. Elle est connue sous le nom de théorème de Green-Ostrogradski ou théorème de flux-divergence.

Divergence en dimension n[modifier | modifier le code]

Cette définition et ses propriétés s'étendent aux champs de vecteurs sur \mathbb{R}^n. Si  A = (A_{x^i})_{1\le i\le n} est un tel champ de vecteurs, on pose


\mathrm{div} A =\sum_{i=1}^n \frac{\part A_{x^i}}{\part x^i}

Pour \omega = dx^1\wedge\dots dx^n, on a encore L_A\omega=(\mathrm{div}A)\omega, d'où l'interprétation en termes de (non)conservation du volume n-dimensionel.



Exemple. Pour le champ linéaire donné par


 A_{x^i}= \sum_{j=1}^n a^i_jx^j

la divergence est la trace de la matrice (a^i_j)_{1\le i,j\le n}.

Propriétés[modifier | modifier le code]

  • \operatorname{div}\bigl(\overrightarrow{\operatorname{rot}}\vec A\bigr) = 0

Cette formule est particulière à la dimension 3. Elle signifie qu'un champ rotationnel est à divergence nulle. Inversement, si un champ de vecteurs \vec B sur un ouvert étoilé de \R^3 est à divergence nulle, il existe un champ \vec A tel que

\vec B=\overrightarrow{\operatorname{rot}}\vec A

(on dit alors que \vec A est un potentiel vecteur). Cette propriété, une fois convenablement interprétée en termes de formes différentielles, est une application directe du lemme de Poincaré.

Attention. Le champ newtonien \tfrac{\vec r}{r^3} est à divergence nulle, mais il n'existe pas de champ de vecteurs \vec A tel que \overrightarrow{\operatorname{rot}}\vec A=\tfrac{\vec r}{r^3}. En effet, si tel était le cas, son flux à travers toute surface fermée serait nul, alors que son flux à travers les sphères centrées à l'origine vaut 4\pi. En fait, ce champ n'est défini que sur l'espace privé de l'origine, qui n'est pas un ouvert étoilé : le lemme de Poincaré ne s'applique pas.

  • \operatorname{div}(f\vec A)=f\operatorname{div}(\vec A)+\nabla f\cdot \vec A

Cette formule, qui est une conséquence directe de la formule de Leibniz, permet de voir l'opérateur divergence comme le transposé (au signe près) de l'opérateur gradient. En effet, d'après le théorème de Stokes, l'intégrale sur \R^n de la divergence d'un champ de vecteurs nul en dehors d'une partie bornée est nulle. Par conséquent, si f est une fonction lisse et \vec A un champ de vecteurs, tous deux nuls en dehors d'une partie bornée (cette condition assurant que les intégrales ont un sens),

\int_{\R^n}f\operatorname{div} \vec A\;\mathrm dx_1\ldots \mathrm dx_n= - \int_{\R^n}\nabla f\cdot \vec A\;\mathrm dx_1\ldots \mathrm dx_n
.

Cette propriété s'interprète de la façon suivante. Soient C^\infty(\R^n) et \operatorname{Vect}(\R^n) respectivement les espaces vectoriels des fonctions lisses et des champs de vecteurs sur \R^n à support compact. On les munit des produits scalaires

\langle f,g\rangle \doteqdot \int_{\R^n}fg A\; \mathrm dx_1\ldots \mathrm dx_n et  \langle \vec A,\vec B \rangle \doteqdot \int_{\R^n}\vec A\cdot\vec B\; \mathrm dx_1\ldots \mathrm dx_n

Alors \langle \nabla f,\vec A\rangle = -\langle f,\operatorname{div}\vec A\rangle
, ce qui permet de voir l'opérateur divergence comme le transposé (au signe près) de l'opérateur gradient.

Cette interprétation de la divergence présente l'avantage de se généraliser aussi bien aux variétés riemanniennes qu'aux tenseurs.

  • \operatorname{div}(\vec A\wedge\vec B)=\vec B\cdot\operatorname{rot}\vec A - \vec A\cdot\operatorname{rot}\vec B

Une application typique de cette formule est le théorème de Poynting en électromagnétisme.

  • \overrightarrow{\operatorname{rot}}\Bigl(\overrightarrow{\operatorname{rot}}\bigl(\vec A\bigr)\Bigr) = \overrightarrow{\operatorname{grad}}\Bigl(\operatorname{div}\bigl(\vec A\bigr)\Bigr)-\Delta \bigl(\vec{A} \bigr)

Ces relations, très utilisées en analyse vectorielle, se comprennent mieux dans le cadre des formes différentielles.

Utilisation en physique[modifier | modifier le code]

Lois de conservation[modifier | modifier le code]

Illustration représentant un volume V intérieur à la surface fermée (S) dont la normale extérieure en tout point est \vec n.

D'une manière générale, la divergence est reliée en physique à l'expression locale de la propriété de conservation d'une grandeur. En considérant une surface fermée quelconque (S), la variation d'une grandeur conservative dans le volume fermé par cette surface est, par définition d'une grandeur conservative, due aux échanges avec l'extérieur (il n'existe pas de sources de création ou d'annihilation d'une grandeur conservative). Le bilan de cette grandeur entre deux instants s'écrit donc uniquement comme la somme du flux de cette grandeur à travers la surface fermée (S) et de la variation temporelle de la grandeur à l'intérieur de la surface (S). Si la grandeur est conservative, ce bilan est nul.

Par exemple, en électromagnétisme, si \vec J est le vecteur densité volumique courant, \rho la densité volumique de charge électrique et V le volume intérieur à la surface (S), la conservation de la charge s'écrit de façon intégrale :

\iint_{(S)} \vec{J}\cdot\vec{n}\,\mathrm dS + \frac{\mathrm d}{\mathrm dt} \iiint_{V} \rho\,\mathrm dV = 0

ou encore, pour une surface (S) fixe :

\iint_{(S)} \vec{J}\cdot\vec{n}\,\mathrm dS + \iiint_{V} \frac{\partial \rho}{\partial t}\,\mathrm dV = 0

avec \vec{n} est un vecteur unitaire normal en tout point à (S).

La formule d'Ostrogradsky permet de réécrire l'équation précédente en termes de divergence :

\iiint_{V} \operatorname{div}(\vec J)\,\mathrm dV + \iiint_{V} \frac{\partial \rho}{\partial t}\,\mathrm dV = 0

Ce qui mène immédiatement à la relation locale de conservation :

\operatorname{div}(\vec J) + \frac{\partial \rho}{\partial t} = 0

Il est ainsi également possible d'exprimer localement, par exemple dans le cadre de la mécanique des fluides, si \rho est la masse volumique en un point et \vec V le champ des vecteurs vitesse :

\operatorname{div}(\rho\,\vec V)+\frac{\partial \rho}{\partial t}=0 (équation de continuité).

D'autres lois de conservation font intervenir la divergence de tenseurs d'ordre 2, comme la conservation de la quantité de mouvement en mécanique des fluides. En relativité générale, on montre aussi la nullité de la divergence du tenseur énergie-impulsion.

Champs radiaux en carré inverse de la distance[modifier | modifier le code]

Lorsqu'une loi d'interaction radiale, due à des sources ponctuelles, varie comme le carré inverse de la distance il est possible d'établir que le flux du champ d'interaction à travers une surface fermée est toujours proportionnel à la quantité de sources présentes à l'intérieur de la surface fermée. Ce type de relation porte généralement en physique le nom de théorème de Gauss. Par exemple, dans le cas du champ électrostatique \vec E, dû aux charges électriques Q_\text{int} présentes à l'intérieur de la surface fermée (S) on a la forme intégrale du théorème de Gauss suivante :

\iint_{(S)} \vec E \cdot \vec n \,\mathrm dS = \frac{Q_\text{int}}{\varepsilon_0}

Grâce au théorème de flux-divergence il est possible d'exprimer une forme locale du théorème de Gauss. L'équation précédente se réécrit :

\iiint_V \operatorname{div} \vec E \, \mathrm dV= \frac{Q_\text{int}}{\varepsilon_0} = \iiint_V \frac{\rho}{\varepsilon_0}\,\mathrm dV

si V est le volume délimité par (S) et \rho la densité volumique de charge. On obtient alors immédiatement :

\operatorname{div} \vec E = \frac{\rho}{\varepsilon_0}

qui est la forme locale du théorème de Gauss. Ce type de relation est également possible pour le champ de gravitation :

\operatorname{div} \vec G  = - 4\pi\,G\,\rho

G est la constante fondamentale de la gravitation, \vec G le champ de gravitation et \rho la masse volumique.

Flux du champ magnétique[modifier | modifier le code]

En électromagnétisme il est possible de montrer, à partir de la loi de Biot et Savart, que la divergence du champ magnétique \vec B est nulle :

\operatorname{div}\vec B = 0

Cette propriété intrinsèque du champ magnétique permet d'établir que le flux du champ magnétique à travers une surface fermée est toujours nul ; on dit que le champ magnétique est à flux conservatif. En effet, si on appelle (S) la surface fermée considérée et V le volume intérieur à cette surface, on a :

\iint_{(S)} \vec B \cdot \vec n \,\mathrm dS = \iiint_V \operatorname{div}(\vec B)\,\mathrm dV = 0

Expression d’un champ de vecteur en dimension 3 dans d'autres systèmes de coordonnées[modifier | modifier le code]

En coordonnées cylindriques[modifier | modifier le code]

\boldsymbol \nabla \cdot \boldsymbol A =
\frac{1}{r}\frac{\partial (r\cdot A_r)}{\partial r} + \frac{1}{r}\frac{\partial A_\theta}{\partial \theta} + \frac{\partial A_z}{\partial z}

En coordonnées sphériques[modifier | modifier le code]

En choisissant pour convention les notations physiques (conformément au standard ISO 31-11), soit (x,y,z) \longrightarrow (r \cos \varphi \sin \theta, r \sin \varphi \sin \theta, r \cos \theta), 0<\theta<\pi, 0<\varphi<2\pi :

\boldsymbol \nabla \cdot \boldsymbol A =
\frac{1}{r^2}\frac{\partial (r^2 A_r)}{\partial 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}

Notes et références[modifier | modifier le code]

Notes[modifier | modifier le code]

  1. J.-P. Pérez et al., Électromagnétisme. Fondements et applications, 3e édition, Masson, Paris, 1997, pages 612-617.
  2. GallotHulinLafontaine, p. 211

Référence[modifier | modifier le code]

(en) Sylvestre Gallot, Dominique Hulin et Jacques Lafontaine, Riemannian Geometry [détail des éditions]

Voir aussi[modifier | modifier le code]

Sur les autres projets Wikimedia :

Articles connexes[modifier | modifier le code]

Bibliographie[modifier | modifier le code]