Équation des ondes — Wikipédia

Cet article ne cite pas suffisamment ses sources (juin 2017).

Si vous disposez d'ouvrages ou d'articles de référence ou si vous connaissez des sites web de qualité traitant du thème abordé ici, merci de compléter l'article en donnant les références utiles à sa vérifiabilité et en les liant à la section « Notes et références ».

En pratique : Quelles sources sont attendues ? Comment ajouter mes sources ?

Pour les articles homonymes, voir ondes (homonymie).

L'équation de d'Alembert ou équation des ondes est une équation aux dérivées partielles en physique qui régit la propagation d'une onde^{[N 1]}. C'est une équation vérifiée par de nombreux phénomènes ondulatoires de la vie courante comme le son ou la lumière.

Énoncé[modifier | modifier le code]

L'équation des ondes se formule^{[N 2]} :

${\displaystyle \Delta {\vec {\mathrm {E} }}={\frac {1}{c^{2}}}{\frac {\partial ^{2}{\vec {\mathrm {E} }}}{\partial t^{2}}}}$

avec :

• ${\displaystyle \Delta }$ l'opérateur laplacien ;
• ${\displaystyle {\vec {\mathrm {E} }}}$ l'onde vectorielle^{[N 3]};
• ${\displaystyle c}$ une constante, vitesse de propagation de ${\displaystyle \mathrm {\vec {E}} }$ dans le milieu considéré ;

L'utilisation du laplacien permet de s'affranchir du choix d'un système de coordonnées.

En coordonnées cartésiennes, l'équation des ondes devient :

${\displaystyle {\frac {\partial ^{2}{\vec {\mathrm {E} }}}{\partial x^{2}}}+{\frac {\partial ^{2}{\vec {\mathrm {E} }}}{\partial y^{2}}}+{\frac {\partial ^{2}{\vec {\mathrm {E} }}}{\partial z^{2}}}={\frac {1}{c^{2}}}{\frac {\partial ^{2}{\vec {\mathrm {E} }}}{\partial t^{2}}}}$

avec :

• ${\displaystyle {\frac {\partial ^{2}{\vec {\mathrm {E} }}}{\partial u^{2}}}}$ l'opérateur de dérivée partielle seconde en ${\displaystyle u}$ appliqué sur ${\displaystyle \mathrm {\vec {E}} }$ ;
• ${\displaystyle x}$, ${\displaystyle y}$, ${\displaystyle z}$ les trois variables cartésiennes de l'espace, et ${\displaystyle t}$ celle du temps.

L'équation des ondes s'applique à des fonctions scalaires ou vectorielles, qu'on formalise en champ vectoriel ou champ scalaire. Le champ ${\displaystyle {\vec {\mathrm {E} }}}$ renseigne à la fois sur l'amplitude de l'onde et sa polarisation. Une équation des ondes vectorielle regroupe trois équations des ondes scalaires.

Histoire[modifier | modifier le code]

L'établissement de l'équation des ondes est venu de l’étude des vibrations d'une corde de violon. Afin de pouvoir modéliser ce comportement, les mathématiciens du XVII^e siècle ont appliqué la deuxième loi de Newton à la corde, d'abord vue comme un ensemble fini de masses ponctuelles reliées par des ressorts (dont le comportement est donné par la loi de Hooke établie en 1660), avant d'augmenter le nombre de masses pour se rapprocher de la corde^[2].

En 1727, Jean Bernoulli reprend l'expérience de la corde de violon et constate que ses vibrations forment une sinusoïde et que la variation de son amplitude en un point forme également une courbe sinusoïdale, mettant ainsi en évidence les modes^[2]. En 1746, Jean Le Rond d'Alembert reprend le modèle des masses ponctuelles liées par des ressorts et établit uniquement à partir des équations que les vibrations de la corde dépendent à la fois de l'espace et du temps.

Exemples en dimension 1[modifier | modifier le code]

Ressort[modifier | modifier le code]

Pour un système masses-ressorts passé à la limite continue, de constante de raideur totale ${\displaystyle k}$, de longueur totale ${\displaystyle L}$ et de masse totale ${\displaystyle m}$, la fonction de déplacement ${\displaystyle u(x,t)}$ vérifie :

${\displaystyle {\frac {\partial ^{2}u}{\partial x^{2}}}={\frac {1}{c^{2}}}{\frac {\partial ^{2}u}{\partial t^{2}}}}$ avec ${\displaystyle c={\sqrt {\frac {kL^{2}}{m}}}}$

