Équation de Poisson

En analyse vectorielle, l'équation de Poisson (ainsi nommée en l'honneur du mathématicien et physicien français Siméon Denis Poisson) est l'équation aux dérivées partielles elliptique du second ordre suivante :

${\displaystyle \displaystyle \Delta \phi =f}$

où ${\displaystyle \displaystyle \Delta }$ est l'opérateur laplacien et ${\displaystyle \displaystyle f}$ est une distribution généralement donnée.

Sur un domaine borné de ${\displaystyle \mathbb {R} ^{N}}$ et de frontière régulière, le problème de trouver ${\displaystyle \displaystyle \phi }$ à partir de ${\displaystyle \displaystyle f}$ et satisfaisant certaines conditions aux limites appropriées est un problème bien posé : la solution existe et est unique.

Ce problème est important en pratique :

• En électrostatique, la formulation classique (voir Équation de Poisson-Boltzmann) exprime le potentiel électrique ${\displaystyle \displaystyle V}$ associé à une distribution connue de charges ${\displaystyle \displaystyle \rho }$ par la relation

${\displaystyle \Delta V=-{\rho \over \varepsilon _{0}}.}$

• En gravitation universelle, le potentiel gravitationnel ${\displaystyle \displaystyle \Phi }$ est relié à la masse volumique ${\displaystyle \displaystyle \rho }$ par la relation

${\displaystyle \displaystyle \Delta \Phi =4\pi \,G\,\rho }$

• En mécanique des fluides, pour des écoulements incompressibles, la pression ${\displaystyle p}$ est reliée au champ de vitesse ${\displaystyle {\boldsymbol {u}}}$ par une équation de Poisson. Par exemple, en 2D, en notant les composantes du champ de vitesse ${\displaystyle {\boldsymbol {u}}=(u_{x},u_{y})}$, la relation s'écrit :

${\displaystyle \Delta p=-\rho \left(\left({\frac {\partial u_{x}}{\partial x}}\right)^{2}+2{\frac {\partial u_{x}}{\partial y}}{\frac {\partial u_{y}}{\partial x}}+\left({\frac {\partial u_{y}}{\partial y}}\right)^{2}\right),}$

où ${\displaystyle \rho }$ représente la masse volumique du fluide.

Conditions aux limites

L'équation de Poisson étant insensible à l’ajout sur ${\displaystyle \displaystyle \phi }$ d’une fonction satisfaisant l’équation de Laplace (ou une simple fonction linéaire par exemple), une condition aux limites est nécessaire pour espérer l'unicité de la solution : par exemple les conditions de Dirichlet, celles de Neumann, ou des conditions mixtes sur des portions de frontière.

Équation de Poisson à deux dimensions

En coordonnées cartésiennes dans ${\displaystyle \mathbb {R} ^{2}}$, considérons un ouvert ${\displaystyle \displaystyle \Omega }$, une fonction ${\displaystyle \displaystyle f}$ continue sur ${\displaystyle \displaystyle \Omega }$ et une fonction ${\displaystyle \displaystyle g}$ continue sur la frontière ${\displaystyle \partial \Omega }$. Le problème consiste à trouver une fonction de deux variables réelles ${\displaystyle \varphi (x,y)}$ définie sur ${\displaystyle \displaystyle \Omega }$ qui vérifie les deux relations :

${\displaystyle {\partial ^{2} \over \partial x^{2}}\varphi (x,y)+{\partial ^{2} \over \partial y^{2}}\varphi (x,y)=f(x,y)}$ sur ${\displaystyle \displaystyle \Omega }$ et ${\displaystyle \varphi =g}$ sur ${\displaystyle \partial \Omega .}$

Cette formulation est un modèle mathématique du problème statique d’une membrane élastique tendue et chargée (une peau de tambour) :

• ${\displaystyle \displaystyle f}$ est la densité de charge (exprimée par exemple en Pa, ceci à un multiple près caractérisant les propriétés d’élasticité de la membrane) ;
• ${\displaystyle \displaystyle g}$ est la cote (élévation verticale) le long de la frontière de fixation de la membrane ;
• la solution ${\displaystyle \varphi (x,y)}$ indique la cote de la membrane dans ${\displaystyle \displaystyle \Omega }$.

