Modélisation de la propagation de la houle
Par SWAN
Simulating Waves Nearshore
Introduction
C’est un modèle numérique spectral de troisième génération conçu pour simuler des vagues
évoluant dans les régions côtières, les lacs et dans les estuaires tenant compte de conditions
données du vent, de la bathymétrie et du courant.
Il a été conçu par l’université de technologie de Delft (Hollande).
Le modèle est basé sur l'équation spectrale de conservation de l’action de la houle qui est
résolue suivant un schéma implicite en différences finies.
La modélisation de l’évolution de la houle du large à la côte, les processus à prendre en
considération sont :
- La propagation.
- Les interactions conservatrices entre ondes.
- Les sources d’apport d’énergie. (Le vent)
- Les sources de dissipation de cette énergie. (Le frottement)
SWAN permet une simulation en mode uni ou bidimensionnel, stationnaire ou non-
stationnaire, en coordonnées Cartésiennes (domaines de faible emprise) ou Sphériques
(pour les petites et grandes échelles).
Elle donne aussi en option trois sortes de grilles pour le maillage du domaine d’étude :
- Une grille régulière.
- Une grille curvilinéaire
- Une grille non-structurée.
SWAN fournit les quantités suivantes comme outputs :
- Le spectre unidimensionnel ou bidimensionnel.
- La hauteur significative de la houle.
- Les différentes périodes : période moyenne, période de pic et période énergétique.
- La direction moyenne de la houle.
- L’étalement directionnel.
Fichier de contrôle SWAN
Le logiciel SWAN n’admet pas une interface graphique. Il est commandé à travers une
fenêtre DOS, et il a donc besoin d’un fichier de contrôle (fichier texte) sous forme de
programme dans un langage précis (instauré par les concepteurs).
- Les informations à entrer dans ce fichier sont principalement :
- Le nom qu’on veut attribuer à notre projet
- La génération de SWAN à utiliser parmi les trois existants.
- Le mode : unidimensionnel ou bidimensionnel.
- Le mode stationnaire ou non : le mode non stationnaire est utilisé lorsqu’on dispose
de l’évolution des caractéristiques (vent, bathymétrie, courant) du milieu sur une
période donnée.
- Le système de coordonnées (cartésiennes ou sphériques)
- Les processus à activer (friction du fond, amplification par le vent, déferlement…) ;
- Les répertoires des fichiers contenant les données d’entrée : bathymétrie, vents,
courant
- Les conditions aux limites
- Les options de sorties
Liste des commandes les plus utilisées
Commandes de démarrage
PROJEC
: Titre du problème à calculer
T
SET : Définit les valeurs de certains paramètres généraux
MODE : Définit le mode stationnaire ou instationnaire aussi
le mode 1D ou mode 2D.
COORD : Définit le système de coordonnées cartésiennes ou sphériques
Commandes pour la description du modèle
Commandes pour la grille de calcul :
CGRID : définit les dimensions de la grille de calcul
READGRID
lit une grille de calcul curviligne ou non structurée
:
Commandes pour les champs de saisie :
INPGRID définit les dimensions du fond, le niveau d'eau, le courant et les grilles de
: frottement
READINP lit les champs de saisie
WIND active l'option vent constant
Commandes pour les conditions initiales et limites :
BOUND définit la forme des spectres à la limite de la grille géographique
BOUNDSPECT définit les spectres (paramétriques) à la limite de la grille géographique
Commandes pour la physique :
GEN1 SWAN fonctionne en mode première génération
GEN2 SWAN fonctionne en mode deuxième génération
GEN3 SWAN fonctionne en mode troisième génération
FRICTIO
active la dissipation par frottement du fond
N
Syntaxe d’un fichier SWAN
PROJECT 'nom' 'nr'
'titre1'
'titre2'
'titre3'
Avec cette commande obligatoire, l'utilisateur définit un certain nombre de chaînes pour
identifier l'exécution SWAN
est le nom du projet, d'une longueur maximale de 16
'nom'
caractères.
'nr' est l'identification de l'exécution
'titre1' est le titre du projet
'titre2' identique à 'title1' .
'titre3' identique à 'title1' .
SET NAUTICAL / CARTESIEN
Convention nautique : la direction d'où viennent le vent ou les vagues, mesurée dans le sens
des aiguilles d'une montre à partir du nord géographique.
Convention cartésienne : c'est-à-dire la direction vers laquelle pointe le vecteur, mesurée
dans le sens inverse des aiguilles d'une montre à partir de l’axe x positif de ce système.
MODE
STATIONARY NONSTATIONARY
MODE
ONEDIMENSIONAL TWODIMENSIONAL
COORDINATES
CARTESIAN
COORDINATES
SPHERICAL
CARTESIA
Toutes les distances sont en m. Les coordonnées sont données par rapport à
N
à x et y choisis par l'utilisateur.
Toutes les coordonnées des emplacements et les tailles des grilles
SPHERICAL
géographiques sont données en degrés ;
Grille de calcul
CGRID
Avec cette commande obligatoire, l'utilisateur définit l'emplacement géographique, la taille, la
résolution et l'orientation de la grille de calcul dans le système de coordonnées du problème
dans le cas d'une grille de calcul uniforme et rectiligne, d'une grille curviligne ou d'un
maillage non structuré.
REGular [xpc] [ypc] [alpc] [xlenc] [ylenc] [mxc] [myc]
CGRID < CURVilinear [mxc] [myc] (EXCeption [xexc] [yexc])
UNSTRUCtured
cette option indique que la grille de calcul doit être considérée comme
REGular
uniforme et rectangulaire.
cette option indique que la grille de calcul doit être considérée comme
CURVilinear
curviligne.
L'utilisateur doit fournir les coordonnées des points de la grille avec la
commande
READGRID COOR .
UNSTRUCture cette option indique que la grille de calcul doit être considérée comme
d non structurée.
L'utilisateur doit fournir les coordonnées des sommets et la numérotation
des triangles avec la table de connectivité associée aux sommets avec la
commande
READGRID NON CONSTRUCT .
localisation géographique de l'origine de la grille de calcul dans le
[xpc]
problème système de coordonnées (x coordonnée, en m)
Par défaut : [xpc] = 0,0 (coordonnées cartésiennes).
Dans le cas de coordonnées sphériques, il n'y a pas de valeur par défaut,
l'utilisateur doit donner une valeur.
localisation géographique de l'origine de la grille de calcul dans le
[ypc]
problème système de coordonnées ( y coordonnée, en m)
Par défaut : [ypc] = 0,0 (coordonnées cartésiennes).
Dans le cas de coordonnées sphériques, il n'y a pas de valeur par défaut,
l'utilisateur doit donner une valeur.
[alpc] direction du positif x axe de la grille de calcul (en degrés, cartésien
convention). En mode 1D, [alpc] doit être égal à la direction [alpinp]
(voir commande INPGRID ).
Par défaut : [alpc] = 0,0.
Longueur de la grille de calcul en x direction (en m). Dans le cas de formes
[xlenc]
sphériques les coordonnées [xlenc] sont en degrés.
[ylenc] Longueur de la grille de calcul en y direction (en m). En mode 1D, [ylenc]
devrait être 0. Dans le cas de coordonnées sphériques, [ylenc] est en
degrés.
nombre de mailles dans la grille de calcul dans x direction pour un
[mxc]
maillage uniforme, rectiligne direction pour une grille curviligne.
nombre de mailles dans la grille de calcul dans y direction pour un
[myc]
maillage uniforme, rectiligne direction pour une grille curviligne.
REGular [xpc] [ypc] [alpc] [xlenc] [ylenc] [mxc] [myc]
| -> ADCirc
|
READgrid UNSTRUCtured < TRIAngle |
| > 'fname'
| EASYmesh |
les informations de grille nécessaires sont lues à partir du fichier fort.14 tel
ADCIRC
qu'utilisé par ADCIRC.
Ce fichier contient également les informations de profondeur qui sont
également lues.
les informations de grille nécessaires sont lues à partir de deux fichiers
TRIANGLE
produits par Triangle.
Les fichiers .node et .ele sont obligatoires. Le nom de base de ces fichiers
doit être indiqué avec le paramètre 'fname' .
EASYMESH les informations de grille nécessaires sont lues à partir de deux fichiers tels
que produits par Easymesh. Les fichiers .n et .e sont obligatoires. Le nom de
base de ces fichiers doit être indiqué avec le paramètre 'fname' .
'nom de nom de base des fichiers requis, c'est-à-dire sans extension. Uniquement
famille' destiné à
Triangle et Easymesh.
INPGRID BOTTOM
REGular [xpinp] [ypinp] [alpinp] [mxinp] [myinp] [dxinp] [dyinp] |
| |
< CURVilinear [stagrx] [stagry] [mxinp] [myinp] >&
| |
| UNSTRUCtured
[xpinp] situation géographique (x coordonnée) de l'origine de la grille d'entrée dans
coordonnées du problème (en m) si des coordonnées cartésiennes sont utilisées
ou en degrés si
des coordonnées sphériques sont utilisées (voir commande COORD ).
Par défaut : [xpinp] = 0. Dans le cas de coordonnées sphériques, il n'y a pas de
valeur par défaut, la
l'utilisateur doit donner une valeur.
[ypinp] situation géographique ( y coordonnée) de l'origine de la grille d'entrée dans
coordonnées du problème (en m) si des coordonnées cartésiennes sont utilisées
ou en degrés si
des coordonnées sphériques sont utilisées (voir commande COORD ).
Par défaut : [ypinp] = 0. Dans le cas de coordonnées sphériques, il n'y a pas de
valeur par défaut, la
l'utilisateur doit donner une valeur.
[alpinp] direction du positif x axe de la grille d'entrée (en degrés, convention cartésienne).
Voir la commande COORD .
Par défaut : [alpinp] = 0.
nombre de mailles dans x direction de la grille d'entrée (ce nombre est un de
[mxinp]
moins
que le nombre de points de la grille dans cette direction !).
[myinp nombre de mailles dans y direction de la grille d'entrée (ce nombre est un de
] moins
que le nombre de points de la grille dans cette direction !).
En mode 1D , [myinp] doit être 0.
[dxinp] taille des mailles en x direction de la grille d'entrée,
en m dans le cas de coordonnées cartésiennes ou
en degrés si des coordonnées sphériques sont utilisées, voir la
commande COORD .
[dyinp] taille des mailles en y direction de la grille d'entrée,
en m dans le cas de coordonnées cartésiennes ou
en degrés si des coordonnées sphériques sont utilisées, voir la
commande COORD .