Risques littoraux majeurs - Webissimo

6 ème journées scientifiques et techniques du CETMEF– Paris– Décembre 2006 – Laborie et Heurtefeux – modélisation et observation de submersionde dunes en Languedoc -Roussillonutilisés lors du calcul courantologique. Comme on s'y attendait, les forces de radiation selon l'axelongitudinal (parallèle à la côte) sont quasiment nulles.IV.3 MODELISATION COURANTOLOGIQUE TENANT COMPTE DES COURANTS DEHOULESAfin de calculer la surélévation du plan d'eau liée à la houle (set-up de houle), le même maillage a étéutilisé dans le cadre d’une simulation numérique des courants à l'aide du code de calcul TELEMAC quirésout les équations de Saint-Venant bidimensionnelles et s'appuie sur une résolution par la méthode deséléments finis.Sur la frontière maritime, la cote de surface libre est imposée à 0,9 m NGF69 ; les forces de radiationcalculées par ARTEMIS sont imposées sur l'ensemble du domaine (de manière physique, logiquement, ellessont non-nulles que du côté maritime). Des tests de sensibilité sur le Strickler et la viscosité ont montré queles résultats n'étaient pas sensibles à ces paramètres du côté maritime. Côté lagune, la sortie est libre.Figure 6 : Profil transversal de la surface libre calculée par TELEMAC (en noir foncé) ; en pointillé, enfiltrant les contraintes de radiation lorsque la hauteur d’eau est inférieur à 1 m.La cote de la surface libre calculée par TELEMAC un peu en amont de la crête est de 2,06 m. Le set-up dehoule avoisinerait donc les 1,10 m. Ce résultat n'est pas en désaccord avec les observations. En effet, le profilmodélisé est celui pris dans le point dépressionnaire ou point faible de la dune, là où se produit la brèche etégalement le transport solide maximal conduisant à la formation des « washover fans ». Suite auxévénements de l'hiver 2003/2004, l'étage supérieur à 2,0 m de la dune a été arasé et les cotes situées entre 1,5et 2,0 m très érodées ; une cote de surface libre de 2,06 au droit de la dune semble aller dans le sens de cesobservations. En effet, elle implique un écoulement dans la brèche dont la hauteur d'eau atteindrait 55 cm etun écoulement transcritique sur la pente côté lagune dont les vitesses seraient comprises entre 1,5 m/s et 2,3m/s. Dans ces conditions, le calcul rapide du critère de Shields pour un sédiment de diamètre moyen 0,35mm conclut effectivement à un phénomène d'érosion.Un calcul dans lequel les forces de radiation ont été filtrées pour des hauteurs d'eau inférieures à 1 m aégalement été mené. La cote de surface libre obtenue est représentée en pointillés sur la figure 6. De fait, seulle premier pic de force de radiation est alors pris en considération. Le set-up de houle calculé par TELEMACest alors de 10 cm et la dune n'est pas franchie par l'écoulement.On en déduit que c'est la précision sur le calcul des forces de radiation liée à la houle dans la zone la plusproche du rivage qui détermine le franchissement du cordon dunaire. Cependant, c'est aussi dans cette zoneque, du fait des hauteurs d'eau les plus faibles, l'imprécision sur la contrainte de radiation est la plus forte.Par ailleurs, l'étude des récifs (cf. [8]) a montré que les courants de houle gouvernent le set-up de houle etque ce set-up décroît lorsque la hauteur d'eau au-dessus des récifs augmente. De fait, deux questions seposent concernant le calcul précédent : se situant dans une zone de hauteur d'eau faibles, quelle fiabilité peutonaccorder aux forces de radiation calculées? Ensuite, le set-up conduisant à une augmentation de la hauteurd'eau au-dessus de la dune et au-dessus de la zone ayant généré les contraintes de radiation, ne peut-ons'attendre à ce qu'elles soient atténuées ? Autrement dit : ne faut-il pas envisager un calcul interactif entre lahoule et le courant pour ne pas courir le risque de surestimer le set-up de houle près du rivage ?