Démonstration

Considérons une chaine de masses ponctuelles ${\displaystyle m_{dx}}$ reliées par des ressorts sans masse de longueur ${\displaystyle h}$ et de raideur ${\displaystyle k_{dx}}$:

Notons ${\displaystyle u(x)}$ le déplacement de la masse ${\displaystyle m_{dx}}$ en ${\displaystyle x}$ par rapport à sa position de repos. La résultante des forces de rappel exercées sur la masse ${\displaystyle m_{dx}}$ au point ${\displaystyle x+h}$ est :

${\displaystyle F_{\mathrm {rappel} }=F_{x+2h}-F_{x}=k_{dx}\left[{u(x+2h,t)-u(x+h,t)}\right]-k_{dx}[u(x+h,t)-u(x,t)]}$

Le théorème de la quantité de mouvement appliqué à la masse ${\displaystyle m_{dx}}$ établit ensuite :

${\displaystyle m_{dx}{\partial ^{2}u \over \partial t^{2}}(x+h,t)=k_{dx}[u(x+2h,t)-u(x+h,t)-u(x+h,t)+u(x,t)]}$

En découpant la chaine en ${\displaystyle N}$ ressorts équidistants, la longueur totale vérifie ${\displaystyle L=N\times h}$, la masse totale ${\displaystyle m=N\times m_{dx}}$, et la raideur totale ${\displaystyle {\frac {1}{k}}={\frac {N}{k_{dx}}}}$ ^{[N 4]}. On obtient alors :

${\displaystyle {\partial ^{2}u \over \partial t^{2}}(x+h,t)={kL^{2} \over m}{u(x+2h,t)-2u(x+h,t)+u(x,t) \over h^{2}}}$

La fonction ${\displaystyle u}$ étant supposée deux fois dérivable, en faisant tendre ${\displaystyle N}$ vers l'infini et donc ${\displaystyle h}$ vers ${\displaystyle 0}$, le taux d'accroissement tend vers la dérivée seconde escomptée.

Corde[modifier | modifier le code]

Pour une corde sans raideur de longueur ${\displaystyle L}$, de masse ${\displaystyle m}$ sous la tension ${\displaystyle T}$, et avec l'hypothèse de petites déformations, l'élongation ${\displaystyle Y}$ vérifie :

${\displaystyle {\frac {\partial ^{2}Y}{\partial x^{2}}}={\frac {1}{c^{2}}}{\frac {\partial ^{2}Y}{\partial t^{2}}}}$ avec ${\displaystyle c={\sqrt {\frac {TL}{m}}}}$

L'énergie ${\displaystyle E=TL}$ d'application de la tension ${\displaystyle T}$ sur la longueur ${\displaystyle L}$ vérifie ${\displaystyle E=mc^{2}}$.

Câble coaxial[modifier | modifier le code]

Pour un câble coaxial de capacité linéique ${\displaystyle \Gamma }$ et d'inductance linéique ${\displaystyle \Lambda }$, l'intensité ${\displaystyle I}$ et la tension ${\displaystyle U}$ vérifient toutes deux :

${\displaystyle {\frac {\partial ^{2}I}{\partial x^{2}}}={\frac {1}{c^{2}}}{\frac {\partial ^{2}I}{\partial t^{2}}}}$ et ${\displaystyle {\frac {\partial ^{2}U}{\partial x^{2}}}={\frac {1}{c^{2}}}{\frac {\partial ^{2}U}{\partial t^{2}}}}$ avec ${\displaystyle c={\sqrt {\frac {1}{\Lambda \Gamma }}}}$

Barreau élastique[modifier | modifier le code]

Pour un barreau élastique de module de Young ${\displaystyle E}$, de volume ${\displaystyle V}$ et de masse ${\displaystyle m}$, l'allongement ${\displaystyle \epsilon }$ vérifie :

${\displaystyle {\frac {\partial ^{2}\epsilon }{\partial x^{2}}}={\frac {1}{c^{2}}}{\frac {\partial ^{2}\epsilon }{\partial t^{2}}}}$ avec ${\displaystyle c={\sqrt {\frac {EV}{m}}}}$

Résolution[modifier | modifier le code]

En dimension 1[modifier | modifier le code]

En dimension 1 d'espace, l'équation des ondes se simplifie en^[3],[4],[5] :

${\displaystyle {\frac {\partial ^{2}U}{\partial x^{2}}}={\frac {1}{c^{2}}}{\frac {\partial ^{2}U}{\partial t^{2}}}}$