Éléments de justification

À une dimension, il s’agit d’une corde élastique chargée qui est fixée en ses deux extrémités.

Sur un petit élément ${\displaystyle \displaystyle [x-\delta x,x+\delta x]}$, considérons l’équilibre statique entre les deux forces de traction ${\displaystyle {\vec {F}}_{1}}$ et ${\displaystyle {\vec {F}}_{2}}$ de la corde (respectivement à gauche et à droite), puis la force de la charge ${\displaystyle {\vec {F}}_{G}}$ induite par une densité de charge linéaire notée ${\displaystyle \displaystyle \rho (x)}$ :

• ${\displaystyle {\vec {F}}_{1}=-(F_{1}/\delta x){\begin{pmatrix}\delta x\\\varphi (x)-\varphi (x-\delta x)\end{pmatrix}},}$
• ${\displaystyle {\vec {F}}_{2}=(F_{2}/\delta x){\begin{pmatrix}\delta x\\\varphi (x+\delta x)-\varphi (x)\end{pmatrix}},}$
• ${\displaystyle {\vec {F}}_{G}={\begin{pmatrix}0\\-2\rho (x)\delta x\end{pmatrix}}.}$

Sans restreindre la généralité, les facteurs ${\displaystyle \displaystyle F_{1}}$ et ${\displaystyle \displaystyle F_{2}}$ ont été divisés par ${\displaystyle \displaystyle \delta x}$ afin de leur conserver une grandeur non différentielle.

La somme vectorielle de ces forces conduit aux égalités :

• ${\displaystyle \displaystyle F_{1}=F_{2}}$ qu’on peut appeler ${\displaystyle 2\displaystyle k}$, un coefficient indépendant de ${\displaystyle \displaystyle x}$ puisque toutes les composantes horizontales se compensent pour se répercuter uniquement sur les points d’attache,
• ${\displaystyle {2k \over \delta x}[\varphi (x+\delta x)-2\varphi (x)+\varphi (x-\delta x)]=2\rho (x)\delta x}$ qui, lorsque ${\displaystyle \delta x}$ tend vers 0, s’écrit ${\displaystyle k\,\varphi ''(x)=\rho (x).}$

Cette dernière relation est bien l’équation de Poisson à une dimension.

Formulation faible et solution

Soit ${\displaystyle \displaystyle \Omega }$ un domaine ouvert et borné de ${\displaystyle \mathbb {R} ^{N}}$ dont la frontière ${\displaystyle \partial \Omega }$ est suffisamment régulière pour satisfaire le théorème de la divergence. Soit ${\displaystyle \mathbf {n} }$ le vecteur normal à ${\displaystyle \partial \Omega }$ et dirigé vers l’extérieur.

Soient ${\displaystyle \displaystyle f}$ une fonction de ${\displaystyle \displaystyle L^{2}(\Omega )}$, puis ${\displaystyle \displaystyle g}$ et ${\displaystyle \alpha >0}$ des fonctions continues définies sur ${\displaystyle \partial \Omega }$.

On cherche une solution ${\displaystyle \displaystyle \phi }$ pour chacun des problèmes suivants :

${\displaystyle \displaystyle -\Delta \phi =f}$ sur ${\displaystyle \displaystyle \Omega }$

satisfaisant l’une des conditions sur ${\displaystyle \partial \Omega }$ :

1. ${\displaystyle \displaystyle \phi =0}$
2. ${\displaystyle \nabla \phi \cdot \mathbf {n} =g}$ et ${\displaystyle \int _{\Omega }\phi \,\mathrm {d} V=0}$ (pour fixer la constante additive d’indétermination)
3. ${\displaystyle \nabla \phi \cdot \mathbf {n} +\alpha \phi =0}$

Pour toute fonction ${\displaystyle \displaystyle \psi }$ régulière, la relation

${\displaystyle {\mathrm {div} }(\psi \,\nabla \phi )=\nabla \phi \cdot \nabla \psi +\psi \Delta \phi }$

et le théorème de la divergence impliquent

