# Théorème de Gauss (physique)

En physique, le théorème de Gauss relie le flux d'un champ de vecteurs sortant d'une surface fermée aux entités à l'origine du champ (charges électriques pour le champ électrique, masses pour le champ gravitationnel, etc.). Un corollaire notable du théorème est que les entités extérieures à la surface ne contribuent pas au flux.
Ce théorème, qui est en fait une application du théorème de la divergence, a été démontré indépendamment par Carl Friedrich Gauss en 1813 (pour la gravitation et dans un cas géométrique particulier), par Siméon Denis Poisson en 1824 (pour l'élasticité), par Mikhaïl Ostrogradski en 1826 (pour les flux de chaleur mais avec une démonstration générale), par George Green en 1828 (pour des cas particuliers) et par Pierre-Frédéric Sarrus en 1828 (pour les corps flottants).

## Théorème de Gauss en électrostatique
### Énoncé

Le flux du champ électrique ${\vec {E}}$ sortant d'une surface fermée $S$ est proportionnel à la charge électrique totale $Q_{\text{int}}$ contenue dans le volume $V$ délimité par cette surface. La constante de proportionnalité est $1/\varepsilon _{0}$, où $\varepsilon _{0}$ est la permittivité diélectrique du vide.
$\,\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{S}{\vec {E}}\cdot \mathrm {d} {\vec {S}}={\frac {Q_{\mathrm {int} }}{\varepsilon _{0}}}$
Le champ électrique ${\vec {E}}$ vérifie l'équation de Maxwell-Gauss
    $\,\operatorname {div} {\vec {E}}={\frac {\rho }{\varepsilon _{0}}}$,
où $\rho$ désigne la densité volumique de charge. L'intégrale sur $V$ de cette équation donne
    $\int \!\!\!\!\!\int \!\!\!\!\!\int _{V}\operatorname {div} {\vec {E}}\;\mathrm {d} V={\frac {1}{\varepsilon _{0}}}\int \!\!\!\!\!\int \!\!\!\!\!\int _{V}\rho \,\mathrm {d} V\equiv {\frac {Q_{\mathrm {int} }}{\varepsilon _{0}}}$,
où $Q_{\mathrm {int} }$ est la charge totale à l'intérieur du volume $V$. Le théorème de la divergence
    $\,\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{S}{\vec {E}}\cdot \mathrm {d} {\vec {S}}=\int \!\!\!\!\!\int \!\!\!\!\!\int _{V}\operatorname {div} {\vec {E}}\;\mathrm {d} V$
nous permet finalement d'écrire
    $\,\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{S}{\vec {E}}\cdot \mathrm {d} {\vec {S}}={\frac {Q_{\mathrm {int} }}{\varepsilon _{0}}}$.

### Utilisation
#### Méthode générale

Il est toujours possible de calculer le champ électrique généré par une distribution de charge en intégrant la loi de Coulomb sur l'espace. Cette méthode étant souvent impraticable, le théorème de Gauss nous permet de simplifier ce calcul dans les cas où la répartition de charge possède des symétries permettant de choisir une surface de Gauss commode.
La méthode générale pour trouver le champ électrique ${\vec {E}}$ en un point ${\vec {r}}$ de l'espace peut être résumée comme suit.
   * Trouver l'orientation du champ ${\vec {E}}$ au point ${\vec {r}}$ par des considérations de symétrie.
   * Choisir une surface de Gauss $\Sigma$
   ** passant par le point ${\vec {r}}$,
   ** ayant les mêmes symétries que la distribution de charges $\rho$,
   ** permettant d'exprimer facilement le flux de $E$,
   ** ne se superposant pas à une distribution surfacique ou linéique de charges (ce qui n'est pas physique).
   * Calculer le flux de $E$ à travers $\Sigma$.
   * Calculer la charge $Q_{\mathrm {int} }$ contenue dans $\Sigma$.
   * Appliquer le théorème de Gauss.
Un cas simple ou le théorème de Gauss permet de faciliter ce calcul est la boule uniformément chargée.

#### Exemple de la boule uniformément chargée

On considère une boule uniformément chargée de rayon $R$ et de charge totale $Q$. On peut la représenter par une distribution de charge
    $\rho ({\vec {r}})={\begin{cases}\rho _{0}&{\text{ si }}|{\vec {r}}|<R\\0&{\text{ sinon}}\end{cases}}\quad$, avec $\quad Q=V_{\text{boule}}\rho _{0}={\frac {4}{3}}\pi R^{3}\rho _{0}$.
En coordonnées sphériques, ${\vec {r}}=(r,\theta ,\phi )$, on remarque que le système est invariant par rotation d'angle $\theta ,\phi$ quelconques, on a donc une symétrie sphérique. Le champ électrique ${\vec {E}}$ généré par cette distribution de charges ne dépendra donc que de la distance au centre de la sphère $r$ donc ${\vec {E}}({\vec {r}})={\vec {E}}(r)$. D'autre part, pour chaque point de l'espace, toutes les coupes de la sphère passant par ce point et le centre de sphère sont équivalentes donc le champ en ce point est radial. Ainsi on a ${\vec {E}}({\vec {r}})=E(r)\cdot {\vec {e}}_{r}$. Cette symétrie nous invite à choisir pour surface de Gauss $\Sigma$ la sphère de rayon $r$.
En appliquant le théorème de Gauss à la surface $\Sigma$, on obtient
    ${\frac {1}{\varepsilon _{0}}}\int \!\!\!\!\!\int \!\!\!\!\!\int _{V}\rho \,{\text{d}}V=\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{\Sigma }{\vec {E}}\cdot {\text{d}}{\vec {S}}=E(r)\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{\Sigma }{\vec {e}}_{r}\cdot {\text{d}}{\vec {S}}=4\pi r^{2}E(r)$.
On a alors deux cas de figure :
    $\int \!\!\!\!\!\int \!\!\!\!\!\int _{V}\rho \,{\text{d}}V={\begin{cases}{\frac {4}{3}}\pi r^{3}\rho _{0}&{\text{ si }}r<R\\{\frac {4}{3}}\pi R^{3}\rho _{0}&{\text{ si }}r>R\end{cases}}$,
donnant finalement une norme pour le champ électrique
    $E(r)={\frac {\rho _{0}}{3\varepsilon _{0}}}\times {\begin{cases}r&{\text{ si }}r<R\\R^{3}/r^{2}&{\text{ si }}r>R\end{cases}}$.
Dans ce cas, on voit que la symétrie de la distribution de charge permet d'exprimer simplement l'intégrale sur le champ ${\vec {E}}$.

## Théorème de Gauss en gravitation
### Énoncé

Le flux du champ gravitationnel ${\vec {\mathcal {G}}}$ sortant d'une surface fermée $S$ est proportionnel à la masse totale $M_{\text{int}}$ contenue dans le volume $V$ délimité par cette surface. La constante de proportionnalité est $-4\pi G$, où $G$ est la constante de gravitation.
$\,\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{S}{\vec {\mathcal {G}}}\cdot \mathrm {d} {\vec {S}}=-4\pi G\,M_{\text{int}}$
Le champ gravitationnel ${\vec {\mathcal {G}}}$ vérifie l'équation locale
    $\operatorname {div} {\vec {\mathcal {G}}}=-4\pi G\rho _{\mathrm {m} }$,
où $\rho _{\mathrm {m} }$ désigne la masse volumique. L'intégrale sur $V$ de cette équation donne
    $\int \!\!\!\!\!\int \!\!\!\!\!\int _{V}\operatorname {div} {\vec {\mathcal {G}}}\;\mathrm {d} V=-4\pi G\int \!\!\!\!\!\int \!\!\!\!\!\int _{V}\rho _{\mathrm {m} }\,\mathrm {d} V\equiv -4\pi G\,M_{\mathrm {int} }$,
où $M_{\mathrm {int} }$ est la masse totale à l'intérieur du volume $V$. Le théorème de la divergence
    $\,\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{S}{\vec {\mathcal {G}}}\cdot \mathrm {d} {\vec {S}}=\int \!\!\!\!\!\int \!\!\!\!\!\int _{V}\operatorname {div} {\vec {\mathcal {G}}}\;\mathrm {d} V$
nous permet finalement d'écrire
    $\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{S}{\vec {\mathcal {G}}}\cdot \mathrm {d} {\vec {S}}=-4\pi G\,M_{\text{int}}$.
Les équations locales pour le champ électrique, $\,\operatorname {div} {\vec {E}}={\frac {\rho }{\varepsilon _{0}}}$ et pour le champ gravitationnel $\,\operatorname {div} {\vec {\mathcal {G}}}=-4\pi G\rho _{\mathrm {m} }$ peuvent respectivement être déduites de la force de Coulomb $\,{\vec {F}}_{\mathrm {C} }={\frac {q_{1}q_{2}{\vec {r}}_{12}}{4\pi \varepsilon _{0}|{\vec {r}}_{12}|^{3}}}\,$et de la force de gravitation $\,{\vec {F}}_{\mathrm {g} }=-G{\frac {m_{1}m_{2}{\vec {r}}_{12}}{|{\vec {r}}_{12}|^{3}}}\,$ qui ont la même forme fonctionnelle.
Par analogie, on assimile :
* Le champ gravitationnel ${\vec {\mathcal {G}}}$ au champ électrique ${\vec {E}}$,
* La masse $m$ à la charge électrique $q$,
* La constante gravitationnelle négative $-G$ à la constante de Coulomb $k=1/4\pi \varepsilon _{0}$.
On obtient alors à partir du théorème de Gauss en électrostatique, le théorème de Gauss en gravitation
$\int \!\!\!\!\!\!\!\subset \!\!\!\supset \!\!\!\!\!\!\!\int _{S}{\vec {\mathcal {G}}}\cdot \mathrm {d} {\vec {S}}=-4\pi G\,M_{\text{int}}$
On remarquera l'apparition du signe $-$ qui provient du fait que la force de gravitation est toujours attractive (la masse étant toujours positive) alors que la force de Coulomb peut être attractive ou répulsive en fonction du signe des charges.

### Exemples d'application
#### Gravité à l'extérieur de la Terre

La Terre pouvant être considérée comme à symétrie sphérique (en première approximation), le champ de gravité ${\vec {g}}$ est partout dirigé vers le centre O de la Terre, et son module $g$ (la gravité) ne dépend que de la distance radiale $r$ (distance à O). Le théorème de Gauss appliqué à la sphère de centre O et de rayon $r$ s'écrit $-4\pi r^{2}g=-4\pi G\,M_{\mathrm {int} }$ où $M_{\mathrm {int} }$ désigne la masse contenue à l'intérieur de la sphère, donc :
    $g=G{\frac {M_{\mathrm {int} }}{r^{2}}}$.
Si l'on se situe à l'extérieur de la Terre ($r>R\,$ où $R$ désigne le rayon terrestre), $M_{\mathrm {int} }$ n'est autre que la masse terrestre $M$ (la masse d'air contenue entre les distances $R$ et $r$ est négligeable) :
    $g={\frac {GM}{r^{2}}}$.
La gravité à l'extérieur de la Terre est donc la même que si toute la masse de la Terre était concentrée au point O.
Si l'on désigne par $g_{0}$ la gravité à la surface de la Terre ($g_{0}=GM/R^{2}$), la gravité à l'extérieur de la Terre s'écrit :
    $g=g_{0}\left({\frac {R}{r}}\right)^{\!2}$.
À l'altitude $z$ on a $\,r=R+z=R\left(1+{\frac {z}{R}}\right)$, la gravité vaut donc :
    $g=g_{0}\left(1+{\frac {z}{R}}\right)^{\!-2}$.
Si l'altitude est petite (comparée au rayon terrestre), on peut remplacer l'expression ci-dessus par un développement limité au premier ordre :
    $g\approx g_{0}\left(1-2{\frac {z}{R}}\right)$.
$R\approx$ 6 371 km. Pour $z=$ 1 km, $2{\frac {z}{R}}\approx$ 0,00031, la gravité diminue donc d'environ 0,03 % par kilomètre d'altitude.

#### Gravité à l'intérieur de la Terre

Comme dans la section précédente, le champ de gravité ${\vec {g}}$ est partout dirigé vers le centre O de la Terre, son module $g$ ne dépend que de la distance radiale $r$, et le théorème de Gauss appliqué à la sphère de centre O et de rayon $r$ donne $\,g=G\,M_{\mathrm {int} }/r^{2}\,$ où $M_{\mathrm {int} }$ désigne la masse contenue à l'intérieur de la sphère.
Quand $r$ est inférieur au rayon terrestre $R$, la masse $M_{\mathrm {int} }$ dépend de $r$. Si l'on suppose dans un premier temps que la masse volumique $\rho$ de la Terre est uniforme, $M_{\mathrm {int} }(r)={\tfrac {4}{3}}\pi r^{3}\rho \,$ donc :
    $g={\frac {4}{3}}\pi \,G\rho \,r$.
La gravité en surface étant égale à $g_{0}={\tfrac {4}{3}}\pi \,G\rho \,R$, la gravité à l'intérieur de la Terre peut s'écrire :
    $g=g_{0}{\frac {r}{R}}$.
En réalité, la masse volumique des matériaux terrestres est une fonction $\rho (r)$ décroissante, en partie en raison de leur composition chimique (le fer du noyau est plus dense que les silicates du manteau) et en partie parce qu'ils sont d'autant plus comprimés qu'ils sont profonds (la pression augmente avec la profondeur, elle est maximale au centre de la Terre). Alors $\,M_{\mathrm {int} }(r)=\int _{0}^{r}\!4\pi \,r_{1}^{2}\,\rho (r_{1})\,\mathrm {d} r_{1}$ donc :
    $g={\frac {4\pi \,G}{r^{2}}}\int _{0}^{r}\!r_{1}^{2}\,\rho (r_{1})\,\mathrm {d} r_{1}$.

#### Théories de la Terre creuse

Les théories de la Terre creuse imaginent qu'il y a dans la Terre une grande cavité centrale, sphérique (de même centre O que la Terre elle-même), et que d'éventuels habitants de la surface (concave) de cette cavité tiendraient debout avec la tête vers O et les pieds au sol, maintenus par une gravité dirigée vers le sous-sol. Le théorème de Gauss montre facilement qu'en fait la gravité serait nulle, que tout l'intérieur de la cavité centrale serait en apesanteur : en prenant comme surface de Gauss une sphère de centre O et de rayon r inférieur ou égal au rayon de la cavité, la même formule que précédemment, $g=G\,M_{\mathrm {int} }/r^{2}$, donne $\,g=0\,$ puisqu'il n'y a aucune masse à l'intérieur de la sphère.
Dans le cycle de Pellucidar d'Edgar Rice Burroughs, il y a en plus un soleil fixe au centre O. Le théorème de Gauss montre de même que la gravité ressentie par des habitants de la surface de la cavité ne serait pas dirigée vers le sous-sol mais vers le soleil fixe.

#### Correction de plateau

La correction de plateau est l'une des corrections apportées à l'accélération de la pesanteur déduite du modèle ellipsoïdal de la Terre pour définir une valeur théorique et la comparer à la valeur mesurée (la différence entre celle-ci et celle-là est l'anomalie de Bouguer). La correction de plateau est la pesanteur due à la présence de roches entre le niveau de la mer et le point de mesure, d'altitude h.
Les roches sont supposées appartenir à une couche d'épaisseur h (un plateau) présentant une différence de masse volumique $\Delta \rho$ par rapport à l'air qui devrait être présent entre les altitudes 0 et h selon le modèle ellipsoïdal (les accidents de la topographie sont pris en compte dans une autre des corrections intervenant dans le calcul de l'anomalie de Bouguer). L'étendue de la couche étant supposée très grande par rapport à h, la correction de plateau $\Delta g_{\mathrm {p} }$ est assimilée à la pesanteur (le module du champ de pesanteur $\Delta {\vec {g}}_{\mathrm {p} }$) créée par une plaque horizontale d'épaisseur h, de masse volumique $\Delta \rho$ et d'extension horizontale infinie dans les deux directions. La symétrie de cette plaque implique que le champ de pesanteur qu'elle crée soit perpendiculaire au plan médiateur de la plaque, de même module $\Delta g_{\mathrm {p} }$ en tous les points de ses deux faces, et de sens différent sur ces deux faces (dirigé vers le bas sur la face du haut, vers le haut sur la face du bas).
On choisit comme surface de Gauss une portion de cylindre de section S, de génératrices perpendiculaires au plan médiateur, et avec l'une de ses deux faces juste au-dessus de la face supérieure de la plaque et l'autre juste au-dessous de la face inférieure. Le flux de $\Delta {\vec {g}}_{\mathrm {p} }$ sortant de cette surface fermée se compose du flux sortant par les parois latérales (nul puisque $\Delta {\vec {g}}_{\mathrm {p} }$ y est parallèle à la surface), du flux sortant par la face du haut ($-S\,\Delta g_{\mathrm {p} }$) et du flux sortant par la face du bas (égal au précédent), le flux total est donc $-2S\,\Delta g_{\mathrm {p} }$. La masse contenue à l'intérieur de la surface de Gauss est égale au produit de la masse volumique $\Delta \rho$ par le volume de plaque compris à l'intérieur, c'est-à-dire $S\,h$. Le théorème de Gauss s'écrit donc $-2S\,\Delta g_{\mathrm {p} }=-4\pi G\,\Delta \rho \,S\,h$, d'où :
    $\Delta g_{\mathrm {p} }=2\pi G\,\Delta \rho \,h$.
La masse volumique de l'air étant négligeable devant celle des roches, on peut remplacer $\Delta \rho$ par $\rho$, la masse volumique de ces roches.