La solution générale de cette équation est alors la somme de deux fonctions indépendantes :

${\displaystyle U(x,t)=F(x-ct)+G(x+ct)}$

${\displaystyle F}$ est une onde nommée progressive, car elle se propage dans le sens des ${\displaystyle x}$ croissants, tandis que ${\displaystyle G}$ est nommée régressive car se propageant dans le sens des ${\displaystyle x}$ décroissants. Lorsqu'on suit des yeux la perturbation, par exemple le haut d'une sinusoïdale, on observe en fait un point de phase constante, c'est-à-dire au point tel que ${\displaystyle \phi =x-ct}$ soit constante, dans le cas de l'onde progressive. Comme le temps avance, ${\displaystyle t}$ croit, et ${\displaystyle x}$ doit alors croître à son tour pour maintenir ${\displaystyle \phi }$ constante. L'onde semble alors avancer dans le sens des ${\displaystyle x}$ croissants.

Il est possible de montrer que ${\displaystyle c}$ représente la vitesse de propagation de l'onde en cherchant une solution progressive de type ${\displaystyle y(x-vt)}$. On obtient alors ${\displaystyle v^{2}=c^{2}}$, le signe de ${\displaystyle v}$ dépendant du sens de propagation de l'onde.

En dimension 3[modifier | modifier le code]

Dans le cas d'une onde scalaire dans un milieu homogène, il convient de travailler en coordonnées sphériques pour résoudre l'équation des ondes :

${\displaystyle {\frac {1}{c^{2}}}{\frac {\partial ^{2}u}{\partial t^{2}}}={\frac {\partial ^{2}u}{\partial r^{2}}}+{\frac {2}{r}}{\frac {\partial u}{\partial r}}.}$

En réécrivant l'équation sous la forme :

${\displaystyle {\frac {1}{c^{2}}}{\frac {\partial ^{2}(ru)}{\partial t^{2}}}-{\frac {\partial ^{2}(ru)}{\partial r^{2}}}=0,}$

il vient, en reprenant les calculs faits sur le problème 1D, que la solution s'écrit sous la forme :

${\displaystyle u(r,t)={\frac {1}{r}}F(r-ct)+{\frac {1}{r}}G(r+ct),}$

où F et G sont des fonctions arbitraires.

Il apparaît ainsi que les solutions sont des ondes sphériques, se propageant ou se rapprochant du point d'origine du repère, considéré comme un point source, où les ondes sont singulières tandis qu'elles s'éloignent avec une amplitude décroissante en ¹⁄_r.

Conservation de l'énergie[modifier | modifier le code]

Si ${\displaystyle u}$ est une solution de l'équation des ondes alors l'énergie

${\displaystyle E(u(t))={\frac {1}{2}}\int _{\mathbb {R} ^{N}}\left|{\frac {\partial u}{\partial t}}(t,x)\right|^{2}\mathrm {d} x+{\frac {c^{2}}{2}}\int _{\mathbb {R} ^{N}}\left|\nabla u(t,x)\right|^{2}\mathrm {d} x}$

est conservée au cours du temps. Ici on a noté ${\displaystyle N}$ la dimension d'espace et

${\displaystyle \left|\nabla u(t,x)\right|^{2}=\sum _{j=1}^{N}\left|{\frac {\partial u}{\partial x_{j}}}(t,x)\right|^{2}.}$

Équation dans un domaine borné avec condition au bord[modifier | modifier le code]

On peut également considérer l'équation des ondes dans un domaine de l'espace ${\displaystyle D}$ :

${\displaystyle \square u(t,x)=0\quad t\in \mathbb {R} ,\quad x\in D}$

avec comme condition aux limites, par exemple :

${\displaystyle u(t,x)=0,\quad t\in \mathbb {R} ,\quad x\in \partial D}$

(condition aux limites de Dirichlet) où ${\displaystyle \partial D}$ est le bord du domaine ${\displaystyle D}$, ou

${\displaystyle \partial _{\nu }u(t,x)=0,\quad t\in \mathbb {R} ,\quad x\in \partial D}$

(condition aux limites de Neumann) où ${\displaystyle \partial _{\nu }}$ est la dérivée normale extérieure au bord ${\displaystyle \partial D}$.

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

Notes[modifier | modifier le code]

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

Voir aussi[modifier | modifier le code]

Onde sur une corde vibrante

• Portail de la physique
• Portail de l'analyse