ATLAS Ouvrir dans l'atlas →

Théorème 6Cohomologie de de Rham

Formule de Stokes pour les chaînes cubiques

Énoncé

Soit UU un ouvert de ℝn{\mathbb R^n}, k⩾1k \geqslant 1 et ω\omega une forme différentielle de degré k−1k - 1 sur UU. Pour toute chaîne cubique z∈Qk(U)z \in Q_k(U), ∫zdω=∫∂zω.\int_z d\omega = \int_{\partial z} \omega. C'est la formule de Stokes. Pour k=1k = 1 et un cube cc, elle s'écrit ∫cdf=f(c(1))−f(c(0))\int_c df = f(c(1)) - f(c(0)).

Démonstration

Le cube unité. Soit WW un ouvert de ℝk\mathbb R^k contenant IkI^k, η\eta une forme de degré k−1k - 1 sur WW, et 1k\mathbf 1_k le cube de WW donné par l'inclusion de IkI^k (Cubes singuliers). Montrons que ∫1kdη=∫∂1kη\int_{\mathbf 1_k} d\eta = \int_{\partial \mathbf 1_k} \eta.

Par Algèbre extérieure, (i), η\eta s'écrit ∑i=1kfi dx1∧⋯∧dxi^∧⋯∧dxk\sum_{i=1}^k f_i\,dx_1 \wedge \cdots \wedge \widehat{dx_i} \wedge \cdots \wedge dx_k, où le chapeau signale le facteur omis, avec des fif_i de classe C∞C^\infty sur WW. Les deux membres étant linéaires en η\eta (Formes différentielles, (i), et Intégrale sur une chaîne, (ii)), il suffit de traiter η=f dx1∧⋯∧dxi^∧⋯∧dxk\eta = f\,dx_1 \wedge \cdots \wedge \widehat{dx_i} \wedge \cdots \wedge dx_k.

Premier membre. Dans dη=∑j∂jf dxj∧dx1∧⋯∧dxi^∧⋯∧dxkd\eta = \sum_j \partial_j f\,dx_j \wedge dx_1 \wedge \cdots \wedge \widehat{dx_i} \wedge \cdots \wedge dx_k, seul le terme j=ij = i subsiste, les autres contenant deux fois le même dxjdx_j. Ramener dxidx_i à sa place coûte i−1i - 1 transpositions : dη=(−1)i−1 ∂if dx1∧⋯∧dxk.d\eta = (-1)^{i-1}\,\partial_i f\,dx_1 \wedge \cdots \wedge dx_k.

L'image réciproque par l'inclusion ne change pas le coefficient, et le nœud Fubini pour les fonctions continues, (iii), permet d'intégrer d'abord en xix_i.

Pour yy fixé, la fonction t↦∫0t∂if(ιi,s(y)) ds−f(ιi,t(y))t \mapsto \int_0^t \partial_i f(\iota_{i,s}(y))\,ds - f(\iota_{i,t}(y)), avec la notation ιi,s(y)\iota_{i,s}(y) du nœud Fubini pour les fonctions continues, (iii), a une dérivée nulle (Théorème fondamental), donc est constante (Signe de la dérivée, (ii)). Ainsi ∫01∂if(ιi,s(y)) ds=f(ιi,1(y))−f(ιi,0(y)),\begin{aligned} &\int_0^1 \partial_i f\big(\iota_{i,s}(y)\big)\,ds \\ &\quad = f\big(\iota_{i,1}(y)\big) - f\big(\iota_{i,0}(y)\big), \end{aligned} d'où ∫1kdη=(−1)i−1∫Ik−1(f∘ιi,1−f∘ιi,0).\int_{\mathbf 1_k} d\eta = (-1)^{i-1}\int_{I^{k-1}} \big(f \circ \iota_{i,1} - f \circ \iota_{i,0}\big).

