Site web | https://www.edf.fr/groupe-edf/inventer-l-avenir-de-l-energie/r-d-un-savoir-faire-mondial/nos-offres/codes-de-calcul/lammi |
---|---|
Extension(s) de fichier | .txt, .prn |
Nombre de dimension | 1 |
Nombre de tronçon | |
Unité(s) de temps | |
Équation | |
Méthode d'analyse numérique | |
Variable(s) aux nœuds | z, h, v, substrat |
Condition(s) de simulation |
Les sorties du logiciel LAMMI pour HABBY sont constitués d’un fichier .txt qui fait référence dans un ordre amont/aval à autant de fichiers .prn qu’il y a de transects.
Chaque transect a une longueur représentative constante pour tous les débits celle-ci est indiquée dans le .txt.
Chaque .prn décrit pour chaque débit simulé un ensemble de verticales fournies de la rive gauche à la rive droite. Le nombre de verticales est variable pour un même transect selon les valeurs de débit.
Dans le fichier .txt le 'chemin' des fichiers .prn n'a pas d'importance seul compte le nom des fichier .prn et le fait que ceux-ci soit présents dans le même répertoire que celui du fichier .txt .
Dans ce répertoire l'utilisateur pourra utilement construire plusieurs fichiers .txt en isolant des .prn constituant un transect, un faciès constitué de plusieurs transects, etc…
Restera alors à créer d'autres .hab dans le même projet en prenant soin de leur indiquer un nom représentatif.
Les calculs d'habitat seront facilités, car les modèles biologiques n'auront besoin d'être sélectionnés qu'une seule fois dans le projet.
Pour pouvoir utiliser dans HABBY les données de sorties hydrauliques et substrats du logiciel LAMMI, il faut convertir ces descriptions en transects, en maillages. Ces maillages seront au choix de l'utilisateur soit de type éléments finis avec des valeurs hydrauliques H,V interpolées aux nœuds (recommandé) soit de type volume fini avec des valeurs hydrauliques H,V réputées constantes pour toute la maille; dans les deux cas la géométrie x,y des mailles est identique. A une différence prés dans le cas de valeur H moyenne aux mailles nulles les mailles 'sèches' n'étant pas conservées dans HABBY.
Le format HABBY de transformation sera directement de type habitat .hab car le substrat est renseigné dans LAMMI sous forme de pourcentages en 8 classes granulométriques.
Deux formats de description du substrat sont possibles aux choix de l'utilisateur soit EDF soit Cemagref (recommandé).
Dans le cas d'une description granulométrique EDF:LAMMI [Code EDF R&D (Cailleux 1954)] celle ci doit être transformée pour répondre à l'une des normes HABBY voisine : le [code Cemagref (Malavoi 1989) ou Code EVHA 2.0 (GINOT 1998)].
Nom de la classe granulométrique | Taille de l’élément (mm) | Code EDF R&D, 8 classes | Code Cemagref, 8 classes |
---|---|---|---|
Argiles | <0.0039 | s1 | s1 |
Limons | 0.0039-0.0625 | s2 | |
Sables fins | 0.0625-0.5000 | s2 | s3 |
Sables grossiers | 0.5-2.0 | ||
Graviers fins | 2-8 | s3 | s4 |
Graviers grossier | 8-16 | ||
Cailloux fins | 16-32 | s4 | s5 |
Cailloux grossiers | 32-64 | ||
Pierres fines | 64-128 | s5 | s6 |
Pierres grossières | 128-256 | s6 | |
Blocs | 256-1024 | s7 | s7 |
Rochers | >1024 | s8 | s8 |
Localement si $\mathit{S_{EDF,k}}$ et $\mathit{S_{Cemagref,k}}$ représentent des % de substrat de la classe $\mathit{S_k}$ pour effectuer les transformations, il faut effectuer les opérations suivantes.
$\mathit{S_{Cemagref,2}=S_{EDF,1}}$, $\mathit{S_{Cemagref,3}=S_{EDF,2}}$, $\mathit{S_{Cemagref,4}=S_{EDF,3}}$, $\mathit{S_{Cemagref,5}=S_{EDF,4}}$, $\mathit{S_{Cemagref,6}=S_{EDF,5}+0.5\times S_{EDF,6}}$, $\mathit{S_{Cemagref,7}=0.5\times S_{EDF,6} +S_{EDF,7}}$, $\mathit{S_{Cemagref,8}=S_{EDF,8}}$
La répartition de $\mathit{S_{EDF,6}}$ dans deux classes $\mathit{S_{Cemagref,6}}$ et $\mathit{S_{Cemagref,7}}$ s'explique, car les classes de substrat du protocole LAMMI ne sont pas identiques à celles du tableau présenté ci dessus.
La figure ci dessus schématise la construction d'un maillage HABBY pour un transect LAMMI comprenant une, deux ou trois verticales.
Pour chaque verticale LAMMI 8 mailles HABBY sont crées et toutes comportent la même information substrat issue de cette verticale.
Pour un noeud i : si l'on note $\mathit{H_{i}}$ la hauteur d 'eau, $\mathit{V_{i}}$ la vitesse moyenne, $\mathit{la_{i}}$ la largeur de representativité d'une verticale :
Dans le cas d'un calcul en éléments finis :
Aux bords $\mathit{(H_{0},V_{0})=(H_{1},V_{1})=(0,0)}$
Dans l'axe d'une verticale LAMMI $\mathit{(H_{2},V_{2})=(H_{3},V_{3})=(H_{7},V_{7})}$
Entre deux verticales LAMMI $\mathit{(H_{5},V_{5})=(H_{6},V_{6})=(\frac{H_{A}\times la_{B}+H_{B}\times la_{A}}{la_{A}+la_{B}},\frac{V_{A}\times la_{B}+V_{B}\times la_{A}}{la_{A}+la_{B}})}$
Dans le cas d'un calcul en volumes finis :
Systématiquement pour les 8 mailles construites autour d'une verticale LAMMI toutes les valeurs moyennes hydrauliques H,V des mailles sont celles de la verticale.
Une fois le maillage réalisé HABBY va supprimer les mailles de hauteur d'eau moyenne nulle, l'aire indiquée dans un calcul d'habitat par exemple correspond à une aire mouillée.
Pour comparer des résultats de calcul d'habitat il faut choisir un calcul en volume fini, puis il faut demander dans HABBY un calcul du substrat en pourcentage ainsi que l'opération s'effectue dans LAMMI.
Soit à calculer la valeur d'habitat pour le substrat dans une maille d'indice i : $\mathit{SI_{i,S}}$.
\[{SI_{i,S}(S_{i,1},S_{i,2},..S_{i,8})}=\frac{\sum_{k=1}^8 S_{i,k}\times SI_S(S_k)}{100} \]
Dans cette équation $\mathit{S_{i,k}}$ représente le % de substrat de la classe $\mathit{S_k}$ dans la maille i, le substrat étant décrit par 8 classes granulométriques k ∈ $\mathit{[1,8]_N}$.
Reste qu'en raison de la transformation du substrat les résultats ne pourrons pas être exactement identiques.
Si l'on effectue une comparaison alors que les calculs sont effectués en éléments finis, la différence sera mécaniquement plus importante.
Car HABBY crée des nœuds supplémentaires en particulier aux bords, par interpolation entre la rive H=0 et la hauteur d'eau H de la première verticale en eau.
Ces nœuds interpolés ont des hauteurs et des vitesses non utilisées dans les calculs d'habitat de LAMMI; ces nouvelles valeurs hydrauliques vont être utilisés par HABBY pour les calculs d'habitat, ce qui va introduire des différences numériques.
S'il souhaite vérifier pour un tronçon les valeurs de hauteur et de vitesse moyenne, l'utilisateur pourra utiliser l'outil hydrosignature d'HABBY. Rappelons que pour obtenir la vitesse moyenne d'un tronçon, cet outil effectue une pondération des vitesses moyennes par les volumes des mailles.
Si les utilisateurs souhaitent utiliser HABBY pour les calculs d'habitats pour bénéficier de la grande bibliothéque de modèles biologiques d'habitat d'HABBY : il est conseillé de décrire sur le terrain le substrat en mode Cemagref, car c'est le format de description des modèles biologiques Français, et donc de n'effectuer les calculs d'habitat que dans HABBY. Lors de l'intégration des données LAMMI dans HABBY il faut donc préciser que le substrat est de type Cemagref, et choisir le mode d'équation en éléments finis ce qui permet notamment selon notre expertise, un calcul plus précis des valeurs d'habitats des jeunes stades en bordure (une maille en bordure comportera deux verticales avec des hauteurs d'eau nulles).