Commit 1360e3ad authored by Theophile Terraz's avatar Theophile Terraz
Browse files

work on Tuto1

No related merge requests found
Showing with 65 additions and 32 deletions
+65 -32
doc/users/Tuto1/img/Geo.png

101 KB | W: | H:

doc/users/Tuto1/img/Geo.png

108 KB | W: | H:

doc/users/Tuto1/img/Geo.png
doc/users/Tuto1/img/Geo.png
doc/users/Tuto1/img/Geo.png
doc/users/Tuto1/img/Geo.png
  • 2-up
  • Swipe
  • Onion skin
doc/users/Tuto1/img/Geo2.png

168 KB | W: | H:

doc/users/Tuto1/img/Geo2.png

163 KB | W: | H:

doc/users/Tuto1/img/Geo2.png
doc/users/Tuto1/img/Geo2.png
doc/users/Tuto1/img/Geo2.png
doc/users/Tuto1/img/Geo2.png
  • 2-up
  • Swipe
  • Onion skin
doc/users/Tuto1/img/Seuil.png

38.3 KB

doc/users/Tuto1/img/dl.png

15.7 KB | W: | H:

doc/users/Tuto1/img/dl.png

16.3 KB | W: | H:

doc/users/Tuto1/img/dl.png
doc/users/Tuto1/img/dl.png
doc/users/Tuto1/img/dl.png
doc/users/Tuto1/img/dl.png
  • 2-up
  • Swipe
  • Onion skin
doc/users/Tuto1/img/fev2002.png

62.5 KB | W: | H:

doc/users/Tuto1/img/fev2002.png

61 KB | W: | H:

doc/users/Tuto1/img/fev2002.png
doc/users/Tuto1/img/fev2002.png
doc/users/Tuto1/img/fev2002.png
doc/users/Tuto1/img/fev2002.png
  • 2-up
  • Swipe
  • Onion skin
doc/users/Tuto1/img/frictions.png

47.9 KB | W: | H:

doc/users/Tuto1/img/frictions.png

49.5 KB | W: | H:

doc/users/Tuto1/img/frictions.png
doc/users/Tuto1/img/frictions.png
doc/users/Tuto1/img/frictions.png
doc/users/Tuto1/img/frictions.png
  • 2-up
  • Swipe
  • Onion skin
doc/users/Tuto1/img/frictions2.png

49.5 KB

doc/users/Tuto1/img/ic.png

93.7 KB | W: | H:

doc/users/Tuto1/img/ic.png

97.2 KB | W: | H:

doc/users/Tuto1/img/ic.png
doc/users/Tuto1/img/ic.png
doc/users/Tuto1/img/ic.png
doc/users/Tuto1/img/ic.png
  • 2-up
  • Swipe
  • Onion skin
doc/users/Tuto1/img/mjpk.png

20.9 KB

doc/users/Tuto1/img/param.png

115 KB

