Propriétés du potentiel newtonien

Propriétés du potentiel newtonien

Cet article fort technique traitant en détail des propriétés d'un potentiel inversement proportionnel à la distance s'adresse avant tout aux étudiants et aux chercheurs en géodésie physique et en géophysique, mais peut aussi intéresser les physiciens théoriciens. Sa lecture demande une tournure d'esprit orientée vers les mathématiques appliquées.

Sommaire

Potentiel gravitationnel et force de gravitation

On appelle le potentiel gravitationnel d'une masse répartie dans le volume B d'un corps (la Terre dans notre cas) un potentiel newtonien. Ce potentiel joue un rôle tellement important en géodésie physique qu'il semble nécessaire d'en établir en détail les propriétés. A cette fin nous supposons que la densité ρ du corps attirant B est en général une fonction continue, mais qu'il peut exister un nombre fini de surfaces dans B sur lesquelles la densité est discontinue. Une telle fonction est appelée continue par morceaux. Nous ne ferons pas de distinction ici entre l'espace ponctuel et l'espace vectoriel associé, et nous écrirons pour les points quelconques P(x1,x2,x3) et Q(y1,y2,y3) aussi x et y, respectivement. Nous désignerons en général les éléments de volume par et les éléments de masse par dM. Nous avons donc la relation :

dM = ρdτ.

Par conséquent, en portant notre attention sur le point Q en particulier, nous utiliserons indistinctement dτ(Q), dy3 ou dy1dy2dy3 pour l'élément de volume en Q, et dM(Q) ou dM(y) pour l'élément de masse en Q. En notation vectorielle nous écrirons donc

V(P) = V( \vec x ) = G \int_B \frac{\rho (y)} {\| \vec y - \vec x \|} \mathrm dy^3

pour le potentiel de gravitation et

X_i(P) = X_i(\vec x) = G \int_B \frac{\rho(\vec y) (y_i-x_i)}{ \| \vec y- \vec x \|^{3}} \mathrm d y^3

pour la ie composante de la force de gravitation. Si nous posons ρ(Q) = 0 en tout point Q extérieur au corps B, c'est-à-dire dans l'espace vide, nous pouvons étendre les intégrations dans les deux expressions précédentes à tout domaine régulier D contenant B. Tout point appartenant au domaine d'intégration est appelé un point intérieur (ou point interne), tous les autres points sont appelés des points extérieurs (ou points externes).

Soit \vec r = \vec y - \vec x le vecteur joignant P à Q. Nous avons alors

ri = yixi pour i = 1,2,3.

La distance d(P,Q) entre P et Q peut s'écrire de l'une quelconque des façons suivantes

d(P,Q) = \| \vec y - \vec x \| = \sqrt{ (y_1 - x_1)^2+(y_2-x_2)^2+(y_3 - x_3)^2  } = r = \sqrt{  r_1^2+r_2^2+r_3^2  } = \sqrt{r_k r_k} ,

