Sillage de Kelvin

Exemple de sillage de Kelvin derrière un bateau.

Le passage d'un bateau ou d'un oiseau aquatique à la surface d'une eau calme produit des ondes de gravité qui forment un motif caractéristique appelé sillage de Kelvin, où deux types d'ondes interviennent, les ondes divergentes, surtout visibles sur les bords du cône de sillage, et les ondes transverses (voir plus bas). C'est Lord Kelvin qui le premier a formalisé et résolu mathématiquement le problème.

Notations et équations utiles

On considère un bateau se déplaçant à la surface d'une eau initialement immobile dans le référentiel terrestre. Le bateau a un vecteur vitesse U dans ce référentiel. Soit un système de coordonnées fixe dans le référentiel terrestre. Soit un système de coordonnées fixe dans le référentiel du bateau et tel que z corresponde à la verticale et que U soit porté par la droite des x. Le passage d'un repère à l'autre se fait donc par la transformation galiléenne suivante :

On considère une onde plane progressive harmonique dont la fonction d'amplitude s'écrit, en notation complexe :

A₀ est l'amplitude maximale (éventuellement complexe), k est le vecteur d'onde (dont la norme est la pulsation spatiale) et ω est la pulsation temporelle de l'onde. En remplaçant r' par sa valeur en fonction de r, et en regroupant les termes, on obtient :

Ainsi l'onde vue depuis le bateau est perçue avec une fréquence temporelle ω-kU différente de celle que perçoit un observateur immobile dans le référentiel terrestre, depuis la plage par exemple (ω). C'est une illustration de l'effet Doppler.

Il existe une relation entre les pulsations spatiale et temporelle, liée intrinsèquement à la nature physique du milieu de propagation, et qui s'appelle la relation de dispersion. Si l'on considère que l'océan est infiniment profond, celle-ci s'écrit :

où g est l'accélération de la pesanteur. La vitesse de phase c des ondes de gravité, dans ce milieu, est donc :

Hypothèse des phases stationnaires

Depuis le bateau, les ondes du sillage semblent immobiles. Donc une condition nécessaire pour que l'onde précédemment décrite contribue au sillage est que l'amplitude de l'onde à r constant soit indépendante du temps, i.e. :

Sur la figure 1, le bateau se trouve au point O et se dirige vers les x négatifs. Le phénomène étant symétrique par rapport à l'axe des x, on restreint le raisonnement aux y positifs. L'instant où le bateau est en O est pris comme origine des temps. Les x positifs correspondent donc aux positions précédemment occupées par le bateau. Soit I le point où il se trouvait au temps -t. On a donc OI=Ut. Le cercle de rayon ct représente la position du front d'onde (à une pulsation ω donnée) né du passage du bateau en I.

Figure 1.

En quel(s) point(s) la perturbation créée en I contribue-t-elle à former le sillage à l'instant zéro ? On s'est donné une pulsation temporelle ω, donnons-nous un vecteur d'onde k (et le vecteur d'onde unitaire correspondant, ). Soit P l'intersection du demi-cercle de centre I et de rayon ct avec la droite passant par I et dirigée par k. C'est la position au temps zéro de la composante de pulsation ω et de vecteur d'onde k de la perturbation créée en I au temps -t.

Supposons que la condition ω=kU soit vérifiée, i.e. en ce point la perturbation crée en I contribue au sillage à l'instant zéro. On a alors :

Erreur math (La conversion en PNG a échoué ; vérifiez l’installation de latex et dvipng (ou dvips + gs + convert)): \|IP\|= c t = \frac \omega k t = \frac \bold k k t \bold U = \bold{k_u} \cdot \overrightarrow{IO}

.

Autrement dit, IP est la projection orthogonale de IO sur la droite passant par I et dirigée selon k. On en déduit que OP est perpendiculaire à IP. Donc P appartient au cercle de diamètre OI, dont le centre, noté O', est à mi-distance entre O et I.

Tout laissant k fixe, faisons varier ω, ce qui revient à faire varier c=g/ω, et donc le rayon ct du cercle en pointillés longs. On voit ainsi que, pour chaque point P du demi-cercle de diamètre OI, on peut trouver une valeur de ω telle que P soit à l'intersection du demi-cercle de rayon ct et du le demi-cercle de diamètre OI. Donc, au temps t₀, il passe une onde stationnaire par rapport au bateau en chaque point du demi-cercle de diamètre OI.