${\displaystyle \int _{\Omega }\nabla \phi \cdot \nabla \psi \,\mathrm {d} V=-\int _{\Omega }\psi \,\Delta \phi \,\mathrm {d} V+\int _{\partial \Omega }\psi \,\nabla \phi \cdot \mathbf {n} \,\mathrm {d} S.}$

Si ${\displaystyle \displaystyle \phi }$ est solution du problème précédent muni de la condition aux limites retenue, alors

1. ${\displaystyle \int _{\Omega }\nabla \phi \cdot \nabla \psi \,\mathrm {d} V=\int _{\Omega }f\,\psi \,\mathrm {d} V}$
2. ${\displaystyle \int _{\Omega }\nabla \phi \cdot \nabla \psi \,\mathrm {d} V=\int _{\Omega }f\,\psi \,\mathrm {d} V+\int _{\partial \Omega }g\,\psi \,\mathrm {d} S}$
3. ${\displaystyle \int _{\Omega }\nabla \phi \cdot \nabla \psi \,\mathrm {d} V+\int _{\partial \Omega }\alpha \,\phi \,\psi \,\mathrm {d} S=\int _{\Omega }f\,\psi \,\mathrm {d} V}$

En notant ${\displaystyle a(\phi ,\,\psi )}$ le membre de gauche et ${\displaystyle \displaystyle b(\psi )}$ celui de droite, la formulation faible consiste à :

• définir un espace vectoriel ${\displaystyle \displaystyle H}$ approprié dans lequel ${\displaystyle \displaystyle a(.,.)}$ et ${\displaystyle \displaystyle b(.)}$ sont définies,
• rechercher ${\displaystyle \displaystyle \phi \in H}$ tel que ${\displaystyle a(\phi ,\,\psi )=b(\psi )}$ pour tout ${\displaystyle \psi \in H}$.

Si elle existe, la solution naturelle de ces formulations se trouve dans l’espace de Sobolev ${\displaystyle \displaystyle H^{1}(\Omega )}$ muni de sa norme ${\displaystyle \|\psi \|_{H^{1}}^{2}=\|\psi \|_{L^{2}}^{2}+\|\nabla \psi \|_{L^{2}}^{2}.}$

En effet, pour chaque problème, ${\displaystyle a(.,.)}$ est une forme bilinéaire symétrique définie sur ${\displaystyle H^{1}(\Omega )\times H^{1}(\Omega )}$, et ${\displaystyle b(.)}$ est une forme linéaire sur ${\displaystyle \displaystyle H^{1}(\Omega )}$.

Proposition — Soit ${\displaystyle \displaystyle \Omega }$ un domaine ouvert et borné de ${\displaystyle \mathbb {R} ^{N}}$ et de frontière ${\displaystyle \partial \Omega }$ régulière (ou régulière par morceaux), ${\displaystyle \displaystyle f}$ dans ${\displaystyle \displaystyle L^{2}(\Omega )}$, puis ${\displaystyle \displaystyle g}$ et ${\displaystyle \alpha >0}$ des fonctions continues définies sur ${\displaystyle \partial \Omega }$.

Alors les trois problèmes précédents possèdent une unique solution ${\displaystyle \displaystyle \phi }$ dans ${\displaystyle \displaystyle H^{1}(\Omega )}$ qui est caractérisée par la formulation faible correspondante mise en œuvre dans les espaces suivants :

1. ${\displaystyle H=H_{0}^{1}(\Omega )}$ qui est l’adhérence dans ${\displaystyle \displaystyle H^{1}(\Omega )}$ des fonctions indéfiniment dérivables et à support compact dans ${\displaystyle \displaystyle \Omega .}$
2. ${\displaystyle H=\left\{\phi \in H^{1}(\Omega )\,|\int _{\Omega }\phi \,\mathrm {d} V=0\right\}.}$
3. ${\displaystyle \displaystyle H=H^{1}(\Omega ).}$

Justification

Si les conditions de continuité et de coercivité des hypothèses du théorème de Lax-Milgram sont satisfaites, ce dernier permet de conclure.

• Continuité

Pour la continuité des deux formes, il s’agit de montrer l’existence de constantes positives notées génériquement ${\displaystyle \displaystyle c}$ telles que

