Groupe de travail du MAPMO

Jeudi 4 novembre 10h-12h

Bilan sur la réunion de mi-novembre à Orléans

Jeudi 28 octobre : pas de réunion, interruption pédagogique

Jeudi 21 octobre - 10h-12h

Discussions sur FullSWOF avec Céline Bonhomme.

Jeudi 14 octobre - 10h-12h

Avancées sur la relation topographie / condition de bord.

Point sur les conditions périodique : voir le fichier.

Jeudi 7 octobre - 10h-12h

Code : bilan sur la question de l'initialisation : FullSWOF prend les valeurs au centre de la cellule, et considère la valeur comme constante sur toute la cellule.

Jeudi 30 septembre - 10h-12h

FullSWOF : question de l'initialisation sur les cellules fictives : la hauteur restait constante.

Jeudi 23 septembre - 10h-12h

Bilan sur le code - pb dû à la variable "verif" pas réinitialisée.

Intérêts de la classe choix : voir le fichier.

Lundi 6 septembre - 14h  / mercredi 9 septembre - 16h

Répétitions soutenance Pierre-Antoine

Jeudi 2 septembre - 10h-12h

Présents : Pierre-Antoine Ksinant, Carine Lucas, Christian Laguerre, Frédéric Darboux

Rapport de stage Pierre-Antoine.

Discussion sur la grille de FullSWOF_2D.

Vendredi 13 aout - 14h-16h

Présents : Pierre-Antoine Ksinant, Carine Lucas, Christian Laguerre, Frédéric Darboux

Bilan sur les avancées de Pierre-Antoine

Jeudi 29 juillet - 10h-12h

Présents : Pierre-Antoine Ksinant, Christian Laguerre, Frédéric Darboux

Recherche du bug sur la création d'eau ...

Jeudi 22 juillet - 14h-16h

Présents : Pierre-Antoine Ksinant, Carine Lucas, Christian Laguerre, Frédéric Darboux

Point sur le plan du rapport de stage de Pierre-Antoine.

Discussion sur les changements du parser et la forme, plus souple, du fichier de paramètres d'entrée.

Jeudi 15 juillet - 10h-12h

Présents : Pierre-Antoine Ksinant, Carine Lucas, Christian Laguerre, Frédéric Darboux

Cas tests pour sillons

Frédéric a proposé une première série de cas tests. Suite à discussion, des modifications et ajouts sont effectués ; Frédéric les intègre.
PDF LATEX

Format de fichier

Actuellement FullSWOF_2D utilise un format ASCII compatible avec Gnuplot pour les entrées, les sorties intermédiaires et les sorties finales. Christian propose d'utiliser le format HDF, qui est standard en climatologie et qui facilite la gestion des données.

Frédéric suggère de rajouter ce format pour pour les entrées, les sorties intermédiaires et les sorties finales, tout en gardant la possibilité d'avoir le format Gnuplot. Cela veut dire que l'utilisateur pourra choisir le format qu'il souhaite.

Forge

Christian compare les caractéristiques et fonctions de la forge CNRS avec SourceForge pour vérifier que la forge CNRS répondra à nos besoins.

Changements dans le code FullSWOF_2D

Christian a modifié le code (variables locales et globales). Pierre-Antoine lui envoie son code actuel ; Christian fait un diff et renvoie à Pierre-Antoine le code mis à jour ; Pierre-Antoine vérifie que les résultats ne sont pas altérés.

Condor

Christian écrit au CRI pour savoir s'il peut installer Condor sur les salles de TP.

GARS

Frédéric et Pierre-Antoine regardent comment paramétrer le bruit avec la librairie FFTW.

Jeudi 10 juin - 10h-12h

Résultats de Pierre-Antoine, comparaison 1D pas de temps fixe / 2D cfl fixe.

Jeudi 3 juin - 10h-12h

Travail avec Pierre-Antoine sur le périodique.

Jeudi 27 mai - 10h-12h

Travail avec Pierre-Antoine sur le périodique.

Jeudi 20 mai - 10h-12h

Le point sur les codes FullSWOF :

• FullSWOF_2D

