open TELEMAC-MASCARET The mathematically superior suite of solvers

Bonjour,

J'utilise la nouvelle version sortie recemment (6.1) de Telemac2D, et j'aimerais lancer une simulation que j'ai faite avec la version precedente (6.0) avec quelques ameliorations si possible.

Cependant, j'ai rencontre un probleme dans la definition d'un traceur: alors que dans la version precedente, je pouvais definir la concentration d'un traceur sur une ou plusieurs frontiere(s), ici j'ai essaye la meme chose mais le resultat que j'obtiens est une concentration nulle en tout point et a tout instant. J'ai vu l'apparition d'un nouveau mot-cle, VALEURS INITIALES DES TRACEURS, et j'ai tente d'entrer des valeurs mais la je me suis retrouve soit encore avec rien du tout, soit avec une concentration constante dans le temps et uniformement repartie.

La seule solution que j'ai trouvee est de placer une source au niveau de la frontiere, mais j'aimerais savoir s'il est toujours possible de definir une concentration en traceur a la frontiere comme precedemment ou s'il y a eu un changement de "philosophie" si je puis dire, et que desormais le seul moyen est d'utiliser une ou plusieurs source(s) ponctuelle(s).

Ci-joint mon fichier cas, parametre avec les sources.

Merci d'avance !
Remi Cousson

Desole, je n'avais pas modifie mon fichier en .txt. C'est maintenant chose faite !

Hello,

No, what you do should indeed still work. The initial values of tracers are only used for initial conditions in the whole domain.

In your case the keyword to be used is :

VALEURS IMPOSEES DES TRACEURS (or PRESCRIBED TRACERS VALUES)

The only case where this is not taken into account is when the flow is exiting your boundary. Of course the tracer boundary condition must be prescribed value, i.e. 5.

Check this and if it does not work send your whole case in a zip file.

With best regards,

Jean-Michel Hervouet

J'ai verifie, et normalement ca devrait mercher, mais ce n'est pas le cas. Le flux est a priori bien dirige car lorsque j'utilise les sources, le traceur se dirige bien dans le bon sens.

Je joins donc le dossier entier, si vous arrivez a trouver de quoi il retourne.

Merci

Hello,

Actually your boundary conditions file gives the tracer as free on all liquid boundaries. When a line in your file cas.conlim starts with 4 5 5 (prescribed discharge and free elevation) the next integer on the line should be 5 and not 4. If you do this you will see tracer entering the domain, the problem was thus not due to the new version (or did it really work like this with version 6.0 ???)

You should change also in your file "cas" :

CORIOLIS : true into CORIOLIS: OUI

(true is not understood with my compiler)

add:

PRECISION POUR LA DIFFUSION DES TRACEURS : 1.E-10

sometimes if the accuracy is not enough, there will be no more iterations of solver and the tracer will stop.

remove

NOMBRE DE SOUS-ITERATIONS : 10

this is very costful and useless here. If problem decrease the time step.

SOLVEUR POUR LA DIFFUSION DES TRACEURS : 7

This is the GMRES solver, put 1 instead (conjugate gradient)

You can also remove the line :

PARAMETER FILE : cas (it is useless and drives you mad when you change the name and forget to change this line)

I hope it will be better like this, with best regards,

Jean-Michel Hervouet

Merci, ca marche parfaitement a present !
Mais je pense que c'etait la meme chose sur la version 6.0, en tout cas je n'ai touche que rarement au fichier de conditions aux limites, je passais plutot par Fudaa, je ne peux pas etre sur. Mais pour ce qui est des valeurs des mots cles sur Fudaa, c'etait a priori les meme que sur la version 6.0.

Encore merci
Cordialement
Remi Cousson

Hi Jean Michel
I downloaded this example and made the modifications that you recommended but I can't seem to get the tracer to enter the flow. Please could you tell me what I should do to get it to work.
Merci et Joyeux Noel

Steve

Hello,

It seems to me that the problem lies in the boundary conditions file. Lines where your velocities are prescribed (starting with 4 5 5) have a free tracer (a 4 further in the lines) whereas it should be a 5.

Hope this helps, and happy X-mas,

Jean-Michel Hervouet

Thanks, Jean-Michel
I've got it to work now