Table of Contents

LAMMI

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

Principes d'utilisation du fichier .txt et des fichiers .prn

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.

Intégration dans HABBY

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.

Comparaison des résultats de calculs d'habitats LAMMI HABBY

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.

Recommandation aux utilisateurs de LAMMI dans HABBY

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).

Cas particulier de hauteurs d'eau nulles provenant de LAMMI, dans le cas d'un choix éléments finis

La figure ci dessus schématise le cas particulier d'un maillage avec une verticale LAMMI de hauteur d'eau nulle $\mathit{H_{A}=0}$, deux solutions sont possibles et c'est la solution 2 qui est retenue par HABBY dans le cas d'un choix éléments finis.