${\displaystyle |a(\phi ,\,\psi )|\leqslant c\,\|\phi \|_{H^{1}}\|\psi \|_{H^{1}},}$

${\displaystyle |b(\,\psi )|\leqslant c\,\|\psi \|_{H^{1}}.}$

Ces constantes existent par définition de la norme de ${\displaystyle \displaystyle H^{1}(\Omega )}$ et par la continuité des opérateurs de trace qui, à une fonction ${\displaystyle \psi \in H^{1}(\Omega )}$ associe une fonction de ${\displaystyle L^{2}(\partial \Omega )}$ définie par la restriction de ${\displaystyle \displaystyle \psi }$ sur ${\displaystyle \partial \Omega }$.

On peut remarquer que la continuité des formes assure simultanément leur définition rigoureuse. Pour le second problème en particulier, ${\displaystyle \displaystyle \Omega }$ borné implique la continuité de l’injection de ${\displaystyle \displaystyle L^{2}(\Omega )}$ dans ${\displaystyle \displaystyle L^{1}(\Omega )}$ pour la norme ${\displaystyle \|.\|_{L^{1}}}$, ce qui justifie la définition de l’espace ${\displaystyle \displaystyle H}$ correspondant.

• Coercivité

Pour la coercivité des ${\displaystyle \displaystyle a(.,.)}$, il s’agit de montrer l’existence d’une constante ${\displaystyle \displaystyle \mu >0}$ indépendante de ${\displaystyle \psi \in H}$ telle que

• ${\displaystyle |a(\psi ,\,\psi )|\geqslant \mu \,\|\phi \|_{H^{1}}^{2}.}$

Cette propriété découle de l’inégalité de Poincaré classique pour la forme ${\displaystyle \displaystyle a_{1}(.,.)}$ et de l’inégalité de Poincaré-Wirtinger pour la forme ${\displaystyle \displaystyle a_{2}(.,.)}$.

La coercivité de la forme ${\displaystyle \displaystyle a_{3}(.,.)}$ peut se montrer par l’absurde. En notant

${\displaystyle d(\psi )=\int _{\partial \Omega }\alpha \,\psi ^{2}\,\mathrm {d} S,}$

supposons qu’il existe une suite ${\displaystyle \psi _{n}\in H^{1}(\Omega )}$ satisfaisant

${\displaystyle \|\psi _{n}\|_{H^{1}}=1}$ et ${\displaystyle a_{3}(\psi _{n},\,\psi _{n})=d(\psi _{n})+\|\nabla \psi _{n}\|_{L^{2}}^{2}}$ tend vers 0.

Par compacité de l’injection canonique de ${\displaystyle \displaystyle H^{1}(\Omega )}$ dans ${\displaystyle \displaystyle L^{2}(\Omega )}$ (lorsque ${\displaystyle \displaystyle \Omega }$ est borné), il existe une sous-suite convergeant vers une fonction ${\displaystyle \displaystyle \psi }$ pour la norme ${\displaystyle \|.\|_{L^{2}}}$. Cette suite est donc une suite de Cauchy dans ${\displaystyle \displaystyle L^{2}(\Omega )}$ et, puisque son gradient tend vers 0 dans ${\displaystyle \displaystyle L^{2}(\Omega )}$, elle est également une suite de Cauchy dans ${\displaystyle \displaystyle H^{1}(\Omega )}$ qui converge vers ${\displaystyle \psi \in H^{1}(\Omega )}$ et qui ne peut être qu’une fonction constante avec ${\displaystyle \|\psi \|_{H^{1}}=1}$. Ainsi, sa trace sur ${\displaystyle \partial \Omega }$ (par continuité) ne peut être que constante non nulle, ce qui contredit ${\displaystyle d(\psi )=0}$.

Résolution

Il y a diverses méthodes pour la résolution numérique. La méthode de relaxation, un algorithme itératif, est un exemple. Les méthodes basées sur les transformées de Fourier sont presque toujours utilisées en gravitation universelle.

Considérations historiques et essais de résolution

L'équation de Poisson est une correction célèbre de l’équation différentielle de Laplace au second degré pour le potentiel :

${\displaystyle \nabla ^{2}\phi =-4\pi \rho \;,}$