Cependant le raisonnement doit être adapté au fait que la vitesse de phase c des ondes de gravité, i.e. la vitesse à laquelle avance la crête des ondes, est différente de leur vitesse de groupe c_g, i.e. la vitesse à laquelle l'énergie de ces ondes est transportée.

La vitesse c_g d'un groupe ondes de gravité en eau profonde, vaut exactement la moitié de leur vitesse de phase c_p individuelle. En effet, la relation de dispersion, ω² = gk, s'écrit 2ωdω = gdk, d'où on tire que

C'est cette vitesse de groupe qu'il faut prendre en compte pour la caractérisation du sillage, et non pas la vitesse de phase. Pour un ω, un k et un t donnés, on a vu qu'on déduisait la position du point P où, à l'instant zéro, l'onde possède la même phase qu'en I à l'instant -t où elle est née. Pour cette onde, le pic d'énergie n'est pas en P mais a voyagé deux fois plus lentement, à la vitesse de groupe. Donc ce pic, à l'instant zéro, se situe au point P', à mi-chemin entre I et P. Par une homothétie de centre I et de rapport 1/2, on en déduit que P' se trouve sur le demi-cercle de diamètre O'I (O' milieu de OI). C'est donc en P' que l'onde créée en I à l'instant -t participe au sillage à l'instant zéro.

Cône de sillage

Inversement, en tout point appartenant à un demi-cercle de diamètre O''J, avec J quelconque sur la demi-droite Ox et O'' le milieu de OJ, une onde participera au sillage à l'instant zéro. Le cône de sillage est l'ensemble de ces points. Il est délimité par deux demi-droites symétriques par rapport à l'axe des x et d'origine O. Soit β₀ l'angle formé par ces demi-droites avec Ox :

ainsi, on en déduit que l’angle β₀ vaut environ 19,5° et l'angle du V formé par le sillage vaut environ 39°.

Ondes transverses et ondes divergentes

Considérons maintenant le problème opposé. Soit P un point de coordonnées (x, y). Où se trouve(nt) le(s) point(s) de la droite Ox d'où est partie, à un instant précédent, une onde qui participe au temps zéro au sillage au point P ? D'après la figure 1, un tel point M vérifie MO=3MP. On a donc :

MO² − 9MP² = 0

(MO − 3MP)(MO + 3MP) = 0

Soit G le barycentre du système {(O, 1), (P, -3)} et G' le barycentre du système {(O, 1), (P, 3)}. On a donc MO-3MP=-2MG et MO+3MP=4MG'. Donc M vérifie MG.MG'=0, i.e. M appartient au cercle de diamètre GG'. Ce cercle a pour centre le point de coordonnées et pour rayon . Un point de coordonnées (u, v) appartient au cercle si et seulement si ses coordonnées vérifient l'équation du cercle :

Le point M de coordonnées (u, v), intersection du cercle avec la droite des x, vérifie donc :

Ce qui donne :

On voit donc que :

• Il n'y a pas de solution pour u si , ce qui correspond à un point P situé hors du cône de sillage (on peut vérifier que ) ;
• il y a une unique solution si  ;
• et il y a deux solutions distinctes si .

Ce résultat signifie qu'aux points situés strictement à l'intérieur du cône de sillage, deux (et seulement deux) ondes de direction de propagation différente participant toutes deux au sillage, nées à des points différents de l'axe Ox, M _{− 1} et M₁, et à des temps différents, se croisent. Les abscisses x _{− 1} et x₁ de M _{− 1} et M₁ sont les deux solutions de l'équation ci-dessus :

Ces deux types d'ondes sont appelées divergentes et transverses. Elles seront caractérisées plus loin. Soit à présent θ l'angle entre MP et MO (ε = +ou- 1). C'est l'angle indiquant la direction de propagation de l'onde. On a :

Il est possible d'inverser cette relation, ce qui donne (le ε disparaît au cours du calcul) :

Articles connexes

• Sillage

• Portail de la physique