open TELEMAC-MASCARET The mathematically superior suite of solvers

Bonjour,

je construis un modèle hydraulique fluviale avec TELEMAC 2D et je n'arrive pas à faire comprendre au logiciel d'utiliser la courbe de tarage comme étant ma condition à la limite avale (exutoire), frontière liquide numéro 2.

J'ai suivi le manuel d'utilisateur:
- mot clé COURBE DE TARAGE réglé sur Q(Z)
- mot clé FICHIER DES COURBES DE TARAGE indiqué mon fichier courbe de tarage
- courbe de tarage construite comme indiqué dans le manuel.
- le Q et Z initial sont dans la gamme de la courbe de tarage

J'ai essayé de définir la frontière liquide comme etant une hauteur imposée, un débit imposé etc., rien n'y fait, la courbe de tarage n'est jamais prise en compte ...

Quand je lance la simulation, la courbe de tarage est bien lue mais pas prise en compte.

Où est mon erreur?
Comment dire a TELEMAC, "utilises la courbe de tarage sur la frontiere liquide numéro 2"?

Merci

JA

A response will be given to you shortly -- can you use English in the futur ?

sure, I will write in english in the future, sorry.
I saw many threads written in French, that's why I used French.

Feel free to reply in English

I have checked my notes and what you wrote -- you seem to have everything I would have though of.

Would it be possible to send us you files (GEO, CAS, CONLIM, etc.) so we can reproduce your test conditions ?

Thank you.

Hello,

The solution is simple (and disapointing), actually the functionality Q(Z) is not yet implemented, only Z(Q) works so far (see subroutine BORD.F and test on STA_DIS_CURVES). A general Q(Z) is very tricky to implement, and would require a choice of a point where to take the elevation. If you really want Q(Z), I would suggest that you implement the function Q.F and put it in your Fortran file. You would have to choose a point where to take the elevation, preferably a point that is slightly upstream the boundary itself, otherwise it may give unstable results. Then having this elevation you deduce the corresponding Q, and put it as a result of the function.

Note that function Q will be called only if the key-word : DEBITS IMPOSES (or PRESCRIBED FLOWRATES) has been used in the parameter file. So you must give an a priori discharge to the corresponding boundary.

Last bad news: in parallel mode the point you choose for taking the elevation will not be in all sub-domains, thus an extra implementation will be required to overcome this, I can send an example if required.

With best regards,

Jean-Michel Hervouet