• licence : relancer la personne du FIST
• rapport de bug suite au mail de Minh :il faut que la personne explique comment elle a trouvé le bug et éventuellement propose une correction.
  Cette information doit être envoyé à tous les utilisateurs, puis validée par Christian qui fait une nouvelle version dans tags.
• bilan sur qui modifie quoi : Minh, il faut qu'il envoie ses modifs !
  Pierre-Antoine code les conditions aux bords périodiques : finalement, il ne va que créer une nouvelle classe, pas d'impact sur le reste du code.
• modifications que Christian va faire : liste de priorités - récupérer ce qu'a fait Olivier dans sa thèse pour faire un manuel et une doc des Benchmarks

• FullSWOF_1D : Carine : version de base, sans contrôle, mise sous la forme du 2D.
  A faire : tester avec les cas tests d'Olivier en fortran pour valider.
  Dépôt également prévu, mais sans le BRGM puisqu'il n'y a pas d'infiltration.

Jeudi 15 avril - 10h-12h

• Finance Methode
• Stage de Pierre-Antoine : le point
• Résultats reçus sur Etape 1 : le point
• Poste modélisation INRA
• Fournisseur sièges du Mapmo = ?
• Survie du wiki ?
• edt 25 mai

Jeudi 1er avril 15h-17h

visite de Pierre-Antoine Ksinant, stagiaire de l'INRA qui va travailler sur FullSWOF_2D.

Mardi 23 mars 14h-15h30

• Frédéric et Carine : exposé à Nimes
• Planning de la prochaine réunion METHODE du 25 mai.

Réunion du jeudi 11 mars - 9h30-12h

• Frédéric : étape 1
• Christian : résultats de R. Joubaud sur paraview. Il faut un script pour avoir les hauteurs d'eau.
• Carine : réunion IHP 25 mai. La salle est réservée, il faut fixer le nombre de personnes pour le buffet => mail général.
• Pas encore de nouvelles de C. Terrasse au CNRS pour le dépôt du code.

Jeudi 4 mars : Pas de réunion car journée interne MAPMO

Réunion du jeudi 25 février - 9h30-12h

• étape 1 (Frédéric) : retour des différentes personnes qui ont testé leurs codes
  • Pablo : OK mais il faut peut-être revoir la topo
  • Ulrich : OK pour les tests qu'il a faits.

Demander à tout le monde de bien envoyer les résultats sous forme de fichier texte, et non pas une image ....

• détermination de la date de la réunion de printemps à Paris : nouveau Doodle restreint à 9 dates.

Jeudi 18 février : Pas de réunion car congés

Réunion du jeudi 11 février - 9h30-12h

Dernières modifications pour le dépôt du dossier de déclaration d'auteurs de FullSWOF_2D :

• Licence CéCILL C et rien d'autre (il semblerait qu'elle soit valable à l'étranger)
• Papiers Ulrich et Marie : OK
• Désaccord avec le BRGM : rédaction d'un historique + rapports de stages pour que le CNRS puisse juger.

Point sur le code avec Ulrich : il va utiliser la dernière version de FullSWOF_2D pour l'étape 1.

Jeudi 4 février : Pas de réunion car journée Math et Géosciences

Réunion du jeudi 28 janvier - 10-12h

Stéphane, Frédéric,  Carine, Christian, Charles (stagiaire 3ème)

Licence FullSWOF _2D

En attente de la réponse du BRGM.

Il faut aussi quelques papiers d'Ulrich et de Marie

Comparaison des codes sur Etape 1

Ulrich va faire tourner le code 1D. Une fois les conditions de simulation fixées et les résultats reçus, Christian fera tourner FullSWOF_2D avec le même paramétrage et comparera les résultats 2D et 1D pour vérifier qu'il n'y a pas d'écart.

Stéphane envoie un mail à Ulrich pour confirmation.

Stage INRA-Mapmo

Le stagiaire disposera aussi d'un bureau au Mapmo pour travailler avec Christian. Il faut lui créer un compte au Mapmo afin de pouvoir rebondir depuis l'INRA (ou ailleurs) afin qu'il puisse accéder au CCSC. L'INRA lui fournira un PC portable.