Second membre. La face (j,ε)(j, \varepsilon) de 1k\mathbf 1_k est prolongée par ιj,ε\iota_{j,\varepsilon}, et ιj,ε∗(dxl)=d(ιj,ε)l\iota_{j,\varepsilon}^*(dx_l) = d(\iota_{j,\varepsilon})_l (Formes différentielles, (iii)) vaut dtldt_l si l<jl < j, 00 si l=jl = j, et dtl−1dt_{l-1} si l>jl > j.

Si j≠ij \neq i, le facteur dxjdx_j figure dans η\eta, dont l'image réciproque est donc nulle. Si j=ij = i, il ne reste que ιi,ε∗η=(f∘ιi,ε) dt1∧⋯∧dtk−1.\iota_{i,\varepsilon}^*\eta = (f \circ \iota_{i,\varepsilon})\,dt_1 \wedge \cdots \wedge dt_{k-1}. Par conséquent ∫∂1kη=∑ε=01(−1)i+ε∫Ik−1f∘ιi,ε=(−1)i+1∫Ik−1(f∘ιi,1−f∘ιi,0),\begin{aligned} \int_{\partial\mathbf 1_k} \eta &= \sum_{\varepsilon = 0}^1 (-1)^{i+\varepsilon}\int_{I^{k-1}} f \circ \iota_{i,\varepsilon} \\ &= (-1)^{i+1}\int_{I^{k-1}} \big(f \circ \iota_{i,1} - f \circ \iota_{i,0}\big), \end{aligned} qui est le premier membre, puisque (−1)i+1=(−1)i−1(-1)^{i+1} = (-1)^{i-1}. Pour k=1k = 1, η=f\eta = f est une fonction, les faces sont les points 00 et 11, et les deux membres valent f(1)−f(0)f(1) - f(0).

Cas général. Soit cc un cube singulier lisse de dimension kk de UU, prolongé par c~:W→U\tilde c : W \to U. Par Cubes singuliers, (iii), c=c~#(1k)c = \tilde c_\#(\mathbf 1_k) et ∂c=c~#(∂1k)\partial c = \tilde c_\#(\partial \mathbf 1_k).

Par Intégrale sur une chaîne, (iv), et la commutation c~∗d=d c~∗\tilde c^* d = d\,\tilde c^* (Formes différentielles, (iii)), ∫cdω=∫1kc~∗(dω)=∫1kd(c~∗ω)=∫∂1kc~∗ω=∫∂cω.\begin{aligned} \int_c d\omega &= \int_{\mathbf 1_k} \tilde c^*(d\omega) = \int_{\mathbf 1_k} d(\tilde c^*\omega) \\ &= \int_{\partial\mathbf 1_k} \tilde c^*\omega = \int_{\partial c} \omega. \end{aligned} Les deux membres de la formule sont additifs en zz, et ∂\partial est un morphisme : elle vaut pour toutes les chaînes.

Aucun choix n'intervient.

□
Remarque

Le cube singulier c(u,v,w)=(u,uv,uvw)c(u, v, w) = (u, uv, uvw) de ℝ3\mathbb R^3 envoie I3I^3 sur le tétraèdre {0⩽z⩽y⩽x⩽1}\{0 \leqslant z \leqslant y \leqslant x \leqslant 1\} ; deux de ses faces sont dégénérées, la face u=0u = 0 écrasée sur l'origine et la face v=0v = 0 sur un segment de l'axe des xx.

Pour ω=x dy∧dz\omega = x\,dy \wedge dz, on a dω=dx∧dy∧dzd\omega = dx \wedge dy \wedge dz, et (i) du nœud Intégrale sur une chaîne donne g=det⁡dc=u2vg = \det dc = u^2v, d'où ∫cdω=∫I3u2v=16\int_c d\omega = \int_{I^3} u^2v = \frac16, qui est le volume du tétraèdre.

Au bord, seules contribuent la face u=1u = 1, partie du plan x=1x = 1, pour +12+\frac12, et la face v=1v = 1, partie du plan y=xy = x, pour −13-\frac13 ; les faces w=0w = 0 et w=1w = 1, contenues dans les plans z=0z = 0 et z=yz = y, où dy∧dzdy \wedge dz a une image réciproque nulle, et les deux faces dégénérées donnent 00. On retrouve bien 12−13=16\frac12 - \frac13 = \frac16 : c'est le flux du champ (x,0,0)(x, 0, 0) à travers la surface du tétraèdre, qui entre par la face oblique et sort par la face x=1x = 1. Ces nombres ont été vérifiés par un calcul formel.