On appelle aussi cette équation : l'équation de la théorie du potentiel publiée en 1813. Si une fonction d’un point donné ρ = 0, nous obtenons l’équation de Laplace :

${\displaystyle \nabla ^{2}\phi =0\;.}$

En 1812, Poisson découvrit que cette équation n’est valide que hors d’un solide. Une preuve rigoureuse pour les masses avec une densité variable fut d’abord donnée par Carl Friedrich Gauss en 1839. Les deux équations ont leurs équivalents en analyse vectorielle. L’étude des champs scalaires φ d’une divergence^[Quoi ?] donne :

${\displaystyle \nabla ^{2}\phi =\rho (x,y,z)\;.}$

Par exemple, une équation de Poisson pour un potentiel électrique en surface Ψ, qui montre sa dépendance de la densité d’une charge électrique ρ_e dans une place particulière :

${\displaystyle \nabla ^{2}\Psi ={\partial ^{2}\Psi \over \partial x^{2}}+{\partial ^{2}\Psi \over \partial y^{2}}+{\partial ^{2}\Psi \over \partial z^{2}}=-{\rho _{e} \over \varepsilon \varepsilon _{0}}\;.}$

La distribution d’une charge dans un fluide est inconnue et nous devons utiliser l’équation de Poisson-Boltzmann :

${\displaystyle \nabla ^{2}\Psi ={n_{0}e \over \varepsilon \varepsilon _{0}}\left(e^{e\Psi (x,y,z) \over k_{B}T}-e^{-e\Psi (x,y,z) \over k_{B}T}\right)\;,}$

ce qui, dans la plupart des cas, ne peut être résolu analytiquement, mais seulement pour des situations particulières. En coordonnées polaires, l’équation de Poisson-Boltzmann s'écrit :

${\displaystyle {1 \over r^{2}}{d \over dr}\left(r^{2}{d\Psi \over dr}\right)={n_{0}e \over \varepsilon \varepsilon _{0}}\left(e^{e\Psi (r) \over k_{B}T}-e^{-e\Psi (r) \over k_{B}T}\right)\;,}$

laquelle ne peut pas non plus être résolue analytiquement. Même si le champ φ n’est pas scalaire, l’équation de Poisson est valide, comme elle peut l’être par exemple dans un espace de Minkowski à quatre dimensions :

${\displaystyle \square \phi _{ik}=\rho (x,y,z,ct)\;.}$

Si ρ(x, y, z) est une fonction continue et si pour r→∞ (ou si un point 'se déplace' à l’infini) une fonction φ va à 0 suffisamment rapidement, une solution à l’équation de Poisson est le potentiel newtonien d’une fonction ρ(x, y, z) :

${\displaystyle \phi _{M}=-{1 \over 4\pi }\int {\rho (x,y,z)dv \over r}\;,}$

où r est une distance entre l’élément avec le volume dv et le point M. L’intégration parcourt la totalité de l’espace. L’intégrale de Poisson en résolvant la fonction de Green pour le Problème de Dirichlet de l’équation de Laplace, si le cercle est le domaine étudié :

${\displaystyle \phi (\xi ,\eta )={1 \over 2\pi }\int _{0}^{2\pi }{R^{2}-\rho ^{2} \over R^{2}+\rho ^{2}-2R\rho \cos(\psi -\chi )}\phi (\chi )d\chi \;,}$

où :

${\displaystyle \xi =\rho \cos \psi \;,\quad \eta =\rho \sin \psi \;.}$

φ(χ) est une fonction prescrite sur une ligne circulaire, qui définit les conditions aux limites de la fonction requise φ de l’équation de Laplace. De la même manière nous définissons la fonction de Green pour le problème de Dirichlet pour l’équation de Laplace ${\displaystyle {\vec {\nabla }}^{2}\phi =0}$ dans l’espace, pour un domaine constitué d’une sphère de rayon R. Cette fois la fonction de Green est:

${\displaystyle G(x,y,z;\xi ,\eta ,\zeta )={1 \over r}-{R \over r_{1}\rho }\;,}$