en sous-entendant implicitement que dans le dernier terme la sommation s'effectue pour toutes les valeurs de l'indice muet k qui se répète (convention de sommation d'Einstein).

Il est souvent avantageux d'écrire les intégrales plus haut sous la forme d'intégrale de Stieltjes, à savoir :

\ V(P) = G \int_M r^{-1} \mathrm dM,
\ X_i(P) = G \int_M r_i r^{-3} \mathrm dM,

M est la masse totale du corps B.

Existence et continuité du potentiel newtonien et de ses dérivées dans l'espace

Nous allons montrer maintenant que le potentiel V(P) et ses dérivées premières existent et sont continus dans tout l'espace pour des distributions de densité continues par morceaux. Comme le nombre de surfaces de discontinuité est fini, le corps B peut être subdivisé en un nombre fini N de corps plus petits Bk de masses Mk telles que

M = \sum_{k=1}^N M_k,

dans chacun desquels la densité est continue. On peut donc écrire

V(P) = G \int_{M} r^{-1} \mathrm{d}M = G \sum_{k=1}^N \int_{M_{k}}r^{-1} \mathrm{d}M.

dans chacun desquels la densité est continue.

Existence et continuité du potentiel newtonien et de la gravité à l'extérieur de la matière

Si le point potentié P se trouve à l'extérieur de B, chacune des intégrales

 \int_{B_k} \rho(y)|y-x|^{-1} \mathrm d^3\mathrm dy

existe, car la fonction ρ(y) | yx | − 1 est continue dans la région d'intégration.

Puisque chaque intégrale existe séparément, leur somme existe aussi, et par conséquent V existe aux points extérieurs. En outre, puisque l'intégrand dans l'expression de V(P) est partout positif de même G, on a

G\frac{M}{r_\text{min}} \leqslant V(P) \leqslant G\frac{M}{r_\text{max}}.

Ici, rmin et rmax désignent respectivement les distances minimum et maximum entre le corps attirant et le point attiré P extérieur au corps. Lorsque P s'éloigne de B, à la fois rmin et rmax grandissent indéfiniment. On voit donc que V possède la limite zéro à l'infini, quelle que soit la direction envisagée. En termes plus précis, on dit que | V | se comporte régulièrement à l'infini et décroît comme d − 1(O,P). En effet, soit Q un point intérieur quelconque. Alors la distance d(O,Q) est bornée et  | d(O,P) - d(P,Q) | \leqslant d(O,Q), de sorte que

d(O,P)V(P) = G \int_{B} \frac{\rho(Q)}{ 1 - \frac{d(O,P)-d(P,Q)}{d(O,P)} } \mathrm d\tau(Q) = G \int_{B} \rho(Q) \mathrm d\tau(Q) + O\left( \frac{1}{d(O,P)} \right).

Par des arguments tout à fait similaires il est possible de prouver que Xi existe dans l'espace extérieur au corps B et que l'intensité du champ de gravité

X = \sqrt{ {X_1}^2 + {X_2}^2 + {X_3}^2}

varie à grande distance comme d − 2(O,P). De plus, il est évident que le potentiel et l'intensité de la gravité sont continus en des points extérieurs. Par la définition-même d'une dérivée, qui dans la théorie du potentiel est le plus souvent envisagée comme une dérivée directionnelle, la dérivation sous le signe intégrale est entièrement justifiée à l'extérieur du corps matériel. Soient h la longueur d'un petit segment de droite et \mathbf x_{i}^{o} un vecteur unitaire porté par l'axe Oxi. Alors, par définition :

\partial_{i}V = G \lim_{h=0} h^{-1} \left[\int_{M} \frac{\mathrm dM}{r(\mathbf x+\zeta h\mathbf x_{i}^{o},\mathbf y)} - \int_{M} \frac{\mathrm dM}{r(\mathbf x,\mathbf y)}\right],

donc

\partial_{i}V = G \lim_{h=0} \int_{M} \left\{\partial_{i}\left[ \frac{1}{r(\mathbf x,\mathbf y)}\right] + \frac12 h \partial_{i}^2\left[ \frac{1}{r(\mathbf x+\zeta h\mathbf x_{i}^{o},\mathbf y)}\right]\right\} \mathrm dM(\mathbf y),

ζ est un certain nombre compris entre 0 et 1. Comme

\partial_{i}r^{-1} = (y_{i}-x_{i}) r^{-3}

\partial_{i}^2r^{-1} = 3 (y_{i}-x_{i})^2 r^{-5} - r^{-3}

on a

|\partial_{i}^2r^{-1}(\mathbf x+\zeta h\mathbf x_{i}^{o},\mathbf y)| < 4 L^{-3}

si L est la distance entre l'élément de masse \mathrm dM(\mathbf y) et le point P' positionné en \mathbf x+\zeta h\mathbf x_{i}^{o}. Supposons que Q' soit le point le plus proche de P' tel que Q' fasse encore partie de B. Alors d(P',Q') est la plus petite valeur possible de L, et

\left| \frac12 h \int_{M} \partial_{i}r^{-1}(\mathbf x+\zeta h\mathbf x_{i}^{o},\mathbf y) \mathrm dM(\mathbf y) \right| < \int_{M} 2 |h| \mathrm dM(\mathbf y) d^{-3}(P',Q') = 2 M |h| d^{-3}(P',Q').

Lorsque h tend vers zéro, le point P' se meut vers P. Pour cette raison, la distance minimum d(P',Q') de P' au corps peut changer. Toutefois, comme P se trouve à l'extérieur de la région B, il existe une valeur minimum de d(P',Q') non nulle lorsque | h | est suffisamment petit. Soit dmin > 0 cette distance minimum. Alors la limite

\lim_{h=0} \left| \frac12 h \int_{M} \partial_{i}r^{-1}(\mathbf x+\zeta h\mathbf x_{i}^{o},\mathbf y) \mathrm dM(\mathbf y) \right| < 2 M |h| d_\text{min}^{-3}

s'annule avec h. Nous pouvons ainsi conclure qu'en tout point extérieur le gradient du potentiel de gravité existe et, à cause des expressions ci-avant fournissant Xi et \partial_{i}r^{-1}, représente bien la gravité :

\partial_{i}V = G \int_{M} \partial_{i}r^{-1} \mathrm dM = X_{i}.

Equation de Laplace

En fait, on peut étendre cette démonstration successivement à des dérivées d'ordres supérieurs et ainsi prouver que dans une région dépourvue de matière (c'est-à-dire dans l'espace vide), le potentiel gravitationnel possède des dérivées partielles continues de tous les ordres et est analytique.

En particulier, si nous considérons les dérivées partielles d'ordre 2, nous trouvons grâce à la relation ci-avant fournissant \partial_{i}^2r^{-1} que

(\partial_{1}^2 + \partial_{2}^2 + \partial_{3}^2) r^{-1} = 3 \left((y_{1}-x_{1})^2 + (y_{2}-x_{2})^2 + (y_{3}-x_{3})^2\right) r^{-5} - 3 r^{-3} = 0.

Nous aboutissons ainsi à la célèbre équation de Laplace

(\partial_{1}^2 + \partial_{2}^2 + \partial_{3}^2) V = G \int_{M} (\partial_{1}^2 + \partial_{2}^2 + \partial_{3}^2) r^{-1} \mathrm dM = 0.

Celle-ci est une équation aux dérivées partielles du second ordre qui doit être satisfaite par le potentiel gravifique dans les régions de l'espace dépourvues de matière. Elle joue un rôle fondamental, non seulement en géodésie, mais dans pratiquement toutes les sciences qui admettent une description mathématique des phénomènes. Nous rappelant la signification physique de l'opérateur laplacien, nous savons que les fonctions qui obéissent à l'équation de Laplace, tel que le potentiel gravifique dans le vide, possèdent la propriété d'être des fonctions moyennes dans le sens que la valeur d'une telle fonction en un point est la moyenne des valeurs de cette fonction dans un voisinage suffisamment restreint de ce point. Une telle fonction est encore appelée fonction harmonique.

Existence du potentiel newtonien et de la gravité à l'intérieur de la matière

Pour établir les propriétés du potentiel gravitique à l'intérieur d'un corps matériel, en particulier à l'intérieur de la Terre, les démonstrations deviennent plus compliquées. Cela est dû au fait que pour des points intérieurs, l'intégrale

V(x) = G \int_B \dfrac{\rho(y)}{|y-x|} \mathrm d^3y

est impropre suite à la singularité qui se produit lorsque | yx | = 0. Dès lors, la dérivation sous le signe d'intégration n'est plus permise sans justification aux points intérieurs. C'est la nécessité d'une telle justification qui complique les démonstrations. En effet, même l'existence et la continuité du potentiel en des points où il existe de la matière ne sont plus des propriétés évidentes, mais demandent une preuve. Supposons donc que P est un point intérieur. Afin de simplifier les démonstrations nous allons passer, sans perte de généralité, du système de coordonnées cartésiennes envisagé jusqu'à présent à un système de coordonnées sphériques (r,θ,λ) admettant P comme origine. On obtient ainsi :

 r_1 = y_1 - x_1 = r \, \sin(\theta) \,\cos(\lambda) ,  \qquad r_2 = y_2 - x_2 = r \, \sin(\theta) \, \sin(\lambda),     \qquad r_3 = y_3 - x_3 = r \, \cos(\theta).

Le vecteur x étant fixé, l'élément de volume en Q en termes des variables r,θ,λ est d3y = d3r = r2sin(θ)drdθdλ. Ayant défini la densité comme étant nulle en dehors de B, nous pouvons écrire

V(x) = G \int_D \dfrac{\rho(x+r)}{r} \mathrm d^3r,

où le domaine d'intégration D possède une frontière sphérique de rayon L, centrée sur le point P. Nous prenons pour L la plus grande distance entre deux points arbitraires appartenant au corps matériel de dimensions finies B. Soit ρmax la valeur maximum de | ρ | dans D. Alors

 V(P) \leqslant 2 \pi G \rho_\text{max} L^2 .

Ainsi, l'intégrale impropre V(x) ci-dessus converge uniformément par rapport aux paramètres x1,x2,x3, ce qui implique que le potentiel V existe en chaque point. D'une manière tout à fait similaire, nous pouvons prouver que les intégrales impropres définissant les composantes Xi(P)(i = 1,2,3) de la gravité existent. En effet, en considérant encore des coordonnées sphériques et en gardant les mêmes notations comme ci-dessus, nous trouvons

|X_i(P) |= G \left| \int_D \frac{r_i \, \rho(x+r)}{r^3} \mathrm d^3r \right| \leqslant 4 \pi \, G \, \rho_\text{max} L.

Ceci prouve la convergence uniforme, et donc l'existence, des intégrales impropres Xi.

Continuité du potentiel newtonien et de la gravité à l'intérieur de la matière

Nous montrons ensuite que le potentiel V est continu en un point P0 de B. Ce n'est pas une restriction essentielle que de supposer que P0 est à l'intérieur de B. En effet, nous avons vu que nous pouvons agrandir B en posant ρ(Q) = 0 dans la région qu'on ajoute. Le raisonnement classique consiste alors à décomposer le domaine B en deux sous-domaines, à savoir Bb et b. La sous-région b est une petite boule centrée sur P0. Nous écrivons

V(P) = V1(P) + V2(P)

moyennant

V_1(P) = G \int_b \dfrac{\rho(Q)}{d(P,Q)} \mathrm d\tau,
V_2(P) = G \int_{B-b} \dfrac{\rho(Q)}{d(P,Q)} \mathrm d\tau.

Or, pour \varepsilon > 0 donné arbitrairement petit, nous pouvons prendre b suffisamment petit pour que

|V_1(P)| < \frac{1}{3} \varepsilon,

indépendamment de la position de P. Pour ce faire, il suffit de prendre le rayon rb de b plus petit que \sqrt{\tfrac{\varepsilon}{6\pi G \rho_\text{max}}}, puisque dans ce cas il vient

|V_1(P)| < 2\pi G \rho_\text{max} {r_b}^2.

Il s'ensuit que pour une telle sphère b on a

| V_1(P) - V_2(P_0)| < \frac{2}{3} \varepsilon .

Alors, avec b fixé, il existe un voisinage de P0 tel que, si p appartient à ce voisinage et Q appartient au sous-volume Bb, nous avons

 \left| \dfrac{1}{d(P,Q)} - \dfrac{1}{d(P_0,Q)} \right| < \dfrac{\epsilon}{3G (B-b) \rho_\text{max} } .

Ainsi, lorsque P est dans ce voisinage, nous trouvons

|V_2(P) - V_2(P_0)| \leqslant G \int_{B-b} |\rho(Q)| \, \left| \dfrac{1}{d(P,Q)} - \dfrac{1}{d(P_0,Q)} \right| \mathrm d\tau(Q) < \dfrac{1}{3} \varepsilon.

En combinant les inégalités pour V1 et V2, nous avons

| V(P) − V(P0) | < ε.

Ainsi V est continu en P0, et dès lors dans tout l'espace.

Comme précédemment, on établira la continuité de la gravité par une partition similaire de B en deux sous-régions. En effet, soit

{X_i}(P) = {X_i}^{(1)}(P) + {X_i}^{(2)}(P),

avec

{X_i}^k (P) = G \int_{b_k} \dfrac{\rho(Q) \, e_i}{d^2(P,Q)} \mathrm d\tau(Q) \quad , k=1,2.

Ici nous utilisons les notations suivantes : b1 = b, b2 = Bb, e est le vecteur unitaire \overrightarrow{PQ} :

|e_i| \leqslant \dfrac{1}{ \sqrt{ {e_1}^2 + {e_2}^2  + {e_3}^2 } } = 1.

Pour tout \varepsilon > 0 fixé arbitrairement petit, nous pouvons prendre le rayon de b plus petit que \tfrac{\varepsilon}{12\pi G \rho_\text{max} }. Alors

\left| {X_i}^{(1)}(P) \right| < \dfrac{2}{3} \varepsilon

indépendamment de la position de P, et

 \left| {X_i}^{(1)} (P) - {X_i}^{(1)}(P_0) \right| < \dfrac{2}{3} \varepsilon .

En outre, avec b fixé mais suffisamment petit et contenant les points P et P0, et avec Q appartenant au volume Bb, nous avons

 \left| \dfrac{1}{d^2(P,Q)} - \dfrac{1}{d^2(P_0,Q)} \right| < \dfrac{ \varepsilon '}{ 3 G (B-b) \rho_\text{max} }

et

 \left| {X_i}^{(2)} (P) - {X_i}^{(2)} (P_0) \right| \leqslant G \int_{B-b} { | \rho(Q)| } \left| \dfrac{1}{d^2(P,Q)} - \dfrac{1}{d^2(P_0,Q)} \right| \mathrm d\tau(Q) < \dfrac{1}{3} \varepsilon ' .

En combinant les inégalités pour Xi(1) et Xi(2), nous établissons la continuité du champ de force gravitique en tout point intérieur P0 et, par conséquent, la continuité de la gravité dans tout l'espace rempli ou non de matière.

Existence et continuité des dérivées partielles du potentiel à l'intérieur de la matière

Toutefois, il n'est pas évident sans investigations complémentaires que les composantes de force Xi sont égales aux dérivées partielles \partial_{i}V (i = 1,2,3) aux points où il y a de la matière. En effet, les conditions usuelles permettant de dériver sous le signe intégrale ne sont pas remplies par les intégrales impropres. Néanmoins, la relation X_{i} = \partial_{i}V subsiste dans le cas présent. Pour le prouver, considérons deux points P et Po intérieurs à B, occupant par rapport à l'origine O les positions \mathbf x et \mathbf x + h \mathbf x_{i}^{o}, respectivement. Soit Q un point attirant occupant la position \mathbf y. Par rapport à P, Q est situé en \mathbf r = \mathbf y - \mathbf x. Considérons maintenant l'expression

I(h) = X_{i}(P) - h^{-1} [V(P) - V(P_{o})] = G \int_{B} \left[r^{-3} r_{i} - h^{-1} (r^{-1} - r_{o}^{-1})\right] \rho(\mathbf r) d^3\mathbf r,

ri = yixi ,

r = d(P,Q) = |\mathbf y - \mathbf x| = \sqrt{r_{1}^2+r_{2}^2+r_{3}^2},

r_{o} = d(P_{o},Q) = |\mathbf y - \mathbf x - h \mathbf x^o_{i} | = \sqrt{r^2+h^2-2hr_{i}}.

Si nous arrivons à démontrer que la limite pour h tendant vers zéro de | I(h) | est plus petite qu'un nombre positif \epsilon arbitrairement petit fixé à l'avance, alors nous pouvons conclure que la dérivée partielle \partial_{i}V de V existe et est continue aux points intérieurs, tout comme pour les points extérieurs, et que sa valeur est Xi. Ainsi, prenons | h | assez petit pour que Po se trouve à l'intérieur du domaine sphérique b de rayon rb, centré sur P. La valeur rb elle-même est choisie assez petite de manière à ce que la boule b soit contenue entièrement dans le volume B. Ainsi, comme P et Po sont extérieurs au volume Bb, il résulte des propriétés du potentiel en des points extérieurs que, quelle que soit la valeur du rayon rb, la distance d(P,Po) = | h | < rb peut être prise suffisamment petite de manière à ce que

I_{2}(h) = | G \int_{B-b} [r^{-3} r_{i} - h^{-1} (r^{-1} - r_{o}^{-1})] \rho(\mathbf y) \mathrm d^3\mathbf y < \frac12 \varepsilon.

En ce qui concerne la contribution I1(h) de l'intégrale prise sur le volume sphérique b, il convient de remarquer tout d'abord que

|r_{i}| \leqslant r

et que

r^{-1} - r_{o}^{-1} = (r_{o} - r)(r r_{o})^{-1}.

Des relations exprimant les inégalités triangulaires appliquées au triangle PPoQ on déduit que

|r_{o} - r| \leqslant |h|.

D'autre part, comme

(r^{-1} - r_{o}^{-1})^2 = r^{-2} + r_{o}^{-2} - 2 r^{-1} r_{o}^{-1} > 0,

il vient

(r r_{o})^{-1} < \frac12 (r^{-2} + r_{o}^{-2}) < r^{-2} + r_{o}^{-2}.

Ainsi, on aboutit à la majoration

| I_{1}(h) | < G \int_{b} (2 r^{-2}  + r_{o}^{-2}) \rho(\mathbf y) \mathrm d^3\mathbf y.

En posant comme de coutume ρ(Q) = 0 pour tout point Q se trouvant en dehors du volume B, nous pouvons écrire

| I_{1}(h) | < 2G \int_{b} \rho r^{-2} \mathrm d\tau + G \int_{b_{o}} \rho r_{o}^{-2} \mathrm d\tau,

bo est un domaine sphérique de rayon 2rb centré sur Po. Ainsi, nous avons

| I_{1}(h) | \leqslant 16\pi G\rho_\text{max}r_{b}

Pour des valeurs de rb suffisamment petites, c'est-à-dire pour

r_{b} < \frac{\varepsilon}{32\pi G\rho_\text{max}}

cette expression est numériquement inférieure à \epsilon/2, et donc

| I(h) | \leqslant | I_{1}(h) | + | I_{2}(h) | < \varepsilon.

De cette façon, nous avons prouvé que les dérivées partielles du premier ordre du potentiel existent et représentent les composantes de la gravité. Cette proposition est tout à fait générale. Elle est valable lorsque le point attiré se trouve à l'extérieur du corps qui attire, mais elle est aussi valable lorsque le point attiré se trouve à l'intérieur de ce corps. Une autre manière de s'exprimer consiste à dire que les dérivées de V du premier ordre s'obtiennent par dérivation sous le signe d'intégration.

Existence des dérivées partielles secondes du potentiel newtonien à l'intérieur de la matière

Il en va tout autrement des dérivées du second ordre à l'intérieur du corps. Le simple fait de supposer la densité continue et bornée ne suffit plus pour garantir l'existence de ces dérivées. Cela se voit clairement si nous dérivons formellement l'expression

\partial_i V(x) = - G \int_B \dfrac{r_i \rho(y)}{r^3} \mathrm d^3y

par rapport à xj sous le signe d'intégration, et essayons de démontrer l'existence de \partial_j \partial_i V comme auparavant pour V et \partial_i V, en agrandissant le domaine B en un domaine sphérique D de rayon L centré sur P(x) et contenant B complètement. En effet, nous trouvons alors formellement

\partial_j \partial_i V(x) = -G \int_B \partial_i \left( \dfrac{r_i}{r^3} \right) \rho(y) \mathrm d^3y = -G \int_B \left( \dfrac{3 r_i r_j}{r^5} - \dfrac{\delta_{ij} }{r^3} \right) \rho(y) \mathrm d^3y

et

\left| \partial_j \partial_i V(x) \right| < \int_D 4G \dfrac{|\rho(r,\theta,\lambda)|}{r^3} r^2 \, \mathrm dr \, \mathrm d\theta \, \mathrm d\lambda .

Condition de Hölder

Ici nous ne pouvons pas simplement remplacer | ρ | par sa borne supérieure ρmax pour pouvoir affirmer que l'intégrale du membre de droite existe, car l'intégrale impropre \scriptstyle\int_{0 \to L} \frac{\mathrm dr}{r} diverge. Nous imposons pour cette raison à la densité ρ(Q) une condition en P qui fut initialement introduite en 1882 par Otto Hölder, à savoir

  • | \rho(Q) - \rho(P) | \leqslant A \; r^{\alpha}
  • pour tous les points Q tels que r = d(P,Q) \leqslant c
  • A, α, c sont des constantes positives.

On peut montrer qu'une condition de Hölder est plus forte que la condition de continuité, mais plus faible que la condition de dérivabilité si α < 1. A fortiori elle est donc plus faible que la condition d'analyticité. S'il existe une région b dans laquelle ρ(Q) obéit à une condition de Hölder en chaque point, avec les mêmes valeurs de A, α et c, alors la fonction ρ(Q) est dite remplir une condition de Hölder uniformément dans b.

Un cas évident pour lequel une condition de Hölder uniforme s'applique est celui d'un domaine b dans lequel la densité est constante, soit ρ(Q) = ρo. Supposons que b est un volume sphérique de rayon rb et montrons que dans cet exemple les dérivées secondes de V existent en effet en un point intérieur P. Sans restreindre la généralité de la démonstration, nous pouvons supposer que le centre de b coïncide avec l'origine O des axes de coordonnées, et nous prenons le vecteur x = \overrightarrow{OP} le long de l'axe Ox3. En termes des coordonnées sphériques (y,θ,λ) , où y = \sqrt{y_k \; y_k} , \theta = \arccos \left( \tfrac{x_k \; y_k}{xy} \right), avec x = \sqrt{ x_k \; x_k}, le potentiel en P engendré par toutes les masses dans b est fourni par[1]

V(P) = 2\pi G \rho_o \left( \dfrac{1}{3} x^2 - {r_b}^2 \right).

Étant une fonction scalaire, le potentiel gravifique est invariant par rapport aux changements du système de coordonnées. Par conséquent, cette formule est générale. Il s'ensuit qu'en chaque point à l'intérieur d'une masse sphérique homogène, nous avons

V = 2\pi G \rho_o \left( \dfrac{1}{3} x_k \; x_k - {r_b}^2 \right),
\partial_i V = 4\pi G \rho_o \dfrac{x_i}{3},
\partial_j \partial_i V = 4 \pi G \rho_o \dfrac{\delta_{ji}}{3}.

Nous constatons, en particulier, que toutes les six dérivées partielles d'ordre 2 de V existent et sont continues. En outre, nous trouvons que

\partial_i \partial_i V = 4\pi G \rho_o.

Nous pouvons maintenant étudier les dérivées partielles du second ordre de V en un point intérieur quelconque P du domaine fini B, pour une distribution de densité générale satisfaisant une condition de Hölder du type ci-dessus. Comme précédemment, soit b un domaine sphérique de rayon rb centré sur P. Nous supposons que le volume b est entièrement contenu dans le volume B. Considérons alors séparément le potentiel V1(P) créé en P par toutes les masses dans b, et le potentiel V2(P) créé en P par toutes les masses restantes. Puisque P est un point extérieur pour l'évaluation de V2 , ce potentiel possède des dérivées continues de tous ordres en P et y est harmonique, c'est-à-dire V2 est une solution de l'équation de Laplace en P. Ainsi, le problème est réduit à une étude de V1.

Si nous écrivons

\rho(Q) = \rho(P) + \left[ \rho(Q) - \rho(P) \right] ,

nous voyons que le potentiel d'une sphère possédant une densité continue en P est la somme de deux potentiels : le potentiel créé par la sphère si elle était remplie de matière de densité constante égale à celle du point P, et le potentiel d'une sphère dont la densité s'annule en P. Comme nous venons juste de démontrer que le potentiel d'une sphère uniforme possède des dérivées continues du second ordre, il nous reste de discuter le cas dans lequel la densité s'annule en P et y remplit une condition de Hölder. En admettant que le rayon rb de la boule b est plus petit qu'une certaine constante positive c donnée, cela signifie que

| \rho_1 | \le A \, r^{\alpha} avec r = d(P,Q) \le r_b \le c ,

ρ1 désigne la fonction ρ(Q) − ρ(P). Celle-ci engendre dans le volume b le potentiel V11 et le champ de force X1j(1) , j = 1,2,3. Définissons la quantité

J_{ij} (P) = -G \int_b \partial_i \rho_1 \partial_j \dfrac{\mathrm d^3r}{r} ,

qui devient successivement

\begin{align}
J_{ij} & = & - G \int_b \rho_1 \partial_i \left( \dfrac{r_j}{r^3} \right) \mathrm d\tau \\
& = & -G \int_b \rho_1 \left(  3 \dfrac{r_i \, r_j}{r^5} - \dfrac{\delta_{ij}}{r^3} \right) \mathrm d\tau
\end{align}.

Cette intégrale impropre, obtenue en dérivant formellement le potentiel gravifique V11 deux fois sous le signe intégrale, converge puisque

| J_{ij} | \le 4 A G \int_b \rho_1 r^{\alpha - 3} \mathrm d^3r \le 16 \pi G A \dfrac{ \rho_1}{ \alpha^{3} } \, {r_b}^{\alpha}

moyennant la condition de Hölder. Considérons maintenant au voisinage de P occupant la position x un point Po occupant la position x + hxio, de manière à ce que d(P,Po) = | h | < rb. Formons l'expression

K(h) = \dfrac{1}{h} \left[ { X_{1j} }^{(1)} (P_o) - { X_{1j} }^{(1)} (P) \right] + J_{ij}

pour h \neq 0 . Cette intégrale est aussi convergente, comme on peut le constater en utilisant le même raisonnement que pour Jij, et parce que Xj existe aux points intérieurs. Nous avons

K(h) = G \int_b \rho_1 
\left[
\dfrac{1}{h} \left( \dfrac{( r_j - h \delta_{ij} )}{ {r_o}^3 } - \dfrac{1}{ {r_j}^3 } \right)
- \left(
3 \dfrac{r_i \, r_j}{r^5} - \dfrac{ \delta_{ij} }{r^3}
\right)
\right] \mathrm d^3r

avec ro2 = r2 + h2 − 2hri. Nous voulons montrer que K(h) tend vers zéro avec h, indépendamment de la position de P. Mais pour arriver à cette fin, nous devons nous débarrasser de h dans le dénominateur du premier terme dans l'intégrale. En remarquant que

 \dfrac{1}{ {r_o}^3 } - \dfrac{1}{r^3} = \dfrac{r^2 - {r_o}^2 }{ r r_o (r + r_o)} \left( \dfrac{1}{ {r_o}^2 }  - \dfrac{1}{r^2} + \dfrac{1}{r r_o} \right)   ,

nous trouvons

 \dfrac{1}{h} \left(  \dfrac{ (r_j - h \delta_{ij}) r_o }{ {K_1}^3(h) } - \dfrac{r_j}{r^3} \right) = - \dfrac{\delta_{ij} }{ {r_o}^3 } + \dfrac{r_j}{h} \left( \dfrac{1}{ {r_o}^3 } - \dfrac{1}{r^3} \right) = - \dfrac{ \delta_{ij} }{ {r_o}^3 }  + r_j \dfrac{2 r_i - h}{r r_o (r + r_o) } \left( \dfrac{1}{ {r_o}^2 } + \dfrac{1}{r^2} + \dfrac{1}{r r_o} \right) .

Pour h = 0, cette dernière expression se réduit à - \tfrac{ \delta_{ij} }{ r^3 } + 3 \tfrac{r_i r_j}{r^5} , rendant ainsi l'intégrand de K(h) égal à zéro et, par conséquent, K(0) = 0. Si nous pouvons prouver que l'expression K(h) est continue par rapport à h en P, nous saurons qu'elle s'annule avec h, et il s'ensuivra que la dérivée de X1j(1) par rapport à xi existe en P et vaut Jij. Traçons donc une petite sphère de volume bε, de rayon rε < rb et de centre P, de sorte que | h | < rε. Comme P et Po sont tous les deux extérieurs au domaine bbε, nous pouvons prendre | h | assez petit pour que la contribution K2(h) de l'intégration sur bbε est plus petite en valeur absolue que ε / 2, où ε est un nombre positif fixé arbitrairement petit. Considérons maintenant la contribution H1(h) de l'intégration sur bε, à savoir

 K_1(h) = G \int_{ b_{\varepsilon} } \left[ 
-\dfrac{ \delta_{ij} }{ {r_o}^3 }
+ r_j \dfrac{ 2 r_i - h }{ r r_o ( r + r_o [r + r_o] ) } \left( \dfrac{1}{ {r_o}^2 } + \dfrac{1}{r^2} + \dfrac{1}{r r_o} \right)
- \left( \dfrac{3 r_i r_j}{r^5} - \dfrac{ \delta_{ij} }{r^3} \right)
\right] 
\rho_1 \mathrm d^3r .

En nous souvenant que rj et rjhδij sont les jiemes composantes (pour i fixé) des vecteurs \overrightarrow{PQ} et \overrightarrow{P_o Q}, respectivement, nous avons

  •  |r_j| \le r ,
  •  | r_j - h \delta_{ij} | \le r_o ,
  •  | 2 r_i - h| \le |r_i| + |r_i -h| \le r + r_o .

Ainsi, l'expression entre crochets est bornée en valeur absolue par  \tfrac{4}{r^3} + \tfrac{1}{r^2 r_o} + \tfrac{1}{r {r_o}^2 } + \tfrac{2}{ {r_o}^3 } . En utilisant alors la condition de Hölder, nous obtenons

K_1(h) \le G A \int_{ b_{\varepsilon} } \left( 4 r^{\alpha - 3} + \dfrac{ r^{\alpha - 2} }{r_o} + \dfrac{ r^{\alpha - 1} }{ {r_o}^2 } + \dfrac{ 2 r^{\alpha} }{ {r_o}^3 } \right) \mathrm d^3r 
\le 16 \pi G A \dfrac{ {r_{\varepsilon}}^{\alpha} }{ \alpha } + G A \int_{ b_{\varepsilon} } 
\left( \dfrac{ r^{\alpha - 2} }{r_o} + \dfrac{ r^{\alpha - 1} }{ {r_o}^2 } + \dfrac{ 2 r^{\alpha} }{ {r_o}^3 } \right) \mathrm d^3r

Nous avons le droit de supposer α < 1, parce qu'une condition de Hölder avec un exposant donné implique toujours une condition avec un exposant positif plus petit. Cependant, dans ce cas, l'intégrand restant dans l'intégrale du membre de droite devient infini dans deux occasions, à savoir pour r = 0 et pour ro = 0. En considérant séparément les cas r < ro et ro < r, nous avons certainement

\int_{ b_{\varepsilon} } \left( \dfrac{ r^{\alpha - 2} }{r_o} + \dfrac{ r^{\alpha - 1} }{ {r_o}^2 } + \dfrac{ 2 r^{\alpha} }{ {r_o}^3 } \right) \mathrm d^3r 
< \int_{ b_{\varepsilon} } 4 r^{\alpha-3} \mathrm d^3r
+ \int_{ {b'}_{\varepsilon} } 4 {r_o}^{\alpha-3} \mathrm d^3r_o ,

b'ε est un domaine sphérique de rayon 2rε centré sur Po. Il s'ensuit que

 |K_1(h)| < 16\pi G ( 2 + 2^{\alpha} ) A \dfrac{ { r_{\varepsilon} }^{\alpha} }{\alpha} .

En choisissant rε suffisamment petit pour que | K1(h) | < ε / 2 , nous avons

| K(h) | < | K1(h) | + | K1(h) | < ε.

Ceci prouve que la fonction K(h) est continue par rapport à h en P, indépendamment de la position de P. Donc, l'existence des dérivées partielles du second ordre de V est établie et, en outre :

\partial_i \partial_j V(P) = - \partial_i X_j(P) = - \partial_j X_i(P).

En particulier, le potentiel V11 pour lequel la densité vérifie une condition Hölder et qui s'annule en P est harmonique, puisque

\partial_i \partial_i V_{11} = J_{ii} = G \int_b \left( \dfrac{3 r^2}{r^5} - \dfrac{3}{r^3} \right) \rho_1 \mathrm d^3r \equiv 0.

Équation de Poisson

Si nous superposons à la distribution de densité ρ1 s'annulant en P une autre distribution de densité ρ(P) constante partout dans le volume sphérique b, nous aboutissons au résultat suivant valable pour une distribution de densité continue dans un domaine sphérique satisfaisant à une condition de Hölder en P : les dérivées de la gravité existent, et

\nabla^2 V(P) = -4 \pi G \rho(P) ,

\nabla^2 désigne l'opérateur laplacien (ou simplement le laplacien) \partial_i \partial_i. Finalement, si nous additionnons les potentiels des distributions extérieures à b, rien n'est contribué au laplacien et la même équation ci-dessus reste valable. Celle-ci est une équation aux dérivées partielles du second ordre, qui renferme l'équation de Laplace comme un cas spécial. Elle est connue comme équation de Poisson. Elle nous permet de trouver le potentiel si nous connaissons la densité et, réciproquement, trouver une distribution de densité si nous connaissons le potentiel. Elle est donc d'une importance capitale en géodésie physique et, plus précisément, pour le problème qui consiste à trouver la forme d'équilibre d'un corps en rotation. C'est pour cette raison que nous avons pris la peine dans cet article d'établir en détail l'équation de Poisson et d'autres propriétés importantes du potentiel newtonien.

Il convient de remarquer que nous avons pas discuté ici la situation prévalant en des points de la frontière du corps. En ces points-frontière, le potentiel d'une distribution de masse volumique, ainsi que ses dérivées premières, sont continues. Par contre, sur une frontière les dérivées secondes du potentiel n'existent pas en général. Il est clair qu'elles ne peuvent pas toutes être continues, car lorsque nous passons d'un point extérieur à un point intérieur au travers d'une frontière où la densité ne s'annule pas, \nabla^2 V subit un saut de − 4πGρ. Une situation semblable est réalisée sur des surfaces à l'intérieur d'un corps sur lesquelles la densité est discontinue. Alors, la valeur de \nabla^2 V saute de − 4πGΔρ lorsqu'on passe à travers la surface de discontinuité de densité, où Δρ désigne le saut de densité à travers cette surface.

Bibliographie

  • C. Denis (1985). The Hydrostatic Figure of the Earth, Geophysical Report 85/02, Institut d'Astrophysique, Université de Liège.
  • C. Denis, E. Majewski, R. Teisseyre & J.B. Zieliński (1989). The Earth's Gravity Field, Chapitre 1, pages 1–77 de l'ouvrage collectif Physics and Evolution of the Earth's Interior, volume 4 : Gravity and Low-Frequency Geodynamics, édité par R. Teisseyre. Elsevier Publications, Amsterdam & PWN-Polish Scientific Publishers, Warszawa.
  • N.M. Günter (1957). Die Potentialtheorie, und ihre Anwendung auf Grundaufgaben der mathematischen Physik, B.G. Teubner, Leipzig.
  • O. Hölder (1882). Beiträge zur Potentialtheorie, Dissertation, Tübingen.
  • O.D. Kellogg (1929). Foundations of Potential Theory, Springer-Verlag, Berlin (Dover reprint, 1953).
  • W.D. MacMillan (1930). The theory of the potential, University of Chicago, Chicago (Dover reprint, 1958).

Voir aussi

Notes

  1. En effet
    V(P) = 2\pi G \rho_o \int_{ 0 \to r_b} \int_{-1 \to +1 } \sqrt{ x^2 + y^2 - 2 \, xyz} \, \mathrm dy \, \mathrm dz,

    qui devient ainsi successivement

    \begin{align}
V(P) & = & 2 \pi G \rho_o \int_{0 \to r_b } \dfrac{y [ \, |x-y| - (x+y) \, ]}{x} \mathrm dy \\
& = & -2 \pi G \rho_o \left[ \int_{0 \to x} \dfrac{2 y^2}{x} \, \mathrm dy + \int_{x \to r_b} 2y \, \mathrm dy \right] \\
& = & -2\pi G \rho_o \left[ \dfrac{2}{3} x^2 + ({r_b}^2 -x^2) \right] \\
& = & 2\pi G \rho_o \left( \dfrac{1}{3} x^2 - {r_b}^2 \right) 
\end{align}.

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Propriétés du potentiel newtonien de Wikipédia en français (auteurs)

Игры ⚽ Поможем написать реферат

Regardez d'autres dictionnaires:

  • Proprietes du potentiel newtonien — Propriétés du potentiel newtonien Statue d Isaac Newton à Trinity College, Cambridge. Cet article fort technique traitant en détail des propriétés d un potentiel inversement proportionnel à la distance s adresse avant tout aux étudiants et aux… …   Wikipédia en Français

  • Propriétés du potentiel newtonien (géodésie physique) — Propriétés du potentiel newtonien Statue d Isaac Newton à Trinity College, Cambridge. Cet article fort technique traitant en détail des propriétés d un potentiel inversement proportionnel à la distance s adresse avant tout aux étudiants et aux… …   Wikipédia en Français

  • Géodésie physique (propriétés du potentiel newtonien) — Propriétés du potentiel newtonien Statue d Isaac Newton à Trinity College, Cambridge. Cet article fort technique traitant en détail des propriétés d un potentiel inversement proportionnel à la distance s adresse avant tout aux étudiants et aux… …   Wikipédia en Français

  • Potentiel newtonien — On appelle potentiel newtonien tout potentiel scalaire « en  ». Cela vient du fait que la force newtonienne de gravitation soit « en  », en effet on a : avec K = GM1M2. D où  …   Wikipédia en Français

  • Potentiel de simple couche — Le plus souvent on conçoit la force gravifique par unité de masse, appelée gravité, comme engendrée par une distribution volumique de masse donnant lieu à un potentiel newtonien où r représente la distance entre un point potentiant Q et un point… …   Wikipédia en Français

  • Potentiel de double couche — Imaginons en un point potentiant quelconque Q un dipôle gravifique formé de deux masses ponctuelles de signes opposés, M + = + MQ et M − = − MQ, situées de part et d autre de Q à une distance h / 2 …   Wikipédia en Français

  • POTENTIEL ET FONCTIONS HARMONIQUES — La théorie du potentiel, directement issue de l’électrostatique, est une source d’inspiration extrêmement riche en analyse. Si, au début du XIXe siècle, on connaissait déjà l’équation de Laplace, la fonction de Green et l’intégrale de Poisson… …   Encyclopédie Universelle

  • Théorie du potentiel gravifique — Propriétés du potentiel newtonien Statue d Isaac Newton à Trinity College, Cambridge. Cet article fort technique traitant en détail des propriétés d un potentiel inversement proportionnel à la distance s adresse avant tout aux étudiants et aux… …   Wikipédia en Français

  • Géodésie physique (potentiel de simple couche) — Potentiel de simple couche Le plus souvent on conçoit la force gravifique par unité de masse, appelée gravité, comme engendrée par une distribution volumique de masse donnant lieu à un potentiel newtonien où r représente la distance entre un… …   Wikipédia en Français

  • Géodésie physique (potentiel de double couche) — Potentiel de double couche Imaginons en un point potentiant quelconque Q un dipôle gravifique formé de deux masses ponctuelles de signes opposés, M + = + MQ et M − = − MQ, situées de part et d autre de Q à une distance h / 2 …   Wikipédia en Français

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”