0 (1, 0, 0) (1, 1, 0) (1, 1, 1) x = 1 +1/2 y = x −1/3 x y z
Le tétraèdre image du cube par c(u,v,w)=(u,uv,uvw)c(u, v, w) = (u, uv, uvw). Les faces x=1x = 1, en bleu, et y=xy = x, en vert, sont les seules qui portent un flux du champ (x,0,0)(x, 0, 0) : +12+\frac12 et −13-\frac13. La face u=0u = 0 du cube est écrasée sur l'origine, la face v=0v = 0 sur le segment rouge de l'axe des xx ; ces faces dégénérées ne contribuent pas.
Remarque

La démonstration ne contient qu'une idée : sur le cube unité, dηd\eta ne fait intervenir que la dérivée de chaque coefficient dans la direction manquante, et le théorème fondamental de l'analyse, appliqué dans cette direction après Fubini, transforme son intégrale en une différence entre deux faces opposées. Le cas d'un cube singulier s'y ramène par image réciproque, parce que dd commute à c~∗\tilde c^*. C'est la présentation de Spivak (Calculus on Manifolds, 1965, chapitre 4), qui intègre sur des cubes plutôt que sur des simplexes, parce que Fubini s'y applique directement.

Remarque

La formule porte le nom de Stokes, qui l'a posée en février 1854 à l'examen du prix Smith de Cambridge ; elle figurait dans une lettre de William Thomson à Stokes du 2 juillet 1850. Volterra (1889) puis Poincaré (1899) ont énoncé la forme générale en dimension quelconque, et Élie Cartan l'a écrite avec les formes différentielles extérieures (Les systèmes différentiels extérieurs et leurs applications géométriques, 1945). Sources : V. J. Katz, The History of Stokes' Theorem, Mathematics Magazine 52 (1979), n° 3.

Remarque

La forme classique de la formule porte sur une variété compacte orientée à bord MM de dimension kk et une forme ω\omega de degré k−1k - 1 : ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M} \omega. Elle est démontrée pour k⩾2k \geqslant 2 au chapitre Intégration sur les variétés de cette théorie (Stokes sur une variété), qui construit les trois outils qu'elle demande ; pour k=1k = 1, c'est le théorème fondamental de l'analyse, l'énoncé ci-dessus pour k=1k = 1.

Le premier est la notion même de variété à bord, modelée sur le demi-espace {xk⩾0}\{x_k \geqslant 0\}, avec l'orientation induite sur le bord (Variétés à bord, Orientation induite) ; la géométrie différentielle (Variétés) ne définit que des variétés sans bord.

Le deuxième est l'intégrale d'une forme de degré maximal sur une variété orientée (Intégrale sur une variété). Définie carte par carte, elle n'est indépendante des cartes que grâce à la formule du changement de variables pour un jacobien positif ; le chapitre la démontre pour l'intégrale de Riemann des fonctions continues à support compact (Changement de variables), sans la théorie de la mesure, puis sur le demi-espace (Changement de variables au bord).

Le troisième est une partition de l'unité. Le nœud Partitions de l'unité les construit sur une variété sans bord ; la même construction, avec des cubes lus dans des cartes à valeurs dans le demi-espace, donne sans choix, au voisinage de tout compact d'une variété à bord, un nombre fini de fonctions dont la somme vaut 11 (Intégrale sur une variété, (i)).

La preuve se ramène alors, carte par carte, au calcul fait ici sur le cube unité : pour une forme à support compact dans un demi-espace, le théorème fondamental appliqué dans la direction xkx_k ne laisse que l'intégrale sur l'hyperplan xk=0x_k = 0 (Stokes sur le demi-espace).

Ouvrir dans l'atlas interactif →