où : ${\displaystyle \rho ={\sqrt {\xi ^{2}+\eta ^{2}+\zeta ^{2}}}}$ est une distance d’un point (ξ, η, ζ) depuis le centre d’une sphère, r une distance entre des points (x, y, z), (ξ, η, ζ), r₁ est une distance entre le point (x, y, z) et le point (Rξ/ρ, Rη/ρ, Rζ/ρ), symétrique au point (ξ, η, ζ). L’intégrale de Poisson est maintenant de la forme:

${\displaystyle \phi (\xi ,\eta ,\zeta )={1 \over 4\pi }\int \!\!\!\int _{S}{R^{2}-\rho ^{2} \over Rr^{3}}\phi ds\;.}$

Notes et références

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Siméon Denis Poisson » (voir la liste des auteurs).

Voir aussi

Bibliographie

• [Poisson 1813] Siméon-Denis Poisson, « Remarques sur une équation qui se présente dans la théorie des attractions des sphéroïdes », Nouveau bulletin des sciences : par la Société philomat(h)ique (de Paris), Paris, J. Klostermann fils, t. III, n^o 75,‎ déc. 1813, p. 388-392 (lire en ligne).
• [Godard et Boer 2020] Roger Godard et John de Boer, « Gauss et le modèle du champ magnétique terrestre », dans Maria Zack et Dirk Schlimm (éd.), Research in history and philosophy of mathematics : the CSHPM 2018 volume [« Recherche en histoire et en philosophie des mathématiques »], Cham, Birkhäuser, coll. « Proceedings of the Canadian Society for History and Philosophy of Mathematics / (Actes de la) Société canadienne d'histoire et de philosophie des mathématiques », janv. 2020, 1^re éd., 1 vol., XIII-172, ill., 15,6 × 23,4 cm (ISBN 978-3-030-31196-4, OCLC 1154674154, DOI 10.1007/978-3-030-31298-5, présentation en ligne, lire en ligne), chap. 8, p. 125-138.
• [Solomentsev 1995] (en) E. D. Solomentsev, « Poisson equation », dans Michiel Hazewinkel (éd.), Encyclopaedia of mathematics : an updated and annotated translation of the Soviet Mathematical encyclopaedia [« Encyclopédie des mathématiques : une traduction mise à jour et annotée de l'Encyclopédie des mathématiques soviétique »], t. IV : Monge-Ampère equation – Rings and algebras [« Équation de Monge-Ampère – Anneaux et algèbre »], Dordrecht, Kluwer Academic, hors coll., janv. 1995, 1^re éd., 1 vol., IV-929, ill., 21 × 29,7 cm (ISBN 1-556-08010-7, EAN 9781556080104, OCLC 36917086, DOI 10.1007/978-1-4899-3791-9, SUDOC 030253195, présentation en ligne, lire en ligne), s.v. Poisson equation [« équation de Poisson »], p. 445.
• [Taillet, Villain et Febvre 2018] Richard Taillet, Loïc Villain et Pascal Febvre, Dictionnaire de physique, Louvain-la-Neuve, De Boeck Supérieur, hors coll. / physique, janv. 2018, 4^e éd. (1^re éd. mai 2008), 1 vol., X-956, ill. et fig., 17 × 24 cm (ISBN 978-2-8073-0744-5, EAN 9782807307445, OCLC 1022951339, BNF 45646901, SUDOC 224228161, présentation en ligne, lire en ligne), s.v. Poisson (équation de), p. 579-580.
• Poisson Equation at EqWorld: The World of Mathematical Equations.
• L.C. Evans, Partial Differential Equations, American Mathematical Society, Providence, 1998. (ISBN 0-8218-0772-2)
• A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002. (ISBN 1-58488-299-9)

Articles connexes

• Opérateur laplacien
• Opérateur laplacien vectoriel
• Équation de Laplace
• Théorie du potentiel

Liens externes

• (en) Eric W. Weisstein, « Poisson Integral », sur MathWorld
• (en) Eric W. Weisstein, « Poisson Kernel », sur MathWorld
• (en) Poisson's Integral for the Unit Disk
• Notices dans des dictionnaires ou encyclopédies généralistes :

  • Britannica
  • Den Store Danske Encyklopædi
  • Nationalencyklopedin

• Portail de l'analyse
• Portail de la physique