ANSWER - Modèle de propagation des ondes longues

Eléments de contexte

Le Groupement d'Intérêt Scientifique "Hydraulique pour l'Ecologie et le Développement durable" (GIS HED2) initie une action collaborative de recensement de solutions analytiques dans le domaine de la mécanique des fluides à surface libre. L'objectif de ce recensement est multiple :

• il permet de disposer d'une base de référence de cas très simples, qui constituent des simplifications des équations de base et qui permettent de caler des codes de calcul
• il peut servir d'accompagnement de programmes de formation pour des techniciens, ingénieurs et chercheurs qui souhaitent disposer d'une information très accessible et concentrée en un site

Pour atteindre cet objectif, chaque fiche comprend un rappel des simplifications réalisées, l'expression des modèles simplifiés, la solution analytique correspondante, la résolution d'un problème simple, une proposition de montage en laboratoire et une proposition de réalisation d'images ou de films en nature pour illustrer ce processus. Pour cela, un ensemble de 2 fiches types sont réalisées : l'une théorique (fiche à venir) explicite la méthode d'obtention de la solution analytique, la seconde, plus pratique concerne le traitement de la solution analytique (qui constitue la présente page).

Hiérarchie des hypothèses simplificatrices

Navier-Stokes

→ fluide incompressible

→ intégration dans une section de calcul (canal rectangulaire) ==> Saint-Venant 1D

→ accélération négligeable

→ frottement négligeable

Expression de l'équation simplifiée

A partir des hypothèses précédentes, considérons un canal infini de forme rectangulaire : largeur $ b $ et profondeur d'eau $ H $.
Soit $ h $ le niveau d'eau, $ u $ la vitesse moyenne de l'écoulement dans la section de surface $ A = bH $ et enfin $ Q=bHu $ le débit, les équations simplifiées prennent la forme suivante :

$ \begin{cases} \dfrac{ \partial h }{ \partial t }+H \dfrac{ \partial u }{ \partial x }=0 } \\ \dfrac{ \partial u }{ \partial t }+g \dfrac{ \partial h }{ \partial x }=0 \end{cases} $ Si l'on dérive la première équation par rapport à $ t $ et la seconde par rapport à $ x $, et en éliminant le terme $ \frac{ \partial^2 u }{ \partial t \partial t } $

Nous obtenons l'équation des ondes suivante: $ \frac{ \partial^2 h }{ \partial t^2 }-gH \frac{ \partial^2 h }{ \partial x^2 } $

Expression de la solution analytique

Cette équation peut se mettre sous la forme: $ (\frac{ \partial h }{ \partial t }-c \frac{ \partial h }{ \partial x }) ( \frac{ \partial h }{ \partial t }+c \frac{ \partial h }{ \partial x })=0 $ Une solution de cette équation pour l'évolution de la surface libre est la suivante:
$ h=acos(kx-\sigma t) $ qui correspond à la propagation de 2 ondes de vitesse $ c= \sqrt{gH} $ dans les deux sens (vers la gauche et vers la droite)

Cas d'application : translation d'une onde sinusoïdale

Onde de marée progressive

Le cas que nous présentons ici est un cas d'école, car il met en scène la propagation d'une onde de marée de période de 12 heures dans un canal long de 10^7 m, ceci pour illustrer plusieurs longueurs d'onde sur une même figure.
Les caractéristiques

Onde de marée stationnaire

Domaine 1D

Conditions initiales

Conditions limites

Animation des résultats

Domaines d'application

houle linéaire

illustration n°1 : essai en canal de laboratoire


Cette vidéo a été réalisée par José Vazkez, professeur d'hydraulique à l'ENGEES

illustration n°2 : clip sur la propagation de vagues en mer

transport de sédiment - évolution des fonds

illustration n°1 : essai en canal de laboratoire

Bibliographie

A compléter