Frédéric vérifie avec le stagiaire si son stage doit se prolonger après sa soutenance.

Le stagiaire modifiera GARS pour :

• pouvoir changer la direction de la pente générale
• avoir aussi des sorties en x y z (option en plus de l'actuel i j )

Il adaptera FullSWOF_2D pour :

• travailler avec des conditions de bord périodique
• arrêter le calcul automatiquement lorsque l'état d'équilibre est atteint.

Il enchaînera la sortie de GARS en entrée de FullSWOF_2D.

Un script devra être écrit pour lancer les calculs paramétriques (est-ce qu'il existe un outil dédié à cette tâche ?). Les calculs paramétriques débuteront en jouant sur :

• direction de la pente
• amplitude de la pente
• direction des sillons
• longueur d'onde des sillons (attention à avoir un nombre entier de sillons dans la topo !)
• amplitude des sillons
• quantité d'eau sur la surface

Dans un deuxième temps, du bruit sera ajouté sur la surface.

La sortie principale sera la répartition des flux (pourcentage dans la direction des sillons et pourcentage dans la direction de la pente générale).

A partir du calcul paramétrique, il formulera une loi phénoménologique permettant, à partir des paramètres, de prédire cette répartition des flux.  Cette loi sera comparée au modèle de Véronique Souchère (et d'Ingrid Takken, éventuellement).

Il poursuivra le travail d'Ulrich (modèle de frottement anisotrope 1D) en l'étendant au 2D. Il comparera les résultats du calcul paramétrique de FullSWOF_2D avec les résultats du modèle de frottement anisotrope 2D.

Stage LIFO de parallélisation de FullSWOF_2D

Il n'y a pas eu de candidat pour ce stage (qui est en fait un travail sur projet). Mais le LIFO jugeant ce travail pertinent pour tester leur librairie de parallélisation, le travail sera effectué après la mi-février, probablement par un doctorant.

FullSWOF_2D sur le bassin de la Loire

Un projet Région va être déposé. le CCSC devrait suffire pour un calcul "temps réel" avec une maille de 10 m de coté. Frédéric se renseigne sur les données topo et occupation du sol.

Date de la prochaine rencontre METHODE

On prévoit une journée à l'IHP fin avril. Stéphane fait un Doodle et réserve la salle. Afin de maximiser notre temps d'échange, on fera un buffet pour le déjeuner.

Nouvel horaire des réunions du GTmapmo

Afin que François puisse être présent, nous les commencerons à 9h30

Réunion du jeudi 21 janvier - 10-12h

Ordre du jour :

quelques questions à discuter sur le code

Réunion du jeudi 14 janvier - 10-12h

Ordre du jour :

Erosion : présentation de Minh.

Réunion du jeudi 10 décembre - 10-12h

Ordre du jour :

• Code Fullswof : présentation de la version de la licence : structure générale, directives pour les développeurs - prévenir le BRGM (Olivier C., Pierre, Marie, Minh)

Pas de réunion le 3 décembre (journée du pôle 3 - contrôle)

Réunion du jeudi 26 novembre - 10-12h

Ordre du jour :

• Fullswof : CR de la réunion du 20/11 avec C. Terrasse
• Exposé Olivier à Cascimodot le 11/12
• Stage Fullswof avec sillons : quelles conditions de bords ? quels calculs ?

Pas de réunion le 19 novembre (journée des doctorants à Blois)

Réunion du jeudi 5 novembre - 10-12h

Ordre du jour :

code Fullswof : licence et développements futurs envisagés par le BRGM dans le cadre du projet RISKERO.

Réunion du jeudi 11 juin - 10-12h

Ordre du jour :

to be defined

Réunion du jeudi 4 juin - 10-12h

Ordre du jour :

to be defined

Réunion du jeudi 28 mai - 10-12h

Ordre du jour :

to be defi

Réunion du jeudi 20 mai - 10-12h

Ordre du jour :

discussions sur les modeles d'érosion (présentation de Carine Lucas)

Réunion du jeudi 14 mai - 10-12h

Ordre du jour :

discussions sur les articles et la presentation d'O. Delestre au congres SMAI 2009.

Réunion du jeudi 7 mai - 10-12h

Ordre du jour :

Presentation Modeles erosion par Marie Martin (stage BRGM).

Réunion du jeudi 16 avril - 10-12h

Ordre du jour :

discussions ruissellement avec Christophe Josserand.

Réunion du jeudi 19 février - 10-12h

Ordre du jour :

Yves Nédélec pour faire le point des contributions du Cemagref

Réunion du jeudi 12 février - 10-12h

Ordre du jour :

Emmanuel Mouche et ses collègues du LSCE (CEA - Gif-sur-Yvette) nous présente leurs travaux sur la modélisation du ruissellement et de l'infiltration.

Pierre Sochala et Marie Rousseau viendront discuter d'un projet de thèse au BRGM

Réunion du jeudi 29 janvier - 10-12h

Ordre du jour :

Hector Gomart nous présentera son "générateur de topographie"

Réunion du jeudi 22 janvier - 10-12h

Christian Laguerre (Mapmo) nous a présenté ses travaux aux LSCE.

Pierre Sochala (BRGM) nous a présenté ses travaux de thèse : Méthodes numériques pour les écoulements souterrains et couplage avec le ruissellement.

Exposé

François et Carine font un bref compte-rendu de la réunion à l'ENPC (21/01/2009). L'approche "ex-Boussinesq" ne semble pas être réalisable (trop forte complexité). Difficulté de comparaison NS-surface libre et Saint-Venant : compatibilité des conditions aux limites (problème difficile, mais intéressant).

Réunion du jeudi 15 janvier - 10-12h

Ordre du jour :

Réunion du jeudi 8 janvier - 10-12h

Frédéric, Carine, Stéphane, Olivier, Christian

FullSWOF

Il faut s'assurer de la convergence des codes pour le tronc commun (Ulrich, Marie)

Générateur de surface

Frédéric présente le code et des réalisations de gegensurf. Actuellement, gegensurf génère correctement la pente générale et la sinusoïde.  La composante bruit est programmée, elle doit à présent être testée. On discute de la façon de paramétrer ; on décide qu'il faut un paramétrage physique (taille des cellules, etc.). Frédéric contacte Ulrich pour s'assurer que les entrées et sorties sont similaires à celle de FullSWOF

Réunion du jeudi 18 décembre - 10-12h

Participants : Nicolas Seguin, Frédéric Darboux, Carine Lucas, Olivier Delestre, François James, Stéphane Cordier

Le point :

• Générateur de surface : pour janvier (Frédéric). Prévoir d'inviter Hector Goma.
• Manip Etape 0 : en cours.
• Manip Etape 1 : OK pour faire les 4 cas extrêmes. Frédéric contacte Michaël Antoine.
• Manip Etape 2 : Frédéric va prendre contact avec Cédric Legout pour planification.
• Articles Etape 0 : un article sur "modèles", un article sur "modèles + expé".
• Biblio : à mettre à jour (Olivier).
• Calcul paramétrique avec conditions de bords périodiques : à voir pour la machine.
• Terrain : sortie à planifier.

Exposé de Nicolas Seguin

Fichier à attacher.

Réunion du jeudi 11 décembre

Pas de réunion, journée CaSciModOT.

Réunion du jeudi 4 décembre - 10-12h

Participants : U. Razafison, O. Delestre, C. Lucas, M. Rousseau, S. Cordier, F. James, B. Cheviron

"Les termes de frictions dans les modèles d'écoulement"

Discussion avec Bruno Cheviron (qui a fait une bibliographie spécifique à ce sujet).

Exposé

Le cas bicouche (Ulrich)

Mise au point sur le modèle "bicouche" : sous l'effet de la pluie (uniquement), les sillons se remplissent jusqu'à une certaine hauteur puis c'est la couche d'eau supérieure qui évolue (pas d'infiltration). On cherche le modèle équivalent qui ne voit pas les sillons.

Parallèle avec le modèle anisotrope.

GegenSurf : calcul de la sinusoïde inclinée (Carine)

On peut (au moins pour une pente inférieure à 10%) prendre l'approximation du modèle cisaillé.

Courbe avec une pente 0 - Courbe avec une pente de 10% : en superposant les deux courbes, on ne voit pas de différence.
paramètres : les courbes sont en cm, longueur d'onde 80 cm, 10 périodes, de crête à crête 20 cm

Réunion du jeudi 27 novembre - 10-12h

Participants : U. Razafison, O. Delestre, C. Lucas, M. Rousseau, F. Darboux

Manip Etape 0 - la suite ?

A la suite de la série de manip de décembre sur l'étape 0, les données produites seront utilisées avec les différents modèles. Il faudrait décider d'un coordinateur afin que l'exploitation des données puissent fournir des résultats cohérents entre eux.

Termes pluie et infiltration (Marie)

Dans le cas "pluie sans infiltration", la différence entre splitting temporel ou par intégration est très faible (de l'ordre du 10^-5). On suggère à Marie de faire maintenant une comparaison (rapide) dans le cas "pluie avec infiltration" pour éviter une mauvaise surprise ultérieure.

La comparaison avec PSEM_2D sur les parcelles du Niger montre des résultats similaires.

Pour le terme de frottement, on verra la semaine prochaine avec Bruno Cheviron. Il devrait nous apporter des éléments supplémentaires pour les situations avec une très faible hauteur d'eau.

Il apparaît qu'après l'arrêt de la pluie, le flux à l'exutoire met très longtemps à s'annuler. C'est d'autant plus bizarre qu'il y a de l'infiltration. On suggère à Marie de tester le cas d'une bassine (fond plat) avec de l'infiltration. La hauteur d'eau devrait s'annuler rapidement. Selon le résultat, il faudra peut-être rediscuter de la transition sec-mouillé.

Marie va continuer à simuler sur les parcelles du Niger. Il faut demander à Guillaume Nord de faire un cas "Etape 0" avec PSEM_2D.

GegenSurf : calcul de la sinusoïde inclinée (Carine)

Il faut regarder si l'erreur dû au cisaillement est tolérable. Sinon, on pourrait utiliser un calcul par approximation

Rédaction Proceedings (Olivier)

Les graphiques sont bientôt finis. Dans la biblio ne doivent apparaître que les références citées dans le texte.

Olivier enverra une dernière version vendredi matin et soumettra le manuscrit vendredi après-midi

Cas extrêmes pour manip Etape 1 Belgique (Ulrich)

Le premier cas simulé (pente et débit faible) ne montre pas de choc (en fin de compte). On serait donc dans un cas complètement fluvial. Faire le calcul du nombre du Froude en chaque point pour vérifier.

Lancer les 3 autres cas extrêmes. Mettre une hauteur d'eau initiale de qq mm afin de diminuer le temps de calcul pour atteindre l'état stationnaire. On vérifiera visuellement (comparaison des images à différents pas de temps) que l'état stationnaire a été atteint.

Modèle WP1 - Frottement anisotrope (Ulrich)

Il faut travailler avec la sinusoïde définie lors de la dernière rencontre Méthode (vérifier que l'on a les paramètres descriptifs des sillons dans le compte-rendu).

Paramétrage de FullSWOFF (Ulrich)

Pour les paramètres spatiaux (carte topographique numérique, etc.), on utilise un repère direct classique :

y ou j

^
|
|
|
|___> x ou i

Dans le manuel, il faudra indiquer clairement que les calculs sont plus efficaces si les valeurs de z diminuent lorsque x et/ou y augmente (pente vers le N, le NE ou l'E).

FullSWOFF prend en entrée un premier fichier de paramètres. En plus des valeurs numériques définissant la taille du système, le pas de temps, etc., il contient les noms des fichiers spatialisés (en ASCII, avec des lignes "i j ...") pour :

• la topo
• les valeurs pour les frottements
• les conditions initiales (h, u et v)
• les valeurs de l'humidité initiale

Il contient aussi le nom du fichier "sauvegarde des images". Celui-ci est une liste des pas de temps auquels il faut sauver des valeurs spatialisés et les valeurs à sauver. Exemple :

3 huv
25 huv
56 huv
89 h
150 h
800 hu

Il serait bon de faire tourne un optimizer pour vérifier qu'elles sont les parties les plus consommatrices en temps de calcul au lieu d'essayer d'optimiser "au doigt mouillé".

Réunion du jeudi 13 novembre - 10h 12h

Participants : F. James, U. Razafison, O. Delestre, C. Lucas, M. Rousseau, F. Darboux

Fête de la Science (Tous)

Liste du matériel demandé : 2 tables, 2 PC portables + rallonges + 1 prise, 1 support poster.
S'il manque des choses, prévenir Philippe Grillot.

Idée de l'expérience : faire varier la pente et étudier la vague … peut-être pas facile à mettre en place !

Faire un poster pour montrer les apports des maths à l'hydrologie, préparer des vidéos et des images à passer en boucle.

Regarder le résultat, sur une parcelle de M. Esteves, de 2 évènements pluvieux :

• une première pluie où il n'y a que de l'infiltration
• un moment sans pluie
• une seconde pluie où ça commence à ruisseler.

Planning : samedi 22 et dimanche 23, 10h-18h.

• Frédéric et Olivier samedi matin
• ?? samedi aprèm
• Carine (ou ??) dimanche matin
• Carine dimanche aprèm

(+ Frédéric à d'autres moments en cas de besoin)

Problèmes sur l'infiltration (Marie)

Il y a un problème sur le code avec infiltration en mettant le coeff d'infiltration à 0 : l'assèchement n'est pas observé.
Essayer le cas Mac Do. sur le code sans infiltration (suite du test d'Olivier) pour voir ce qu'il se passe.

Tant que ce problème n'est pas résolu, on ne peut pas s'attaquer à la comparaison "infiltration par splitting" - "infiltration comme terme source".

Faux 2d (Ulrich)

Sillons perpendiculaires à la pente (Etape 1) : simulations pour dimensionner les expériences de Belgique.
Pour l'instant, il y a des problèmes.

Il serait nécessaire de mettre en place un module de débuggage qui, en particulier, vérifie la conservation de la masse !

Générateur de topo ("GegenSurf") - (Frédéric)

On revient encore et toujours sur l'équation pour obtenir des sillons pivotés (et pas cisaillés...) à faire (Carine).

Il faudra que GegenSurf renvoie un fichier topo.dat au format (i j z(ij)). Les autres paramètres (Nx,Ny, dx, dy) sont connus a priori (c'est à dire avant l'appel de GegenSurf)

Code 2d (Ulrich principalement)

"Le point" sur les entrées / sorties du code, ce qui se fait en pré-traitement, ce qui est entré par l'utilisateur, ce dont a besoin le code pour tourner ...

Nom proposé : FullSWOF, pour "Full Shallow Water equations for Overland Flow" (vérification faite, ce n'est pas une insulte ;-), c'est même plutôt flatteur (http://www.abbreviations.com/SWOF).

FullSWOF se concentre sur une seule tâche : le simulation des écoulements. Le reste (pré- et post-traitement) est fait par d'autres programmes (tel que GegenSurf). FullSWOF prend en entrée un fichier "parametres.dat" incluant les parametres à utiliser, dont le nom des fichiers d'initialisation de la grille ("topo.dat", "infilt.dat", "friction.dat", "depth.dat", etc.). Tout les fichiers d'initialisation de la grille ont le même format : un fichier ASCII au format (i, j, x(ij)). Il faudra préciser une fois pour toutes le sens du repère dans FullSWOF

Ulrich enverra un exemple de fichiers afin que l'on puisse s'assurer que l'on est bien d'accord sur les formats.

Rédaction (Olivier)

Ajout des graphiques en cours, nouvelle version du proceeding dans les jours à venir.

Réunion du jeudi 6 novembre - 10h 12h

Participants : S. Cordier, U. Razafison, O. Delestre, C. Lucas, M. Rousseau, F. Darboux

Fête de la Science (Tous)

Préparer, à l'échelle du Mapmo un planning pour le week-end du 22-23 pour répartir les membres de methode dans tous les créneaux.

Termes d'infiltration et de pluie (Marie)

On peut prendre en compte les termes d'infiltration par la technique du splitting temporel ou par intégration dans les équations de SV.

Il faut comparer ces deux approches sur un cas test (un plan incliné 2D) et comparer les flux calculés à l'exutoire en régime stationnaire. Le faire pour :

• pluie seule (pas d'infiltration, donc) et pour plusieurs dx
• pluie seule et voir quel est le débit le plus réaliste
• pluie et infiltration.

Si le splitting donne des résultats satisfaisants, on préfèrera cette technique (pour avoir une gestion "indépendante" de l'infiltration)

Avant de comparer avec PSEM_2D, il faut savoir comment la pluie et l'infiltration sont implémentées dans ce modèle (et faire le même choix !).

Cas tests

Placettes Niger (Marie)

Il faut intégrer dans les benchmarks du code C++ les placettes du Niger

Cas tests synthétiques (Ulrich)

1. pente avec sillons perpendiculaires ("faux 2D") avec débit amont imposé. A été réalisée ; cela se remplit et déborde comme prévu
2. pente avec sillons perpendiculaires ("faux 2D") avec pluie
3. vallée + sillons + pluie
4. pente avec sillons perpendiculaires et bruit ("faux 2D") avec débit amont imposé
5. pente avec sillons perpendiculaires et bruit ("faux 2D") avec pluie
6. vallée + sillons + bruit + pluie

Modèle frottement anisotrope (Ulrich)

La formulation du coefficient de frottement a été modifiée pour intégrer la hauteur d'eau stockée dans les sillons.

A rédiger.

Expérimentation Etape 1 (Frédéric)

Frédéric est en contact avec un laboratoire en Belgique qui a un dispositif expérimental correspondant à l'étape 1 en stationnaire : il est possible de mesurer les hauteurs d'eau le long de la pente. Il serait intéressant de disposer de quelques mesures pour des cas contrastés, en particulier avec un régime fluvial et avec un régime torrentiel.

Frédéric va se renseigner sur les conditions expérimentales (largeur du système, gamme de pente, gamme de débit). Il demandera à Olivier de réaliser des simulations numériques afin de préciser les conditions expérimentales pertinentes (dimensionnement des cas tests).

Générateur de surface (Frédéric)

Pour l'étape 2, il sera important de disposer d'un outil commun pour générer des surfaces (calculs paramétriques en particulier). Ce générateur devra être capable de produire des grilles de hauteur à mailles carrées avec :

• une pente générale (valeur de la pente et direction) ;
• des sillons sinusoïdaux (longueur d'onde, amplitude, phase = décalage latéral). Il faut que cela soit une rotation de la sinusoïde et pas un cisaillement ;
• un bruit simulant la présence des mottes. Des caractérisations de ce type de bruit ont déjà été faites en science du sol (pente du spectre, fréquence de coupure).

Frédéric avait mis au point un générateur de surface au cours de sa thèse. Il faudra le perfectionner un peu avant de le mettre à la disposition des membres du projet Méthode. Stéphane fera une interface web.

Modèle "bi-couche" (limitation de hauteur dans SV)

A écrire (Stéphane, Ulrich, Carine)

Réunion du jeudi 16 octobre. 9h30-11h

Participants : S. Cordier, F. James, U. Razafison, O. Delestre, C. Lucas, M. Rousseau, F. Darboux, A. Elhajj

Developpement du code C++ Saint-Venant2D (Marie, Ulrich, Olivier)

Architecture

• L'architecture actuelle du code, sous forme de librairies, a été présentée (Ulrich).
• Discussion sur les modifications et les améliorations futures à apporter pour les utilisateurs.

Gestion des versions

• Le code est désormais sous svn (Nicolas).

Les entrées

• Le codage des entrées par un parser ASCII est mis en place (Ulrich).

Réunion du jeudi 9 octobre. 10h-12h (réunion d'équipe, interne)

Participants : S. Cordier, F. James (- - -), U. Razafison, O. Delestre, C. Lucas, M. Rousseau, F. Darboux, A. Elhajj

Avancement du code 2D en C++ (Marion, Ulrich)

Débogage

• Les problèmes rencontrés sont résolus (cause = seuils fixés différemment sur le code fortran).
• La comparaison avec les solutions de référence est en cours. On refera les ruptures de barrage oblique sur fonds sec et mouillé.

Ajout pluie et infiltration

• Après le codage en C++ de la pluie et de l'infiltration, refaire les tests sur la parcelle du Niger.
• Implémenter aussi la pluie en non-splitté.

Méthode et ordre

• La modification est faite pour avoir la même méthode pour les 2 ordres.

Entrées (Ulrich)

• Le codage des entrées et des sorties en XML imposerait l'installation d'un logiciel spécifique. On utilisera donc un codage ASCII.
• Nécessité de discuter de l'organisation des fichiers.
• Le parser doit être mis en place.

Gestion des versions

• A mettre en place (Nicolas)

Documentation

• On cherche toujours à pouvoir mettre de la biblio dans Doxygen...

2D Tôle ondulée (Ulrich)

• Tester les cas de J. Peugeot (2007)
• Tester le McDonald périodique

Note : on appelera "Etape 3" une parcelle (homogène) et "Etape 4" un bassin versant (ensemble de parcelles).

Thèse Olivier

• La rédaction du chapitre sur les frottements est en cours.
• Il prendra la forme d'un proceedings d'HYP2008. Date limite = 30/11/08

Frottement et anisotropie

• Aucun élément nouveau n'a été trouvé dans le livre "Geophysical Fluid Dynamics".
• Il faut contacter Hervouet pour savoir d'où viennent ses coefficients de frottements dans Navier-Stokes. Est-ce une analogie par rapport à Saint-Venant ou est-ce référencé ailleurs ? (Carine, fait le 13/10 : réponse à lire)
• Stéphane Cordier fait au tableau le développement des équations de Saint-Venant pour un frottement anisotrope. (Retranscription par Ulrich)

Réunion du jeudi 2 octobre. 10h-12h

Le point sur le code 2D en C++

Problème aux transitions sec-mouillé (qui n'existait pas avec le code fortran), en particulier sur le calcul de la vitesse

• essayer sans l'optimisation du code à la compilation
• regarder sur un seul pas de temps, avec faible résolution
• regarder le calcul des flux

(--> Marie Rousseau)

Entrées

Elles sont maintenant centralisées dans la classe mère.La prochaine étape est de lire des fichiers de paramétrage. Il faut pour cela inclure un parser. Ulrich propose que les fichiers d'entrées soient en XML. Il lui est demandé de montre à quoi ressemblerait un fichier d'entrée XML. Plus-value par rapport à de l'ASCII "de base" ?

(--> Ulrich)

Documentation

Il faut trouver comment inclure les références bibliographiques dans Doxygen

Gestion des version par CVS

Nicolas, le responsable informatique du Mapmo, va mettre en place un serveur CVS.

Sorties graphiques

Il serait possible d'inclure des routines openGL dans le code C++ pour gérer l'affichage des sorties graphiques.

Coefficients de frottement

Les coefficients de Darcy-Weisbach et de Manning se différencie essentiellement par la valeur de l'exposant sur la hauteur d'eau. D'autres valeurs de l'exposant pourrait être plus appropriées selon le type d'écoulement.

Stages 2009

• modélisation, aspects théoriques (profil mathématicien)
• modélisation, aspects numériques (profil méca flu)
• création d'une interface graphique (profil informaticien)

Fête de la science

On se décide pour deux gouttières (l'une avec un écoulement critique, l'autre sous-critique) arrivant dans une bassine ("lac"). L'INRA s'occupe de la partie matérielle, le Mapmo de la partie numérique (prévision du choc).

Réunion du jeudi 25 septembre. 10h-12h

• Questionnements sur la forme du frottement de fond.
• Passage de Navier-Stokes à Saint-Venant avec les formulations de Chézy et Manning (Olivier, Ulrich, Carine).
  De Navier-Stokes à Saint-Venant avec les formules de Chézy et Manning. (Carine)
• Réflexions sur les tests numériques à faire.

Réunion du jeudi 18 septembre. 10h-12h

• Passage de Navier-Stokes à Saint-Venant : les différentes étapes et choix d'ordres de grandeur.
• Cas avec fond plat : De Navier-Stokes à Saint-Venant. (Carine)
• Cas plus complexe avec fond variable : calculs.