...@@ -110,13 +110,14 @@ Auteurs : & Pierre-Antoine Rouby & pierre-antoine.rouby@inrae.fr\tabularnewline ...@@ -110,13 +110,14 @@ Auteurs : & Pierre-Antoine Rouby & pierre-antoine.rouby@inrae.fr\tabularnewline
Pamhyr2 peut être téléchargé ici : \url{https://gitlab.irstea.fr/theophile.terraz/pamhyr}. Pamhyr2 peut être téléchargé ici : \url{https://gitlab.irstea.fr/theophile.terraz/pamhyr}.
Si vous arrivez sur une fenêtre de connexion, cliquez sur "parcourir la liste des dépôts publics" et cherchez l'URL donnée plus haut.
Utilisez le bouton de téléchargement GNU Linux ou Windows en fonction de votre système.
\begin{center} \begin{center}
\includegraphics[width=15cm]{img/dl.png} \includegraphics[width=15cm]{img/dl.png}
\par\end{center} \par\end{center}
Décompressez l'archive et lancez Pamhyr2.
Utilisez le bouton de téléchargement GNU Linux ou Windows en fonction de votre système. Sous Windows, lancez le programme d'installation. Sous Linux, décompressez l'archive et lancez Pamhyr2.
\section{Créer une première étude} \section{Créer une première étude}
...@@ -233,10 +234,9 @@ Vous arrivez sur la fen ...@@ -233,10 +234,9 @@ Vous arrivez sur la fen
Utilisez le bouton \includegraphics[width=0.5cm]{../../../src/View/ui/ressources/add.png} en haut à gauche de la fenêtre pour ajouter une condition limite pour le liquide. Utilisez le bouton \includegraphics[width=0.5cm]{../../../src/View/ui/ressources/add.png} en haut à gauche de la fenêtre pour ajouter une condition limite pour le liquide.
Nous pouvons cliquer sur la nouvelle ligne pour sélectionner la ligne en entier, et double-cliquer pour sélectionner une seule cellule. Nous pouvons cliquer sur la nouvelle ligne pour sélectionner la ligne en entier, et double-cliquer pour sélectionner une seule cellule.
Sélectionnez la cellule \textit{Nom} pour donner un nom à la condition limite.
Ici, nous définirons le débit mesuré lors de la crue de février 2002. Ici, nous définirons le débit mesuré lors de la crue de février 2002.
Vous pouvez nommer cette condition limite "crue2002". Sélectionnez la cellule \textit{Nom} pour donner un nom à la condition limite par exemple "crue2002".
Sélectionner la cellule \textit{Type} et utiliser la combo box pour mettre une loi \textit{Q(t)} : débit en fonction du temps. Sélectionner la cellule \textit{Type} et utiliser la combo box pour mettre une loi \textit{Q(t)} : débit en fonction du temps (hydrogramme).
Sélectionnez la cellule \textit{N½ud} et attribuez cette condition au n½ud amont. Sélectionnez la cellule \textit{N½ud} et attribuez cette condition au n½ud amont.
Les noms des n½uds sont rappelés dans le panneau de droite, avec une vue du réseau. Les noms des n½uds sont rappelés dans le panneau de droite, avec une vue du réseau.
Sélectionnez maintenant la ligne entière et cliquez sur le bouton d'édition \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/edit.png"}. Sélectionnez maintenant la ligne entière et cliquez sur le bouton d'édition \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/edit.png"}.
...@@ -249,12 +249,19 @@ Vous pouvez maintenant voir la courbe de d ...@@ -249,12 +249,19 @@ Vous pouvez maintenant voir la courbe de d
\includegraphics[width=15cm]{img/fev2002.png} \includegraphics[width=15cm]{img/fev2002.png}
\par\end{center} \par\end{center}
Fermer cette fenêtre. Revenez sur la fenêtre \textit{Conditions aux limites}. Fermez cette fenêtre. Revenez sur la fenêtre \textit{Conditions aux limites}.
Ajouter une nouvelle ligne, lui donner un nom, lui donner le type \textit{Z(T)} et l'associer au n½ud aval du réseau. Ajoutez une nouvelle ligne, donnez lui un nom, donnez lui le type \textit{Z(T)} (limnigramme) et associez là au n½ud aval du réseau.
Ouvrez la fenêtre d'édition des conditions aux limites (\includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/edit.png"}). Ouvrez la fenêtre d'édition des conditions aux limites (\includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/edit.png"}).
Ajoutez deux lignes. Dans la première, entrez le temps : 0.00.00 et Z : 15.000. Ajoutez deux lignes et rentrez les valeurs suivantes :
Dans la seconde, saisissez le temps : 1.00.00 et Z : 15.000. \begin{center}
Cela crée une côte constante de l'eau en aval. \begin{tabular}{|c|c|}
\hline
temps & cote \\ \hline
0.00.00 & 15.000\\ \hline
1.00.00 & 15.000\\ \hline
\end{tabular}\\
\par\end{center}
Cela crée une cote constante de l'eau en aval.
Pour le calcul, le solveur extrapolera continûment l'élévation de l'eau si le temps de la simulation dépasse le dernier point de la courbe. Pour le calcul, le solveur extrapolera continûment l'élévation de l'eau si le temps de la simulation dépasse le dernier point de la courbe.
Vous pouvez fermer les fenêtres \textit{Éditer les conditions aux limites} et \textit{Conditions aux limites}. Vous pouvez fermer les fenêtres \textit{Éditer les conditions aux limites} et \textit{Conditions aux limites}.
...@@ -273,7 +280,7 @@ Vous pouvez fermer les fen ...@@ -273,7 +280,7 @@ Vous pouvez fermer les fen
Dans la fenêtre principale, cliquez sur \texttt{[Hydraulique] => [Conditions initiales]} ou sur le raccourci \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/boundary_condition.png"}. Dans la fenêtre principale, cliquez sur \texttt{[Hydraulique] => [Conditions initiales]} ou sur le raccourci \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/boundary_condition.png"}.
Pour démarrer, le code de calcul hydraulique a besoin de connaître le débit et la cote de la surface libre de la rivière en tout point à l'instant initial. Pour démarrer, le code de calcul hydraulique a besoin de connaître le débit et la cote de la surface libre de la rivière en tout point à l'instant initial.
Si vous ne connaissez pas ces conditions initiales, vous pouvez utiliser les boutons \texttt{[Générer une profondeur uniforme]}, \texttt{[Générer un débit uniforme]} ou \texttt{[Générer une cote uniforme]} pour laisser Pamhyr2 estimer une condition initiale à l'aide de la formule de Manning-Strickler. Si vous ne connaissez pas ces conditions initiales, vous pouvez utiliser les boutons \texttt{[Générer une profondeur uniforme]}, \texttt{[Générer un débit uniforme]} ou \texttt{[Générer une cote uniforme]} pour laisser Pamhyr2 estimer une condition initiale.
Cliquez sur \texttt{[Générer un débit uniforme]} et saisissez un débit de 4 m$^3$/s dans la fenêtre contextuelle, et cochez \texttt{[Générer une profondeur]} pour générer une condition initiale de hauteur d'eau basée sur la formule de Manning-Strickler pour le débit donné. Cliquez sur \texttt{[Générer un débit uniforme]} et saisissez un débit de 4 m$^3$/s dans la fenêtre contextuelle, et cochez \texttt{[Générer une profondeur]} pour générer une condition initiale de hauteur d'eau basée sur la formule de Manning-Strickler pour le débit donné.
Vous devriez voir : Vous devriez voir :
...@@ -285,11 +292,11 @@ Vous devriez voir : ...@@ -285,11 +292,11 @@ Vous devriez voir :
% Vous pouvez également utiliser \texttt{[Générer une cote uniforme]} et saisir une cote de 21 m à l'amont et à l'aval avec un débit nul. % Vous pouvez également utiliser \texttt{[Générer une cote uniforme]} et saisir une cote de 21 m à l'amont et à l'aval avec un débit nul.
% Cela revient à créer une bassine et à laisser le solveur la vider pour trouver un écoulement initial satisfaisant. % Cela revient à créer une bassine et à laisser le solveur la vider pour trouver un écoulement initial satisfaisant.
Lors des prochaines simulations, vous pouvez utiliser l'état final de la simulation précédente comme état initial. Lors des prochaines simulations, vous pourrez utiliser l'état final de la simulation précédente comme état initial.
Pour cela, cliquez sur \includegraphics[width=0.5cm]{../../../src/View/ui/ressources/import.png} et retrouvez le résultat sous forme de fichier \textit{.BIN}. Pour cela, cliquez sur \includegraphics[width=0.5cm]{../../../src/View/ui/ressources/import.png} et retrouvez le résultat sous forme de fichier \textit{.BIN}.
Ce fichier de résultats se trouve normalement dans le sous-dossier \textit{\_PAMHYR\_/Hogneau/default-mage}. Ce fichier de résultats se trouve normalement dans le sous-dossier \texttt{\_PAMHYR\_/Hogneau/default-mage}.
Fermer la fenêtre \textit{Conditions initiales} Fermer la fenêtre \textit{Conditions initiales}.
\section{Éditer les coefficients de frottement} \section{Éditer les coefficients de frottement}
...@@ -310,13 +317,13 @@ Vous pouvez utiliser le bouton \includegraphics[width=0.5cm]{../../../src/View/u ...@@ -310,13 +317,13 @@ Vous pouvez utiliser le bouton \includegraphics[width=0.5cm]{../../../src/View/u
Fermez la fenêtre \textit{Strickler}. Fermez la fenêtre \textit{Strickler}.
Dans la fenêtre \textit{Éditer les frottements}, ajoutez quatre lignes avec le bouton \includegraphics[width=0.5cm]{../../../src/View/ui/ressources/add.png} pour créer quatre zones de frottement. Dans la fenêtre \textit{Éditer les frottements}, ajoutez quatre lignes avec le bouton \includegraphics[width=0.5cm]{../../../src/View/ui/ressources/add.png} pour créer quatre zones de frottement.
Chaque zone est définie par un PK de \textit{début} et de \textit{fin} associés à un couple de coefficients de Strickler de \textit{début} et de \textit{fin}. Chaque zone est définie par un PK de \textit{début} et de \textit{fin} associés à un couple de coefficients de Strickler.
Les couples de coefficients de Strickler à l'intérieur d'une zone sont interpolés à partir des couples \textit{début} et \textit{fin}. % Les couples de coefficients de Strickler à l'intérieur d'une zone sont interpolés à partir des couples \textit{début} et \textit{fin}.
Dans notre cas, nous utiliserons des coefficients uniformes par zone. % Dans notre cas, nous utiliserons des coefficients uniformes par zone.
Définissez les zones comme suit : Définissez les zones comme suit :
\begin{center} \begin{center}
\includegraphics[width=15cm]{img/frictions2.png} \includegraphics[width=15cm]{img/frictions.png}
\par\end{center} \par\end{center}
La zone sélectionnée est surlignée en bleu. La zone sélectionnée est surlignée en bleu.
...@@ -325,10 +332,10 @@ Vous pouvez maintenant fermer la fen ...@@ -325,10 +332,10 @@ Vous pouvez maintenant fermer la fen
\section{Modélisation des ouvrages hydrauliques} \section{Modélisation des ouvrages hydrauliques}
Parfois, il peut y avoir des sections en travers dans lesquelles les équations de Saint-Venant ne peuvent pas être utilisées pour modéliser l'écoulement de l'eau. Parfois, il peut y avoir des sections en travers dans lesquelles les équations de Saint-Venant ne peuvent pas être utilisées pour modéliser l'écoulement de l'eau.
Dans ce cas, nous devons définir une autre loi pour relier la côte de l'eau et le débit. Dans ce cas, nous devons définir une autre loi pour relier la cote de l'eau et le débit.
C'est le cas, par exemple, au niveau des ponts lorsque la hauteur d'eau est trop élevée et que l'écoulement passe en charge. C'est le cas, par exemple, au niveau des ponts lorsque la hauteur d'eau est trop élevée et que l'écoulement passe en charge.
Pamhyr2 permet de définir différents ouvrages hydrauliques avec des lois paramétrables. Pamhyr2 permet de définir différents ouvrages hydrauliques avec des lois paramétrables.
Dans notre cas, nous allons devoir modéliser deux ponts par des ouvrages hydrauliques. Dans notre cas, nous allons devoir modéliser deux ponts et un seuil par des ouvrages hydrauliques.
Dans la fenêtre principale, cliquez sur \texttt{[Hydraulique] => [Ouvrages hydrauliques]} pour ouvrir la fenêtre des ouvrages hydrauliques. Dans la fenêtre principale, cliquez sur \texttt{[Hydraulique] => [Ouvrages hydrauliques]} pour ouvrir la fenêtre des ouvrages hydrauliques.
Cliquez trois fois sur le bouton \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/add.png"} pour créer trois ouvrages hydrauliques. Cliquez trois fois sur le bouton \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/add.png"} pour créer trois ouvrages hydrauliques.
% Cliquez deux fois sur le bouton \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/add.png"} pour créer deux ouvrages hydrauliques. % Cliquez deux fois sur le bouton \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/add.png"} pour créer deux ouvrages hydrauliques.
...@@ -339,7 +346,7 @@ D ...@@ -339,7 +346,7 @@ D
\includegraphics[width=15cm]{img/hs.png} \includegraphics[width=15cm]{img/hs.png}
\par\end{center} \par\end{center}
Sélectionnez le seuil aval et cliquez sur \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/edit.png"} pour modifier les lois de cet ouvrage. Les ouvrages hydrauliques sont composés d'ouvrages hydrauliques élémentaires. Vous pouvez combiner les lois de plusieurs ouvrages hydrauliques élémentaires pour créer votre ouvrage. Pour ce seuil, nous n'avons besoin que d'un ouvrage hydraulique élémentaire de type déversoir. Cliquez sur \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/add.png"} pour ajouter un nouvel ouvrage hydraulique élémentaire, donnez-lui le type \textit{seuil} et configurez-le comme suit : Sélectionnez le seuil aval et cliquez sur \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/edit.png"} pour modifier les lois de cet ouvrage. Les ouvrages hydrauliques sont composés d'ouvrages hydrauliques élémentaires. Vous pouvez combiner les lois de plusieurs ouvrages hydrauliques élémentaires pour créer votre ouvrage. Pour ce seuil, nous n'avons besoin que d'un ouvrage hydraulique élémentaire de type déversoir. Cliquez sur \includegraphics[width=0.5cm]{"../../../src/View/ui/ressources/add.png"} pour ajouter un nouvel ouvrage hydraulique élémentaire, donnez-lui le type \textit{Déversoir rectangulaire} et configurez-le comme suit :
\begin{center} \begin{center}
\includegraphics[width=15cm]{img/Seuil.png} \includegraphics[width=15cm]{img/Seuil.png}
...@@ -361,10 +368,30 @@ Cr ...@@ -361,10 +368,30 @@ Cr
Revenez à la fenêtre \textit{Ouvrages hydrauliques} et appliquez la même procédure pour le pont de Thivencelle : Revenez à la fenêtre \textit{Ouvrages hydrauliques} et appliquez la même procédure pour le pont de Thivencelle :
\begin{center} % \begin{center}
\includegraphics[width=15cm]{img/SeuilThivencelle.png} % \includegraphics[width=15cm]{img/SeuilThivencelle.png}
%
% \includegraphics[width=15cm]{img/OrificeThivencelle.png}
% \par\end{center}
\includegraphics[width=15cm]{img/OrificeThivencelle.png} \begin{center}
\begin{tabular}{|c|c|}
\hline
\multicolumn{2}{|c|}{Orifice Thivencelle} \\ \hline
Largeur & 9.07 \\ \hline
cote & 17.28\\ \hline
cote de mise en charge & 20.74\\ \hline
coefficient de débit & 0.4\\ \hline
cote de mise en charge maximale & 9999.999\\ \hline
\end{tabular}
\hspace{1cm}
\begin{tabular}{|c|c|}
\hline
\multicolumn{2}{|c|}{Seuil Thivencelle} \\ \hline
Largeur & 13.0 \\ \hline
cote & 21.74\\ \hline
coefficient de débit & 0.4\\ \hline
\end{tabular}\\
\par\end{center} \par\end{center}
...@@ -374,13 +401,19 @@ Si vous r ...@@ -374,13 +401,19 @@ Si vous r
\section{Paramètres du solveur} \section{Paramètres du solveur}
Dans la fenêtre principale, cliquez sur \texttt{[Exécuter] => [Parameters numériques des solveurs]}. Dans la fenêtre principale, cliquez sur \texttt{[Exécuter] => [Parameters numériques des solveurs]}.
Dans la fenêtre \textit{Paramètres du solveur}, sélectionnez l'onglet \textit{Mage v8}. Dans la fenêtre \textit{Paramètres du solveur}, sélectionnez l'onglet \texttt{[Mage v8]}.
Gardez les paramètres du solveur par défaut, sauf le pas de temps minimum, que vous mettrez à 0.1. Ces paramètres pilotent le comportement du solveur numérique.
En effet, durant la montée du pic de crue, le solveur a besoin de réduire le pas de temps suffisament pour permettre la convergence des itérations. la valeur 999:99:00:00 du temps final indique au solveur de s'arrêter lorsqu'il a atteint un régime permanent.
% Pour accélérer les calculs, nous pouvons également dégrader la précision, à l'aide des facteurs de réduction de la précision. Vous pouvez changer la fréquence d'écriture des résultats dans la ligne \textit{Pas de temps d'écriture dans le fichier .BIN}.
% Les précisions internes du solveur sont de $10{-9}$. Une valeur inférieure à 1 seconde indique que la valeur de la ligne \textit{Pas de temps d'écriture dans le fichier .TRA} s sera prise à la place.
% Cette précision est multipliée par le facteur de réduction de la précision : un facteur de 1000 ramènera donc la précision à $10{-5}$. % En effet, durant la montée du pic de crue, le solveur a besoin de réduire le pas de temps suffisament pour permettre la convergence des itérations.
% Pour utiliser ce facteur de réduction de la précision, il faut donner un \textit{nombre d'itérations à précision maximum} inférieur au \textit{nombre maximum d'itérations} : le solveur va d'abord tenter de converger avec un certain nombre d'itératons à la précision maximum avant de basculer sur une précision dégradée pour le reste des itérations. Pour accélérer les calculs et pour aider le solveur à démarer, nous allons l'autoriser à dégrader la précision, à l'aide des facteurs de réduction de la précision.
Les précisions internes du solveur sont de 10$^{-9}$.
Cette précision est multipliée par le facteur de réduction de la précision : un facteur de 1000 ramènera donc la précision à 10$^{-5}$.
Pour utiliser ce facteur de réduction de la précision, il faut donner un \textit{nombre d'itérations à précision maximum} inférieur au \textit{nombre maximum d'itérations} : le solveur va d'abord tenter de converger avec un certain nombre d'itératons à la précision maximum avant de basculer sur une précision dégradée pour le reste des itérations.
Dans notre cas, c'est nécessaire pour lancer le solveur à partir de la condition initiale calculée par Pamhyr2.
rentrez 1000 dans les trois lignes \textit{facteurs de réduction de la précision}, rentrez 99 pour le \textit{Nombre d'itérations} et 5 pour le \textit{nombre d'itérations à la précision maximum}.
Gardez les autres paramètres du solveur par défaut.
Fermer la fenêtre \textit{Paramètres du solveur}. Fermer la fenêtre \textit{Paramètres du solveur}.
\section{Lancer la simulation} \section{Lancer la simulation}
......
No preview for this file type
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment